Body size interacts with the structure of the central nervous system: A multi-center in vivo neuroimaging study

Clinical research emphasizes the implementation of rigorous and reproducible study designs that rely on between-group matching or controlling for sources of biological variation such as subject's sex and age. However, corrections for body size (i.e. height and weight) are mostly lacking in clinical neuroimaging designs. This study investigates the importance of body size parameters in their relationship with spinal cord (SC) and brain magnetic resonance imaging (MRI) metrics. Data were derived from a cosmopolitan population of 267 healthy human adults (age 30.1±6.6 years old, 125 females). We show that body height correlated strongly or moderately with brain gray matter (GM) volume, cortical GM volume, total cerebellar volume, brainstem volume, and cross-sectional area (CSA) of cervical SC white matter (CSA-WM; 0.44≤r≤0.62). In comparison, age correlated weakly with cortical GM volume, precentral GM volume, and cortical thickness (-0.21≥r≥-0.27). Body weight correlated weakly with magnetization transfer ratio in the SC WM, dorsal columns, and lateral corticospinal tracts (-0.20≥r≥-0.23). Body weight further correlated weakly with the mean diffusivity derived from diffusion tensor imaging (DTI) in SC WM (r=-0.20) and dorsal columns (-0.21), but only in males. CSA-WM correlated strongly or moderately with brain volumes (0.39≤r≤0.64), and weakly with precentral gyrus thickness and DTI-based fractional anisotropy in SC dorsal columns and SC lateral corticospinal tracts (-0.22≥r≥-0.25). Linear mixture of sex and age explained 26±10% of data variance in brain volumetry and SC CSA. The amount of explained variance increased at 33±11% when body height was added into the mixture model. Age itself explained only 2±2% of such variance. In conclusion, body size is a significant biological variable. Along with sex and age, body size should therefore be included as a mandatory variable in the design of clinical neuroimaging studies examining SC and brain structure.

Breath-hold BOLD fMRI without CO2 sampling enables estimation of venous cerebral blood volume: potential use in normalization of stimulus-evoked BOLD fMRI data

BOLD fMRI signal has been used in conjunction with vasodilatory stimulation as a marker of cerebrovascular reactivity (CVR): the relative change in cerebral blood flow (CBF) arising from a unit change in the vasodilatory stimulus. Using numerical simulations, we demonstrate that the variability in the relative BOLD signal change induced by vasodilation is strongly influenced by the variability in deoxyhemoglobin-containing cerebral blood volume (CBV), as this source of variability is likely to be more prominent than that of CVR. It may, therefore, be more appropriate to describe the relative BOLD signal change induced by an isometabolic vasodilation as a proxy of deoxygenated CBV (CBVdHb) rather than CVR. With this in mind, a new method was implemented to map a marker of CBVdHb, termed BOLD-CBV, based on the normalization of voxel-wise BOLD signal variation by an estimate of the intravascular venous BOLD signal from voxels filled with venous blood. The intravascular venous BOLD signal variation, recorded during repeated breath-holding, was extracted from the superior sagittal sinus in a cohort of 27 healthy volunteers and used as a regressor across the whole brain, yielding maps of BOLD-CBV. In the same cohort, we demonstrated the potential use of BOLD-CBV for the normalization of stimulus-evoked BOLD fMRI by comparing group-level BOLD fMRI responses to a visuomotor learning task with and without the inclusion of voxel-wise vascular covariates of BOLD-CBV and the BOLD signal change per mmHg variation in end-tidal carbon dioxide (BOLD-CVR). The empirical measure of BOLD-CBV accounted for more between-subject variability in the motor task-induced BOLD responses than BOLD-CVR estimated from end-tidal carbon dioxide recordings. The new method can potentially increase the power of group fMRI studies by including a measure of vascular characteristics and has the strong practical advantage of not requiring experimental measurement of end-tidal carbon dioxide, unlike traditional methods to estimate BOLD-CVR. It also more closely represents a specific physiological characteristic of brain vasculature than BOLD-CVR, namely blood volume.

Changes in brain perfusion with training-related visuomotor improvement in MS

Multiple sclerosis (MS) is a chronic inflammatory disease of the central nervous system. A better understanding of the mechanisms supporting brain plasticity in MS would help to develop targeted interventions to promote recovery. A total of 29 MS patients and 19 healthy volunteers underwent clinical assessment and multi-modal MRI acquisition [fMRI during serial reaction time task (SRT), DWI, T1w structural scans and ASL of resting perfusion] at baseline and after 4-weeks of SRT training. Reduction of functional hyperactivation was observed in MS patients following the training, shown by the stronger reduction of the BOLD response during task execution compared to healthy volunteers. The functional reorganization was accompanied by a positive correlation between improvements in task accuracy and the change in resting perfusion after 4 weeks' training in right angular and supramarginal gyri in MS patients. No longitudinal changes in WM and GM measures and no correlation between task performance improvements and brain structure were observed in MS patients. Our results highlight a potential role for CBF as an early marker of plasticity, in terms of functional (cortical reorganization) and behavioral (performance improvement) changes in MS patients that may help to guide future interventions that exploit preserved plasticity mechanisms.

Pathophysiology of multiple sclerosis damage and repair: Linking cerebral hypoperfusion to the development of irreversible tissue loss in multiple sclerosis using magnetic resonance imaging

BACKGROUND AND PURPOSE

Reduced cerebral perfusion has been observed in multiple sclerosis (MS) and may contribute to tissue loss both acutely and chronically. Here, we test the hypothesis that hypoperfusion occurs in MS and relates to the presence of irreversible tissue damage.

METHODS

In 91 patients with relapsing MS and 26 healthy controls (HC), gray matter (GM) cerebral blood flow (CBF) was assessed using pulsed arterial spin labeling. GM volume, T1 hypointense and T2 hyperintense lesion volumes (T1LV and T2LV, respectively), and the proportion of T2-hyperintense lesion volume that appears hypointense on T1-weighted magnetic resonance imaging (T1LV/T2LV) were quantified. GM CBF and GM volume were evaluated globally, as well as regionally, using an atlas-based approach.

RESULTS

Global GM CBF was lower in patients (56.9 ± 12.3 mL/100 g/min) than in HC (67.7 ± 10.0 mL/100 g/min; p < 0.001), a difference that was widespread across brain regions. Although total GM volume was comparable between groups, significant reductions were observed in a subset of subcortical structures. GM CBF negatively correlated with T1LV (r = -0.43, p = 0.0002) and T1LV/T2LV (r = -0.37, p = 0.0004), but not with T2LV.

CONCLUSIONS

GM hypoperfusion occurs in MS and is associated with irreversible white matter damage, thus suggesting that cerebral hypoperfusion may actively contribute and possibly precede neurodegeneration by hampering tissue repair abilities in MS.

Cerebral blood flow patterns in preterm and term neonates assessed with pseudo-continuous arterial spin labeling perfusion MRI

In preterm (PT) infants, regional cerebral blood flow (CBF) disturbances may predispose to abnormal brain maturation even without overt brain injury. Therefore, it would be informative to determine the spatial distribution of grey matter (GM) CBF in PT and full-term (FT) newborns at term-equivalent age (TEA) and to assess the relationship between the features of the CBF pattern and both prematurity and prematurity-related brain lesions. In this prospective study, we obtained measures of CBF in 66 PT (51 without and 15 with prematurity-related brain lesions) and 38 FT newborns through pseudo-continuous arterial spin labeling (pCASL) MRI acquired at TEA. The pattern of GM CBF was characterized by combining an atlas-based automated segmentation of structural MRI with spatial normalization and hierarchical clustering. The effects of gestational age (GA) at birth and brain injury on the CBF pattern were investigated. We identified 4 physiologically-derived clusters of brain regions that were labeled Fronto-Temporal, Parieto-Occipital, Insular-Deep GM (DGM) and Sensorimotor, from the least to the most perfused. We demonstrated that GM perfusion was associated with GA at birth in the Fronto-Temporal and Sensorimotor clusters, positively and negatively, respectively. Moreover, the presence of periventricular leukomalacia was associated with significantly increased Fronto-Temporal GM perfusion and decreased Insular-DGM perfusion, while the presence of germinal matrix hemorrhage appeared to mildly decrease the Insular-DGM perfusion. Prematurity and prematurity-related brain injury heterogeneously affect brain perfusion. ASL MRI may, therefore, have strong potential as a noninvasive tool for the accurate stratification of individuals at risk of domain-specific impairment.

Reduced brain oxygen metabolism in patients with multiple sclerosis: Evidence from dual-calibrated functional MRI

Cerebral energy deficiency is increasingly recognised as an important feature of multiple sclerosis (MS). Until now, we have lacked non-invasive imaging methods to quantify energy utilisation and mitochondrial function in the human brain. Here, we used novel dual-calibrated functional magnetic resonance imaging (dc-fMRI) to map grey-matter (GM) deoxy-haemoglobin sensitive cerebral blood volume (CBVdHb), cerebral blood flow (CBF), oxygen extraction fraction (OEF), and cerebral metabolic rate of oxygen consumption (CMRO2) in patients with MS (PwMS) and age/sex matched controls. By integrating a flow-diffusion model of oxygen transport, we evaluated the effective oxygen diffusivity of the capillary network (DC) and the partial pressure of oxygen at the mitochondria (PmO2). Significant between-group differences were observed as decreased CBF (p = 0.010), CMRO2 (p < 0.001) and DC (p = 0.002), and increased PmO2 (p = 0.043) in patients compared to controls. No significant differences were observed for CBVdHb (p = 0.389), OEF (p = 0.358), or GM volume (p = 0.302). Regional analysis showed widespread reductions in CMRO2 and DC for PwMS. Our findings may be indicative of reduced oxygen demand or utilisation in the MS brain and mitochondrial dysfunction. Our results suggest changes in brain physiology may precede MRI-detectable GM loss and may contribute to disease progression and neurodegeneration.

A flow-diffusion model of oxygen transport for quantitative mapping of cerebral metabolic rate of oxygen (CMRO2) with single gas calibrated fMRI

One promising approach for mapping CMRO₂ is dual-calibrated functional MRI (dc-fMRI). This method exploits the Fick Principle to combine estimates of CBF from ASL, and OEF derived from BOLD-ASL measurements during arterial O₂ and CO₂ modulations. Multiple gas modulations are required to decouple OEF and deoxyhemoglobin-sensitive blood volume. We propose an alternative single gas calibrated fMRI framework, integrating a model of oxygen transport, that links blood volume and CBF to OEF and creates a mapping between the maximum BOLD signal, CBF and OEF (and CMRO₂). Simulations demonstrated the method's viability within physiological ranges of mitochondrial oxygen pressure, P_mO₂, and mean capillary transit time. A dc-fMRI experiment, performed on 20 healthy subjects using O₂ and CO₂ challenges, was used to validate the approach. The validation conveyed expected estimates of model parameters (e.g., low P_mO₂), with spatially uniform OEF maps (grey matter, GM, OEF spatial standard deviation ≈ 0.13). GM OEF estimates obtained with hypercapnia calibrated fMRI correlated with dc-fMRI (r = 0.65, p = 2·10^-3). For 12 subjects, OEF measured with dc-fMRI and the single gas calibration method were correlated with whole-brain OEF derived from phase measures in the superior sagittal sinus (r = 0.58, p = 0.048; r = 0.64, p = 0.025 respectively). Simplified calibrated fMRI using hypercapnia holds promise for clinical application.

Cerebrovascular Variants and the Role of the Selfish Brain in Young-Onset Hypertension

BACKGROUND

Variants in the posterior anatomy of the cerebral circulation are associated with hypertension and lower cerebral blood flow in midlife (age ≈55 years); however, whether these variants are a result of aging or long-term exposure to high blood pressure is unclear. Additionally, the role these variants play in early onset of hypertension (<40 years) and poor cerebral perfusion in this population is unknown.

METHODS

We retrospectively examined whether specific cerebrovascular variants (vertebral artery hypoplasia and absent/hypoplastic posterior communicating arteries (an incomplete posterior circle of Willis) measured via magnetic resonance angiography) were associated with a diagnosis of hypertension in 220 young adults (<40 years; n=164 primary hypertensive [mean age±SD, 32±6 years] and n=56 [30±6 years] normotensive adults). Whether cerebrovascular variants were associated with lower cerebral blood flow (phase-contrast angiography) was measured in the hypertensive group only (n=146).

RESULTS

Binary logistic regression (adjusted for age, sex, and body mass index) showed that vertebral artery hypoplasia with an incomplete posterior circle of Willis was associated with hypertension diagnosis (P<0.001, odds ratio; 11.79 [95% CI, 3.34-41.58]). Vertebral artery hypoplasia plus an incomplete circle of Willis was associated with lower cerebral blood flow in young adults with hypertension (P=0.0172).

CONCLUSIONS

Vertebral artery hypoplasia plus an incomplete posterior circle of Willis independently predicts hypertension in young adults suggesting that this variant is not acquired with aging into midlife. Importantly this variant combination was associated with lower cerebral perfusion, which may have long-term consequences on cerebrovascular health in young adults with hypertension.

A comparison of GABA-ergic (propofol) and non-GABA-ergic (dexmedetomidine) sedation on visual and motor cortical oscillations, using magnetoencephalography

Studying changes in cortical oscillations can help elucidate the mechanistic link between receptor physiology and the clinical effects of anaesthetic drugs. Propofol, a GABA-ergic drug produces divergent effects on visual cortical activity: increasing induced gamma-band responses (GBR) while decreasing evoked responses. Dexmedetomidine, an α2- adrenergic agonist, differs from GABA-ergic sedatives both mechanistically and clinically as it allows easy arousability from deep sedation with less cognitive side-effects. Here we use magnetoencephalography (MEG) to characterize and compare the effects of GABA-ergic (propofol) and non-GABA-ergic (dexmedetomidine) sedation, on visual and motor cortical oscillations. Sixteen male participants received target-controlled infusions of propofol and dexmedetomidine, producing mild-sedation, in a placebo-controlled, cross-over study. MEG data was collected during a combined visuomotor task. The key findings were that propofol significantly enhanced visual stimulus induced GBR (44% increase in amplitude) while dexmedetomidine decreased it (40%). Propofol also decreased the amplitudes of the Mv100 (visual M100) (27%) and Mv150 (52%) visual evoked fields (VEF), whilst dexmedetomidine had no effect on these. During the motor task, neither drug had any significant effect on movement related gamma synchrony (MRGS), movement related beta de-synchronisation (MRBD) or Mm100 (movement-related M100) movement-related evoked fields (MEF), although dexmedetomidine slowed the Mm300. Dexmedetomidine increased (92%) post-movement beta synchronisation/rebound (PMBR) power while propofol reduced it (70%, statistically non- significant). Overall, dexmedetomidine and propofol, at equi-sedative doses, produce contrasting effects on visual induced GBR, VEF, PMBR and MEF. These findings provide a mechanistic link between the known receptor physiology of these sedative drugs with their known clinical effects and may be used to explore mechanisms of other anaesthetic drugs on human consciousness.

Contribution of Metabolomics to Multiple Sclerosis Diagnosis, Prognosis and Treatment

Metabolomics-based technologies map in vivo biochemical changes that may be used as early indicators of pathological abnormalities prior to the development of clinical symptoms in neurological conditions. Metabolomics may also reveal biochemical pathways implicated in tissue dysfunction and damage and thus assist in the development of novel targeted therapeutics for neuroinflammation and neurodegeneration. Metabolomics holds promise as a non-invasive, high-throughput and cost-effective tool for early diagnosis, follow-up and monitoring of treatment response in multiple sclerosis (MS), in combination with clinical and imaging measures. In this review, we offer evidence in support of the potential of metabolomics as a biomarker and drug discovery tool in MS. We also use pathway analysis of metabolites that are described as potential biomarkers in the literature of MS biofluids to identify the most promising molecules and upstream regulators, and show novel, still unexplored metabolic pathways, whose investigation may open novel avenues of research.

Open-access quantitative MRI data of the spinal cord and reproducibility across participants, sites and manufacturers

In a companion paper by Cohen-Adad et al. we introduce the spine generic quantitative MRI protocol that provides valuable metrics for assessing spinal cord macrostructural and microstructural integrity. This protocol was used to acquire a single subject dataset across 19 centers and a multi-subject dataset across 42 centers (for a total of 260 participants), spanning the three main MRI manufacturers: GE, Philips and Siemens. Both datasets are publicly available via git-annex. Data were analysed using the Spinal Cord Toolbox to produce normative values as well as inter/intra-site and inter/intra-manufacturer statistics. Reproducibility for the spine generic protocol was high across sites and manufacturers, with an average inter-site coefficient of variation of less than 5% for all the metrics. Full documentation and results can be found at https://spine-generic.rtfd.io/ . The datasets and analysis pipeline will help pave the way towards accessible and reproducible quantitative MRI in the spinal cord.

Hippocampal Blood Flow Is Increased After 20 min of Moderate-Intensity Exercise

Long-term exercise interventions have been shown to be a potent trigger for both neurogenesis and vascular plasticity. However, little is known about the underlying temporal dynamics and specifically when exercise-induced vascular adaptations first occur, which is vital for therapeutic applications. In this study, we investigated whether a single session of moderate-intensity exercise was sufficient to induce changes in the cerebral vasculature. We employed arterial spin labeling magnetic resonance imaging to measure global and regional cerebral blood flow (CBF) before and after 20 min of cycling. The blood vessels’ ability to dilate, measured by cerebrovascular reactivity (CVR) to CO₂ inhalation, was measured at baseline and 25-min postexercise. Our data showed that CBF was selectively increased by 10–12% in the hippocampus 15, 40, and 60 min after exercise cessation, whereas CVR to CO₂ was unchanged in all regions. The absence of a corresponding change in hippocampal CVR suggests that the immediate and transient hippocampal adaptations observed after exercise are not driven by a mechanical vascular change and more likely represents an adaptive metabolic change, providing a framework for exploring the therapeutic potential of exercise-induced plasticity (neural, vascular, or both) in clinical and aged populations.

Predictors of training-related improvement in visuomotor performance in patients with multiple sclerosis: A behavioural and MRI study

BACKGROUND

The development of tailored recovery-oriented strategies in multiple sclerosis requires early identification of an individual's potential for functional recovery.

OBJECTIVE

To identify predictors of visuomotor performance improvements, a proxy of functional recovery, using a predictive statistical model that combines demographic, clinical and magnetic resonance imaging (MRI) data.

METHODS

Right-handed multiple sclerosis patients underwent baseline disability assessment and MRI of the brain structure, function and vascular health. They subsequently undertook 4 weeks of right upper limb visuomotor practice. Changes in performance with practice were our outcome measure. We identified predictors of improvement in a training set of patients using lasso regression; we calculated the best performing model in a validation set and applied this model to a test set.

RESULTS

Patients improved their visuomotor performance with practice. Younger age, better visuomotor abilities, less severe disease burden and concurrent use of preventive treatments predicted improvements. Neuroimaging localised outcome-relevant sensory motor regions, the microstructure and activity of which correlated with performance improvements.

CONCLUSION

Initial characteristics, including age, disease duration, visuo-spatial abilities, hand dexterity, self-evaluated disease impact and the presence of disease-modifying treatments, can predict functional recovery in individual patients, potentially improving their clinical management and stratification in clinical trials. MRI is a correlate of outcome, potentially supporting individual prognosis.

Distinct effects of prematurity on MRI metrics of brain functional connectivity, activity, and structure: Univariate and multivariate analyses

Premature birth affects the developmental trajectory of the brain during a period of intense maturation with possible lifelong consequences. To better understand the effect of prematurity on brain structure and function, we performed blood‐oxygen‐level dependent (BOLD) and anatomical magnetic resonance imaging (MRI) at 40 weeks of postmenstrual age on 88 newborns with variable gestational age (GA) at birth and no evident radiological alterations. We extracted measures of resting‐state functional connectivity and activity in a set of 90 cortical and subcortical brain regions through the evaluation of BOLD correlations between regions and of fractional amplitude of low‐frequency fluctuation (fALFF) within regions, respectively. Anatomical information was acquired through the assessment of regional volumes. We performed univariate analyses on each metric to examine the association with GA at birth, the spatial distribution of the effects, and the consistency across metrics. Moreover, a data‐driven multivariate analysis (i.e., Machine Learning) framework exploited the high dimensionality of the data to assess the sensitivity of each metric to the effect of premature birth. Prematurity was associated with bidirectional alterations of functional connectivity and regional volume and, to a lesser extent, of fALFF. Notably, the effects of prematurity on functional connectivity were spatially diffuse, mainly within cortical regions, whereas effects on regional volume and fALFF were more focal, involving subcortical structures. While the two analytical approaches delivered consistent results, the multivariate analysis was more sensitive in capturing the complex pattern of prematurity effects. Future studies might apply multivariate frameworks to identify premature infants at risk of a negative neurodevelopmental outcome.

Extended-amygdala intrinsic functional connectivity networks: A population study

Pre-clinical and human neuroimaging research implicates the extended-amygdala (ExtA) (including the bed nucleus of the stria terminalis [BST] and central nucleus of the amygdala [CeA]) in networks mediating negative emotional states associated with stress and substance-use behaviours. The extent to which individual ExtA structures form a functionally integrated unit is controversial. We utilised a large sample (n > 1,000 healthy young adult humans) to compare the intrinsic functional connectivity networks (ICNs) of the BST and CeA using task-free functional magnetic resonance imaging (fMRI) data from the Human Connectome Project. We assessed whether inter-individual differences within these ICNs were related to two principal components representing negative disposition and alcohol use. Building on recent primate evidence, we tested whether within BST-CeA intrinsic functional connectivity (iFC) was heritable and further examined co-heritability with our principal components. We demonstrate the BST and CeA to have discrete, but largely overlapping ICNs similar to previous findings. We found no evidence that within BST-CeA iFC was heritable; however, post hoc analyses found significant BST iFC heritability with the broader superficial and centromedial amygdala regions. There were no significant correlations or co-heritability associations with our principal components either across the ICNs or for specific BST-Amygdala iFC. Possible differences in phenotype associations across task-free, task-based, and clinical fMRI are discussed, along with suggestions for more causal investigative paradigms that make use of the now well-established ExtA ICNs.

Evaluation of denoising strategies for task-based functional connectivity: Equalizing residual motion artifacts between rest and cognitively demanding tasks

In‐scanner head motion represents a major confounding factor in functional connectivity studies and it raises particular concerns when motion correlates with the effect of interest. One such instance regards research focused on functional connectivity modulations induced by sustained cognitively demanding tasks. Indeed, cognitive engagement is generally associated with substantially lower in‐scanner movement compared with unconstrained, or minimally constrained, conditions. Consequently, the reliability of condition‐dependent changes in functional connectivity relies on effective denoising strategies. In this study, we evaluated the ability of common denoising pipelines to minimize and balance residual motion‐related artifacts between resting‐state and task conditions. Denoising pipelines—including realignment/tissue‐based regression, PCA/ICA‐based methods (aCompCor and ICA‐AROMA, respectively), global signal regression, and censoring of motion‐contaminated volumes—were evaluated according to a set of benchmarks designed to assess either residual artifacts or network identifiability. We found a marked heterogeneity in pipeline performance, with many approaches showing a differential efficacy between rest and task conditions. The most effective approaches included aCompCor, optimized to increase the noise prediction power of the extracted confounding signals, and global signal regression, although both strategies performed poorly in mitigating the spurious distance‐dependent association between motion and connectivity. Censoring was the only approach that substantially reduced distance‐dependent artifacts, yet this came at the great cost of reduced network identifiability. The implications of these findings for best practice in denoising task‐based functional connectivity data, and more generally for resting‐state data, are discussed.

Cardiorespiratory fitness is associated with increased middle cerebral arterial compliance and decreased cerebral blood flow in young healthy adults: A pulsed ASL MRI study

Cardiorespiratory fitness is thought to have beneficial effects on systemic vascular health, in part, by decreasing arterial stiffness. However, in the absence of non-invasive methods, it remains unknown whether this effect extends to the cerebrovasculature. The present study uses a novel pulsed arterial spin labelling (pASL) technique to explore the relationship between cardiorespiratory fitness and arterial compliance of the middle cerebral arteries (MCAC). Other markers of cerebrovascular health, including resting cerebral blood flow (CBF) and cerebrovascular reactivity to CO2 (CVRCO2) were also investigated. Eleven healthy males aged 21 ± 2 years with varying levels of cardiorespiratory fitness (maximal oxygen uptake (V · O2MAX) 38-76 ml/min/kg) underwent MRI scanning at 3 Tesla. Higher V · O2MAX was associated with greater MCAC (R2 = 0.64, p < 0.01) and lower resting grey matter CBF (R2 = 0.75, p < 0.01). However, V · O2MAX was not predictive of global grey matter BOLD-based CVR (R2 = 0.47, p = 0.17) or CBF-based CVR (R2 = 0.19, p = 0.21). The current experiment builds upon the established benefits of exercise on arterial compliance in the systemic vasculature, by showing that increased cardiorespiratory fitness is associated with greater cerebral arterial compliance in early adulthood.

Graded fMRI Neurofeedback Training of Motor Imagery in Middle Cerebral Artery Stroke Patients: A Preregistered Proof-of-Concept Study

Ischemic stroke of the middle cerebral artery (MCA), a major brain vessel that supplies the primary motor and premotor cortex, is one of the most common causes for severe upper limb impairment. Currently available motor rehabilitation training largely lacks satisfying efficacy with over 70% of stroke survivors showing residual upper limb dysfunction. Motor imagery-based functional magnetic resonance imaging neurofeedback (fMRI-NF) has been suggested as a potential therapeutic technique to improve motor impairment in stroke survivors. In this preregistered proof-of-concept study (https://osf.io/y69jc/), we translated graded fMRI-NF training, a new paradigm that we have previously studied in healthy participants, to first-time MCA stroke survivors with residual mild to severe impairment of upper limb motor function. Neurofeedback was provided from the supplementary motor area (SMA) targeting two different neurofeedback target levels (low and high). We hypothesized that MCA stroke survivors will show (1) sustained SMA-region of interest (ROI) activation and (2) a difference in SMA-ROI activation between low and high neurofeedback conditions during graded fMRI-NF training. At the group level, we found only anecdotal evidence for these preregistered hypotheses. At the individual level, we found anecdotal to moderate evidence for the absence of the hypothesized graded effect for most subjects. These null findings are relevant for future attempts to employ fMRI-NF training in stroke survivors. The study introduces a Bayesian sequential sampling plan, which incorporates prior knowledge, yielding higher sensitivity. The sampling plan was preregistered together with a priori hypotheses and all planned analysis before data collection to address potential publication/researcher biases. Unforeseen difficulties in the translation of our paradigm to a clinical setting required some deviations from the preregistered protocol. We explicitly detail these changes, discuss the accompanied additional challenges that can arise in clinical neurofeedback studies, and formulate recommendations for how these can be addressed. Taken together, this work provides new insights about the feasibility of motor imagery-based graded fMRI-NF training in MCA stroke survivors and serves as a first example for comprehensive study preregistration of an (fMRI) neurofeedback experiment.

Altered cerebrovascular response to acute exercise in patients with Huntington's disease

The objective of this study was to determine whether a single session of exercise was sufficient to induce cerebral adaptations in individuals with Huntington’s disease and to explore the time dynamics of any acute cerebrovascular response. In this case–control study, we employed arterial-spin labelling MRI in 19 Huntington’s disease gene-positive participants (32–65 years, 13 males) and 19 controls (29–63 years, 10 males) matched for age, gender, body mass index and self-reported activity levels, to measure global and regional perfusion in response to 20 min of moderate-intensity cycling. Cerebral perfusion was measured at baseline and 15, 40 and 60 min after exercise cessation. Relative to baseline, we found that cerebral perfusion increased in patients with Huntington’s disease yet was unchanged in control participants in the precentral gyrus (P = 0.016), middle frontal gyrus (P = 0.046) and hippocampus (P = 0.048) 40 min after exercise cessation (+15 to +32.5% change in Huntington’s disease participants, −7.7 to 0.8% change in controls). The length of the disease‐causing trinucleotide repeat expansion in the huntingtin gene predicted the change in the precentral gyrus (P = 0.03) and the intensity of the exercise intervention predicted hippocampal perfusion change in Huntington’s disease participants (P < 0.001). In both groups, exercise increased hippocampal blood flow 60 min after exercise cessation (P = 0.039). These findings demonstrate the utility of acute exercise as a clinically sensitive experimental paradigm to modulate the cerebrovasculature. Twenty minutes of aerobic exercise induced transient cerebrovascular adaptations in the hippocampus and cortex selectively in Huntington’s disease participants and likely represents latent neuropathology not evident at rest.

Brain Network Modularity During a Sustained Working-Memory Task

Spontaneous oscillations of the blood oxygenation level-dependent (BOLD) signal are spatially synchronized within specific brain networks and are thought to reflect synchronized brain activity. Networks are modulated by the performance of a task, even if the exact features and degree of such modulations are still elusive. The presence of networks showing anticorrelated fluctuations lend initially to suppose that a competitive relationship between the default mode network (DMN) and task positive networks (TPNs) supports the efficiency of brain processing. However, more recent results indicate that cooperative and competitive dynamics between networks coexist during task performance. In this study, we used graph analysis to assess the functional relevance of the topological reorganization of brain networks ensuing the execution of a steady state working-memory (WM) task. Our results indicate that the performance of an auditory WM task is associated with a switching between different topological configurations of several regions of specific networks, including frontoparietal, ventral attention, and dorsal attention areas, suggesting segregation of ventral attention regions in the presence of increased overall integration. However, the correct execution of the task requires integration between components belonging to all the involved networks.

Most Small Cerebral Cortical Veins Demonstrate Significant Flow Pulsatility: A Human Phase Contrast MRI Study at 7T

Phase contrast MRI (pcMRI) has been used to investigate flow pulsatility in cerebral arteries, larger cerebral veins, and the cerebrospinal fluid (CSF). Such measurements of intracranial pulsatility and compliance are beginning to inform understanding of the pathophysiology of conditions including normal pressure hydrocephalus, multiple sclerosis, and dementias. We demonstrate the presence of flow pulsatility in small cerebral cortical veins, for the first time using pcMRI at 7 T, with the aim of improving our understanding of the hemodynamics of this little-studied vascular compartment. A method for establishing where venous flow is pulsatile is introduced, revealing significant pulsatility in 116 out of 146 veins, across eight healthy participants, assessed in parietal and frontal regions. Distributions of pulsatility index (PI) and pulse waveform delay were characterized, indicating a small, but statistically significant (p < 0.05), delay of 59 ± 41 ms in cortical veins with respect to the superior sagittal sinus, but no differences between veins draining different arterial supply territories. Measurements of pulsatility in smaller cortical veins, a hitherto unstudied compartment closer to the capillary bed, could lead to a better understanding of intracranial compliance and cerebrovascular (patho)physiology.

Cerebral Metabolic Changes During Visuomotor Adaptation Assessed Using Quantitative fMRI

The brain retains a lifelong ability to adapt through learning and in response to injury or disease-related damage, a process known as functional neuroplasticity. The neural energetics underlying functional brain plasticity have not been thoroughly investigated experimentally in the healthy human brain. A better understanding of the blood flow and metabolic changes that accompany motor skill acquisition, and which facilitate plasticity, is needed before subsequent translation to treatment interventions for recovery of function in disease. The aim of the current study was to characterize cerebral blood flow (CBF) and oxygen consumption (relative CMRO2) responses, using calibrated fMRI conducted in 20 healthy participants, during performance of a serial reaction time task which induces rapid motor adaptation. Regions of interest (ROIs) were defined from areas showing task-induced BOLD and CBF responses that decreased over time. BOLD, CBF and relative CMRO2 responses were calculated for each block of the task. Motor and somatosensory cortices and the cerebellum showed statistically significant positive responses to the task compared to baseline, but with decreasing amplitudes of BOLD, CBF, and CMRO2 response as the task progressed. In the cerebellum, there was a sustained positive BOLD response in the absence of a significant CMRO2 increase from baseline, for all but the first task blocks. This suggests that the brain may continue to elevate the supply energy even after CMRO2 has returned to near baseline levels. Relying on BOLD fMRI data alone in studies of plasticity may not reveal the nature of underlying metabolic responses and their changes over time. Calibrated fMRI approaches may offer a more complete picture of the energetic changes supporting plasticity and learning.

Assessment of the Effects of Aerobic Fitness on Cerebrovascular Function in Young Adults Using Multiple Inversion Time Arterial Spin Labeling MRI

This cross-sectional study investigated the effects of aerobic fitness on cerebrovascular function in the healthy brain. Gray matter cerebral blood flow (CBF) and cerebrovascular reactivity (CVR) were quantified in a sample of young adults within a normal fitness range. Based on existing Transcranial Doppler ultrasound and fMRI evidence, we predicted a positive relationship between fitness and resting gray matter CBF and CVR. Exploratory hypotheses that higher V . O2peak would be associated with higher GM volume and cognitive performance were also investigated. 20 adults underwent a V . O2peak test and a battery of cognitive tests. All subjects also underwent an MRI scan where multiple inversion time (MTI) pulsed arterial spin labeling (PASL) was used to quantify resting CBF and CVR to 5% CO2. Region of interest analysis showed a non-significant inverse correlation between whole-brain gray matter CBF and V . O2peak; r = -0.4, p = 0.08, corrected p (p') = 0.16 and a significant positive correlation between V . O2peak and whole-brain averaged gray matter CVR; r = 0.62, p = 0.003, p' = 0.006. Voxel-wise analysis revealed a significant inverse association between V . O2peak and resting CBF in the left and right thalamus, brainstem, right lateral occipital cortex, left intra-calcarine cortex and cerebellum. The results of this study suggest that aerobic fitness is associated with lower baseline CBF and greater CVR in young adults.

Specific absorption rate implications of within-scan patient head motion for ultra-high field MRI

PURPOSE

This study investigates the implications of all degrees of freedom of within-scan patient head motion on patient safety.

METHODS

Electromagnetic simulations were performed by displacing and/or rotating a virtual body model inside an 8-channel transmit array to simulate 6 degrees of freedom of motion. Rotations of up to 20° and displacements of up to 20 mm including off-axis axial/coronal translations were investigated, yielding 104 head positions. Quadrature excitation, RF shimming, and multi-spoke parallel-transmit excitation pulses were designed for axial slice-selection at 7T, for seven slices across the head. Variation of whole-head specific absorption rate (SAR) and 10-g averaged local SAR of the designed pulses, as well as the change in the maximum eigenvalue (worst-case pulse) were investigated by comparing off-center positions to the central position.

RESULTS

In their respective worst-cases, patient motion increased the eigenvalue-based local SAR by 42%, whole-head SAR by 60%, and the 10-g averaged local SAR by 210%. Local SAR was observed to be more sensitive to displacements along right-left and anterior-posterior directions than displacement in the superior-inferior direction and rotation.

CONCLUSION

This is the first study to investigate the effect of all 6 degrees of freedom of motion on safety of practical pulses. Although the results agree with the literature for overlapping cases, the results demonstrate higher increases (up to 3.1-fold) in local SAR for off-axis displacement in the axial plane, which had received less attention in the literature. This increase in local SAR could potentially affect the local SAR compliance of subjects, unless realistic within-scan patient motion is taken into account during pulse design.

Retrograde blood flow in the internal jugular veins of humans with hypertension may have implications for cerebral arterial blood flow

OBJECTIVES

To use multi-parametric magnetic resonance imaging (MRI) to test the hypothesis that hypertensives would have higher retrograde venous blood flow (RVBF) in the internal jugular veins (IJV) vs. normotensives, and that this would inversely correlate with arterial inflow and gray matter, white matter, and cerebrospinal fluid volumes.

METHODS

Following local institutional review board approval and written consent, a prospective observational 3-T MRI study of 42 hypertensive patients (53 ± 2 years, BMI 28.2 ± 0.6 kg/m², ambulatory daytime systolic BP 148 ± 2 mmHg, ambulatory daytime diastolic BP 101 ± 2 mmHg) and 35 normotensive patients (48 ± 2 years, BMI 25.2 ± 0.8 kg/m², ambulatory daytime systolic BP 119 ± 3 mmHg, ambulatory daytime diastolic BP 90 ± 2 mmHg) was performed. Phase contrast imaging calculated percentage retrograde venous blood flow (%RVBF), brain segmentation estimated regional brain volumes from 3D T1-weighted images, and pseudo-continuous arterial spin labeling measured regional cerebral blood perfusion. Statistical analysis included two-sample equal variance Student's T tests, two-way analysis of variance with Tukey's post hoc correction, and permutation-based two-group general linear modeling (p < 0.05).

RESULTS

In the left IJV, %RVBF was higher in hypertensives (6.1 ± 1.5%) vs. normotensives (1.1 ± 0.3%, p = 0.003). In hypertensives, there was an inverse relationship of %RVBF (permutation-based general linear modeling) to cerebral blood flow in several brain regions, including the left occipital pole and the cerebellar vermis (p < 0.01). Percentage retrograde flow in the left IJV correlated inversely with the total matter volume (gray plus white matter volume) in hypertensives (r = - 0.49, p = 0.004).

CONCLUSION

RVBF in the left IJV is greater in hypertensives vs. normotensives and is linked to regional hypoperfusion and brain total matter volume.

KEY POINTS

• Hypertensive humans have higher retrograde cerebral venous blood flow, associated with regional brain hypoperfusion and lower tissue volume, compared with controls. • Cerebral retrograde venous blood flow may add further stress to already hypoperfused tissue in hypertensive patients. • The amount of retrograde venous blood flow in hypertensive patients may predict which patients might be at higher risk of developing cerebral pathologies.

Whole brain 31 P MRSI at 7T with a dual-tuned receive array

PURPOSE

The design and performance of a novel head coil setup for 31 P spectroscopy at ultra-high field strengths (7T) is presented. The described system supports measurements at both the 1 H and 31 P resonance frequencies.

METHODS

The novel coil consists of 2, actively detunable, coaxial birdcage coils to give homogeneous transmit, combined with a double resonant 30 channel receive array. This allows for anatomical imaging combined with 31 P acquisitions over the whole head, without changing coils or disturbing the subject. A phosphate buffer phantom and 3 healthy volunteers were scanned with a pulse acquire CSI sequence using both the novel array coil and a conventional transceiver birdcage. Four different methods of combining the array channels were compared at 3 different levels of SNR.

RESULTS

The novel coil setup delivers significantly increased 31 P SNR in the peripheral regions of the brain, reaching up to factor 8, while maintaining comparable performance relative to the birdcage in the center.

CONCLUSIONS

The new system offers the potential to acquire whole brain 31 P MRSI with superior signal relative to the standard options.

Mild Propofol Sedation Reduces Frontal Lobe and Thalamic Cerebral Blood Flow: An Arterial Spin Labeling Study

Mechanisms of anesthetic drug-induced sedation and unconsciousness are still incompletely understood. Functional neuroimaging modalities provide a window to study brain function changes during anesthesia allowing us to explore the sequence of neuro-physiological changes associated with anesthesia. Cerebral perfusion change under an assumption of intact neurovascular coupling is an indicator of change in large-scale neural activity. In this experiment, we have investigated resting state cerebral blood flow (CBF) changes in the human brain during mild sedation, with propofol. Arterial spin labeling (ASL) provides a non-invasive, reliable, and robust means of measuring cerebral blood flow (CBF) and can therefore be used to investigate central drug effects. Mild propofol sedation-related CBF changes were studied at rest (n = 15), in a 3 T MR scanner using a PICORE-QUIPSS II ASL technique. CBF was reduced in bilateral paracingulate cortex, premotor cortex, Broca's areas, right superior frontal gyrus and also the thalamus. This cerebral perfusion study demonstrates that propofol induces suppression of key cortical (frontal lobe) and subcortical (thalamus) regions during mild sedation.

Traumatic imagery following glucocorticoid administration in earthquake-related post-traumatic stress disorder: A preliminary functional magnetic resonance imaging study

OBJECTIVE

Post-traumatic stress disorder involves excessive retrieval of traumatic memories. Glucocorticoids impair declarative memory retrieval. This preliminary study examined the effect of acute hydrocortisone administration on brain activation in individuals with earthquake-related post-traumatic stress disorder compared with earthquake-exposed healthy individuals, during retrieval of traumatic memories.

METHOD

Participants exposed to earthquakes with (n = 11) and without post-traumatic stress disorder (n = 11) underwent two functional magnetic resonance imaging scans, 1-week apart, in a double-blind, placebo-controlled, counter-balanced design. On one occasion, they received oral hydrocortisone (20 mg), and on the other, placebo, 1 hour before scanning. Symptom provocation involved script-driven imagery (traumatic and neutral scripts) and measures of self-reported anxiety.

RESULTS

Arterial spin labelling showed that both post-traumatic stress disorder and trauma-exposed controls had significantly reduced cerebral blood flow in response to retrieval of traumatic versus neutral memories in the right hippocampus, parahippocampal gyrus, calcarine sulcus, middle and superior temporal gyrus, posterior cingulate, Heschl's gyrus, inferior parietal lobule, angular gyrus, middle occipital gyrus, supramarginal gyrus, lingual gyrus and cuneus, and the left prefrontal cortex. Hydrocortisone resulted in non-significant trends of increasing subjective distress and reduced regional cerebral blood flow in the left inferior frontal gyrus, left anterior cingulate gyrus, middle temporal gyrus, cerebellum, postcentral gyrus and right frontal pole, during the trauma script.

CONCLUSION

Findings do not fit with some aspects of the accepted neurocircuitry model of post-traumatic stress disorder, i.e., failure of the medial prefrontal cortex to quieten hyperresponsive amygdala activity, and the potential therapeutic benefits of hydrocortisone. They do, however, provide further evidence that exposure to earthquake trauma, regardless of whether post-traumatic stress disorder eventuates, impacts brain activity and highlights the importance of inclusion of trauma-exposed comparisons in studies of post-traumatic stress disorder.

Venous contribution to sodium MRI in the human brain

PURPOSE

Sodium MRI shows great promise as a marker for cerebral metabolic dysfunction in stroke, brain tumor, and neurodegenerative pathologies. However, cerebral blood vessels, whose volume and function are perturbed in these pathologies, have elevated sodium concentrations relative to surrounding tissue. This study aims to assess whether this fluid compartment could bias measurements of tissue sodium using MRI.

METHODS

Density-weighted and B₁ corrected sodium MRI of the brain was acquired in 9 healthy participants at 4.7T. Veins were identified using co-registered ¹ H T 2 ∗ -weighted images and venous partial volume estimates were calculated by down-sampling the finer spatial resolution venous maps from the T 2 ∗ -weighted images to the coarser spatial resolution of the sodium data. Linear regressions of venous partial volume estimates and sodium signal were performed for regions of interest including just gray matter, just white matter, and all brain tissue.

RESULTS

Linear regression demonstrated a significant venous sodium contribution above the underlying tissue signal. The apparent venous sodium concentrations derived from regression were 65.8 ± 4.5 mM (all brain tissue), 71.0 ± 7.4 mM (gray matter), and 55.0 ± 4.7 mM (white matter).

CONCLUSION

Although the partial vein linear regression did not yield the expected sodium concentration in blood (~87 mM), likely the result of point spread function smearing, this regression highlights that blood compartments may bias brain tissue sodium signals across neurological conditions where blood volumes may differ.

Regional Cerebral Activation Accompanies Sympathoexcitation in Women With Polycystic Ovary Syndrome

CONTEXT

Polycystic ovary syndrome (PCOS) is associated with increased sympathetic nervous system activation, but the cerebral pathways involved are unclear.

OBJECTIVE

To compare cerebral [blood oxygen level-dependent (BOLD) functional MRI], pressor [blood pressure (BP), heart rate (HR], and muscle sympathetic nerve activity (MSNA) responses to isometric forearm contraction (IFC) in women with PCOS and matched control subjects.

DESIGN

Case-control study.

SETTING

Referral center.

PARTICIPANTS

Patients with PCOS (n = 20; mean ± SD data: age, 29.8 ± 4.8 years; body mass index (BMI), 26.1 ± 4.9 kg/ m2) and 20 age- and BMI-matched control subjects (age, 29.7 ± 5.0 years; BMI, 26.1 ± 4.8 kg/ m2).

MAIN OUTCOME MEASURES

BP, HR, catecholamine, and MSNA responses to 30% IFC. BOLD signal change was modeled for BP response to 30% IFC.

RESULTS

Although HR and BP increased to a similar extent in both groups after IFC, MSNA burst frequency increased by 68% in the PCOS group compared with 11.9% in control subjects (n = 7 in both groups; P = 0.002). Brain activation indexed by the BOLD signal in response to IFC was significantly greater in the PCOS group (n = 15) compared with controls (n = 15) in the right orbitofrontal cortex (P < 0.0001). Adjustment for insulin sensitivity, but not hyperandrogenism, abolished these between-group differences.

CONCLUSION

Our study confirms enhanced sympathoexcitation in women with PCOS and demonstrates increased regional brain activation in response to IFC. The right orbitofrontal cortex BOLD signal change in women with PCOS is associated with insulin sensitivity. Additional studies are warranted to clarify whether this may offer a novel target for cardiovascular risk reduction.

Comparing MRI metrics to quantify white matter microstructural damage in multiple sclerosis

Quantifying white matter damage in vivo is becoming increasingly important for investigating the effects of neuroprotective and repair strategies in multiple sclerosis (MS). While various approaches are available, the relationship between MRI‐based metrics of white matter microstructure in the disease, that is, to what extent the metrics provide complementary versus redundant information, remains largely unexplored. We obtained four microstructural metrics from 123 MS patients: fractional anisotropy (FA), radial diffusivity (RD), myelin water fraction (MWF), and magnetisation transfer ratio (MTR). Coregistration of maps of these four indices allowed quantification of microstructural damage through voxel‐wise damage scores relative to healthy tissue, as assessed in a group of 27 controls. We considered three white matter tissue‐states, which were expected to vary in microstructural damage: normal appearing white matter (NAWM), T2‐weighted hyperintense lesional tissue without T1‐weighted hypointensity (T2L), and T1‐weighted hypointense lesional tissue with corresponding T2‐weighted hyperintensity (T1L). All MRI indices suggested significant damage in all three tissue‐states, the greatest damage being in T1L. The correlations between indices ranged from r = 0.18 to r = 0.87. MWF was most sensitive when differentiating T2L from NAWM, while MTR was most sensitive when differentiating T1L from NAWM and from T2L. Combining the four metrics into one, through a principal component analysis, did not yield a measure more sensitive to damage than any single measure. Our findings suggest that the metrics are (at least partially) correlated with each other, but sensitive to the different aspects of pathology. Leveraging these differences could be beneficial in clinical trials testing the effects of therapeutic interventions.

Calcium channel blockade with nimodipine reverses MRI evidence of cerebral oedema following acute hypoxia

Acute cerebral hypoxia causes rapid calcium shifts leading to neuronal damage and death. Calcium channel antagonists improve outcomes in some clinical conditions, but mechanisms remain unclear. In 18 healthy participants we: (i) quantified with multiparametric MRI the effect of hypoxia on the thalamus, a region particularly sensitive to hypoxia, and on the whole brain in general; (ii) investigated how calcium channel antagonism with the drug nimodipine affects the brain response to hypoxia. Hypoxia resulted in a significant decrease in apparent diffusion coefficient (ADC), a measure particularly sensitive to cell swelling, in a widespread network of regions across the brain, and the thalamus in particular. In hypoxia, nimodipine significantly increased ADC in the same brain regions, normalizing ADC towards normoxia baseline. There was positive correlation between blood nimodipine levels and ADC change. In the thalamus, there was a significant decrease in the amplitude of low frequency fluctuations (ALFF) in resting state functional MRI and an apparent increase of grey matter volume in hypoxia, with the ALFF partially normalized towards normoxia baseline with nimodipine. This study provides further evidence that the brain response to acute hypoxia is mediated by calcium, and importantly that manipulation of intracellular calcium flux following hypoxia may reduce cerebral cytotoxic oedema.

The BOLD response in primary motor cortex and supplementary motor area during kinesthetic motor imagery based graded fMRI neurofeedback

There is increasing interest in exploring the use of functional MRI neurofeedback (fMRI-NF) as a therapeutic technique for a range of neurological conditions such as stroke and Parkinson's disease (PD). One main therapeutic potential of fMRI-NF is to enhance volitional control of damaged or dysfunctional neural nodes and networks via a closed-loop feedback model using mental imagery as the catalyst of self-regulation. The choice of target node/network and direction of regulation (increase or decrease activity) are central design considerations in fMRI-NF studies. Whilst it remains unclear whether the primary motor cortex (M1) can be activated during motor imagery, the supplementary motor area (SMA) has been robustly activated during motor imagery. Such differences in the regulation potential between primary and supplementary motor cortex are important because these areas can be differentially affected by a stroke or PD, and the choice of fMRI-NF target and grade of self-regulation of activity likely have substantial influence on the clinical effects and cost effectiveness of NF-based interventions. In this study we therefore investigated firstly whether healthy subjects would be able to achieve self-regulation of the hand-representation areas of M1 and the SMA using fMRI-NF training. There was a significant decrease in M1 neural activity during fMRI-NF, whereas SMA neural activity was increased, albeit not with the predicated graded effect. This study has important implications for fMRI-NF protocols that employ motor imagery to modulate activity in specific target regions of the brain and to determine how they may be tailored for neurorehabilitation.

Dual-calibrated fMRI measurement of absolute cerebral metabolic rate of oxygen consumption and effective oxygen diffusivity

Dual-calibrated fMRI is a multi-parametric technique that allows for the quantification of the resting oxygen extraction fraction (OEF), the absolute rate of cerebral metabolic oxygen consumption (CMRO2), cerebral vascular reactivity (CVR) and baseline perfusion (CBF). It combines measurements of arterial spin labelling (ASL) and blood oxygenation level dependent (BOLD) signal changes during hypercapnic and hyperoxic gas challenges. Here we propose an extension to this methodology that permits the simultaneous quantification of the effective oxygen diffusivity of the capillary network (DC). The effective oxygen diffusivity has the scope to be an informative biomarker and useful adjunct to CMRO2, potentially providing a non-invasive metric of microvascular health, which is known to be disturbed in a range of neurological diseases. We demonstrate the new method in a cohort of healthy volunteers (n = 19) both at rest and during visual stimulation. The effective oxygen diffusivity was found to be highly correlated with CMRO2 during rest and activation, consistent with previous PET observations of a strong correlation between metabolic oxygen demand and effective diffusivity. The increase in effective diffusivity during functional activation was found to be consistent with previously reported increases in capillary blood volume, supporting the notion that measured oxygen diffusivity is sensitive to microvascular physiology.

Can Blood Oxygenation Level Dependent Functional Magnetic Resonance Imaging Be Used Accurately to Compare Older and Younger Populations? A Mini Literature Review

A wealth of research has investigated the aging brain using blood oxygenation level dependent functional MRI [Blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI)]. However, many studies do not consider the aging of the cerebrovascular system, which can influence the BOLD signal independently from neural activity, limiting what can be inferred when comparing age groups. Here, we discuss the ways in which the aging neurovascular system can impact BOLD fMRI, the consequences for age-group comparisons and possible strategies for mitigation. While BOLD fMRI is a valuable tool in this context, this review highlights the importance of consideration of vascular confounds.

Modeling of dynamic cerebrovascular reactivity to spontaneous and externally induced CO2 fluctuations in the human brain using BOLD-fMRI

In this work, we investigate the regional characteristics of the dynamic interactions between arterial CO₂ and BOLD (dynamic cerebrovascular reactivity - dCVR) during normal breathing and hypercapnic, externally induced step CO₂ challenges. To obtain dCVR curves at each voxel, we use a custom set of basis functions based on the Laguerre and gamma basis sets. This allows us to obtain robust dCVR estimates both in larger regions of interest (ROIs), as well as in individual voxels. We also implement classification schemes to identify brain regions with similar dCVR characteristics. Our results reveal considerable variability of dCVR across different brain regions, as well as during different experimental conditions (normal breathing and hypercapnic challenges), suggesting a differential response of cerebral vasculature to spontaneous CO₂ fluctuations and larger, externally induced CO₂ changes that are possibly associated with the underlying differences in mean arterial CO₂ levels. The clustering results suggest that anatomically distinct brain regions are characterized by different dCVR curves that in some cases do not exhibit the standard, positive valued curves that have been previously reported. They also reveal a consistent set of dCVR cluster shapes for resting and forcing conditions, which exhibit different distribution patterns across brain voxels.

Scale-invariant rearrangement of resting state networks in the human brain under sustained stimulation

Brain activity at rest is characterized by widely distributed and spatially specific patterns of synchronized low-frequency blood-oxygenation level-dependent (BOLD) fluctuations, which correspond to physiologically relevant brain networks. This network behaviour is known to persist also during task execution, yet the details underlying task-associated modulations of within- and between-network connectivity are largely unknown. In this study we exploited a multi-parametric and multi-scale approach to investigate how low-frequency fluctuations adapt to a sustained n-back working memory task. We found that the transition from the resting state to the task state involves a behaviourally relevant and scale-invariant modulation of synchronization patterns within both task-positive and default mode networks. Specifically, decreases of connectivity within networks are accompanied by increases of connectivity between networks. In spite of large and widespread changes of connectivity strength, the overall topology of brain networks is remarkably preserved. We show that these findings are strongly influenced by connectivity at rest, suggesting that the absolute change of connectivity (i.e., disregarding the baseline) may not be the most suitable metric to study dynamic modulations of functional connectivity. Our results indicate that a task can evoke scale-invariant, distributed changes of BOLD fluctuations, further confirming that low frequency BOLD oscillations show a specialized response and are tightly bound to task-evoked activation.

Investigation Of Interaction Between Physiological Signals And fMRI Dynamic Functional Connectivity Using Independent Component Analysis

The blood oxygen level dependent (BOLD) fMRI signal is influenced not only by neuronal activity but also by fluctuations in physiological signals, including respiration, arterial CO2 and heart rate/ heart rate variability (HR/HRV). Even spontaneous physiological signal fluctuations have been shown to influence the BOLD fMRI signal in a regionally specific manner. Consequently, estimates of functional connectivity between different brain regions, performed when the subject is at rest, may be confounded by the effects of physiological signal fluctuations. In addition, resting functional connectivity has been shown to vary with respect to time (dynamic functional connectivity - DFC), with the sources of this variation not fully elucidated. The effect of physiological factors on dynamic (time-varying) resting-state functional connectivity has not been studied extensively, to our knowledge. In our previous study, we investigated the effect of heart rate (HR) and end-tidal CO2 (PETCO2) on the time-varying network degree of three well-described RSNs (DMN, SMN and Visual Network) using mask-based and seed-based analysis, and we identified brain-heart interactions which were more pronounced in specific frequency bands. Here, we extend this work, by estimating DFC and its corresponding network degree for the RSNs, employing a data-driven approach to extract the RSNs (low-and high-dimensional Independent Component Analysis (ICA)), which we subsequently correlate with the characteristics of simultaneously collected physiological signals. The results confirm that physiological signals have a modulatory effect on resting-state, fMRI-based DFC.

Assessing the repeatability of absolute CMRO2, OEF and haemodynamic measurements from calibrated fMRI

As energy metabolism in the brain is largely oxidative, the measurement of cerebral metabolic rate of oxygen consumption (CMRO2) is a desirable biomarker for quantifying brain activity and tissue viability. Currently, PET techniques based on oxygen isotopes are the gold standard for obtaining whole brain CMRO2 maps. Among MRI techniques that have been developed as an alternative are dual calibrated fMRI (dcFMRI) methods, which exploit simultaneous measurements of BOLD and ASL signals during a hypercapnic-hyperoxic experiment to modulate brain blood flow and oxygenation. In this study we quantified the repeatability of a dcFMRI approach developed in our lab, evaluating its limits and informing its application in studies aimed at characterising the metabolic state of human brain tissue over time. Our analysis focussed on the estimates of oxygen extraction fraction (OEF), cerebral blood flow (CBF), CBF-related cerebrovascular reactivity (CVR) and CMRO2 based on a forward model that describes analytically the acquired dual echo GRE signal. Indices of within- and between-session repeatability are calculated from two different datasets both at a bulk grey matter and at a voxel-wise resolution and finally compared with similar indices obtained from previous MRI and PET measurements. Within- and between-session values of intra-subject coefficient of variation (CVintra) calculated from bulk grey matter estimates 6.7 ± 6.6% (mean ± std.) and 10.5 ± 9.7% for OEF, 6.9 ± 6% and 5.5 ± 4.7% for CBF, 12 ± 9.7% and 12.3 ± 10% for CMRO2. Coefficient of variation (CV) and intraclass correlation coefficient (ICC) maps showed the spatial distribution of the repeatability metrics, informing on the feasibility limits of the method. In conclusion, results show an overall consistency of the estimated physiological parameters with literature reports and a satisfactory level of repeatability considering the higher spatial sensitivity compared to other MRI methods, with varied performance depending on the specific parameter under analysis, on the spatial resolution considered and on the study design.

Calibrated fMRI for mapping absolute CMRO2: Practicalities and prospects

Functional magnetic resonance imaging (fMRI) is an essential workhorse of modern neuroscience, providing valuable insight into the functional organisation of the brain. The physiological mechanisms underlying the blood oxygenation level dependent (BOLD) effect are complex and preclude a straightforward interpretation of the signal. However, by employing appropriate calibration of the BOLD signal, quantitative measurements can be made of important physiological parameters including the absolute rate of cerebral metabolic oxygen consumption or oxygen metabolism (CMRO₂) and oxygen extraction (OEF). The ability to map such fundamental parameters has the potential to greatly expand the utility of fMRI and to broaden its scope of application in clinical research and clinical practice. In this review article we discuss some of the practical issues related to the calibrated-fMRI approach to the measurement of CMRO₂. We give an overview of the necessary precautions to ensure high quality data acquisition, and explore some of the pitfalls and challenges that must be considered as it is applied and interpreted in a widening array of diseases and research questions.

Neurovascular Coupling During Visual Stimulation in Multiple Sclerosis: A MEG-fMRI Study

•

A reduced electrophysiological response to a visual stimulus in MS, characterized by reduced gamma power (30–80 Hz), with MEG.

•

A reduced hemodynamic response to a visual stimulus in MS, characterized by reduced BOLD and CBF responses, with fMRI.

•

The coupling between gamma power and BOLD/CBF was not significantly impaired in the MS group.

The process of neurovascular coupling ensures that increases in neuronal activity are fed by increases in cerebral blood flow. Evidence suggests that neurovascular coupling may be impaired in Multiple Sclerosis (MS) due to a combination of brain hypoperfusion, altered cerebrovascular reactivity and oxygen metabolism, and altered levels of vasoactive compounds.

Here, we tested the hypothesis that neurovascular coupling is impaired in MS. We characterized neurovascular coupling as the relationship between changes in neuronal oscillatory power within the gamma frequency band (30–80 Hz), as measured by magnetoencephalography (MEG), and associated hemodynamic changes (blood oxygenation level dependent, BOLD, and cerebral blood flow, CBF) as measured by functional MRI. We characterized these responses in the visual cortex in 13 MS patients and in 10 matched healthy controls using a reversing checkerboard stimulus at five visual contrasts.

There were no significant group differences in visual acuity, P100 latencies, occipital gray matter (GM) volumes and baseline CBF. However, in the MS patients we found a significant reduction in peak gamma power, BOLD and CBF responses. There were no significant differences in neurovascular coupling between groups, in the visual cortex.

Our results suggest that neuronal and vascular responses are altered in MS. Gamma power reduction could be an indicator of GM dysfunction, possibly mediated by GABAergic changes. Altered hemodynamic responses confirm previous reports of a vascular dysfunction in MS. Despite altered neuronal and vascular responses, neurovascular coupling appears to be preserved in MS, at least within the range of damage and disability studied here.

Estimation of voxel-wise dynamic cerebrovascular reactivity curves from resting-state fMRI data

In this work, we investigate the linear dynamic interactions between fluctuations in arterial CO2 that occur during normal breathing, and the BOLD fMRI signal. We cast this problem within a systems-theoretic framework, where we employ functional expansions for the estimation of the impulse responses in large regions of interest, as well as in individual voxels. We also implement classification schemes in order to identify different brain regions with similar cerebrovascular reactivity characteristics. Our results reveal that it is feasible to obtain reliable estimates of cerebrovascular reactivity curves from resting-state data and that these curves exhibit considerable variability across different brain regions that may be related to the underlying anatomy.

Task-Related Modulations of BOLD Low-Frequency Fluctuations within the Default Mode Network

Spontaneous low-frequency Blood-Oxygenation Level-Dependent (BOLD) signals acquired during resting state are characterized by spatial patterns of synchronous fluctuations, ultimately leading to the identification of robust brain networks. The resting-state brain networks, including the Default Mode Network (DMN), are demonstrated to persist during sustained task execution, but the exact features of task-related changes of network properties are still not well characterized. In this work we sought to examine in a group of 20 healthy volunteers (age 33 ± 6 years, 8 F/12 M) the relationship between changes of spectral and spatiotemporal features of one prominent resting-state network, namely the DMN, during the continuous execution of a working memory n-back task. We found that task execution impacted on both functional connectivity and amplitude of BOLD fluctuations within large parts of the DMN, but these changes correlated between each other only in a small area of the posterior cingulate. We conclude that combined analysis of multiple parameters related to connectivity, and their changes during the transition from resting state to continuous task execution, can contribute to a better understanding of how brain networks rearrange themselves in response to a task.

Graded Hypercapnia-Calibrated BOLD: Beyond the Iso-metabolic Hypercapnic Assumption

Calibrated BOLD is a promising technique that overcomes the sensitivity of conventional fMRI to the cerebrovascular state; measuring either the basal level, or the task-induced response of cerebral metabolic rate of oxygen consumption (CMRO₂). The calibrated BOLD method is susceptible to errors in the measurement of the calibration parameter M, the theoretical BOLD signal change that would occur if all deoxygenated hemoglobin were removed. The original and most popular method for measuring M uses hypercapnia (an increase in arterial CO₂), making the assumption that it does not affect CMRO₂. This assumption has since been challenged and recent studies have used a corrective term, based on literature values of a reduction in basal CMRO₂ with hypercapnia. This is not ideal, as this value may vary across subjects and regions of the brain, and will depend on the level of hypercapnia achieved. Here we propose a new approach, using a graded hypercapnia design and the assumption that CMRO₂ changes linearly with hypercapnia level, such that we can measure M without assuming prior knowledge of the scale of CMRO₂ change. Through use of a graded hypercapnia gas challenge, we are able to remove the bias caused by a reduction in basal CMRO₂ during hypercapnia, whilst simultaneously calculating the dose-wise CMRO₂ change with hypercapnia. When compared with assuming no change in CMRO₂, this approach resulted in significantly lower M-values in both visual and motor cortices, arising from significant dose-dependent hypercapnia reductions in basal CMRO₂ of 1.5 ± 0.6%/mmHg (visual) and 1.8 ± 0.7%/mmHg (motor), where mmHg is the unit change in end-tidal CO₂ level. Variability in the basal CMRO₂ response to hypercapnia, due to experimental differences and inter-subject variability, is accounted for in this approach, unlike previous correction approaches, which use literature values. By incorporating measurement of, and correction for, the reduction in basal CMRO₂ during hypercapnia in the measurement of M-values, application of our approach will correct for an overestimation in both CMRO₂ task-response values and absolute CMRO₂.

Opioid suppression of conditioned anticipatory brain responses to breathlessness

Opioid painkillers are a promising treatment for chronic breathlessness, but are associated with potentially fatal side effects. In the treatment of breathlessness, their mechanisms of action are unclear. A better understanding might help to identify safer alternatives. Learned associations between previously neutral stimuli (e.g. stairs) and repeated breathlessness induce an anticipatory threat response that may worsen breathlessness, contributing to the downward spiral of decline seen in clinical populations. As opioids are known to influence associative learning, we hypothesized that they may interfere with the brain processes underlying a conditioned anticipatory response to breathlessness in relevant brain areas, including the amygdala and the hippocampus.

Healthy volunteers viewed visual cues (neutral stimuli) immediately before induction of experimental breathlessness with inspiratory resistive loading. Thus, an association was formed between the cue and breathlessness. Subsequently, this paradigm was repeated in two identical neuroimaging sessions with intravenous infusions of either low-dose remifentanil (0.7 ng/ml target-controlled infusion) or saline (randomised).

During saline infusion, breathlessness anticipation activated the right anterior insula and the adjacent operculum. Breathlessness was associated with activity in a network including the insula, operculum, dorsolateral prefrontal cortex, anterior cingulate cortex and the primary sensory and motor cortices.

Remifentanil reduced breathlessness unpleasantness but not breathlessness intensity. Remifentanil depressed anticipatory activity in the amygdala and the hippocampus that correlated with reductions in breathlessness unpleasantness. During breathlessness, remifentanil decreased activity in the anterior insula, anterior cingulate cortex and sensory motor cortices. Remifentanil-induced reduction in breathlessness unpleasantness was associated with increased activity in the rostral anterior cingulate cortex and nucleus accumbens, components of the endogenous opioid system known to decrease the perception of aversive stimuli.

These findings suggest that in addition to effects on brainstem respiratory control, opioids palliate breathlessness through an interplay of altered associative learning mechanisms. These mechanisms provide potential targets for novel ways to develop and assess treatments for chronic breathlessness.

•

The mechanisms of how low-dose opioids relieve breathlessness are unknown.

•

We tested whether low-dose opioids affect conditioned anticipation and perception of breathlessness.

•

Low-dose opioids reduced unpleasantness, but not intensity of breathlessness.

•

Reduced breathlessness unpleasantness was associated with activation of the endogenous opioid system.

•

Breathlessness relief was predicted by decreased anticipatory activity in amygdala/hippocampus.