Differentiating BOLD and non-BOLD signals in fMRI time series from anesthetized rats using multi-echo EPI at 11.7 T

Sleep loss, caffeine, sleep aids and sedation modify brain abnormalities of mild traumatic brain injury

Little evidence exists about how mild traumatic brain injury (mTBI) is affected by commonly encountered exposures of sleep loss, sleep aids, and caffeine that might be potential therapeutic opportunities. In addition, while propofol sedation is administered in severe TBI, its potential utility in mild TBI is unclear. Each of these exposures is known to have pronounced effects on cerebral metabolism and blood flow and neurochemistry. We hypothesized that they each interact with cerebral metabolic dynamics post-injury and change the subclinical characteristics of mTBI. MTBI in rats was produced by head rotational acceleration injury that mimics the biomechanics of human mTBI. Three mTBIs spaced 48 h apart were used to increase the likelihood that vulnerabilities induced by repeated mTBI would be manifested without clinically relevant structural damage. After the third mTBI, rats were immediately sleep deprived or administered caffeine or suvorexant (an orexin antagonist and sleep aid) for the next 24 h or administered propofol for 5 h. Resting state functional magnetic resonance imaging (rs-fMRI) and diffusion tensor imaging (DTI) were performed 24 h after the third mTBI and again after 30 days to determine changes to the brain mTBI phenotype. Multi-modal analyses on brain regions of interest included measures of functional connectivity and regional homogeneity from rs-fMRI, and mean diffusivity (MD) and fractional anisotropy (FA) from DTI. Each intervention changed the mTBI profile of subclinical effects that presumably underlie healing, compensation, damage, and plasticity. Sleep loss during the acute post-injury period resulted in dramatic changes to functional connectivity. Caffeine, propofol sedation and suvorexant were especially noteworthy for differential effects on microstructure in gray and white matter regions after mTBI. The present results indicate that commonplace exposures and short-term sedation alter the subclinical manifestations of repeated mTBI and therefore likely play roles in symptomatology and vulnerability to damage by repeated mTBI.

Decomposing Neural Representational Patterns of Discriminatory and Hedonic Information during Somatosensory Stimulation

The ability to interrogate specific representations in the brain, determining how, and where, difference sources of information are instantiated can provide invaluable insight into neural functioning. Pattern component modeling (PCM) is a recent analytic technique for human neuroimaging that allows the decomposition of representational patterns in brain into contributing subcomponents. In the current study, we present a novel PCM variant that tracks the contribution of prespecified representational patterns to brain representation across areas, thus allowing hypothesis-guided employment of the technique. We apply this technique to investigate the contributions of hedonic and nonhedonic information to the neural representation of tactile experience. We applied aversive pressure (AP) and appetitive brush (AB) to stimulate distinct peripheral nerve pathways for tactile information (C-/CT-fibers, respectively) while patients underwent functional magnetic resonance imaging (fMRI) scanning. We performed representational similarity analyses (RSAs) with pattern component modeling to dissociate how discriminatory versus hedonic tactile information contributes to population code representations in the human brain. Results demonstrated that information about appetitive and aversive tactile sensation is represented separately from nonhedonic tactile information across cortical structures. This also demonstrates the potential of new hypothesis-guided PCM variants to help delineate how information is instantiated in the brain.

Functional Connectivity of the Brain Across Rodents and Humans

Resting-state functional magnetic resonance imaging (rs-fMRI), which measures the spontaneous fluctuations in the blood oxygen level-dependent (BOLD) signal, is increasingly utilized for the investigation of the brain's physiological and pathological functional activity. Rodents, as a typical animal model in neuroscience, play an important role in the studies that examine the neuronal processes that underpin the spontaneous fluctuations in the BOLD signal and the functional connectivity that results. Translating this knowledge from rodents to humans requires a basic knowledge of the similarities and differences across species in terms of both the BOLD signal fluctuations and the resulting functional connectivity. This review begins by examining similarities and differences in anatomical features, acquisition parameters, and preprocessing techniques, as factors that contribute to functional connectivity. Homologous functional networks are compared across species, and aspects of the BOLD fluctuations such as the topography of the global signal and the relationship between structural and functional connectivity are examined. Time-varying features of functional connectivity, obtained by sliding windowed approaches, quasi-periodic patterns, and coactivation patterns, are compared across species. Applications demonstrating the use of rs-fMRI as a translational tool for cross-species analysis are discussed, with an emphasis on neurological and psychiatric disorders. Finally, open questions are presented to encapsulate the future direction of the field.

An isotropic EPI database and analytical pipelines for rat brain resting-state fMRI

Resting-state functional magnetic resonance imaging (fMRI) has drastically expanded the scope of brain research by advancing our knowledge about the topologies, dynamics, and interspecies translatability of functional brain networks. Several databases have been developed and shared in accordance with recent key initiatives in the rodent fMRI community to enhance the transparency, reproducibility, and interpretability of data acquired at various sites. Despite these pioneering efforts, one notable challenge preventing efficient standardization in the field is the customary choice of anisotropic echo planar imaging (EPI) schemes with limited spatial coverage. Imaging with anisotropic resolution and/or reduced brain coverage has significant shortcomings including reduced registration accuracy and increased deviation in brain feature detection. Here we proposed a high-spatial-resolution (0.4 mm), isotropic, whole-brain EPI protocol for the rat brain using a horizontal slicing scheme that can maintain a functionally relevant repetition time (TR), avoid high gradient duty cycles, and offer unequivocal whole-brain coverage. Using this protocol, we acquired resting-state EPI fMRI data from 87 healthy rats under the widely used dexmedetomidine sedation supplemented with low-dose isoflurane on a 9.4 T MRI system. We developed an EPI template that closely approximates the Paxinos and Watson's rat brain coordinate system and demonstrated its ability to improve the accuracy of group-level approaches and streamline fMRI data pre-processing. Using this database, we employed a multi-scale dictionary-learning approach to identify reliable spatiotemporal features representing rat brain intrinsic activity. Subsequently, we performed k-means clustering on those features to obtain spatially discrete, functional regions of interest (ROIs). Using Euclidean-based hierarchical clustering and modularity-based partitioning, we identified the topological organizations of the rat brain. Additionally, the identified group-level FC network appeared robust across strains and sexes. The "triple-network" commonly adapted in human fMRI were resembled in the rat brain. Through this work, we disseminate raw and pre-processed isotropic EPI data, a rat brain EPI template, as well as identified functional ROIs and networks in standardized rat brain coordinates. We also make our analytical pipelines and scripts publicly available, with the hope of facilitating rat brain resting-state fMRI study standardization.

Contribution of animal models toward understanding resting state functional connectivity

Functional connectivity, which reflects the spatial and temporal organization of intrinsic activity throughout the brain, is one of the most studied measures in human neuroimaging research. The noninvasive acquisition of resting state functional magnetic resonance imaging (rs-fMRI) allows the characterization of features designated as functional networks, functional connectivity gradients, and time-varying activity patterns that provide insight into the intrinsic functional organization of the brain and potential alterations related to brain dysfunction. Functional connectivity, hence, captures dimensions of the brain's activity that have enormous potential for both clinical and preclinical research. However, the mechanisms underlying functional connectivity have yet to be fully characterized, hindering interpretation of rs-fMRI studies. As in other branches of neuroscience, the identification of the neurophysiological processes that contribute to functional connectivity largely depends on research conducted on laboratory animals, which provide a platform where specific, multi-dimensional investigations that involve invasive measurements can be carried out. These highly controlled experiments facilitate the interpretation of the temporal correlations observed across the brain. Indeed, information obtained from animal experimentation to date is the basis for our current understanding of the underlying basis for functional brain connectivity. This review presents a compendium of some of the most critical advances in the field based on the efforts made by the animal neuroimaging community.

In vivo evidence of functional disconnection between brainstem monoaminergic nuclei and brain networks in multiple sclerosis

BACKGROUND

brainstem monoaminergic (dopaminergic, noradrenergic, and serotoninergic) nuclei (BrMn) contain a variety of ascending neurons that diffusely project to the whole brain, crucially regulating normal brain function. BrMn are directly affected in multiple sclerosis (MS) by inflammation and neurodegeneration. Moreover, inflammation reduces the synthesis of monoamines. Aberrant monoaminergic neurotransmission contributes to the pathogenesis of MS and explains some clinical features of MS. We used resting-state functional MRI (RS-fMRI) to characterize abnormal patterns of BrMn functional connectivity (FC) in MS.

METHODS

BrMn FC was studied with multi-echo RS-fMRI in n = 68 relapsing-remitting MS patients and n = 39 healthy controls (HC), by performing a seed-based analysis, after producing standard space seed masks of the BrMn. FC was assessed between ventral tegmental area (VTA), locus coeruleus (LC), median raphe (MR), dorsal raphe (DR), and the rest of the brain and compared between MS patients and HC. Between-group comparisons were carried out only within the main effect observed in HC, setting p<0.05 family-wise-error corrected (_FWE).

RESULTS

in HC, VTA displayed FC with the core regions of the default-mode network. As compared to HC, MS patients showed altered FC between VTA and posterior cingulate cortex (p<0.05_FWE). LC displayed FC with core regions of the executive-control network with a reduced functional connection between LC and right prefrontal cortex in MS patients (p<0.05_FWE). Raphe nuclei was functionally connected with cerebellar cortex, with a significantly lower FC between these nuclei and cerebellum in MS patients, as compared to HC (p<0.05_FWE).

CONCLUSIONS

our study demonstrated in MS patients a functional disconnection between BrMn and cortical/subcortical efferent targets of central brain networks, possibly due to a loss or a dysregulation of BrMn neurons. This adds new information about how monoaminergic systems contribute to MS pathogenesis and suggests new potential therapeutic targets.

Echo Time Dependency of Local Activity Metrics of Resting-State Functional MRI

Local activity metrics of resting-state functional MRI (RS-fMRI), such as the amplitude of low-frequency fluctuation (ALFF), fractional ALFF (fALFF), regional homogeneity (ReHo), and degree centrality (DC), are widely used to detect brain abnormalities based on signal fluctuations. Although signal changes with echo time (TE) have been widely studied, the effect of TE on local activity metrics has not been investigated. RS-fMRI datasets from 12 healthy subjects with eyes open (EO) and eyes closed (EC) were obtained with a four-echo gradient-echo-planar imaging pulse sequence with the following parameters: repetition time/TE1/TE2/TE3/TE4 = 2,000/13/30.93/48.86/66.79 ms. Six representative regions were selected for simulating the spatial feature of TE dependency of local activity metrics. Moreover, whole-brain local activity metrics were calculated from each echo dataset and compared between EO and EC conditions. Dice overlap coefficient (DOC) was then employed to calculate the overlap between the T maps. We found that all the local activity metrics displayed different TE dependency characteristics, while their overall change patterns were similar: an initial large change followed by a slow variation. The T maps for local activity metrics also varied greatly with TE. For ALFF, fALFF, ReHo, and DC, the DOCs for voxels in four TE datasets were 6.87, 0.73, 5.08, and 0.93%, respectively. Collectively, these findings demonstrate that local metrics are greatly dependent on TE. Therefore, TE should be carefully considered for the optimization of data acquisition and multi-center data analysis in RS-fMRI.

Cognitive fatigue in multiple sclerosis is associated with alterations in the functional connectivity of monoamine circuits

Fatigue is a highly prevalent and debilitating symptom in multiple sclerosis, but currently the available treatment options have limited efficacy. The development of innovative and efficacious targeted treatments for fatigue in multiple sclerosis has been marred by the limited knowledge of the underlying mechanisms. One of the hypotheses postulates that multiple sclerosis pathology might cause reduced monoaminergic release in the central nervous system with consequences on motivation, mood and attention. Here, we applied the recently developed Receptor-Enriched Analysis of Functional Connectivity by Targets method to investigate whether patients with high and low fatigue differ in the functional connectivity (FC) of the monoamine circuits in the brain. We recruited 55 patients with multiple sclerosis, which were then classified as highly fatigued or mildly fatigued based on their scores on the cognitive sub-scale of the Modified Fatigue Impact scale. We acquired resting-state functional MRI scans and derived individual maps of connectivity associated with the distribution of the dopamine, noradrenaline and serotonin transporters as measured by positron emission tomography. We found that patients with high fatigue present decreased noradrenaline transporter (NAT)-enriched connectivity in several frontal and prefrontal areas when compared to those with lower fatigue. The NAT-enriched FC predicted negatively individual cognitive fatigue scores. Our findings support the idea that alterations in the catecholaminergic functional circuits underlie fatigue in multiple sclerosis and identify the NAT as a putative therapeutic target directed to pathophysiology.

Systematic Review: Anesthetic Protocols and Management as Confounders in Rodent Blood Oxygen Level Dependent Functional Magnetic Resonance Imaging (BOLD fMRI)-Part B: Effects of Anesthetic Agents, Doses and Timing

Simple Summary

To understand brain function in rats and mice functional magnetic resonance imaging of the brain is used. With this type of “brain scan” regional changes in blood flow and oxygen consumption are measured as an indirect surrogate for activity of brain regions. Animals are often anesthetized for the experiments to prevent stress and blurred images due to movement. However, anesthesia may alter the measurements, as blood flow within the brain is differently affected by different anesthetics, and anesthetics also directly affect brain function. Consequently, results obtained under one anesthetic protocol may not be comparable with those obtained under another, and/or not representative for awake animals and humans. We have systematically searched the existing literature for studies analyzing the effects of different anesthesia methods or studies that compared anesthetized and awake animals. Most studies reported that anesthetic agents, doses and timing had an effect on functional magnetic resonance imaging results. To obtain results which promote our understanding of brain function, it is therefore essential that a standard for anesthetic protocols for functional magnetic resonance is defined and their impact is well characterized.

Abstract

In rodent models the use of functional magnetic resonance imaging (fMRI) under anesthesia is common. The anesthetic protocol might influence fMRI readouts either directly or via changes in physiological parameters. As long as those factors cannot be objectively quantified, the scientific validity of fMRI in rodents is impaired. In the present systematic review, literature analyzing in rats and mice the influence of anesthesia regimes and concurrent physiological functions on blood oxygen level dependent (BOLD) fMRI results was investigated. Studies from four databases that were searched were selected following pre-defined criteria. Two separate articles publish the results; the herewith presented article includes the analyses of 83 studies. Most studies found differences in BOLD fMRI readouts with different anesthesia drugs and dose rates, time points of imaging or when awake status was compared to anesthetized animals. To obtain scientifically valid, reproducible results from rodent fMRI studies, stable levels of anesthesia with agents suitable for the model under investigation as well as known and objectively quantifiable effects on readouts are, thus, mandatory. Further studies should establish dose ranges for standardized anesthetic protocols and determine time windows for imaging during which influence of anesthesia on readout is objectively quantifiable.

Rapid Precision Functional Mapping of Individuals Using Multi-Echo fMRI

Resting-state functional magnetic resonance imaging (fMRI) is widely used in cognitive and clinical neuroscience, but long-duration scans are currently needed to reliably characterize individual differences in functional connectivity (FC) and brain network topology. In this report, we demonstrate that multi-echo fMRI can improve the reliability of FC-based measurements. In four densely sampled individual humans, just 10 min of multi-echo data yielded better test-retest reliability than 30 min of single-echo data in independent datasets. This effect is pronounced in clinically important brain regions, including the subgenual cingulate, basal ganglia, and cerebellum, and is linked to three biophysical signal mechanisms (thermal noise, regional variability in the rate of T₂* decay, and S₀-dependent artifacts) with spatially distinct influences. Together, these findings establish the potential utility of multi-echo fMRI for rapid precision mapping using experimentally and clinically tractable scan times and will facilitate longitudinal neuroimaging of clinical populations.

Lynch et al. demonstrate that the test-retest reliability of resting-state connectivity measurements can be improved using multi-echo fMRI. This effect is pronounced in clinically important brain regions and could help facilitate precision mapping of functional brain networks in healthy people and patient populations.

Sustained perturbation in functional connectivity induced by cold pain

BACKGROUND

Functional connectivity (FC) perturbations have been reported in multiple chronic pain phenotypes, but the nature of reported changes varies between cohorts and may relate to the consequences of living with chronic-pain related comorbidities, such as anxiety and depression. Healthy volunteer studies provide opportunities to study the effects of tonic noxious stimulation independently of these sequelae. Connectivity changes in task negative and positive networks, for example, the default mode and salience networks (DMN/SN), respectively, have been described, but how these and other connectivity networks, for example, those governing descending pain control are affected by the presence of tonic, noxious stimulation in healthy, pain-free individuals, remains unknown.

METHOD

In 20 healthy volunteers, we assessed FC prior to, during, and following tonic cold painful stimulation in the ventromedial prefrontal cortex (vmPFC), rostral anterior insula (rAI), subgenual anterior cingulate cortex (ACC) and periaqueductal grey (PAG). We also recorded subjectively reported pain using a computerised visual analogue scale.

RESULTS

We saw DMN FC changes during painful stimulation and that inter-network connectivity between the rAI with the vmPFC increased during pain, whereas PAG-precuneus FC decreased. Pain-induced FC alterations persisted following noxious stimulation. FC changes related to the magnitude of individuals' subjectively reported pain.

CONCLUSIONS

We demonstrate FC changes during and following tonic cold-pain in healthy participants. Similarities between our findings and reports of patients with chronic pain suggest that some FC changes observed in these patients may relate to the presence of an ongoing afferent nociceptive drive.

SIGNIFICANCE

How pain-related resting state networks are affected by tonic cold-pain remains unknown. We investigated functional connectivity alterations during and following tonic cold pain in healthy volunteers. Cold pain perturbed the functional connectivity of the ventro-medial prefrontal cortex, anterior insula, and the periacquaductal grey area. These connectivity changes were associated with the magnitude of individuals' reported pain. We suggest that some connectivity changes described in chronic pain patients may be due to an ongoing afferent peripheral drive.

Reduced fronto-striatal volume in attention-deficit/hyperactivity disorder in two cohorts across the lifespan

Attention-Deficit/Hyperactivity Disorder (ADHD) has been associated with altered brain anatomy in neuroimaging studies. However, small and heterogeneous study samples, and the use of region-of-interest and tissue-specific analyses have limited the consistency and replicability of these effects. We used a data-driven multivariate approach to investigate neuroanatomical features associated with ADHD in two independent cohorts: the Dutch NeuroIMAGE cohort (n = 890, 17.2 years) and the Brazilian IMpACT cohort (n = 180, 44.2 years). Using independent component analysis of whole-brain morphometry images, 375 neuroanatomical components were assessed for association with ADHD. In both discovery (corrected-p = 0.0085) and replication (p = 0.032) cohorts, ADHD was associated with reduced volume in frontal lobes, striatum, and their interconnecting white-matter. Current results provide further evidence for the role of the fronto-striatal circuit in ADHD in children, and for the first time show its relevance to ADHD in adults. The fact that the cohorts are from different continents and comprise different age ranges highlights the robustness of the findings.

Unravelling the effects of methylphenidate on the dopaminergic and noradrenergic functional circuits

Functional magnetic resonance imaging (fMRI) can be combined with drugs to investigate the system-level functional responses in the brain to such challenges. However, most psychoactive agents act on multiple neurotransmitters, limiting the ability of fMRI to identify functional effects related to actions on discrete pharmacological targets. We recently introduced a multimodal approach, REACT (Receptor-Enriched Analysis of functional Connectivity by Targets), which offers the opportunity to disentangle effects of drugs on different neurotransmitters and clarify the biological mechanisms driving clinical efficacy and side effects of a compound. Here, we focus on methylphenidate (MPH), which binds to the dopamine transporter (DAT) and the norepinephrine transporter (NET), to unravel its effects on dopaminergic and noradrenergic functional circuits in the healthy brain at rest. We then explored the relationship between these target-enriched resting state functional connectivity (FC) maps and inter-individual variability in behavioural responses to a reinforcement-learning task encompassing a novelty manipulation to disentangle the molecular systems underlying specific cognitive/behavioural effects. Our main analysis showed a significant MPH-induced FC increase in sensorimotor areas in the functional circuit associated with DAT. In our exploratory analysis, we found that MPH-induced regional variations in the DAT and NET-enriched FC maps were significantly correlated with some of the inter-individual differences on key behavioural responses associated with the reinforcement-learning task. Our findings show that main MPH-related FC changes at rest can be understood through the distribution of DAT in the brain. Furthermore, they suggest that when compounds have mixed pharmacological profiles, REACT may be able to capture regional functional effects that are underpinned by the same cognitive mechanism but are related to engagement of distinct molecular targets.

Local and global consequences of reward-evoked striatal dopamine release

The neurotransmitter dopamine is required for the reinforcement of actions by rewarding stimuli¹. Neuroscientists have tried to define dopamine’s functions in concise conceptual terms², but the practical significance of dopamine release depends on its diverse brain-wide consequences. Although the molecular and cellular effects of dopaminergic signaling have been extensively studied³, its impact on larger-scale neural activity profiles is less understood. Here we combine dynamic dopamine-sensitive molecular imaging⁴ and functional magnetic resonance imaging (fMRI) to determine how striatal dopamine release shapes local and global responses to rewarding stimulation in the rodent brain. We find that dopamine consistently alters the duration but not the magnitude of stimulus responses across much of striatum, via quantifiable postsynaptic effects that vary across subregions. Striatal dopamine release also potentiates a network of distal responses we delineate using neurochemically-dependent functional connectivity analyses. Hot spots of dopaminergic drive notably include cortical regions associated with both limbic and motor function. Our results thus reveal distinct neuromodulatory actions of striatal dopamine that extend well beyond its sites of peak release, and that result in enhanced activation of remote neural populations necessary for performance of motivated actions. Our findings also suggest brain-wide biomarkers of dopaminergic function and could provide a basis for improved interpretation of neuroimaging results relevant to learning and addiction.

Separable neural mechanisms for the pleiotropic association of copy number variants with neuropsychiatric traits

22q11.2, 15q13.3, and 1q21.1 microdeletions attract considerable interest by conferring high risk for a range of neuropsychiatric disorders, including schizophrenia and autism. A fundamental open question is whether divergent or convergent neural mechanisms mediate this genetic pleiotropic association with the same behavioral phenotypes. We use a combination of rodent microdeletion models with high-field neuroimaging to perform a comparative whole-brain characterization of functional and structural mechanisms linked to high-risk states. Resting-state functional and structural magnetic resonance imaging data were acquired on mice carrying heterozygous microdeletions in 22q11.2 (N = 12), 15q13.3 (N = 11), and 1q21.1 (N = 11) loci. We performed network-based statistic, graph, and morphometric analyses. The three microdeletions did not share significant systems-level features. Instead, morphometric analyses revealed microcephaly in 1q21.1 and macrocephaly in 15q13.3 deletions, whereas cerebellar volume was specifically reduced in 22q11.2 deletion. In function, 22q11.2 deletion mice showed widespread cortical hypoconnectivity, accompanied by opposing hyperconnectivity in dopaminergic pathways, which was confirmed by graph analysis. 1q21.1 exhibited distinct changes in posterior midbrain morphology and function, especially in periaqueductal gray, whereas 15q13.3 demonstrated alterations in auditory/striatal system. The combination of cortical hypoconnectivity and dopaminergic hyperconnectivity and reduced cerebellum in 22q11.2 deletion mirrors key neurodevelopmental features of schizophrenia, whereas changes in midbrain and auditory/striatal morphology and topology in 1q21.1 and 15q13.3 rather indicate focal processes possibly linked to the emergence of abnormal salience perception and hallucinations. In addition to insights into pathophysiological processes in these microdeletions, our results establish the general point that microdeletions might increase risk for overlapping neuropsychiatric phenotypes through separable neural mechanisms.

Linking Pain Sensation to the Autonomic Nervous System: The Role of the Anterior Cingulate and Periaqueductal Gray Resting-State Networks

There are bi-directional interactions between the autonomic nervous system (ANS) and pain. This is likely underpinned by a substantial overlap between brain areas of the central autonomic network and areas involved in pain processing and modulation. To date, however, relatively little is known about the neuronal substrates of the ANS-pain association. Here, we acquired resting state fMRI scans in 21 healthy subjects at rest and during tonic noxious cold stimulation. As indicators of autonomic function, we examined how heart rate variability (HRV) frequency measures were influenced by tonic noxious stimulation and how these variables related to participants’ pain perception and to brain functional connectivity in regions known to play a role in both ANS regulation and pain perception, namely the right dorsal anterior cingulate cortex (dACC) and periaqueductal gray (PAG). Our findings support a role of the cardiac ANS in brain connectivity during pain, linking functional connections of the dACC and PAG with measurements of low frequency (LF)-HRV. In particular, we identified a three-way relationship between the ANS, cortical brain networks known to underpin pain processing, and participants’ subjectively reported pain experiences. LF-HRV both at rest and during pain correlated with functional connectivity between the seed regions and other cortical areas including the right dorsolateral prefrontal cortex (dlPFC), left anterior insula (AI), and the precuneus. Our findings link cardiovascular autonomic parameters to brain activity changes involved in the elaboration of nociceptive information, thus beginning to elucidate underlying brain mechanisms associated with the reciprocal relationship between autonomic and pain-related systems.

Animal Functional Magnetic Resonance Imaging: Trends and Path Toward Standardization

Animal whole-brain functional magnetic resonance imaging (fMRI) provides a non-invasive window into brain activity. A collection of associated methods aims to replicate observations made in humans and to identify the mechanisms underlying the distributed neuronal activity in the healthy and disordered brain. Animal fMRI studies have developed rapidly over the past years, fueled by the development of resting-state fMRI connectivity and genetically encoded neuromodulatory tools. Yet, comparisons between sites remain hampered by lack of standardization. Recently, we highlighted that mouse resting-state functional connectivity converges across centers, although large discrepancies in sensitivity and specificity remained. Here, we explore past and present trends within the animal fMRI community and highlight critical aspects in study design, data acquisition, and post-processing operations, that may affect the results and influence the comparability between studies. We also suggest practices aimed to promote the adoption of standards within the community and improve between-lab reproducibility. The implementation of standardized animal neuroimaging protocols will facilitate animal population imaging efforts as well as meta-analysis and replication studies, the gold standards in evidence-based science.

Global brain signal in awake rats

Although often used as a nuisance in resting-state functional magnetic resonance imaging (rsfMRI), the global brain signal in humans and anesthetized animals has important neural basis. However, our knowledge of the global signal in awake rodents is sparse. To bridge this gap, we systematically analyzed rsfMRI data acquired with a conventional single-echo (SE) echo planar imaging (EPI) sequence in awake rats. The spatial pattern of rsfMRI frames during peaks of the global signal exhibited prominent co-activations in the thalamo-cortical and hippocampo-cortical networks, as well as in the basal forebrain, hinting that these neural networks might contribute to the global brain signal in awake rodents. To validate this concept, we acquired rsfMRI data using a multi-echo (ME) EPI sequence and removed non-neural components in the rsfMRI signal. Consistent co-activation patterns were obtained in extensively de-noised ME-rsfMRI data, corroborating the finding from SE-rsfMRI data. Furthermore, during rsfMRI experiments, we simultaneously recorded neural spiking activities in the hippocampus using GCaMP-based fiber photometry. The hippocampal calcium activity exhibited significant correspondence with the global rsfMRI signal. These data collectively suggest that the global rsfMRI signal contains significant neural components that involve coordinated activities in the thalamo-cortical and hippocampo-cortical networks. These results provide important insight into the neural substrate of the global brain signal in awake rodents.

Dependence of resting-state fMRI fluctuation amplitudes on cerebral cortical orientation relative to the direction of B0 and anatomical axes

Functional magnetic resonance imaging (fMRI) is now capable of sub-millimetre scale measurements over the entire human brain, however with such high resolutions each voxel is influenced by the local fine-scale details of the cerebral cortical vascular anatomy. The cortical vasculature is structured with the pial vessels lying tangentially along the grey matter surface, intracortical diving arterioles and ascending venules running perpendicularly to the surface, and a randomly oriented capillary network within the parenchyma. It is well-known that the amplitude of the blood-oxygenation level dependent (BOLD) signal emanating from a vessel depends on its orientation relative to the B0-field. Thus the vascular geometric hierarchy will impart an orientation dependence to the BOLD signal amplitudes and amplitude differences due to orientation differences constitute a bias for interpreting neuronal activity. Here, we demonstrate a clear effect of cortical orientation to B0 in the resting-state BOLD-fMRI amplitude (quantified as the coefficient of temporal signal variation) for 1.1 mm isotropic data at 7T and 2 mm isotropic at 3T. The maximum bias, i.e. the fluctuation amplitude difference between regions where cortex is perpendicular to vs. parallel to B0, is about +70% at the pial surface at 7T and +11% at 3T. The B0 orientation bias declines with cortical depth, becomes progressively smaller closer to the white matter surface, but then increases again to a local maximum within the white matter just beneath the cortical grey matter, suggesting a distinct tangential network of white matter vessels that also generate a BOLD orientation effect. We further found significant (negative) biases with the cortex orientation to the anterior-posterior anatomical axis of the head: a maximum negative bias of about -30% at the pial surface at 7T and about -13% at 3T. The amount of signal variance explained by the low frequency drift, motion and the respiratory cycle also showed a cortical orientation dependence; only the cardiac cycle induced signal variance was independent of cortical orientation, suggesting that the cardiac induced component of the image time-series fluctuations is not related to a significant change in susceptibility. Although these orientation effects represent a signal bias, and are likely to be a nuisance in high-resolution analyses, they may help characterize the vascular influences on candidate fMRI acquisitions and, thereby, may be exploited to improve the neuronal specificity of fMRI.

Mind the gap: Congruence between present and future motivational states shapes prospective decisions

Poor estimation of one's future actions has been associated with the influence of reward over executive control processes during prospection. However, the neural mechanisms underlying this reward-control trade-off remain poorly understood. In the present study, we take advantage of projection bias, underestimating how motivations will change in the future, to examine brain and behavior changes during prospection about future decisions. To manipulate motivation, we altered satiety (hungry vs. satiated) and asked human participants (N = 25) to place bids on snack foods while undergoing fMRI scanning across two sessions. While hungry, participants bid for the right to consume snacks in both a future congruent motivational state (hungry) and a future incongruent motivational state (satiated). In a second session, while satiated, participants placed bids for the right to immediately consume the items. Imagination of a congruent future state was associated with brain activity in regions implicated in prospection. Imagination of an incongruent future state was related to brain activity in areas related to cognitive control. Projection bias, the difference between bids during incongruent prospection (hungry to satiated, session one) and realization (satiated, session two), was negatively related to thalamic and insular engagement. Bias was positively related to engagement of the ventral striatum, a region involved in reward processing. These results suggest that the relative activation between reward and control systems is influenced by the congruence of present and future motivational states, and shapes bias in predictions about future behavior.

Gastric stimulation drives fast BOLD responses of neural origin

Functional magnetic resonance imaging (fMRI) is commonly thought to be too slow to capture any neural dynamics faster than 0.1 Hz. However, recent findings demonstrate the feasibility of detecting fMRI activity at higher frequencies beyond 0.2 Hz. The origin, reliability, and generalizability of fast fMRI responses are still under debate and await confirmation through animal experiments with fMRI and invasive electrophysiology. Here, we acquired single-echo and multi-echo fMRI, as well as local field potentials, from anesthetized rat brains given gastric electrical stimulation modulated at 0.2, 0.4 and 0.8 Hz. Such gastric stimuli could drive widespread fMRI responses at corresponding frequencies from the somatosensory and cingulate cortices. Such fast fMRI responses were linearly dependent on echo times and thus indicative of blood oxygenation level dependent nature (BOLD). Local field potentials recorded during the same gastric stimuli revealed transient and phase-locked broadband neural responses, preceding the fMRI responses by as short as 0.5 s. Taken together, these results suggest that gastric stimulation can drive widespread and rapid fMRI responses of BOLD and neural origin, lending support to the feasibility of using fMRI to detect rapid changes in neural activity up to 0.8 Hz under visceral stimulation.

Receptor-Enriched Analysis of functional connectivity by targets (REACT): A novel, multimodal analytical approach informed by PET to study the pharmacodynamic response of the brain under MDMA

One of the main limitations of pharmacological fMRI is its inability to provide a molecular insight into the main effect of compounds, leaving an open question about the relationship between drug effects and haemodynamic response. The aim of this study is to investigate the acute effects of 3,4-methylenedioxymethamphetamine (MDMA) on functional connectivity (FC) using a novel multimodal method (Receptor-Enriched Analysis of functional Connectivity by Targets - REACT). This approach enriches the resting state (rs-)fMRI analysis with the molecular information about the distribution density of serotonin receptors in the brain, given the serotonergic action of MDMA.

Twenty healthy subjects participated in this double-blind, placebo-controlled, crossover study. A high-resolution in vivo atlas of four serotonin receptors (5-HT_1A, 5-HT_1B, 5-HT_2A, and 5-HT₄) and its transporter (5-HTT) was used as a template in a two-step multivariate regression analysis to estimate the spatial maps reflecting the whole-brain connectivity behaviour related to each target under placebo and MDMA.

Results showed that the networks exhibiting significant changes after MDMA administration are the ones informed by the 5-HTT and 5-HT_1A distribution density maps, which are the main targets of this compound. Changes in the 5-HT_1A-enriched functional maps were also associated with the pharmacokinetic levels of MDMA and MDMA-induced FC changes in the 5-HT_2A-enriched maps correlated with the spiritual experience subscale of the Altered States of Consciousness Questionnaire.

By enriching the rs-fMRI analysis with molecular data of voxel-wise distribution of the serotonin receptors across the brain, we showed that MDMA effects on FC can be understood through the distribution of its main targets. This result supports the ability of this method to characterise the specificity of the functional response of the brain to MDMA binding to serotonergic receptors, paving the way to the definition of a new fingerprint in the characterization of new compounds and potentially to a further understanding to the response to treatment.

•

MDMA connectivity effects understood through the distribution of 5-HT_1A and 5-HTT.

•

Direct link between PK levels of MDMA and 5-HT_1A-enriched functional connectivity maps.

•

Ability to link receptor targets to functional mechanisms underlying behaviour.

•

Mapping pharmacodynamic effects onto the drug's molecular targets.

Differential impact of reward and punishment on functional connectivity after skill learning

Reward and punishment shape behavior, but the mechanisms underlying their effect on skill learning are not well understood. Here, we tested whether the functional connectivity of premotor cortex (PMC), a region known to be critical for learning of sequencing skills, is altered after training when reward or punishment is given during training. Resting-state fMRI was collected in two experiments before and after participants trained on either a serial reaction time task (SRTT; n = 36) or force-tracking task (FTT; n = 36) with reward, punishment, or control feedback. In each experiment, training-related change in PMC functional connectivity was compared across feedback groups. In both tasks, we found that reward and punishment differentially affected PMC functional connectivity. On the SRTT, participants trained with reward showed an increase in functional connectivity between PMC and cerebellum as well as PMC and striatum, while participants trained with punishment showed an increase in functional connectivity between PMC and medial temporal lobe connectivity. After training on the FTT, subjects trained with control and reward showed increases in PMC connectivity with parietal and temporal cortices after training, while subjects trained with punishment showed increased PMC connectivity with ventral striatum. While the results from the two experiments overlapped in some areas, including ventral pallidum, temporal lobe, and cerebellum, these regions showed diverging patterns of results across the two tasks for the different feedback conditions. These findings suggest that reward and punishment strongly influence spontaneous brain activity after training, and that the regions implicated depend on the task learned.

Improving the Assessment of Breath-Holding Induced Cerebral Vascular Reactivity Using a Multiband Multi-echo ASL/BOLD Sequence

Breath holding (BH) is a viable vasodilatory stimulus for calculating functional MRI-derived cerebral vascular reactivity (CVR). The BH technique suffers from reduced repeatability compared with gas inhalation techniques; however, extra equipment is needed to perform gas inhalation techniques, and this equipment is not available at all institutions. This study aimed to determine the sensitivity and repeatability of BH activation and CVR using a multiband multi-echo simultaneous arterial spin labelling/blood oxygenation level dependent (ASL/BOLD) sequence. Whole-brain images were acquired in 14 volunteers. Ten subjects returned for repeat imaging. Each subject performed four cycles of 16 s BH on expiration interleaved with paced breathing. Following standard preprocessing, the echoes were combined using a T2*-weighted approach. BOLD and ASL BH activation was computed, and CVR was then determined as the percent signal change related to the activation. The "M" parameter from the Davis Model was also computed by incorporating the ASL signal. Our results showed higher BH activation strength, volume, and repeatability for the combined multi-echo (MEC) data compared with the single-echo data. MEC CVR also had higher repeatability, sensitivity, specificity, and reliability compared with the single-echo BOLD data. These data support the usefulness of an MBME ASL/BOLD acquisition for BH CVR and M measurements.

Anterior cingulate cortex connectivity is associated with suppression of behaviour in a rat model of chronic pain

A cardinal feature of persistent pain that follows injury is a general suppression of behaviour, in which motivation is inhibited in a way that promotes energy conservation and recuperation. Across species, the anterior cingulate cortex is associated with the motivational aspects of phasic pain, but whether it mediates motivational functions in persistent pain is less clear. Using burrowing behaviour as an marker of non-specific motivated behaviour in rodents, we studied the suppression of burrowing following painful confirmatory factor analysis or control injection into the right knee joint of 30 rats (14 with pain) and examined associated neural connectivity with ultra-high-field resting state functional magnetic resonance imaging. We found that connectivity between anterior cingulate cortex and subcortical structures including hypothalamic/preoptic nuclei and the bed nucleus of the stria terminalis correlated with the reduction in burrowing behaviour observed following the pain manipulation. In summary, the findings implicate anterior cingulate cortex connectivity as a correlate of the motivational aspect of persistent pain in rodents.

Multiband multi-echo simultaneous ASL/BOLD for task-induced functional MRI

Typical simultaneous blood oxygenation-level dependent (BOLD) and arterial spin labeling (ASL) sequences acquire two echoes, one perfusion-sensitive and one BOLD-sensitive. However, for ASL, spatial resolution and brain coverage are limited due to the T1 decay of the labeled blood. This study applies a sequence combining a multiband acquisition with four echoes for simultaneous BOLD and pseudo-continuous ASL (pCASL) echo planar imaging (MBME ASL/BOLD) for block-design task-fMRI. A multiband acceleration of four was employed to increase brain coverage and reduce slice-timing effects on the ASL signal. Multi-echo independent component analysis (MEICA) was implemented to automatically denoise the BOLD signal by regressing non-BOLD components. This technique led to increased temporal signal-to-noise ratio (tSNR) and BOLD sensitivity. The MEICA technique was also modified to denoise the ASL signal by regressing artifact and BOLD signals from the first echo time-series. The MBME ASL/BOLD sequence was applied to a finger-tapping task functional MRI (fMRI) experiment. Signal characteristics and activation were evaluated using single echo BOLD, combined ME BOLD, combined ME BOLD after MEICA denoising, perfusion-weighted (PW), and perfusion-weighted after MEICA denoising time-series. The PW data was extracted using both surround subtraction and high-pass filtering followed by demodulation. In addition, the CBF/BOLD response ratio and CBF/BOLD coupling were analyzed. Results showed that the MEICA denoising procedure significantly improved the BOLD signal, leading to increased BOLD sensitivity, tSNR, and activation statistics compared to conventional single echo BOLD data. At the same time, the denoised PW data showed increased tSNR and activation statistics compared to the non-denoised PW data. CBF/BOLD coupling was also increased using the denoised ASL and BOLD data. Our preliminary data suggest that the MBME ASL/BOLD sequence can be employed to collect whole-brain task-fMRI with improved data quality for both BOLD and PW time series, thus improving the results of block-design task fMRI.

Investigating the spatiotemporal characteristics of the deoxyhemoglobin-related and deoxyhemoglobin-unrelated functional hemodynamic response across cortical layers in awake marmosets

Blood oxygenation level dependent (BOLD) functional magnetic resonance imaging (fMRI) has become a major tool to map neural activity. However, the spatiotemporal characteristics of the BOLD functional hemodynamic response across the cortical layers remain poorly understood. While human fMRI studies suffer from low spatiotemporal resolution, the use of anesthesia in animal models introduces confounding factors. Additionally, inflow contributions to the fMRI signal become non-negligible when short repetition times (TRs) are used. In the present work, we mapped the BOLD fMRI response to somatosensory stimulation in awake marmosets. To address the above technical concerns, we used a dual-echo gradient-recalled echo planar imaging (GR-EPI) sequence to separate the deoxyhemoglobin-related response (absolute T2* differences) from the deoxyhemoglobin-unrelated response (relative S0 changes). We employed a spatial saturation pulse to saturate incoming arterial spins and reduce inflow effects. Functional GR-EPI images were obtained from a single coronal slice with two different echo times (13.5 and 40.5ms) and TR=0.2s. BOLD, T2*, and S0 images were calculated and their functional responses were detected in both hemispheres of primary somatosensory cortex, from which five laminar regions (L1+2, L3, L4, L5, and L6) were derived. The spatiotemporal distribution of the BOLD response across the cortical layers was heterogeneous, with the middle layers having the highest BOLD amplitudes and shortest onset times. ΔT2* also showed a similar trend. However, functional S0 changes were detected only in L1+2, with a fast onset time. Because inflow effects were minimized, the source of S0 functional changes in L1+2 could be attributed to a reduction of cerebrospinal fluid volume fraction due to the functional increase in cerebral blood volume and to unmodeled T2* changes in the extra- and intra-venous compartments. Caution should be exercised when interpreting laminar BOLD fMRI changes in superficial layers as surrogates of underlying neural activity.

A Resilient, Non-neuronal Source of the Spatiotemporal Lag Structure Detected by BOLD Signal-Based Blood Flow Tracking

Recent evidence has suggested that blood oxygenation level-dependent (BOLD) signals convey information about brain circulation via low frequency oscillation of systemic origin (sLFO) that travels through the vascular structure ("lag mapping"). Prompted by its promising application in both physiology and pathology, we examined this signal component using multiple approaches. A total of 30 healthy volunteers were recruited to perform two reproducibility experiments at 3 Tesla using multiband echo planar imaging. The first experiment investigated the effect of denoising and the second was designed to study the effect of subject behavior on lag mapping. The lag map's intersession test-retest reproducibility and image contrast were both diminished by removal of either the neuronal or the non-neuronal (e.g., cardiac, respiratory) components by independent component analysis-based denoising, suggesting that the neurovascular coupling also comprises a part of the BOLD lag structure. The lag maps were, at the same time, robust against local perfusion increases due to visuomotor task and global changes in perfusion induced by breath-holding at the same level as the intrasession reliability. The lag structure was preserved after time-locked averaging to the visuomotor task and breath-holding events, while any preceding signal changes were canceled out for the visuomotor task, consistent with the passive effect of neurovascular coupling in the venous side of the vasculature. These findings support the current assumption that lag mapping primarily reflects vascular structure despite the presence of sLFO perturbation of neuronal or non-neuronal origin and, thus, emphasize the vascular origin of the lag map, encouraging application of BOLD-based blood flow tracking.

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.

Comparing resting state fMRI de-noising approaches using multi- and single-echo acquisitions

Artifact removal in resting state fMRI (rfMRI) data remains a serious challenge, with even subtle head motion undermining reliability and reproducibility. Here we compared some of the most popular single-echo de-noising methods—regression of Motion parameters, White matter and Cerebrospinal fluid signals (MWC method), FMRIB’s ICA-based X-noiseifier (FIX) and ICA-based Automatic Removal Of Motion Artifacts (ICA-AROMA)—with a multi-echo approach (ME-ICA) that exploits the linear dependency of BOLD on the echo time. Data were acquired using a clinical scanner and included 30 young, healthy participants (minimal head motion) and 30 Attention Deficit Hyperactivity Disorder patients (greater head motion). De-noising effectiveness was assessed in terms of data quality after each cleanup procedure, ability to uncouple BOLD signal and motion and preservation of default mode network (DMN) functional connectivity. Most cleaning methods showed a positive impact on data quality. However, based on the investigated metrics, ME-ICA was the most robust. It minimized the impact of motion on FC even for high motion participants and preserved DMN functional connectivity structure. The high-quality results obtained using ME-ICA suggest that using a multi-echo EPI sequence, reliable rfMRI data can be obtained in a clinical setting.

Multiband multi-echo imaging of simultaneous oxygenation and flow timeseries for resting state connectivity

A novel sequence has been introduced that combines multiband imaging with a multi-echo acquisition for simultaneous high spatial resolution pseudo-continuous arterial spin labeling (ASL) and blood-oxygenation-level dependent (BOLD) echo-planar imaging (MBME ASL/BOLD). Resting-state connectivity in healthy adult subjects was assessed using this sequence. Four echoes were acquired with a multiband acceleration of four, in order to increase spatial resolution, shorten repetition time, and reduce slice-timing effects on the ASL signal. In addition, by acquiring four echoes, advanced multi-echo independent component analysis (ME-ICA) denoising could be employed to increase the signal-to-noise ratio (SNR) and BOLD sensitivity. Seed-based and dual-regression approaches were utilized to analyze functional connectivity. Cerebral blood flow (CBF) and BOLD coupling was also evaluated by correlating the perfusion-weighted timeseries with the BOLD timeseries. These metrics were compared between single echo (E2), multi-echo combined (MEC), multi-echo combined and denoised (MECDN), and perfusion-weighted (PW) timeseries. Temporal SNR increased for the MECDN data compared to the MEC and E2 data. Connectivity also increased, in terms of correlation strength and network size, for the MECDN compared to the MEC and E2 datasets. CBF and BOLD coupling was increased in major resting-state networks, and that correlation was strongest for the MECDN datasets. These results indicate our novel MBME ASL/BOLD sequence, which collects simultaneous high-resolution ASL/BOLD data, could be a powerful tool for detecting functional connectivity and dynamic neurovascular coupling during the resting state. The collection of more than two echoes facilitates the use of ME-ICA denoising to greatly improve the quality of resting state functional connectivity MRI.

Brain network reorganization differs in response to stress in rats genetically predisposed to depression and stress-resilient rats

Treatment-resistant depression (TRD) remains a pressing clinical problem. Optimizing treatment requires better definition of the specificity of the involved brain circuits. The rat strain bred for negative cognitive state (NC) represents a genetic animal model of TRD with high face, construct and predictive validity. Vice versa, the positive cognitive state (PC) strain represents a stress-resilient phenotype. Although NC rats show depressive-like behavior, some symptoms such as anhedonia require an external trigger, i.e. a stressful event, which is similar to humans when stressful event induces a depressive episode in genetically predisposed individuals (gene-environment interaction). We aimed to distinguish neurobiological predisposition from the depressogenic pathology at the level of brain-network reorganization. For this purpose, resting-state functional magnetic resonance imaging time series were acquired at 9.4 Tesla scanner in NC (N=11) and PC (N=7) rats before and after stressful event. We used a graph theory analytical approach to calculate the brain-network global and local properties. There was no difference in the global characteristics between the strains. At the local level, the response in the risk strain was characterized with an increased internodal role and reduced local clustering and efficiency of the anterior cingulate cortex (ACC) and prelimbic cortex compared to the stress-resilient strain. We suggest that the increased internodal role of these prefrontal regions could be due to the enhancement of some of their long-range connections, given their connectivity with the amygdala and other default-mode-like network hubs, which could create a bias to attend to negative information characteristic for depression.

Scene-Selectivity and Retinotopy in Medial Parietal Cortex

Functional imaging studies in human reliably identify a trio of scene-selective regions, one on each of the lateral [occipital place area (OPA)], ventral [parahippocampal place area (PPA)], and medial [retrosplenial complex (RSC)] cortical surfaces. Recently, we demonstrated differential retinotopic biases for the contralateral lower and upper visual fields within OPA and PPA, respectively. Here, using functional magnetic resonance imaging, we combine detailed mapping of both population receptive fields (pRF) and category-selectivity, with independently acquired resting-state functional connectivity analyses, to examine scene and retinotopic processing within medial parietal cortex. We identified a medial scene-selective region, which was contained largely within the posterior and ventral bank of the parieto-occipital sulcus (POS). While this region is typically referred to as RSC, the spatial extent of our scene-selective region typically did not extend into retrosplenial cortex, and thus we adopt the term medial place area (MPA) to refer to this visually defined scene-selective region. Intriguingly MPA co-localized with a region identified solely on the basis of retinotopic sensitivity using pRF analyses. We found that MPA demonstrates a significant contralateral visual field bias, coupled with large pRF sizes. Unlike OPA and PPA, MPA did not show a consistent bias to a single visual quadrant. MPA also co-localized with a region identified by strong differential functional connectivity with PPA and the human face-selective fusiform face area (FFA), commensurate with its functional selectivity. Functional connectivity with OPA was much weaker than with PPA, and similar to that with face-selective occipital face area (OFA), suggesting a closer link with ventral than lateral cortex. Consistent with prior research, we also observed differential functional connectivity in medial parietal cortex for anterior over posterior PPA, as well as a region on the lateral surface, the caudal inferior parietal lobule (cIPL). However, the differential connectivity in medial parietal cortex was found principally anterior of MPA. We suggest that there is posterior–anterior gradient within medial parietal cortex, with posterior regions in the POS showing retinotopically based scene-selectivity and more anterior regions showing connectivity that may be more reflective of abstract, navigationally pertinent and possibly mnemonic representations.

Connectome analysis for pre-operative brain mapping in neurosurgery

Object: Brain mapping has entered a new era focusing on complex network connectivity. Central to this is the search for the connectome or the brains ‘wiring diagram’. Graph theory analysis of the connectome allows understanding of the importance of regions to network function, and the consequences of their impairment or excision. Our goal was to apply connectome analysis in patients with brain tumours to characterise overall network topology and individual patterns of connectivity alterations.

Methods: Resting-state functional MRI data were acquired using multi-echo, echo planar imaging pre-operatively from five participants each with a right temporal–parietal–occipital glioblastoma. Complex networks analysis was initiated by parcellating the brain into anatomically regions amongst which connections were identified by retaining the most significant correlations between the respective wavelet decomposed time-series.

Results: Key characteristics of complex networks described in healthy controls were preserved in these patients, including ubiquitous small world organization. An exponentially truncated power law fit to the degree distribution predicted findings of general network robustness to injury but with a core of hubs exhibiting disproportionate vulnerability. Tumours produced a consistent reduction in local and long-range connectivity with distinct patterns of connection loss depending on lesion location.

Conclusions: Connectome analysis is a feasible and novel approach to brain mapping in individual patients with brain tumours. Applications to pre-surgical planning include identifying regions critical to network function that should be preserved and visualising connections at risk from tumour resection. In the future one could use such data to model functional plasticity and recovery of cognitive deficits.

Improving effect size estimation and statistical power with multi-echo fMRI and its impact on understanding the neural systems supporting mentalizing

Functional magnetic resonance imaging (fMRI) research is routinely criticized for being statistically underpowered due to characteristically small sample sizes and much larger sample sizes are being increasingly recommended. Additionally, various sources of artifact inherent in fMRI data can have detrimental impact on effect size estimates and statistical power. Here we show how specific removal of non-BOLD artifacts can improve effect size estimation and statistical power in task-fMRI contexts, with particular application to the social-cognitive domain of mentalizing/theory of mind. Non-BOLD variability identification and removal is achieved in a biophysical and statistically principled manner by combining multi-echo fMRI acquisition and independent components analysis (ME-ICA). Without smoothing, group-level effect size estimates on two different mentalizing tasks were enhanced by ME-ICA at a median rate of 24% in regions canonically associated with mentalizing, while much more substantial boosts (40-149%) were observed in non-canonical cerebellar areas. Effect size boosting occurs via reduction of non-BOLD noise at the subject-level and consequent reductions in between-subject variance at the group-level. Smoothing can attenuate ME-ICA-related effect size improvements in certain circumstances. Power simulations demonstrate that ME-ICA-related effect size enhancements enable much higher-powered studies at traditional sample sizes. Cerebellar effects observed after applying ME-ICA may be unobservable with conventional imaging at traditional sample sizes. Thus, ME-ICA allows for principled design-agnostic non-BOLD artifact removal that can substantially improve effect size estimates and statistical power in task-fMRI contexts. ME-ICA could mitigate some issues regarding statistical power in fMRI studies and enable novel discovery of aspects of brain organization that are currently under-appreciated and not well understood.

Subtle in-scanner motion biases automated measurement of brain anatomy from in vivo MRI

While the potential for small amounts of motion in functional magnetic resonance imaging (fMRI) scans to bias the results of functional neuroimaging studies is well appreciated, the impact of in-scanner motion on morphological analysis of structural MRI is relatively under-studied. Even among "good quality" structural scans, there may be systematic effects of motion on measures of brain morphometry. In the present study, the subjects' tendency to move during fMRI scans, acquired in the same scanning sessions as their structural scans, yielded a reliable, continuous estimate of in-scanner motion. Using this approach within a sample of 127 children, adolescents, and young adults, significant relationships were found between this measure and estimates of cortical gray matter volume and mean curvature, as well as trend-level relationships with cortical thickness. Specifically, cortical volume and thickness decreased with greater motion, and mean curvature increased. These effects of subtle motion were anatomically heterogeneous, were present across different automated imaging pipelines, showed convergent validity with effects of frank motion assessed in a separate sample of 274 scans, and could be demonstrated in both pediatric and adult populations. Thus, using different motion assays in two large non-overlapping sets of structural MRI scans, convergent evidence showed that in-scanner motion-even at levels which do not manifest in visible motion artifact-can lead to systematic and regionally specific biases in anatomical estimation. These findings have special relevance to structural neuroimaging in developmental and clinical datasets, and inform ongoing efforts to optimize neuroanatomical analysis of existing and future structural MRI datasets in non-sedated humans. Hum Brain Mapp 37:2385-2397, 2016. © 2016 Wiley Periodicals, Inc.

Robust resting state fMRI processing for studies on typical brain development based on multi-echo EPI acquisition

Several methodological challenges affect the study of typical brain development based on resting state blood oxygenation level dependent (BOLD) functional MRI (fMRI). One such challenge is mitigating artifacts such as those from head motion, known to be more substantial in younger subjects than older subjects. Other challenges include controlling for potential age-dependence in cerebrospinal fluid (CSF) volume affecting anatomical-functional coregistration; in vascular density affecting BOLD contrast-to-noise; and in CSF pulsation creating time series artifacts. Historically, these confounds have been approached through incorporating artifact-specific temporal and/or spatial filtering into preprocessing pipelines. However, such paths often come with new confounds or limitations. In this study we take the approach of a bottom-up revision of fMRI methodology based on acquisition of multi-echo fMRI and comprehensive utilization of the information in the TE-domain to enhance several aspects of fMRI analysis in the context of a developmental study. We show in a cohort of 25 healthy subjects, aged 9 to 43 years, that the analysis of multi-echo fMRI data eliminates a number of arbitrary processing steps such as bandpass filtering and spatial smoothing, while enabling procedures such as [Formula: see text] mapping, BOLD contrast normalization and signal dropout recovery, precise anatomical-functional coregistration based on [Formula: see text] measurements, automatic denoising through removing subject motion, scanner-related signal drifts and physiology, as well as statistical inference for seed-based connectivity. These enhancements are of both theoretical significance and practical benefit in the study of typical brain development.

Large-scale functional connectivity networks in the rodent brain

Resting-state functional Magnetic Resonance Imaging (rsfMRI) of the human brain has revealed multiple large-scale neural networks within a hierarchical and complex structure of coordinated functional activity. These distributed neuroanatomical systems provide a sensitive window on brain function and its disruption in a variety of neuropathological conditions. The study of macroscale intrinsic connectivity networks in preclinical species, where genetic and environmental conditions can be controlled and manipulated with high specificity, offers the opportunity to elucidate the biological determinants of these alterations. While rsfMRI methods are now widely used in human connectivity research, these approaches have only relatively recently been back-translated into laboratory animals. Here we review recent progress in the study of functional connectivity in rodent species, emphasising the ability of this approach to resolve large-scale brain networks that recapitulate neuroanatomical features of known functional systems in the human brain. These include, but are not limited to, a distributed set of regions identified in rats and mice that may represent a putative evolutionary precursor of the human default mode network (DMN). The impact and control of potential experimental and methodological confounds are also critically discussed. Finally, we highlight the enormous potential and some initial application of connectivity mapping in transgenic models as a tool to investigate the neuropathological underpinnings of the large-scale connectional alterations associated with human neuropsychiatric and neurological conditions. We conclude by discussing the translational potential of these methods in basic and applied neuroscience.

Enhanced identification of BOLD-like components with multi-echo simultaneous multi-slice (MESMS) fMRI and multi-echo ICA

The recent introduction of simultaneous multi-slice (SMS) acquisitions has enabled the acquisition of blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) data with significantly higher temporal sampling rates. In a parallel development, the use of multi-echo fMRI acquisitions in conjunction with a multi-echo independent component analysis (ME-ICA) approach has been introduced as a means to automatically distinguish functionally-related BOLD signal components from signal artifacts, with significant gains in sensitivity, statistical power, and specificity. In this work, we examine the gains that can be achieved with a combined approach in which data obtained with a multi-echo simultaneous multi-slice (MESMS) acquisition are analyzed with ME-ICA. We find that ME-ICA identifies significantly more BOLD-like components in the MESMS data as compared to data acquired with a conventional multi-echo single-slice acquisition. We demonstrate that the improved performance of MESMS derives from both an increase in the number of temporal samples and the enhanced ability to filter out high-frequency artifacts.