High resolution atlas of the venous brain vasculature from 7 T quantitative susceptibility maps

Widespread Autonomic Physiological Coupling Across the Brain-Body Axis

The brain is closely attuned to visceral signals from the body's internal environment, as evidenced by the numerous associations between neural, hemodynamic, and peripheral physiological signals. We show that these brain-body co-fluctuations can be captured by a single spatiotemporal pattern. Across several independent samples, as well as single-echo and multi-echo fMRI data acquisition sequences, we identify widespread co-fluctuations in the low-frequency range (0.01 - 0.1 Hz) between resting-state global fMRI signals, neural activity, and a host of autonomic signals spanning cardiovascular, pulmonary, exocrine and smooth muscle systems. The same brain-body co-fluctuations observed at rest are elicited by arousal induced by cued deep breathing and intermittent sensory stimuli, as well as spontaneous phasic EEG events during sleep. Further, we show that the spatial structure of global fMRI signals is maintained under experimental suppression of end-tidal carbon dioxide (PETCO2) variations, suggesting that respiratory-driven fluctuations in arterial CO2 accompanying arousal cannot explain the origin of these signals in the brain. These findings establish the global fMRI signal as a significant component of the arousal response governed by the autonomic nervous system.

Oxygen and the Spark of Human Brain Evolution: Complex Interactions of Metabolism and Cortical Expansion across Development and Evolution

Scientific theories on the functioning and dysfunction of the human brain require an understanding of its development-before and after birth and through maturation to adulthood-and its evolution. Here we bring together several accounts of human brain evolution by focusing on the central role of oxygen and brain metabolism. We argue that evolutionary expansion of human transmodal association cortices exceeded the capacity of oxygen delivery by the vascular system, which led these brain tissues to rely on nonoxidative glycolysis for additional energy supply. We draw a link between the resulting lower oxygen tension and its effect on cytoarchitecture, which we posit as a key driver of genetic developmental programs for the human brain-favoring lower intracortical myelination and the presence of biosynthetic materials for synapse turnover. Across biological and temporal scales, this protracted capacity for neural plasticity sets the conditions for cognitive flexibility and ongoing learning, supporting complex group dynamics and intergenerational learning that in turn enabled improved nutrition to fuel the metabolic costs of further cortical expansion. Our proposed model delineates explicit mechanistic links among metabolism, molecular and cellular brain heterogeneity, and behavior, which may lead toward a clearer understanding of brain development and its disorders.

The molecular genetic landscape of human brain size variation

Human brain size changes dynamically through early development, peaks in adolescence, and varies up to 2-fold among adults. However, the molecular genetic underpinnings of interindividual variation in brain size remain unknown. Here, we leveraged postmortem brain RNA sequencing and measurements of brain weight (BW) in 2,531 individuals across three independent datasets to identify 928 genome-wide significant associations with BW. Genes associated with higher or lower BW showed distinct neurodevelopmental trajectories and spatial patterns that mapped onto functional and cellular axes of brain organization. Expression of BW genes was predictive of interspecies differences in brain size, and bioinformatic annotation revealed enrichment for neurogenesis and cell-cell communication. Genome-wide, transcriptome-wide, and phenome-wide association analyses linked BW gene sets to neuroimaging measurements of brain size and brain-related clinical traits. Cumulatively, these results represent a major step toward delineating the molecular pathways underlying human brain size variation in health and disease.

Modeling venous bias in resting state functional MRI metrics

Resting-state (rs) functional magnetic resonance imaging (fMRI) is used to detect low-frequency fluctuations in the blood oxygen-level dependent (BOLD) signal across brain regions. Correlations between temporal BOLD signal fluctuations are commonly used to infer functional connectivity. However, because BOLD is based on the dilution of deoxyhemoglobin, it is sensitive to veins of all sizes, and its amplitude is biased by draining veins. These biases affect local BOLD signal location and amplitude, and may also influence BOLD-derived connectivity measures, but the magnitude of this venous bias and its relation to vein size and proximity is unknown. Here, veins were identified using high-resolution quantitative susceptibility maps and utilized in a biophysical model to investigate systematic venous biases on common local rsfMRI-derived measures. Specifically, we studied the impact of vein diameter and distance to veins on the amplitude of low-frequency fluctuations (ALFF), fractional ALFF (fALFF), Hurst exponent (HE), regional homogeneity (ReHo), and eigenvector centrality values in the grey matter. Values were higher across all distances in smaller veins, and decreased with increasing vein diameter. Additionally, rsfMRI values associated with larger veins decrease with increasing distance from the veins. ALFF and ReHo were the most biased by veins, while HE and fALFF exhibited the smallest bias. Across all metrics, the amplitude of the bias was limited in voxel-wise data, confirming that venous structure is not the dominant source of contrast in these rsfMRI metrics. Finally, the models presented can be used to correct this venous bias in rsfMRI metrics.

Cerebral oxygen metabolism from MRI susceptibility

This article provides an overview of MRI methods exploiting magnetic susceptibility properties of blood to assess cerebral oxygen metabolism, including the tissue oxygen extraction fraction (OEF) and the cerebral metabolic rate of oxygen (CMRO2). The first section is devoted to describing blood magnetic susceptibility and its effect on the MRI signal. Blood circulating in the vasculature can have diamagnetic (oxyhemoglobin) or paramagnetic properties (deoxyhemoglobin). The overall balance between oxygenated and deoxygenated hemoglobin determines the induced magnetic field which, in turn, modulates the transverse relaxation decay of the MRI signal via additional phase accumulation. The following sections of this review then illustrate the principles underpinning susceptibility-based techniques for quantifying OEF and CMRO2. Here, it is detailed whether these techniques provide global (OxFlow) or local (Quantitative Susceptibility Mapping - QSM, calibrated BOLD - cBOLD, quantitative BOLD - qBOLD, QSM+qBOLD) measurements of OEF or CMRO2, and what signal components (magnitude or phase) and tissue pools they consider (intravascular or extravascular). Validations studies and potential limitations of each method are also described. The latter include (but are not limited to) challenges in the experimental setup, the accuracy of signal modeling, and assumptions on the measured signal. The last section outlines the clinical uses of these techniques in healthy aging and neurodegenerative diseases and contextualizes these reports relative to results from gold-standard PET.

Vessel Distance Mapping: A novel methodology for assessing vascular-induced cognitive resilience

The association between cerebral blood supply and cognition has been widely discussed in the recent literature. One focus of this discussion has been the anatomical variability of the circle of Willis, with morphological differences being present in more than half of the general population. While previous studies have attempted to classify these differences and explore their contribution to hippocampal blood supply and cognition, results have been controversial. To disentangle these previously inconsistent findings, we introduce Vessel Distance Mapping (VDM) as a novel methodology for evaluating blood supply, which allows for obtaining vessel pattern metrics with respect to the surrounding structures, extending the previously established binary classification into a continuous spectrum. To accomplish this, we manually segmented hippocampal vessels obtained from high-resolution 7T time-of-flight MR angiographic imaging in older adults with and without cerebral small vessel disease, generating vessel distance maps by computing the distances of each voxel to its nearest vessel. Greater values of VDM-metrics, which reflected higher vessel distances, were associated with poorer cognitive outcomes in subjects affected by vascular pathology, while this relation was not observed in healthy controls. Therefore, a mixed contribution of vessel pattern and vessel density is proposed to confer cognitive resilience, consistent with previous research findings. In conclusion, VDM provides a novel platform, based on a statistically robust and quantitative method of vascular mapping, for addressing a variety of clinical research questions.

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

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

Assessment of perivascular space filtering methods using a three-dimensional computational model

Growing interest surrounds the assessment of perivascular spaces (PVS) on magnetic resonance imaging (MRI) and their validation as a clinical biomarker of adverse brain health. Nonetheless, the limits of validity of current state-of-the-art segmentation methods are still unclear. Here, we propose an open-source three-dimensional computational framework comprising 3D digital reference objects and evaluate the performance of three PVS filtering methods under various spatiotemporal imaging considerations (including sampling, motion artefacts, and Rician noise). Specifically, we study the performance of the Frangi, Jerman and RORPO filters in enhancing PVS-like structures to facilitate segmentation. Our findings were three-fold. First, as long as voxels are isotropic, RORPO outperforms the other two filters, regardless of imaging quality. Unlike the Frangi and Jerman filters, RORPO's performance does not deteriorate as PVS volume increases. Second, the performance of all "vesselness" filters is heavily influenced by imaging quality, with sampling and motion artefacts being the most damaging for these types of analyses. Third, none of the filters can distinguish PVS from other hyperintense structures (e.g. white matter hyperintensities, stroke lesions, or lacunes) effectively, the area under precision-recall curve dropped substantially (Frangi: from 94.21 [IQR 91.60, 96.16] to 43.76 [IQR 25.19, 63.38]; Jerman: from 94.51 [IQR 91.90, 95.37] to 58.00 [IQR 35.68, 64.87]; RORPO: from 98.72 [IQR 95.37, 98.96] to 71.87 [IQR 57.21, 76.63] without and with other hyperintense structures, respectively). The use of our computational model enables comparing segmentation methods and identifying their advantages and disadvantages, thereby providing means for testing and optimising pipelines for ongoing and future studies.

Imaging of the pial arterial vasculature of the human brain in vivo using high-resolution 7T time-of-flight angiography

The pial arterial vasculature of the human brain is the only blood supply to the neocortex, but quantitative data on the morphology and topology of these mesoscopic arteries (diameter 50–300 µm) remains scarce. Because it is commonly assumed that blood flow velocities in these vessels are prohibitively slow, non-invasive time-of-flight magnetic resonance angiography (TOF-MRA)—which is well suited to high 3D imaging resolutions—has not been applied to imaging the pial arteries. Here, we provide a theoretical framework that outlines how TOF-MRA can visualize small pial arteries in vivo, by employing extremely small voxels at the size of individual vessels. We then provide evidence for this theory by imaging the pial arteries at 140 µm isotropic resolution using a 7 Tesla (T) magnetic resonance imaging (MRI) scanner and prospective motion correction, and show that pial arteries one voxel width in diameter can be detected. We conclude that imaging pial arteries is not limited by slow blood flow, but instead by achievable image resolution. This study represents the first targeted, comprehensive account of imaging pial arteries in vivo in the human brain. This ultra-high-resolution angiography will enable the characterization of pial vascular anatomy across the brain to investigate patterns of blood supply and relationships between vascular and functional architecture.

Simultaneous multifactor Bayesian analysis (SiMBA) of PET time activity curve data

Positron emission tomography (PET) is an in vivo imaging method essential for studying the neurochemical pathophysiology of psychiatric and neurological disease. However, its high cost and exposure of participants to radiation make it unfeasible to employ large sample sizes. The major shortcoming of PET imaging is therefore its lack of power for studying clinically-relevant research questions. Here, we introduce a new method for performing PET quantification and analysis called SiMBA, which helps to alleviate these issues by improving the efficiency of PET analysis by exploiting similarities between both individuals and regions within individuals. In simulated [¹¹C]WAY100635 data, SiMBA greatly improves both statistical power and the consistency of effect size estimation without affecting the false positive rate. This approach makes use of hierarchical, multifactor, multivariate Bayesian modelling to effectively borrow strength across the whole dataset to improve stability and robustness to measurement error. In so doing, parameter identifiability and estimation are improved, without sacrificing model interpretability. This comes at the cost of increased computational overhead, however this is practically negligible relative to the time taken to collect PET data. This method has the potential to make it possible to test clinically-relevant hypotheses which could never be studied before given the practical constraints. Furthermore, because this method does not require any additional information over and above that required for traditional analysis, it makes it possible to re-examine data which has already previously been collected at great expense. In the absence of dramatic advancements in PET image data quality, radiotracer development, or data sharing, PET imaging has been fundamentally limited in the scope of research hypotheses which could be studied. This method, especially combined with the recent steps taken by the PET imaging community to embrace data sharing, will make it possible to greatly improve the research possibilities and clinical relevance of PET neuroimaging.

Medullary vein architecture modulates the white matter BOLD cerebrovascular reactivity signal response to CO2: Observations from high-resolution T2* weighted imaging at 7T

Brain stress testing using blood oxygenation level-dependent (BOLD) MRI to evaluate changes in cerebrovascular reactivity (CVR) is of growing interest for evaluating white matter integrity. However, even under healthy conditions, the white matter BOLD-CVR response differs notably from that observed in the gray matter. In addition to actual arterial vascular control, the venous draining topology may influence the WM-CVR response leading to signal delays and dispersions. These types of alterations in hemodynamic parameters are sometimes linked with pathology, but may also arise from differences in normal venous architecture. In this work, high-resolution T2*weighted anatomical images combined with BOLD imaging during a hypercapnic breathing protocol were acquired using a 7 tesla MRI system. Hemodynamic parameters including base CVR, hemodynamic lag, lag-corrected CVR, response onset and signal dispersion, and finally ΔCVR (corrected CVR minus base CVR) were calculated in 8 subjects. Parameter maps were spatially normalized and correlated against an MNI-registered white matter medullary vein atlas. Moderate correlations (Pearson's rho) were observed between medullary vessel frequency (MVF) and ΔCVR (0.52; 0.58 for total WM), MVF and hemodynamic lag (0.42; 0.54 for total WM), MVF and signal dispersion (0.44; 0.53 for total WM), and finally MVF and signal onset (0.43; 0.52 for total WM). Results indicate that, when assessed in the context of the WM venous architecture, changes in the response shape may only be partially reflective of the actual vascular reactivity response occurring further upstream by control vessels. This finding may have implications when attributing diseases mechanisms and/or progression to presumed impaired WM BOLD-CVR.

Abnormalities in normal-appearing white matter from which multiple sclerosis lesions arise

Normal-appearing white matter is far from normal in multiple sclerosis; little is known about the precise pathology or spatial pattern of this alteration and its relation to subsequent lesion formation. This study was undertaken to evaluate normal-appearing white matter abnormalities in brain areas where multiple sclerosis lesions subsequently form, and to investigate the spatial distribution of normal-appearing white matter abnormalities in persons with multiple sclerosis. Brain MRIs of pre-lesion normal-appearing white matter were analysed in participants with new T2 lesions, pooled from three clinical trials: SYNERGY (NCT01864148; n = 85 with relapsing multiple sclerosis) was the test data set; ASCEND (NCT01416181; n = 154 with secondary progressive multiple sclerosis) and ADVANCE (NCT00906399; n = 261 with relapsing-remitting multiple sclerosis) were used as validation data sets. Focal normal-appearing white matter tissue state was analysed prior to lesion formation in areas where new T2 lesions later formed (pre-lesion normal-appearing white matter) using normalized magnetization transfer ratio and T2-weighted (nT2) intensities, and compared with overall normal-appearing white matter and spatially matched contralateral normal-appearing white matter. Each outcome was analysed using linear mixed-effects models. Follow-up time (as a categorical variable), patient-level characteristics (including treatment group) and other baseline variables were treated as fixed effects. In SYNERGY, nT2 intensity was significantly higher, and normalized magnetization transfer ratio was lower in pre-lesion normal-appearing white matter versus overall and contralateral normal-appearing white matter at all time points up to 24 weeks before new T2 lesion onset. In ASCEND and ADVANCE (for which normalized magnetization transfer ratio was not available), nT2 intensity in pre-lesion normal-appearing white matter was significantly higher compared to both overall and contralateral normal-appearing white matter at all pre-lesion time points extending up to 2 years prior to lesion formation. In all trials, nT2 intensity in the contralateral normal-appearing white matter was also significantly higher at all pre-lesion time points compared to overall normal-appearing white matter. Brain atlases of normal-appearing white matter abnormalities were generated using measures of voxel-wise differences in normalized magnetization transfer ratio of normal-appearing white matter in persons with multiple sclerosis compared to scanner-matched healthy controls. We observed that overall spatial distribution of normal-appearing white matter abnormalities in persons with multiple sclerosis largely recapitulated the anatomical distribution of probabilities of T2 hyperintense lesions. Overall, these findings suggest that intrinsic spatial properties and/or longstanding precursory abnormalities of normal-appearing white matter tissue may contribute to the risk of autoimmune acute demyelination in multiple sclerosis.

Physiological and motion signatures in static and time-varying functional connectivity and their subject identifiability

Human brain connectivity yields significant potential as a noninvasive biomarker. Several studies have used fMRI-based connectivity fingerprinting to characterize individual patterns of brain activity. However, it is not clear whether these patterns mainly reflect neural activity or the effect of physiological and motion processes. To answer this question, we capitalize on a large data sample from the Human Connectome Project and rigorously investigate the contribution of the aforementioned processes on functional connectivity (FC) and time-varying FC, as well as their contribution to subject identifiability. We find that head motion, as well as heart rate and breathing fluctuations, induce artifactual connectivity within distinct resting-state networks and that they correlate with recurrent patterns in time-varying FC. Even though the spatiotemporal signatures of these processes yield above-chance levels in subject identifiability, removing their effects at the preprocessing stage improves identifiability, suggesting a neural component underpinning the inter-individual differences in connectivity.

Investigating the effect of flow compensation and quantitative susceptibility mapping method on the accuracy of venous susceptibility measurement

Quantitative susceptibility mapping (QSM) is a promising non-invasive method for obtaining information relating to oxygen metabolism. However, the optimal acquisition sequence and QSM reconstruction method for reliable venous susceptibility measurements are unknown. Full flow compensation is generally recommended to correct for the influence of venous blood flow, although the effect of flow compensation on the accuracy of venous susceptibility values has not been systematically evaluated. In this study, we investigated the effect of different acquisition sequences, including different flow compensation schemes, and different QSM reconstruction methods on venous susceptibilities. Ten healthy subjects were scanned with five or six distinct QSM sequence designs using monopolar readout gradients and different flow compensation schemes. All data sets were processed using six different QSM pipelines and venous blood susceptibility was evaluated in whole-brain segmentations of the venous vasculature and single veins. The quality of vein segmentations and the accuracy of venous susceptibility values were analyzed and compared between all combinations of sequences and reconstruction methods. The influence of the QSM reconstruction method on average venous susceptibility values was found to be 2.7-11.6 times greater than the influence of the acquisition sequence, including flow compensation. The majority of the investigated QSM reconstruction methods tended to underestimate venous susceptibility values in the vein segmentations that were obtained. In summary, we found that multi-echo gradient-echo acquisition sequences without full flow compensation yielded venous susceptibility values comparable to sequences with full flow compensation. However, the QSM reconstruction method had a great influence on susceptibility values and thus needs to be selected carefully for accurate venous QSM.

Evolution of Human Brain Atlases in Terms of Content, Applications, Functionality, and Availability

Human brain atlases have been evolving tremendously, propelled recently by brain big projects, and driven by sophisticated imaging techniques, advanced brain mapping methods, vast data, analytical strategies, and powerful computing. We overview here this evolution in four categories: content, applications, functionality, and availability, in contrast to other works limited mostly to content. Four atlas generations are distinguished: early cortical maps, print stereotactic atlases, early digital atlases, and advanced brain atlas platforms, and 5 avenues in electronic atlases spanning the last two generations. Content-wise, new electronic atlases are categorized into eight groups considering their scope, parcellation, modality, plurality, scale, ethnicity, abnormality, and a mixture of them. Atlas content developments in these groups are heading in 23 various directions. Application-wise, we overview atlases in neuroeducation, research, and clinics, including stereotactic and functional neurosurgery, neuroradiology, neurology, and stroke. Functionality-wise, tools and functionalities are addressed for atlas creation, navigation, individualization, enabling operations, and application-specific. Availability is discussed in media and platforms, ranging from mobile solutions to leading-edge supercomputers, with three accessibility levels. The major application-wise shift has been from research to clinical practice, particularly in stereotactic and functional neurosurgery, although clinical applications are still lagging behind the atlas content progress. Atlas functionality also has been relatively neglected until recently, as the management of brain data explosion requires powerful tools. We suggest that the future human brain atlas-related research and development activities shall be founded on and benefit from a standard framework containing the core virtual brain model cum the brain atlas platform general architecture.

Automatic collateral circulation scoring in ischemic stroke using 4D CT angiography with low-rank and sparse matrix decomposition

PURPOSE

Sufficient collateral blood supply is crucial for favorable outcomes with endovascular treatment. The current practice of collateral scoring relies on visual inspection and thus can suffer from inter and intra-rater inconsistency. We present a robust and automatic method to score cerebral collateral blood supply to aid ischemic stroke treatment decision making. The developed method is based on 4D dynamic CT angiography (CTA) and the ASPECTS scoring protocol.

METHODS

The proposed method, ACCESS (Automatic Collateral Circulation Evaluation in iSchemic Stroke), estimates a target patient's unfilled cerebrovasculature in contrast-enhanced CTA using the lack of contrast agent due to clotting. To do so, the fast robust matrix completion algorithm with in-face extended Frank-Wolfe optimization is applied on a cohort of healthy subjects and a target patient, to model the patient's unfilled vessels and the estimated full vasculature as sparse and low-rank components, respectively. The collateral score is computed as the ratio of the unfilled vessels to the full vasculature, mimicking existing clinical protocols.

RESULTS

ACCESS was tested with 46 stroke patients and obtained an overall accuracy of 84.78%. The optimal threshold selection was evaluated using a receiver operating characteristics curve with the leave-one-out approach, and a mean area under the curve of 85.39% was obtained.

CONCLUSION

ACCESS automates collateral scoring to mitigate the shortcomings of the standard clinical practice. It is a robust approach, which resembles how radiologists score clinical scans, and can be used to help radiologists in clinical decisions of stroke treatment.

Denoising High-Field Multi-Dimensional MRI With Local Complex PCA

Modern high field and ultra high field magnetic resonance imaging (MRI) experiments routinely collect multi-dimensional data with high spatial resolution, whether multi-parametric structural, diffusion or functional MRI. While diffusion and functional imaging have benefited from recent advances in multi-dimensional signal analysis and denoising, structural MRI has remained untouched. In this work, we propose a denoising technique for multi-parametric quantitative MRI, combining a highly popular denoising method from diffusion imaging, over-complete local PCA, with a reconstruction of the complex-valued MR signal in order to define stable estimates of the noise in the decomposition. With this approach, we show signal to noise ratio (SNR) improvements in high resolution MRI without compromising the spatial accuracy or generating spurious perceptual boundaries.

A critical assessment of data quality and venous effects in sub-millimeter fMRI

Advances in hardware, pulse sequences, and reconstruction techniques have made it possible to perform functional magnetic resonance imaging (fMRI) at sub-millimeter resolution while maintaining high spatial coverage and acceptable signal-to-noise ratio. Here, we examine whether sub-millimeter fMRI can be used as a routine method for obtaining accurate measurements of fine-scale local neural activity. We conducted fMRI in human visual cortex during a simple event-related visual experiment (7 T, gradient-echo EPI, 0.8-mm isotropic voxels, 2.2-s sampling rate, 84 slices), and developed analysis and visualization tools to assess the quality of the data. Our results fall along three lines of inquiry. First, we find that the acquired fMRI images, combined with appropriate surface-based processing, provide reliable and accurate measurements of fine-scale blood oxygenation level dependent (BOLD) activity patterns. Second, we show that the highly folded structure of cortex causes substantial biases on spatial resolution and data visualization. Third, we examine the well-recognized issue of venous contributions to fMRI signals. In a systematic assessment of large sections of cortex measured at a fine scale, we show that time-averaged T2*-weighted EPI intensity is a simple, robust marker of venous effects. These venous effects are unevenly distributed across cortex, are more pronounced in gyri and outer cortical depths, and are, to a certain degree, in consistent locations across subjects relative to cortical folding. Furthermore, we show that these venous effects are strongly correlated with BOLD responses evoked by the experiment. We conclude that sub-millimeter fMRI can provide robust information about fine-scale BOLD activity patterns, but special care must be exercised in visualizing and interpreting these patterns, especially with regards to the confounding influence of the brain's vasculature. To help translate these methodological findings to neuroscience research, we provide practical suggestions for both high-resolution and standard-resolution fMRI studies.