The development and future of perfusion fMRI for dynamic imaging of human brain activity

Coupling of cerebral blood flow and functional connectivity is decreased in healthy aging

Aging leads to cerebral perfusion and functional connectivity changes that have been assessed using various neuroimaging techniques. In addition, a link between these two parameters has been demonstrated in healthy young adults. In this work, we employed arterial spin labeling (ASL) fMRI to measure global and voxel-wise differences in cerebral blood flow (CBF) and intrinsic connectivity contrast (ICC) in the resting state in a group of cognitively normal elderly subjects and a group of cognitively normal young subjects, in order to assess the effects of aging on CBF-ICC coupling, which had not been previously evaluated. Our results showed age-related global and regional CBF decreases in prefrontal mesial areas, lateral frontal regions, insular cortex, lateral parietal areas, precuneus and occipital regions. Subcortically, perfusion was reduced in the medial thalamus and caudate nucleus. ICC was also found reduced with age in prefrontal cortical areas and insular cortex, affecting key nodes of the default mode and salience networks. Areas of ICC and CBF decrease partially overlapped, however, the CBF reduction was more extensive and encompassed more areas. This dissociation was accompanied by a decrease in CBF-ICC coupling. These results suggest that aging leads to a disruption in the relationship between CBF and intrinsic functional connectivity that could be due to neurovascular dysregulation.

Advances in simultaneous PET/MR for imaging neuroreceptor function

Hybrid imaging using PET/MRI has emerged as a platform for elucidating novel neurobiology, molecular and functional changes in disease, and responses to physiological or pharmacological interventions. For the central nervous system, PET/MRI has provided insights into biochemical processes, linking selective molecular targets and distributed brain function. This review highlights several examples that leverage the strengths of simultaneous PET/MRI, which includes measuring the perturbation of multi-modal imaging signals on dynamic timescales during pharmacological challenges, physiological interventions or behavioral tasks. We discuss important considerations for the experimental design of dynamic PET/MRI studies and data analysis approaches for comparing and quantifying simultaneous PET/MRI data. The primary focus of this review is on functional PET/MRI studies of neurotransmitter and receptor systems, with an emphasis on the dopamine, opioid, serotonin and glutamate systems as molecular neuromodulators. In this context, we provide an overview of studies that employ interventions to alter the activity of neuroreceptors or the release of neurotransmitters. Overall, we emphasize how the synergistic use of simultaneous PET/MRI with appropriate study design and interventions has the potential to expand our knowledge about the molecular and functional dynamics of the living human brain. Finally, we give an outlook on the future opportunities for simultaneous PET/MRI.

Improved sensitivity and temporal resolution in perfusion FMRI using velocity selective inversion ASL

PURPOSE

This work aims to investigate the utility of velocity selective inversion pulses for perfusion weighted functional MRI.

METHODS

Tracer kinetic properties of velocity selective inversion (VSI) pulses as an input function for an arterial spin labeling (ASL) experiment were characterized in a group of healthy participants. Numerical simulations were conducted to search for a robust set of timing parameters for FMRI time series acquisition with maximal signal to noise ratio efficiency. The performance of three VSI pulse sequences with different timing parameters was compared with a pseudocontinuous ASL sequence in a simple FMRI experiment conducted on healthy participants.

RESULTS

The fit to the tracer kinetic model yielded arterial CBV of 1.24% ± 0.52% and 0.45 ± 0.11% and perfusion rates of 60.8 ± 32.3 and 34.4 ± 5.4 mL/min/100 g for gray and white matter, respectively. Bolus arrival times were estimated as 75.7 ± 21 ms and 349 ± 78 ms for gray and white matter, respectively. The FMRI experiments showed that VSI pulses yield comparable sensitivity to PCASL with similar timing parameters (TR = 4 s). However, VSI pulses could be used at a faster acquisition speed (TR = 3 s) and were more sensitive to neuronal activity than PCASL pulses, as evidenced by the 31% higher Z scores obtained on average in the active regions.

CONCLUSION

VSI pulses can be very beneficial for perfusion weighted functional MRI because of their tracer kinetic characteristics, which allow a faster acquisition rate while maintaining an efficient labeling input function.

Brain Regional Blood Flow and Working Memory Performance Predict Change in Blood Pressure Over 2 Years

Hypertension is a presumptive risk factor for premature cognitive decline. However, lowering blood pressure (BP) does not uniformly reverse cognitive decline, suggesting that high BP per se may not cause cognitive decline. We hypothesized that essential hypertension has initial effects on the brain that, over time, manifest as cognitive dysfunction in conjunction with both brain vascular abnormalities and systemic BP elevation. Accordingly, we tested whether neuropsychological function and brain blood flow responses to cognitive challenges among prehypertensive individuals would predict subsequent progression of BP. Midlife adults (n=154; mean age, 49; 45% men) with prehypertensive BP underwent neuropsychological testing and assessment of regional cerebral blood flow (rCBF) response to cognitive challenges. Neuropsychological performance measures were derived for verbal and logical memory (memory), executive function, working memory, mental efficiency, and attention. A pseudo-continuous arterial spin labeling magnetic resonance imaging sequence compared rCBF responses with control and active phases of cognitive challenges. Brain areas previously associated with BP were grouped into composites for frontoparietal, frontostriatal, and insular-subcortical rCBF areas. Multiple regression models tested whether BP after 2 years was predicted by initial BP, initial neuropsychological scores, and initial rCBF responses to cognitive challenge. The neuropsychological composite of working memory (standardized beta, -0.276; se=0.116; P=0.02) and the frontostriatal rCBF response to cognitive challenge (standardized beta, 0.234; se=0.108; P=0.03) significantly predicted follow-up BP. Initial BP failed to significantly predict subsequent cognitive performance or rCBF. Changes in brain function may precede or co-occur with progression of BP toward hypertensive levels in midlife.

Multi-echo fMRI: A review of applications in fMRI denoising and analysis of BOLD signals

In recent years the field of fMRI research has enjoyed expanded technical abilities related to resolution, as well as use across many fields of brain research. At the same time, the field has also dealt with uncertainty related to many known and unknown effects of artifact in fMRI data. In this review we discuss an emerging fMRI technology, called multi-echo (ME)-fMRI, which focuses on improving the fidelity and interpretability of fMRI. Where the essential problem of standard single-echo fMRI is the indeterminacy of sources of signals, whether BOLD or artifact, this is not the case for ME-fMRI. By acquiring multiple echo images per slice, the ME approach allows T₂* decay to be modeled at every voxel at every time point. Since BOLD signals arise by changes in T₂* over time, an fMRI experiment sampling the T₂* signal decay can be analyzed to distinguish BOLD from artifact signal constituents. While the ME approach has a long history of use in theoretical and validation studies, modern MRI systems enable whole-brain multi-echo fMRI at high resolution. This review covers recent multi-echo fMRI acquisition methods, and the analysis steps for this data to make fMRI at once more principled, straightforward, and powerful. After a brief overview of history and theory, T₂* modeling and applications will be discussed. These applications include T₂* mapping and combining echoes from ME data to increase BOLD contrast and mitigate dropout artifacts. Next, the modeling of fMRI signal changes to detect signal origins in BOLD-related T₂* versus artifact-related S₀ changes will be reviewed. A focus is on the use of ME-fMRI data to extract and classify components from spatial ICA, called multi-echo ICA (ME-ICA). After describing how ME-fMRI and ME-ICA lead to a general model for analysis of fMRI signals, applications in animal and human imaging will be discussed. Applications include removing motion artifacts in resting state data at subject and group level. New imaging methods such as multi-band multi-echo fMRI and imaging at 7T are demonstrated throughout the review, and a practical analysis pipeline is described. The review culminates with evidence from recent studies of major boosts in statistical power from using multi-echo fMRI for detecting activation and connectivity in healthy individuals and patients with neuropsychiatric disease. In conclusion, the review shows evidence that the multi-echo approach expands the range of experiments that is practicable using fMRI. These findings suggest a compelling future role of the multi-echo approach in subject-level and clinical fMRI.

Sex differences in brain and behavior in adolescence: Findings from the Philadelphia Neurodevelopmental Cohort

Sex differences in brain and behavior were investigated across the lifespan. Parameters include neurobehavioral measures linkable to neuroanatomic and neurophysiologic indicators of brain structure and function. Sexual differentiation of behavior has been related to organizational factors during sensitive periods of development, with adolescence and puberty gaining increased attention. Adolescence is a critical developmental period where transition to adulthood is impacted by multiple factors that can enhance vulnerability to brain dysfunction. Here we highlight sex differences in neurobehavioral measures in adolescence that are linked to brain function. We summarize neuroimaging studies examining brain structure, connectivity and perfusion, underscoring the relationship to sex differences in behavioral measures and commenting on hormonal findings. We focus on relevant data from the Philadelphia Neurodevelopmental Cohort (PNC), a community-based sample of nearly 10,000 clinically and neurocognitively phenotyped youths age 8-21 of whom 1600 have received multimodal neuroimaging. These data indicate early and pervasive sexual differentiation in neurocognitive measures that is linkable to brain parameters. We conclude by describing possible clinical implications.

Multi-regional investigation of the relationship between functional MRI blood oxygenation level dependent (BOLD) activation and GABA concentration

Several recent studies have reported an inter-individual correlation between regional GABA concentration, as measured by MRS, and the amplitude of the functional blood oxygenation level dependent (BOLD) response in the same region. In this study, we set out to investigate whether this coupling generalizes across cortex. In 18 healthy participants, we performed edited MRS measurements of GABA and BOLD-fMRI experiments using regionally related activation paradigms. Regions and tasks were the: occipital cortex with a visual grating stimulus; auditory cortex with a white noise stimulus; sensorimotor cortex with a finger-tapping task; frontal eye field with a saccade task; and dorsolateral prefrontal cortex with a working memory task. In contrast to the prior literature, no correlation between GABA concentration and BOLD activation was detected in any region. The origin of this discrepancy is not clear. Subtle differences in study design or insufficient power may cause differing results; these and other potential reasons for the discrepant results are discussed. This negative result, although it should be interpreted with caution, has a larger sample size than prior positive results, and suggests that the relationship between GABA and the BOLD response may be more complex than previously thought.

The pediatric template of brain perfusion

Magnetic resonance imaging (MRI) captures the dynamics of brain development with multiple modalities that quantify both structure and function. These measurements may yield valuable insights into the neural patterns that mark healthy maturation or that identify early risk for psychiatric disorder. The Pediatric Template of Brain Perfusion (PTBP) is a free and public neuroimaging resource that will help accelerate the understanding of childhood brain development as seen through the lens of multiple modality neuroimaging and in relation to cognitive and environmental factors. The PTBP uses cross-sectional and longitudinal MRI to quantify cortex, white matter, resting state functional connectivity and brain perfusion, as measured by Arterial Spin Labeling (ASL), in 120 children 7-18 years of age. We describe the PTBP and show, as a demonstration of validity, that global summary measurements capture the trajectories that demarcate critical turning points in brain maturation. This novel resource will allow a more detailed understanding of the network-level, structural and functional landmarks that are obtained during normal adolescent brain development.

Separating slow BOLD from non-BOLD baseline drifts using multi-echo fMRI

The functional magnetic resonance (fMRI) baseline is known to drift over the course of an experiment and is often attributed to hardware instability. These ultraslow fMRI fluctuations are inseparable from blood oxygenation level dependent (BOLD) changes in standard single echo fMRI and they are therefore typically removed before further analysis in both resting-state and task paradigms. However, some part of these fluctuations may be of neuronal origin, as neural activity can indeed fluctuate at the scale of several minutes or even longer, such as after the administration of drugs or during the ultradian rhythms. Here, we show that it is possible to separate the slow BOLD and non-BOLD drifts automatically using multi-echo fMRI and multi-echo independent components analysis (ME-ICA) denoising by demonstrating the detection of a visual signal evoked from a flickering checkerboard with slowly changing contrast.

Data collection and analysis strategies for phMRI

Although functional MRI traditionally has been applied mainly to study changes in task-induced brain function, evolving acquisition methodologies and improved knowledge of signal mechanisms have increased the utility of this method for studying responses to pharmacological stimuli, a technique often dubbed "phMRI". The proliferation of higher magnetic field strengths and the use of exogenous contrast agent have boosted detection power, a critical factor for successful phMRI due to the restricted ability to average multiple stimuli within subjects. Receptor-based models of neurovascular coupling, including explicit pharmacological models incorporating receptor densities and affinities and data-driven models that incorporate weak biophysical constraints, have demonstrated compelling descriptions of phMRI signal induced by dopaminergic stimuli. This report describes phMRI acquisition and analysis methodologies, with an emphasis on data-driven analyses. As an example application, statistically efficient data-driven regressors were used to describe the biphasic response to the mu-opioid agonist remifentanil, and antagonism using dopaminergic and GABAergic ligands revealed modulation of the mesolimbic pathway. Results illustrate the power of phMRI as well as our incomplete understanding of mechanisms underlying the signal. Future directions are discussed for phMRI acquisitions in human studies, for evolving analysis methodologies, and for interpretative studies using the new generation of simultaneous PET/MRI scanners. This article is part of the Special Issue Section entitled 'Neuroimaging in Neuropharmacology'.

Impact of puberty on the evolution of cerebral perfusion during adolescence

Puberty is the defining biological process of adolescent development, yet its effects on fundamental properties of brain physiology such as cerebral blood flow (CBF) have never been investigated. Capitalizing on a sample of 922 youths ages 8-22 y imaged using arterial spin labeled MRI as part of the Philadelphia Neurodevelopmental Cohort, we studied normative developmental differences in cerebral perfusion in males and females, as well as specific associations between puberty and CBF. Males and females had conspicuously divergent nonlinear trajectories in CBF evolution with development as modeled by penalized splines. Seventeen brain regions, including hubs of the executive and default mode networks, showed a robust nonlinear age-by-sex interaction that surpassed Bonferroni correction. Notably, within these regions the decline in CBF was similar between males and females in early puberty and only diverged in midpuberty, with CBF actually increasing in females. Taken together, these results delineate sex-specific growth curves for CBF during youth and for the first time to our knowledge link such differential patterns of development to the effects of puberty.

Parameterisation of multi-scale continuum perfusion models from discrete vascular networks

Experimental data and advanced imaging techniques are increasingly enabling the extraction of detailed vascular anatomy from biological tissues. Incorporation of anatomical data within perfusion models is non-trivial, due to heterogeneous vessel density and disparate radii scales. Furthermore, previous idealised networks have assumed a spatially repeating motif or periodic canonical cell, thereby allowing for a flow solution via homogenisation. However, such periodicity is not observed throughout anatomical networks. In this study, we apply various spatial averaging methods to discrete vascular geometries in order to parameterise a continuum model of perfusion. Specifically, a multi-compartment Darcy model was used to provide vascular scale separation for the fluid flow. Permeability tensor fields were derived from both synthetic and anatomically realistic networks using (1) porosity-scaled isotropic, (2) Huyghe and Van Campen, and (3) projected-PCA methods. The Darcy pressure fields were compared via a root-mean-square error metric to an averaged Poiseuille pressure solution over the same domain. The method of Huyghe and Van Campen performed better than the other two methods in all simulations, even for relatively coarse networks. Furthermore, inter-compartment volumetric flux fields, determined using the spatially averaged discrete flux per unit pressure difference, were shown to be accurate across a range of pressure boundary conditions. This work justifies the application of continuum flow models to characterise perfusion resulting from flow in an underlying vascular network.

Subject-specific functional localizers increase sensitivity and functional resolution of multi-subject analyses

One important goal of cognitive neuroscience is to discover and explain properties common to all human brains. The traditional solution for comparing functional activations across brains in fMRI is to align each individual brain to a template brain in a Cartesian coordinate system (e.g., the Montreal Neurological Institute template). However, inter-individual anatomical variability leads to decreases in sensitivity (ability to detect a significant activation when it is present) and functional resolution (ability to discriminate spatially adjacent but functionally different neural responses) in group analyses. Subject-specific functional localizers have been previously argued to increase the sensitivity and functional resolution of fMRI analyses in the presence of inter-subject variability in the locations of functional activations (e.g., Brett et al., 2002; Fedorenko and Kanwisher, 2009, 2011; Fedorenko et al., 2010; Kanwisher et al., 1997; Saxe et al., 2006). In the current paper we quantify this dependence of sensitivity and functional resolution on functional variability across subjects in order to illustrate the highly detrimental effects of this variability on traditional group analyses. We show that analyses that use subject-specific functional localizers usually outperform traditional group-based methods in both sensitivity and functional resolution, even when the same total amount of data is used for each analysis. We further discuss how the subject-specific functional localization approach, which has traditionally only been considered in the context of ROI-based analyses, can be extended to whole-brain voxel-based analyses. We conclude that subject-specific functional localizers are particularly well suited for investigating questions of functional specialization in the brain. An SPM toolbox that can perform all of the analyses described in this paper is publicly available, and the analyses can be applied retroactively to any dataset, provided that multiple runs were acquired per subject, even if no explicit "localizer" task was included.