Comprehensive Brain MRI Segmentation in High Risk Preterm Newborns

Most extremely preterm newborns exhibit cerebral atrophy/growth disturbances and white matter signal abnormalities on MRI at term-equivalent age. MRI brain volumes could serve as biomarkers for evaluating the effects of neonatal intensive care and predicting neurodevelopmental outcomes. This requires detailed, accurate, and reliable brain MRI segmentation methods. We describe our efforts to develop such methods in high risk newborns using a combination of manual and automated segmentation tools. After intensive efforts to accurately define structural boundaries, two trained raters independently performed manual segmentation of nine subcortical structures using axial T2-weighted MRI scans from 20 randomly selected extremely preterm infants. All scans were re-segmented by both raters to assess reliability. High intra-rater reliability was achieved, as assessed by repeatability and intra-class correlation coefficients (ICC range: 0.97 to 0.99) for all manually segmented regions. Inter-rater reliability was slightly lower (ICC range: 0.93 to 0.99). A semi-automated segmentation approach was developed that combined the parametric strengths of the Hidden Markov Random Field Expectation Maximization algorithm with non-parametric Parzen window classifier resulting in accurate white matter, gray matter, and CSF segmentation. Final manual correction of misclassification errors improved accuracy (similarity index range: 0.87 to 0.89) and facilitated objective quantification of white matter signal abnormalities. The semi-automated and manual methods were seamlessly integrated to generate full brain segmentation within two hours. This comprehensive approach can facilitate the evaluation of large cohorts to rigorously evaluate the utility of regional brain volumes as biomarkers of neonatal care and surrogate endpoints for neurodevelopmental outcomes

A four-dimensional probabilistic atlas of the human brain

The authors describe the development of a four-dimensional atlas and reference system that includes both macroscopic and microscopic information on structure and function of the human brain in persons between the ages of 18 and 90 years. Given the presumed large but previously unquantified degree of structural and functional variance among normal persons in the human population, the basis for this atlas and reference system is probabilistic. Through the efforts of the International Consortium for Brain Mapping (ICBM), 7,000 subjects will be included in the initial phase of database and atlas development. For each subject, detailed demographic, clinical, behavioral, and imaging information is being collected. In addition, 5,800 subjects will contribute DNA for the purpose of determining genotype-phenotype-behavioral correlations. The process of developing the strategies, algorithms, data collection methods, validation approaches, database structures, and distribution of results is described in this report. Examples of applications of the approach are described for the normal brain in both adults and children as well as in patients with schizophrenia. This project should provide new insights into the relationship between microscopic and macroscopic structure and function in the human brain and should have important implications in basic neuroscience, clinical diagnostics, and cerebral disorders

MRI quantification of multiple sclerosis pathology

Background: Multiple sclerosis (MS) is a chronic neuroinflammatory and neurodegenerative disease and a common cause of neurologic disability. MS pathology is characterized by demyelination, neuroaxonal loss and atrophy. Magnetic Resonance Imaging (MRI) is an essential tool in diagnosing and monitoring MS, but its clinical value could be even further expanded by more advanced and quantitative MRI methods, which may also provide additional pathophysiological insights. Purpose: The overall aim of this thesis was to quantify MS pathology using volumetric brain MRI, ultra-high field brain and cervical spinal cord MRI as well as a newly developed rapid myelin imaging technique in relation to cognitive and physical MS disability. Study I, a prospective 17-year longitudinal study of 37 MS participants with 23 age/sex- matched healthy controls for comparison at the last follow-up. Longitudinal volumetric brain 1.5 Tesla MRI during the second half of the study showed that lesion accumulation and corpus callosum atrophy were the most strongly associated neuroanatomical correlates of cognitive disability, with the lesion fraction being an independent predictor of cognitive performance 8.5 years later. Study II, a prospective cross-sectional study of 35 MS participants and 11 age-matched healthy controls using 3 and 7 Tesla MRI. The study demonstrated involvement of both grey and white matter in MS, not only the brain but also the cervical spinal cord, associated with MS disability. Lesions appeared in proximity to the cerebrospinal fluid (CSF), with special predilection to the periventricular and grey matter surrounding the central canal in secondary progressive MS. Study III, a prospective in vivo (71 MS participants and 21 age/sex-matched healthy controls) and ex vivo (brain tissue from 3 MS donors) study at 3 Tesla, showed that a new clinically approved and feasible rapid myelin imaging technique correlates well with myelin stainings and produces robust in vivo myelin quantification that is related to both current and future cognitive and physical disability in MS. Study IV, an in-depth topographical analysis based on Study III, mapped the distribution of demyelination, both in vivo and ex vivo, in the periventricular and perilesional regions of the brain. A gradient of demyelination with predominance near the CSF spaces was seen. Measures of clinical disability were consistently and more strongly associated with the myelin content in normal-appearing tissue compared to the intralesional myelin content. Conclusions: Lesions and atrophy contribute to cognitive and physical disability in MS but to a varying degree, likely dependent on the relative involvement of white vs. grey matter. Both focal lesions/demyelination as well as diffuse demyelination in normal-appearing white matter shows an apparent gradient from the CSF, which differ between relapsing-remitting and progressive MS subtypes/phases. The growing utility and clinical availability of advanced and quantitative MRI techniques holds promise for improved monitoring of MS pathology and likely represents a vital tool for assessing the efficacy of potential remyelinating/reparative therapies in MS

Developing advanced mathematical models for detecting abnormalities in 2D/3D medical structures.

Detecting abnormalities in two-dimensional (2D) and three-dimensional (3D) medical structures is among the most interesting and challenging research areas in the medical imaging field. Obtaining the desired accurate automated quantification of abnormalities in medical structures is still very challenging. This is due to a large and constantly growing number of different objects of interest and associated abnormalities, large variations of their appearances and shapes in images, different medical imaging modalities, and associated changes of signal homogeneity and noise for each object. The main objective of this dissertation is to address these problems and to provide proper mathematical models and techniques that are capable of analyzing low and high resolution medical data and providing an accurate, automated analysis of the abnormalities in medical structures in terms of their area/volume, shape, and associated abnormal functionality. This dissertation presents different preliminary mathematical models and techniques that are applied in three case studies: (i) detecting abnormal tissue in the left ventricle (LV) wall of the heart from delayed contrast-enhanced cardiac magnetic resonance images (MRI), (ii) detecting local cardiac diseases based on estimating the functional strain metric from cardiac cine MRI, and (iii) identifying the abnormalities in the corpus callosum (CC) brain structure—the largest fiber bundle that connects the two hemispheres in the brain—for subjects that suffer from developmental brain disorders. For detecting the abnormal tissue in the heart, a graph-cut mathematical optimization model with a cost function that accounts for the object’s visual appearance and shape is used to segment the the inner cavity. The model is further integrated with a geometric model (i.e., a fast marching level set model) to segment the outer border of the myocardial wall (the LV). Then the abnormal tissue in the myocardium wall (also called dead tissue, pathological tissue, or infarct area) is identified based on a joint Markov-Gibbs random field (MGRF) model of the image and its region (segmentation) map that accounts for the pixel intensities and the spatial interactions between the pixels. Experiments with real in-vivo data and comparative results with ground truth (identified by a radiologist) and other approaches showed that the proposed framework can accurately detect the pathological tissue and can provide useful metrics for radiologists and clinicians. To estimate the strain from cardiac cine MRI, a novel method based on tracking the LV wall geometry is proposed. To achieve this goal, a partial differential equation (PDE) method is applied to track the LV wall points by solving the Laplace equation between the LV contours of each two successive image frames over the cardiac cycle. The main advantage of the proposed tracking method over traditional texture-based methods is its ability to track the movement and rotation of the LV wall based on tracking the geometric features of the inner, mid-, and outer walls of the LV. This overcomes noise sources that come from scanner and heart motion. To identify the abnormalities in the CC from brain MRI, the CCs are aligned using a rigid registration model and are segmented using a shape-appearance model. Then, they are mapped to a simple unified space for analysis. This work introduces a novel cylindrical mapping model, which is conformal (i.e., one to one transformation and bijective), that enables accurate 3D shape analysis of the CC in the cylindrical domain. The framework can detect abnormalities in all divisions of the CC (i.e., splenium, rostrum, genu and body). In addition, it offers a whole 3D analysis of the CC abnormalities instead of only area-based analysis as done by previous groups. The initial classification results based on the centerline length and CC thickness suggest that the proposed CC shape analysis is a promising supplement to the current techniques for diagnosing dyslexia. The proposed techniques in this dissertation have been successfully tested on complex synthetic and MR images and can be used to advantage in many of today’s clinical applications of computer-assisted medical diagnostics and intervention

Morphological correlates to cognitive dysfunction in schizophrenia as studied with Bayesian regression

BACKGROUND: Relationships between cognitive deficits and brain morphological changes observed in schizophrenia are alternately explained by less gray matter in the brain cerebral cortex, by alterations in neural circuitry involving the basal ganglia, and by alteration in cerebellar structures and related neural circuitry. This work explored a model encompassing all of these possibilities to identify the strongest morphological relationships to cognitive skill in schizophrenia. METHODS: Seventy-one patients with schizophrenia and sixty-five healthy control subjects were characterized by neuropsychological tests covering six functional domains. Measures of sixteen brain morphological structures were taken using semi-automatic and fully manual tracing of MRI images, with the full set of measures completed on thirty of the patients and twenty controls. Group differences were calculated. A Bayesian decision-theoretic method identified those morphological features, which best explained neuropsychological test scores in the context of a multivariate response linear model with interactions. RESULTS: Patients performed significantly worse on all neuropsychological tests except some regarding executive function. The most prominent morphological observations were enlarged ventricles, reduced posterior superior vermis gray matter volumes, and increased putamen gray matter volumes in the patients. The Bayesian method associated putamen volumes with verbal learning, vigilance, and (to a lesser extent) executive function, while caudate volumes were associated with working memory. Vermis regions were associated with vigilance, executive function, and, less strongly, visuo-motor speed. Ventricular volume was strongly associated with visuo-motor speed, vocabulary, and executive function. Those neuropsychological tests, which were strongly associated to ventricular volume, showed only weak association to diagnosis, possibly because ventricular volume was regarded a proxy for diagnosis. Diagnosis was strongly associated with the other neuropsychological tests, implying that the morphological associations for these tasks reflected morphological effects and not merely group volumetric differences. Interaction effects were rarely associated, indicating that volumetric relationships to neuropsychological performance were similar for both patients and controls. CONCLUSION: The association of subcortical and cerebellar structures to verbal learning, vigilance, and working memory supports the importance of neural connectivity to these functions. The finding that a morphological indicator of diagnosis (ventricular volume) provided more explanatory power than diagnosis itself for visuo-motor speed, vocabulary, and executive function suggests that volumetric abnormalities in the disease are more important for cognition than non-morphological features

A proposal for a coordinated effort for the determination of brainwide neuroanatomical connectivity in model organisms at a mesoscopic scale

In this era of complete genomes, our knowledge of neuroanatomical circuitry remains surprisingly sparse. Such knowledge is however critical both for basic and clinical research into brain function. Here we advocate for a concerted effort to fill this gap, through systematic, experimental mapping of neural circuits at a mesoscopic scale of resolution suitable for comprehensive, brain-wide coverage, using injections of tracers or viral vectors. We detail the scientific and medical rationale and briefly review existing knowledge and experimental techniques. We define a set of desiderata, including brain-wide coverage; validated and extensible experimental techniques suitable for standardization and automation; centralized, open access data repository; compatibility with existing resources, and tractability with current informatics technology. We discuss a hypothetical but tractable plan for mouse, additional efforts for the macaque, and technique development for human. We estimate that the mouse connectivity project could be completed within five years with a comparatively modest budget.Comment: 41 page

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

Traumatic brain injury: a comparison of diffusion and volumetric magnetic resonance imaging measures

Cognitive impairment after traumatic brain injury remains hard to predict. This is partly because axonal injury, which is of fundamental importance, is difficult to measure clinically. Advances in MRI allow axonal injury to be detected after traumatic brain injury, but the most sensitive approach is unclear. Here, we compare the performance of diffusion tensor imaging, neurite orientation dispersion and density-imaging and volumetric measures of brain atrophy in the identification of white-matter abnormalities after traumatic brain injury. Thirty patients with moderate-severe traumatic brain injury in the chronic phase and 20 age-matched controls had T1-weighted and diffusion MRI. Neuropsychological tests of processing speed, executive functioning and memory were used to detect cognitive impairment. Extensive abnormalities in neurite density index and orientation dispersion index were observed, with distinct spatial patterns. Fractional anisotropy and mean diffusivity also indicated widespread abnormalities of white-matter structure. Neurite density index was significantly correlated with processing speed. Slower processing speed was also related to higher mean diffusivity in the corticospinal tracts. Lower white-matter volumes were seen after brain injury with greater effect sizes compared to diffusion metrics; however, volume was not sensitive to changes in cognitive performance. Volume was the most sensitive at detecting change between groups but was not specific for determining relationships with cognition. Abnormalities in fractional anisotropy and mean diffusivity were the most sensitive diffusion measures; however, neurite density index and orientation dispersion index may be more spatially specific. Lower neurite density index may be a useful metric for examining slower processing speed