Computerized Analysis of Magnetic Resonance Images to Study Cerebral Anatomy in Developing Neonates

The study of cerebral anatomy in developing neonates is of great importance for the understanding of brain development during the early period of life. This dissertation therefore focuses on three challenges in the modelling of cerebral anatomy in neonates during brain development. The methods that have been developed all use Magnetic Resonance Images (MRI) as source data. To facilitate study of vascular development in the neonatal period, a set of image analysis algorithms are developed to automatically extract and model cerebral vessel trees. The whole process consists of cerebral vessel tracking from automatically placed seed points, vessel tree generation, and vasculature registration and matching. These algorithms have been tested on clinical Time-of- Flight (TOF) MR angiographic datasets. To facilitate study of the neonatal cortex a complete cerebral cortex segmentation and reconstruction pipeline has been developed. Segmentation of the neonatal cortex is not effectively done by existing algorithms designed for the adult brain because the contrast between grey and white matter is reversed. This causes pixels containing tissue mixtures to be incorrectly labelled by conventional methods. The neonatal cortical segmentation method that has been developed is based on a novel expectation-maximization (EM) method with explicit correction for mislabelled partial volume voxels. Based on the resulting cortical segmentation, an implicit surface evolution technique is adopted for the reconstruction of the cortex in neonates. The performance of the method is investigated by performing a detailed landmark study. To facilitate study of cortical development, a cortical surface registration algorithm for aligning the cortical surface is developed. The method first inflates extracted cortical surfaces and then performs a non-rigid surface registration using free-form deformations (FFDs) to remove residual alignment. Validation experiments using data labelled by an expert observer demonstrate that the method can capture local changes and follow the growth of specific sulcus

Multimodality evaluation of the pediatric brain: DTI and its competitors

The development of the human brain, from the fetal period until childhood, happens in a series of intertwined neurogenetical and histogenetical events that are influenced by environment. Neuronal proliferation and migration, cell aggregation, axonal ingrowth and outgrowth, dendritic arborisation, synaptic pruning and myelinisation contribute to the ‘plasticity of the developing brain'. These events taken together contribute to the establishment of adult-like neuroarchitecture required for normal brain function. With the advances in technology today, mostly due to the development of non-invasive neuroimaging tools, it is possible to analyze these structural events not only in anatomical space but also longitudinally in time. In this review we have highlighted current ‘state of the art' neuroimaging tools. Development of the new MRI acquisition sequences (DTI, CHARMED and phase imaging) provides valuable insight into the changes of the microstructural environment of the cortex and white matter. Development of MRI imaging tools dedicated for analysis of the acquired images (i) TBSS and ROI fiber tractography, (ii) new tissue segmentation techniques and (iii) morphometric analysis of the cortical mantle (cortical thickness and convolutions) allows the researchers to map the longitudinal changes in the macrostructure of the developing brain that go hand-in-hand with the acquisition of cognitive skills during childhood. Finally, the latest and the newest technologies, like connectom analysis and resting state fMRI connectivity analysis, today, for the first time provide the opportunity to study the developing brain through the prism of maturation of the systems and networks beyond individual anatomical areas. Combining these methods in the future and modeling the hierarchical organization of the brain might ultimately help to understand the mechanisms underlying complex brain structure function relationships of normal development and of developmental disorder

Quantification of cortical folding using MR image data

The cerebral cortex is a thin layer of tissue lining the brain where neural circuits perform important high level functions including sensory perception, motor control and language processing. In the third trimester the fetal cortex folds rapidly from a smooth sheet into a highly convoluted arrangement of gyri and sulci. Premature birth is a high risk factor for poor neurodevelopmental outcome and has been associated with abnormal cortical development, however the nature of the disruption to developmental processes is not fully understood. Recent developments in magnetic resonance imaging have allowed the acquisition of high quality brain images of preterms and also fetuses in-utero. The aim of this thesis is to develop techniques which quantify folding from these images in order to better understand cortical development in these two populations. A framework is presented that quantifies global and regional folding using curvature-based measures. This methodology was applied to fetuses over a wide gestational age range (21.7 to 38.9 weeks) for a large number of subjects (N = 80) extending our understanding of how the cortex folds through this critical developmental period. The changing relationship between the folding measures and gestational age was modelled with a Gompertz function which allowed an accurate prediction of physiological age. A spectral-based method is outlined for constructing a spatio-temporal surface atlas (a sequence of mean cortical surface meshes for weekly intervals). A key advantage of this method is the ability to do group-wise atlasing without bias to the anatomy of an initial reference subject. Mean surface templates were constructed for both fetuses and preterms allowing a preliminary comparison of mean cortical shape over the postmenstrual age range 28-36 weeks. Displacement patterns were revealed which intensified with increasing prematurity, however more work is needed to evaluate the reliability of these findings.Open Acces

Mapping the Early Cortical Folding Process in the Preterm Newborn Brain

In the developing human brain, the cortical sulci formation is a complex process starting from 14 weeks of gestation onward. The potential influence of underlying mechanisms (genetic, epigenetic, mechanical or environmental) is still poorly understood, because reliable quantification in vivo of the early folding is lacking. In this study, we investigate the sulcal emergence noninvasively in 35 preterm newborns, by applying dedicated postprocessing tools to magnetic resonance images acquired shortly after birth over a developmental period critical for the human cortex maturation (26-36 weeks of age). Through the original three-dimensional reconstruction of the interface between developing cortex and white matter and correlation with volumetric measurements, we document early sulcation in vivo, and quantify changes with age, gender, and the presence of small white matter lesions. We observe a trend towards lower cortical surface, smaller cortex, and white matter volumes, but equivalent sulcation in females compared with males. By precisely mapping the sulci, we highlight interindividual variability in time appearance and interhemispherical asymmetries, with a larger right superior temporal sulcus than the left. Thus, such an approach, included in a longitudinal follow-up, may provide early indicators on the structural basis of cortical functional specialization and abnormalities induced by genetic and environmental factor

Automatic segmentation of MR brain images with a convolutional neural network

Automatic segmentation in MR brain images is important for quantitative analysis in large-scale studies with images acquired at all ages. This paper presents a method for the automatic segmentation of MR brain images into a number of tissue classes using a convolutional neural network. To ensure that the method obtains accurate segmentation details as well as spatial consistency, the network uses multiple patch sizes and multiple convolution kernel sizes to acquire multi-scale information about each voxel. The method is not dependent on explicit features, but learns to recognise the information that is important for the classification based on training data. The method requires a single anatomical MR image only. The segmentation method is applied to five different data sets: coronal T2-weighted images of preterm infants acquired at 30 weeks postmenstrual age (PMA) and 40 weeks PMA, axial T2- weighted images of preterm infants acquired at 40 weeks PMA, axial T1-weighted images of ageing adults acquired at an average age of 70 years, and T1-weighted images of young adults acquired at an average age of 23 years. The method obtained the following average Dice coefficients over all segmented tissue classes for each data set, respectively: 0.87, 0.82, 0.84, 0.86 and 0.91. The results demonstrate that the method obtains accurate segmentations in all five sets, and hence demonstrates its robustness to differences in age and acquisition protocol

Advanced neuroimaging techniques to study the development of the cerebral cortex, subplate and thalamus in preterm infants at 3 Tesla

Preterm infants are at increased risk of neurodevelopmental delay, cognitive dysfunction, and behavioural disturbances. Recent studies of older preterm children with cognitive impairments implicate morphological and functional cortical abnormalities. However elucidation of the preterm cortical abnormalities has been challenging due to specific neonatal features. Using 3 Tesla neonatal MR images and Expectation Maximisation/Markov Random Field segmentation with incorporation of a novel knowledge based technique for removal of mislabelled partial volume voxels, neonatal 3D cortical extraction was possible from 25 to 48 weeks gestation. This enabled the study of the true cortical scaling exponent, cortical thickness, regional volumes and curvature measurements. It showed a relative excess of the cortical surface area for its volume which corresponded with a change in the intrinsic curvature and fissuration up to 36 weeks gestation, after which, the relative growth of the surface area and volume were proportional leading to dominant changes in the extrinsic curvature and cortical folding. Thus the curvature measurements showed an important mechanistic property of convolution. By term equivalent age, the cortex was thicker and there were changes in cortical curvature although there were no differences in the cortical surface area of preterm infants compared to term born controls. There were specific frontal and parietal deficits in the cortical volume. Diffusion MR showed that although the early cortical anisotropy diminished to noise levels by 35 weeks, the mean diffusivity reduced during the entire third trimester due to changes in the radial diffusivity. Regional variations in the mean diffusivity occurred during development with frontal abnormalities persisting at term equivalent age. Subplate and thalamic quantification showed important development features during the third trimester, however in the absence of overt lesions no associations with cortical measures were found. Thus this thesis provides interesting and novel insights into the macroscopic and microscopic development of the cortex.Imperial Users onl

The Developing Human Connectome Project: a minimal processing pipeline for neonatal cortical surface reconstruction

The Developing Human Connectome Project (dHCP) seeks to create the first 4-dimensional connectome of early life. Understanding this connectome in detail may provide insights into normal as well as abnormal patterns of brain development. Following established best practices adopted by the WU-MINN Human Connectome Project (HCP), and pioneered by FreeSurfer, the project utilises cortical surface-based processing pipelines. In this paper, we propose a fully automated processing pipeline for the structural Magnetic Resonance Imaging (MRI) of the developing neonatal brain. This proposed pipeline consists of a refined framework for cortical and sub-cortical volume segmentation, cortical surface extraction, and cortical surface inflation, which has been specifically designed to address considerable differences between adult and neonatal brains, as imaged using MRI. Using the proposed pipeline our results demonstrate that images collected from 465 subjects ranging from 28 to 45 weeks post-menstrual age (PMA) can be processed fully automatically; generating cortical surface models that are topologically correct, and correspond well with manual evaluations of tissue boundaries in 85% of cases. Results improve on state-of-the-art neonatal tissue segmentation models and significant errors were found in only 2% of cases, where these corresponded to subjects with high motion. Downstream, these surfaces will enhance comparisons of functional and diffusion MRI datasets, supporting the modelling of emerging patterns of brain connectivity

Characterisation of the Haemodynamic Response Function (HRF) in the neonatal brain using functional MRI

Background: Preterm birth is associated with a marked increase in the risk of later neurodevelopmental impairment. With the incidence rising, novel tools are needed to provide an improved understanding of the underlying pathology and better prognostic information. Functional Magnetic Resonance Imaging (fMRI) with Blood Oxygen Level Dependent (BOLD) contrast has the potential to add greatly to the knowledge gained through traditional MRI techniques. However, it has been rarely used with neonatal subjects due to difficulties in application and inconsistent results. Central to this is uncertainity regarding the effects of early brain development on the Haemodynamic Response Function (HRF), knowledge of which is fundamental to fMRI methodology and analysis. Hypotheses: (1) Well localised and positive BOLD functional responses can be identified in the neonatal brain. (2) The morphology of the neonatal HRF differs significantly during early human development. (3) The application of an age-appropriate HRF will improve the identification of functional responses in neonatal fMRI studies. Methods: To test these hypotheses, a systematic fMRI study of neonatal subjects was carried out using a custom made somatosensory stimulus, and an adapted study design and analysis pipeline. The neonatal HRF was then characterised using an event related study design. The potential future application of the findings was then tested in a series of small experiments. Results: Well localised and positive BOLD functional responses were identified in neonatal subjects, with a maturational tendency towards an increasingly complex pattern of activation. A positive amplitude HRF was identified in neonatal subjects, with a maturational trend of a decreasing time-to-peak and increasing positive peak amplitude. Application of the empirical HRF significantly improved the precision of analysis in further fMRI studies. Conclusions: fMRI can be used to study functional activity in the neonatal brain, and may provide vital new information about both development and pathology

Examining the development of brain structure in utero with fetal MRI, acquired as part of the Developing Human Connectome Project

The human brain is an incredibly complex organ, and the study of it traverses many scales across space and time. The development of the brain is a protracted process that begins embryonically but continues into adulthood. Although neural circuits have the capacity to adapt and are modulated throughout life, the major structural foundations are laid in utero during the fetal period, through a series of rapid but precisely timed, dynamic processes. These include neuronal proliferation, migration, differentiation, axonal pathfinding, and myelination, to name a few. The fetal origins of disease hypothesis proposed that a variety of non-communicable diseases emerging in childhood and adulthood could be traced back to a series of risk factors effecting neurodevelopment in utero (Barker 1995). Since this publication, many studies have shown that the structural scaffolding of the brain is vulnerable to external environmental influences and the perinatal developmental window is a crucial determinant of neurological health later in life. However, there remain many fundamental gaps in our understanding of it. The study of human brain development is riddled with biophysical, ethical, and technical challenges. The Developing Human Connectome Project (dHCP) was designed to tackle these specific challenges and produce high quality open-access perinatal MRI data, to enable researchers to investigate normal and abnormal neurodevelopment (Edwards et al., 2022). This thesis will focus on investigating the diffusion-weighted and anatomical (T2) imaging data acquired in the fetal period, between the second to third trimester (22 – 37 gestational weeks). The limitations of fetal MR data are ill-defined due to a lack of literature and therefore this thesis aims to explore the data through a series of critical and strategic analysis approaches that are mindful of the biophysical challenges associated with fetal imaging. A variety of analysis approaches are optimised to quantify structural brain development in utero, exploring avenues to relate the changes in MR signal to possible neurobiological correlates. In doing so, the work in this thesis aims to improve mechanistic understanding about how the human brain develops in utero, providing the clinical and medical imaging community with a normative reference point. To this aim, this thesis investigates fetal neurodevelopment with advanced in utero MRI methods, with a particular emphasis on diffusion MRI. Initially, the first chapter outlines a descriptive, average trajectory of diffusion metrics in different white matter fiber bundles across the second to third trimester. This work identified unique polynomial trajectories in diffusion metrics that characterise white matter development (Wilson et al., 2021). Guided by previous literature on the sensitivity of DWI to cellular processes, I formulated a hypothesis about the biophysical correlates of diffusion signal components that might underpin this trend in transitioning microstructure. This hypothesis accounted for the high sensitivity of the diffusion signal to a multitude of simultaneously occurring processes, such as the dissipating radial glial scaffold, commencement of pre-myelination and arborization of dendritic trees. In the next chapter, the methods were adapted to address this hypothesis by introducing another dimension, and charting changes in diffusion properties along developing fiber pathways. With this approach it was possible to identify compartment-specific microstructural maturation, refining the spatial and temporal specificity (Wilson et al., 2023). The results reveal that the dynamic fluctuations in the components of the diffusion signal correlate with observations from previous histological work. Overall, this work allowed me to consolidate my interpretation of the changing diffusion signal from the first chapter. It also serves to improve understanding about how diffusion signal properties are affected by processes in transient compartments of the fetal brain. The third chapter of this thesis addresses the hypothesis that cortical gyrification is influenced by both underlying fiber connectivity and cytoarchitecture. Using the same fetal imaging dataset, I analyse the tissue microstructural change underlying the formation of cortical folds. I investigate correlations between macrostructural surface features (curvature, sulcal depth) and tissue microstructural measures (diffusion tensor metrics, and multi-shell multi-tissue decomposition) in the subplate and cortical plate across gestational age, exploring this relationship both at the population level and within subjects. This study provides empirical evidence to support the hypotheses that microstructural properties in the subplate and cortical plate are altered with the development of sulci. The final chapter explores the data without anatomical priors, using a data-driven method to extract components that represent coordinated structural maturation. This analysis aims to examine if brain regions with coherent patterns of growth over the fetal period converge on neonatal functional networks. I extract spatially independent features from the anatomical imaging data and quantify the spatial overlap with pre-defined neonatal resting state networks. I hypothesised that coherent spatial patterns of anatomical development over the fetal period would map onto the functional networks observed in the neonatal period. Overall, this thesis provides new insight about the developmental contrast over the second to third trimester of human development, and the biophysical correlates affecting T2 and diffusion MR signal. The results highlight the utility of fetal MRI to research critical mechanisms of structural brain maturation in utero, including white matter development and cortical gyrification, bridging scales from neurobiological processes to whole brain macrostructure. their gendered constructions relating to women