Widespread, depth-dependent cortical microstructure alterations in pediatric focal epilepsy

OBJECTIVE

Tissue abnormalities in focal epilepsy may extend beyond the presumed focus. The underlying pathophysiology of these broader changes is unclear, and it is not known whether they result from ongoing disease processes or treatment-related side effects, or whether they emerge earlier. Few studies have focused on the period of onset for most focal epilepsies, childhood. Fewer still have utilized quantitative magnetic resonance imaging (MRI), which may provide a more sensitive and interpretable measure of tissue microstructural change. Here, we aimed to determine common spatial modes of changes in cortical architecture in children with heterogeneous drug-resistant focal epilepsy and, secondarily, whether changes were related to disease severity.

METHODS

To assess cortical microstructure, quantitative T1 and T2 relaxometry (qT1 and qT2) was measured in 43 children with drug-resistant focal epilepsy (age range = 4-18 years) and 46 typically developing children (age range = 2-18 years). We assessed depth-dependent qT1 and qT2 values across the neocortex, as well as their gradient of change across cortical depths. We also determined whether global changes seen in group analyses were driven by focal pathologies in individual patients. Finally, as a proof-of-concept, we trained a classifier using qT1 and qT2 gradient maps from patients with radiologically defined abnormalities (MRI positive) and healthy controls, and tested whether this could classify patients without reported radiological abnormalities (MRI negative).

RESULTS

We uncovered depth-dependent qT1 and qT2 increases in widespread cortical areas in patients, likely representing microstructural alterations in myelin or gliosis. Changes did not correlate with disease severity measures, suggesting they may represent antecedent neurobiological alterations. Using a classifier trained with MRI-positive patients and controls, sensitivity was 71.4% at 89.4% specificity on held-out MRI-negative patients.

SIGNIFICANCE

These findings suggest the presence of a potential imaging endophenotype of focal epilepsy, detectable irrespective of radiologically identified abnormalities.

Multi-Channel 4D Parametrized Atlas of Macro- and Microstructural Neonatal Brain Development

Structural (also known as anatomical) and diffusion MRI provide complimentary anatomical and microstructural characterization of early brain maturation. However, the existing models of the developing brain in time include only either structural or diffusion MRI channels. Furthermore, there is a lack of tools for combined analysis of structural and diffusion MRI in the same reference space. In this work, we propose a methodology to generate a multi-channel (MC) continuous spatio-temporal parametrized atlas of the brain development that combines multiple MRI-derived parameters in the same anatomical space during 37-44 weeks of postmenstrual age range. We co-align structural and diffusion MRI of 170 normal term subjects from the developing Human Connectomme Project using MC registration driven by both T2-weighted and orientation distribution functions channels and fit the Gompertz model to the signals and spatial transformations in time. The resulting atlas consists of 14 spatio-temporal microstructural indices and two parcellation maps delineating white matter tracts and neonatal transient structures. In order to demonstrate applicability of the atlas for quantitative region-specific studies, a comparison analysis of 140 term and 40 preterm subjects scanned at the term-equivalent age is performed using different MRI-derived microstructural indices in the atlas reference space for multiple white matter regions, including the transient compartments. The atlas and software will be available after publication of the article.

The developing Human Connectome Project (dHCP) automated resting-state functional processing framework for newborn infants

•

An automated and robust pipeline to minimally pre-process highly confounded neonatal fMRI data.

•

Includes integrated dynamic distortion and slice-to-volume motion correction.

•

A robust multimodal registration approach which includes custom neonatal templates.

•

Incorporates an automated and self-reporting QC framework to quantify data quality and identify issues for further inspection.

•

Data analysis of 538 infants imaged at 26–45 weeks post-menstrual age.

The developing Human Connectome Project (dHCP) aims to create a detailed 4-dimensional connectome of early life spanning 20–45 weeks post-menstrual age. This is being achieved through the acquisition of multi-modal MRI data from over 1000 in- and ex-utero subjects combined with the development of optimised pre-processing pipelines. In this paper we present an automated and robust pipeline to minimally pre-process highly confounded neonatal resting-state fMRI data, robustly, with low failure rates and high quality-assurance. The pipeline has been designed to specifically address the challenges that neonatal data presents including low and variable contrast and high levels of head motion. We provide a detailed description and evaluation of the pipeline which includes integrated slice-to-volume motion correction and dynamic susceptibility distortion correction, a robust multimodal registration approach, bespoke ICA-based denoising, and an automated QC framework. We assess these components on a large cohort of dHCP subjects and demonstrate that processing refinements integrated into the pipeline provide substantial reduction in movement related distortions, resulting in significant improvements in SNR, and detection of high quality RSNs from neonates.

Modelling brain development to detect white matter injury in term and preterm born neonates

Premature birth occurs during a period of rapid brain growth. In this context, interpreting clinical neuroimaging can be complicated by the typical changes in brain contrast, size and gyrification occurring in the background to any pathology. To model and describe this evolving background in brain shape and contrast, we used a Bayesian regression technique, Gaussian process regression, adapted to multiple correlated outputs. Using MRI, we simultaneously estimated brain tissue intensity on T1- and T2-weighted scans as well as local tissue shape in a large cohort of 408 neonates scanned cross-sectionally across the perinatal period. The resulting model provided a continuous estimate of brain shape and intensity, appropriate to age at scan, degree of prematurity and sex. Next, we investigated the clinical utility of this model to detect focal white matter injury. In individual neonates, we calculated deviations of a neonate's observed MRI from that predicted by the model to detect punctate white matter lesions with very good accuracy (area under the curve > 0.95). To investigate longitudinal consistency of the model, we calculated model deviations in 46 neonates who were scanned on a second occasion. These infants' voxelwise deviations from the model could be used to identify them from the other 408 images in 83% (T2-weighted) and 76% (T1-weighted) of cases, indicating an anatomical fingerprint. Our approach provides accurate estimates of non-linear changes in brain tissue intensity and shape with clear potential for radiological use.

Multimodal surface matching with higher-order smoothness constraints

In brain imaging, accurate alignment of cortical surfaces is fundamental to the statistical sensitivity and spatial localisation of group studies, and cortical surface-based alignment has generally been accepted to be superior to volume-based approaches at aligning cortical areas. However, human subjects have considerable variation in cortical folding, and in the location of functional areas relative to these folds. This makes alignment of cortical areas a challenging problem. The Multimodal Surface Matching (MSM) tool is a flexible, spherical registration approach that enables accurate registration of surfaces based on a variety of different features. Using MSM, we have previously shown that driving cross-subject surface alignment, using areal features, such as resting state-networks and myelin maps, improves group task fMRI statistics and map sharpness. However, the initial implementation of MSM's regularisation function did not penalize all forms of surface distortion evenly. In some cases, this allowed peak distortions to exceed neurobiologically plausible limits, unless regularisation strength was increased to a level which prevented the algorithm from fully maximizing surface alignment. Here we propose and implement a new regularisation penalty, derived from physically relevant equations of strain (deformation) energy, and demonstrate that its use leads to improved and more robust alignment of multimodal imaging data. In addition, since spherical warps incorporate projection distortions that are unavoidable when mapping from a convoluted cortical surface to the sphere, we also propose constraints that enforce smooth deformation of cortical anatomies. We test the impact of this approach for longitudinal modelling of cortical development for neonates (born between 31 and 43 weeks of post-menstrual age) and demonstrate that the proposed method increases the biological interpretability of the distortion fields and improves the statistical significance of population-based analysis relative to other spherical methods.