Studying cyto and myeloarchitecture of the human cortex at ultra-high field with quantitative imaging: R1, R2* and magnetic susceptibility

Reproducibility of tailored and universal nonselective excitation pulses at 7 T for human cardiac MRI: A 3-year and an interday study

PURPOSE

Ultrahigh-field (UHF; ≥7 T) MRI is challenging due to spatially heterogeneous B1 + profiles. This longitudinal study evaluates the reproducibility of three parallel-transmission excitation strategies to enable UHF cardiac MRI: vendor-supplied radiofrequency (RF) shim, subject-tailored kT-points pulses (TPs), and universal kT-points pulses (UPs).

METHODS

Six healthy subjects underwent 7 T MRI scans performed by different MR operators using a 32-element parallel-transmission body array at four time points over 3 years. A single UP was computed and applied to all subjects. TPs were computed individually for each scan and organized into four configurations. Each configuration was applied to all scans from each subject to analyze intrasubject variability. Reproducibility was assessed by comparing the coefficient of variation (CV) of simulated flip angles (FAs) within the heart volume across scan sessions.

RESULTS

TPs designed for a specific scan session yielded lower CVs (2-fold reduction) than UP. Applying TPs to other scan sessions of the same subject, however, resulted in approximately 40% higher CVs and lower FA uniformity compared with the UP. On average, the UP consistently achieved the most reproducible results across inter-year, inter-day, and same-operator studies, with CVs of approximately 12%.

CONCLUSION

Although TPs showed advantages when tailored for a specific target volume, they struggled with long-term consistency and required lengthy calibration. The precomputed UP kT-points pulses proved to be the most consistent across all scans acquired in the 3 years by different operators, minimizing CV-data dispersion and maintaining FA uniformity.

In vivo human neurite exchange imaging (NEXI) at 500 mT/m diffusion gradients

Evaluating tissue microstructure and membrane integrity in the living human brain through diffusion-water exchange imaging is challenging due to requirements for a high signal-to-noise ratio and short diffusion times dictated by relatively fast exchange processes. The goal of this work was to demonstrate the feasibility of in vivo imaging of tissue micro-geometries and water exchange within the brain gray matter using the state-of-the-art Connectome 2.0 scanner equipped with an ultra-high-performance gradient system (maximum gradient strength=500 mT/m, maximum slew rate=600 T/m/s). We performed diffusion MRI measurements in 15 healthy volunteers at multiple diffusion times (13-30 ms) and b -values up to 17.5 ms/μm2. The anisotropic Kärger model was applied to estimate the exchange time between intra-neurite and extracellular water in gray matter. The estimated exchange time across the cortical ribbon was around (median±interquartile range) 13±8 ms on Connectome 2.0, substantially faster than that measured using an imaging protocol compatible with Connectome 1.0-alike systems on the same cohort. Our investigation suggested that the NEXI exchange time estimation using a Connectome 1.0 compatible protocol was more prone to residual noise floor biases due to the small time-dependent signal contrasts across diffusion times when the exchange is fast (≤20 ms). Furthermore, spatial variation of exchange time was observed across the cortex, where the motor cortex, somatosensory cortex and visual cortex exhibit longer exchange times compared to other cortical regions. Non-linear fitting for the anisotropic Kärger model was accelerated 100 times using a GPU-based pipeline compared to the conventional CPU-based approach. This study highlighted the importance of the chosen diffusion times and measures to address Rician noise in dMRI data, which can have a substantial impact on the estimated NEXI exchange time and require extra attention when comparing NEXI results between various hardware setups.

Rapid submillimeter QSM and R2* mapping using interleaved multishot 3D-EPI at 7 and 3 Tesla

PURPOSE

To explore the high signal-to-noise ratio (SNR) efficiency of interleaved multishot 3D-EPI with standard image reconstruction for fast and robust high-resolution whole-brain quantitative susceptibility (QSM) andR 2 ∗ $$ {R}_2^{\ast } $$ mapping at 7 and 3T.

METHODS

Single- and multi-TE segmented 3D-EPI is combined with conventional CAIPIRINHA undersampling for up to 72-fold effective gradient echo (GRE) imaging acceleration. Across multiple averages, scan parameters are varied (e.g., dual-polarity frequency-encoding) to additionally correct forB 0 $$ {\mathrm{B}}_0 $$ -induced artifacts, geometric distortions and motion retrospectively. A comparison to established GRE protocols is made. Resolutions range from 1.4 mm isotropic (1 multi-TE average in 36 s) up to 0.4 mm isotropic (2 single-TE averages in approximately 6 min) with whole-head coverage.

RESULTS

Only 1-4 averages are needed for sufficient SNR with 3D-EPI, depending on resolution and field strength. Fast scanning and small voxels together with retrospective corrections result in substantially reduced image artifacts, which improves susceptibility andR 2 ∗ $$ {R}_2^{\ast } $$ mapping. Additionally, much finer details are obtained in susceptibility-weighted image projections through significantly reduced partial voluming.

CONCLUSION

Using interleaved multishot 3D-EPI, single-TE and multi-TE data can readily be acquired 10 times faster than with conventional, accelerated GRE imaging. Even 0.4 mm isotropic whole-head QSM within 6 min becomes feasible at 7T. At 3T, motion-robust 0.8 mm isotropic whole-brain QSM andR 2 ∗ $$ {R}_2^{\ast } $$ mapping with no apparent distortion in less than 7 min becomes clinically feasible. Stronger gradient systems may allow for even higher effective acceleration rates through larger EPI factors while maintaining optimal contrast.

Quantitative MRI at 7-Tesla reveals novel frontocortical myeloarchitecture anomalies in major depressive disorder

Whereas meta-analytical data highlight abnormal frontocortical macrostructure (thickness/surface area/volume) in Major Depressive Disorder (MDD), the underlying microstructural processes remain uncharted, due to the use of conventional MRI scanners and acquisition techniques. We uniquely combined Ultra-High Field MRI at 7.0 Tesla with Quantitative Imaging to map intracortical myelin (proxied by longitudinal relaxation time T1) and iron concentration (proxied by transverse relaxation time T2*), microstructural processes deemed particularly germane to cortical macrostructure. Informed by meta-analytical evidence, we focused specifically on orbitofrontal and rostral anterior cingulate cortices among adult MDD patients (N = 48) and matched healthy controls (HC; N = 10). Analyses probed the association of MDD diagnosis and clinical profile (severity, medication use, comorbid anxiety disorders, childhood trauma) with aforementioned microstructural properties. MDD diagnosis (p's < 0.05, Cohen's D = 0.55-0.66) and symptom severity (p's < 0.01, r = 0.271-0.267) both related to decreased intracortical myelination (higher T1 values) within the lateral orbitofrontal cortex, a region tightly coupled to processing negative affect and feelings of sadness in MDD. No relations were found with local iron concentrations. These findings allow uniquely fine-grained insights on frontocortical microstructure in MDD, and cautiously point to intracortical demyelination as a possible driver of macroscale cortical disintegrity in MDD.

High-resolution myelin-water fraction and quantitative relaxation mapping using 3D ViSTa-MR fingerprinting

PURPOSE

This study aims to develop a high-resolution whole-brain multi-parametric quantitative MRI approach for simultaneous mapping of myelin-water fraction (MWF), T1, T2, and proton-density (PD), all within a clinically feasible scan time.

METHODS

We developed 3D visualization of short transverse relaxation time component (ViSTa)-MRF, which combined ViSTa technique with MR fingerprinting (MRF), to achieve high-fidelity whole-brain MWF and T1/T2/PD mapping on a clinical 3T scanner. To achieve fast acquisition and memory-efficient reconstruction, the ViSTa-MRF sequence leverages an optimized 3D tiny-golden-angle-shuffling spiral-projection acquisition and joint spatial-temporal subspace reconstruction with optimized preconditioning algorithm. With the proposed ViSTa-MRF approach, high-fidelity direct MWF mapping was achieved without a need for multicompartment fitting that could introduce bias and/or noise from additional assumptions or priors.

RESULTS

The in vivo results demonstrate the effectiveness of the proposed acquisition and reconstruction framework to provide fast multi-parametric mapping with high SNR and good quality. The in vivo results of 1 mm- and 0.66 mm-isotropic resolution datasets indicate that the MWF values measured by the proposed method are consistent with standard ViSTa results that are 30× slower with lower SNR. Furthermore, we applied the proposed method to enable 5-min whole-brain 1 mm-iso assessment of MWF and T1/T2/PD mappings for infant brain development and for post-mortem brain samples.

CONCLUSIONS

In this work, we have developed a 3D ViSTa-MRF technique that enables the acquisition of whole-brain MWF, quantitative T1, T2, and PD maps at 1 and 0.66 mm isotropic resolution in 5 and 15 min, respectively. This advancement allows for quantitative investigations of myelination changes in the brain.

Analyzing the cortical fine structure as revealed by ex-vivo anatomical MRI

The human cortex has a rich fiber structure as revealed by myelin-staining of histological slices. Myelin also contributes to the image contrast in Magnetic Resonance Imaging (MRI). Recent advances in Magnetic Resonance (MR) scanner and imaging technology allowed the acquisition of an ex-vivo data set at an isotropic resolution of 100 µm. This study focused on a computational analysis of this data set with the aim of bridging between histological knowledge and MRI-based results. This work highlights: (1) the design and implementation of a processing chain that extracts intracortical features from a high-resolution MR image; (2) a demonstration of the correspondence between MRI-based cortical intensity profiles and the myelo-architectonic layering of the cortex; (3) the characterization and classification of four basic myelo-architectonic profile types; (4) the distinction of cortical regions based on myelo-architectonic features; and (5) the segmentation of cortical modules in the entorhinal cortex.

Normative quantitative relaxation atlases for characterization of cortical regions using magnetic resonance fingerprinting

Quantitative magnetic resonance (MR) has been used to study cyto- and myelo-architecture of the human brain non-invasively. However, analyzing brain cortex using high-resolution quantitative MR acquisition can be challenging to perform using 3T clinical scanners. MR fingerprinting (MRF) is a highly efficient and clinically feasible quantitative MR technique that simultaneously provides T1 and T2 relaxation maps. Using 3D MRF from 40 healthy subjects (mean age = 25.6 ± 4.3 years) scanned on 3T magnetic resonance imaging, we generated whole-brain gyral-based normative MR relaxation atlases and investigated cortical-region-based T1 and T2 variations. Gender and age dependency of T1 and T2 variations were additionally analyzed. The coefficient of variation of T1 and T2 for each cortical-region was 3.5% and 7.3%, respectively, supporting low variability of MRF measurements across subjects. Significant differences in T1 and T2 were identified among 34 brain regions (P < 0.001), lower in the precentral, postcentral, paracentral lobule, transverse temporal, lateral occipital, and cingulate areas, which contain sensorimotor, auditory, visual, and limbic functions. Significant correlations were identified between age and T1 and T2 values. This study established whole-brain MRF T1 and T2 atlases of healthy subjects using a clinical 3T scanner, which can provide a quantitative and region-specific baseline for future brain studies and pathology detection.

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.

Rapid estimation of 2D relative B 1 + -maps from localizers in the human heart at 7T using deep learning

PURPOSE

Subject-tailored parallel transmission pulses for ultra-high fields body applications are typically calculated based on subject-specific B 1 + $$ {\mathrm{B}}_1^{+} $$ -maps of all transmit channels, which require lengthy adjustment times. This study investigates the feasibility of using deep learning to estimate complex, channel-wise, relative 2D B 1 + $$ {\mathrm{B}}_1^{+} $$ -maps from a single gradient echo localizer to overcome long calibration times.

METHODS

126 channel-wise, complex, relative 2D B 1 + $$ {\mathrm{B}}_1^{+} $$ -maps of the human heart from 44 subjects were acquired at 7T using a Cartesian, cardiac gradient-echo sequence obtained under breath-hold to create a library for network training and cross-validation. The deep learning predicted maps were qualitatively compared to the ground truth. Phase-only B 1 + $$ {\mathrm{B}}_1^{+} $$ -shimming was subsequently performed on the estimated B 1 + $$ {\mathrm{B}}_1^{+} $$ -maps for a region of interest covering the heart. The proposed network was applied at 7T to 3 unseen test subjects.

RESULTS

The deep learning-based B 1 + $$ {\mathrm{B}}_1^{+} $$ -maps, derived in approximately 0.2 seconds, match the ground truth for the magnitude and phase. The static, phase-only pulse design performs best when maximizing the mean transmission efficiency. In-vivo application of the proposed network to unseen subjects demonstrates the feasibility of this approach: the network yields predicted B 1 + $$ {\mathrm{B}}_1^{+} $$ -maps comparable to the acquired ground truth and anatomical scans reflect the resulting B 1 + $$ {\mathrm{B}}_1^{+} $$ -pattern using the deep learning-based maps.

CONCLUSION

The feasibility of estimating 2D relative B 1 + $$ {\mathrm{B}}_1^{+} $$ -maps from initial localizer scans of the human heart at 7T using deep learning is successfully demonstrated. Because the technique requires only sub-seconds to derive channel-wise B 1 + $$ {\mathrm{B}}_1^{+} $$ -maps, it offers high potential for advancing clinical body imaging at ultra-high fields.

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.

Towards automated in vivo parcellation of the human cerebral cortex using supervised classification of magnetic resonance fingerprinting residuals

BACKGROUND

Accurate parcellation of the cerebral cortex in an individual is a guide to its underlying organization. The most promising in vivo quantitative magnetic resonance (MR)-based microstructural cortical mapping methods are yet to achieve a level of parcellation accuracy comparable to quantitative histology.

METHODS

We scanned 6 participants using a 3D echo-planar imaging MR fingerprinting (EPI-MRF) sequence on a 7T Siemens scanner. After projecting MRF signals to the individual-specific inflated model of the cortical surface, normalized autocorrelations of MRF residuals of vertices of 8 microstructurally distinct areas (BA1, BA2, BA4a, BA6, BA44, BA45, BA17, and BA18) from 3 cortical regions were used as feature vector inputs into linear support vector machine (SVM), radial basis function SVM (RBF-SVM), random forest, and k-nearest neighbors supervised classification algorithms. The algorithms' prediction performance was compared using: (i) features from each vertex or (ii) features from neighboring vertices.

RESULTS

The neighborhood-based RBF-SVM classifier achieved the highest prediction score of 0.85 for classification of MRF residuals in the central region from a held-out participant.

CONCLUSIONS

We developed an automated method of cortical parcellation using a combination of MR fingerprinting residual analysis and machine learning classification. Our findings provide the basis for employing unsupervised learning algorithms for whole-cortex structural parcellation in individuals.

Mesoscopic in vivo human T2* dataset acquired using quantitative MRI at 7 Tesla

Mesoscopic (0.1-0.5 mm) interrogation of the living human brain is critical for advancing neuroscience and bridging the resolution gap with animal models. Despite the variety of MRI contrasts measured in recent years at the mesoscopic scale, in vivo quantitative imaging of T₂^* has not been performed. Here we provide a dataset containing empirical T₂^* measurements acquired at 0.35 × 0.35 × 0.35 mm³ voxel resolution using 7 Tesla MRI. To demonstrate unique features and high quality of this dataset, we generate flat map visualizations that reveal fine-scale cortical substructures such as layers and vessels, and we report quantitative depth-dependent T₂^* (as well as R₂^*) values in primary visual cortex and auditory cortex that are highly consistent across subjects. This dataset is freely available at https://doi.org/10.17605/OSF.IO/N5BJ7, and may prove useful for anatomical investigations of the human brain, as well as for improving our understanding of the basis of the T₂^*-weighted (f)MRI signal.

The interindividual variability of multimodal brain connectivity maintains spatial heterogeneity and relates to tissue microstructure

Humans differ from each other in a wide range of biometrics, but to what extent brain connectivity varies between individuals remains largely unknown. By combining diffusion-weighted imaging (DWI) and magnetoencephalography (MEG), this study characterizes the inter-subject variability (ISV) of multimodal brain connectivity. Structural connectivity is characterized by higher ISV in association cortices including the core multiple-demand network and lower ISV in the sensorimotor cortex. MEG ISV exhibits frequency-dependent signatures, and the extent of MEG ISV is consistent with that of structural connectivity ISV in selective macroscopic cortical clusters. Across the cortex, the ISVs of structural connectivity and beta-band MEG functional connectivity are negatively associated with cortical myelin content indexed by the quantitative T1 relaxation rate measured by high-resolution 7 T MRI. Furthermore, MEG ISV from alpha to gamma bands relates to the hindrance and restriction of the white-matter tissue estimated by DWI microstructural models. Our findings depict the inter-relationship between the ISV of brain connectivity from multiple modalities, and highlight the role of tissue microstructure underpinning the ISV.

A whole-brain 3D myeloarchitectonic atlas: Mapping the Vogt-Vogt legacy to the cortical surface

Building precise and detailed parcellations of anatomically and functionally distinct brain areas has been a major focus in Neuroscience. Pioneer anatomists parcellated the cortical manifold based on extensive histological studies of post-mortem brain, harnessing local variations in cortical cyto- and myeloarchitecture to define areal boundaries. Compared to the cytoarchitectonic field, where multiple neuroimaging studies have recently translated this old legacy data into useful analytical resources, myeloarchitectonics, which parcellate the cortex based on the organization of myelinated fibers, has received less attention. Here, we present the neocortical surface-based myeloarchitectonic atlas based on the histology-derived maps of the Vogt-Vogt school and its 2D translation by Nieuwenhuys. In addition to a myeloarchitectonic parcellation, our package includes intracortical laminar profiles of myelin content based on Vogt-Vogt-Hopf original publications. Histology-derived myelin density mapped on our atlas demonstrated a close overlap with in vivo quantitative MRI markers for myelin and relates to cytoarchitectural features. Complementing the existing battery of approaches for digital cartography, the whole-brain myeloarchitectonic atlas offers an opportunity to validate imaging surrogate markers of myelin in both health and disease.

Empirical transmit field bias correction of T1w/T2w myelin maps

T1-weighted divided by T2-weighted (T1w/T2w) myelin maps were initially developed for neuroanatomical analyses such as identifying cortical areas, but they are increasingly used in statistical comparisons across individuals and groups with other variables of interest. Existing T1w/T2w myelin maps contain radiofrequency transmit field (B1+) biases, which may be correlated with these variables of interest, leading to potentially spurious results. Here we propose two empirical methods for correcting these transmit field biases using either explicit measures of the transmit field or alternatively a 'pseudo-transmit' approach that is highly correlated with the transmit field at 3T. We find that the resulting corrected T1w/T2w myelin maps are both better neuroanatomical measures (e.g., for use in cross-species comparisons), and more appropriate for statistical comparisons of relative T1w/T2w differences across individuals and groups (e.g., sex, age, or body-mass-index) within a consistently acquired study at 3T. We recommend that investigators who use the T1w/T2w approach for mapping cortical myelin use these B1+ transmit field corrected myelin maps going forward.

The prediction of resilience to alcohol consumption in youths: insular and subcallosal cingulate myeloarchitecture

BACKGROUND

The prediction of alcohol consumption in youths and particularly biomarkers of resilience, is critical for early intervention to reduce the risk of subsequent harmful alcohol use.

METHODS

At baseline, the longitudinal relaxation rate (R1), indexing grey matter myelination (i.e. myeloarchitecture), was assessed in 86 adolescents/young adults (mean age = 21.76, range: 15.75-26.67 years). The Alcohol Use Disorder Identification Test (AUDIT) was assessed at baseline, 1- and 2-year follow-ups (12- and 24-months post-baseline). We used a whole brain data-driven approach controlled for age, gender, impulsivity and other substance and behavioural addiction measures, such as problematic cannabis use, drug use-related problems, internet gaming, pornography use, binge eating, and levels of externalization, to predict the change in AUDIT scores from R1.

RESULTS

Greater baseline bilateral anterior insular and subcallosal cingulate R1 (cluster-corrected family-wise error p < 0.05) predict a lower risk for harmful alcohol use (measured as a reduction in AUDIT scores) at 2-year follow-up. Control analyses show that other grey matter measures (local volume or fractional anisotropy) did not reveal such an association. An atlas-based machine learning approach further confirms the findings.

CONCLUSIONS

The insula is critically involved in predictive coding of autonomic function relevant to subjective alcohol cue/craving states and risky decision-making processes. The subcallosal cingulate is an essential node underlying emotion regulation and involved in negative emotionality addiction theories. Our findings highlight insular and cingulate myeloarchitecture as a potential protective biomarker that predicts resilience to alcohol misuse in youths, providing novel identifiers for early intervention.

In vivo irreversible and reversible transverse relaxation rates in human cerebral cortex via line scans at 7 T with 250 micron resolution perpendicular to the cortical surface

Understanding how and why MR signals and their associated relaxation rates vary with cortical depth could ultimately enable the noninvasive investigation of the laminar architecture of cerebral cortex in the living human brain. However, cortical gray matter is typically only a few millimeters thick, making it challenging to sample many cortical depths with the voxel sizes commonly used in MRI studies. Line-scan techniques provide a way to overcome this challenge and here we implemented a novel line-scan GESSE pulse sequence that allowed us to measure irreversible and reversible transverse relaxation rates-R₂ and R₂´, respectively-with extremely high resolution (250 μm) in the radial direction, perpendicular to the cortical surface. Eight healthy human subjects were scanned at 7 T using this sequence, with primary visual cortex (V1) targeted in three subjects and primary motor (M1) and somatosensory cortex (S1) targeted in the other five. In all three cortical areas, a peak in R₂ values near the central depths was seen consistently across subjects-an observation that has not been made before, to our knowledge. On the other hand, no consistent pattern was apparent for R₂´ values as a function of cortical depth. The intracortical R₂ peak reported here is unlikely to be explained by myelin content or by deoxyhemoglobin in the microvasculature; however, this peak is in accord with the laminar distribution of non-heme iron in these cortical areas, known from prior histology studies. Obtaining information about tissue microstructure via measurements of transverse relaxation (and other quantitative MR contrast mechanisms) at the extremely high radial resolutions achievable through the use of line-scan techniques could therefore bring us closer to being able to perform "in vivo histology" of the cerebral cortex.

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.

A unified 3D map of microscopic architecture and MRI of the human brain

We present the first three-dimensional (3D) concordance maps of cyto- and fiber architecture of the human brain, combining histology, immunohistochemistry, and 7-T quantitative magnetic resonance imaging (MRI), in two individual specimens. These 3D maps each integrate data from approximately 800 microscopy sections per brain, showing neuronal and glial cell bodies, nerve fibers, and interneuronal populations, as well as ultrahigh-field quantitative MRI, all coaligned at the 200-μm scale to the stacked blockface images obtained during sectioning. These unprecedented 3D multimodal datasets are shared without any restrictions and provide a unique resource for the joint study of cell and fiber architecture of the brain, detailed anatomical atlasing, or modeling of the microscopic underpinnings of MRI contrasts.

3D Echo Planar Time-resolved Imaging (3D-EPTI) for ultrafast multi-parametric quantitative MRI

Multi-parametric quantitative MRI has shown great potential to improve the sensitivity and specificity of clinical diagnosis and to enhance our understanding of complex brain processes, but suffers from long scan time especially at high spatial resolution. To address this longstanding challenge, we introduce a novel approach, termed 3D Echo Planar Time-resolved Imaging (3D-EPTI), which significantly increases the acceleration capacity of MRI sampling, and provides high acquisition efficiency for multi-parametric MRI. This is achieved by exploiting the spatiotemporal correlation of MRI data at multiple timescales through new encoding strategies within and between efficient continuous readouts. Specifically, an optimized spatiotemporal CAIPI encoding within the readouts combined with a radial-block sampling strategy across the readouts enables an acceleration rate of 800 fold in the k-t space. A subspace reconstruction was employed to resolve thousands of high-quality multi-contrast images. We have demonstrated the ability of 3D-EPTI to provide robust and repeatable whole-brain simultaneous T1, T2, T2*, PD and B1+ mapping at high isotropic resolution within minutes (e.g., 1-mm isotropic resolution in 3 minutes), and to enable submillimeter multi-parametric imaging to study detailed brain structures.

Bundle myelin fraction (BMF) mapping of different white matter connections using microstructure informed tractography

To date, we have scarce information about the relative myelination level of different fiber bundles in the human brain. Indirect evidence comes from postmortem histology data but histological stainings are unable to follow a specific bundle and determine its intrinsic myelination. In this context, quantitative MRI, and diffusion MRI tractography may offer a viable solution by providing, respectively, voxel-wise myelin sensitive maps and the pathways of the major tracts of the brain. Then, "tractometry" can be used to combine these two pieces of information by averaging tissue features (obtained from any voxel-wise map) along the streamlines recovered with diffusion tractography. Although this method has been widely used in the literature, in cases of voxels containing multiple fiber populations (each with different levels of myelination), tractometry provides biased results because the same value will be attributed to any bundle passing through the voxel. To overcome this bias, we propose a new method - named "myelin streamline decomposition" (MySD) - which extends convex optimization modeling for microstructure informed tractography (COMMIT) allowing the actual value measured by a microstructural map to be deconvolved on each individual streamline, thereby recovering unique bundle-specific myelin fractions (BMFs). We demonstrate the advantage of our method with respect to tractometry in well-studied bundles and compare the cortical projection of the obtained bundle-wise myelin values of both methods. We also prove the stability of our approach across different subjects and different MRI sensitive myelin mapping approaches. This work provides a proof-of-concept of in vivo investigations of entire neuronal pathways that, to date, are not possible.

Shared genetic influences on resting‐state functional networks of the brain

The amplitude of activation in brain resting state networks (RSNs), measured with resting‐state functional magnetic resonance imaging, is heritable and genetically correlated across RSNs, indicating pleiotropy. Recent univariate genome‐wide association studies (GWASs) explored the genetic underpinnings of individual variation in RSN activity. Yet univariate genomic analyses do not describe the pleiotropic nature of RSNs. In this study, we used a novel multivariate method called genomic structural equation modeling to model latent factors that capture the shared genomic influence on RSNs and to identify single nucleotide polymorphisms (SNPs) and genes driving this pleiotropy. Using summary statistics from GWAS of 21 RSNs reported in UK Biobank (N = 31,688), the genomic latent factor analysis was first conducted in a discovery sample (N = 21,081), and then tested in an independent sample from the same cohort (N = 10,607). In the discovery sample, we show that the genetic organization of RSNs can be best explained by two distinct but correlated genetic factors that divide multimodal association networks and sensory networks. Eleven of the 17 factor loadings were replicated in the independent sample. With the multivariate GWAS, we found and replicated nine independent SNPs associated with the joint architecture of RSNs. Further, by combining the discovery and replication samples, we discovered additional SNP and gene associations with the two factors of RSN amplitude. We conclude that modeling the genetic effects on brain function in a multivariate way is a powerful approach to learn more about the biological mechanisms involved in brain function.

Charting human subcortical maturation across the adult lifespan with in vivo 7 T MRI

The human subcortex comprises hundreds of unique structures. Subcortical functioning is crucial for behavior, and disrupted function is observed in common neurodegenerative diseases. Despite their importance, human subcortical structures continue to be difficult to study in vivo. Here we provide a detailed account of 17 prominent subcortical structures and ventricles, describing their approximate iron and myelin contents, morphometry, and their age-related changes across the normal adult lifespan. The results provide compelling insights into the heterogeneity and intricate age-related alterations of these structures. They also show that the locations of many structures shift across the lifespan, which is of direct relevance for the use of standard magnetic resonance imaging atlases. The results further our understanding of subcortical morphometry and neuroimaging properties, and of normal aging processes which ultimately can improve our understanding of neurodegeneration.

Magnetic resonance fingerprinting residual signals can disassociate human grey matter regions

The importance of accurate structural discrimination of the human grey matter regions has motivated the development of observer-independent reproducible methods that account for inter-individual architectonic variations. We introduce a non-invasive statistical residual analysis framework, employing unique tissue-specific magnetic resonance fingerprinting (MRF) signals after adjusting for the effect of T₁ and T₂^* MR relaxometry parameters (here termed MRF residuals). A 7 T Siemens MR scanner was used to acquire MRF signals, quantitative transmit magnetic field (B₁⁺) maps and T₁-weighted anatomical images of eleven cortical areas (5L, 5M, 5Ci, 7A, 7P, 7PC, hIP3, BA2, BA4a, BA4p and BA6) from six female participants. MRF residual signal for each voxel was calculated as the difference between the actual and best matching MRF signal evolutions from a precomputed MRF dictionary covering a range of T₁, T₂^* and B₁⁺ values. To compare MRF residuals between regions of interest, normalised autocorrelation was used as a shape-based statistical signal characterisation method and the Euclidean distance between autocorrelation profiles of residuals was used to measure the interareal dissimilarity. In the eleven cortical areas in both cerebral hemispheres of six participants, the proposed MRF residual analysis consistently showed interareal dissimilarity profiles that concorded with histological studies, indicating that MRF residuals potentially contain tissue microstructural information. MRF residual signals provide additional area-specific information that is complementary to the MR relaxometry-based (T₁, T₂^*) information used previously for distinguishing microstructural differences between human cerebral cortex regions in vivo. The proposed approach led to more accurate identification of structural variations across cortical areas of interest.

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

Ultra-high field MRI across the depth of the cortex has the potential to provide anatomically precise biomarkers and mechanistic insights into neurodegenerative disease like Huntington's disease that show layer-selective vulnerability. Here we compare multi-parametric mapping (MPM) measures across cortical depths for a 7T 500 μm whole brain acquisition to (a) layer-specific cell measures from the von Economo histology atlas, (b) layer-specific gene expression, using the Allen Human Brain atlas and (c) white matter connections using high-fidelity diffusion tractography, at a 1.3 mm isotropic voxel resolution, from a 300mT/m Connectom MRI system. We show that R2*, but not R1, across cortical depths is highly correlated with layer-specific cell number and layer-specific gene expression. R1- and R2*-weighted connectivity strength of cortico-striatal and intra-hemispheric cortical white matter connections was highly correlated with grey matter R1 and R2* across cortical depths. Limitations of the layer-specific relationships demonstrated are at least in part related to the high cross-correlations of von Economo atlas cell counts and layer-specific gene expression across cortical layers. These findings demonstrate the potential and limitations of combining 7T MPMs, gene expression and white matter connections to provide an anatomically precise framework for tracking neurodegenerative disease.

Calibration-free pTx of the human heart at 7T via 3D universal pulses

PURPOSE

MRI at ultra-high fields in the human body is highly challenging and requires lengthy calibration times to compensate for spatially heterogeneous B 1 + profiles. This study investigates the feasibility of using pre-computed universal pulses for calibration-free homogeneous 3D flip angle distribution in the human heart at 7T.

METHODS

Twenty-two channel-wise 3D B 1 + data sets were acquired under free-breathing in 19 subjects to generate a library for an offline universal pulse (UP) design (group 1: 12 males [M] and 7 females [F], 21-66 years, 19.8-28.3 kg/m² ). Three of these subjects (2M/1F, 21-33 years, 20.8-23.6 kg/m² ) were re-scanned on different days. A 4kT-points UP optimized for the 22 channel-wise 3D B 1 + data sets in group 1 (UP22-4kT) is proposed and applied at 7T in 9 new and unseen subjects (group 2: 4M/5F, 25-56 years, 19.5-35.3 kg/m² ). Multiple tailored and universal static and dynamic parallel-transmit (pTx) pulses were designed and evaluated for different permutations of the B 1 + data sets in group 1 and 2.

RESULTS

The proposed UP22-4kT provides low B 1 + variation in all subjects, seen and unseen, without severe signal drops. Experimental data at 7T acquired with UP22-4kT shows comparable image quality as data acquired with tailored-4kT pulses and demonstrates successful calibration-free pTx of the human heart.

CONCLUSION

UP22-4kT allows for calibration-free homogeneous flip angle distributions across the human heart at 7T. Large inter-subject variations because of sex, age, and body mass index are well tolerated. The proposed universal pulse removes the need for lengthy (10-15 min) calibration scans and therefore has the potential to bring body imaging at 7T closer to the clinical application.

Minimal specifications for non-human primate MRI: Challenges in standardizing and harmonizing data collection

Recent methodological advances in MRI have enabled substantial growth in neuroimaging studies of non-human primates (NHPs), while open data-sharing through the PRIME-DE initiative has increased the availability of NHP MRI data and the need for robust multi-subject multi-center analyses. Streamlined acquisition and analysis protocols would accelerate and improve these efforts. However, consensus on minimal standards for data acquisition protocols and analysis pipelines for NHP imaging remains to be established, particularly for multi-center studies. Here, we draw parallels between NHP and human neuroimaging and provide minimal guidelines for harmonizing and standardizing data acquisition. We advocate robust translation of widely used open-access toolkits that are well established for analyzing human data. We also encourage the use of validated, automated pre-processing tools for analyzing NHP data sets. These guidelines aim to refine methodological and analytical strategies for small and large-scale NHP neuroimaging data. This will improve reproducibility of results, and accelerate the convergence between NHP and human neuroimaging strategies which will ultimately benefit fundamental and translational brain science.

QSM reconstruction challenge 2.0: A realistic in silico head phantom for MRI data simulation and evaluation of susceptibility mapping procedures

PURPOSE

To create a realistic in silico head phantom for the second QSM reconstruction challenge and for future evaluations of processing algorithms for QSM.

METHODS

We created a digital whole-head tissue property phantom by segmenting and postprocessing high-resolution (0.64 mm isotropic), multiparametric MRI data acquired at 7 T from a healthy volunteer. We simulated the steady-state magnetization at 7 T using a Bloch simulator and mimicked a Cartesian sampling scheme through Fourier-based processing. Computer code for generating the phantom and performing the MR simulation was designed to facilitate flexible modifications of the phantom in the future, such as the inclusion of pathologies as well as the simulation of a wide range of acquisition protocols. Specifically, the following parameters and effects were implemented: TR and TE, voxel size, background fields, and RF phase biases. Diffusion-weighted imaging phantom data are provided, allowing future investigations of tissue-microstructure effects in phase and QSM algorithms.

RESULTS

The brain part of the phantom featured realistic morphology with spatial variations in relaxation and susceptibility values similar to the in vivo setting. We demonstrated some of the phantom's properties, including the possibility of generating phase data with nonlinear evolution over TE due to partial-volume effects or complex distributions of frequency shifts within the voxel.

CONCLUSION

The presented phantom and computer programs are publicly available and may serve as a ground truth in future assessments of the faithfulness of quantitative susceptibility reconstruction algorithms.

Improved cortical surface reconstruction using sub-millimeter resolution MPRAGE by image denoising

Automatic cerebral cortical surface reconstruction is a useful tool for cortical anatomy quantification, analysis and visualization. Recently, the Human Connectome Project and several studies have shown the advantages of using T₁-weighted magnetic resonance (MR) images with sub-millimeter isotropic spatial resolution instead of the standard 1-mm isotropic resolution for improved accuracy of cortical surface positioning and thickness estimation. Nonetheless, sub-millimeter resolution images are noisy by nature and require averaging multiple repetitions to increase the signal-to-noise ratio for precisely delineating the cortical boundary. The prolonged acquisition time and potential motion artifacts pose significant barriers to the wide adoption of cortical surface reconstruction at sub-millimeter resolution for a broad range of neuroscientific and clinical applications. We address this challenge by evaluating the cortical surface reconstruction resulting from denoised single-repetition sub-millimeter T₁-weighted images. We systematically characterized the effects of image denoising on empirical data acquired at 0.6 mm isotropic resolution using three classical denoising methods, including denoising convolutional neural network (DnCNN), block-matching and 4-dimensional filtering (BM4D) and adaptive optimized non-local means (AONLM). The denoised single-repetition images were found to be highly similar to 6-repetition averaged images, with a low whole-brain averaged mean absolute difference of ~0.016, high whole-brain averaged peak signal-to-noise ratio of ~33.5 dB and structural similarity index of ~0.92, and minimal gray matter–white matter contrast loss (2% to 9%). The whole-brain mean absolute discrepancies in gray matter–white matter surface placement, gray matter–cerebrospinal fluid surface placement and cortical thickness estimation were lower than 165 μm, 155 μm and 145 μm—sufficiently accurate for most applications. These discrepancies were approximately one third to half of those from 1-mm isotropic resolution data. The denoising performance was equivalent to averaging ~2.5 repetitions of the data in terms of image similarity, and 1.6–2.2 repetitions in terms of the cortical surface placement accuracy. The scan-rescan variability of the cortical surface positioning and thickness estimation was lower than 170 μm. Our unique dataset and systematic characterization support the use of denoising methods for improved cortical surface reconstruction at sub-millimeter resolution.

Variable flip angle echo planar time-resolved imaging (vFA-EPTI) for fast high-resolution gradient echo myelin water imaging

Myelin water imaging techniques based on multi-compartment relaxometry have been developed as an important tool to measure myelin concentration in vivo, but are limited by the long scan time of multi-contrast multi-echo acquisition. In this work, a fast imaging technique, termed variable flip angle Echo Planar Time-Resolved Imaging (vFA-EPTI), is developed to acquire multi-echo and multi-flip-angle gradient-echo data with significantly reduced acquisition time, providing rich information for multi-compartment analysis of gradient-echo myelin water imaging (GRE-MWI). The proposed vFA-EPTI method achieved 26 folds acceleration with good accuracy by utilizing an efficient continuous readout, optimized spatiotemporal encoding across echoes and flip angles, as well as a joint subspace reconstruction. An approach to estimate off-resonance field changes between different flip-angle acquisitions was also developed to ensure high-quality joint reconstruction across flip angles. The accuracy of myelin water fraction (MWF) estimate under high acceleration was first validated by a retrospective undersampling experiment using a lengthy fully-sampled data as reference. Prospective experiments were then performed where whole-brain MWF and multi-compartment quantitative maps were obtained in 5 min at 1.5 mm isotropic resolution and 24 min at 1 mm isotropic resolution at 3T. Additionally, ultra-high resolution data at 600 μm isotropic resolution were acquired at 7T, which show detailed structures within the cortex such as the line of Gennari, demonstrating the ability of the proposed method for submillimeter GRE-MWI that can be used to study cortical myeloarchitecture in vivo.

An in-vivo study of BOLD laminar responses as a function of echo time and static magnetic field strength

Layer specific functional MRI requires high spatial resolution data. To compensate the associated poor signal to noise ratio it is common to integrate the signal from voxels at a given cortical depth. If the region is sufficiently large then physiological noise will be the dominant noise source. In this work, activation profiles in response to the same visual stimulus are compared at 1.5 T, 3 T and 7 T using a multi-echo, gradient echo (GE) FLASH sequence, with a 0.75 mm isotropic voxel size and the cortical integration approach. The results show that after integrating over a cortical volume of 40, 60 and 100 mm³ (at 7 T, 3 T, and 1.5 T, respectively), the signal is in the physiological noise dominated regime. The activation profiles obtained are similar for equivalent echo times. BOLD-like noise is found to be the dominant source of physiological noise. Consequently, the functional contrast to noise ratio is not strongly echo-time or field-strength dependent. We conclude that laminar GE-BOLD fMRI at lower field strengths is feasible but that larger patches of cortex will need to be examined, and that the acquisition efficiency is reduced.

Improving in vivo human cerebral cortical surface reconstruction using data-driven super-resolution

Accurate and automated reconstruction of the in vivo human cerebral cortical surface from anatomical magnetic resonance (MR) images facilitates the quantitative analysis of cortical structure. Anatomical MR images with sub-millimeter isotropic spatial resolution improve the accuracy of cortical surface and thickness estimation compared to the standard 1-millimeter isotropic resolution. Nonetheless, sub-millimeter resolution acquisitions require averaging multiple repetitions to achieve sufficient signal-to-noise ratio and are therefore long and potentially vulnerable to subject motion. We address this challenge by synthesizing sub-millimeter resolution images from standard 1-millimeter isotropic resolution images using a data-driven supervised machine learning-based super-resolution approach achieved via a deep convolutional neural network. We systematically characterize our approach using a large-scale simulated dataset and demonstrate its efficacy in empirical data. The super-resolution data provide improved cortical surfaces similar to those obtained from native sub-millimeter resolution data. The whole-brain mean absolute discrepancy in cortical surface positioning and thickness estimation is below 100 μm at the single-subject level and below 50 μm at the group level for the simulated data, and below 200 μm at the single-subject level and below 100 μm at the group level for the empirical data, making the accuracy of cortical surfaces derived from super-resolution sufficient for most applications.

Characterization of the cortical myeloarchitecture with inhomogeneous magnetization transfer imaging (ihMT)

BACKGROUND

Myelin specific imaging techniques to characterize white matter in demyelinating diseases such as multiple sclerosis (MS) have become an area of increasing focus. Gray matter myelination is an important marker of cortical microstructure, and its impairment is relevant in progressive MS. However, its assessment is challenging due to its thin layers. While myelin water imaging and ultra-short TE imaging have not yet been implemented to assess cortical myeloarchitecture, magnetization transfer (MT) shows promise. A recent development of the MT technique, ihMT, has demonstrated greater myelin sensitivity/specificity. Here we implemented a 3D ihMT acquisition and analysis to characterize cortical gray matter myeloarchitecture.

METHODS

20 young healthy volunteers were imaged with a 3D ihMTRAGE sequence and quantitative metrics of ihMT (ihMTsat), and dual frequency-offset MT (dual MTsat) were calculated. Cortical surface-based analysis of ihMTsat and dual MTsat were performed and compared. We also compared the cortical ihMTsat map to a cortical surface-based map of T₁-weighted images (T₁w), defined as a proxy of myelin content.

RESULTS

Cortical ihMTsat and dual MTsat maps were in qualitative agreement with previous work and the cortical T₁w map, showing higher values in primary cortices and lower values in the insula. IhMTsat and dual MTsat were significantly correlated but with important regional differences. The ratio ihMTsat/dual MTsat highlighted higher ihMTsat values in the primary cortices and sulci.

CONCLUSION

ihMTsat, a quantitative metric of ihMT, can be reliably measured in cortical gray matter and shows unique contrast between cortical regions.

Effects of MP2RAGE B1+ sensitivity on inter-site T1 reproducibility and hippocampal morphometry at 7T

Most neuroanatomical studies are based on T₁-weighted MR images, whose intensity profiles are not solely determined by the tissue's longitudinal relaxation times (T₁), but also affected by varying non-T₁ contributions, hampering data reproducibility. In contrast, quantitative imaging using the MP2RAGE sequence, for example, allows direct characterization of the brain based on the tissue property of interest. Combined with 7 Tesla (7T) MRI, this offers unique opportunities to obtain robust high-resolution brain data characterized by a high reproducibility, sensitivity and specificity. However, specific MP2RAGE parameter choices - e.g., to emphasize intracortical myelin-dependent contrast variations - can substantially impact image quality and cortical analyses through remnants of B₁⁺-related intensity variations, as illustrated in our previous work. To follow up on this: we (1) validate this protocol effect using a dataset acquired with a particularly B₁⁺ insensitive set of MP2RAGE parameters combined with parallel transmission excitation; and (2) extend our analyses to evaluate the effects on hippocampal morphometry. The latter remained unexplored initially, but can provide important insights related to generalizability and reproducibility of neurodegenerative research using 7T MRI. We confirm that B₁⁺ inhomogeneities have a considerably variable effect on cortical T₁ estimates, as well as on hippocampal morphometry depending on the MP2RAGE setup. While T₁ differed substantially across datasets initially, we show the inter-site T₁ comparability improves after correcting for the spatially varying B₁⁺ field using a separately acquired Sa2RAGE B₁⁺ map. Finally, removal of B₁⁺ residuals affects hippocampal volumetry and boundary definitions, particularly near structures characterized by strong intensity changes (e.g. cerebral spinal fluid). Taken together, we show that the choice of MP2RAGE parameters can impact T₁ comparability across sites and present evidence that hippocampal segmentation results are modulated by B₁⁺ inhomogeneities. This calls for careful (1) consideration of sequence parameters when setting acquisition protocols, as well as (2) acquisition of a B₁⁺ map to correct MP2RAGE data for potential B₁⁺ variations to allow comparison across datasets.

Probing in vivo cortical myeloarchitecture in humans via line-scan diffusion acquisitions at 7 T with 250-500 micron radial resolution

PURPOSE

The goal of this study was to measure diffusion signals within the cerebral cortex using the line-scan technique to achieve extremely high resolution in the radial direction (ie, perpendicular to the cortical surface) and to demonstrate the utility of these measurements for investigating laminar architecture in the living human brain.

METHODS

Line-scan diffusion data with 250-500 micron radial resolution were acquired at 7 T on 8 healthy volunteers, with each line prescribed perpendicularly to primary somatosensory cortex (S1) and primary motor cortex (M1). Apparent diffusion coefficients, fractional anisotropy values, and radiality indices were measured as a function of cortical depth.

RESULTS

In the deep layers of S1, we found evidence for high anisotropy and predominantly tangential diffusion, with low anisotropy observed in superficial S1. In M1, moderate anisotropy and predominantly radial diffusion was seen at almost all cortical depths. These patterns were consistent across subjects and were conspicuous without averaging data across different locations on the cortical sheet.

CONCLUSION

Our results are in accord with the myeloarchitecture of S1 and M1, known from prior histology studies: in S1, dense bands of tangential myelinated fibers run through the deep layers but not the superficial ones, and in M1, radial myelinated fibers are prominent at most cortical depths. This work therefore provides support for the idea that high-resolution diffusion signals, measured with the line-scan technique and receiving a boost in SNR at 7 T, may serve as a sensitive probe of in vivo laminar architecture.

B0-field dependence of MRI T1 relaxation in human brain

Tissue longitudinal relaxation characterized by recovery time T1 or rate R1 is a fundamental MRI contrast mechanism that is increasingly being used to study the brain's myelination patterns in both health and disease. Nevertheless, the quantitative relationship between T1 and myelination, and its dependence on B0 field strength, is still not well known. It has been theorized that in much of brain tissue, T1 field-dependence is driven by that of macromolecular protons (MP) through a mechanism called magnetization transfer (MT). Despite the explanatory power of this theory and substantial support from in-vitro experiments at low fields (<3 ​T), in-vivo evidence across clinically relevant field strengths is lacking. In this study, T1-weighted MRI was acquired in a group of eight healthy volunteers at four clinically relevant field strengths (0.55, 1.5, 3 and 7 ​T) using the same pulse sequence at a single site, and jointly analyzed based on the two-pool model of MT. MP fraction and free-water pool T1 were obtained in several brain structures at 3 and 7 ​T, which allowed distinguishing between contributions from macromolecular content and iron to tissue T1. Based on this, the T1 of MP in white matter, indirectly determined by assuming a field independent T1 of free water, was shown to increase approximately linearly with B0. This study advances our understanding of the T1 contrast mechanism and its relation to brain myelin content across the wide range of currently available MRI strengths, and it has the potential to inform design of T1 mapping methods for improved reproducibility in the human brain.

Structural Associations of Cortical Contrast and Thickness in First Episode Psychosis

There is growing evidence that psychosis is characterized by brain network abnormalities. Analyzing morphological abnormalities with T1-weighted structural MRI may be limited in discovering the extent of deviations in cortical associations. We assess whether structural associations of either cortical white-gray contrast (WGC) or cortical thickness (CT) allow for a better understanding of brain structural relationships in first episode of psychosis (FEP) patients. Principal component and structural covariance analyses were applied to WGC and CT derived from T1-weighted MRI for 116 patients and 88 controls, to explore sets of brain regions that showed group differences, and associations with symptom severity and cognitive ability in patients. We focused on 2 principal components: one encompassed primary somatomotor regions, which showed trend-like group differences in WGC, and the second included heteromodal cortices. Patients' component scores were related to general psychopathology for WGC, but not CT. Structural covariance analyses with WGC revealed group differences in pairwise correlations across widespread brain regions, mirroring areas derived from PCA. More group differences were uncovered with WGC compared with CT. WGC holds potential as a proxy measure of myelin from commonly acquired T1-weighted MRI and may be sensitive in detecting systems-level aberrations in early psychosis, and relationships with clinical/cognitive profiles.

Reduced Venous Oxygen Saturation Associates With Increased Dependence of Patients With Cerebral Autosomal Dominant Arteriopathy With Subcortical Infarcts and Leukoencephalopathy: A 7.0-T Magnetic Resonance Imaging Study

Background and Purpose- Accumulating evidence has demonstrated hemodynamic abnormalities and cerebral hypoperfusion in patients with cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL). Increased venous susceptibility assessed by susceptibility weighted imaging and mapping has been shown to indicate a decrease in venous oxygen saturation. This study aimed to investigate whether altered venous oxygen saturation is related to clinical phenotypes of CADASIL patients. Methods- Using 7.0-T susceptibility weighted imaging and mapping, we compared venous susceptibility of cortical veins between 41 CADASIL patients and 43 age- and sex-matched healthy controls. The magnetic resonance imaging lesion load, mini-mental state examination score, Barthel Index, and modified Rankin Scale were examined in the patient group, and the correlations between venous susceptibility and clinical characteristics were analyzed. Results- Venous susceptibility increased with age (r=0.508, P=0.001) and was higher in CADASIL patients than in healthy controls (t=-4.673; P<0.001). We found a positive association between venous susceptibility and the age-related white matter change scores (r=0.364; P=0.019), number of lacunar infarctions (r=0.520; P<0.001), number of cerebral microbleeds (ρ=0.445; P=0.004), and small-vessel disease scores (ρ=0.465; P=0.002) in CADASIL patients. Moreover, increased venous susceptibility was associated with higher modified Rankin Scale scores in CADASIL patients after adjustment for age- and small-vessel disease scores (odds ratio=3.178; 95% CI, 1.101-9.179; P=0.033). Conclusions- Our findings indicate that extensive cerebral hypoperfusion may induce central nervous system impairment in CADASIL, and susceptibility weighted imaging and mapping could be used clinically to assess the condition of CADASIL patients.

Biophysically motivated efficient estimation of the spatially isotropic R 2 * component from a single gradient-recalled echo measurement

Purpose

To propose and validate an efficient method, based on a biophysically motivated signal model, for removing the orientation‐dependent part of R2* using a single gradient‐recalled echo (GRE) measurement.

Methods

The proposed method utilized a temporal second‐order approximation of the hollow‐cylinder‐fiber model, in which the parameter describing the linear signal decay corresponded to the orientation‐independent part of R2*. The estimated parameters were compared to the classical, mono‐exponential decay model for R2* in a sample of an ex vivo human optic chiasm (OC). The OC was measured at 16 distinct orientations relative to the external magnetic field using GRE at 7T. To show that the proposed signal model can remove the orientation dependence of R2*, it was compared to the established phenomenological method for separating R2* into orientation‐dependent and ‐independent parts.

Results

Using the phenomenological method on the classical signal model, the well‐known separation of R2* into orientation‐dependent and ‐independent parts was verified. For the proposed model, no significant orientation dependence in the linear signal decay parameter was observed.

Conclusions

Since the proposed second‐order model features orientation‐dependent and ‐independent components at distinct temporal orders, it can be used to remove the orientation dependence of R2* using only a single GRE measurement.

A comparison of in vivo MRI based cortical myelin mapping using T1w/T2w and R1 mapping at 3T

In this manuscript, we compare two commonly used methods to perform cortical mapping based on myelination of the human neocortex. T1w/T2w and R1 maps with matched total acquisition times were obtained from a young cohort in randomized order and using a test–retest design. Both methodologies showed cortical myelin maps that enhanced similar anatomical features, namely primary sensory regions known to be myelin rich. T1w/T2w maps showed increased robustness to movement artifacts in comparison to R1 maps, while the test re-test reproducibility of both methods was comparable. Based on Brodmann parcellation, both methods showed comparable variability within each region. Having parcellated cortical myelin maps into VDG11b areas of 4a, 4p, 3a, 3b, 1, 2, V2, and MT, both methods behave identically with R1 showing an increased variability between subjects. In combination with the test re-test evaluation, we concluded that this increased variability between subjects reflects relevant tissue variability. A high level of correlation was found between the R1 and T1w/T2w regions with regions of higher deviations being co-localized with those where the transmit RF field deviated most from its nominal value. We conclude that R1 mapping strategies might be preferable when studying different population cohorts where cortical properties are expected to be altered while T1w/T2w mapping will have advantages when performing cortical based segmentation.

hMRI - A toolbox for quantitative MRI in neuroscience and clinical research

Neuroscience and clinical researchers are increasingly interested in quantitative magnetic resonance imaging (qMRI) due to its sensitivity to micro-structural properties of brain tissue such as axon, myelin, iron and water concentration. We introduce the hMRI-toolbox, an open-source, easy-to-use tool available on GitHub, for qMRI data handling and processing, presented together with a tutorial and example dataset. This toolbox allows the estimation of high-quality multi-parameter qMRI maps (longitudinal and effective transverse relaxation rates R1 and R2⋆, proton density PD and magnetisation transfer MT saturation) that can be used for quantitative parameter analysis and accurate delineation of subcortical brain structures. The qMRI maps generated by the toolbox are key input parameters for biophysical models designed to estimate tissue microstructure properties such as the MR g-ratio and to derive standard and novel MRI biomarkers. Thus, the current version of the toolbox is a first step towards in vivo histology using MRI (hMRI) and is being extended further in this direction. Embedded in the Statistical Parametric Mapping (SPM) framework, it benefits from the extensive range of established SPM tools for high-accuracy spatial registration and statistical inferences and can be readily combined with existing SPM toolboxes for estimating diffusion MRI parameter maps. From a user's perspective, the hMRI-toolbox is an efficient, robust and simple framework for investigating qMRI data in neuroscience and clinical research.

The myeloarchitecture of impulsivity: premature responding in youth is associated with decreased myelination of ventral putamen

Impulsivity has been suggested as a neurocognitive endophenotype conferring risk across a number of neuropsychiatric conditions, including substance and behavioural addictions, eating disorders, and attention deficit/hyperactivity disorder. We used a paradigm with interspecies translation validity (the four-choice serial reaction time task, 4CSRTT) to assess 'waiting' impulsivity in a youth sample (N = 99, aged 16-26 years). We collected magnetization prepared two rapid acquisition gradient echo (MP2RAGE) scans, which enabled us to measure R1, the longitudinal relaxation rate, a parameter closely related to tissue myelin content, as well as quantify grey matter volume. We also assessed inhibitory control (commission errors) on a Go/NoGo task and measured decisional impulsivity (delay discounting) using the Monetary Choice Questionnaire (MCQ). We found R1 of the bilateral ventral putamen was negatively correlated with premature responding, the index of waiting impulsivity on the 4CSRTT. Heightened impulsivity in youth was significantly and specifically associated with lower levels of myelination in the ventral putamen. Impulsivity was not associated with grey matter volume. The association with myelination was specific to waiting impulsivity: R1 was not associated with decisional impulsivity on the MCQ or inhibitory control on the Go/NoGo task. We report that heightened waiting impulsivity, measured as premature responding on the 4CSRTT, is specifically associated with lower levels of ventral putaminal myelination, measured using R1. This may represent a neural signature of vulnerability to diseases associated with excessive impulsivity and demonstrates the added explanatory power of quantifying the mesoscopic organization of the human brain, over and above macroscopic volumetric measurements.

Laminar signal extraction over extended cortical areas by means of a spatial GLM

There is converging evidence that distinct neuronal processes leave distinguishable footprints in the laminar BOLD response. However, even though the achievable spatial resolution in functional MRI has much improved over the years, it is still challenging to separate signals arising from different cortical layers. In this work, we propose a new method to extract laminar signals. We use a spatial General Linear Model in combination with the equivolume principle of cortical layers to unmix laminar signals instead of interpolating through and integrating over a cortical area: thus reducing partial volume effects. Not only do we provide a mathematical framework for extracting laminar signals with a spatial GLM, we also illustrate that the best case scenarios of existing methods can be seen as special cases within the same framework. By means of simulation, we show that this approach has a sharper point spread function, providing better signal localisation. We further assess the partial volume contamination in cortical profiles from high resolution human ex vivo and in vivo structural data, and provide a full account of the benefits and potential caveats. We eschew here any attempt to validate the spatial GLM on the basis of fMRI data as a generally accepted ground-truth pattern of laminar activation does not currently exist. This approach is flexible in terms of the number of layers and their respective thickness, and naturally integrates spatial regularisation along the cortex, while preserving laminar specificity. Care must be taken, however, as this procedure of unmixing is susceptible to sources of noise in the data or inaccuracies in the laminar segmentation.

Age-related deficits in intracortical myelination in young adults with bipolar disorder type I

BACKGROUND

Previous studies have implicated white-matter-related changes in the pathophysiology of bipolar disorder. However, most of what is known is derived from in vivo subcortical white-matter imaging or postmortem studies. In this study, we investigated whole-brain intracortical myelin (ICM) content in people with bipolar disorder type I and controls.

METHODS

Between Sept. 1, 2014, and Jan. 31, 2017, we used a 3 T General Electric scanner to collect T1-weighted images in 45 people with bipolar disorder type I and 60 controls aged 17 to 45 years using an optimized sequence that was sensitive to ICM content. We analyzed images using a surfacebased approach. We used general linear models with quadratic age terms to examine the signal trajectory of ICM across the age range.

RESULTS

In healthy controls, the T1-weighted signal followed an inverted-U trajectory over age; in people with bipolar disorder type I, the association between ICM and age followed a flat trajectory (p < 0.05, Bonferroni corrected). Exploratory analyses showed that ICM signal intensity was associated with duration of illness, age of onset, and anticonvulsant and antipsychotic use in people with bipolar disorder type I (p < 0.05, uncorrected).

LIMITATIONS

Because of the cross-sectional nature of the study, we were unable to comment on whether the effects were due to dysmyelination or demyelination in bipolar disorder.

CONCLUSION

This foundational study is, to our knowledge, the first to show global age-related deficits in ICM maturation throughout the cortex in bipolar disorder. Considering the impact of myelination on the maintenance of neural synchrony and the integrity of neural connections, this work may help us better understand the cognitive and behavioural deficits seen in bipolar disorder.

Brain tissue segmentation based on MP2RAGE multi-contrast images in 7 T MRI

We proposed a method for segmentation of brain tissues-gray matter, white matter, and cerebrospinal fluid-using multi-contrast images, including a T1 map and a uniform T1-weighted image, from a magnetization-prepared 2 rapid acquisition gradient echoes (MP2RAGE) sequence at 7 Tesla. The proposed method was evaluated with respect to the processing time and the similarity of the segmented masks of brain tissues with those obtained using FSL, FreeSurfer, and SPM12. The processing time of the proposed method (28 ± 0 s) was significantly shorter than those of FSL and SPM12 (444 ± 4 s and 159 ± 2 s for FSL and SPM12, respectively). In the similarity assessment, the tissue mask of the brain obtained by the proposed method showed higher consistency with those obtained using FSL than with those obtained using SPM12. The proposed method misclassified the subcortical structures and large vessels since it is based on the intensities of multi-contrast images obtained using MP2RAGE, which uses a similar segmentation approach as FSL but is not based on a template image or a parcellated brain atlas, which are used for FreeSurfer and SPM12, respectively. However, the proposed method showed good segmentation in the cerebellum and white matter in the medial part of the brain in comparison with the other methods. Thus, because the proposed method using different contrast images of MP2RAGE sequence showed the shortest processing time and similar segmentation ability as the other methods, it may be useful for both neuroimaging research and clinical diagnosis.

The effect of motion correction interpolation on quantitative T1 mapping with MRI

Quantitative magnetic resonance imaging (qMRI) is a technique for mapping the physical properties of the underlying tissue using several MR images with different contrasts. To overcome subject motion between the acquired images, it is necessary to register the images to a common reference frame. A drawback of registration is the use of interpolation and resampling techniques, which can introduce artifacts into the interpolated data. These artifacts could have unfavorable effects on the accuracy of the estimated tissue's physical properties. Here, we quantified the error of interpolation and resampling on T1-weighted images and studied its effects on the mapping of the longitudinal relaxation time (T1) using variable flip angles. We simulated T1-weighted images and calculated the transformation error resulting from interpolation and resampling. We found that the error is a function of the image contrast (i.e., flip angle) and of the translation and rotation of the image. Furthermore, we found that the error in the T1-weighted images has a substantial effect on the T1 estimation, of the order of 10% of the signal in the brain's gray and white matter. Hence, minimizing the registration error can enable more accurate in vivo modeling of brain microstructure.