Quantitative functional MRI: concepts, issues and future challenges

A systematic review of preclinical and clinical transcranial ultrasound neuromodulation and opportunities for functional connectomics

BACKGROUND

Low-intensity transcranial ultrasound has surged forward as a non-invasive and disruptive tool for neuromodulation with applications in basic neuroscience research and the treatment of neurological and psychiatric conditions.

OBJECTIVE

To provide a comprehensive overview and update of preclinical and clinical transcranial low intensity ultrasound for neuromodulation and emphasize the emerging role of functional brain mapping to guide, better understand, and predict responses.

METHODS

A systematic review was conducted by searching the Web of Science and Scopus databases for studies on transcranial ultrasound neuromodulation, both in humans and animals.

RESULTS

187 relevant studies were identified and reviewed, including 116 preclinical and 71 clinical reports with subjects belonging to diverse cohorts. Milestones of ultrasound neuromodulation are described within an overview of the broader landscape. General neural readouts and outcome measures are discussed, potential confounds are noted, and the emerging use of functional magnetic resonance imaging is highlighted.

CONCLUSION

Ultrasound neuromodulation has emerged as a powerful tool to study and treat a range of conditions and its combination with various neural readouts has significantly advanced this platform. In particular, the use of functional magnetic resonance imaging has yielded exciting inferences into ultrasound neuromodulation and has the potential to advance our understanding of brain function, neuromodulatory mechanisms, and ultimately clinical outcomes. It is anticipated that these preclinical and clinical trials are the first of many; that transcranial low intensity focused ultrasound, particularly in combination with functional magnetic resonance imaging, has the potential to enhance treatment for a spectrum of neurological conditions.

Mitigation of Fetal Irradiation Injury from Mid-Gestation Total Body Radiation with Mitochondrial-Targeted GS-Nitroxide JP4-039

Victims of a radiation terrorist event will include pregnant women and unborn fetuses. Mitochondrial dysfunction and oxidative stress are key pathogenic factors of fetal irradiation injury. The goal of this preclinical study is to investigate the efficacy of mitigating fetal irradiation injury by maternal administration of the mitochondrial-targeted gramicidin S (GS)- nitroxide radiation mitigator, JP4-039. Pregnant female C57BL/6NTac mice received 3 Gy total body ionizing irradiation (TBI) at mid-gestation embryonic day 13.5 (E13.5). Using novel time- and-motion-resolved 4D in utero magnetic resonance imaging (4D-uMRI), we found TBI caused extensive injury to the fetal brain that included cerebral hemorrhage, loss of cerebral tissue, and hydrocephalus with excessive accumulation of cerebrospinal fluid (CSF). Histopathology of the fetal mouse brain showed broken cerebral vessels and elevated apoptosis. Further use of novel 4D Oxy-wavelet MRI capable of probing in vivo mitochondrial function in intact brain revealed significant reduction of mitochondrial function in the fetal brain after 3Gy TBI. This was validated by ex vivo Oroboros mitochondrial respirometry. Maternal administration JP4-039 one day after TBI (E14.5), which can pass through the placental barrier, significantly reduced fetal brain radiation injury and improved fetal brain mitochondrial respiration. This also preserved cerebral brain tissue integrity and reduced cerebral hemorrhage and cell death. As JP4-039 administration did not change litter sizes or fetus viability, together these findings indicate JP4-039 can be deployed as a safe and effective mitigator of fetal radiation injury from mid-gestational in utero ionizing radiation exposure.

One Sentence Summary

Mitochondrial-targeted gramicidin S (GS)-nitroxide JP4-039 is safe and effective radiation mitigator for mid-gestational fetal irradiation injury.

Specificity of Quantitative Functional Brain Mapping with Arterial Spin-Labeling for Preoperative Assessment

BACKGROUND AND PURPOSE

Arterial spin-labeling is a noninvasive MR imaging technique allowing direct and quantitative measurement of brain perfusion. Arterial spin-labeling is well-established in clinics for investigating the overall cerebral perfusion, but it is still occasionally employed during tasks. The typical contrast for functional MR imaging is blood oxygen level-dependent (BOLD) imaging, whose specificity could be biased in neurologic patients due to altered neurovascular coupling. This work aimed to validate the use of functional ASL as a noninvasive tool for presurgical functional brain mapping. This is achieved by comparing the spatial accuracy of functional ASL with transcranial magnetic stimulation as the criterion standard.

MATERIALS AND METHODS

Twenty-eight healthy participants executed a motor task and received a somatosensory stimulation, while BOLD imaging and arterial spin-labeling were acquired simultaneously. Transcranial magnetic stimulation was subsequently used to define hand somatotopy.

RESULTS

Functional ASL was found more adjacent to transcranial magnetic stimulation than BOLD imaging, with a significant shift along the inferior-to-superior direction. With respect to BOLD imaging, functional ASL was localized significantly more laterally, anteriorly, and inferiorly during motor tasks and pneumatic stimulation.

CONCLUSIONS

Our results confirm the specificity of functional ASL in targeting the regional neuronal excitability. Functional ASL could be considered as a valid supplementary technique to BOLD imaging for presurgical mapping when spatial accuracy is crucial for delineating eloquent cortex.

Imagining the future of optical microscopy: everything, everywhere, all at once

The optical microscope has revolutionized biology since at least the 17th Century. Since then, it has progressed from a largely observational tool to a powerful bioanalytical platform. However, realizing its full potential to study live specimens is hindered by a daunting array of technical challenges. Here, we delve into the current state of live imaging to explore the barriers that must be overcome and the possibilities that lie ahead. We venture to envision a future where we can visualize and study everything, everywhere, all at once - from the intricate inner workings of a single cell to the dynamic interplay across entire organisms, and a world where scientists could access the necessary microscopy technologies anywhere.

Correspondence between BOLD fMRI task response and cerebrovascular reactivity across the cerebral cortex

BOLD sensitivity to baseline perfusion and blood volume is a well-acknowledged fMRI confound. Vascular correction techniques based on cerebrovascular reactivity (CVR) might reduce variance due to baseline cerebral blood volume, however this is predicated on an invariant linear relationship between CVR and BOLD signal magnitude. Cognitive paradigms have relatively low signal, high variance and involve spatially heterogenous cortical regions; it is therefore unclear whether the BOLD response magnitude to complex paradigms can be predicted by CVR. The feasibility of predicting BOLD signal magnitude from CVR was explored in the present work across two experiments using different CVR approaches. The first utilized a large database containing breath-hold BOLD responses and 3 different cognitive tasks. The second experiment, in an independent sample, calculated CVR using the delivery of a fixed concentration of carbon dioxide and a different cognitive task. An atlas-based regression approach was implemented for both experiments to evaluate the shared variance between task-invoked BOLD responses and CVR across the cerebral cortex. Both experiments found significant relationships between CVR and task-based BOLD magnitude, with activation in the right cuneus (R ² = 0.64) and paracentral gyrus (R ² = 0.71), and the left pars opercularis (R ² = 0.67), superior frontal gyrus (R ² = 0.62) and inferior parietal cortex (R ² = 0.63) strongly predicted by CVR. The parietal regions bilaterally were highly consistent, with linear regressions significant in these regions for all four tasks. Group analyses showed that CVR correction increased BOLD sensitivity. Overall, this work suggests that BOLD signal response magnitudes to cognitive tasks are predicted by CVR across different regions of the cerebral cortex, providing support for the use of correction based on baseline vascular physiology.

No difference in cerebral perfusion between the wild-type and the 5XFAD mouse model of Alzheimer's disease

Neuroimaging with [2,2-dimethyl-3-[(2R,3E)-3-oxidoiminobutan-2-yl]azanidylpropyl]-[(2R,3E)-3-hydroxyiminobutan-2-yl]azanide;oxo(⁹⁹Tc)technetium-99(3+) ([^99mTc]HMPAO) single photon emission computed tomography (SPECT) is used in Alzheimer's disease (AD) to evaluate regional cerebral blood flow (rCBF). Hypoperfusion in select temporoparietal regions has been observed in human AD. However, it is unknown whether AD hypoperfusion signatures are also present in the 5XFAD mouse model. The current study was undertaken to compare baseline brain perfusion between 5XFAD and wild-type (WT) mice using [99mTc]HMPAO SPECT and determine whether hypoperfusion is recapitulated in 5XFAD mice. 5XFAD and WT mice underwent a 45 min SPECT scan, 20 min after [^99mTc]HMPAO administration. Whole brain and regional standardized uptake values (SUV) and regional relative standardized uptake values (SUVR) with whole brain reference were compared between groups. Brain perfusion was similar between WT and 5XFAD brains. Whole brain [^99mTc]HMPAO retention revealed no significant difference in SUV (5XFAD, 0.372 ± 0.762; WT, 0.640 ± 0.955; p = 0.536). Similarly, regional analysis revealed no significant differences in [^99mTc]HMPAO metrics between groups (SUV: 0.357 ≤ p ≤ 0.640; SUVR: 0.595 ≤ p ≤ 0.936). These results suggest apparent discrepancies in rCBF between human AD and the 5XFAD model. Establishing baseline perfusion patterns in 5XFAD mice is essential to inform pre-clinical diagnostic and therapeutic drug discovery programs.

Application of DTI and fMRI in moyamoya disease

Moyamoya disease (MMD) is a chronic and progressive cerebrovascular stenosis or occlusive disease that occurs near Willis blood vessels. Diffusion tensor imaging (DTI) and functional magnetic resonance imaging (fMRI) are used to detect the microstructure of white matter and the function of gray matter, respectively. The damage of these structures will lead to the change of cognitive level in patients with moyamoya disease. In this paper, the principles of DTI and fMRI, their applications and challenges in moyamoya disease are reviewed.

Mapping oxidative metabolism in the human brain with calibrated fMRI in health and disease

Conventional functional MRI (fMRI) with blood-oxygenation level dependent (BOLD) contrast is an important tool for mapping human brain activity non-invasively. Recent interest in quantitative fMRI has renewed the importance of oxidative neuroenergetics as reflected by cerebral metabolic rate of oxygen consumption (CMR_O2) to support brain function. Dynamic CMR_O2 mapping by calibrated fMRI require multi-modal measurements of BOLD signal along with cerebral blood flow (CBF) and/or volume (CBV). In human subjects this "calibration" is typically performed using a gas mixture containing small amounts of carbon dioxide and/or oxygen-enriched medical air, which are thought to produce changes in CBF (and CBV) and BOLD signal with minimal or no CMR_O2 changes. However non-human studies have demonstrated that the "calibration" can also be achieved without gases, revealing good agreement between CMR_O2 changes and underlying neuronal activity (e.g., multi-unit activity and local field potential). Given the simpler set-up of gas-free calibrated fMRI, there is evidence of recent clinical applications for this less intrusive direction. This up-to-date review emphasizes technological advances for such translational gas-free calibrated fMRI experiments, also covering historical progression of the calibrated fMRI field that is impacting neurological and neurodegenerative investigations of the human brain.

Reliability of task-specific neuronal activation assessed with functional PET, ASL and BOLD imaging

Mapping the neuronal response during cognitive processing is of crucial importance to gain new insights into human brain function. BOLD imaging and ASL are established MRI methods in this endeavor. Recently, the novel approach of functional PET (fPET) was introduced, enabling absolute quantification of glucose metabolism at rest and during task execution in a single measurement. Here, we report test-retest reliability of fPET in direct comparison to BOLD imaging and ASL. Twenty healthy subjects underwent two PET/MRI measurements, providing estimates of glucose metabolism, cerebral blood flow (CBF) and blood oxygenation. A cognitive task was employed with different levels of difficulty requiring visual-motor coordination. Task-specific neuronal activation was robustly detected with all three imaging approaches. The highest reliability was obtained for glucose metabolism at rest. Although this dropped during task performance it was still comparable to that of CBF. In contrast, BOLD imaging yielded high performance only for qualitative spatial overlap of task effects but not for quantitative comparison. Hence, the combined assessment of fPET and ASL offers reliable and simultaneous absolute quantification of glucose metabolism and CBF at rest and task.

Separating vascular and neuronal effects of age on fMRI BOLD signals

Accurate identification of brain function is necessary to understand the neurobiology of cognitive ageing, and thereby promote well-being across the lifespan. A common tool used to investigate neurocognitive ageing is functional magnetic resonance imaging (fMRI). However, although fMRI data are often interpreted in terms of neuronal activity, the blood oxygenation level-dependent (BOLD) signal measured by fMRI includes contributions of both vascular and neuronal factors, which change differentially with age. While some studies investigate vascular ageing factors, the results of these studies are not well known within the field of neurocognitive ageing and therefore vascular confounds in neurocognitive fMRI studies are common. Despite over 10 000 BOLD-fMRI papers on ageing, fewer than 20 have applied techniques to correct for vascular effects. However, neurovascular ageing is not only a confound in fMRI, but an important feature in its own right, to be assessed alongside measures of neuronal ageing. We review current approaches to dissociate neuronal and vascular components of BOLD-fMRI of regional activity and functional connectivity. We highlight emerging evidence that vascular mechanisms in the brain do not simply control blood flow to support the metabolic needs of neurons, but form complex neurovascular interactions that influence neuronal function in health and disease. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

The thermodynamics of thinking: connections between neural activity, energy metabolism and blood flow

Several current functional neuroimaging methods are sensitive to cerebral metabolism and cerebral blood flow (CBF) rather than the underlying neural activity itself. Empirically, the connections between metabolism, flow and neural activity are complex and somewhat counterintuitive: CBF and glycolysis increase more than seems to be needed to provide oxygen and pyruvate for oxidative metabolism, and the oxygen extraction fraction is relatively low in the brain and decreases when oxygen metabolism increases. This work lays a foundation for the idea that this unexpected pattern of physiological changes is consistent with basic thermodynamic considerations related to metabolism. In the context of this thermodynamic framework, the apparent mismatches in metabolic rates and CBF are related to preserving the entropy change of oxidative metabolism, specifically the O₂/CO₂ ratio in the mitochondria. However, the mechanism supporting this CBF response is likely not owing to feedback from a hypothetical O₂ sensor in tissue, but rather is consistent with feed-forward control by signals from both excitatory and inhibitory neural activity. Quantitative predictions of the thermodynamic framework, based on models of O₂ and CO₂ transport and possible neural drivers of CBF control, are in good agreement with a wide range of experimental data, including responses to neural activation, hypercapnia, hypoxia and high-altitude acclimatization. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

Oxygen extraction fraction mapping with multi-parametric quantitative BOLD MRI: Reduced transverse relaxation bias using 3D-GraSE imaging

Magnetic resonance imaging (MRI)-based quantification of the blood-oxygenation-level-dependent (BOLD) effect allows oxygen extraction fraction (OEF) mapping. The multi-parametric quantitative BOLD (mq-BOLD) technique facilitates relative OEF (rOEF) measurements with whole brain coverage in clinically applicable scan times. Mq-BOLD requires three separate scans of cerebral blood volume and transverse relaxation rates measured by gradient-echo (1/T₂*) and spin-echo (1/T₂). Although the current method is of clinical merit in patients with stroke, glioma and internal carotid artery stenosis (ICAS), there are relaxation measurement artefacts that impede the sensitivity of mq-BOLD and artificially elevate reported rOEF values.

We posited that T₂-related biases caused by slice refocusing imperfections during rapid 2D-GraSE (Gradient and Spin Echo) imaging can be reduced by applying 3D-GraSE imaging sequences, because the latter requires no slice selective pulses. The removal of T₂-related biases would decrease overestimated rOEF values measured by mq-BOLD. We characterized effects of T₂-related bias in mq-BOLD by comparing the initially employed 2D-GraSE and two proposed 3D-GraSE sequences to multiple single spin-echo reference measurements, both in vitro and in vivo. A phantom and 25 participants, including young and elderly healthy controls as well as ICAS-patients, were scanned. We additionally proposed a procedure to reliably identify and exclude artefact affected voxels. In the phantom, 3D-GraSE derived T₂ values had 57% lower deviation from the reference. For in vivo scans, the formerly overestimated rOEF was reduced by −27% (p < 0.001). We obtained rOEF = 0.51, which is much closer to literature values from positron emission tomography (PET) measurements. Furthermore, increased sensitivity to a focal rOEF elevation in an ICAS-patient was demonstrated.

In summary, the application of 3D-GraSE improves the mq-BOLD-based rOEF quantification while maintaining clinically feasible scan times. Thus, mq-BOLD with non-slice selective T₂ imaging is highly promising to improve clinical diagnostics of cerebrovascular diseases such as ICAS.

Cortical laminar resting-state signal fluctuations scale with the hypercapnic blood oxygenation level-dependent response

Calibrated functional magnetic resonance imaging can remove unwanted sources of signal variability in the blood oxygenation level‐dependent (BOLD) response. This is achieved by scaling, using information from a perfusion‐sensitive scan during a purely vascular challenge, typically induced by a gas manipulation or a breath‐hold task. In this work, we seek for a validation of the use of the resting‐state fluctuation amplitude (RSFA) as a scaling factor to remove vascular contributions from the BOLD response. Given the peculiarity of depth‐dependent vascularization in gray matter, BOLD and vascular space occupancy (VASO) data were acquired at submillimeter resolution and averaged across cortical laminae. RSFA from the primary motor cortex was, thus, compared to the amplitude of hypercapnia‐induced signal changes (tSD_hc) and with the M factor of the Davis model on a laminar level. High linear correlations were observed for RSFA and tSD_hc (R² = 0.92 ± 0.06) and somewhat reduced for RSFA and M (R² = 0.62 ± 0.19). Laminar profiles of RSFA‐normalized BOLD signal changes yielded good agreement with corresponding VASO profiles. Overall, this suggests that RSFA contains strong vascular components and is also modulated by baseline quantities contained in the M factor. We conclude that RSFA may replace the scaling factor tSD_hc for normalizing the laminar BOLD response.

The Neural Correlates of Consciousness and Attention: Two Sister Processes of the Brain

During the last three decades our understanding of the brain processes underlying consciousness and attention has significantly improved, mainly because of the advances in functional neuroimaging techniques. Still, caution is needed for the correct interpretation of these empirical findings, as both research and theoretical proposals are hampered by a number of conceptual difficulties. We review some of the most significant theoretical issues concerning the concepts of consciousness and attention in the neuroscientific literature, and put forward the implications of these reflections for a coherent model of the neural correlates of these brain functions. Even though consciousness and attention have an overlapping pattern of neural activity, they should be considered as essentially separate brain processes. The contents of phenomenal consciousness are supposed to be associated with the activity of multiple synchronized networks in the temporo-parietal-occipital areas. Only subsequently, attention, supported by fronto-parietal networks, enters the process of consciousness to provide focal awareness of specific features of reality.

Test-retest reliability of cerebral blood flow for assessing brain function at rest and during a vigilance task

Arterial spin labeled (ASL) perfusion magnetic resonance imaging (MRI) is increasingly used to assess regional brain activity and cerebrovascular function in both healthy and clinical populations. ASL perfusion imaging provides a quantitative measure of regional brain activity by determining absolute cerebral blood flow (CBF) values at a resting state or during task performance. However, the comparative reliability of these ASL measures is not well characterized. It is also unclear whether the test-retest reliability of absolute CBF or task-induced CBF change measures would be comparable to the reliability of task performance. In this study, fifteen healthy participants were scanned three times in a strictly controlled in-laboratory study while at rest and during performing a simple and reliable psychomotor vigilance test (PVT). The reliability of absolute CBF and task-induced CBF changes was evaluated using the intraclass correlation coefficient (ICC) and compared to that of task performance. Absolute CBF showed excellent test-retest reliability across the three scans for both resting and PVT scans. The reliability of regional absolute CBF was comparable to that of behavioral measures of PVT performance, and was slightly higher during PVT scans as compared with resting scans. Task-induced regional CBF changes demonstrated only poor to moderate reliability across three scans. These findings suggest that absolute CBF measures are more reliable than task-induced CBF changes for characterizing regional brain function, especially for longitudinal and clinical studies.

Toolbox for In Vivo Imaging of Host-Parasite Interactions at Multiple Scales

Animal models have for long been pivotal for parasitology research. Over the last few years, techniques such as intravital, optoacoustic and magnetic resonance imaging, optical projection tomography, and selective plane illumination microscopy developed promising potential for gaining insights into host-pathogen interactions by allowing different visualization forms in vivo and ex vivo. Advances including increased resolution, penetration depth, and acquisition speed, together with more complex image analysis methods, facilitate tackling biological problems previously impossible to study and/or quantify. Here we discuss advances and challenges in the in vivo imaging toolbox, which hold promising potential for the field of parasitology.

Dependence of the MR signal on the magnetic susceptibility of blood studied with models based on real microvascular networks

PURPOSE

The primary goal of this study was to estimate the value of β , the exponent in the power law relating changes of the transverse relaxation rate and intra-extravascular local magnetic susceptibility differences as Δ R 2 ∗ ∝ ( Δ χ ) β . The secondary objective was to evaluate any differences that might exist in the value of β obtained using a deoxyhemoglobin-weighted Δ χ distribution versus a constant Δ χ distribution assumed in earlier computations. The third objective was to estimate the value of β that is relevant for methods based on susceptibility contrast agents with a concentration of Δ χ higher than that used for BOLD fMRI calculations.

METHODS

Our recently developed model of real microvascular anatomical networks is used to extend the original simplified Monte-Carlo simulations to compute β from the first principles.

RESULTS

Our results show that β = 1 for most BOLD fMRI measurements of real vascular networks, as opposed to earlier predictions of β = 1 .5 using uniform Δ χ distributions. For perfusion or fMRI methods based on contrast agents, which generate larger values for Δ χ , β = 1 for B 0 ≤ 9.4 T, whereas at 14 T β can drop below 1 and the variation across subjects is large, indicating that a lower concentration of contrast agent with a lower value of Δ χ is desired for experiments at high B₀ .

CONCLUSION

These results improve our understanding of the relationship between R₂ ^* and the underlying microvascular properties. The findings will help to infer the cerebral metabolic rate of oxygen and cerebral blood volume from BOLD and perfusion MRI, respectively.

Cerebral blood volume changes during the BOLD post-stimulus undershoot measured with a combined normoxia/hyperoxia method

Cerebral blood flow (CBF) and blood oxygenation level dependent (BOLD) signal measurements make it possible to estimate steady-state changes in the cerebral metabolic rate of oxygen (CMRO2) with a calibrated BOLD method. However, extending this approach to measure the dynamics of CMRO2 requires an additional assumption: that deoxygenated cerebral blood volume (CBVdHb) follows CBF in a predictable way. A test-case for this assumption is the BOLD post-stimulus undershoot, for which one proposed explanation is a strong uncoupling of flow and blood volume with an elevated level of CBVdHb during the post-stimulus period compared to baseline due to slow blood volume recovery (Balloon Model). A challenge in testing this model is that CBVdHb differs from total blood volume, which can be measured with other techniques. In this study, the basic hypothesis of elevated CBVdHb during the undershoot was tested, based on the idea that the BOLD signal change when a subject switches from breathing a normoxic gas to breathing a hyperoxic gas is proportional to the absolute CBVdHb. In 19 subjects (8F), dual-echo BOLD responses were measured in primary visual cortex during a flickering radial checkerboard stimulus in normoxia, and the identical experiment was repeated in hyperoxia (50% O2/balance N2). The BOLD signal differences between normoxia and hyperoxia for the pre-stimulus baseline, stimulus, and post-stimulus periods were compared using an equivalent BOLD signal calculated from measured R2* changes to eliminate signal drifts. Relative to the pre-stimulus baseline, the average BOLD signal change from normoxia to hyperoxia was negative during the undershoot period (p = 0.0251), consistent with a reduction of CBVdHb and contrary to the prediction of the Balloon Model. Based on these results, the BOLD post-stimulus undershoot does not represent a case of strong uncoupling of CBVdHb and CBF, supporting the extension of current calibrated BOLD methods to estimate the dynamics of CMRO2.

Current technical approaches to brain energy metabolism

Neuroscience is a technology-driven discipline and brain energy metabolism is no exception. Once satisfied with mapping metabolic pathways at organ level, we are now looking to learn what it is exactly that metabolic enzymes and transporters do and when, where do they reside, how are they regulated, and how do they relate to the specific functions of neurons, glial cells, and their subcellular domains and organelles, in different areas of the brain. Moreover, we aim to quantify the fluxes of metabolites within and between cells. Energy metabolism is not just a necessity for proper cell function and viability but plays specific roles in higher brain functions such as memory processing and behavior, whose mechanisms need to be understood at all hierarchical levels, from isolated proteins to whole subjects, in both health and disease. To this aim, the field takes advantage of diverse disciplines including anatomy, histology, physiology, biochemistry, bioenergetics, cellular biology, molecular biology, developmental biology, neurology, and mathematical modeling. This article presents a well-referenced synopsis of the technical side of brain energy metabolism research. Detail and jargon are avoided whenever possible and emphasis is given to comparative strengths, limitations, and weaknesses, information that is often not available in regular articles.

Brain-derived neurotrophic factor Val66Met genotype and ovarian steroids interactively modulate working memory-related hippocampal function in women: a multimodal neuroimaging study

Preclinical evidence suggests that the actions of ovarian steroid hormones and brain-derived neurotrophic factor (BDNF) are highly convergent on brain function. Studies in humanized mice document an interaction between estrus cycle-related changes in estradiol secretion and BDNF Val⁶⁶Met genotype on measures of hippocampal function and anxiety-like behavior. We believe our multimodal imaging data provide the first demonstration in women that the effects of the BDNF Val/Met polymorphism on hippocampal function are selectively modulated by estradiol. In a 6-month pharmacological hormone manipulation protocol, healthy, regularly menstruating, asymptomatic women completed positron emission tomography (PET) and functional magnetic resonance imaging (fMRI) scans while performing the n-back working memory task during three hormone conditions: ovarian suppression induced by the gonadotropin-releasing hormone agonist, leuprolide acetate; leuprolide plus estradiol; and leuprolide plus progesterone. For each of the three hormone conditions, a discovery data set was obtained with oxygen-15 water regional cerebral blood flow PET in 39 healthy women genotyped for BDNF Val⁶⁶Met, and a confirmatory data set was obtained with fMRI in 27 women. Our results, in close agreement across the two imaging platforms, demonstrate an ovarian hormone-by-BDNF interaction on working memory-related hippocampal function (PET: F_2,37=9.11, P=0.00026 uncorrected, P=0.05, familywise error corrected with small volume correction; fMRI: F_2,25=5.43, P=0.01, uncorrected) that reflects differential hippocampal recruitment in Met carriers but only in the presence of estradiol. These findings have clinical relevance for understanding the neurobiological basis of individual differences in the cognitive and behavioral effects of ovarian steroids in women, and may provide a neurogenetic framework for understanding neuropsychiatric disorders related to reproductive hormones as well as illnesses with sex differences in disease expression and course.

Advances in Imaging Brain Metabolism

Metabolism is central to neuroimaging because it can reveal pathways by which neuronal and glial cells use nutrients to fuel their growth and function. We focus on advanced magnetic resonance imaging (MRI) and magnetic resonance spectroscopy (MRS) methods used in brain metabolic studies. ¹⁷O-MRS and ³¹P-MRS, respectively, provide rates of oxygen use and ATP synthesis inside mitochondria, whereas ¹⁹F-MRS enables measurement of cytosolic glucose metabolism. Calibrated functional MRI (fMRI), an advanced form of fMRI that uses contrast generated by deoxyhemoglobin, provides maps of oxygen use that track neuronal firing across brain regions. ¹³C-MRS is the only noninvasive method of measuring both glutamatergic neurotransmission and cell-specific energetics with signaling and nonsignaling purposes. Novel MRI contrasts, arising from endogenous diamagnetic agents and exogenous paramagnetic agents, permit pH imaging of glioma. Overall, these magnetic resonance methods for imaging brain metabolism demonstrate translational potential to better understand brain disorders and guide diagnosis and treatment.

The roadmap for estimation of cell-type-specific neuronal activity from non-invasive measurements

The computational properties of the human brain arise from an intricate interplay between billions of neurons connected in complex networks. However, our ability to study these networks in healthy human brain is limited by the necessity to use non-invasive technologies. This is in contrast to animal models where a rich, detailed view of cellular-level brain function with cell-type-specific molecular identity has become available due to recent advances in microscopic optical imaging and genetics. Thus, a central challenge facing neuroscience today is leveraging these mechanistic insights from animal studies to accurately draw physiological inferences from non-invasive signals in humans. On the essential path towards this goal is the development of a detailed 'bottom-up' forward model bridging neuronal activity at the level of cell-type-specific populations to non-invasive imaging signals. The general idea is that specific neuronal cell types have identifiable signatures in the way they drive changes in cerebral blood flow, cerebral metabolic rate of O2 (measurable with quantitative functional Magnetic Resonance Imaging), and electrical currents/potentials (measurable with magneto/electroencephalography). This forward model would then provide the 'ground truth' for the development of new tools for tackling the inverse problem-estimation of neuronal activity from multimodal non-invasive imaging data.This article is part of the themed issue 'Interpreting BOLD: a dialogue between cognitive and cellular neuroscience'.

Validation and optimization of hypercapnic-calibrated fMRI from oxygen-sensitive two-photon microscopy

Hypercapnic-calibrated fMRI allows the estimation of the relative changes in the cerebral metabolic rate of oxygen (rCMRO2) from combined BOLD and arterial spin labelling measurements during a functional task, and promises to permit more quantitative analyses of brain activity patterns. The estimation relies on a macroscopic model of the BOLD effect that balances oxygen delivery and consumption to predict haemoglobin oxygenation and the BOLD signal. The accuracy of calibrated fMRI approaches has not been firmly established, which is limiting their broader adoption. We use our recently developed microscopic vascular anatomical network model in mice as a ground truth simulator to test the accuracy of macroscopic, lumped-parameter BOLD models. In particular, we investigate the original Davis model and a more recent heuristic simplification. We find that these macroscopic models are inaccurate using the originally defined parameters, but that the accuracy can be significantly improved by redefining the model parameters to take on new values. In particular, we find that the parameter α that relates cerebral blood-volume changes to cerebral blood-flow changes is significantly smaller than typically assumed and that the optimal value changes with magnetic field strength. The results are encouraging in that they support the use of simple BOLD models to quantify BOLD signals, but further work is needed to understand the physiological interpretation of the redefined model parameters.This article is part of the themed issue 'Interpreting BOLD: a dialogue between cognitive and cellular neuroscience'.

Identification of neurovascular changes associated with cerebral amyloid angiopathy from subject-specific hemodynamic response functions

Cerebral amyloid angiopathy (CAA) is a small-vessel disease preferentially affecting posterior brain regions. Recent evidence has demonstrated the efficacy of functional MRI in detecting CAA-related neurovascular injury, however, it is unknown whether such perturbations are associated with changes in the hemodynamic response function (HRF). Here we estimated HRFs from two different brain regions from block design activation data, in light of recent findings demonstrating how block designs can accurately reflect HRF parameter estimates while maximizing signal detection. Patients with a diagnosis of probable CAA and healthy controls performed motor and visual stimulation tasks. Time-to-peak (TTP), full-width at half-maximum (FWHM), and area under the curve (AUC) of the estimated HRFs were compared between groups and to MRI features associated with CAA including cerebral microbleed (CMB) count. Motor HRFs in CAA patients showed significantly wider FWHM ( P = 0.006) and delayed TTP ( P = 0.03) compared to controls. In the patient group, visual HRF FWHM was positively associated with CMB count ( P = 0.03). These findings indicate that hemodynamic abnormalities in patients with CAA may be reflected in HRFs estimated from block designs across different brain regions. Moreover, visual FWHM may be linked to structural MR indications associated with CAA.

Molecular imaging with engineered physiology

In vivo imaging techniques are powerful tools for evaluating biological systems. Relating image signals to precise molecular phenomena can be challenging, however, due to limitations of the existing optical, magnetic and radioactive imaging probe mechanisms. Here we demonstrate a concept for molecular imaging which bypasses the need for conventional imaging agents by perturbing the endogenous multimodal contrast provided by the vasculature. Variants of the calcitonin gene-related peptide artificially activate vasodilation pathways in rat brain and induce contrast changes that are readily measured by optical and magnetic resonance imaging. CGRP-based agents induce effects at nanomolar concentrations in deep tissue and can be engineered into switchable analyte-dependent forms and genetically encoded reporters suitable for molecular imaging or cell tracking. Such artificially engineered physiological changes, therefore, provide a highly versatile means for sensitive analysis of molecular events in living organisms.

Physiological basis of vascular autocalibration (VasA): Comparison to hypercapnia calibration methods

Purpose

The statistical power of functional MRI (fMRI) group studies is significantly hampered by high intersubject spatial and magnitude variance. We recently presented a vascular autocalibration method (VasA) to account for vascularization differences between subjects and hence improve the sensitivity in group studies. Here, we validate the novel calibration method by means of direct comparisons of VasA with more established measures of baseline venous blood volume (and indirectly vascular reactivity), the M‐value.

Methods

Seven healthy volunteers participated in two 7 T (T) fMRI experiments to compare M‐values with VasA estimates: (i) a hypercapnia experiment to estimate voxelwise M‐value maps, and (ii) an fMRI experiment using visual stimulation to estimate voxelwise VasA maps.

Results

We show that VasA and M‐value calibration maps show the same spatial profile, providing strong evidence that VasA is driven by local variations in vascular reactivity as reflected in the M‐value.

Conclusion

The agreement of vascular reactivity maps obtained with VasA when compared with M‐value maps confirms empirically the hypothesis that the VasA method is an adequate tool to account for variations in fMRI response amplitudes caused by vascular reactivity differences in healthy volunteers. VasA can therefore directly account for them and increase the statistical power of group studies. The VasA toolbox is available as a statistical parametric mapping (SPM) toolbox, facilitating its general application. Magn Reson Med 78:1168–1173, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine.

EEG-fMRI Bayesian framework for neural activity estimation: a simulation study

OBJECTIVE

Due to the complementary nature of electroencephalography (EEG) and functional magnetic resonance imaging (fMRI), and given the possibility of simultaneous acquisition, the joint data analysis can afford a better understanding of the underlying neural activity estimation. In this simulation study we want to show the benefit of the joint EEG-fMRI neural activity estimation in a Bayesian framework.

APPROACH

We built a dynamic Bayesian framework in order to perform joint EEG-fMRI neural activity time course estimation. The neural activity is originated by a given brain area and detected by means of both measurement techniques. We have chosen a resting state neural activity situation to address the worst case in terms of the signal-to-noise ratio. To infer information by EEG and fMRI concurrently we used a tool belonging to the sequential Monte Carlo (SMC) methods: the particle filter (PF).

MAIN RESULTS

First, despite a high computational cost, we showed the feasibility of such an approach. Second, we obtained an improvement in neural activity reconstruction when using both EEG and fMRI measurements.

SIGNIFICANCE

The proposed simulation shows the improvements in neural activity reconstruction with EEG-fMRI simultaneous data. The application of such an approach to real data allows a better comprehension of the neural dynamics.

Impact of abnormal cerebrovascular reactivity on BOLD fMRI: a preliminary investigation of moyamoya disease

Blood oxygen level‐dependent (BOLD) functional magnetic resonance imaging (fMRI) studies of patients with cerebrovascular disease have largely ignored the confounds associated with abnormal cerebral blood flow, vascular reactivity and neurovascular coupling. We studied BOLD fMRI activation and cerebrovascular reactivity in moyamoya disease. To characterize the impact of remote vascular demands on BOLD fMRI measurements, we varied the vascular territories engaged by manipulating the experimental task performed by the participants. Vascular territories affected by disease were identified using BOLD cerebrovascular reactivity. Preliminary evidence from two patients pre‐ and postrevascularization surgery and four controls indicates that neurovascular coupling in affected brain regions can be modulated by the task‐related vascular demands in unaffected regions. In one patient studied, we observed that brain regions with improved cerebrovascular reactivity after surgery demonstrated normalized neurovascular coupling, that is the degree to which neurovascular coupling was modulated by task‐related vascular demands was decreased. We propose that variations in task‐dependent neurovascular coupling in patients with moyamoya disease are likely related to vascular steal. While preliminary, our findings are a proof of concept of the limitations of BOLD fMRI in cerebrovascular disease and suggest a role for assessment of cerebrovascular reactivity to improve interpretation of task‐related BOLD fMRI activation.

Behavior, neuropsychology and fMRI

Cognitive neuroscientists in the late 20th century began the task of identifying the part(s) of the brain concerned with normal behavior as manifest in the psychological capacities as affective powers, reasoning, behaving purposively and the pursuit of goals, following introduction of the 'functional magnetic resonance imaging' (fMRI) method for identifying brain activity. For this research program to be successful two questions require satisfactory answers. First, as the fMRI method can currently only be used on stationary subjects, to what extent can neuropsychological tests applicable to such stationary subjects be correlated with normal behavior. Second, to what extent can correlations between the various neuropsychological tests on the one hand, and sites of brain activity determined with fMRI on the other, be regarded as established. The extent to which these questions have yet received satisfactory answers is reviewed, and suggestions made both for improving correlations of neuropsychological tests with behavior as well as with the results of fMRI-based observations.

Cell type specificity of neurovascular coupling in cerebral cortex

Identification of the cellular players and molecular messengers that communicate neuronal activity to the vasculature driving cerebral hemodynamics is important for (1) the basic understanding of cerebrovascular regulation and (2) interpretation of functional Magnetic Resonance Imaging (fMRI) signals. Using a combination of optogenetic stimulation and 2-photon imaging in mice, we demonstrate that selective activation of cortical excitation and inhibition elicits distinct vascular responses and identify the vasoconstrictive mechanism as Neuropeptide Y (NPY) acting on Y1 receptors. The latter implies that task-related negative Blood Oxygenation Level Dependent (BOLD) fMRI signals in the cerebral cortex under normal physiological conditions may be mainly driven by the NPY-positive inhibitory neurons. Further, the NPY-Y1 pathway may offer a potential therapeutic target in cerebrovascular disease.

Building a Science of Individual Differences from fMRI

To date, fMRI research has been concerned primarily with evincing generic principles of brain function through averaging data from multiple subjects. Given rapid developments in both hardware and analysis tools, the field is now poised to study fMRI-derived measures in individual subjects, and to relate these to psychological traits or genetic variations. We discuss issues of validity, reliability and statistical assessment that arise when the focus shifts to individual subjects and that are applicable also to other imaging modalities. We emphasize that individual assessment of neural function with fMRI presents specific challenges and necessitates careful consideration of anatomical and vascular between-subject variability as well as sources of within-subject variability.

New horizons in neurometabolic and neurovascular coupling from calibrated fMRI

Neurovascular coupling relates changes in neuronal activity to constriction/dilation of microvessels. However neurometabolic coupling, which is less well known, relates alterations in neuronal activity with metabolic demands. The link between the blood oxygenation level dependent (BOLD) signal and neural activity opened doors for functional MRI (fMRI) to be a powerful neuroimaging tool in the neurosciences. But due to the complex makeup of BOLD contrast, researchers began to investigate the relationship between BOLD signal and blood flow and/or volume changes during functional brain activation, which together provided the tools to measure oxygen consumption on the basis of the biophysical model of BOLD. This field is called calibrated fMRI, thereby allowed probing of both neurometabolic and neurovascular couplings for a variety of health conditions in animals and humans. Calibrated fMRI may provide brain disorder biomarkers that could be used for monitoring effective therapies.

Alcohol-induced impairment of inhibitory control is linked to attenuated brain responses in right fronto-temporal cortex

BACKGROUND

A self-enhancing loop between impaired inhibitory control under alcohol and alcohol consumption has been proposed as a possible mechanism underlying dysfunctional drinking in susceptible people. However, the neural underpinnings of alcohol-induced impairment of inhibitory control are widely unknown.

METHODS

We measured inhibitory control in 50 young adults with a stop-signal task during functional magnetic resonance imaging. In a single-blind placebo-controlled cross-over design, all participants performed the stop-signal task once under alcohol with a breath alcohol concentration of .6 g/kg and once under placebo. In addition, alcohol consumption was assessed with a free-access alcohol self-administration paradigm in the same participants.

RESULTS

Inhibitory control was robustly decreased under alcohol compared with placebo, indicated by longer stop-signal reaction times. On the neural level, impaired inhibitory control under alcohol was associated with attenuated brain responses in the right fronto-temporal portion of the inhibition network that supports the attentional capture of infrequent stop-signals and subsequent updating of action plans from response execution to inhibition. Furthermore, the extent of alcohol-induced impairment of inhibitory control predicted free-access alcohol consumption.

CONCLUSIONS

We suggest that during inhibitory control alcohol affects cognitive processes preceding actual motor inhibition. Under alcohol, decreased brain responses in right fronto-temporal areas might slow down the attentional capture of infrequent stop-signals and subsequent updating of action plans, which leads to impaired inhibitory control. In turn, pronounced alcohol-induced impairment of inhibitory control might enhance alcohol consumption in young adults, which might promote future alcohol problems.

Dynamic functional imaging of brain glucose utilization using fPET-FDG

Glucose is the principal source of energy for the brain and yet the dynamic response of glucose utilization to changes in brain activity is still not fully understood. Positron emission tomography (PET) allows quantitative measurement of glucose metabolism using 2-[(18)F]-fluorodeoxyglucose (FDG). However, FDG PET in its current form provides an integral (or average) of glucose consumption over tens of minutes and lacks the temporal information to capture physiological alterations associated with changes in brain activity induced by tasks or drug challenges. Traditionally, changes in glucose utilization are inferred by comparing two separate scans, which significantly limits the utility of the method. We report a novel method to track changes in FDG metabolism dynamically, with higher temporal resolution than exists to date and within a single session. Using a constant infusion of FDG, we demonstrate that our technique (termed fPET-FDG) can be used in an analysis pipeline similar to fMRI to define within-session differential metabolic responses. We use visual stimulation to demonstrate the feasibility of this method. This new method has a great potential to be used in research protocols and clinical settings since fPET-FDG imaging can be performed with most PET scanners and data acquisition and analysis are straightforward. fPET-FDG is a highly complementary technique to MRI and provides a rich new way to observe functional changes in brain metabolism.

Variability of the coupling of blood flow and oxygen metabolism responses in the brain: a problem for interpreting BOLD studies but potentially a new window on the underlying neural activity

Recent studies from our group and others using quantitative fMRI methods have found that variations of the coupling ratio of blood flow (CBF) and oxygen metabolism (CMRO2) responses to a stimulus have a strong effect on the BOLD response. Across a number of studies an empirical pattern is emerging in the way CBF and CMRO2 changes are coupled to neural activation: if the stimulus is modulated to create a stronger response (e.g., increasing stimulus contrast), CBF is modulated more than CMRO2; on the other hand, if the brain state is altered such that the response to the same stimulus is increased (e.g., modulating attention, adaptation, or excitability), CMRO2 is modulated more than CBF. Because CBF and CMRO2 changes conflict in producing BOLD signal changes, this finding has an important implication for conventional BOLD-fMRI studies: the BOLD response exaggerates the effects of stimulus variation but is only weakly sensitive to modulations of the brain state that alter the response to a standard stimulus. A speculative hypothesis is that variability of the coupling ratio of the CBF and CMRO2 responses reflects different proportions of inhibitory and excitatory evoked activity, potentially providing a new window on neural activity in the human brain.

Obese patients after gastric bypass surgery have lower brain-hedonic responses to food than after gastric banding

OBJECTIVES

Roux-en-Y gastric bypass (RYGB) has greater efficacy for weight loss in obese patients than gastric banding (BAND) surgery. We hypothesise that this may result from different effects on food hedonics via physiological changes secondary to distinct gut anatomy manipulations.

DESIGN

We used functional MRI, eating behaviour and hormonal phenotyping to compare body mass index (BMI)-matched unoperated controls and patients after RYGB and BAND surgery for obesity.

RESULTS

Obese patients after RYGB had lower brain-hedonic responses to food than patients after BAND surgery. RYGB patients had lower activation than BAND patients in brain reward systems, particularly to high-calorie foods, including the orbitofrontal cortex, amygdala, caudate nucleus, nucleus accumbens and hippocampus. This was associated with lower palatability and appeal of high-calorie foods and healthier eating behaviour, including less fat intake, in RYGB compared with BAND patients and/or BMI-matched unoperated controls. These differences were not explicable by differences in hunger or psychological traits between the surgical groups, but anorexigenic plasma gut hormones (GLP-1 and PYY), plasma bile acids and symptoms of dumping syndrome were increased in RYGB patients.

CONCLUSIONS

The identification of these differences in food hedonic responses as a result of altered gut anatomy/physiology provides a novel explanation for the more favourable long-term weight loss seen after RYGB than after BAND surgery, highlighting the importance of the gut-brain axis in the control of reward-based eating behaviour.

Measurement of OEF and absolute CMRO2: MRI-based methods using interleaved and combined hypercapnia and hyperoxia

Blood oxygenation level dependent (BOLD) functional magnetic resonance imaging (fMRI) is most commonly used in a semi-quantitative manner to infer changes in brain activity. Despite the basis of the image contrast lying in the cerebral venous blood oxygenation level, quantification of absolute cerebral metabolic rate of oxygen consumption (CMRO2) has only recently been demonstrated. Here we examine two approaches to the calibration of fMRI signal to measure absolute CMRO2 using hypercapnic and hyperoxic respiratory challenges. The first approach is to apply hypercapnia and hyperoxia separately but interleaved in time and the second is a combined approach in which we apply hyperoxic challenges simultaneously with different levels of hypercapnia. Eleven healthy volunteers were studied at 3T using a dual gradient-echo spiral readout pulsed arterial spin labelling (ASL) imaging sequence. Respiratory challenges were conducted using an automated system of dynamic end-tidal forcing. A generalised BOLD signal model was applied, within a Bayesian estimation framework, that aims to explain the effects of modulation of CBF and arterial oxygen content to estimate venous deoxyhaemoglobin concentration ([dHb]0). Using CBF measurements combined with the estimated oxygen extraction fraction (OEF), absolute CMRO2 was calculated. The interleaved approach to hypercapnia and hyperoxia, as well as yielding estimates of CMRO2 and OEF demonstrated a significant increase in regional CBF, venous oxygen saturation (SvO2) (a decrease in OEF) and absolute CMRO2 in visual cortex in response to a continuous (20 min) visual task, demonstrating the potential for the method in measuring long term changes in CMRO2. The combined approach to oxygen and carbon dioxide modulation, as well as taking less time to acquire data, yielded whole brain grey matter estimates of CMRO2 and OEF of 184±45 μmol/100 g/min and 0.42±0.12 respectively, along with additional estimates of the vascular parameters α=0.33±0.06, the exponent relating relative increases in CBF to CBV, and β=1.35±0.13, the exponent relating deoxyhaemoglobin concentration to the relaxation rate R2*. Maps of cerebrovascular and cerebral metabolic parameters were also calculated. We show that combined modulation of oxygen and carbon dioxide can offer an experimentally more efficient approach to estimating OEF and absolute CMRO2 along with the additional vascular parameters that form an important part of the commonly used calibrated fMRI signal model.

The physics of functional magnetic resonance imaging (fMRI)

Functional magnetic resonance imaging (fMRI) is a methodology for detecting dynamic patterns of activity in the working human brain. Although the initial discoveries that led to fMRI are only about 20 years old, this new field has revolutionized the study of brain function. The ability to detect changes in brain activity has a biophysical basis in the magnetic properties of deoxyhemoglobin, and a physiological basis in the way blood flow increases more than oxygen metabolism when local neural activity increases. These effects translate to a subtle increase in the local magnetic resonance signal, the blood oxygenation level dependent (BOLD) effect, when neural activity increases. With current techniques, this pattern of activation can be measured with resolution approaching 1 mm(3) spatially and 1 s temporally. This review focuses on the physical basis of the BOLD effect, the imaging methods used to measure it, the possible origins of the physiological effects that produce a mismatch of blood flow and oxygen metabolism during neural activation, and the mathematical models that have been developed to understand the measured signals. An overarching theme is the growing field of quantitative fMRI, in which other MRI methods are combined with BOLD methods and analyzed within a theoretical modeling framework to derive quantitative estimates of oxygen metabolism and other physiological variables. That goal is the current challenge for fMRI: to move fMRI from a mapping tool to a quantitative probe of brain physiology.

Quantitative basis for neuroimaging of cortical laminae with calibrated functional MRI

Layer-specific neurophysiologic, hemodynamic, and metabolic measurements are needed to interpret high-resolution functional magnetic resonance imaging (fMRI) data in the cerebral cortex. We examined how neurovascular and neurometabolic couplings vary vertically in the rat's somatosensory cortex. During sensory stimulation we measured dynamic layer-specific responses of local field potential (LFP) and multiunit activity (MUA) as well as blood oxygenation level-dependent (BOLD) signal and cerebral blood volume (CBV) and blood flow (CBF), which in turn were used to calculate changes in oxidative metabolism (CMR(O2)) with calibrated fMRI. Both BOLD signal and CBV decreased from superficial to deep laminae, but these responses were not well correlated with either layer-specific LFP or MUA. However, CBF changes were quite stable across laminae, similar to LFP. However, changes in CMR(O2) and MUA varied across cortex in a correlated manner and both were reduced in superficial lamina. These results lay the framework for quantitative neuroimaging across cortical laminae with calibrated fMRI methods.