Relating quantitative 7T MRI across cortical depths to cytoarchitectonics, gene expression and connectomics

Unraveling the major role of vascular (R2') contributions to R2* signal relaxation at ultra-high-field MRI: A comprehensive analysis with quantitative gradient recalled echo in mouse brain

PURPOSE

Ultra-high-field (UHF) R2* relaxometry is often used for in vivo analysis of biological tissue microstructure without accounting for vascular contributions to R2* signal, that is, the BOLD signal component, and magnetic field inhomogeneities. These effects are especially important at UHF as their contribution to R2* scales linearly with magnetic field. Our study aims to report on the results of separate contributions of R2t* (tissue-specific sub-component) and R2' (vascular BOLD sub-component), corrected for the adverse effects of magnetic field inhomogeneities, to the total R2* signal at in vivo UHF MRI of mouse brain.

METHODS

Four healthy, 8-week-old C57BL/6J mice were imaged in vivo with multi-gradient echo MRI at 9.4 T and analyzed using the quantitative gradient recalled echo (qGRE) approach. A segmentation protocol was established using the Dorr Mouse Brain Atlas and ANTs Syn registration to warp template brain region labels to subject qGRE maps.

RESULTS

By separating R2' contribution from R2* signal, we have established normative R2t* data in mouse brain. Our findings revealed significant contributions of R2' to R2*, with approximately 42% of the R2* signal arising from vascular contributions, thus suggesting the R2t* as a more accurate metric for quantifying tissue microstructural information and its changes in neurodegenerative diseases.

CONCLUSION

qGRE approach allows efficient separation of tissue microstructure-specific (R2t*), vascular BOLD (R2'), and background gradients contributions to the total R2* relaxation at UHF MRI. Due to low concentration of non-heme iron in mouse brain, major contribution to R2t* results from tissue cellular components.

Neuroimaging Techniques in Huntington's Disease: A Critical Review

BACKGROUND

Huntington's disease (HD) is a hereditary neurodegenerative disorder characterized by cognitive, neuropsychiatric and motor symptoms caused by a CAG trinucleotide repeat expansion in the huntingtin gene. Imaging techniques are crucial for understanding HD pathophysiology and monitoring disease progression.

OBJECTIVES

This review is targeted at general neurologists and movement disorders specialists with an interest in HD and aims to bring complex imaging, including new experimental techniques, closer to the practicing clinician.

METHODS

We provide a summary of findings from conventional structural, diffusion and functional imaging in HD studies, together with an update on emerging novel techniques, including multiparametric mapping, multi-shell diffusion techniques, ultra-high field 7-Tesla MRI, positron emission tomography and magnetoencephalography.

RESULTS

Conventional imaging techniques have deepened our understanding of neuropathological progression in HD, from striatal atrophy to widespread cortical and white matter changes. The integration of novel imaging techniques reviewed has further improved our ability to interrogate, quantify and visualize disease-specific alterations with high precision.

CONCLUSIONS

Novel imaging techniques have promising roles to further our understanding of HD pathology and as imaging markers for clinical trials, disease staging and therapeutic monitoring. Additionally, the synergistic potential of combining imaging modalities with molecular and genetic data, along with wet biomarkers and clinical data, will help provide a complete and comprehensive view of HD pathology and progression.

Closed-loop fMRI at the mesoscopic scale of columns and layers: Can we do it and why would we want to?

Technological advances in fMRI including ultra-high magnetic fields (≥ 7 T) and acquisition methods that increase spatial specificity have paved the way for studies of the human cortex at the scale of layers and columns. This mesoscopic scale promises an improved mechanistic understanding of human cortical function so far only accessible to invasive animal neurophysiology. In recent years, an increasing number of studies have applied such methods to better understand the cortical function in perception and cognition. This future perspective article asks whether closed-loop fMRI studies could equally benefit from these methods to achieve layer and columnar specificity. We outline potential applications and discuss the conceptual and concrete challenges, including data acquisition and volitional control of mesoscopic brain activity. We anticipate an important role of fMRI with mesoscopic resolution for closed-loop fMRI and neurofeedback, yielding new insights into brain function and potentially clinical applications.This article is part of the theme issue 'Neurofeedback: new territories and neurocognitive mechanisms of endogenous neuromodulation'.

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 intra-cortical 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 notably high homology with those in humans, suggesting cortical myelin shares a similar developmental program across 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.

Segmentation of supragranular and infragranular layers in ultra-high-resolution 7T ex vivo MRI of the human cerebral cortex

Accurate labeling of specific layers in the human cerebral cortex is crucial for advancing our understanding of neurodevelopmental and neurodegenerative disorders. Building on recent advancements in ultra-high-resolution ex vivo MRI, we present a novel semi-supervised segmentation model capable of identifying supragranular and infragranular layers in ex vivo MRI with unprecedented precision. On a dataset consisting of 17 whole-hemisphere ex vivo scans at 120 $\mu $m, we propose a Multi-resolution U-Nets framework that integrates global and local structural information, achieving reliable segmentation maps of the entire hemisphere, with Dice scores over 0.8 for supra- and infragranular layers. This enables surface modeling, atlas construction, anomaly detection in disease states, and cross-modality validation while also paving the way for finer layer segmentation. Our approach offers a powerful tool for comprehensive neuroanatomical investigations and holds promise for advancing our mechanistic understanding of progression of neurodegenerative diseases.

7T MP2RAGE for cortical myelin segmentation: Impact of aging

BACKGROUND

Myelin and iron are major contributors to the cortical MR signal. The aim of this study was to investigate 1. Can MP2RAGE-derived contrasts at 7T in combination with k-means clustering be used to distinguish between heavily and sparsely myelinated layers in cortical gray matter (GM)? 2. Does this approach provide meaningful biological information?

METHODS

The following contrasts were generated from the 7T MP2RAGE images from 45 healthy controls (age: 19-75, f/m = 23/22) from the ATAG data repository: 1. T1 weighted image (UNI). 2. T1 relaxation image (T1map). 3. INVC/T1map ratio (RATIO). K-means clustering identified 6 clusters/tissue maps (csf, csf/gm-transition, wm, wm/gm transition, heavily myelinated cortical GM (dGM), sparsely myelinated cortical GM (sGM)). These tissue maps were then processed with SPM/DARTEL (volume-based analyses) and Freesurfer (surface-based analyses) and dGM and sGM volume/thickness of young adults (n = 27, 19-27 years) compared to those of older adults (n = 18, 42-75 years) at p<0.001 uncorrected.

RESULTS

The resulting maps showed good agreement with histological maps in the literature. Volume- and surface analyses found age-related dGM loss/thinning in the mid-posterior cingulate and parahippocampal/entorhinal gyrus and age-related sGM losses in lateral, mesial and orbitofrontal frontal, insular cortex and superior temporal gyrus.

CONCLUSION

The MP2RAGE derived UNI, T1map and RATIO contrasts can be used to identify dGM and sGM. Considering the close relationship between cortical myelo- and cytoarchitecture, the findings reported here indicate that this new technique might provide new insights into the nature of cortical GM loss in physiological and pathological conditions.

High-resolution quantitative and functional MRI indicate lower myelination of thin and thick stripes in human secondary visual cortex

The characterization of cortical myelination is essential for the study of structure-function relationships in the human brain. However, knowledge about cortical myelination is largely based on post-mortem histology, which generally renders direct comparison to function impossible. The repeating pattern of pale-thin-pale-thick stripes of cytochrome oxidase (CO) activity in the primate secondary visual cortex (V2) is a prominent columnar system, in which histology also indicates different myelination of thin/thick versus pale stripes. We used quantitative magnetic resonance imaging (qMRI) in conjunction with functional magnetic resonance imaging (fMRI) at ultra-high field strength (7 T) to localize and study myelination of stripes in four human participants at sub-millimeter resolution in vivo. Thin and thick stripes were functionally localized by exploiting their sensitivity to color and binocular disparity, respectively. Resulting functional activation maps showed robust stripe patterns in V2 which enabled further comparison of quantitative relaxation parameters between stripe types. Thereby, we found lower longitudinal relaxation rates (R1) of thin and thick stripes compared to surrounding gray matter in the order of 1-2%, indicating higher myelination of pale stripes. No consistent differences were found for effective transverse relaxation rates (R2*). The study demonstrates the feasibility to investigate structure-function relationships in living humans within one cortical area at the level of columnar systems using qMRI.

The Potential of Myelin-Sensitive Imaging: Redefining Spatiotemporal Patterns of Myeloarchitecture

Recent advances in magnetic resonance imaging (MRI) have paved the way for approximation of myelin content in vivo. In this review, our main goal was to determine how to best capitalize on myelin-sensitive imaging. First, we briefly overview the theoretical and empirical basis for the myelin sensitivity of different MRI markers and, in doing so, highlight how multimodal imaging approaches are important for enhancing specificity to myelin. Then, we discuss recent studies that have probed the nonuniform distribution of myelin across cortical layers and along white matter tracts. These approaches, collectively known as myelin profiling, have provided detailed depictions of myeloarchitecture in both the postmortem and living human brain. Notably, MRI-based profiling studies have recently focused on investigating whether it can capture interindividual variability in myelin characteristics as well as trajectories across the lifespan. Finally, another line of recent evidence emphasizes the contribution of region-specific myelination to large-scale organization, demonstrating the impact of myelination on global brain networks. In conclusion, we suggest that combining well-validated MRI markers with profiling techniques holds strong potential to elucidate individual differences in myeloarchitecture, which has important implications for understanding brain function and disease.

Quantitative MRI maps of human neocortex explored using cell type-specific gene expression analysis

Quantitative magnetic resonance imaging (qMRI) allows extraction of reproducible and robust parameter maps. However, the connection to underlying biological substrates remains murky, especially in the complex, densely packed cortex. We investigated associations in human neocortex between qMRI parameters and neocortical cell types by comparing the spatial distribution of the qMRI parameters longitudinal relaxation rate (${R_{1}}$), effective transverse relaxation rate (${R_{2}}^{\ast }$), and magnetization transfer saturation (MTsat) to gene expression from the Allen Human Brain Atlas, then combining this with lists of genes enriched in specific cell types found in the human brain. As qMRI parameters are magnetic field strength-dependent, the analysis was performed on MRI data at 3T and 7T. All qMRI parameters significantly covaried with genes enriched in GABA- and glutamatergic neurons, i.e. they were associated with cytoarchitecture. The qMRI parameters also significantly covaried with the distribution of genes enriched in astrocytes (${R_{2}}^{\ast }$ at 3T, ${R_{1}}$ at 7T), endothelial cells (${R_{1}}$ and MTsat at 3T), microglia (${R_{1}}$ and MTsat at 3T, ${R_{1}}$ at 7T), and oligodendrocytes and oligodendrocyte precursor cells (${R_{1}}$ at 7T). These results advance the potential use of qMRI parameters as biomarkers for specific cell types.

Error quantification in multi-parameter mapping facilitates robust estimation and enhanced group level sensitivity

Multi-Parameter Mapping (MPM) is a comprehensive quantitative neuroimaging protocol that enables estimation of four physical parameters (longitudinal and effective transverse relaxation rates R₁ and R₂^*, proton density PD, and magnetization transfer saturation MT_sat) that are sensitive to microstructural tissue properties such as iron and myelin content. Their capability to reveal microstructural brain differences, however, is tightly bound to controlling random noise and artefacts (e.g. caused by head motion) in the signal. Here, we introduced a method to estimate the local error of PD, R₁, and MT_sat maps that captures both noise and artefacts on a routine basis without requiring additional data. To investigate the method's sensitivity to random noise, we calculated the model-based signal-to-noise ratio (mSNR) and showed in measurements and simulations that it correlated linearly with an experimental raw-image-based SNR map. We found that the mSNR varied with MPM protocols, magnetic field strength (3T vs. 7T) and MPM parameters: it halved from PD to R₁ and decreased from PD to MT_sat by a factor of 3-4. Exploring the artefact-sensitivity of the error maps, we generated robust MPM parameters using two successive acquisitions of each contrast and the acquisition-specific errors to down-weight erroneous regions. The resulting robust MPM parameters showed reduced variability at the group level as compared to their single-repeat or averaged counterparts. The error and mSNR maps may better inform power-calculations by accounting for local data quality variations across measurements. Code to compute the mSNR maps and robustly combined MPM maps is available in the open-source hMRI toolbox.

Neuroimaging at 7 Tesla: a pictorial narrative review

Neuroimaging using the 7-Tesla (7T) human magnetic resonance (MR) system is rapidly gaining popularity after being approved for clinical use in the European Union and the USA. This trend is the same for functional MR imaging (MRI). The primary advantages of 7T over lower magnetic fields are its higher signal-to-noise and contrast-to-noise ratios, which provide high-resolution acquisitions and better contrast, making it easier to detect lesions and structural changes in brain disorders. Another advantage is the capability to measure a greater number of neurochemicals by virtue of the increased spectral resolution. Many structural and functional studies using 7T have been conducted to visualize details in the white matter and layers of the cortex and hippocampus, the subnucleus or regions of the putamen, the globus pallidus, thalamus and substantia nigra, and in small structures, such as the subthalamic nucleus, habenula, perforating arteries, and the perivascular space, that are difficult to observe at lower magnetic field strengths. The target disorders for 7T neuroimaging range from tumoral diseases to vascular, neurodegenerative, and psychiatric disorders, including Alzheimer's disease, Parkinson's disease, multiple sclerosis, epilepsy, major depressive disorder, and schizophrenia. MR spectroscopy has also been used for research because of its increased chemical shift that separates overlapping peaks and resolves neurochemicals more effectively at 7T than a lower magnetic field. This paper presents a narrative review of these topics and an illustrative presentation of images obtained at 7T. We expect 7T neuroimaging to provide a new imaging biomarker of various brain disorders.

Combining navigator and optical prospective motion correction for high-quality 500 μm resolution quantitative multi-parameter mapping at 7T

PURPOSE

High-resolution quantitative multi-parameter mapping shows promise for non-invasively characterizing human brain microstructure but is limited by physiological artifacts. We implemented corrections for rigid head movement and respiration-related B0-fluctuations and evaluated them in healthy volunteers and dementia patients.

METHODS

Camera-based optical prospective motion correction (PMC) and FID navigator correction were implemented in a gradient and RF-spoiled multi-echo 3D gradient echo sequence for mapping proton density (PD), longitudinal relaxation rate (R1) and effective transverse relaxation rate (R2*). We studied their effectiveness separately and in concert in young volunteers and then evaluated the navigator correction (NAVcor) with PMC in a group of elderly volunteers and dementia patients. We used spatial homogeneity within white matter (WM) and gray matter (GM) and scan-rescan measures as quality metrics.

RESULTS

NAVcor and PMC reduced artifacts and improved the homogeneity and reproducibility of parameter maps. In elderly participants, NAVcor improved scan-rescan reproducibility of parameter maps (coefficient of variation decreased by 14.7% and 11.9% within WM and GM respectively). Spurious inhomogeneities within WM were reduced more in the elderly than in the young cohort (by 9% vs. 2%). PMC increased regional GM/WM contrast and was especially important in the elderly cohort, which moved twice as much as the young cohort. We did not find a significant interaction between the two corrections.

CONCLUSION

Navigator correction and PMC significantly improved the quality of PD, R1, and R2* maps, particularly in less compliant elderly volunteers and dementia patients.

Rhes protein transits from neuron to neuron and facilitates mutant huntingtin spreading in the brain

Rhes (RASD2) is a thyroid hormone-induced gene that regulates striatal motor activity and promotes neurodegeneration in Huntington disease (HD) and tauopathy. Rhes moves and transports the HD protein, polyglutamine-expanded huntingtin (mHTT), via tunneling nanotube (TNT)-like membranous protrusions between cultured neurons. However, similar intercellular Rhes transportation in the intact brain was unknown. Here, we report that Rhes induces TNT-like protrusions in the striatal medium spiny neurons (MSNs) and transported between dopamine-1 receptor (D1R)-MSNs and D2R-MSNs of intact striatum and organotypic brain slices. Notably, mHTT is robustly transported within the striatum and from the striatum to the cortical areas in the brain, and Rhes deletion diminishes such transport. Moreover, Rhes moves to the cortical regions following restricted expression in the MSNs of the striatum. Thus, Rhes is a first striatum-enriched protein demonstrated to move and transport mHTT between neurons and brain regions, providing new insights into interneuronal protein transport in the brain.

Cortical Cross-Frequency Coupling Is Affected by in utero Exposure to Antidepressant Medication

Up to five percent of human infants are exposed to maternal antidepressant medication by serotonin reuptake inhibitors (SRI) during pregnancy, yet the SRI effects on infants’ early neurodevelopment are not fully understood. Here, we studied how maternal SRI medication affects cortical frequency-specific and cross-frequency interactions estimated, respectively, by phase-phase correlations (PPC) and phase-amplitude coupling (PAC) in electroencephalographic (EEG) recordings. We examined the cortical activity in infants after fetal exposure to SRIs relative to a control group of infants without medical history of any kind. Our findings show that the sleep-related dynamics of PPC networks are selectively affected by in utero SRI exposure, however, those alterations do not correlate to later neurocognitive development as tested by neuropsychological evaluation at two years of age. In turn, phase-amplitude coupling was found to be suppressed in SRI infants across multiple distributed cortical regions and these effects were linked to their neurocognitive outcomes. Our results are compatible with the overall notion that in utero drug exposures may cause subtle, yet measurable changes in the brain structure and function. Our present findings are based on the measures of local and inter-areal neuronal interactions in the cortex which can be readily used across species, as well as between different scales of inspection: from the whole animals to in vitro preparations. Therefore, this work opens a framework to explore the cellular and molecular mechanisms underlying neurodevelopmental SRI effects at all translational levels.

Striatal Induction and Spread of the Huntington's Disease Protein: A Novel Rhes Route

The CAG/CAA expansion encoding polyQ huntingtin (mutant huntingtin [mHTT]) causes Huntington's disease (HD), which is characterized by atrophy and loss of striatal medium spiny neurons (MSNs), which are preceded by neuropathological alterations in the cortex. Previous studies have shown that mHTT can spread in the brain, but the mechanisms involved in the stereotyped degeneration and dysfunction of the neurons from the striatum to the cortex remain unclear. In this study, we found that the mHTT expression initially restricted in the striatum later spread to the cortical regions in mouse brains. Such transmission was diminished in mice that lacked the striatal-enriched protein Ras-homolog enriched in the striatum (Rhes). Rhes restricted to MSNs was also found in the cortical layers of the brain, indicating a new transmission route for the Rhes protein to the brain. Mechanistically, Rhes promotes such transmission via a direct cell-to-cell contact mediated by tunneling nanotubes (TNTs), the membranous protrusions that enable the transfer of mHTT, Rhes, and other vesicular cargoes. These transmission patterns suggest that Rhes and mHTT are likely co-transported in the brain using TNT-like cell-to-cell contacts. On the basis of these new results, a perspective is presented in this review: Rhes may ignite the mHTT transmission from the striatum that may coincide with HD onset and disease progression through an anatomically connected striato-cortical retrograde route.

Simulating Local Deformations in the Human Cortex Due to Blood Flow-Induced Changes in Mechanical Tissue Properties: Impact on Functional Magnetic Resonance Imaging

Investigating human brain tissue is challenging due to the complexity and the manifold interactions between structures across different scales. Increasing evidence suggests that brain function and microstructural features including biomechanical features are related. More importantly, the relationship between tissue mechanics and its influence on brain imaging results remains poorly understood. As an important example, the study of the brain tissue response to blood flow could have important theoretical and experimental consequences for functional magnetic resonance imaging (fMRI) at high spatial resolutions. Computational simulations, using realistic mechanical models can predict and characterize the brain tissue behavior and give us insights into the consequent potential biases or limitations of in vivo, high-resolution fMRI. In this manuscript, we used a two dimensional biomechanical simulation of an exemplary human gyrus to investigate the relationship between mechanical tissue properties and the respective changes induced by focal blood flow changes. The model is based on the changes in the brain’s stiffness and volume due to the vasodilation evoked by neural activity. Modeling an exemplary gyrus from a brain atlas we assessed the influence of different potential mechanisms: (i) a local increase in tissue stiffness (at the level of a single anatomical layer), (ii) an increase in local volume, and (iii) a combination of both effects. Our simulation results showed considerable tissue displacement because of these temporary changes in mechanical properties. We found that the local volume increase causes more deformation and consequently higher displacement of the gyrus. These displacements introduced considerable artifacts in our simulated fMRI measurements. Our results underline the necessity to consider and characterize the tissue displacement which could be responsible for fMRI artifacts.