Effect of scanning duration and sample size on reliability in resting state fMRI dynamic causal modeling analysis

Despite its widespread use, resting-state functional magnetic resonance imaging (rsfMRI) has been criticized for low test-retest reliability. To improve reliability, researchers have recommended using extended scanning durations, increased sample size, and advanced brain connectivity techniques. However, longer scanning runs and larger sample sizes may come with practical challenges and burdens, especially in rare populations. Here we tested if an advanced brain connectivity technique, dynamic causal modeling (DCM), can improve reliability of fMRI effective connectivity (EC) metrics to acceptable levels without extremely long run durations or extremely large samples. Specifically, we employed DCM for EC analysis on rsfMRI data from the Human Connectome Project. To avoid bias, we assessed four distinct DCMs and gradually increased sample sizes in a randomized manner across ten permutations. We employed pseudo true positive and pseudo false positive rates to assess the efficacy of shorter run durations (3.6, 7.2, 10.8, 14.4 min) in replicating the outcomes of the longest scanning duration (28.8 min) when the sample size was fixed at the largest (n = 160 subjects). Similarly, we assessed the efficacy of smaller sample sizes (n = 10, 20, …, 150 subjects) in replicating the outcomes of the largest sample (n = 160 subjects) when the scanning duration was fixed at the longest (28.8 min). Our results revealed that the pseudo false positive rate was below 0.05 for all the analyses. After the scanning duration reached 10.8 min, which yielded a pseudo true positive rate of 92%, further extensions in run time showed no improvements in pseudo true positive rate. Expanding the sample size led to enhanced pseudo true positive rate outcomes, with a plateau at n = 70 subjects for the targeted top one-half of the largest ECs in the reference sample, regardless of whether the longest run duration (28.8 min) or the viable run duration (10.8 min) was employed. Encouragingly, smaller sample sizes exhibited pseudo true positive rates of approximately 80% for n = 20, and 90% for n = 40 subjects. These data suggest that advanced DCM analysis may be a viable option to attain reliable metrics of EC when larger sample sizes or run times are not feasible.

Spatial Bayesian GLM on the cortical surface produces reliable task activations in individuals and groups

The general linear model (GLM) is a widely popular and convenient tool for estimating the functional brain response and identifying areas of significant activation during a task or stimulus. However, the classical GLM is based on a massive univariate approach that does not explicitly leverage the similarity of activation patterns among neighboring brain locations. As a result, it tends to produce noisy estimates and be underpowered to detect significant activations, particularly in individual subjects and small groups. A recently proposed alternative, a cortical surface-based spatial Bayesian GLM, leverages spatial dependencies among neighboring cortical vertices to produce more accurate estimates and areas of functional activation. The spatial Bayesian GLM can be applied to individual and group-level analysis. In this study, we assess the reliability and power of individual and group-average measures of task activation produced via the surface-based spatial Bayesian GLM. We analyze motor task data from 45 subjects in the Human Connectome Project (HCP) and HCP Retest datasets. We also extend the model to multi-run analysis and employ subject-specific cortical surfaces rather than surfaces inflated to a sphere for more accurate distance-based modeling. Results show that the surface-based spatial Bayesian GLM produces highly reliable activations in individual subjects and is powerful enough to detect trait-like functional topologies. Additionally, spatial Bayesian modeling enhances reliability of group-level analysis even in moderately sized samples (n = 45). Notably, the power of the spatial Bayesian GLM to detect activations above a scientifically meaningful effect size is nearly invariant to sample size, exhibiting high power even in small samples (n = 10). The spatial Bayesian GLM is computationally efficient in individuals and groups and is convenient to implement with the open-source BayesfMRI R package.

Validating dynamicity in resting state fMRI with activation-informed temporal segmentation

Confirming the presence (or absence) of dynamic functional connectivity (dFC) states during rest is an important open question in the field of cognitive neuroscience. The prevailing dFC framework aims to identify dynamics directly from connectivity estimates with a sliding window approach, however this method suffers from several drawbacks including sensitivity to window size and poor test-retest reliability. We hypothesize that time-varying changes in functional connectivity are mirrored by significant temporal changes in functional activation, and that this coupling can be leveraged to study dFC without the need for a predefined sliding window. Here, we introduce a data-driven dFC framework, which involves informed segmentation of fMRI time series at candidate FC state transition points estimated from changes in whole-brain functional activation, rather than a fixed-length sliding window. We show our approach reliably identifies true cognitive state change points when applied on block-design working memory task data and outperforms the standard sliding window approach in both accuracy and computational efficiency in this context. When applied to data from four resting state fMRI scanning sessions, our method consistently recovers five reliable FC states, and subject-specific features derived from these states show significant correlation with behavioral phenotypes of interest (cognitive ability, personality). Overall, these results suggest abrupt whole-brain changes in activation can be used as a marker for changes in connectivity states and provides new evidence for the existence of time-varying FC in rest.

The bright side and dark side of daydreaming predict creativity together through brain functional connectivity

Daydreaming and creativity have similar cognitive processes and neural basis. However, few empirical studies have examined the relationship between daydreaming and creativity using cognitive neuroscience methods. The present study explored the relationship between different types of daydreaming and creativity and their common neural basis. The behavioral results revealed that positive constructive daydreaming is positively related to creativity, while poor attentional control is negatively related to it. Machine learning framework was adopted to examine the predictive effect of daydreaming-related brain functional connectivity (FC) on creativity. The results demonstrated that task FCs related to positive constructive daydreaming and task FCs related to poor attentional control both predicted an individual's creativity score successfully. In addition, task FCs combining the positive constructive daydreaming and poor attentional control also had significant predictive effect on creativity score. Furthermore, predictive analysis based on resting-state FCs showed similar patterns. Both of the subscale-related FCs and combined FCs had significant predictive effect on creativity score. Further analysis showed the task and the resting-state FCs both mainly located in the default mode network, central executive network, salience network, and attention network. These results showed that daydreaming was closely related to creativity, as they shared common FC basis.

A systematic evaluation of source reconstruction of resting MEG of the human brain with a new high-resolution atlas: Performance, precision, and parcellation

Noninvasive functional neuroimaging of the human brain can give crucial insight into the mechanisms that underpin healthy cognition and neurological disorders. Magnetoencephalography (MEG) measures extracranial magnetic fields originating from neuronal activity with high temporal resolution, but requires source reconstruction to make neuroanatomical inferences from these signals. Many source reconstruction algorithms are available, and have been widely evaluated in the context of localizing task-evoked activities. However, no consensus yet exists on the optimum algorithm for resting-state data. Here, we evaluated the performance of six commonly-used source reconstruction algorithms based on minimum-norm and beamforming estimates. Using human resting-state MEG, we compared the algorithms using quantitative metrics, including resolution properties of inverse solutions and explained variance in sensor-level data. Next, we proposed a data-driven approach to reduce the atlas from the Human Connectome Project's multi-modal parcellation of the human cortex based on metrics such as MEG signal-to-noise-ratio and resting-state functional connectivity gradients. This procedure produced a reduced cortical atlas with 230 regions, optimized to match the spatial resolution and the rank of MEG data from the current generation of MEG scanners. Our results show that there is no "one size fits all" algorithm, and make recommendations on the appropriate algorithms depending on the data and aimed analyses. Our comprehensive comparisons and recommendations can serve as a guide for choosing appropriate methodologies in future studies of resting-state MEG.

Exploring MEG brain fingerprints: Evaluation, pitfalls, and interpretations

Individual characterization of subjects based on their functional connectome (FC), termed "FC fingerprinting", has become a highly sought-after goal in contemporary neuroscience research. Recent functional magnetic resonance imaging (fMRI) studies have demonstrated unique characterization and accurate identification of individuals as an accomplished task. However, FC fingerprinting in magnetoencephalography (MEG) data is still widely unexplored. Here, we study resting-state MEG data from the Human Connectome Project to assess the MEG FC fingerprinting and its relationship with several factors including amplitude- and phase-coupling functional connectivity measures, spatial leakage correction, frequency bands, and behavioral significance. To this end, we first employ two identification scoring methods, differential identifiability and success rate, to provide quantitative fingerprint scores for each FC measurement. Secondly, we explore the edgewise and nodal MEG fingerprinting patterns across the different frequency bands (delta, theta, alpha, beta, and gamma). Finally, we investigate the cross-modality fingerprinting patterns obtained from MEG and fMRI recordings from the same subjects. We assess the behavioral significance of FC across connectivity measures and imaging modalities using partial least square correlation analyses. Our results suggest that fingerprinting performance is heavily dependent on the functional connectivity measure, frequency band, identification scoring method, and spatial leakage correction. We report higher MEG fingerprinting performances in phase-coupling methods, central frequency bands (alpha and beta), and in the visual, frontoparietal, dorsal-attention, and default-mode networks. Furthermore, cross-modality comparisons reveal a certain degree of spatial concordance in fingerprinting patterns between the MEG and fMRI data, especially in the visual system. Finally, the multivariate correlation analyses show that MEG connectomes have strong behavioral significance, which however depends on the considered connectivity measure and temporal scale. This comprehensive, albeit preliminary investigation of MEG connectome test-retest identifiability offers a first characterization of MEG fingerprinting in relation to different methodological and electrophysiological factors and contributes to the understanding of fingerprinting cross-modal relationships. We hope that this first investigation will contribute to setting the grounds for MEG connectome identification.

Boost in Test-Retest Reliability in Resting State fMRI with Predictive Modeling

Recent studies found low test-retest reliability in functional magnetic resonance imaging (fMRI), raising serious concerns among researchers, but these studies mostly focused on the reliability of individual fMRI features (e.g., individual connections in resting state connectivity maps). Meanwhile, neuroimaging researchers increasingly employ multivariate predictive models that aggregate information across a large number of features to predict outcomes of interest, but the test-retest reliability of predicted outcomes of these models has not previously been systematically studied. Here we apply 10 predictive modeling methods to resting state connectivity maps from the Human Connectome Project dataset to predict 61 outcome variables. Compared with mean reliability of individual resting state connections, we find mean reliability of the predicted outcomes of predictive models is substantially higher for all 10 modeling methods assessed. Moreover, improvement was consistently observed across all scanning and processing choices (i.e., scan lengths, censoring thresholds, volume- vs. surface-based processing). For the most reliable methods, the reliability of predicted outcomes was mostly, though not exclusively, in the "good" range (above 0.60). Finally, we identified three mechanisms that help to explain why predicted outcomes of predictive models have higher reliability than individual imaging features. We conclude that researchers can potentially achieve higher test-retest reliability by making greater use of predictive models.

Bayesian estimation of maximum entropy model for individualized energy landscape analysis of brain state dynamics

The pairwise maximum entropy model (MEM) for resting state functional MRI (rsfMRI) has been used to generate energy landscape of brain states and to explore nonlinear brain state dynamics. Researches using MEM, however, has mostly been restricted to fixed‐effect group‐level analyses, using concatenated time series across individuals, due to the need for large samples in the parameter estimation of MEM. To mitigate the small sample problem in analyzing energy landscapes for individuals, we propose a Bayesian estimation of individual MEM using variational Bayes approximation (BMEM). We evaluated the performances of BMEM with respect to sample sizes and prior information using simulation. BMEM showed advantages over conventional maximum likelihood estimation in reliably estimating model parameters for individuals with small sample data, particularly utilizing the empirical priors derived from group data. We then analyzed individual rsfMRI of the Human Connectome Project to show the usefulness of MEM in differentiating individuals and in exploring neural correlates for human behavior. MEM and its energy landscape properties showed high subject specificity comparable to that of functional connectivity. Canonical correlation analysis identified canonical variables for MEM highly associated with cognitive scores. Inter‐individual variations of cognitive scores were also reflected in energy landscape properties such as energies, occupation times, and basin sizes at local minima. We conclude that BMEM provides an efficient method to characterize dynamic properties of individuals using energy landscape analysis of individual brain states.

Association of Prenatal Acetaminophen Exposure Measured in Meconium With Risk of Attention-Deficit/Hyperactivity Disorder Mediated by Frontoparietal Network Brain Connectivity

Importance

Despite evidence of an association between prenatal acetaminophen exposure and attention-deficit/hyperactivity disorder (ADHD) in offspring, the drug is not contraindicated during pregnancy, possibly because prior studies have relied on maternal self-report, failed to quantify acetaminophen dose, and lacked mechanistic insight.

Objective

To examine the association between prenatal acetaminophen exposure measured in meconium (hereinafter referred to as meconium acetaminophen) and ADHD in children aged 6 to 7 years, along with the potential for mediation by functional brain connectivity.

Design, Setting, and Participants

This prospective birth cohort study from the Centre Hospitalier Université de Sherbrooke in Sherbrooke, Québec, Canada, included 394 eligible children, of whom 345 had meconium samples collected at delivery and information on ADHD diagnosis. Mothers were enrolled from September 25, 2007, to September 10, 2009, at their first prenatal care visit or delivery and were followed up when children were aged 6 to 7 years. When children were aged 9 to 11 years, resting-state brain connectivity was assessed with magnetic resonance imaging. Data for the present study were collected from September 25, 2007, to January 18, 2020, and analyzed from January 7, 2019, to January 22, 2020.

Exposures

Acetaminophen levels measured in meconium.

Main Outcomes and Measures

Physician diagnosis of ADHD was determined at follow-up when children were aged 6 to 7 years or from medical records. Resting-state brain connectivity was assessed with magnetic resonance imaging; attention problems and hyperactivity were assessed with the Behavioral Assessment System for Children Parent Report Scale. Associations between meconium acetaminophen levels and outcomes were estimated with linear and logistic regressions weighted on the inverse probability of treatment to account for potential confounders. Causal mediation analysis was used to test for mediation of the association between prenatal acetaminophen exposure and hyperactivity by resting-state brain connectivity.

Results

Among the 345 children included in the analysis (177 boys [51.3%]; mean [SD] age, 6.58 [0.54] years), acetaminophen was detected in 199 meconium samples (57.7%), and ADHD was diagnosed in 33 children (9.6%). Compared with no acetaminophen, detection of acetaminophen in meconium was associated with increased odds of ADHD (odds ratio [OR], 2.43; 95% CI, 1.41-4.21). A dose-response association was detected; each doubling of exposure increased the odds of ADHD by 10% (OR, 1.10; 95% CI, 1.02-1.19). Children with acetaminophen detected in meconium showed increased negative connectivity between frontoparietal and default mode network nodes to clusters in the sensorimotor cortices, which mediated an indirect effect on increased child hyperactivity (14%; 95% CI, 1%-26%).

Conclusions and Relevance

Together with the multitude of other cohort studies showing adverse neurodevelopment associated with prenatal acetaminophen exposure, this work suggests caution should be used in administering acetaminophen during pregnancy. Research into alternative pain management strategies for pregnant women could be beneficial.

Resting state functional network switching rate is differently altered in bipolar disorder and major depressive disorder

The clinical misdiagnosis ratio of bipolar disorder (BD) patients to major depressive disorder (MDD) patients is high. Recent findings hypothesize that the ability to flexibly recruit functional neural networks is differently altered in BD and MDD patients. This study aimed to explore distinct aberrance of network flexibility during dynamic networks configuration in BD and MDD patients. Resting state functional magnetic resonance imaging of 40 BD patients, 61 MDD patients, and 61 matched healthy controls were recruited. Dynamic functional connectivity matrices for each subject were constructed with a sliding window method. Then, network switching rate of each node was calculated and compared among the three groups. BD and MDD patients shared decreased network switching rate of regions including left precuneus, bilateral parahippocampal gyrus, and bilateral dorsal medial prefrontal cortex. Apart from these regions, MDD patients presented specially decreased network switching rate in the bilateral anterior insula, left amygdala, and left striatum. Taken together, BD and MDD patients shared decreased network switching rate of key hubs in default mode network and MDD patients presented specially decreased switching rate in salience network and striatum. We found shared and distinct aberrance of network flexibility which revealed altered adaptive functions during dynamic networks configuration of BD and MDD.

Alterations and test-retest reliability of functional connectivity network measures in cerebral small vessel disease

While structural network analysis consolidated the hypothesis of cerebral small vessel disease (SVD) being a disconnection syndrome, little is known about functional changes on the level of brain networks. In patients with genetically defined SVD (CADASIL, n = 41) and sporadic SVD (n = 46), we independently tested the hypothesis that functional networks change with SVD burden and mediate the effect of disease burden on cognitive performance, in particular slowing of processing speed. We further determined test-retest reliability of functional network measures in sporadic SVD patients participating in a high-frequency (monthly) serial imaging study (RUN DMC-InTENse, median: 8 MRIs per participant). Functional networks for the whole brain and major subsystems (i.e., default mode network, DMN; fronto-parietal task control network, FPCN; visual network, VN; hand somatosensory-motor network, HSMN) were constructed based on resting-state multi-band functional MRI. In CADASIL, global efficiency (a graph metric capturing network integration) of the DMN was lower in patients with high disease burden (standardized beta = -.44; p [corrected] = .035) and mediated the negative effect of disease burden on processing speed (indirect path: std. beta = -.20, p = .047; direct path: std. beta = -.19, p = .25; total effect: std. beta = -.39, p = .02). The corresponding analyses in sporadic SVD showed no effect. Intraclass correlations in the high-frequency serial MRI dataset of the sporadic SVD patients revealed poor test-retest reliability and analysis of individual variability suggested an influence of age, but not disease burden, on global efficiency. In conclusion, our results suggest that changes in functional connectivity networks mediate the effect of SVD-related brain damage on cognitive deficits. However, limited reliability of functional network measures, possibly due to age-related comorbidities, impedes the analysis in elderly SVD patients.

Optimising network modelling methods for fMRI

A major goal of neuroimaging studies is to develop predictive models to analyze the relationship between whole brain functional connectivity patterns and behavioural traits. However, there is no single widely-accepted standard pipeline for analyzing functional connectivity. The common procedure for designing functional connectivity based predictive models entails three main steps: parcellating the brain, estimating the interaction between defined parcels, and lastly, using these integrated associations between brain parcels as features fed to a classifier for predicting non-imaging variables e.g., behavioural traits, demographics, emotional measures, etc. There are also additional considerations when using correlation-based measures of functional connectivity, resulting in three supplementary steps: utilising Riemannian geometry tangent space parameterization to preserve the geometry of functional connectivity; penalizing the connectivity estimates with shrinkage approaches to handle challenges related to short time-series (and noisy) data; and removing confounding variables from brain-behaviour data. These six steps are contingent on each-other, and to optimise a general framework one should ideally examine these various methods simultaneously. In this paper, we investigated strengths and short-comings, both independently and jointly, of the following measures: parcellation techniques of four kinds (categorized further depending upon number of parcels), five measures of functional connectivity, the decision of staying in the ambient space of connectivity matrices or in tangent space, the choice of applying shrinkage estimators, six alternative techniques for handling confounds and finally four novel classifiers/predictors. For performance evaluation, we have selected two of the largest datasets, UK Biobank and the Human Connectome Project resting state fMRI data, and have run more than 9000 different pipeline variants on a total of ∼14000 individuals to determine the optimum pipeline. For independent performance validation, we have run some best-performing pipeline variants on ABIDE and ACPI datasets (∼1000 subjects) to evaluate the generalisability of proposed network modelling methods.

Image acquisition and quality assurance in the Boston Adolescent Neuroimaging of Depression and Anxiety study

The Connectomes Related to Human Diseases (CRHD) initiative was developed with the Human Connectome Project (HCP) to provide high-resolution, open-access, multi-modal MRI data to better understand the neural correlates of human disease. Here, we present an introduction to a CRHD project, the Boston Adolescent Neuroimaging of Depression and Anxiety (BANDA) study, which is collecting multimodal neuroimaging, clinical, and neuropsychological data from 225 adolescents (ages 14-17), 150 of whom are expected to have a diagnosis of depression and/or anxiety. Our transdiagnostic recruitment approach samples the full spectrum of depressed/anxious symptoms and their comorbidity, consistent with NIMH Research Domain Criteria (RDoC). We focused on an age range that is critical for brain development and for the onset of mental illness. This project sought to harmonize imaging sequences, hardware, and functional tasks with other HCP studies, although some changes were made to canonical HCP methods to accommodate our study population and questions. We present a thorough overview of our imaging sequences, hardware, and scanning protocol. We detail similarities and differences between this study and other HCP studies. We evaluate structural-, diffusion-, and functional-image-quality measures that may be influenced by clinical factors (e.g., disorder, symptomatology). Signal-to-noise and motion estimates from the first 140 adolescents suggest minimal influence of clinical factors on image quality. We anticipate enrollment of an additional 85 participants, most of whom are expected to have a diagnosis of anxiety and/or depression. Clinical and neuropsychological data from the first 140 participants are currently freely available through the National Institute of Mental Health Data Archive (NDA).

Brain networks, dimensionality, and global signal averaging in resting-state fMRI: Hierarchical network structure results in low-dimensional spatiotemporal dynamics

One of the most controversial practices in resting-state fMRI functional connectivity studies is whether or not to regress out the global average brain signal (GS) during artifact removal. Some groups have argued that it is absolutely essential to regress out the GS in order to fully remove head motion, respiration, and other global imaging artifacts. Others have argued that removing the GS distorts the resulting correlation matrices and inappropriately alters the results of group comparisons and relationships to behavior. At the core of this argument is the assessment of dimensionality in terms of the number of brain networks with uncorrelated time series. If the dimensionality is high, then the distortions due to GS removal could be effectively negligible. In the current paper, we examine the dimensionality of resting-state fMRI data using principal component analyses (PCA) and network clustering analyses. In two independent datasets (Set 1: N = 62, Set 2: N = 32), scree plots of the eigenvalues level off at or prior to 10 principal components, with prominent elbows at 3 and 7 components. While network clustering analyses have previously demonstrated that numerous networks can be distinguished with high thresholding of the voxel-wise correlation matrices, lower thresholding reveals a lower-dimensional hierarchical structure, with the first prominent branch at 2 networks (corresponding to the previously described "task-positive"/"task-negative" distinction) and further stable subdivisions at 4, 7 and 17. Since inter-correlated time series within these larger branches do not cancel to zero when averaged, the hierarchical nature of the correlation structure results in low effective dimensionality. Consistent with this, partial correlation analyses revealed that network-specific variance remains present in the GS at each level of the hierarchy, accounting for at least 14-18% of the overall GS variance in each dataset. These results demonstrate that GS regression is expected to remove substantial portions of network-specific brain signals along with artifacts, not simply whole-brain signals corresponding to arousal levels. We highlight alternative means of controlling for residual global artifacts when not removing the GS.

Generalizability of machine learning for classification of schizophrenia based on resting-state functional MRI data

Machine learning has increasingly been applied to classification of schizophrenia in neuroimaging research. However, direct replication studies and studies seeking to investigate generalizability are scarce. To address these issues, we assessed within-site and between-site generalizability of a machine learning classification framework which achieved excellent performance in a previous study using two independent resting-state functional magnetic resonance imaging data sets collected from different sites and scanners. We established within-site generalizability of the classification framework in the main data set using cross-validation. Then, we trained a model in the main data set and investigated between-site generalization in the validated data set using external validation. Finally, recognizing the poor between-site generalization performance, we updated the unsupervised algorithm to investigate if transfer learning using additional unlabeled data were able to improve between-site classification performance. Cross-validation showed that the published classification procedure achieved an accuracy of 0.73 using majority voting across all selected components. External validation found a classification accuracy of 0.55 (not significant) and 0.70 (significant) using the direct and transfer learning procedures, respectively. The failure of direct generalization from one site to another demonstrates the limitation of within-site cross-validation and points toward the need to incorporate efforts to facilitate application of machine learning across multiple data sets. The improvement in performance with transfer learning highlights the importance of taking into account the properties of data when constructing predictive models across samples and sites. Our findings suggest that machine learning classification result based on a single study should be interpreted cautiously.

Effective degrees of freedom of the Pearson's correlation coefficient under autocorrelation

The dependence between pairs of time series is commonly quantified by Pearson's correlation. However, if the time series are themselves dependent (i.e. exhibit temporal autocorrelation), the effective degrees of freedom (EDF) are reduced, the standard error of the sample correlation coefficient is biased, and Fisher's transformation fails to stabilise the variance. Since fMRI time series are notoriously autocorrelated, the issue of biased standard errors - before or after Fisher's transformation - becomes vital in individual-level analysis of resting-state functional connectivity (rsFC) and must be addressed anytime a standardised Z-score is computed. We find that the severity of autocorrelation is highly dependent on spatial characteristics of brain regions, such as the size of regions of interest and the spatial location of those regions. We further show that the available EDF estimators make restrictive assumptions that are not supported by the data, resulting in biased rsFC inferences that lead to distorted topological descriptions of the connectome on the individual level. We propose a practical "xDF" method that accounts not only for distinct autocorrelation in each time series, but instantaneous and lagged cross-correlation. We find the xDF correction varies substantially over node pairs, indicating the limitations of global EDF corrections used previously. In addition to extensive synthetic and real data validations, we investigate the impact of this correction on rsFC measures in data from the Young Adult Human Connectome Project, showing that accounting for autocorrelation dramatically changes fundamental graph theoretical measures relative to no correction.

Classification of schizophrenia by intersubject correlation in functional connectome

Functional connectomes have been suggested as fingerprinting for individual identification. Accordingly, we hypothesized that subjects in the same phenotypic group have similar functional connectome features, which could help to discriminate schizophrenia (SCH) patients from healthy controls (HCs) and from depression patients. To this end, we included resting-state functional magnetic resonance imaging data of SCH, depression patients, and HCs from three centers. We first investigated the characteristics of connectome similarity between individuals, and found higher similarity between subjects belonging to the same group (i.e., SCH-SCH) than different groups (i.e., HC-SCH). These findings suggest that the average connectome within group (termed as group-specific functional connectome [GFC]) may help in individual classification. Consistently, significant accuracy (75-77%) and area under curve (81-86%) were found in discriminating SCH from HC or depression patients by GFC-based leave-one-out cross-validation. Cross-center classification further suggests a good generalizability of the GFC classification. We additionally included normal aging data (255 young and 242 old subjects with different scanning sequences) to show factors could be improved for better classification performance, and the findings emphasized the importance of increasing sample size but not temporal resolution during scanning. In conclusion, our findings suggest that the average functional connectome across subjects contained group-specific biological features and may be helpful in clinical diagnosis for schizophrenia.

Machine learning identifies unaffected first-degree relatives with functional network patterns and cognitive impairment similar to those of schizophrenia patients

Schizophrenia (SCZ) patients and their unaffected first-degree relatives (FDRs) share similar functional neuroanatomy. However, it remains largely unknown to what extent unaffected FDRs with functional neuroanatomy patterns similar to patients can be identified at an individual level. In this study, we used a multivariate pattern classification method to learn informative large-scale functional networks (FNs) and build classifiers to distinguish 32 patients from 30 healthy controls and to classify 34 FDRs as with or without FNs similar to patients. Four informative FNs-the cerebellum, default mode network (DMN), ventral frontotemporal network, and posterior DMN with parahippocampal gyrus-were identified based on a training cohort and pattern classifiers built upon these FNs achieved a correct classification rate of 83.9% (sensitivity 87.5%, specificity 80.0%, and area under the receiver operating characteristic curve [AUC] 0.914) estimated based on leave-one-out cross-validation for the training cohort and a correct classification rate of 77.5% (sensitivity 72.5%, specificity 82.5%, and AUC 0.811) for an independent validation cohort. The classification scores of the FDRs and patients were negatively correlated with their measures of cognitive function. FDRs identified by the classifiers as having SCZ patterns were similar to the patients, but significantly different from the controls and FDRs with normal patterns in terms of their cognitive measures. These results demonstrate that the pattern classifiers built upon the informative FNs can serve as biomarkers for quantifying brain alterations in SCZ and help to identify FDRs with FN patterns and cognitive impairment similar to those of SCZ patients.

General functional connectivity: Shared features of resting-state and task fMRI drive reliable and heritable individual differences in functional brain networks

Intrinsic connectivity, measured using resting-state fMRI, has emerged as a fundamental tool in the study of the human brain. However, due to practical limitations, many studies do not collect enough resting-state data to generate reliable measures of intrinsic connectivity necessary for studying individual differences. Here we present general functional connectivity (GFC) as a method for leveraging shared features across resting-state and task fMRI and demonstrate in the Human Connectome Project and the Dunedin Study that GFC offers better test-retest reliability than intrinsic connectivity estimated from the same amount of resting-state data alone. Furthermore, at equivalent scan lengths, GFC displayed higher estimates of heritability than resting-state functional connectivity. We also found that predictions of cognitive ability from GFC generalized across datasets, performing as well or better than resting-state or task data alone. Collectively, our work suggests that GFC can improve the reliability of intrinsic connectivity estimates in existing datasets and, subsequently, the opportunity to identify meaningful correlates of individual differences in behavior. Given that task and resting-state data are often collected together, many researchers can immediately derive more reliable measures of intrinsic connectivity through the adoption of GFC rather than solely using resting-state data. Moreover, by better capturing heritable variation in intrinsic connectivity, GFC represents a novel endophenotype with broad applications in clinical neuroscience and biomarker discovery.

On the analysis of rapidly sampled fMRI data

Recent advances in parallel imaging and simultaneous multi-slice techniques have permitted whole-brain fMRI acquisitions at sub-second sampling intervals, without significantly sacrificing the spatial coverage and resolution. Apart from probing brain function at finer temporal scales, faster sampling rates may potentially lead to enhanced functional sensitivity, owing possibly to both cleaner neural representations (due to less aliased physiological noise) and additional statistical benefits (due to more degrees of freedom for a fixed scan duration). Accompanying these intriguing aspects of fast acquisitions, however, confusion has also arisen regarding (1) how to preprocess/analyze these fast fMRI data, and (2) what exactly is the extent of benefits with fast acquisitions, i.e., how fast is fast enough for a specific research aim? The first question is motivated by the altered spectral distribution and noise characteristics at short sampling intervals, while the second question seeks to reconcile the complicated trade-offs between the functional contrast-to-noise ratio and the effective degrees of freedom. Although there have been recent efforts to empirically approach different aspects of these two questions, in this work we discuss, from a theoretical perspective accompanied by some illustrative, proof-of-concept experimental in vivo human fMRI data, a few considerations that are rarely mentioned, yet are important for both preprocessing and optimizing statistical inferences for studies that employ acquisitions with sub-second sampling intervals. Several summary recommendations include concerns regarding advisability of relying on low-pass filtering to de-noise physiological contributions, employment of statistical models with sufficient complexity to account for the substantially increased serial correlation, and cautions regarding using rapid sampling to enhance functional sensitivity given that different analysis models may associate with distinct trade-offs between contrast-to-noise ratios and the effective degrees of freedom. As an example, we demonstrate that as TR shortens, the intrinsic differences in how noise is accommodated in general linear models and Pearson correlation analyses (assuming Gaussian distributed stochastic signals and noise) can result in quite different outcomes, either gaining or losing statistical power.

Efficacy of different dynamic functional connectivity methods to capture cognitively relevant information

Given the dynamic nature of the human brain, there has been an increasing interest in investigating short-term temporal changes in functional connectivity, also known as dynamic functional connectivity (dFC), i.e., the time-varying inter-regional statistical dependence of blood oxygenation level-dependent (BOLD) signal within the constraints of a single scan. Numerous methodologies have been proposed to characterize dFC during rest and task, but few studies have compared them in terms of their efficacy to capture behavioral and clinically relevant dynamics. This is mostly due to lack of a well-defined ground truth, especially for rest scans. In this study, with a multitask dataset (rest, memory, video, and math) serving as ground truth, we investigated the efficacy of several dFC estimation techniques at capturing cognitively relevant dFC modulation induced by external tasks. We evaluated two framewise methods (dFC estimates for a single time point): dynamic conditional correlation (DCC) and jackknife correlation (JC); and five window-based methods: sliding window correlation (SWC), sliding window correlation with L1-regularization (SWC_L1), a combination of DCC and SWC called moving average DCC (DCC_MA), multiplication of temporal derivatives (MTD), and a variant of jackknife correlation called delete-d jackknife correlation (dJC). The efficacy is defined as each dFC metric's ability to successfully subdivide multitask scans into cognitively homogenous segments (even if those segments are not temporally continuous). We found that all window-based dFC methods performed well for commonly used window lengths (WL ≥ 30sec), with sliding window methods (SWC, SWC_L1) as well as the hybrid DCC_MA approach performing slightly better. For shorter window lengths (WL ≤ 15sec), DCC_MA and dJC produced the best results. Neither framewise method (i.e., DCC and JC) led to dFC estimates with high accuracy.

Tracking and validation techniques for topographically organized tractography

Topographic regularity of axonal connections is commonly understood as the preservation of spatial relationships between nearby neurons and is a fundamental structural property of the brain. In particular the retinotopic mapping of the visual pathway can even be quantitatively computed. Inspired from this previously untapped anatomical knowledge, we propose a novel tractography method that preserves both topographic and geometric regularity. We make use of parameterized curves with Frenet-Serret frame and introduce a highly flexible mechanism for controlling geometric regularity. At the same time, we incorporate a novel local data support term in order to account for topographic organization. Unifying geometry with topographic regularity, we develop a Bayesian framework for generating highly organized streamlines that accurately follow neuroanatomy. We additionally propose two novel validation techniques to quantify topographic regularity. In our experiments, we studied the results of our approach with respect to connectivity, reproducibility and topographic regularity aspects. We present both qualitative and quantitative comparisons of our technique against three algorithms from MRtrix3. We show that our method successfully generates highly organized fiber tracks while capturing bundle anatomy that are geometrically challenging for other approaches.

Sparse Graphical Models for Functional Connectivity Networks: Best Methods and the Autocorrelation Issue

Sparse graphical models are frequently used to explore both static and dynamic functional brain networks from neuroimaging data. However, the practical performance of the models has not been studied in detail for brain networks. In this work, we have two objectives. First, we compare several sparse graphical model estimation procedures and several selection criteria under various experimental settings, such as different dimensions, sample sizes, types of data, and sparsity levels of the true model structures. We discuss in detail the superiority and deficiency of each combination. Second, in the same simulation study, we show the impact of autocorrelation and whitening on the estimation of functional brain networks. We apply the methods to a resting-state functional magnetic resonance imaging (fMRI) data set. Our results show that the best sparse graphical model, in terms of detection of true connections and having few false-positive connections, is the smoothly clipped absolute deviation (SCAD) estimating method in combination with the Bayesian information criterion (BIC) and cross-validation (CV) selection method. In addition, the presence of autocorrelation in the data adversely affects the estimation of networks but can be helped by using the CV selection method. These results question the validity of a number of fMRI studies where inferior graphical model techniques have been used to estimate brain networks.

On the Reliability of Individual Brain Activity Networks

There is intense interest in fMRI research on whole-brain functional connectivity, and however, two fundamental issues are still unresolved: the impact of spatiotemporal data resolution (spatial parcellation and temporal sampling) and the impact of the network construction method on the reliability of functional brain networks. In particular, the impact of spatiotemporal data resolution on the resulting connectivity findings has not been sufficiently investigated. In fact, a number of studies have already observed that functional networks often give different conclusions across different parcellation scales. If the interpretations from functional networks are inconsistent across spatiotemporal scales, then the whole validity of the functional network paradigm is called into question. This paper investigates the consistency of resting state network structure when using different temporal sampling or spatial parcellation, or different methods for constructing the networks. To pursue this, we develop a novel network comparison framework based on persistent homology from a topological data analysis. We use the new network comparison tools to characterize the spatial and temporal scales under which consistent functional networks can be constructed. The methods are illustrated on Human Connectome Project data, showing that the DISCOH₂ network construction method outperforms other approaches at most data spatiotemporal resolutions.

Reliability correction for functional connectivity: Theory and implementation

Network properties can be estimated using functional connectivity MRI (fcMRI). However, regional variation of the fMRI signal causes systematic biases in network estimates including correlation attenuation in regions of low measurement reliability. Here we computed the spatial distribution of fcMRI reliability using longitudinal fcMRI datasets and demonstrated how pre-estimated reliability maps can correct for correlation attenuation. As a test case of reliability-based attenuation correction we estimated properties of the default network, where reliability was significantly lower than average in the medial temporal lobe and higher in the posterior medial cortex, heterogeneity that impacts estimation of the network. Accounting for this bias using attenuation correction revealed that the medial temporal lobe's contribution to the default network is typically underestimated. To render this approach useful to a greater number of datasets, we demonstrate that test-retest reliability maps derived from repeated runs within a single scanning session can be used as a surrogate for multi-session reliability mapping. Using data segments with different scan lengths between 1 and 30 min, we found that test-retest reliability of connectivity estimates increases with scan length while the spatial distribution of reliability is relatively stable even at short scan lengths. Finally, analyses of tertiary data revealed that reliability distribution is influenced by age, neuropsychiatric status and scanner type, suggesting that reliability correction may be especially important when studying between-group differences. Collectively, these results illustrate that reliability-based attenuation correction is an easily implemented strategy that mitigates certain features of fMRI signal nonuniformity.

Reproducibility and robustness of graph measures of the associative-semantic network

Graph analysis is a promising tool to quantify brain connectivity. However, an essential requirement is that the graph measures are reproducible and robust. We have studied the reproducibility and robustness of various graph measures in group based and in individual binary and weighted networks derived from a task fMRI experiment during explicit associative-semantic processing of words and pictures. The nodes of the network were defined using an independent study and the connectivity was based on the partial correlation of the time series between any pair of nodes. The results showed that in case of binary networks, global graph measures exhibit a good reproducibility and robustness for networks which are not too sparse and these figures of merit depend on the graph measure and on the density of the network. Furthermore, group based binary networks should be derived from groups of sufficient size and the lower the density the more subjects are required to obtain robust values. Local graph measures are very variable in terms of reproducibility and should be interpreted with care. For weighted networks, we found good reproducibility (average test-retest variability <5% and ICC values >0.4) when using subject specific networks and this will allow us to relate network properties to individual subject information.

Multimodal imaging in autism: an early review of comprehensive neural circuit characterization

There is accumulating evidence that the neurobiology of autism spectrum disorders (ASD) is linked to atypical neural communication and connectivity. This body of work emphasizes the need to characterize the function of multiple regions that comprise neural circuits rather than focusing on singular regions as contributing to deficits in ASD. Multimodal neuroimaging - the formal combination of multiple functional and structural measures of the brain - is extremely promising as an approach to understanding neural deficits in ASD. This review provides an overview of the multimodal imaging approach, and then provides a snapshot of how multimodal imaging has been applied in the study of ASD to date. This body of work is separated into two categories: one concerning whole brain connectomics and the other focused on characterizing neural circuits implicated as altered in ASD. We end this review by highlighting emerging themes from the existing body of literature, and new resources that will likely influence future multimodal imaging studies.

Human Connectome Project informatics: quality control, database services, and data visualization

The Human Connectome Project (HCP) has developed protocols, standard operating and quality control procedures, and a suite of informatics tools to enable high throughput data collection, data sharing, automated data processing and analysis, and data mining and visualization. Quality control procedures include methods to maintain data collection consistency over time, to measure head motion, and to establish quantitative modality-specific overall quality assessments. Database services developed as customizations of the XNAT imaging informatics platform support both internal daily operations and open access data sharing. The Connectome Workbench visualization environment enables user interaction with HCP data and is increasingly integrated with the HCP's database services. Here we describe the current state of these procedures and tools and their application in the ongoing HCP study.