Fast diffusion MRI based on sparse acquisition and reconstruction for long-term population imaging

Diffusion weighted magnetic resonance imaging (dMRI) is a unique MRI modality to probe the diffusive molecular transport in biological tissue. Due to its noninvasiveness and its ability to investigate the living human brain at submillimeter scale, dMRI is frequently performed in clinical and biomedical research to study the brain’s complex microstructural architecture. Over the last decades large prospective cohort studies have been set up with the aim to gain new insights into the development and progression of brain diseases across the life span and to discover biomarkers for disease prediction and potentially prevention. To allow for diverse brain imaging using different MRI modalities, stringent scan time limits are typically imposed in population imaging. Nevertheless, population studies aim to apply advanced and thereby time consuming dMRI protocols that deliver high quality data with great potential for future analysis. To allow for time-efficient but also versatile diffusion imaging, this thesis contributes to the investigation of accelerating diffusion spectrum imaging (DSI), an advanced dMRI technique that acquires imaging data with high intra-voxel resolution of tissue microstructure. Combining state-of-the-art parallel imaging and the theory of compressed sensing (CS) enables the acceleration of spatial encoding and diffusion encoding in dMRI. In this way, the otherwise long acquisition times in DSI can be reduced significantly. In this thesis, first, suitable q-space sampling strategies and basis functions are explored that fulfill the requirements of CS theory for accurate sparse DSI reconstruction. Novel 3D q-space sample distributions are investigated for CS-DSI. Moreover, conventional CS-DSI based on the discrete Fourier transform is compared for the first time to CS-DSI based on the continuous SHORE (simple harmonic oscillator based reconstruction and estimation) basis functions. Based on these findings, a CS-DSI protocol is proposed for application in a prospective cohort study, the Rhineland Study. A pilot study was designed and conducted to evaluate the CS-DSI protocol in comparison with state-of-the-art 3-shell dMRI and dedicated protocols for diffusion tensor imaging (DTI) and for the combined hindered and restricted model of diffusion (CHARMED). Population imaging requires processing techniques preferably with low computational cost to process and analyze the acquired big data within a reasonable time frame. Therefore, a pipeline for automated processing of CS-DSI acquisitions was implemented including both in-house developed and existing state-of-the-art processing tools. The last contribution of this thesis is a novel method for automatic detection and imputation of signal dropout due to fast bulk motion during the diffusion encoding in dMRI. Subject motion is a common source of artifacts, especially when conducting clinical or population studies with children, the elderly or patients. Related artifacts degrade image quality and adversely affect data analysis. It is, thus, highly desired to detect and then exclude or potentially impute defective measurements prior to dMRI analysis. Our proposed method applies dMRI signal modeling in the SHORE basis and determines outliers based on the weighted model residuals. Signal imputation reconstructs corrupted and therefore discarded measurements from the sparse set of inliers. This approach allows for fast and robust correction of imaging artifacts in dMRI which is essential to estimate accurate and precise model parameters that reflect the diffusive transport of water molecules and the underlying microstructural environment in brain tissue.Die diffusionsgewichtete Magnetresonanztomographie (dMRT) ist ein einzigartiges MRTBildgebungsverfahren, um die Diffusionsbewegung von Wassermolekülen in biologischem Gewebe zu messen. Aufgrund der Möglichkeit Schichtbilder nicht invasiv aufzunehmen und das lebende menschliche Gehirn im Submillimeter-Bereich zu untersuchen, ist die dMRT ein häufig verwendetes Bildgebungsverfahren in klinischen und biomedizinischen Studien zur Erforschung der komplexen mikrostrukturellen Architektur des Gehirns. In den letzten Jahrzehnten wurden große prospektive Kohortenstudien angelegt, um neue Einblicke in die Entwicklung und den Verlauf von Gehirnkrankheiten über die Lebenspanne zu erhalten und um Biomarker zur Krankheitserkennung und -vorbeugung zu bestimmen. Um durch die Verwendung unterschiedlicher MRT-Verfahren verschiedenartige Schichtbildaufnahmen des Gehirns zu ermöglich, müssen Scanzeiten typischerweise stark begrenzt werden. Dennoch streben Populationsstudien die Anwendung von fortschrittlichen und daher zeitintensiven dMRT-Protokollen an, um Bilddaten in hoher Qualität und mit großem Potential für zukünftige Analysen zu akquirieren. Um eine zeiteffizente und gleichzeitig vielseitige Diffusionsbildgebung zu ermöglichen, leistet diese Dissertation Beiträge zur Untersuchung von Beschleunigungsverfahren für die Bildgebung mittels diffusion spectrum imaging (DSI). DSI ist ein fortschrittliches dMRT-Verfahren, das Bilddaten mit hoher intra-voxel Auflösung der Gewebestruktur erhebt. Werden modernste Verfahren zur parallelen MRT-Bildgebung mit der compressed sensing (CS) Theorie kombiniert, ermöglicht dies eine Beschleunigung der räumliche Kodierung und der Diffusionskodierung in der dMRT. Dadurch können die ansonsten langen Aufnahmezeiten für DSI erheblich reduziert werden. In dieser Arbeit werden zuerst geeigenete Strategien zur Abtastung des q-space sowie Basisfunktionen untersucht, welche die Anforderungen der CS-Theorie für eine korrekte Signalrekonstruktion der dünnbesetzten DSI-Daten erfüllen. Neue 3D-Verteilungen von Messpunkten im q-space werden für die Verwendung in CS-DSI untersucht. Außerdem wird konventionell auf der diskreten Fourier-Transformation basierendes CS-DSI zum ersten Mal mit einem CS-DSI Verfahren verglichen, welches kontinuierliche SHORE (simple harmonic oscillator based reconstruction and estimation) Basisfunktionen verwendet. Aufbauend auf diesen Ergebnissen wird ein CS-DSI-Protokoll zur Anwendung in einer prospektiven Kohortenstudie, der Rheinland Studie, vorgestellt. Eine Pilotstudie wurde entworfen und durchgeführt, um das CS-DSI-Protokoll im Vergleich mit modernster 3-shell-dMRT und mit dedizierten Protokollen für diffusion tensor imaging (DTI) und für das combined hindered and restricted model of diffusion (CHARMED) zu evaluieren. Populationsbildgebung erfordert Prozessierungsverfahren mit möglichst geringem Rechenaufwand, um große akquirierte Datenmengen in einem angemessenen Zeitrahmen zu verarbeiten und zu analysieren. Dafür wurde eine Pipeline zur automatisierten Verarbeitung von CS-DSI-Daten implementiert, welche sowohl eigenentwickelte als auch bereits existierende moderene Verarbeitungsprogramme enthält. Der letzte Beitrag dieser Arbeit ist eine neue Methode zur automatischen Detektion und Imputation von Signalabfall, welcher durch schnelle Bewegungen während der Diffusionskodierung in der dMRT entsteht. Bewegungen der Probanden während der dMRT-Aufnahme sind eine häufige Ursache für Bildfehler, vor allem in klinischen oder Populationsstudien mit Kindern, alten Menschen oder Patienten. Diese Artefakte vermindern die Datenqualität und haben einen negativen Einfluss auf die Datenanalyse. Daher ist es das Ziel, fehlerhafte Messungen vor der dMRI-Analyse zu erkennen und dann auszuschließen oder wenn möglich zu ersetzen. Die vorgestellte Methode verwendet die SHORE-Basis zur dMRT-Signalmodellierung und bestimmt Ausreißer mit Hilfe von gewichteten Modellresidualen. Die Datenimputation rekonstruiert die unbrauchbaren und daher verworfenen Messungen mit Hilfe der verbleibenden, dünnbesetzten Menge an Messungen. Dieser Ansatz ermöglicht eine schnelle und robuste Korrektur von Bildartefakten in der dMRT, welche erforderlich ist, um korrekte und präzise Modellparameter zu schätzen, die die Diffusionsbewegung von Wassermolekülen und die zugrundeliegende Mikrostruktur des Gehirngewebes reflektieren

A 3D Printed Axon-Mimetic Diffusion MRI Phantom

Diffusion MRI is used to non-invasively characterize the microstructure of the brain. However, the accuracy of the characterization is difficult to verify because no other non-invasive imaging modality provides the same information. This thesis presents a novel 3D printed axon-mimetic (3AM) diffusion MRI phantom, a synthetic object designed to mimic the brain\u27s microstructure. The phantoms were characterized using microscopy, synchrotron micro-computed tomography, and diffusion MRI, and found to have sufficiently axon-mimetic properties to be useful as diffusion MRI phantoms. A set of phantoms designed to have anatomically realistic and complex fibre structures was used to test the response of diffusion MRI models of white matter to fibre orientation dispersion. All tested models were found to respond to orientation dispersion, but some robust metrics were identified. The studies in this thesis demonstrate that 3AM phantoms are a novel, flexible, and inexpensive tool for validating diffusion MRI models of white matter

Super Resolution of HARDI images Using Compressed Sensing Techniques

Effective techniques of inferring the condition of neural tracts in the brain is invaluable for clinicians and researchers towards investigation of neurological disorders in patients. It was not until the advent of diffusion Magnetic Resonance Imaging (dMRI), a noninvasive imaging method used to detect the diffusion of water molecules, that scientists have been able to assess the characteristics of cerebral diffusion in vivo. Among different dMRI methods, High Angular Resolution Diffusion Imaging (HARDI) is well known for striking a balance between ability to distinguish crossing neural fibre tracts while requiring a modest number of diffusion measurements (which is directly related to acquisition time). HARDI data provides insight into the directional properties of water diffusion in cerebral matter as a function of spatial coordinates. Ideally, one would be interested in having this information available at fine spatial resolution while minimizing the probing along different spatial orientations (so as to minimize the acquisition time). Unfortunately, availability of such datasets in reasonable acquisition times are hindered by limitations in current hardware and scanner protocols. On the other hand, post processing techniques prove promising in increasing the effective spatial resolution, allowing more detailed depictions of cerebral matter, while keeping the number of diffusion measurements within a feasible range. In light of the preceding developments, the main purpose of this research is to look into super resolution of HARDI data, using the modern theory of compressed sensing. The method proposed in this thesis allows an accurate approximation of HARDI signals at a higher spatial resolution compared to data obtained with a typical scanner. At the same time, ideas for reducing the number of diffusion measurements in the angular domain to improve the acquisition time are explored. Accordingly, the novel method of applying two distinct compressed sensing approaches in both spatial and angular domain, and combining them into a single framework for performing super resolution forms the main contribution provided by this thesis

A novel mechanism of contrast in MRI: pseudo super-diffusion of water molecules unveils microstructural details in biological tissues

The goal of this work is to investigate the properties of the contrast provided by Anomalous Diffusion (AD) γ-imaging technique and to test its potential in detecting tissue microstructure. The collateral purpose is to implement this technique by optimizing data acquisition and data processing, with the long term perspective of adoption in massive in vitro, in vivo and clinical studies. The AD γ-imaging technique is a particular kind of Diffusion Weighted- Magnetic Resonance Imaging (DW-MRI). It represents a refinement of conventionally used DW-MRI methods, sharing with them the advantage of being non invasive, since it uses water as an endogenous contrast agent. Besides, it is more suitable to the study of complex tissues, because it is based on a theoretical model that overcomes the simplistic Gaussian assumption. While the Gaussian assumption predicates the linearity between the average molecular displacement of water and the diffusing time, as in case of diffusion in isotropic, homogeneous and infinite environments, a number of experiments performed in vitro and in vivo on both animals and humans showed an anomalous behavior of water molecules, with a non linear relation between the distance travelled and the elapsed time. In particular, the γ-parameter quantifies water pseudo super-diffusion, a peculiarity due to the fact that water diffusion occurs in multi-compartments and it is probed by means of MRI. In fact, a restricted diffusion is rather predicted for water diffusing in biological tissues. Recently, the trick that allows to make the traditional DW-MRI acquisition sequence suitable for pseudo super-diffusion quantification has been unveiled, and in short it consists in performing DW experiments varying the diffusion gradient strengths, at a constant diffusive time. The γ-parameter is extracted by fitting DW-data to a stretched-exponential function. Finally, probing water diffusion in different directions allows to reconstruct a γ-tensor, with scalar invariants that quantify the entity of AD and its anisotropy in a given volume element. In vitro results on inert materials revealed that γ correlates with internal gradients arising from magnetic susceptibility differences (Δ) between neighboring compartments, and that it reflects the multi-compartmentalization of the space explored by diffusing molecules. Furthermore, values of γ compatible with a description of super-diffusive motion were found. This anomaly can be explained considering that the presence of Δ induce an additional attenuation to the signal, simulating a pseudo super-diffusion. Finally, In vivo results on human brain showed that γ is more effective in discriminating among different brain regions compared to conventional DWMRI parameters. These studies suggest that the contrast provided by AD γ-imaging is influenced by an interplay of two factors, Δ -effects on one hand, multicompartmentalization on the other hand, through which γ could reflect tissue microstructure. With the aim to shed some light on this issue I performed AD γ-imaging in excised mouse spinal cord (MSC) at 9.4 T and healthy human brain at 3.0 T. The adoption of MSC was motivated by its current use in studies of demyelination due to an induced pathology that mimics Multiple Sclerosis alterations, and by its simplified geometry. I acquired DW-data with parameters optimized for the particular system chosen: the MSC was scanned along 3 orthogonal directions, thus an apparent γ was derived; for the in vivo studies I used more directions and I extracted a γ-tensor. I found that γ and its anisotropy reflected the microstructure of spinal cord tracts (such as the axon diameters and the axonal density). I investigated both in MSC and human brain the relation between γ and the rate of relaxation (R2*), a parameter well-known to reflect Δ, and found significant linear correlations. Because of this γ was able to differentiate white matter regions on the basis of their spatial orientation, and gray matter regions on the basis of their intrinsic iron content in human brain imaged at 3.0 T. These results suggest that AD γ-imaging could be an alternative or complementary technique to DW-MRI in the field of neuroscience. Indeed it could be useful for the assessment of the bulk susceptibility inhomogeneity, which reflects iron deposition, the hallmark of several neurodegenerative diseases. The part of this thesis work concerning the in vivo experiment in human brain gave rise to a paper published on NeuroImage, a relevant scientific journal in the field of MRI applied to brain investigation

Development of High Angular Resolution Diffusion Imaging Analysis Paradigms for the Investigation of Neuropathology

Diffusion weighted magnetic resonance imaging (DW-MRI), provides unique insight into the microstructure of neural white matter tissue, allowing researchers to more fully investigate white matter disorders. The abundance of clinical research projects incorporating DW-MRI into their acquisition protocols speaks to the value this information lends to the study of neurological disease. However, the most widespread DW-MRI technique, diffusion tensor imaging (DTI), possesses serious limitations which restrict its utility in regions of complex white matter. Fueled by advances in DW-MRI acquisition protocols and technologies, a group of exciting new DW-MRI models, developed to address these concerns, are now becoming available to clinical researchers. The emergence of these new imaging techniques, categorized as high angular resolution diffusion imaging (HARDI), has generated the need for sophisticated computational neuroanatomic techniques able to account for the high dimensionality and structure of HARDI data. The goal of this thesis is the development of such techniques utilizing prominent HARDI data models. Specifically, methodologies for spatial normalization, population atlas building and structural connectivity have been developed and validated. These methods form the core of a comprehensive analysis paradigm allowing the investigation of local white matter microarcitecture, as well as, systemic properties of neuronal connectivity. The application of this framework to the study of schizophrenia and the autism spectrum disorders demonstrate its sensitivity sublte differences in white matter organization, as well as, its applicability to large population DW-MRI studies

Homogeneity based segmentation and enhancement of Diffusion Tensor Images : a white matter processing framework

In diffusion magnetic resonance imaging (DMRI) the Brownian motion of the water molecules, within biological tissue, is measured through a series of images. In diffusion tensor imaging (DTI) this diffusion is represented using tensors. DTI describes, in a non-invasive way, the local anisotropy pattern enabling the reconstruction of the nervous fibers - dubbed tractography. DMRI constitutes a powerful tool to analyse the structure of the white matter within a voxel, but also to investigate the anatomy of the brain and its connectivity. DMRI has been proved useful to characterize brain disorders, to analyse the differences on white matter and consequences in brain function. These procedures usually involve the virtual dissection of white matters tracts of interest. The manual isolation of these bundles requires a great deal of neuroanatomical knowledge and can take up to several hours of work. This thesis focuses on the development of techniques able to automatically perform the identification of white matter structures. To segment such structures in a tensor field, the similarity of diffusion tensors must be assessed for partitioning data into regions, which are homogeneous in terms of tensor characteristics. This concept of tensor homogeneity is explored in order to achieve new methods for segmenting, filtering and enhancing diffusion images. First, this thesis presents a novel approach to semi-automatically define the similarity measures that better suit the data. Following, a multi-resolution watershed framework is presented, where the tensor field’s homogeneity is used to automatically achieve a hierarchical representation of white matter structures in the brain, allowing the simultaneous segmentation of different structures with different sizes. The stochastic process of water diffusion within tissues can be modeled, inferring the homogeneity characteristics of the diffusion field. This thesis presents an accelerated convolution method of diffusion images, where these models enable the contextual processing of diffusion images for noise reduction, regularization and enhancement of structures. These new methods are analysed and compared on the basis of their accuracy, robustness, speed and usability - key points for their application in a clinical setting. The described methods enrich the visualization and exploration of white matter structures, fostering the understanding of the human brain

What's new and what's next in diffusion MRI preprocessing

Diffusion MRI (dMRI) provides invaluable information for the study of tissue microstructure and brain connectivity, but suffers from a range of imaging artifacts that greatly challenge the analysis of results and their interpretability if not appropriately accounted for. This review will cover dMRI artifacts and preprocessing steps, some of which have not typically been considered in existing pipelines or reviews, or have only gained attention in recent years: brain/skull extraction, B-matrix incompatibilities w.r.t the imaging data, signal drift, Gibbs ringing, noise distribution bias, denoising, between- and within-volumes motion, eddy currents, outliers, susceptibility distortions, EPI Nyquist ghosts, gradient deviations, bias fields, and spatial normalization. The focus will be on “what’s new” since the notable advances prior to and brought by the Human Connectome Project (HCP), as presented in the predecessing issue on “Mapping the Connectome” in 2013. In addition to the development of novel strategies for dMRI preprocessing, exciting progress has been made in the availability of open source tools and reproducible pipelines, databases and simulation tools for the evaluation of preprocessing steps, and automated quality control frameworks, amongst others. Finally, this review will consider practical considerations and our view on “what’s next” in dMRI preprocessing

Imaging brain microstructure with diffusion MRI: practicality and applications

This article gives an overview of microstructure imaging of the brain with diffusion MRI and reviews the state of the art. The microstructure-imaging paradigm aims to estimate and map microscopic properties of tissue using a model that links these properties to the voxel scale MR signal. Imaging techniques of this type are just starting to make the transition from the technical research domain to wide application in biomedical studies. We focus here on the practicalities of both implementing such techniques and using them in applications. Specifically, the article summarizes the relevant aspects of brain microanatomy and the range of diffusion-weighted MR measurements that provide sensitivity to them. It then reviews the evolution of mathematical and computational models that relate the diffusion MR signal to brain tissue microstructure, as well as the expanding areas of application. Next we focus on practicalities of designing a working microstructure imaging technique: model selection, experiment design, parameter estimation, validation, and the pipeline of development of this class of technique. The article concludes with some future perspectives on opportunities in this topic and expectations on how the field will evolve in the short-to-medium term

Improved Quantification of Connectivity in Human Brain Mapping

Diffusion magnetic resonance imaging (dMRI) is an advanced MRI methodology that can be used to probe the microstructure of biological tissue. dMRI can provide orientation information by modeling the process of water diffusion in white matter. This thesis presents contributions in three areas of diffusion imaging technology: diffusion reconstruction, quantification, and validation of derived metrics. It presents a novel reconstruction method by combining generalized q-sampling imaging, spherical harmonic basis functions and constrained spherical deconvolution methods to estimate the fiber orientation distribution function (ODF). This method provides improved spatial localization of brain nuclei and fiber tract separation. A novel diffusion anisotropy metric is presented that provides anatomically interpretable measurements of tracts that are robust in crossing areas of the brain. The metric, directional Axonal Volume (dAV) provides an estimate of directional water content of the tract based on the (ODF) and proton density map. dAV is a directionally sensitive metric and can separate anisotropic water content for each fiber population, providing a quantification in milliliters of water. A method is provided to map voxel-based dAV onto tracts that is not confounded by crossing areas and follows the tract morphology. This work introduces a novel textile based hollow fiber anisotropic phantom (TABIP) for validation of reconstruction and quantification methods. This provides a ground truth reference for axonal scale water tubular structures arranged in various anatomical configurations, crossing and mixing patterns. Analysis shows that: 1) the textile tracts are identifiable with scans used in human imaging and produced tracts and voxel metrics in the range of human tissue; 2) the current methods could resolve crossing at 90o and 45o but not 30o; 3) dAV/NODDI model closely matches (r=0.95) the number of fibers whereas conventional metrics poorly match (i.e., FA r=0.32). This work represents a new accurate quantification of axonal water content through diffusion imaging. dAV shows promise as a new anatomically interpretable metric of axonal connectivity that is not confounded by factors such as axon dispersion, crossing and local isotropic water content. This will provide better anatomical mapping of white matter and potentially improve the detection of axonal tract pathology