Comparison of heritability estimates on resting state fMRI connectivity phenotypes using the ENIGMA analysis pipeline

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

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

Shared genetics linking sociability with the brain's default mode network

The brain's default mode network (DMN) plays a role in social cognition, with altered DMN function being associated with social impairments across various neuropsychiatric disorders. In the present study, we examined the genetic relationship between sociability and DMN-related resting-state functional magnetic resonance imaging (rs-fMRI) traits. To this end, we used genome-wide association summary statistics for sociability and 31 activity and 64 connectivity DMN-related rs-fMRI traits (N=34,691-342,461). First, we examined global and local genetic correlations between sociability and the rs-fMRI traits. Second, to assess putatively causal relationships between the traits, we conducted bi-directional Mendelian randomisation (MR) analyses. Finally, we prioritised genes influencing both sociability and rs-fMRI traits by combining three methods: gene-expression eQTL MR analyses, the CELLECT framework using single-nucleus RNA-seq data, and network propagation in the context of a protein-protein interaction network. Significant local genetic correlations were found between sociability and two rs-fMRI traits, one representing spontaneous activity within the temporal cortex, the other representing connectivity between the frontal/cingulate and angular/temporal cortices. Sociability affected 12 rs-fMRI traits when allowing for weakly correlated genetic instruments. Combing all three methods for gene prioritisation, we defined 17 highly prioritised genes, with DRD2 and LINGO1 showing the most robust evidence across all analyses. By integrating genetic and transcriptomics data, our gene prioritisation strategy may serve as a blueprint for future studies. The prioritised genes could be explored as potential biomarkers for social dysfunction in the context of neuropsychiatric disorders and as drug target genes.

Modeling brain sex in the limbic system as phenotype for female-prevalent mental disorders

BACKGROUND

Sex differences exist in the prevalence and clinical manifestation of several mental disorders, suggesting that sex-specific brain phenotypes may play key roles. Previous research used machine learning models to classify sex from imaging data of the whole brain and studied the association of class probabilities with mental health, potentially overlooking regional specific characteristics.

METHODS

We here investigated if a regionally constrained model of brain volumetric imaging data may provide estimates that are more sensitive to mental health than whole brain-based estimates. Given its known role in emotional processing and mood disorders, we focused on the limbic system. Using two different cohorts of healthy subjects, the Human Connectome Project and the Queensland Twin IMaging, we investigated sex differences and heritability of brain volumes of limbic structures compared to non-limbic structures, and subsequently applied regionally constrained machine learning models trained solely on limbic or non-limbic features. To investigate the biological underpinnings of such models, we assessed the heritability of the obtained sex class probability estimates, and we investigated the association with major depression diagnosis in an independent clinical sample. All analyses were performed both with and without controlling for estimated total intracranial volume (eTIV).

RESULTS

Limbic structures show greater sex differences and are more heritable compared to non-limbic structures in both analyses, with and without eTIV control. Consequently, machine learning models performed well at classifying sex based solely on limbic structures and achieved performance as high as those on non-limbic or whole brain data, despite the much smaller number of features in the limbic system. The resulting class probabilities were heritable, suggesting potentially meaningful underlying biological information. Applied to an independent population with major depressive disorder, we found that depression is associated with male-female class probabilities, with largest effects obtained using the limbic model. This association was significant for models not controlling for eTIV whereas in those controlling for eTIV the associations did not pass significance correction.

CONCLUSIONS

Overall, our results highlight the potential utility of regionally constrained models of brain sex to better understand the link between sex differences in the brain and mental disorders.

Heritability of cognitive and emotion processing during functional MRI in a twin sample

Despite compelling evidence that brain structure is heritable, the evidence for the heritability of task-evoked brain function is less robust. Findings from previous studies are inconsistent possibly reflecting small samples and methodological variations. In a large national twin sample, we systematically evaluated heritability of task-evoked brain activity derived from functional magnetic resonance imaging. We used established standardised tasks to engage brain regions involved in cognitive and emotional functions. Heritability was evaluated across a conscious and nonconscious Facial Expressions of Emotion Task (FEET), selective attention Oddball Task, N-back task of working memory maintenance, and a Go-NoGo cognitive control task in a sample of Australian adult twins (N ranged from 136 to 226 participants depending on the task and pairs). Two methods for quantifying associations of heritability and brain activity were utilised; a multivariate independent component analysis (ICA) approach and a univariate brain region-of-interest (ROI) approach. Using ICA, we observed that a significant proportion of task-evoked brain activity was heritable, with estimates ranging from 23% to 26% for activity elicited by nonconscious facial emotion stimuli, 27% to 34% for N-back working memory maintenance and sustained attention, and 32% to 33% for selective attention in the Oddball task. Using the ROI approach, we found that activity of regions specifically implicated in emotion processing and selective attention showed significant heritability for three ROIs, including estimates of 33%-34% for the left and right amygdala in the nonconscious processing of sad faces and 29% in the medial superior prefrontal cortex for the Oddball task. Although both approaches show similar levels of heritability for the Nonconscious Faces and Oddball tasks, ICA results displayed a more extensive network of heritable brain function, including additional regions beyond the ROI analysis. Furthermore, multivariate twin modelling of both ICA networks and ROI activation suggested a mix of common genetic and unique environmental factors that contribute to the associations between networks/regions. Together, the results indicate a complex relationship between genetic factors and environmental interactions that ultimately give rise to neural activation underlying cognition and emotion.

A systematic review exploring the association between the human gut microbiota and brain connectivity in health and disease

A body of pre-clinical evidence shows how the gut microbiota influence brain functioning, including brain connectivity. Linking measures of brain connectivity to the gut microbiota can provide important mechanistic insights into the bi-directional gut-brain communication. In this systematic review, we therefore synthesized the available literature assessing this association, evaluating the degree of consistency in microbiota-connectivity associations. Following the PRISMA guidelines, a PubMed search was conducted, including studies published up to September 1, 2022. We identified 16 studies that met the inclusion criteria. Several bacterial genera, including Prevotella, Bacteroides, Ruminococcus, Blautia, and Collinsella were most frequently reported in association with brain connectivity. Additionally, connectivity of the salience (specifically the insula and anterior cingulate cortex), default mode, and frontoparietal networks were most frequently associated with the gut microbiota, both in terms of microbial diversity and composition. There was no discernible pattern in the association between microbiota and brain connectivity. Altogether, based on our synthesis, there is evidence for an association between the gut microbiota and brain connectivity. However, many findings were poorly replicated across studies, and the specificity of the association is yet unclear. The current studies show substantial inter-study heterogeneity in methodology and reporting, limiting the robustness and reproducibility of the findings and emphasizing the need to harmonize methodological approaches. To enhance comparability and replicability, future research should focus on further standardizing processing pipelines and employing data-driven multivariate analysis strategies.

Genetic and environmental influence on resting state networks in young male and female adults: a cartographer mapping study

We propose a unique, minimal assumption, approach based on variance analyses (compared with standard approaches) to investigate genetic influence on individual differences on the functional connectivity of the brain using 65 monozygotic and 65 dizygotic healthy young adult twin pairs' low-frequency oscillation resting state functional Magnetic Resonance Imaging (fMRI) data from the Human Connectome Project. Overall, we found high number of genetically-influenced functional (GIF) connections involving posterior to posterior brain regions (occipital/temporal/parietal) implicated in low-level processes such as vision, perception, motion, categorization, dorsal/ventral stream visuospatial, and long-term memory processes, as well as high number across midline brain regions (cingulate) implicated in attentional processes, and emotional responses to pain. We found low number of GIF connections involving anterior to anterior/posterior brain regions (frontofrontal > frontoparietal, frontotemporal, frontooccipital) implicated in high-level processes such as working memory, reasoning, emotional judgment, language, and action planning. We found very low number of GIF connections involving subcortical/noncortical networks such as basal ganglia, thalamus, brainstem, and cerebellum. In terms of sex-specific individual differences, individual differences in males were more genetically influenced while individual differences in females were more environmentally influenced in terms of the interplay of interactions of Task positive networks (brain regions involved in various task-oriented processes and attending to and interacting with environment), extended Default Mode Network (a central brain hub for various processes such as internal monitoring, rumination, and evaluation of self and others), primary sensorimotor systems (vision, audition, somatosensory, and motor systems), and subcortical/noncortical networks. There were >8.5-19.1 times more GIF connections in males than females. These preliminary (young adult cohort-specific) findings suggest that individual differences in the resting state brain may be more genetically influenced in males and more environmentally influenced in females; furthermore, standard approaches may suggest that it is more substantially nonadditive genetics, rather than additive genetics, which contribute to the differences in sex-specific individual differences based on this young adult (male and female) specific cohort. Finally, considering the preliminary cohort-specific results, based on standard approaches, environmental influences on individual differences may be substantially greater than that of genetics, for either sex, frontally and brain-wide. [Correction added on 10 May 2023, after first online publication: added: functional Magnetic Resonance Imaging. Added: individual differences in, twice. Added statement between furthermore … based on standard approaches.].

Cerebral blood flow and cardiovascular risk effects on resting brain regional homogeneity

Regional homogeneity (ReHo) is a measure of local functional brain connectivity that has been reported to be altered in a wide range of neuropsychiatric disorders. Computed from brain resting-state functional MRI time series, ReHo is also sensitive to fluctuations in cerebral blood flow (CBF) that in turn may be influenced by cerebrovascular health. We accessed cerebrovascular health with Framingham cardiovascular risk score (FCVRS). We hypothesize that ReHo signal may be influenced by regional CBF; and that these associations can be summarized as FCVRS→CBF→ReHo. We used three independent samples to test this hypothesis. A test-retest sample of N = 30 healthy volunteers was used for test-retest evaluation of CBF effects on ReHo. Amish Connectome Project (ACP) sample (N = 204, healthy individuals) was used to evaluate association between FCVRS and ReHo and testing if the association diminishes given CBF. The UKBB sample (N = 6,285, healthy participants) was used to replicate the effects of FCVRS on ReHo. We observed strong CBF→ReHo links (p<2.5 × 10^-3) using a three-point longitudinal sample. In ACP sample, marginal and partial correlations analyses demonstrated that both CBF and FCVRS were significantly correlated with the whole-brain average (p<10^-6) and regional ReHo values, with the strongest correlations observed in frontal, parietal, and temporal areas. Yet, the association between ReHo and FCVRS became insignificant once the effect of CBF was accounted for. In contrast, CBF→ReHo remained significantly linked after adjusting for FCVRS and demographic covariates (p<10^-6). Analysis in N = 6,285 replicated the FCVRS→ReHo effect (p = 2.7 × 10^-27). In summary, ReHo alterations in health and neuropsychiatric illnesses may be partially driven by region-specific variability in CBF, which is, in turn, influenced by cardiovascular factors.

Interactive Effects of HIV Infection and Cannabis Use on Insula Subregion Functional Connectivity

Chronic inflammation in the central nervous system is one mechanism through which human immunodeficiency virus (HIV) may lead to progressive cognitive decline. Given cannabis's (CB's) anti-inflammatory properties, use prevalence among people living with HIV (PLWH), and evidence implicating the insula in both, we examined independent and interactive effects of HIV and CB on insular circuitry, cognition, and immune function. We assessed resting-state functional connectivity (rsFC) of three insula subregions among 106 participants across four groups (co-occurring: HIV+/CB+; HIV-only: HIV+/CB-; CB-only: HIV-/CB+; controls: HIV-/CB-). Participants completed a neurocognitive battery assessing functioning across multiple domains and self-reported somatic complaints. Blood samples quantified immune function (T-cell counts) and inflammation (tumor necrosis factor alpha [TNF-α]). We observed interactive HIV × CB effects on rsFC strength between two anterior insula (aI) subregions and sensorimotor cortices such that, CB appeared to normalize altered rsFC among non-using PLWH. Specifically, compared to controls, HIV-only and CB-only groups displayed decreased dorsal anterior insula (DI) - postcentral gyrus rsFC and increased ventral anterior insula (VI) - supplementary motor area rsFC, whereas the co-occurring group displayed DI and VI rsFC more akin to that of controls. Altered DI - postcentral rsFC correlated with decreased processing speed and somatic complaints, but did not significantly correlate with inflammation (TNF-α). These outcomes implicate insula - sensorimotor neurocircuitries in HIV and CB and are consistent with prior work suggesting that CB use may normalize insula functioning among PLWH.

Dynamic trajectories of connectome state transitions are heritable

The brain’s functional connectome is dynamic, constantly reconfiguring in an individual-specific manner. However, which characteristics of such reconfigurations are subject to genetic effects, and to what extent, is largely unknown. Here, we identified heritable dynamic features, quantified their heritability, and determined their association with cognitive phenotypes. In resting-state fMRI, we obtained multivariate features, each describing a temporal or spatial characteristic of connectome dynamics jointly over a set of connectome states. We found strong evidence for heritability of temporal features, particularly, Fractional Occupancy (FO) and Transition Probability (TP), representing the duration spent in each connectivity configuration and the frequency of shifting between configurations, respectively. These effects were robust against methodological choices of number of states and global signal regression. Genetic effects explained a substantial proportion of phenotypic variance of these features (h² = 0.39, 95% CI = [.24,.54] for FO; h² = 0. 43, 95% CI = [.29,.57] for TP). Moreover, these temporal phenotypes were associated with cognitive performance. Contrarily, we found no robust evidence for heritability of spatial features of the dynamic states (i.e., states’ Modularity and connectivity pattern). Genetic effects may therefore primarily contribute to how the connectome transitions across states, rather than the precise spatial instantiation of the states in individuals. In sum, genetic effects impact the dynamic trajectory of state transitions (captured by FO and TP), and such temporal features may act as endophenotypes for cognitive abilities.

Unsupervised classification reveals consistency and degeneracy in neural network patterns of emotion

In the present study, we used an unsupervised classification algorithm to reveal both consistency and degeneracy in neural network connectivity during anger and anxiety. Degeneracy refers to the ability of different biological pathways to produce the same outcomes. Previous research is suggestive of degeneracy in emotion, but little research has explicitly examined whether degenerate functional connectivity patterns exist for emotion categories such as anger and anxiety. Twenty-four subjects underwent functional magnetic resonance imaging (fMRI) while listening to unpleasant music and self-generating experiences of anger and anxiety. A data-driven model building algorithm with unsupervised classification (subgrouping Group Iterative Multiple Model Estimation) identified patterns of connectivity among 11 intrinsic networks that were associated with anger vs anxiety. As predicted, degenerate functional connectivity patterns existed within these overarching consistent patterns. Degenerate patterns were not attributable to differences in emotional experience or other individual-level factors. These findings are consistent with the constructionist account that emotions emerge from flexible functional neuronal assemblies and that emotion categories such as anger and anxiety each describe populations of highly variable instances.

Effects of copy number variations on brain structure and risk for psychiatric illness: Large-scale studies from the ENIGMA working groups on CNVs

The Enhancing NeuroImaging Genetics through Meta-Analysis copy number variant (ENIGMA-CNV) and 22q11.2 Deletion Syndrome Working Groups (22q-ENIGMA WGs) were created to gain insight into the involvement of genetic factors in human brain development and related cognitive, psychiatric and behavioral manifestations. To that end, the ENIGMA-CNV WG has collated CNV and magnetic resonance imaging (MRI) data from ~49,000 individuals across 38 global research sites, yielding one of the largest studies to date on the effects of CNVs on brain structures in the general population. The 22q-ENIGMA WG includes 12 international research centers that assessed over 533 individuals with a confirmed 22q11.2 deletion syndrome, 40 with 22q11.2 duplications, and 333 typically developing controls, creating the largest-ever 22q11.2 CNV neuroimaging data set. In this review, we outline the ENIGMA infrastructure and procedures for multi-site analysis of CNVs and MRI data. So far, ENIGMA has identified effects of the 22q11.2, 16p11.2 distal, 15q11.2, and 1q21.1 distal CNVs on subcortical and cortical brain structures. Each CNV is associated with differences in cognitive, neurodevelopmental and neuropsychiatric traits, with characteristic patterns of brain structural abnormalities. Evidence of gene-dosage effects on distinct brain regions also emerged, providing further insight into genotype-phenotype relationships. Taken together, these results offer a more comprehensive picture of molecular mechanisms involved in typical and atypical brain development. This "genotype-first" approach also contributes to our understanding of the etiopathogenesis of brain disorders. Finally, we outline future directions to better understand effects of CNVs on brain structure and behavior.

An overview of the first 5 years of the ENIGMA obsessive-compulsive disorder working group: The power of worldwide collaboration

Neuroimaging has played an important part in advancing our understanding of the neurobiology of obsessive-compulsive disorder (OCD). At the same time, neuroimaging studies of OCD have had notable limitations, including reliance on relatively small samples. International collaborative efforts to increase statistical power by combining samples from across sites have been bolstered by the ENIGMA consortium; this provides specific technical expertise for conducting multi-site analyses, as well as access to a collaborative community of neuroimaging scientists. In this article, we outline the background to, development of, and initial findings from ENIGMA's OCD working group, which currently consists of 47 samples from 34 institutes in 15 countries on 5 continents, with a total sample of 2,323 OCD patients and 2,325 healthy controls. Initial work has focused on studies of cortical thickness and subcortical volumes, structural connectivity, and brain lateralization in children, adolescents and adults with OCD, also including the study on the commonalities and distinctions across different neurodevelopment disorders. Additional work is ongoing, employing machine learning techniques. Findings to date have contributed to the development of neurobiological models of OCD, have provided an important model of global scientific collaboration, and have had a number of clinical implications. Importantly, our work has shed new light on questions about whether structural and functional alterations found in OCD reflect neurodevelopmental changes, effects of the disease process, or medication impacts. We conclude with a summary of ongoing work by ENIGMA-OCD, and a consideration of future directions for neuroimaging research on OCD within and beyond ENIGMA.

Where the genome meets the connectome: Understanding how genes shape human brain connectivity

The integration of modern neuroimaging methods with genetically informative designs and data can shed light on the molecular mechanisms underlying the structural and functional organization of the human connectome. Here, we review studies that have investigated the genetic basis of human brain network structure and function through three complementary frameworks: (1) the quantification of phenotypic heritability through classical twin designs; (2) the identification of specific DNA variants linked to phenotypic variation through association and related studies; and (3) the analysis of correlations between spatial variations in imaging phenotypes and gene expression profiles through the integration of neuroimaging and transcriptional atlas data. We consider the basic foundations, strengths, limitations, and discoveries associated with each approach. We present converging evidence to indicate that anatomical connectivity is under stronger genetic influence than functional connectivity and that genetic influences are not uniformly distributed throughout the brain, with phenotypic variation in certain regions and connections being under stronger genetic control than others. We also consider how the combination of imaging and genetics can be used to understand the ways in which genes may drive brain dysfunction in different clinical disorders.

Intergenerational transmission of the patterns of functional and structural brain networks

There is clear evidence of intergenerational transmission of life values, cognitive traits, psychiatric disorders, and even aspects of daily decision making. To investigate biological substrates of this phenomenon, the brain has received increasing attention as a measurable biomarker and potential target for intervention. However, no previous study has quantitatively and comprehensively investigated the effects of intergenerational transmission on functional and structural brain networks. Here, by employing an unusually large cohort dataset (N = 84 parent-child dyads; 45 sons, 39 daughters, 81 mothers, and 3 fathers), we show that patterns of functional and structural brain networks are preserved over a generation. We also demonstrate that several demographic factors and behavioral/physiological phenotypes have a relationship with brain similarity. Collectively, our results provide a comprehensive picture of neurobiological substrates of intergenerational transmission and demonstrate the usability of our dataset for investigating the neurobiological substrates of intergenerational transmission.

Genetic Influence on Nociceptive Processing in the Human Brain-A Twin Study

Nociceptive processing in the human brain is complex and involves several brain structures and varies across individuals. Determining the structures that contribute to interindividual differences in nociceptive processing is likely to improve our understanding of why some individuals feel more pain than others. Here, we found specific parts of the cerebral response to nociception that are under genetic influence by employing a classic twin-design. We found genetic influences on nociceptive processing in the midcingulate cortex and bilateral posterior insula. In addition to brain activations, we found genetic contributions to large-scale functional connectivity (FC) during nociceptive processing. We conclude that additive genetics influence specific brain regions involved in nociceptive processing. The genetic influence on FC during nociceptive processing is not limited to core nociceptive brain regions, such as the dorsal posterior insula and somatosensory areas, but also involves cognitive and affective brain circuitry. These findings improve our understanding of human pain perception and increases chances to find new treatments for clinical pain.

Inferring the heritability of large-scale functional networks with a multivariate ACE modeling approach

Recent evidence suggests that the human functional connectome is stable at different timescales and is unique. These characteristics posit the functional connectome not only as an individual marker but also as a powerful discriminatory measure characterized by high intersubject variability. Among distinct sources of intersubject variability, the long-term sources include functional patterns that emerge from genetic factors. Here, we sought to investigate the contribution of additive genetic factors to the variability of functional networks by determining the heritability of the connectivity strength in a multivariate fashion. First, we reproduced and extended the connectome fingerprinting analysis to the identification of twin pairs. Then, we estimated the heritability of functional networks by a multivariate ACE modeling approach with bootstrapping. Twin pairs were identified above chance level using connectome fingerprinting, with monozygotic twin identification accuracy equal to 57.2% on average for whole-brain connectome. Additionally, we found that a visual (0.37), the medial frontal (0.31), and the motor (0.30) functional networks were the most influenced by additive genetic factors. Our findings suggest that genetic factors not only partially determine intersubject variability of the functional connectome, such that twins can be identified using connectome fingerprinting, but also differentially influence connectivity strength in large-scale functional networks.

The functional connectome is a unique representation of the functional organization of the human brain. As such, it has been extensively used as an individual marker, a “fingerprint,” because of its high intersubject variability. Here, we sought to investigate the influence of genetic factors on intersubject variability of functional networks. Therefore, we extended the connectome fingerprinting analysis to the identification of twin pairs, and we estimated the heritability of functional networks by a multivariate ACE modeling approach with bootstrapping. We found that genetic factors not only partially determine intersubject variability of the functional connectome, such that monozygotic twin identification accuracy achieved 57.2% on average using whole-brain connectome in the fingerprinting analysis, but also differentially influence connectivity strength in large-scale functional networks.

Extended-amygdala intrinsic functional connectivity networks: A population study

Pre-clinical and human neuroimaging research implicates the extended-amygdala (ExtA) (including the bed nucleus of the stria terminalis [BST] and central nucleus of the amygdala [CeA]) in networks mediating negative emotional states associated with stress and substance-use behaviours. The extent to which individual ExtA structures form a functionally integrated unit is controversial. We utilised a large sample (n > 1,000 healthy young adult humans) to compare the intrinsic functional connectivity networks (ICNs) of the BST and CeA using task-free functional magnetic resonance imaging (fMRI) data from the Human Connectome Project. We assessed whether inter-individual differences within these ICNs were related to two principal components representing negative disposition and alcohol use. Building on recent primate evidence, we tested whether within BST-CeA intrinsic functional connectivity (iFC) was heritable and further examined co-heritability with our principal components. We demonstrate the BST and CeA to have discrete, but largely overlapping ICNs similar to previous findings. We found no evidence that within BST-CeA iFC was heritable; however, post hoc analyses found significant BST iFC heritability with the broader superficial and centromedial amygdala regions. There were no significant correlations or co-heritability associations with our principal components either across the ICNs or for specific BST-Amygdala iFC. Possible differences in phenotype associations across task-free, task-based, and clinical fMRI are discussed, along with suggestions for more causal investigative paradigms that make use of the now well-established ExtA ICNs.

Reliability modelling of resting-state functional connectivity

Resting-state functional magnetic resonance imaging (rs-fMRI) has an inherently low signal-to-noise ratio largely due to thermal and physiological noise that attenuates the functional connectivity (FC) estimates. Such attenuation limits the reliability of FC and may bias its association with other traits. Low reliability also limits heritability estimates. Classical test theory can be used to obtain a true correlation estimate free of random measurement error from parallel tests, such as split-half sessions of a rs-fMRI scan. We applied a measurement model to split-half FC estimates from the resting-state fMRI data of 1003 participants from the Human Connectome Project (HCP) to examine the benefit of reliability modelling of FC in association with traits from various domains. We evaluated the efficiency of the measurement model on extracting a stable and reliable component of FC and its association with several traits for various sample sizes and scan durations. In addition, we aimed to replicate our previous findings of increased heritability estimates when using a measurement model in a longitudinal adolescent twin cohort. The split-half measurement model improved test-retest reliability of FC on average with +0.33 points (from +0.49 to +0.82), improved strength of associations between FC and various traits on average 1.2-fold (range 1.09-1.35), and increased heritability estimates on average with +20% points (from 39% to 59%) for the full HCP dataset. On average, about half of the variance in split-session FC estimates was attributed to the stable and reliable component of FC. Shorter scan durations showed greater benefit of reliability modelling (up to 1.6-fold improvement), with an additional gain for smaller sample sizes (up to 1.8-fold improvement). Reliability modelling of FC based on a split-half using a measurement model can benefit genetic and behavioral studies by extracting a stable and reliable component of FC that is free from random measurement error and improves genetic and behavioral associations.

Local versus long-range connectivity patterns of auditory disturbance in schizophrenia

Auditory hallucinations are a debilitating symptom of schizophrenia. Effective treatment is limited because the underlying neural mechanisms remain unknown. Our study investigates how local and long-range functional connectivity is associated with auditory perceptual disturbances (APD) in schizophrenia. APD was assessed using the Auditory Perceptual Trait and State Scale. Resting state fMRI data were collected for N=99 patients with schizophrenia. Local functional connectivity was estimated using regional homogeneity (ReHo) analysis; long-range connectivity was estimated using resting state functional connectivity (rsFC) analysis. Mediation analyses tested whether local (ReHo) connectivity significantly mediated associations between long-distance rsFC and APD. Severity of APD was significantly associated with reduced ReHo in left and right putamen, left temporoparietal junction (TPJ), and right hippocampus-pallidum. Higher APD was also associated with reduced rsFC between the right putamen and the contralateral putamen and auditory cortex. Local and long-distance connectivity measures together explained 40.3% of variance in APD (P < 0.001), with the strongest predictor being the left TPJ ReHo (P < 0.001). Additionally, TPJ ReHo significantly mediated the relationship between right putamen - left putamen rsFC and APD (Sobel test, P = 0.001). Our findings suggest that both local and long-range functional connectivity deficits contribute to APD, emphasizing the role of striatum and auditory cortex. Considering the translational impact of these circuit-based findings within the context of prior clinical trials to treat auditory hallucinations, we propose a model in which correction of both local and long-distance functional connectivity deficits may be necessary to treat auditory hallucinations.

Heritability of Functional Connectivity in Resting State: Assessment of the Dynamic Mean, Dynamic Variance, and Static Connectivity across Networks

Recent efforts to evaluate the heritability of the brain's functional connectome have predominantly focused on static connectivity. However, evaluating connectivity changes across time can provide valuable insight about the inherent dynamic nature of brain function. Here, the heritability of Human Connectome Project resting-state fMRI data was examined to determine whether there is a genetic basis for dynamic fluctuations in functional connectivity. The dynamic connectivity variance, in addition to the dynamic mean and standard static connectivity, was evaluated. Heritability was estimated using Accelerated Permutation Inference for the ACE (APACE), which models the additive genetic (h2), common environmental (c2), and unique environmental (e2) variance. Heritability was moderate (mean h2: dynamic mean = 0.35, dynamic variance = 0.45, and static = 0.37) and tended to be greater for dynamic variance compared to either dynamic mean or static connectivity. Further, heritability of dynamic variance was reliable across both sessions for several network connections, particularly between higher-order cognitive and visual networks. For both dynamic mean and static connectivity, similar patterns of heritability were found across networks. The findings support the notion that dynamic connectivity is genetically influenced. The flexibility of network connections, not just their strength, is a heritable endophenotype that may predispose trait behavior.

Temporal-thalamic and cingulo-opercular connectivity in people with schizophrenia

A growing body of research has suggested that people with schizophrenia (SZ) exhibit altered patterns of functional and anatomical brain connectivity. For example, many previous resting state functional connectivity (rsFC) studies have shown that, compared to healthy controls (HC), people with SZ demonstrate hyperconnectivity between subregions of the thalamus and sensory cortices, as well as hypoconnectivity between subregions of the thalamus and prefrontal cortex. In addition to thalamic findings, hypoconnectivity between cingulo-opercular brain regions thought to be involved in salience detection has also been commonly reported in people with SZ. However, previous studies have largely relied on seed-based analyses. Seed-based approaches require researchers to define a single a priori brain region, which is then used to create a rsFC map across the entire brain. While useful for testing specific hypotheses, these analyses are limited in that only a subset of connections across the brain are explored. In the current manuscript, we leverage novel network statistical techniques in order to detect latent functional connectivity networks with organized topology that successfully differentiate people with SZ from HCs. Importantly, these techniques do not require a priori seed selection and allow for whole brain investigation, representing a comprehensive, data-driven approach to determining differential connectivity between diagnostic groups. Across two samples, (Sample 1: 35 SZ, 44 HC; Sample 2: 65 SZ, 79 HC), we found evidence for differential rsFC within a network including temporal and thalamic regions. Connectivity in this network was greater for people with SZ compared to HCs. In the second sample, we also found evidence for hypoconnectivity within a cingulo-opercular network of brain regions in people with SZ compared to HCs. In summary, our results replicate and extend previous studies suggesting hyperconnectivity between the thalamus and sensory cortices and hypoconnectivity between cingulo-opercular regions in people with SZ using data-driven statistical and graph theoretical techniques.

Test-retest reliability of the human functional connectome over consecutive days: identifying highly reliable portions and assessing the impact of methodological choices

Countless studies have advanced our understanding of the human brain and its organization by using functional magnetic resonance imaging (fMRI) to derive network representations of human brain function. However, we do not know to what extent these "functional connectomes" are reliable over time. In a large public sample of healthy participants (N = 833) scanned on two consecutive days, we assessed the test-retest reliability of fMRI functional connectivity and the consequences on reliability of three common sources of variation in analysis workflows: atlas choice, global signal regression, and thresholding. By adopting the intraclass correlation coefficient as a metric, we demonstrate that only a small portion of the functional connectome is characterized by good (6-8%) to excellent (0.08-0.14%) reliability. Connectivity between prefrontal, parietal, and temporal areas is especially reliable, but also average connectivity within known networks has good reliability. In general, while unreliable edges are weak, reliable edges are not necessarily strong. Methodologically, reliability of edges varies between atlases, global signal regression decreases reliability for networks and most edges (but increases it for some), and thresholding based on connection strength reduces reliability. Focusing on the reliable portion of the connectome could help quantify brain trait-like features and investigate individual differences using functional neuroimaging.

ENIGMA MDD: seven years of global neuroimaging studies of major depression through worldwide data sharing

A key objective in the field of translational psychiatry over the past few decades has been to identify the brain correlates of major depressive disorder (MDD). Identifying measurable indicators of brain processes associated with MDD could facilitate the detection of individuals at risk, and the development of novel treatments, the monitoring of treatment effects, and predicting who might benefit most from treatments that target specific brain mechanisms. However, despite intensive neuroimaging research towards this effort, underpowered studies and a lack of reproducible findings have hindered progress. Here, we discuss the work of the ENIGMA Major Depressive Disorder (MDD) Consortium, which was established to address issues of poor replication, unreliable results, and overestimation of effect sizes in previous studies. The ENIGMA MDD Consortium currently includes data from 45 MDD study cohorts from 14 countries across six continents. The primary aim of ENIGMA MDD is to identify structural and functional brain alterations associated with MDD that can be reliably detected and replicated across cohorts worldwide. A secondary goal is to investigate how demographic, genetic, clinical, psychological, and environmental factors affect these associations. In this review, we summarize findings of the ENIGMA MDD disease working group to date and discuss future directions. We also highlight the challenges and benefits of large-scale data sharing for mental health research.

ENIGMA and global neuroscience: A decade of large-scale studies of the brain in health and disease across more than 40 countries

This review summarizes the last decade of work by the ENIGMA (Enhancing NeuroImaging Genetics through Meta Analysis) Consortium, a global alliance of over 1400 scientists across 43 countries, studying the human brain in health and disease. Building on large-scale genetic studies that discovered the first robustly replicated genetic loci associated with brain metrics, ENIGMA has diversified into over 50 working groups (WGs), pooling worldwide data and expertise to answer fundamental questions in neuroscience, psychiatry, neurology, and genetics. Most ENIGMA WGs focus on specific psychiatric and neurological conditions, other WGs study normal variation due to sex and gender differences, or development and aging; still other WGs develop methodological pipelines and tools to facilitate harmonized analyses of "big data" (i.e., genetic and epigenetic data, multimodal MRI, and electroencephalography data). These international efforts have yielded the largest neuroimaging studies to date in schizophrenia, bipolar disorder, major depressive disorder, post-traumatic stress disorder, substance use disorders, obsessive-compulsive disorder, attention-deficit/hyperactivity disorder, autism spectrum disorders, epilepsy, and 22q11.2 deletion syndrome. More recent ENIGMA WGs have formed to study anxiety disorders, suicidal thoughts and behavior, sleep and insomnia, eating disorders, irritability, brain injury, antisocial personality and conduct disorder, and dissociative identity disorder. Here, we summarize the first decade of ENIGMA's activities and ongoing projects, and describe the successes and challenges encountered along the way. We highlight the advantages of collaborative large-scale coordinated data analyses for testing reproducibility and robustness of findings, offering the opportunity to identify brain systems involved in clinical syndromes across diverse samples and associated genetic, environmental, demographic, cognitive, and psychosocial factors.

Genetic influence on ageing-related changes in resting-state brain functional networks in healthy adults: A systematic review

This systematic review examines the genetic and epigenetic factors associated with resting-state functional connectivity (RSFC) in healthy human adult brains across the lifespan, with a focus on genes associated with Alzheimer's disease (AD). There were 58 studies included. The key findings are: (i) genetic factors have a low to moderate contribution; (ii) the apolipoprotein E ε2/3/4 polymorphism was the most studied genetic variant, with the APOE-ε4 allele most consistently associated with deficits of the default mode network, but there were insufficient studies to determine the relationships with other AD candidate risk genes; (iii) a single genome-wide association study identified several variants related to RSFC; (iv) two epigenetic independent studies showed a positive relationship between blood DNA methylation of the SLC6A4 promoter and RSFC measures. Thus, there is emerging evidence that genetic and epigenetic variation influence the brain's functional organisation and connectivity over the adult lifespan. However, more studies are required to elucidate the roles genetic and epigenetic factors play in RSFC measures across the adult lifespan.

The Role of Hippocampal Functional Connectivity on Multisystem Subclinical Abnormalities in Schizophrenia

OBJECTIVE

Schizophrenia is associated with excess medical mortality: patients have an average life expectancy one to two decades shorter than the general population. This study investigates the relationship between aberrant hippocampal resting-state functional connectivity in schizophrenia and cumulative subclinical effects of chronic stress on metabolic, cardiovascular, and immune function using the allostatic load index.

METHODS

Cumulative stress was estimated using allostatic load total score (range, 0-13) in 46 patients with schizophrenia and 31 controls matched for age and sex (patients: age = 36.1 [13.7] years, sex = 32/14 male/female; controls: age = 35.5 [14.1], sex = 21/10 male/female). Hippocampal functional connectivity was assessed using resting-state functional magnetic resonance imaging; hippocampal structural connectivity was assessed using fornix fractional anisotropy. Linear regression analysis was used a) to test the hypothesis that aberrant hippocampal resting-state functional connectivity in schizophrenia (identified in analysis of schizophrenia - control differences) is associated with elevated allostatic load scores in patients and b) to determine the association between fornix fractional anisotropy with allostatic load.

RESULTS

In patients, higher allostatic load was significantly associated with reduced resting functional connectivity between the left hippocampus and right cingulate cortex and left precentral gyrus, but higher connectivity between the right hippocampus and left cerebellum lobe VI (corrected p values <. 05). In controls, reductions in both hippocampal structural connectivity and hippocampal-cingulate functional connectivity were associated with higher allostatic load scores.

CONCLUSIONS

These findings support basic neuroscience evidence that cumulative stress and hippocampal function are closely connected and suggest that abnormal hippocampal functional communication in schizophrenia may be related to elevated multisystem subclinical medical issues in patients as indexed by allostatic load.

Effects of ketamine and midazolam on resting state connectivity and comparison with ENIGMA connectivity deficit patterns in schizophrenia

Subanesthetic administration of ketamine is a pharmacological model to elicit positive and negative symptoms of psychosis in healthy volunteers. We used resting‐state pharmacological functional MRI (rsPhfMRI) to identify cerebral networks affected by ketamine and compared them to the functional connectivity (FC) in schizophrenia. Ketamine can produce sedation and we contrasted its effects with the effects of the anxiolytic drug midazolam. Thirty healthy male volunteers (age = 19–37 years) underwent a randomized, three‐way, cross‐over study consisting of three imaging sessions, with 48 hr between sessions. A session consisted of a control period followed by infusion of placebo or ketamine or midazolam. The ENIGMA rsfMRI pipeline was used to derive two long‐distance (seed‐based and dual‐regression) and one local (regional homogeneity, ReHo) FC measures. Ketamine induced significant reductions in the connectivity of the salience network (Cohen's d: 1.13 ± 0.28, p = 4.0 × 10⁻³), auditory network (d: 0.67 ± 0.26, p = .04) and default mode network (DMN, d: 0.63 ± 0.26, p = .05). Midazolam significantly reduced connectivity in the DMN (d: 0.77 ± 0.27, p = .03). The effect sizes for ketamine for resting networks showed a positive correlation (r = .59, p = .07) with the effect sizes for schizophrenia‐related deficits derived from ENIGMA's study of 261 patients and 327 controls. Effect sizes for midazolam were not correlated with the schizophrenia pattern (r = −.17, p = .65). The subtraction of ketamine and midazolam patterns showed a significant positive correlation with the pattern of schizophrenia deficits (r = .68, p = .03). RsPhfMRI reliably detected the shared and divergent pharmacological actions of ketamine and midazolam on cerebral networks. The pattern of disconnectivity produced by ketamine was positively correlated with the pattern of connectivity deficits observed in schizophrenia, suggesting a brain functional basis for previously poorly understood effects of the drug.

Genetic and environmental influences on functional connectivity within and between canonical cortical resting-state networks throughout adolescent development in boys and girls

The human brain is active during rest and hierarchically organized into intrinsic functional networks. These functional networks are largely established early in development, with reports of a shift from a local to more distributed organization during childhood and adolescence. It remains unknown to what extent genetic and environmental influences on functional connectivity change throughout adolescent development. We measured functional connectivity within and between eight cortical networks in a longitudinal resting-state fMRI study of adolescent twins and their older siblings on two occasions (mean ages 13 and 18 years). We modelled the reliability for these inherently noisy and head-motion sensitive measurements by analyzing data from split-half sessions. Functional connectivity between resting-state networks decreased with age whereas functional connectivity within resting-state networks generally increased with age, independent of general cognitive functioning. Sex effects were sparse, with stronger functional connectivity in the default mode network for girls compared to boys, and stronger functional connectivity in the salience network for boys compared to girls. Heritability explained up to 53% of the variation in functional connectivity within and between resting-state networks, and common environment explained up to 33%. Genetic influences on functional connectivity remained stable during adolescent development. In conclusion, longitudinal age-related changes in functional connectivity within and between cortical resting-state networks are subtle but wide-spread throughout adolescence. Genes play a considerable role in explaining individual variation in functional connectivity with mostly stable influences throughout adolescence.

Functional network connectivity impairments and core cognitive deficits in schizophrenia

Cognitive deficits contribute to functional disability in patients with schizophrenia and may be related to altered functional networks that serve cognition. We evaluated the integrity of major functional networks and assessed their role in supporting two cognitive functions affected in schizophrenia: processing speed (PS) and working memory (WM). Resting-state functional magnetic resonance imaging (rsfMRI) data, N = 261 patients and 327 controls, were aggregated from three independent cohorts and evaluated using Enhancing NeuroImaging Genetics through Meta Analysis rsfMRI analysis pipeline. Meta- and mega-analyses were used to evaluate patient-control differences in functional connectivity (FC) measures. Canonical correlation analysis was used to study the association between cognitive deficits and FC measures. Patients showed consistent patterns of cognitive and resting-state FC (rsFC) deficits across three cohorts. Patient-control differences in rsFC calculated using seed-based and dual-regression approaches were consistent (Cohen's d: 0.31 ± 0.09 and 0.29 ± 0.08, p < 10^-4 ). RsFC measures explained 12-17% of the individual variations in PS and WM in the full sample and in patients and controls separately, with the strongest correlations found in salience, auditory, somatosensory, and default-mode networks. The pattern of association between rsFC (within-network) and PS (r = .45, p = .07) and WM (r = .36, p = .16), and rsFC (between-network) and PS (r = .52, p = 8.4 × 10^-3 ) and WM (r = .47, p = .02), derived from multiple networks was related to effect size of patient-control differences in the functional networks. No association was detected between rsFC and current medication dose or psychosis ratings. Patients demonstrated significant reduction in several FC networks that may partially underlie some of the core neurocognitive deficits in schizophrenia. The strength of connectivity-cognition relationships in different networks was strongly associated with network's vulnerability to schizophrenia.

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.

Comparison of heritability estimates on resting state fMRI connectivity phenotypes using the ENIGMA analysis pipeline

We measured and compared heritability estimates for measures of functional brain connectivity extracted using the Enhancing Neuroimaging Genetics through Meta-Analysis (ENIGMA) rsfMRI analysis pipeline in two cohorts: the genetics of brain structure (GOBS) cohort and the HCP (the Human Connectome Project) cohort. These two cohorts were assessed using conventional (GOBS) and advanced (HCP) rsfMRI protocols, offering a test case for harmonization of rsfMRI phenotypes, and to determine measures that show consistent heritability for in-depth genome-wide analysis. The GOBS cohort consisted of 334 Mexican-American individuals (124M/210F, average age = 47.9 ± 13.2 years) from 29 extended pedigrees (average family size = 9 people; range 5-32). The GOBS rsfMRI data was collected using a 7.5-min acquisition sequence (spatial resolution = 1.72 × 1.72 × 3 mm³ ). The HCP cohort consisted of 518 twins and family members (240M/278F; average age = 28.7 ± 3.7 years). rsfMRI data was collected using 28.8-min sequence (spatial resolution = 2 × 2 × 2 mm³ ). We used the single-modality ENIGMA rsfMRI preprocessing pipeline to estimate heritability values for measures from eight major functional networks, using (1) seed-based connectivity and (2) dual regression approaches. We observed significant heritability (h² = 0.2-0.4, p < .05) for functional connections from seven networks across both cohorts, with a significant positive correlation between heritability estimates across two cohorts. The similarity in heritability estimates for resting state connectivity measurements suggests that the additive genetic contribution to functional connectivity is robustly detectable across populations and imaging acquisition parameters. The overarching genetic influence, and means to consistently detect it, provides an opportunity to define a common genetic search space for future gene discovery studies.