Recent meta-analyses combining direct genome-wide association studies (GWAS) with those of family history (GWAX) have indicated very low SNP heritability of Alzheimer's disease (AD). These low ...estimates may call into question the prospects of continued progress in genetic discovery for AD within the spectrum of common variants. We highlight dramatic downward biases in previous methods, and we validate a novel method for the estimation of SNP heritability via integration of GWAS and GWAX summary data. We apply our method to investigate the genetic architecture of AD using GWAX from UK Biobank and direct case-control GWAS from the International Genomics of Alzheimer's Project (IGAP). We estimate the liability scale common variant SNP heritability of Clinical AD outside of APOE region at ~7-11%, and we project the corresponding estimate for AD pathology to be up to approximately 23%. We estimate that nearly 90% of common variant SNP heritability of Clinical AD exists outside the APOE region. Rare variants not tagged in standard GWAS may account for additional variance. Our results indicate that, while GWAX for AD in UK Biobank may result in greater attenuation of genetic effects beyond that conventionally assumed, it does not introduce appreciable contamination of signal by genetically distinct traits relative to direct case-control GWAS in IGAP. Genetic risk for AD represents a strong effect of APOE superimposed upon a highly polygenic background.
Estimation of direct and indirect (i.e. parental and/or sibling) genetic effects on phenotypes is becoming increasingly important. We compare several multivariate methods that utilize summary results ...statistics from genome-wide association studies to determine how well they estimate direct and indirect genetic effects. Using data from the UK Biobank, we contrast point estimates and standard errors at individual loci compared to those obtained using individual level data. We show that Genomic structural equation modelling (SEM) outperforms the other methods in accurately estimating conditional genetic effects and their standard errors. We apply Genomic SEM to fertility data in the UK Biobank and partition the genetic effect into female and male fertility and a sibling specific effect. We identify a novel locus for fertility and genetic correlations between fertility and educational attainment, risk taking behaviour, autism and subjective well-being. We recommend Genomic SEM be used to partition genetic effects into direct and indirect components when using summary results from genome-wide association studies.
Gene-environment correlations affect associations between genetic variants and complex traits in genome-wide association studies (GWASs). Here we showed in up to 43,516 British siblings that ...educational attainment polygenic scores capture gene-environment correlations, and that migration extends these gene-environment correlations beyond the family to broader geographic regions. We then ran GWASs on 56 complex traits in up to 254,387 British individuals. Controlling for geographic regions significantly decreased the heritability for socioeconomic status (SES)-related traits, most strongly for educational attainment and income. For most traits, controlling for regions significantly reduced genetic correlations with educational attainment and income, most significantly for body mass index/body fat, sedentary behavior and substance use, consistent with gene-environment correlations related to regional socio-economic differences. The effects of controlling for birthplace and current address suggest both passive and active sources of gene-environment correlations. Our results show that the geographic clustering of DNA and SES introduces gene-environment correlations that affect GWAS results.
Understanding how parents' cognitive and non-cognitive skills influence offspring education is essential for educational, family and economic policy. We use genetics (GWAS-by-subtraction) to assess a ...latent, broad non-cognitive skills dimension. To index parental effects controlling for genetic transmission, we estimate indirect parental genetic effects of polygenic scores on childhood and adulthood educational outcomes, using siblings (N = 47,459), adoptees (N = 6407), and parent-offspring trios (N = 2534) in three UK and Dutch cohorts. We find that parental cognitive and non-cognitive skills affect offspring education through their environment: on average across cohorts and designs, indirect genetic effects explain 36-40% of population polygenic score associations. However, indirect genetic effects are lower for achievement in the Dutch cohort, and for the adoption design. We identify potential causes of higher sibling- and trio-based estimates: prenatal indirect genetic effects, population stratification, and assortative mating. Our phenotype-agnostic, genetically sensitive approach has established overall environmental effects of parents' skills, facilitating future mechanistic work.
Human DNA polymorphisms vary across geographic regions, with the most commonly observed variation reflecting distant ancestry differences. Here we investigate the geographic clustering of common ...genetic variants that influence complex traits in a sample of ~450,000 individuals from Great Britain. Of 33 traits analysed, 21 showed significant geographic clustering at the genetic level after controlling for ancestry, probably reflecting migration driven by socioeconomic status (SES). Alleles associated with educational attainment (EA) showed the most clustering, with EA-decreasing alleles clustering in lower SES areas such as coal mining areas. Individuals who leave coal mining areas carry more EA-increasing alleles on average than those in the rest of Great Britain. The level of geographic clustering is correlated with genetic associations between complex traits and regional measures of SES, health and cultural outcomes. Our results are consistent with the hypothesis that social stratification leaves visible marks in geographic arrangements of common allele frequencies and gene-environment correlations.
Functional genomic methods are needed that consider multiple genetically correlated traits. Here we develop and validate Transcriptome-wide Structural Equation Modeling (T-SEM), a multivariate method ...for studying the effects of tissue-specific gene expression across genetically overlapping traits. T-SEM allows for modeling effects on broad dimensions spanning constellations of traits, while safeguarding against false positives that can arise when effects of gene expression are specific to a subset of traits. We apply T-SEM to investigate the biological mechanisms shared across seven distinct cognitive traits (N = 11,263-331,679), as indexed by a general dimension of genetic sharing (g). We identify 184 genes whose tissue-specific expression is associated with g, including 10 genes not identified in univariate analysis for the individual cognitive traits for any tissue type, and three genes whose expression explained a significant portion of the genetic sharing across g and different subclusters of psychiatric disorders. We go on to apply Stratified Genomic SEM to identify enrichment for g within 28 functional categories. This includes categories indexing the intersection of protein-truncating variant intolerant (PI) genes and specific neuronal cell types, which we also find to be enriched for the genetic covariance between g and a psychotic disorders factor.
We introduce two novel methods for multivariate genome-wide-association meta-analysis (GWAMA) of related traits that correct for sample overlap. A broad range of simulation scenarios supports the ...added value of our multivariate methods relative to univariate GWAMA. We applied the novel methods to life satisfaction, positive affect, neuroticism, and depressive symptoms, collectively referred to as the well-being spectrum (N
= 2,370,390), and found 304 significant independent signals. Our multivariate approaches resulted in a 26% increase in the number of independent signals relative to the four univariate GWAMAs and in an ~57% increase in the predictive power of polygenic risk scores. Supporting transcriptome- and methylome-wide analyses (TWAS and MWAS, respectively) uncovered an additional 17 and 75 independent loci, respectively. Bioinformatic analyses, based on gene expression in brain tissues and cells, showed that genes differentially expressed in the subiculum and GABAergic interneurons are enriched in their effect on the well-being spectrum.
The association between circadian rhythms and diseases has been well established, while the association with mental health is less explored. Given the heritable nature of circadian rhythms, this ...study aimed to investigate the relationship between genes underlying circadian rhythms and mental health outcomes, as well as a possible gene-environment correlation for circadian rhythms. Polygenic scores (PGSs) represent the genetic predisposition to develop a certain trait or disease. In a sample from the Netherlands Twin Register (N = 14,021), PGSs were calculated for two circadian rhythm measures: morningness and relative amplitude (RA). The PGSs were used to predict mental health outcomes such as subjective happiness, quality of life, and depressive symptoms. In addition, we performed the same prediction analysis in a within-family design in a subset of dizygotic twins. The PGS for morningness significantly predicted morningness (R2 = 1.55%) and depressive symptoms (R2 = 0.22%). The PGS for RA significantly predicted general health (R2 = 0.12%) and depressive symptoms (R2 = 0.20%). Item analysis of the depressive symptoms showed that 4 out of 14 items were significantly associated with the PGSs. Overall, the results showed that people with a genetic predisposition of being a morning person or with a high RA are likely to have fewer depressive symptoms. The four associated depressive symptoms described symptoms related to decision-making, energy, and feeling worthless or inferior, rather than sleep. Based on our findings future research should include a substantial role for circadian rhythms in depression research and should further explore the gene-environment correlation in circadian rhythms.
Here we present a method of genome-wide inferred study (GWIS) that provides an approximation of genome-wide association study (GWAS) summary statistics for a variable that is a function of phenotypes ...for which GWAS summary statistics, phenotypic means, and covariances are available. A GWIS can be performed regardless of sample overlap between the GWAS of the phenotypes on which the function depends. Because a GWIS provides association estimates and their standard errors for each SNP, a GWIS can form the basis for polygenic risk scoring, LD score regression, Mendelian randomization studies, biological annotation, and other analyses. GWISs can also be used to boost power of a GWAS meta-analysis where cohorts have not measured all constituent phenotypes in the function. We demonstrate the accuracy of a BMI GWIS by performing power simulations and type I error simulations under varying circumstances, and we apply a GWIS by reconstructing a body mass index (BMI) GWAS based on a weight GWAS and a height GWAS. Furthermore, we apply a GWIS to further our understanding of the underlying genetic structure of bipolar disorder and schizophrenia and their relation to educational attainment. Our analyses suggest that the previously reported genetic correlation between schizophrenia and educational attainment is probably induced by the observed genetic correlation between schizophrenia and bipolar disorder and the previously reported genetic correlation between bipolar disorder and educational attainment.
In recent years, multiple eQTL (expression quantitative trait loci) catalogs have become available that can help understand the functionality of complex trait-related single nucleotide polymorphisms ...(SNPs). In eQTL catalogs, gene expression is often strongly associated with multiple SNPs, which may reflect either one or multiple independent associations. Conditional eQTL analysis allows a distinction between dependent and independent eQTLs. We performed conditional eQTL analysis in 4,896 peripheral blood microarray gene expression samples. Our analysis showed that 35% of genes with a cis eQTL have at least two independent cis eQTLs; for several genes up to 13 independent cis eQTLs were identified. Also, 12% (671) of the independent cis eQTLs identified in conditional analyses were not significant in unconditional analyses. The number of GWAS catalog SNPs identified as eQTL in the conditional analyses increases with 24% as compared to unconditional analyses. We provide an online conditional cis eQTL mapping catalog for whole blood (https://eqtl.onderzoek.io/), which can be used to lookup eQTLs more accurately than in standard unconditional whole blood eQTL databases.