Medical Conditions in the First Years of Life Associated with Future Diagnosis of ASD in Children

This study examines medical conditions diagnosed prior to the diagnosis of autism spectrum disorder (ASD). Using a matched case control design with 3911 ASD cases and 38,609 controls, we found that 38 out of 79 medical conditions were associated with increased ASD risk. Developmental delay, mental health, and neurology conditions had the strongest associations (ORs 2.0–23.3). Moderately strong associations were observed for nutrition, genetic, ear nose and throat, and sleep conditions (ORs 2.1–3.2). Using machine learning methods, we clustered children based on their medical conditions prior to ASD diagnosis and demonstrated ASD risk stratification. Our findings provide new evidence indicating that children with ASD have a disproportionate burden of certain medical conditions preceding ASD diagnosis. Electronic supplementary material The online version of this article (doi:10.1007/s10803-017-3130-4) contains supplementary material, which is available to authorized users.


Introduction
Autism spectrum disorders (ASD) are neurodevelopmental disorders defined by impaired social interaction and repetitive behaviors. The most recent estimates of ASD prevalence in the United States range from 1 to 2% (Kogan et al. 2009;CDC 2009;Blumberget al. 2013). Experienced clinicians can diagnose ASD in children as young as age 2 (Charman and Baird 2002). However, the average age of ASD diagnosis is >3 years in the general population and even older in minority populations, suggesting that many children are not diagnosed as quickly as possible (Mandell et al. 2005). Earlier diagnosis of ASD is critical because earlier treatment can reduce the degree of impairment and improve function (Lord et al. 2005). A possible first step toward earlier ASD diagnosis is risk stratification of children into groups who have lower and higher risks of ASD so that higher risk children can be monitored more closely. Clinical factors that precede an ASD diagnosis and are strongly associated with risk of ASD could be immensely useful for risk stratification.
Many children with ASD have also been diagnosed with co-occurring medical conditions during their later 1 3 childhood and adulthood. The medical conditions most often reported to co-occur in children with ASD include immune conditions (Comi et al. 1999;Croen et al. 2005;Rosen et al. 2006;Sweeten et al. 2003;Vargas et al. 2005), seizure (Bauman 2010;Hara 2007;Rapin 1996;Tuchman et al. 1991;Volkmar and Nelson 1990), sleep disorders (Bauman 2010;Ferber 1996;Goldman et al. 2011;Johnson et al. 2009;Krakowiak et al. 2008;Richdale 1999), gastrointestinal disorders (Bauman 2010;Buie et al. 2010), metabolic disorders (Bauman 2010), hormone dysfunction (Bauman 2010), mitochondrial disorders (Oliveira et al. 2005;Rossignol and Frye 2012), nutritional deficits (Bauman 2010), and mood disorders (Bauman 2010). Most of these studies have focused on medical conditions that occur after the diagnosis of ASD. Much less is known about medical conditions that occur in the time prior to ASD diagnosis, i.e. the pre-diagnostic period. Furthermore, previous studies examining the time prior to ASD diagnosis rely on parent reports of child histories collected after ASD diagnosis that could be influenced by recall bias. Large-scale epidemiologic studies with prospectively collected medical data are needed to quantify meaningful and reliable estimates of the burden of co-occurring medical conditions among children with ASD during this time period before ASD diagnosis.
This study considered the diagnoses of medical conditions in the period preceding the first ASD diagnosis to achieve two objectives. First, we sought to determine whether certain medical conditions occurring in early childhood were associated with a higher risk of future ASD diagnosis. Second, we used supervised clustering methods developed for machine learning to identify risk-stratified clusters of children based on their documented medical conditions in early childhood before ASD diagnosis.

Study Population
The study population was drawn from the population of all children who were Kaiser Permanente (KP) members in Northern California, Georgia, and the Northwest born from January 1, 2000 through December 31, 2009. In order to capture all early medical conditions, children included in the study were required to have at least 9 months of KP membership in both their first and second years of life. Additionally, all children must have utilized health services at least once after their second birthday. We collected data from electronic medical records (EMR) and administrative claims for these children up to June 2012. KP Northern California (KPNC) is the largest and oldest group-model, pre-paid, integrated health care delivery organization in the United States, providing health care to over 3.8 million people; KP Northwest (KPNW) serves over 500,000 members in Oregon and Washington; KP Georgia (KPGA) has over 250,000 members. Table 1 shows the demographic characteristics of the study population.

Case and Control Definition
ASD cases were defined as children with either (a) an ASD diagnosis from an ASD specialist or (b) two or more ASD diagnoses from non-specialists, separated in time. Specialists were doctors with known expertise in pediatric ASD (many of whom were psychologists or developmental pediatricians); non-specialists were the remaining doctors. These criteria were based on an ASD case validation exercise completed across all sites (Coleman et al. 2015). Controls were required to never have had an ASD diagnosis as of June 2012. Controls were matched to cases at a 10:1 ratio on site, year and month of birth, and length of membership in the 12 months before and after the case's first ASD diagnosis.

Medical Conditions
We restricted our analysis to medical conditions recorded in the EMR during the "pre-diagnostic period," defined for ASD cases as the time period prior to the first ASD diagnosis, and for matched controls as the time period prior to the age of first ASD diagnosis of the matched ASD case. Medical conditions were identified by ICD-9 codes. Over 1000 ICD-9 codes were grouped into 79 medical conditions (e.g., constipation) within 19 domains (e.g., gastrointestinal). First, based on clinician input, each condition was required to have been recorded a minimum number of times in the EMR to be considered valid (Appendix A: ESM Online). Preliminary grouping was guided by previous studies of ASD-related conditions and independently reviewed and refined by two clinicians. Additional ICD-9 codes related to each condition were identified using phenotype-wide association study methodology (Doshi-Velez et al. 2014) and reviewed by one clinician; codes deemed relevant were included for that condition. Finally, rare conditions (prevalence <0.5% among ASD cases) were regrouped into the "other" category within each domain.

Statistical Analyses
We tabulated the count and observed prevalence of each condition among cases and among controls during the prediagnostic period. Conditional logistic regression analyses were carried out to estimate the odds ratio (OR) while accounting for the matched case-control design. The OR can be interpreted as a relative risk of ASD diagnosis in this study, because ASD meets the rare outcome criterion (Jewell 2003). Adjusted models accounted for the matching by site and age and adjusting for sex, maternal race, maternal education, and household income. Wald-based 95% confidence intervals (CI) were constructed for each OR using a level of 0.05. To account for multiple comparisons, we used the Bonferroni cutoff, which considers only p-values <0.00063 to be statistically significant based on testing 79 conditions.
To identify clusters of multiple conditions that were strongly associated with subsequent ASD diagnosis, we used the conditional inference tree (CIT) method (Hothorn et al. 2006). CIT is a supervised clustering method that creates clusters defined by specific combinations of conditions and detects interactions between conditions without having to pre-specify those interactions. As a supervised algorithm, the splits for each cluster are based on differentiating ASD cases from controls, resulting in a set of clusters that should distinguish condition patterns of varying association with subsequent diagnosis of ASD. To account for multiple comparisons and avoid overfitting, we used the Bonferroni cutoff for statistical significance as the stopping criteria. CIT was implemented using the party package in R (R Core Team 2014). The categorical variable created by CIT was then used as a covariate in our conditional logistic regression model, and we computed crude and adjusted ORs using the lowest risk condition cluster as the reference.
We conducted two clustering analyses of medical conditions. The first clustering analysis identified co-occurring medical condition clusters among medical conditions documented during the pre-diagnostic period. This analysis was aimed at understanding the overall burden of multiple medical conditions that occur prior to ASD diagnosis. The second clustering analysis focused specifically on medical conditions present before age 3. By focusing on conditions that occur in early life, this analysis was aimed at identifying comorbidity clusters that might have potential to improve early ASD risk stratification.

Results
Demographic differences between cases and controls are presented in Table 1. We observed a large and significant difference in the gender of ASD cases and controls: 82% of ASD cases were males, compared with 51% of the control group. There were also small but statistically significant differences in the distributions of race/ethnicity, median income, and education. A higher percentage of ASD cases were white or Asian and a lower percentage of ASD cases were Hispanic or other compared to the control group. ASD cases also lived in areas with slightly higher median incomes and slightly higher levels of education compared to the control group. There were no statistically significant differences in the matching variables age and site. Among ASD cases, the average age of ASD diagnosis was 3.99 years. Table 2 lists each medical condition, its prevalence among cases and among controls during the pre-diagnostic period, and the estimated odds ratio (OR) for the crude and adjusted associations with ASD risk. We found that 38 of the 79 medical conditions had statistically significant associations with ASD risk after multiple testing adjustment by the Bonferroni cutoff. Within the developmental delay domain, language delays were the most frequently diagnosed and most strongly associated with ASD risk; learning & cognitive disorders, global delays, motor delays also had strong associations. All mental health conditions were statistically significant; disruptive impulse conduct disorders, attentiondeficit/hyperactivity disorders (ADHD) and anxiety disorders were the most prevalent among ASD cases during the pre-diagnostic period and had stronger associations with ASD risk than depression, adjustment disorders, tic disorder and other mental health conditions, each of which had low prevalence among ASD cases (<1%) during the pre-diagnostic period. Most neurologic conditions showed strong associations with risk of ASD. Cerebral palsy, epilepsy & recurrent seizures, and other disorders of the central nervous system had the strongest associations. Weaker but statistically significant associations were observed for other congenital anomalies of the nervous system, abnormalities of skull and face, hereditary and degenerative diseases of the central nervous system, and paroxysmal disorders.
Conditions in the domains of nutrition; genetics; sleep; and ear, nose & throat demonstrated moderate increases in ASD risk. Statistically significant increases in ASD risk of low to moderate magnitude were observed for conditions in the domains of ophthalmology, pulmonary, metabolic, musculoskeletal, and gastrointestinal (GI) conditions. Conditions in the domains of autoimmune, infections, hematology, allergy, injury, cardiovascular, endocrine, and asthma had weak associations, most of which were not strong enough to reach our Bonferroni cutoff for statistical significance. Table 3 shows the results from our first clustering analysis with condition clusters selected by conditional inference tree regression. The reference cluster (children without developmental delay; mental health; ear, nose & throat; nutrition; or ophthalmology) has the lowest risk of ASD diagnosis and includes 66.81% of the controls and 19.66% of the ASD cases. The condition clusters are listed from lowest to highest OR in comparison to the reference group, illustrating the increasing relative risk of ASD diagnosis. Note that the clusters are mutually exclusive such that every child belongs to only one cluster, in contrast to the results displayed in Table 2 where conditions are studied individually.

Associations of Condition Clusters Using Supervised Clustering
Developmental delay and mental health condition clusters were associated with the highest relative risks of ASD. Children with language delay were categorized into four distinct condition clusters; the relative risks for ASD were very high in each cluster but also showed substantial variability (ORs ranged from 28 to 101). The cluster of language delay with global delay had the highest relative risk. Among children without developmental delays, three mental health condition clusters emerged with 9-38-fold elevated risks of ASD compared to children in the lowest risk category. The disruptive impulse conduct disorder cluster had the highest relative risk for ASD and the next highest relative risk cluster included children with ADHD but without disruptive impulse conduct disorder. Beyond the mental health and developmental delay domains, condition clusters with moderate to low increases in ASD risk were defined by combinations of ophthalmology; ear, nose & throat; and nutrition conditions. Table 4 displays the results from our second clustering analysis focused on early ASD predictors by further restricting the pre-diagnostic period conditions to those identified prior to age 3. The reference cluster (children without developmental delay; ear, nose & throat; nutrition; ophthalmology; musculoskeletal; or constipation) has the lowest risk of ASD diagnosis and includes 65.31% of the controls and 29.30% of the ASD cases. Compared to the reference cluster, developmental delay condition clusters have the strongest associations with ASD risk. The three clusters with the highest relative risks of ASD all included language delay, with substantial heterogeneity in the relative risk depending on the co-occurrence of global delay or other disorders of the central nervous system. The mental health conditions do not appear in any cluster definitions. Outside of the developmental delay domain, the conditions of seizures; musculoskeletal; ear, nose & throat; nutrition; ophthalmology and gastrointestinal conditions defined   eight condition clusters, each associated with low to moderate increases in ASD risk. Latency time in months between diagnosis of all conditions in the cluster and subsequent diagnosis of ASD is reported in the last column of Tables 3 and 4, computed only among ASD cases. When a cluster includes multiple conditions diagnosed on different dates, the diagnosis date for the last condition identified is used, since diagnosis of all conditions are required to be identified in that cluster.
We found that the mean latency time was greater than 12 months for every cluster identified in Tables 3 and 4, with mean latency time as high as 50 months. In general, condition clusters involving developmental delays or mental health conditions had shorter latency times than other condition clusters; however, the cluster of children diagnosed with a language delay and a mental health condition who had no diagnosed global delay had the longest mean latency time in Table 3 (mean = 32.7 months, SD = 25.5). Table 4 Associations of early life child medical condition clusters with risk of ASD up to age 12, comparing medical conditions diagnosed before age 3 and before the first ASD diagnosis to medical conditions diagnosed during the same time period among matched controls a Condition clusters selected by conditional inference tree regression identifying combinations of conditions along the main branches (in bold) of whether a developmental delay had been diagnosed b The adjusted conditional logistic regression model includes matching by site and age and adjusting by sex, maternal race, and census tract levels of education and household income c Latency time between having all conditions in cluster and being diagnosed with ASD among ASD cases

Discussion
By leveraging comprehensive prospectively recorded EMR data on a large and diverse population of children in an integrated health care delivery system, this study generated robust estimates of the burden of childhood medical conditions diagnosed early in life among children subsequently diagnosed with ASD. Our manuscript presents three key findings. First, our study provides new evidence indicating that children who are diagnosed with ASD experience a higher burden of 38 medical conditions in the years preceding ASD diagnosis, including GI conditions, developmental delays, mental health disorders, musculoskeletal abnormalities, neurological conditions, and sleep disorders. Second, our highly-powered study found no evidence that children who are diagnosed with ASD experience a higher burden of the other 41 medical conditions in the years preceding ASD diagnosis, even though the conditions were all hypothesized to be more common in children with ASD based on previous studies. Third, we demonstrated how combinations of medical comorbidities could aid in risk stratification for ASD prior to ASD diagnosis using a supervised clustering analysis based on machine-learning methods. We identified 14 risk clusters such that ASD risk increases across these strata, with ORs ranging from 1.5 to 100.

Medical Comorbidities Documented Before ASD Diagnosis
The associations of 38 medical conditions with future ASD diagnosis reported here strengthen the evidence presented in previous studies of the higher prevalence of gastrointestinal conditions, developmental delays, mental health disorders, musculoskeletal abnormalities, neurological conditions, and sleep disorders among children with ASD. In contrast to the few previous studies that relied on parent reports of medical conditions in the years before ASD diagnosis-which is highly subject to recall bias-the key innovation of our study is that we relied on diagnoses made by clinicians and recorded in the child's medical record prospectively before the ASD diagnosis was made. Our findings provide new evidence indicating that children with ASD have a disproportionate burden of certain medical conditions in the years preceding ASD diagnosis. A high prevalence of GI conditions among children with ASD has been reported in several studies with no comparison to control children (Molloy and Manning-Courtney 2003;Nikolov et al. 2009;Ming et al. 2008). However, a recent consensus report evaluating the diagnosis and treatment of GI conditions among individuals with ASD concluded that: "The prevalence of gastrointestinal abnormalities in individuals with ASDs is incompletely understood," and went on to highlight the methodological limitations of previous studies and the possible influence of referral bias (Buie et al. 2010). For example, one recent study using descriptive reports from structured interviews estimated the prevalence of GI symptoms among children with ASD to be 70% compared to a prevalence of only 28% among control children (Valicenti-McDermott et al. 2006). In stark contrast, a small retrospective case-control study of the medical records of 96 ASD children and 449 control children reported that children with ASD did not have higher odds of a history of GI disorders before ASD diagnosis compared to children without ASD (OR 1.0, 95% CI 0.5-2.2) (Black et al. 2002). This study is the most similar in design to ours, although much smaller, and the CI bounds are consistent with our results (adjusted OR estimates ranging from 1.0 to 2.2 across GI conditions). Overall, our findings support our hypothesis that children with ASD are more likely than control children to have a documented comorbid GI condition in the time before ASD diagnosis, but the difference is quite modest, with effects smaller than we expected. Our findings suggest that descriptive reporting of gastrointestinal conditions may substantially over-estimate the difference in prevalence between children with ASD and controls.
We found that developmental delays and mental health conditions showed the strongest associations with future ASD diagnosis. The strong association between language delay and ASD risk is not surprising given that language delay is a component that characterizes ASD diagnosis (Rapin 1991;American Psychiatric Association 2013). Notably, we found a moderately strong association of motor delay with future ASD diagnosis. Motor delays are not considered prodromal to ASD development, but some studies have reported increased motor impairment among children already diagnosed with ASD as compared to control children (Noterdaeme et al. 2002;Mari et al. 2003). To our knowledge, our finding of an association between ASD and a motor delay diagnosis that precedes ASD diagnosis has not previously been reported.
For 41 other medical conditions, we did not find statistically significant evidence of an association with ASD risk in this highly-powered study. We found that children who are later diagnosed with ASD were not more likely than controls to be diagnosed with allergy, asthma, infections, autoimmune conditions, hematology, genitourinary, accidental injury, endocrine disorders or cardiovascular disorders in the years preceding ASD diagnosis. In contrast, some previous studies have reported that ASD is associated with comorbid medical conditions of allergies (Yaghmaie et al. 2013;Theije et al. 2014), asthma (Kotey et al. 2014), infections (Atladottir et al. 2010), autoimmune conditions (Gesundheit et al. 2013), and endocrine disorders (Chen et al. 2009).
We previously reported that both allergies and autoimmune diseases were diagnosed more often among children already diagnosed with autism than among controls (Zerbo et al. 2015). One possible explanation for why results during the time period prior to ASD diagnosis differ from results when studying children already diagnosed with ASD is that some medical conditions may have a later age of onset than the age of ASD diagnosis. Alternatively, increased health care utilization among children already diagnosed with ASD could lead to increased opportunity for diagnoses of other medical conditions relative to children without ASD who have fewer encounters with the health care system. Further research is needed to elucidate how medical conditions that co-occur with ASD may differ before and after ASD diagnosis and the reasons driving those differences.

ASD Risk Stratification Based on Medical Comorbidities
Our study also demonstrated the ability of medical condition clusters to differentiate relative risks of ASD. We identified 14 risk strata and illustrated a gradient of increasing ASD risk across these strata, with ORs ranging from 1.5 to 100 in comparison to the lowest risk group of children. Risk stratification of children into groups who have lower and higher risks of ASD could be a first step toward earlier ASD diagnosis by identifying groups of children for whom increased monitoring or earlier screening protocols could be developed. Earlier diagnosis of ASD is critical because earlier treatment can reduce the degree of impairment and improve function (Lord et al. 2005). Our analysis of the latency time between diagnosis of all conditions in each cluster and subsequent diagnosis of ASD found that the mean latency time was greater than 12 months for every cluster, and as high as 50 months for some, illustrating that there is a substantial window of time between identification of higher risk children and the diagnosis of ASD. Not surprisingly, the high risk clusters were dominated by developmental delays and mental health conditions. However, there were interesting features regarding how these delays clustered together in different combinations. Interestingly, both cluster analyses identified several distinct clusters within the high risk group of children with language delays, where language delays in combination with certain other comorbid conditions showed heterogeneous relative risks for ASD that were all higher than the relative risk in the cluster of language delay alone. These results underscore the importance of considering combinations of conditions. Global delay sometimes incorporates language delay because global developmental delay is defined as significant delay in two or more areas of development (American Psychiatric Association 2013). Children diagnosed with a global delay who do not have a separate language delay diagnosis may have a language delay as a component. Thus, it is difficult to interpret the higher relative risk for the cluster of global delay with language delay compared to the cluster of global delay without language delay. Our findings of an association between motor delay and ASD are consistent with other recent studies suggesting an increased prevalence of motor impairment among those with ASD as compared to controls (Ming et al. 2007). Notably, our cluster analyses showed an association of motor delay with ASD risk even among children who had no co-occurring language delay diagnosis.
Mental health conditions in the pre-diagnostic period showed very strong associations with subsequent ASD diagnosis in the main analysis. The all-age cluster analysis showed that mental health diagnoses were important in both the absence and the presence of developmental delays. These mental health condition clusters were not identified in the age 3 cluster analysis which was restricted to medical conditions identified before age 3. This finding is not surprising given that typical age-of-onset is 4-11 years for ADHD and 5-15 years for conduct disorders (Kessler et al. 2007). Increased monitoring or screening of children with mental health conditions may be important considering the evidence of (i) the strength and timing of these associations with mental health conditions, (ii) more than 10% of ASD cases were diagnosed with a mental health condition prior to their first ASD diagnosis, and (iii) the cluster of children diagnosed with a language delay and a mental health condition who had no diagnosed global delay had the longest mean latency time in Table 3 (mean = 32.7 months, SD = 25.5). These findings are supported by a recent study of children with ASD, ages 2-17, in which approximately 20% of ASD cases were first diagnosed with ADHD, and on average those children were not diagnosed with ASD until 3 years after their ADHD diagnosis (Miodovnik et al. 2015).
Outside of the domains of mental health and developmental delay, the conditions of ophthalmology; ear, nose and throat; and nutrition defined other condition clusters. Individually, each condition was also associated with a moderate increase in ASD risk. The conditions of constipation, gastroesophageal reflux disease, musculoskeletal abnormalities, and seizures contributed to defining the clusters only when studying conditions diagnosed before age 3. Several conditions reached Bonferroni significance in individual analyses but not cluster analyses, including metabolic, gastrointestinal, genetic, pulmonary, sleep, and neurologic conditions; these conditions tended to have either low prevalence or a weak association with ASD risk (OR < 2).
Several recent studies of ASD have utilized a clustering approach to identify risk factors, but the focus of those studies differs from ours. Doshi-Velez et al. only examined children diagnosed with ASD and sought to cluster underlying patterns that may represent heterogeneous etiologic pathways for ASD (Doshi-Velez et al. 2014). Brennan et al. focused on symptoms identified in screening protocols and used hierarchical cluster analyses to detect underlying clinically distinct subgroups (Brennan et al. 2015). Chawarska et al. studied younger siblings of children with ASD, a high risk population, and used a tree-based method to identify clusters of behavioral features associated with subsequent ASD diagnosis or atypical development (Chawarska et al. 2014). In contrast, our study included children with ASD and age-matched controls and identified condition clusters that indicated increased risk of ASD in the pre-diagnostic period.

Strengths and Limitations
Our study has a number of strengths and limitations. Our large and diverse study population across three health systems in different parts of the U.S. contributes to the generalizability of these results. Physician documentation of conditions was recorded prospectively across a broad range of medical conditions prior to ASD diagnosis, preventing recall bias and diagnostic bias. However, the ascertainment of medical diagnoses relies on the documentation in each record, so differential utilization of health services that leads to a differential reporting of conditions by case-control status could bias our results. In addition, there is potential for some misclassification in ASD phenotypic assessment since not all diagnoses were made by experts.
A statistical strength of this study is the Bonferroni control for multiple testing in individual and cluster analyses. A concern in tree-based supervised clustering is model overfitting by creating too many clusters that are not meaningful and cannot be replicated (Strobl et al. 2009). Pruning, a cross-validation procedure, can combat overfitting, but a strength of our study is using a Bonferroni test for the stopping criteria which provides a strong degree of control for overfitting without loss of predictive performance (Hothorn et al. 2006). The main statistical limitation of our analyses is that the matched case-control design prohibits us from calculating an estimate and confidence interval for the absolute risk, thus limiting our statistical inference to ORs (Agresti and Kateri 2011). We cannot directly evaluate the predictive characteristics of this clustering model in terms of the absolute risk of ASD for each cluster; this is a topic for future work. This future work could extend to quantifying the sensitivity and specificity of different screening protocols that might be targeted to particular high risk clusters. Although we have demonstrated risk stratification with a gradient of increasing relative risk across clusters, the clinical utility of integrating the risk stratification with screening protocols would depend on a number of factors including the strength of the association of each risk cluster with ASD and the proportions of cases and controls in each risk cluster.

Conclusion
By taking advantage of the comprehensive EMR data available for a large and diverse population of children in an integrated health care delivery system, this study has generated robust estimates of the burden of medical conditions among children with ASD during the pre-diagnostic period. This work can increase awareness of the diversity of conditions affecting children with ASD prior to their first ASD diagnosis. Further study of conditions during the prediagnostic period may increase our ability to identify children at risk of having ASD, and improved risk stratification may lead to earlier diagnoses and interventions.
Funding This study was funded by Autism Speaks Grant# 7788.
Author Contributions LAC conceptualized and designed the study, supervised all aspects of the study, coordinated data collection across sites; SEA designed the statistical analysis plan, interpreted the results, and led the writing of the manuscript; VY participated in the design of the study, performed data analyses, and helped to draft the manuscript; YQ led the data acquisition, integrated the data across sites, and performed data analyses; MD participated in the design of the study, and contributed to the interpretation of the results; FL participated in the design of the study, supervised and coordinated data collection at one site, and contributed to the interpretation of the results; PC contributed to the data acquisition and performed data analyses; RD participated in the design of the study, supervised and coordinated data collection at one site, and contributed to the interpretation of the results; All authors critically reviewed and edited the manuscript for important intellectual content, and all authors approved the final manuscript.

Compliance with Ethical Standards
Conflict of interest Each author declares that he/she has no conflict of interest.
Ethical Approval Medical research involving human subjects includes research on identifiable human material or identifiable data. All procedures performed in studies involving human subjects were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.
Informed Consent This retrospective study of medical records data did not involve any participant contact.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.