Novel diffuse white matter abnormality biomarker at term-equivalent age enhances prediction of long-term motor development in very preterm children

Our objective was to evaluate the independent prognostic value of a novel MRI biomarker—objectively diagnosed diffuse white matter abnormality volume (DWMA; diffuse excessive high signal intensity)—for prediction of motor outcomes in very preterm infants. We prospectively enrolled a geographically-based cohort of very preterm infants without severe brain injury and born before 32 weeks gestational age. Structural brain MRI was obtained at term-equivalent age and DWMA volume was objectively quantified using a published validated algorithm. These results were compared with visually classified DWMA. We used multivariable linear regression to assess the value of DWMA volume, independent of known predictors, to predict motor development as assessed using the Bayley Scales of Infant & Toddler Development, Third Edition at 3 years of age. The mean (SD) gestational age of the cohort was 28.3 (2.4) weeks. In multivariable analyses, controlling for gestational age, sex, and abnormality on structural MRI, DWMA volume was an independent prognostic biomarker of Bayley Motor scores (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}β= −12.59 [95% CI −18.70, −6.48] R2 = 0.41). Conversely, visually classified DWMA was not predictive of motor development. In conclusion, objectively quantified DWMA is an independent prognostic biomarker of long-term motor development in very preterm infants and warrants further study.

quantitative cerebral morphometric analyses representing the most clinically feasible approach. Of these morphometric measures, objective assessment of diffuse excessive high signal intensity (DEHSI) abnormality is of import due to its high prevalence in preterm infants 14,[16][17][18][19] , association with in-vivo microstructural 18,20 , metabolic 21 , and postmortem pathology 22 , and early evidence suggesting a correlation with neurodevelopmental impairments (NDI) 17,[23][24][25][26] . While most studies that diagnosed DEHSI visually/qualitatively have not reported a significant association with NDI 18,19,[27][28][29][30][31][32] , when quantified objectively, DEHSI appears to significantly predict cognitive and language development in extremely preterm infants 24,25 . To better reflect its pathologic nature, we will henceforth use the label diffuse white matter abnormality [DWMA] in place of DEHSI. The goal of this study was to examine the prognostic value of objectively quantified DWMA volume at term-equivalent age for prediction of motor development in a prospective cohort of very preterm infants. We hypothesized that objectively quantified DWMA volume would be an independent predictor of motor development at 3 years of age.

Methods
Population. All very preterm infants born at 31 weeks completed gestation or earlier and admitted to any of the four level III neonatal intensive care units (NICUs) in Columbus, Ohio from November 2014 to March 2016 were eligible for inclusion 25 . We prospectively enrolled 110 very preterm infants from a consecutively eligible sample of infants during this period. The four NICUs were Nationwide Children's Hospital (NCH), Ohio State University Medical Center, Riverside Hospital, and Mount Carmel St. Ann's Hospital. These NICUs care for approximately 80% of all very preterm infants in the Columbus, Ohio region. We excluded any infants with congenital or chromosomal anomalies that affected the central nervous system and likely result in a poor outcome. Data collection occurred between January 2015 and July 2018. The NCH Institutional Review Board approved the study at NCH and the other study sites through established reciprocity agreements. Written informed consent was obtained from a parent or guardian of each very preterm infant after they were given sufficient time to determine if they wished to participate. All methods/research activities were carried out in accordance with the NCH Institutional Review Board guidelines and regulations. All study infants were invited for routine developmental follow-up in the NCH High-Risk Follow-up Clinic up to 3 years corrected age.
Magnetic resonance imaging acquisition. We performed brain structural MRI scans on all 110 study infants at NCH on a 3T Siemens Skyra MRI scanner with at 32-channel pediatric head coil between the ages of 39 and 44 weeks post-menstrual age (PMA). Most infants from NCH were typically imaged while they were still inpatients, while all infants cared for at the other three NICUs were imaged as outpatients after being discharged. All inpatient MRI scans were attended by a skilled neonatal nurse and a neonatologist. Heart rate and oxygen saturation of all infants were monitored continuously during all scans. We performed all imaging without sedation by feeding the infants 30 min prior to the scan, applying silicone earplugs and swaddling the infants in a blanket and a vacuum immobilization device (MedVac, CFI Medical Solutions, Fenton, MI) to promote natural sleep. There were no adverse events. The following structural MRI sequence parameters were used for all infants: axial T2-weighted: echo time 147, repetition time 9,500 ms, echo train length 16 Image post-processing. We applied our previously published algorithm to objectively detect and quantify DWMA on T2-weighted MRI ( Fig. 1; For in-depth methods and additional examples of DWMA segmentation see He et al.) 25 . To summarize, first we conducted bias field correction (removal of signal intensity inhomogeneity caused mainly by the radiofrequency coils) and intensity normalization (reducing the variations in signal intensity and contrast across slices and across subjects). Next, we conducted brain tissue segmentation using a neonatal probabilistic brain atlas as a guide and defined DWMA to be any voxels with signal intensity values greater than ∝ standard deviations above the mean for all cerebral tissues (white and gray matter). We refer to ∝ as our cut-off threshold. For this study, we examined a cut-off thresholds of 2.0. However, this threshold was too restrictive and defined only very small regions as DWMA; therefore, we chose a lower threshold of 1.8 SD. We controlled for partial volume artifacts by only labeling voxels with high gray and white matter membership probability (≥ 95%) as cerebral tissues. We manually removed the few isolated false positive voxels detected by the algorithm. Total DWMA volume was calculated as the product of a single voxel volume (determined by the imaging resolution) and the total number of voxels in the detected DWMA region. This algorithm was first validated on simulated preterm infant brains with manually drawn DWMA that represented the ground truth by demonstrating that our DWMA algorithm exhibited strong agreement with this ground truth (both qualitatively and quantitatively) 25 . We limited DWMA detection to the centrum semiovale only because we have found this to be the most predictive white matter region and it is not confounded by the normal high signal intensity of the periventricular crossroads [24][25][26] . We defined the centrum semiovale as the central white matter in the two slices immediately above the lateral ventricles on axial view. We calculated a normalized DWMA volume by dividing DWMA volume by total cerebral white matter volume. All analyses were performed masked to clinical and outcome data. MRI imaging assessment. All brain structural MRI readings were performed by pediatric neuroradiologists who used a standardized scoring system graded for degree of brain injury/maturation, and the objective quantitative biometric measurements were performed separately by a trained expert, per Kidokoro et al. 33 . This Scientific RepoRtS | (2020) 10:15920 | https://doi.org/10.1038/s41598-020-72632-0 www.nature.com/scientificreports/ approached yielded a global brain abnormality score, which was categorized as normal (total score, 0-3), mild (total score, 4-7), moderate (total score, [8][9][10][11], or severe abnormality (total score ≥ 12). A single reader (NAP) with greater than 10 years of experience interpreting neonatal MRI scans performed visual qualitative classification of DWMA, masked to clinical and outcome data. The DWMA score was based on severity and extent as described by Kidokoro et al. 18 . Infants were assigned grade 0 if there was no DWMA or if high signal intensity was present only in the periventricular crossroads, grade 1 if DMWA was only visible in one region, grade 2 if DWMA was visible in two regions, and grade 3 if three or more regions were involved in addition to the normal signal intensity observed in the crossroads. While DWMA was observed in all white matter regions, the centrum semiovale was the most commonly identified region with qualitatively defined DWMA. The reader also assessed whether the margins of the posterior crossroads were invisible (defined as invisible posterior crossroads). The same reader reevaluated 20 randomly chosen MRI scans three weeks later and used kappa ( κ ) statistics to assess intra-rater agreement for DWMA grade. Of the 20 subjects, complete agreement was seen in 60% of cases (expected agreement 31.0%) for a κ of 0.42. This represents a fair to moderate agreement strength 34,35 . Neurodevelopmental assessment. Participating infants underwent a comprehensive neurodevelopmental evaluation at a median age of 36.1 (IQR: 35.3-37.5) months in the NCH High-Risk Follow-up Clinic. We assessed overall motor development using the standardized Bayley Scales of Infant and Toddler Development, Third Edition (Bayley-III). A Motor composite score (composite of fine and gross motor development scores) that was 3 SD below the normative mean was assigned to children who could not complete the test due to difficulty resulting from likely severe disability. The composite score for the Bayley-III Motor subscale is scaled to metric with a mean of 100 (SD 15) and range of 40-160. Examiners performed the standardized Amiel-Tison neurologic exam 36 , which included evaluation of tone, reflexes, posture, and strength; gross motor function was classified using the Gross Motor Function Classification System 37 . Cerebral palsy was defined as abnormal muscle tone in at least one extremity and abnormal control of movement and posture. All assessments were www.nature.com/scientificreports/ performed by assessors who were masked to the quantitative DWMA diagnosis but not masked to clinical information.

Statistical analyses.
In univariate analyses, we examined the relationship between the normalized DWMA volume and the Bayley-III Motor composite score using linear regression. To evaluate the independent prognostic value of DWMA volume, we performed multivariable regression by adding known perinatal predictors of Bayley score, including sex, gestational age, and global brain abnormality score. In addition, we also added center/NICU and PMA at MRI to the multivariable model to control for their potential confounding effects.
We also tested a model that substituted global brain abnormality with a composite variable that included sMRI injury variables known to be strong predictors of motor impairment: cystic white matter abnormalities, hemorrhage (intraventricular, parenchymal, and/or cerebellar), and punctate white matter lesions 7 . The internal validity of our final model was tested by estimating a bias-corrected confidence interval derived from a bootstrap procedure involving 10,000 resamples 38 . In secondary analyses, to assess prediction accuracy for CP, we used Fisher's exact test to evaluate prognostic properties, including sensitivity, specificity, positive and negative likelihood ratios, for: (1) objectively quantified severe DWMA (normalized DWMA volume dichotomized at > 90th percentile [pre-specified cut-off]), (2) global brain abnormality (moderate or greater), and (3) visually-classified severe DWMA (grade 3). We used logistic regression to evaluate the relationship between DWMA and CP. Last, we used Pearson's correlation and multivariable linear regression to assess the relationship between (1) DWMA volume and global brain abnormality score and (2) DWMA volume and visually defined DWMA. We used the traditional two-sided P value of < 0.05 to indicate statistical significance. All analyses were performed using STATA 16.0 (Stata Corp., College Station, TX).

Results
Of the original cohort of 110 very preterm infants, we excluded one infant due to excessive motion artifacts and excluded all 11 infants with severe brain injury since this interfered with accurate DWMA segmentation (e.g. severe ventriculomegaly resulting in loss of centrum semiovale white matter). Structural MRI was performed at a mean (SD) PMA of 40.3 (0.5) weeks. By 3 years of age, 77 infants (79%) returned for Bayley motor testing. The baseline characteristics for infants who returned for follow-up were not significantly different from those who did not ( Table 1). The mean (SD) Bayley-III raw Gross Motor, raw Fine Motor, and Composite Motor scores were 60.1 (5.6), 44.0 (5.0), and 90.8 (12.2), respectively. Cerebral palsy was diagnosed in six infants (7.3%). Four infants were diagnosed with spastic diplegia, one had spastic left hemiplegia, and one had spastic quadriplegia. The latter infant was classified as GMFCS level 4 while the other five infants were classified as GMFCS level 1. Based on the global brain abnormality score, five infants were classified as having moderate injury (6.5%), 19 had mild injury (24.7%), and 53 had no injury (68.8%) on their sMRI at term-equivalent age. Moderate injury was noted on sMRI in three of the six infants (50%) with CP. As stated above, all infants with severe injury were excluded from the study. Visually/qualitatively classified DWMA was diagnosed as severe (grade 3) in 11 infants (13.4%), moderate (grade 2) in 22 infants (26.8%), and no/mild (grade 0/1) in 49 infants (59.8%). Only one infant was diagnosed with invisible posterior crossroads.
In univariate analyses, DWMA volume was significantly predictive of Bayley-III Motor scores, explaining 26% of the variance in motor development ( Table 2; Fig. 2). This association remained significant even when raw DWMA volume was tested in the regression analyses, suggesting that the normalization by total white matter volume did not have a significant effect on the association with Bayley Motor scores. In multivariable analyses, controlling for other known predictors of Bayley scores, including sex, gestational age, and global brain abnormality, normalized DWMA volume ( β = −12.59 [95% CI −18.70, −6.48]) remained a significant predictor of  Table 2). Replacing global brain abnormality score with cystic abnormalities, hemorrhage, and punctate white matter lesion variables actually reduced the model adjusted R 2 (38.9%) and enhanced the predictive power of DWMA ( β = −14.33). The bootstrap bias-corrected confidence intervals were comparable ( β 95% CI −18.60, −4.31), supporting the internal validity of the final model. To confirm that the significant predictive relationship we observed between normalized volume of DWMA and Motor scores (t = −4.11; p < 0.001) was not a function of white matter volume loss, we replaced normalized DWMA volume with DWMA volume that was uncorrected for white matter volume. This replacement did not result in a meaningful difference in the model parameters (t = −3.92; p < 0.001). When we control for the effects of different head sizes by including total white matter in the multivariable model (with uncorrected DWMA volume), the model remains very comparable (t = −3.89; p < 0.001); when controlled for total intracranial volume, the results again remain comparable (t = −3.95; p < 0.001). The total explained variance in Bayley Motor scores was 39% and 40%, respectively (compared to 41% for normalized DWMA volume). These analyses suggest that DWMA is independent of other white matter pathology.
Visual, qualitative diagnosis of DWMA was not significantly predictive of Motor scores in univariate analyses (P = 0.23). Inclusion of known predictors and confounders in the model did not substantially change this  Table 2). Global brain abnormality score was also a significant predictor of Motor scores, even after controlling for other predictors. However, it was not as predictive as DWMA volume. The addition of DWMA volume increased the variance in Bayley Motor scores by another 13.2% (41.1 vs. 27.9%; Table 2). Finally, normalized DWMA volume was significantly correlated with global brain abnormality score (r = 0.30; p = 0.003) but not with visually defined qualitative DWMA (r = 0.09; p = 0.23) in univariate analyses. In multivariable linear regression models, controlling for gestational age, sex, PMA at MRI, and center, we observed a significant relationship between objectively defined DWMA volume and global brain abnormality score ( β =0.032; 95% CI 0.003, 0.062; p = 0.032). In similarly controlled multivariable linear regression models, we did not observe a significant relationship between objectively defined DWMA volume and visually defined DWMA ( β=0.044; 95% CI -0.064, 0.152; p = 0.418).
In secondary logistic regression analyses, a 10% increase in DWMA volume was associated with an odds ratio of 31.64 (95% CI 3.96, 253.03) of developing CP (p < 0.001). A one point increase in global brain abnormality was associated with an odds ratio of 2.25 (95% CI 1.44, 3.51) of developing CP (p < 0.001). The relationship between DWMA volume and CP remained significant in a multivariable model after controlling for global brain abnormality, gestational age, and center (OR 12.64, 95% CI 1.41, 113.43). Conversely, we did not find a significant relationship between qualitatively diagnosed DWMA and CP (OR 2.29; 95% CI 0.81, 6.46; p = 0.118). Objectively quantified severe DWMA (P < 0.001) and global brain abnormality on structural MRI (P = 0.004) were both significantly predictive of CP, while visually-classified severe DWMA (grade 3) was a poor predictor of CP (P = 1.00) ( Table 3).

Discussion
We demonstrated that objectively quantified DWMA is a significant and independent prognostic biomarker of motor development at 3 years of age in very preterm infants. Normalized DWMA volume remained a prominent predictor of standardized Motor scores even after controlling for other known predictors of motor development such as sMRI abnormalities, gestational age, and sex. This is notable because we excluded all infants with severe brain injury, which is the most prominent risk factor for the development of CP and motor impairments. Of the six infants that developed CP in this cohort, five had mild CP. Two of the three infants with a normal sMRI developed spastic diplegia. A recent CP registry study found that infants with these CP subtypes were twice as likely to have normal sMRI scans 39 . Our results suggest that DWMA is pathologic and deserves further testing in larger studies to externally validate its prognostic value for the early detection of minor motor impairments and CP.
In two previous independent cohorts, we have shown that DWMA volume quantified using our objective algorithm is predictive of cognitive and language development at 2 years 25,26 . For motor development, it is difficult to compare our objective DWMA biomarker results with other studies because no prior study has attempted to predict Motor outcomes using objectively quantified DWMA. However, similar to our findings, multiple studies have examined the link between visually classified DWMA and motor development and found no correlation [17][18][19]23,[27][28][29][30]32 . There are likely several confounding factors that contributed to a lack of association. First, visually diagnosed DWMA is subjective and exhibits suboptimal retest reliability, as demonstrated in our study and several prior studies 14,31,40,41 . This may be due to the signal inhomogeneity and the occurrence of developmental crossroads that are routinely present on all MRI scans at term-equivalent age. This is especially true for periventricular white matter regions, which potentially explains why our previous study showed lower predictive value for this white matter region as compared to the centrum semiovale 24 . Therefore, we limited assessment of DWMA to only the centrum semiovale white matter. Lower diagnostic reliability increases measurement error, which can reduce the likelihood of finding a significant association 42 . Our intra-rater reliability for visual diagnosis of DWMA was only fair to moderate ( κ 0.42). This reliability was comparable to our prior study where a pediatric neuroradiologist diagnosed DWMA ( κ 0.46) 14 . Second, for a given sample size, a qualitative diagnosis is categorical and will therefore exhibit lower study power than a quantitative diagnosis (continuous measure) 43,44 . Lastly, even when qualitative DWMA diagnostic reliability is excellent 18 , DWMA diagnosis may still be inaccurate because there is no gold standard test to confirm true DWMA pathology.
Our DWMA segmentation algorithm is automated, easy to use, and can generate DWMA and whole-brain tissue volumes within 5 min. For about half of the MRI scans, it requires no further manual correction and the other half it incorrectly labels between 2 to 8 voxels, most commonly in the interhemispheric fissure or peripheral white matter. While we currently remove these mislabeled voxels manually, we have recently developed a fully-automated approach using machine learning to overcome this limitation 45 . Such a tool could be used Table 3. Prognostic test properties for objectively-diagnosed severe DWMA, abnormal structural MRI, and visually-classified severe diffuse white matter abnormality (DWMA) for predicting cerebral palsy in very preterm infants. Statistically significant results are given in italics (P < 0.05)

Predictors
Sensitivity (95% CI) Specificity (95% CI) Although it's predictive accuracy has been touted as approximately 90% 7 , when more robust evidence is considered, its sensitivity is closer to 70% and its predictive probability is only 35% [11][12][13][14]46 . This leaves a substantial gap in our ability to accurately detect CP at term-equivalent age in order to take advantage of the early critical window for neuroplasticity in the first two years. This is the period during which proven (re)habilitative interventions could restore motor function and thus improve quality of life. Prognostic tests such as the general movements assessment and the Hammersmith Infant Neurological Exam are being increasingly performed in many centers at 3 months corrected age 7 , however more research is still needed to determine their accuracy in combination with injury on sMRI in predicting CP and minor motor abnormalities 47 . Other advanced quantitative MRI measures such as diffusion MRI and structural and functional sensorimotor tract connectivity 48 could potentially fill this gap, as highlighted in a recent systematic review 15 . For example, several studies have reported abnormal microstructural properties of DWMA using diffusion MRI 18,20 . However, when these measures were compared to quantitative DWMA volumes determined via signal intensity, the addition of microstructural measures did not improve outcome prediction 24 . The incremental predictive value of other promising advanced MRI biomarkers, independent of sMRI, remains to be validated in larger, population-based cohort studies.
Our study has several limitations. Our follow-up rate was 79%, which may have introduced ascertainment bias. However, the baseline characteristics of infants with and without follow-up testing were comparable, suggesting a low risk of bias. Only six infants developed CP and thus our secondary CP prognostic analyses will need to be validated in larger studies. Also, we were unable to determine the predictive ability of DWMA volume over and above general movements assessment or early standardized motor exam because these tests were not part of the research study or clinical care during study enrollment. This limitation is being addressed in our newer and larger cohort study. Our Bayley assessors were blinded to DWMA result but not blinded to clinical or structural MRI information; thus, this may have biased our findings. However, this bias would reduce the independent association of DWMA volume with Bayley Motor scores by potentially strengthening the association of clinical and global brain abnormality scores with Motor scores. Strengths of our study include a geographically-based cohort, objective quantification of DWMA on sMRI that can be readily translated clinically, comparison with the current standard method of visual classification, standardized assessments of motor outcomes up to 3 years of age when minor motor abnormalities are more evident, and independent validation of DWMA volume as a new prognostic biomarker, over and above existing predictors.
In conclusion, in this multicenter prospective cohort study, we were able to demonstrate for the first time that objectively quantified DWMA is an independent predictor of motor development in very preterm infants. We also validated prior research showing that visually classified DWMA is not predictive of neurodevelopmental outcomes and is therefore suboptimal for use in clinical practice. Additional studies are needed to externally validate the use of DWMA volume as an early prognostic biomarker for cerebral palsy and minor motor impairments and to enable clinical translation of our DWMA algorithm. If externally validated, our findings could be applied to improve risk stratification at hospital discharge for targeted, aggressive early intervention therapies.

Data availability
All data, software, and code from this study are being submitted for publication and can be accessed from the lead author in the meantime.