Diffeomorphic Metric Mapping of High Angular Resolution Diffusion Imaging based on Riemannian Structure of Orientation Distribution Functions

In this paper, we propose a novel large deformation diffeomorphic registration algorithm to align high angular resolution diffusion images (HARDI) characterized by orientation distribution functions (ODFs). Our proposed algorithm seeks an optimal diffeomorphism of large deformation between two ODF fields in a spatial volume domain and at the same time, locally reorients an ODF in a manner such that it remains consistent with the surrounding anatomical structure. To this end, we first review the Riemannian manifold of ODFs. We then define the reorientation of an ODF when an affine transformation is applied and subsequently, define the diffeomorphic group action to be applied on the ODF based on this reorientation. We incorporate the Riemannian metric of ODFs for quantifying the similarity of two HARDI images into a variational problem defined under the large deformation diffeomorphic metric mapping (LDDMM) framework. We finally derive the gradient of the cost function in both Riemannian spaces of diffeomorphisms and the ODFs, and present its numerical implementation. Both synthetic and real brain HARDI data are used to illustrate the performance of our registration algorithm

Investigating brain connectivity heritability in a twin study using diffusion imaging data

Heritability of brain anatomical connectivity has been studied with diffusion-weighted imaging (DWI) mainly by modeling each voxel's diffusion pattern as a tensor (e.g., to compute fractional anisotropy), but this method cannot accurately represent the many crossing connections present in the brain. We hypothesized that different brain networks (i.e., their component fibers) might have different heritability and we investigated brain connectivity using High Angular Resolution Diffusion Imaging (HARDI) in a cohort of twins comprising 328 subjects that included 70 pairs of monozygotic and 91 pairs of dizygotic twins. Water diffusion was modeled in each voxel with a Fiber Orientation Distribution (FOD) function to study heritability for multiple fiber orientations in each voxel. Precision was estimated in a test-retest experiment on a sub-cohort of 39 subjects. This was taken into account when computing heritability of FOD peaks using an ACE model on the monozygotic and dizygotic twins. Our results confirmed the overall heritability of the major white matter tracts but also identified differences in heritability between connectivity networks. Inter-hemispheric connections tended to be more heritable than intra-hemispheric and cortico-spinal connections. The highly heritable tracts were found to connect particular cortical regions, such as medial frontal cortices, postcentral, paracentral gyri, and the right hippocampus

Methods for improved mapping of brain lesion connectivity

Recent advances over the past two decades in neuroimaging methods have enabled us to map the connectivity of the brain. In parallel, pathophysiological models of brain disease have shifted from an emphasis on understanding pathology in specific brain regions to characterizing disruptions to interconnected neural networks. Nevertheless, these recent methods for mapping brain connectivity are still under development. Every step of the mapping process becomes a potential source for additional error due to noise or artifacts that could impact final analyses. Segmentation, parcellation, registration, and tractography are some of the steps where this occurs. Moreover, mapping the connectivity in a brain lesion is even more susceptible to errors in these steps. In this body of work, I describe multiple new methods for improving the accuracy of mapping lesion connectivity by reducing errors at the tractography stage which is the most error prone stage. First, we develop an approach for directly normalizing streamlines into a template space that avoids performing tractography in the normalized template space, reducing the error of connectomes constructed in the template space with respect to the ground truth native space connectome. Second, we develop a rapid approach for performing shortest path tractography and constructing shortest path probability weighted connectomes which increases the connection specificity relative to local streamline tracking approaches. We then demonstrate how our shortest path tractography approach can be used construct a disconnectome, a connectivity map of the proportion of connections lost due to intersecting a lesion. We then develop a fast, greedy graph-theoretic algorithm that extracts the maximally disconnected subgraph containing brain regions with the greatest shared loss of connectivity. Finally, we demonstrate how combining methods from diffusion based image inpainting and optimal estimation can be used to restore or inpaint corrupted fiber diffusion models in lesioned white matter tissue, enabling tractography and the study of lesion connectivity and modeling of microstructural measures in the patient’s native space

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

Effect of different spatial normalization approaches on tractography and structural brain networks

To facilitate the comparison of white matter morphologic connectivity across target populations, it is invaluable to map the data to a standardized neuroanatomical space. Here, we evaluated direct streamline normalization (DSN), where the warping was applied directly to the streamlines, with two publically available approaches that spatially normalize the diffusion data and then reconstruct the streamlines. Prior work has shown that streamlines generated after normalization from reoriented diffusion data do not reliably match the streamlines generated in native space. To test the impact of these different normalization methods on quantitative tractography measures, we compared the reproducibility of the resulting normalized connectivity matrices and network metrics with those originally obtained in native space. The two methods that reconstruct streamlines after normalization led to significant differences in network metrics with large to huge standardized effect sizes, reflecting a dramatic alteration of the same subject’s native connectivity. In contrast, after normalizing with DSN we found no significant difference in network metrics compared with native space with only very small-to-small standardized effect sizes. DSN readily outperformed the other methods at preserving native space connectivity and introduced novel opportunities to define connectome networks without relying on gray matter parcellations. Direct streamline normalization (DSN) directly warps the streamlines into any template space by using the transformations output from Advanced Normalization Tools (ANTs). DSN overcomes the limitations of diffusion weighted images (DWI) spatial normalization. It allows DWIs to be acquired with any desired sampling scheme. Fiber orientation distributions (FODs) or orientation distribution functions (ODFs) can also be reconstructed using any desired method and streamlines generated using any algorithm. Most importantly, it avoids the problem of generating tracts from FODs or ODFs that have become distorted because of spatial normalization. Our results show that DSN has minimal influence on tractography measures such as tract count and structure and does not significantly alter structural networks with only very small to small effect sizes. We have developed a framework in Python that works with most diffusion software platforms. It is available at http://github.com/clintg6/DSN

How our brains are wired: Are the applications of diffusion imaging useful given the current limitations?

Diffusion imaging (DI) enables researchers to study white matter (WM) pathways in the human brain in-vivo by labelling water molecules and measuring their diffusion into different directions. Connectivity patterns are inferred assuming that water diffuses rather along than across fibre bundles. This paper introduces the concept of DI, addresses suitable applications and evaluates gains versus limitations. Common applications are (1) generating WM atlases, (2) mapping connectional models of functionally subdivided brain regions, (3) linking disorders to connectivity abnormalities, (4) verifying WM pathways from animal studies, (5) linking personality traits to particular connectivity patterns, (6) measuring structural changes resulting from experience or ageing and (7) presurgical planning. Despite limitations like the moderate spatial resolution, or – more fundamentally – the lack of a gold standard and the kissing/crossing problem, DI can be regarded as a useful tool if researchers choose methods carefully and consider the known limitations

Genetic Contributions to the Midsagittal Area of the Corpus Callosum

The degree to which genes and environment determine variations in brain structure and function is fundamentally important to understanding normal and disease-related patterns of neural organization and activity. We studied genetic contributions to the midsagittal area of the corpus callosum (CC) in pedigreed baboons (68 males, 112 females) to replicate findings of high genetic contribution to that area of the CC reported in humans, and to determine if the heritability of the CC midsagittal area in adults was modulated by fetal development rate. Measurements of callosal area were obtained from high-resolution MRI scans. Heritability was estimated from pedigree-based maximum likelihood estimation of genetic and non-genetic variance components as implemented in Sequential Oligogenic Linkage Analysis Routines (SOLAR). Our analyses revealed significant heritability for the total area of the CC and all of its subdivisions, with h2 = .46 for the total CC, and h 2 = .54, .37, .62, .56, and .29 for genu, anterior midbody, medial midbody, posterior midbody and splenium, respectively. Genetic correlation analysis demonstrated that the individual subdivisions shared between 41% and 98% of genetic variability. Combined with previous research reporting high heritability of other brain structures in baboons, these results reveal a consistent pattern of high heritability for brain morphometric measures in baboons