Partitioning heritability analyses unveil the genetic architecture of human brain multidimensional functional connectivity patterns

Genomic structural equation modeling reveals latent phenotypes in the human cortex with distinct genetic architecture

Genetic contributions to human cortical structure manifest pervasive pleiotropy. This pleiotropy may be harnessed to identify unique genetically-informed parcellations of the cortex that are neurobiologically distinct from functional, cytoarchitectural, or other cortical parcellation schemes. We investigated genetic pleiotropy by applying genomic structural equation modeling (SEM) to map the genetic architecture of cortical surface area (SA) and cortical thickness (CT) for 34 brain regions recently reported in the ENIGMA cortical GWAS. Genomic SEM uses the empirical genetic covariance estimated from GWAS summary statistics with LD score regression (LDSC) to discover factors underlying genetic covariance, which we are denoting genetically informed brain networks (GIBNs). Genomic SEM can fit a multivariate GWAS from summary statistics for each of the GIBNs, which can subsequently be used for LD score regression (LDSC). We found the best-fitting model of cortical SA identified 6 GIBNs and CT identified 4 GIBNs, although sensitivity analyses indicated that other structures were plausible. The multivariate GWASs of the GIBNs identified 74 genome-wide significant (GWS) loci (p < 5 × 10-8), including many previously implicated in neuroimaging phenotypes, behavioral traits, and psychiatric conditions. LDSC of GIBN GWASs found that SA-derived GIBNs had a positive genetic correlation with bipolar disorder (BPD), and cannabis use disorder, indicating genetic predisposition to a larger SA in the specific GIBN is associated with greater genetic risk of these disorders. A negative genetic correlation was observed between attention deficit hyperactivity disorder (ADHD) and major depressive disorder (MDD). CT GIBNs displayed a negative genetic correlation with alcohol dependence. Even though we observed model instability in our application of genomic SEM to high-dimensional data, jointly modeling the genetic architecture of complex traits and investigating multivariate genetic links across neuroimaging phenotypes offers new insights into the genetics of cortical structure and relationships to psychopathology.

Common and distinct functional brain network abnormalities in adolescent, early-middle adult, and late adult major depressive disorders

BACKGROUND

The age-related heterogeneity in major depressive disorder (MDD) has received significant attention. However, the neural mechanisms underlying such heterogeneity still need further investigation. This study aimed to explore the common and distinct functional brain abnormalities across different age groups of MDD patients from a large-sample, multicenter analysis.

METHODS

The analyzed sample consisted of a total of 1238 individuals including 617 MDD patients (108 adolescents, 12-17 years old; 411 early-middle adults, 18-54 years old; and 98 late adults, > = 55 years old) and 621 demographically matched healthy controls (60 adolescents, 449 early-middle adults, and 112 late adults). MDD-related abnormalities in brain functional connectivity (FC) patterns were investigated in each age group separately and using the whole pooled sample, respectively.

RESULTS

We found shared FC reductions among the sensorimotor, visual, and auditory networks across all three age groups of MDD patients. Furthermore, adolescent patients uniquely exhibited increased sensorimotor-subcortical FC; early-middle adult patients uniquely exhibited decreased visual-subcortical FC; and late adult patients uniquely exhibited wide FC reductions within the subcortical, default-mode, cingulo-opercular, and attention networks. Analysis of covariance models using the whole pooled sample further revealed: (1) significant main effects of age group on FCs within most brain networks, suggesting that they are decreased with aging; and (2) a significant age group × MDD diagnosis interaction on FC within the default-mode network, which may be reflective of an accelerated aging-related decline in default-mode FCs.

CONCLUSIONS

To summarize, these findings may deepen our understanding of the age-related biological and clinical heterogeneity in MDD.

Improved estimation of general cognitive ability and its neural correlates with a large battery of cognitive tasks

Elucidating the neural mechanisms of general cognitive ability (GCA) is an important mission of cognitive neuroscience. Recent large-sample cohort studies measured GCA through multiple cognitive tasks and explored its neural basis, but they did not investigate how task number, factor models, and neural data type affect the estimation of GCA and its neural correlates. To address these issues, we tested 1,605 Chinese young adults with 19 cognitive tasks and Raven's Advanced Progressive Matrices (RAPM) and collected resting state and n-back task fMRI data from a subsample of 683 individuals. Results showed that GCA could be reliably estimated by multiple tasks. Increasing task number enhances both reliability and validity of GCA estimates and reliably strengthens their correlations with brain data. The Spearman model and hierarchical bifactor model yield similar GCA estimates. The bifactor model has better model fit and stronger correlation with RAPM but explains less variance and shows weaker correlations with brain data than does the Spearman model. Notably, the n-back task-based functional connectivity patterns outperform resting-state fMRI in predicting GCA. These results suggest that GCA derived from a multitude of cognitive tasks serves as a valid measure of general intelligence and that its neural correlates could be better characterized by task fMRI than resting-state fMRI data.

Genomic Structural Equation Modeling Reveals Latent Phenotypes in the Human Cortex with Distinct Genetic Architecture

Genetic contributions to human cortical structure manifest pervasive pleiotropy. This pleiotropy may be harnessed to identify unique genetically-informed parcellations of the cortex that are neurobiologically distinct from functional, cytoarchitectural, or other cortical parcellation schemes. We investigated genetic pleiotropy by applying genomic structural equation modeling (SEM) to map the genetic architecture of cortical surface area (SA) and cortical thickness (CT) for the 34 brain regions recently reported in the ENIGMA cortical GWAS. Genomic SEM uses the empirical genetic covariance estimated from GWAS summary statistics with LD score regression (LDSC) to discover factors underlying genetic covariance, which we are denoting genetically informed brain networks (GIBNs). Genomic SEM can fit a multivariate GWAS from summary statistics for each of the GIBNs, which can subsequently be used for LD score regression (LDSC). We found the best-fitting model of cortical SA identified 6 GIBNs and CT identified 4 GIBNs. The multivariate GWASs of these GIBNs identified 74 genome-wide significant (GWS) loci (p<5×10-8), including many previously implicated in neuroimaging phenotypes, behavioral traits, and psychiatric conditions. LDSC of GIBN GWASs found that SA-derived GIBNs had a positive genetic correlation with bipolar disorder (BPD), and cannabis use disorder, indicating genetic predisposition to a larger SA in the specific GIBN is associated with greater genetic risk of these disorders. A negative genetic correlation was observed with attention deficit hyperactivity disorder (ADHD), major depressive disorder (MDD), and insomnia, indicating genetic predisposition to a larger SA in the specific GIBN is associated with lower genetic risk of these disorders. CT GIBNs displayed a negative genetic correlation with alcohol dependence. Jointly modeling the genetic architecture of complex traits and investigating multivariate genetic links across phenotypes offers a new vantage point for mapping the cortex into genetically informed networks.

Structural architecture and brain network efficiency link polygenic scores to intelligence

Intelligence is highly heritable. Genome-wide association studies (GWAS) have shown that thousands of alleles contribute to variation in intelligence with small effect sizes. Polygenic scores (PGS), which combine these effects into one genetic summary measure, are increasingly used to investigate polygenic effects in independent samples. Whereas PGS explain a considerable amount of variance in intelligence, it is largely unknown how brain structure and function mediate this relationship. Here, we show that individuals with higher PGS for educational attainment and intelligence had higher scores on cognitive tests, larger surface area, and more efficient fiber connectivity derived by graph theory. Fiber network efficiency as well as the surface of brain areas partly located in parieto-frontal regions were found to mediate the relationship between PGS and cognitive performance. These findings are a crucial step forward in decoding the neurogenetic underpinnings of intelligence, as they identify specific regional networks that link polygenic predisposition to intelligence.

TTLL11 gene is associated with sustained attention performance and brain networks: A genome-wide association study of a healthy Chinese sample

Genetic studies on attention have mainly focused on children with attention-deficit/hyperactivity disorder (ADHD), so little systematic research has been conducted on genetic correlates of attention performance and their potential brain mechanisms among healthy individuals. The current study included a genome-wide association study (GWAS, N = 1145 healthy young adults) aimed to identify genes associated with sustained attention and an imaging genetics study (an independent sample of 483 healthy young adults) to examine any identified genes' influences on brain function. The GWAS found that TTLL11 showed genome-wide significant associations with sustained attention, with rs13298112 as the most significant SNP and the GG homozygotes showing more impulsive but also more focused responses than the A allele carriers. A retrospective examination of previously published ADHD GWAS results confirmed an un-reported, small but statistically significant effect of TTLL11 on ADHD. The imaging genetics study replicated this association and showed that the TTLL11 gene was associated with resting state activity and connectivity of the somatomoter network, and can be predicted by dorsal attention network connectivity. Specifically, the GG homozygotes showed lower brain activity, weaker brain network connectivity, and non-significant brain-attention association compared to the A allele carriers. Expression database showed that expression of this gene is enriched in the brain and that the G allele is associated with lower expression level than the A allele. These results suggest that TTLL11 may play a major role in healthy individuals' attention performance and may also contribute to the etiology of ADHD.

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.

A review of SNP heritability estimation methods

Over the past decade, statistical methods have been developed to estimate single nucleotide polymorphism (SNP) heritability, which measures the proportion of phenotypic variance explained by all measured SNPs in the data. Estimates of SNP heritability measure the degree to which the available genetic variants influence phenotypes and improve our understanding of the genetic architecture of complex phenotypes. In this article, we review the recently developed and commonly used SNP heritability estimation methods for continuous and binary phenotypes from the perspective of model assumptions and parameter optimization. We primarily focus on their capacity to handle multiple phenotypes and longitudinal measurements, their ability for SNP heritability partition and their use of individual-level data versus summary statistics. State-of-the-art statistical methods that are scalable to the UK Biobank dataset are also elucidated in detail.

Shared genetic influences on resting‐state functional networks of the brain

The amplitude of activation in brain resting state networks (RSNs), measured with resting‐state functional magnetic resonance imaging, is heritable and genetically correlated across RSNs, indicating pleiotropy. Recent univariate genome‐wide association studies (GWASs) explored the genetic underpinnings of individual variation in RSN activity. Yet univariate genomic analyses do not describe the pleiotropic nature of RSNs. In this study, we used a novel multivariate method called genomic structural equation modeling to model latent factors that capture the shared genomic influence on RSNs and to identify single nucleotide polymorphisms (SNPs) and genes driving this pleiotropy. Using summary statistics from GWAS of 21 RSNs reported in UK Biobank (N = 31,688), the genomic latent factor analysis was first conducted in a discovery sample (N = 21,081), and then tested in an independent sample from the same cohort (N = 10,607). In the discovery sample, we show that the genetic organization of RSNs can be best explained by two distinct but correlated genetic factors that divide multimodal association networks and sensory networks. Eleven of the 17 factor loadings were replicated in the independent sample. With the multivariate GWAS, we found and replicated nine independent SNPs associated with the joint architecture of RSNs. Further, by combining the discovery and replication samples, we discovered additional SNP and gene associations with the two factors of RSN amplitude. We conclude that modeling the genetic effects on brain function in a multivariate way is a powerful approach to learn more about the biological mechanisms involved in brain function.

Twin studies to GWAS: there and back again

The field of human behavioral genetics has come full circle. It began by using twin/family studies to estimate the relative importance of genetic and environmental influences. As large-scale genotyping became cost-effective, genome-wide association studies (GWASs) yielded insights about the nature of genetic influences and new methods that use GWAS data to estimate heritability and genetic correlations invigorated the field. Yet these newer GWAS methods have not replaced twin/family studies. In this review, we discuss the strengths and weaknesses of the two approaches with respect to characterizing genetic and environmental influences, measurement of behavioral phenotypes, and evaluation of causal models, with a particular focus on cognitive neuroscience. This discussion highlights how twin/family studies and GWAS complement and mutually reinforce one another.

Inter-annual and spatial climatic variability have led to a balance between local fluctuating selection and wide-range directional selection in a perennial grass species

BACKGROUND AND AIMS

The persistence of a plant population under a specific local climatic regime requires phenotypic adaptation with underlying particular combinations of alleles at adaptive loci. The level of allele diversity at adaptive loci within a natural plant population conditions its potential to evolve, notably towards adaptation to a change in climate. Investigating the environmental factors that contribute to the maintenance of adaptive diversity in populations is thus worthwhile. Within-population allele diversity at adaptive loci can be partly driven by the mean climate at the population site but also by its temporal variability.

METHODS

The effects of climate temporal mean and variability on within-population allele diversity at putatively adaptive quantitative trait loci (QTLs) were evaluated using 385 natural populations of Lolium perenne (perennial ryegrass) collected right across Europe. For seven adaptive traits related to reproductive phenology and vegetative potential growth seasonality, the average within-population allele diversity at major QTLs (HeA) was computed.

KEY RESULTS

Significant relationships were found between HeA of these traits and the temporal mean and variability of the local climate. These relationships were consistent with functional ecology theory.

CONCLUSIONS

Results indicated that temporal variability of local climate has likely led to fluctuating directional selection, which has contributed to the maintenance of allele diversity at adaptive loci and thus potential for further adaptation.

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.

The coupling of BOLD signal variability and degree centrality underlies cognitive functions and psychiatric diseases

Brain signal variability has been consistently linked to functional integration; however, whether this coupling is associated with cognitive functions and/or psychiatric diseases has not been clarified. Using multiple multimodality datasets, including resting-state functional magnetic resonance imaging (rsfMRI) data from the Human Connectome Project (HCP: N = 927) and a Beijing sample (N = 416) and cerebral blood flow (CBF) and rsfMRI data from a Hangzhou sample (N = 29), we found that, compared with the existing variability measure (i.e., SD_BOLD), the mean-scaled (standardized) fractional standard deviation of the BOLD signal (mfSD_BOLD) maintained very high test-retest reliability, showed greater cross-site reliability and was less affected by head motion. We also found strong reproducible couplings between the mfSD_BOLD and functional integration measured by the degree centrality (DC), both cross-voxel and cross-subject, which were robust to scanning and preprocessing parameters. Moreover, both mfSD_BOLD and DC were correlated with CBF, suggesting a common physiological basis for both measures. Critically, the degree of coupling between mfSD_BOLD and long-range DC was positively correlated with individuals' cognitive total composite scores. Brain regions with greater mismatches between mfSD_BOLD and long-range DC were more vulnerable to brain diseases. Our results suggest that BOLD signal variability could serve as a meaningful index of local function that underlies functional integration in the human brain and that a strong coupling between BOLD signal variability and functional integration may serve as a hallmark of balanced brain networks that are associated with optimal brain functions.

Partitioning heritability analyses unveil the genetic architecture of human brain multidimensional functional connectivity patterns

Resting-state functional connectivity profiles have been increasingly shown to be important endophenotypes that are tightly linked to human cognitive functions and psychiatric diseases, yet the genetic architecture of this multidimensional trait is barely understood. Using a unique sample of 1,704 unrelated, young and healthy Chinese Han individuals, we revealed a significant heritability of functional connectivity patterns in the whole brain and several subnetworks. We further proposed a partitioned heritability analysis for multidimensional functional connectivity patterns, which revealed the common and unique enrichment patterns of the genetic contributions to brain connectivity patterns for several gene sets linked to brain functions, including the genes expressed preferentially in the central nervous system and those associated with intelligence, educational attainment, attention-deficit/hyperactivity disorder, and schizophrenia. These results for the first time reveal the genetic architecture of multidimensional brain connectivity patterns across different networks and advance our understanding of the complex relationship between gene sets, neural networks, and behaviors.