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

7T MP2RAGE for cortical myelin segmentation: Impact of aging

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

The prosubiculum in the human hippocampus: A rostrocaudal, feature-driven, and systematic approach

The hippocampal subfield prosubiculum (ProS), is a conserved neuroanatomic region in mouse, monkey, and human. This area lies between CA1 and subiculum (Sub) and particularly lacks consensus on its boundaries; reports have varied on the description of its features and location. In this report, we review, refine, and evaluate four cytoarchitectural features that differentiate ProS from its neighboring subfields: (1) small neurons, (2) lightly stained neurons, (3) superficial clustered neurons, and (4) a cell sparse zone. ProS was delineated in all cases (n = 10). ProS was examined for its cytoarchitectonic features and location rostrocaudally, from the anterior head through the body in the hippocampus. The most common feature was small pyramidal neurons, which were intermingled with larger pyramidal neurons in ProS. We quantitatively measured ProS pyramidal neurons, which showed (average, width at pyramidal base = 14.31 µm, n = 400 per subfield). CA1 neurons averaged 15.57 µm and Sub neurons averaged 15.63 µm, both were significantly different than ProS (Kruskal-Wallis test, p < .0001). The other three features observed were lightly stained neurons, clustered neurons, and a cell sparse zone. Taken together, these findings suggest that ProS is an independent subfield, likely with distinct functional contributions to the broader interconnected hippocampal network. Our results suggest that ProS is a cytoarchitecturally varied subfield, both for features and among individuals. This diverse architecture in features and individuals for ProS could explain the long-standing complexity regarding the identification of this subfield.

Evaluating the effect of denoising submillimeter auditory fMRI data with NORDIC

Functional magnetic resonance imaging (fMRI) has emerged as an essential tool for exploring human brain function. Submillimeter fMRI, in particular, has emerged as a tool to study mesoscopic computations. The inherently low signal-to-noise ratio (SNR) at submillimeter resolutions warrants the use of denoising approaches tailored at reducing thermal noise - the dominant contributing noise component in high resolution fMRI. NORDIC PCA is one of such approaches, and has been benchmarked against other approaches in several applications. Here, we investigate the effects that two versions of NORDIC denoising have on auditory submillimeter data. As investigating auditory functional responses poses unique challenges, we anticipated that the benefit of this technique would be especially pronounced. Our results show that NORDIC denoising improves the detection sensitivity and the reliability of estimates in submillimeter auditory fMRI data. These effects can be explained by the reduction of the noise-induced signal variability. However, we also observed a reduction in the average response amplitude (percent signal), which may suggest that a small amount of signal was also removed. We conclude that, while evaluating the effects of the signal reduction induced by NORDIC may be necessary for each application, using NORDIC in high resolution auditory fMRI studies may be advantageous because of the large reduction in variability of the estimated responses.

Resting state fMRI connectivity is sensitive to laminar connectional architecture in the human brain

Previous invasive studies indicate that human neocortical graymatter contains cytoarchitectonically distinct layers, with notable differences in their structural connectivity with the rest of the brain. Given recent improvements in the spatial resolution of anatomical and functional magnetic resonance imaging (fMRI), we hypothesize that resting state functional connectivity (FC) derived from fMRI is sensitive to layer-specific thalamo-cortical and cortico-cortical microcircuits. Using sub-millimeter resting state fMRI data obtained at 7 T, we found that: (1) FC between the entire thalamus and cortical layers I and VI was significantly stronger than between the thalamus and other layers. Furthermore, FC between somatosensory thalamus (ventral posterolateral nucleus, VPL) and layers IV, VI of the primary somatosensory cortex were stronger than with other layers; (2) Inter-hemispheric cortico-cortical FC between homologous regions in superficial layers (layers I-III) was stronger compared to deep layers (layers V-VI). These findings are in agreement with structural connections inferred from previous invasive studies that showed that: (i) M-type neurons in the entire thalamus project to layer-I; (ii) Pyramidal neurons in layer-VI target all thalamic nuclei, (iii) C-type neurons in the VPL project to layer-IV and receive inputs from layer-VI of the primary somatosensory cortex, and (iv) 80% of collosal projecting neurons between homologous cortical regions connect superficial layers. Our results demonstrate for the first time that resting state fMRI is sensitive to structural connections between cortical layers (previously inferred through invasive studies), specifically in thalamo-cortical and cortico-cortical networks.

A sub+cortical fMRI-based surface parcellation

Both cortical and subcortical structures are organized into a large number of distinct areas reflecting functional and cytoarchitectonic differences. Mapping these areas is of fundamental importance to neuroscience. A central obstacle to this task is the inaccuracy associated with bringing results from individuals into a common space. The vast individual differences in morphology pose a serious problem for volumetric registration. Surface‐based approaches fare substantially better, but have thus far been used only for cortical parcellation, leaving subcortical parcellation in volumetric space. We extend the surface‐based approach to include also the subcortical deep gray‐matter structures, thus achieving a uniform representation across both cortex and subcortex, suitable for use with surface‐based metrics that span these structures, for example, white/gray contrast. Using data from the Enhanced Nathan Klein Institute—Rockland Sample, limited to individuals between 19 and 69 years of age, we generate a functional parcellation of both the cortical and subcortical surfaces. To assess this extended parcellation, we show that (a) our parcellation provides greater homogeneity of functional connectivity patterns than do arbitrary parcellations matching in the number and size of parcels; (b) our parcels align with known cortical and subcortical architecture; and (c) our extended functional parcellation provides an improved fit to the complexity of life‐span (6–85 years) changes in white/gray contrast data compared to arbitrary parcellations matching in the number and size of parcels, supporting its use with surface‐based measures. We provide our extended functional parcellation for the use of the neuroimaging community.

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.

Surface-based analysis increases the specificity of cortical activation patterns and connectivity results

Spatial smoothing of functional magnetic resonance imaging (fMRI) data can be performed on volumetric images and on the extracted surface of the brain. Smoothing on the unfolded cortex should theoretically improve the ability to separate signals between brain areas that are near together in the folded cortex but are more distant in the unfolded cortex. However, surface-based method approaches (SBA) are currently not utilized as standard procedure in the preprocessing of neuroimaging data. Recent improvements in the quality of cortical surface modeling and improvements in its usability nevertheless advocate this method. In the current study, we evaluated the benefits of an up-to-date surface-based smoothing in comparison to volume-based smoothing. We focused on the effect of signal contamination between different functional systems using the primary motor and primary somatosensory cortex as an example. We were particularly interested in how this signal contamination influences the results of activity and connectivity analyses for these brain regions. We addressed this question by performing fMRI on 19 subjects during a tactile stimulation paradigm and by using simulated BOLD responses. We demonstrated that volume-based smoothing causes contamination of the primary motor cortex by somatosensory cortical responses, leading to false positive motor activation. These false positive motor activations were not found by using surface-based smoothing for reasonable kernel sizes. Accordingly, volume-based smoothing caused an exaggeration of connectivity estimates between these regions. In conclusion, this study showed that surface-based smoothing decreases signal contamination considerably between neighboring functional brain regions and improves the validity of activity and connectivity results.

Resting-State Functional Network Organization Is Stable Across Adolescent Development for Typical and Psychosis Spectrum Youth

BACKGROUND

Resting-state functional neuroimaging captures large-scale network organization; whether this organization is intact or disrupted during adolescent development across the psychosis spectrum is unresolved. We investigated the integrity of network organization in psychosis spectrum youth and those with first episode psychosis (FEP) from late childhood through adulthood.

METHODS

We analyzed data from Philadelphia Neurodevelopmental Cohort (PNC; typically developing = 450, psychosis spectrum = 273, 8-22 years), a longitudinal cohort of typically developing youth (LUNA; N = 208, 1-3 visits, 10-25 years), and a sample of FEP (N = 39) and matched controls (N = 34). We extracted individual time series and calculated correlations from brain regions and averaged them for 4 age groups: late childhood, early adolescence, late adolescence, adulthood. Using multiple analytic approaches, we assessed network stability across 4 age groups, compared stability between controls and psychosis spectrum youth, and compared group-level network organization of FEP to controls. We explored whether variability in cognition or clinical symptomatology was related to network organization.

RESULTS

Network organization was stable across the 4 age groups in the PNC and LUNA typically developing youth and PNC psychosis spectrum youth. Psychosis spectrum and typically developing youth had similar functional network organization during all age ranges. Network organization was intact in PNC youth who met full criteria for psychosis and in FEP. Variability in cognitive functioning or clinical symptomatology was not related to network organization in psychosis spectrum youth or FEP.

DISCUSSION

These findings provide rigorous evidence supporting intact functional network organization in psychosis risk and psychosis from late childhood through adulthood.

A likelihood ratio approach for functional localization in fMRI

BACKGROUND

To increase power when analyzing fMRI data, researchers often define functional regions of interest (fROIs). It is crucial that this fROI is defined with an optimal balance between both false positives and false negatives to ensure maximal spatial accuracy and to avoid potentially biased results in the main fMRI experiment. Additionally, since the fROI is defined in each subject separately, the used method should attune to the general level of activation of the individual.

NEW METHOD

We investigate the benefits of the maximized likelihood ratio (mLR) method. This method is based on the likelihood paradigm where likelihood ratios are used to reflect relative statistical evidence in favor of an a priori defined practically relevant alternative hypothesis as compared to the null hypothesis of no activation.

RESULTS

Through both simulations and real data, we show that the mLR method provides cumulative evidence for voxels that are active with an effect size that is larger than the one a priori defined in the alternative. Furthermore, an optimal balance between Type I and Type II errors is achieved when the alternative is an underestimation of the true effect size.

COMPARISON WITH EXISTING METHODS

The mLR method is compared with false discovery rate corrected null hypothesis significance testing and regular likelihood ratio testing. It performs as good as or outperformed both methods in both detection of practically relevant voxels and the trade- off between false positives and false negatives.

CONCLUSIONS

The mLR method provides fROIs that are both spatially accurate and practically relevant.

Educational fMRI: From the Lab to the Classroom

Functional MRI (fMRI) findings hold many potential applications for education, and yet, the translation of fMRI findings to education has not flowed. Here, we address the types of fMRI that could better support applications of neuroscience to the classroom. This 'educational fMRI' comprises eight main challenges: (1) collecting artifact-free fMRI data in school-aged participants and in vulnerable young populations, (2) investigating heterogenous cohorts with wide variability in learning abilities and disabilities, (3) studying the brain under natural and ecological conditions, given that many practical topics of interest for education can be addressed only in ecological contexts, (4) depicting complex age-dependent associations of brain and behaviour with multi-modal imaging, (5) assessing changes in brain function related to developmental trajectories and instructional intervention with longitudinal designs, (6) providing system-level mechanistic explanations of brain function, so that useful individualized predictions about learning can be generated, (7) reporting negative findings, so that resources are not wasted on developing ineffective interventions, and (8) sharing data and creating large-scale longitudinal data repositories to ensure transparency and reproducibility of fMRI findings for education. These issues are of paramount importance to the development of optimal fMRI practices for educational applications.

Myelin and Modeling: Bootstrapping Cortical Microcircuits

Histological studies of myelin-stained sectioned cadaver brain and in vivo myelin-weighted magnetic resonance imaging (MRI) show that the cerebral cortex is organized into cortical areas with generally well-defined boundaries, which have consistent internal patterns of myelination. The process of myelination is largely driven by neural experience, in which the axonal passage of action potentials stimulates neighboring oligodendrocytes to perform their task. This bootstrapping process, such that the traffic of action potentials facilitates increased traffic, suggests the hypothesis that the specific pattern of myelination (myeloarchitecture) in each cortical area reveals the principal cortical microcircuits required for the function of that area. If this idea is correct, the observable sequential maturation of specific brain areas can provide evidence for models of the stages of cognitive development.

Size and shape matter: The impact of voxel geometry on the identification of small nuclei

How, and to what extent do size and shape of a voxel measured with magnetic resonance imaging (MRI) affect the ability to visualize small brain nuclei? Despite general consensus that voxel geometry affects volumetric properties of regions of interest, particularly those of small brain nuclei, no quantitative data on the influence of voxel size and shape on labeling accuracy is available. Using simulations, we investigated the selective influence of voxel geometry by reconstructing simulated ellipsoid structures with voxels varying in shape and size. For each reconstructed ellipsoid, we calculated differences in volume and similarity between the labeled volume and the predefined dimensions of the ellipsoid. Probability functions were derived from one or two individual raters and a simulated ground truth for reference. As expected, larger voxels (i.e., coarser resolution) and increasing anisotropy results in increased deviations of both volume and shape measures, which is of particular relevance for small brain structures. Our findings clearly illustrate the anatomical inaccuracies introduced by the application of large and/or anisotropic voxels. To ensure deviations occur within the acceptable range (Dice coefficient scores; DCS > 0.75, corresponding to < 57% volume deviation), the volume of isotropic voxels should not exceed 5% of the total volume of the region of interest. When high accuracy is required (DCS > 0.90, corresponding to a < 19% volume deviation), the volumes of isotropic voxels should not exceed 0.08%, of the total volume. Finally, when large anisotropic factors (>3) are used, and the ellipsoid is orthogonal to the slice axes, having its long axis in the imaging plane, the voxel volume should not exceed 0.005% of the total volume. This allows sufficient compensation of anisotropy effects, in order to reach accuracy in the acceptable range (DCS > 0.75, corresponding to >57% volume deviation).

Addressing challenges of high spatial resolution UHF fMRI for group analysis of higher-order cognitive tasks: An inter-sensory task directing attention between visual and somatosensory domains

Functional MRI at ultra‐high field (UHF, ≥7 T) provides significant increases in BOLD contrast‐to‐noise ratio (CNR) compared with conventional field strength (3 T), and has been exploited for reduced field‐of‐view, high spatial resolution mapping of primary sensory areas. Applying these high spatial resolution methods to investigate whole brain functional responses to higher‐order cognitive tasks leads to a number of challenges, in particular how to perform robust group‐level statistical analyses. This study addresses these challenges using an inter‐sensory cognitive task which modulates top‐down attention at graded levels between the visual and somatosensory domains. At the individual level, highly focal functional activation to the task and task difficulty (modulated by attention levels) were detectable due to the high CNR at UHF. However, to assess group level effects, both anatomical and functional variability must be considered during analysis. We demonstrate the importance of surface over volume normalisation and the requirement of no spatial smoothing when assessing highly focal activity. Using novel group analysis on anatomically parcellated brain regions, we show that in higher cognitive areas (parietal and dorsal‐lateral‐prefrontal cortex) fMRI responses to graded attention levels were modulated quadratically, whilst in visual cortex and VIP, responses were modulated linearly. These group fMRI responses were not seen clearly using conventional second‐level GLM analyses, illustrating the limitations of a conventional approach when investigating such focal responses in higher cognitive regions which are more anatomically variable. The approaches demonstrated here complement other advanced analysis methods such as multivariate pattern analysis, allowing UHF to be fully exploited in cognitive neuroscience.

LISA improves statistical analysis for fMRI

One of the principal goals in functional magnetic resonance imaging (fMRI) is the detection of local activation in the human brain. However, lack of statistical power and inflated false positive rates have recently been identified as major problems in this regard. Here, we propose a non-parametric and threshold-free framework called LISA to address this demand. It uses a non-linear filter for incorporating spatial context without sacrificing spatial precision. Multiple comparison correction is achieved by controlling the false discovery rate in the filtered maps. Compared to widely used other methods, it shows a boost in statistical power and allows to find small activation areas that have previously evaded detection. The spatial sensitivity of LISA makes it especially suitable for the analysis of high-resolution fMRI data acquired at ultrahigh field (≥7 Tesla).

Visualizing the Human Subcortex Using Ultra-high Field Magnetic Resonance Imaging

With the recent increased availability of ultra-high field (UHF) magnetic resonance imaging (MRI), substantial progress has been made in visualizing the human brain, which can now be done in extraordinary detail. This review provides an extensive overview of the use of UHF MRI in visualizing the human subcortex for both healthy and patient populations. The high inter-subject variability in size and location of subcortical structures limits the usability of atlases in the midbrain. Fortunately, the combined results of this review indicate that a large number of subcortical areas can be visualized in individual space using UHF MRI. Current limitations and potential solutions of UHF MRI for visualizing the subcortex are also discussed.

The impact of traditional neuroimaging methods on the spatial localization of cortical areas

Localizing human brain functions is a long-standing goal in systems neuroscience. Toward this goal, neuroimaging studies have traditionally used volume-based smoothing, registered data to volume-based standard spaces, and reported results relative to volume-based parcellations. A novel 360-area surface-based cortical parcellation was recently generated using multimodal data from the Human Connectome Project, and a volume-based version of this parcellation has frequently been requested for use with traditional volume-based analyses. However, given the major methodological differences between traditional volumetric and Human Connectome Project-style processing, the utility and interpretability of such an altered parcellation must first be established. By starting from automatically generated individual-subject parcellations and processing them with different methodological approaches, we show that traditional processing steps, especially volume-based smoothing and registration, substantially degrade cortical area localization compared with surface-based approaches. We also show that surface-based registration using features closely tied to cortical areas, rather than to folding patterns alone, improves the alignment of areas, and that the benefits of high-resolution acquisitions are largely unexploited by traditional volume-based methods. Quantitatively, we show that the most common version of the traditional approach has spatial localization that is only 35% as good as the best surface-based method as assessed using two objective measures (peak areal probabilities and "captured area fraction" for maximum probability maps). Finally, we demonstrate that substantial challenges exist when attempting to accurately represent volume-based group analysis results on the surface, which has important implications for the interpretability of studies, both past and future, that use these volume-based methods.

Statistical power comparisons at 3T and 7T with a GO / NOGO task

The field of cognitive neuroscience is weighing evidence about whether to move from standard field strength to ultra-high field (UHF). The present study contributes to the evidence by comparing a cognitive neuroscience paradigm at 3 Tesla (3T) and 7 Tesla (7T). The goal was to test and demonstrate the practical effects of field strength on a standard GO/NOGO task using accessible preprocessing and analysis tools. Two independent matched healthy samples (N = 31 each) were analyzed at 3T and 7T. Results show gains at 7T in statistical strength, the detection of smaller effects and group-level power. With an increased availability of UHF scanners, these gains may be exploited by cognitive neuroscientists and other neuroimaging researchers to develop more efficient or comprehensive experimental designs and, given the same sample size, achieve greater statistical power at 7T.

Intrinsic functional connectivity of the central extended amygdala

The central extended amygdala (EAc)-including the bed nucleus of the stria terminalis (BST) and central nucleus of the amygdala (Ce)-plays a critical role in triggering fear and anxiety and is implicated in the development of a range of debilitating neuropsychiatric disorders. Although it is widely believed that these disorders reflect the coordinated activity of distributed neural circuits, the functional architecture of the EAc network and the degree to which the BST and the Ce show distinct patterns of functional connectivity is unclear. Here, we used a novel combination of imaging approaches to trace the connectivity of the BST and the Ce in 130 healthy, racially diverse, community-dwelling adults. Multiband imaging, high-precision registration techniques, and spatially unsmoothed data maximized anatomical specificity. Using newly developed seed regions, whole-brain regression analyses revealed robust functional connectivity between the BST and Ce via the sublenticular extended amygdala, the ribbon of subcortical gray matter encompassing the ventral amygdalofugal pathway. Both regions displayed coupling with the ventromedial prefrontal cortex (vmPFC), midcingulate cortex (MCC), insula, and anterior hippocampus. The BST showed stronger connectivity with the thalamus, striatum, periaqueductal gray, and several prefrontal territories. The only regions showing stronger functional connectivity with the Ce were neighboring regions of the dorsal amygdala, amygdalohippocampal area, and anterior hippocampus. These observations provide a baseline against which to compare a range of special populations, inform our understanding of the role of the EAc in normal and pathological fear and anxiety, and showcase image registration techniques that are likely to be useful for researchers working with "deidentified" neuroimaging data.

Cortical parcellation based on structural connectivity: A case for generative models

One of the major challenges in systems neuroscience is to identify brain networks and unravel their significance for brain function -this has led to the concept of the 'connectome'. Connectomes are currently extensively studied in large-scale international efforts at multiple scales, and follow different definitions with respect to their connections as well as their elements. Perhaps the most promising avenue for defining the elements of connectomes originates from the notion that individual brain areas maintain distinct (long-range) connection profiles. These connectivity patterns determine the areas' functional properties and also allow for their anatomical delineation and mapping. This rationale has motivated the concept of connectivity-based cortex parcellation. In the past ten years, non-invasive mapping of human brain connectivity has led to immense advances in the development of parcellation techniques and their applications. Unfortunately, many of these approaches primarily aim for confirmation of well-known, existing architectonic maps and, to that end, unsuitably incorporate prior knowledge and frequently build on circular argumentation. Often, current approaches also tend to disregard the specific apertures of connectivity measurements, as well as the anatomical specificities of cortical areas, such as spatial compactness, regional heterogeneity, inter-subject variability, the multi-scaling nature of connectivity information, and potential hierarchical organisation. From a methodological perspective, however, a useful framework that regards all of these aspects in an unbiased way is technically demanding. In this commentary, we first outline the concept of connectivity-based cortex parcellation and discuss its prospects and limitations in particular with respect to structural connectivity. To improve reliability and efficiency, we then strongly advocate for connectivity-based cortex parcellation as a modelling approach; that is, an approximation of the data based on (model) parameter inference. As such, a parcellation algorithm can be formally tested for robustness -the precision of its predictions can be quantified and statistics about potential generalization of the results can be derived. Such a framework also allows the question of model constraints to be reformulated in terms of hypothesis testing through model selection and offers a formative way to integrate anatomical knowledge in terms of prior distributions.

Comparing functional MRI protocols for small, iron-rich basal ganglia nuclei such as the subthalamic nucleus at 7 T and 3 T

The basal ganglia (BG) form a network of subcortical nuclei. Functional magnetic resonance imaging (fMRI) in the BG could provide insight in its functioning and the underlying mechanisms of Deep Brain Stimulation (DBS). However, fMRI of the BG with high specificity is challenging, because the nuclei are small and variable in their anatomical location. High resolution fMRI at field strengths of 7 Tesla (T) could help resolve these challenges to some extent. A set of MR protocols was developed for functional imaging of the BG nuclei at 3 T and 7 T. The protocols were validated using a stop-signal reaction task (Logan et al. []: J Exp Psychol: Human Percept Perform 10:276-291). Compared with sub-millimeter 7 T fMRI protocols aimed at cortex, a reduction of echo time and spatial resolution was strictly necessary to obtain robust Blood Oxygen Level Dependent (BOLD) sensitivity in the BG. An fMRI protocol at 3 T with identical resolution to the 7 T showed no robust BOLD sensitivity in any of the BG nuclei. The results suggest that the subthalamic nucleus, as well as the substantia nigra, red nucleus, and the internal and external parts of the globus pallidus show increased activation in failed stop trials compared with successful stop and go trials. Hum Brain Mapp 38:3226-3248, 2017. © 2017 Wiley Periodicals, Inc.

Introducing Alternative-Based Thresholding for Defining Functional Regions of Interest in fMRI

In fMRI research, one often aims to examine activation in specific functional regions of interest (fROIs). Current statistical methods tend to localize fROIs inconsistently, focusing on avoiding detection of false activation. Not missing true activation is however equally important in this context. In this study, we explored the potential of an alternative-based thresholding (ABT) procedure, where evidence against the null hypothesis of no effect and evidence against a prespecified alternative hypothesis is measured to control both false positives and false negatives directly. The procedure was validated in the context of localizer tasks on simulated brain images and using a real data set of 100 runs per subject. Voxels categorized as active with ABT can be confidently included in the definition of the fROI, while inactive voxels can be confidently excluded. Additionally, the ABT method complements classic null hypothesis significance testing with valuable information by making a distinction between voxels that show evidence against both the null and alternative and voxels for which the alternative hypothesis cannot be rejected despite lack of evidence against the null.

Uses, misuses, new uses and fundamental limitations of magnetic resonance imaging in cognitive science

When blood oxygenation level-dependent (BOLD) contrast functional magnetic resonance imaging (fMRI) was discovered in the early 1990s, it provoked an explosion of interest in exploring human cognition, using brain mapping techniques based on MRI. Standards for data acquisition and analysis were rapidly put in place, in order to assist comparison of results across laboratories. Recently, MRI data acquisition capabilities have improved dramatically, inviting a rethink of strategies for relating functional brain activity at the systems level with its neuronal substrates and functional connections. This paper reviews the established capabilities of BOLD contrast fMRI, the perceived weaknesses of major methods of analysis, and current results that may provide insights into improved brain modelling. These results have inspired the use of in vivo myeloarchitecture for localizing brain activity, individual subject analysis without spatial smoothing and mapping of changes in cerebral blood volume instead of BOLD activation changes. The apparent fundamental limitations of all methods based on nuclear magnetic resonance are also discussed.This article is part of the themed issue 'Interpreting BOLD: a dialogue between cognitive and cellular neuroscience'.

The technological future of 7 T MRI hardware

In this article we present our projections of future hardware developments on 7 T human MRI systems. These include compact cryogen-light magnets, improved gradient performance, integrated RF-receive and direct current shimming coil arrays, new RF technology with adaptive impedance matching, patient-specific specific absorption rate estimation and monitoring, and increased integration of physiological monitoring systems. Copyright © 2015 John Wiley & Sons, Ltd.

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.

The Human Connectome Project's neuroimaging approach

Noninvasive human neuroimaging has yielded many discoveries about the brain. Numerous methodological advances have also occurred, though inertia has slowed their adoption. This paper presents an integrated approach to data acquisition, analysis and sharing that builds upon recent advances, particularly from the Human Connectome Project (HCP). The 'HCP-style' paradigm has seven core tenets: (i) collect multimodal imaging data from many subjects; (ii) acquire data at high spatial and temporal resolution; (iii) preprocess data to minimize distortions, blurring and temporal artifacts; (iv) represent data using the natural geometry of cortical and subcortical structures; (v) accurately align corresponding brain areas across subjects and studies; (vi) analyze data using neurobiologically accurate brain parcellations; and (vii) share published data via user-friendly databases. We illustrate the HCP-style paradigm using existing HCP data sets and provide guidance for future research. Widespread adoption of this paradigm should accelerate progress in understanding the brain in health and disease.

Task-Related Edge Density (TED)-A New Method for Revealing Dynamic Network Formation in fMRI Data of the Human Brain

The formation of transient networks in response to external stimuli or as a reflection of internal cognitive processes is a hallmark of human brain function. However, its identification in fMRI data of the human brain is notoriously difficult. Here we propose a new method of fMRI data analysis that tackles this problem by considering large-scale, task-related synchronisation networks. Networks consist of nodes and edges connecting them, where nodes correspond to voxels in fMRI data, and the weight of an edge is determined via task-related changes in dynamic synchronisation between their respective times series. Based on these definitions, we developed a new data analysis algorithm that identifies edges that show differing levels of synchrony between two distinct task conditions and that occur in dense packs with similar characteristics. Hence, we call this approach "Task-related Edge Density" (TED). TED proved to be a very strong marker for dynamic network formation that easily lends itself to statistical analysis using large scale statistical inference. A major advantage of TED compared to other methods is that it does not depend on any specific hemodynamic response model, and it also does not require a presegmentation of the data for dimensionality reduction as it can handle large networks consisting of tens of thousands of voxels. We applied TED to fMRI data of a fingertapping and an emotion processing task provided by the Human Connectome Project. TED revealed network-based involvement of a large number of brain areas that evaded detection using traditional GLM-based analysis. We show that our proposed method provides an entirely new window into the immense complexity of human brain function.

A map of the human neocortex showing the estimated overall myelin content of the individual architectonic areas based on the studies of Adolf Hopf

During the period extending from 1910 to 1970, Oscar and Cécile Vogt and their numerous collaborators published a large number of myeloarchitectonic studies on the cortex of the various lobes of the human cerebrum. In a previous publication [Nieuwenhuys et al (Brain Struct Funct 220:2551–2573, 2015; Erratum in Brain Struct Funct 220: 3753–3755, 2015)], we used the data provided by the Vogt–Vogt school for the composition of a myeloarchitectonic map of the entire human neocortex. Because these data were derived from many different brains, a standard brain had to be introduced to which all data available could be transferred. As such the Colin 27 structural scan, aligned to the MNI305 template was selected. The resultant map includes 180 myeloarchitectonic areas, 64 frontal, 30 parietal, 6 insular, 17 occipital and 63 temporal. Here we present a supplementary map in which the overall density of the myelinated fibers in the individual architectonic areas is indicated, based on a meta-analysis of data provided by Adolf Hopf, a prominent collaborator of the Vogts. This map shows that the primary sensory and motor regions are densely myelinated and that, in general, myelination decreases stepwise with the distance from these primary regions. The map also reveals the presence of a number of heavily myelinated formations, situated beyond the primary sensory and motor domains, each consisting of two or more myeloarchitectonic areas. These formations were provisionally designated as the orbitofrontal, intraparietal, posterolateral temporal, and basal temporal dark clusters. Recently published MRI-based in vivo myelin content mappings show, with regard to the primary sensory and motor regions, a striking concordance with our map. As regards the heavily myelinated clusters shown by our map, scrutiny of the current literature revealed that correlates of all of these clusters have been identified in in vivo structural MRI studies and appear to correspond either entirely or largely to known cytoarchitectonic entities. Moreover, functional neuroimaging studies indicate that all of these clusters are involved in vision-related cognitive functions.

Resting-state functional connectivity of the human habenula in healthy individuals: Associations with subclinical depression

INTRODUCTION

The habenula (Hb) is postulated to play a critical role in reward and aversion processing across species, including humans, and has been increasingly implicated in depression. However, technical constraints have limited in vivo investigation of the human Hb, and its function remains poorly characterized. We sought to overcome these challenges by examining the whole-brain resting-state functional connectivity of the Hb and its possible relationship to depressive symptomatology using the high-resolution WU-Minn Human Connectome Project (HCP) dataset.

METHODS

Anatomical and resting-state functional MRI data from 50 healthy subjects with low or high subclinical depression scores (n = 25 each) were analyzed. Using novel semi-automated segmentation and optimization techniques, we generated individual-specific Hb seeds and calculated whole-brain functional connectivity for the entire cohort and the contrast of high vs. low depression groups.

RESULTS

In the entire cohort, the Hb exhibited significant connectivity with key brainstem structures (i.e., ventral tegmental area, substantia nigra, pons) as well as the anterior and posterior cingulate cortices, precuneus, thalamus, and sensorimotor cortex. Multiple regions showed differential Hb connectivity based on subclinical depression scores, including the amygdala, insula, and prefrontal, mid-cingulate, and entorhinal cortices.

CONCLUSIONS

Hb connectivity findings converged on areas associated with salience processing, sensorimotor systems, and the default mode network. We also detected substantial Hb-brainstem connectivity, consistent with prior histological and animal research. High and low subclinical depression groups exhibited differences in Hb connectivity with multiple regions previously linked to depression, suggesting the relationship between these structures as a potential target for future research and treatment. Hum Brain Mapp 37:2369-2384, 2016. © 2016 Wiley Periodicals, Inc.

Neurobiological origin of spurious brain morphological changes: A quantitative MRI study

The high gray‐white matter contrast and spatial resolution provided by T1‐weighted magnetic resonance imaging (MRI) has made it a widely used imaging protocol for computational anatomy studies of the brain. While the image intensity in T1‐weighted images is predominantly driven by T1, other MRI parameters affect the image contrast, and hence brain morphological measures derived from the data. Because MRI parameters are correlates of different histological properties of brain tissue, this mixed contribution hampers the neurobiological interpretation of morphometry findings, an issue which remains largely ignored in the community. We acquired quantitative maps of the MRI parameters that determine signal intensities in T1‐weighted images (R₁ (=1/T1), R₂*, and PD) in a large cohort of healthy subjects (n = 120, aged 18–87 years). Synthetic T1‐weighted images were calculated from these quantitative maps and used to extract morphometry features—gray matter volume and cortical thickness. We observed significant variations in morphometry measures obtained from synthetic images derived from different subsets of MRI parameters. We also detected a modulation of these variations by age. Our findings highlight the impact of microstructural properties of brain tissue—myelination, iron, and water content—on automated measures of brain morphology and show that microstructural tissue changes might lead to the detection of spurious morphological changes in computational anatomy studies. They motivate a review of previous morphological results obtained from standard anatomical MRI images and highlight the value of quantitative MRI data for the inference of microscopic tissue changes in the healthy and diseased brain. Hum Brain Mapp 37:1801–1815, 2016. © 2016 The Authors. Human Brain Mapping Published by Wiley Periodicals, Inc.