Video Abstract

Video Abstract

Close modal
BACKGROUND:

Children born preterm are at high risk for autism spectrum disorder (ASD). However, there is still a lack of appropriate developmental markers. In this study, we aim to examine whether early mental performance trajectory is related to ASD outcome in the preterm population.

METHODS:

The population-based cohort included 414 very preterm survivors born between 2008 and 2014. After excluding children with severe neurosensory impairment, 319 children with available records of developmental quotients before age 2 years were enrolled. The trajectory of mental performance evaluated by using the Bayley Scales of Infant Development across 6, 12, and 24 months of age was analyzed with group-based trajectory modeling. At 5 years of age, the ASD diagnosis was established by using the Autism Diagnostic Observation Schedule and the Autism Diagnostic Interview–Revised.

RESULTS:

There were 29 children with ASD and 290 children without ASD. The mental performances from age 6 to 24 months could be classified into 3 trajectory patterns: low declining, high declining, and high stable, which corresponded to ASD prevalence at age 5 years of 35%, 9%, and 3%, respectively. ASD odds was 15 times higher in the low-declining group than in the high-stable group (odds ratio 15; 95% confidence interval 3.8–59; P < .001). Through the analysis of multinomial logistic regression, we found that male infants with longer exposure to oxygen therapy whose mothers had lower maternal education levels tended to follow the low-declining trajectory.

CONCLUSIONS:

The early-life mental trajectory patterns, by using the Bayley Scales of Infant Development, may lead to identification of vulnerable children born preterm for early ASD diagnosis and targeted intervention.

What’s Known on This Subject:

Early developmental trajectory is an indicator for autism spectrum disorder (ASD) in the general population. Preterm infants are at high risk of ASD. However, appropriate developmental markers at toddler age are still lacking.

What This Study Adds:

In this population-based cohort, using group-based trajectory modeling, we found there are 3 patterns of mental performance trajectory for children born preterm from age 6 to 24 months, which is related to different susceptibility to ASD at age 5 years.

Autism spectrum disorder (ASD) represents a spectrum of social communication deficits characterized by occurrence of restricted, repetitive behavioral patterns starting in early childhood.1  Early diagnosis for early intervention is considered a crucial factor for good outcomes for ASD.2  As the prevalence of ASD increases, it has raised concerns related to early detection through developmental surveillance in early life.3,4  Children born preterm have recently been recognized to be at high risk for ASD.5,6  However, several ASD screening tools, which are commonly administered in the general population, have failed to reveal accuracy in the preterm population at 18 to 36 months of age.710  Hence, there is still a lack of appropriate developmental markers for early detection of ASD risk in toddlers who are born preterm.

In the term population, the developmental signature before 36 months of age is considered an early ASD marker.1114  An association between the neurodevelopmental trajectory and ASD risk has been validated in a general population cohort; however, this study did not include children born preterm.15  Whether the early-life developmental trajectory contributes to the ASD risk in preterm populations remains unclear.

The Bayley Scales of Infant Development (BSID), which measures mental, language, and motor performance, is a widely used tool in the preterm population because of their high risk of neurodevelopmental impairment.16,17  However, the relationship between the early neurodevelopmental trajectory patterns evaluated by using the BSID and the ASD diagnosis in the preterm population is unknown. Group-based trajectory modeling (GBTM) has been applied to capture the heterogeneity of trajectory within study populations and also to facilitate a causal inference.18  Here, we used GBTM to examine the relationship between the longitudinal neurodevelopmental trajectory from age 6 to 24 months (assessed with the BSID) and the universal ASD outcome (measured with the Autism Diagnostic Interview–Revised [ADI-R] and the Autism Diagnostic Observation Schedule [ADOS]) at 5 years of age in the preterm population in south Taiwan.1925  We hypothesized that the mental performance trajectory patterns and the associated neonatal risk factors are useful to characterize the early developmental markers of ASD in the preterm population.

In total, 414 very preterm infants (birth weight <1500 g; gestational age <32 weeks) who survived to discharge from 2008 to 2014 from the 4 NICUs in medical centers of Tainan City in southern Taiwan were enrolled. The demographics, perinatal and neonatal risk factors, and morbidities were collected during the infants’ hospitalization in the NICUs (Demographic Data and Risk Factors section of the Supplemental Information). After discharge, the very preterm survivors were longitudinally followed-up for neurodevelopmental status in a single outpatient clinic at the university hospital.

Among the 414 preterm survivors, 360 (87%) could be followed-up to age 5 years (Supplemental Fig 3). The demographics of the children who missed the follow-up are shown in Supplemental Table 4. After excluding the preterm children who had cerebral palsy (n = 24), hearing impairment (n = 5), or congenital malformation (n = 7) for the confounding effect of concomitant neurologic disorders, 324 children were used for analysis (Criteria of Exclusion section of the Supplemental Information).26,27  This study was approved by the university hospital’s institutional review board. For each subject, informed consent was obtained from the parents during the hospitalization and at the follow-up visit.

At the follow-ups at 6, 12, and 24 months’ corrected age, child psychologists evaluated the children’s neurodevelopment using the Bayley Scales of Infant Development, Second Edition (BSID-II) for children born before 2011 and the Bayley Scales of Infant and Toddler Development, Third Edition (BSID-III) for those born after 2011. The BSID-II produces 2 composite scores: a mental developmental index (MDI) to assess cognition and language and a psychomotor developmental index for motor skills. The BSID-III has 3 composite scores: cognition, language, and motor function (BSID Scores section of the Supplemental Information). The cognition and language composite scores of the BSID-III were converted into the predicted MDI according to the algorithm proposed by Moore et al28  for analysis. The predicted MDI converted from the BSID-III cognition and language scores, along with the MDI in the BSID-II, was used for a trajectory analysis to represent the mental performance of these preterm children.

At age 5 years, cognitive function was evaluated by using the Wechsler Preschool and Primary Scale of Intelligence, Revised for children born before 2012 and the Wechsler Preschool and Primary Scale of Intelligence, Fourth Edition for those born after 2012. ASD diagnosis was made according to the Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition, ADI-R, and ADOS, First Edition.1,21,22  The team for the ASD evaluation was supervised by C.-H.C., an ASD specialist who received research training and research certification in the United States.24,25  The psychologists performing the ADOS and the ADI-R were blinded to their previous neurodevelopmental scores. The diagnosis of ASD was established after a thorough discussion of the results of the diagnostic assessments by the team.

Module 3 of the ADOS was performed.22  The cutoffs for ASD in the domains of communication and social interaction and the communication and social interaction total score were 2, 4, and 7, respectively. However, when interactions could not be elicited by tasks in module 3, examiners would use module 2 or 1 for appropriate evaluations. The ASD symptom severity was represented by the ADOS severity scores (including social-affect deficits, restricted, repetitive behaviors, and a total score), and the ADI-R scores included qualitative abnormalities in reciprocal social interaction (domain A); qualitative abnormalities in communication (domain B); restricted, repetitive, and stereotyped behavioral patterns (domain C); and abnormality in development evident at or before age 36 months (domain D) (The ADOS Severity Scores and the ADI-R Score section of the Supplemental Information).21,22,2931 

Mental performance represented by the MDI from the BSID-II and the cognition and language scores from the BSID-III at age 6, 12, and 24 months were analyzed by using GBTM to classify homogeneous trajectory clusters without a priori assignment of ASD.18,32  Fit criterion assessment plots were employed to select the best model fit for further analysis.33  The Bayesian information criterion and the Akaike information criterion were used for model selection. Association of trajectory groups and ASD was analyzed by using a logistic regression with a Firth approximation.34  A stepwise multinomial logistic regression was applied to select the important risk factors for the mental performance trajectories at an entry P value of .1 and a stay P value of .15.35  Bootstrapping with 10 repeated samples was then used to compare the subsets of the selected variables by using the averaged Akaike information criterion and the Bayesian information criterion.36  Risk factors were used to predict the attribution of trajectory groups by the multinomial logistic regression and multiclass receiver operating characteristic curves.37  The association between the mental developmental trajectory groups with ASD was analyzed by using a multinomial logistic regression adjusted for the selected risk factors in the best fit model.38  The 5-year-old cognition and ASD symptoms among the trajectory groups were compared by using the Kruskal-Wallis test and the Dunn’s post hoc comparison. A statistical analysis was performed by using GraphPad Prism 5, SPSS software version 17.0 (SPSS Inc, Chicago, IL), SAS software version 9.2 (SAS Institute, Inc, Cary, NC), and the R project version 3.6.1.

Among the 324 children born preterm followed-up at 5 years of age, 30 were diagnosed with ASD and 294 did not have ASD. Serial BSID data from age 6 to 24 months were available in 319 children (BSID-II: n = 132; BSID-III: n = 187), including 29 children with ASD and 290 without ASD. The ADI-R and ADOS scores of the 29 children with ASD are shown in Supplemental Figs 4 and 5. Scores in domain D of the ADI-R revealed that 86% (25 of 29) of the children with ASD had ASD-related developmental problems at or before 36 months of age.

GBTM was used to examine the mental performance trajectory of these children across the period of 6 to 24 months. According to the Bayesian information criterion and the Akaike information criterion, a quadratic shape revealed a better fit (Supplemental Table 5). However, the model fits for the 3- and 4-group classifications were similar. For the purpose of the appropriateness of the estimated sample size, we selected the 3-group trajectory model for further analysis (Supplemental Table 6).

The 3-trajectory classification indicated the patterns of low-declining (n = 23; 7%), high-declining (n = 198; 62%), and high-stable (n = 98; 31%) mental performance dynamics (Fig 1). The high-stable group maintained normal and stable mental development across time, whereas both the low-declining and high-declining groups showed quadratic changes in mental performance with time. The 3 trajectory groups started at different intercepts at 6 months of age (high stable [93.8] > high declining [89.6] > low declining [72.7]). By 12 months of age, the low-declining group still had the worst mental performance (high stable [94.8] = high declining [90.3] > low declining [79.6]), but the mental scores of the high-declining and high-stable groups became indistinguishable. It was between age 12 and 24 months that the mental abilities of the high-declining and high-stable groups started to diverge. By 24 months of age, the high-stable group performed significantly better than the other 2 groups (high stable [96.9] > high declining [79.7] = low declining [59.1]).

FIGURE 1

The trajectory patterns of mental performance from 6 to 24 months of age. The mental performance trajectory was classified into 3 trajectory groups: low-declining group (n = 23 [7%]), high-declining group (n = 198 [62%]), and high-stable group (n = 98 [31%]). The dotted lines represent 95% CIs.

FIGURE 1

The trajectory patterns of mental performance from 6 to 24 months of age. The mental performance trajectory was classified into 3 trajectory groups: low-declining group (n = 23 [7%]), high-declining group (n = 198 [62%]), and high-stable group (n = 98 [31%]). The dotted lines represent 95% CIs.

Close modal

The ASD prevalence rates revealed a gradient of changes across the 3 trajectory groups: as high as 35% in the low-declining group, 9% in the high-declining group, and 3% in the high-stable group (Fig 2). Compared with the high-stable group, the low-declining group had 15 times higher odds of ASD (odds ratio 15; 95% confidence interval [CI] 3.8–59; P < .001). The high-declining group also had a higher prevalence of ASD, although the odds did not reveal statistical differences (odds ratio 2.9; 95% CI 0.9–9.3; P = .08).

FIGURE 2

The ASD prevalence at age 5 years among the 3 mental performance trajectory groups across 6 to 24 months of age. ASD was most prevalent in the low-declining group (35%), followed by the high-declining group (9%) and the high-stable group (3%). The low-declining group had significantly higher odds of ASD than the high-stable group (odds ratio 15; 95% CI 3.8–59; P < .001).

FIGURE 2

The ASD prevalence at age 5 years among the 3 mental performance trajectory groups across 6 to 24 months of age. ASD was most prevalent in the low-declining group (35%), followed by the high-declining group (9%) and the high-stable group (3%). The low-declining group had significantly higher odds of ASD than the high-stable group (odds ratio 15; 95% CI 3.8–59; P < .001).

Close modal

Among the 70 children whose 6-month mental scores were <85, 18 children exhibited improved mental scores >85, whereas the other 52 children did not by 24 months of age. None of the children with improved mental scores developed ASD at 5 years of age, in contrast to the 23% (12 of 52) children without improvement who did develop ASD.

From the low-declining to high-declining group and up to the high-stable group, the children who showed inferior mental performance tended to comprise more boys and higher rates of lower gestational age, smaller birth weight, small for gestation age, and lower maternal education (Table 1). These children with lower mental performance were significantly more prone to have a longer duration on oxygen therapy, a higher rate of chronic lung disease, and a greater length of hospitalization. A multivariate regression with bootstrapping using stepwise selection revealed that 3 risk factors, sex, maternal educational status, and duration on oxygen therapy, significantly affected the patterns of mental trajectory that led to different rates of ASD at age 5 years (Table 2). As days on oxygen extended, the likelihood for the low-declining trajectory increased, whereas that for the high-stable trajectory decreased (Supplemental Fig 6). The area under the receiver operating characteristic curve was 0.74 by using the 3 risk factors to predict the development of the low-declining trajectory (Supplemental Fig 7).

TABLE 1

Differences in Demographic Data and Neonatal Morbidities Among the 3 Mental Performance Trajectory Groups

Low-Declining Group (n = 23)High-Declining Group (n = 198)High-Stable Group (n = 98)
Demographics    
 Male sex,an (%) 18 (78) 112 (57) 38 (39) 
 Gestational age,a wk, mean (SD) 27 (3) 28 (2) 29 (2) 
 Birth wt,a g, mean (SD) 922 (235) 1076 (254) 1141 (225) 
 Maternal age, y, mean (SD) 31 (5) 32 (5) 32 (4) 
 Paternal age, y, mean (SD) 36 (7) 34 (5) 35 (4) 
 Maternal education, university and graduate school,an (%) 5 (22) 60 (30) 43 (44) 
 Paternal education, university and graduate school, n (%) 8 (35) 79 (40) 43 (44) 
Neonatal morbidities    
 Small for gestational age,an (%) 10 (44) 38 (19) 17 (17) 
 Intraventricular hemorrhage, grades 1 and 2, n (%) 2 (9) 34 (17) 24 (25) 
 Intraventricular hemorrhage, grades 3 and 4, n (%) 1 (4) 8 (4) 6 (6) 
 Cystic periventricular leukomalacia, n (%) 1 (4) 3 (2) 0 (0) 
 Hemodynamically significant patent ductus arteriosus,bn (%) 8 (35) 87 (44) 33 (34) 
 Blood culture–proven sepsis, n (%) 5 (22) 26 (13) 19 (19) 
 Necrotizing enterocolitis,cn (%) 2 (9) 12 (6) 4 (4) 
 Chronic lung disease,an (%) 16 (70) 133 (67) 49 (50) 
 Retinopathy of prematurity,dn (%) 12 (52) 81 (41) 42 (43) 
 Steroids for chronic lung disease, n (%) 2 (9) 15 (8) 5 (5) 
 NSAIDs for patent ductus arteriosus, n (%) 7 (30) 87 (44) 31 (32) 
 Duration on oxygen therapy,a d, mean (SD) 63 (43) 44 (26) 35 (24) 
 Length of hospitalization,a d, mean (SD) 85 (36) 65 (25) 60 (22) 
Low-Declining Group (n = 23)High-Declining Group (n = 198)High-Stable Group (n = 98)
Demographics    
 Male sex,an (%) 18 (78) 112 (57) 38 (39) 
 Gestational age,a wk, mean (SD) 27 (3) 28 (2) 29 (2) 
 Birth wt,a g, mean (SD) 922 (235) 1076 (254) 1141 (225) 
 Maternal age, y, mean (SD) 31 (5) 32 (5) 32 (4) 
 Paternal age, y, mean (SD) 36 (7) 34 (5) 35 (4) 
 Maternal education, university and graduate school,an (%) 5 (22) 60 (30) 43 (44) 
 Paternal education, university and graduate school, n (%) 8 (35) 79 (40) 43 (44) 
Neonatal morbidities    
 Small for gestational age,an (%) 10 (44) 38 (19) 17 (17) 
 Intraventricular hemorrhage, grades 1 and 2, n (%) 2 (9) 34 (17) 24 (25) 
 Intraventricular hemorrhage, grades 3 and 4, n (%) 1 (4) 8 (4) 6 (6) 
 Cystic periventricular leukomalacia, n (%) 1 (4) 3 (2) 0 (0) 
 Hemodynamically significant patent ductus arteriosus,bn (%) 8 (35) 87 (44) 33 (34) 
 Blood culture–proven sepsis, n (%) 5 (22) 26 (13) 19 (19) 
 Necrotizing enterocolitis,cn (%) 2 (9) 12 (6) 4 (4) 
 Chronic lung disease,an (%) 16 (70) 133 (67) 49 (50) 
 Retinopathy of prematurity,dn (%) 12 (52) 81 (41) 42 (43) 
 Steroids for chronic lung disease, n (%) 2 (9) 15 (8) 5 (5) 
 NSAIDs for patent ductus arteriosus, n (%) 7 (30) 87 (44) 31 (32) 
 Duration on oxygen therapy,a d, mean (SD) 63 (43) 44 (26) 35 (24) 
 Length of hospitalization,a d, mean (SD) 85 (36) 65 (25) 60 (22) 

NSAID, nonsteroidal antiinflammatory drug.

a

Differences with statistical significances among the 3 trajectory groups.

b

Patent ductus arteriosus requiring medical or surgical treatments.

c

Including stage IIA and above.

d

Including stage I and above.

TABLE 2

Multivariate Analysis of Demographics and Neonatal Morbidities on Trajectory Membership (the High-Stable Group Served as the Reference)

Risk FactorsLow-Declining GroupHigh-Declining Group
Odds Ratio95% CIPOdds Ratio95% CIP
Male sex 7.3 2.2–24.3 .001 2.1 1.3–3.6 .004 
Maternal education, university and graduate school 0.3 0.1–0.9 .03 0.5 0.3–0.8 .008 
Duration on oxygen therapy, d 1.04 1.02–1.05 <.001 1.02 1.004–1.03 .006 
Risk FactorsLow-Declining GroupHigh-Declining Group
Odds Ratio95% CIPOdds Ratio95% CIP
Male sex 7.3 2.2–24.3 .001 2.1 1.3–3.6 .004 
Maternal education, university and graduate school 0.3 0.1–0.9 .03 0.5 0.3–0.8 .008 
Duration on oxygen therapy, d 1.04 1.02–1.05 <.001 1.02 1.004–1.03 .006 

A multinomial logistic regression model was used to elucidate the roles of the 3 risk factors associated with the mental performance trajectory from age 6 to 24 months and the occurrence of ASD at age 5 years (Table 3). After adjustment, maternal educational status and the duration on oxygen therapy became less significant to the odds of ASD, revealing that their mediations to ASD were dependent on the effects of the mental performance trajectory. In contrast, sex remained significant, suggesting its direct influence on ASD odds.

TABLE 3

Mediation Effects of Risk Factors on the Association Between Trajectory Classification and ASD

EstimateSEWald χ2 TestPOdds Ratio95% CI
Factors predicting ASD       
 Low-declining trajectorya 2.5 0.8 9.8 .002 11.8 2.5–55.2 
 High-declining trajectorya 0.9 0.7 2.0 .2 2.5 0.7–9.0 
 Male sex 1.0 0.5 4.3 .04 2.8 1.1–7.3 
 Maternal education, university and graduate school −0.2 0.4 0.1 .7 0.8 0.4–2.0 
 Duration on oxygen therapy 0.007 0.007 0.9 .3 1.007 0.993–1.021 
Overall model −4.1 0.8 28.4 <.001 — — 
EstimateSEWald χ2 TestPOdds Ratio95% CI
Factors predicting ASD       
 Low-declining trajectorya 2.5 0.8 9.8 .002 11.8 2.5–55.2 
 High-declining trajectorya 0.9 0.7 2.0 .2 2.5 0.7–9.0 
 Male sex 1.0 0.5 4.3 .04 2.8 1.1–7.3 
 Maternal education, university and graduate school −0.2 0.4 0.1 .7 0.8 0.4–2.0 
 Duration on oxygen therapy 0.007 0.007 0.9 .3 1.007 0.993–1.021 
Overall model −4.1 0.8 28.4 <.001 — — 

—, not applicable.

a

Compared with the high-stable group.

At age 5 years, the 3 mental trajectory groups had significant differences in mental capabilities and autistic behavior (Supplemental Table 7). Mental abilities, including general IQ, verbal communication skill, and visual spatial performance, varied among the 3 trajectory groups. In addition, the ADOS severity scores and ADI-R scores for ASD revealed that the low-declining group had significantly higher scores in every domain, indicating more severe autistic behaviors related to social communication and restricted, repetitive behaviors than in the other 2 groups.

We examined the ASD risk at 5 years of age in children born preterm with different mental trajectory patterns from age 6 to 24 months using BSID examinations. We found 3 distinct trajectory patterns, and the low-declining pattern was associated the highest ASD odds. The vulnerable trajectory pattern for ASD tended to be from male infants, with higher rates of lower maternal education levels and prolonged use of oxygen therapy. Our findings suggest that an early mental developmental trajectory is useful to characterize the differential risk for ASD outcome and that demographic and neonatal risk factors contribute to the trajectory, leading to higher ASD vulnerability.

The evolution of mental ability is better defined by longitudinal assessments than a cross-sectional evaluation.39  In our study, we classified the mental performance trajectory of children born preterm from age 6 to 24 months using GBTM and then associated the trajectory pattern to the ASD odds at 5 years of age. Because the heterogeneity of mental development existed in both the children with ASD and the children without ASD, our results provided a different perspective on ASD susceptibility based on the dynamic patterns of neurodevelopment in preterm infants.

The BSID has been applied extensively in the preterm population.16,17  However, corresponding the BSID results to future ASD susceptibility has not been reported. Because the validity of ASD screening tools has not yet been confirmed for the preterm population at toddler age, the mental performance trajectory provided by the BSID may give alerts to ASD concerns. Clinicians and parents may put emphasis on the behavioral aspect in the child born preterm who shows a consistently low or a declining mental trajectory.

The trajectory of early language and nonverbal cognitive development has been recognized as a risk indicator in term children.11  Children with ASD exhibit slower acquisition of both verbal and nonverbal skills, leading to a dramatic decline in the composite IQ score during the second year of life.11,12,14,15,40  However, the hypothesis has not been tested in children born preterm. Our finding of the association between mental performance trajectories and future ASD odds suggests that the early mental developmental trajectory may serve as a marker for ASD in preterm populations.

Early interventions to enhance person-environment fit is crucial to good outcomes in children with ASD.2,41,42  The mental trajectory findings suggested that ASD can be suspected by 24 months of age. In addition, most children with ASD had evident developmental problems at or before 36 months of age according to our ADI-R scores. Taken together, our study indicates the possibility of ASD diagnosis at a younger age in preterm populations.

A study on term siblings with ASD revealed that the developmental trajectory of children with ASD began to deviate from 14 to 24 months of age.12  The cognitive trajectory in a general population also identified a high risk of ASD in children with declining abilities after 12 months of age.15,43  Our study revealed that preterm children with an ASD risk of 9% in the high-declining group were also characterized by a decline in mental performance between 12 and 24 months. The mental ability of the low-declining group, who showed the highest ASD prevalence of 35%, initiated from a significantly lower intercept at 6 months of age and continued to exhibit the worst scores at age 12 and 24 months. Hence, our study revealed that the differential time window for deviations in mental developmental carried different ASD risks in the preterm population: the children with atypical mental development in the first year of life had the highest ASD risk, whereas the children with a delay in mental development in the second year had an intermediate ASD risk.

Demographic and neonatal risk factors also interacted to affect the development of the trajectory pattern. We found that intrinsic factors and neonatal morbidities affected preterm infants’ mental trajectory development. Using multinomial selection, we identified that male sex, prolonged oxygen exposure, and lower maternal education level had influential effects on the preterm infants, causing them to follow the low-declining mental trajectory.

Sex is known to play a significant role in outcomes in preterm populations.4449  We recognized that male sex was predictive of the trajectory of poorer mental performance and higher ASD odds. In addition, sex not only influenced the attribution of the mental developmental trajectory, which in turn brought about a distinct ASD risk, but also imposed a direct effect on ASD irrespective of cognitive function. These findings suggest the dual role of sex on ASD risk in children born preterm.

Ventilator support and chronic lung disease tend to increase cognitive impairment and ASD risk in preterm infants.27,5053  Preterm infants frequently require oxygen supplementation for respiratory insufficiency, but they are also vulnerable to oxidative stress because of immature antioxidant defenses.54,55  In our study, we identified the effects of chronic lung disease and longer duration of oxygen therapy on mental performance trajectories and ASD risk, but the underlying mechanism affecting immature brain development remains unclear.55 

The significance of maternal education is well-established in developmental psychopathology.5659  Our study revealed that higher maternal education exerted a protective effect on children’s mental trajectory development against ASD risk. Maternal educational status may represent socioeconomical class, the quality of parenting skills, and also opportunities for cognitive stimulation at home.60  The influence of maternal education highlights the importance of nurturing environments, apart from the physical issues, such as morbidities, associated with preterm birth.

In our cohort, 87% of children were managed for neurodevelopment from age 6 months to age 5 years. However, it is possible that the missing data in the remaining 13% of children could have affected the results.18,32  Detailed family history and parent-child interaction patterns are needed for analysis of the genetic and environmental influences on ASD risk in children born preterm. Because extrauterine immature brain development is vulnerable to adverse exposures during critical care, further studies incorporating the clinical big data, as well as sophisticated neuroimaging analyses, may assist in predicting which trajectory an infant might follow for ASD outcome. Further multicenter studies recruiting more preterm infants for longitudinal neurodevelopmental follow-up and ASD diagnosis are necessary to validate our findings.

The current study revealed that early-life mental performance trajectory patterns by using the BSID for children born preterm were related to different susceptibility to ASD at age 5 years. The mental trajectory patterns may lead to early identification of vulnerable children for early ASD diagnosis and targeted intervention.

The corresponding author, Dr Chao-Ching Huang, had full access to the data set used and analyzed during the current study. This article was edited by the Foreign Language Center, National Cheng Kung University (Tainan, Taiwan). We thank the Premature Baby Foundation of Taiwan and all team members in charge of data collection and assessment of the children. None of these individuals were compensated for their contributions.

Dr Chen designed and conducted the study, assisted in the statistical analysis, and drafted the manuscript; Dr S.-T. Wang performed the formal statistical analysis and drafted the manuscript; Drs L.-W. Wang and Kao conducted the study, assisted in the interpretation of data, and revised the manuscript for important intellectual content; Dr Chu conducted the study, interpreted the data, provided input on the data analysis, and revised the manuscript for important intellectual content; Dr Wu conducted the study, interpreted the data, and revised the manuscript for important intellectual content; Dr Chiang coordinated and supervised the conduction of the study, curated the data, and made critical manuscript revisions; Dr Huang conceptualized and designed the study, curated the data, and revised the manuscript; and all authors approved the final manuscript as submitted and agree to be accountable for all aspects of the work.

FUNDING: Supported by grants from the Taiwan Ministry of Science and Technology (MOST 108-2321-B006-023-MY2 and MOST-105-2410-H-004-071-MY3), the joint grant of National Cheng Kung University and E-Da Hospital (NCKUEDA10306), and the grant of National Cheng Kung University Hospital (NCKUH-10809002). The funding sources had no roles in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication. The content is solely the responsibility of the authors and does not necessarily represent the official views of the funding agencies.

ADI-R

Autism Diagnostic Interview–Revised

ADOS

Autism Diagnostic Observation Schedule

ASD

autism spectrum disorder

BSID

Bayley Scales of Infant Development

BSID-II

Bayley Scales of Infant Development, Second Edition

BSID-III

Bayley Scales of Infant and Toddler Development, Third Edition

CI

confidence interval

GBTM

group-based trajectory modeling

MDI

mental developmental index

1
American Psychiatric Association
.
Diagnostic and Statistical Manual of Mental Disorders
, 5th ed.
Washington, DC
:
American Psychiatric Association
;
2013
2
Lai
MC
,
Anagnostou
E
,
Wiznitzer
M
,
Allison
C
,
Baron-Cohen
S
.
Evidence-based support for autistic people across the lifespan: maximising potential, minimising barriers, and optimising the person-environment fit
.
Lancet Neurol
.
2020
;
19
(
5
):
434
451
3
Baio
J
,
Wiggins
L
,
Christensen
DL
, et al
.
Prevalence of autism spectrum disorder among children aged 8 years - autism and developmental disabilities monitoring network, 11 sites, United States, 2014
.
MMWR Surveill Summ
.
2018
;
67
(
6
):
1
23
4
Zwaigenbaum
L
,
Bauman
ML
,
Fein
D
, et al
.
Early screening of autism spectrum disorder: recommendations for practice and research
.
Pediatrics
.
2015
;
136
(
suppl 1
):
S41
S59
5
Agrawal
S
,
Rao
SC
,
Bulsara
MK
,
Patole
SK
.
Prevalence of autism spectrum disorder in preterm infants: a meta-analysis
.
Pediatrics
.
2018
;
142
(
3
):
e20180134
6
Xie
S
,
Heuvelman
H
,
Magnusson
C
, et al
.
Prevalence of autism spectrum disorders with and without intellectual disability by gestational age at birth in the Stockholm youth cohort: a register linkage study
.
Paediatr Perinat Epidemiol
.
2017
;
31
(
6
):
586
594
7
Kim
SH
,
Joseph
RM
,
Frazier
JA
, et al;
Extremely Low Gestational Age Newborn (ELGAN) Study Investigators
.
Predictive validity of the modified checklist for autism in toddlers (M-CHAT) born very preterm
.
J Pediatr
.
2016
;
178
:
101
107.e2
8
Hrdlicka
M
,
Dudova
I
.
Screening preterm children for autism at 2 years of age
.
J Pediatr
.
2015
;
167
(
1
):
212
9
Gray
PH
.
M-CHAT autism screening may be inaccurate among toddlers born very preterm
.
J Pediatr
.
2017
;
182
:
401
404
10
Boone
KM
,
Brown
AK
,
Keim
SA
.
Screening accuracy of the brief infant toddler social-emotional assessment to identify autism spectrum disorder in toddlers born at less than 30 weeks’ gestation
.
Child Psychiatry Hum Dev
.
2018
;
49
(
4
):
493
504
11
Zwaigenbaum
L
,
Bauman
ML
,
Stone
WL
, et al
.
Early identification of autism spectrum disorder: recommendations for practice and research
.
Pediatrics
.
2015
;
136
(
suppl 1
):
S10
S40
12
Landa
R
,
Garrett-Mayer
E
.
Development in infants with autism spectrum disorders: a prospective study
.
J Child Psychol Psychiatry
.
2006
;
47
(
6
):
629
638
13
Landa
RJ
,
Gross
AL
,
Stuart
EA
,
Faherty
A
.
Developmental trajectories in children with and without autism spectrum disorders: the first 3 years
.
Child Dev
.
2013
;
84
(
2
):
429
442
14
Landa
RJ
,
Gross
AL
,
Stuart
EA
,
Bauman
M
.
Latent class analysis of early developmental trajectory in baby siblings of children with autism
.
J Child Psychol Psychiatry
.
2012
;
53
(
9
):
986
996
15
Nishimura
T
,
Takei
N
,
Tsuchiya
KJ
.
Neurodevelopmental trajectory during infancy and diagnosis of autism spectrum disorder as an outcome at 32 months of age
.
Epidemiology
.
2019
;
30
(
suppl 1
):
S9
S14
16
Bode
MM
,
DʼEugenio
DB
,
Mettelman
BB
,
Gross
SJ
.
Predictive validity of the Bayley, Third Edition at 2 years for intelligence quotient at 4 years in preterm infants
.
J Dev Behav Pediatr
.
2014
;
35
(
9
):
570
575
17
Yu
YT
,
Hsieh
WS
,
Hsu
CH
, et al
.
A psychometric study of the Bayley Scales of Infant and Toddler Development - 3rd Edition for term and preterm Taiwanese infants
.
Res Dev Disabil
.
2013
;
34
(
11
):
3875
3883
18
Nagin
DS
,
Odgers
CL
.
Group-based trajectory modeling in clinical research
.
Annu Rev Clin Psychol
.
2010
;
6
:
109
138
19
Tu
YF
,
Wang
LW
,
Wang
ST
,
Yeh
TF
,
Huang
CC
.
Postnatal steroids and febrile seizure susceptibility in preterm children
.
Pediatrics
.
2016
;
137
(
4
):
e20153404
20
Wang
LW
,
Lin
YC
,
Wang
ST
,
Huang
CC
;
on behalf of the Taiwan Premature Infant Developmental Collaborative Study Group
.
Identifying risk factors shared by bronchopulmonary dysplasia, severe retinopathy, and cystic periventricular leukomalacia in very preterm infants for targeted intervention
.
Neonatology
.
2018
;
114
(
1
):
17
24
21
Lord
C
,
Rutter
M
,
Le Couteur
A
.
Autism Diagnostic Interview-Revised: a revised version of a diagnostic interview for caregivers of individuals with possible pervasive developmental disorders
.
J Autism Dev Disord
.
1994
;
24
(
5
):
659
685
22
Lord
C
,
Rutter
M
,
Goode
S
, et al
.
Autism diagnostic observation schedule: a standardized observation of communicative and social behavior
.
J Autism Dev Disord
.
1989
;
19
(
2
):
185
212
23
Chen
LW
,
Wang
ST
,
Wang
LW
, et al
.
Behavioral characteristics of autism spectrum disorder in very preterm birth children
.
Mol Autism
.
2019
;
10
:
32
24
Chu
CL
,
Chiang
CH
,
Wu
CC
,
Hou
YM
,
Liu
JH
.
Service system and cognitive outcomes for young children with autism spectrum disorders in a rural area of Taiwan
.
Autism
.
2017
;
21
(
5
):
581
591
25
Chiang
CH
,
Wu
CC
,
Hou
YM
,
Chu
CL
,
Liu
JH
,
Soong
WT
.
Development of T-STAT for early autism screening
.
J Autism Dev Disord
.
2013
;
43
(
5
):
1028
1037
26
Palisano
R
,
Rosenbaum
P
,
Walter
S
,
Russell
D
,
Wood
E
,
Galuppi
B
.
Development and reliability of a system to classify gross motor function in children with cerebral palsy
.
Dev Med Child Neurol
.
1997
;
39
(
4
):
214
223
27
Adams-Chapman
I
,
Heyne
RJ
,
DeMauro
SB
, et al;
Follow-Up Study of the Eunice Kennedy Shriver National Institute of Child Health and Human Development Neonatal Research Network
.
Neurodevelopmental impairment among extremely preterm infants in the Neonatal Research Network
.
Pediatrics
.
2018
;
141
(
5
):
e20173091
28
Moore
T
,
Johnson
S
,
Haider
S
,
Hennessy
E
,
Marlow
N
.
Relationship between test scores using the second and third editions of the Bayley Scales in extremely preterm children
.
J Pediatr
.
2012
;
160
(
4
):
553
558
29
Gotham
K
,
Risi
S
,
Pickles
A
,
Lord
C
.
The Autism Diagnostic Observation Schedule: revised algorithms for improved diagnostic validity
.
J Autism Dev Disord
.
2007
;
37
(
4
):
613
627
30
Gotham
K
,
Pickles
A
,
Lord
C
.
Standardizing ADOS scores for a measure of severity in autism spectrum disorders
.
J Autism Dev Disord
.
2009
;
39
(
5
):
693
705
31
Hus
V
,
Gotham
K
,
Lord
C
.
Standardizing ADOS domain scores: separating severity of social affect and restricted and repetitive behaviors
.
J Autism Dev Disord
.
2014
;
44
(
10
):
2400
2412
32
Nagin
DS
.
Group-based trajectory modeling: an overview
.
Ann Nutr Metab
.
2014
;
65
(
2–3
):
205
210
33
Klijn
SL
,
Weijenberg
MP
,
Lemmens
P
,
van den Brandt
PA
,
Lima Passos
V
.
Introducing the fit-criteria assessment plot - a visualisation tool to assist class enumeration in group-based trajectory modelling
.
Stat Methods Med Res
.
2017
;
26
(
5
):
2424
2436
34
Puhr
R
,
Heinze
G
,
Nold
M
,
Lusa
L
,
Geroldinger
A
.
Firth’s logistic regression with rare events: accurate effect estimates and predictions?
Stat Med
.
2017
;
36
(
14
):
2302
2317
35
Tutz
G
,
Pößnecker
W
,
Uhlmann
L
.
Variable selection in general multinomial logit models
.
Comput Stat Data Anal
.
2015
;
82
:
207
222
36
Cherrie
J.
Variable screening for multinomial logistic regression on very large data sets as applied to direct response modeling.
SAS Global Forum;
2007
37
Wang
J
,
Wei
R
,
Jia
W
.
A quick tour of multiROC. Available at: https://mran.microsoft.com/snapshot/2018-02-12/web/packages/multiROC/vignettes/my-vignette.html. Accessed April 1, 2020
38
Hayes
AF
,
Rockwood
NJ
.
Regression-based statistical mediation and moderation analysis in clinical research: observations, recommendations, and implementation
.
Behav Res Ther
.
2017
;
98
:
39
57
39
Ployhart
RE
,
Vandenberg
RJ
.
Longitudinal research: the theory, design, and analysis of change
.
J Manage
.
2010
;
36
(
1
):
94
120
40
Bryson
SE
,
Zwaigenbaum
L
,
Brian
J
, et al
.
A prospective case series of high-risk infants who developed autism
.
J Autism Dev Disord
.
2007
;
37
(
1
):
12
24
41
Nevill
RE
,
Lecavalier
L
,
Stratis
EA
.
Meta-analysis of parent-mediated interventions for young children with autism spectrum disorder
.
Autism
.
2018
;
22
(
2
):
84
98
42
Gengoux
GW
,
Abrams
DA
,
Schuck
R
, et al
.
A pivotal response treatment package for children with autism spectrum disorder: an RCT
.
Pediatrics
.
2019
;
144
(
3
):
e20190178
43
Nishimura
T
,
Takei
N
,
Tsuchiya
KJ
,
Asano
R
,
Mori
N
.
Identification of neurodevelopmental trajectories in infancy and of risk factors affecting deviant development: a longitudinal birth cohort study
.
Int J Epidemiol
.
2016
;
45
(
2
):
543
553
44
Garfinkle
J
,
Yoon
EW
,
Alvaro
R
, et al;
Canadian Neonatal Network Investigators
.
Trends in sex-specific differences in outcomes in extreme preterms: progress or natural barriers?
Arch Dis Child Fetal Neonatal Ed
.
2020
;
105
(
2
):
158
163
45
Joseph
RM
,
O’Shea
TM
,
Allred
EN
, et al
.
Prevalence and associated features of autism spectrum disorder in extremely low gestational age newborns at age 10 years
.
Autism Res
.
2017
;
10
(
2
):
224
232
46
Pinto-Martin
JA
,
Levy
SE
,
Feldman
JF
,
Lorenz
JM
,
Paneth
N
,
Whitaker
AH
.
Prevalence of autism spectrum disorder in adolescents born weighing <2000 grams
.
Pediatrics
.
2011
;
128
(
5
):
883
891
47
Johnson
S
,
Hollis
C
,
Kochhar
P
,
Hennessy
E
,
Wolke
D
,
Marlow
N
.
Autism spectrum disorders in extremely preterm children
.
J Pediatr
.
2010
;
156
(
4
):
525
531.e2
48
Kuban
KC
,
Joseph
RM
,
O’Shea
TM
, et al;
Extremely Low Gestational Age Newborn (ELGAN) Study Investigators
.
Girls and boys born before 28 weeks gestation: risks of cognitive, behavioral, and neurologic outcomes at age 10 years
.
J Pediatr
.
2016
;
173
:
69
75.e1
49
Ferri
SL
,
Abel
T
,
Brodkin
ES
.
Sex differences in autism spectrum disorder: a review
.
Curr Psychiatry Rep
.
2018
;
20
(
2
):
9
50
Twilhaar
ES
,
Wade
RM
,
de Kieviet
JF
,
van Goudoever
JB
,
van Elburg
RM
,
Oosterlaan
J
.
Cognitive outcomes of children born extremely or very preterm since the 1990s and associated risk factors: a meta-analysis and meta-regression
.
JAMA Pediatr
.
2018
;
172
(
4
):
361
367
51
Hack
M
,
Taylor
HG
,
Schluchter
M
,
Andreias
L
,
Drotar
D
,
Klein
N
.
Behavioral outcomes of extremely low birth weight children at age 8 years
.
J Dev Behav Pediatr
.
2009
;
30
(
2
):
122
130
52
Stålnacke
SR
,
Tessma
M
,
Böhm
B
,
Herlenius
E
.
Cognitive development trajectories in preterm children with very low birth weight longitudinally followed until 11 years of age
.
Front Physiol
.
2019
;
10
:
307
53
Kuzniewicz
MW
,
Wi
S
,
Qian
Y
,
Walsh
EM
,
Armstrong
MA
,
Croen
LA
.
Prevalence and neonatal factors associated with autism spectrum disorders in preterm infants
.
J Pediatr
.
2014
;
164
(
1
):
20
25
54
Tarnow-Mordi
W
,
Stenson
B
,
Kirby
A
, et al;
BOOST-II Australia and United Kingdom Collaborative Groups
.
Outcomes of two trials of oxygen-saturation targets in preterm infants
.
N Engl J Med
.
2016
;
374
(
8
):
749
760
55
Saugstad
OD
,
Oei
JL
,
Lakshminrusimha
S
,
Vento
M
.
Oxygen therapy of the newborn from molecular understanding to clinical practice
.
Pediatr Res
.
2019
;
85
(
1
):
20
29
56
Linsell
L
,
Johnson
S
,
Wolke
D
, et al
.
Cognitive trajectories from infancy to early adulthood following birth before 26 weeks of gestation: a prospective, population-based cohort study
.
Arch Dis Child
.
2018
;
103
(
4
):
363
370
57
McManus
BM
,
Poehlmann
J
.
Maternal depression and perceived social support as predictors of cognitive function trajectories during the first 3 years of life for preterm infants in Wisconsin
.
Child Care Health Dev
.
2012
;
38
(
3
):
425
434
58
Yaari
M
,
Mankuta
D
,
Harel-Gadassi
A
, et al
.
Early developmental trajectories of preterm infants
.
Res Dev Disabil
.
2018
;
81
:
12
23
59
Luu
TM
,
Vohr
BR
,
Allan
W
,
Schneider
KC
,
Ment
LR
.
Evidence for catch-up in cognition and receptive vocabulary among adolescents born very preterm
.
Pediatrics
.
2011
;
128
(
2
):
313
322
60
Harding
JF
.
Increases in maternal education and low-income children’s cognitive and behavioral outcomes
.
Dev Psychol
.
2015
;
51
(
5
):
583
599

Competing Interests

POTENTIAL CONFLICT OF INTEREST: The authors have indicated they have no potential conflicts of interest to disclose.

FINANCIAL DISCLOSURE: The authors have indicated they have no financial relationships relevant to this article to disclose.

Supplementary data