Resting-state functional connectivity in children cooled for neonatal encephalopathy.

Therapeutic hypothermia improves outcomes following neonatal hypoxic-ischaemic encephalopathy, reducing cases of death and severe disability such as cerebral palsy compared with normothermia management. However, when cooled children reach early school-age, they have cognitive and motor impairments which are associated with underlying alterations to brain structure and white matter connectivity. It is unknown whether these differences in structural connectivity are associated with differences in functional connectivity between cooled children and healthy controls. Resting-state functional MRI has been used to characterize static and dynamic functional connectivity in children, both with typical development and those with neurodevelopmental disorders. Previous studies of resting-state brain networks in children with hypoxic-ischaemic encephalopathy have focussed on the neonatal period. In this study, we used resting-state fMRI to investigate static and dynamic functional connectivity in children aged 6-8 years who were cooled for neonatal hypoxic-ischaemic without cerebral palsy [n = 22, median age (interquartile range) 7.08 (6.85-7.52) years] and healthy controls matched for age, sex and socioeconomic status [n = 20, median age (interquartile range) 6.75 (6.48-7.25) years]. Using group independent component analysis, we identified 31 intrinsic functional connectivity networks consistent with those previously reported in children and adults. We found no case-control differences in the spatial maps of these intrinsic connectivity networks. We constructed subject-specific static functional connectivity networks by measuring pairwise Pearson correlations between component time courses and found no case-control differences in functional connectivity after false discovery rate correction. To study the time-varying organization of resting-state networks, we used sliding window correlations and deep clustering to investigate dynamic functional connectivity characteristics. We found k = 4 repetitively occurring functional connectivity states, which exhibited no case-control differences in dwell time, fractional occupancy or state functional connectivity matrices. In this small cohort, the spatiotemporal characteristics of resting-state brain networks in cooled children without severe disability were too subtle to be differentiated from healthy controls at early school-age, despite underlying differences in brain structure and white matter connectivity, possibly reflecting a level of recovery of healthy resting-state brain function. To our knowledge, this is the first study to investigate resting-state functional connectivity in children with hypoxic-ischaemic encephalopathy beyond the neonatal period and the first to investigate dynamic functional connectivity in any children with hypoxic-ischaemic encephalopathy.

Cerebellar growth, volume and diffusivity in children cooled for neonatal encephalopathy without cerebral palsy.

Children cooled for HIE and who did not develop cerebral palsy (CP) still underperform at early school age in motor and cognitive domains and have altered supra-tentorial brain volumes and white matter connectivity. We obtained T1-weighted and diffusion-weighted MRI, motor (MABC-2) and cognitive (WISC-IV) scores from children aged 6-8 years who were cooled for HIE secondary to perinatal asphyxia without CP (cases), and controls matched for age, sex, and socioeconomic status. In 35 case children, we measured cerebellar growth from infancy (age 4-15 days after birth) to childhood. In childhood, cerebellar volumes were measured in 26 cases and 23 controls. Diffusion properties (mean diffusivity, MD and fractional anisotropy, FA) were calculated in 24 cases and 19 controls, in 9 cerebellar regions. Cases with FSIQ ≤ 85 had reduced growth of cerebellar width compared to those with FSIQ > 85 (p = 0.0005). Regional cerebellar volumes were smaller in cases compared to controls (p < 0.05); these differences were not significant when normalised to total brain volume. There were no case-control differences in MD or FA. Interposed nucleus volume was more strongly associated with IQ in cases than in controls (p = 0.0196). Other associations with developmental outcome did not differ between cases and controls.

Mammillary body abnormalities and cognitive outcomes in children cooled for neonatal encephalopathy.

AIM: To evaluate mammillary body abnormalities in school-age children without cerebral palsy treated with therapeutic hypothermia for neonatal hypoxic-ischaemic encephalopathy (cases) and matched controls, and associations with cognitive outcome, hippocampal volume, and diffusivity in the mammillothalamic tract (MTT) and fornix. METHOD: Mammillary body abnormalities were scored from T1-weighted magnetic resonance imaging (MRI) in 32 cases and 35 controls (median age [interquartile range] 7 years [6 years 7 months-7 years 7 months] and 7 years 4 months [6 years 7 months-7 years 7 months] respectively). Cognition was assessed using the Wechsler Intelligence Scale for Children, Fourth Edition. Hippocampal volume (normalized by total brain volume) was measured from T1-weighted MRI. Radial diffusivity and fractional anisotropy were measured in the MTT and fornix, from diffusion-weighted MRI using deterministic tractography. RESULTS: More cases than controls had mammillary body abnormalities (34% vs 0%; p < 0.001). Cases with abnormal mammillary bodies had lower processing speed (p = 0.016) and full-scale IQ (p = 0.028) than cases without abnormal mammillary bodies, and lower scores than controls in all cognitive domains (p < 0.05). Cases with abnormal mammillary bodies had smaller hippocampi (left p = 0.016; right p = 0.004) and increased radial diffusivity in the right MTT (p = 0.004) compared with cases without mammillary body abnormalities. INTERPRETATION: Cooled children with mammillary body abnormalities at school-age have reduced cognitive scores, smaller hippocampi, and altered MTT microstructure compared with those without mammillary body abnormalities, and matched controls. WHAT THIS PAPER ADDS: Cooled children are at higher risk of mammillary body abnormalities than controls. Abnormal mammillary bodies are associated with reduced cognitive scores and smaller hippocampi. Abnormal mammillary bodies are associated with altered mammillothalamic tract diffusivity.

Development of the corpus callosum and cognition after neonatal encephalopathy.

OBJECTIVE: Neonatal imaging studies report corpus callosum abnormalities after neonatal hypoxic-ischaemic encephalopathy (HIE), but corpus callosum development and relation to cognition in childhood are unknown. Using magnetic resonance imaging (MRI), we examined the relationship between corpus callosum size, microstructure and cognitive and motor outcomes at early school-age children cooled for HIE (cases) without cerebral palsy compared to healthy, matched controls. A secondary aim was to examine the impact of HIE-related neonatal brain injury on corpus callosum size, microstructure and growth. METHODS: Participants aged 6-8 years underwent MRI, the Movement Assessment Battery for Children Second Edition and Wechsler Intelligence Scale for Children Fourth Edition. Cross-sectional area, volume, fractional anisotropy and radial diffusivity of the corpus callosum and five subdivisions were measured. Multivariable regression was used to assess associations between total motor score, full-scale IQ (FSIQ) and imaging metrics. RESULTS: Adjusting for age, sex and intracranial volume, cases (N = 40) compared to controls (N = 39) demonstrated reduced whole corpus callosum area (ß = -26.9, 95% confidence interval [CI] = -53.17, -0.58), volume (ß = -138.5, 95% CI = -267.54, -9.56), fractional anisotropy and increased radial diffusivity (P < 0.05) within segments II-V. In cases, segment V area (ß = 0.18, 95% CI = 0.004, 0.35), volume (ß = 0.04, 95% CI = 0.001, 0.079), whole corpus callosum fractional anisotropy (ß = 13.8 95% CI = 0.6, 27.1) and radial diffusivity (ß = -11.3, 95% CI = -22.22, -0.42) were associated with FSIQ. Growth of the corpus callosum was restricted in cases with a FSIQ ≤85, and volume was reduced in cases with mild neonatal multifocal injury compared to white matter injury alone. INTERPRETATION: Following neonatal HIE, morphological and microstructural changes in the corpus callosum are associated with reduced cognitive function at early school age.

Factors associated with MRI success in children cooled for neonatal encephalopathy and controls.

OBJECTIVE: To investigate if an association exists between motion artefacts on brain MRI and comprehension, co-ordination, or hyperactivity scores in children aged 6-8 years, cooled for neonatal encephalopathy (cases) and controls. METHODS: Case children (n = 50) without cerebral palsy were matched with 43 controls for age, sex, and socioeconomic status. Children underwent T1-weighted (T1w), diffusion-weighted image (DWI) brain MRI and cognitive, behavioural, and motor skills assessment. Stepwise multivariable logistic regression assessed associations between unsuccessful MRI and comprehension (including Weschler Intelligence Scale for Children (WISC-IV) verbal comprehension, working memory, processing speed and full-scale IQ), co-ordination (including Movement Assessment Battery for Children (MABC-2) balance, manual dexterity, aiming and catching, and total scores) and hyperactivity (including Strengths and Difficulties Questionnaire (SDQ) hyperactivity and total difficulties scores). RESULTS: Cases had lower odds of completing both T1w and DWIs (OR: 0.31, 95% CI 0.11-0.89). After adjusting for case-status and sex, lower MABC-2 balance score predicted unsuccessful T1w MRI (OR: 0.81, 95% CI 0.67-0.97, p = 0.022). Processing speed was negatively correlated with relative motion on DWI (r = -0.25, p = 0.026) and SDQ total difficulties score was lower for children with successful MRIs (p = 0.049). CONCLUSIONS: Motion artefacts on brain MRI in early school-age children are related to the developmental profile. IMPACT: Children who had moderate/severe neonatal encephalopathy are less likely to have successful MRI scans than matched controls. Motion artefact on MRI is associated with lower MABC-2 balance scores in both children who received therapeutic hypothermia for neonatal encephalopathy and matched controls, after controlling for case-status and sex. Exclusion of children with motion artefacts on brain MRI can introduce sampling bias, which impacts the utility of neuroimaging to understand the brain-behaviour relationship in children with functional impairments.

Brain volumes and functional outcomes in children without cerebral palsy after therapeutic hypothermia for neonatal hypoxic-ischaemic encephalopathy.

AIM: To investigate whether brain volumes were reduced in children aged 6 to 8 years without cerebral palsy, who underwent therapeutic hypothermia for neonatal hypoxic-ischaemic encephalopathy (patients), and matched controls, and to examine the relation between subcortical volumes and functional outcome. METHOD: We measured regional brain volumes in 31 patients and 32 controls (median age 7 years and 7 years 2 months respectively) from T1-weighted magnetic resonance imaging (MRI). We assessed cognition using the Wechsler Intelligence Scales for Children, Fourth Edition and motor ability using the Movement Assessment Battery for Children, Second Edition (MABC-2). RESULTS: Patients had lower volume of whole-brain grey matter, white matter, pallidi, hippocampi, and thalami than controls (false discovery rate-corrected p < 0.05). Differences in subcortical grey-matter volumes were not independent of total brain volume (TBV). In patients, hippocampal and thalamic volumes correlated with full-scale IQ (hippocampi, r = 0.477, p = 0.010; thalami, r = 0.452, p = 0.016) and MABC-2 total score (hippocampi, r = 0.526, p = 0.004; thalami, r = 0.505, p = 0.006) independent of age, sex, and TBV. No significant correlations were found in controls. In patients, cortical injury on neonatal MRI was associated with reduced volumes of hippocampi (p = 0.001), thalami (p = 0.002), grey matter (p = 0.015), and white matter (p = 0.013). INTERPRETATION: Children who underwent therapeutic hypothermia have reduced whole-brain grey and white-matter volumes, with associations between hippocampal and thalamic volumes and functional outcomes.

Communication skills in children aged 6-8 years, without cerebral palsy cooled for neonatal hypoxic-ischemic encephalopathy.

We assessed communication skills of 48 children without cerebral palsy (CP) treated with therapeutic hypothermia (TH) for neonatal hypoxic-ischemic encephalopathy (HIE) (cases) compared to 42 controls at early school-age and examined their association with white matter diffusion properties in both groups and 18-month Bayley-III developmental assessments in cases. Parents completed a Children's Communication Checklist (CCC-2) yielding a General Communication Composite (GCC), structural and pragmatic language scores and autistic-type behavior score. GCC ≤ 54 and thresholds of structural and pragmatic language score differences defined language impairment. Using tract-based spatial statistics (TBSS), fractional anisotropy (FA) was compared between 31 cases and 35 controls. Compared to controls, cases had lower GCC (p = 0.02), structural (p = 0.03) and pragmatic language score (p = 0.04) and higher language impairments (p = 0.03). GCC correlated with FA in the mid-body of the corpus callosum, the cingulum and the superior longitudinal fasciculus (p < 0.05) in cases. Bayley-III Language Composite correlated with GCC (r = 0.34, p = 0.017), structural (r = 0.34, p = 0.02) and pragmatic (r = 0.32, p = 0.03) language scores and autistic-type behaviors (r = 0.36, p = 0.01).

Using deep clustering to improve fMRI dynamic functional connectivity analysis.

Dynamic functional connectivity (dFC) analysis of resting-state fMRI data is commonly performed by calculating sliding-window correlations (SWC), followed by k-means clustering in order to assign each window to a given state. Studies using synthetic data have shown that k-means performance is highly dependent on sliding window parameters and signal-to-noise ratio. Additionally, sources of heterogeneity between subjects may affect the accuracy of group-level clustering, thus affecting measurements of dFC state temporal properties such as dwell time and fractional occupancy. This may result in spurious conclusions regarding differences between groups (e.g. when comparing a clinical population to healthy controls). Therefore, is it important to quantify the ability of k-means to estimate dFC state temporal properties when applied to cohorts of multiple subjects, and to explore ways in which clustering performance can be maximised. Here, we explore the use of dimensionality reduction methods prior to clustering in order to map high-dimensional data to a lower dimensional space, providing salient features to the subsequent clustering step. We assess the use of deep autoencoders for dimensionality reduction prior to applying k-means clustering to the encoded data. We compare this deep clustering method to dimensionality reduction using principle component analysis (PCA), uniform manifold approximation and projection (UMAP), as well as applying k-means to the original feature space using either L1 or L2 distance. We provide extensive quantitative evaluation of clustering performance using synthetic datasets, representing data from multiple heterogeneous subjects. In synthetic data we find that deep clustering gives the best performance, while other approaches are often insufficient to capture temporal properties of dFC states. We then demonstrate the application of each method to real-world data from human subjects and show that the choice of dimensionality reduction method has a significant effect on group-level measurements of state temporal properties.

An Age-Specific Atlas for Delineation of White Matter Pathways in Children Aged 6-8 Years.

Introduction: Diffusion magnetic resonance imaging (MRI) allows noninvasive assessment of white matter connectivity in typical development and of changes due to brain injury or pathology. Probabilistic white matter atlases allow diffusion metrics to be measured in specific white matter pathways, and are a critical component in spatial normalization for group analysis. However, given the known developmental changes in white matter it may be suboptimal to use an adult template when assessing data acquired from children. Methods: By averaging subject-specific fiber bundles from 28 children aged from 6 to 8 years, we created an age-specific probabilistic white matter atlas for 12 major white matter tracts. Using both the newly developed and Johns Hopkins adult atlases, we compared the atlas with subject-specific fiber bundles in two independent validation cohorts, assessing accuracy in terms of volumetric overlap and measured diffusion metrics. Results: Our age-specific atlas gave better overall performance than the adult atlas, achieving higher volumetric overlap with subject-specific fiber tracking and higher correlation of fractional anisotropy (FA) measurements with those measured from subject-specific fiber bundles. Specifically, estimates of FA values for corticospinal tract, uncinate fasciculus, forceps minor, cingulate gyrus part of the cingulum, and anterior thalamic radiation were all significantly more accurate when estimated with an age-specific atlas. Discussion: The age-specific atlas allows delineation of white matter tracts in children aged 6-8 years, without the need for tractography, more accurately than when normalizing to an adult atlas. To our knowledge, this is the first publicly available probabilistic atlas of white matter tracts for this age group.

Motor function and white matter connectivity in children cooled for neonatal encephalopathy.

Therapeutic hypothermia reduces the incidence of severe motor disability, such as cerebral palsy, following neonatal hypoxic-ischaemic encephalopathy. However, cooled children without cerebral palsy at school-age demonstrate motor deficits and altered white matter connectivity. In this study, we used diffusion-weighted imaging to investigate the relationship between white matter connectivity and motor performance, measured using the Movement Assessment Battery for Children-2, in children aged 6-8 years treated with therapeutic hypothermia for neonatal hypoxic-ischaemic encephalopathy at birth, who did not develop cerebral palsy (cases), and matched typically developing controls. Correlations between total motor scores and diffusion properties in major white matter tracts were assessed in 33 cases and 36 controls. In cases, significant correlations (FDR-corrected P < 0.05) were found in the anterior thalamic radiation bilaterally (left: r = 0.513; right: r = 0.488), the cingulate gyrus part of the left cingulum (r = 0.588), the hippocampal part of the left cingulum (r = 0.541), and the inferior fronto-occipital fasciculus bilaterally (left: r = 0.445; right: r = 0.494). No significant correlations were found in controls. We then constructed structural connectivity networks, for 22 cases and 32 controls, in which nodes represent brain regions and edges were determined by probabilistic tractography and weighted by fractional anisotropy. Analysis of whole-brain network metrics revealed correlations (FDR-corrected P < 0.05), in cases, between total motor scores and average node strength (r = 0.571), local efficiency (r = 0.664), global efficiency (r = 0.677), clustering coefficient (r = 0.608), and characteristic path length (r = -0.652). No significant correlations were found in controls. We then investigated edge-level association with motor function using the network-based statistic. This revealed subnetworks which exhibited group differences in the association between motor outcome and edge weights, for total motor scores (P = 0.0109) as well as for balance (P = 0.0245) and manual dexterity (P = 0.0233) domain scores. All three of these subnetworks comprised numerous frontal lobe regions known to be associated with motor function, including the superior frontal gyrus and middle frontal gyrus. The subnetwork associated with total motor scores was highly left-lateralised. These findings demonstrate an association between impaired motor function and brain organisation in school-age children treated with therapeutic hypothermia for neonatal hypoxic-ischaemic encephalopathy.

Disrupted brain connectivity in children treated with therapeutic hypothermia for neonatal encephalopathy.

Therapeutic hypothermia following neonatal encephalopathy due to birth asphyxia reduces death and cerebral palsy. However, school-age children without cerebral palsy treated with therapeutic hypothermia for neonatal encephalopathy still have reduced performance on cognitive and motor tests, attention difficulties, slower reaction times and reduced visuo-spatial processing abilities compared to typically developing controls. We acquired diffusion-weighted imaging data from school-age children without cerebral palsy treated with therapeutic hypothermia for neonatal encephalopathy at birth, and a matched control group. Voxelwise analysis (33 cases, 36 controls) confirmed reduced fractional anisotropy in widespread areas of white matter in cases, particularly in the fornix, corpus callosum, anterior and posterior limbs of the internal capsule bilaterally and cingulum bilaterally. In structural brain networks constructed using probabilistic tractography (22 cases, 32 controls), graph-theoretic measures of strength, local and global efficiency, clustering coefficient and characteristic path length were found to correlate with IQ in cases but not controls. Network-based statistic analysis implicated brain regions involved in visuo-spatial processing and attention, aligning with previous behavioural findings. These included the precuneus, thalamus, left superior parietal gyrus and left inferior temporal gyrus. Our findings demonstrate that, despite the manifest successes of therapeutic hypothermia, brain development is impaired in these children.