A decade of test-retest reliability of functional connectivity: A systematic review and meta-analysis

Statistical power in network neuroscience

Network neuroscience has emerged as a leading method to study brain connectivity. The success of these investigations is dependent not only on approaches to accurately map connectivity but also on the ability to detect real effects in the data - that is, statistical power. We review the state of statistical power in the field and discuss sample size, effect size, measurement error, and network topology as key factors that influence the power of brain connectivity investigations. We use the term 'differential power' to describe how power can vary between nodes, edges, and graph metrics, leaving traces in both positive and negative connectome findings. We conclude with strategies for working with, rather than around, power in connectivity studies.

Uncovering individual differences in fine-scale dynamics of functional connectivity

Functional connectivity (FC) profiles contain subject-specific features that are conserved across time and have potential to capture brain-behavior relationships. Most prior work has focused on spatial features (nodes and systems) of these FC fingerprints, computed over entire imaging sessions. We propose a method for temporally filtering FC, which allows selecting specific moments in time while also maintaining the spatial pattern of node-based activity. To this end, we leverage a recently proposed decomposition of FC into edge time series (eTS). We systematically analyze functional magnetic resonance imaging frames to define features that enhance identifiability across multiple fingerprinting metrics, similarity metrics, and data sets. Results show that these metrics characteristically vary with eTS cofluctuation amplitude, similarity of frames within a run, transition velocity, and expression of functional systems. We further show that data-driven optimization of features that maximize fingerprinting metrics isolates multiple spatial patterns of system expression at specific moments in time. Selecting just 10% of the data can yield stronger fingerprints than are obtained from the full data set. Our findings support the idea that FC fingerprints are differentially expressed across time and suggest that multiple distinct fingerprints can be identified when spatial and temporal characteristics are considered simultaneously.

Reliability of Decision-Making and Reinforcement Learning Computational Parameters

Computational models can offer mechanistic insight into cognition and therefore have the potential to transform our understanding of psychiatric disorders and their treatment. For translational efforts to be successful, it is imperative that computational measures capture individual characteristics reliably. Here we examine the reliability of reinforcement learning and economic models derived from two commonly used tasks. Healthy individuals (N = 50) completed a restless four-armed bandit and a calibrated gambling task twice, two weeks apart. Reward and punishment learning rates from the reinforcement learning model showed good reliability and reward and punishment sensitivity from the same model had fair reliability; while risk aversion and loss aversion parameters from a prospect theory model exhibited good and excellent reliability, respectively. Both models were further able to predict future behaviour above chance within individuals. This prediction was better when based on participants' own model parameters than other participants' parameter estimates. These results suggest that reinforcement learning, and particularly prospect theory parameters, as derived from a restless four-armed bandit and a calibrated gambling task, can be measured reliably to assess learning and decision-making mechanisms. Overall, these findings indicate the translational potential of clinically-relevant computational parameters for precision psychiatry.

Are inter-subject functional correlations consistent across different movies?

Movie fMRI has been increasingly used in investigations of human brain function. Inter-subject functional correlation (ISFC), which evaluates stimulus-dependent inter-regional synchrony between brains exposed to the same stimulus, is emerging as an influencing measure for movie fMRI data analyses. Before the wide application of ISFC analyses, it will be useful to investigate the degree to which they are similar and different across different movies. Based on the four movie fMRI runs of 178 subjects included in the "human connectome project (HCP) S1200 Release", we evaluated ISFCs throughout the brain and analyzed their consistency across different movies using intra-class correlation (ICC). We also investigated the generalizability of ISFC-based predictive models, which is closely related to their consistency, with sex classification and grip strength prediction used as test cases. The results showed that the intensity of ISFCs was generally weak (0.047). Except a few within-network ones (e.g., ICC of ISFC in the PON was 0.402), ISFCs throughout the brain exhibited low consistency, as indicated by a mean ICC of 0.130. The accuracies for inter-run predictions (60.7-72.8% for sex classification, and R = 0.122-0.275 for grip strength prediction) were much lower than those for intra-run predictions (73.2-83.0% for sex classification, and R = 0.325-0.403 for grip strength prediction), and this indicates poor generalizability of predictive models based on ISFCs. According to these findings, ISFC analyses capture aspects of brain function that are specific to each individual movie, and this specificity should be taken into account (in some cases might be especially useful) in future naturalistic studies.

Differences in intrinsic functional networks in patients with essential tremor who had good and poor long-term responses after thalamotomy performed using MR-guided ultrasound

OBJECTIVE

Thalamotomy at the nucleus ventralis intermedius using MR-guided focused ultrasound has been an effective treatment method for essential tremor (ET). However, this is not true for all cases, even for successful ablation. How the brain differs in patients with ET between those with long-term good and poor outcomes is not clear. To analyze the functional connectivity difference between patients in whom thalamotomy was effective and those in whom thalamotomy was ineffective and its prognostic role in ET treatment, the authors evaluated preoperative resting-state functional MRI in thalamotomy-treated patients.

METHODS

Preoperative resting-state functional MRI data in 85 patients with ET, who were experiencing tremor relief at the time of treatment and were followed up for a minimum of 6 months after the procedure, were collected for the study. The authors conducted a graph independent component analysis of the functional connectivity matrices of tremor-related networks. The patients were divided into thalamotomy-effective and thalamotomy-ineffective groups (thalamotomy-effective group, ≥ 50% motor symptom reduction; thalamotomy-ineffective group, < 50% motor symptom reduction at 6 months after treatment) and the authors compared network components between groups.

RESULTS

Seventy-two (84.7%) of the 85 patients showed ≥ 50% tremor reduction from baseline at 6 months after thalamotomy. The network analysis shows significant suppression of functional network components with connections between the areas of the cerebellum and the basal ganglia and thalamus, but enhancement of those between the premotor cortex and supplementary motor area in the noneffective group compared to the effective group.

CONCLUSIONS

The present study demonstrates that patients in the noneffective group have suppressed functional subnetworks in the cerebellum and subcortex regions and have enhanced functional subnetworks among motor-sensory cortical networks compared to the thalamotomy-effective group. Therefore, the authors suggest that the functional connectivity pattern might be a possible predictive factor for outcomes of MR-guided focused ultrasound thalamotomy.

Functional Connectivity Between Basal Forebrain and Superficial Amygdala Negatively Correlates with Social Fearfulness

Social anxiety is characterized by an intense fear of evaluation from others and/or withdrawal from social situations. Extreme social anxiety can lead to social anxiety disorder. There remains an urgent need to investigate the neural substrates of subclinical social anxiety for early diagnosis and intervention to reduce the risk to develop social anxiety disorder. Twenty-nine young adults were recruited (10 males/19 females; mean age (SD) = 20.34 (2.29)). Trait-like social anxiety was assessed by Liebowitz Social Anxiety Scale. Functional magnetic resonance imaging was used with an emotional face-matching paradigm to probe brain activation in response to emotional stimuli including angry, fearful, and happy faces, with shape-matching as a control condition. Behavioral results showed positive correlations between Liebowitz Social Anxiety Scale scores and the reaction time in both angry and fearful conditions. The activation of superficial amygdala and the deactivation of basal forebrain in response to angry condition showed positive correlations with the level of social anxiety. In addition, the resting-state functional connectivity between these two regions was negatively correlated with the level of social anxiety. These results may help to understand the individual difference and corresponding neural underpinnings of social anxiety in the subclinical population, and might provide some insight to develop strategies for early diagnosis and interventions of social anxiety to reduce the risk of deterioration from subclinical to clinical level of social anxiety.

Automatic Preprocessing Pipeline for White Matter Functional Analyses of Large-Scale Databases

Recently, increasing evidence suggests that fMRI signals in white matter (WM), conventionally ignored as nuisance, are robustly detectable using appropriate processing methods and are related to neural activity, while changes in WM with aging and degeneration are also well documented. These findings suggest variations in patterns of BOLD signals in WM should be investigated. However, existing fMRI analysis tools, which were designed for processing gray matter signals, are not well suited for large-scale processing of WM signals in fMRI data. We developed an automatic pipeline for high-performance preprocessing of fMRI images with emphasis on quantifying changes in BOLD signals in WM in an aging population. At the image processing level, the pipeline integrated existing software modules with fine parameter tunings and modifications to better extract weaker WM signals. The preprocessing results primarily included whole-brain time-courses, functional connectivity, maps and tissue masks in a common space. At the job execution level, this pipeline exploited a local XNAT to store datasets and results, while using DAX tool to automatic distribute batch jobs that run on high-performance computing clusters. Through the pipeline, 5,034 fMRI/T1 scans were preprocessed. The intraclass correlation coefficient (ICC) of test-retest experiment based on the preprocessed data is 0.52 - 0.86 (N=1000), indicating a high reliability of our pipeline, comparable to previously reported ICC in gray matter experiments. This preprocessing pipeline highly facilitates our future analyses on WM functional alterations in aging and may be of benefit to a larger community interested in WM fMRI studies.

Within-subject reliability of concurrent TMS-fMRI during a single session

Concurrent transcranial magnetic stimulation with functional MRI (concurrent TMS-fMRI) allows real-time causative probing of brain connectivity. However, technical challenges, safety, and tolerability may limit the number of trials employed during a concurrent TMS-fMRI experiment. We leveraged an existing data set with 100 trials of active TMS compared to a sub-threshold control condition to assess the reliability of the evoked BOLD response during concurrent TMS-fMRI. This data will permit an analysis of the minimum number of trials that should be employed in a concurrent TMS-fMRI protocol in order to achieve reliable spatial changes in activity. Single-subject maps of brain activity were created by splitting the trials within the same experimental session into groups of 50, 40, 30, 25, 20, 15, or 10 trials, correlations (R) between t-maps derived from paired subsets of trials within the same individual were calculated as reliability. R was moderate-high for 50 trials (mean R = .695) and decreased as the number of trials decreased. Consistent with previous findings of high individual variability in the spatial patterns of evoked neuronal changes following a TMS pulse, the spatial pattern of Rs differed across participants, but regional R was correlated with the magnitude of TMS-evoked activity. These results demonstrate concurrent TMS-fMRI produces a reliable pattern of activity at the individual level at higher trial numbers, particularly within localized regions. The spatial pattern of reliability is individually idiosyncratic and related to the individual pattern of evoked changes.

The functional connectivity between left insula and left medial superior frontal gyrus underlying the relationship between rumination and procrastination

Procrastination is regarded as a prevalent problematic behavior that impairs people's physical and mental health. Although previous studies have indicated that trait rumination is robustly positively correlated with procrastination, it remains unknown about the neural substrates underlying the relationship between trait rumination and procrastination. To address this issue, we used voxel-based morphometry (VBM) and resting-state functional connectivity (RSFC) approaches to explore the neural basis of the relationship between trait rumination and procrastination. Our behavior results found that trait rumination was significantly positively correlated to procrastination, while the VBM analysis showed that trait rumination was negatively correlated with gray matter volume of the insula. Furthermore, the RSFC results revealed a negative association of the left insula-lmSFG (left medial superior frontal gyrus) functional connectivity with trait rumination. More importantly, the mediation analysis showed that trait rumination could completely mediate the relationship between left insula-lmSFG functional connectivity and procrastination. These results suggest that the left insula-lmSFG functional connectivity involved in emotion regulation modulates the association between trait rumination and procrastination, which provides neural evidence for the relationship between trait rumination and procrastination.

Executive control network resting state fMRI functional and effective connectivity and delay discounting in cocaine dependent subjects compared to healthy controls

Resting state functional magnetic resonance imaging (fMRI) has been used to study functional connectivity of brain networks in addictions. However, most studies to-date have focused on the default mode network (DMN) with fewer studies assessing the executive control network (ECN) and salience network (SN), despite well-documented cognitive executive behavioral deficits in addictions. The present study assessed the functional and effective connectivity of the ECN, DMN, and SN in cocaine dependent subjects (CD) (n = 22) compared to healthy control subjects (HC) (n = 22) matched on age and education. This study also investigated the relationship between impulsivity measured by delay discounting and functional and effective connectivity of the ECN, DMN, and SN. The Left ECN (LECN), Right ECN (RECN), DMN, and SN functional networks were identified using FSL MELODIC independent component analysis. Functional connectivity differences between CD and HC were assessed using FSL Dual Regression analysis and FSLNets. Effective connectivity differences between CD and HC were measured using the Parametric Empirical Bayes module of Dynamic Causal Modeling. The relationship between delay discounting and functional and effective connectivity were examined using regression analyses. Dynamic causal modeling (DCM) analysis showed strong evidence (posterior probability > 0.95) for CD to have greater effective connectivity than HC in the RECN to LECN pathway when tobacco use was included as a factor in the model. DCM analysis showed strong evidence for a positive association between delay discounting and effective connectivity for the RECN to LECN pathway and for the DMN to DMN self-connection. There was strong evidence for a negative association between delay discounting and effective connectivity for the DMN to RECN pathway and for the SN to DMN pathway. Results also showed strong evidence for a negative association between delay discounting and effective connectivity for the RECN to SN pathway in CD but a positive association in HC. These novel findings provide preliminary support that RECN effective connectivity may differ between CD and HC after controlling for tobacco use. RECN effective connectivity may also relate to tobacco use and impulsivity as measured by delay discounting.

Does the current state of biomarker discovery in autism reflect the limits of reductionism in precision medicine? Suggestions for an integrative approach that considers dynamic mechanisms between brain, body, and the social environment

Over the past decade, precision medicine has become one of the most influential approaches in biomedical research to improve early detection, diagnosis, and prognosis of clinical conditions and develop mechanism-based therapies tailored to individual characteristics using biomarkers. This perspective article first reviews the origins and concept of precision medicine approaches to autism and summarises recent findings from the first "generation" of biomarker studies. Multi-disciplinary research initiatives created substantially larger, comprehensively characterised cohorts, shifted the focus from group-comparisons to individual variability and subgroups, increased methodological rigour and advanced analytic innovations. However, although several candidate markers with probabilistic value have been identified, separate efforts to divide autism by molecular, brain structural/functional or cognitive markers have not identified a validated diagnostic subgroup. Conversely, studies of specific monogenic subgroups revealed substantial variability in biology and behaviour. The second part discusses both conceptual and methodological factors in these findings. It is argued that the predominant reductionist approach, which seeks to parse complex issues into simpler, more tractable units, let us to neglect the interactions between brain and body, and divorce individuals from their social environment. The third part draws on insights from systems biology, developmental psychology and neurodiversity approaches to outline an integrative approach that considers the dynamic interaction between biological (brain, body) and social mechanisms (stress, stigma) to understanding the origins of autistic features in particular conditions and contexts. This requires 1) closer collaboration with autistic people to increase face validity of concepts and methodologies; (2) development of measures/technologies that enable repeat assessment of social and biological factors in different (naturalistic) conditions and contexts, (3) new analytic methods to study (simulate) these interactions (including emergent properties), and (4) cross-condition designs to understand which mechanisms are transdiagnostic or specific for particular autistic sub-populations. Tailored support may entail both creating more favourable conditions in the social environment and interventions for some autistic people to increase well-being.

Dwelling in prolonged grief: Resting state functional connectivity during oxytocin and placebo administration

Clinical theories of adaptation in bereavement highlight a need for flexible shifting between mental states. However, prolonged motivational salience of the deceased partner may be a complicating factor, particularly when coupled with perseverative thinking about the loss. We investigated how prolonged grief symptoms might relate to resting state functional brain network connectivity in a sample of older adults (n = 38) who experienced the death of a partner 6-36 months prior, and whether intranasal oxytocin (as a neuropeptide involved in pair-bonding) had differential effects in participants with higher prolonged grief symptoms. Higher scores on the Inventory of Complicated Grief (ICG) were associated with lower anticorrelation (i.e., higher functional connectivity) between the defaultretrosplenial - cingulo-operculardACC network pair. Intranasal oxytocin increased functional connectivity in the same defaultretrosplenial - cingulo-operculardACC circuit but ICG scores did not moderate effects of oxytocin, contrary to our prediction. Higher ICG scores were associated with longer dwell time in a dynamic functional connectivity state featuring positive correlations among default, frontoparietal, and cingulo-opercular networks, across both placebo and oxytocin sessions. Dwell time was not significantly affected by oxytocin, and higher prolonged grief symptoms were not associated with more variability in dynamic functional connectivity states over the scan. Results offer preliminary evidence that prolonged grief symptoms in older adults are associated with patterns of static and time-varying functional network connectivity and may specifically involve a default network-salience-related circuit that is sensitive to oxytocin.

Assessing the Repeatability of Multi-Frequency Multi-Layer Brain Network Topologies Across Alternative Researcher's Choice Paths

There is a growing interest in the neuroscience community on the advantages of multilayer functional brain networks. Researchers usually treated different frequencies separately at distinct functional brain networks. However, there is strong evidence that these networks share complementary information while their interdependencies could reveal novel findings. For this purpose, neuroscientists adopt multilayer networks, which can be described mathematically as an extension of trivial single-layer networks. Multilayer networks have become popular in neuroscience due to their advantage to integrate different sources of information. Here, Ι will focus on the multi-frequency multilayer functional connectivity analysis on resting-state fMRI (rs-fMRI) recordings. However, constructing a multilayer network depends on selecting multiple pre-processing steps that can affect the final network topology. Here, I analyzed the rs-fMRI dataset from a single human performing scanning over a period of 18 months (84 scans in total), and the rs-fMRI dataset containing 25 subjects with 3 repeat scans. I focused on assessing the reproducibility of multi-frequency multilayer topologies exploring the effect of two filtering methods for extracting frequencies from BOLD activity, three connectivity estimators, with or without a topological filtering scheme, and two spatial scales. Finally, I untangled specific combinations of researchers' choices that yield consistently brain networks with repeatable topologies, giving me the chance to recommend best practices over consistent topologies.

Functional underpinnings of feedback-enhanced test-potentiated encoding

The testing effect describes the finding that retrieval practice enhances memory performance compared to restudy practice. Prior evidence demonstrates that this effect can be boosted by providing feedback after retrieval attempts (i.e. test-potentiated encoding [TPE]). The present fMRI study investigated the neural processes during successful memory retrieval underlying this beneficial effect of correct answer feedback compared with restudy and whether additional performance feedback leads to further benefits. Twenty-seven participants learned cue-target pairs by (i) restudying, (ii) standard TPE including a restudy opportunity, or (iii) TPE including a restudy opportunity immediately after a positive or negative performance feedback. One day later, a cued retrieval recognition test was performed inside the MRI scanner. Behavioral results confirmed the testing effect and that adding explicit performance feedback-enhanced memory relative to restudy and standard TPE. Stronger functional engagement while retrieving items previously restudied was found in lateral prefrontal cortex and superior parietal lobe. By contrast, lateral temporo-parietal areas were more strongly recruited while retrieving items previously tested. Performance feedback increased the hippocampal activation and resulted in stronger functional coupling between hippocampus, supramarginal gyrus, and ventral striatum with lateral temporo-parietal cortex. Our results unveil the main functional dynamics and connectivity nodes underlying memory benefits from additional performance feedback.

Patterns of connectome variability in autism across five functional activation tasks: findings from the LEAP project

BACKGROUND

Autism spectrum disorder (autism) is a complex neurodevelopmental condition with pronounced behavioral, cognitive, and neural heterogeneities across individuals. Here, our goal was to characterize heterogeneity in autism by identifying patterns of neural diversity as reflected in BOLD fMRI in the way individuals with autism engage with a varied array of cognitive tasks.

METHODS

All analyses were based on the EU-AIMS/AIMS-2-TRIALS multisite Longitudinal European Autism Project (LEAP) with participants with autism (n = 282) and typically developing (TD) controls (n = 221) between 6 and 30 years of age. We employed a novel task potency approach which combines the unique aspects of both resting state fMRI and task-fMRI to quantify task-induced variations in the functional connectome. Normative modelling was used to map atypicality of features on an individual basis with respect to their distribution in neurotypical control participants. We applied robust out-of-sample canonical correlation analysis (CCA) to relate connectome data to behavioral data.

RESULTS

Deviation from the normative ranges of global functional connectivity was greater for individuals with autism compared to TD in each fMRI task paradigm (all tasks p < 0.001). The similarity across individuals of the deviation pattern was significantly increased in autistic relative to TD individuals (p < 0.002). The CCA identified significant and robust brain-behavior covariation between functional connectivity atypicality and autism-related behavioral features.

CONCLUSIONS

Individuals with autism engage with tasks in a globally atypical way, but the particular spatial pattern of this atypicality is nevertheless similar across tasks. Atypicalities in the tasks originate mostly from prefrontal cortex and default mode network regions, but also speech and auditory networks. We show how sophisticated modeling methods such as task potency and normative modeling can be used toward unravelling complex heterogeneous conditions like autism.

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

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

A graph theory neuroimaging approach to distinguish the depression of bipolar disorder from major depressive disorder in adolescents and young adults

BACKGROUND

Markers to differentiate depressions of bipolar disorder (BD-Dep) from depressions of major depressive disorder (MDD-Dep), and for more targeted treatments, are critically needed to decrease current high rates of misdiagnosis that can lead to ineffective or potentially deleterious treatments. Distinguishing, and specifically treating the depressions, during the adolescent/young adult epoch is especially important to decrease illness progression and improve prognosis, and suicide, as it is the epoch when suicide thoughts and behaviors often emerge. With differences in functional connectivity patterns reported when BD-Dep and MDD-Dep have been studied separately, this study used a graph theory approach aimed to identify functional connectivity differences in their direct comparison.

METHODS

Functional magnetic resonance imaging whole-brain functional connectivity (Intrinsic Connectivity Distribution, ICD) measures were compared across adolescents/young adults with BD-Dep (n = 28), MDD-Dep (n = 20) and HC (n = 111). Follow-up seed-based connectivity was conducted on regions of significant ICD differences. Relationships with demographic and clinical measures were assessed.

RESULTS

Compared to the HC group, both the BD-Dep and MDD-Dep groups exhibited left-sided frontal, insular, and medial temporal ICD increases. The BD-Dep group had additional right-sided ICD increases in frontal, basal ganglia, and fusiform areas. In seed-based analyses, the BD-Dep group exhibited increased interhemispheric functional connectivity between frontal areas not seen in the MDD-Dep group.

LIMITATIONS

Modest sample size; medications not studied systematically.

CONCLUSIONS

This study supports bilateral and interhemispheric functional dysconnectivity as features of BD-Dep that may differentiate it from MDD-Dep in adolescents/young adults and serve as a target for early diagnosis and treatment strategies.

Test-retest reliability of modular-relevant analysis in brain functional network

Introduction

The human brain could be modeled as a complex network via functional magnetic resonance imaging (fMRI), and the architecture of these brain functional networks can be studied from multiple spatial scales with different graph theory tools. Detecting modules is an important mesoscale network measuring approach that has provided crucial insights for uncovering how brain organizes itself among different functional subsystems. Despite its successful application in a wide range of brain network studies, the lack of comprehensive reliability assessment prevents its potential extension to clinical trials.

Methods

To fill this gap, this paper, using resting-state test-retest fMRI data, systematically explored the reliabilities of five popular network metrics derived from modular structure. Considering the repeatability of network partition depends heavily on network size and module detection algorithm, we constructed three types of brain functional networks for each subject by using a set of coarse-to-fine brain atlases and adopted four methods for single-subject module detection and twelve methods for group-level module detection.

Results

The results reported moderate-to-good reliability in modularity, intra- and inter-modular functional connectivities, within-modular degree and participation coefficient at both individual and group levels, indicating modular-relevant network metrics can provide robust evaluation results. Further analysis identified the significant influence of module detection algorithm and node definition approach on reliabilities of network partitions and its derived network analysis results.

Discussion

This paper provides important guidance for choosing reliable modular-relevant network metrics and analysis strategies in future studies.

Test-retest reliability of prepulse inhibition (PPI) and PPI correlation with working memory

OBJECTIVE

Sensorimotor gating is experimentally operationalized by the prepulse inhibition (PPI) of the startle response (SR). Previous studies suggest high test-retest reliability of PPI and potential correlation with working memory (WM). Here, we aimed to validate and extend the test-retest reliability of PPI in healthy humans and its correlation with WM performance.

METHODS

We applied an acoustic startle PPI paradigm with four different prepulse intensities (64, 68, 72 and 76 dB) and two different WM tasks [n-back, change detection task (CDT)] in a group of 26 healthy adults (final sample size n = 23). To assess test-retest reliability, we performed all tests on two separate days ~27 days (range: 21-32 days) apart.

RESULTS

We were able to confirm high test-retest reliability of the PPI with a mean intraclass correlation (ICC) of > 0.80 and significant positive correlation of PPI with n-back but not with CDT performance. Detailed analysis showed that PPI across all prepulse intensities significantly correlated with both the 2-back and 0-back conditions, suggesting regulation by cross-conditional processes (e.g. attention). However, when removing the 0-back component from the 2-back data, we found a specific and significant correlation with WM for the 76-dB PPI condition.

CONCLUSION

With the present study, we were able to confirm the high test-retest reliability of the PPI in humans and could validate and expand on its correlation with WM performance.

Cross-cohort replicability and generalizability of connectivity-based psychometric prediction patterns

An increasing number of studies have investigated the relationships between inter-individual variability in brain regions' connectivity and behavioral phenotypes, making use of large population neuroimaging datasets. However, the replicability of brain-behavior associations identified by these approaches remains an open question. In this study, we examined the cross-dataset replicability of brain-behavior association patterns for fluid cognition and openness predictions using a previously developed region-wise approach, as well as using a standard whole-brain approach. Overall, we found moderate similarity in patterns for fluid cognition predictions across cohorts, especially in the Human Connectome Project Young Adult, Human Connectome Project Aging, and Enhanced Nathan Kline Institute Rockland Sample cohorts, but low similarity in patterns for openness predictions. In addition, we assessed the generalizability of prediction models in cross-dataset predictions, by training the model in one dataset and testing in another. Making use of the region-wise prediction approach, we showed that first, a moderate extent of generalizability could be achieved with fluid cognition prediction, and that, second, a set of common brain regions related to fluid cognition across cohorts could be identified. Nevertheless, the moderate replicability and generalizability could only be achieved in specific contexts. Thus, we argue that replicability and generalizability in connectivity-based prediction remain limited and deserve greater attention in future studies.

Striatal resting-state connectivity after long-term diacetylmorphine treatment in opioid-dependent patients

New treatment approaches for opioid-dependent patients include injectable opioid agonist treatment with diacetylmorphine. While evidence has shown beneficial clinical effects of diacetylmorphine, it is still not clear how long-term diacetylmorphine treatment affects the brain and whether functional brain changes are accompanied by clinical improvements. Therefore, this prospective case-control study focuses on long-term effects of diacetylmorphine on resting-state functional connectivity. We included opioid-dependent patients (N = 22, age range 33-58, 16 males) treated with diacetylmorphine and healthy controls (N = 9, age range 27-55, 5 males) that underwent two MRI assessments approximately nine years apart. For the patients, the assessments took part shortly after the diacetylmorphine intake to be able to explore changes in resting-state functional connectivity in brain regions related to the stage of binge and intoxication (caudate, putamen, nucleus accumbens). A cluster in the right superior frontal gyrus was detected, showing over nine years an increase in functional connectivity originating from the left caudate and the left accumbens in patients but not in healthy controls. These connectivity changes in patients were related to the duration of the diacetylmorphine treatment at the follow-up, indicating smaller increases in functional connectivity with longer treatment duration (r = 0.63, P < 0.01). These results suggest that long-term diacetylmorphine treatment in opioid-dependent patients increases fronto-striatal connections, an effect that is linked to the duration of the treatment duration. Future research needs to further address the wide-ranging effects of diacetylmorphine on brain functioning and deepen the understanding of their clinical relevance.

How Repeatable Is the Ergogenic Effect of Caffeine? Limited Reproducibility of Acute Caffeine (3 mg.kg-1) Ingestion on Muscular Strength, Power, and Muscular Endurance

This study aimed to determine the effect of 3 mg.kg−1 acute caffeine ingestion on muscular strength, power and strength endurance and the repeatability of potential ergogenic effects across multiple trials. Twenty-two university standard male rugby union players (20 ± 2 years) completed the study. Using a double-blind, randomized, and counterbalanced within-subject experimental design. Participants completed six experimental trials (three caffeine and three placebo) where force time characteristic of the Isometric Mid-Thigh Pull (IMTP), Countermovement Jump (CMJ) and Drop Jumps (DJ) were assessed followed by assessments of Chest Press (CP), Shoulder Press (SP), Squats (SQ), and Deadlifts (DL) Repetitions Until Failure (RTF at 70% 1 RM). ANOVA indicated that caffeine improved both the CMJ and DJ (p < 0.044) and increased RTF in all RTF assessments (p < 0.002). When individual caffeine trials were compared to corresponding placebo trials, effect sizes ranged from trivial-large favoring caffeine irrespective of a main effect of treatment being identified in the ANOVA. These results demonstrate for the first time that the performance enhancing effects of caffeine may not be repeatable between days, where our data uniquely indicates that this is in part attributable to between sessions variation in caffeine’s ergogenic potential.

Functional Connectome-Based Predictive Modeling in Autism

Autism is a heterogeneous neurodevelopmental condition, and functional magnetic resonance imaging-based studies have helped advance our understanding of its effects on brain network activity. We review how predictive modeling, using measures of functional connectivity and symptoms, has helped reveal key insights into this condition. We discuss how different prediction frameworks can further our understanding of the brain-based features that underlie complex autism symptomatology and consider how predictive models may be used in clinical settings. Throughout, we highlight aspects of study interpretation, such as data decay and sampling biases, that require consideration within the context of this condition. We close by suggesting exciting future directions for predictive modeling in autism.

Spatial variation of changes in test-retest reliability of functional connectivity after global signal regression: The effect of considering hemodynamic delay

Global signal regression (GSR) is a controversial analysis method, since its removal of signal has been observed to reduce the reliability of functional connectivity estimates. Here, we used test-retest reliability to characterize potential differences in spatial patterns between conventional, static GSR (sGSR) and a novel dynamic form of GSR (dGSR). In contrast with sGSR, dGSR models the global signal at a time delay to correct for blood arrival time. Thus, dGSR accounts for greater variation in global signal, removes blood-flow-related nuisance signal, and leaves higher quality neuronal signal remaining. We used intraclass correlation coefficients (ICCs) to estimate the reliability of functional connectivity in 462 healthy controls from the Human Connectome Project. We tested across two factors: denoising method used (control, sGSR, and dGSR), and interacquisition interval (between days, or within session while varying phase encoding direction). Reliability was estimated regionally to identify topographic patterns for each condition. sGSR and dGSR provided global reductions in reliability compared with the non-GSR control. Test-retest reliability was highest in the frontoparietal and default mode regions, and lowest in sensorimotor cortex for all conditions. dGSR provides more effective denoising in regions where both strategies greatly reduce reliability. Both GSR methods substantially reduced test-retest reliability, which was most evident in brain regions that had low reliability prior to denoising. These findings suggest that reliability of interregional correlation is likely inflated by the global signal, which is thought to primarily reflect dynamic blood flow.

A comparison of multiband and multiband multiecho gradient-echo EPI for task fMRI at 3 T

A multiband (MB) echo-planar imaging (EPI) sequence is compared to a multiband multiecho (MBME) EPI protocol to investigate differences in sensitivity for task functional magnetic resonance imaging (fMRI) at 3 T. Multiecho sampling improves sensitivity in areas where single-echo-EPI suffers from dropouts. However, It requires in-plane acceleration to reduce the echo train length, limiting the slice acceleration factor and the temporal and spatial resolution Data were acquired for both protocols in two sessions 24 h apart using an adapted color-word interference Stroop task. Besides protocol comparison statistically, we performed test-retest reliability across sessions for different protocols and denoising methods. We evaluated the sensitivity of two different echo-combination strategies for MBME-EPI. We examined the performance of three different data denoising approaches: "Standard," "AROMA," and "FIX" for MB and MBME, and assessed whether a specific method is preferable. We consider using an appropriate autoregressive model order within the general linear model framework to correct TR differences between the protocols. The comparison between protocols and denoising methods showed at group level significantly higher mean z-scores and the number of active voxels for MBME in the motor, subcortical and medial frontal cortices. When comparing different echo combinations, our results suggest that a contrast-to-noise ratio weighted echo combination improves sensitivity in MBME compared to simple echo-summation. This study indicates that MBME can be a preferred protocol in task fMRI at spatial resolution (≥2 mm), primarily in medial prefrontal and subcortical areas.

Benchmarking functional connectivity by the structure and geometry of the human brain

The brain's structural connectivity supports the propagation of electrical impulses, manifesting as patterns of coactivation, termed functional connectivity. Functional connectivity emerges from the underlying sparse structural connections, particularly through polysynaptic communication. As a result, functional connections between brain regions without direct structural links are numerous, but their organization is not completely understood. Here we investigate the organization of functional connections without direct structural links. We develop a simple, data-driven method to benchmark functional connections with respect to their underlying structural and geometric embedding. We then use this method to reweigh and reexpress functional connectivity. We find evidence of unexpectedly strong functional connectivity among distal brain regions and within the default mode network. We also find unexpectedly strong functional connectivity at the apex of the unimodal-transmodal hierarchy. Our results suggest that both phenomena-functional modules and functional hierarchies-emerge from functional interactions that transcend the underlying structure and geometry. These findings also potentially explain recent reports that structural and functional connectivity gradually diverge in transmodal cortex. Collectively, we show how structural connectivity and geometry can be used as a natural frame of reference with which to study functional connectivity patterns in the brain.

The Location Reliability of the Resting-State fMRI FC of Emotional Regions Towards rTMS Therapy

Resting-state magnetic resonance imaging (RS-fMRI) studies indicated that the repetitive transcranial magnetic stimulation (rTMS) exerts antidepression effect through the functional connectivity (FC) of the DLPFC with the subgenual anterior cingulate cortex (sgACC), pregneual ACC (pgACC), or nucleus accumbens (NAc). It is proposed that the FC-guided individualized precise stimulation on the DLPFC would be more effective. The current study systematically investigated the reliability of the RS-fMRI FC location as well as the FC strength with multiple potential factors. We aimed to provide a stable stimulation target for future FC-guided TMS therapy for affective related disorders. Twenty-one subjects under RS-fMRI conditions with the first two times (V1, V2) scanned on a GE 3 T scanner and the third visit (V3) on a Siemens 3 T scanner. Then the FC strength and location reliability were assessed by using intra-class correlation (ICC) and intra-individual distance, respectively. The factors included deep seed ROIs (midline (mid-) sgACC, left pgACC, mid-pgACC, and left NAc), eyes closed (EC) vs eyes open (EO), frequency bands, FC algorithm (Pearson vs Spearman), scanning length (half a session vs whole session), and location method (FC peak vs center of gravity (COG)). The reliability of the voxel-wise FC strength was low to moderate. The intra-individual distances of the COG were 3.8-7.3 mm across all factors, much smaller than that of FC peak (approximately 30 mm). The COG of seed-based FC might be a potential rTMS stimulation target. Anyway, all potential stimulation targets should be tested in future rTMS treatment studies.

Reliability and Sensitivity to Alterered Hemodynamics Measured with Resting-state fMRI Metrics: Comparison with 123I-IMP SPECT

Blood oxygenation level-dependent (BOLD) contrast is sensitive to local hemodynamic changes and thus is applicable to imaging perfusion or vascular reactivity. However, knowledge about its measurement characteristics compared to reference standard perfusion imaging is limited. This study longitudinally evaluated perfusion in patients with steno-occlusive disease using resting-state functional MRI (rsfMRI) acquired before and within nine days of anterior circulation revascularization in patients with large cerebral artery steno-occlusive diseases. The reliability and sensitivity to longitudinal changes of rsfMRI temporal correlation (Rc) and time delay (TDc) relative to the cerebellar signal were examined voxel-wise in comparison with single-photon emission CT (SPECT) cerebral blood flow (CBF) using the within-subject standard deviation (Sw) and intraclass correlation coefficients (ICCs). For statistical comparisons, the standard deviation (SD) of longitudinal changes within the cerebellum, the number of voxels with significant changes in the left middle cerebral artery territory ipsilateral to surgery, and their average changes relative to the cerebellar SD were evaluated. The test-retest reliability of the fMRI metrics was also similarly evaluated using the human connectome project (HCP) healthy young adult dataset. The test-retest time interval was 31 ± 18 days. Test-retest reliability was significantly higher for SPECT (cerebellar SD: -2.59 ± 0.20) than for fMRI metrics (cerebellar SD: Rc, -2.34 ± 0.24, p = 0.04; TDc, -2.19 ± 0.21, p = 0.003). Sensitivity to postoperative changes, which was evaluated as the number of voxels, was significantly higher for fMRI TDc (8.78 ± 0.72) than for Rc (7.42 ± 1.48, p = 0.03) or SPECT CBF (6.88 ± 0.67, p < 0.001). The ratio between the average Rc, TDc, and SPECT CBF changes within the left MCA target region and cerebellar SD was also significantly higher for fMRI TDc (1.21 ± 0.79) than Rc (0.48 ± 0.94, p = 0.006) or SPECT CBF (0.23 ± 0.57, p = 0.001). The measurement variability of time delay was also larger than that of temporal correlation in HCP data within the cerebellum (t = -8.7, p < 0.001) or in the whole-brain (t = -27.4, p < 0.001) gray matter. These data suggest that fMRI time delay is more sensitive to the hemodynamic changes than SPECT CBF, although the reliability is lower. The implication for fMRI connectivity studies is that temporal correlation can be significantly decreased due to altered hemodynamics, even in cases with normal CBF.

Validity of age estimation methods and reproducibility of bone/dental maturity indices for chronological age estimation: a systematic review and meta-analysis of validation studies

Several approaches have been developed to estimate age, an important aspect of forensics and orthodontics, using different measures and radiological examinations. Here, through meta-analysis, we determined the validity of age estimation methods and reproducibility of bone/dental maturity indices used for age estimation. The PubMed and Google Scholar databases were searched to December 31, 2021 for human cross-sectional studies meeting pre-defined PICOS criteria that simultaneously assessed the reproducibility and validity. Meta-estimates of validity (mean error: estimated age-chronological age) and intra- and inter-observer reproducibility (Cohen's kappa, intraclass correlation coefficient) and their predictive intervals (PI) were calculated using mixed-effect models when heterogeneity was high (I2 > 50%). The literature search identified 433 studies, and 23 met the inclusion criteria. The mean error meta-estimate (mixed effects model) was 0.08 years (95% CI - 0.12; 0.29) in males and 0.09 (95% CI - 0.12; 0.30) in females. The PI of each method spanned zero; of nine reported estimation methods, Cameriere's had the smallest (- 0.82; 0.47) and Haavikko's the largest (- 7.24; 4.57) PI. The reproducibility meta-estimate (fixed effects model) was 0.98 (95% CI 0.97; 1.00) for intra- and 0.99 (95% CI 0.98; 1.00) for inter-observer agreement. All methods were valid but with different levels of precision. The intra- and inter-observer reproducibility was high and homogeneous across studies.

Functional network alterations in young brain tumor patients with radiotherapy-induced memory impairments and vascular injury

Background

Cognitive impairment and cerebral microbleeds (CMBs) are long-term side-effects of cranial radiation therapy (RT). Previously we showed that memory function is disrupted in young patients and that the rate of cognitive decline correlates with CMB development. However, vascular injury alone cannot explain RT-induced cognitive decline. Here we use resting-state functional MRI (rsfMRI) to further investigate the complex mechanisms underlying memory impairment after RT.

Methods

Nineteen young patients previously treated with or without focal or whole-brain RT for a brain tumor underwent cognitive testing followed by 7T rsfMRI and susceptibility-weighted imaging for CMB detection. Global brain modularity and efficiency, and rsfMRI signal variability within the dorsal attention, salience, and frontoparietal networks were computed. We evaluated whether MR metrics could distinguish age- and sex-matched controls (N = 19) from patients and differentiate patients based on RT exposure and aggressiveness. We also related MR metrics with memory performance, CMB burden, and risk factors for cognitive decline after RT.

Results

Compared to controls, patients exhibited widespread hyperconnectivity, similar modularity, and significantly increased efficiency (p < 0.001) and network variability (p < 0.001). The most abnormal values were detected in patients treated with high dose whole-brain RT, having supratentorial tumors, and who did not undergo RT but had hydrocephalus. MR metrics and memory performance were correlated (R = 0.34–0.53), though MR metrics were more strongly related to risk factors for cognitive worsening and CMB burden with evidence of functional recovery.

Conclusions

MR metrics describing brain connectivity and variability represent promising candidate imaging biomarkers for monitoring of long-term cognitive side-effects after RT.

Choice of Voxel-based Morphometry processing pipeline drives variability in the location of neuroanatomical brain markers

Fundamental and clinical neuroscience has benefited tremendously from the development of automated computational analyses. In excess of 600 human neuroimaging papers using Voxel-based Morphometry (VBM) are now published every year and a number of different automated processing pipelines are used, although it remains to be systematically assessed whether they come up with the same answers. Here we examined variability between four commonly used VBM pipelines in two large brain structural datasets. Spatial similarity and between-pipeline reproducibility of the processed gray matter brain maps were generally low between pipelines. Examination of sex-differences and age-related changes revealed considerable differences between the pipelines in terms of the specific regions identified. Machine learning-based multivariate analyses allowed accurate predictions of sex and age, however accuracy differed between pipelines. Our findings suggest that the choice of pipeline alone leads to considerable variability in brain structural markers which poses a serious challenge for reproducibility and interpretation.

Four common processing pipelines tested on two Voxel-based Morphometry (VBM) datasets yield considerable variations in results, raising issues on the interpretability and robustness of VBM results.

Reproducibility of the Blood and Urine Exposome: A Systematic Literature Review and Meta-Analysis

Endogenous and exogenous metabolite concentrations may be susceptible to variation over time. This variability can lead to misclassification of exposure levels and in turn to biased results. To assess the reproducibility of metabolites, the intraclass correlation coefficient (ICC) is computed. A literature search in three databases from 2000 to May 2021 was conducted to identify studies reporting ICCs for blood and urine metabolites. This review includes 192 studies, of which 31 studies are included in the meta-analyses. The ICCs of 359 single metabolites are reported, and the ICCs of 10 metabolites were meta-analyzed. The reproducibility of the single metabolites ranges from poor to excellent and is highly compound-dependent. The reproducibility of bisphenol A (BPA), mono-ethyl phthalate (MEP), mono-n-butyl phthalate (MnBP), mono-2-ethylhexyl phthalate (MEHP), mono(2-ethyl-5-hydroxyhexyl) phthalate (MEHHP), mono-benzyl phthalate (MBzP), mono-(2-ethyl-5-oxohexyl) phthalate (MEOHP), methylparaben, and propylparaben is poor to moderate (ICC median: 0.32; range: 0.15-0.49), and for 25-hydroxyvitamin D [25(OH)D], it is excellent (ICC: 0.95; 95% CI, 0.90-0.99). Pharmacokinetics, mainly the half-life of elimination and exposure patterns, can explain reproducibility. This review describes the reproducibility of the blood and urine exposome, provides a vast dataset of ICC estimates, and hence constitutes a valuable resource for future reproducibility and clinical epidemiologic studies.

Resting-state functional connectivity does not predict individual differences in the effects of emotion on memory

Emotion-laden events and objects are typically better remembered than neutral ones. This is usually explained by stronger functional coupling in the brain evoked by emotional content. However, most research on this issue has focused on functional connectivity evoked during or after learning. The effect of an individual’s functional connectivity at rest is unknown. Our pre-registered study addresses this issue by analysing a large database, the Cambridge Centre for Ageing and Neuroscience, which includes resting-state data and emotional memory scores from 303 participants aged 18–87 years. We applied regularised regression to select the relevant connections and replicated previous findings that whole-brain resting-state functional connectivity can predict age and intelligence in younger adults. However, whole-brain functional connectivity predicted neither an emotional enhancement effect (i.e., the degree to which emotionally positive or negative events are remembered better than neutral events) nor a positivity bias effect (i.e., the degree to which emotionally positive events are remembered better than negative events), failing to support our pre-registered hypotheses. These results imply a small or no association between individual differences in functional connectivity at rest and emotional memory, and support recent notions that resting-state functional connectivity is not always useful in predicting individual differences in behavioural measures.

Functional neuroimaging in psychiatry and the case for failing better

Psychiatric disorders encompass complex aberrations of cognition and affect and are among the most debilitating and poorly understood of any medical condition. Current treatments rely primarily on interventions that target brain function (drugs) or learning processes (psychotherapy). A mechanistic understanding of how these interventions mediate their therapeutic effects remains elusive. From the early 1990s, non-invasive functional neuroimaging, coupled with parallel developments in the cognitive neurosciences, seemed to signal a new era of neurobiologically grounded diagnosis and treatment in psychiatry. Yet, despite three decades of intense neuroimaging research, we still lack a neurobiological account for any psychiatric condition. Likewise, functional neuroimaging plays no role in clinical decision making. Here, we offer a critical commentary on this impasse and suggest how the field might fare better and deliver impactful neurobiological insights.

The effect of general anesthesia on the test–retest reliability of resting-state fMRI metrics and optimization of scan length

Resting-state functional magnetic resonance imaging (rs-fMRI) has been known as a powerful tool in neuroscience. However, exploring the test–retest reliability of the metrics derived from the rs-fMRI BOLD signal is essential, particularly in the studies of patients with neurological disorders. Here, two factors, namely, the effect of anesthesia and scan length, have been estimated on the reliability of rs-fMRI measurements. A total of nine patients with drug-resistant epilepsy (DRE) requiring interstitial thermal therapy (LITT) were scanned in two states. The first scan was performed in an awake state before surgery on the same patient. The second scan was performed 2 weeks later under general anesthesia necessary for LITT surgery. At each state, two rs-fMRI sessions were obtained that each one lasted 15 min, and the effect of scan length was evaluated. Voxel-wise rs-fMRI metrics, including the amplitude of low-frequency fluctuation (ALFF), the fractional amplitude of low-frequency fluctuation (fALFF), functional connectivity (FC), and regional homogeneity (ReHo), were measured. Intraclass correlation coefficient (ICC) was calculated to estimate the reliability of the measurements in two states of awake and under anesthesia. Overall, it appeared that the reliability of rs-fMRI metrics improved under anesthesia. From the 15-min data, we found mean ICC values in awake state including 0.81, 0.51, 0.65, and 0.84 for ALFF, fALFF, FC, and ReHo, respectively, as well as 0.80, 0.59, 0.83, and 0.88 for ALFF, fALFF, FC, and ReHo, respectively, under anesthesia. Additionally, our findings revealed that reliability increases as the function of scan length. We showed that the optimized scan length to achieve less variability of rs-fMRI measurements was 3.1–7.5 min shorter in an anesthetized, compared to a wakeful state.

Reliability and subject specificity of personalized whole-brain dynamical models

Dynamical whole-brain models were developed to link structural (SC) and functional connectivity (FC) together into one framework. Nowadays, they are used to investigate the dynamical regimes of the brain and how these relate to behavioral, clinical and demographic traits. However, there is no comprehensive investigation on how reliable and subject specific the modeling results are given the variability of the empirical FC. In this study, we show that the parameters of these models can be fitted with a "poor" to "good" reliability depending on the exact implementation of the modeling paradigm. We find, as a general rule of thumb, that enhanced model personalization leads to increasingly reliable model parameters. In addition, we observe no clear effect of the model complexity evaluated by separately sampling results for linear, phase oscillator and neural mass network models. In fact, the most complex neural mass model often yields modeling results with "poor" reliability comparable to the simple linear model, but demonstrates an enhanced subject specificity of the model similarity maps. Subsequently, we show that the FC simulated by these models can outperform the empirical FC in terms of both reliability and subject specificity. For the structure-function relationship, simulated FC of individual subjects may be identified from the correlations with the empirical SC with an accuracy up to 70%, but not vice versa for non-linear models. We sample all our findings for 8 distinct brain parcellations and 6 modeling conditions and show that the parcellation-induced effect is much more pronounced for the modeling results than for the empirical data. In sum, this study provides an exploratory account on the reliability and subject specificity of dynamical whole-brain models and may be relevant for their further development and application. In particular, our findings suggest that the application of the dynamical whole-brain modeling should be tightly connected with an estimate of the reliability of the results.

Temporal stability of the hemodynamic response function across the majority of human cerebral cortex

The hemodynamic response function (HRF) measured with functional magnetic resonance imaging is generated by vascular and metabolic responses evoked by brief (<4 s) stimuli. It is known that the human HRF varies across cortex, between subjects, with stimulus paradigms, and even between different measurements in the same cortical location. However, our results demonstrate that strong HRFs are remarkably repeatable across sessions separated by time intervals up to 3 months. In this study, a multisensory stimulus was used to activate and measure the HRF across the majority of cortex (>70%, with lesser reliability observed in some areas of prefrontal cortex). HRFs were measured with high spatial resolution (2‐mm voxels) in central gray matter to minimize variations caused by partial‐volume effects. HRF amplitudes and temporal dynamics were highly repeatable across four sessions in 20 subjects. Positive and negative HRFs were consistently observed across sessions and subjects. Negative HRFs were generally weaker and, thus, more variable than positive HRFs. Statistical measurements showed that across‐session variability is highly correlated to the variability across events within a session; these measurements also indicated a normal distribution of variability across cortex. The overall repeatability of the HRFs over long time scales generally supports the long‐term use of event‐related functional magnetic resonance imaging protocols.

Intrinsic brain abnormalities in female major depressive disorder patients with childhood trauma: A resting-state functional magnetic resonance imaging study

Objective

Childhood trauma is a strong predictor of major depressive disorder (MDD). Women are more likely to develop MDD than men. However, the neural basis of female MDD patients with childhood trauma remains unclear. We aimed to identify the specific brain regions that are associated with female MDD patients with childhood trauma.

Methods

We recruited 16 female MDD patients with childhood trauma, 16 female MDD patients without childhood trauma, and 20 age- and education level-matched healthy controls. All participants underwent resting-state functional magnetic resonance imaging (MRI). Regional brain activity was evaluated as the amplitude of low-frequency fluctuation (ALFF). Furthermore, functional connectivity (FC) analyses were performed on areas with altered ALFF to explore alterations in FC patterns.

Results

There was increased ALFF in the left middle frontal gyrus (MFG) and the right postcentral gyrus (PoCG) in MDD with childhood trauma compared with MDD without childhood trauma. The areas with significant ALFF discrepancies were selected as seeds for the FC analyses. There was increased FC between the left MFG and the bilateral putamen gyrus. Moreover, ALFF values were correlated with childhood trauma severity.

Conclusion

Our findings revealed abnormal intrinsic brain activity and FC patterns in female MDD patients with childhood trauma, which provides new possibilities for exploring the pathophysiology of this disorder in women.

A cognitive neurogenetic approach to uncovering the structure of executive functions

One central mission of cognitive neuroscience is to understand the ontology of complex cognitive functions. We addressed this question with a cognitive neurogenetic approach using a large-scale dataset of executive functions (EFs), whole-brain resting-state functional connectivity, and genetic polymorphisms. We found that the bifactor model with common and shifting-specific components not only was parsimonious but also showed maximal dissociations among the EF components at behavioral, neural, and genetic levels. In particular, the genes with enhanced expression in the middle frontal gyrus (MFG) and the subcallosal cingulate gyrus (SCG) showed enrichment for the common and shifting-specific component, respectively. Finally, High-dimensional mediation models further revealed that the functional connectivity patterns significantly mediated the genetic effect on the common EF component. Our study not only reveals insights into the ontology of EFs and their neurogenetic basis, but also provides useful tools to uncover the structure of complex constructs of human cognition.

Covariance shrinkage can assess and improve functional connectomes

Connectomes derived from resting-state functional MRI scans have significantly benefited from the development of dedicated fMRI motion correction and denoising algorithms. But they are based on empirical correlations that can produce unreliable results in high dimension low sample size settings. A family of statistical estimators, the covariance shrinkage methods, could mitigate this issue. Unfortunately, these methods have rarely been used to correct functional connectomes and no extensive experiment has been conducted so far to compare the shrinkage methods available for this task. In this work, we propose to fix this issue by processing a benchmark dataset made of a thousand high-resolution resting-state fMRI scans provided by the Human Connectome Project to compare the ability of five prominent covariance shrinkage methods to produce reliable functional connectomes at different spatial resolutions and scans duration: the pioneer linear covariance shrinkage method introduced by Ledoit and Wolf, the Oracle Approximating Shrinkage, the QuEST method, the NERCOME method, and a recent analytical approximation of the QuEST approach. Our experiments establish that all covariance shrinkage methods significantly improve functional connectomes derived from short fMRI scans. The Oracle Approximating Shrinkage and the QuEST method produced the best results. Lastly, we present shrinkage intensity charts that can be used for designing and analyzing fMRI studies. These charts indicate that sparse connectomes are difficult to estimate from short fMRI scans, and they describe a range of settings where dynamic functional connectivity should not be computed.

Human electromagnetic and haemodynamic networks systematically converge in unimodal cortex and diverge in transmodal cortex

Whole-brain neural communication is typically estimated from statistical associations among electromagnetic or haemodynamic time-series. The relationship between functional network architectures recovered from these 2 types of neural activity remains unknown. Here, we map electromagnetic networks (measured using magnetoencephalography (MEG)) to haemodynamic networks (measured using functional magnetic resonance imaging (fMRI)). We find that the relationship between the 2 modalities is regionally heterogeneous and systematically follows the cortical hierarchy, with close correspondence in unimodal cortex and poor correspondence in transmodal cortex. Comparison with the BigBrain histological atlas reveals that electromagnetic-haemodynamic coupling is driven by laminar differentiation and neuron density, suggesting that the mapping between the 2 modalities can be explained by cytoarchitectural variation. Importantly, haemodynamic connectivity cannot be explained by electromagnetic activity in a single frequency band, but rather arises from the mixing of multiple neurophysiological rhythms. Correspondence between the two is largely driven by MEG functional connectivity at the beta (15 to 29 Hz) frequency band. Collectively, these findings demonstrate highly organized but only partly overlapping patterns of connectivity in MEG and fMRI functional networks, opening fundamentally new avenues for studying the relationship between cortical microarchitecture and multimodal connectivity patterns.

Computational approaches to predicting treatment response to obesity using neuroimaging

Obesity is a worldwide disease associated with multiple severe adverse consequences and comorbid conditions. While an increased body weight is the defining feature in obesity, etiologies, clinical phenotypes and treatment responses vary between patients. These variations can be observed within individual treatment options which comprise lifestyle interventions, pharmacological treatment, and bariatric surgery. Bariatric surgery can be regarded as the most effective treatment method. However, long-term weight regain is comparably frequent even for this treatment and its application is not without risk. A prognostic tool that would help predict the effectivity of the individual treatment methods in the long term would be essential in a personalized medicine approach. In line with this objective, an increasing number of studies have combined neuroimaging and computational modeling to predict treatment outcome in obesity. In our review, we begin by outlining the central nervous mechanisms measured with neuroimaging in these studies. The mechanisms are primarily related to reward-processing and include "incentive salience" and psychobehavioral control. We then present the diverse neuroimaging methods and computational prediction techniques applied. The studies included in this review provide consistent support for the importance of incentive salience and psychobehavioral control for treatment outcome in obesity. Nevertheless, further studies comprising larger sample sizes and rigorous validation processes are necessary to answer the question of whether or not the approach is sufficiently accurate for clinical real-world application.

Resting-state functional connectivity identifies individuals and predicts age in 8-to-26-month-olds

Resting-state functional connectivity (rsFC) measured with fMRI has been used to characterize functional brain maturation in typically and atypically developing children and adults. However, its reliability and utility for predicting development in infants and toddlers is less well understood. Here, we use fMRI data from the Baby Connectome Project study to measure the reliability and uniqueness of rsFC in infants and toddlers and predict age in this sample (8-to-26 months old; n = 170). We observed medium reliability for within-session infant rsFC in our sample, and found that individual infant and toddler's connectomes were sufficiently distinct for successful functional connectome fingerprinting. Next, we trained and tested support vector regression models to predict age-at-scan with rsFC. Models successfully predicted novel infants' age within ± 3.6 months error and a prediction R2 = .51. To characterize the anatomy of predictive networks, we grouped connections into 11 infant-specific resting-state functional networks defined in a data-driven manner. We found that connections between regions of the same network-i.e. within-network connections-predicted age significantly better than between-network connections. Looking ahead, these findings can help characterize changes in functional brain organization in infancy and toddlerhood and inform work predicting developmental outcome measures in this age range.

Adolescent depression and resting-state fMRI brain networks: a scoping review of longitudinal studies

The neurobiological factors associated with the emergence of major depressive disorder (MDD) in adolescence are still unclear. Previous cross-sectional studies have documented aberrant connectivity in resting-state functional magnetic resonance imaging (rs-fMRI) networks. However, whether these findings precede MDD onset has not been established. This scoping review mapped key methodological aspects and main findings of longitudinal rs-fMRI studies of MDD in adolescence. Three sets of neuroimaging methods to analyze rs-fMRI data were identified: seed-based analysis, independent component analysis, and network-based approaches. Main findings involved aberrant connectivity within and between the default mode network (DMN), the cognitive control network (CCN), and the salience network (SN). Accordingly, we utilized Menon's (2011) triple-network model for neuropsychiatric disorders to summarize key results. Adolescent MDD was associated with hyperconnectivity within the SN and between DMN and SN, as well as hypoconectivity within the CCN. These findings suggested that dysfunctional connectivity among the three main large-scale brain networks preceded MDD onset. However, there was high heterogeneity in neuroimaging methods and sampling procedures, which may limit comparisons between studies. Future studies should consider some level of harmonization for clinical instruments and neuroimaging methods.

Reproducibility and repeatability of magnetic resonance imaging in dementia

PURPOSE

Individualised predictive models of cognitive decline require disease-monitoring markers that are repeatable. For wide-spread adoption, such markers also need to be reproducible at different locations. This study assessed the repeatability and reproducibility of MRI markers derived from a dementia protocol.

METHODS

Six participants were scanned at three different sites with a 3T MRI scanner. The protocol employed: T1-weighted (T1w) imaging, resting state functional MRI (rsfMRI), arterial spin labelling (ASL), diffusion-weighted imaging (DWI), T2-weighted fluid attenuation inversion recovery (FLAIR), T2-weighted (T2w) imaging, and susceptibility weighted imaging (SWI). Participants were scanned repeatedly, up to six times over a maximum period of five years. One participant was also scanned a further three times on sequential days on one scanner. Fifteen derived metrics were computed from the seven different modalities.

RESULTS

Reproducibility (coefficient of variation; CoV, across sites) was best for T1w derived grey matter, white matter and hippocampal volume (CoV < 1.5%), compared to rsfMRI and SWI derived metrics (CoV, 19% and 21%). For a given metric, long-term repeatability (CoV across time) was comparable to reproducibility, with short-term repeatability considerably better.

CONCLUSIONS

Reproducibility and repeatability were assessed for a suite of markers calculated from a dementia MRI protocol. In general, structural markers were less variable than functional MRI markers. Variability over time on the same scanner was comparable to variability measured across different scanners. Overall, the results support the viability of multi-site longitudinal studies for monitoring cognitive decline.

NeoRS: A Neonatal Resting State fMRI Data Preprocessing Pipeline

Resting state functional MRI (rsfMRI) has been shown to be a promising tool to study intrinsic brain functional connectivity and assess its integrity in cerebral development. In neonates, where functional MRI is limited to very few paradigms, rsfMRI was shown to be a relevant tool to explore regional interactions of brain networks. However, to identify the resting state networks, data needs to be carefully processed to reduce artifacts compromising the interpretation of results. Because of the non-collaborative nature of the neonates, the differences in brain size and the reversed contrast compared to adults due to myelination, neonates can’t be processed with the existing adult pipelines, as they are not adapted. Therefore, we developed NeoRS, a rsfMRI pipeline for neonates. The pipeline relies on popular neuroimaging tools (FSL, AFNI, and SPM) and is optimized for the neonatal brain. The main processing steps include image registration to an atlas, skull stripping, tissue segmentation, slice timing and head motion correction and regression of confounds which compromise functional data interpretation. To address the specificity of neonatal brain imaging, particular attention was given to registration including neonatal atlas type and parameters, such as brain size variations, and contrast differences compared to adults. Furthermore, head motion was scrutinized, and motion management optimized, as it is a major issue when processing neonatal rsfMRI data. The pipeline includes quality control using visual assessment checkpoints. To assess the effectiveness of NeoRS processing steps we used the neonatal data from the Baby Connectome Project dataset including a total of 10 neonates. NeoRS was designed to work on both multi-band and single-band acquisitions and is applicable on smaller datasets. NeoRS also includes popular functional connectivity analysis features such as seed-to-seed or seed-to-voxel correlations. Language, default mode, dorsal attention, visual, ventral attention, motor and fronto-parietal networks were evaluated. Topology found the different analyzed networks were in agreement with previously published studies in the neonate. NeoRS is coded in Matlab and allows parallel computing to reduce computational times; it is open-source and available on GitHub (https://github.com/venguix/NeoRS). NeoRS allows robust image processing of the neonatal rsfMRI data that can be readily customized to different datasets.

Functional connectomes become more longitudinally self-stable, but not more distinct from others, across early childhood

Functional connectomes, as measured with functional magnetic resonance imaging (fMRI), are highly individualized, and evidence suggests this individualization may increase across childhood. A connectome can become more individualized either by increasing self-stability or decreasing between-subject-similarity. Here we used a longitudinal early childhood dataset to investigate age associations with connectome self-stability, between-subject-similarity, and developmental individualization, defined as an individual's self-stability across a 12-month interval relative to their between-subject-similarity. fMRI data were collected during an 18-minute passive viewing scan from 73 typically developing children aged 4-7 years, at baseline and 12-month follow-up. We found that young children had highly individualized connectomes, with sufficient self-stability across 12-months for 98% identification accuracy. Linear models showed a significant relationship between age and developmental individualization across the whole brain and in most networks. This association appeared to be largely driven by an increase in self-stability with age, with only weak evidence for relationships between age and similarity across participants. Together our findings suggest that children's connectomes become more individualized across early childhood, and that this effect is driven by increasing self-stability rather than decreasing between-subject-similarity.

A graph theory-based whole brain approach to assess mood state differences in adolescents and young adults with bipolar disorder

OBJECTIVES

Identifying hubs of brain dysfunction in adolescents and young adults with Bipolar I Disorder (BDAYA ) could provide targets for early detection, prevention, and treatment. Previous neuroimaging studies across mood states of BDAYA are scarce and often examined limited brain regions potentially prohibiting detection of other important regions. We used a data-driven whole-brain Intrinsic Connectivity Distribution (ICD) approach to investigate dysconnectivity hubs across mood states in BDAYA .

METHODS

Functional magnetic resonance imaging whole-brain ICD data were investigated for differences across four groups: BDAYA -depressed (n = 22), BDAYA -euthymic (n = 45), BDAYA -elevated (n = 24), and healthy controls (HC, n = 111). Clusters of ICD differences were assessed for regional dysconnectivity and mood symptom relationships. Analyses were also performed for BDAYA overall (vs. HC) ICD differences persisting across mood states.

RESULTS

ICD was higher in the BDAYA- depressed group than other groups in bilateral ventral/rostral/dorsal prefrontal cortex (PFC) and right lenticular nucleus (LN) (pcorrected  <0.05). In BDAYA -depressed, functional connectivity (FC) was increased between these regions with their contralateral homologues and PFC-medial temporal FC was more negative (p < 0.005). PFC-related findings correlated with depression scores (p < 0.05). The overall BDAYA group showed ICD increases in more ventral left PFC and right cerebellum, present across euthymia and acute mood states.

CONCLUSIONS

This ICD approach supports a PFC hub of inter- and intra-hemispheric frontotemporal dysconnectivity in BDAYA with potential trait features and disturbances of higher magnitude during depression. Hubs were also revealed in LN and cerebellum, less common foci of BD research. The hubs are potential targets for early interventions to detect, prevent, and treat BD.

Future Directions for Chemosensory Connectomes: Best Practices and Specific Challenges

Ecological chemosensory stimuli almost always evoke responses in more than one sensory system. Moreover, any sensory processing takes place along a hierarchy of brain regions. So far, the field of chemosensory neuroimaging is dominated by studies that examine the role of brain regions in isolation. However, to completely understand neural processing of chemosensation, we must also examine interactions between regions. In general, the use of connectivity methods has increased in the neuroimaging field, providing important insights to physical sensory processing, such as vision, audition, and touch. A similar trend has been observed in chemosensory neuroimaging, however, these established techniques have largely not been rigorously applied to imaging studies on the chemical senses, leaving network insights overlooked. In this article, we first highlight some recent work in chemosensory connectomics and we summarize different connectomics techniques. Then, we outline specific challenges for chemosensory connectome neuroimaging studies. Finally, we review best practices from the general connectomics and neuroimaging fields. We recommend future studies to develop or use the following methods we perceive as key to improve chemosensory connectomics: (1) optimized study designs, (2) reporting guidelines, (3) consensus on brain parcellations, (4) consortium research, and (5) data sharing.

Identification of Pathogenetic Brain Regions via Neuroimaging Data for Diagnosis of Autism Spectrum Disorders

Autism spectrum disorder (ASD) is a kind of neurodevelopmental disorder that often occurs in children and has a hidden onset. Patients usually have lagged development of communication ability and social behavior and thus suffer an unhealthy physical and mental state. Evidence has indicated that diseases related to ASD have commonalities in brain imaging characteristics. This study aims to study the pathogenesis of ASD based on brain imaging data to locate the ASD-related brain regions. Specifically, we collected the functional magnetic resonance image data of 479 patients with ASD and 478 normal subjects matched in age and gender and used a machine-learning framework named random support vector machine cluster to extract distinctive brain regions from the preprocessed data. According to the experimental results, compared with other existing approaches, the method used in this study can more accurately distinguish patients from normal individuals based on brain imaging data. At the same time, this study found that the development of ASD was highly correlated with certain brain regions, e.g., lingual gyrus, superior frontal gyrus, medial gyrus, insular lobe, and olfactory cortex. This study explores the effectiveness of a novel machine-learning approach in the study of ASD brain imaging and provides a reference brain area for the medical research and clinical treatment of ASD.