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.

Inhibitory TMS over Visual Area V5/MT Disrupts Visual Speech Recognition

During face-to-face communication, the perception and recognition of facial movements can facilitate individuals' understanding of what is said. Facial movements are a form of complex biological motion. Separate neural pathways are thought to processing (1) simple, nonbiological motion with an obligatory waypoint in the motion-sensitive visual middle temporal area (V5/MT); and (2) complex biological motion. Here, we present findings that challenge this dichotomy. Neuronavigated offline transcranial magnetic stimulation (TMS) over V5/MT on 24 participants (17 females and 7 males) led to increased response times in the recognition of simple, nonbiological motion as well as visual speech recognition compared with TMS over the vertex, an active control region. TMS of area V5/MT also reduced practice effects on response times, that are typically observed in both visual speech and motion recognition tasks over time. Our findings provide the first indication that area V5/MT causally influences the recognition of visual speech.SIGNIFICANCE STATEMENT In everyday face-to-face communication, speech comprehension is often facilitated by viewing a speaker's facial movements. Several brain areas contribute to the recognition of visual speech. One area of interest is the motion-sensitive visual medial temporal area (V5/MT), which has been associated with the perception of simple, nonbiological motion such as moving dots, as well as more complex, biological motion such as visual speech. Here, we demonstrate using noninvasive brain stimulation that area V5/MT is causally relevant in recognizing visual speech. This finding provides new insights into the neural mechanisms that support the perception of human communication signals, which will help guide future research in typically developed individuals and populations with communication difficulties.

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.

Identification of Laminar Composition in Cerebral Cortex Using Low-Resolution Magnetic Resonance Images and Trust Region Optimization Algorithm

Pathological changes in the cortical lamina can cause several mental disorders. Visualization of these changes in vivo would enhance their diagnostics. Recently a framework for visualizing cortical structures by magnetic resonance imaging (MRI) has emerged. This is based on mathematical modeling of multi-component T1 relaxation at the sub-voxel level. This work proposes a new approach for their estimation. The approach is validated using simulated data. Sixteen MRI experiments were carried out on healthy volunteers. A modified echo-planar imaging (EPI) sequence was used to acquire 105 individual volumes. Data simulating the images were created, serving as the ground truth. The model was fitted to the data using a modified Trust Region algorithm. In single voxel experiments, the estimation accuracy of the T1 relaxation times depended on the number of optimization starting points and the level of noise. A single starting point resulted in a mean percentage error (MPE) of 6.1%, while 100 starting points resulted in a perfect fit. The MPE was <5% for the signal-to-noise ratio (SNR) ≥ 38 dB. Concerning multiple voxel experiments, the MPE was <5% for all components. Estimation of T1 relaxation times can be achieved using the modified algorithm with MPE < 5%.

The Myelin Content of the Human Precentral Hand Knob Reflects Interindividual Differences in Manual Motor Control at the Physiological and Behavioral Level

The primary motor cortex hand area (M1_HAND) and adjacent dorsal premotor cortex (PMd) form the so-called motor hand knob in the precentral gyrus. M1_HAND and PMd are critical for dexterous hand use and are densely interconnected via corticocortical axons, lacking a sharp demarcating border. In 24 young right-handed volunteers, we performed multimodal mapping to delineate the relationship between structure and function in the right motor hand knob. Quantitative structural magnetic resonance imaging (MRI) at 3 tesla yielded regional R1 maps as a proxy of cortical myelin content. Participants also underwent functional MRI (fMRI). We mapped task-related activation and temporal precision, while they performed a visuomotor synchronization task requiring visually cued abduction movements with the left index or little finger. We also performed sulcus-aligned transcranial magnetic stimulation of the motor hand knob to localize the optimal site (hotspot) for evoking a motor evoked potential (MEP) in two intrinsic hand muscles. Individual motor hotspot locations varied along the rostrocaudal axis. The more rostral the motor hotspot location in the precentral crown, the longer were corticomotor MEP latencies. “Hotspot rostrality” was associated with the regional myelin content in the precentral hand knob. Cortical myelin content also correlated positively with task-related activation of the precentral crown and temporal precision during the visuomotor synchronization task. Together, our results suggest a link among cortical myelination, the spatial cortical representation, and temporal precision of finger movements. We hypothesize that the myelination of cortical axons facilitates neuronal integration in PMd and M1_HAND and, hereby, promotes the precise timing of movements.

SIGNIFICANCE STATEMENT Here we used magnetic resonance imaging and transcranial magnetic stimulation of the precentral motor hand knob to test for a link among cortical myelin content, functional corticomotor representations, and manual motor control. A higher myelin content of the precentral motor hand knob was associated with more rostral corticomotor presentations, with stronger task-related activation and a higher precision of movement timing during a visuomotor synchronization task. We propose that a high precentral myelin content enables fast and precise neuronal integration in M1 (primary motor cortex) and dorsal premotor cortex, resulting in higher temporal precision during dexterous hand use. Our results identify the degree of myelination as an important structural feature of the neocortex that is tightly linked to the function and behavior supported by the cortical area.

Quantitative T1 mapping using multi-slice multi-shot inversion recovery EPI

An efficient multi-slice inversion–recovery EPI (MS-IR-EPI) sequence for fast, high spatial resolution, quantitative T₁ mapping is presented, using a segmented simultaneous multi-slice acquisition, combined with slice order shifting across multiple acquisitions. The segmented acquisition minimises the effective TE and readout duration compared to a single-shot EPI scheme, reducing geometric distortions to provide high quality T₁ maps with a narrow point-spread function. The precision and repeatability of MS-IR-EPI T₁ measurements are assessed using both T₁-calibrated and T₂-calibrated ISMRM/NIST phantom spheres at 3 and 7 T and compared with single slice IR and MP2RAGE methods. Magnetization transfer (MT) effects of the spectrally-selective fat-suppression (FS) pulses required for in vivo imaging are shown to shorten the measured in-vivo T₁ values. We model the effect of these fat suppression pulses on T₁ measurements and show that the model can remove their MT contribution from the measured T₁, thus providing accurate T₁ quantification. High spatial resolution T₁ maps of the human brain generated with MS-IR-EPI at 7 T are compared with those generated with the widely implemented MP2RAGE sequence. Our MS-IR-EPI sequence provides high SNR per unit time and sharper T₁ maps than MP2RAGE, demonstrating the potential for ultra-high resolution T₁ mapping and the improved discrimination of functionally relevant cortical areas in the human brain.

Cross-species neuroscience: closing the explanatory gap

Neuroscience has seen substantial development in non-invasive methods available for investigating the living human brain. However, these tools are limited to coarse macroscopic measures of neural activity that aggregate the diverse responses of thousands of cells. To access neural activity at the cellular and circuit level, researchers instead rely on invasive recordings in animals. Recent advances in invasive methods now permit large-scale recording and circuit-level manipulations with exquisite spatio-temporal precision. Yet, there has been limited progress in relating these microcircuit measures to complex cognition and behaviour observed in humans. Contemporary neuroscience thus faces an explanatory gap between macroscopic descriptions of the human brain and microscopic descriptions in animal models. To close the explanatory gap, we propose adopting a cross-species approach. Despite dramatic differences in the size of mammalian brains, this approach is broadly justified by preserved homology. Here, we outline a three-armed approach for effective cross-species investigation that highlights the need to translate different measures of neural activity into a common space. We discuss how a cross-species approach has the potential to transform basic neuroscience while also benefiting neuropsychiatric drug development where clinical translation has, to date, seen minimal success. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

Quantitative third-harmonic generation imaging of mouse visual cortex areas reveals correlations between functional maps and structural substrates

The structure of brain regions is assumed to correlate with their function, but there are very few instances in which the relationship has been demonstrated in the live brain. This is due to the difficulty of simultaneously measuring functional and structural properties of brain areas, particularly at cellular resolution. Here, we performed label-free, third-harmonic generation (THG) microscopy to obtain a key structural signature of cortical areas, their effective attenuation lengths (EAL), in the vertical columns of functionally defined primary visual cortex and five adjacent visual areas in awake mice. EALs measured by THG microscopy in the cortex and white matter showed remarkable correspondence with the functional retinotopic sign map of each area. Structural features such as cytoarchitecture, myeloarchitecture and blood vessel architecture were correlated with areal EAL values, suggesting that EAL is a function of these structural features as an optical property of these areas. These results demonstrate for the first time a strong relationship between structural substrates of visual cortical areas and their functional representation maps in vivo. This study may also help in understanding the coupling between structure and function in other animal models as well as in humans.

Functional Segregation of the Middle Temporal Visual Motion Area Revealed With Coactivation-Based Parcellation

Traditionally, the visual motion area (MT) is considered as a brain region specialized for visual motion perception. However, accumulating evidence showed that MT is also related to various functions, suggesting that it is a complex functional area and different functional subregions might exist in this area. To delineate functional subregions of this area, left and right masks of MT were defined using meta-analysis in the BrainMap database, and coactivation-based parcellation was then performed on these two masks. Two dorsal subregions (Cl1 and Cl2) and one ventral subregion (Cl3) of left MT, as well as two dorsal-anterior subregions (Cl1 and Cl2), one ventral-anterior subregion (Cl3), and an additional posterior subregion (Cl4) of right MT were identified. In addition to vision motion, distinct and specific functions were identified in different subregions characterized by task-dependent functional connectivity mapping and forward/reverse inference on associated functions. These results not only were in accordance with the previous findings of a hemispheric asymmetry of MT, but also strongly confirmed the existence of subregions in this region with distinct and specific functions. Furthermore, our results extend the special role of visual motion perception on this area and might facilitate future cognitive study.

In vivo Probabilistic Structural Atlas of the Inferior and Superior Colliculi, Medial and Lateral Geniculate Nuclei and Superior Olivary Complex in Humans Based on 7 Tesla MRI

Despite extensive neuroimaging research of primary sensory cortices involved in auditory and visual functions, subcortical structures within these domains, such as the inferior and superior colliculi, the medial and lateral geniculate nuclei and the superior olivary complex, are currently understudied with magnetic resonance imaging (MRI) in living humans. This is because a precise localization of these nuclei is hampered by the limited contrast and sensitivity of conventional neuroimaging methods for deep brain nuclei. In this work, we used 7 Tesla multi-modal (T₂-weighted and diffusion fractional anisotropy) 1.1 mm isotropic resolution MRI to achieve high sensitivity and contrast for single-subject brainstem and thalamic nuclei delineation. After precise coregistration to stereotactic space, we generated an in vivo human probabilistic atlas of auditory (medial geniculate nucleus, inferior colliculus, and superior olivary complex) and visual (lateral geniculate nucleus and superior colliculus) subcortical nuclei. We foresee the use of this atlas as a tool to precisely identify the location and shape of auditory/visual deep nuclei in research as well as clinical human studies.

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.

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.

Parcellating Cerebral Cortex: How Invasive Animal Studies Inform Noninvasive Mapmaking in Humans

The cerebral cortex in mammals contains a mosaic of cortical areas that differ in function, architecture, connectivity, and/or topographic organization. A combination of local connectivity (within-area microcircuitry) and long-distance (between-area) connectivity enables each area to perform a unique set of computations. Some areas also have characteristic within-area mesoscale organization, reflecting specialized representations of distinct types of information. Cortical areas interact with one another to form functional networks that mediate behavior, and each area may be a part of multiple, partially overlapping networks. Given their importance to the understanding of brain organization, mapping cortical areas across species is a major objective of systems neuroscience and has been a century-long challenge. Here, we review recent progress in multi-modal mapping of mouse and nonhuman primate cortex, mainly using invasive experimental methods. These studies also provide a neuroanatomical foundation for mapping human cerebral cortex using noninvasive neuroimaging, including a new map of human cortical areas that we generated using a semiautomated analysis of high-quality, multimodal neuroimaging data. We contrast our semiautomated approach to human multimodal cortical mapping with various extant fully automated human brain parcellations that are based on only a single imaging modality and offer suggestions on how to best advance the noninvasive brain parcellation field. We discuss the limitations as well as the strengths of current noninvasive methods of mapping brain function, architecture, connectivity, and topography and of current approaches to mapping the brain's functional networks.

Using diffusion MRI to discriminate areas of cortical grey matter

Cortical area parcellation is a challenging problem that is often approached by combining structural imaging (e.g., quantitative T1, diffusion-based connectivity) with functional imaging (e.g., task activations, topological mapping, resting state correlations). Diffusion MRI (dMRI) has been widely adopted to analyse white matter microstructure, but scarcely used to distinguish grey matter regions because of the reduced anisotropy there. Nevertheless, differences in the texture of the cortical 'fabric' have long been mapped by histologists to distinguish cortical areas. Reliable area-specific contrast in the dMRI signal has previously been demonstrated in selected occipital and sensorimotor areas. We expand upon these findings by testing several diffusion-based feature sets in a series of classification tasks. Using Human Connectome Project (HCP) 3T datasets and a supervised learning approach, we demonstrate that diffusion MRI is sensitive to architectonic differences between a large number of different cortical areas defined in the HCP parcellation. By employing a surface-based cortical imaging pipeline, which defines diffusion features relative to local cortical surface orientation, we show that we can differentiate areas from their neighbours with higher accuracy than when using only fractional anisotropy or mean diffusivity. The results suggest that grey matter diffusion may provide a new, independent source of information for dividing up the cortex.

•

Diffusion MRI provides architectonic contrast between a variety of cortical areas.

•

dMRI provides complementary information in areas of low myelin content.

•

High order rotational invariants provide more reliable parcellation results.

High resolution data analysis strategies for mesoscale human functional MRI at 7 and 9.4T

The advent of ultra-high field functional magnetic resonance imaging (fMRI) has greatly facilitated submillimeter resolution acquisitions (voxel volume below (1 mm³)), allowing the investigation of cortical columns and cortical depth dependent (i.e. laminar) structures in the human brain. Advanced data analysis techniques are essential to exploit the information in high resolution functional measures. In this article, we use recent, exemplary 9.4 T human functional and anatomical data to review the advantages and disadvantages of (1) pooling high resolution data across regions of interest for cortical depth profile analysis, (2) pooling across cortical depths for mapping patches of cortex while discarding depth-dependent (i.e. columnar) effects, and (3) isotropic sampling without pooling to assess individual voxel’s responses. A set of cortical depth meshes may be a solution to sampling information tangentially while keeping correspondence across depths. For quantitative analysis of the spatial organization in fine-grained structures, a cortical grid approach is advantageous. We further extend this general framework by combining it with a previously introduced cortical layer volume-preserving (equi-volume) approach. This framework can readily accommodate the research questions which allow for spatial smoothing within or across layers. We demonstrate and discuss that equi-volume sampling yields a slight advantage over equidistant sampling given the current limitations of fMRI voxel size, participant motion, coregistration and segmentation. Our 9.4 T human anatomical and functional data indicate the advantage over lower fields including 7 T and demonstrate the practical applicability of T₂^* and T₂-weighted fMRI acquisitions.

•

High resolution regular cortical grids are advantageous for local applications.

•

Equi-volume sampling is slightly advantageous over equidistant sampling in-vivo.

•

Isotropic submillimeter cortical sampling without spatial pooling requires high SNR.

•

9.4 T human T2 and T2* BOLD fMRI are practically feasible and provide high SNR.

•

9.4 T T2*-weighted 0.35 mm iso. res. anatomical images for laminar contrast in vivo.

The impact of ultra-high field MRI on cognitive and computational neuroimaging

The ability to measure functional brain responses non-invasively with ultra high field MRI (7 T and above) represents a unique opportunity in advancing our understanding of the human brain. Compared to lower fields (3 T and below), ultra high field MRI has an increased sensitivity, which can be used to acquire functional images with greater spatial resolution, and greater specificity of the blood oxygen level dependent (BOLD) signal to the underlying neuronal responses. Together, increased resolution and specificity enable investigating brain functions at a submillimeter scale, which so far could only be done with invasive techniques. At this mesoscopic spatial scale, perception, cognition and behavior can be probed at the level of fundamental units of neural computations, such as cortical columns, cortical layers, and subcortical nuclei. This represents a unique and distinctive advantage that differentiates ultra high from lower field imaging and that can foster a tighter link between fMRI and computational modeling of neural networks. So far, functional brain mapping at submillimeter scale has focused on the processing of sensory information and on well-known systems for which extensive information is available from invasive recordings in animals. It remains an open challenge to extend this methodology to uniquely human functions and, more generally, to systems for which animal models may be problematic. To succeed, the possibility to acquire high-resolution functional data with large spatial coverage, the availability of computational models of neural processing as well as accurate biophysical modeling of neurovascular coupling at mesoscopic scale all appear necessary.

In vivo evidence of functional and anatomical stripe-based subdivisions in human V2 and V3

Visual cortex contains a hierarchy of visual areas. The earliest cortical area (V1) contains neurons responding to colour, form and motion. Later areas specialize on processing of specific features. The second visual area (V2) in non-human primates contains a stripe-based anatomical organization, initially defined using cytochrome-oxidase staining of post-mortem tissue. Neurons in these stripes have been proposed to serve distinct functional specializations, e.g. processing of color, form and motion. These stripes represent an intermediate stage in visual hierarchy and serve a key role in the increasing functional specialization of visual areas. Using sub-millimeter high-field functional and anatomical MRI (7T), we provide in vivo evidence for stripe-based subdivisions in humans. Using functional MRI, we contrasted responses elicited by stimuli alternating at slow and fast temporal frequencies. We revealed stripe-based subdivisions in V2 ending at the V1/V2 border. The human stripes reach into V3. Using anatomical MRI optimized for myelin contrast within gray matter, we also observe a stripe pattern. Stripe subdivisions preferentially responding to fast temporal frequencies are more myelinated. As such, functional and anatomical measures provide independent and converging evidence for functional organization into striped-based subdivisions in human V2 and V3.

Relationships between cortical myeloarchitecture and electrophysiological networks

The human brain relies upon the dynamic formation and dissolution of a hierarchy of functional networks to support ongoing cognition. However, how functional connectivities underlying such networks are supported by cortical microstructure remains poorly understood. Recent animal work has demonstrated that electrical activity promotes myelination. Inspired by this, we test a hypothesis that gray-matter myelin is related to electrophysiological connectivity. Using ultra-high field MRI and the principle of structural covariance, we derive a structural network showing how myelin density differs across cortical regions and how separate regions can exhibit similar myeloarchitecture. Building upon recent evidence that neural oscillations mediate connectivity, we use magnetoencephalography to elucidate networks that represent the major electrophysiological pathways of communication in the brain. Finally, we show that a significant relationship exists between our functional and structural networks; this relationship differs as a function of neural oscillatory frequency and becomes stronger when integrating oscillations over frequency bands. Our study sheds light on the way in which cortical microstructure supports functional networks. Further, it paves the way for future investigations of the gray-matter structure/function relationship and its breakdown in pathology.

Histological Basis of Laminar MRI Patterns in High Resolution Images of Fixed Human Auditory Cortex

Functional magnetic resonance imaging (fMRI) studies of the auditory region of the temporal lobe would benefit from the availability of image contrast that allowed direct identification of the primary auditory cortex, as this region cannot be accurately located using gyral landmarks alone. Previous work has suggested that the primary area can be identified in magnetic resonance (MR) images because of its relatively high myelin content. However, MR images are also affected by the iron content of the tissue and in this study we sought to confirm that different MR image contrasts did correlate with the myelin content in the gray matter and were not primarily affected by iron content as is the case in the primary visual and somatosensory areas. By imaging blocks of fixed post-mortem cortex in a 7 T scanner and then sectioning them for histological staining we sought to assess the relative contribution of myelin and iron to the gray matter contrast in the auditory region. Evaluating the image contrast in T2*-weighted images and quantitative R2* maps showed a reasonably high correlation between the myelin density of the gray matter and the intensity of the MR images. The correlation with T₁-weighted phase sensitive inversion recovery (PSIR) images was better than with the previous two image types, and there were clearly differentiated borders between adjacent cortical areas in these images. A significant amount of iron was present in the auditory region, but did not seem to contribute to the laminar pattern of the cortical gray matter in MR images. Similar levels of iron were present in the gray and white matter and although iron was present in fibers within the gray matter, these fibers were fairly uniformly distributed across the cortex. Thus, we conclude that T₁- and T2*-weighted imaging sequences do demonstrate the relatively high myelin levels that are characteristic of the deep layers in primary auditory cortex and allow it and some of the surrounding areas to be reliably distinguished.

Recent applications of UHF-MRI in the study of human brain function and structure: a review

The increased availability of ultra-high-field (UHF) MRI has led to its application in a wide range of neuroimaging studies, which are showing promise in transforming fundamental approaches to human neuroscience. This review presents recent work on structural and functional brain imaging, at 7 T and higher field strengths. After a short outline of the effects of high field strength on MR images, the rapidly expanding literature on UHF applications of blood-oxygenation-level-dependent-based functional MRI is reviewed. Structural imaging is then discussed, divided into sections on imaging weighted by relaxation time, including quantitative relaxation time mapping, phase imaging and quantitative susceptibility mapping, angiography, diffusion-weighted imaging, and finally magnetization-transfer imaging. The final section discusses studies using the high spatial resolution available at UHF to identify explicit links between structure and function. Copyright © 2015 John Wiley & Sons, Ltd.

A method for the dynamic correction of B0-related distortions in single-echo EPI at 7T

We propose a method to calculate field maps from the phase of each EPI in an fMRI time series. These field maps can be used to correct the corresponding magnitude images for distortion caused by inhomogeneity in the static magnetic field. In contrast to conventional static distortion correction, in which one ‘snapshot’ field map is applied to all subsequent fMRI time points, our method also captures dynamic changes to B₀ which arise due to motion and respiration. The approach is based on the assumption that the non-B₀-related contribution to the phase measured by each radio-frequency coil, which is dominated by the coil sensitivity, is stable over time and can therefore be removed to yield a field map from EPI.

Our solution addresses imaging with multi-channel coils at ultra-high field (7 T), where phase offsets vary rapidly in space, phase processing is non-trivial and distortions are comparatively large. We propose using dual-echo gradient echo reference scan for the phase offset calculation, which yields estimates with high signal-to-noise ratio. An extrapolation method is proposed which yields reliable estimates for phase offsets even where motion is large and a tailored phase unwrapping procedure for EPI is suggested which gives robust results in regions with disconnected tissue or strong signal decay.

Phase offsets are shown to be stable during long measurements (40 min) and for large head motions. The dynamic distortion correction proposed here is found to work accurately in the presence of large motion (up to 8.1°), whereas a conventional method based on single field map fails to correct or even introduces distortions (up to 11.2 mm). Finally, we show that dynamic unwarping increases the temporal stability of EPI in the presence of motion.

Our approach can be applied to any EPI measurements without the need for sequence modification.

Individual Differences in the Alignment of Structural and Functional Markers of the V5/MT Complex in Primates

Extrastriate visual area V5/MT in primates is defined both structurally by myeloarchitecture and functionally by distinct responses to visual motion. Myelination is directly identifiable from postmortem histology but also indirectly by image contrast with structural magnetic resonance imaging (sMRI). First, we compared the identification of V5/MT using both sMRI and histology in Rhesus macaques. A section-by-section comparison of histological slices with in vivo and postmortem sMRI for the same block of cortical tissue showed precise correspondence in localizing heavy myelination for V5/MT and neighboring MST. Thus, sMRI in macaques accurately locates histologically defined myelin within areas known to be motion selective. Second, we investigated the functionally homologous human motion complex (hMT+) using high-resolution in vivo imaging. Humans showed considerable intersubject variability in hMT+ location, when defined with myelin-weighted sMRI signals to reveal structure. When comparing sMRI markers to functional MRI in response to moving stimuli, a region of high myelin signal was generally located within the hMT+ complex. However, there were considerable differences in the alignment of structural and functional markers between individuals. Our results suggest that variation in area identification for hMT+ based on structural and functional markers reflects individual differences in human regional brain architecture.

Lines of Baillarger in vivo and ex vivo: Myelin contrast across lamina at 7T MRI and histology

The human cerebral cortex is characterized by a number of features that are not uniformly distributed, such as the presence of multiple cytoarchitectonic elements and of myelinated layers running tangentially to the cortex surface. The presence and absence of these features are the basis of the parcellation of the cerebral cortex in several areas. A number of areas show myelin increases localized within the cortex, e.g., the stria of Gennari located in layer IV of the primary visual cortex. Sub-millimeter MRI can resolve myelin variations across the human cortex and may allow in vivo parcellation of these brain areas. Here, we image within-area myelination. We modified a T1-weighted (T1-w) MPRAGE sequence to enhance myelin visualization within the cortex. First, we acquired images from an ex vivo sample, and compared MRI laminar profiles from calcarine (corresponding to primary visual cortex) and extra-calcarine areas with histology sections from the same locations. Laminar profiles between myelin stained sections and the T1-w images were similar both in calcarine as well as extra-calcarine cortex. In calcarine cortex, the profile reveals the stria of Gennari. In extra-calcarine cortex, a similar profile exists which we suggest corresponds to the lines of Baillarger. Next, we adapted the same sequence to image within-area myelination in vivo. Also in in vivo data, we discriminated similar laminar profiles in calcarine and extra-calcarine cortex, extending into parietal and frontal lobes. We argue that this myelin pattern outside the calcarine cortex represents the lines of Baillarger.

Functional and anatomical properties of human visual cortical fields

Human visual cortical fields (VCFs) vary in size and anatomical location across individual subjects. Here, we used functional magnetic resonance imaging (fMRI) with retinotopic stimulation to identify VCFs on the cortical surface. We found that aligning and averaging VCF activations across the two hemispheres provided clear delineation of multiple retinotopic fields in visual cortex. The results show that VCFs have consistent locations and extents in different subjects that provide stable and accurate landmarks for functional and anatomical mapping. Interhemispheric comparisons revealed minor differences in polar angle and eccentricity tuning in comparable VCFs in the left and right hemisphere, and somewhat greater intersubject variability in the right than left hemisphere. We then used the functional boundaries to characterize the anatomical properties of VCFs, including fractional anisotropy (FA), magnetization transfer ratio (MTR) and the ratio of T1W and T2W images and found significant anatomical differences between VCFs and between hemispheres.

Physiological measurements using ultra-high field fMRI: a review

Functional MRI (fMRI) has grown to be the neuroimaging technique of choice for investigating brain function. This topical review provides an outline of fMRI methods and applications, with a particular emphasis on the recent advances provided by ultra-high field (UHF) scanners to allow functional mapping with greater sensitivity and improved spatial specificity. A short outline of the origin of the blood oxygenation level dependent (BOLD) contrast is provided, followed by a review of BOLD fMRI methods based on gradient-echo (GE) and spin-echo (SE) contrast. Phase based fMRI measures, as well as perfusion contrast obtained with the technique of arterial spin labelling (ASL), are also discussed. An overview of 7 T based functional neuroimaging is provided, outlining the potential advances to be made and technical challenges to be addressed.

Comparing like with like: the power of knowing where you are

Magnetic resonance imaging can now provide human brain images of structure, function, and connectivity with isotropic voxels smaller than one millimeter, and thus much smaller than the cortical thickness. This resolution, achievable in a scan time of less than 1 h, enables visualization of myeloarchitectural layer structure, intracortical variations in functional activity--recorded in changes in blood oxygenation level dependent signal or cerebral blood volume CBV--and intracortical axonal orientational structure via diffusion-weighted magnetic resonance imaging. While recent improvements in radiofrequency receiver coils now enable excellent image data to be obtained at 3T, scanning at the ultra-high field of 7 T offers further gains in signal-to-noise ratio and speed of image acquisition, with a structural image resolution of about 300 μm. These improvements throw into sharp question the strategies that have become conventional for the analysis of functional imaging data, especially the practice of spatial smoothing of raw functional data before further analysis. Creation of a native cortical map for each human subject that provides a reliable individual parcellation into cortical areas related to Brodmann Areas enables a strikingly different approach to functional image analysis. This proposed approach involves surface registration of the cortices of groups of subjects using maps of the longitudinal relaxation time T1 as an index of myelination, and methods for inferring statistical significance that do not entail spatial smoothing. The outcome should be a far more precise comparison of like-with-like cortical areas across subjects, with the potential to greatly increase experimental power, to discriminate activity in neighboring cortical areas, and to enable correlation of function and connectivity with specific cytoarchitecture. Such analyses should enable a far more convincing modeling of brain mechanisms than current graph-based methods that require gross over-simplification of brain activity patterns in order to be computationally tractable.

High-Resolution Mapping of Myeloarchitecture In Vivo: Localization of Auditory Areas in the Human Brain

The precise delineation of auditory areas in vivo remains problematic. Histological analysis of postmortem tissue indicates that the relation of areal borders to macroanatomical landmarks is variable across subjects. Furthermore, functional parcellation schemes based on measures of, for example, frequency preference (tonotopy) remain controversial. Here, we propose a 7 Tesla magnetic resonance imaging method that enables the anatomical delineation of auditory cortical areas in vivo and in individual brains, through the high-resolution visualization (0.6 × 0.6 × 0.6 mm(3)) of intracortical anatomical contrast related to myelin. The approach combines the acquisition and analysis of images with multiple MR contrasts (T1, T2*, and proton density). Compared with previous methods, the proposed solution is feasible at high fields and time efficient, which allows collecting myelin-related and functional images within the same measurement session. Our results show that a data-driven analysis of cortical depth-dependent profiles of anatomical contrast allows identifying a most densely myelinated cortical region on the medial Heschl's gyrus. Analyses of functional responses show that this region includes neuronal populations with typical primary functional properties (single tonotopic gradient and narrow frequency tuning), thus indicating that it may correspond to the human homolog of monkey A1.

High-resolution imaging of magnetisation transfer and nuclear Overhauser effect in the human visual cortex at 7 T

The aim of this study was to optimise a pulse sequence for high-resolution imaging sensitive to the effects of conventional macromolecular magnetisation transfer (MT(m)) and nuclear Overhauser enhancement (NOE), and to use it to investigate variations in these parameters across the cerebral cortex. A high-spatial-resolution magnetisation transfer-prepared turbo field echo (MT-TFE) sequence was designed to have high sensitivity to MT(m) and NOE effects, whilst being robust to B0 and B1 inhomogeneities, and producing a good point spread function across the cortex. This was achieved by optimising the saturation and imaging components of the sequence using simulations based on the Bloch equations, including exchange and an image simulator. This was used to study variations in these parameters across the cortex. Using the sequence designed to be sensitive to NOE and MT(m), a variation in signals corresponding to a variation in MT(m) and NOE across the cortex, consistent with a reduction in myelination from the white matter surface to the pial surface of the cortex, was observed. In regions in which the stria was visible on T2*-weighted images, it could also be detected in signals sensitive to MT(m) and NOE. There was greater variation in signals sensitive to NOE, suggesting that the NOE signal is more sensitive to myelination. A sequence has been designed to image variations in MT(m) and NOE at high spatial resolution and has been used to investigate variations in contrast in these parameters across the cortex.

BigBrain: an ultrahigh-resolution 3D human brain model

Reference brains are indispensable tools in human brain mapping, enabling integration of multimodal data into an anatomically realistic standard space. Available reference brains, however, are restricted to the macroscopic scale and do not provide information on the functionally important microscopic dimension. We created an ultrahigh-resolution three-dimensional (3D) model of a human brain at nearly cellular resolution of 20 micrometers, based on the reconstruction of 7404 histological sections. "BigBrain" is a free, publicly available tool that provides considerable neuroanatomical insight into the human brain, thereby allowing the extraction of microscopic data for modeling and simulation. BigBrain enables testing of hypotheses on optimal path lengths between interconnected cortical regions or on spatial organization of genetic patterning, redefining the traditional neuroanatomy maps such as those of Brodmann and von Economo.

Trends and properties of human cerebral cortex: correlations with cortical myelin content

"In vivo Brodmann mapping" or non-invasive cortical parcellation using MRI, especially by measuring cortical myelination, has recently become a popular research topic, though myeloarchitectonic cortical parcellation in humans previously languished in favor of cytoarchitecture. We review recent in vivo myelin mapping studies and discuss some of the different methods for estimating myelin content. We discuss some ways in which myelin maps may improve surface registration and be useful for cross-modal and cross-species comparisons, including some preliminary cross-species results. Next, we consider neurobiological aspects of why some parts of cortex are more myelinated than others. Myelin content is inversely correlated with intracortical circuit complexity - in general, more myelin content means simpler and perhaps less dynamic intracortical circuits. Using existing PET data and functional network parcellations, we examine metabolic differences in the differently myelinated cortical functional networks. Lightly myelinated cognitive association networks tend to have higher aerobic glycolysis than heavily myelinated early sensory-motor ones, perhaps reflecting greater ongoing dynamic anabolic cortical processes. This finding is consistent with the hypothesis that intracortical myelination may stabilize intracortical circuits and inhibit synaptic plasticity. Finally, we discuss the future of the in vivo myeloarchitectural field and cortical parcellation--"in vivo Brodmann mapping"--in general.

What can we learn from T2* maps of the cortex?

Studies have shown that T2* contrast can reveal features of cortical anatomy. However, understanding the relationship between T2* contrast and the underlying cyto- and myelo-architecture is not an easy task, given the number of confounds, such as myelin, iron, blood vessels and structure orientation. Moreover, it is difficult to obtain reliable T2* measurements in the cortex due to its thin and folded geometry and the presence of artifacts. This review addresses issues associated with T2* mapping in the human cortex. After describing the theory behind T2* relaxation, a list of practical steps is proposed to reliably acquire and process T2* data and then map these values within the cortex using surface-based analysis. The last section addresses the question: "What can we gain from T2* cortical mapping?", with particular emphasis on Brodmann mapping.