On the origin of respiratory artifacts in BOLD-EPI of the human brain

Robust data-driven segmentation of pulsatile cerebral vessels using functional magnetic resonance imaging

Functional magnetic resonance imaging (fMRI) captures rich physiological and neuronal information, offering insight into neurofluid dynamics, vascular health and waste clearance. Accurate cerebral vessel segmentation could greatly facilitate fluid dynamics research in fMRI. However, existing vessel identification methods, such as magnetic resonance angiography or deep-learning-based segmentation on structural MRI, cannot reliably locate cerebral vessels in fMRI space due to misregistration from inherent fMRI distortions. To address this challenge, we developed a data-driven, automatic segmentation of cerebral vessels directly within fMRI space. This approach identified large cerebral arteries and the superior sagittal sinus (SSS) by leveraging these vessels' distinct pulsatile signal patterns during the cardiac cycle. The method was validated in a local dataset by comparing it to ground truth cerebral artery and SSS segmentations. Using the Human Connectome Project (HCP) ageing dataset, the method's reproducibility was tested on 422 participants aged 36-90, each with four repeated fMRI scans. The method demonstrated high reproducibility, with an intraclass correlation coefficient > 0.7 in both cerebral artery and SSS segmentation volumes. This study demonstrates that large cerebral arteries and SSS can be reproducibly and automatically segmented in fMRI datasets, facilitating reliable fluid dynamics investigation in these regions.

Robust data-driven segmentation of pulsatile cerebral vessels using functional magnetic resonance imaging

Functional magnetic resonance imaging (fMRI) captures rich physiological and neuronal information that can offer insights into neurofluid dynamics, vascular health, and waste clearance function. The availability of cerebral vessel segmentation could facilitate fluid dynamics research in fMRI. However, without magnetic resonance angiography scans, cerebral vessel segmentation is challenging and time-consuming. This study leverages cardiac-induced pulsatile fMRI signal to develop a data-driven, automatic segmentation of large cerebral arteries and the superior sagittal sinus (SSS). The method was validated in a local dataset by comparing it to ground truth cerebral artery and SSS segmentations. Using the Human Connectome Project (HCP) aging dataset, the method's reproducibility was tested on 422 participants aged 36 to 100 years, each with four repeated fMRI scans. The method demonstrated high reproducibility, with an intraclass correlation coefficient > 0.7 in both cerebral artery and SSS segmentation volumes. This study demonstrates that the large cerebral arteries and SSS can be reproducibly and automatically segmented in fMRI datasets, facilitating the investigation of fluid dynamics in these regions.

"Listen to your heart": A novel interoceptive strategy for real-time fMRI neurofeedback training of anterior insula activity

Real-time fMRI (rt-fMRI) neurofeedback (NF) training is a novel non-invasive technique for volitional brain modulation. Given the important role of the anterior insula (AI) in human cognitive and affective processes, it has become one of the most investigated regions in rt-fMRI studies. Most rt-fMRI insula studies employed emotional recall/imagery as the regulation strategy, which may be less effective for psychiatric disorders characterized by altered emotional processing. The present study thus aimed to examine the feasibility of a novel interoceptive strategy based on heartbeat detection in rt-fMRI guided AI regulation and its associated behavioral changes using a randomized double-blind, sham feedback-controlled between-subject design. 66 participants were recruited and randomly assigned to receive either NF from the left AI (LAI) or sham feedback from a control region while using the interoceptive strategy. N = 57 participants were included in the final data analyses. Empathic and interoceptive pre-post training changes were collected as behavioral measures of NF training effects. Results showed that participants in the NF group exhibited stronger LAI activity than the control group with LAI activity being positively correlated with interoceptive accuracy following NF training, although there were no significant increases of LAI activity over training sessions. Importantly, ability of LAI regulation could be maintained in a transfer session without feedback. Successful LAI regulation was associated with strengthened functional connectivity of the LAI with cognitive control, memory and learning, and salience/interoceptive networks. The present study demonstrated for the first time the efficacy of a novel regulation strategy based on interoceptive processing in up-regulating LAI activity. Our findings also provide proof of concept for the translational potential of this strategy in rt-fMRI AI regulation of psychiatric disorders characterized by altered emotional processing.

Differences in Olivo-Cerebellar Circuit and Cerebellar Network Connectivity in Essential Tremor: a Resting State fMRI Study

The olivo-cerebellar circuit is thought to play a crucial role in the pathophysiology of essential tremor (ET). Whether olivo-cerebellar circuit dysfunction is also present at rest, in the absence of clinical tremor and linked voluntary movement, remains unclear. Assessing this network in detail with fMRI is challenging, considering the brainstem is close to major arteries and pulsatile cerebrospinal fluid-filled spaces obscuring signals of interest. Here, we used methods tailored to the analysis of infratentorial structures. We hypothesize that the olivo-cerebellar circuit shows altered intra-network connectivity at rest and decreased functional coupling with other parts of the motor network in ET. In 17 ET patients and 19 healthy controls, we investigated using resting state fMRI intracerebellar functional and effective connectivity on a dedicated cerebellar atlas. With independent component analysis, we investigated data-driven cerebellar motor network activations during rest. Finally, whole-brain connectivity of cerebellar motor structures was investigated using identified components. In ET, olivo-cerebellar pathways show decreased functional connectivity compared with healthy controls. Effective connectivity analysis showed an increased inhibitory influence of the dentate nucleus towards the inferior olive. Cerebellar independent component analyses showed motor resting state networks are less strongly connected to the cerebral cortex compared to controls. Our results indicate the olivo-cerebellar circuit to be affected at rest. Also, the cerebellum is "disconnected" from the rest of the motor network. Aberrant activity, generated within the olivo-cerebellar circuit could, during action, spread towards other parts of the motor circuit and potentially underlie the characteristic tremor of this patient group.

Quantifying numerical and spatial reliability of hippocampal and amygdala subdivisions in FreeSurfer

On-going, large-scale neuroimaging initiatives can aid in uncovering neurobiological causes and correlates of poor mental health, disease pathology, and many other important conditions. As projects grow in scale with hundreds, even thousands, of individual participants and scans collected, quantification of brain structures by automated algorithms is becoming the only truly tractable approach. Here, we assessed the spatial and numerical reliability for newly deployed automated segmentation of hippocampal subfields and amygdala nuclei in FreeSurfer 7. In a sample of participants with repeated structural imaging scans (N = 928), we found numerical reliability (as assessed by intraclass correlations, ICCs) was reasonable. Approximately 95% of hippocampal subfields had "excellent" numerical reliability (ICCs ≥ 0.90), while only 67% of amygdala subnuclei met this same threshold. In terms of spatial reliability, 58% of hippocampal subfields and 44% of amygdala subnuclei had Dice coefficients ≥ 0.70. Notably, multiple regions had poor numerical and/or spatial reliability. We also examined correlations between spatial reliability and person-level factors (e.g., participant age; T1 image quality). Both sex and image scan quality were related to variations in spatial reliability metrics. Examined collectively, our work suggests caution should be exercised for a few hippocampal subfields and amygdala nuclei with more variable reliability.

Unsupervised physiological noise correction of functional magnetic resonance imaging data using phase and magnitude information (PREPAIR)

Of the sources of noise affecting blood oxygen level-dependent functional magnetic resonance imaging (fMRI), respiration and cardiac fluctuations are responsible for the largest part of the variance, particularly at high and ultrahigh field. Existing approaches to removing physiological noise either use external recordings, which can be unwieldy and unreliable, or attempt to identify physiological noise from the magnitude fMRI data. Data-driven approaches are limited by sensitivity, temporal aliasing, and the need for user interaction. In the light of the sensitivity of the phase of the MR signal to local changes in the field stemming from physiological processes, we have developed an unsupervised physiological noise correction method using the information carried in the phase and the magnitude of echo-planar imaging data. Our technique, Physiological Regressor Estimation from Phase and mAgnItude, sub-tR (PREPAIR) derives time series signals sampled at the slice TR from both phase and magnitude images. It allows physiological noise to be captured without aliasing, and efficiently removes other sources of signal fluctuations not related to physiology, prior to regressor estimation. We demonstrate that the physiological signal time courses identified with PREPAIR agree well with those from external devices and retrieve challenging cardiac dynamics. The removal of physiological noise was as effective as that achieved with the most used approach based on external recordings, RETROICOR. In comparison with widely used recording-free physiological noise correction tools-PESTICA and FIX, both performed in unsupervised mode-PREPAIR removed significantly more respiratory and cardiac noise than PESTICA, and achieved a larger increase in temporal signal-to-noise-ratio at both 3 and 7 T.

Neural underpinning of a respiration-associated resting-state fMRI network

Respiration can induce motion and CO2 fluctuation during resting-state fMRI (rsfMRI) scans, which will lead to non-neural artifacts in the rsfMRI signal. In the meantime, as a crucial physiologic process, respiration can directly drive neural activity change in the brain, and may thereby modulate the rsfMRI signal. Nonetheless, this potential neural component in the respiration-fMRI relationship is largely unexplored. To elucidate this issue, here we simultaneously recorded the electrophysiology, rsfMRI, and respiration signals in rats. Our data show that respiration is indeed associated with neural activity changes, evidenced by a phase-locking relationship between slow respiration variations and the gamma-band power of the electrophysiological signal recorded in the anterior cingulate cortex. Intriguingly, slow respiration variations are also linked to a characteristic rsfMRI network, which is mediated by gamma-band neural activity. In addition, this respiration-related brain network disappears when brain-wide neural activity is silenced at an isoelectrical state, while the respiration is maintained, further confirming the necessary role of neural activity in this network. Taken together, this study identifies a respiration-related brain network underpinned by neural activity, which represents a novel component in the respiration-rsfMRI relationship that is distinct from respiration-related rsfMRI artifacts. It opens a new avenue for investigating the interactions between respiration, neural activity, and resting-state brain networks in both healthy and diseased conditions.

fMRI studies evaluating central respiratory control in humans

A plethora of neural centers in the central nervous system control the fundamental respiratory pattern. This control is ensured by neurons that act as pacemakers, modulating activity through chemical control driven by changes in the O2/CO2 balance. Most of the respiratory neural centers are located in the brainstem, but difficult to localize on magnetic resonance imaging (MRI) due to their small size, lack of visually-detectable borders with neighboring areas, and significant physiological noise hampering detection of its activity with functional MRI (fMRI). Yet, several approaches make it possible to study the normal response to different abnormal stimuli or conditions such as CO2 inhalation, induced hypercapnia, volitional apnea, induced hypoxia etc. This review provides a comprehensive overview of the majority of available studies on central respiratory control in humans.

WHOCARES: WHOle-brain CArdiac signal REgression from highly accelerated simultaneous multi-Slice fMRI acquisitions

Objective. To spatio-temporally resolve cardiac signals in functional magnetic resonance imaging (fMRI) time-series of the human brain using neither external physiological measurements nor ad hoc modelling assumptions.Approach. Cardiac pulsation is a physiological confound of fMRI time-series that introduces spurious signal fluctuations in proximity to blood vessels. fMRI alone is not sufficiently fast to resolve cardiac pulsation. Depending on the ratio between the instantaneous heart-rate and the acquisition sampling frequency (1/TR, with TR being the repetition time), the cardiac signal may alias into the frequency band of neural activation so that its removal through spectral filtering techniques is generally not possible. In this paper, we show that it is feasible to temporally and spatially resolve cardiac signals throughout the brain even when cardiac aliasing occurs by combining fMRI hyper-sampling with simultaneous multislice (SMS) imaging. The technique, which we name WHOle-brain CArdiac signal REgression from highly accelerated simultaneous multi-Slice fMRI acquisitions (WHOCARES), was developed on 695 healthy subjects selected from the Human Connectome Project and its performance validated against the RETROICOR, HAPPY and the pulse oxymeter signal regression methods.Main results.WHOCARES is capable of retrieving voxel-wise cardiac signal regressors. This is achieved without employing external physiological recordings nor through ad hoc modelling assumptions. The performance of WHOCARES was, on average, superior to RETROICOR, HAPPY and the pulse oxymeter regression methods.Significance.WHOCARES holds basis for the reliable mapping of cardiac activity in fMRI time-series. WHOCARES can be employed for the retrospective removal of cardiac noise in publicly available fMRI datasets where physiological recordings are not available. WHOCARES is freely available athttps://github.com/gferrazzi/WHOCARES.

Neural synchronization predicts marital satisfaction

Marital attachment plays an important role in maintaining intimate personal relationships and sustaining psychological well-being. Mate-selection theories suggest that people are more likely to marry someone with a similar personality and social status, yet evidence for the association between personality-based couple similarity measures and marital satisfaction has been inconsistent. A more direct and useful approach for understanding fundamental processes underlying marital satisfaction is to probe similarity of dynamic brain responses to maritally and socially relevant communicative cues, which may better reflect how married couples process information in real time and make sense of their mates and themselves. Here, we investigate shared neural representations based on intersubject synchronization (ISS) of brain responses during free viewing of marital life-related, and nonmarital, object-related movies. Compared to randomly selected pairs of couples, married couples showed significantly higher levels of ISS during viewing of marital movies and ISS between married couples predicted higher levels of marital satisfaction. ISS in the default mode network emerged as a strong predictor of marital satisfaction and canonical correlation analysis revealed a specific relation between ISS in this network and shared communication and egalitarian components of martial satisfaction. Our findings demonstrate that brain similarities that reflect real-time mental responses to subjective perceptions, thoughts, and feelings about interpersonal and social interactions are strong predictors of marital satisfaction, reflecting shared values and beliefs. Our study advances foundational knowledge of the neurobiological basis of human pair bonding.

Parsing neural circuits of fear learning and extinction across basic and clinical neuroscience: Towards better translation

Over the past decades, studies of fear learning and extinction have advanced our understanding of the neurobiology of threat and safety learning. Animal studies can provide mechanistic/causal insights into human brain regions and their functional connectivity involved in fear learning and extinction. Findings in humans, conversely, may further enrich our understanding of neural circuits in animals by providing macroscopic insights at the level of brain-wide networks. Nevertheless, there is still much room for improvement in translation between basic and clinical research on fear learning and extinction. Through the lens of neural circuits, in this article, we aim to review the current knowledge of fear learning and extinction in both animals and humans, and to propose strategies to fill in the current knowledge gap for the purpose of enhancing clinical benefits.

Cardiovascular Pulsatility Increases in Visual Cortex Before Blood Oxygen Level Dependent Response During Stimulus

The physiological pulsations that drive tissue fluid homeostasis are not well characterized during brain activation. Therefore, we used fast magnetic resonance encephalography (MREG) fMRI to measure full band (0–5 Hz) blood oxygen level-dependent (BOLD_FB) signals during a dynamic visual task in 23 subjects. This revealed brain activity in the very low frequency (BOLD_VLF) as well as in cardiac and respiratory bands. The cardiovascular hemodynamic envelope (CHe) signal correlated significantly with the visual BOLD_VLF response, considered as an independent signal source in the V1-V2 visual cortices. The CHe preceded the canonical BOLD_VLF response by an average of 1.3 (± 2.2) s. Physiologically, the observed CHe signal could mark increased regional cardiovascular pulsatility following vasodilation.

Metabolism modulates network synchrony in the aging brain

Brain aging is associated with hypometabolism and global changes in functional connectivity. Using functional MRI (fMRI), we show that network synchrony, a collective property of brain activity, decreases with age. Applying quantitative methods from statistical physics, we provide a generative (Ising) model for these changes as a function of the average communication strength between brain regions. We find that older brains are closer to a critical point of this communication strength, in which even small changes in metabolism lead to abrupt changes in network synchrony. Finally, by experimentally modulating metabolic activity in younger adults, we show how metabolism alone-independent of other changes associated with aging-can provide a plausible candidate mechanism for marked reorganization of brain network topology.

Spectral analysis of physiological brain pulsations affecting the BOLD signal

Physiological pulsations have been shown to affect the global blood oxygen level dependent (BOLD) signal in human brain. While these pulsations have previously been regarded as noise, recent studies show their potential as biomarkers of brain pathology. We used the extended 5 Hz spectral range of magnetic resonance encephalography (MREG) data to investigate spatial and frequency distributions of physiological BOLD signal sources. Amplitude spectra of the global image signals revealed cardiorespiratory envelope modulation (CREM) peaks, in addition to the previously known very low frequency (VLF) and cardiorespiratory pulsations. We then proceeded to extend the amplitude of low frequency fluctuations (ALFF) method to each of these pulsations. The respiratory pulsations were spatially dominating over most brain structures. The VLF pulsations overcame the respiratory pulsations in frontal and parietal gray matter, whereas cardiac and CREM pulsations had this effect in central cerebrospinal fluid (CSF) spaces and major blood vessels. A quasi‐periodic pattern (QPP) analysis showed that the CREM pulsations propagated as waves, with a spatiotemporal pattern differing from that of respiratory pulsations, indicating them to be distinct intracranial physiological phenomenon. In conclusion, the respiration has a dominant effect on the global BOLD signal and directly modulates cardiovascular brain pulsations.

Estimates of locus coeruleus function with functional magnetic resonance imaging are influenced by localization approaches and the use of multi-echo data

The locus coeruleus (LC) plays a central role in regulating human cognition, arousal, and autonomic states. Efforts to characterize the LC’s function in humans using functional magnetic resonance imaging have been hampered by its small size and location near a large source of noise, the fourth ventricle. We tested whether the ability to characterize LC function is improved by employing neuromelanin-T1 weighted images (nmT1) for LC localization and multi-echo functional magnetic resonance imaging (ME-fMRI) for estimating intrinsic functional connectivity (iFC). Analyses indicated that, relative to a probabilistic atlas, utilizing nmT1 images to individually localize the LC increases the specificity of seed time series and clusters in the iFC maps. When combined with independent components analysis (ME-ICA), ME-fMRI data provided significant improvements in the temporal signal to noise ratio and DVARS relative to denoised single echo data (1E-fMRI). The effects of acquiring nmT1 images and ME-fMRI data did not appear to only reflect increases in power: iFC maps for each approach overlapped only moderately. This is consistent with findings that ME-fMRI offers substantial advantages over 1E-fMRI acquisition and denoising. It also suggests that individually identifying LC with nmT1 scans is likely to reduce the influence of other nearby brainstem regions on estimates of LC function.

Task-related activity in human visual cortex

The brain exhibits widespread endogenous responses in the absence of visual stimuli, even at the earliest stages of visual cortical processing. Such responses have been studied in monkeys using optical imaging with a limited field of view over visual cortex. Here, we used functional MRI (fMRI) in human participants to study the link between arousal and endogenous responses in visual cortex. The response that we observed was tightly entrained to task timing, was spatially extensive, and was independent of visual stimulation. We found that this response follows dynamics similar to that of pupil size and heart rate, suggesting that task-related activity is related to arousal. Finally, we found that higher reward increased response amplitude while decreasing its trial-to-trial variability (i.e., the noise). Computational simulations suggest that increased temporal precision underlies both of these observations. Our findings are consistent with optical imaging studies in monkeys and support the notion that arousal increases precision of neural activity.

The brain exhibits widespread endogenous responses in the absence of visual stimuli, even at the earliest stages of visual cortical processing. This fMRI study characterizes a widespread hemodynamic response in early visual cortex that is not related to visual input but instead reflects a participant’s engagement in a task, is modulated by expected monetary reward, and may reflect neural quenching.

Reconstruction of respiratory variation signals from fMRI data

Functional MRI signals can be heavily influenced by systemic physiological processes in addition to local neural activity. For example, widespread hemodynamic fluctuations across the brain have been found to correlate with natural, low-frequency variations in the depth and rate of breathing over time. Acquiring peripheral measures of respiration during fMRI scanning not only allows for modeling such effects in fMRI analysis, but also provides valuable information for interrogating brain-body physiology. However, physiological recordings are frequently unavailable or have insufficient quality. Here, we propose a computational technique for reconstructing continuous low-frequency respiration volume (RV) fluctuations from fMRI data alone. We evaluate the performance of this approach across different fMRI preprocessing strategies. Further, we demonstrate that the predicted RV signals can account for similar patterns of temporal variation in resting-state fMRI data compared to measured RV fluctuations. These findings indicate that fluctuations in respiration volume can be extracted from fMRI alone, in the common scenario of missing or corrupted respiration recordings. The results have implications for enriching a large volume of existing fMRI datasets through retrospective addition of respiratory variations information.

Characterization of Cerebellar Atrophy and Resting State Functional Connectivity Patterns in Sporadic Adult-Onset Ataxia of Unknown Etiology (SAOA)

Sporadic adult-onset ataxia of unknown etiology (SAOA) is a non-genetic neurodegenerative disorder of the cerebellum of unknown cause which manifests with progressive ataxia without severe autonomic failure. Although SAOA is associated with cerebellar degeneration, little is known about the specific cerebellar atrophy pattern in SAOA. Thirty-seven SAOA patients and 49 healthy controls (HCs) were included at two centers. We investigated the structural and functional characteristics of SAOA brains using voxel-based morphometry (VBM) and resting-state functional imaging (rs-fMRI). In order to examine the functional consequence of structural cerebellar alterations, the amplitude of low-frequency fluctuation (ALFF) and degree centrality (DC) were analyzed, and then assessed their relation with disease severity, disease duration, and age of onset within these regions. Group differences were investigated using two-sample t tests, controlling for age, gender, site, and the total intracranial volume. The VBM analysis revealed a significant, mostly bilateral reduction of local gray matter (GM) volume in lobules I-V, V, VI, IX, X, and vermis VIII a/b in SAOA patients, compared with HCs. The GM volume loss in these regions was significantly associated with disease severity, disease duration, and age of onset. The disease-related atrophy regions did not show any functional alternations compared with HCs but were functionally characterized by high ALFF and poor DC compared with intact cerebellar regions. Our data revealed volume reduction in SAOA in cerebellar regions that are known to be involved in motor and somatosensory processing, corresponding with the clinical phenotype of SAOA. Our data suggest that the atrophy occurs in those cerebellar regions which are characterized by high ALFF and poor DC. Further studies have to show if these findings are specific for SAOA, and if they can be used to predict disease progression.

Functional neuroimaging in the acute phase of Takotsubo syndrome: volumetric and functional changes of the right insular cortex

Background

A brain–heart interaction has been proposed in Takotsubo syndrome (TTS). Structural changes in the limbic system and hypoconnectivity between certain brain areas in the chronic phase of the disease have been reported, but little is known concerning functional neuroimaging in the acute phase. We hypothesized anatomical and functional changes in the central nervous system and investigated whole-brain volumetric and functional connectivity alterations in the acute phase TTS patients compared to controls.

Methods

Anatomical and resting-state functional magnetic resonance imaging were performed in postmenopausal females: thirteen in the acute TTS phase and thirteen healthy controls without evidence of coronary artery disease. Voxel-based morphometry and graph theoretical analysis were applied to identify anatomical and functional differences between patients and controls.

Results

Significantly lower gray matter volumes were found in TTS patients in the right middle frontal gyrus (p = 0.004) and right subcallosal cortex (p = 0.009) compared to healthy controls. When lower threshold was applied, volumetric changes were noted in the right insular cortex (p = 0.0113), the right paracingulate cortex (p = 0.012), left amygdala (p = 0.018), left central opercular cortex (p = 0.017), right (p = 0.013) and left thalamus (p = 0.017), and left cerebral cortex (p = 0.017). Graph analysis revealed significantly (p < 0.01) lower functional connectivity in TTS patients compared to healthy controls, particularly in the connections originating from the right insular cortex, temporal lobes, and precuneus.

Conclusion

In the acute phase of TTS volumetric changes in frontal regions and the central autonomic network (i.e. insula, anterior cingulate cortex, and amygdala) were noted. In particular, the right insula, associated with sympathetic autonomic tone, had both volumetric and functional changes.

Graphic abstract

Multimodal Parcellations and Extensive Behavioral Profiling Tackling the Hippocampus Gradient

The hippocampus displays a complex organization and function that is perturbed in many neuropathologies. Histological work revealed a complex arrangement of subfields along the medial-lateral and the ventral-dorsal dimension, which contrasts with the anterior-posterior functional differentiation. The variety of maps has raised the need for an integrative multimodal view. We applied connectivity-based parcellation to 1) intrinsic connectivity 2) task-based connectivity, and 3) structural covariance, as complementary windows into structural and functional differentiation of the hippocampus. Strikingly, while functional properties (i.e., intrinsic and task-based) revealed similar partitions dominated by an anterior-posterior organization, structural covariance exhibited a hybrid pattern reflecting both functional and cytoarchitectonic subdivision. Capitalizing on the consistency of functional parcellations, we defined robust functional maps at different levels of partitions, which are openly available for the scientific community. Our functional maps demonstrated a head-body and tail partition, subdivided along the anterior-posterior and medial-lateral axis. Behavioral profiling of these fine partitions based on activation data indicated an emotion-cognition gradient along the anterior-posterior axis and additionally suggested a self-world-centric gradient supporting the role of the hippocampus in the construction of abstract representations for spatial navigation and episodic memory.

Modulating functional connectivity between medial frontopolar cortex and amygdala by inhibitory and excitatory transcranial magnetic stimulation

The prefrontal-limbic network in the human brain plays a major role in social cognition, especially cognitive control of emotion. The medial frontopolar cortex (mFP; Brodmann Area 10) and the amygdala are part of this network and display correlated neuronal activity in time, as measured by functional magnetic resonance imaging (fMRI). This functional connectivity is dynamic, sensitive to training, and affected in mental disorders. However, the effects of neurostimulation on functional connectivity within this network have not yet been systematically investigated. Here, we investigate the effects of both low- and high-frequency repetitive transcranial magnetic stimulation (rTMS) to the right mFP on functional connectivity between mFP and amygdala, as measured with resting state fMRI (rsfMRI). Three groups of healthy participants received either low-frequency rTMS (1 Hz; N = 18), sham TMS (1 Hz, subthreshold; N = 18) or high-frequency rTMS (20 Hz; N = 19). rsfMRI was acquired before and after (separate days). We hypothesized a modulation of functional connectivity in opposite directions compared to sham TMS through adjustment of the stimulation frequency. Groups differed in functional connectivity between mFP and amygdala after stimulation compared to before stimulation (low-frequency: decrease, high-frequency: increase). Motion or induced changes in neuronal activity were excluded as confounders. Results show that rTMS is effective for increasing and decreasing functional coherence between prefrontal and limbic regions. This finding is relevant for social and affective neuroscience as well as novel treatment approaches in psychiatry.

On the Use of Linear-Modelling-based Algorithms for Physiological Noise Correction in fMRI Studies of the Central Breathing Control

A full characterization of the physiological behavior of human central chemoreceptors through fMRI is crucial to understand the pathophysiology of central abnormal breathing patterns. In this scenario, physiological noise and activity of interest may be naturally correlated. Here, we examined the adequacy of linear-modelling-based retrospective physiological noise correction for studies of the central breathing control. We focused on the relationship between a nonlinear model of BOLD response, hypothesized to describe neuronal specific activity, and noise modelled by correction algorithms. Analyses were performed on fMRI acquisitions from healthy subjects during a breath hold task. A general linear model including static nonlinearities in the response to end-tidal CO₂ was applied to data preprocessed both with and without physiological noise correction. Relations between physiological noise and PETCO₂ were explored both with linear and nonlinear measures. Lastly, parametric maps of noise spatial distribution were extracted. Our results evidenced that correction algorithms based on linear modelling remove components that are both linearly and nonlinearly related to end-tidal CO₂, whereas uncorrected data showed spurious activations in regions outside gray matter. Thus, despite a correction step is fundamental, these algorithms are shown to be over-conservative approaches to noise correction and need to be adapted to the specific purpose.

A connectivity-based parcellation improved functional representation of the human cerebellum

The cerebellum is traditionally well known for its role in motor learning and coordination. Recently, it is recognized that the function of the cerebellum is highly diverse and extends to non-motor domains, such as working memory, emotion and language. The diversity of the cerebellum can be appreciated by examining its extensive connectivity to the cerebral regions selective for both motor and cognitive functions. Importantly, the pattern of cerebro-cerebellar connectivity is specific and distinct to different cerebellar subregions. Therefore, to understand the cerebellum and the various functions it involves, it is essential to identify and differentiate its subdivisions. However, most studies are still referring the cerebellum as one brain structure or by its gross anatomical subdivisions, which does not necessarily reflect the functional mapping of the cerebellum. We here employed a data-driven method to generate a functional connectivity-based parcellation of the cerebellum. Our results demonstrated that functional connectivity-based atlas is superior to existing atlases in regards to cluster homogeneity, accuracy of functional connectivity representation and individual identification. Furthermore, our functional atlas improves statistical results of task fMRI analyses, as compared to the standard voxel-based approach and existing atlases. Our detailed functional parcellation provides a valuable tool for elucidating the functional diversity and connectivity of the cerebellum as well as its network relationships with the whole brain.

The resting-state fMRI arterial signal predicts differential blood transit time through the brain

Previous studies have found that aperiodic, systemic low-frequency oscillations (sLFOs) are present in blood-oxygen-level-dependent (BOLD) data. These signals are in the same low frequency band as the "resting state" signal; however, they are distinct signals which represent non-neuronal, physiological oscillations. The same sLFOs are found in the periphery (i.e. finger tips) as changes in oxy/deoxy-hemoglobin concentration using concurrent near-infrared spectroscopy. Together, this evidence points toward an extra-cerebral origin of these sLFOs. If this is the case, it is expected that these sLFO signals would be found in the carotid arteries with time delays that precede the signals found in the brain. To test this hypothesis, we employed the publicly available MyConnectome dataset (a two-year longitudinal study of a single subject) to extract the sLFOs in the internal carotid arteries (ICAs) with the help of the T1/T2-weighted images. Significant, but negative, correlations were found between the LFO BOLD signals from the ICAs and (1) the global signal (GS), (2) the superior sagittal sinus, and (3) the jugulars. We found the consistent time delays between the sLFO signals from ICAs, GS and veins which coincide with the blood transit time through the cerebral vascular tree.

Retrospective Correction of Physiological Noise: Impact on Sensitivity, Specificity, and Reproducibility of Resting-State Functional Connectivity in a Reading Network Model

It is well accepted that physiological noise (PN) obscures the detection of neural fluctuations in resting-state functional connectivity (rsFC) magnetic resonance imaging. However, a clear consensus for an optimal PN correction (PNC) methodology and how it can impact the rsFC signal characteristics is still lacking. In this study, we probe the impact of three PNC methods: RETROICOR: (Glover et al., 2000 ), ANATICOR: (Jo et al., 2010 ), and RVTMBPM: (Bianciardi et al., 2009 ). Using a reading network model, we systematically explore the effects of PNC optimization on sensitivity, specificity, and reproducibility of rsFC signals. In terms of specificity, ANATICOR was found to be effective in removing local white matter (WM) fluctuations and also resulted in aggressive removal of expected cortical-to-subcortical functional connections. The ability of RETROICOR to remove PN was equivalent to removal of simulated random PN such that it artificially inflated the connection strength, thereby decreasing sensitivity. RVTMBPM maintained specificity and sensitivity by balanced removal of vasodilatory PN and local WM nuisance edges. Another aspect of this work was exploring the effects of PNC on identifying reading group differences. Most PNC methods accounted for between-subject PN variability resulting in reduced intersession reproducibility. This effect facilitated the detection of the most consistent group differences. RVTMBPM was most effective in detecting significant group differences due to its inherent sensitivity to removing spatially structured and temporally repeating PN arising from dense vasculature. Finally, results suggest that combining all three PNC resulted in "overcorrection" by removing signal along with noise.

A three-wave longitudinal study of subcortical-cortical resting-state connectivity in adolescence: Testing age- and puberty-related changes

Adolescence is the transitional period between childhood and adulthood, characterized by substantial changes in reward‐driven behavior. Although reward‐driven behavior is supported by subcortical‐medial prefrontal cortex (PFC) connectivity, the development of these circuits is not well understood. Particularly, while puberty has been hypothesized to accelerate organization and activation of functional neural circuits, the relationship between age, sex, pubertal change, and functional connectivity has hardly been studied. Here, we present an analysis of resting‐state functional connectivity between subcortical structures and the medial PFC, in 661 scans of 273 participants between 8 and 29 years, using a three‐wave longitudinal design. Generalized additive mixed model procedures were used to assess the effects of age, sex, and self‐reported pubertal status on connectivity between subcortical structures (nucleus accumbens, caudate, putamen, hippocampus, and amygdala) and cortical medial structures (dorsal anterior cingulate, ventral anterior cingulate, subcallosal cortex, frontal medial cortex). We observed an age‐related strengthening of subcortico‐subcortical and cortico‐cortical connectivity. Subcortical–cortical connectivity, such as, between the nucleus accumbens—frontal medial cortex, and the caudate—dorsal anterior cingulate cortex, however, weakened across age. Model‐based comparisons revealed that for specific connections pubertal development described developmental change better than chronological age. This was particularly the case for changes in subcortical–cortical connectivity and distinctively for boys and girls. Together, these findings indicate changes in functional network strengthening with pubertal development. These changes in functional connectivity may maximize the neural efficiency of interregional communication and set the stage for further inquiry of biological factors driving adolescent functional connectivity changes.

Influence of respiration-induced B0 variations in real-time phase-contrast echo planar imaging of the cervical cerebrospinal fluid

PURPOSE

Respiration induces temporal variations of the main magnetic field B₀ along the spinal cord. These variations are typically not compensated for in velocity quantifications using phase-contrast MRI. The goal of this study was to analyze errors caused by respiration-induced B₀ variations in real-time phase-contrast echo planar imaging (PCEPI) of cervical cerebrospinal fluid (CSF) velocity measurements and to evaluate this effect for various sequence parameters using numerical simulations.

METHODS

Real-time B₀ measurements with double gradient echo sequence and PCEPI measurements were acquired in the cervical CSF of 10 healthy subjects. Dynamic phase offsets attributed to respiration-induced B₀ variations were analyzed by quantifying amplitudes and comparing the temporal behavior with respiratory signals. In experiments and simulations, the influence of the echo time (TE) and the delay between PCEPI images (Δt) with respect to respiration on the dynamic phase offsets were investigated.

RESULTS

A good agreement was found between phase offsets extracted from both acquisition types. Furthermore, respiratory signals qualitatively matched the temporal behavior of the measured phase offsets showing a dependency on subject-dependent local B₀ distribution and respiration physiology. Simulations revealed residual background phases in PCEPI velocity quantification varying with TE and Δt.

CONCLUSION

Respiration-induced B₀ variations result in dynamic background phases in real-time PCEPI velocity quantifications of the CSF in the cervical spine. The current work underlines that these background phases need to be corrected to avoid confounding effects.

Role of the striatum in incidental learning of sound categories

Humans are born as “universal listeners.” However, over the first year, infants’ perception is shaped by native speech categories. How do these categories naturally emerge without explicit training or overt feedback? Using fMRI, we examined the neural basis of incidental sound category learning as participants played a videogame in which sound category exemplars had functional utility in guiding videogame success. Even without explicit categorization of the sounds, participants learned functionally relevant sound categories that generalized to novel exemplars when exemplars had an organized distributional structure. Critically, the striatum was engaged and functionally connected to the auditory cortex during game play, and this activity and connectivity predicted the learning outcome. These findings elucidate the neural mechanism by which humans incidentally learn “real-world” categories.

Humans are born as “universal listeners” without a bias toward any particular language. However, over the first year of life, infants’ perception is shaped by learning native speech categories. Acoustically different sounds—such as the same word produced by different speakers—come to be treated as functionally equivalent. In natural environments, these categories often emerge incidentally without overt categorization or explicit feedback. However, the neural substrates of category learning have been investigated almost exclusively using overt categorization tasks with explicit feedback about categorization decisions. Here, we examined whether the striatum, previously implicated in category learning, contributes to incidental acquisition of sound categories. In the fMRI scanner, participants played a videogame in which sound category exemplars aligned with game actions and events, allowing sound categories to incidentally support successful game play. An experimental group heard nonspeech sound exemplars drawn from coherent category spaces, whereas a control group heard acoustically similar sounds drawn from a less structured space. Although the groups exhibited similar in-game performance, generalization of sound category learning and activation of the posterior striatum were significantly greater in the experimental than control group. Moreover, the experimental group showed brain–behavior relationships related to the generalization of all categories, while in the control group these relationships were restricted to the categories with structured sound distributions. Together, these results demonstrate that the striatum, through its interactions with the left superior temporal sulcus, contributes to incidental acquisition of sound category representations emerging from naturalistic learning environments.

Ultra-high field MRI reveals mood-related circuit disturbances in depression: a comparison between 3-Tesla and 7-Tesla

Ultra-high field 7-Tesla (7 T) MRI has the potential to advance our understanding of neuropsychiatric disorders, including major depressive disorder (MDD). To date, few studies have quantified the advantage of resting state functional MRI (fMRI) at 7 T compared to 3-Tesla (3 T). We conducted a series of experiments that demonstrate the improvement in temporal signal-to-noise ratio (TSNR) of a multi-echo multi-band fMRI protocol with ultra-high field 7 T MRI, compared to a similar protocol using 3 T MRI in healthy controls (HC). We also directly tested the enhancement in ultra-high field 7 T fMRI signal power by examining the ventral tegmental area (VTA), a small midbrain structure that is critical to the expected neuropathology of MDD but difficult to discern with standard 3 T MRI. We demonstrate up to 300% improvement in TSNR and resting state functional connectivity coefficients provided by ultra-high field 7 T fMRI compared to 3 T, indicating enhanced power for detection of functional neural architecture. A multi-echo based acquisition protocol and signal denoising pipeline afforded greater gain in signal power compared to classic acquisition and denoising pipelines. Furthermore, ultra-high field fMRI revealed mood-related neurocircuit disturbances in patients with MDD compared to HC, which were not detectable with 3 T fMRI. Ultra-high field 7 T fMRI may provide an effective tool for studying functional neural architecture relevant to MDD and other neuropsychiatric disorders.

Evaluation of confound regression strategies for the mitigation of micromovement artifact in studies of dynamic resting-state functional connectivity and multilayer network modularity

Dynamic functional connectivity reflects the spatiotemporal organization of spontaneous brain activity in health and disease. Dynamic functional connectivity may be susceptible to artifacts induced by participant motion. This report provides a systematic evaluation of 12 commonly used participant-level confound regression strategies designed to mitigate the effects of micromovements in a sample of 393 youths (ages 8-22 years). Each strategy was evaluated according to a number of benchmarks, including (a) the residual association between participant motion and edge dispersion, (b) distance-dependent effects of motion on edge dispersion, (c) the degree to which functional subnetworks could be identified by multilayer modularity maximization, and (d) measures of module reconfiguration, including node flexibility and node promiscuity. Results indicate variability in the effectiveness of the evaluated pipelines across benchmarks. Methods that included global signal regression were the most consistently effective de-noising strategies.

Beware detrending: Optimal preprocessing pipeline for low-frequency fluctuation analysis

Resting‐state functional magnetic resonance imaging (rs‐fMRI) offers the possibility to assess brain function independent of explicit tasks and individual performance. This absence of explicit stimuli in rs‐fMRI makes analyses more susceptible to nonneural signal fluctuations than task‐based fMRI. Data preprocessing is a critical procedure to minimise contamination by artefacts related to motion and physiology. We herein investigate the effects of different preprocessing strategies on the amplitude of low‐frequency fluctuations (ALFFs) and its fractional counterpart, fractional ALFF (fALFF). Sixteen artefact reduction schemes based on nuisance regression are applied to data from 82 subjects acquired at 1.5 T, 30 subjects at 3 T, and 23 subjects at 7 T, respectively. In addition, we examine test–retest variance and effects of bias correction. In total, 569 data sets are included in this study. Our results show that full artefact reduction reduced test–retest variance by up to 50%. Polynomial detrending of rs‐fMRI data has a positive effect on group‐level t‐values for ALFF but, importantly, a negative effect for fALFF. We show that the normalisation process intrinsic to fALFF calculation causes the observed reduction and introduce a novel measure for low‐frequency fluctuations denoted as high‐frequency ALFF (hfALFF). We demonstrate that hfALFF values are not affected by the negative detrending effects seen in fALFF data. Still, highest grey matter (GM) group‐level t‐values were obtained for fALFF data without detrending, even when compared to an exploratory detrending approach based on autocorrelation measures. From our results, we recommend the use of full nuisance regression including polynomial detrending in ALFF data, but to refrain from using polynomial detrending in fALFF data. Such optimised preprocessing increases GM group‐level t‐values by up to 60%.

Limitations of Resting-State Functional MR Imaging in the Setting of Focal Brain Lesions

Methods of image acquisition and analysis for resting-state functional MR imaging (rsfMR imaging) are still evolving. Neurovascular uncoupling and susceptibility artifact are important confounds of rsfMR imaging in the setting of focal brain lesions such as brain tumors. This article reviews the detection of these confounds using rsfMR imaging metrics in the setting of focal brain lesions. In the near future, with the wide range of ongoing research in rsfMR imaging, these issues likely will be overcome and will open new windows into brain function and connectivity.

Reprint of 'Noise contributions to the fMRI signal: An Overview'

The ability to discriminate signal from noise plays a key role in the analysis and interpretation of functional magnetic resonance imaging (fMRI) measures of brain activity. Over the past two decades, a number of major sources of noise have been identified, including system-related instabilities, subject motion, and physiological fluctuations. This article reviews the characteristics of the various noise sources as well as the mechanisms through which they affect the fMRI signal. Approaches for distinguishing signal from noise and the associated challenges are also reviewed. These challenges reflect the fact that some noise sources, such as respiratory activity, are generated by the same underlying brain networks that give rise to functional signals that are of interest.

Similar neural responses predict friendship

Human social networks are overwhelmingly homophilous: individuals tend to befriend others who are similar to them in terms of a range of physical attributes (e.g., age, gender). Do similarities among friends reflect deeper similarities in how we perceive, interpret, and respond to the world? To test whether friendship, and more generally, social network proximity, is associated with increased similarity of real-time mental responding, we used functional magnetic resonance imaging to scan subjects’ brains during free viewing of naturalistic movies. Here we show evidence for neural homophily: neural responses when viewing audiovisual movies are exceptionally similar among friends, and that similarity decreases with increasing distance in a real-world social network. These results suggest that we are exceptionally similar to our friends in how we perceive and respond to the world around us, which has implications for interpersonal influence and attraction.

Though we are often friends with people similar to ourselves, it is unclear if neural responses to perceptual stimuli are also similar. Here, authors show that the similarity of neural responses evoked by a range of videos was highest for close friends and decreased with increasing social distance.

Serial correlations in single-subject fMRI with sub-second TR

When performing statistical analysis of single-subject fMRI data, serial correlations need to be taken into account to allow for valid inference. Otherwise, the variability in the parameter estimates might be under-estimated resulting in increased false-positive rates. Serial correlations in fMRI data are commonly characterized in terms of a first-order autoregressive (AR) process and then removed via pre-whitening. The required noise model for the pre-whitening depends on a number of parameters, particularly the repetition time (TR). Here we investigate how the sub-second temporal resolution provided by simultaneous multislice (SMS) imaging changes the noise structure in fMRI time series. We fit a higher-order AR model and then estimate the optimal AR model order for a sequence with a TR of less than 600 ms providing whole brain coverage. We show that physiological noise modelling successfully reduces the required AR model order, but remaining serial correlations necessitate an advanced noise model. We conclude that commonly used noise models, such as the AR(1) model, are inadequate for modelling serial correlations in fMRI using sub-second TRs. Rather, physiological noise modelling in combination with advanced pre-whitening schemes enable valid inference in single-subject analysis using fast fMRI sequences.

Dual-echo ASL contributes to decrypting the link between functional connectivity and cerebral blow flow

Arterial spin labeling (ASL) MRI with a dual-echo readout module (DE-ASL) enables noninvasive simultaneous acquisition of cerebral blood flow (CBF)-weighted images and blood oxygenation level dependent (BOLD) contrast. Up to date, resting-state functional connectivity (FC) studies based on CBF fluctuations have been very limited, while the BOLD is still the method most frequently used. The purposes of this technical report were (i) to assess the potentiality of the DE-ASL sequence for the quantification of resting-state FC and brain organization, with respect to the conventional BOLD (cvBOLD) and (ii) to investigate the relationship between a series of complex network measures and the CBF information. Thirteen volunteers were scanned on a 3 T scanner acquiring a pseudocontinuous multislice DE-ASL sequence, from which the concomitant BOLD (ccBOLD) simultaneously to the ASL can be extracted. In the proposed comparison, the brain FC and graph-theoretical analysis were used for quantifying the connectivity strength between pairs of regions and for assessing the network model properties in all the sequences. The main finding was that the ccBOLD part of the DE-ASL sequence provided highly comparable connectivity results compared to cvBOLD. As expected, because of its different nature, ASL sequence showed different patterns of brain connectivity and graph indices compared to BOLD sequences. To conclude, the resting-state FC can be reliably detected using DE-ASL, simultaneously to CBF quantifications, whereas a single fMRI experiment precludes the quantitative measurement of BOLD signal changes. Hum Brain Mapp 38:5831-5844, 2017. © 2017 Wiley Periodicals, Inc.

Methods for cleaning the BOLD fMRI signal

Blood oxygen-level-dependent functional magnetic resonance imaging (BOLD fMRI) has rapidly become a popular technique for the investigation of brain function in healthy individuals, patients as well as in animal studies. However, the BOLD signal arises from a complex mixture of neuronal, metabolic and vascular processes, being therefore an indirect measure of neuronal activity, which is further severely corrupted by multiple non-neuronal fluctuations of instrumental, physiological or subject-specific origin. This review aims to provide a comprehensive summary of existing methods for cleaning the BOLD fMRI signal. The description is given from a methodological point of view, focusing on the operation of the different techniques in addition to pointing out the advantages and limitations in their application. Since motion-related and physiological noise fluctuations are two of the main noise components of the signal, techniques targeting their removal are primarily addressed, including both data-driven approaches and using external recordings. Data-driven approaches, which are less specific in the assumed model and can simultaneously reduce multiple noise fluctuations, are mainly based on data decomposition techniques such as principal and independent component analysis. Importantly, the usefulness of strategies that benefit from the information available in the phase component of the signal, or in multiple signal echoes is also highlighted. The use of global signal regression for denoising is also addressed. Finally, practical recommendations regarding the optimization of the preprocessing pipeline for the purpose of denoising and future venues of research are indicated. Through the review, we summarize the importance of signal denoising as an essential step in the analysis pipeline of task-based and resting state fMRI studies.

Impact of automated ICA-based denoising of fMRI data in acute stroke patients

Different strategies have been developed using Independent Component Analysis (ICA) to automatically de-noise fMRI data, either focusing on removing only certain components (e.g. motion-ICA-AROMA, Pruim et al., 2015a) or using more complex classifiers to remove multiple types of noise components (e.g. FIX, Salimi-Khorshidi et al., 2014 Griffanti et al., 2014). However, denoising data obtained in an acute setting might prove challenging: the presence of multiple noise sources may not allow focused strategies to clean the data enough and the heterogeneity in the data may be so great to critically undermine complex approaches. The purpose of this study was to explore what automated ICA based approach would better cope with these limitations when cleaning fMRI data obtained from acute stroke patients. The performance of a focused classifier (ICA-AROMA) and a complex classifier (FIX) approaches were compared using data obtained from twenty consecutive acute lacunar stroke patients using metrics determining RSN identification, RSN reproducibility, changes in the BOLD variance, differences in the estimation of functional connectivity and loss of temporal degrees of freedom. The use of generic-trained FIX resulted in misclassification of components and significant loss of signal (< 80%), and was not explored further. Both ICA-AROMA and patient-trained FIX based denoising approaches resulted in significantly improved RSN reproducibility (p < 0.001), localized reduction in BOLD variance consistent with noise removal, and significant changes in functional connectivity (p < 0.001). Patient-trained FIX resulted in higher RSN identifiability (p < 0.001) and wider changes both in the BOLD variance and in functional connectivity compared to ICA-AROMA. The success of ICA-AROMA suggests that by focusing on selected components the full automation can deliver meaningful data for analysis even in population with multiple sources of noise. However, the time invested to train FIX with appropriate patient data proved valuable, particularly in improving the signal-to-noise ratio.

The global signal in fMRI: Nuisance or Information?

The global signal is widely used as a regressor or normalization factor for removing the effects of global variations in the analysis of functional magnetic resonance imaging (fMRI) studies. However, there is considerable controversy over its use because of the potential bias that can be introduced when it is applied to the analysis of both task-related and resting-state fMRI studies. In this paper we take a closer look at the global signal, examining in detail the various sources that can contribute to the signal. For the most part, the global signal has been treated as a nuisance term, but there is growing evidence that it may also contain valuable information. We also examine the various ways that the global signal has been used in the analysis of fMRI data, including global signal regression, global signal subtraction, and global signal normalization. Furthermore, we describe new ways for understanding the effects of global signal regression and its relation to the other approaches.

Analysis and correction of field fluctuations in fMRI data using field monitoring

This work investigates the role of magnetic field fluctuations as a confound in fMRI. In standard fMRI experiments with single-shot EPI acquisition at 3 Tesla the uniform and gradient components of the magnetic field were recorded with NMR field sensors. By principal component analysis it is found that differences of field evolution between the EPI readouts are explainable by few components relating to slow and within-shot field dynamics of hardware and physiological origin. The impact of fluctuating field components is studied by selective data correction and assessment of its influence on image fluctuation and SFNR. Physiological field fluctuations, attributed to breathing, were found to be small relative to those of hardware origin. The dominant confounds were hardware-related and attributable to magnet drift and thermal changes. In raw image time series, field fluctuation caused significant SFNR loss, reflected by a 67% gain upon correction. Large part of this correction can be accomplished by traditional image realignment, which addresses slow and spatially uniform field changes. With realignment, explicit field correction increased the SFNR on the order of 6%. In conclusion, field fluctuations are a relevant confound in fMRI and can be addressed effectively by retrospective data correction. Based on the physics involved it is anticipated that the advantage of full field correction increases with field strength, with non-Cartesian readouts, and upon phase-sensitive BOLD analysis.

Noise contributions to the fMRI signal: An overview

The ability to discriminate signal from noise plays a key role in the analysis and interpretation of functional magnetic resonance imaging (fMRI) measures of brain activity. Over the past two decades, a number of major sources of noise have been identified, including system-related instabilities, subject motion, and physiological fluctuations. This article reviews the characteristics of the various noise sources as well as the mechanisms through which they affect the fMRI signal. Approaches for distinguishing signal from noise and the associated challenges are also reviewed. These challenges reflect the fact that some noise sources, such as respiratory activity, are generated by the same underlying brain networks that give rise to functional signals that are of interest.

Noise Reduction in Arterial Spin Labeling Based Functional Connectivity Using Nuisance Variables

Arterial Spin Labeling (ASL) perfusion image series have recently been utilized for functional connectivity (FC) analysis in healthy volunteers and children with autism spectrum disorders (ASD). Noise reduction by using nuisance variables has been shown to be necessary to minimize potential confounding effects of head motion and physiological signals on BOLD based FC analysis. The purpose of the present study is to systematically evaluate the effectiveness of different noise reduction strategies (NRS) using nuisance variables to improve perfusion based FC analysis in two cohorts of healthy adults using state of the art 3D background-suppressed (BS) GRASE pseudo-continuous ASL (pCASL) and dual-echo 2D-EPI pCASL sequences. Five different NRS were performed in healthy volunteers to compare their performance. We then compared seed-based FC analysis using 3D BS GRASE pCASL in a cohort of 12 children with ASD (3f/9m, age 12.8 ± 1.3 years) and 13 typically developing (TD) children (1f/12m; age 13.9 ± 3 years) in conjunction with NRS. Regression of different combinations of nuisance variables affected FC analysis from a seed in the posterior cingulate cortex (PCC) to other areas of the default mode network (DMN) in both BOLD and pCASL data sets. Consistent with existing literature on BOLD-FC, we observed improved spatial specificity after physiological noise reduction and improved long-range connectivity using head movement related regressors. Furthermore, 3D BS GRASE pCASL shows much higher temporal SNR compared to dual-echo 2D-EPI pCASL and similar effects of noise reduction as those observed for BOLD. Seed-based FC analysis using 3D BS GRASE pCASL in children with ASD and TD children showed that noise reduction including physiological and motion related signals as nuisance variables is crucial for identifying altered long-range connectivity from PCC to frontal brain areas associated with ASD. This is the first study that systematically evaluated the effects of different NRS on ASL based FC analysis. 3D BS GRASE pCASL is the preferred ASL sequence for FC analysis due to its superior temporal SNR. Removing physiological noise and motion parameters is critical for detecting altered FC in neurodevelopmental disorders such as ASD.

Resting-state test-retest reliability of a priori defined canonical networks over different preprocessing steps

Resting-state functional connectivity analysis has become a widely used method for the investigation of human brain connectivity and pathology. The measurement of neuronal activity by functional MRI, however, is impeded by various nuisance signals that reduce the stability of functional connectivity. Several methods exist to address this predicament, but little consensus has yet been reached on the most appropriate approach. Given the crucial importance of reliability for the development of clinical applications, we here investigated the effect of various confound removal approaches on the test-retest reliability of functional-connectivity estimates in two previously defined functional brain networks. Our results showed that gray matter masking improved the reliability of connectivity estimates, whereas denoising based on principal components analysis reduced it. We additionally observed that refraining from using any correction for global signals provided the best test-retest reliability, but failed to reproduce anti-correlations between what have been previously described as antagonistic networks. This suggests that improved reliability can come at the expense of potentially poorer biological validity. Consistent with this, we observed that reliability was proportional to the retained variance, which presumably included structured noise, such as reliable nuisance signals (for instance, noise induced by cardiac processes). We conclude that compromises are necessary between maximizing test-retest reliability and removing variance that may be attributable to non-neuronal sources.

Abnormal functional connectivity and cortical integrity influence dominant hand motor disability in multiple sclerosis: a multimodal analysis

Functional reorganization and structural damage occur in the brains of people with multiple sclerosis (MS) throughout the disease course. However, the relationship between resting-state functional connectivity (FC) reorganization in the sensorimotor network and motor disability in MS is not well understood. This study used resting-state fMRI, T1-weighted and T2-weighted, and magnetization transfer (MT) imaging to investigate the relationship between abnormal FC in the sensorimotor network and upper limb motor disability in people with MS, as well as the impact of disease-related structural abnormalities within this network. Specifically, the differences in FC of the left hemisphere hand motor region between MS participants with preserved (n = 17) and impaired (n = 26) right hand function, compared with healthy controls (n = 20) was investigated. Differences in brain atrophy and MT ratio measured at the global and regional levels were also investigated between the three groups. Motor preserved MS participants had stronger FC in structurally intact visual information processing regions relative to motor impaired MS participants. Motor impaired MS participants showed weaker FC in the sensorimotor and somatosensory association cortices and more severe structural damage throughout the brain compared with the other groups. Logistic regression analysis showed that regional MTR predicted motor disability beyond the impact of global atrophy whereas regional grey matter volume did not. More importantly, as the first multimodal analysis combining resting-state fMRI, T1-weighted, T2-weighted and MTR images in MS, we demonstrate how a combination of structural and functional changes may contribute to motor impairment or preservation in MS. Hum Brain Mapp 37:4262-4275, 2016. © 2016 Wiley Periodicals, Inc.

Subgenual Cingulate-Amygdala Functional Disconnection and Vulnerability to Melancholic Depression

The syndromic heterogeneity of major depressive disorder (MDD) hinders understanding of the etiology of predisposing vulnerability traits and underscores the importance of identifying neurobiologically valid phenotypes. Distinctive fMRI biomarkers of vulnerability to MDD subtypes are currently lacking. This study investigated whether remitted melancholic MDD patients, who are at an elevated lifetime risk for depressive episodes, demonstrate distinctive patterns of resting-state connectivity with the subgenual cingulate cortex (SCC), known to be of core pathophysiological importance for severe and familial forms of MDD. We hypothesized that patterns of disrupted SCC connectivity would be a distinguishing feature of melancholia. A total of 63 medication-free remitted MDD (rMDD) patients (33 melancholic and 30 nonmelancholic) and 39 never-depressed healthy controls (HC) underwent resting-state fMRI scanning. SCC connectivity was investigated with closely connected bilateral a priori regions of interest (ROIs) relevant to MDD (anterior temporal, ventromedial prefrontal, dorsomedial prefrontal cortices, amygdala, hippocampus, septal region, and hypothalamus). Decreased (less positive) SCC connectivity with the right parahippocampal gyrus and left amygdala distinguished melancholic rMDD patients from the nonmelancholic rMDD and HC groups (cluster-based familywise error-corrected p⩽0.007 over individual a priori ROIs corresponding to approximate Bonferroni-corrected p⩽0.05 across all seven a priori ROIs). No areas demonstrating increased (more positive) connectivity were observed. Abnormally decreased connectivity of the SCC with the amygdala and parahippocampal gyrus distinguished melancholic from nonmelancholic rMDD. These results provide the first resting-state neural signature distinctive of melancholic rMDD and may reflect a subtype-specific primary vulnerability factor given a lack of association with the number of previous episodes.

Progressive Bidirectional Age-Related Changes in Default Mode Network Effective Connectivity across Six Decades

The default mode network (DMN) is a set of regions that is tonically engaged during the resting state and exhibits task-related deactivation that is readily reproducible across a wide range of paradigms and modalities. The DMN has been implicated in numerous disorders of cognition and, in particular, in disorders exhibiting age-related cognitive decline. Despite these observations, investigations of the DMN in normal aging are scant. Here, we used blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) acquired during rest to investigate age-related changes in functional connectivity of the DMN in 120 healthy normal volunteers comprising six, 20-subject, decade cohorts (from 20–29 to 70–79). Structural equation modeling (SEM) was used to assess age-related changes in inter-regional connectivity within the DMN. SEM was applied both using a previously published, meta-analytically derived, node-and-edge model, and using exploratory modeling searching for connections that optimized model fit improvement. Although the two models were highly similar (only 3 of 13 paths differed), the sample demonstrated significantly better fit with the exploratory model. For this reason, the exploratory model was used to assess age-related changes across the decade cohorts. Progressive, highly significant changes in path weights were found in 8 (of 13) paths: four rising, and four falling (most changes were significant by the third or fourth decade). In all cases, rising paths and falling paths projected in pairs onto the same nodes, suggesting compensatory increases associated with age-related decreases. This study demonstrates that age-related changes in DMN physiology (inter-regional connectivity) are bidirectional, progressive, of early onset and part of normal aging.