Identifying and characterizing systematic temporally-lagged BOLD artifacts

Cognitive Dynamics Estimation: A whole-brain spatial regression paradigm for extracting the temporal dynamics of cognitive processes

Functional MRI (fMRI) has been instrumental in understanding how cognitive processes are spatially mapped in the brain, yielding insights about brain regions and functions. However, in case the orthogonality of behavioral or stimulus timing is not guaranteed, the estimated brain maps fail to dissociate each cognitive process, and the resultant maps become unstable. Also, the brain mapping exercise can not provide temporal information on the cognitive process. Here we propose a qualitatively different approach to fMRI analysis, named Cognitive Dynamics Estimation (CDE), that estimates how multiple cognitive processes change over time even when behavior or stimulus logs are unavailable. This method transposes the conventional brain mapping; the brain activity pattern at each time point is subject to regression analysis with data-driven maps of cognitive processes as regressors, resulting in the time series of cognitive processes. The estimated time series captured the fluctuation of intensity and timing of cognitive processes on a trial-by-trial basis, which conventional analysis could not capture. Notably, the estimated time series predicted participants' cognitive ability to perform each psychological task. As an addition to our fMRI analytic toolkit, these results suggest the potential for CDE to elucidate underexplored cognitive phenomena, especially in the temporal domain.

Edge-centric analysis of time-varying functional brain networks with applications in autism spectrum disorder

The interaction between brain regions changes over time, which can be characterized using time-varying functional connectivity (tvFC). The common approach to estimate tvFC uses sliding windows and offers limited temporal resolution. An alternative method is to use the recently proposed edge-centric approach, which enables the tracking of moment-to-moment changes in co-fluctuation patterns between pairs of brain regions. Here, we first examined the dynamic features of edge time series and compared them to those in the sliding window tvFC (sw-tvFC). Then, we used edge time series to compare subjects with autism spectrum disorder (ASD) and healthy controls (CN). Our results indicate that relative to sw-tvFC, edge time series captured rapid and bursty network-level fluctuations that synchronize across subjects during movie-watching. The results from the second part of the study suggested that the magnitude of peak amplitude in the collective co-fluctuations of brain regions (estimated as root sum square (RSS) of edge time series) is similar in CN and ASD. However, the trough-to-trough duration in RSS signal is greater in ASD, compared to CN. Furthermore, an edge-wise comparison of high-amplitude co-fluctuations showed that the within-network edges exhibited greater magnitude fluctuations in CN. Our findings suggest that high-amplitude co-fluctuations captured by edge time series provide details about the disruption of functional brain dynamics that could potentially be used in developing new biomarkers of mental disorders.

Decreased resting-state functional connectivity of the habenula-cerebellar in a major depressive disorder

BACKGROUND

In animal experiments, the habenula and septal nuclei are known as the key brain areas of depression. However, there are few magnetic resonance imaging (MRI) studies on the functional connectivity between these areas and the subcortical areas in humans with major depression. We aimed to investigate the difference in resting-state functional connectivity (RSFC) among the major regions of interest (ROI) in the subcortical areas, including both the habenula and septal nuclei.

METHODS

We performed the seed-to-voxel analysis to investigate the RSFC between both the habenula and septal nucleus, as well as other subcortical regions. Furthermore, ROI-to-ROI analysis was performed among the combinations of ROI pairs in the subcortical areas.

RESULTS

The seed-to-voxel analysis showed a lower RSFC between the left habenula and the cerebellum in major depressive disorder (MDD) than in healthy controls (HCs). As a result of ROI-to-ROI analysis in subcortical areas, a total of 31 pairs of FCs in the MDD group showed a lower RSFC than in the HCs group.

CONCLUSION

This study revealed a lower RSFC between the left habenula and cerebellum in patients with MDD and reduced RSFC among numerous subcortical areas. These new findings on the neural circuitry of MDD might contribute to an in-depth understanding of depression.

Multi-Band Brain Network Analysis for Functional Neuroimaging Biomarker Identification

The functional connectomic profile is one of the non-invasive imaging biomarkers in the computer-assisted diagnostic system for many neuro-diseases. However, the diagnostic power of functional connectivity is challenged by mixed frequency-specific neuronal oscillations in the brain, which makes the single Functional Connectivity Network (FCN) often underpowered to capture the disease-related functional patterns. To address this challenge, we propose a novel functional connectivity analysis framework to conduct joint feature learning and personalized disease diagnosis, in a semi-supervised manner, aiming at focusing on putative multi-band functional connectivity biomarkers from functional neuroimaging data. Specifically, we first decompose the Blood Oxygenation Level Dependent (BOLD) signals into multiple frequency bands by the discrete wavelet transform, and then cast the alignment of all fully-connected FCNs derived from multiple frequency bands into a parameter-free multi-band fusion model. The proposed fusion model fuses all fully-connected FCNs to obtain a sparsely-connected FCN (sparse FCN for short) for each individual subject, as well as lets each sparse FCN be close to its neighbored sparse FCNs and be far away from its furthest sparse FCNs. Furthermore, we employ the l1 -SVM to conduct joint brain region selection and disease diagnosis. Finally, we evaluate the effectiveness of our proposed framework on various neuro-diseases, i.e., Fronto-Temporal Dementia (FTD), Obsessive-Compulsive Disorder (OCD), and Alzheimer's Disease (AD), and the experimental results demonstrate that our framework shows more reasonable results, compared to state-of-the-art methods, in terms of classification performance and the selected brain regions. The source code can be visited by the url https://github.com/reynard-hu/mbbna.

You read my mind: fMRI markers of threatening appraisals in people with persistent psychotic experiences

Anomalous perceptual experiences are relatively common in the general population. Evidence indicates that the key to distinguishing individuals with persistent psychotic experiences (PEs) with a need for care from those without is how they appraise their anomalous experiences. Here, we aimed to characterise the neural circuits underlying threatening and non-threatening appraisals in people with and without a need for care for PEs, respectively. A total of 48 participants, consisting of patients with psychosis spectrum disorder (clinical group, n = 16), non-need-for-care participants with PEs (non-clinical group, n = 16), and no-PE healthy control participants (n = 16), underwent functional magnetic resonance imaging while completing the Telepath task, designed to induce an anomalous perceptual experience. Appraisals of the anomalous perceptual experiences were examined, as well as functional brain responses during this window, for significant group differences. We also examined whether activation co-varied with the subjective threat appraisals reported in-task by participants. The clinical group reported elevated subjective threat appraisals compared to both the non-clinical and no-PE control groups, with no differences between the two non-clinical groups. This pattern of results was accompanied by reduced activation in the superior and inferior frontal gyri in the clinical group as compared to the non-clinical and control groups. Precuneus activation scaled with threat appraisals reported in-task. Resilience in the context of persistent anomalous experiences may be explained by intact functioning of fronto-parietal regions, and may correspond to the ability to contextualise and flexibly evaluate psychotic experiences.

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.

Co-activation pattern alterations in autism spectrum disorder-A volume-wise hierarchical clustering fMRI study

INTRODUCTION

There has been a growing effort to characterize the time-varying functional connectivity of resting state (RS) fMRI brain networks (RSNs). Although voxel-wise connectivity studies have examined different sliding window lengths, nonsequential volume-wise approaches have been less common.

METHODS

Inspired by earlier co-activation pattern (CAP) studies, we applied hierarchical clustering (HC) to classify the image volumes of the RS-fMRI data on 28 adolescents with autism spectrum disorder (ASD) and their 27 typically developing (TD) controls. We compared the distribution of the ASD and TD groups' volumes in CAPs as well as their voxel-wise means. For simplification purposes, we conducted a group independent component analysis to extract 14 major RSNs. The RSNs' average z-scores enabled us to meaningfully regroup the RSNs and estimate the percentage of voxels within each RSN for which there was a significant group difference. These results were jointly interpreted to find global group-specific patterns.

RESULTS

We found similar brain state proportions in 58 CAPs (clustering interval from 2 to 30). However, in many CAPs, the voxel-wise means differed significantly within a matrix of 14 RSNs. The rest-activated default mode-positive and default mode-negative brain state properties vary considerably in both groups over time. This division was seen clearly when the volumes were partitioned into two CAPs and then further examined along the HC dendrogram of the diversifying brain CAPs. The ASD group network activations followed a more heterogeneous distribution and some networks maintained higher baselines; throughout the brain deactivation state, the ASD participants had reduced deactivation in 12/14 networks. During default mode-negative CAPs, the ASD group showed simultaneous visual network and either dorsal attention or default mode network overactivation.

CONCLUSION

Nonsequential volume gathering into CAPs and the comparison of voxel-wise signal changes provide a complementary perspective to connectivity and an alternative to sliding window analysis.

Functional MRI of arousals in nonrapid eye movement sleep

Arousals commonly occur during human sleep and have been associated with several sleep disorders. Arousals are characterized as an abrupt electroencephalography (EEG) frequency change to higher frequencies during sleep. However, the human brain regions involved in arousal are not yet clear. Simultaneous EEG and functional magnetic resonance imaging (fMRI) data were recorded during the early portion of the sleep period in healthy young adults. Arousals were identified based on the EEG data, and fMRI signal changes associated with 83 arousals from 19 subjects were analyzed. Subcortical regions, including the midbrain, thalamus, basal ganglia, and cerebellum, were activated with arousal. Cortices, including the temporal gyrus, occipital gyrus, and frontal gyrus, were deactivated with arousal. The activations associated with arousal in the subcortical regions were consistent with previous findings of subcortical involvement in behavioral arousal and consciousness. Cortical deactivations may serve as a mechanism to direct incoming sensory stimuli to specific brain regions, thereby monitoring environmental perturbations during sleep.

Associations between Neighborhood SES and Functional Brain Network Development

Higher socioeconomic status (SES) in childhood is associated with stronger cognitive abilities, higher academic achievement, and lower incidence of mental illness later in development. While prior work has mapped the associations between neighborhood SES and brain structure, little is known about the relationship between SES and intrinsic neural dynamics. Here, we capitalize upon a large cross-sectional community-based sample (Philadelphia Neurodevelopmental Cohort, ages 8-22 years, n = 1012) to examine associations between age, SES, and functional brain network topology. We characterize this topology using a local measure of network segregation known as the clustering coefficient and find that it accounts for a greater degree of SES-associated variance than mesoscale segregation captured by modularity. High-SES youth displayed stronger positive associations between age and clustering than low-SES youth, and this effect was most pronounced for regions in the limbic, somatomotor, and ventral attention systems. The moderating effect of SES on positive associations between age and clustering was strongest for connections of intermediate length and was consistent with a stronger negative relationship between age and local connectivity in these regions in low-SES youth. Our findings suggest that, in late childhood and adolescence, neighborhood SES is associated with variation in the development of functional network structure in the human brain.

A Hilbert-based method for processing respiratory timeseries

In this technical note, we introduce a new method for estimating changes in respiratory volume per unit time (RVT) from respiratory bellows recordings. By using techniques from the electrophysiological literature, in particular the Hilbert transform, we show how we can better characterise breathing rhythms, with the goal of improving physiological noise correction in functional magnetic resonance imaging (fMRI). Specifically, our approach leads to a representation with higher time resolution and better captures atypical breathing events than current peak-based RVT estimators. Finally, we demonstrate that this leads to an increase in the amount of respiration-related variance removed from fMRI data when used as part of a typical preprocessing pipeline. Our implementation is publicly available as part of the PhysIO package, which is distributed as part of the open-source TAPAS toolbox (https://translationalneuromodeling.org/tapas).

Dynamic time warping outperforms Pearson correlation in detecting atypical functional connectivity in autism spectrum disorders

Resting state fMRI (rsfMRI) is frequently used to study brain function, including in clinical populations. Similarity of blood-oxygen-level-dependent (BOLD) fluctuations during rsfMRI between brain regions is thought to reflect intrinsic functional connectivity (FC), potentially due to history of coactivation. To quantify similarity, studies have almost exclusively relied on Pearson correlation, which assumes linearity and can therefore underestimate FC if the hemodynamic response function differs regionally or there is BOLD signal lag between timeseries. Here we show in three cohorts of children, adolescents and adults, with and without autism spectrum disorders (ASDs), that measuring the similarity of BOLD signal fluctuations using non-linear dynamic time warping (DTW) is more robust to global signal regression (GSR), has higher test-retest reliability and is more sensitive to task-related changes in FC. Additionally, when comparing FC between individuals with ASDs and typical controls, more group differences are detected using DTW. DTW estimates are also more related to ASD symptom severity and executive function, while Pearson correlation estimates of FC are more strongly associated with respiration during rsfMRI. Together these findings suggest that non-linear methods such as DTW improve estimation of resting state FC, particularly when studying clinical populations whose hemodynamics or neurovascular coupling may be altered compared to typical controls.

Dissimilarity of functional connectivity uncovers the influence of participant's motion in functional magnetic resonance imaging studies

Head motion is a major confounding factor impairing the quality of functional magnetic resonance imaging (fMRI) data. In particular, head motion can reduce analytical efficiency, and its effects are still present even after preprocessing. To examine the validity of motion removal and to evaluate the remaining effects of motion on the quality of the preprocessed fMRI data, a new metric of group quality control (QC), dissimilarity of functional connectivity, is introduced. Here, we investigate the association between head motion, represented by mean framewise displacement, and dissimilarity of functional connectivity by applying four preprocessing methods in two independent resting‐state fMRI datasets: one consisting of healthy participants (N = 167) scanned in a 3T GE‐Discovery 750 with longer TR (2.5 s), and the other of chronic back pain patients (N = 143) in a 3T Siemens Magnetom Prisma scanner with shorter TR (0.555 s). We found that dissimilarity of functional connectivity uncovers the influence of participant's motion, and this relationship is independent of population, scanner, and preprocessing method. The association between motion and dissimilarity of functional connectivity, and how the removal of high‐motion participants affects this association, is a new strategy for group‐level QC following preprocessing.

Accurate prediction of individual subject identity and task, but not autism diagnosis, from functional connectomes

Despite enthusiasm about the potential for using fMRI-based functional connectomes in the development of biomarkers for autism spectrum disorder (ASD), the literature is full of negative findings-failures to distinguish ASD functional connectomes from those of typically developing controls (TD)-and positive findings that are inconsistent across studies. Here, we report on a new study designed to either better differentiate ASD from TD functional connectomes-or, alternatively, to refine our understanding of the factors underlying the current state of affairs. We scanned individuals with ASD and controls both at rest and while watching videos with social content. Using multiband fMRI across repeat sessions, we improved both data quantity and scanning duration by collecting up to 2 hr of data per individual. This is about 50 times the typical number of temporal samples per individual in ASD fcMRI studies. We obtained functional connectomes that were discriminable, allowing for near-perfect individual identification regardless of diagnosis, and equally reliable in both groups. However, contrary to what one might expect, we did not consistently or robustly observe in the ASD group either reductions in similarity to TD functional connectivity (FC) patterns or shared atypical FC patterns. Accordingly, FC-based predictions of diagnosis group achieved accuracy levels around chance. However, using the same approaches to predict scan type (rest vs. video) achieved near-perfect accuracy. Our findings suggest that neither the limitations of resting state as a "task," data resolution, data quantity, or scan duration can be considered solely responsible for failures to differentiate ASD from TD functional connectomes.

Combining Physiological and Neuroimaging Measures to Predict Affect Processing Induced by Affectively Valent Image Stimuli

The importance of affect processing to human behavior has long driven researchers to pursue its measurement. In this study, we compared the relative fidelity of measurements of neural activation and physiology (i.e., heart rate change) in detecting affective valence induction across a broad continuum of conveyed affective valence. We combined intra-subject neural activation based multivariate predictions of affective valence with measures of heart rate (HR) deceleration to predict predefined normative affect rating scores for stimuli drawn from the International Affective Picture System (IAPS) in a population (n = 50) of healthy adults. In sum, we found that patterns of neural activation and HR deceleration significantly, and uniquely, explain the variance in normative valent scores associated with IAPS stimuli; however, we also found that patterns of neural activation explain a significantly greater proportion of that variance. These traits persisted across a range of stimulus sets, differing by the polar-extremity of their positively and negatively valent subsets, which represent the positively and negatively valent polar-extremity of stimulus sets reported in the literature. Overall, these findings support the acquisition of heart rate deceleration concurrently with fMRI to provide convergent validation of induced affect processing in the dimension of affective valence.

Temporal fluctuations in the brain's modular architecture during movie-watching

Brain networks are flexible and reconfigure over time to support ongoing cognitive processes. However, tracking statistically meaningful reconfigurations across time has proven difficult. This has to do largely with issues related to sampling variability, making instantaneous estimation of network organization difficult, along with increased reliance on task-free (cognitively unconstrained) experimental paradigms, limiting the ability to interpret the origin of changes in network structure over time. Here, we address these challenges using time-varying network analysis in conjunction with a naturalistic viewing paradigm. Specifically, we developed a measure of inter-subject network similarity and used this measure as a coincidence filter to identify synchronous fluctuations in network organization across individuals. Applied to movie-watching data, we found that periods of high inter-subject similarity coincided with reductions in network modularity and increased connectivity between cognitive systems. In contrast, low inter-subject similarity was associated with increased system segregation and more rest-like architectures. We then used a data-driven approach to uncover clusters of functional connections that follow similar trajectories over time and are more strongly correlated during movie-watching than at rest. Finally, we show that synchronous fluctuations in network architecture over time can be linked to a subset of features in the movie. Our findings link dynamic fluctuations in network integration and segregation to patterns of inter-subject similarity, and suggest that moment-to-moment fluctuations in functional connectivity reflect shared cognitive processing across individuals.

Transient Arousal Modulations Contribute to Resting-State Functional Connectivity Changes Associated with Head Motion Parameters

Correlations of resting-state functional magnetic resonance imaging (rsfMRI) signals are being widely used for assessing the functional brain connectivity in health and disease. However, an association was recently observed between rsfMRI connectivity modulations and the head motion parameters and regarded as a causal relationship, which has raised serious concerns about the validity of many rsfMRI findings. Here, we studied the origin of this rsfMRI-motion association and its relationship to arousal modulations. By using a template-matching method to locate arousal-related fMRI changes, we showed that the effects of high motion time points on rsfMRI connectivity are largely due to their significant overlap with arousal-affected time points. The finding suggests that the association between rsfMRI connectivity and the head motion parameters arises from their comodulations at transient arousal modulations, and this information is critical not only for proper interpretation of motion-associated rsfMRI connectivity changes, but also for controlling the potential confounding effects of arousal modulation on rsfMRI metrics.

Disambiguating the role of blood flow and global signal with partial information decomposition

Global signal (GS) is an ubiquitous construct in resting state functional magnetic resonance imaging (rs-fMRI), associated to nuisance, but containing by definition most of the neuronal signal. Global signal regression (GSR) effectively removes the impact of physiological noise and other artifacts, but at the same time it alters correlational patterns in unpredicted ways. Performing GSR taking into account the underlying physiology (mainly the blood arrival time) has been proven to be beneficial. From these observations we aimed to: 1) characterize the effect of GSR on network-level functional connectivity in a large dataset; 2) assess the complementary role of global signal and vessels; and 3) use the framework of partial information decomposition to further look into the joint dynamics of the global signal and vessels, and their respective influence on the dynamics of cortical areas. We observe that GSR affects intrinsic connectivity networks in the connectome in a non-uniform way. Furthermore, by estimating the predictive information of blood flow and the global signal using partial information decomposition, we observe that both signals are present in different amounts across intrinsic connectivity networks. Simulations showed that differences in blood arrival time can largely explain this phenomenon, while using hemodynamic and calcium mouse recordings we were able to confirm the presence of vascular effects, as calcium recordings lack hemodynamic information. With these results we confirm network-specific effects of GSR and the importance of taking blood flow into account for improving de-noising methods. Additionally, and beyond the mere issue of data denoising, we quantify the diverse and complementary effect of global and vessel BOLD signals on the dynamics of cortical areas.

Resting-state "physiological networks"

Slow changes in systemic brain physiology can elicit large fluctuations in fMRI time series, which manifest as structured spatial patterns of temporal correlations between distant brain regions. Here, we investigated whether such "physiological networks"-sets of segregated brain regions that exhibit similar responses following slow changes in systemic physiology-resemble patterns associated with large-scale networks typically attributed to remotely synchronized neuronal activity. By analyzing a large group of subjects from the 3T Human Connectome Project (HCP) database, we demonstrate brain-wide and noticeably heterogenous dynamics tightly coupled to either respiratory variation or heart rate changes. We show, using synthesized data generated from physiological recordings across subjects, that these physiologically-coupled fluctuations alone can produce networks that strongly resemble previously reported resting-state networks, suggesting that, in some cases, the "physiological networks" seem to mimic the neuronal networks. Further, we show that such physiologically-relevant connectivity estimates appear to dominate the overall connectivity observations in multiple HCP subjects, and that this apparent "physiological connectivity" cannot be removed by the use of a single nuisance regressor for the entire brain (such as global signal regression) due to the clear regional heterogeneity of the physiologically-coupled responses. Our results challenge previous notions that physiological confounds are either localized to large veins or globally coherent across the cortex, therefore emphasizing the necessity to consider potential physiological contributions in fMRI-based functional connectivity studies. The rich spatiotemporal patterns carried by such "physiological" dynamics also suggest great potential for clinical biomarkers that are complementary to large-scale neuronal networks.

State and trait characteristics of anterior insula time-varying functional connectivity

The human anterior insula (aINS) is a topographically organized brain region, in which ventral portions contribute to socio-emotional function through limbic and autonomic connections, whereas the dorsal aINS contributes to cognitive processes through frontal and parietal connections. Open questions remain, however, regarding how aINS connectivity varies over time. We implemented a novel approach combining seed-to-whole-brain sliding-window functional connectivity MRI and k-means clustering to assess time-varying functional connectivity of aINS subregions. We studied three independent large samples of healthy participants and longitudinal datasets to assess inter- and intra-subject stability, and related aINS time-varying functional connectivity profiles to dispositional empathy. We identified four robust aINS time-varying functional connectivity modes that displayed both "state" and "trait" characteristics: while modes featuring connectivity to sensory regions were modulated by eye closure, modes featuring connectivity to higher cognitive and emotional processing regions were stable over time and related to empathy measures.

Nonreplication of functional connectivity differences in autism spectrum disorder across multiple sites and denoising strategies

A rapidly growing number of studies on autism spectrum disorder (ASD) have used resting‐state fMRI to identify alterations of functional connectivity, with the hope of identifying clinical biomarkers or underlying neural mechanisms. However, results have been largely inconsistent across studies, and there remains a pressing need to determine the primary factors influencing replicability. Here, we used resting‐state fMRI data from the Autism Brain Imaging Data Exchange to investigate two potential factors: denoising strategy and data site (which differ in terms of sample, data acquisition, etc.). We examined the similarity of both group‐averaged functional connectomes and group‐level differences (ASD vs. control) across 33 denoising pipelines and four independently‐acquired datasets. The group‐averaged connectomes were highly consistent across pipelines (r = 0.92 ± 0.06) and sites (r = 0.88 ± 0.02). However, the group differences, while still consistent within site across pipelines (r = 0.76 ± 0.12), were highly inconsistent across sites regardless of choice of denoising strategies (r = 0.07 ± 0.04), suggesting lack of replication may be strongly influenced by site and/or cohort differences. Across‐site similarity remained low even when considering the data at a large‐scale network level or when considering only the most significant edges. We further show through an extensive literature survey that the parameters chosen in the current study (i.e., sample size, age range, preprocessing methods) are quite representative of the published literature. These results highlight the importance of examining replicability in future studies of ASD, and, more generally, call for extra caution when interpreting alterations in functional connectivity across groups of individuals.

Arousal Contributions to Resting-State fMRI Connectivity and Dynamics

Resting-state functional magnetic resonance imaging (rsfMRI) is being widely used for charting brain connectivity and dynamics in healthy and diseased brains. However, the resting state paradigm allows an unconstrained fluctuation of brain arousal, which may have profound effects on resting-state fMRI signals and associated connectivity/dynamic metrics. Here, we review current understandings of the relationship between resting-state fMRI and brain arousal, in particular the effect of a recently discovered event of arousal modulation on resting-state fMRI. We further discuss potential implications of arousal-related fMRI modulation with a focus on its potential role in mediating spurious correlations between resting-state connectivity/dynamics with physiology and behavior. Multiple hypotheses are formulated based on existing evidence and remain to be tested by future studies.

Characteristics of respiratory measures in young adults scanned at rest, including systematic changes and "missed" deep breaths

Breathing rate and depth influence the concentration of carbon dioxide in the blood, altering cerebral blood flow and thus functional magnetic resonance imaging (fMRI) signals. Such respiratory fluctuations can have substantial influence in studies of fMRI signal covariance in subjects at rest, the so-called "resting state functional connectivity" technique. If respiration is monitored during fMRI scanning, it is typically done using a belt about the subject's abdomen to record abdominal circumference. Several measures have been derived from these belt records, including the windowed envelope of the waveform (ENV), the windowed variance in the waveform (respiration variation, RV), and a measure of the amplitude of each breath divided by the cycle time of the breath (respiration volume per time, RVT). Any attempt to gauge respiratory contributions to fMRI signals requires a respiratory measure, but little is known about how these measures compare to each other, or how they perform beyond the small studies in which they were initially proposed. Here, we examine the properties of these measures in hundreds of healthy young adults scanned for an hour each at rest, a subset of the Human Connectome Project chosen for having high-quality physiological records. We find: 1) ENV, RV, and RVT are all correlated, and ENV and RV are more highly correlated to each other than to RVT; 2) respiratory events like deep breaths exhibit characteristic heart rate elevations, fMRI signal changes, head motions, and image quality abnormalities time-locked to large deflections in the belt traces; 3) all measures can "miss" deep breaths; 4) RVT "misses" deep breaths more than ENV or RV; 5) all respiratory measures change systematically over the course of a 14.4-min scan. We discuss the implications of these findings for the literature and ways to move forward in modeling respiratory influences on fMRI scans.

Axial variation of deoxyhemoglobin density as a source of the low-frequency time lag structure in blood oxygenation level-dependent signals

Perfusion-related information is reportedly embedded in the low-frequency component of a blood oxygen level-dependent (BOLD) functional magnetic resonance imaging (fMRI) signal. The blood-propagation pattern through the cerebral vascular tree is detected as an interregional lag variation of spontaneous low-frequency oscillations (sLFOs). Mapping of this lag, or phase, has been implicitly treated as a projection of the vascular tree structure onto real space. While accumulating evidence supports the biological significance of this signal component, the physiological basis of the “perfusion lag structure,” a requirement for an integrative resting-state fMRI-signal model, is lacking. In this study, we conducted analyses furthering the hypothesis that the sLFO is not only largely of systemic origin, but also essentially intrinsic to blood, and hence behaves as a virtual tracer. By summing the small fluctuations of instantaneous phase differences between adjacent vascular regions, a velocity response to respiratory challenges was detected. Regarding the relationship to neurovascular coupling, the removal of the whole lag structure, which can be considered as an optimized global-signal regression, resulted in a reduction of inter-individual variance while preserving the fMRI response. Examination of the T2* and S₀, or non-BOLD, components of the fMRI signal revealed that the lag structure is deoxyhemoglobin dependent, while paradoxically presenting a signal-magnitude reduction in the venous side of the cerebral vasculature. These findings provide insight into the origin of BOLD sLFOs, suggesting that they are highly intrinsic to the circulating blood.

Identification of physiological response functions to correct for fluctuations in resting-state fMRI related to heart rate and respiration

Functional magnetic resonance imaging (fMRI) is widely viewed as the gold standard for studying brain function due to its high spatial resolution and non-invasive nature. However, it is well established that changes in breathing patterns and heart rate strongly influence the blood oxygen-level dependent (BOLD) fMRI signal and this, in turn, can have considerable effects on fMRI studies, particularly resting-state studies. The dynamic effects of physiological processes are often quantified by using convolution models along with simultaneously recorded physiological data. In this context, physiological response function (PRF) curves (cardiac and respiratory response functions), which are convolved with the corresponding physiological fluctuations, are commonly employed. While it has often been suggested that the PRF curves may be region- or subject-specific, it is still an open question whether this is the case. In the present study, we propose a novel framework for the robust estimation of PRF curves and use this framework to rigorously examine the implications of using population-, subject-, session- and scan-specific PRF curves. The proposed framework was tested on resting-state fMRI and physiological data from the Human Connectome Project. Our results suggest that PRF curves vary significantly across subjects and, to a lesser extent, across sessions from the same subject. These differences can be partly attributed to physiological variables such as the mean and variance of the heart rate during the scan. The proposed methodological framework can be used to obtain robust scan-specific PRF curves from data records with duration longer than 5 min, exhibiting significantly improved performance compared to previously defined canonical cardiac and respiration response functions. Besides removing physiological confounds from the BOLD signal, accurate modeling of subject- (or session-/scan-) specific PRF curves is of importance in studies that involve populations with altered vascular responses, such as aging subjects.

Global signal regression strengthens association between resting-state functional connectivity and behavior

Global signal regression (GSR) is one of the most debated preprocessing strategies for resting-state functional MRI. GSR effectively removes global artifacts driven by motion and respiration, but also discards globally distributed neural information and introduces negative correlations between certain brain regions. The vast majority of previous studies have focused on the effectiveness of GSR in removing imaging artifacts, as well as its potential biases. Given the growing interest in functional connectivity fingerprinting, here we considered the utilitarian question of whether GSR strengthens or weakens associations between resting-state functional connectivity (RSFC) and multiple behavioral measures across cognition, personality and emotion. By applying the variance component model to the Brain Genomics Superstruct Project (GSP), we found that behavioral variance explained by whole-brain RSFC increased by an average of 47% across 23 behavioral measures after GSR. In the Human Connectome Project (HCP), we found that behavioral variance explained by whole-brain RSFC increased by an average of 40% across 58 behavioral measures, when GSR was applied after ICA-FIX de-noising. To ensure generalizability, we repeated our analyses using kernel regression. GSR improved behavioral prediction accuracies by an average of 64% and 12% in the GSP and HCP datasets respectively. Importantly, the results were consistent across methods. A behavioral measure with greater RSFC-explained variance (using the variance component model) also exhibited greater prediction accuracy (using kernel regression). A behavioral measure with greater improvement in behavioral variance explained after GSR (using the variance component model) also enjoyed greater improvement in prediction accuracy after GSR (using kernel regression). Furthermore, GSR appeared to benefit task performance measures more than self-reported measures. Since GSR was more effective at removing motion-related and respiratory-related artifacts, GSR-related increases in variance explained and prediction accuracies were unlikely the result of motion-related or respiratory-related artifacts. However, it is worth emphasizing that the current study focused on whole-brain RSFC, so it remains unclear whether GSR improves RSFC-behavioral associations for specific connections or networks. Overall, our results suggest that at least in the case for young healthy adults, GSR strengthens the associations between RSFC and most (although not all) behavioral measures. Code for the variance component model and ridge regression can be found here: https://github.com/ThomasYeoLab/CBIG/tree/master/stable_projects/preprocessing/Li2019_GSR.

On time delay estimation and sampling error in resting-state fMRI

Accumulating evidence indicates that resting-state functional magnetic resonance imaging (rsfMRI) signals correspond to propagating electrophysiological infra-slow activity (<0.1 Hz). Thus, pairwise correlations (zero-lag functional connectivity (FC)) and temporal delays among regional rsfMRI signals provide useful, complementary descriptions of spatiotemporal structure in infra-slow activity. However, the slow nature of fMRI signals implies that practical scan durations cannot provide sufficient independent temporal samples to stabilize either of these measures. Here, we examine factors affecting sampling variability in both time delay estimation (TDE) and FC. Although both TDE and FC accuracy are highly sensitive to data quantity, we use surrogate fMRI time series to study how the former is additionally related to the magnitude of a given pairwise correlation and, to a lesser extent, the temporal sampling rate. These contingencies are further explored in real data comprising 30-min rsfMRI scans, where sampling error (i.e., limited accuracy owing to insufficient data quantity) emerges as a significant but underappreciated challenge to FC and, even more so, to TDE. Exclusion of high-motion epochs exacerbates sampling error; thus, both sides of the bias-variance (or data quality-quantity) tradeoff associated with data exclusion should be considered when analyzing rsfMRI data. Finally, we present strategies for TDE in motion-corrupted data, for characterizing sampling error in TDE and FC, and for mitigating the influence of sampling error on lag-based analyses.

High-accuracy individual identification using a "thin slice" of the functional connectome

Connectome fingerprinting-a method that uses many thousands of functional connections in aggregate to identify individuals-holds promise for individualized neuroimaging. A better characterization of the features underlying successful fingerprinting performance-how many and which functional connections are necessary and/or sufficient for high accuracy-will further inform our understanding of uniqueness in brain functioning. Thus, here we examine the limits of high-accuracy individual identification from functional connectomes. Using ∼3,300 scans from the Human Connectome Project in a split-half design and an independent replication sample, we find that a remarkably small "thin slice" of the connectome-as few as 40 out of 64,620 functional connections-was sufficient to uniquely identify individuals. Yet, we find that no specific connections or even specific networks were necessary for identification, as even small random samples of the connectome were sufficient. These results have important conceptual and practical implications for the manifestation and detection of uniqueness in the brain.

Functional connectivity-based subtypes of individuals with and without autism spectrum disorder

Autism spectrum disorder (ASD) is a heterogeneous neurodevelopmental disorder, characterized by impairments in social communication and restricted, repetitive behaviors. Neuroimaging studies have shown complex patterns and functional connectivity (FC) in ASD, with no clear consensus on brain-behavior relationships or shared patterns of FC with typically developing controls. Here, we used a dimensional approach to characterize two distinct clusters of FC patterns across both ASD participants and controls using k-means clustering. Using multivariate statistical analyses, a categorical approach was taken to characterize differences in FC between subtypes and between diagnostic groups. One subtype was defined by increased FC within resting-state networks and decreased FC across networks compared with the other subtype. A separate FC pattern distinguished ASD from controls, particularly within default mode, cingulo-opercular, sensorimotor, and occipital networks. There was no significant interaction between subtypes and diagnostic groups. Finally, a dimensional analysis of FC patterns with behavioral measures of IQ, social responsiveness, and ASD severity showed unique brain-behavior relations in each subtype and a continuum of brain-behavior relations from ASD to controls within one subtype. These results demonstrate that distinct clusters of FC patterns exist across ASD and controls, and that FC subtypes can reveal unique information about brain-behavior relationships.

Physical characteristics not psychological state or trait characteristics predict motion during resting state fMRI

Head motion (HM) during fMRI acquisition can significantly affect measures of brain activity or connectivity even after correction with preprocessing methods. Moreover, any systematic relationship between HM and variables of interest can introduce systematic bias. There is a large and growing interest in identifying neural biomarkers for psychiatric disorders using resting state fMRI (rsfMRI). However, the relationship between HM and different psychiatric symptoms domains is not well understood. The aim of this investigation was to determine whether psychiatric symptoms and other characteristics of the individual predict HM during rsfMRI. A sample of n = 464 participants (174 male) from the Tulsa1000, a naturalistic longitudinal study recruiting subjects with different levels of severity in mood/anxiety/substance use disorders based on the dimensional NIMH Research Domain Criteria framework was used for this study. Based on a machine learning (ML) pipeline with nested cross-validation to avoid overfitting, the stacked model with 15 anthropometric (like body mass index, BMI) and demographic (age and sex) variables identifies BMI and weight as the most important variables and explained 10.9 percent of the HM variance (95% CI: 9.9–11.8). In comparison ML models with 105 self-report measures for state and trait psychological characteristics identified nicotine and alcohol use variables as well as impulsivity inhibitory control variables but explain only 5 percent of HM variance (95% CI: 3.5–6.4). A combined ML model using all 120 variables did not perform significantly better than the model using only 15 physical variables (combined model 95% confidence interval: 10.2–12.4). Taken together, after considering physical variables, state or trait psychological characteristics do not provide additional power to predict motion during rsfMRI.