Bayesian generative learning of brain and spinal cord templates from neuroimaging datasets

In the field of neuroimaging, Bayesian modelling techniques have been largely adopted and recognised as powerful tools for the purpose of extracting quantitative anatomical and functional information from medical scans. Nevertheless the potential of Bayesian inference has not yet been fully exploited, as many available tools rely on point estimation techniques, such as maximum likelihood estimation, rather than on full Bayesian inference. The aim of this thesis is to explore the value of approximate learning schemes, for instance variational Bayes, to perform inference from brain and spinal cord MRI data. The applications that will be explored in this work mainly concern image segmentation and atlas construction, with a particular emphasis on the problem of shape and intensity prior learning, from large training data sets of structural MR scans. The resulting computational tools are intended to enable integrated brain and spinal cord morphometric analyses, as opposed to the approach that is most commonly adopted in neuroimaging, which consists in optimising separate tools for brain and spine morphometrics

Quantifying neurodegeneration of the cervical cord and brain in degenerative cervical myelopathy : A multicentre study using quantitative magnetic resonance imaging

BACKGROUND AND PURPOSE Simultaneous assessment of neurodegeneration in both the cervical cord and brain across multiple centres can enhance the effectiveness of clinical trials. Thus, this study aims to simultaneously assess microstructural changes in the cervical cord and brain above the stenosis in degenerative cervical myelopathy (DCM) using quantitative magnetic resonance imaging (MRI) in a multicentre study. METHODS We applied voxelwise analysis with a probabilistic brain/spinal cord template embedded in statistical parametric mappin (SPM-BSC) to process multi parametric mapping (MPM) including effective transverse relaxation rate (R2*), longitudinal relaxation rate (R1), and magnetization transfer (MT), which are indirectly sensitive to iron and myelin content. Regression analysis was conducted to establish associations between neurodegeneration and clinical impairment. Thirty-eight DCM patients (mean age ± SD = 58.45 ± 11.47 years) and 38 healthy controls (mean age ± SD = 41.18 ± 12.75 years) were recruited at University Hospital Balgrist, Switzerland and Toronto Western Hospital, Canada. RESULTS Remote atrophy was observed in the cervical cord (p = 0.002) and in the left thalamus (0.026) of the DCM group. R1 was decreased in the periaqueductal grey matter (p = 0.014), thalamus (p = 0.001), corpus callosum (p = 0.0001), and cranial corticospinal tract (p = 0.03). R2* was increased in the primary somatosensory cortices (p = 0.008). Sensory impairments were associated with increased iron-sensitive R2* in the thalamus and periaqueductal grey matter in DCM. CONCLUSIONS Simultaneous assessment of the spinal cord and brain revealed DCM-induced demyelination, iron deposition, and atrophy. The extent of remote neurodegeneration was associated with sensory impairment, highlighting the intricate and expansive nature of microstructural neurodegeneration in DCM, reaching beyond the stenosis level

Computational modelling of imaging markers to support the diagnosis and monitoring of multiple sclerosis

Multiple sclerosis is a leading cause of neurological disability in young adults which affects more than 2.5 million people worldwide. An important substrate of disability accrual is the loss of neurons and connections between them (neurodegeneration) which can be captured by serial brain imaging, especially in the cerebral grey matter. In this thesis in four separate subprojects, I aimed to assess the strength of imaging-derived grey matter volume as a biomarker in the diagnosis, predicting the evolution of multiple sclerosis, and developing a staging system to stratify patients. In total, I retrospectively studied 1701 subjects, of whom 1548 had longitudinal brain imaging data. I used advanced computational models to investigate cross-sectional and longitudinal datasets. In the cross-sectional study, I demonstrated that grey matter volumes could distinguish multiple sclerosis from another demyelinating disorder (neuromyelitis optica) with an accuracy of 74%. In longitudinal studies, I showed that over time the deep grey matter nuclei had the fastest rate of volume loss (up to 1.66% annual loss) across the brain regions in multiple sclerosis. The volume of the deep grey matter was the strongest predictor of disability progression. I found that multiple sclerosis affects different brain areas with a specific temporal order (or sequence) that starts with the deep grey matter nuclei, posterior cingulate cortex, precuneus, and cerebellum. Finally, with multivariate mechanistic and causal modelling, I showed that brain volume loss causes disability and cognitive worsening which can be delayed with a potential neuroprotective treatment (simvastatin). This work provides conclusive evidence that grey matter volume loss affects some brain regions more severely, can predict future disability progression, can be used as an outcome measure in phase II clinical trials, and causes clinical and cognitive worsening. This thesis also provides a subject staging system based on which patients can be scored during multiple sclerosis

Deep learning in medical image registration: introduction and survey

Image registration (IR) is a process that deforms images to align them with respect to a reference space, making it easier for medical practitioners to examine various medical images in a standardized reference frame, such as having the same rotation and scale. This document introduces image registration using a simple numeric example. It provides a definition of image registration along with a space-oriented symbolic representation. This review covers various aspects of image transformations, including affine, deformable, invertible, and bidirectional transformations, as well as medical image registration algorithms such as Voxelmorph, Demons, SyN, Iterative Closest Point, and SynthMorph. It also explores atlas-based registration and multistage image registration techniques, including coarse-fine and pyramid approaches. Furthermore, this survey paper discusses medical image registration taxonomies, datasets, evaluation measures, such as correlation-based metrics, segmentation-based metrics, processing time, and model size. It also explores applications in image-guided surgery, motion tracking, and tumor diagnosis. Finally, the document addresses future research directions, including the further development of transformers

Generative Models for Preprocessing of Hospital Brain Scans

I will in this thesis present novel computational methods for processing routine clinical brain scans. Such scans were originally acquired for qualitative assessment by trained radiologists, and present a number of difficulties for computational models, such as those within common neuroimaging analysis software. The overarching objective of this work is to enable efficient and fully automated analysis of large neuroimaging datasets, of the type currently present in many hospitals worldwide. The methods presented are based on probabilistic, generative models of the observed imaging data, and therefore rely on informative priors and realistic forward models. The first part of the thesis will present a model for image quality improvement, whose key component is a novel prior for multimodal datasets. I will demonstrate its effectiveness for super-resolving thick-sliced clinical MR scans and for denoising CT images and MR-based, multi-parametric mapping acquisitions. I will then show how the same prior can be used for within-subject, intermodal image registration, for more robustly registering large numbers of clinical scans. The second part of the thesis focusses on improved, automatic segmentation and spatial normalisation of routine clinical brain scans. I propose two extensions to a widely used segmentation technique. First, a method for this model to handle missing data, which allows me to predict entirely missing modalities from one, or a few, MR contrasts. Second, a principled way of combining the strengths of probabilistic, generative models with the unprecedented discriminative capability of deep learning. By introducing a convolutional neural network as a Markov random field prior, I can model nonlinear class interactions and learn these using backpropagation. I show that this model is robust to sequence and scanner variability. Finally, I show examples of fitting a population-level, generative model to various neuroimaging data, which can model, e.g., CT scans with haemorrhagic lesions

Statistical Computing on Non-Linear Spaces for Computational Anatomy

International audienceComputational anatomy is an emerging discipline that aims at analyzing and modeling the individual anatomy of organs and their biological variability across a population. However, understanding and modeling the shape of organs is made difficult by the absence of physical models for comparing different subjects, the complexity of shapes, and the high number of degrees of freedom implied. Moreover, the geometric nature of the anatomical features usually extracted raises the need for statistics on objects like curves, surfaces and deformations that do not belong to standard Euclidean spaces. We explain in this chapter how the Riemannian structure can provide a powerful framework to build generic statistical computing tools. We show that few computational tools derive for each Riemannian metric can be used in practice as the basic atoms to build more complex generic algorithms such as interpolation, filtering and anisotropic diffusion on fields of geometric features. This computational framework is illustrated with the analysis of the shape of the scoliotic spine and the modeling of the brain variability from sulcal lines where the results suggest new anatomical findings

Automated Lesion Detection-Based Symptome Mapping In Patients With Sensory Deficits

The human brain is the most complex structure known. With its high number of cells, number of connections and number of pathways it is the source of every thought in the world. It consumes 25% of our oxygen and suffers very fast from a disruption of its supply. An acute event, like a stroke, results in rapid dysfunction referable to the affected area. A few minutes without oxygen and neuronal cells die and subsequently degenerate. Changes in the brains incoming blood flow alternate the anatomy and physiology of the brain. All stroke events leave behind a brain tissue lesion. To rapidly react and improve the prediction of outcome in stroke patients, accurate lesion detection and reliable lesion-based function correlation would be very helpful. With a number of neuroimaging and clinical data of cerebral injured patients this study aims to investigate correlations of structural lesion locations with sensory functions

Deformable Image Registration in the Analysis of Multiple Sclerosis

In medical image analysis, image registration is the task of finding corresponding features in two or more images, and using them to solve for the transformation that best aligns the images. Knowing the alignment allows information, such as landmarks and functional metrics, to be easily transferred between images, and allows them to be analyzed together. This dissertation focuses on the development of deformable image registration techniques for the analysis of multiple sclerosis (MS), a neurodegenerative disease that damages the myelin sheath of nervous tissue. MS is known to affect the entire central nervous system (CNS), and can result in the loss of sensorimotor control, cognition, and vision. Hence, the four primary contributions of this dissertation are on the development and application of deformable image registration in the three areas of the CNS that are most currently studied for MS -- the spinal cord, the retina, and the brain. First, for spinal cord magnetic resonance imaging (MRI), an approach is presented that uses deformable registration to provide atlas priors for automatic topology-preserving segmentation of the spinal cord and cerebrospinal fluid. The method shows high accuracy and robustness when compared to manual raters, and allows spinal cord atrophy to be analyzed on large datasets without manual segmentations. Second, for spinal cord diffusion tensor imaging, a pipeline is presented that uses deformable registration to correct for susceptibility distortions in the images. The pipeline allows for accurate computation of spinal cord diffusion metrics, which are shown to be significantly correlated with clinical measures of sensorimotor function and disability levels. Third, for optical coherence tomography (OCT) of the retina, a deformable registration technique is presented that constrains the transformation to follow the OCT acquisition geometry. 3D voxel-based analysis using the algorithm found significant differences between healthy and MS cohorts in regions of the retina that is consistent with previous findings using 2D analysis. Lastly, for brain MRI, a multi-channel registration framework is presented that can use distance transforms and image synthesis to improve registration accuracy. Together, these techniques have enabled several types of analysis that were previously unavailable for the study of MS