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.

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.

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.

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.

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.

Genome-wide gene-environment interactions in neuroticism: an exploratory study across 25 environments.

Gene-environment interactions (GxE) are often suggested to play an important role in the aetiology of psychiatric phenotypes, yet so far, only a handful of genome-wide environment interaction studies (GWEIS) of psychiatric phenotypes have been conducted. Representing the most comprehensive effort of its kind to date, we used data from the UK Biobank to perform a series of GWEIS for neuroticism across 25 broadly conceptualised environmental risk factors (trauma, social support, drug use, physical health). We investigated interactions on the level of SNPs, genes, and gene-sets, and computed interaction-based polygenic risk scores (PRS) to predict neuroticism in an independent sample subset (N = 10,000). We found that the predictive ability of the interaction-based PRSs did not significantly improve beyond that of a traditional PRS based on SNP main effects from GWAS, but detected one variant and two gene-sets showing significant interaction signal after correction for the number of analysed environments. This study illustrates the possibilities and limitations of a comprehensive GWEIS in currently available sample sizes.

Systematic assessment of variability in the proteome of iPSC derivatives.

The use of induced pluripotent stem cells (iPSC) to model human complex diseases is gaining popularity as it allows investigation of human cells that are otherwise sparsely available. However, due to its laborious and cost intensive nature, iPSC research is often plagued by limited sample size and putative large variability between clones, decreasing statistical power for detecting experimental effects. Here, we investigate the source and magnitude of variability in the proteome of parallel differentiated astrocytes using mass spectrometry. We compare three possible sources of variability: inter-donor variability, inter- and intra-clonal variability, at different stages of maturation. We show that the interclonal variability is significantly smaller than the inter-donor variability, and that including more donors has a much larger influence on statistical power than adding more clones per donor. Our results provide insight into the sources of variability at protein level between iPSC samples derived in parallel and will aid in optimizing iPSC studies.

Author Correction: Genetic mapping of cell type specificity for complex traits.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Genome-wide meta-analysis of brain volume identifies genomic loci and genes shared with intelligence.

The phenotypic correlation between human intelligence and brain volume (BV) is considerable (r ≈ 0.40), and has been shown to be due to shared genetic factors. To further examine specific genetic factors driving this correlation, we present genomic analyses of the genetic overlap between intelligence and BV using genome-wide association study (GWAS) results. First, we conduct a large BV GWAS meta-analysis (N = 47,316 individuals), followed by functional annotation and gene-mapping. We identify 18 genomic loci (14 not previously associated), implicating 343 genes (270 not previously associated) and 18 biological pathways for BV. Second, we use an existing GWAS for intelligence (N = 269,867 individuals), and estimate the genetic correlation (rg) between BV and intelligence to be 0.24. We show that the rg is partly attributable to physical overlap of GWAS hits in 5 genomic loci. We identify 92 shared genes between BV and intelligence, which are mainly involved in signaling pathways regulating cell growth. Out of these 92, we prioritize 32 that are most likely to have functional impact. These results provide information on the genetics of BV and provide biological insight into BV's shared genetic etiology with intelligence.

Recognizing peripheral ecosystems in marine protected areas: A case study of golden jellyfish lakes in Raja Ampat, Indonesia.

Peripheral marine ecosystems can harbor endemic diversity and attract tourism attention, yet are generally not included in conservation management plans due to their remoteness or inland positioning. A case study in Raja Ampat of seven landlocked marine lakes containing golden jellyfish (Mastigias spp.) was conducted to address the lack of fundamental insights into evolutionary, ecological and social contexts of these ecosystems. An interdisciplinary approach was taken towards identifying the jellyfish lakes as distinct management units in order to incorporate them into existing Marine Protected Areas. Mastigias papua populations showed strong genetic (ϕST: 0.30-0.86) and morphological (F = 28.62, p-value = 0.001) structure among lakes, with putative new subspecies. Risks arising from rapid increase in tourism to Raja Ampat (30-fold since 2007) warrant restrictions on jellyfish lake use. Recommendations are provided for adaptive management and science-based conservation policies for jellyfish lakes across Indonesia.

Genetic mapping of cell type specificity for complex traits.

Single-cell RNA sequencing (scRNA-seq) data allows to create cell type specific transcriptome profiles. Such profiles can be aligned with genome-wide association studies (GWASs) to implicate cell type specificity of the traits. Current methods typically rely only on a small subset of available scRNA-seq datasets, and integrating multiple datasets is hampered by complex batch effects. Here we collated 43 publicly available scRNA-seq datasets. We propose a 3-step workflow with conditional analyses within and between datasets, circumventing batch effects, to uncover associations of traits with cell types. Applying this method to 26 traits, we identify independent associations of multiple cell types. These results lead to starting points for follow-up functional studies aimed at gaining a mechanistic understanding of these traits. The proposed framework as well as the curated scRNA-seq datasets are made available via an online platform, FUMA, to facilitate rapid evaluation of cell type specificity by other researchers.

Genome-wide analysis of insomnia in 1,331,010 individuals identifies new risk loci and functional pathways.

Insomnia is the second most prevalent mental disorder, with no sufficient treatment available. Despite substantial heritability, insight into the associated genes and neurobiological pathways remains limited. Here, we use a large genetic association sample (n = 1,331,010) to detect novel loci and gain insight into the pathways, tissue and cell types involved in insomnia complaints. We identify 202 loci implicating 956 genes through positional, expression quantitative trait loci, and chromatin mapping. The meta-analysis explained 2.6% of the variance. We show gene set enrichments for the axonal part of neurons, cortical and subcortical tissues, and specific cell types, including striatal, hypothalamic, and claustrum neurons. We found considerable genetic correlations with psychiatric traits and sleep duration, and modest correlations with other sleep-related traits. Mendelian randomization identified the causal effects of insomnia on depression, diabetes, and cardiovascular disease, and the protective effects of educational attainment and intracranial volume. Our findings highlight key brain areas and cell types implicated in insomnia, and provide new treatment targets.

Genome-wide association study identifies 30 loci associated with bipolar disorder.

Bipolar disorder is a highly heritable psychiatric disorder. We performed a genome-wide association study (GWAS) including 20,352 cases and 31,358 controls of European descent, with follow-up analysis of 822 variants with P < 1 × 10-4 in an additional 9,412 cases and 137,760 controls. Eight of the 19 variants that were genome-wide significant (P < 5 × 10-8) in the discovery GWAS were not genome-wide significant in the combined analysis, consistent with small effect sizes and limited power but also with genetic heterogeneity. In the combined analysis, 30 loci were genome-wide significant, including 20 newly identified loci. The significant loci contain genes encoding ion channels, neurotransmitter transporters and synaptic components. Pathway analysis revealed nine significantly enriched gene sets, including regulation of insulin secretion and endocannabinoid signaling. Bipolar I disorder is strongly genetically correlated with schizophrenia, driven by psychosis, whereas bipolar II disorder is more strongly correlated with major depressive disorder. These findings address key clinical questions and provide potential biological mechanisms for bipolar disorder.

Exome Chip Meta-analysis Fine Maps Causal Variants and Elucidates the Genetic Architecture of Rare Coding Variants in Smoking and Alcohol Use.

BACKGROUND: Smoking and alcohol use have been associated with common genetic variants in multiple loci. Rare variants within these loci hold promise in the identification of biological mechanisms in substance use. Exome arrays and genotype imputation can now efficiently genotype rare nonsynonymous and loss of function variants. Such variants are expected to have deleterious functional consequences and to contribute to disease risk. METHODS: We analyzed ∼250,000 rare variants from 16 independent studies genotyped with exome arrays and augmented this dataset with imputed data from the UK Biobank. Associations were tested for five phenotypes: cigarettes per day, pack-years, smoking initiation, age of smoking initiation, and alcoholic drinks per week. We conducted stratified heritability analyses, single-variant tests, and gene-based burden tests of nonsynonymous/loss-of-function coding variants. We performed a novel fine-mapping analysis to winnow the number of putative causal variants within associated loci. RESULTS: Meta-analytic sample sizes ranged from 152,348 to 433,216, depending on the phenotype. Rare coding variation explained 1.1% to 2.2% of phenotypic variance, reflecting 11% to 18% of the total single nucleotide polymorphism heritability of these phenotypes. We identified 171 genome-wide associated loci across all phenotypes. Fine mapping identified putative causal variants with double base-pair resolution at 24 of these loci, and between three and 10 variants for 65 loci. Twenty loci contained rare coding variants in the 95% credible intervals. CONCLUSIONS: Rare coding variation significantly contributes to the heritability of smoking and alcohol use. Fine-mapping genome-wide association study loci identifies specific variants contributing to the biological etiology of substance use behavior.

Conditional and interaction gene-set analysis reveals novel functional pathways for blood pressure.

Gene-set analysis provides insight into which functional and biological properties of genes are aetiologically relevant for a particular phenotype. But genes have multiple properties, and these properties are often correlated across genes. This can cause confounding in a gene-set analysis, because one property may be statistically associated even if biologically irrelevant to the phenotype, by being correlated with gene properties that are relevant. To address this issue we present a novel conditional and interaction gene-set analysis approach, which attains considerable functional refinement of its conclusions compared to traditional gene-set analysis. We applied our approach to blood pressure phenotypes in the UK Biobank data (N = 360,243), the results of which we report here. We confirm and further refine several associations with multiple processes involved in heart and blood vessel formation but also identify novel interactions, among others with cardiovascular tissues involved in regulatory pathways of blood pressure homoeostasis.

Meta-analysis of genome-wide association studies for neuroticism in 449,484 individuals identifies novel genetic loci and pathways.

Neuroticism is an important risk factor for psychiatric traits, including depression1, anxiety2,3, and schizophrenia4-6. At the time of analysis, previous genome-wide association studies7-12 (GWAS) reported 16 genomic loci associated to neuroticism10-12. Here we conducted a large GWAS meta-analysis (n = 449,484) of neuroticism and identified 136 independent genome-wide significant loci (124 new at the time of analysis), which implicate 599 genes. Functional follow-up analyses showed enrichment in several brain regions and involvement of specific cell types, including dopaminergic neuroblasts (P = 3.49 × 10-8), medium spiny neurons (P = 4.23 × 10-8), and serotonergic neurons (P = 1.37 × 10-7). Gene set analyses implicated three specific pathways: neurogenesis (P = 4.43 × 10-9), behavioral response to cocaine processes (P = 1.84 × 10-7), and axon part (P = 5.26 × 10-8). We show that neuroticism's genetic signal partly originates in two genetically distinguishable subclusters13 ('depressed affect' and 'worry'), suggesting distinct causal mechanisms for subtypes of individuals. Mendelian randomization analysis showed unidirectional and bidirectional effects between neuroticism and multiple psychiatric traits. These results enhance neurobiological understanding of neuroticism and provide specific leads for functional follow-up experiments.