Identifying relevant biomarkers of brain injury from structural MRI: Validation using automated approaches in children with unilateral cerebral palsy

Abstract

Previous studies have proposed that the early elucidation of brain injury from structural Magnetic Resonance Images (sMRI) is critical for the clinical assessment of children with cerebral palsy (CP). Although distinct aetiologies, including cortical maldevelopments, white and grey matter lesions and ventricular enlargement, have been categorised, these injuries are commonly only assessed in a qualitative fashion. As a result, sMRI remains relatively underexploited for clinical assessments, despite its widespread use. In this study, several automated and validated techniques to automatically quantify these three classes of injury were generated in a large cohort of children (n = 139) aged 5–17, including 95 children diagnosed with unilateral CP. Using a feature selection approach on a training data set (n = 97) to find severity of injury biomarkers predictive of clinical function (motor, cognitive, communicative and visual function), cortical shape and regional lesion burden were most often chosen associated with clinical function. Validating the best models on the unseen test data (n = 42), correlation values ranged between 0.545 and 0.795 (p<0.008), indicating significant associations with clinical function. The measured prevalence of injury, including ventricular enlargement (70%), white and grey matter lesions (55%) and cortical malformations (30%), were similar to the prevalence observed in other cohorts of children with unilateral CP. These findings support the early characterisation of injury from sMRI into previously defined aetiologies as part of standard clinical assessment. Furthermore, the strong and significant association between quantifications of injury observed on structural MRI and multiple clinical scores accord with empirically established structure-function relationships.

About the Authors:

Alex M. Pagnozzi

Roles Conceptualization, Data curation, Formal analysis, Methodology, Software, Validation, Visualization, Writing – original draft

* E-mail: [email protected]

Affiliations CSIRO Health and Biosecurity, The Australian e-Health Research Centre, Brisbane, Australia, The School of Information Technology and Electrical Engineering, The University of Queensland, Brisbane, Australia

ORCID http://orcid.org/0000-0001-7931-1204

Nicholas Dowson

Roles Formal analysis, Methodology, Project administration, Resources, Software, Supervision, Validation, Writing – review & editing

Affiliation: CSIRO Health and Biosecurity, The Australian e-Health Research Centre, Brisbane, Australia

James Doecke

Roles Formal analysis, Methodology, Validation, Writing – review & editing

Affiliation: CSIRO Health and Biosecurity, The Australian e-Health Research Centre, Brisbane, Australia

Simona Fiori

Roles Conceptualization, Writing – review & editing

Affiliation: Stella Maris Scientific Institute, Pisa, Italy

Andrew P. Bradley

Roles Formal analysis, Methodology, Software, Writing – review & editing

Affiliation: The School of Information Technology and Electrical Engineering, The University of Queensland, Brisbane, Australia

Roslyn N. Boyd

Roles Conceptualization, Investigation, Writing – review & editing

Affiliation: Queensland Cerebral Palsy and Rehabilitation Research Centre, School of Medicine and Science, Centre for Children’s Health Research, The University of Queensland, Brisbane, Australia

Stephen Rose

Roles Conceptualization, Funding acquisition, Methodology, Resources, Supervision, Writing – review & editing

Affiliation: CSIRO Health and Biosecurity, The Australian e-Health Research Centre, Brisbane, Australia

Introduction

Cerebral palsy (CP) is an umbrella term that covers a range of injury occurring in the developing brain, pre- or perinatally [1]. Although impairments to motor function are the most characteristic symptom of CP [2], impairments to cognition, vision and communication have also been observed as a result of these injuries [3,4]. Magnetic resonance imaging (MRI) is strongly recommended to elucidate the presumed timing and aetiology of brain insults that cause CP [5] to help facilitate therapeutic selection [6], and is the current standard for assessing injury in clinical practice. Current characterisation of injuries related to CP are broadly grouped into three classes based on aetiological patterns [7]: brain maldevelopments (of which cortical malformations are the main type) occurring from disturbances in the first and second trimesters; periventricular white matter injury (leading to white matter lesions, and potentially the secondary ventricular enlargement due to primary tissue loss) occurring from disturbances in the early third trimester; and cortical/deep grey matter injury occurring from disturbances in the late third trimester (leading to grey matter lesions). The use of these characterisations are common, and many studies have investigated the prevalence of these timing-related aetiologies of injury within specific cohorts of children with CP [7–9]. However, these classifications have only been qualitatively assessed, disregarding the location and severity of lesions observed from the structural MRI (sMRI).

Quantitative measures of injury, accounting for injury severity and anatomical location, have the potential to better quantify the relationship between brain lesions and functional outcomes [10]. Furthermore, in future, such characterisations of injury may allow for valid and reliable predictions of patient impairment from MRIs acquired very early in life, which has implications for the selection of treatment strategies. However, this requires lesion volumes to be segmented in three dimensional images, which is too laborious to perform manually on large cohorts of data. Automated techniques are necessary to perform these segmentations in a repeatable and time-efficient manner. Such techniques need to be tailored to the specific challenges present in the sMRI of children with CP, including the potentially severe morphological alterations and the heterogeneous appearance of lesions, which are illustrated in Fig 1. However, no study has utilised the comprehensive quantifications of injury that automated approaches allow, nor have they examined the utility of using these quantifications to predict clinical function. As a result, it remains unclear whether quantitative measures of brain injury can usefully augment the current classifications of brain injury in the clinical setting to produce finer estimates of patient function. To this end, this study aims to quantitatively characterise the prevalence and extent of different manifestations of injury in a cohort of 139 children; including 95 children diagnosed with unilateral CP and 44 children with healthy development (CHD), using three previously developed automated approaches. These approaches identify the three main aetiologies of injury based on current classifications [11]; including a tailored segmentation strategy and cortical shape analysis pipeline to detect cortical malformations [12], a lesion-as-outlier segmentation strategy using T1- and T2-weighted MRIs for the detection of white and grey matter lesions [13], and a statistical shape model (SSM) of healthy ventricular shape to detect the secondary enlargement of ventricles [14]. The quantitative biomarkers derived from these approaches were used to characterise the prevalence of injury in this cohort of children with CP, and to construct regression models, in order to identify those biomarkers that are significantly associated with clinical scores of motor, cognitive, communicative and visual function. Such models can improve our understanding of the relationship between the extent and topography of brain lesions and clinical function in children with CP, and can provide valid and reliable estimates of patient function, which can help to tailor treatment strategies, potentially leading to improved gains in functional outcomes for children with CP.

[Figure omitted. See PDF.]

Fig 1. Illustrations of severe injury (top row) and subtle injury (bottom row) in cases of cortical malformations (first column), white/grey matter lesions (middle column) and ventricular enlargement (last column).

(A) WM atrophy with corresponding ventricular enlargement is presented with abnormal sulcal depth predominantly in the frontal and occipital lobes. (B) WM lesions resulting from periventricular leukomalacia are shown as local regions of high intensity. (C) The consequences of a periventricular haemorrhagic infarction leading to a severe loss of WM and secondary enlargement of the lateral ventricles is shown, particularly on the left hemisphere. (D) Bilateral perisylvian polymicrogyria is shown, which are visible as excessive numbers of small gyri. (E) Gliosis in a small location in the posterior limb of the internal capsule is shown. (F) Periventricular cystic lesion in the right ventricular cella media is shown, leading to a small region of ventricular enlargement on the lateral side of the ventricle.

https://doi.org/10.1371/journal.pone.0181605.g001

Methods

Data

A total of 139 patients were included in this study: 95 patients with clinically diagnosed unilateral CP (50 male, 44 female, mean age 11.4, age range 5–17), and 44 CHD children (15 male, 29 female, mean age 10.7, age range 7–16) were included. Imaging data was made available from the Queensland Cerebral Palsy and Rehabilitation Research Centre (QCPRRC), as well as from the Stella Maris Institute, Pisa. Study participants included children who were recruited as part of ongoing studies of children with CP [15,16]. Diagnoses of CP were made based on clinical assessment by experienced clinicians in the field of CP. For both studies, ethical approval was granted and informed parental consent was obtained for all participants.

Imaging acquisition

All 139 patients underwent T1 Magnetization Prepared Rapid Gradient Echo (MPRAGE) scanning from one of two different scanners, and one of three different scanning parameters, including a 3T Siemens’ scanner with scanning parameters (TR = 1900 ms, TE = 2.32 ms, flip angle = 9 degrees, slice thickness = 0.9 mm), and a 1.5T GE scanner with two different scanning parameters (TR = 12.36 ms, TE = 5.17 ms, flip angle = 13 degrees, slice thickness = 1 mm) and (TR = 124.29 ms, TE = 4.37 ms, flip angle = 10 degrees, slice thickness = 1 mm). The developed image processing algorithm is robust to imaging sequences, and has been validated previously on this dataset [17]. The influence on the Expectation Maximisation (EM)/ Markov Random Field (MRF) algorithm is minimal as it adaptively models the intensity distribution of each tissue class.

A subset of patients (n = 125) also underwent either T2 Turbo Inversion Recovery Magnitude (TIRM) (TR = 7000 ms, TE = 79 ms, flip angle = 120 degrees, slice thickness = 4 mm) or T2 Half-Fourier Acquisition Single-Shot Turbo Spin-Echo (HASTE) (TR = 1500 ms, TE = 81 ms, flip angle = 150 degrees, slice thickness = 4 mm), acquired using a 3T Siemens’ scanner.

Clinical scores

Motor and cognitive impairments are the most common forms of disabilities in children with CP [3], although accompanying impairments to speech and vision have also been shown in up to 70% of patients [4]. Therefore, measures of patient motor, cognitive, visual and communicative function are all considered important in clinical assessment. As part of the ongoing studies on children with unilateral CP [15,16], six available clinical scores covering this potential range of functional impairments were Used in this study. The Assisting Hand Assessment (AHA) is a reliable measure of how well the assisting hand is used by the patient in bimanual tasks [18], and in this study is used to quantify the manual capabilities of children with CP. Measures of patient cognition, although difficult to quantify, is measured in this study using two parent reported questionnaires scoring their child’s behavioural and emotional function in daily life; the Behaviour Rating Inventory of Executive Function (BRIEF) [19] and Strengths and Difficulties Questionnaire (SDQ) [20]. To quantify visual acuity in this study, the Test of Visual Perception Skills (TVPS) measure [21] was utilised to assess patients’ ability to discriminate and memorise visual cues. Communicative ability was quantified in this study using the vocabulary (VOC) and word reasoning (WR) subtests of the Wechsler Preschool and Primary Scale of Intelligence (WPPSI-III) [22], which assess the patients’ ability to express and comprehend language, respectively.

Demographic information

The demographics of the children with CP and the CHD children in this cohort are provided in Table 1 below. There were slightly less males compared with females in the CHD group (p = 0.042), and no difference in the ages between CHD and CP groups (p = 0.548).

[Figure omitted. See PDF.]

Table 1. Demographic characteristics for the CHD and CP cohorts.

For the CHD cohort, their clinical scores and injury severity score were not obtained.

https://doi.org/10.1371/journal.pone.0181605.t001

Image pre-processing

Several image pre-processing steps were performed on the T1-MPRAGE, T2-TIRM or T2-HASTE images prior to the extraction of biomarkers of injury. Firstly, image alignment to the Colin 27 Average Brain Atlas [24] using an affine block matching registration algorithm [25]. Image bias was corrected for using the N4 algorithm [26]. Intensity normalisation and image de-noising, using anisotropic diffusion [27] with modified curvature diffusion equation [28], was performed with the Insight Toolkit (ITK) in order to minimise the effect of Rician-distributed noise in the MR image, while attempting to preserve high resolution features within the image. Skull stripping of the T1-MPRAGE, T2-TIRM and T2-HASTE images was performed using an in-house algorithm developed in MATLAB (Mathworks, Natick, MA). In this approach, intradural CSF was identified using thresholding and morphological operations, following which brain tissues were isolated. By identifying CSF internal to the skull boundary, this approach is capable of accurately segmenting the brain in cases of large lesions, which may be present in children with CP. Tissue probability maps (TPMs) from the Colin 27 Atlas were registered to the T1-weighted MRIs of each patient using the fast free-form deformation registration algorithm [29] to assist the downstream lesion segmentation.

Image biomarkers

Image-processing techniques were used to identify and quantify three types of injury observed in the MRIs of children with CP [11]. These methods and their performance are detailed in the following sections. An overview of the pipeline of these automated methods are illustrated in Fig 2.

[Figure omitted. See PDF.]

Fig 2. Illustration of the overall processing pipeline.

The pre-processing steps are shown in red, registration steps in yellow, the utilised atlases in green and blue, and the brain masking and tissue segmentation approaches in orange which facilitate the detection of ventricular enlargement (shown in light green), cortical shape measures (in brown) and lesion burden (in purple).

https://doi.org/10.1371/journal.pone.0181605.g002

Identifying cortical malformations

Primary (i.e. brain maldevelopments) or secondary (i.e. due to hypoxic insult) cortical abnormalities, which appear as a heterogeneous range of altered cortical shapes, were identified in this study using a brain tissue segmentation algorithm that accommodates severe injury [17], and three shape measures computed from the cortical grey matter segmentation (cortical thickness, curvature, sulcal depth) in order to quantify shape abnormalities [12] in each cortical region defined by the Automated Anatomical Labelling (AAL) atlas. An illustration of these cortical measures in three patients with observed cortical alterations are shown in Fig 3. Note for these figures, the cortical measures, which were obtained only on the cortical surface, were smoothed across the entire cortical grey matter segmentation. Biomarkers obtained from this approach are an absolute z-score from healthy cortical shape measures measured from the corresponding cortical region among the 44 CHD children. Unlike an absolute measure of cortical shape (such as cortical thickness in millimetres), a z-score relative to healthy cortical shape aligns with the subsequent volumes of lesion burden and ventricular enlargement volume, where larger values represent greater injury severity. This also facilitates the interpretation of the regression model coefficients.

[Figure omitted. See PDF.]

Fig 3. An illustration of three cortical abnormalities.

The first column shows the measured cortical thickness (in mm), the second columns shows sulcal depth (in mm) and the third column shows the curvature (in mm^-2) of these cortices. Regions of injury were observed to have higher sulcal depths, and reduced cortical thickness and curvature.

https://doi.org/10.1371/journal.pone.0181605.g003

Identifying white and grey matter lesions

White and grey matter lesions were identified using a lesion segmentation algorithm tailored to CP data, and with multiple lesion classes to identify both types of lesions [13]. Biomarkers from this approach are a lesion volume (in mL) in the different brain regions, identified using the AAL grey matter atlas and the ICBM DTI-81 white matter parcellation atlas (International Consortium for Brain Mapping, CA) respectively. This approach achieved a sensitivity of 94% and a specificity of 93% for both white and grey matter lesions [13]. Although it was observed that the specificity of WM lesion segmentations alone were comparatively lower due to more frequent variations in WM intensity, while the sensitivity of GM lesion segmentations alone were comparatively lower due to reduced contrast observed for GM lesions. Illustrations of white and grey matter lesion segmentations are shown in Fig 4.

[Figure omitted. See PDF.]

Fig 4. Illustration of lesion segmentation.

(A) GM lesion segmentation, (B), WM lesion segmentation and (C) internal capsule lesion segmentation.

https://doi.org/10.1371/journal.pone.0181605.g004

Identifying ventricular enlargement

Ventricular enlargement was identified using a SSM of healthy lateral ventricles to extract volumes of enlargement, and compute their impingement on nearby deep grey matter anatomies [14]. Biomarkers from this approach include a volume of ventricular enlargement (in mL) in the deep grey matter anatomies, as determined by the AAL grey matter atlas. Note that this is not the entire lateral ventricular volume, solely the volume of the region that is thought to be enlarged, compared to typically developing children. Three cases illustrating the segmented regions are shown in Fig 5.

[Figure omitted. See PDF.]

Fig 5. Illustration of three cases of localised enlarged ventricles due to injury, and the segmentation of the enlarged volumes from the nearest healthy ventricle shape.

https://doi.org/10.1371/journal.pone.0181605.g005

Statistical methodology

The Least Absolute Shrinkage and Selection Operator (LASSO) method [30] was used to define optimal sets of relevant biomarkers of brain injury associated with patient outcomes. Data was partitioned, with models trained on 75% of the data and their performance validated on the independent 25% test set. In the LASSO approach, predictor variables that are not strongly associated with the clinical outcomes are removed from the model. The LASSO was run with default lambda and alpha parameters that define the weight, and type of the regularisation penalty respectively. Patient age and gender were included in all regression models, as it has been previously shown that there are age [31] and gender [32] impacts on both motor and cognitive outcomes. To account for the variance in different scanner sequences used in this study, scanner sequence was also included as a covariate in each model.

The regression coefficients of each model, which encode the relative weights of individual biomarkers impacting clinical outcome, were enforced to be negative. This was done in order to only extract correlations where observed injury led to reductions in patient outcome. In total, six models were constructed, one for each of the six clinical scores (AHA, BRIEF, SDQ, TVPS, WR and VOC). Multiple comparisons of the models were corrected for using a Bonferroni correction. An analysis of variance analysis (ANOVA) was performed to compare the complete models from all biomarkers combined, and the biomarkers of each type of injury individually, to ascertain if it is beneficial to look at all kinds of injury observed in the MRIs concurrently. All statistical analysis were conducted using the R statistical software Version 3.2 [33].

Results

Characterisation of injury in the cohort

The prevalence of injury, as determined using the automated analyses, in the cohort of children with CP is illustrated in Fig 6. Prevalence is shown as a Venn diagram, to illustrate that each individual patients may have a combination of the three classes of injury. Binary characterisation of injury was classified as the presence of any enlarged ventricles or white/grey matter lesions (greater than 0 mL anywhere in the brain), or a regional cortical shape with a computed z-score > 2.5 compared to corresponding healthy cortical shape (for any cortical shape and region). Of the 139 children in our cohort, 109 were classified as having some form of injury, with 15% of the cohort having all three types of injury identified. Ventricular enlargement was the most common form of injury (observed in 68% of children), followed by white/grey matter lesions (55%) and cortical abnormalities (30%). As expected from the known aetiology of injury, ventricular enlargement was observed in over 60% of patients with white/grey matter lesions. False positive characterisations of ventricular enlargement, caused by slight alignment errors of the ventricular meshes, explains why 14 children with typical brain development were incorrectly identified as possessing some form of brain injury.

[Figure omitted. See PDF.]

Fig 6. A Venn diagram characterising the observed prevalence of ventricular enlargement, cortical malformations and white/grey matter lesions observed in our cohort using the described automated techniques.

https://doi.org/10.1371/journal.pone.0181605.g006

Observed structure-function relationships

The image biomarkers retained from LASSO, and their respective regression coefficients for each of the six models, are given in Table 2 below. In these models, the cortical shape biomarkers (labelled ‘Cortical thickness’, ‘Curvature’ or ‘Sulcal depth’) represent the reduction in clinical outcome for a 1 unit increase in z-score from healthy shape in that cortical region, while the lesion biomarkers (labelled ‘Lesion’) and ventricular enlargement biomarkers (labelled ‘Ventricle enlargement’) represent the decrease in clinical score from an increase in 1mL of lesion, or enlarged ventricle respectively, in that anatomy. In this table, we note that the p-values of each biomarker are uncorrected, and simply reflect the strength of that feature within the chosen model. However, the adjusted R-squared of the models, which measure strength of the correlation between the outcome and the model predictions, were compared against a Bonferroni corrected alpha value (0.05/6 tests, 0.008).

[Figure omitted. See PDF.]

Table 2. The retained anatomical regions, and corresponding regression coefficients (including standard errors) of the six regression models modelled on the 75% training set.

For each model, the multiple R-squared is provided. Model features that are significant (p < 0.05) are bolded.

https://doi.org/10.1371/journal.pone.0181605.t002

Analysis of predictive biomarkers

The performance of the models on the unseen test set, shown in Table 3, demonstrates how well these observed relationships could predict functional outcomes in children with observed brain injury. Five of the six data-driven models were statistically significant (p < 0.008) in the 25% test set. Although the correlations for the manually chosen biomarkers were generally smaller than the corresponding data-driven models, four of these six models were statistically significant (p < 0.008).

[Figure omitted. See PDF.]

Table 3. Pearson’s R correlation between the predicted outcomes in the test set using the trained regression models and the clinical scores of the test set, for both the data-driven and manually chosen models.

https://doi.org/10.1371/journal.pone.0181605.t003

Assessment of independence of biomarker-type

To assess the independence of the different categories of biomarkers, an ANOVA was performed to compare the complete models obtained from the data-driven variable selection, and models constructed using only the biomarkers from this model of a specific injury class (i.e. cortical biomarkers and lesion biomarkers separately). The summary of these analyses is provided in Supplementary Table 1. Model comparisons to the ventricle enlargement biomarkers alone were omitted as no ventricular enlargement biomarkers were retained by LASSO in any model. All models containing only cortical shape or lesion burden biomarkers were found to be significantly different from the complete regression models (p < 0.05), except for the cortical biomarker only model for TVPS (p = 0.109). This suggests that the two different sets of biomarkers explain different portions of variance in the clinical score of children with CP.

Discussion

The prevalence of brain injury observed in the present study using the automated approaches was found to be comparable to studies on other cohorts of children with CP [8,9,11]. Leveraging the quantifications that these automated methods allow, models of brain injury to patient function were generated. Five of the six models were observed to obtain significant correlations in this test set, demonstrating the clinical utility of these predictive models in estimating patient function in order to provide early and effective therapeutic interventions. Only the Brief model did not achieve significance, which may be due to reduced executive outcomes manifesting primarily as the loss of white matter connectivity, which would require diffusion MRI and tractography to identify [34]. On its own sMRI does not assess neural connectivity and the inclusion of diffusion data may result in greater association of these outcomes. However, unlike sMRI, there is not yet a routine and widely accepted method for diffusion MRI acquisition and processing of data [35]. These findings takes steps towards establishing standardised biomarkers used in the assessment of children with CP, with implications on estimating patient function, allowing effective interventions to be performed early in life, and to characterise structural brain changes as a result of therapy, improving the understanding of neuroplastic mechanisms.

Although the model features were retained from a data-driven process, the significant regions (p<0.05) conform to known structure-function relationships of the brain. For instance, the supplementary motor area has known associations with intended voluntary action [36], the primary somatosensory cortex has been associated with the mental rehearsal of motor acts [37], the anterior cingulate has known motor regions regulating the interactions between cognitive and motor control [38], supporting the observed importance of these regions in the AHA model (Table 2). Additional significant regions in the AHA model include the lenticular nucleus, which is known to control a variety of movements [39], and is connected to the external capsule [40], and the cerebral peduncle, which contains WM fibres of the corticospinal tract and has previously been shown to be predictive of motor deficits [41]. For the Brief and SDQ models of cognitive function, several of the retained regions (inferior temporal gyrus, primary motor cortex, insular cortex, cingulate cortex, premotor cortex and corona radiata) have known associations with cognition and executive functioning. For instance, the inferior temporal gyrus, which was significantly retained in both models, has a known role in verbal fluency and cognition, and is commonly affected in AD [42]. The observation in the Brief model that girls had higher reported executive functioning than boys may be due to females maturing earlier than males [43], or arising from the currently unexplained phenomenon that males with CP tend to exhibit more severe impairment [44]. Other significant features retained in the SDQ model include the primary motor and premotor cortices, which have both been shown to be responsible for motor learning and cognitive actions [45,46], the cingulate cortex that has been associated with the regulation of varied mental and emotional activity [47], and the corona radiata which has been associated with mental calculation and information processing [48]. The insular cortex, which was retained in both models for cognition (SDQ) and vision (TVPS), has known associations with multiple sensory areas, influencing awareness of somatosensation and goal directed cognition [49], and the processing of visual information [50]. We hypothesise that the presence of the occipital gyrus in the word reasoning model may arise from the known roles in the auditory and spatial information and stimuli [51], including sound localisation, which are required to perform well in the WPPSI-III subtests for communication. Age was retained in both models for vocabulary (VOC) and word reasoning (WR), potentially reflecting the known improvements in articulation [52] and phonological skills [53] occurring during childhood development.

Ventricular enlargement biomarkers were observed to be removed by LASSO’s variable selection for all models of clinical function. This may due to ventricular enlargement being a secondary injury caused by primary periventricular white matter tissue loss [54,55]. Consistent with this aetiology, most children with ventricular enlargement in this cohort had white or grey matter tissue lesions. As a result, both ventricular enlargement and lesion biomarkers explain similar portions of variance in these children’s functional outcomes. Since the ventricular enlargement biomarkers were removed by LASSO, this indicates that they explain relatively less variance in the functional outcomes than the lesion biomarkers. However, it was observed that both cortical morphology and lesion burden variables were retained by LASSO for all six models. The independence of these biomarkers is further demonstrated in the ANOVA supplementary table, identifying that the two sets of factors explain independent portions of variance in five of the six models. These findings highlight the benefit of quantifying both biomarkers for the assessment of MRIs.

Although the developed automated approaches were designed to be robust to the presence of severe brain injury, there are a number of technical limitations with these methods. Firstly, despite most of these approaches requiring only an affine alignment of atlases, which is mostly driven by the alignment of the skull, there are registration errors introduced in the deformable registration of the tissue probability maps used in the lesion segmentation pipeline. The brain masking segmentation step occasionally incorporates dura into the brain mask and miss subtle cortical sulci. In these cases, manual editing of the segmentations is required, as instances of either error will affect the cortical shape measures computed in the cortical analysis pipeline. The presence of slight correspondence errors between the ventricular shape models led to a number of false positive classifications of ventricular enlargement. Due to the risk of over-fitting on the limited available training data, no interaction terms were included for account for potential feature covariance, nor was WM connectivity information incorporated in the current study. The inclusion of diffusion data the model will be important in the future, because structural MRI alone cannot account for the variable influence of neuroplastic mechanisms that lead to altered structure-function relationships to compensate for the presence of injury [56,57]. These altered relationships confound the relationships between injury and impairment that the regression models attempt to elucidate, introducing unexplained variance in the clinical scores and potentially reducing the multiple R-squared of these models. Performing WM tract labelling, which has previously been achieved using label fusion [58], regional diffusion measures such as Fractional Anisotropy (FA) and Mean Diffusivity (MD) can be computed and incorporated into the regression models presented here. Diffusion measures such as FA and MD in key regions, such as the corpus callosum, have been associated with neurocognitive performance [59,60]. In the future, similar approaches can leverage longitudinal information when it becomes available, such as MRI taken early in life (around 2 years of age) and developmental trajectories through early childhood, in order to provide estimates of long-term function, which have been performed in the context of Alzheimer’s disease progression [61].

The main strength of this study is that the developed injury segmentation approaches only requiring well established T1-weighted sequences. More sophisticated diffusion and functional sequences are not as widely available, and have longer scanning times which complicate the imaging of young children. The development of automated techniques that only require these rapid and established sequences help facilitate potential translation to clinical practice. Secondly, this study encapsulates the full range of primary and secondary injuries characterised previously [7], as well as patient age and gender, allowing for comprehensive characterisation of brain injury from sMRI. In future, image analysis in the CP setting needs to move towards more automated quantification of injury, as it allows quantitative relationships between brain injury and outcome to be elucidated for children with both unilateral and bilateral CP. These associations have the potential to build upon the understanding of the brains structure-function relationships, which is an important prerequisite for understanding reorganisation and plasticity. Furthermore, there is potential in using these associations to produce estimates of patient function, which have important clinical implications for dictating the type and intensity of intervention that may be required, early in life in order to optimise functional outcomes.

Conclusion

In this study, several validated, automated approaches for identifying three different aetiologies of brain injury were applied to a cohort of children diagnosed with unilateral CP. A similar prevalence of cortical malformations, white and grey matter lesions and secondary ventricular enlargement were observed to previous studies on populations of children with CP. Significant correlations were observed between the biomarkers of injury and multiple scores of patient motor, cognitive, visual and communicative function. These structure-function relationships generalised to unseen data, with correlations between 0.482 and 0.795, and retained predictor variables that conformed to known roles of brain structures. There was significant independence between cortical morphology and tissue lesions in explaining functional outcomes, highlighting the benefit of quantifying both of these types of injury for clinical assessment. Furthermore, the strong and significant association between quantifications of injury obtained from on structural MRI and multiple clinical scores agree with empirically established structure-function relationships. These findings support the early elucidation and characterisation of brain injury from sMRI in the clinical assessment of children with CP. In future, the methods detailed in this paper could facilitate the longitudinal assessment of patients during atypical development and rehabilitation, and will assist in quantifying structural changes due to therapeutic effects. As a result, the selection of optimal treatment strategies for individual cases can be performed, to ultimately improve outcomes for children with CP.

Supporting information

[Figure omitted. See PDF.]

S1 Table. Analysis of independence between cortical morphology and lesion burden.

ANOVA model comparisons between the complete regression models, and the models constructed with the cortical shape and lesion burden biomarkers only.

https://doi.org/10.1371/journal.pone.0181605.s001

(PDF)

Acknowledgments

Alex M. Pagnozzi is supported by the Australian Postgraduate Award (APA) from The University of Queensland, and the Commonwealth Scientific Industrial and Research Organisation (CSIRO). Roslyn N. Boyd is supported by a Foundation for Children Grant, NHMRC Research Fellowship (2015–2020) and a NHMRC Project Grant COMBIT (1003887). Roslyn N. Boyd and Stephen Rose are supported by the by the Smart Futures Co-Investment Program Grant. Andrew Bradley is supported by an ARC Future fellowship (FT110100623). The funding bodies have not contributed to the study design, the collection, management, analysis and interpretation of data, the writing of final reports or the decision to submit findings for publication. No other authors have potential conflicts of interest to declare.

Citation: Pagnozzi AM, Dowson N, Doecke J, Fiori S, Bradley AP, Boyd RN, et al. (2017) Identifying relevant biomarkers of brain injury from structural MRI: Validation using automated approaches in children with unilateral cerebral palsy. PLoS ONE 12(8): e0181605. https://doi.org/10.1371/journal.pone.0181605

References

1. Rosenbaum P, Nigel P, Leviton A, Goldstein M, Bax M. A report: the definition and classification of cerebral palsy April 2006. Dev Med Child Neurol. 2007;49: 8–14.

2. Rosenbaum PL, Walter SD, Hanna SE, Palisano RJ, Russell DJ, Raina P, et al. Prognosis for Gross Motor Function in Cerebral Palsy. JAMA. American Medical Association; 2002;288: 1357. pmid:12234229

3. Himmelmann K, Beckung E, Hagberg G, Uvebrant P. Gross and fine motor function and accompanying impairments in cerebral palsy. Dev Med Child Neurol. 2006;48: 417–23. pmid:16700930

4. Odding E, Roebrock ME, Stam HJ, Roebroeck ME, Stam HJ. The epidemiology of cerebral palsy: Incidence, impairments and risk factors. Disabil Rehabil. 2006;28: 183–191. pmid:16467053

5. Korzeniewski SJ, Birbeck G, DeLano MC, Potchen MJ, Paneth N. A Systematic Review of Neuroimaging for Cerebral Palsy. J Child Neurol. 2008;23: 216–227. pmid:18263759

6. Accardo J, Kammann H, Hoon AH Jr. Neuroimaging in cerebral palsy. J Pediatr. 2004;145: S19–S27. pmid:15292883

7. Krägeloh-Mann I, Horber V, Krägeloh-Mann I, Horber V. The role of magnetic resonance imaging in elucidating the pathogenesis of cerebral palsy: a systematic review. Dev Med Child Neurol. 2007;49: 144–151. pmid:17254004

8. Legault G, Shevell MI, Dagenais L. Predicting Comorbidities With Neuroimaging in Children With Cerebral Palsy. Pediatr Neurol. 2011;45: 229–232. pmid:21907883

9. Bax M, Tydeman C, Flodmark O. Clinical and MRI correlates of cerebral palsy: the European Cerebral Palsy Study. JAMA. American Medical Association; 2006;296: 1602–8. pmid:17018805

10. Arnfield E, Guzzetta A, Boyd R. Relationship between brain structure on magnetic resonance imaging and motor outcomes in children with cerebral palsy: A systematic review. Res Dev Disabil. 2013;34: 2234–2250. pmid:23643774

11. Krägeloh-Mann I, Horber V. The role of magnetic resonance imaging in elucidating the pathogenesis of cerebral palsy: a systematic review. Dev Med Child Neurol. 2007;49: 144–151. pmid:17254004

12. Pagnozzi AM, Dowson N, Fiori S, Doecke J, Bradley AP, Boyd RN, et al. Alterations in regional shape on ipsilateral and contralateral cortex contrast in children with unilateral cerebral palsy and are predictive of multiple outcomes. Hum Brain Mapp. 2016;37: 3588–3603. pmid:27259165

13. Pagnozzi AM, Dowson N, Doecke J, Fiori S, Bradley AP, Boyd RN, et al. Automated, quantitative measures of grey and white matter lesion burden correlates with motor and cognitive function in children with unilateral cerebral palsy. NeuroImage Clin. Elsevier; 2016;11: 751–759. pmid:27330975

14. Pagnozzi AM, Shen K, Doecke JD, Boyd RN, Bradley AP, Rose S, et al. Using ventricular modeling to robustly probe significant deep gray matter pathologies: Application to cerebral palsy. Hum Brain Mapp. 2016;37: 3795–3809. pmid:27257958

15. Boyd RN, Ziviani J, Sakzewski L, Miller L, Bowden J, Cunnington R, et al. COMBIT: protocol of a randomised comparison trial of COMbined modified constraint induced movement therapy and bimanual intensive training with distributed model of standard upper limb rehabilitation in children with congenital hemiplegia. BMC Neurol. 2013;13: 68. pmid:23809257

16. Boyd RN, Mitchell LE, James ST, Ziviani J, Sakzewski L, Smith A, et al. Move it to improve it (Mitii): study protocol of a randomised controlled trial of a novel web-based multimodal training program for children and adolescents with cerebral palsy. BMJ Open. 2013;3: e002853. pmid:23578686

17. Pagnozzi AM, Dowson N, Bourgeat P, Bradley AP, Boyd RN, Rose S. Expectation-maximization with image-weighted Markov Random Fields to handle severe pathology. Digital Image Computing: Techniques and Applications (DICTA). Adelaide, SA; 2015. pp. 1–6.

18. Krumlinde-Sundholm L, Holmefur M, Kottorp A, Eliasson AC. The Assisting Hand Assessment: current evidence of validity, reliability, and responsiveness to change. Dev Med Child Neurol. 2007;49: 259–64. pmid:17376135

19. Gioia GA, Isquith PK, Retzlaff PD, Espy KA. Confirmatory factor analysis of the Behavior Rating Inventory of Executive Function (BRIEF) in a clinical sample. Child Neuropsychol. Routledge; 2002;8: 249–57. pmid:12759822

20. Bourdon KH, Goodman R, Rae DS, Simpson G, Koretz DS. The Strengths and Difficulties Questionnaire: U.S. normative data and psychometric properties. J Am Acad Child Adolesc Psychiatry. 2005;44: 557–64. pmid:15908838

21. Frostig M, Lefever DW, Whittlesey JRB. A developmental test of visual perception for evaluating normal and neurologically handicapped children. Percept Mot Skills. Ammons Scientific; 1961;12: 383–394.

22. Wechsler D. Wechsler preschool and primary scale of intelligence. Psychological Corporation New York; 1967.

23. Fiori S, Guzzetta A, Pannek K, Ware RS, Rossi G, Klingels K, et al. Validity of semi-quantitative scale for brain MRI in unilateral cerebral palsy due to periventricular white matter lesions: Relationship with hand sensorimotor function and structural connectivity. NeuroImage Clin. 2015;8: 104–109. pmid:26106533

24. The McConnell Brain Imaging Centre. Colin 27 Average Brain, Stereotaxic Registration Model, original 1998 version [Internet]. 2012 [cited 24 Jul 2014]. Available: http://www.bic.mni.mcgill.ca/ServicesAtlases/Colin27

25. Rivest-Hénault D, Dowson N, Greer PB, Fripp J, Dowling J. Robust inverse-consistent affine CT–MR registration in MRI-assisted and MRI-alone prostate radiation therapy. Med Image Anal. 2015;23: 56–69. pmid:25966468

26. Tustison NJ, Avants BB, Cook PA, Zheng Y, Egan A, Yushkevich PA, et al. N4ITK: improved N3 bias correction. IEEE Trans Med Imaging. 2010;29: 1310–20. pmid:20378467

27. Perona P, Malik J. Scale-space and edge detection using anisotropic diffusion. IEEE Trans Pattern Anal Mach Intell. 1990;12: 629–639.

28. Yezzi A. Modified curvature motion for image smoothing and enhancement. IEEE Trans Image Process. 1998;7: 345–352. pmid:18276254

29. Modat M, Ridgway GR, Taylor ZA, Lehmann M, Barnes J, Hawkes DJ, et al. Fast free-form deformation using graphics processing units. Comput Methods Programs Biomed. 2010;98: 278–84. pmid:19818524

30. Tibshirani R. Regression Shrinkage and Selection via the Lasso. J R Stat Soc Ser B. 1996;58: 267–288.

31. HANNA SE, ROSENBAUM PL, BARTLETT DJ, PALISANO RJ, WALTER SD, AVERY L, et al. Stability and decline in gross motor function among children and youth with cerebral palsy aged 2 to 21 years. Dev Med Child Neurol. Blackwell Publishing Ltd; 2009;51: 295–302. pmid:19391185

32. Romeo DMM, Cioni M, Battaglia LR, Palermo F, Mazzone D. Spectrum of gross motor and cognitive functions in children with cerebral palsy: Gender differences. Eur J Paediatr Neurol. 2011;15: 53–58. pmid:20542713

33. The R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing [Internet]. Vienna, Austria; 2008. Available: http://www.r-project.org.

34. Rose S, Guzzetta A, Pannek K, Boyd R. MRI Structural Connectivity, Disruption of Primary Sensorimotor Pathways, and Hand Function in Cerebral Palsy. Brain Connect. Mary Ann Liebert, Inc. 140 Huguenot Street, 3rd Floor New Rochelle, NY 10801 USA; 2011;1: 309–316. pmid:22432420

35. Jones DK, Knösche TR, Turner R. White matter integrity, fiber count, and other fallacies: The do’s and don’ts of diffusion MRI. Neuroimage. 2013;73: 239–254. pmid:22846632

36. Goldberg G. Supplementary motor area structure and function: Review and hypotheses. Behav Brain Sci. Cambridge University Press; 2010;8: 567.

37. Porro CA, Francescato MP, Cettolo V, Diamond ME, Baraldi P, Zuiani C, et al. Primary motor and sensory cortex activation during motor performance and motor imagery: a functional magnetic resonance imaging study. J Neurosci. 1996;16: 7688–98. Available: http://www.ncbi.nlm.nih.gov/pubmed/8922425 pmid:8922425

38. Paus T. Primate Anterior Cingulate Cortex: Where Motor Control, Drive and Cognition Interface. Nat Rev Neurosci. 2001;2: 417–24. pmid:11389475

39. Middleton FA, Strick PL. Basal ganglia and cerebellar loops: motor and cognitive circuits. Brain Res Rev. 2000;31: 236–250. pmid:10719151

40. Henry RG, Berman JI, Nagarajan SS, Mukherjee P, Berger MS. Subcortical pathways serving cortical language sites: initial experience with diffusion tensor imaging fiber tracking combined with intraoperative language mapping. Neuroimage. 2004;21: 616–22. pmid:14980564

41. Laundre BJ, Jellison BJ, Badie B, Alexander AL, Field AS. Diffusion Tensor Imaging of the Corticospinal Tract before and after Mass Resection as Correlated with Clinical Motor Findings: Preliminary Data. AJNR Am J Neuroradiol. 2005;26: 791–796. Available: http://www.ajnr.org/content/26/4/791.short pmid:15814922

42. Scheff SW, Price DA, Schmitt FA, Scheff MA, Mufson EJ. Synaptic loss in the inferior temporal gyrus in mild cognitive impairment and Alzheimer’s disease. J Alzheimers Dis. IOS Press; 2011;24: 547–57. pmid:21297265

43. Waber D. Sex differences in cognition: a function of maturation rate? Science (80-). American Association for the Advancement of Science; 1976;192: 572–574.

44. Jarvis S, Glinianaia S V, Arnaud C, Fauconnier J, Johnson A, McManus V, et al. Case gender and severity in cerebral palsy varies with intrauterine growth. Arch Dis Child. 2005;90: 474–9. pmid:15851428

45. Sanes JN, Donoghue JP. Plasticity and primary motor cortex. Annu Rev Neurosci. 2000;23: 393–415. pmid:10845069

46. Rizzolatti G, Fogassi L, Gallese V. Motor and cognitive functions of the ventral premotor cortex. Curr Opin Neurobiol. 2002;12: 149–154. pmid:12015230

47. Bush G, Luu P, Posner MI. Cognitive and emotional influences in anterior cingulate cortex. Trends Cogn Sci. 2000;4: 215–222. pmid:10827444

48. Sasson E, Doniger GM, Pasternak O, Tarrasch R, Assaf Y. Structural correlates of cognitive domains in normal aging with diffusion tensor imaging. Brain Struct Funct. 2012;217: 503–15. pmid:21909706

49. Chang LJ, Yarkoni T, Khaw MW, Sanfey AG. Decoding the role of the insula in human cognition: functional parcellation and large-scale reverse inference. Cereb Cortex. Oxford University Press; 2013;23: 739–49. pmid:22437053

50. Nagai M, Kishi K, Kato S. Insular cortex and neuropsychiatric disorders: a review of recent literature. Eur Psychiatry. 2007;22: 387–94. pmid:17416488

51. Renier LA, Anurova I, De Volder AG, Carlson S, VanMeter J, Rauschecker JP. Preserved functional specialization for spatial processing in the middle occipital gyrus of the early blind. Neuron. 2010;68: 138–48. pmid:20920797

52. Kilminster MGE, Laird EM. Articulation Development in Children Aged Three to Nine Years. Aust J Hum Commun Disord. Taylor & Francis; 2014;6: 23–30.

53. Snowling M, Hulme C. The development of phonological skills. Philos Trans R Soc Lond B Biol Sci. 1994;346: 21–7. pmid:7886149

54. Melhem ER, Hoon AH, Ferrucci JT, Quinn CB, Reinhardt EM, Demetrides SW, et al. Periventricular Leukomalacia: Relationship between Lateral Ventricular Volume on Brain MR Images and Severity of Cognitive and Motor Impairment1. Radiology. 2000;214: 199–204. Available: http://radiology.rsna.org/content/214/1/199 pmid:10644124

55. Palmer FB. Strategies for the early diagnosis of cerebral palsy. J Pediatr. 2004;145: S8–S11. pmid:15292881

56. Thiel A, Herholz K, Koyuncu A, Ghaemi M, Kracht LW, Habedank B, et al. Plasticity of language networks in patients with brain tumors: a positron emission tomography activation study. Ann Neurol. 2001;50: 620–9. Available: http://www.ncbi.nlm.nih.gov/pubmed/11706968 pmid:11706968

57. Carmichael ST. Plasticity of Cortical Projections after Stroke. Neurosci. 2003;9: 64–75. pmid:12580341

58. Jin Y, Shi Y, Zhan L, Gutman BA, de Zubicaray GI, McMahon KL, et al. Automatic clustering of white matter fibers in brain diffusion MRI with an application to genetics. Neuroimage. Elsevier Inc.; 2014;100: 75–90. pmid:24821529

59. Dennis EL, Jin Y, Villalon-Reina JE, Zhan L, Kernan CL, Babikian T, et al. White matter disruption in moderate/severe pediatric traumatic brain injury: advanced tract-based analyses. NeuroImage Clin. 2015;7: 493–505. pmid:25737958

60. Dennis EL, Ellis MU, Marion SD, Jin Y, Moran L, Olsen A, et al. Callosal Function in Pediatric Traumatic Brain Injury Linked to Disrupted White Matter Integrity. J Neurosci. Society for Neuroscience; 2015;35: 10202–11. pmid:26180196

61. Huang L, Jin Y, Gao Y, Thung K-H, Shen D. Longitudinal clinical score prediction in Alzheimer’s disease with soft-split sparse regression based random forest. Neurobiol Aging. 2016;46: 180–191. pmid:27500865