Physiological and motion signatures in static and time-varying functional connectivity and their subject identifiability

Functional MRI of the Brainstem for Assessing Its Autonomic Functions: From Imaging Parameters and Analysis to Functional Atlas

BACKGROUND

The brainstem is a crucial component of the central autonomic nervous (CAN) system. Functional MRI (fMRI) of the brainstem remains challenging due to a range of factors, including diverse imaging protocols, analysis, and interpretation.

PURPOSE

To develop an fMRI protocol for establishing a functional atlas in the brainstem.

STUDY TYPE

Prospective cross-sectional study.

SUBJECTS

Ten healthy subjects (four males, six females).

FIELD STRENGTH/SEQUENCE

Using a 3.0 Tesla MR scanner, we acquired T1-weighted images and three different fMRI scans using fMRI protocols of the optimized functional Imaging of Brainstem (FIBS), the Human Connectome Project (HCP), and the Adolescent Brain Cognitive Development (ABCD) project.

ASSESSMENT

The temporal signal-to-noise-ratio (TSNR) of fMRI data was compared between the FIBS, HCP, and ABCD protocols. Additionally, the main normalization algorithms (i.e., FSL-FNIRT, SPM-DARTEL, and ANTS-SyN) were compared to identify the best approach to normalize brainstem data using root-mean-square (RMS) error computed based on manually defined reference points. Finally, a functional autonomic brainstem atlas that maps brainstem regions involved in the CAN system was defined using meta-analysis and data-driven approaches.

STATISTICAL TESTS

ANOVA was used to compare the performance of different imaging and preprocessing pipelines with multiple comparison corrections (P ≤ 0.05). Dice coefficient estimated ROI overlap, with 50% overlap between ROIs identified in each approach considered significant.

RESULTS

The optimized FIBS protocol showed significantly higher brainstem TSNR than the HCP and ABCD protocols (P ≤ 0.05). Furthermore, FSL-FNIRT RMS error (2.1 ± 1.22 mm; P ≤ 0.001) exceeded SPM (1.5 ± 0.75 mm; P ≤ 0.01) and ANTs (1.1 ± 0.54 mm). Finally, a set of 12 final brainstem ROIs with dice coefficient ≥0.50, as a step toward the development of a functional brainstem atlas.

DATA CONCLUSION

The FIBS protocol yielded more robust brainstem CAN results and outperformed both the HCP and ABCD protocols.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 1.

Cortical timescales and the modular organization of structural and functional brain networks

Recent years have seen growing interest in characterizing the properties of regional brain dynamics and their relationship to other features of brain structure and function. In particular, multiple studies have observed regional differences in the "timescale" over which activity fluctuates during periods of quiet rest. In the cerebral cortex, these timescales have been associated with both local circuit properties as well as patterns of inter-regional connectivity, including the extent to which each region exhibits widespread connectivity to other brain areas. In the current study, we build on prior observations of an association between connectivity and dynamics in the cerebral cortex by investigating the relationship between BOLD fMRI timescales and the modular organization of structural and functional brain networks. We characterize network community structure across multiple scales and find that longer timescales are associated with greater within-community functional connectivity and diverse structural connectivity. We also replicate prior observations of a positive correlation between timescales and structural connectivity degree. Finally, we find evidence for preferential functional connectivity between cortical areas with similar timescales. We replicate these findings in an independent dataset. These results contribute to our understanding of functional brain organization and structure-function relationships in the human brain, and support the notion that regional differences in cortical dynamics may in part reflect the topological role of each region within macroscale brain networks.

On the variability of dynamic functional connectivity assessment methods

BACKGROUND

Dynamic functional connectivity (dFC) has become an important measure for understanding brain function and as a potential biomarker. However, various methodologies have been developed for assessing dFC, and it is unclear how the choice of method affects the results. In this work, we aimed to study the results variability of commonly used dFC methods.

METHODS

We implemented 7 dFC assessment methods in Python and used them to analyze the functional magnetic resonance imaging data of 395 subjects from the Human Connectome Project. We measured the similarity of dFC results yielded by different methods using several metrics to quantify overall, temporal, spatial, and intersubject similarity.

RESULTS

Our results showed a range of weak to strong similarity between the results of different methods, indicating considerable overall variability. Somewhat surprisingly, the observed variability in dFC estimates was found to be comparable to the expected functional connectivity variation over time, emphasizing the impact of methodological choices on the final results. Our findings revealed 3 distinct groups of methods with significant intergroup variability, each exhibiting distinct assumptions and advantages.

CONCLUSIONS

Overall, our findings shed light on the impact of dFC assessment analytical flexibility and highlight the need for multianalysis approaches and careful method selection to capture the full range of dFC variation. They also emphasize the importance of distinguishing neural-driven dFC variations from physiological confounds and developing validation frameworks under a known ground truth. To facilitate such investigations, we provide an open-source Python toolbox, PydFC, which facilitates multianalysis dFC assessment, with the goal of enhancing the reliability and interpretability of dFC studies.

Discovering individual fingerprints in resting-state functional connectivity using deep neural networks

Non-negligible idiosyncrasy due to interindividual differences is an ongoing issue in resting-state functional MRI (rfMRI) analysis. We show that a deep neural network (DNN) can be employed for individual identification by learning important features from the time-varying functional connectivity (FC) of rfMRI in the Human Connectome Project. We employed the trained DNN to identify individuals from an independent dataset acquired at our institution. The results revealed that the DNN could successfully identify 300 individuals with an error rate of 2.9% using 15 s time-window and 870 individuals with an error rate of 6.7%. A trained DNN with nonlinear hidden layers led to the proposal of the "fingerprint of FC" (fpFC) as representative edges of individual FC. The fpFCs for individuals exhibited commonly important and individual-specific edges across time-window lengths (from 5 min to 15 s). Furthermore, the utility of our model for another group of subjects was validated, supporting the feasibility of our technique in the context of transfer learning. In conclusion, our study offers an insight into the discovery of the intrinsic mode of the human brain using whole-brain resting-state FC and DNNs.

Intermediately synchronised brain states optimise trade-off between subject specificity and predictive capacity

Functional connectivity (FC) refers to the statistical dependencies between activity of distinct brain areas. To study temporal fluctuations in FC within the duration of a functional magnetic resonance imaging (fMRI) scanning session, researchers have proposed the computation of an edge time series (ETS) and their derivatives. Evidence suggests that FC is driven by a few time points of high-amplitude co-fluctuation (HACF) in the ETS, which may also contribute disproportionately to interindividual differences. However, it remains unclear to what degree different time points actually contribute to brain-behaviour associations. Here, we systematically evaluate this question by assessing the predictive utility of FC estimates at different levels of co-fluctuation using machine learning (ML) approaches. We demonstrate that time points of lower and intermediate co-fluctuation levels provide overall highest subject specificity as well as highest predictive capacity of individual-level phenotypes.

Deep-learning-enabled brain hemodynamic mapping using resting-state fMRI

Cerebrovascular disease is a leading cause of death globally. Prevention and early intervention are known to be the most effective forms of its management. Non-invasive imaging methods hold great promises for early stratification, but at present lack the sensitivity for personalized prognosis. Resting-state functional magnetic resonance imaging (rs-fMRI), a powerful tool previously used for mapping neural activity, is available in most hospitals. Here we show that rs-fMRI can be used to map cerebral hemodynamic function and delineate impairment. By exploiting time variations in breathing pattern during rs-fMRI, deep learning enables reproducible mapping of cerebrovascular reactivity (CVR) and bolus arrival time (BAT) of the human brain using resting-state CO2 fluctuations as a natural "contrast media". The deep-learning network is trained with CVR and BAT maps obtained with a reference method of CO2-inhalation MRI, which includes data from young and older healthy subjects and patients with Moyamoya disease and brain tumors. We demonstrate the performance of deep-learning cerebrovascular mapping in the detection of vascular abnormalities, evaluation of revascularization effects, and vascular alterations in normal aging. In addition, cerebrovascular maps obtained with the proposed method exhibit excellent reproducibility in both healthy volunteers and stroke patients. Deep-learning resting-state vascular imaging has the potential to become a useful tool in clinical cerebrovascular imaging.

Brain Connectivity Correlates of Breathing and Cardiac Irregularities in SUDEP: A Resting-State fMRI Study

Sudden unexpected death in epilepsy (SUDEP) is the leading cause of premature mortality among people with epilepsy. Evidence from witnessed and monitored SUDEP cases indicate seizure-induced cardiovascular and respiratory failures; yet, the underlying mechanisms remain obscure. SUDEP occurs often during the night and early morning hours, suggesting that sleep or circadian rhythm-induced changes in physiology contribute to the fatal event. Resting-state fMRI studies have found altered functional connectivity between brain structures involved in cardiorespiratory regulation in later SUDEP cases and in individuals at high-risk of SUDEP. However, those connectivity findings have not been related to changes in cardiovascular or respiratory patterns. Here, we compared fMRI patterns of brain connectivity associated with regular and irregular cardiorespiratory rhythms in SUDEP cases with those of living epilepsy patients of varying SUDEP risk, and healthy controls. We analysed resting-state fMRI data from 98 patients with epilepsy (9 who subsequently succumbed to SUDEP, 43 categorized as low SUDEP risk (no tonic-clonic seizures (TCS) in the year preceding the fMRI scan), and 46 as high SUDEP risk (>3 TCS in the year preceding the scan)) and 25 healthy controls. The global signal amplitude (GSA), defined as the moving standard deviation of the fMRI global signal, was used to identify periods with regular ('low state') and irregular ('high state') cardiorespiratory rhythms. Correlation maps were derived from seeds in twelve regions with a key role in autonomic or respiratory regulation, for the low and high states. Following principal component analysis, component weights were compared between the groups. We found widespread alterations in connectivity of precuneus/posterior cingulate cortex in epilepsy compared to controls, in the low state (regular cardiorespiratory activity). In the low state, and to a lesser degree in the high state, reduced anterior insula connectivity (mainly with anterior and posterior cingulate cortex) in epilepsy appeared, relative to healthy controls. For SUDEP cases, the insula connectivity differences were inversely related to the interval between the fMRI scan and death. The findings suggest that anterior insula connectivity measures may provide a biomarker of SUDEP risk. The neural correlates of autonomic brain structures associated with different cardiorespiratory rhythms may shed light on the mechanisms underlying terminal apnea observed in SUDEP.

Comparing the test-retest reliability of resting-state functional magnetic resonance imaging metrics across single band and multiband acquisitions in the context of healthy aging

The identification of meaningful functional magnetic resonance imaging (fMRI) biomarkers requires measures that reliably capture brain performance across different subjects and over multiple scanning sessions. Recent developments in fMRI acquisition, such as the introduction of multiband (MB) protocols and in-plane acceleration, allow for increased scanning speed and improved temporal resolution. However, they may also lead to reduced temporal signal to noise ratio and increased signal leakage between simultaneously excited slices. These methods have been adopted in several scanning modalities including diffusion weighted imaging and fMRI. To our knowledge, no study has formally compared the reliability of the same resting-state fMRI (rs-fMRI) metrics (amplitude of low-frequency fluctuations; seed-to-voxel and region of interest [ROI]-to-ROI connectivity) across conventional single-band fMRI and different MB acquisitions, with and without in-plane acceleration, across three sessions. In this study, 24 healthy older adults were scanned over three visits, on weeks 0, 1, and 4, and, on each occasion, underwent a conventional single band rs-fMRI scan and three different rs-fMRI scans with MB factors 4 and 6, with and without in-plane acceleration. Across all three rs-fMRI metrics, the reliability scores were highest with MB factor 4 with no in-plane acceleration for cortical areas and with conventional single band for subcortical areas. Recommendations for future research studies are discussed.

Altered Resting State Networks Before and After Temporal Lobe Epilepsy Surgery

OBJECTIVES

To explore the resting state networks (RSNs) alterations in patients with unilateral mesial temporal lobe epilepsy (mTLE) before and after successful surgery.

METHODS

Resting-state functional MRI and T1-weighted structural MRI were obtained in 37 mTLE patients who achieved seizure freedom after anterior temporal lobectomy. Patients were scanned before surgery and at two years after surgery. Twenty-eight age- and sex-matched healthy controls were scanned once. Functional connectivity (FC) changes within and between ten common RSNs before and after surgery, and FC changes between hippocampus and RSNs were explored.

RESULTS

Before surgery, decreased FC was found within visual network and basal ganglia network, while after surgery, FC within basal ganglia network further decreased but FC within sensorimotor network and dorsal attention network increased. Before surgery, between-network FC related to basal ganglia network, visual network and dorsal attention network decreased, while between-network FC related to default mode network increased. After surgery, between-network FC related to visual network and dorsal attention network significantly increased. In addition, before surgery, ipsilateral hippocampus showed decreased FC with visual network, basal ganglia network, sensorimotor network, default mode network and frontoparietal network, while contralateral rostral hippocampus showed increased FC with salience network. After surgery, no obvious FC changes were found between contralateral hippocampus and these RSNs.

CONCLUSION

MTLE patients showed significant RSNs alterations before and after surgery. Basal ganglia network showed progressive decline in functional connectivity. Successful surgery may lead to RSNs reorganization. These results provide preliminary evidence for postoperative functional remodeling at whole-brain-network level.

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.

Extracting electrophysiological correlates of functional magnetic resonance imaging data using the canonical polyadic decomposition

The relation between electrophysiology and BOLD-fMRI requires further elucidation. One approach for studying this relation is to find time-frequency features from electrophysiology that explain the variance of BOLD time-series. Convolution of these features with a canonical hemodynamic response function (HRF) is often required to model neurovascular coupling mechanisms and thus account for time shifts between electrophysiological and BOLD-fMRI data. We propose a framework for extracting the spatial distribution of these time-frequency features while also estimating more flexible, region-specific HRFs. The core component of this method is the decomposition of a tensor containing impulse response functions using the Canonical Polyadic Decomposition. The outputs of this decomposition provide insight into the relation between electrophysiology and BOLD-fMRI and can be used to construct estimates of BOLD time-series. We demonstrated the performance of this method on simulated data while also examining the effects of simulated measurement noise and physiological confounds. Afterwards, we validated our method on publicly available task-based and resting-state EEG-fMRI data. We adjusted our method to accommodate the multisubject nature of these datasets, enabling the investigation of inter-subject variability with regards to EEG-to-BOLD neurovascular coupling mechanisms. We thus also demonstrate how EEG features for modelling the BOLD signal differ across subjects.

Heritability and cross-species comparisons of human cortical functional organization asymmetry

The human cerebral cortex is symmetrically organized along large-scale axes but also presents inter-hemispheric differences in structure and function. The quantified contralateral homologous difference, that is asymmetry, is a key feature of the human brain left-right axis supporting functional processes, such as language. Here, we assessed whether the asymmetry of cortical functional organization is heritable and phylogenetically conserved between humans and macaques. Our findings indicate asymmetric organization along an axis describing a functional trajectory from perceptual/action to abstract cognition. Whereas language network showed leftward asymmetric organization, frontoparietal network showed rightward asymmetric organization in humans. These asymmetries were heritable in humans and showed a similar spatial distribution with macaques, in the case of intra-hemispheric asymmetry of functional hierarchy. This suggests (phylo)genetic conservation. However, both language and frontoparietal networks showed a qualitatively larger asymmetry in humans relative to macaques. Overall, our findings suggest a genetic basis for asymmetry in intrinsic functional organization, linked to higher order cognitive functions uniquely developed in humans.

Resting-state functional MRI signal fluctuation amplitudes are correlated with brain amyloid- deposition in patients with mild cognitive impairment

Mounting evidence suggests that amyloid-β (Aβ) and vascular etiologies are intertwined in the pathogenesis of Alzheimer's disease (AD). Blood-oxygen-level-dependent (BOLD) signals, measured by resting-state functional MRI (rs-fMRI), are associated with neuronal activity and cerebrovascular hemodynamics. Nevertheless, it is unclear if BOLD fluctuations are associated with Aβ deposition in individuals at high risk of AD. Thirty-three patients with amnestic mild cognitive impairment underwent rs-fMRI and AV45 PET. The AV45 standardized uptake value ratio (AV45-SUVR) was calculated using cerebral white matter as reference, to assess Aβ deposition. The whole-brain normalized amplitudes of low-frequency fluctuations (sALFF) of local BOLD signals were calculated in the frequency band of 0.01-0.08 Hz. Stepwise increasing physiological/vascular signal regressions on the rs-fMRI data examined whether sALFF-AV45 correlations were driven by vascular hemodynamics, neuronal activities, or both. We found that sALFF and AV45-SUVR were negatively correlated in regions of default-mode and visual networks (precuneus, angular, lingual and fusiform gyri). Regions with higher sALFF had less Aβ accumulation. Correlated cluster sizes in MNI space (r ≈ -0.47) were reduced from 3018 mm³ to 1072 mm³ with stronger cardiovascular regression. These preliminary findings imply that local brain blood fluctuations due to vascular hemodynamics or neuronal activity can affect Aβ homeostasis.

Hemodynamic and metabolic correspondence of resting-state voxel-based physiological metrics in healthy adults

Voxel-based physiological (VBP) variables derived from blood oxygen level dependent (BOLD) fMRI time-course variations include: amplitude of low frequency fluctuations (ALFF), fractional amplitude of low frequency fluctuations (fALFF) and regional homogeneity (ReHo). Although these BOLD-derived variables can detect between-group (e.g. disease vs control) spatial pattern differences, physiological interpretations are not well established. The primary objective of this study was to quantify spatial correspondences between BOLD VBP variables and PET measurements of cerebral metabolic rate and hemodynamics, being well-validated physiological standards. To this end, quantitative, whole-brain PET images of metabolic rate of glucose (MRGlu; 18FDG) and oxygen (MRO2; 15OO), blood flow (BF; H215O) and blood volume (BV; C15O) were obtained in 16 healthy controls. In the same subjects, BOLD time-courses were obtained for computation of ALFF, fALFF and ReHo images. PET variables were compared pair-wise with BOLD variables. In group-averaged, across-region analyses, ALFF corresponded significantly only with BV (R = 0.64; p < 0.0001). fALFF corresponded most strongly with MRGlu (R = 0.79; p < 0.0001), but also significantly (p < 0.0001) with MRO2 (R = 0.68), BF (R = 0.68) and BV (R=0.68). ReHo performed similarly to fALFF, with significant strong correspondence (p < 0.0001) with MRGlu (R = 0.78), MRO2 (R = 0.54), and, but less strongly with BF (R = 0.50) and BV (R=0.50). Mutual information analyses further clarified these physiological interpretations. When conditioned by BV, ALFF retained no significant MRGlu, MRO2 or BF information. When conditioned by MRGlu, fALFF and ReHo retained no significant MRO2, BF or BV information. Of concern, however, the strength of PET-BOLD correspondences varied markedly by brain region, which calls for future investigation on physiological interpretations at a regional and per-subject basis.

The effect of acquisition duration on cerebral blood flow-based resting-state functional connectivity

Resting-state functional connectivity (rs-FC) is widely used to examine the functional architecture of the brain, and the blood-oxygenation-level-dependent (BOLD) signal is often utilized for determining rs-FC. However, the BOLD signal is susceptible to various factors that have less influence on the cerebral blood flow (CBF). Therefore, CBF could comprise an alternative for determining rs-FC. Since acquisition duration is one of the essential parameters for obtaining reliable rs-FC, we investigated the effect of acquisition duration on CBF-based rs-FC to examine the reliability of CBF-based rs-FC. Nineteen participants underwent CBF scanning for a total duration of 50 min. Variance of CBF-based rs-FC within the whole brain and 13 large-scale brain networks at various acquisition durations was compared to that with a 50-min duration using the Levene's test. Variance of CBF-based rs-FC at any durations did not differ from that at a 50-min duration (p > .05). Regarding variance of rs-FC within each large-scale brain network, the acquisition duration required to obtain reliable estimates of CBF-based rs-FC was shorter than 10 min and varied across large-scale brain networks. Altogether, an acquisition duration of at least 10 min is required to obtain reliable CBF-based rs-FC. These results indicate that CBF-based resting-state functional magnetic resonance imaging (rs-fMRI) with more than 10 min of total acquisition duration could be an alternative method to BOLD-based rs-fMRI to obtain reliable rs-FC.

Aberrant Resting-State Functional Brain Connectivity of Insular Subregions in Obstructive Sleep Apnea

The insular cortex is a cortical regulatory area involved in dyspnea, cognition, emotion, and sensorimotor function. Previous studies reported that obstructive sleep apnea (OSA) shows insular tissue damage and abnormal functional connections for the whole insula. The insula can be divided into different subregions with distinct functional profiles, including the ventral anterior insula (vAI) participating in affective processing, dorsal anterior insula (dAI) involved in cognitive processing, and posterior insula (PI) involved in the processing of sensorimotor information. However, the functional connectivity (FC) of these insular subregions in OSA has yet to be established. Hence, the purpose of this study was to explore the resting-state FC of the insular subregions with other brain areas and its relationship with clinical symptoms of OSA. Resting-state functional magnetic resonance imaging data from 83 male OSA patients and 84 healthy controls were analyzed by whole-brain voxel-based FC using spherical seeds from six insular subregions, namely, the bilateral vAI, dAI, and PI, to identify abnormalities in the insular subregions network and related brain regions. Ultimately, the Pearson correlation analysis was carried out between the dysfunction results and the neuropsychological tests. Compared with the healthy control group, the OSA patients exhibited disturbed FC from the dAI to areas relevant to cognition, such as the bilateral cerebellum posterior lobe, superior frontal gyrus, right middle frontal gyrus and middle temporal gyrus; decreased FC from the vAI to areas linked with emotion, such as the bilateral fusiform gyrus, superior parietal lobule, precuneus and cerebellum posterior lobe; and abnormal FC from the PI to the brain regions involved in sensorimotor such as the bilateral precentral gyrus, right superior/middle temporal gyrus and left superior frontal gyrus. The linear regression result showed that the apnea-hypopnea index was positively correlated with the increased FC between the right PI and the right precuneus (after Bonferroni correlation, P < 0.001) In conclusion, the abnormal FC between insular subregions and other brain regions were related to cognitive, emotional and sensorimotor networks in OSA patients. These results may provide a new imaging perspective for further understanding of OSA-related cognitive and affective disorders.

Advancing motion denoising of multiband resting-state functional connectivity fMRI data

Simultaneous multi-slice (multiband) accelerated functional magnetic resonance imaging (fMRI) provides dramatically improved temporal and spatial resolution for resting-state functional connectivity (RSFC) studies of the human brain in health and disease. However, multiband acceleration also poses unique challenges for denoising of subject motion induced data artifacts, the presence of which is a major confound in RSFC research that substantively diminishes reliability and reproducibility. We comprehensively evaluated existing and novel approaches to volume censoring-based motion denoising in the Human Connectome Project (HCP) dataset. We show that assumptions underlying common metrics for evaluating motion denoising pipelines, especially those based on quality control-functional connectivity (QC-FC) correlations and differences between high- and low-motion participants, are problematic, and appear to be inappropriate in their current widespread use as indicators of comparative pipeline performance and as targets for investigators to use when tuning pipelines for their own datasets. We further develop two new quantitative metrics that are instead agnostic to QC-FC correlations and other measures that rely upon the null assumption that no true relationships exist between trait measures of subject motion and functional connectivity, and demonstrate their use as benchmarks for comparing volume censoring methods. Finally, we develop and validate quantitative methods for determining dataset-specific optimal volume censoring parameters prior to the final analysis of a dataset, and provide straightforward recommendations and code for all investigators to apply this optimized approach to their own RSFC datasets.

Gradual Restraint Habituation for Awake Functional Magnetic Resonance Imaging Combined With a Sparse Imaging Paradigm Reduces Motion Artifacts and Stress Levels in Rodents

Functional magnetic resonance imaging, as a non-invasive technique, offers unique opportunities to assess brain function and connectivity under a broad range of applications, ranging from passive sensory stimulation to high-level cognitive abilities, in awake animals. This approach is confounded, however, by the fact that physical restraint and loud unpredictable acoustic noise must inevitably accompany fMRI recordings. These factors induce marked stress in rodents, and stress-related elevations of corticosterone levels are known to alter information processing and cognition in the rodent. Here, we propose a habituation strategy that spans specific stages of adaptation to restraint, MRI noise, and confinement stress in awake rats and circumvents the need for surgical head restraint. This habituation protocol results in stress levels during awake fMRI that do not differ from pre-handling levels and enables stable image acquisition with very low motion artifacts. For this, rats were gradually trained over a period of three weeks and eighteen training sessions. Stress levels were assessed by analysis of fecal corticosterone metabolite levels and breathing rates. We observed significant drops in stress levels to below pre-handling levels at the end of the habituation procedure. During fMRI in awake rats, after the conclusion of habituation and using a non-invasive head-fixation device, breathing was stable and head motion artifacts were minimal. A task-based fMRI experiment, using acoustic stimulation, conducted 2 days after the end of habituation, resulted in precise whole brain mapping of BOLD signals in the brain, with clear delineation of the expected auditory-related structures. The active discrimination by the animals of the acoustic stimuli from the backdrop of scanner noise was corroborated by significant increases in BOLD signals in the thalamus and reticular formation. Taken together, these data show that effective habituation to awake fMRI can be achieved by gradual and incremental acclimatization to the experimental conditions. Subsequent BOLD recordings, even during superimposed acoustic stimulation, reflect low stress-levels, low motion and a corresponding high-quality image acquisition. Furthermore, BOLD signals obtained during fMRI indicate that effective habituation facilitates selective attention to sensory stimuli that can in turn support the discrimination of cognitive processes in the absence of stress confounds.

Inter-electrode correlations measured with EEG predict individual differences in cognitive ability

Human brains share a broadly similar functional organization with consequential individual variation. This duality in brain function has primarily been observed when using techniques that consider the spatial organization of the brain, such as MRI. Here, we ask whether these common and unique signals of cognition are also present in temporally sensitive but spatially insensitive neural signals. To address this question, we compiled electroencephalogram (EEG) data from individuals of both sexes while they performed multiple working memory tasks at two different data-collection sites (n = 171 and 165). Results revealed that trial-averaged EEG activity exhibited inter-electrode correlations that were stable within individuals and unique across individuals. Furthermore, models based on these inter-electrode correlations generalized across datasets to predict participants' working memory capacity and general fluid intelligence. Thus, inter-electrode correlation patterns measured with EEG provide a signature of working memory and fluid intelligence in humans and a new framework for characterizing individual differences in cognitive abilities.

A multi-measure approach for assessing the performance of fMRI preprocessing strategies in resting-state functional connectivity

It is well established that head motion and physiological processes (e.g. cardiac and breathing activity) should be taken into consideration when analyzing and interpreting results in fMRI studies. However, even though recent studies aimed to evaluate the performance of different preprocessing pipelines there is still no consensus on the optimal strategy. This is partly due to the fact that the quality control (QC) metrics used to evaluate differences in performance across pipelines have often yielded contradictory results. Furthermore, preprocessing techniques based on physiological recordings or data decomposition techniques (e.g. aCompCor) have not been comprehensively examined. Here, to address the aforementioned issues, we propose a framework that summarizes the scores from eight previously proposed and novel QC metrics to a reduced set of two QC metrics that reflect the signal-to-noise ratio and the reduction in motion artifacts and biases in the preprocessed fMRI data. Using this framework, we evaluate the performance of three commonly used practices on the quality of data: 1) Removal of nuisance regressors from fMRI data, 2) discarding motion-contaminated volumes (i.e., scrubbing) before regression, and 3) low-pass filtering the data and the nuisance regressors before their removal. Using resting-state fMRI data from the Human Connectome Project, we show that the scores of the examined QC metrics improve the most when the global signal (GS) and about 17% of principal components from white matter (WM) are removed from the data. Finally, we observe a small further improvement with low-pass filtering at 0.20 Hz and milder variants of WM denoising, but not with scrubbing.

Physiological noise modeling in fMRI based on the pulsatile component of photoplethysmograph

The blood oxygenation level-dependent (BOLD) contrast mechanism allows the noninvasive monitoring of changes in deoxyhemoglobin content. As such, it is commonly used in functional magnetic resonance imaging (fMRI) to study brain activity since levels of deoxyhemoglobin are indirectly related to local neuronal activity through neurovascular coupling mechanisms. However, the BOLD signal is severely affected by physiological processes as well as motion. Due to this, several noise correction techniques have been developed to correct for the associated confounds. The present study focuses on cardiac pulsatility fMRI confounds, aiming to refine model-based techniques that utilize the photoplethysmograph (PPG) signal. Specifically, we propose a new technique based on convolution filtering, termed cardiac pulsatility model (CPM) and compare its performance with the cardiac-related RETROICOR (Card-RETROICOR), which is a technique commonly used to model fMRI fluctuations due to cardiac pulsatility. Further, we investigate whether variations in the amplitude of the PPG pulses (PPG-Amp) covary with variations in amplitude of pulse-related fMRI fluctuations, as well as with the systemic low frequency oscillations (SLFOs) component of the fMRI global signal (GS - defined as the mean signal across all gray matter voxels). Capitalizing on 3T fMRI data from the Human Connectome Project, CPM was found to explain a significantly larger fraction of the fMRI signal variance compared to Card-RETROICOR, particularly for subjects with larger heart rate variability during the scan. The amplitude of the fMRI pulse-related fluctuations did not covary with PPG-Amp; however, PPG-Amp explained significant variance in the GS that was not attributed to variations in heart rate or breathing patterns. Our results suggest that the proposed approach can model high-frequency fluctuations due to pulsation as well as low-frequency physiological fluctuations more accurately compared to model-based techniques commonly employed in fMRI studies.

Physiological and motion signatures in static and time-varying functional connectivity and their subject identifiability

Human brain connectivity yields significant potential as a noninvasive biomarker. Several studies have used fMRI-based connectivity fingerprinting to characterize individual patterns of brain activity. However, it is not clear whether these patterns mainly reflect neural activity or the effect of physiological and motion processes. To answer this question, we capitalize on a large data sample from the Human Connectome Project and rigorously investigate the contribution of the aforementioned processes on functional connectivity (FC) and time-varying FC, as well as their contribution to subject identifiability. We find that head motion, as well as heart rate and breathing fluctuations, induce artifactual connectivity within distinct resting-state networks and that they correlate with recurrent patterns in time-varying FC. Even though the spatiotemporal signatures of these processes yield above-chance levels in subject identifiability, removing their effects at the preprocessing stage improves identifiability, suggesting a neural component underpinning the inter-individual differences in connectivity.