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

Migraine aura: retracting particle-like waves in weakly susceptible cortex

Cortical spreading depression (SD) has been suggested to underlie migraine aura. Despite a precise match in speed, the spatio-temporal patterns of SD and aura symptoms on the cortical surface ordinarily differ in aspects of size and shape. We show that this mismatch is reconciled by utilizing that both pattern types bifurcate from an instability point of generic reaction-diffusion models. To classify these spatio-temporal pattern we suggest a susceptibility scale having the value [sigma]=1 at the instability point. We predict that human cortex is only weakly susceptible to SD ([sigma]<1), and support this prediction by directly matching visual aura symptoms with anatomical landmarks using fMRI retinotopic mapping. We discuss the increased dynamical repertoire of cortical tissue close to [sigma]=1, in particular, the resulting implications on migraine pharmacology that is hitherto tested in the regime ([sigma]>>1), and potentially silent aura occurring below a second bifurcation point at [sigma]=0 on the susceptible scale

Imaging Informatics and the Human Brain Project: the Role of Structure, Review

(no abstract

A geometric network model of intrinsic grey-matter connectivity of the human brain

Network science provides a general framework for analysing the large-scale brain networks that naturally arise from modern neuroimaging studies, and a key goal in theoretical neuro- science is to understand the extent to which these neural architectures influence the dynamical processes they sustain. To date, brain network modelling has largely been conducted at the macroscale level (i.e. white-matter tracts), despite growing evidence of the role that local grey matter architecture plays in a variety of brain disorders. Here, we present a new model of intrinsic grey matter connectivity of the human connectome. Importantly, the new model incorporates detailed information on cortical geometry to construct ‘shortcuts’ through the thickness of the cortex, thus enabling spatially distant brain regions, as measured along the cortical surface, to communicate. Our study indicates that structures based on human brain surface information differ significantly, both in terms of their topological network characteristics and activity propagation properties, when compared against a variety of alternative geometries and generative algorithms. In particular, this might help explain histological patterns of grey matter connectivity, highlighting that observed connection distances may have arisen to maximise information processing ability, and that such gains are consistent with (and enhanced by) the presence of short-cut connections

Surface-Based tools for Characterizing the Human Brain Cortical Morphology

Tesis por compendio de publicacionesThe cortex of the human brain is highly convoluted. These characteristic convolutions present advantages over lissencephalic brains. For instance, gyrification allows an expansion of cortical surface area without significantly increasing the cranial volume, thus facilitating the pass of the head through the birth channel. Studying the human brain’s cortical morphology and the processes leading to the cortical folds has been critical for an increased understanding of the pathological processes driving psychiatric disorders such as schizophrenia, bipolar disorders, autism, or major depression. Furthermore, charting the normal developmental changes in cortical morphology during adolescence or aging can be of great importance for detecting deviances that may be precursors for pathology. However, the exact mechanisms that push cortical folding remain largely unknown. The accurate characterization of the neurodevelopment processes is challenging. Multiple mechanisms co-occur at a molecular or cellular level and can only be studied through the analysis of ex-vivo samples, usually of animal models. Magnetic Resonance Imaging can partially fill the breach, allowing the portrayal of the macroscopic processes surfacing on in-vivo samples. Different metrics have been defined to measure cortical structure to describe the brain’s morphological changes and infer the associated microstructural events. Metrics such as cortical thickness, surface area, or cortical volume help establish a relation between the measured voxels on a magnetic resonance image and the underlying biological processes. However, the existing methods present limitations or room for improvement. Methods extracting the lines representing the gyral and sulcal morphology tend to over- or underestimate the total length. These lines can provide important information about how sulcal and gyral regions function differently due to their distinctive ontogenesis. Nevertheless, some methods label every small fold on the cortical surface as a sulcal fundus, thus losing the perspective of lines that travel through the deeper zones of a sulcal basin. On the other hand, some methods are too restrictive, labeling sulcal fundi only for a bunch of primary folds. To overcome this issue, we have proposed a Laplacian-collapse-based algorithm that can delineate the lines traversing the top regions of the gyri and the fundi of the sulci avoiding anastomotic sulci. For this, the cortex, represented as a 3D surface, is segmented into gyral and sulcal surfaces attending to the curvature and depth at every point of the mesh. Each resulting surface is spatially filtered, smoothing the boundaries. Then, a Laplacian-collapse-based algorithm is applied to obtain a thinned representation of the morphology of each structure. These thin curves are processed to detect where the extremities or endpoints lie. Finally, sulcal fundi and gyral crown lines are obtained by eroding the surfaces while preserving the structure topology and connectivity between the endpoints. The assessment of the presented algorithm showed that the labeled sulcal lines were close to the proposed ground truth length values while crossing through the deeper (and more curved) regions. The tool also obtained reproducibility scores better or similar to those of previous algorithms. A second limitation of the existing metrics concerns the measurement of sulcal width. This metric, understood as the physical distance between the points on opposite sulcal banks, can come in handy in detecting cortical flattening or complementing the information provided by cortical thickness, gyrification index, or such features. Nevertheless, existing methods only provided averaged measurements for different predefined sulcal regions, greatly restricting the possibilities of sulcal width and ignoring the intra-region variability. Regarding this, we developed a method that estimates the distance from each sulcal point in the cortex to its corresponding opposite, thus providing a per-vertex map of the physical sulcal distances. For this, the cortical surface is sampled at different depth levels, detecting the points where the sulcal banks change. The points corresponding to each sulcal wall are matched with the closest point on a different one. The distance between those points is the sulcal width. The algorithm was validated against a simulated sulcus that resembles a simple fold. Then the tool was used on a real dataset and compared against two widely-used sulcal width estimation methods, averaging the proposed algorithm’s values into the same region definition those reference tools use. The resulting values were similar for the proposed and the reference methods, thus demonstrating the algorithm’s accuracy. Finally, both algorithms were tested on a real aging population dataset to prove the methods’ potential in a use-case scenario. The main idea was to elucidate fine-grained morphological changes in the human cortex with aging by conducting three analyses: a comparison of the age-dependencies of cortical thickness in gyral and sulcal lines, an analysis of how the sulcal and gyral length changes with age, and a vertex-wise study of sulcal width and cortical thickness. These analyses showed a general flattening of the cortex with aging, with interesting findings such as a differential age-dependency of thickness thinning in the sulcal and gyral regions. By demonstrating that our method can detect this difference, our results can pave the way for future in vivo studies focusing on macro- and microscopic changes specific to gyri or sulci. Our method can generate new brain-based biomarkers specific to sulci and gyri, and these can be used on large samples to establish normative models to which patients can be compared. In parallel, the vertex-wise analyses show that sulcal width is very sensitive to changes during aging, independent of cortical thickness. This corroborates the concept of sulcal width as a metric that explains, in the least, the unique variance of morphology not fully captured by existing metrics. Our method allows for sulcal width vertex-wise analyses that were not possible previously, potentially changing our understanding of how changes in sulcal width shape cortical morphology. In conclusion, this thesis presents two new tools, open source and publicly available, for estimating cortical surface-based morphometrics. The methods have been validated and assessed against existing algorithms. They have also been tested on a real dataset, providing new, exciting insights into cortical morphology and showing their potential for defining innovative biomarkers.Programa de Doctorado en Ciencia y Tecnología Biomédica por la Universidad Carlos III de MadridPresidente: Juan Domingo Gispert López.- Secretario: Norberto Malpica González de Vega.- Vocal: Gemma Cristina Monté Rubi

Automated Sulcal Depth Measurement on Cortical Surface Reflecting Geometrical Properties of Sulci

Sulcal depth that is one of the quantitative measures of cerebral cortex has been widely used as an important marker for brain morphological studies. Several studies have employed Euclidean (EUD) or geodesic (GED) algorithms to measure sulcal depth, which have limitations that ignore sulcal geometry in highly convoluted regions and result in under or overestimated depth. In this study, we proposed an automated measurement for sulcal depth on cortical surface reflecting geometrical properties of sulci, which named the adaptive distance transform (ADT). We first defined the volume region of cerebrospinal fluid between the 3D convex hull and the cortical surface, and constructed local coordinates for that restricted region. Dijkstra’s algorithm was then used to compute the shortest paths from the convex hull to the vertices of the cortical surface based on the local coordinates, which may be the most proper approach for defining sulcal depth. We applied our algorithm to both a clinical dataset including patients with mild Alzheimer’s disease (AD) and 25 normal controls and a simulated dataset whose shape was similar to a single sulcus. The mean sulcal depth in the mild AD group was significantly lower than controls (p = 0.007, normal [mean±SD]: 7.29±0.23 mm, AD: 7.11±0.29) and the area under the receiver operating characteristic curve was relatively high, showing the value of 0.818. Results from clinical dataset that were consistent with former studies using EUD or GED demonstrated that ADT was sensitive to cortical atrophy. The robustness against inter-individual variability of ADT was highlighted through simulation dataset. ADT showed a low and constant normalized difference between the depth of the simulated data and the calculated depth, whereas EUD and GED had high and variable differences. We suggest that ADT is more robust than EUD or GED and might be a useful alternative algorithm for measuring sulcal depth

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