{'fraction_non_alphanumeric': 0.03379276303299541, 'fraction_numerical': 0.020617382389793826, 'mean_word_length': 5.006385404789054, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 1, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Developmental Trajectories in Very Preterm Born Children Up to 8 Years: A Longitudinal Cohort Study

Aim: Long-term outcome data in preterm children is often limited to cross-sectional measurement of neurodevelopmental impairment (NDI) at the corrected age of 24-36 months. However, impairments may only become overt during childhood or resolve with time, and individual trajectories in outcome over time may vary. The primary aim of this study was to describe NDI in very preterm born children at three subsequent ages of 2, 5, and 8 years of age. As a secondary aim, a longitudinal analysis was performed on the individual longitudinal trajectories in NDI from 2 to 8 years of age. Methods: Single-center prospective cohort study including children born between 1990 and 2011 below 30 weeks' gestation and followed into 2019. The outcome measurement was NDI assessed at 2, 5, and 8 years of age. NDI is a composite score that includes cognitive, neurological, visual, and auditory functions, in which problems were categorized as none, mild, moderate, or severe. Cognitive function measured as total DQ/IQ score was assessed by standardized psychometric tests. Neurological, visual, and auditory functions were assessed by the neonatologist. Results: In total, 921 children were eligible for follow-up, of whom 726 (79%) children were assessed. No NDI was seen in 54, 54, and 62%, mild NDI was seen in 31, 36, and 30%, and moderate-to-severe NDI was seen in 15, 9.2, and 8.6% of the children at 2, 5, and 8 years, respectively. From 2 to 8 years, 63% of the children remained in the same NDI category, 20% of the children improved to a better NDI category, and 17% deteriorated toward a worse NDI category. No differences were found in baseline characteristics of infants that improved or deteriorated. Extreme prematurity, male gender and low parental education were associated with worse NDI status at all time points. Although we observed considerable individual variation over time in NDI status, the course of the trajectories in NDI were not associated with gestation, gender, and parental education. Conclusions: Continued follow-up until school life is essential in order to provide optimal and individually focused referrals and care when needed.


INTRODUCTION
The number of preterm deliveries below 30 weeks' gestation has increased over the last decades, with increasing survival rates of preterm children (1,2). However, improved survival rates still raise the concern of adverse long-term outcome in the increasing number of surviving preterm children. Preterm born children are known to have a higher risk of physical disabilities as well as cognitive problems later in life (3)(4)(5). Knowledge on neurodevelopmental outcomes of children born at these early gestational ages (GA) is crucial for clinicians and families as this may influence antenatal counseling, resuscitation polices, and NICU guidelines (6)(7)(8).
As the impact of developmental impairment may be different at different stages of development, there is increasing interest in studying development as a dynamic process (9,10). Currently available outcome data are often limited to cross-sectional measurements in toddlerhood, but longitudinal follow-up of children is important. Early suboptimal functioning may form an important signal for later problems or an indication for early intervention, and impairments may persist over childhood into adolescence and adulthood (11)(12)(13)(14)(15)(16). Moreover, there may be considerable variation in individual trajectories that is not detectable in cross-sectional studies (17).
Studies evaluating developmental trajectories in preterm children often have focused on specific components of development, like cognitive, behavioral, or social problems (5,10,13,(18)(19)(20). However, a composite outcome score combining different domains provides a general insight in the amount and kind of disabilities of preterm born children. As developmental problems can arise over a broad spectrum of outcome measures, evaluation of developmental trajectories using a composite outcome might provide additional information. A frequently used indication of adverse long-term outcome is the composite measure of neurodevelopmental impairment (NDI), a score that takes cognitive, neurological, visual, and auditory function into account (8,(21)(22)(23)(24). This outcome measure focuses on severe impairments and provides important prognostic information for clinicians and parents (21).
Since three decades, preterm children born below 30 weeks' gestation are eligible for an extensive follow-up program in our perinatal center. This includes outpatient clinic visits to the neonatologist and psychologist at the corrected age of 2 years and the uncorrected age of 5 and 8 years, making an NDI assessment possible at three subsequent ages. The data collected over a period of more than 20 years provides unique information on the development of very preterm children. Therefore, the primary aim of this study was to describe NDI in very preterm born children who were evaluated at three subsequent ages of 2, 5, and 8 years of age. As a secondary aim, a longitudinal analysis was performed on the individual longitudinal trajectories in NDI from 2 to 8 years of age.

Patient Population
This cohort study included all children born between 1990 and 2011 and followed into 2019, with a gestational age below 30 weeks, who were admitted within 24 h after birth to the neonatal intensive care unit (NICU) of Máxima Medical Centre (MMC). The NICU of MMC serves a 1.6 million population including antenatal and postnatal transfer from six other hospitals in the region. Children from parents living outside the adherence area of MMC and referrals from other NICUs were excluded. The ethical review board of MMC approved the study in accordance with the Dutch law on medical research with humans (WMO).

Data Collection
Data from the outpatient clinic visits were collected prospectively. Neonatal data were retrieved from the individual medical records. Individual characteristics and medical data included gender (male or female); birth weight in grams; gestational age in days (based on ultrasound findings or on the first day of last menstrual period if ultrasound data was not available); small for gestational age [defined as birth weight below the 10th percentile (25)]; multiplicity (dichotomized as single or multiple birth); mode of delivery (dichotomized as vaginal or by caesarean section); complete course of antenatal corticosteroids (defined as two doses of betamethasone given 24 h apart before the start of labor); Apgar score at 5 min postpartum; inborn or outborn NICU; rate of artificial ventilation > 12 h; days of endotracheal intubation on any mode of ventilation; surgical treatment of a persistent ductus arteriosus; intraventricular hemorrhage grade 3 or 4 based on ultrasound (26); cystic periventricular leukomalacia grade 3 (27); severe brain injury (defined as intraventricular hemorrhage grade 3 or 4 or cystic periventricular leukomalacia grade 3); laparotomy for necrotizing enterocolitis or single intestinal perforation; surgical treatment or laser therapy for retinopathy of prematurity; and total days of NICU admission. Socio-economic status was assessed using scores defined by the Netherlands Institute for Social and Cultural Research (The Hague, Netherlands) based on postal code at birth, with an average score of 0 and a positive score reflecting a higher than average status and a negative score reflecting a lower than average status (28). For children not seen for follow-up, reasons for no show were identified.

Follow-Up
All preterm children below 30 weeks' gestation were eligible for our follow-up program. This consisted of outpatient clinic visits to the neonatologist and psychologist at the corrected age of 2 years and the uncorrected age of 5 and 8 years. The neonatologist assessed the child's health and evaluated the neurological, visual, and auditory functions. Neurological outcome was scored as normal, mildly abnormal, or unilateral/bilateral CP, according to the GMFCS classification (29). The psychologist evaluated the child's cognitive function, emotional, and behavioral development.  (30). In addition, the psychologist collected information on educational status of the parents, which was classified as low, middle, or high according to the CBS classification (31). This variable was dichotomized describing whether there was a low education or middle-to-high education. If one of the parents was classified as middle-to-high educated, parental education was classified as middle-to-high.

Neurodevelopmental Outcome
The outcome measure was neurodevelopmental impairment (NDI), a composite score based on cognition, neurological assessment, and presence of visual and/or hearing impairment ( Table 1). NDI was categorized as none, mild, moderate, or severe. NDI was classified as mild if cognitive scores showed a developmental quotient (DQ) or intelligence quotient (IQ) between 70 and 84 (−2 to −1 SD); vision or hearing loss without an aid or with good correction, or abnormal neurological tests in the absence of a neurological syndrome (e.g., posture, coordination, and tone dysregulation disorders). NDI was scored as moderate if cognitive DQ/IQ scores were between 55 and 69 (−3 to −2 SD); limited vision or hearing and the use of aids or the presence of a unilateral cerebral palsy. NDI was scored as severe if cognitive DQ/IQ scores were below 55 (>-3 SD), or blindness, deafness, or bilateral cerebral palsy were present. NDI score was based on the worst determinant in either one of the four categories. If one category was missing, NDI was classified as missing. NDI was determined for examinations at 2, 5, and 8 years of age.

Statistical Analyses
Children with and without follow-up were compared using the Student's T-test or Mann-Whitney U-test for continuous variables, depending on distribution of the data, and using

Study Population and Loss-to-follow-Up
Within the study period (1990-2011), 1,107 children born below < 30 weeks' GA were admitted to the NICU. Of these children 186 (17%) died, leaving 921 children eligible for follow-up at the outpatient clinic (Figure 1). Of these children, 726 (79%) were seen for follow-up. In total, 693 (75%), 658 (71%), and 579 (63%) children had follow-up at 2, 5, and 8 years, respectively. Reasons for the total group of 195 loss-to-follow-up children are shown in Table 2. During a limited period of time, difficulties in availability of staff resulted in a group of 102 relatively low risk children who were not invited for follow-up. Table 3 shows the baseline characteristics, separately for children with and without follow-up. Children with follow-up were more immature at birth, compared to children without follow-up. Socio-economic status was higher in children seen for follow-up. Their NICU admission was significantly more often complicated by PDA and ROP, but less often complicated by a laparotomy. The length of NICU stay in children with follow-up was longer, compared to children without follow-up.  When they got older, more children were seen in a clinic for rehabilitation medicine and dropped-out on follow-up. Including these children in the category moderate-to-severe NDI, the percentage at 8 years increased up to 16%. In further analysis, the original data was used categorizing this subgroup as missing. Separate presentation of NDI rates for EP and VP infants showed decreased "no NDI" and increased "mild NDI" rates in EP infants compared to VP infants, but similar moderate-to-severe NDI rates ( Table 4). In Appendix 1, classifications for the separate components of NDI are presented for each follow-up age.

NDI From 2 to 8 Years of Age
In the 554 children with three follow-up contacts, NDI could be calculated at all time points for 495 children. No NDI during the complete trajectory at 2, 5, and 8 years of age was seen for 179 (36%) children and both no-or-mild NDI during the complete trajectory was seen for 427 (81%) children. Moderateto-severe NDI during the complete trajectory was seen for 21 (4.2%) children. In these 495 children, from 2 to 8 years 314 (63%) children remained in the same NDI category, 101 (20%) children improved toward a better NDI category, and 80 (17%) children deteriorated toward a worse NDI category (Figure 2). Of all 293 children with normal NDI at 2 years, 223 (76%) remained in the normal NDI category at 8 years of age. For mild impaired infants 43% (66/152) and for moderate-to-severe impaired infants 50% (25/50) remained in the same NDI category. No differences were found in the characteristics of infants that remained in the same category, improved or deteriorated from 2 to 8 years ( Table 5).

Individual Longitudinal Trajectories in NDI
In clinical work individuals are more important than (sub)groups. Therefore Figure 3 presents the horizontal line plot for NDI at 2, 5, and 8 years of age, showing individual patterns

DISCUSSION
In this study, neurodevelopmental impairment at 2, 5, and 8 years was evaluated in very and extremely preterm children born below 30 weeks gestation. In addition the course of the individual longitudinal trajectories over time was studied. We observed individual variation over time in NDI status in 37% of the children, with 17% showing a change to a more worrisome category, but 20% showing an improvement. However, 63% of the children remained in the same category over time. Longitudinal analysis showed a clear association of gestation, gender, and parental education with the severity of NDI at all time points. No differences were found in the characteristics between children that improved and deteriorated, and the course of the trajectories in NDI was not affected by gestation, gender, and parental education. Compared to other studies we observed higher rates for children without or with mild NDI. At the age of 5 and 8 years, respectively, 54 and 62% of the surviving children showed a normal neurodevelopment, and 36 and 30% of the surviving children showed a mild neurodevelopmental impairment. In EP children, normal development rates were 47 and 55%, and mild NDI rates were 45 and 36% at 5 and 8 years, respectively. The Swedish EXPRESS study found rates of 36 and 30% for children without and with mild NDI at 6.5 years in children born below 27 weeks' GA (32). The EPICure study from the UK showed a rate of 75% for children with none-to-mild NDI at 6 years and a rate of 53% for children with none-to-mild NDI in 53% at 11 years, in children born below 26 weeks' GA (3,33). Unfortunately, international comparisons are hampered by differences in age of follow-up, definition of neurodevelopmental impairment and study population (21). For example, the EPICURE and EXPRESS studies included substantially more immature children, born at 22-24 weeks' gestation, whereas in our sample the youngest children were born at 25 weeks' gestation. FIGURE 2 | Shifts in NDI from 2 to 8 years of age. This figure shows NDI rate at 2 vs. 8 years of age for infants with NDI calculation at all three follow-up contacts (N = 495). The numbers are presented as N (%), with the % calculated relatively to the full group of N = 495 infants. The row sums show the total number of infants at 2 years of age for normal, mild, and moderate-to-severe NDI. The column sums show the total number of infants at 8 years of age for normal, mild, and moderate-to-severe NDI. The dark gray boxes represent all infants that deteriorated toward a worse NDI category from 2 to 8 years, the light gray boxed represent all infants that improved toward a better NDI category from 2 to 8 years, and the white boxes represent all infants that remained in the same NDI category.
Mild neurodevelopmental problems were seen in 31, 36, and 30% of the infants at 2, 5, and 8 years of age. However, mild deficits in multiple domains might be of the same severity as one moderate-to-severe deficit in a single domain. At 2, 5, and 8 years of age, 12.4, 25.3, and 18.6% of the infants with mild NDI had mild problems in more than one domain. Apparently, at a later age more multiple mild deficits become overt. Multiple deficits across domains may have combined long-term effects, which unfortunately is not reflected by the NDI definition. The significance of milder forms of neurocognitive deficits might need additional research (34).
The moderate-to-severe disability rate in the current study initially appeared to be 8.8% at 8 years of age. However, it was found that 13% of the children lost for follow-up at 8 years of age did not attend follow-up because they were already in treatment in rehabilitation medicine. Including these children as having moderate-to-severe disability resulted in a disability rate of 16%, which is slightly higher than the severe disability rate of 13% reported in both the EPICure and EXPRESS studies (3,32). On the other hand, children in our study were also lostto-follow-up because they did not experience any problems. The real moderate-to-severe disability rate probably is somewhere between 8.8 and 16%. This emphasizes the importance of presenting impairment rates in the context of reasons for lossto-follow-up. Despite the abundancy of cross-sectional follow-up studies a paucity exists in longitudinal follow-up. This study showed that approximately two third of the children assessed at 2 years of age were classified in the same NDI category at 8 years of age, and that 16% of all children became worse at 8 years of age. Similar results have been reported before in the EXPRESS study, reporting that 47% of all children remained in the similar NDI category and 32% of all children deteriorated toward a worse NDI category from 2 to 6.5 years of age (32). Although overall NDI rates remained comparable over time, these results demonstrate considerable individual variation over time. Indeed this also shows the importance of continuing follow-up until school life for individual and specific referrals and advise.
EP/VP status, gender, and parental education were found to be associated with severity of NDI at all time points. These results were in line with previously published studies, reporting genderdifferences in neurodevelopmental outcomes in the favor of girls (35)(36)(37)(38). Moreover, these results enhance the formerly reported association between gestational age and neurodevelopmental outcome as well as the association between parental education and neurodevelopmental outcome (11,36,(38)(39)(40).
Although EP/VP status, gender and parental education were found to be associated with NDI, these associations remained stable over time and the course of the trajectories was not affected by these factors. Children with these characteristics therefore seem to have the same developmental growth potential as children without these characteristics (18). Moreover, no differences were found in characteristics between infants that improved or deteriorated from 2 to 8 years. Nevertheless, considerable individual differences were seen in trajectories. This indicates the importance of other factors that might influence development over time, for example early childhood interventions such as an extensive physiotherapy program or special education assistance. Moreover, socio-environmental factors such as the quality of the parent-child relationship are important throughout development (16).
Although extensive evaluation of separate domains is important, the added value of a composite outcome is that it provides an overall impression of the outcomes after very preterm birth. Problems after preterm birth occur in a range of developmental domains and therefore it is important not to focus on single domains of development. Looking separately at the specific domains in this study, the majority of the children did not have any impairment. However, combining the different domains into the NDI composite outcome showed no NDI during the complete trajectory for a minority of 36% of all children. Apparently, the majority of the very preterm children do experience some clear problems at some time during childhood. Moreover, the combined outcome measure used in this study is the longer term outcome most frequently used for comparisons both within and between countries (21). International comparisons can guide clinical decision-making and provide prognostic information for families.
In this study, cognitive scores were corrected for prematurity at age 2, but not at age 5 and 8. The current Dutch national guideline on follow-up and most international guidelines recommends the use of corrected scores for preterm children up to 2-3 years. However, in very preterm children at age 5, a significant difference between corrected and uncorrected IQ was found, with corrected scores of course being higher than uncorrected scores (41). For future research, consistent reporting of cognitive outcome based on corrected scores is recommended (42).
The overall follow-up rate in this study was 79%, which is comparable to follow-up rates of other studies, showing rates varying from 71 to 92% at different ages (3,5,32). Moreover, more than 60% of the children completed follow-up at all time points during the longitudinal follow-up program, which demonstrates a high follow-up rate compared to other longitudinal studies such as the recently published EPICure2 study (follow-up rate 19%) (33). Our results might represent the worst-case scenario as medium risk children have not always been invited for follow-up during the study period because of limited resources as shown in Table 2. These children without follow-up were children with an appropriate birth weight, without severe brain injury and with uncomplicated NICU admission. This explains why infants seen for followup were more immature at birth and had an increased length of stay in the NICU compared to infants without followup. On the other hand, a significant difference was found in socio-economic status between infants with and without follow-up, with a higher SES score in the children that did have follow-up. This finding is similar to findings in previous studies, showing that drop-out was more likely to occur in families with social disadvantages, while preterm children from socially disadvantaged families may have poorer neurodevelopment (3,32,43). However, considering the high follow-up rate in this study, limited influence on the presented results is expected.

Strengths and Limitations
Strengths of our paper included the size of the cohort and the high follow-up rate with most children assessed at three ages. Moreover, the longitudinal nature of this study provides important information reading the developmental course of the children. However, this study has several limitations. First, different tests measuring cognitive performance had to be used, both at different ages to be developmentally appropriate, but also over time in order to use the most recent population norms. Use of different tests intending to measure the same constructs at different ages is inevitable when performing longterm longitudinal studies as development continues at high pace and differentiates strongly during infancy and toddlerhood as well as preschool age (5,19,44). In addition, tests need to be re-evaluated and updated over time to allow ecologically valid assessments (e.g., think of the appearance and use of phones in the nineties and zero's, causing the need for revision of the images used in cognitive test). As all tests were standardized however, with a mean of 100, results could be compared. Second, defining NDI based on four determinants (cognitive, neurological, auditory, and visual function) has its limitations to delineate a child's development. Additional domains such as behavioral problems could not be taken into account but may also impair children over time. Third, ideally the GMFCS classification would have been used for classifying the severity of problems in the neurological domain. However, this system was not routinely used in 1990. In order to distinguish between moderate and severe neurological problems, uni-and bi-lateral paresis was used as a proxy for GMFCS 1-2 and GMFCS 3-5, respectively. Last, in this retrospective study, no specific information on interventions was available. Improvement during the trajectories could potentially be the result of adequate interventions after detection of NDI at early age, resulting in improved NDI at early age. Future research may elaborate on the effect of interventions on the individual trajectories.
In conclusion, this study evaluated neurodevelopmental impairment at three different ages up to the age of 8 in very preterm children, next to the course of the longitudinal trajectories in these outcomes. A clear association was found of gestation, gender, and parental education with the severity of NDI at all time points. Although we observed considerable individual variation over time in NDI status, the course of the trajectories in NDI were not associated with gestation, gender, and parental education. These results point to the importance of other (unknown) influences on developmental trajectories. Continued follow-up until school life for extremely preterm born children is essential in order to provide optimal individually focused referrals and care when needed.

DATA AVAILABILITY STATEMENT
The datasets presented in this article are not readily available because of privacy regulations according to the General Data Protection Regulation (GDPR). Requests to access the datasets should be directed to Pauline E. van Beek, pauline.van.beek@mmc.nl.

ETHICS STATEMENT
The studies involving human participants were reviewed and approved by Medical Ethics Committee Máxima MC. Written informed consent from the participants' legal guardian/next of kin was not required to participate in this study in accordance with the national legislation and the institutional requirements.
psychologist Titia Katgert, for their contribution in examining the children's neurodevelopment at (pre)school age. The authors thank Anne Verheijen, Guusje Thijssen, Anne van Och, René Blom, and Jasmijn van Erp for their help in creating the database for research use.