Structure-function coupling in highly sampled individual brains.

Structural connectivity (SC) between distant regions of the brain support synchronized function known as functional connectivity (FC) and give rise to the large-scale brain networks that enable cognition and behavior. Understanding how SC enables FC is important to understand how injuries to SC may alter brain function and cognition. Previous work evaluating whole-brain SC-FC relationships showed that SC explained FC well in unimodal visual and motor areas, but only weakly in association areas, suggesting a unimodal-heteromodal gradient organization of SC-FC coupling. However, this work was conducted in group-averaged SC/FC data. Thus, it could not account for inter-individual variability in the locations of cortical areas and white matter tracts. We evaluated the correspondence of SC and FC within three highly sampled healthy participants. For each participant, we collected 78 min of diffusion-weighted MRI for SC and 360 min of resting state fMRI for FC. We found that FC was best explained by SC in visual and motor systems, as well as in anterior and posterior cingulate regions. A unimodal-to-heteromodal gradient could not fully explain SC-FC coupling. We conclude that the SC-FC coupling of the anterior-posterior cingulate circuit is more similar to unimodal areas than to heteromodal areas.

Cingulo-opercular control network and disused motor circuits joined in standby mode.

Whole-brain resting-state functional MRI (rs-fMRI) during 2 wk of upper-limb casting revealed that disused motor regions became more strongly connected to the cingulo-opercular network (CON), an executive control network that includes regions of the dorsal anterior cingulate cortex (dACC) and insula. Disuse-driven increases in functional connectivity (FC) were specific to the CON and somatomotor networks and did not involve any other networks, such as the salience, frontoparietal, or default mode networks. Censoring and modeling analyses showed that FC increases during casting were mediated by large, spontaneous activity pulses that appeared in the disused motor regions and CON control regions. During limb constraint, disused motor circuits appear to enter a standby mode characterized by spontaneous activity pulses and strengthened connectivity to CON executive control regions.

Comparison of Resting State Functional Connectivity in Persons With and Without HIV: A Cross-sectional Study.

BACKGROUND: This study examined the effects of human immunodeficiency virus (HIV) on resting state functional connectivity (RSFC) in a large cohort of people with HIV (PWH) and healthy controls without HIV (PWoH). Within PWH analyses focused on the effects of viral suppression and cognitive impairment on RSFC. METHODS: A total of 316 PWH on stable combination antiretroviral therapy and 209 demographically matched PWoH were scanned at a single institution. Effects of the virus were examined by grouping PWH by detectable (viral load > 20 copies/mL; VLD) and undetectable (VLU) viral loads and as being cognitively impaired (CI) (Global Deficit Score ≥ 0.5) or cognitively normal (CN). Regression analysis, object oriented data analysis, and spring embedded graph models were applied to RSFC measures from 298 established brain regions of interest comprising 13 brain networks to examine group differences. RESULTS: No significant RSFC differences were observed between PWH and PWoH. Within PWH, there were no significant differences in RSFC between VLD and VLU subgroups and CI and CN subgroups. CONCLUSIONS: There were no significant effects of HIV on RSFC in our relatively large cohort of PWH and PWoH. Future studies could increase the sample size and combine with other imaging modalities.

White matter hyperintensity longitudinal morphometric analysis in association with Alzheimer disease.

INTRODUCTION: Vascular damage in Alzheimer's disease (AD) has shown conflicting findings particularly when analyzing longitudinal data. We introduce white matter hyperintensity (WMH) longitudinal morphometric analysis (WLMA) that quantifies WMH expansion as the distance from lesion voxels to a region of interest boundary. METHODS: WMH segmentation maps were derived from 270 longitudinal fluid-attenuated inversion recovery (FLAIR) ADNI images. WLMA was performed on five data-driven WMH patterns with distinct spatial distributions. Amyloid accumulation was evaluated with WMH expansion across the five WMH patterns. RESULTS: The preclinical group had significantly greater expansion in the posterior ventricular WM compared to controls. Amyloid significantly associated with frontal WMH expansion primarily within AD individuals. WLMA outperformed WMH volume changes for classifying AD from controls primarily in periventricular and posterior WMH. DISCUSSION: These data support the concept that localized WMH expansion continues to proliferate with amyloid accumulation throughout the entirety of the disease in distinct spatial locations.

BOLD cofluctuation 'events' are predicted from static functional connectivity.

Recent work identified single time points ("events") of high regional cofluctuation in functional Magnetic Resonance Imaging (fMRI) which contain more large-scale brain network information than other, low cofluctuation time points. This suggested that events might be a discrete, temporally sparse signal which drives functional connectivity (FC) over the timeseries. However, a different, not yet explored possibility is that network information differences between time points are driven by sampling variability on a constant, static, noisy signal. Using a combination of real and simulated data, we examined the relationship between cofluctuation and network structure and asked if this relationship was unique, or if it could arise from sampling variability alone. First, we show that events are not discrete - there is a gradually increasing relationship between network structure and cofluctuation; ∼50% of samples show very strong network structure. Second, using simulations we show that this relationship is predicted from sampling variability on static FC. Finally, we show that randomly selected points can capture network structure about as well as events, largely because of their temporal spacing. Together, these results suggest that, while events exhibit particularly strong representations of static FC, there is little evidence that events are unique timepoints that drive FC structure. Instead, a parsimonious explanation for the data is that events arise from a single static, but noisy, FC structure.

Accuracy and reliability of diffusion imaging models.

Diffusion imaging aims to non-invasively characterize the anatomy and integrity of the brain's white matter fibers. We evaluated the accuracy and reliability of commonly used diffusion imaging methods as a function of data quantity and analysis method, using both simulations and highly sampled individual-specific data (927-1442 diffusion weighted images [DWIs] per individual). Diffusion imaging methods that allow for crossing fibers (FSL's BedpostX [BPX], DSI Studio's Constant Solid Angle Q-Ball Imaging [CSA-QBI], MRtrix3's Constrained Spherical Deconvolution [CSD]) estimated excess fibers when insufficient data were present and/or when the data did not match the model priors. To reduce such overfitting, we developed a novel Bayesian Multi-tensor Model-selection (BaMM) method and applied it to the popular ball-and-stick model used in BedpostX within the FSL software package. BaMM was robust to overfitting and showed high reliability and the relatively best crossing-fiber accuracy with increasing amounts of diffusion data. Thus, sufficient data and an overfitting resistant analysis method enhance precision diffusion imaging. For potential clinical applications of diffusion imaging, such as neurosurgical planning and deep brain stimulation (DBS), the quantities of data required to achieve diffusion imaging reliability are lower than those needed for functional MRI.

Trait-like variants in human functional brain networks.

Resting-state functional magnetic resonance imaging (fMRI) has provided converging descriptions of group-level functional brain organization. Recent work has revealed that functional networks identified in individuals contain local features that differ from the group-level description. We define these features as network variants. Building on these studies, we ask whether distributions of network variants reflect stable, trait-like differences in brain organization. Across several datasets of highly-sampled individuals we show that 1) variants are highly stable within individuals, 2) variants are found in characteristic locations and associate with characteristic functional networks across large groups, 3) task-evoked signals in variants demonstrate a link to functional variation, and 4) individuals cluster into subgroups on the basis of variant characteristics that are related to differences in behavior. These results suggest that distributions of network variants may reflect stable, trait-like, functionally relevant individual differences in functional brain organization.

Probabilistic mapping of human functional brain networks identifies regions of high group consensus.

Many recent developments surrounding the functional network organization of the human brain have focused on data that have been averaged across groups of individuals. While such group-level approaches have shed considerable light on the brain's large-scale distributed systems, they conceal individual differences in network organization, which recent work has demonstrated to be common and widespread. This individual variability produces noise in group analyses, which may average together regions that are part of different functional systems across participants, limiting interpretability. However, cost and feasibility constraints may limit the possibility for individual-level mapping within studies. Here our goal was to leverage information about individual-level brain organization to probabilistically map common functional systems and identify locations of high inter-subject consensus for use in group analyses. We probabilistically mapped 14 functional networks in multiple datasets with relatively high amounts of data. All networks show "core" (high-probability) regions, but differ from one another in the extent of their higher-variability components. These patterns replicate well across four datasets with different participants and scanning parameters. We produced a set of high-probability regions of interest (ROIs) from these probabilistic maps; these and the probabilistic maps are made publicly available, together with a tool for querying the network membership probabilities associated with any given cortical location. These quantitative estimates and public tools may allow researchers to apply information about inter-subject consensus to their own fMRI studies, improving inferences about systems and their functional specializations.

A Critical, Event-Related Appraisal of Denoising in Resting-State fMRI Studies.

This article advances two parallel lines of argument about resting-state functional magnetic resonance imaging (fMRI) signals, one empirical and one conceptual. The empirical line creates a four-part organization of the text: (1) head motion and respiration commonly cause distinct, major, unwanted influences (artifacts) in fMRI signals; (2) head motion and respiratory changes are, confoundingly, both related to psychological and clinical and biological variables of interest; (3) many fMRI denoising strategies fail to identify and remove one or the other kind of artifact; and (4) unremoved artifact, due to correlations of artifacts with variables of interest, renders studies susceptible to identifying variance of noninterest as variance of interest. Arising from these empirical observations is a conceptual argument: that an event-related approach to task-free scans, targeting common behaviors during scanning, enables fundamental distinctions among the kinds of signals present in the data, information which is vital to understanding the effects of denoising procedures. This event-related perspective permits statements like "Event X is associated with signals A, B, and C, each with particular spatial, temporal, and signal decay properties". Denoising approaches can then be tailored, via performance in known events, to permit or suppress certain kinds of signals based on their desirability.

Removal of high frequency contamination from motion estimates in single-band fMRI saves data without biasing functional connectivity.

Denoising fMRI data requires assessment of frame-to-frame head motion and removal of the biases motion introduces. This is usually done through analysis of the parameters calculated during retrospective head motion correction (i.e., 'motion' parameters). However, it is increasingly recognized that respiration introduces factitious head motion via perturbations of the main (B0) field. This effect appears as higher-frequency fluctuations in the motion parameters (>0.1 â€‹Hz, here referred to as 'HF-motion'), primarily in the phase-encoding direction. This periodicity can sometimes be obscured in standard single-band fMRI (TR 2.0-2.5 â€‹s) due to aliasing. Here we examined (1) how prevalent HF-motion effects are in seven single-band datasets with TR from 2.0 to 2.5 â€‹s and (2) how HF-motion affects functional connectivity. We demonstrate that HF-motion is more common in older adults, those with higher body mass index, and those with lower cardiorespiratory fitness. We propose a low-pass filtering approach to remove the contamination of high frequency effects from motion summary measures, such as framewise displacement (FD). We demonstrate that in most datasets this filtering approach saves a substantial amount of data from FD-based frame censoring, while at the same time reducing motion biases in functional connectivity measures. These findings suggest that filtering motion parameters is an effective way to improve the fidelity of head motion estimates, even in single band datasets. Particularly large data savings may accrue in datasets acquired in older and less fit participants.

Correction of respiratory artifacts in MRI head motion estimates.

Head motion represents one of the greatest technical obstacles in magnetic resonance imaging (MRI) of the human brain. Accurate detection of artifacts induced by head motion requires precise estimation of movement. However, head motion estimates may be corrupted by artifacts due to magnetic main field fluctuations generated by body motion. In the current report, we examine head motion estimation in multiband resting state functional connectivity MRI (rs-fcMRI) data from the Adolescent Brain and Cognitive Development (ABCD) Study and comparison 'single-shot' datasets. We show that respirations contaminate movement estimates in functional MRI and that respiration generates apparent head motion not associated with functional MRI quality reductions. We have developed a novel approach using a band-stop filter that accurately removes these respiratory effects from motion estimates. Subsequently, we demonstrate that utilizing a band-stop filter improves post-processing fMRI data quality. Lastly, we demonstrate the real-time implementation of motion estimate filtering in our FIRMM (Framewise Integrated Real-Time MRI Monitoring) software package.

Reward-related regions form a preferentially coupled system at rest.

Neuroimaging studies have implicated a set of striatal and orbitofrontal cortex (OFC) regions that are commonly activated during reward processing tasks. Resting-state functional connectivity (RSFC) studies have demonstrated that the human brain is organized into several functional systems that show strong temporal coherence in the absence of goal-directed tasks. Here we use seed-based and graph-theory RSFC approaches to characterize the systems-level organization of putative reward regions of at rest. Peaks of connectivity from seed-based RSFC patterns for the nucleus accumbens (NAcc) and orbitofrontal cortex (OFC) were used to identify candidate reward regions which were merged with a previously used set of regions (Power et al., 2011). Graph-theory was then used to determine system-level membership for all regions. Several regions previously implicated in reward-processing (NAcc, lateral and medial OFC, and ventromedial prefrontal cortex) comprised a distinct, preferentially coupled system. This RSFC system is stable across a range of connectivity thresholds and shares strong overlap with meta-analyses of task-based reward studies. This reward system shares between-system connectivity with systems implicated in cognitive control and self-regulation, including the fronto-parietal, cingulo-opercular, and default systems. Differences may exist in the pathways through which control systems interact with reward system components. Whereas NAcc is functionally connected to cingulo-opercular and default systems, OFC regions show stronger connectivity with the fronto-parietal system. We propose that future work may be able to interrogate group or individual differences in connectivity profiles using the regions delineated in this work to explore potential relationships to appetitive behaviors, self-regulation failure, and addiction.

Individual Variability of the System-Level Organization of the Human Brain.

Recent functional magnetic resonance imaging-based resting-state functional connectivity analyses of group average data have characterized large-scale systems that represent a high level in the organizational hierarchy of the human brain. These systems are likely to vary spatially across individuals, even after anatomical alignment, but the characteristics of this variance are unknown. Here, we characterized large-scale brain systems across two independent datasets of young adults. In these individuals, we were able to identify brain systems that were similar to those described in the group average, and we observed that individuals had consistent topological arrangement of the system features present in the group average. However, the size of system features varied across individuals in systematic ways, such that expansion of one feature of a given system predicted expansion of other parts of the system. Individual-specific systems also contained unique topological features not present in group average systems; some of these features were consistent across a minority of individuals. These effects were observed even after controlling for data quality and for the accuracy of anatomical registration. The variability characterized here has important implications for cognitive neuroscience investigations, which often assume the functional equivalence of aligned brain regions across individuals.

On the Stability of BOLD fMRI Correlations.

Measurement of correlations between brain regions (functional connectivity) using blood oxygen level dependent (BOLD) fMRI has proven to be a powerful tool for studying the functional organization of the brain. Recently, dynamic functional connectivity has emerged as a major topic in the resting-state BOLD fMRI literature. Here, using simulations and multiple sets of empirical observations, we confirm that imposed task states can alter the correlation structure of BOLD activity. However, we find that observations of "dynamic" BOLD correlations during the resting state are largely explained by sampling variability. Beyond sampling variability, the largest part of observed "dynamics" during rest is attributable to head motion. An additional component of dynamic variability during rest is attributable to fluctuating sleep state. Thus, aside from the preceding explanatory factors, a single correlation structure-as opposed to a sequence of distinct correlation structures-may adequately describe the resting state as measured by BOLD fMRI. These results suggest that resting-state BOLD correlations do not primarily reflect moment-to-moment changes in cognitive content. Rather, resting-state BOLD correlations may predominantly reflect processes concerned with the maintenance of the long-term stability of the brain's functional organization.

Individual-specific features of brain systems identified with resting state functional correlations.

Recent work has made important advances in describing the large-scale systems-level organization of human cortex by analyzing functional magnetic resonance imaging (fMRI) data averaged across groups of subjects. However, new findings have emerged suggesting that individuals' cortical systems are topologically complex, containing small but reliable features that cannot be observed in group-averaged datasets, due in part to variability in the position of such features along the cortical sheet. This previous work has reported only specific examples of these individual-specific system features; to date, such features have not been comprehensively described. Here we used fMRI to identify cortical system features in individual subjects within three large cross-subject datasets and one highly sampled within-subject dataset. We observed system features that have not been previously characterized, but 1) were reliably detected across many scanning sessions within a single individual, and 2) could be matched across many individuals. In total, we identified forty-three system features that did not match group-average systems, but that replicated across three independent datasets. We described the size and spatial distribution of each non-group feature. We further observed that some individuals were missing specific system features, suggesting individual differences in the system membership of cortical regions. Finally, we found that individual-specific system features could be used to increase subject-to-subject similarity. Together, this work identifies individual-specific features of human brain systems, thus providing a catalog of previously unobserved brain system features and laying the foundation for detailed examinations of brain connectivity in individuals.

Postoperative seizure freedom does not normalize altered connectivity in temporal lobe epilepsy.

OBJECTIVES: Specific changes in the functional connectivity of brain networks occur in patients with epilepsy. Yet whether such changes reflect a stable disease effect or one that is a function of active seizure burden remains unclear. Here, we longitudinally assessed the connectivity of canonical cognitive functional networks in patients with intractable temporal lobe epilepsy (TLE), both before and after patients underwent epilepsy surgery and achieved seizure freedom. METHODS: Seventeen patients with intractable TLE who underwent epilepsy surgery with Engel class I outcome and 17 matched healthy controls took part in the study. The functional connectivity of a set of cognitive functional networks derived from typical cognitive tasks was assessed in patients, preoperatively and postoperatively, as well as in controls, using stringent methods of artifact reduction. RESULTS: Preoperatively, functional networks in TLE patients differed significantly from healthy controls, with differences that largely, but not exclusively, involved the default mode and temporal/auditory subnetworks. However, undergoing epilepsy surgery and achieving seizure freedom did not lead to significant changes in network connectivity, with postoperative functional network abnormalities closely mirroring the preoperative state. SIGNIFICANCE: This result argues for a stable chronic effect of the disease on brain connectivity, with changes that are largely "burned in" by the time a patient with intractable TLE undergoes epilepsy surgery, which typically occurs years after the initial diagnosis. The result has potential implications for the treatment of intractable epilepsy, suggesting that delaying surgical intervention that may achieve seizure freedom may lead to functional network changes that are no longer reversible by the time of epilepsy surgery.

Generation and Evaluation of a Cortical Area Parcellation from Resting-State Correlations.

The cortical surface is organized into a large number of cortical areas; however, these areas have not been comprehensively mapped in the human. Abrupt transitions in resting-state functional connectivity (RSFC) patterns can noninvasively identify locations of putative borders between cortical areas (RSFC-boundary mapping; Cohen et al. 2008). Here we describe a technique for using RSFC-boundary maps to define parcels that represent putative cortical areas. These parcels had highly homogenous RSFC patterns, indicating that they contained one unique RSFC signal; furthermore, the parcels were much more homogenous than a null model matched for parcel size when tested in two separate datasets. Several alternative parcellation schemes were tested this way, and no other parcellation was as homogenous as or had as large a difference compared with its null model. The boundary map-derived parcellation contained parcels that overlapped with architectonic mapping of areas 17, 2, 3, and 4. These parcels had a network structure similar to the known network structure of the brain, and their connectivity patterns were reliable across individual subjects. These observations suggest that RSFC-boundary map-derived parcels provide information about the location and extent of human cortical areas. A parcellation generated using this method is available at http://www.nil.wustl.edu/labs/petersen/Resources.html.

Multivariate pattern classification of pediatric Tourette syndrome using functional connectivity MRI.

Tourette syndrome (TS) is a developmental neuropsychiatric disorder characterized by motor and vocal tics. Individuals with TS would benefit greatly from advances in prediction of symptom timecourse and treatment effectiveness. As a first step, we applied a multivariate method - support vector machine (SVM) classification - to test whether patterns in brain network activity, measured with resting state functional connectivity (RSFC) MRI, could predict diagnostic group membership for individuals. RSFC data from 42 children with TS (8-15 yrs) and 42 unaffected controls (age, IQ, in-scanner movement matched) were included. While univariate tests identified no significant group differences, SVM classified group membership with ~70% accuracy (p < .001). We also report a novel adaptation of SVM binary classification that, in addition to an overall accuracy rate for the SVM, provides a confidence measure for the accurate classification of each individual. Our results support the contention that multivariate methods can better capture the complexity of some brain disorders, and hold promise for predicting prognosis and treatment outcome for individuals with TS.

Two common and distinct forms of variation in human functional brain networks.

The cortex has a characteristic layout with specialized functional areas forming distributed large-scale networks. However, substantial work shows striking variation in this organization across people, which relates to differences in behavior. While most previous work treats individual differences as linked to boundary shifts between the borders of regions, here we show that cortical 'variants' also occur at a distance from their typical position, forming ectopic intrusions. Both 'border' and 'ectopic' variants are common across individuals, but differ in their location, network associations, properties of subgroups of individuals, activations during tasks, and prediction of behavioral phenotypes. Border variants also track significantly more with shared genetics than ectopic variants, suggesting a closer link between ectopic variants and environmental influences. This work argues that these two dissociable forms of variation-border shifts and ectopic intrusions-must be separately accounted for in the analysis of individual differences in cortical systems across people.

Network-level enrichment provides a framework for biological interpretation of machine learning results.

Machine learning algorithms are increasingly being utilized to identify brain connectivity biomarkers linked to behavioral and clinical outcomes. However, research often prioritizes prediction accuracy at the expense of biological interpretability, and inconsistent implementation of ML methods may hinder model accuracy. To address this, our paper introduces a network-level enrichment approach, which integrates brain system organization in the context of connectome-wide statistical analysis to reveal network-level links between brain connectivity and behavior. To demonstrate the efficacy of this approach, we used linear support vector regression (LSVR) models to examine the relationship between resting-state functional connectivity networks and chronological age. We compared network-level associations based on raw LSVR weights to those produced from the forward and inverse models. Results indicated that not accounting for shared family variance inflated prediction performance, the k-best feature selection via Pearson correlation reduced accuracy and reliability, and raw LSVR model weights produced network-level associations that deviated from the significant brain systems identified by forward and inverse models. Our findings offer crucial insights for applying machine learning to neuroimaging data, emphasizing the value of network enrichment for biological interpretation.