In vivo myeloarchitectonic analysis of human striate and extrastriate cortex using magnetic resonance imaging

Toward in vivo histology: a comparison of quantitative susceptibility mapping (QSM) with magnitude-, phase-, and R2*-imaging at ultra-high magnetic field strength

Quantitative magnetic susceptibility mapping (QSM) has recently been introduced to provide a novel quantitative and local MRI contrast. However, the anatomical contrast represented by in vivo susceptibility maps has not yet been compared systematically and comprehensively with gradient (recalled) echo (GRE) magnitude, frequency, and R(2)(*) images. Therefore, this study compares high-resolution quantitative susceptibility maps with conventional GRE imaging approaches (magnitude, frequency, R(2)(*)) in healthy individuals at 7 T with respect to anatomic tissue contrast. Volumes-of-interest were analyzed in deep and cortical gray matter (GM) as well as in white matter (WM) on R(2)(*) and susceptibility maps. High-resolution magnetic susceptibility maps of the human brain exhibited superb contrast that allowed the identification of substructures of the thalamus, midbrain and basal ganglia, as well as of the cerebral cortex. These were consistent with histology but not generally visible on magnitude, frequency or R(2)(*)-maps. Common target structures for deep brain stimulation, including substantia nigra pars reticulata, ventral intermediate nucleus, subthalamic nucleus, and the substructure of the internal globus pallidus, were clearly distinguishable from surrounding tissue on magnetic susceptibility maps. The laminar substructure of the cortical GM differed depending on the anatomical region, i.e., a cortical layer with increased magnetic susceptibility, corresponding to the Stria of Gennari, was found in the GM of the primary visual cortex, V1, whereas a layer with reduced magnetic susceptibility was observed in the GM of the temporal cortex. Both magnetic susceptibility and R(2)(*) values differed substantially in cortical GM depending on the anatomic regions. Regression analysis between magnetic susceptibility and R(2)(*) values of WM and GM structures suggested that variations in myelin content cause the overall contrast between gray and white matter on susceptibility maps and that both R(2)(*) and susceptibility values provide linear measures for iron content in GM. In conclusion, quantitative magnetic susceptibility mapping provides a non-invasive and spatially specific contrast that opens the door to the assessment of diseases characterized by variation in iron and/or myelin concentrations. Its ability to reflect anatomy of deep GM structures with superb delineation may be useful for neurosurgical applications.

Mapping the human cortical surface by combining quantitative T(1) with retinotopy

We combined quantitative relaxation rate (R₁= 1/T₁) mapping—to measure local myelination—with fMRI-based retinotopy. Gray–white and pial surfaces were reconstructed and used to sample R₁ at different cortical depths. Like myelination, R₁ decreased from deeper to superficial layers. R₁ decreased passing from V1 and MT, to immediately surrounding areas, then to the angular gyrus. High R₁ was correlated across the cortex with convex local curvature so the data was first “de-curved”. By overlaying R₁ and retinotopic maps, we found that many visual area borders were associated with significant R₁ increases including V1, V3A, MT, V6, V6A, V8/VO1, FST, and VIP. Surprisingly, retinotopic MT occupied only the posterior portion of an oval-shaped lateral occipital R₁ maximum. R₁ maps were reproducible within individuals and comparable between subjects without intensity normalization, enabling multi-center studies of development, aging, and disease progression, and structure/function mapping in other modalities.

Improved prediction of Alzheimer's disease with longitudinal white matter/gray matter contrast changes

Brain morphometry measures derived from magnetic resonance imaging (MRI) are important biomarkers for Alzheimer's disease (AD). The objective of the present study was to test whether we could improve morphometry-based detection and prediction of disease state by use of white matter/gray matter (WM/GM) signal intensity contrast obtained from conventional MRI scans. We hypothesized that including WM/GM contrast change along with measures of atrophy in the entorhinal cortex and the hippocampi would yield better classification of AD patients, and more accurate prediction of early disease progression. T1 -weighted MRI scans from two sessions approximately 2 years apart from 78 participants with AD (Clinical Dementia Rating (CDR) = 0.5-2) and 71 age-matched controls were used to calculate annual change rates. Results showed that WM/GM contrast decay was larger in AD compared with controls in the medial temporal lobes. For the discrimination between AD and controls, entorhinal WM/GM contrast decay contributed significantly when included together with decrease in entorhinal cortical thickness and hippocampal volume, and increased the area under the curve to 0.79 compared with 0.75 when using the two morphometric variables only. Independent effects of WM/GM contrast decay and improved classification were also observed for the CDR-based subgroups, including participants converting from either a non-AD status to very mild AD, or from very mild to mild AD. Thus, WM/GM contrast decay increased diagnostic accuracy beyond what was obtained by well-validated morphometric measures alone. The findings suggest that signal intensity properties constitute a sensitive biomarker for cerebral degeneration in AD.

Generalized multiple-layer appearance of the cerebral cortex with 3D FLAIR 7.0-T MR imaging

To examine the multiple-layer appearance of the cerebral cortex with fluid-attenuated inversion recovery (FLAIR) magnetic resonance (MR) imaging at 7.0 T, whole-brain volumetric three-dimensional (3D) magnetization prepared FLAIR images were acquired in 12 volunteers (0.8 3 0.8 3 0.8-mm spatial resolution). Signal intensity profiles were evaluated for the anterior frontal (Brodmann area [BA] 10), posterior frontal (BA 6), parietal (BA 7), precentral (BA 4), postcentral (BA 3), occipital (BA 18), and calcarine (BA 17) regions. Variance of the normalized profile was used as the metric for the multiple-layer appearance. Wilcoxon signed-rank tests were performed to compare variances of the profiles between all areas. All cortical areas showed multiple-layered appearances, with a prominent hyperintense band at the external surface of the cortex, a hypointense band deeper in the cortex, and a hyperintense third band. The ranking from least- to most-pronounced layer appearance was as follows: postcentral (variance, 0.04), posterior frontal (variance, 0.05), calcarine (variance, 0.05), precentral (variance, 0.06), parietal (variance, 0.08), anterior frontal (variance, 0.10), and occipital (variance, 0.11). Each region was significantly different from at least one other region. In conclusion, a multiple-layer appearance of the cerebral cortex was found for all cortical regions with high-spatial-resolution 3D FLAIR MR imaging at 7.0 T.

T₂* mapping and B₀ orientation-dependence at 7 T reveal cyto- and myeloarchitecture organization of the human cortex

Ultra-high field MRI (≥ 7 T) has recently shown great sensitivity to depict patterns of tissue microarchitecture. Moreover, recent studies have demonstrated a dependency between T₂* and orientation of white matter fibers with respect to the main magnetic field B₀. In this study we probed the potential of T₂* mapping at 7 T to provide new markers of cortical architecture. We acquired multi-echo measurements at 7 T and mapped T₂* over the entire cortex of eight healthy individuals using surface-based analysis. B₀ dependence was tested by computing the angle θ(z) between the normal of the surface and the direction of B₀, then fitting T₂*(θ(z)) using model from the literature. Average T₂* in the cortex was 32.20 +/- 1.35 ms. Patterns of lower T₂* were detected in the sensorimotor, visual and auditory cortices, likely reflecting higher myelin content. Significantly lower T₂* was detected in the left hemisphere of the auditory region (p<0.005), suggesting higher myelin content, in accordance with previous investigations. B₀ orientation dependence was detected in some areas of the cortex, the strongest being in the primary motor cortex (∆R₂*=4.10 Hz). This study demonstrates that quantitative T₂* measures at 7 T MRI can reveal patterns of cytoarchitectural organization of the human cortex in vivo and that B₀ orientation dependence can probe the coherency and orientation of gray matter fibers in the cortex, shedding light into the potential use of this type of contrast to characterize cyto-/myeloarchitecture and to understand the pathophysiology of diseases associated with changes in iron and/or myelin concentration.

A history of randomized task designs in fMRI

In the early days of fMRI, data were acquired using methods adapted mainly from PET imaging. Sets of similar stimuli were presented in extended blocks, with stimulus conditions changing from block to block. While this method provided optimum statistical power, it also presented a variety of potential confounds, including changes in attention, alertness, expectancy, and practice effects within and between blocks. Event-related paradigms using unpredictable, randomized stimulus sequences had been used in EEG studies for over 50 years before the development of fMRI, and provided a means to overcome these issues. However, the temporal dispersion of BOLD fMRI activity resulted in responses to successive stimuli adding together, making it difficult to perform rapid event-related paradigms using fMRI. Here we describe the background and history of methods developed to overcome this limitation, allowing rapid, randomized stimulus sequences to be used with fMRI. The advantages and disadvantages of this technique and how these methods have been applied in a variety of experimental settings are discussed.

Entorhinal verrucae geometry is coincident and correlates with Alzheimer's lesions: a combined neuropathology and high-resolution ex vivo MRI analysis

Entorhinal cortex displays a distinctive organization in layer II and forms small elevations on its surface called entorhinal verrucae. In Alzheimer's disease, the verrucae disappear due to neurofibrillary tangle formation and neuronal death. Isosurface models were reconstructed from high-resolution ex vivo MRI volumes scanned at 7.0 T and individual verruca were measured quantitatively for height, width, volume, and surface area on control and mild Alzheimer's cases. Mean verruca height was 0.13 ± 0.04 mm for our cognitively normal (controls) sample set whereas for mild AD samples mean height was 0.11 mm ± 0.05 mm (p < 0.001) in entorhinal cortex (n = 10 cases). These quantitative methods were validated by a significant correlation of verrucae height and volume with qualitative verrucae ratings (n = 36 cases). Entorhinal surfaces were significantly different from other cortical heights such as, cingulate, frontal, occipital, parietal and temporal cortices. Colocalization of verrucae with entorhinal islands was confirmed in ex vivo MRI and, moreover, verrucae ratings were negatively correlated to Braak and Braak pathological stage. This study characterizes novel methods to measure individual entorhinal verruca size, and shows that verrucae size correlates to Alzheimer's pathology. Taken together, these results suggest that verrucae may have the potential to serve as an early and specific morphological marker for mild cognitive impairment and Alzheimer's disease.

Identification of cortical lamination in awake monkeys by high resolution magnetic resonance imaging

Brodmann divided the neocortex into 47 different cortical areas based on histological differences in laminar myeloarchitectonic and cytoarchitectonic defined structure. The ability to do so in vivo with anatomical magnetic resonance (MR) methods in awake subjects would be extremely advantageous for many functional studies. However, due to the limitations of spatial resolution and contrast, this has been difficult to achieve in awake subjects. Here, we report that by using a combination of MR microscopy and novel contrast effects, cortical layers can be delineated in the visual cortex of awake subjects (nonhuman primates) at 4.7 T. We obtained data from 30-min acquisitions at voxel size of 62.5 × 62.5 × 1000 μm(3) (4 nl). Both the phase and magnitude components of the T(2)*-weighted image were used to generate laminar profiles which are believed to reflect variations in myelin and local cell density content across cortical depth. Based on this, we were able to identify six layers characteristic of the striate cortex (V1). These were the stripe of Kaes-Bechterew (in layer II/III), the stripe of Gennari (in layer IV), the inner band of Baillarger (in layer V), as well as three sub-layers within layer IV (IVa, IVb, and IVc). Furthermore, we found that the laminar structure of two extrastriate visual cortex (V2, V4) can also be detected. Following the tradition of Brodmann, this significant improvement in cortical laminar visualization should make it possible to discriminate cortical regions in awake subjects corresponding to differences in myeloarchitecture and cytoarchitecture.

Novel frontiers in ultra-structural and molecular MRI of the brain

PURPOSE OF REVIEW

Recent developments in the MRI of the brain continue to expand its use in basic and clinical neuroscience. This review highlights some areas of recent progress.

RECENT FINDINGS

Higher magnetic field strengths and improved signal detectors have allowed improved visualization of the various properties of the brain, facilitating the anatomical definition of function-specific areas and their connections. For example, by sensitizing the MRI signal to the magnetic susceptibility of tissue, it is starting to become possible to reveal the laminar structure of the cortex and identify millimeter-scale fiber bundles. Using exogenous contrast agents, and innovative ways to manipulate contrast, it is becoming possible to highlight specific fiber tracts and cell populations. These techniques are bringing us closer to understanding the evolutionary blueprint of the brain, improving the detection and characterization of disease, and help to guide treatment.

SUMMARY

Recent MRI techniques are leading to more detailed and more specific contrast in the study of the brain.

Mapping human cortical areas in vivo based on myelin content as revealed by T1- and T2-weighted MRI

Noninvasively mapping the layout of cortical areas in humans is a continuing challenge for neuroscience. We present a new method of mapping cortical areas based on myelin content as revealed by T1-weighted (T1w) and T2-weighted (T2w) MRI. The method is generalizable across different 3T scanners and pulse sequences. We use the ratio of T1w/T2w image intensities to eliminate the MR-related image intensity bias and enhance the contrast to noise ratio for myelin. Data from each subject were mapped to the cortical surface and aligned across individuals using surface-based registration. The spatial gradient of the group average myelin map provides an observer-independent measure of sharp transitions in myelin content across the surface--i.e., putative cortical areal borders. We found excellent agreement between the gradients of the myelin maps and the gradients of published probabilistic cytoarchitectonically defined cortical areas that were registered to the same surface-based atlas. For other cortical regions, we used published anatomical and functional information to make putative identifications of dozens of cortical areas or candidate areas. In general, primary and early unimodal association cortices are heavily myelinated and higher, multimodal, association cortices are more lightly myelinated, but there are notable exceptions in the literature that are confirmed by our results. The overall pattern in the myelin maps also has important correlations with the developmental onset of subcortical white matter myelination, evolutionary cortical areal expansion in humans compared with macaques, postnatal cortical expansion in humans, and maps of neuronal density in non-human primates.

Fractional anisotropy of cerebral white matter and thickness of cortical gray matter across the lifespan

We examined age trajectories of fractional anisotropy (FA) of cerebral white matter (WM) and thickness of cortical gray matter (GM) in 1031 healthy human subjects (aged 11-90 years). Whole-brain FA and GM thickness values followed quadratic trajectories with age but the relationship between them was linear, indicating that a putative biological mechanism may explain the non-linearity of their age trajectories. Inclusion of the FA values into the quadratic model of the whole-brain and regional GM thickness changes with age made the effect of the age(2) term no longer significant for the whole-brain GM thickness and greatly reduced its significance for regional GM thickness measurements. The phylogenetic order of cerebral myelination helped to further explain the intersubject variability in GM thickness. FA values for the early maturing WM were significantly better (p=10(-6)) at explaining variability in GM thickness in maturing (aged 11-20) subjects than FA values for the late maturing WM. The opposite trend was observed for aging subjects (aged 40-90) where FA values for the late maturing WM were better (p=10(-16)) at explaining the variability in GM thickness. We concluded that the non-linearity of the age trajectory for GM thickness, measured from T1-weighted MRI, was partially explained by the heterogeneity and the heterochronicity of the age-related changes in the microintegrity of cerebral WM. We consider these findings as the evidence that the measurements of age-related changes in GM thickness and FA are driven, in part, by a common biological mechanism, presumed to be related to changes in cerebral myelination.

Temporal dynamics and spatial specificity of arterial and venous blood volume changes during visual stimulation: implication for BOLD quantification

Determination of compartment-specific cerebral blood volume (CBV) changes is important for understanding neurovascular physiology and quantifying blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI). In isoflurane-anesthetized cats, we measured the spatiotemporal responses of arterial CBV (CBV(a)) and total CBV (CBV(t)) induced by a 40-second visual stimulation, using magnetization transfer (MT)-varied BOLD and contrast-agent fMRI techniques at 9.4 T. To determine the venous CBV (CBV(v)) change, we calculated the difference between CBV(t) and CBV(a) changes. The dynamic response of CBV(a) was an order of magnitude faster than that of CBV(v), while the magnitude of change under steady-state conditions was similar between the two. Following stimulation offset, ΔCBV(a) showed small poststimulus undershoots, while ΔCBV(v) slowly returned to baseline. The largest CBV(a) and CBV(t) response occurred after 10 seconds of simulation in cortical layer 4, which we identified as the stripe of Gennari by T(1)-weighted MRI. The CBV(v) response, however, was not specific across the cortical layers during the entire stimulation period. Our data indicate that rapid, more-specific arterial vasodilation is followed by slow, less-specific venous dilation. Our finding implies that the contribution of CBV(v) changes to BOLD signals is significant for long, but not short, stimulation periods.

Microstructural Parcellation of the Human Cerebral Cortex - From Brodmann's Post-Mortem Map to in vivo Mapping with High-Field Magnetic Resonance Imaging

The year 2009 marked the 100th anniversary of the publication of the famous brain map of Korbinian Brodmann. Although a “classic” guide to microanatomical parcellation of the cerebral cortex, it is – from today's state-of-the-art neuroimaging perspective – problematic to use Brodmann's map as a structural guide to functional units in the cortex. In this article we discuss some of the reasons, especially the problematic compatibility of the “post-mortem world” of microstructural brain maps with the “in vivo world” of neuroimaging. We conclude with some prospects for the future of in vivo structural brain mapping: a new approach which has the enormous potential to make direct correlations between microstructure and function in living human brains: “in vivo Brodmann mapping” with high-field magnetic resonance imaging.

Differentiating maturational and aging-related changes of the cerebral cortex by use of thickness and signal intensity

Cortical thickness decreases from childhood throughout life, as estimated by magnetic resonance imaging (MRI). This monotone trajectory does not reflect the fundamentally different neurobiological processes underlying morphometric changes in development versus aging. We hypothesized that intracortical gray matter (GM) and subjacent white matter (WM) T1-weighted signal intensity would distinguish developmental and age-related changes in the cortex better than thickness. Intracortical GM and subjacent WM signal intensity and cortical thickness was measured across the brain surface in a healthy life span sample (n=429, 8-85 years). We also computed the relaxation rate of T2* (R2*) from multiecho sequences and mapped intracortical GM and subjacent WM values to the surface to delineate age-related variability in R2* and to adjust the T1 signal intensity for possible confounds of accumulated iron. While monotone age-related reductions in thickness were found, both intracortical GM and subcortical WM signal intensity showed inverted U patterns with peaks from eight to approximately 30 years of age. The spatial pattern of intracortical neurodevelopment followed a posterior-anterior gradient, with earliest maturation of occipital visual cortices and most protracted in superior frontal regions. From 50s and 60s, substantial signal reductions were observed in several regions, including the insula, cingulate, and inferior temporal gyrus. R2* showed similar patterns but peaked much later than the T1-weighted signal intensity measures. The results are presented as animations yielding detailed depictions of the dynamic regional variability in cortical neurodevelopment and aging and demonstrate that cortical thickness and T1-weighted signal intensity are sensitive to different cortical maturational and aging-related processes.

Layer-specific variation of iron content in cerebral cortex as a source of MRI contrast

Recent advances in high-field MRI have dramatically improved the visualization of human brain anatomy in vivo. Most notably, in cortical gray matter, strong contrast variations have been observed that appear to reflect the local laminar architecture. This contrast has been attributed to subtle variations in the magnetic properties of brain tissue, possibly reflecting varying iron and myelin content. To establish the origin of this contrast, MRI data from postmortem brain samples were compared with electron microscopy and histological staining for iron and myelin. The results show that iron is distributed over laminae in a pattern that is suggestive of each region's myeloarchitecture and forms the dominant source of the observed MRI contrast.

In vivo detection of individual glomeruli in the rodent olfactory bulb using manganese enhanced MRI

MRI contrast based on relaxation times, proton density, or signal phase have been applied to delineate neural structures in the brain. However, neural units such as cortical layers and columns have been difficult to identify using these methods. Manganese ion delivered either systemically or injected directly has been shown to accumulate specifically within cellular areas of the brain enabling the differentiation of layers within the hippocampus, cortex, cerebellum, and olfactory bulb in vivo. Here we show the ability to detect individual olfactory glomeruli using manganese enhanced MRI (MEMRI). Glomeruli are anatomically distinct structures ( approximately 150 microm in diameter) on the surface of the olfactory bulb that represent the first processing units for olfactory sensory information. Following systemic delivery of MnCl(2) we used 3D-MRI with 50 microm isotropic resolution to detect discrete spots of increased signal intensity between 100 and 200 microm in diameter in the glomerular layer of the rat olfactory bulb. Inflow effects of arterial blood and susceptibility effects of venous blood were suppressed and were evaluated by comparing the location of vessels in the bulb to areas of manganese enhancement using iron oxide to increase vessel contrast. These potential vascular effects did not explain the contrast detected. Nissl staining of individual glomeruli were also compared to MEMRI images from the same animals clearly demonstrating that many of the manganese enhanced regions corresponded to individual olfactory glomeruli. Thus, MEMRI can be used as a non-invasive means to detect olfactory glomeruli for longitudinal studies looking at neural plasticity during olfactory development or possible degeneration associated with disease.

Visualizing the entire cortical myelination pattern in marmosets with magnetic resonance imaging

Myeloarchitecture, the pattern of myelin density across the cerebral cortex, has long been visualized in histological sections to identify distinct anatomical areas of the cortex. In humans, two-dimensional (2D) magnetic resonance imaging (MRI) has been used to visualize myeloarchitecture in select areas of the cortex, such as the stripe of Gennari in the primary visual cortex and Heschl's gyrus in the primary auditory cortex. Here, we investigated the use of MRI contrast based on longitudinal relaxation time (T(1)) to visualize myeloarchitecture in vivo over the entire cortex of the common marmoset (Callithrix jacchus), a small non-human primate that is becoming increasingly important in neuroscience and neurobiology research. Using quantitative T(1) mapping, we found that T(1) at 7T in a cortical region with a high myelin content was 15% shorter than T(1) in a region with a low myelin content. To maximize this T(1) contrast for imaging cortical myelination patterns, we optimized a magnetization-prepared rapidly acquired gradient echo (MP-RAGE) sequence. In whole-brain, 3D T(1)-weighted images made in vivo with the sequence, we identified six major cortical areas with high myelination and confirmed the results with histological sections stained for myelin. We also identified several subtle features of myeloarchitecture, showing the sensitivity of our technique. The ability to image myeloarchitecture over the entire cortex may prove useful in studies of longitudinal changes of the topography of the cortex associated with development and neuronal plasticity, as well as for guiding and confirming the location of functional measurements.

Plasticity and stability of visual field maps in adult primary visual cortex

It is important to understand the balance between cortical plasticity and stability in various systems and across spatial scales in the adult brain. Here we review studies of adult plasticity in primary visual cortex (V1), which has a key role in distributing visual information. There are claims of plasticity at multiple spatial scales in adult V1, but a number of inconsistencies in the supporting data raise questions about the extent and nature of such plasticity. Our understanding of the extent of plasticity in V1 is further limited by a lack of quantitative models to guide the interpretation of the data. These problems limit efforts to translate research findings about adult cortical plasticity into significant clinical, educational and policy applications.

MRI of cellular layers in mouse brain in vivo

Noninvasive imaging of the brain of animal models demands the detection of increasingly smaller structures by in vivo MRI. The purpose of this work was to elucidate the spatial resolution and structural contrast that can be obtained for studying the brain of C57BL/6J mice by optimized T2-weighted fast spin-echo MRI at 9.4 T. As a prerequisite for high-resolution imaging in vivo, motion artifacts were abolished by combining volatile anesthetics and positive pressure ventilation with a specially designed animal bed for fixation. Multiple substructures in the cortex, olfactory bulb, hippocampus, and cerebellum were resolved at 30 to 40 microm in-plane resolution and 200 to 300 microm section thickness as well as for relatively long echo times of 65 to 82 ms. In particular, the approach resulted in the differentiation of up to five cortical layers. In the olfactory bulb the images unraveled the mitral cell layer which has a thickness of mostly single cells. In the hippocampus at least five substructures could be separated. The molecular layer, Purkinje layer, and granular layer of the cerebellum could be clearly differentiated from the white matter. In conclusion, even without the use of a contrast agent, suitable adjustments of a widely available T2-weighted MRI sequence at high field allow for structural MRI of living mice at near single-cell layer resolution.

Predicting the location of entorhinal cortex from MRI

Entorhinal cortex (EC) is a medial temporal lobe area critical to memory formation and spatial navigation that is among the earliest parts of the brain affected by Alzheimer's disease (AD). Accurate localization of EC would thus greatly facilitate early detection and diagnosis of AD. In this study, we used ultra-high resolution ex vivo MRI to directly visualize the architectonic features that define EC rostrocaudally and mediolaterally, then applied surface-based registration techniques to quantify the variability of EC with respect to cortical geometry, and made predictions of its location on in vivo scans. The results indicate that EC can be localized quite accurately based on cortical folding patterns, within 3 mm in vivo, a significant step forward in our ability to detect the earliest effects of AD when clinical intervention is most likely to be effective.

Locating the functional and anatomical boundaries of human primary visual cortex

The primary visual cortex (V1) can be delineated both functionally by its topographic map of the visual field and anatomically by its distinct pattern of laminar myelination. Although it is commonly assumed that the specialized anatomy V1 exhibits corresponds in location with functionally defined V1, demonstrating this in human has not been possible thus far due to the difficulty of determining the location of V1 both functionally and anatomically in the same individual. In this study we use MRI to measure the anatomical and functional V1 boundaries in the same individual and demonstrate close agreement between them. Functional V1 location was measured by parcellating occipital cortex of 10 living humans into visual cortical areas based on the topographic map of the visual field measured using functional MRI. Anatomical V1 location was estimated for these same subjects using a surface-based probabilistic atlas derived from high-resolution structural MRI of the stria of Gennari in 10 intact ex vivo human hemispheres. To ensure that the atlas prediction was correct, it was validated against V1 location measured using an observer-independent cortical parcellation based on the laminar pattern of cell density in serial brain sections from 10 separate individuals. The close agreement between the independent anatomically and functionally derived V1 boundaries indicates that the whole extent of V1 can be accurately predicted based on cortical surface reconstructions computed from structural MRI scans, eliminating the need for functional localizers of V1. In addition, that the primary cortical folds predict the location of functional V1 suggests that the mechanism giving rise to V1 location is tied to the development of the cortical folds.

Magnetic resonance imaging of neural circuits

A major goal of modern MRI research is to be able to image neural circuits in the central nervous system. Critical to this mission is the ability to describe a number of important parameters associated with neural circuits. These parameters include neural architecture, functional activation of neural circuits, anatomical and functional connectivity of neural circuits, and factors that might alter neural circuits, such as trafficking of immune cells and brain precursor cells in the brain. Remarkably, a variety of work in human and animal brains has demonstrated that all these features of neural circuits can be visualized with MRI. In this Article we provide a brief summary of the new directions in neural imaging research, which should prove useful in future analyses of normal and pathological human brains and in studies of animal models of neurological and psychiatric disorders. At present, few MRI data characterizing the neural circuits in the heart are available, but in this Article we discuss the applicable present developments and the prospects for the future.

The intrinsic shape of human and macaque primary visual cortex

Previous studies have reported considerable variability in primary visual cortex (V1) shape in both humans and macaques. Here, we demonstrate that much of this variability is due to the pattern of cortical folds particular to an individual and that V1 shape is similar among individual humans and macaques as well as between these 2 species. Human V1 was imaged ex vivo using high-resolution (200 microm) magnetic resonance imaging at 7 T. Macaque V1 was identified in published histological serial section data. Manual tracings of the stria of Gennari were used to construct a V1 surface, which was computationally flattened with minimal metric distortion of the cortical surface. Accurate flattening allowed investigation of intrinsic geometric features of cortex, which are largely independent of the highly variable cortical folds. The intrinsic shape of V1 was found to be similar across human subjects using both nonparametric boundary matching and a simple elliptical shape model fit to the data and is very close to that of the macaque monkey. This result agrees with predictions derived from current models of V1 topography. In addition, V1 shape similarity suggests that similar developmental mechanisms are responsible for establishing V1 shape in these 2 species.

Cortical folding patterns and predicting cytoarchitecture

The human cerebral cortex is made up of a mosaic of structural areas, frequently referred to as Brodmann areas (BAs). Despite the widespread use of cortical folding patterns to perform ad hoc estimations of the locations of the BAs, little is understood regarding 1) how variable the position of a given BA is with respect to the folds, 2) whether the location of some BAs is more variable than others, and 3) whether the variability is related to the level of a BA in a putative cortical hierarchy. We use whole-brain histology of 10 postmortem human brains and surface-based analysis to test how well the folds predict the locations of the BAs. We show that higher order cortical areas exhibit more variability than primary and secondary areas and that the folds are much better predictors of the BAs than had been previously thought. These results further highlight the significance of cortical folding patterns and suggest a common mechanism for the development of the folds and the cytoarchitectonic fields.

Complete pattern of ocular dominance columns in human primary visual cortex

The occipital lobes were obtained after death from six adult subjects with monocular visual loss. Flat-mounts were processed for cytochrome oxidase (CO) to reveal metabolic activity in the primary (V1) and secondary (V2) visual cortices. Mean V1 surface area was 2643 mm2 (range, 1986-3477 mm2). Ocular dominance columns were present in all cases, having a mean width of 863 microm. There were 78-126 column pairs along the V1 perimeter. Human column patterns were highly variable, but in at least one person they resembled a scaled-up version of macaque columns. CO patches in the upper layers were centered on ocular dominance columns in layer 4C, with one exception. In this individual, the columns in a local area resembled those present in the squirrel monkey, and no evidence was found for column/patch alignment. In every subject, the blind spot of the contralateral eye was conspicuous as an oval region without ocular dominance columns. It provided a precise landmark for delineating the central 15 degrees of the visual field. A mean of 53.1% of striate cortex was devoted to the representation of the central 15 degrees. This fraction was less than the proportion of striate cortex allocated to the representation of the central 15 degrees in the macaque. Within the central 15 degrees, each eye occupied an equal territory. Beyond this eccentricity, the contralateral eye predominated, occupying 63% of the cortex. In one subject, monocular visual loss began at age 4 months, causing shrinkage of ocular dominance columns. In V2, which had a larger surface area than V1, CO stripes were present but could not be classified as thick or thin.

Accurate prediction of V1 location from cortical folds in a surface coordinate system

Previous studies demonstrated substantial variability of the location of primary visual cortex (V1) in stereotaxic coordinates when linear volume-based registration is used to match volumetric image intensities [Amunts, K., Malikovic, A., Mohlberg, H., Schormann, T., and Zilles, K. (2000). Brodmann's areas 17 and 18 brought into stereotaxic space-where and how variable? Neuroimage, 11(1):66-84]. However, other qualitative reports of V1 location [Smith, G. (1904). The morphology of the occipital region of the cerebral hemisphere in man and the apes. Anatomischer Anzeiger, 24:436-451; Stensaas, S.S., Eddington, D.K., and Dobelle, W.H. (1974). The topography and variability of the primary visual cortex in man. J Neurosurg, 40(6):747-755; Rademacher, J., Caviness, V.S., Steinmetz, H., and Galaburda, A.M. (1993). Topographical variation of the human primary cortices: implications for neuroimaging, brain mapping, and neurobiology. Cereb Cortex, 3(4):313-329] suggested a consistent relationship between V1 and the surrounding cortical folds. Here, the relationship between folds and the location of V1 is quantified using surface-based analysis to generate a probabilistic atlas of human V1. High-resolution (about 200 microm) magnetic resonance imaging (MRI) at 7 T of ex vivo human cerebral hemispheres allowed identification of the full area via the stria of Gennari: a myeloarchitectonic feature specific to V1. Separate, whole-brain scans were acquired using MRI at 1.5 T to allow segmentation and mesh reconstruction of the cortical gray matter. For each individual, V1 was manually identified in the high-resolution volume and projected onto the cortical surface. Surface-based intersubject registration [Fischl, B., Sereno, M.I., Tootell, R.B., and Dale, A.M. (1999b). High-resolution intersubject averaging and a coordinate system for the cortical surface. Hum Brain Mapp, 8(4):272-84] was performed to align the primary cortical folds of individual hemispheres to those of a reference template representing the average folding pattern. An atlas of V1 location was constructed by computing the probability of V1 inclusion for each cortical location in the template space. This probabilistic atlas of V1 exhibits low prediction error compared to previous V1 probabilistic atlases built in volumetric coordinates. The increased predictability observed under surface-based registration suggests that the location of V1 is more accurately predicted by the cortical folds than by the shape of the brain embedded in the volume of the skull. In addition, the high quality of this atlas provides direct evidence that surface-based intersubject registration methods are superior to volume-based methods at superimposing functional areas of cortex and therefore are better suited to support multisubject averaging for functional imaging experiments targeting the cerebral cortex.

High-field MRI of brain cortical substructure based on signal phase

The ability to detect brain anatomy and pathophysiology with MRI is limited by the contrast-to-noise ratio (CNR), which depends on the contrast mechanism used and the spatial resolution. In this work, we show that in MRI of the human brain, large improvements in contrast to noise in high-resolution images are possible by exploiting the MRI signal phase at high magnetic field strength. Using gradient-echo MRI at 7.0 tesla and a multichannel detector, a nominal voxel size of 0.24 x 0.24 x 1.0 mm3 (58 nl) was achieved. At this resolution, a strong phase contrast was observed both between as well as within gray matter (GM) and white matter (WM). In gradient-echo phase images obtained on normal volunteers at this high resolution, the CNR between GM and WM ranged from 3:1 to 20:1 over the cortex. This CNR is an almost 10-fold improvement over conventional MRI techniques that do not use image phase, and it is an approximately 100-fold improvement when including the gains in resolution from high-field and multichannel detection. Within WM, phase contrast appeared to be associated with the major fiber bundles, whereas contrast within GM was suggestive of the underlying layer structure. The observed phase contrast is attributed to local variations in magnetic susceptibility, which, at least in part, appeared to originate from iron stores. The ability to detect cortical substructure from MRI phase contrast at high field is expected to greatly enhance the study of human brain anatomy in vivo.

Timing of V1/V2 and V5+ activations during coherent motion of dots: an MEG study

In order to study the temporal activation course of visual areas V1 and V5 in response to a motion stimulus, a random dots kinematogram paradigm was applied to eight subjects while magnetic fields were recorded using magnetoencephalography (MEG). Sources generating the registered magnetic fields were localized with Magnetic Field Tomography (MFT). Anatomical identification of cytoarchitectonically defined areas V1/V2 and V5 was achieved by means of probabilistic cytoarchitectonic maps. We found that the areas V1/V2 and V5+ (V5 and other adjacent motion sensitive areas) exhibited two main activations peaks at 100-130 ms and at 140-200 ms after motion onset. The first peak found for V1/V2, which corresponds to the visual evoked field (VEF) M1, always preceded the peak found in V5+. Additionally, the V5+ peak was correlated significantly and positively with the second V1/V2 peak. This result supports the idea that the M1 component is generated not only by the visual area V1/V2 (as it is usually proposed), but also by V5+. It reflects a forward connection between both structures, and a feedback projection to V1/V2, which provokes a second activation in V1/V2 around 200 ms. This second V1/V2 activation (corresponding to motion VEF M2) appeared earlier than the second V5+ activation but both peaked simultaneously. This result supports the hypothesis that both areas also generate the M2 component, which reflects a feedback input from V5+ to V1/V2 and a crosstalk between both structures. Our study indicates that during visual motion analysis, V1/V2 and V5+ are activated repeatedly through forward and feedback connections and both contribute to m-VEFs M1 and M2.

Observer-independent analysis of high-resolution MR images of the human cerebral cortex: in vivo delineation of cortical areas

Using high-resolution MRI, it is now possible to examine the living human cortex down to a resolution of less than 300 mum. Thus, in vivo imaging is now approaching the resolution that has been successfully used in histological analysis of the cerebral cortex for many years, e.g., low-magnification light microscopy. This allows unprecedented views of cortical microstructure that reflect defined histological features, specifically, individual cortical layers. As in histological brain mapping, it is possible to use the changes in the cortical lamination patterns to define individual cortical areas. This allows in vivo neuroanatomical maps to be generated for individual subjects and precise correlation of the results from functional imaging studies in these subjects with their own microanatomical information. To this end, we adapted the well-established observer-independent cytoarchitectonic mapping techniques for defining cortical borders based on changes in cortical lamination for in vivo parcellation of high-resolution structural MR images.

Improving whole brain structural MRI at 4.7 Tesla using 4 irregularly shaped receiver coils

Both higher magnetic field strengths (> or =3 T) and multiple receiver "array coils" can provide increased signal-to-noise ratio (SNR) for MRI. This increase in SNR can be used to obtain images with higher resolution, enabling better visualisation of structures within the human brain. However, high field strength systems also suffer from increased B(1) non-uniformity and increased power deposition, reaching specific absorption rate (SAR) limits more quickly. For these problems to be mitigated, a careful choice of both the pulse sequence design and transmit RF coil is required. This paper describes the use of a prototype array coil consisting of 4 irregularly shaped coils within a standard configuration for neuroimaging at 4.7 T (a head transmit/receive volume coil to minimise SAR and a head gradient insert for maximum gradient performance). With a fast spin echo (FSE) pulse sequence optimised for 4.7 T, this provides dramatically increased quality and resolution over a large brain volume. Using the array coil, a SNR improvement relative to the volume coil of 1-1.5 times in central brain areas and 2-3 times in cortical regions was obtained. Array coil images with a resolution of 352 x 352 x 2000 mum had a SNR of 16.0 to 26.2 in central regions and 19.9 to 34.8 in cortical areas. Such images easily demonstrate cortical myeloarchitecture, while still covering most of the brain in a approximately 12 min scan.

High-resolution MRI: in vivo histology?

For centuries scientists have been fascinated with the question of how the brain works. Investigators have looked at both where different functions are localized and how the anatomical microstructure varies across the brain surface. Here we discuss how advances in magnetic resonance imaging (MRI) have allowed in vivo visualization of the fine structure of the brain that was previously only visible in post-mortem brains. We present data showing the correspondence between definitions of the primary visual cortex defined anatomically using very high-resolution MRI and functionally using functional MRI. We consider how this technology can be applied to allow the investigation of brains that differ from normal, and what this ever-evolving technology may be able to reveal about in vivo brain structure in the next few years.

Methodological issues relating to in vivo cortical myelography using MRI

The relationship between neocortical structure and function is a key area of research in neuroscience. Most studies of neural function, whether using neurophysiology or neuroimaging methods, are interpreted with relation to the underlying cortical myelo- and cytoarchitecture. For functional neuroimaging studies this often means using cytoarchitectonic maps based on the study of a limited number of brains, despite evidence for substantial interindividual variation. Improvements in MR technology, resulting in wider availability of high-field MRI systems, have led to an increase in the achievable resolution in MR scans. Several groups have reported the in vivo detection of myelination patterns within the cortex, consistent with those observed in postmortem tissue. This leads to the possibility of predefining areas for fMRI analysis based on the cortical architecture. To do this it is essential to know, in a quantitative way, how reliably myeloarchitectonic areas and boundaries can be detected using MRI. Here we investigate the striate cortex, known to be coincident with V1, to assess the detectability of the stria of Gennari across V1 and across subjects. Under optimal conditions, 80% of the stria of Gennari was visualized using our methodology, although there was considerable variability in the level of detection across subjects. We discuss the limitations of the methodology and propose ways to improve the detection level of cortical myeloarchitecture more generally.

High-resolution MRI reflects myeloarchitecture and cytoarchitecture of human cerebral cortex

Maps of cytoarchitectonically defined cortical areas have proven to be a valuable tool for anatomic localization of activated brain regions revealed by functional imaging studies. However, architectonic data require observations in a sample of postmortem brains. They can only be used reliably for comparison with functional data as probabilistic maps after spatial normalization to a common reference space. The complete architectonic analysis of an individual living brain has not been achievable to date, because the relationship remains unclear between laminar gray value changes of cerebral cortex in magnetic resonance (MR) images and those of cyto- and myeloarchitectonic histologic sections. We examined intensity profiles through the cortex in five imaging modalities: in vivo T1 and postmortem T2 MRI, one cell body stain, and two myelin stains. After visualizing the dissimilarities in the shapes of these profiles using a canonical analysis, differences between the profiles from the different image modalities were compared quantitatively. Subsequently, the profiles extracted from the in vivo T1-weighted images were estimated from profiles extracted from cyto- and myeloarchitectonic sections using linear combinations. We could verify statistically the mixed nature of the cortical T1 signal obtained in vivo: The MR intensity profiles were significantly more similar to myeloarchitectonic than to cytoarchitectonic profiles, but a weighted sum of both fitted the T1 profiles best.

Localization of the Human Cortical Visual Area MT Based on Computer Aided Histological Analysis

We describe human area MT histologically based on the observer independent analysis of cortical myeloarchiteture, multiple complementary staining techniques and 3-D reconstruction. The topography of an architectonic field that presented constant structural characteristics across specimens was studied in relation to the sulcal geography of the occipito-temporal region. Objective and semi-automated analysis of local microstructure revealed a distinct cortical architecture and matched topographically the localization of MT derived from functional imaging. MT was localized by the histotopographic method in relation to definite macroscopic landmarks. This study demonstrates a new set of distinguishing architectonic features of human MT that permit localization on structural grounds and suggests that the characteristic laminar structure of this area may be related to its unique pattern of connections and to its role in visual perception.

Changes in connectivity profiles define functionally distinct regions in human medial frontal cortex

A fundamental issue in neuroscience is the relation between structure and function. However, gross landmarks do not correspond well to microstructural borders and cytoarchitecture cannot be visualized in a living brain used for functional studies. Here, we used diffusion-weighted and functional MRI to test structure-function relations directly. Distinct neocortical regions were defined as volumes having similar connectivity profiles and borders identified where connectivity changed. Without using prior information, we found an abrupt profile change where the border between supplementary motor area (SMA) and pre-SMA is expected. Consistent with this anatomical assignment, putative SMA and pre-SMA connected to motor and prefrontal regions, respectively. Excellent spatial correlations were found between volumes defined by using connectivity alone and volumes activated during tasks designed to involve SMA or pre-SMA selectively. This finding demonstrates a strong relationship between structure and function in medial frontal cortex and offers a strategy for testing such correspondences elsewhere in the brain.

A myelo-architectonic method for the structural classification of cortical areas

We describe an automatic and reproducible method to analyze the histological design of the cerebral cortex as applied to brain sections stained to reveal myelinated fibers. The technique provides an evaluation of the distribution of myelination across the width of the cortical mantle in accordance with a model of its curvature and its intrinsic geometry. The profile lines along which the density of staining is measured are generated from the solution of a partial differential equation (PDE) that models the intermediate layers of the cortex. Cortical profiles are classified according to significant components that emerge from wavelet analysis. Intensity profiles belonging to each distinct class are normalized and averaged to produce area-specific templates of cortical myelo-architecture.

Anatomy of the visual system

This brief systemic overview presents the anatomic details of the orbit with respect to imaging modalities CT and MR. The structures of the four orbital compartments, intrakonal and extrakonal space, globe and optic nerve are demonstrated in detail on different CT and MR views (axial, coronal, in soft tissue and bone window, T1-weighted, T2-weighted) with corresponding diagrams. The intracranial visual pathway is explained in detail with emphasis to the striate cortex and extrastriate visual association cortex, presented with diagrams and high-resolution MR.

In vivo identification of human cortical areas using high-resolution MRI: an approach to cerebral structure-function correlation

Understanding the relationship between the structural and functional organization of the human brain is one of the most important goals of neuroscience. Individual variability in brain structure means that it is essential to obtain this information from the same subject. To date, this has been almost impossible. Even though noninvasive functional imaging techniques such as functional MRI (fMRI) are now commonplace, there is no complementary noninvasive structural technique. We present an in vivo method of examining the detailed neuroanatomy of any individual, which can then be correlated with that individual's own functional results. This method utilizes high-resolution structural MRI to identify distinct cortical regions based on cortical lamination structure. We demonstrate that the observed MR lamination patterns relate to myeloarchitecture through a correlation of histology with MRI. In vivo high-resolution MRI studies identify striate cortex, as well as visual area V5, in four individuals, as defined by using fMRI. The anatomical identification of a cortical area (V5MT) outside of striate cortex is a significant advance, proving it possible to identify extra-striate cortical areas and demonstrating that in vivo structural mapping of the human cerebral cortex is possible.

Imaging cortical anatomy by high-resolution MR at 3.0T: detection of the stripe of Gennari in visual area 17

The brain can be parcellated into numerous anatomical and functional subunits. The classic work by Brodmann (Vergleichende Lokalisationslehre der Grosshirnrinde in ihren Prinzipien dargestellt auf Grund des Zellenbaues. Leipzig: Barth; 1909) identified areas of the cerebral cortex based on histological differences. An alternative to his cytoarchitectonic approach is the myeloarchitectonic approach. MRI has excellent white/gray matter contrast in the brain due to the presence of myelin, and thus seems uniquely suited for in vivo studies of cortical myeloarchitecture. Here it is demonstrated that the stripe or stria of Gennari can be consistently detected in human occipital cortex. T(1)-weighted images obtained at 3T from six of 10 normal volunteers, with resolutions of 350 x 350 x 600 mu clearly demonstrate this myelin-rich intracortical layer. It is concluded that the striate cortex (area 17 of Brodmann) of the human brain can be delineated in vivo on T(1)-weighted images, potentially enabling detection of specific cortical boundaries within individual brains.

Ultra high-resolution fMRI in monkeys with implanted RF coils

Spatiotemporally resolved functional MRI (fMRI) in animals can reveal how wide-spread neural networks are organized and accompanying electrophysiological recordings can show how small neural assemblies contribute to this organization. Here we present a novel technique that yields high-resolution structural and functional images of the monkey brain with small, tissue-compatible, intraosteally implantable radiofrequency coils. Voxel sizes as small as 0.0113 microl with high signal-to-noise and contrast-to-noise ratios were obtained, revealing both structural and functional cortical architecture in great detail. Up to a certain point, contrast sensitivity increased with decreasing voxel size, probably because of the decreased partial volume effects. Spatial specificity was demonstrated by the lamina-specific activation in experiments comparing responses to moving and flickering stimuli. The implications of this technique for combined fMRI/electrophysiology experiments and its limitations in terms of spatial coverage are discussed.