On the interpretation of transcriptome-wide association studies.

Transcriptome-wide association studies (TWAS) aim to detect relationships between gene expression and a phenotype, and are commonly used for secondary analysis of genome-wide association study (GWAS) results. Results from TWAS analyses are often interpreted as indicating a genetic relationship between gene expression and a phenotype, but this interpretation is not consistent with the null hypothesis that is evaluated in the traditional TWAS framework. In this study we provide a mathematical outline of this TWAS framework, and elucidate what interpretations are warranted given the null hypothesis it actually tests. We then use both simulations and real data analysis to assess the implications of misinterpreting TWAS results as indicative of a genetic relationship between gene expression and the phenotype. Our simulation results show considerably inflated type 1 error rates for TWAS when interpreted this way, with 41% of significant TWAS associations detected in the real data analysis found to have insufficient statistical evidence to infer such a relationship. This demonstrates that in current implementations, TWAS cannot reliably be used to investigate genetic relationships between gene expression and a phenotype, but that local genetic correlation analysis can serve as a potential alternative.

The statistical properties of gene-set analysis.

The rapid increase in loci discovered in genome-wide association studies has created a need to understand the biological implications of these results. Gene-set analysis provides a means of gaining such understanding, but the statistical properties of gene-set analysis are not well understood, which compromises our ability to interpret its results. In this Analysis article, we provide an extensive statistical evaluation of the core structure that is inherent to all gene- set analyses and we examine current implementations in available tools. We show which factors affect valid and successful detection of gene sets and which provide a solid foundation for performing and interpreting gene-set analysis.

Meta-analysis of up to 622,409 individuals identifies 40 novel smoking behaviour associated genetic loci.

Smoking is a major heritable and modifiable risk factor for many diseases, including cancer, common respiratory disorders and cardiovascular diseases. Fourteen genetic loci have previously been associated with smoking behaviour-related traits. We tested up to 235,116 single nucleotide variants (SNVs) on the exome-array for association with smoking initiation, cigarettes per day, pack-years, and smoking cessation in a fixed effects meta-analysis of up to 61 studies (up to 346,813 participants). In a subset of 112,811 participants, a further one million SNVs were also genotyped and tested for association with the four smoking behaviour traits. SNV-trait associations with P < 5 × 10-8 in either analysis were taken forward for replication in up to 275,596 independent participants from UK Biobank. Lastly, a meta-analysis of the discovery and replication studies was performed. Sixteen SNVs were associated with at least one of the smoking behaviour traits (P < 5 × 10-8) in the discovery samples. Ten novel SNVs, including rs12616219 near TMEM182, were followed-up and five of them (rs462779 in REV3L, rs12780116 in CNNM2, rs1190736 in GPR101, rs11539157 in PJA1, and rs12616219 near TMEM182) replicated at a Bonferroni significance threshold (P < 4.5 × 10-3) with consistent direction of effect. A further 35 SNVs were associated with smoking behaviour traits in the discovery plus replication meta-analysis (up to 622,409 participants) including a rare SNV, rs150493199, in CCDC141 and two low-frequency SNVs in CEP350 and HDGFRP2. Functional follow-up implied that decreased expression of REV3L may lower the probability of smoking initiation. The novel loci will facilitate understanding the genetic aetiology of smoking behaviour and may lead to the identification of potential drug targets for smoking prevention and/or cessation.

Synaptic and brain-expressed gene sets relate to the shared genetic risk across five psychiatric disorders.

BACKGROUND: Mounting evidence shows genetic overlap between multiple psychiatric disorders. However, the biological underpinnings of shared risk for psychiatric disorders are not yet fully uncovered. The identification of underlying biological mechanisms is crucial for the progress in the treatment of these disorders. METHODS: We applied gene-set analysis including 7372 gene sets, and 53 tissue-type specific gene-expression profiles to identify sets of genes that are involved in the etiology of multiple psychiatric disorders. We included genome-wide meta-association data of the five psychiatric disorders schizophrenia, bipolar disorder, major depressive disorder, autism spectrum disorder, and attention-deficit/hyperactivity disorder. The total dataset contained 159 219 cases and 262 481 controls. RESULTS: We identified 19 gene sets that were significantly associated with the five psychiatric disorders combined, of which we excluded five sets because their associations were likely driven by schizophrenia only. Conditional analyses showed independent effects of several gene sets that in particular relate to the synapse. In addition, we found independent effects of gene expression levels in the cerebellum and frontal cortex. CONCLUSIONS: We obtained novel evidence for shared biological mechanisms that act across psychiatric disorders and we showed that several gene sets that have been related to individual disorders play a role in a broader range of psychiatric disorders.

Common genetic variants associated with cognitive performance identified using the proxy-phenotype method.

We identify common genetic variants associated with cognitive performance using a two-stage approach, which we call the proxy-phenotype method. First, we conduct a genome-wide association study of educational attainment in a large sample (n = 106,736), which produces a set of 69 education-associated SNPs. Second, using independent samples (n = 24,189), we measure the association of these education-associated SNPs with cognitive performance. Three SNPs (rs1487441, rs7923609, and rs2721173) are significantly associated with cognitive performance after correction for multiple hypothesis testing. In an independent sample of older Americans (n = 8,652), we also show that a polygenic score derived from the education-associated SNPs is associated with memory and absence of dementia. Convergent evidence from a set of bioinformatics analyses implicates four specific genes (KNCMA1, NRXN1, POU2F3, and SCRT). All of these genes are associated with a particular neurotransmitter pathway involved in synaptic plasticity, the main cellular mechanism for learning and memory.

Involvement of astrocyte and oligodendrocyte gene sets in migraine.

BACKGROUND: Migraine is a common episodic brain disorder characterized by recurrent attacks of severe unilateral headache and additional neurological symptoms. Two main migraine types can be distinguished based on the presence of aura symptoms that can accompany the headache: migraine with aura and migraine without aura. Multiple genetic and environmental factors confer disease susceptibility. Recent genome-wide association studies (GWAS) indicate that migraine susceptibility genes are involved in various pathways, including neurotransmission, which have already been implicated in genetic studies of monogenic familial hemiplegic migraine, a subtype of migraine with aura. METHODS: To further explore the genetic background of migraine, we performed a gene set analysis of migraine GWAS data of 4954 clinic-based patients with migraine, as well as 13,390 controls. Curated sets of synaptic genes and sets of genes predominantly expressed in three glial cell types (astrocytes, microglia and oligodendrocytes) were investigated. DISCUSSION: Our results show that gene sets containing astrocyte- and oligodendrocyte-related genes are associated with migraine, which is especially true for gene sets involved in protein modification and signal transduction. Observed differences between migraine with aura and migraine without aura indicate that both migraine types, at least in part, seem to have a different genetic background.

MAGMA: generalized gene-set analysis of GWAS data.

By aggregating data for complex traits in a biologically meaningful way, gene and gene-set analysis constitute a valuable addition to single-marker analysis. However, although various methods for gene and gene-set analysis currently exist, they generally suffer from a number of issues. Statistical power for most methods is strongly affected by linkage disequilibrium between markers, multi-marker associations are often hard to detect, and the reliance on permutation to compute p-values tends to make the analysis computationally very expensive. To address these issues we have developed MAGMA, a novel tool for gene and gene-set analysis. The gene analysis is based on a multiple regression model, to provide better statistical performance. The gene-set analysis is built as a separate layer around the gene analysis for additional flexibility. This gene-set analysis also uses a regression structure to allow generalization to analysis of continuous properties of genes and simultaneous analysis of multiple gene sets and other gene properties. Simulations and an analysis of Crohn's Disease data are used to evaluate the performance of MAGMA and to compare it to a number of other gene and gene-set analysis tools. The results show that MAGMA has significantly more power than other tools for both the gene and the gene-set analysis, identifying more genes and gene sets associated with Crohn's Disease while maintaining a correct type 1 error rate. Moreover, the MAGMA analysis of the Crohn's Disease data was found to be considerably faster as well.

Shifts in dominance of benthic communities along a gradient of water temperature and turbidity in tropical coastal ecosystems.

Tropical coastal benthic communities will change in species composition and relative dominance due to global (e.g., increasing water temperature) and local (e.g., increasing terrestrial influence due to land-based activity) stressors. This study aimed to gain insight into possible trajectories of coastal benthic assemblages in Raja Ampat, Indonesia, by studying coral reefs at varying distances from human activities and marine lakes with high turbidity in three temperature categories (<31 °C, 31-32 °C, and >32 °C). The benthic community diversity and relative coverage of major benthic groups were quantified via replicate photo transects. The composition of benthic assemblages varied significantly among the reef and marine lake habitats. The marine lakes <31 °C contained hard coral, crustose coralline algae (CCA), and turf algae with coverages similar to those found in the coral reefs (17.4-18.8% hard coral, 3.5-26.3% CCA, and 15-15.5% turf algae, respectively), while the higher temperature marine lakes (31-32 °C and >32 °C) did not harbor hard coral or CCA. Benthic composition in the reefs was significantly influenced by geographic distance among sites but not by human activity or depth. Benthic composition in the marine lakes appeared to be structured by temperature, salinity, and degree of connection to the adjacent sea. Our results suggest that beyond a certain temperature (>31 °C), benthic communities shift away from coral dominance, but new outcomes of assemblages can be highly distinct, with a possible varied dominance of macroalgae, benthic cyanobacterial mats, or filter feeders such as bivalves and tubeworms. This study illustrates the possible use of marine lake model systems to gain insight into shifts in the benthic community structure of tropical coastal ecosystems if hard corals are no longer dominant.

Refining the scope of genetic influences on alcohol misuse through environmental stratification and gene-environment interaction.

BACKGROUND: Gene-environment interaction (G × E) is likely an important influence shaping individual differences in alcohol misuse (AM), yet it has not been extensively studied in molecular genetic research. In this study, we use a series of genome-wide gene-environment interaction (GWEIS) and in silico annotation methods with the aim of improving gene identification and biological understanding of AM. METHODS: We carried out GWEIS for four AM phenotypes in the large UK Biobank sample (N = 360,314), with trauma exposure and socioeconomic status (SES) as moderators of the genetic effects. Exploratory analyses compared stratified genome-wide association (GWAS) and GWEIS modeling approaches. We applied functional annotation, gene- and gene-set enrichment, and polygenic score analyses to interpret the GWEIS results. RESULTS: GWEIS models showed few genetic variants with significant interaction effects across gene-environment pairs. Enrichment analyses identified moderation by SES of the genes NOXA1, DLGAP1, and UBE2L3 on drinking quantity and the gene IFIT1B on drinking frequency. Except for DLGAP1, these genes have not previously been linked to AM. The most robust results (GWEIS interaction p = 4.59e-09) were seen for SES moderating the effects of variants linked to immune-related genes on a pattern of drinking with versus without meals. CONCLUSIONS: Our results highlight several genes and a potential mechanism of immune system functioning behind the moderating effect of SES on the genetic influences on AM. Although GWEIS seems to be a preferred approach over stratified GWAS, modeling G × E effects at the molecular level remains a challenge even in large samples. Understanding these effects will require substantial effort and more in-depth phenotypic measurement.

The Genetic Architecture of the Human Corpus Callosum and its Subregions.

The corpus callosum (CC) is the largest set of white matter fibers connecting the two hemispheres of the brain. In humans, it is essential for coordinating sensorimotor responses, performing associative/executive functions, and representing information in multiple dimensions. Understanding which genetic variants underpin corpus callosum morphometry, and their shared influence on cortical structure and susceptibility to neuropsychiatric disorders, can provide molecular insights into the CC's role in mediating cortical development and its contribution to neuropsychiatric disease. To characterize the morphometry of the midsagittal corpus callosum, we developed a publicly available artificial intelligence based tool to extract, parcellate, and calculate its total and regional area and thickness. Using the UK Biobank (UKB) and the Adolescent Brain Cognitive Development study (ABCD), we extracted measures of midsagittal corpus callosum morphometry and performed a genome-wide association study (GWAS) meta-analysis of European participants (combined N = 46,685). We then examined evidence for generalization to the non-European participants of the UKB and ABCD cohorts (combined N = 7,040). Post-GWAS analyses implicate prenatal intracellular organization and cell growth patterns, and high heritability in regions of open chromatin, suggesting transcriptional activity regulation in early development. Results suggest programmed cell death mediated by the immune system drives the thinning of the posterior body and isthmus. Global and local genetic overlap, along with causal genetic liability, between the corpus callosum, cerebral cortex, and neuropsychiatric disorders such as attention-deficit/hyperactivity and bipolar disorders were identified. These results provide insight into variability of corpus callosum development, its genetic influence on the cerebral cortex, and biological mechanisms related to neuropsychiatric dysfunction.

Improved functional mapping of complex trait heritability with GSA-MiXeR implicates biologically specific gene sets.

While genome-wide association studies are increasingly successful in discovering genomic loci associated with complex human traits and disorders, the biological interpretation of these findings remains challenging. Here we developed the GSA-MiXeR analytical tool for gene set analysis (GSA), which fits a model for the heritability of individual genes, accounting for linkage disequilibrium across variants and allowing the quantification of partitioned heritability and fold enrichment for small gene sets. We validated the method using extensive simulations and sensitivity analyses. When applied to a diverse selection of complex traits and disorders, including schizophrenia, GSA-MiXeR prioritizes gene sets with greater biological specificity compared to standard GSA approaches, implicating voltage-gated calcium channel function and dopaminergic signaling for schizophrenia. Such biologically relevant gene sets, often with fewer than ten genes, are more likely to provide insights into the pathobiology of complex diseases and highlight potential drug targets.

Rare variant aggregation in 148,508 exomes identifies genes associated with proxy dementia.

Proxy phenotypes allow for the utilization of genetic data from large population cohorts to analyze late-onset diseases by using parental diagnoses as a proxy for genetic disease risk. Proxy phenotypes based on parental diagnosis status have been used in previous studies to identify common variants associated with Alzheimer's disease. As of yet, proxy phenotypes have not been used to identify genes associated with Alzheimer's disease through rare variants. Here we show that a proxy Alzheimer's disease/dementia phenotype can capture known Alzheimer's disease risk genes through rare variant aggregation. We generated a proxy Alzheimer's disease/dementia phenotype for 148,508 unrelated individuals of European ancestry in the UK biobank in order to perform exome-wide rare variant aggregation analyses to identify genes associated with proxy Alzheimer's disease/dementia. We identified four genes significantly associated with the proxy phenotype, three of which were significantly associated with proxy Alzheimer's disease/dementia in an independent replication cohort consisting of 197,506 unrelated individuals of European ancestry in the UK biobank. All three of the replicated genes have been previously associated with clinically diagnosed Alzheimer's disease (SORL1, TREM2, and TOMM40/APOE). We show that proxy Alzheimer's disease/dementia can be used to identify genes associated with Alzheimer's disease through rare variant aggregation.

The Genetic Architectures of Functional and Structural Connectivity Properties within Cerebral Resting-State Networks.

Functional connectivity within resting-state networks (RSN-FC) is vital for cognitive functioning. RSN-FC is heritable and partially translates to the anatomic architecture of white matter, but the genetic component of structural connections of RSNs (RSN-SC) and their potential genetic overlap with RSN-FC remain unknown. Here, we perform genome-wide association studies (N discovery = 24,336; N replication = 3412) and annotation on RSN-SC and RSN-FC. We identify genes for visual network-SC that are involved in axon guidance and synaptic functioning. Genetic variation in RSN-FC impacts biological processes relevant to brain disorders that previously were only phenotypically associated with RSN-FC alterations. Correlations of the genetic components of RSNs are mostly observed within the functional domain, whereas less overlap is observed within the structural domain and between the functional and structural domains. This study advances the understanding of the complex functional organization of the brain and its structural underpinnings from a genetics viewpoint.

Sponge diversification in marine lakes: Implications for phylogeography and population genomic studies on sponges.

The relative influence of geography, currents, and environment on gene flow within sessile marine species remains an open question. Detecting subtle genetic differentiation at small scales is challenging in benthic populations due to large effective population sizes, general lack of resolution in genetic markers, and because barriers to dispersal often remain elusive. Marine lakes can circumvent confounding factors by providing discrete and replicated ecosystems. Using high-resolution double digest restriction-site-associated DNA sequencing (4826 Single Nucleotide Polymorphisms, SNPs), we genotyped populations of the sponge Suberites diversicolor (n = 125) to test the relative importance of spatial scales (1-1400 km), local environmental conditions, and permeability of seascape barriers in shaping population genomic structure. With the SNP dataset, we show strong intralineage population structure, even at scales <10 km (average F ST = 0.63), which was not detected previously using single markers. Most variation was explained by differentiation between populations (AMOVA: 48.8%) with signatures of population size declines and bottlenecks per lake. Although the populations were strongly structured, we did not detect significant effects of geographic distance, local environments, or degree of connection to the sea on population structure, suggesting mechanisms such as founder events with subsequent priority effects may be at play. We show that the inclusion of morphologically cryptic lineages that can be detected with the COI marker can reduce the obtained SNP set by around 90%. Future work on sponge genomics should confirm that only one lineage is included. Our results call for a reassessment of poorly dispersing benthic organisms that were previously assumed to be highly connected based on low-resolution markers.

Understanding the assumptions underlying Mendelian randomization.

With the rapidly increasing availability of large genetic data sets in recent years, Mendelian Randomization (MR) has quickly gained popularity as a novel secondary analysis method. Leveraging genetic variants as instrumental variables, MR can be used to estimate the causal effects of one phenotype on another even when experimental research is not feasible, and therefore has the potential to be highly informative. It is dependent on strong assumptions however, often producing biased results if these are not met. It is therefore imperative that these assumptions are well-understood by researchers aiming to use MR, in order to evaluate their validity in the context of their analyses and data. The aim of this perspective is therefore to further elucidate these assumptions and the role they play in MR, as well as how different kinds of data can be used to further support them.

An integrated framework for local genetic correlation analysis.

Genetic correlation (rg) analysis is used to identify phenotypes that may have a shared genetic basis. Traditionally, rg is studied globally, considering only the average of the shared signal across the genome, although this approach may fail when the rg is confined to particular genomic regions or in opposing directions at different loci. Current tools for local rg analysis are restricted to analysis of two phenotypes. Here we introduce LAVA, an integrated framework for local rg analysis that, in addition to testing the standard bivariate local rgs between two phenotypes, can evaluate local heritabilities and analyze conditional genetic relations between several phenotypes using partial correlation and multiple regression. Applied to 25 behavioral and health phenotypes, we show considerable heterogeneity in the bivariate local rgs across the genome, which is often masked by the global rg patterns, and demonstrate how our conditional approaches can elucidate more complex, multivariate genetic relations.

Disentangling Genetic Risks for Metabolic Syndrome.

A quarter of the world's population is estimated to meet the criteria for metabolic syndrome (MetS), a cluster of cardiometabolic risk factors that promote development of coronary artery disease and type 2 diabetes, leading to increased risk of premature death and significant health costs. In this study we investigate whether the genetics associated with MetS components mirror their phenotypic clustering. A multivariate approach that leverages genetic correlations of fasting glucose, HDL cholesterol, systolic blood pressure, triglycerides, and waist circumference was used, which revealed that these genetic correlations are best captured by a genetic one factor model. The common genetic factor genome-wide association study (GWAS) detects 235 associated loci, 174 more than the largest GWAS on MetS to date. Of these loci, 53 (22.5%) overlap with loci identified for two or more MetS components, indicating that MetS is a complex, heterogeneous disorder. Associated loci harbor genes that show increased expression in the brain, especially in GABAergic and dopaminergic neurons. A polygenic risk score drafted from the MetS factor GWAS predicts 5.9% of the variance in MetS. These results provide mechanistic insights into the genetics of MetS and suggestions for drug targets, especially fenofibrate, which has the promise of tackling multiple MetS components.

Exploring the genetic overlap between twelve psychiatric disorders.

The widespread comorbidity among psychiatric disorders demonstrated in epidemiological studies1-5 is mirrored by non-zero, positive genetic correlations from large-scale genetic studies6-10. To identify shared biological processes underpinning this observed phenotypic and genetic covariance and enhance molecular characterization of general psychiatric disorder liability11-13, we used several strategies aimed at uncovering pleiotropic, that is, cross-trait-associated, single-nucleotide polymorphisms (SNPs), genes and biological pathways. We conducted cross-trait meta-analysis on 12 psychiatric disorders to identify pleiotropic SNPs. The meta-analytic signal was driven by schizophrenia, hampering interpretation and joint biological characterization of the cross-trait meta-analytic signal. Subsequent pairwise comparisons of psychiatric disorders identified substantial pleiotropic overlap, but mainly among pairs of psychiatric disorders, and mainly at less stringent P-value thresholds. Only annotations related to evolutionarily conserved genomic regions were significant for multiple (9 out of 12) psychiatric disorders. Overall, identification of shared biological mechanisms remains challenging due to variation in power and genetic architecture between psychiatric disorders.

Genome-wide association study of cerebellar volume provides insights into heritable mechanisms underlying brain development and mental health.

Cerebellar volume is highly heritable and associated with neurodevelopmental and neurodegenerative disorders. Understanding the genetic architecture of cerebellar volume may improve our insight into these disorders. This study aims to investigate the convergence of cerebellar volume genetic associations in close detail. A genome-wide associations study for cerebellar volume was performed in a discovery sample of 27,486 individuals from UK Biobank, resulting in 30 genome-wide significant loci and a SNP heritability of 39.82%. We pinpoint the likely causal variants and those that have effects on amino acid sequence or cerebellar gene-expression. Additionally, 85 genome-wide significant genes were detected and tested for convergence onto biological pathways, cerebellar cell types, human evolutionary genes or developmental stages. Local genetic correlations between cerebellar volume and neurodevelopmental and neurodegenerative disorders reveal shared loci with Parkinson's disease, Alzheimer's disease and schizophrenia. These results provide insights into the heritable mechanisms that contribute to developing a brain structure important for cognitive functioning and mental health.

Sensitive period-regulating genetic pathways and exposure to adversity shape risk for depression.

Animal and human studies have documented the existence of developmental windows (or sensitive periods) when experience can have lasting effects on brain structure or function, behavior, and disease. Although sensitive periods for depression likely arise through a complex interplay of genes and experience, this possibility has not yet been explored in humans. We examined the effect of genetic pathways regulating sensitive periods, alone and in interaction with common childhood adversities, on depression risk. Guided by a translational approach, we: (1) performed association analyses of three gene sets (60 genes) shown in animal studies to regulate sensitive periods using summary data from a genome-wide association study of depression (n = 807,553); (2) evaluated the developmental expression patterns of these genes using data from BrainSpan (n = 31), a transcriptional atlas of postmortem brain samples; and (3) tested gene-by-development interplay (dGxE) by analyzing the combined effect of common variants in sensitive period genes and time-varying exposure to two types of childhood adversity within a population-based birth cohort (n = 6254). The gene set regulating sensitive period opening associated with increased depression risk. Notably, 6 of the 15 genes in this set showed developmentally regulated gene-level expression. We also identified a statistical interaction between caregiver physical or emotional abuse during ages 1-5 years and genetic risk for depression conferred by the opening genes. Genes involved in regulating sensitive periods are differentially expressed across the life course and may be implicated in depression vulnerability. Our findings about gene-by-development interplay motivate further research in large, more diverse samples to further unravel the complexity of depression etiology through a sensitive period lens.