Meta-analysis of genome-wide DNA methylation identifies shared associations across neurodegenerative disorders

BACKGROUND

People with neurodegenerative disorders show diverse clinical syndromes, genetic heterogeneity, and distinct brain pathological changes, but studies report overlap between these features. DNA methylation (DNAm) provides a way to explore this overlap and heterogeneity as it is determined by the combined effects of genetic variation and the environment. In this study, we aim to identify shared blood DNAm differences between controls and people with Alzheimer's disease, amyotrophic lateral sclerosis, and Parkinson's disease.

RESULTS

We use a mixed-linear model method (MOMENT) that accounts for the effect of (un)known confounders, to test for the association of each DNAm site with each disorder. While only three probes are found to be genome-wide significant in each MOMENT association analysis of amyotrophic lateral sclerosis and Parkinson's disease (and none with Alzheimer's disease), a fixed-effects meta-analysis of the three disorders results in 12 genome-wide significant differentially methylated positions. Predicted immune cell-type proportions are disrupted across all neurodegenerative disorders. Protein inflammatory markers are correlated with profile sum-scores derived from disease-associated immune cell-type proportions in a healthy aging cohort. In contrast, they are not correlated with MOMENT DNAm-derived profile sum-scores, calculated using effect sizes of the 12 differentially methylated positions as weights.

CONCLUSIONS

We identify shared differentially methylated positions in whole blood between neurodegenerative disorders that point to shared pathogenic mechanisms. These shared differentially methylated positions may reflect causes or consequences of disease, but they are unlikely to reflect cell-type proportion differences.

Quantifying genetic heterogeneity between continental populations for human height and body mass index

Genome-wide association studies (GWAS) in samples of European ancestry have identified thousands of genetic variants associated with complex traits in humans. However, it remains largely unclear whether these associations can be used in non-European populations. Here, we seek to quantify the proportion of genetic variation for a complex trait shared between continental populations. We estimated the between-population correlation of genetic effects at all SNPs (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r_{g}$$\end{document}rg) or genome-wide significant SNPs (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r_{{g\left( {GWS} \right)}}$$\end{document}rgGWS) for height and body mass index (BMI) in samples of European (EUR; \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n = 49,839$$\end{document}n=49,839) and African (AFR; \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n = 17,426$$\end{document}n=17,426) ancestry. The \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hat{r}_{g}$$\end{document}r^g between EUR and AFR was 0.75 (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text{s}}.{\text{e}}. = 0.035$$\end{document}s.e.=0.035) for height and 0.68 (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text{s}}.{\text{e}}. = 0.062$$\end{document}s.e.=0.062) for BMI, and the corresponding \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hat{r}_{{g\left( {GWS} \right)}}$$\end{document}r^gGWS was 0.82 (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text{s}}.{\text{e}}. = 0.030$$\end{document}s.e.=0.030) for height and 0.87 (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text{s}}.{\text{e}}. = 0.064$$\end{document}s.e.=0.064) for BMI, suggesting that a large proportion of GWAS findings discovered in Europeans are likely applicable to non-Europeans for height and BMI. There was no evidence that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hat{r}_{g}$$\end{document}r^g differs in SNP groups with different levels of between-population difference in allele frequency or linkage disequilibrium, which, however, can be due to the lack of power.

GWAS of peptic ulcer disease implicates Helicobacter pylori infection, other gastrointestinal disorders and depression

Genetic factors are recognized to contribute to peptic ulcer disease (PUD) and other gastrointestinal diseases, such as gastro-oesophageal reflux disease (GORD), irritable bowel syndrome (IBS) and inflammatory bowel disease (IBD). Here, genome-wide association study (GWAS) analyses based on 456,327 UK Biobank (UKB) individuals identify 8 independent and significant loci for PUD at, or near, genes MUC1, MUC6, FUT2, PSCA, ABO, CDX2, GAST and CCKBR. There are previously established roles in susceptibility to Helicobacter pylori infection, response to counteract infection-related damage, gastric acid secretion or gastrointestinal motility for these genes. Only two associations have been previously reported for duodenal ulcer, here replicated trans-ancestrally. The results highlight the role of host genetic susceptibility to infection. Post-GWAS analyses for PUD, GORD, IBS and IBD add insights into relationships between these gastrointestinal diseases and their relationships with depression, a commonly comorbid disorder.

Creating and validating a DNA methylation-based proxy for interleukin-6

Background

Studies evaluating the relationship between chronic inflammation and cognitive functioning have produced heterogeneous results. A potential reason for this is the variability of inflammatory mediators which could lead to misclassifications of individuals’ persisting levels of inflammation. DNA methylation (DNAm) has shown utility in indexing environmental exposures and could be leveraged to provide proxy signatures of chronic inflammation.

Method

We conducted an elastic net regression of interleukin-6 (IL-6) in a cohort of 875 older adults (Lothian Birth Cohort 1936; mean age: 70 years) to develop a DNAm-based predictor. The predictor was tested in an independent cohort (Generation Scotland; N = 7028 [417 with measured IL-6], mean age: 51 years).

Results

A weighted score from 35 CpG sites optimally predicted IL-6 in the independent test set (Generation Scotland; R² = 4.4%, p = 2.1 × 10⁻⁵). In the independent test cohort, both measured IL-6 and the DNAm proxy increased with age (serum IL-6: n = 417, β = 0.02, SE = 0.004, p = 1.3 × 10⁻⁷; DNAm IL-6 score: N = 7028, β = 0.02, SE = 0.0009, p < 2 × 10⁻¹⁶). Serum IL-6 did not associate with cognitive ability (n = 417, β = −0.06, SE = 0.05, p = .19); however, an inverse association was identified between the DNAm score and cognitive functioning (N = 7028, β = −0.16, SE = 0.02, p_FDR < 2 × 10⁻¹⁶).

Conclusions

These results suggest methylation-based predictors can be used as proxies for inflammatory markers, potentially allowing for further insight into the relationship between inflammation and pertinent health outcomes.

Phenotypic covariance across the entire spectrum of relatedness for 86 billion pairs of individuals

Attributing the similarity between individuals to genetic and non-genetic factors is central to genetic analyses. In this paper we use the genomic relationship ([Formula: see text]) among 417,060 individuals to investigate the phenotypic covariance between pairs of individuals for 32 traits across the spectrum of relatedness, from unrelated pairs through to identical twins. We find linear relationships between phenotypic covariance and [Formula: see text] that agree with the SNP-based heritability ([Formula: see text]) in unrelated pairs ([Formula: see text]), and with pedigree-estimated heritability in close relatives ([Formula: see text]). The covariance increases faster than [Formula: see text] in distant relatives ([Formula: see text]), and we attribute this to imperfect linkage disequilibrium between causal variants and the common variants used to construct [Formula: see text]. We also examine the effect of assortative mating on heritability estimates from different experimental designs. We find that full-sib identity-by-descent regression estimates for height (0.66 s.e. 0.07) are consistent with estimates from close relatives (0.82 s.e. 0.04) after accounting for the effect of assortative mating.

Analysis of common genetic variation and rare CNVs in the Australian Autism Biobank

BACKGROUND

Autism spectrum disorder (ASD) is a complex neurodevelopmental condition whose biological basis is yet to be elucidated. The Australian Autism Biobank (AAB) is an initiative of the Cooperative Research Centre for Living with Autism (Autism CRC) to establish an Australian resource of biospecimens, phenotypes and genomic data for research on autism.

METHODS

Genome-wide single-nucleotide polymorphism genotypes were available for 2,477 individuals (after quality control) from 546 families (436 complete), including 886 participants aged 2 to 17 years with diagnosed (n = 871) or suspected (n = 15) ASD, 218 siblings without ASD, 1,256 parents, and 117 unrelated children without an ASD diagnosis. The genetic data were used to confirm familial relationships and assign ancestry, which was majority European (n = 1,964 European individuals). We generated polygenic scores (PGS) for ASD, IQ, chronotype and height in the subset of Europeans, and in 3,490 unrelated ancestry-matched participants from the UK Biobank. We tested for group differences for each PGS, and performed prediction analyses for related phenotypes in the AAB. We called copy-number variants (CNVs) in all participants, and intersected these with high-confidence ASD- and intellectual disability (ID)-associated CNVs and genes from the public domain.

RESULTS

The ASD (p = 6.1e-13), sibling (p = 4.9e-3) and unrelated (p = 3.0e-3) groups had significantly higher ASD PGS than UK Biobank controls, whereas this was not the case for height-a control trait. The IQ PGS was a significant predictor of measured IQ in undiagnosed children (r = 0.24, p = 2.1e-3) and parents (r = 0.17, p = 8.0e-7; 4.0% of variance), but not the ASD group. Chronotype PGS predicted sleep disturbances within the ASD group (r = 0.13, p = 1.9e-3; 1.3% of variance). In the CNV analysis, we identified 13 individuals with CNVs overlapping ASD/ID-associated CNVs, and 12 with CNVs overlapping ASD/ID/developmental delay-associated genes identified on the basis of de novo variants.

LIMITATIONS

This dataset is modest in size, and the publicly-available genome-wide-association-study (GWAS) summary statistics used to calculate PGS for ASD and other traits are relatively underpowered.

CONCLUSIONS

We report on common genetic variation and rare CNVs within the AAB. Prediction analyses using currently available GWAS summary statistics are largely consistent with expected relationships based on published studies. As the size of publicly-available GWAS summary statistics grows, the phenotypic depth of the AAB dataset will provide many opportunities for analyses of autism profiles and co-occurring conditions, including when integrated with other omics datasets generated from AAB biospecimens (blood, urine, stool, hair).

Publisher Correction: Genome-wide analyses of behavioural traits are subject to bias by misreports and longitudinal changes

A Correction to this paper has been published: https://doi.org/10.1038/s41467-021-21294-1.

Multi-omic and multi-species meta-analyses of nicotine consumption

Cross-species translational approaches to human genomic analyses are lacking. The present study uses an integrative framework to investigate how genes associated with nicotine use in model organisms contribute to the genetic architecture of human tobacco consumption. First, we created a model organism geneset by collecting results from five animal models of nicotine exposure (RNA expression changes in brain) and then tested the relevance of these genes and flanking genetic variation using genetic data from human cigarettes per day (UK BioBank N = 123,844; all European Ancestry). We tested three hypotheses: (1) DNA variation in, or around, the 'model organism geneset' will contribute to the heritability to human tobacco consumption, (2) that the model organism genes will be enriched for genes associated with human tobacco consumption, and (3) that a polygenic score based off our model organism geneset will predict tobacco consumption in the AddHealth sample (N = 1667; all European Ancestry). Our results suggested that: (1) model organism genes accounted for ~5-36% of the observed SNP-heritability in human tobacco consumption (enrichment: 1.60-31.45), (2) model organism genes, but not negative control genes, were enriched for the gene-based associations (MAGMA, H-MAGMA, SMultiXcan) for human cigarettes per day, and (3) polygenic scores based on our model organism geneset predicted cigarettes per day in an independent sample. Altogether, these findings highlight the advantages of using multiple species evidence to isolate genetic factors to better understand the etiological complexity of tobacco and other nicotine consumption.

Genome-wide analyses of behavioural traits are subject to bias by misreports and longitudinal changes

Genome-wide association studies (GWAS) have discovered numerous genetic variants associated with human behavioural traits. However, behavioural traits are subject to misreports and longitudinal changes (MLC) which can cause biases in GWAS and follow-up analyses. Here, we demonstrate that individuals with higher disease burden in the UK Biobank (n = 455,607) are more likely to misreport or reduce their alcohol consumption levels, and propose a correction procedure to mitigate the MLC-induced biases. The alcohol consumption GWAS signals removed by the MLC corrections are enriched in metabolic/cardiovascular traits. Almost all the previously reported negative estimates of genetic correlations between alcohol consumption and common diseases become positive/non-significant after the MLC corrections. We also observe MLC biases for smoking and physical activities in the UK Biobank. Our findings provide a plausible explanation of the controversy about the effects of alcohol consumption on health outcomes and a caution for future analyses of self-reported behavioural traits in biobank data.

Variation in VKORC1 Is Associated with Vascular Dementia

BACKGROUND

The genetic variant rs9923231 (VKORC1) is associated with differences in the coagulation of blood and consequentially with sensitivity to the drug warfarin. Variation in VKORC1 has been linked in a gene-based test to dementia/Alzheimer's disease in the parents of participants, with suggestive evidence for an association for rs9923231 (p = 1.8×10-7), which was included in the genome-wide significant KAT8 locus.

OBJECTIVE

Our study aimed to investigate whether the relationship between rs9923231 and dementia persists only for certain dementia sub-types, and if those taking warfarin are at greater risk.

METHODS

We used logistic regression and data from 238,195 participants from UK Biobank to examine the relationship between VKORC1, risk of dementia, and the interplay with warfarin use.

RESULTS

Parental history of dementia, APOE variant, atrial fibrillation, diabetes, hypertension, and hypercholesterolemia all had strong associations with vascular dementia (p < 4.6×10-6). The T-allele in rs9923231 was linked to a lower warfarin dose (βperT - allele = -0.29, p < 2×10-16) and risk of vascular dementia (OR = 1.17, p = 0.010), but not other dementia sub-types. However, the risk of vascular dementia was not affected by warfarin use in carriers of the T-allele.

CONCLUSION

Our study reports for the first time an association between rs9923231 and vascular dementia, but further research is warranted to explore potential mechanisms and specify the relationship between rs9923231 and features of vascular dementia.

From Basic Science to Clinical Application of Polygenic Risk Scores: A Primer

IMPORTANCE

Polygenic risk scores (PRS) are predictors of the genetic susceptibilities of individuals to diseases. All individuals have DNA risk variants for all common diseases, but genetic susceptibility differences between people reflect the cumulative burden of these. Polygenic risk scores for an individual are calculated as weighted counts of thousands of risk variants that they carry, where the risk variants and their weights have been identified in genome-wide association studies. Here, we review the underlying basic science of PRS, providing a foundation for understanding the potential clinical utility and limitations of PRS.

OBSERVATIONS

Polygenic risk scores can be calculated for a wide range of diseases from a saliva or blood sample using genotyping technologies that are inexpensive. While genotyping only needs to be done once for each individual in their lifetime, the PRS can be recalculated as identification of risk variants improves. On their own, PRS will never be able to establish or definitively predict future diagnoses of common complex conditions because genetic factors only contribute part of the risk, and PRS will only ever capture part of the genetic contributions. Nonetheless, just as clinical medicine uses a multitude of other predictive measures, PRS either on their own or as part of multivariable predictive algorithms could play a role.

CONCLUSIONS AND RELEVANCE

Utility of PRS in clinical medicine and ethical issues related to their use should be evaluated in the context of realistic expectations of what PRS can and cannot deliver. For different diseases, PRS could have utility in community settings (stratification to better triage people into established screening programs) or could contribute to clinical decision-making for those presenting with symptoms but where formal diagnosis is unclear. In principle, PRS could contribute to treatment choices, but more data are needed to allow development of PRS in this context.

Evolutionary Genomics at the Human-Environment Interface in Africa

We report on the first meeting of SMBE in Africa. SMBE Malawi was initiated to bring together African and international researchers who use genetics or genomics to study natural systems impacted by human activities. The goals of this conference were 1) to reach a world-class standard of science with a large number of contributions from Africa, 2) to initiate exchange between African and international researchers, and 3) to identify challenges and opportunities for evolutionary genomics research in Africa. As repored, we think that we have achieved these goals and make suggestions on the way forward for African evolutionary genomics research.

A unified framework for association and prediction from vertex-wise grey-matter structure

The recent availability of large‐scale neuroimaging cohorts facilitates deeper characterisation of the relationship between phenotypic and brain architecture variation in humans. Here, we investigate the association (previously coined morphometricity) of a phenotype with all 652,283 vertex‐wise measures of cortical and subcortical morphology in a large data set from the UK Biobank (UKB; N = 9,497 for discovery, N = 4,323 for replication) and the Human Connectome Project (N = 1,110). We used a linear mixed model with the brain measures of individuals fitted as random effects with covariance relationships estimated from the imaging data. We tested 167 behavioural, cognitive, psychiatric or lifestyle phenotypes and found significant morphometricity for 58 phenotypes (spanning substance use, blood assay results, education or income level, diet, depression, and cognition domains), 23 of which replicated in the UKB replication set or the HCP. We then extended the model for a bivariate analysis to estimate grey‐matter correlation between phenotypes, which revealed that body size (i.e., height, weight, BMI, waist and hip circumference, body fat percentage) could account for a substantial proportion of the morphometricity (confirmed using a conditional analysis), providing possible insight into previous MRI case–control results for psychiatric disorders where case status is associated with body mass index. Our LMM framework also allowed to predict some of the associated phenotypes from the vertex‐wise measures, in two independent samples. Finally, we demonstrated additional new applications of our approach (a) region of interest (ROI) analysis that retain the vertex‐wise complexity; (b) comparison of the information retained by different MRI processings.

Risk prediction of late-onset Alzheimer's disease implies an oligogenic architecture

Genetic association studies have identified 44 common genome-wide significant risk loci for late-onset Alzheimer's disease (LOAD). However, LOAD genetic architecture and prediction are unclear. Here we estimate the optimal P-threshold (Poptimal) of a genetic risk score (GRS) for prediction of LOAD in three independent datasets comprising 676 cases and 35,675 family history proxy cases. We show that the discriminative ability of GRS in LOAD prediction is maximised when selecting a small number of SNPs. Both simulation results and direct estimation indicate that the number of causal common SNPs for LOAD may be less than 100, suggesting LOAD is more oligogenic than polygenic. The best GRS explains approximately 75% of SNP-heritability, and individuals in the top decile of GRS have ten-fold increased odds when compared to those in the bottom decile. In addition, 14 variants are identified that contribute to both LOAD risk and age at onset of LOAD.

Using prior information from humans to prioritize genes and gene-associated variants for complex traits in livestock

Genome-Wide Association Studies (GWAS) in large human cohorts have identified thousands of loci associated with complex traits and diseases. For identifying the genes and gene-associated variants that underlie complex traits in livestock, especially where sample sizes are limiting, it may help to integrate the results of GWAS for equivalent traits in humans as prior information. In this study, we sought to investigate the usefulness of results from a GWAS on human height as prior information for identifying the genes and gene-associated variants that affect stature in cattle, using GWAS summary data on samples sizes of 700,000 and 58,265 for humans and cattle, respectively. Using Fisher's exact test, we observed a significant proportion of cattle stature-associated genes (30/77) that are also associated with human height (odds ratio = 5.1, p = 3.1e-10). Result of randomized sampling tests showed that cattle orthologs of human height-associated genes, hereafter referred to as candidate genes (C-genes), were more enriched for cattle stature GWAS signals than random samples of genes in the cattle genome (p = 0.01). Randomly sampled SNPs within the C-genes also tend to explain more genetic variance for cattle stature (up to 13.2%) than randomly sampled SNPs within random cattle genes (p = 0.09). The most significant SNPs from a cattle GWAS for stature within the C-genes did not explain more genetic variance for cattle stature than the most significant SNPs within random cattle genes (p = 0.87). Altogether, our findings support previous studies that suggest a similarity in the genetic regulation of height across mammalian species. However, with the availability of a powerful GWAS for stature that combined data from 8 cattle breeds, prior information from human-height GWAS does not seem to provide any additional benefit with respect to the identification of genes and gene-associated variants that affect stature in cattle.

Theoretical and empirical quantification of the accuracy of polygenic scores in ancestry divergent populations

Polygenic scores (PGS) have been widely used to predict disease risk using variants identified from genome-wide association studies (GWAS). To date, most GWAS have been conducted in populations of European ancestry, which limits the use of GWAS-derived PGS in non-European ancestry populations. Here, we derive a theoretical model of the relative accuracy (RA) of PGS across ancestries. We show through extensive simulations that the RA of PGS based on genome-wide significant SNPs can be predicted accurately from modelling linkage disequilibrium (LD), minor allele frequencies (MAF), cross-population correlations of causal SNP effects and heritability. We find that LD and MAF differences between ancestries can explain between 70 and 80% of the loss of RA of European-based PGS in African ancestry for traits like body mass index and type 2 diabetes. Our results suggest that causal variants underlying common genetic variation identified in European ancestry GWAS are mostly shared across continents.

Blood DNA methylation sites predict death risk in a longitudinal study of 12, 300 individuals

DNA methylation has fundamental roles in gene programming and aging that may help predict mortality. However, no large-scale study has investigated whether site-specific DNA methylation predicts all-cause mortality. We used the Illumina-HumanMethylation450-BeadChip to identify blood DNA methylation sites associated with all-cause mortality for 12, 300 participants in 12 Cohorts of the Heart and Aging Research in Genetic Epidemiology (CHARGE) Consortium. Over an average 10-year follow-up, there were 2,561 deaths across the cohorts. Nine sites mapping to three intergenic and six gene-specific regions were associated with mortality (P < 9.3x10^-7) independently of age and other mortality predictors. Six sites (cg14866069, cg23666362, cg20045320, cg07839457, cg07677157, cg09615688)—mapping respectively to BMPR1B, MIR1973, IFITM3, NLRC5, and two intergenic regions—were associated with reduced mortality risk. The remaining three sites (cg17086398, cg12619262, cg18424841)—mapping respectively to SERINC2, CHST12, and an intergenic region—were associated with increased mortality risk. DNA methylation at each site predicted 5%–15% of all deaths. We also assessed the causal association of those sites to age-related chronic diseases by using Mendelian randomization, identifying weak causal relationship between cg18424841 and cg09615688 with coronary heart disease. Of the nine sites, three (cg20045320, cg07839457, cg07677157) were associated with lower incidence of heart disease risk and two (cg20045320, cg07839457) with smoking and inflammation in prior CHARGE analyses. Methylation of cg20045320, cg07839457, and cg17086398 was associated with decreased expression of nearby genes (IFITM3, IRF, NLRC5, MT1, MT2, MARCKSL1) linked to immune responses and cardiometabolic diseases. These sites may serve as useful clinical tools for mortality risk assessment and preventative care.

Multi-method genome- and epigenome-wide studies of inflammatory protein levels in healthy older adults

BACKGROUND

The molecular factors which control circulating levels of inflammatory proteins are not well understood. Furthermore, association studies between molecular probes and human traits are often performed by linear model-based methods which may fail to account for complex structure and interrelationships within molecular datasets.

METHODS

In this study, we perform genome- and epigenome-wide association studies (GWAS/EWAS) on the levels of 70 plasma-derived inflammatory protein biomarkers in healthy older adults (Lothian Birth Cohort 1936; n = 876; Olink® inflammation panel). We employ a Bayesian framework (BayesR+) which can account for issues pertaining to data structure and unknown confounding variables (with sensitivity analyses using ordinary least squares- (OLS) and mixed model-based approaches).

RESULTS

We identified 13 SNPs associated with 13 proteins (n = 1 SNP each) concordant across OLS and Bayesian methods. We identified 3 CpG sites spread across 3 proteins (n = 1 CpG each) that were concordant across OLS, mixed-model and Bayesian analyses. Tagged genetic variants accounted for up to 45% of variance in protein levels (for MCP2, 36% of variance alone attributable to 1 polymorphism). Methylation data accounted for up to 46% of variation in protein levels (for CXCL10). Up to 66% of variation in protein levels (for VEGFA) was explained using genetic and epigenetic data combined. We demonstrated putative causal relationships between CD6 and IL18R1 with inflammatory bowel disease and between IL12B and Crohn's disease.

CONCLUSIONS

Our data may aid understanding of the molecular regulation of the circulating inflammatory proteome as well as causal relationships between inflammatory mediators and disease.

Genome-wide association study identifies 143 loci associated with 25 hydroxyvitamin D concentration

Vitamin D deficiency is a candidate risk factor for a range of adverse health outcomes. In a genome-wide association study of 25 hydroxyvitamin D (25OHD) concentration in 417,580 Europeans we identify 143 independent loci in 112 1-Mb regions, providing insights into the physiology of vitamin D and implicating genes involved in lipid and lipoprotein metabolism, dermal tissue properties, and the sulphonation and glucuronidation of 25OHD. Mendelian randomization models find no robust evidence that 25OHD concentration has causal effects on candidate phenotypes (e.g. BMI, psychiatric disorders), but many phenotypes have (direct or indirect) causal effects on 25OHD concentration, clarifying the epidemiological relationship between 25OHD status and the health outcomes examined in this study.

Tissue-specific sex differences in human gene expression

Despite extensive sex differences in human complex traits and disease, the male and female genomes differ only in the sex chromosomes. This implies that most sex-differentiated traits are the result of differences in the expression of genes that are common to both sexes. While sex differences in gene expression have been observed in a range of different tissues, the biological mechanisms for tissue-specific sex differences (TSSDs) in gene expression are not well understood. A total of 30 640 autosomal and 1021 X-linked transcripts were tested for heterogeneity in sex difference effect sizes in n = 617 individuals across 40 tissue types in Genotype-Tissue Expression (GTEx). This identified 65 autosomal and 66 X-linked TSSD transcripts (corresponding to unique genes) at a stringent significance threshold. Results for X-linked TSSD transcripts showed mainly concordant direction of sex differences across tissues and replicate previous findings. Autosomal TSSD transcripts had mainly discordant direction of sex differences across tissues. The top cis-expression quantitative trait loci (eQTLs) across tissues for autosomal TSSD transcripts are located a similar distance away from the nearest androgen and estrogen binding motifs and the nearest enhancer, as compared to cis-eQTLs for transcripts with stable sex differences in gene expression across tissue types. Enhancer regions that overlap top cis-eQTLs for TSSD transcripts, however, were found to be more dispersed across tissues. These observations suggest that androgen and estrogen regulatory elements in a cis region may play a common role in sex differences in gene expression, but TSSD in gene expression may additionally be due to causal variants located in tissue-specific enhancer regions.

Analysis of DNA methylation associates the cystine-glutamate antiporter SLC7A11 with risk of Parkinson's disease

An improved understanding of etiological mechanisms in Parkinson's disease (PD) is urgently needed because the number of affected individuals is projected to increase rapidly as populations age. We present results from a blood-based methylome-wide association study of PD involving meta-analysis of 229 K CpG probes in 1,132 cases and 999 controls from two independent cohorts. We identify two previously unreported epigenome-wide significant associations with PD, including cg06690548 on chromosome 4. We demonstrate that cg06690548 hypermethylation in PD is associated with down-regulation of the SLC7A11 gene and show this is consistent with an environmental exposure, as opposed to medications or genetic factors with effects on DNA methylation or gene expression. These findings are notable because SLC7A11 codes for a cysteine-glutamate anti-porter regulating levels of the antioxidant glutathione, and it is a known target of the environmental neurotoxin β-methylamino-L-alanine (BMAA). Our study identifies the SLC7A11 gene as a plausible biological target in PD.

Significant out-of-sample classification from methylation profile scoring for amyotrophic lateral sclerosis

We conducted DNA methylation association analyses using Illumina 450K data from whole blood for an Australian amyotrophic lateral sclerosis (ALS) case–control cohort (782 cases and 613 controls). Analyses used mixed linear models as implemented in the OSCA software. We found a significantly higher proportion of neutrophils in cases compared to controls which replicated in an independent cohort from the Netherlands (1159 cases and 637 controls). The OSCA MOMENT linear mixed model has been shown in simulations to best account for confounders. When combined in a methylation profile score, the 25 most-associated probes identified by MOMENT significantly classified case–control status in the Netherlands sample (area under the curve, AUC = 0.65, CI_95% = [0.62–0.68], p = 8.3 × 10⁻²²). The maximum AUC achieved was 0.69 (CI_95% = [0.66–0.71], p = 4.3 × 10⁻³⁴) when cell-type proportion was included in the predictor.

Exploring the role of genetic confounding in the association between maternal and offspring body mass index: evidence from three birth cohorts

BACKGROUND

Maternal pre-pregnancy body mass index (BMI) is positively associated with offspring birth weight (BW) and BMI in childhood and adulthood. Each of these associations could be due to causal intrauterine effects, or confounding (genetic or environmental), or some combination of these. Here we estimate the extent to which the association between maternal BMI and offspring body size is explained by offspring genotype, as a first step towards establishing the importance of genetic confounding.

METHODS

We examined the associations of maternal pre-pregnancy BMI with offspring BW and BMI at 1, 5, 10 and 15 years, in three European birth cohorts (n ≤11 498). Bivariate Genomic-relatedness-based Restricted Maximum Likelihood implemented in the GCTA software (GCTA-GREML) was used to estimate the extent to which phenotypic covariance was explained by offspring genotype as captured by common imputed single nucleotide polymorphisms (SNPs). We merged individual participant data from all cohorts, enabling calculation of pooled estimates.

RESULTS

Phenotypic covariance (equivalent here to Pearson's correlation coefficient) between maternal BMI and offspring phenotype was 0.15 [95% confidence interval (CI): 0.13, 0.17] for offspring BW, increasing to 0.29 (95% CI: 0.26, 0.31) for offspring 15 year BMI. Covariance explained by offspring genotype was negligible for BW [-0.04 (95% CI: -0.09, 0.01)], but increased to 0.12 (95% CI: 0.04, 0.21) at 15 years, which is equivalent to 43% (95% CI: 15%, 72%) of the phenotypic covariance. Sensitivity analyses using weight, BMI and ponderal index as the offspring phenotype at all ages showed similar results.

CONCLUSIONS

Offspring genotype explains a substantial fraction of the covariance between maternal BMI and offspring adolescent BMI. This is consistent with a potentially important role for genetic confounding as a driver of the maternal BMI-offspring BMI association.

Genome-wide association and Mendelian randomisation analysis provide insights into the pathogenesis of heart failure

Heart failure (HF) is a leading cause of morbidity and mortality worldwide. A small proportion of HF cases are attributable to monogenic cardiomyopathies and existing genome-wide association studies (GWAS) have yielded only limited insights, leaving the observed heritability of HF largely unexplained. We report results from a GWAS meta-analysis of HF comprising 47,309 cases and 930,014 controls. Twelve independent variants at 11 genomic loci are associated with HF, all of which demonstrate one or more associations with coronary artery disease (CAD), atrial fibrillation, or reduced left ventricular function, suggesting shared genetic aetiology. Functional analysis of non-CAD-associated loci implicate genes involved in cardiac development (MYOZ1, SYNPO2L), protein homoeostasis (BAG3), and cellular senescence (CDKN1A). Mendelian randomisation analysis supports causal roles for several HF risk factors, and demonstrates CAD-independent effects for atrial fibrillation, body mass index, and hypertension. These findings extend our knowledge of the pathways underlying HF and may inform new therapeutic strategies.

An epigenome-wide association study of sex-specific chronological ageing

BACKGROUND

Advanced age is associated with cognitive and physical decline and is a major risk factor for a multitude of disorders. There is also a gap in life expectancy between males and females. DNA methylation differences have been shown to be associated with both age and sex. Here, we investigate age-by-sex differences in blood-based DNA methylation in an unrelated cohort of 2586 individuals between the ages of 18 and 87 years, with replication in a further 4450 individuals between the ages of 18 and 93 years.

METHODS

Linear regression models were applied, with stringent genome-wide significance thresholds (p < 3.6 × 10-8) used in both the discovery and replication data. A second, highly conservative mixed linear model method that better controls the false-positive rate was also applied, using the same genome-wide significance thresholds.

RESULTS

Using the linear regression method, 52 autosomal and 597 X-linked CpG sites, mapping to 251 unique genes, replicated with concordant effect size directions in the age-by-sex interaction analysis. The site with the greatest difference mapped to GAGE10, an X-linked gene. Here, DNA methylation levels remained stable across the male adult age range (DNA methylation by age r = 0.02) but decreased across female adult age range (DNA methylation by age r = - 0.61). One site (cg23722529) with a significant age-by-sex interaction also had a quantitative trait locus (rs17321482) that is a genome-wide significant variant for prostate cancer. The mixed linear model method identified 11 CpG sites associated with the age-by-sex interaction.

CONCLUSION

The majority of differences in age-associated DNA methylation trajectories between sexes are present on the X chromosome. Several of these differences occur within genes that have been implicated in sexually dimorphic traits.

Identification of novel risk loci, causal insights, and heritable risk for Parkinson's disease: a meta-analysis of genome-wide association studies

BACKGROUND

Genome-wide association studies (GWAS) in Parkinson's disease have increased the scope of biological knowledge about the disease over the past decade. We aimed to use the largest aggregate of GWAS data to identify novel risk loci and gain further insight into the causes of Parkinson's disease.

METHODS

We did a meta-analysis of 17 datasets from Parkinson's disease GWAS available from European ancestry samples to nominate novel loci for disease risk. These datasets incorporated all available data. We then used these data to estimate heritable risk and develop predictive models of this heritability. We also used large gene expression and methylation resources to examine possible functional consequences as well as tissue, cell type, and biological pathway enrichments for the identified risk factors. Additionally, we examined shared genetic risk between Parkinson's disease and other phenotypes of interest via genetic correlations followed by Mendelian randomisation.

FINDINGS

Between Oct 1, 2017, and Aug 9, 2018, we analysed 7·8 million single nucleotide polymorphisms in 37 688 cases, 18 618 UK Biobank proxy-cases (ie, individuals who do not have Parkinson's disease but have a first degree relative that does), and 1·4 million controls. We identified 90 independent genome-wide significant risk signals across 78 genomic regions, including 38 novel independent risk signals in 37 loci. These 90 variants explained 16-36% of the heritable risk of Parkinson's disease depending on prevalence. Integrating methylation and expression data within a Mendelian randomisation framework identified putatively associated genes at 70 risk signals underlying GWAS loci for follow-up functional studies. Tissue-specific expression enrichment analyses suggested Parkinson's disease loci were heavily brain-enriched, with specific neuronal cell types being implicated from single cell data. We found significant genetic correlations with brain volumes (false discovery rate-adjusted p=0·0035 for intracranial volume, p=0·024 for putamen volume), smoking status (p=0·024), and educational attainment (p=0·038). Mendelian randomisation between cognitive performance and Parkinson's disease risk showed a robust association (p=8·00 × 10^-7).

INTERPRETATION

These data provide the most comprehensive survey of genetic risk within Parkinson's disease to date, to the best of our knowledge, by revealing many additional Parkinson's disease risk loci, providing a biological context for these risk factors, and showing that a considerable genetic component of this disease remains unidentified. These associations derived from European ancestry datasets will need to be followed-up with more diverse data.

FUNDING

The National Institute on Aging at the National Institutes of Health (USA), The Michael J Fox Foundation, and The Parkinson's Foundation (see appendix for full list of funding sources).

No Evidence for Social Genetic Effects or Genetic Similarity Among Friends Beyond that Due to Population Stratification: A Reappraisal of Domingue et al (2018)

Using data from 5500 adolescents from the National Longitudinal Study of Adolescent to Adult Health, Domingue et al. (Proc Natl Acad Sci 25:256., 2018) claimed to show that friends are genetically more similar to one another than randomly selected peers, beyond the confounding effects of population stratification by ancestry. The authors also claimed to show 'social-genetic' effects, whereby individuals' educational attainment (EA) is influenced by their friends' genes. We argue that neither claim is justified by the data. Mathematically we show that (1) the genetic similarity reported between friends is far larger than theoretically possible if it was caused by phenotypic assortment as the authors claim; uncontrolled population stratification is a likely reason for the genetic similarity they observed, and (2) significant association between individuals' EA and their friends' polygenic scores for EA is a necessary consequence of EA similarity among friends, and does not provide evidence for social-genetic effects. Going forward, we urge caution in the analysis and interpretation of data at the intersection of human genetics and the social sciences.

Genetic correlates of social stratification in Great Britain

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.

Improved precision of epigenetic clock estimates across tissues and its implication for biological ageing

BACKGROUND

DNA methylation changes with age. Chronological age predictors built from DNA methylation are termed 'epigenetic clocks'. The deviation of predicted age from the actual age ('age acceleration residual', AAR) has been reported to be associated with death. However, it is currently unclear how a better prediction of chronological age affects such association.

METHODS

In this study, we build multiple predictors based on training DNA methylation samples selected from 13,661 samples (13,402 from blood and 259 from saliva). We use the Lothian Birth Cohorts of 1921 (LBC1921) and 1936 (LBC1936) to examine whether the association between AAR (from these predictors) and death is affected by (1) improving prediction accuracy of an age predictor as its training sample size increases (from 335 to 12,710) and (2) additionally correcting for confounders (i.e., cellular compositions). In addition, we investigated the performance of our predictor in non-blood tissues.

RESULTS

We found that in principle, a near-perfect age predictor could be developed when the training sample size is sufficiently large. The association between AAR and mortality attenuates as prediction accuracy increases. AAR from our best predictor (based on Elastic Net, https://github.com/qzhang314/DNAm-based-age-predictor ) exhibits no association with mortality in both LBC1921 (hazard ratio = 1.08, 95% CI 0.91-1.27) and LBC1936 (hazard ratio = 1.00, 95% CI 0.79-1.28). Predictors based on small sample size are prone to confounding by cellular compositions relative to those from large sample size. We observed comparable performance of our predictor in non-blood tissues with a multi-tissue-based predictor.

CONCLUSIONS

This study indicates that the epigenetic clock can be improved by increasing the training sample size and that its association with mortality attenuates with increased prediction of chronological age.

Genotype-by-environment interactions inferred from genetic effects on phenotypic variability in the UK Biobank

Genotype-by-environment interaction (GEI) is a fundamental component in understanding complex trait variation. However, it remains challenging to identify genetic variants with GEI effects in humans largely because of the small effect sizes and the difficulty of monitoring environmental fluctuations. Here, we demonstrate that GEI can be inferred from genetic variants associated with phenotypic variability in a large sample without the need of measuring environmental factors. We performed a genome-wide variance quantitative trait locus (vQTL) analysis of ~5.6 million variants on 348,501 unrelated individuals of European ancestry for 13 quantitative traits in the UK Biobank and identified 75 significant vQTLs with P < 2.0 × 10-9 for 9 traits, especially for those related to obesity. Direct GEI analysis with five environmental factors showed that the vQTLs were strongly enriched with GEI effects. Our results indicate pervasive GEI effects for obesity-related traits and demonstrate the detection of GEI without environmental data.

Genome and epigenome wide studies of neurological protein biomarkers in the Lothian Birth Cohort 1936

Although plasma proteins may serve as markers of neurological disease risk, the molecular mechanisms responsible for inter-individual variation in plasma protein levels are poorly understood. Therefore, we conduct genome- and epigenome-wide association studies on the levels of 92 neurological proteins to identify genetic and epigenetic loci associated with their plasma concentrations (n = 750 healthy older adults). We identify 41 independent genome-wide significant (P < 5.4 × 10^-10) loci for 33 proteins and 26 epigenome-wide significant (P < 3.9 × 10^-10) sites associated with the levels of 9 proteins. Using this information, we identify biological pathways in which putative neurological biomarkers are implicated (neurological, immunological and extracellular matrix metabolic pathways). We also observe causal relationships (by Mendelian randomisation analysis) between changes in gene expression (DRAXIN, MDGA1 and KYNU), or DNA methylation profiles (MATN3, MDGA1 and NEP), and altered plasma protein levels. Together, this may help inform causal relationships between biomarkers and neurological diseases.

The effect of X-linked dosage compensation on complex trait variation

Quantitative genetics theory predicts that X-chromosome dosage compensation (DC) will have a detectable effect on the amount of genetic and therefore phenotypic trait variances at associated loci in males and females. Here, we systematically examine the role of DC in humans in 20 complex traits in a sample of more than 450,000 individuals from the UK Biobank and 1600 gene expression traits from a sample of 2000 individuals as well as across-tissue gene expression from the GTEx resource. We find approximately twice as much X-linked genetic variation across the UK Biobank traits in males (mean h²_SNP = 0.63%) compared to females (mean h²_SNP = 0.30%), confirming the predicted DC effect. Our DC estimates for complex traits and gene expression are consistent with a small proportion of genes escaping X-inactivation in a trait- and tissue-dependent manner. Finally, we highlight examples of biologically relevant X-linked heterogeneity between the sexes that bias DC estimates if unaccounted for.

Dosage compensation (DC) on the X chromosome has predictable effects on genetic and phenotypic trait variance. Here, the authors use information for 20 quantitative traits in the UK Biobank and across-tissue gene expression to compare X-linked heritability and the effects of trait-associated SNPs between the sexes.

"Arte et Labore"-A Blackburn Rovers fan's legacy in human complex trait genetics

Through his own research contributions on the modelling and genetic analysis of quantitative traits and through his former students and postdocs, Robin Thompson has indirectly left a major legacy in human genetics. In this short note, we highlight examples of the long-lasting relevance and impact of Robin's work in human genetics. A lone early study of marker-assisted selection developed many of the tools and approaches later exploited (often after reinvention) by the human genetics community in GWAS studies and for prediction. Furthermore, a particularly clear example of the pervasive impact of Robin's work is that REML has become the default method to estimate variance components and that genetic predictions exploiting linkage disequilibrium in the population are starting to become used in precision medicine applications.

Comprehensive Multiple eQTL Detection and Its Application to GWAS Interpretation

Expression QTL (eQTL) detection has emerged as an important tool for unraveling the relationship between genetic risk factors and disease or clinical phenotypes. Most studies are predicated on the assumption that only a single causal variant explains the association signal in each interval. This greatly simplifies the statistical modeling, but is liable to biases in scenarios where multiple local causal-variants are responsible. Here, our primary goal was to address the prevalence of secondary cis-eQTL signals regulating peripheral blood gene expression locally, utilizing two large human cohort studies, each >2500 samples with accompanying whole genome genotypes. The CAGE (Consortium for the Architecture of Gene Expression) dataset is a compendium of Illumina microarray studies, and the Framingham Heart Study is a two-generation Affymetrix dataset. We also describe Bayesian colocalization analysis of the extent of sharing of cis-eQTL detected in both studies as well as with the BIOS RNAseq dataset. Stepwise conditional modeling demonstrates that multiple eQTL signals are present for ∼40% of over 3500 eGenes in both microarray datasets, and that the number of loci with additional signals reduces by approximately two-thirds with each conditioning step. Although <20% of the peak signals across platforms fine map to the same credible interval, the colocalization analysis finds that as many as 50-60% of the primary eQTL are actually shared. Subsequently, colocalization of eQTL signals with GWAS hits detected 1349 genes whose expression in peripheral blood is associated with 591 human phenotype traits or diseases, including enrichment for genes with regulatory functions. At least 10%, and possibly as many as 40%, of eQTL-trait colocalized signals are due to nonprimary cis-eQTL peaks, but just one-quarter of these colocalization signals replicated across the gene expression datasets. Our results are provided as a web-based resource for visualization of multi-site regulation of gene expression and its association with human complex traits and disease states.

Parkinson's disease age at onset genome-wide association study: Defining heritability, genetic loci, and α-synuclein mechanisms

BACKGROUND

Increasing evidence supports an extensive and complex genetic contribution to PD. Previous genome-wide association studies (GWAS) have shed light on the genetic basis of risk for this disease. However, the genetic determinants of PD age at onset are largely unknown.

OBJECTIVES

To identify the genetic determinants of PD age at onset.

METHODS

Using genetic data of 28,568 PD cases, we performed a genome-wide association study based on PD age at onset.

RESULTS

We estimated that the heritability of PD age at onset attributed to common genetic variation was ∼0.11, lower than the overall heritability of risk for PD (∼0.27), likely, in part, because of the subjective nature of this measure. We found two genome-wide significant association signals, one at SNCA and the other a protein-coding variant in TMEM175, both of which are known PD risk loci and a Bonferroni-corrected significant effect at other known PD risk loci, GBA, INPP5F/BAG3, FAM47E/SCARB2, and MCCC1. Notably, SNCA, TMEM175, SCARB2, BAG3, and GBA have all been shown to be implicated in α-synuclein aggregation pathways. Remarkably, other well-established PD risk loci, such as GCH1 and MAPT, did not show a significant effect on age at onset of PD.

CONCLUSIONS

Overall, we have performed the largest age at onset of PD genome-wide association studies to date, and our results show that not all PD risk loci influence age at onset with significant differences between risk alleles for age at onset. This provides a compelling picture, both within the context of functional characterization of disease-linked genetic variability and in defining differences between risk alleles for age at onset, or frank risk for disease. © 2019 International Parkinson and Movement Disorder Society.

Publisher Correction: Multi-trait analysis of genome-wide association summary statistics using MTAG

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

OSCA: a tool for omic-data-based complex trait analysis

The rapid increase of omic data has greatly facilitated the investigation of associations between omic profiles such as DNA methylation (DNAm) and complex traits in large cohorts. Here, we propose a mixed-linear-model-based method called MOMENT that tests for association between a DNAm probe and trait with all other distal probes fitted in multiple random-effect components to account for unobserved confounders. We demonstrate by simulations that MOMENT shows a lower false positive rate and more robustness than existing methods. MOMENT has been implemented in a versatile software package called OSCA together with a number of other implementations for omic-data-based analyses.

Meta-analysis of genome-wide association studies for height and body mass index in ∼700000 individuals of European ancestry

Recent genome-wide association studies (GWAS) of height and body mass index (BMI) in ∼250000 European participants have led to the discovery of ∼700 and ∼100 nearly independent single nucleotide polymorphisms (SNPs) associated with these traits, respectively. Here we combine summary statistics from those two studies with GWAS of height and BMI performed in ∼450000 UK Biobank participants of European ancestry. Overall, our combined GWAS meta-analysis reaches N ∼700000 individuals and substantially increases the number of GWAS signals associated with these traits. We identified 3290 and 941 near-independent SNPs associated with height and BMI, respectively (at a revised genome-wide significance threshold of P < 1 × 10-8), including 1185 height-associated SNPs and 751 BMI-associated SNPs located within loci not previously identified by these two GWAS. The near-independent genome-wide significant SNPs explain ∼24.6% of the variance of height and ∼6.0% of the variance of BMI in an independent sample from the Health and Retirement Study (HRS). Correlations between polygenic scores based upon these SNPs with actual height and BMI in HRS participants were ∼0.44 and ∼0.22, respectively. From analyses of integrating GWAS and expression quantitative trait loci (eQTL) data by summary-data-based Mendelian randomization, we identified an enrichment of eQTLs among lead height and BMI signals, prioritizing 610 and 138 genes, respectively. Our study demonstrates that, as previously predicted, increasing GWAS sample sizes continues to deliver, by the discovery of new loci, increasing prediction accuracy and providing additional data to achieve deeper insight into complex trait biology. All summary statistics are made available for follow-up studies.

Bi-allelic Loss-of-Function CACNA1B Mutations in Progressive Epilepsy-Dyskinesia

The occurrence of non-epileptic hyperkinetic movements in the context of developmental epileptic encephalopathies is an increasingly recognized phenomenon. Identification of causative mutations provides an important insight into common pathogenic mechanisms that cause both seizures and abnormal motor control. We report bi-allelic loss-of-function CACNA1B variants in six children from three unrelated families whose affected members present with a complex and progressive neurological syndrome. All affected individuals presented with epileptic encephalopathy, severe neurodevelopmental delay (often with regression), and a hyperkinetic movement disorder. Additional neurological features included postnatal microcephaly and hypotonia. Five children died in childhood or adolescence (mean age of death: 9 years), mainly as a result of secondary respiratory complications. CACNA1B encodes the pore-forming subunit of the pre-synaptic neuronal voltage-gated calcium channel Ca_v2.2/N-type, crucial for SNARE-mediated neurotransmission, particularly in the early postnatal period. Bi-allelic loss-of-function variants in CACNA1B are predicted to cause disruption of Ca²⁺ influx, leading to impaired synaptic neurotransmission. The resultant effect on neuronal function is likely to be important in the development of involuntary movements and epilepsy. Overall, our findings provide further evidence for the key role of Ca_v2.2 in normal human neurodevelopment.

Genome-wide association study of medication-use and associated disease in the UK Biobank

Genome-wide association studies (GWASs) of medication use may contribute to understanding of disease etiology, could generate new leads relevant for drug discovery and can be used to quantify future risk of medication taking. Here, we conduct GWASs of self-reported medication use from 23 medication categories in approximately 320,000 individuals from the UK Biobank. A total of 505 independent genetic loci that meet stringent criteria (P < 10⁻⁸/23) for statistical significance are identified. We investigate the implications of these GWAS findings in relation to biological mechanism, potential drug target identification and genetic risk stratification of disease. Amongst the medication-associated genes are 16 known therapeutic-effect target genes for medications from 9 categories. Two of the medication classes studied are for disorders that have not previously been subject to large GWAS (hypothyroidism and gastro-oesophageal reflux disease).

An understanding of the genetic variants associated with medication use may shed light on the underlying biological pathways of disease, and aid in drug development. Here, Wu and colleagues conduct a GWAS for self-reported medication-use in the UK Biobank, finding more than 500 independent variants and many promising leads for future work.

Complex Trait Prediction from Genome Data: Contrasting EBV in Livestock to PRS in Humans: Genomic Prediction

In this Review, we focus on the similarity of the concepts underlying prediction of estimated breeding values (EBVs) in livestock and polygenic risk scores (PRS) in humans. Our research spans both fields and so we recognize factors that are very obvious for those in one field, but less so for those in the other. Differences in family size between species is the wedge that drives the different viewpoints and approaches. Large family size achievable in nonhuman species accompanied by selection generates a smaller effective population size, increased linkage disequilibrium and a higher average genetic relationship between individuals within a population. In human genetic analyses, we select individuals unrelated in the classical sense (coefficient of relationship <0.05) to estimate heritability captured by common SNPs. In livestock data, all animals within a breed are to some extent "related," and so it is not possible to select unrelated individuals and retain a data set of sufficient size to analyze. These differences directly or indirectly impact the way data analyses are undertaken. In livestock, genetic segregation variance exposed through samplings of parental genomes within families is directly observable and taken for granted. In humans, this genomic variation is under-recognized for its contribution to variation in polygenic risk of common disease, in both those with and without family history of disease. We explore the equation that predicts the expected proportion of variance explained using PRS, and quantify how GWAS sample size is the key factor for maximizing accuracy of prediction in both humans and livestock. Last, we bring together the concepts discussed to address some frequently asked questions.

From R.A. Fisher's 1918 Paper to GWAS a Century Later

The genetics and evolution of complex traits, including quantitative traits and disease, have been hotly debated ever since Darwin. A century ago, a paper from R.A. Fisher reconciled Mendelian and biometrical genetics in a landmark contribution that is now accepted as the main foundation stone of the field of quantitative genetics. Here, we give our perspective on Fisher's 1918 paper in the context of how and why it is relevant in today's genome era. We mostly focus on human trait variation, in part because Fisher did so too, but the conclusions are general and extend to other natural populations, and to populations undergoing artificial selection.

Common Disease Is More Complex Than Implied by the Core Gene Omnigenic Model

The evidence that most adult-onset common diseases have a polygenic genetic architecture fully consistent with robust biological systems supported by multiple back-up mechanisms is now overwhelming. In this context, we consider the recent "omnigenic" or "core genes" model. A key assumption of the model is that there is a relatively small number of core genes relevant to any disease. While intuitively appealing, this model may underestimate the biological complexity of common disease, and therefore, the goal to discover core genes should not guide experimental design. We consider other implications of polygenicity, concluding that a focus on patient stratification is needed to achieve the goals of precision medicine.

Meta-analysis of genome-wide association studies for body fat distribution in 694 649 individuals of European ancestry

More than one in three adults worldwide is either overweight or obese. Epidemiological studies indicate that the location and distribution of excess fat, rather than general adiposity, are more informative for predicting risk of obesity sequelae, including cardiometabolic disease and cancer. We performed a genome-wide association study meta-analysis of body fat distribution, measured by waist-to-hip ratio (WHR) adjusted for body mass index (WHRadjBMI), and identified 463 signals in 346 loci. Heritability and variant effects were generally stronger in women than men, and we found approximately one-third of all signals to be sexually dimorphic. The 5% of individuals carrying the most WHRadjBMI-increasing alleles were 1.62 times more likely than the bottom 5% to have a WHR above the thresholds used for metabolic syndrome. These data, made publicly available, will inform the biology of body fat distribution and its relationship with disease.

Assortative mating on complex traits revisited: Double first cousins and the X-chromosome

Mate choice through direct assortment on heritable traits, assortative mating (AM), is predicted in theory to inflate the genetic variance in a population and the correlation between relatives. Here, we revisit the theory of AM, first established in the landmark 1918 paper from RA Fisher, and provide new theory and analytical results. In particular, we shed light on inconsistencies in the literature regarding the correlation between double first cousins under AM and provide a solution. We derive new theory for AM due to X-chromosome loci. We show in the latter case that the inflation of genetic variance induced under AM is twice as large in females compared to males. These two theoretical contributions are verified and illustrated through simulations. We also provide a more general unified framework for the correlation between relatives in a non-inbred population.