Imaging laminar structures in the gray matter with diffusion MRI

Direct segmentation of cortical cytoarchitectonic domains using ultra-high-resolution whole-brain diffusion MRI

We assess the potential of detecting cortical laminar patterns and areal borders by directly clustering voxel values of microstructural parameters derived from high-resolution mean apparent propagator (MAP) magnetic resonance imaging (MRI), as an alternative to conventional template-warping-based cortical parcellation methods. We acquired MAP-MRI data with 200μm resolution in a fixed macaque monkey brain. To improve the sensitivity to cortical layers, we processed the data with a local anisotropic Gaussian filter determined voxel-wise by the plane tangent to the cortical surface. We directly clustered all cortical voxels using only the MAP-derived microstructural imaging biomarkers, with no information regarding their relative spatial location or dominant diffusion orientations. MAP-based 3D cytoarchitectonic segmentation revealed laminar patterns similar to those observed in the corresponding histological images. Moreover, transition regions between these laminar patterns agreed more accurately with histology than the borders between cortical areas estimated using conventional atlas/template-warping cortical parcellation. By cross-tabulating all cortical labels in the atlas- and MAP-based segmentations, we automatically matched the corresponding MAP-derived clusters (i.e., cytoarchitectonic domains) across the left and right hemispheres. Our results demonstrate that high-resolution MAP-MRI biomarkers can effectively delineate three-dimensional cortical cytoarchitectonic domains in single individuals. Their intrinsic tissue microstructural contrasts enable the construction of whole-brain mesoscopic cortical atlases.

Selective Medial Septum Lesions in Healthy Rats Induce Longitudinal Changes in Microstructure of Limbic Regions, Behavioral Alterations, and Increased Susceptibility to Status Epilepticus

Septo-hippocampal pathway, crucial for physiological functions and involved in epilepsy. Clinical monitoring during epileptogenesis is complicated. We aim to evaluate tissue changes after lesioning the medial septum (MS) of normal rats and assess how the depletion of specific neuronal populations alters the animals' behavior and susceptibility to establishing a pilocarpine-induced status epilepticus. Male Sprague-Dawley rats were injected into the MS with vehicle or saporins (to deplete GABAergic or cholinergic neurons; n = 16 per group). Thirty-two animals were used for diffusion tensor imaging (DTI); scanned before surgery and 14 and 49 days post-injection. Fractional anisotropy and apparent diffusion coefficient were evaluated in the fimbria, dorsal hippocampus, ventral hippocampus, dorso-medial thalamus, and amygdala. Between scans 2 and 3, animals were submitted to diverse behavioral tasks. Stainings were used to analyze tissue alterations. Twenty-four different animals received pilocarpine to evaluate the latency and severity of the status epilepticus 2 weeks after surgery. Additionally, eight different animals were only used to evaluate the neuronal damage inflicted on the MS 1 week after the molecular surgery. Progressive changes in DTI parameters in both white and gray matter structures of the four evaluated groups were observed. Behaviorally, the GAT1-saporin injection impacted spatial memory formation, while 192-IgG-saporin triggered anxiety-like behaviors. Histologically, the GABAergic toxin also induced aberrant mossy fiber sprouting, tissue damage, and neuronal death. Regarding the pilocarpine-induced status epilepticus, this agent provoked an increased mortality rate. Selective septo-hippocampal modulation impacts the integrity of limbic regions crucial for certain behavioral skills and could represent a precursor for epilepsy development.

Charting Cortical-Layer Specific Area Boundaries Using Gibbs' Ringing Attenuated T1w/T2w-FLAIR Myelin MRI

Cortical areas have traditionally been defined by their distinctive layer cyto- and/or myelo- architecture using postmortem histology. Recent studies have delineated many areas by measuring overall cortical myelin content and its spatial gradients using the T1w/T2w ratio MRI in living primates, including humans. While T1w/T2w studies of areal transitions might benefit from using the layer profile of this myelin-related contrast, a significant confound is Gibbs' ringing artefact, which produces signal fluctuations resembling cortical layers. Here, we address these issues with a novel approach using cortical layer thickness-adjusted T1w/T2w-FLAIR imaging, which effectively cancels out Gibbs' ringing artefacts while enhancing intracortical myelin contrast. Whole-brain MRI measures were mapped onto twelve equivolumetric layers, and layer-specific sharp myeloarchitectonic transitions were identified using spatial gradients resulting in a putative 182 area/subarea partition of the macaque cerebral cortex. The myelin maps exhibit unexpectedly high homology with humans suggesting cortical myelin shares the same developmental program across the species. Comparison with histological Gallyas myelin stains explains over 80% of the variance in the laminar T1w/T2w-FLAIR profiles, substantiating the validity of the method. Altogether, our approach provides a novel, noninvasive means for precision mapping layer myeloarchitecture in the primate cerebral cortex, advancing the pioneering work of classical neuroanatomists.

Tutorial: a guide to diffusion MRI and structural connectomics

Diffusion magnetic resonance imaging (dMRI) is a versatile imaging technique that has gained popularity thanks to its sensitive ability to measure displacement of water molecules within a living tissue on a micrometer scale. Although dMRI has been around since the early 1990s, its applications are constantly evolving, primarily regarding the inference of structural connectomics from nerve fiber trajectories. However, these applications require expertise in image processing and statistics, and it can be difficult for a newcomer to choose an appropriate pipeline to fit their research needs, not least because dMRI is such a flexible methodology that dozens of acquisition and analysis pipelines have been developed over the years. This introductory guide is designed for graduate students and researchers in the neuroscience community who are interested in integrating this new methodology regardless of their background in neuroimaging and computational tools. The guide provides a brief overview of the basic dMRI methodologies but focuses on its applications in neuroplasticity and connectomics. The guide starts with dMRI experimental designs and a complete step-by-step pipeline for structural connectomics. The following section covers the basics of dMRI, including parameters and clinical applications (apparent diffusion coefficient, mean diffusivity, fractional anisotropy and microscopic fractional anisotropy), as well as different approaches and models. The final section focuses on structural connectomics, covering subjects from fiber tracking (techniques, evaluation and limitations) to structural networks (constructing, analyzing and visualizing a network).

β-Alanine supplementation improves fractional anisotropy scores in the hippocampus and amygdala in 60-80-year-old men and women

Recently, β-alanine (BA) supplementation was shown to improve cognitive function in older adults with decreased cognitive function. Mechanisms supporting these improvements have not been well defined. This study examined the effects of 10-weeks of BA supplementation on changes in circulating brain inflammatory markers, brain derived neurotrophic factor (BDNF), and brain morphology. Twenty participants were initially randomized into BA (2.4 g·d-1) or placebo (PL) groups. At each testing session, participants provided a resting blood sample and completed the Montreal cognitive assessment (MoCA) test and magnetic resonance imaging, which included diffusion tensor imaging to assess brain tissue integrity. Only participants that scored at or below normal for the MoCA assessment were analyzed (6 BA and 4 PL). The Mann-Whitney U test was used to examine Δ (POST-PRE) differences between the groups. No differences in Δ scores were noted in any blood marker (BDNF, CRP, TNF-α and GFAP). Changes in fractional anisotropy scores were significantly greater for BA than PL in the right hippocampus (p = 0.033) and the left amygdala (p = 0.05). No other differences were noted. The results provide a potential mechanism of how BA supplementation may improve cognitive function as reflected by improved tissue integrity within the hippocampus and amygdala.

Advanced Brain Imaging in Central Nervous System Demyelinating Diseases

In recent decades, advances in neuroimaging have profoundly transformed our comprehension of central nervous system demyelinating diseases. Remarkable technological progress has enabled the integration of cutting-edge acquisition and postprocessing techniques, proving instrumental in characterizing subtle focal changes, diffuse microstructural alterations, and macroscopic pathologic processes. This review delves into state-of-the-art modalities applied to multiple sclerosis, neuromyelitis optica spectrum disorders, and myelin oligodendrocyte glycoprotein antibody-associated disease. Furthermore, it explores how this dynamic landscape holds significant promise for the development of effective and personalized clinical management strategies, encompassing support for differential diagnosis, prognosis, monitoring treatment response, and patient stratification.

Circumscribing Laser Cuts Attenuate Seizure Propagation in a Mouse Model of Focal Epilepsy

In partial onset epilepsy, seizures arise focally in the brain and often propagate. Patients frequently become refractory to medical management, leaving neurosurgery, which can cause neurologic deficits, as a primary treatment. In the cortex, focal seizures spread through horizontal connections in layers II/III, suggesting that severing these connections can block seizures while preserving function. Focal neocortical epilepsy is induced in mice, sub-surface cuts are created surrounding the seizure focus using tightly-focused femtosecond laser pulses, and electrophysiological recordings are acquired at multiple locations for 3-12 months. Cuts reduced seizure frequency in most animals by 87%, and only 5% of remaining seizures propagated to the distant electrodes, compared to 80% in control animals. These cuts produced a modest decrease in cortical blood flow that recovered and left a ≈20-µm wide scar with minimal collateral damage. When placed over the motor cortex, cuts do not cause notable deficits in a skilled reaching task, suggesting they hold promise as a novel neurosurgical approach for intractable focal cortical epilepsy.

N-butyl-β-carboline-3-carboxylate (β-CCB) systemic administration promotes remyelination in the cuprizone demyelinating model in mice

Demyelination is generated in several nervous system illnesses. Developing strategies for effective clinical treatments requires the discovery of promyelinating drugs. Increased GABAergic signaling through γ-aminobutyric acid type A receptor (GABAAR) activation in oligodendrocytes has been proposed as a promyelinating condition. GABAAR expressed in oligodendroglia is strongly potentiated by n-butyl-β-carboline-3-carboxylate (β-CCB) compared to that in neurons. Here, mice were subjected to 0.3% cuprizone (CPZ) added in the food to induce central nervous system demyelination, a well-known model for multiple sclerosis. Then β-CCB (1 mg/Kg) was systemically administered to analyze the remyelination status in white and gray matter areas. Myelin content was evaluated using Black-Gold II (BGII) staining, immunofluorescence (IF), and magnetic resonance imaging (MRI). Evidence indicates that β-CCB treatment of CPZ-demyelinated animals promoted remyelination in several white matter structures, such as the fimbria, corpus callosum, internal capsule, and cerebellar peduncles. Moreover, using IF, it was observed that CPZ intake induced an increase in NG2+ and a decrease in CC1+ cell populations, alterations that were importantly retrieved by β-CCB treatment. Thus, the promyelinating character of β-CCB was confirmed in a generalized demyelination model, strengthening the idea that it has clinical potential as a therapeutic drug.

Associations between maternal pre-pregnancy BMI and infant striatal mean diffusivity

BACKGROUND

It is well-established that parental obesity is a strong risk factor for offspring obesity. Further, a converging body of evidence now suggests that maternal weight profiles may affect the developing offspring's brain in a manner that confers future obesity risk. Here, we investigated how pre-pregnancy maternal weight status influences the reward-related striatal areas of the offspring's brain during in utero development.

METHODS

We used diffusion tensor imaging to quantify the microstructure of the striatal brain regions of interest in neonates (N = 116 [66 males, 50 females], mean gestational weeks at birth [39.88], SD = 1.14; at scan [43.56], SD = 1.05). Linear regression was used to test the associations between maternal pre-pregnancy body mass index (BMI) and infant striatal mean diffusivity.

RESULTS

High maternal pre-pregnancy BMI was associated with higher mean MD values in the infant's left caudate nucleus. Results remained unchanged after the adjustment for covariates.

CONCLUSIONS

In utero exposure to maternal adiposity might have a growth-impairing impact on the mean diffusivity of the infant's left caudate nucleus. Considering the involvement of the caudate nucleus in regulating eating behavior and food-related reward processing later in life, this finding calls for further investigations to define the prognostic relevance of early-life caudate nucleus development and weight trajectories of the offspring.

Clustering the cortical laminae: in vivo parcellation

The laminar microstructure of the cerebral cortex has distinct anatomical characteristics of the development, function, connectivity, and even various pathologies of the brain. In recent years, multiple neuroimaging studies have utilized magnetic resonance imaging (MRI) relaxometry to visualize and explore this intricate microstructure, successfully delineating the cortical laminar components. Despite this progress, T1 is still primarily considered a direct measure of myeloarchitecture (myelin content), rather than a probe of tissue cytoarchitecture (cellular composition). This study aims to offer a robust, whole-brain validation of T1 imaging as a practical and effective tool for exploring the laminar composition of the cortex. To do so, we cluster complex microstructural cortical datasets of both human (N = 30) and macaque (N = 1) brains using an adaptation of an algorithm for clustering cell omics profiles. The resulting cluster patterns are then compared to established atlases of cytoarchitectonic features, exhibiting significant correspondence in both species. Lastly, we demonstrate the expanded applicability of T1 imaging by exploring some of the cytoarchitectonic features behind various unique skillsets, such as musicality and athleticism.

Microstructural dynamics of motor learning and sleep-dependent consolidation: A diffusion imaging study

Memory consolidation can benefit from post-learning sleep, eventually leading to long-term microstructural brain modifications to accommodate new memory representations. Non-invasive diffusion-weighted magnetic resonance imaging (DWI) allows the observation of (micro)structural brain remodeling after time-limited motor learning. Here, we combine conventional diffusion tensor imaging (DTI) and neurite orientation dispersion and density imaging (NODDI) that allows modeling dendritic and axonal complexity in gray matter to investigate with improved specificity the microstructural brain mechanisms underlying time- and sleep-dependent motor memory consolidation dynamics. Sixty-one young healthy adults underwent four DWI sessions, two sequential motor trainings, and a night of total sleep deprivation or regular sleep distributed over five days. We observed rapid-motor-learning-related remodeling in occipitoparietal, temporal, and motor-related subcortical regions, reflecting temporary dynamics in learning-related neuronal brain plasticity processes. Sleep-related consolidation seems not to exert a detectable impact on diffusion parameters, at least on the timescale of a few days.

Changes in Cortical Microstructure of the Human Brain Resulting from Long-Term Motor Learning

The mechanisms subserving motor skill acquisition and learning in the intact human brain are not fully understood. Previous studies in animals have demonstrated a causal relationship between motor learning and structural rearrangements of synaptic connections, raising the question of whether neurite-specific changes are also observable in humans. Here, we use advanced diffusion magnetic resonance imaging (MRI), sensitive to dendritic and axonal processes, to investigate neuroplasticity in response to long-term motor learning. We recruited healthy male and female human participants (age range 19-29) who learned a challenging dynamic balancing task (DBT) over four consecutive weeks. Diffusion MRI signals were fitted using Neurite Orientation Dispersion and Density Imaging (NODDI), a theory-driven biophysical model of diffusion, yielding measures of tissue volume, neurite density and the organizational complexity of neurites. While NODDI indices were unchanged and reliable during the control period, neurite orientation dispersion increased significantly during the learning period mainly in primary sensorimotor, prefrontal, premotor, supplementary, and cingulate motor areas. Importantly, reorganization of cortical microstructure during the learning phase predicted concurrent behavioral changes, whereas there was no relationship between microstructural changes during the control phase and learning. Changes in neurite complexity were independent of alterations in tissue density, cortical thickness, and intracortical myelin. Our results are in line with the notion that structural modulation of neurites is a key mechanism supporting complex motor learning in humans.SIGNIFICANCE STATEMENT The structural correlates of motor learning in the human brain are not fully understood. Results from animal studies suggest that synaptic remodeling (e.g., reorganization of dendritic spines) in sensorimotor-related brain areas is a crucial mechanism for the formation of motor memory. Using state-of-the-art diffusion magnetic resonance imaging (MRI), we found a behaviorally relevant increase in the organizational complexity of neocortical microstructure, mainly in primary sensorimotor, prefrontal, premotor, supplementary, and cingulate motor regions, following training of a challenging dynamic balancing task (DBT). Follow-up analyses suggested structural modulation of synapses as a plausible mechanism driving this increase, while colocalized changes in cortical thickness, tissue density, and intracortical myelin could not be detected. These results advance our knowledge about the neurobiological basis of motor learning in humans.

Mapping the individual human cortex using multidimensional MRI and unsupervised learning

Human evolution has seen the development of higher-order cognitive and social capabilities in conjunction with the unique laminar cytoarchitecture of the human cortex. Moreover, early-life cortical maldevelopment has been associated with various neurodevelopmental diseases. Despite these connections, there is currently no noninvasive technique available for imaging the detailed cortical laminar structure. This study aims to address this scientific and clinical gap by introducing an approach for imaging human cortical lamina. This method combines diffusion-relaxation multidimensional MRI with a tailored unsupervised machine learning approach that introduces enhanced microstructural sensitivity. This new imaging method simultaneously encodes the microstructure, the local chemical composition and importantly their correlation within complex and heterogenous tissue. To validate our approach, we compared the intra-cortical layers obtained using our ex vivo MRI-based method with those derived from Nissl staining of postmortem human brain specimens. The integration of unsupervised learning with diffusion-relaxation correlation MRI generated maps that demonstrate sensitivity to areal differences in cytoarchitectonic features observed in histology. Significantly, our observations revealed layer-specific diffusion-relaxation signatures, showing reductions in both relaxation times and diffusivities at the deeper cortical levels. These findings suggest a radial decrease in myelin content and changes in cell size and anisotropy, reflecting variations in both cytoarchitecture and myeloarchitecture. Additionally, we demonstrated that 1D relaxation and high-order diffusion MRI scalar indices, even when aggregated and used jointly in a multimodal fashion, cannot disentangle the cortical layers. Looking ahead, our technique holds the potential to open new avenues of research in human neurodevelopment and the vast array of disorders caused by disruptions in neurodevelopment.

Acute cycling exercise and hippocampal subfield function and microstructure in healthy older adults

Aging is associated with deterioration in dentate gyrus (DG) and CA3, both crucial hippocampal subfields for age susceptible memory processes such as mnemonic discrimination (MD). Meanwhile, a single aerobic exercise session alters DG/CA3 function and neural activity in both rats and younger adults and can elicit short-term microstructural alterations in the hippocampus of older adults. However, our understanding of the effects of acute exercise on hippocampal subfield integrity via function and microstructure in older adults is limited. Thus, a within subject-design was employed to determine if 20-min of moderate to vigorous aerobic exercise alters bilateral hippocampal subfield function and microstructure using high-resolution functional magnetic resonance imaging (fMRI) during an MD task (n = 35) and high angular resolution multi-shell diffusion imaging (n = 31), in healthy older adults, compared to seated rest. Following the exercise condition, participants exhibited poorer MD performance, particularly when their perception of effort was higher. Exercise was also related to lower MD-related activity within the DG/CA3 but not CA1 subfield. Finally, after controlling for whole brain gray matter diffusion, exercise was associated with lower neurite density index (NDI) within the DG/CA3. However, exercise-related differences in DG/CA3 activity and NDI were not associated with differences in MD performance. Our results suggest moderate to vigorous aerobic exercise may temporarily inhibit MD performance, and suppress DG/CA3 MD-related activity and NDI, potentially through neuroinflammatory/glial processes. However, additional studies are needed to confirm whether these short-term changes in behavior and hippocampal subfield neurophysiology are beneficial and how they might relate to long-term exercise habits.

Patient-specific models link neurotransmitter receptor mechanisms with motor and visuospatial axes of Parkinson's disease

Parkinson's disease involves multiple neurotransmitter systems beyond the classical dopaminergic circuit, but their influence on structural and functional alterations is not well understood. Here, we use patient-specific causal brain modeling to identify latent neurotransmitter receptor-mediated mechanisms contributing to Parkinson's disease progression. Combining the spatial distribution of 15 receptors from post-mortem autoradiography with 6 neuroimaging-derived pathological factors, we detect a diverse set of receptors influencing gray matter atrophy, functional activity dysregulation, microstructural degeneration, and dendrite and dopaminergic transporter loss. Inter-individual variability in receptor mechanisms correlates with symptom severity along two distinct axes, representing motor and psychomotor symptoms with large GABAergic and glutamatergic contributions, and cholinergically-dominant visuospatial, psychiatric and memory dysfunction. Our work demonstrates that receptor architecture helps explain multi-factorial brain re-organization, and suggests that distinct, co-existing receptor-mediated processes underlie Parkinson's disease.

Structural and functional changes of the cerebellum in temporal lobe epilepsy

Aims

This study aimed to comprehensively explore the cerebellar structural and functional changes in temporal lobe epilepsy (TLE) and its association with clinical information.

Methods

The SUIT toolbox was utilized to perform cerebellar volume and diffusion analysis. In addition, we extracted the average diffusion values of cerebellar peduncle tracts to investigate microstructure alterations. Seed-based whole-brain analysis was used to investigate cerebellar-cerebral functional connectivity (FC). Subgroup analyses were performed to identify the cerebellar participation in TLE with/without hippocampal sclerosis (HS)/focal-to-bilateral tonic-clonic seizure (FBTCS) and TLE with different lateralization.

Results

TLE showed widespread gray matter atrophy in bilateral crusII, VIIb, VIIIb, left crusI, and left VIIIa. Both voxel and tract analysis observed diffusion abnormalities in cerebellar afferent peduncles. Reduced FC between the right crus II and the left parahippocampal cortex was found in TLE. Additionally, TLE showed increased FCs between left lobules VI-VIII and cortical nodes of the dorsal attention and visual networks. Across all patients, decreased FC was associated with poorer cognitive function, while increased FCs appeared to reflect compensatory effects. The cerebellar structural changes were mainly observed in HS and FBTCS subgroups and were regardless of seizure lateralization, while cerebellar-cerebral FC alterations were similar in all subgroups.

Conclusion

TLE exhibited microstructural changes in the cerebellum, mainly related to HS and FBTCS. In addition, altered cerebellar-cerebral functional connectivity is associated with common cognitive alterations in TLE.

High-resolution multi-shot diffusion imaging of structural networks in healthy neurocognitive aging

Healthy neurocognitive aging has been associated with the microstructural degradation of white matter pathways that connect distributed gray matter regions, assessed by diffusion-weighted imaging (DWI). However, the relatively low spatial resolution of standard DWI has limited the examination of age-related differences in the properties of smaller, tightly curved white matter fibers, as well as the relatively more complex microstructure of gray matter. Here, we capitalize on high-resolution multi-shot DWI, which allows spatial resolutions < 1 mm3 to be achieved on clinical 3T MRI scanners. We assessed whether traditional diffusion tensor-based measures of gray matter microstructure and graph theoretical measures of white matter structural connectivity assessed by standard (1.5 mm3 voxels, 3.375 μl volume) and high-resolution (1 mm3 voxels, 1μl volume) DWI were differentially related to age and cognitive performance in 61 healthy adults 18-78 years of age. Cognitive performance was assessed using an extensive battery comprising 12 separate tests of fluid (speed-dependent) cognition. Results indicated that the high-resolution data had larger correlations between age and gray matter mean diffusivity, but smaller correlations between age and structural connectivity. Moreover, parallel mediation models including both standard and high-resolution measures revealed that only the high-resolution measures mediated age-related differences in fluid cognition. These results lay the groundwork for future studies planning to apply high-resolution DWI methodology to further assess the mechanisms of both healthy aging and cognitive impairment.

High-fidelity mesoscale in-vivo diffusion MRI through gSlider-BUDA and circular EPI with S-LORAKS reconstruction

PURPOSE

To develop a high-fidelity diffusion MRI acquisition and reconstruction framework with reduced echo-train-length for less T2* image blurring compared to typical highly accelerated echo-planar imaging (EPI) acquisitions at sub-millimeter isotropic resolution.

METHODS

We first proposed a circular-EPI trajectory with partial Fourier sampling on both the readout and phase-encoding directions to minimize the echo-train-length and echo time. We then utilized this trajectory in an interleaved two-shot EPI acquisition with reversed phase-encoding polarity, to aid in the correction of off-resonance-induced image distortions and provide complementary k-space coverage in the missing partial Fourier regions. Using model-based reconstruction with structured low-rank constraint and smooth phase prior, we corrected the shot-to-shot phase variations across the two shots and recover the missing k-space data. Finally, we combined the proposed acquisition/reconstruction framework with an SNR-efficient RF-encoded simultaneous multi-slab technique, termed gSlider, to achieve high-fidelity 720 µm and 500 µm isotropic resolution in-vivo diffusion MRI.

RESULTS

Both simulation and in-vivo results demonstrate the effectiveness of the proposed acquisition and reconstruction framework to provide distortion-corrected diffusion imaging at the mesoscale with markedly reduced T2*-blurring. The in-vivo results of 720 µm and 500 µm datasets show high-fidelity diffusion images with reduced image blurring and echo time using the proposed approaches.

CONCLUSIONS

The proposed method provides high-quality distortion-corrected diffusion-weighted images with ∼40% reduction in the echo-train-length and T2* blurring at 500µm-isotropic-resolution compared to standard multi-shot EPI.

Expanding connectomics to the laminar level: A perspective

Despite great progress in uncovering the complex connectivity patterns of the human brain over the last two decades, the field of connectomics still experiences a bias in its viewpoint of the cerebral cortex. Due to a lack of information regarding exact end points of fiber tracts inside cortical gray matter, the cortex is commonly reduced to a single homogenous unit. Concurrently, substantial developments have been made over the past decade in the use of relaxometry and particularly inversion recovery imaging for exploring the laminar microstructure of cortical gray matter. In recent years, these developments have culminated in an automated framework for cortical laminar composition analysis and visualization, followed by studies of cortical dyslamination in epilepsy patients and age-related differences in laminar composition in healthy subjects. This perspective summarizes the developments and remaining challenges of multi-T1 weighted imaging of cortical laminar substructure, the current limitations in structural connectomics, and the recent progress in integrating these fields into a new model-based subfield termed 'laminar connectomics'. In the coming years, we predict an increased use of similar generalizable, data-driven models in connectomics with the purpose of integrating multimodal MRI datasets and providing a more nuanced and detailed characterization of brain connectivity.

Music impacts brain cortical microstructural maturation in very preterm infants: A longitudinal diffusion MR imaging study

Preterm birth disrupts important neurodevelopmental processes occurring from mid-fetal to term-age. Musicotherapy, by enriching infants' sensory input, might enhance brain maturation during this critical period of activity-dependent plasticity. To study the impact of music on preterm infants' brain structural changes, we recruited 54 very preterm infants randomized to receive or not a daily music intervention, that have undergone a longitudinal multi-shell diffusion MRI acquisition, before the intervention (at 33 weeks' gestational age) and after it (at term-equivalent-age). Using whole-brain fixel-based (FBA) and NODDI analysis (n = 40), we showed a longitudinal increase of fiber cross-section (FC) and fiber density (FD) in all major cerebral white matter fibers. Regarding cortical grey matter, FD decreased while FC and orientation dispersion index (ODI) increased, reflecting intracortical multidirectional complexification and intracortical myelination. The music intervention resulted in a significantly higher longitudinal increase of FC and ODI in cortical paralimbic regions, namely the insulo-orbito-temporopolar complex, precuneus/posterior cingulate gyrus, as well as the auditory association cortex. Our results support a longitudinal early brain macro and microstructural maturation of white and cortical grey matter in preterm infants. The music intervention led to an increased intracortical complexity in regions important for socio-emotional development, known to be impaired in preterm infants.

Integrating multimodality magnetic resonance imaging to the Allen Mouse Brain Common Coordinate Framework

High-resolution magnetic resonance imaging (MRI) affords unique image contrasts to nondestructively probe the tissue microstructure; validation of MRI findings with conventional histology is essential to better understand the MRI contrasts. However, the dramatic difference in the spatial resolution and image contrast of these two techniques impedes accurate comparison between MRI metrics and traditional histology. To better validate various MRI metrics, we acquired whole mouse brain multigradient recalled-echo and multishell diffusion MRI datasets at 25-μm isotropic resolution. The recently developed Allen Mouse Brain Common Coordinate Framework (CCFv3) provides opportunities to integrate multimodal and multiscale datasets of the whole mouse brain in a common three-dimensional (3D) space. The T2*, quantitative susceptibility mapping, diffusion tensor imaging, and neurite orientation dispersion and density imaging parameters were compared with both serial two-photon tomography images and 3D Nissl staining images in the CCFv3 at the same spatial resolution. The correlation between MRI and Nissl staining strongly depends on different metrics and different regions of the brain. Integrating different imaging modalities to the same space may substantially improve our understanding of the complexity of the brain at different scales.

Acquisition and Analysis of Excised Neocortex from Pediatric Patients with Focal Cortical Dysplasia Using Mesoscale Diffusion MRI

Non-invasive classification of focal cortical dysplasia (FCD) subtypes remains challenging from a radiology perspective. Quantitative imaging biomarkers (QIBs) have the potential to distinguish subtypes that lack pathognomonic features and might help in defining the extent of abnormal connectivity associated with each FCD subtype. A key motivation of diagnostic imaging is to improve the localization of a "lesion" that can guide the surgical resection of affected tissue, which is thought to cause seizures. Conversely, surgical resections to eliminate or reduce seizures provided unique opportunities to develop magnetic resonance imaging (MRI)-based QIBs by affording long scan times to evaluate multiple contrast mechanisms at the mesoscale (0.5 mm isotropic voxel dimensions). Using ex vivo hybrid diffusion tensor imaging on a 9.4 T MRI scanner, the grey to white matter ratio of scalar indices was lower in the resected middle temporal gyrus (MTG) of two neuropathologically confirmed cases of FCD compared to non-diseased control postmortem fixed temporal lobes. In contrast, fractional anisotropy was increased within FCD and also adjacent white matter tracts. Connectivity (streamlines/mm3) in the MTG was higher in FCD, suggesting that an altered connectivity at the lesion locus can potentially provide a tangible QIB to distinguish and characterize FCD abnormalities. However, as illustrated here, a major challenge for a robust tractographical comparison lies in the considerable differences in the ex vivo processing of bioptic and postmortem samples. Mesoscale diffusion MRI has the potential to better define and characterize epileptic tissues obtained from surgical resection to advance our understanding of disease etiology and treatment.

Column-based cortical depth analysis of the diffusion anisotropy and radiality in submillimeter whole-brain diffusion tensor imaging of the human cortical gray matter in vivo

High-resolution diffusion tensor imaging (DTI) can noninvasively probe the microstructure of cortical gray matter in vivo. In this study, 0.9-mm isotropic whole-brain DTI data were acquired in healthy subjects with an efficient multi-band multi-shot echo-planar imaging sequence. A column-based analysis that samples the fractional anisotropy (FA) and radiality index (RI) along radially oriented cortical columns was then performed to quantitatively analyze the FA and RI dependence on the cortical depth, cortical region, cortical curvature, and cortical thickness across the whole brain, which has not been simultaneously and systematically investigated in previous studies. The results showed characteristic FA and RI vs. cortical depth profiles, with an FA local maximum and minimum (or two inflection points) and a single RI maximum at intermediate cortical depths in most cortical regions, except for the postcentral gyrus where no FA peaks and a lower RI were observed. These results were consistent between repeated scans from the same subjects and across different subjects. They were also dependent on the cortical curvature and cortical thickness in that the characteristic FA and RI peaks were more pronounced i) at the banks than at the crown of gyri or at the fundus of sulci and ii) as the cortical thickness increases. This methodology can help characterize variations in microstructure along the cortical depth and across the whole brain in vivo, potentially providing quantitative biomarkers for neurological disorders.

High-fidelity, high-spatial-resolution diffusion magnetic resonance imaging of ex vivo whole human brain at ultra-high gradient strength with structured low-rank echo-planar imaging ghost correction

Diffusion magnetic resonance imaging (dMRI) of whole ex vivo human brain specimens enables three-dimensional (3D) mapping of structural connectivity at the mesoscopic scale, providing detailed evaluation of fiber architecture and tissue microstructure at a spatial resolution that is difficult to access in vivo. To account for the short T2 and low diffusivity of fixed tissue, ex vivo dMRI is often acquired using strong diffusion-sensitizing gradients and multishot/segmented 3D echo-planar imaging (EPI) sequences to achieve high spatial resolution. However, the combination of strong diffusion-sensitizing gradients and multishot/segmented EPI readout can result in pronounced ghosting artifacts incurred by nonlinear spatiotemporal variations in the magnetic field produced by eddy currents. Such ghosting artifacts cannot be corrected with conventional correction solutions and pose a significant roadblock to leveraging human MRI scanners with ultrahigh gradients for ex vivo whole-brain dMRI. Here, we show that ghosting-correction approaches that correct for either polarity-related ghosting or shot-to-shot variations in a separate manner are suboptimal for 3D multishot diffusion-weighted EPI experiments in fixed human brain specimens using strong diffusion-sensitizing gradients on the 3-T Connectom MRI scanner, resulting in orientationally biased dMRI estimates. We apply a recently developed advanced k-space reconstruction method based on structured low-rank matrix (SLM) modeling that handles both polarity-related ghosting and shot-to-shot variation simultaneously, to mitigate artifacts in high-angular resolution multishot dMRI data acquired in several fixed human brain specimens at 0.7-0.8-mm isotropic spatial resolution using b-values up to 10,000 s/mm2 and gradient strengths up to 280 mT/m. We demonstrate the improved mapping of diffusion tensor imaging and fiber orientation distribution functions in key neuroanatomical areas distributed across the whole brain using SLM-based EPI ghost correction compared with alternative techniques.

Reorganization of thalamocortical connections in congenitally blind humans

Cross-modal plasticity in blind individuals has been reported over the past decades showing that nonvisual information is carried and processed by "visual" brain structures. However, despite multiple efforts, the structural underpinnings of cross-modal plasticity in congenitally blind individuals remain unclear. We mapped thalamocortical connectivity and assessed the integrity of white matter of 10 congenitally blind individuals and 10 sighted controls. We hypothesized an aberrant thalamocortical pattern of connectivity taking place in the absence of visual stimuli from birth as a potential mechanism of cross-modal plasticity. In addition to the impaired microstructure of visual white matter bundles, we observed structural connectivity changes between the thalamus and occipital and temporal cortices. Specifically, the thalamic territory dedicated to connections with the occipital cortex was smaller and displayed weaker connectivity in congenitally blind individuals, whereas those connecting with the temporal cortex showed greater volume and increased connectivity. The abnormal pattern of thalamocortical connectivity included the lateral and medial geniculate nuclei and the pulvinar nucleus. For the first time in humans, a remapping of structural thalamocortical connections involving both unimodal and multimodal thalamic nuclei has been demonstrated, shedding light on the possible mechanisms of cross-modal plasticity in humans. The present findings may help understand the functional adaptations commonly observed in congenitally blind individuals.

Microstructural Neurodegeneration of the Entorhinal-Hippocampus Pathway along the Alzheimer's Disease Continuum

BACKGROUND

Conventional neuroimaging biomarkers for the neurodegeneration of Alzheimer's disease (AD) are not sensitive enough to detect neurodegenerative alterations during the preclinical stage of AD individuals.

OBJECTIVE

We examined whether neurodegeneration of the entorhinal-hippocampal pathway could be detected along the AD continuum using ultra-high-field diffusion tensor imaging and tractography for ex vivo brain tissues.

METHODS

Postmortem brain specimens from a cognitively unimpaired individual without AD pathological changes (non-AD), a cognitively unimpaired individual with AD pathological changes (preclinical AD), and a demented individual with AD pathological changes (AD dementia) were scanned with an 11.7T diffusion magnetic resonance imaging. Fractional anisotropy (FA) values of the entorhinal layer II and number of perforant path fibers counted by tractography were compared among the AD continuum. Following the imaging analyses, the status of myelinated fibers and neuronal cells were verified by subsequent serial histological examinations.

RESULTS

At 250μm (zipped to 125μm) isotropic resolution, the entorhinal layer II islands and the perforant path fibers could be identified in non-AD and preclinical AD, but not in AD dementia, followed by histological verification. The FA value of the entorhinal layer II was the highest among the entorhinal laminae in non-AD and preclinical AD, whereas the FA values in the entorhinal laminae were homogeneously low in AD dementia. The FA values and number of perforant path fibers decreased along the AD continuum (non-AD>preclinical AD > AD dementia).

CONCLUSION

We successfully detected neurodegenerative alterations of the entorhinal-hippocampal pathway at the preclinical stage of the AD continuum.

COnstrained Reference frame diffusion TEnsor Correlation Spectroscopic (CORTECS) MRI: A practical framework for high-resolution diffusion tensor distribution imaging

High-resolution imaging studies have consistently shown that in cortical tissue water diffuses preferentially along radial and tangential orientations with respect to the cortical surface, in agreement with histology. These dominant orientations do not change significantly even if the relative contributions from microscopic water pools to the net voxel signal vary across experiments that use different diffusion times, b-values, TEs, and TRs. With this in mind, we propose a practical new framework for imaging non-parametric diffusion tensor distributions (DTDs) by constraining the microscopic diffusion tensors of the DTD to be diagonalized using the same orthonormal reference frame of the mesoscopic voxel. In each voxel, the constrained DTD (cDTD) is completely determined by the correlation spectrum of the microscopic principal diffusivities associated with the axes of the voxel reference frame. Consequently, all cDTDs are inherently limited to the domain of positive definite tensors and can be reconstructed efficiently using Inverse Laplace Transform methods. Moreover, the cDTD reconstruction can be performed using only data acquired efficiently with single diffusion encoding, although it also supports datasets with multiple diffusion encoding. In tissues with a well-defined architecture, such as the cortex, we can further constrain the cDTD to contain only cylindrically symmetric diffusion tensors and measure the 2D correlation spectra of principal diffusivities along the radial and tangential orientation with respect to the cortical surface. To demonstrate this framework, we perform numerical simulations and analyze high-resolution dMRI data from a fixed macaque monkey brain. We estimate 2D cDTDs in the cortex and derive, in each voxel, the marginal distributions of the microscopic principal diffusivities, the corresponding distributions of the microscopic fractional anisotropies and mean diffusivities along with their 2D correlation spectra to quantify the cDTD shape-size characteristics. Signal components corresponding to specific bands in these cDTD-derived spectra show high specificity to cortical laminar structures observed with histology. Our framework drastically simplifies the measurement of non-parametric DTDs in high-resolution datasets with mesoscopic voxel sizes much smaller than the radius of curvature of the underlying anatomy, e.g., cortical surface, and can be applied retrospectively to analyze existing diffusion MRI data from fixed cortical tissues.

High-resolution cortical MAP-MRI reveals areal borders and laminar substructures observed with histological staining

The variations in cellular composition and tissue architecture measured with histology provide the biological basis for partitioning the brain into distinct cytoarchitectonic areas and for characterizing neuropathological tissue alterations. Clearly, there is an urgent need to develop whole-brain neuroradiological methods that can assess cortical cyto- and myeloarchitectonic features non-invasively. Mean apparent propagator (MAP) MRI is a clinically feasible diffusion MRI method that quantifies efficiently and comprehensively the net microscopic displacements of water molecules diffusing in tissues. We investigate the sensitivity of high-resolution MAP-MRI to detecting areal and laminar variations in cortical cytoarchitecture and compare our results with observations from corresponding histological sections in the entire brain of a rhesus macaque monkey. High-resolution images of MAP-derived parameters, in particular the propagator anisotropy (PA), non-gaussianity (NG), and the return-to-axis probability (RTAP) reveal cortical area-specific lamination patterns in good agreement with the corresponding histological stained sections. In a few regions, the MAP parameters provide superior contrast to the five histological stains used in this study, delineating more clearly boundaries and transition regions between cortical areas and laminar substructures. Throughout the cortex, various MAP parameters can be used to delineate transition regions between specific cortical areas observed with histology and to refine areal boundaries estimated using atlas registration-based cortical parcellation. Using surface-based analysis of MAP parameters we quantify the cortical depth dependence of diffusion propagators in multiple regions-of-interest in a consistent and rigorous manner that is largely independent of the cortical folding geometry. The ability to assess cortical cytoarchitectonic features efficiently and non-invasively, its clinical feasibility, and translatability make high-resolution MAP-MRI a promising 3D imaging tool for studying whole-brain cortical organization, characterizing abnormal cortical development, improving early diagnosis of neurodegenerative diseases, identifying targets for biopsies, and complementing neuropathological investigations.

Diffusion MRI anisotropy in the cerebral cortex is determined by unmyelinated tissue features

Diffusion magnetic resonance imaging (dMRI) is commonly used to assess the tissue and cellular substructure of the human brain. In the white matter, myelinated axons are the principal neural elements that shape dMRI through the restriction of water diffusion; however, in the gray matter the relative contributions of myelinated axons and other tissue features to dMRI are poorly understood. Here we investigate the determinants of diffusion in the cerebral cortex. Specifically, we ask whether myelinated axons significantly shape dMRI fractional anisotropy (dMRI-FA), a measure commonly used to characterize tissue properties in humans. We compared ultra-high resolution ex vivo dMRI data from the brain of a marmoset monkey with both myelin- and Nissl-stained histological sections obtained from the same brain after scanning. We found that the dMRI-FA did not match the spatial distribution of myelin in the gray matter. Instead dMRI-FA was more closely related to the anisotropy of stained tissue features, most prominently those revealed by Nissl staining and to a lesser extent those revealed by myelin staining. Our results suggest that unmyelinated neurites such as large caliber apical dendrites are the primary features shaping dMRI measures in the cerebral cortex.

Hippocampal dentate gyrus integrity revealed with ultrahigh resolution diffusion imaging predicts memory performance in older adults

Medial temporal lobe (MTL) atrophy is a core feature of age-related cognitive decline and Alzheimer's disease (AD). While regional volumes and thickness are often used as a proxy for neurodegeneration, they lack the sensitivity to serve as an accurate diagnostic test and indicate advanced neurodegeneration. Here, we used a submillimeter resolution diffusion weighted MRI sequence (ZOOMit) to quantify microstructural properties of hippocampal subfields in older adults (63-98 years old) using tensor derived measures: fractional anisotropy (FA) and mean diffusivity (MD). We demonstrate that the high-resolution sequence, and not a standard resolution sequence, identifies dissociable profiles for CA1, dentate gyrus (DG), and the collateral sulcus. Using ZOOMit, we show that advanced age is associated with increased MD of the CA1 and DG as well as decreased FA of the DG. Increased MD of the DG, reflecting decreased cellular density, mediated the relationship between age and word list recall. Further, increased MD in the DG, but not DG volume, was linked to worse spatial pattern separation. Our results demonstrate that ultrahigh-resolution diffusion imaging enables the detection of microstructural differences in hippocampal subfield integrity and will lead to novel insights into the mechanisms of age-related memory loss.

Inter- and intra-individual variation in brain structural-cognition relationships in aging

The sources of inter- and intra-individual variability in age-related cognitive decline remain poorly understood. We examined the association between 20-year trajectories of cognitive decline and multimodal brain structure and morphology in older age. We used the Whitehall II Study, an extensively characterised cohort with 3T brain magnetic resonance images acquired at older age (mean age = 69.52 ± 4.9) and 5 repeated cognitive performance assessments between mid-life (mean age = 53.2 ±4.9 years) and late-life (mean age = 67.7 ± 4.9). Using non-negative matrix factorization, we identified 10 brain components integrating cortical thickness, surface area, fractional anisotropy, and mean and radial diffusivities. We observed two latent variables describing distinct brain-cognition associations. The first describes variations in 5 structural components associated with low mid-life performance across multiple cognitive domains, decline in reasoning, but maintenance of fluency abilities. The second describes variations in 6 structural components associated with low mid-life performance in fluency and memory, but retention of multiple abilities. Expression of latent variables predicts future cognition 3.2 years later (mean age = 70.87 ± 4.9). This data-driven approach highlights brain-cognition relationships wherein individuals degrees of cognitive decline and maintenance across diverse cognitive functions are both positively and negatively associated with markers of cortical structure.

High B-value diffusion tensor imaging for early detection of hippocampal microstructural alteration in a mouse model of multiple sclerosis

Several studies have highlighted the value of diffusion tensor imaging (DTI) with strong diffusion weighting to reveal white matter microstructural lesions, but data in gray matter (GM) remains scarce. Herein, the effects of b-values combined with different numbers of diffusion-encoding directions (NDIRs) on DTI metrics to capture the normal hippocampal microstructure and its early alterations were investigated in a mouse model of multiple sclerosis (experimental autoimmune encephalomyelitis [EAE]). Two initial DTI datasets (B2700-43Dir acquired with b = 2700 s.mm⁻² and NDIR = 43; B1000-22Dir acquired with b = 1000 s.mm⁻² and NDIR = 22) were collected from 18 normal and 18 EAE mice at 4.7 T. Three additional datasets (B2700-22Dir, B2700-12Dir and B1000-12Dir) were extracted from the initial datasets. In healthy mice, we found a significant influence of b-values and NDIR on all DTI metrics. Confronting unsupervised hippocampal layers classification to the true anatomical classification highlighted the remarkable discrimination of the molecular layer with B2700-43Dir compared with the other datasets. Only DTI from the B2700 datasets captured the dendritic loss occurring in the molecular layer of EAE mice. Our findings stress the needs for both high b-values and sufficient NDIR to achieve a GM DTI with more biologically meaningful correlations, though DTI-metrics should be interpreted with caution in these settings.

Cerebral cortex layer segmentation using diffusion magnetic resonance imaging in vivo with applications to laminar connections and working memory analysis

Understanding the laminar brain structure is of great help in further developing our knowledge of the functions of the brain. However, since most layer segmentation methods are invasive, it is difficult to apply them to the human brain in vivo. To systematically explore the human brain's laminar structure noninvasively, the K-means clustering algorithm was used to automatically segment the left hemisphere into two layers, the superficial and deep layers, using a 7 Tesla (T) diffusion magnetic resonance imaging (dMRI)open dataset. The obtained layer thickness was then compared with the layer thickness of the BigBrain reference dataset, which segmented the neocortex into six layers based on the von Economo atlas. The results show a significant correlation not only between our automatically segmented superficial layer thickness and the thickness of layers 1-3 from the reference histological data, but also between our automatically segmented deep layer thickness and the thickness of layers 4-6 from the reference histological data. Second, we constructed the laminar connections between two pairs of unidirectional connected regions, which is consistent with prior research. Finally, we conducted the laminar analysis of the working memory, which was challenging to do in the past, and explained the conclusions of the functional analysis. Our work successfully demonstrates that it is possible to segment the human cortex noninvasively into layers using dMRI data and further explores the mechanisms of the human brain.

In vivo measurements of lamination patterns in the human cortex

The laminar composition of the cerebral cortex is tightly connected to the development and connectivity of the brain, as well as to function and pathology. Although most of the research on the cortical layers is done with the aid of ex vivo histology, there have been recent attempts to use magnetic resonance imaging (MRI) with potential in vivo applications. However, the high-resolution MRI technology and protocols required for such studies are neither common nor practical. In this article, we present a clinically feasible method for assessing the laminar properties of the human cortex using standard pulse sequence available on any common MRI scanner. Using a series of low-resolution inversion recovery (IR) MRI scans allows us to calculate multiple T1 relaxation time constants for each voxel. Based on the whole-brain T1 -distribution, we identify six different gray matter T1 populations and their variation across the cortex. Based on this, we show age-related differences in these population and demonstrate that this method is able to capture the difference in laminar composition across varying brain areas. We also provide comparison to ex vivo high-resolution MRI scans. We show that this method is feasible for the estimation of layer variability across large population cohorts, which can lead to research into the links between the cortical layers and function, behavior and pathologies that was heretofore unexplorable.

Using Advanced Diffusion-Weighted Imaging to Predict Cell Counts in Gray Matter: Potential and Pitfalls

Recent advances in diffusion imaging have given it the potential to non-invasively detect explicit neurobiological properties, beyond what was previously possible with conventional structural imaging. However, there is very little known about what cytoarchitectural properties these metrics, especially those derived from newer multi-shell models like Neurite Orientation Dispersion and Density Imaging (NODDI) correspond to. While these diffusion metrics do not promise any inherent cell type specificity, different brain cells have varying morphologies, which could influence the diffusion signal in distinct ways. This relationship is currently not well-characterized. Understanding the possible cytoarchitectural signatures of diffusion measures could allow them to estimate important neurobiological properties like cell counts, potentially resulting in a powerful clinical diagnostic tool. Here, using advanced diffusion imaging (NODDI) in the mouse brain, we demonstrate that different regions have unique relationships between cell counts and diffusion metrics. We take advantage of this exclusivity to introduce a framework to predict cell counts of different types of cells from the diffusion metrics alone, in a region-specific manner. We also outline the challenges of reliably developing such a model and discuss the precautions the field must take when trying to tie together medical imaging modalities and histology.

Higher-order multi-shell diffusion measures complement tensor metrics and volume in gray matter when predicting age and cognition

Recent advances in diffusion-weighted imaging have enabled us to probe the microstructure of even gray matter non-invasively. However, these advanced multi-shell protocols are often not included in large-scale studies as they significantly increase scan time. In this study, we investigated whether one set of multi-shell diffusion metrics commonly used in gray matter (as derived from Neurite Orientation Dispersion and Density Imaging, NODDI) provide enough additional information over typical tensor and volume metrics to justify the increased acquisition time, using the cognitive aging framework in the human hippocampus as a testbed. We first demonstrated that NODDI metrics are robust and reliable by replicating previous findings from our lab in a larger population of 79 younger (20.41 ± 1.89 years, 46 females) and 75 older (73.56 ± 6.26 years, 45 females) adults, showing that these metrics in the hippocampal subfields are sensitive to age and memory performance. We then asked how these subfield specific hippocampal NODDI metrics compared with standard tensor metrics and volume in predicting age and memory ability. We discovered that both NODDI and tensor measures separately predicted age and cognition in comparable capacities. However, integrating these modalities together considerably increased the predictive power of our logistic models, indicating that NODDI and tensor measures may be capturing independent microstructural information. We use these findings to encourage neuroimaging data collection consortiums to include a multi-shell diffusion sequence in their protocols since existing NODDI measures (and potential future multi-shell measures) may be able to capture microstructural variance that is missed by traditional approaches, even in studies exclusively examining gray matter.

Neurite dispersion and density mediates the relationship between cardiorespiratory fitness and cognition in healthy younger adults

Growing evidence suggests physical activity and cardiorespiratory fitness are associated with better cognition across the lifespan. However, the neurobiological underpinnings relating fitness and cognition remain unclear, particularly in healthy younger adults. Using a well-established and popular multi-compartment diffusion modeling approach, called Neurite Orientation and Dispersion and Density Imaging (NODDI), we investigated the relationship between physical fitness (measured via a 2-min walk test), cognition (fluid and crystallized), and gray and white matter microstructure, in a large sample (n = 816) of healthy younger adults (ages 22-35 years) from the human connectome project (HCP). Concurrent with previous literature, we found that fitness was positively associated with both fluid and crystallized cognition. Furthermore, we found that physical fitness was negatively associated with white matter orientation dispersion index (ODIWM) around the cerebellar peduncle and was negatively associated with widespread cortical and subcortical gray matter neurite density index (NDIGM). Lower ODIWM of the cerebral peduncle was associated with better fluid cognitive performance, while lower NDIGM was associated with better crystallized cognition. Finally, we found that while ODIWM partially mediated the relationship between fitness and fluid cognition, NDIGM partially mediated the relationship between fitness and crystallized cognition. This study is the first to explore the relationship between physical fitness and white and gray matter microstructure measures using NODDI. Our findings suggest that in addition to improved cognitive performance, higher physical fitness may be associated with lower white matter tract dispersion and lower neurite density in the cortical and subcortical gray matter of healthy younger adults.

Resolution and b value dependent Structural Connectome in ex vivo Mouse Brain

Diffusion magnetic resonance imaging has been widely used in both clinical and preclinical studies to characterize tissue microstructure and structural connectivity. The diffusion MRI protocol for the Human Connectome Project (HCP) has been developed and optimized to obtain high-quality, high-resolution diffusion MRI (dMRI) datasets. However, such efforts have not been fully explored in preclinical studies, especially for rodents. In this study, high quality dMRI datasets of mouse brains were acquired at 9.4T system from two vendors. In particular, we acquired a high-spatial resolution dMRI dataset (25 μm isotropic with 126 diffusion encoding directions), which we believe to be the highest spatial resolution yet obtained; and a high-angular resolution dMRI dataset (50 μm isotropic with 384 diffusion encoding directions), which we believe to be the highest angular resolution compared to the dMRI datasets at the microscopic resolution. We systematically investigated the effects of three important parameters that affect the final outcome of the connectome: b value (1000s/mm² to 8000 s/mm²), angular resolution (10 to 126), and spatial resolution (25 μm to 200 μm). The stability of tractography and connectome increase with the angular resolution, where more than 50 angles is necessary to achieve consistent results. The connectome and quantitative parameters derived from graph theory exhibit a linear relationship to the b value (R² > 0.99); a single-shell acquisition with b value of 3000 s/mm² shows comparable results to the multi-shell high angular resolution dataset. The dice coefficient decreases and both false positive rate and false negative rate gradually increase with coarser spatial resolution. Our study provides guidelines and foundations for exploration of tradeoffs among acquisition parameters for the structural connectome in ex vivo mouse brain.

Signature laminar distributions of pathology in frontotemporal lobar degeneration

Frontotemporal lobar degeneration (FTLD) with either tau (FTLD-tau) or TDP-43 (FTLD-TDP) inclusions are distinct proteinopathies that frequently cause similar frontotemporal dementia (FTD) clinical syndromes. FTD syndromes often display macroscopic signatures of neurodegeneration at the level of regions and networks, but it is unclear if subregional laminar pathology display patterns unique to proteinopathy or clinical syndrome. We hypothesized that FTLD-tau and FTLD-TDP accumulate pathology in relatively distinct cortical layers independent of clinical syndrome, with greater involvement of lower layers in FTLD-tau. The current study examined 170 patients with either FTLD-tau (n = 73) or FTLD-TDP (n = 97) spanning dementia and motor phenotypes in the FTD spectrum. We digitally measured the percent area occupied by tau and TDP-43 pathology in upper layers (I-III), lower layers (IV-VI), and juxtacortical white matter (WM) from isocortical regions in both hemispheres where available. Linear mixed-effects models compared ratios of upper to lower layer pathology between FTLD groups and investigated relationships with regions, WM pathology, and global cognitive impairment while adjusting for demographics. We found lower ratios of layer pathology in FTLD-tau and higher ratios of layer pathology in FTLD-TDP, reflecting lower layer-predominant tau pathology and upper layer-predominant TDP-43 pathology, respectively (p < 0.001). FTLD-tau displayed lower ratios of layer pathology related to greater WM tau pathology (p = 0.002) and to earlier involved/severe pathology regions (p = 0.007). In contrast, FTLD-TDP displayed higher ratios of layer pathology not related to either WM pathology or regional severity. Greater cognitive impairment was associated with higher ratios of layer pathology in FTLD-tau (p = 0.018), but was not related to ratios of layer pathology in FTLD-TDP. Lower layer-predominant tau pathology and upper layer-predominant TDP-43 pathology are proteinopathy-specific, regardless of clinical syndromes or regional networks that define these syndromes. Thus, patterns of laminar change may provide a useful anatomical framework for investigating how degeneration of select cells and corresponding laminar circuits influence large-scale networks and clinical symptomology in FTLD.

Plasma Neurofilament Light Concentration Is Associated with Diffusion-Tensor MRI-Based Measures of Neurodegeneration in Early Parkinson's Disease

BACKGROUND

Neurofilament light is a marker of axonal degeneration, whose measurement from peripheral blood was recently made possible by new assays.

OBJECTIVE

We aimed to determine whether plasma neurofilament light chain (NfL) concentration reflects brain white matter integrity in patients with early Parkinson's disease (PD).

METHODS

137 early PD patients and 51 healthy controls were included. Plasma NfL levels were measured using ultrasensitive single molecule array. 3T MRI including diffusion tensor imaging was acquired for voxelwise analysis of association between NfL and both fractional anisotropy (FA) and mean diffusivity (MD) in white matter tracts and subcortical nuclei.

RESULTS

A pattern of brain microstructural changes consistent with neurodegeneration was associated with increased plasma NfL in most of the frontal lobe and right internal capsule, with decreased FA and increased MD. The same clusters were also associated with poorer global cognition. A significant cluster in the left putamen was associated with increased NfL, with a significantly greater effect in PD than controls.

CONCLUSION

Plasma NfL may be associated with brain microstructure, as measured using diffusion tensor imaging, in patients with early PD. Higher plasma NfL was associated with a frontal pattern of neurodegeneration that also correlates with cognitive performance in our cohort. This may support a future role for plasma NfL as an accessible biomarker for neurodegeneration and cognitive dysfunction in PD.

Slab boundary artifact correction in multislab imaging using convolutional-neural-network-enabled inversion for slab profile encoding

PURPOSE

This study aims to propose a novel algorithm for slab boundary artifact correction in both single-band multislab imaging and simultaneous multislab (SMSlab) imaging.

THEORY AND METHODS

In image domain, the formation of slab boundary artifacts can be regarded as modulating the artifact-free images using the slab profiles and introducing aliasing along the slice direction. Slab boundary artifact correction is the inverse problem of this process. An iterative algorithm based on convolutional neural networks (CNNs) is proposed to solve the problem, termed CNN-enabled inversion for slab profile encoding (CPEN). Diffusion-weighted SMSlab images and reference images without slab boundary artifacts were acquired in 7 healthy subjects for training. Images of 5 healthy subjects were acquired for testing, including single-band multislab and SMSlab images with 1.3-mm or 1-mm isotropic resolution. CNN-enabled inversion for slab profile encoding was compared with a previously reported method (i.e., nonlinear inversion for slab profile encoding [NPEN]).

RESULTS

CNN-enabled inversion for slab profile encoding reduces the slab boundary artifacts in both single-band multislab and SMSlab images. It also suppresses the slab boundary artifacts in the diffusion metric maps. Compared with NPEN, CPEN shows fewer residual artifacts in different acquisition protocols and more significant improvements in quantitative assessment, and it also accelerates the computation by more than 35 times.

CONCLUSION

CNN-enabled inversion for slab profile encoding can reduce the slab boundary artifacts in multislab acquisitions. It shows better slab boundary artifact correction capacity, higher robustness, and computation efficiency when compared with NPEN. It has the potential to improve the accuracy of multislab acquisitions in high-resolution DWI and functional MRI.

On the generalizability of diffusion MRI signal representations across acquisition parameters, sequences and tissue types: Chronicles of the MEMENTO challenge

Diffusion MRI (dMRI) has become an invaluable tool to assess the microstructural organization of brain tissue. Depending on the specific acquisition settings, the dMRI signal encodes specific properties of the underlying diffusion process. In the last two decades, several signal representations have been proposed to fit the dMRI signal and decode such properties. Most methods, however, are tested and developed on a limited amount of data, and their applicability to other acquisition schemes remains unknown. With this work, we aimed to shed light on the generalizability of existing dMRI signal representations to different diffusion encoding parameters and brain tissue types. To this end, we organized a community challenge - named MEMENTO, making available the same datasets for fair comparisons across algorithms and techniques. We considered two state-of-the-art diffusion datasets, including single-diffusion-encoding (SDE) spin-echo data from a human brain with over 3820 unique diffusion weightings (the MASSIVE dataset), and double (oscillating) diffusion encoding data (DDE/DODE) of a mouse brain including over 2520 unique data points. A subset of the data sampled in 5 different voxels was openly distributed, and the challenge participants were asked to predict the remaining part of the data. After one year, eight participant teams submitted a total of 80 signal fits. For each submission, we evaluated the mean squared error, the variance of the prediction error and the Bayesian information criteria. The received submissions predicted either multi-shell SDE data (37%) or DODE data (22%), followed by cartesian SDE data (19%) and DDE (18%). Most submissions predicted the signals measured with SDE remarkably well, with the exception of low and very strong diffusion weightings. The prediction of DDE and DODE data seemed more challenging, likely because none of the submissions explicitly accounted for diffusion time and frequency. Next to the choice of the model, decisions on fit procedure and hyperparameters play a major role in the prediction performance, highlighting the importance of optimizing and reporting such choices. This work is a community effort to highlight strength and limitations of the field at representing dMRI acquired with trending encoding schemes, gaining insights into how different models generalize to different tissue types and fiber configurations over a large range of diffusion encodings.

Myelin densities in retinotopically defined dorsal visual areas of the macaque

The visuotopic organization of dorsal visual cortex rostral to area V2 in primates has been a longstanding source of controversy. Using sub-millimeter phase-encoded retinotopic fMRI mapping, we recently provided evidence for a surprisingly similar visuotopic organization in dorsal visual cortex of macaques compared to previously published maps in New world monkeys (Zhu and Vanduffel, Proc Natl Acad Sci USA 116:2306–2311, 2019). Although individual quadrant representations could be robustly delineated in that study, their grouping into hemifield representations remains a major challenge. Here, we combined in-vivo high-resolution myelin density mapping based on MR imaging (400 µm isotropic resolution) with fine-grained retinotopic fMRI to quantitatively compare myelin densities across retinotopically defined visual areas in macaques. Complementing previously documented differences in populational receptive-field (pRF) size and visual field signs, myelin densities of both quadrants of the dorsolateral posterior area (DLP) and area V3A are significantly different compared to dorsal and ventral area V3. Moreover, no differences in myelin density were observed between the two matching quadrants belonging to areas DLP, V3A, V1, V2 and V4, respectively. This was not the case, however, for the dorsal and ventral quadrants of area V3, which showed significant differences in MR-defined myelin densities, corroborating evidence of previous myelin staining studies. Interestingly, the pRF sizes and visual field signs of both quadrant representations in V3 are not different. Although myelin density correlates with curvature and anticorrelates with cortical thickness when measured across the entire cortex, exactly as in humans, the myelin density results in the visual areas cannot be explained by variability in cortical thickness and curvature between these areas. The present myelin density results largely support our previous model to group the two quadrants of DLP and V3A, rather than grouping DLP- with V3v into a single area VLP, or V3d with V3A+ into DM.

On the potential for mapping apparent neural soma density via a clinically viable diffusion MRI protocol

•

B-tensor encoding enables estimation of spherical cellular structures in the brain.

•

Spherical compartments may provide markers for apparent neural soma density.

•

Model parameters can be estimated in a fast and robust way using deep learning.

•

Practical acquisition times are achievable on widely available clinical scanners.

Diffusion MRI is a valuable tool for probing tissue microstructure in the brain noninvasively. Today, model-based techniques are widely available and used for white matter characterisation where their development is relatively mature. Conversely, tissue modelling in grey matter is more challenging, and no generally accepted models exist. With advances in measurement technology and modelling efforts, a clinically viable technique that reveals salient features of grey matter microstructure, such as the density of quasi-spherical cell bodies and quasi-cylindrical cell projections, is an exciting prospect. As a step towards capturing the microscopic architecture of grey matter in clinically feasible settings, this work uses a biophysical model that is designed to disentangle the diffusion signatures of spherical and cylindrical structures in the presence of orientation heterogeneity, and takes advantage of B-tensor encoding measurements, which provide additional sensitivity compared to standard single diffusion encoding sequences. For the fast and robust estimation of microstructural parameters, we leverage recent advances in machine learning and replace conventional fitting techniques with an artificial neural network that fits complex biophysical models within seconds. Our results demonstrate apparent markers of spherical and cylindrical geometries in healthy human subjects, and in particular an increased volume fraction of spherical compartments in grey matter compared to white matter. We evaluate the extent to which spherical and cylindrical geometries may be interpreted as correlates of neural soma and neural projections, respectively, and quantify parameter estimation errors in the presence of various departures from the modelling assumptions. While further work is necessary to translate the ideas presented in this work to the clinic, we suggest that biomarkers focussing on quasi-spherical cellular geometries may be valuable for the enhanced assessment of neurodevelopmental disorders and neurodegenerative diseases.

Microstructural MRI Correlates of Cognitive Impairment in Multiple Sclerosis: The Role of Deep Gray Matter

Although cognitive impairment (CI) is frequently observed in people with multiple sclerosis (pwMS), its pathogenesis is still controversial. Conflicting results emerged concerning the role of microstructural gray matter (GM) damage especially when involving the deep GM structures. In this study, we aimed at evaluating whether differences in cortical and deep GM structures between apparently cognitively normal (ACN) and CI pwMS (36 subjects in total) are present, using an extensive set of diffusion MRI (dMRI) indices and conventional morphometry measures. The results revealed increased anisotropy and restriction over several deep GM structures in CI compared with ACN pwMS, while no changes in volume were present in the same areas. Conversely, reduced anisotropy/restriction values were detected in cortical regions, mostly the pericalcarine cortex and precuneus, combined with reduced thickness of the superior frontal gyrus and insula. Most of the dMRI metrics but none of the morphometric indices correlated with the Symbol Digit Modality Test. These results suggest that deep GM microstructural damage can be a strong anatomical substrate of CI in pwMS and might allow identifying pwMS at higher risk of developing CI.

Validating layer-specific VASO across species

Cerebral blood volume (CBV) has been shown to be a robust and important physiological parameter for quantitative interpretation of functional (f)MRI, capable of delivering highly localized mapping of neural activity. Indeed, with recent advances in ultra-high-field (≥7T) MRI hardware and associated sequence libraries, it has become possible to capture non-invasive CBV weighted fMRI signals across cortical layers. One of the most widely used approaches to achieve this (in humans) is through vascular-space-occupancy (VASO) fMRI. Unfortunately, the exact contrast mechanisms of layer-dependent VASO fMRI have not been validated for human fMRI and thus interpretation of such data is confounded. Here we validate the signal source of layer-dependent SS-SI VASO fMRI using multi-modal imaging in a rat model in response to neuronal activation (somatosensory cortex) and respiratory challenge (hypercapnia). In particular VASO derived CBV measures are directly compared to concurrent measures of total haemoglobin changes from high resolution intrinsic optical imaging spectroscopy (OIS). Quantified cortical layer profiling is demonstrated to be in agreement between VASO and contrast enhanced fMRI (using monocrystalline iron oxide nanoparticles, MION). Responses show high spatial localisation to layers of cortical processing independent of confounding large draining veins which can hamper BOLD fMRI studies, (depending on slice positioning). Thus, a cross species comparison is enabled using VASO as a common measure. We find increased VASO based CBV reactivity (3.1 ± 1.2 fold increase) in humans compared to rats. Together, our findings confirm that the VASO contrast is indeed a reliable estimate of layer-specific CBV changes. This validation study increases the neuronal interpretability of human layer-dependent VASO fMRI as an appropriate method in neuroscience application studies, in which the presence of large draining intracortical and pial veins limits neuroscientific inference with BOLD fMRI.

Automated cerebral cortex segmentation based solely on diffusion tensor imaging for investigating cortical anisotropy

To extract Diffusion Tensor Imaging (DTI) parameters from the human cortex, the inner and outer boundaries of the cortex are usually defined on 3D-T1-weighted images and then applied to the co-registered DTI. However, this analysis requires the acquisition of an additional high-resolution structural image that may not be practical in various imaging studies. Here an automatic cortical boundary segmentation method was developed to work directly only on the native DTI images by using fractional anisotropy (FA) maps and mean diffusion weighted images (DWI), the latter with acceptable gray-white matter image contrast. This new method was compared to the conventional cortical segmentations generated from high-resolution T1 structural images in 5 participants. In addition, the proposed method was applied to 15 healthy young adults (10 cross-sectional, 5 test-retest) to measure FA, MD, and radiality of the primary eigenvector across the cortex on whole-brain 1.5 mm isotropic images acquired in 3.5 min at 3T. The proposed method generated reasonable segmentations of the cortical boundaries for all individuals and large proportions of the proposed method segmentations (more than 85%) were within ±1 mm from those generated with the conventional approach on higher resolution T1 structural images. Both FA (~0.15) and MD (~0.77 × 10^-3 mm²/s) extracted halfway between the cortical boundaries were relatively stable across the cortex, although focal regions such as the posterior bank of the central sulcus, anterior insula, and medial temporal lobe showed higher FA. The primary eigenvectors were primarily oriented radially to the middle cortical surface, but there were tangential orientations in the sulcal fundi as well as in the posterior bank of the central sulcus. The proposed method demonstrates the feasibility and accuracy of cortical analysis in native DTI space while avoiding the acquisition of other imaging contrasts like 3D T1-weighted scans.

Exercise Training-Related Changes in Cortical Gray Matter Diffusivity and Cognitive Function in Mild Cognitive Impairment and Healthy Older Adults

Individuals with Mild Cognitive Impairment (MCI) are at an elevated risk of dementia and exhibit deficits in cognition and cortical gray matter (GM) volume, thickness, and microstructure. Meanwhile, exercise training appears to preserve brain function and macrostructure may help delay or prevent the onset of dementia in individuals with MCI. Yet, our understanding of the neurophysiological effects of exercise training in individuals with MCI remains limited. Recent work suggests that the measures of gray matter microstructure using diffusion imaging may be sensitive to early cognitive and neurophysiological changes in the aging brain. Therefore, this study is aimed to determine the effects of exercise training in cognition and cortical gray matter microstructure in individuals with MCI vs. cognitively healthy older adults. Fifteen MCI participants and 17 cognitively intact controls (HC) volunteered for a 12-week supervised walking intervention. Following the intervention, MCI and HC saw improvements in cardiorespiratory fitness, performance on Trial 1 of the Rey Auditory Verbal Learning Test (RAVLT), a measure of verbal memory, and the Controlled Oral Word Association Test (COWAT), a measure of verbal fluency. After controlling for age, a voxel-wise analysis of cortical gray matter diffusivity showed individuals with MCI exhibited greater increases in mean diffusivity (MD) in the left insular cortex than HC. This increase in MD was positively associated with improvements in COWAT performance. Additionally, after controlling for age, the voxel-wise analysis indicated a main effect of Time with both groups experiencing an increase in left insular and left and right cerebellar MD. Increases in left insular diffusivity were similarly found to be positively associated with improvements in COWAT performance in both groups, while increases in cerebellar MD were related to gains in episodic memory performance. These findings suggest that exercise training may be related to improvements in neural circuits that govern verbal fluency performance in older adults through the microstructural remodeling of cortical gray matter. Furthermore, changes in left insular cortex microstructure may be particularly relevant to improvements in verbal fluency among individuals diagnosed with MCI.

Distortion-free, high-isotropic-resolution diffusion MRI with gSlider BUDA-EPI and multicoil dynamic B0 shimming

PURPOSE

We combine SNR-efficient acquisition and model-based reconstruction strategies with newly available hardware instrumentation to achieve distortion-free in vivo diffusion MRI of the brain at submillimeter-isotropic resolution with high fidelity and sensitivity on a clinical 3T scanner.

METHODS

We propose blip-up/down acquisition (BUDA) for multishot EPI using interleaved blip-up/blip-down phase encoding and incorporate B₀ forward-modeling into structured low-rank reconstruction to enable distortion-free and navigator-free diffusion MRI. We further combine BUDA-EPI with an SNR-efficient simultaneous multislab acquisition (generalized slice-dithered enhanced resolution ["gSlider"]), to achieve high-isotropic-resolution diffusion MRI. To validate gSlider BUDA-EPI, whole-brain diffusion data at 860-μm and 780-μm data sets were acquired. Finally, to improve the conditioning and minimize noise penalty in BUDA reconstruction at very high resolutions where B₀ inhomogeneity can have a detrimental effect, the level of B₀ inhomogeneity was reduced by incorporating slab-by-slab dynamic shimming with a 32-channel AC/DC coil into the acquisition. Whole-brain 600-μm diffusion data were then acquired with this combined approach of gSlider BUDA-EPI with dynamic shimming.

RESULTS

The results of 860-μm and 780-μm datasets show high geometry fidelity with gSlider BUDA-EPI. With dynamic shimming, the BUDA reconstruction's noise penalty was further alleviated. This enables whole-brain 600-μm isotropic resolution diffusion imaging with high image quality.

CONCLUSIONS

The gSlider BUDA-EPI method enables high-quality, distortion-free diffusion imaging across the whole brain at submillimeter resolution, where the use of multicoil dynamic B₀ shimming further improves reconstruction performance, which can be particularly useful at very high resolutions.