Efficient Bayesian mixed-model analysis increases association power in large cohorts

Parent-of-Origin inference for biobanks

Identical genetic variations can have different phenotypic effects depending on their parent of origin. Yet, studies focusing on parent-of-origin effects have been limited in terms of sample size due to the lack of parental genomes or known genealogies. We propose a probabilistic approach to infer the parent-of-origin of individual alleles that does not require parental genomes nor prior knowledge of genealogy. Our model uses Identity-By-Descent sharing with second- and third-degree relatives to assign alleles to parental groups and leverages chromosome X data in males to distinguish maternal from paternal groups. We combine this with robust haplotype inference and haploid imputation to infer the parent-of-origin for 26,393 UK Biobank individuals. We screen 99 phenotypes for parent-of-origin effects and replicate the discoveries of 6 GWAS studies, confirming signals on body mass index, type 2 diabetes, standing height and multiple blood biomarkers, including the known maternal effect at the MEG3/DLK1 locus on platelet phenotypes. We also report a novel maternal effect at the TERT gene on telomere length, thereby providing new insights on the heritability of this phenotype. All our summary statistics are publicly available to help the community to better characterize the molecular mechanisms leading to parent-of-origin effects and their implications for human health.

Mosaic patterns of selection in genomic regions associated with diverse human traits

Natural selection shapes the genetic architecture of many human traits. However, the prevalence of different modes of selection on genomic regions associated with variation in traits remains poorly understood. To address this, we developed an efficient computational framework to calculate positive and negative enrichment of different evolutionary measures among regions associated with complex traits. We applied the framework to summary statistics from >900 genome-wide association studies (GWASs) and 11 evolutionary measures of sequence constraint, population differentiation, and allele age while accounting for linkage disequilibrium, allele frequency, and other potential confounders. We demonstrate that this framework yields consistent results across GWASs with variable sample sizes, numbers of trait-associated SNPs, and analytical approaches. The resulting evolutionary atlas maps diverse signatures of selection on genomic regions associated with complex human traits on an unprecedented scale. We detected positive enrichment for sequence conservation among trait-associated regions for the majority of traits (>77% of 290 high power GWASs), which included reproductive traits. Many traits also exhibited substantial positive enrichment for population differentiation, especially among hair, skin, and pigmentation traits. In contrast, we detected widespread negative enrichment for signatures of balancing selection (51% of GWASs) and absence of enrichment for evolutionary signals in regions associated with late-onset Alzheimer's disease. These results support a pervasive role for negative selection on regions of the human genome that contribute to variation in complex traits, but also demonstrate that diverse modes of evolution are likely to have shaped trait-associated loci. This atlas of evolutionary signatures across the diversity of available GWASs will enable exploration of the relationship between the genetic architecture and evolutionary processes in the human genome.

Rare and common genetic determinants of metabolic individuality and their effects on human health

Garrod's concept of 'chemical individuality' has contributed to comprehension of the molecular origins of human diseases. Untargeted high-throughput metabolomic technologies provide an in-depth snapshot of human metabolism at scale. We studied the genetic architecture of the human plasma metabolome using 913 metabolites assayed in 19,994 individuals and identified 2,599 variant-metabolite associations (P < 1.25 × 10^-11) within 330 genomic regions, with rare variants (minor allele frequency ≤ 1%) explaining 9.4% of associations. Jointly modeling metabolites in each region, we identified 423 regional, co-regulated, variant-metabolite clusters called genetically influenced metabotypes. We assigned causal genes for 62.4% of these genetically influenced metabotypes, providing new insights into fundamental metabolite physiology and clinical relevance, including metabolite-guided discovery of potential adverse drug effects (DPYD and SRD5A2). We show strong enrichment of inborn errors of metabolism-causing genes, with examples of metabolite associations and clinical phenotypes of non-pathogenic variant carriers matching characteristics of the inborn errors of metabolism. Systematic, phenotypic follow-up of metabolite-specific genetic scores revealed multiple potential etiological relationships.

Genetic and psychosocial influence on the association between early childhood infections and later psychiatric disorders

To evaluate the influence of extensive genetic and psychosocial confounding on the association between early childhood infection and five major psychiatric disorders METHODS: A case-cohort study including participants from the Danish iPSYCH2012 sample, a case-cohort sample where all cases born between May 1, 1981, and December 31, 2005, diagnosed with attention-deficit/hyperactivity disorder (ADHD), autism spectrum disorder (ASD), bipolar affective disorder (BIP), Major Depressive Disorder (MDD) or schizophrenia (SCZ), were identified and pooled with a representative sample (subcohort) of the Danish population. We used Cox proportional hazards regression customized to the case-cohort setup to calculate hazard ratios of outcome with 95% confidence intervals (CIs), following exposure to early childhood infection before the age of 5 years for ADHD and ASD, and before the age of 10 years for BIP, MDD, and SCZ. To evaluate psychosocial confounding we included sex, calendar period, sibling infections, urbanicity, parental socio-economic status, parental mental health information, and polygenic risk scores for all five disorders, as covariates. To estimate how liability for psychiatric disorders measured through the PRS influenced the risk of early childhood infection, we calculated odds ratios (ORs) with 95% CIs, using logistic regression RESULTS: Early childhood infection was associated with ADHD, ASD, MDD, and SCZ with number of childhood infections increasing the hazard. The HR was still significant in the model with full adjustments after 1 infection for ADHD (HR 1.29, 95% CI: 1.19-1.41), ASD (HR 1.28, 95% CI: 1.18-1.40), MDD (HR 1.23, 95% CI: 1.14-1.33), and SCZ (HR 1.21, 95% CI: 1.07-1.36), but not for BIP (HR1.17, 95% CI: 0.96-1.42). Probands exposed to sibling infections, but not own infection had an absolute risk of ADHD, BIP, MDD, and SCZ that closely approached the absolute risk for individuals exposed to own infections. We found evidence of gene-environment correlation with higher PRS of MDD and to some extent SCZ increasing the risk of infections and higher PRS of BIP associated with significantly decreased risk CONCLUSION: Early childhood infection is significantly associated with ADHD, ASD, MDD, and SCZ and not explained by genetic or psychosocial confounding. Although we found evidence of gene-environment correlation, it had minor impact on the results.

Assessing the contribution of genetic nurture to refractive error

Parents pass on both their genes and environment to offspring, prompting debate about the relative importance of nature versus nurture in the inheritance of complex traits. Advances in molecular genetics now make it possible to quantify an individual's genetic predisposition to a trait via his or her 'polygenic score'. However, part of the risk captured by an individual's polygenic score may actually be attributed to the genotype of their parents. In the most well-studied example of this indirect 'genetic nurture' effect, about half the genetic contribution to educational attainment was found to be attributed to parental alleles, even if those alleles were not inherited by the child. Refractive errors, such as myopia, are a common cause of visual impairment and pose high economic and quality-of-life costs. Despite strong evidence that refractive errors are highly heritable, the extent to which genetic risk is conferred directly via transmitted risk alleles or indirectly via the environment that parents create for their children is entirely unknown. Here, an instrumental variable analysis in 1944 pairs of adult siblings from the United Kingdom was used to quantify the proportion of the genetic risk ('single nucleotide polymorphism (SNP) heritability') of refractive error contributed by genetic nurture. We found no evidence of a contribution from genetic nurture: non-within-family SNP-heritability estimate = 0.213 (95% confidence interval 0.134-0.310) and within-family SNP-heritability estimate = 0.250 (0.152-0.372). Our findings imply the genetic contribution to refractive error is principally an intrinsic effect from alleles transmitted from parents to offspring.

The causal associations of circulating amino acids with blood pressure: a Mendelian randomization study

BACKGROUND

Circulating levels of amino acids were associated with blood pressure (BP) in observational studies. However, the causation of such associations has been hypothesized but is difficult to prove in human studies. Here, we aimed to use two-sample Mendelian randomization analyses to evaluate the potential causal associations of circulating levels of amino acids with BP and risk of hypertension.

METHODS

We generated genetic instruments for circulating levels of nine amino acids by conducting meta-analyses of genome-wide association study (GWAS) in UK Biobank participants with metabolomic data (n = 98,317) and another published metabolomics GWAS (n = 24,925). Data on the associations of the genetic variants with BP and hypertension were obtained in the UK Biobank participants without metabolomic data (n = 286,390). The causal effects were estimated using inverse-variance weighted method.

RESULTS

Significant evidence consistently supported the causal effects of increased branched-chain amino acids (BCAAs, i.e., leucine, isoleucine, and valine) levels on higher BP and risk of hypertension (all P < 0.006 after Bonferroni correction except for P_{leucine-on-diastolicBP} = 0.008). For example, per standard deviation higher of genetically predicted isoleucine levels were associated with 2.71 ± 0.78 mmHg higher systolic BP and 1.24 ± 0.34 mmHg higher diastolic BP, as well as with 7% higher risk of hypertension (odds ratio: 1.07, [95% CI: 1.04-1.10]). In addition, per standard deviation higher of genetically predicted glycine level was associated with lower systolic BP (- 0.70 ± 0.17 mmHg, P = 4.04 × 10^-5) and a lower risk of hypertension (0.99 [0.98-0.99], P = 6.46 × 10^-5). In the reverse direction, genetically predicted higher systolic BP was associated with lower circulating levels of glycine (- 0.025±0.008, P = 0.001).

CONCLUSIONS

This study provides evidence for causal impacts of genetically predicted circulating BCAAs and glycine levels on BP. Meanwhile, genetically predicted higher BP was associated with lower glycine levels. Further investigations are warranted to clarify the underlying mechanisms.

Influences of rare copy-number variation on human complex traits

The human genome contains hundreds of thousands of regions harboring copy-number variants (CNV). However, the phenotypic effects of most such polymorphisms are unknown because only larger CNVs have been ascertainable from SNP-array data generated by large biobanks. We developed a computational approach leveraging haplotype sharing in biobank cohorts to more sensitively detect CNVs. Applied to UK Biobank, this approach accounted for approximately half of all rare gene inactivation events produced by genomic structural variation. This CNV call set enabled a detailed analysis of associations between CNVs and 56 quantitative traits, identifying 269 independent associations (p < 5 × 10-8) likely to be causally driven by CNVs. Putative target genes were identifiable for nearly half of the loci, enabling insights into dosage sensitivity of these genes and uncovering several gene-trait relationships. These results demonstrate the ability of haplotype-informed analysis to provide insights into the genetic basis of human complex traits.

Genome-Wide Association Analysis of Over 170,000 Individuals from the UK Biobank Identifies Seven Loci Associated with Dietary Approaches to Stop Hypertension (DASH) Diet

Diet is a modifiable risk factor for common chronic diseases and mental health disorders, and its effects are under partial genetic control. To estimate the impact of diet on individual health, most epidemiological and genetic studies have focused on individual aspects of dietary intake. However, analysing individual food groups in isolation does not capture the complexity of the whole diet pattern. Dietary indices enable a holistic estimation of diet and account for the intercorrelations between food and nutrients. In this study we performed the first ever genome-wide association study (GWA) including 173,701 individuals from the UK Biobank to identify genetic variants associated with the Dietary Approaches to Stop Hypertension (DASH) diet. DASH was calculated using the 24 h-recall questionnaire collected by UK Biobank. The GWA was performed using a linear mixed model implemented in BOLT-LMM. We identified seven independent single-nucleotide polymorphisms (SNPs) associated with DASH. Significant genetic correlations were observed between DASH and several educational traits with a significant enrichment for genes involved in the AMP-dependent protein kinase (AMPK) activation that controls the appetite by regulating the signalling in the hypothalamus. The colocalization analysis implicates genes involved in body mass index (BMI)/obesity and neuroticism (ARPP21, RP11-62H7.2, MFHAS1, RHEBL1). The Mendelian randomisation analysis suggested that increased DASH score, which reflect a healthy diet style, is causal of lower glucose, and insulin levels. These findings further our knowledge of the pathways underlying the relationship between diet and health outcomes. They may have significant implications for global public health and provide future dietary recommendations for the prevention of common chronic diseases.

The HUNT study: A population-based cohort for genetic research

The Trøndelag Health Study (HUNT) is a population-based cohort of ∼229,000 individuals recruited in four waves beginning in 1984 in Trøndelag County, Norway. Approximately 88,000 of these individuals have available genetic data from array genotyping. HUNT participants were recruited during four community-based recruitment waves and provided information on health-related behaviors, self-reported diagnoses, family history of disease, and underwent physical examinations. Linkage via the Norwegian personal identification number integrates digitized health care information from doctor visits and national health registries including death, cancer and prescription registries. Genome-wide association studies of HUNT participants have provided insights into the mechanism of cardiovascular, metabolic, osteoporotic, and liver-related diseases, among others. Unique features of this cohort that facilitate research include nearly 40 years of longitudinal follow-up in a motivated and well-educated population, family data, comprehensive phenotyping, and broad availability of DNA, RNA, urine, fecal, plasma, and serum samples.

Constructing an atlas of associations between polygenic scores from across the human phenome and circulating metabolic biomarkers

Background

Polygenic scores (PGS) are becoming an increasingly popular approach to predict complex disease risk, although they also hold the potential to develop insight into the molecular profiles of patients with an elevated genetic predisposition to disease.

Methods

We sought to construct an atlas of associations between 125 different PGS derived using results from genome-wide association studies and 249 circulating metabolites in up to 83,004 participants from the UK Biobank.

Results

As an exemplar to demonstrate the value of this atlas, we conducted a hypothesis-free evaluation of all associations with glycoprotein acetyls (GlycA), an inflammatory biomarker. Using bidirectional Mendelian randomization, we find that the associations highlighted likely reflect the effect of risk factors, such as adiposity or liability towards smoking, on systemic inflammation as opposed to the converse direction. Moreover, we repeated all analyses in our atlas within age strata to investigate potential sources of collider bias, such as medication usage. This was exemplified by comparing associations between lipoprotein lipid profiles and the coronary artery disease PGS in the youngest and oldest age strata, which had differing proportions of individuals undergoing statin therapy. Lastly, we generated all PGS-metabolite associations stratified by sex and separately after excluding 13 established lipid-associated loci to further evaluate the robustness of findings.

Conclusions

We envisage that the atlas of results constructed in our study will motivate future hypothesis generation and help prioritize and deprioritize circulating metabolic traits for in-depth investigations. All results can be visualized and downloaded at http://mrcieu.mrsoftware.org/metabolites_PRS_atlas.

Funding

This work is supported by funding from the Wellcome Trust, the British Heart Foundation, and the Medical Research Council Integrative Epidemiology Unit.

Constructing an atlas of associations between polygenic scores from across the human phenome and circulating metabolic biomarkers

BACKGROUND

Polygenic scores (PGS) are becoming an increasingly popular approach to predict complex disease risk, although they also hold the potential to develop insight into the molecular profiles of patients with an elevated genetic predisposition to disease.

METHODS

We sought to construct an atlas of associations between 125 different PGS derived using results from genome-wide association studies and 249 circulating metabolites in up to 83,004 participants from the UK Biobank.

RESULTS

As an exemplar to demonstrate the value of this atlas, we conducted a hypothesis-free evaluation of all associations with glycoprotein acetyls (GlycA), an inflammatory biomarker. Using bidirectional Mendelian randomization, we find that the associations highlighted likely reflect the effect of risk factors, such as adiposity or liability towards smoking, on systemic inflammation as opposed to the converse direction. Moreover, we repeated all analyses in our atlas within age strata to investigate potential sources of collider bias, such as medication usage. This was exemplified by comparing associations between lipoprotein lipid profiles and the coronary artery disease PGS in the youngest and oldest age strata, which had differing proportions of individuals undergoing statin therapy. Lastly, we generated all PGS-metabolite associations stratified by sex and separately after excluding 13 established lipid-associated loci to further evaluate the robustness of findings.

CONCLUSIONS

We envisage that the atlas of results constructed in our study will motivate future hypothesis generation and help prioritize and deprioritize circulating metabolic traits for in-depth investigations. All results can be visualized and downloaded at http://mrcieu.mrsoftware.org/metabolites_PRS_atlas.

FUNDING

This work is supported by funding from the Wellcome Trust, the British Heart Foundation, and the Medical Research Council Integrative Epidemiology Unit.

Assessment of the bi-directional relationship between blood mitochondrial DNA copy number and type 2 diabetes mellitus: a multivariable-adjusted regression and Mendelian randomisation study

AIMS/HYPOTHESIS

Mitochondrial dysfunction, which can be approximated by blood mitochondrial DNA copy number (mtDNA-CN), has been implicated in the pathogenesis of type 2 diabetes mellitus. Thus far, however, insights from prospective cohort studies and Mendelian randomisation (MR) analyses on this relationship are limited. We assessed the association between blood mtDNA-CN and incident type 2 diabetes using multivariable-adjusted regression analyses, and the associations between blood mtDNA-CN and type 2 diabetes and BMI using bi-directional MR.

METHODS

Multivariable-adjusted Cox proportional hazard models were used to estimate the association between blood mtDNA-CN and incident type 2 diabetes in 285,967 unrelated European individuals from UK Biobank free of type 2 diabetes at baseline. Additionally, a cross-sectional analysis was performed to investigate the association between blood mtDNA-CN and BMI. We also assessed the potentially causal relationship between blood mtDNA-CN and type 2 diabetes (N=898,130 from DIAGRAM, N=215,654 from FinnGen) and BMI (N=681,275 from GIANT) using bi-directional two-sample MR.

RESULTS

During a median follow-up of 11.87 years, 15,111 participants developed type 2 diabetes. Participants with a higher level of blood mtDNA-CN are at lower risk of developing type 2 diabetes (HR 0.90 [95% CI 0.89, 0.92]). After additional adjustment for BMI and other confounders, these results attenuated moderately and remained present. The multivariable-adjusted cross-sectional analyses showed that higher blood mtDNA-CN was associated with lower BMI (-0.12 [95% CI -0.14, -0.10]) kg/m2. In the bi-directional MR analyses, we found no evidence for causal associations between blood mtDNA-CN and type 2 diabetes, and blood mtDNA-CN and BMI in either direction.

CONCLUSIONS/INTERPRETATION

The results from the present study indicate that the observed association between low blood mtDNA-CN and higher risk of type 2 diabetes is likely not causal.

COX7A2L genetic variants determine cardiorespiratory fitness in mice and human

Mitochondrial respiratory complexes form superassembled structures called supercomplexes. COX7A2L is a supercomplex-specific assembly factor in mammals, although its implication for supercomplex formation and cellular metabolism remains controversial. Here we identify a role for COX7A2L for mitochondrial supercomplex formation in humans. By using human cis-expression quantitative trait loci data, we highlight genetic variants in the COX7A2L gene that affect its skeletal muscle expression specifically. The most significant cis-expression quantitative trait locus is a 10-bp insertion in the COX7A2L 3' untranslated region that increases messenger RNA stability and expression. Human myotubes harboring this insertion have more supercomplexes and increased respiration. Notably, increased COX7A2L expression in the muscle is associated with lower body fat and improved cardiorespiratory fitness in humans. Accordingly, specific reconstitution of Cox7a2l expression in C57BL/6J mice leads to higher maximal oxygen consumption, increased lean mass and increased energy expenditure. Furthermore, Cox7a2l expression in mice is induced specifically in the muscle upon exercise. These findings elucidate the genetic basis of mitochondrial supercomplex formation and function in humans and show that COX7A2L plays an important role in cardiorespiratory fitness, which could have broad therapeutic implications in reducing cardiovascular mortality.

Including diverse and admixed populations in genetic epidemiology research

The inclusion of ancestrally diverse participants in genetic studies can lead to new discoveries and is important to ensure equitable health care benefit from research advances. Here, members of the Ethical, Legal, Social, Implications (ELSI) committee of the International Genetic Epidemiology Society (IGES) offer perspectives on methods and analysis tools for the conduct of inclusive genetic epidemiology research, with a focus on admixed and ancestrally diverse populations in support of reproducible research practices. We emphasize the importance of distinguishing socially defined population categorizations from genetic ancestry in the design, analysis, reporting, and interpretation of genetic epidemiology research findings. Finally, we discuss the current state of genomic resources used in genetic association studies, functional interpretation, and clinical and public health translation of genomic findings with respect to diverse populations.

Quantitative trait locus analysis for endophenotypes reveals genetic substrates of core symptom domains and neurocognitive function in autism spectrum disorder

Autism spectrum disorder (ASD) represents a heterogeneous group of neurodevelopmental disorders and is largely attributable to genetic risk factors. Phenotypic and genetic heterogeneity of ASD have been well-recognized; however, genetic substrates for endophenotypes that constitute phenotypic heterogeneity are not yet known. In the present study, we compiled data from the Autism Genetic Resource Exchange, which contains the demographic and detailed phenotype information of 11,961 individuals. Notably, the whole-genome sequencing data available from MSSNG and iHART for 3833 individuals in this dataset was used to perform an endophenotype-wide association study. Using a linear mixed model, genome-wide association analyses were performed for 29 endophenotype scores and 0.58 million common variants with variant allele frequency ≥ 5%. We discovered significant associations between 9 genetic variants and 6 endophenotype scores comprising neurocognitive development and severity scores for core symptoms of ASD at a significance threshold of p < 5 × 10^-7. Of note, the Stereotyped Behaviors and Restricted Interests total score in Autism Diagnostic Observation Schedule Module 3 was significantly associated with multiple variants in the VPS13B gene, a causal gene for Cohen syndrome and a candidate gene for syndromic ASD. Our findings yielded loci with small effect sizes due to the moderate sample size and, thus, require validation in another cohort. Nonetheless, our endophenotype-wide association analysis extends previous candidate gene discovery in the context of genotype and endophenotype association. As a result, these candidate genes may be responsible for specific traits that constitute core symptoms and neurocognitive function of ASD rather than the disorder itself.

Efficient permutation-based genome-wide association studies for normal and skewed phenotypic distributions

MOTIVATION

Genome-wide association studies (GWAS) are an integral tool for studying the architecture of complex genotype and phenotype relationships. Linear mixed models (LMMs) are commonly used to detect associations between genetic markers and a trait of interest, while at the same time allowing to account for population structure and cryptic relatedness. Assumptions of LMMs include a normal distribution of the residuals and that the genetic markers are independent and identically distributed-both assumptions are often violated in real data. Permutation-based methods can help to overcome some of these limitations and provide more realistic thresholds for the discovery of true associations. Still, in practice, they are rarely implemented due to the high computational complexity.

RESULTS

We propose permGWAS, an efficient LMM reformulation based on 4D tensors that can provide permutation-based significance thresholds. We show that our method outperforms current state-of-the-art LMMs with respect to runtime and that permutation-based thresholds have lower false discovery rates for skewed phenotypes compared to the commonly used Bonferroni threshold. Furthermore, using permGWAS we re-analyzed more than 500 Arabidopsis thaliana phenotypes with 100 permutations each in less than 8 days on a single GPU. Our re-analyses suggest that applying a permutation-based threshold can improve and refine the interpretation of GWAS results.

AVAILABILITY AND IMPLEMENTATION

permGWAS is open-source and publicly available on GitHub for download: https://github.com/grimmlab/permGWAS.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Efficient and accurate frailty model approach for genome-wide survival association analysis in large-scale biobanks

With decades of electronic health records linked to genetic data, large biobanks provide unprecedented opportunities for systematically understanding the genetics of the natural history of complex diseases. Genome-wide survival association analysis can identify genetic variants associated with ages of onset, disease progression and lifespan. We propose an efficient and accurate frailty model approach for genome-wide survival association analysis of censored time-to-event (TTE) phenotypes by accounting for both population structure and relatedness. Our method utilizes state-of-the-art optimization strategies to reduce the computational cost. The saddlepoint approximation is used to allow for analysis of heavily censored phenotypes (>90%) and low frequency variants (down to minor allele count 20). We demonstrate the performance of our method through extensive simulation studies and analysis of five TTE phenotypes, including lifespan, with heavy censoring rates (90.9% to 99.8%) on ~400,000 UK Biobank participants with white British ancestry and ~180,000 individuals in FinnGen. We further analyzed 871 TTE phenotypes in the UK Biobank and presented the genome-wide scale phenome-wide association results with the PheWeb browser.

Incorporating family disease history and controlling case-control imbalance for population-based genetic association studies

MOTIVATION

In the genome-wide association analysis of population-based biobanks, most diseases have low prevalence, which results in low detection power. One approach to tackle the problem is using family disease history, yet existing methods are unable to address type I error inflation induced by increased correlation of phenotypes among closely related samples, as well as unbalanced phenotypic distribution.

RESULTS

We propose a new method for genetic association test with family disease history, mixed-model-based Test with Adjusted Phenotype and Empirical saddlepoint approximation, which controls for increased phenotype correlation by adopting a two-variance-component mixed model, accounts for case-control imbalance by using empirical saddlepoint approximation, and is flexible to incorporate any existing adjusted phenotypes, such as phenotypes from the LT-FH method. We show through simulation studies and analysis of UK Biobank data of white British samples and the Korean Genome and Epidemiology Study of Korean samples that the proposed method is robust and yields better calibration compared to existing methods while gaining power for detection of variant-phenotype associations.

AVAILABILITY AND IMPLEMENTATION

The summary statistics and code generated in this study are available at https://github.com/styvon/TAPE.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Cold shock domain-containing protein E1 is a posttranscriptional regulator of the LDL receptor

The low-density lipoprotein receptor (LDLR) controls cellular delivery of cholesterol and clears LDL from the bloodstream, protecting against atherosclerotic heart disease, the leading cause of death in the United States. We therefore sought to identify regulators of the LDLR beyond the targets of current therapies and known causes of familial hypercholesterolemia. We found that cold shock domain-containing protein E1 (CSDE1) enhanced hepatic LDLR messenger RNA (mRNA) decay via its 3' untranslated region and regulated atherogenic lipoproteins in vivo. Using parallel phenotypic genome-wide CRISPR interference screens in a tissue culture model, we identified 40 specific regulators of the LDLR that were not previously identified by observational human genetic studies. Among these, we demonstrated that, in HepG2 cells, CSDE1 regulated the LDLR at least as strongly as statins and proprotein convertase subtilisin/kexin type 9 (PCSK9) inhibitors. In addition, we showed that hepatic gene silencing of Csde1 treated diet-induced dyslipidemia in mice to a similar degree as Pcsk9 silencing. These results suggest the therapeutic potential of targeting CSDE1 to manipulate the posttranscriptional regulation of the LDLR mRNA for the prevention of cardiovascular disease. Our approach of modeling a clinically relevant phenotype in a forward genetic screen, followed by mechanistic pharmacologic dissection and in vivo validation, may serve as a generalizable template for the identification of therapeutic targets in other human disease states.

Genetic risk factors have a substantial impact on healthy life years

The impact of genetic variation on overall disease burden has not been comprehensively evaluated. We introduce an approach to estimate the effect of genetic risk factors on disability-adjusted life years (DALYs; 'lost healthy life years'). We use genetic information from 735,748 individuals and consider 80 diseases. Rare variants had the highest effect on DALYs at the individual level. Among common variants, rs3798220 (LPA) had the strongest individual-level effect, with 1.18 DALYs from carrying 1 versus 0 copies. Being in the top 10% versus the bottom 90% of a polygenic score for multisite chronic pain had an effect of 3.63 DALYs. Some common variants had a population-level effect comparable to modifiable risk factors such as high sodium intake and low physical activity. Attributable DALYs vary between males and females for some genetic exposures. Genetic risk factors can explain a sizable number of healthy life years lost both at the individual and population level.

Integrative analysis of metabolite GWAS illuminates the molecular basis of pleiotropy and genetic correlation

Pleiotropy and genetic correlation are widespread features in genome-wide association studies (GWAS), but they are often difficult to interpret at the molecular level. Here, we perform GWAS of 16 metabolites clustered at the intersection of amino acid catabolism, glycolysis, and ketone body metabolism in a subset of UK Biobank. We utilize the well-documented biochemistry jointly impacting these metabolites to analyze pleiotropic effects in the context of their pathways. Among the 213 lead GWAS hits, we find a strong enrichment for genes encoding pathway-relevant enzymes and transporters. We demonstrate that the effect directions of variants acting on biology between metabolite pairs often contrast with those of upstream or downstream variants as well as the polygenic background. Thus, we find that these outlier variants often reflect biology local to the traits. Finally, we explore the implications for interpreting disease GWAS, underscoring the potential of unifying biochemistry with dense metabolomics data to understand the molecular basis of pleiotropy in complex traits and diseases.

Genome-wide association study of the human brain functional connectome reveals strong vascular component underlying global network efficiency

Complex brain networks play a central role in integrating activity across the human brain, and such networks can be identified in the absence of any external stimulus. We performed 10 genome-wide association studies of resting state network measures of intrinsic brain activity in up to 36,150 participants of European ancestry in the UK Biobank. We found that the heritability of global network efficiency was largely explained by blood oxygen level-dependent (BOLD) resting state fluctuation amplitudes (RSFA), which are thought to reflect the vascular component of the BOLD signal. RSFA itself had a significant genetic component and we identified 24 genomic loci associated with RSFA, 157 genes whose predicted expression correlated with it, and 3 proteins in the dorsolateral prefrontal cortex and 4 in plasma. We observed correlations with cardiovascular traits, and single-cell RNA specificity analyses revealed enrichment of vascular related cells. Our analyses also revealed a potential role of lipid transport, store-operated calcium channel activity, and inositol 1,4,5-trisphosphate binding in resting-state BOLD fluctuations. We conclude that that the heritability of global network efficiency is largely explained by the vascular component of the BOLD response as ascertained by RSFA, which itself has a significant genetic component.

Brain-specific genes contribute to chronic but not to acute back pain

Introduction

Back pain is the leading cause of disability worldwide. Although most back pain cases are acute, 20% of acute pain patients experience chronic back pain symptoms. It is unclear whether acute pain and chronic pain have similar or distinct underlying genetic mechanisms.

Objectives

To characterize the molecular and cellular pathways contributing to acute and chronic pain states.

Methods

Cross-sectional observational genome-wide association study.

Results

A total of 375,158 individuals from the UK Biobank cohort were included in the discovery of genome-wide association study. Of those, 70,633 (19%) and 32,209 (9%) individuals met the definition of chronic and acute back pain, respectively. A total of 355 single nucleotide polymorphism grouped into 13 loci reached the genome-wide significance threshold (5x10^-8) for chronic back pain, but none for acute. Of these, 7 loci were replicated in the Nord-Trøndelag Health Study (HUNT) cohort (19,760 chronic low back pain cases and 28,674 pain-free controls). Single nucleotide polymorphism heritability was 4.6% (P=1.4x10^-78) for chronic back pain and 0.81% (P=1.4x10-8) for acute back pain. Similar differences in heritability estimates between acute and chronic back pain were found in the HUNT cohort: 3.4% (P=0.0011) and 0.6% (P=0.851), respectively. Pathway analyses, tissue-specific heritability enrichment analyses, and epigenetic characterization suggest a substantial genetic contribution to chronic but not acute back pain from the loci predominantly expressed in the central nervous system.

Conclusion

Chronic back pain is substantially more heritable than acute back pain. This heritability is mostly attributed to genes expressed in the brain.

Detection and characterization of male sex chromosome abnormalities in the UK Biobank study

PURPOSE

The study aimed to systematically ascertain male sex chromosome abnormalities, 47,XXY (Klinefelter syndrome [KS]) and 47,XYY, and characterize their risks of adverse health outcomes.

METHODS

We analyzed genotyping array or exome sequence data in 207,067 men of European ancestry aged 40 to 70 years from the UK Biobank and related these to extensive routine health record data.

RESULTS

Only 49 of 213 (23%) of men whom we identified with KS and only 1 of 143 (0.7%) with 47,XYY had a diagnosis of abnormal karyotype on their medical records or self-report. We observed expected associations for KS with reproductive dysfunction (late puberty: risk ratio [RR] = 2.7; childlessness: RR = 4.2; testosterone concentration: RR = -3.8 nmol/L, all P < 2 × 10-8), whereas XYY men appeared to have normal reproductive function. Despite this difference, we identified several higher disease risks shared across both KS and 47,XYY, including type 2 diabetes (RR = 3.0 and 2.6, respectively), venous thrombosis (RR = 6.4 and 7.4, respectively), pulmonary embolism (RR = 3.3 and 3.7, respectively), and chronic obstructive pulmonary disease (RR = 4.4 and 4.6, respectively) (all P < 7 × 10-6).

CONCLUSION

KS and 47,XYY were mostly unrecognized but conferred substantially higher risks for metabolic, vascular, and respiratory diseases, which were only partially explained by higher levels of body mass index, deprivation, and smoking.

Genetic analyses of the electrocardiographic QT interval and its components identify additional loci and pathways

The QT interval is an electrocardiographic measure representing the sum of ventricular depolarization and repolarization, estimated by QRS duration and JT interval, respectively. QT interval abnormalities are associated with potentially fatal ventricular arrhythmia. Using genome-wide multi-ancestry analyses (>250,000 individuals) we identify 177, 156 and 121 independent loci for QT, JT and QRS, respectively, including a male-specific X-chromosome locus. Using gene-based rare-variant methods, we identify associations with Mendelian disease genes. Enrichments are observed in established pathways for QT and JT, and previously unreported genes indicated in insulin-receptor signalling and cardiac energy metabolism. In contrast for QRS, connective tissue components and processes for cell growth and extracellular matrix interactions are significantly enriched. We demonstrate polygenic risk score associations with atrial fibrillation, conduction disease and sudden cardiac death. Prioritization of druggable genes highlight potential therapeutic targets for arrhythmia. Together, these results substantially advance our understanding of the genetic architecture of ventricular depolarization and repolarization.

Genome-Wide Association Study Identifies New Genetic Determinants of Cardiorespiratory Fitness: The Trøndelag Health Study

PURPOSE

Low cardiorespiratory fitness (CRF) is a major risk factor for cardiovascular disease (CVD) and a stronger predictor of CVD morbidity and mortality than established risk factors. The genetic component of CRF, quantified as peak oxygen uptake (V̇O 2peak ), is estimated to be ~60%. Unfortunately, current studies on genetic markers for CRF have been limited by small sample sizes and using estimated CRF. To overcome these limitations, we performed a large-scale systematic screening for genetic variants associated with V̇O 2peak .

METHODS

A genome-wide association study was performed with BOLT-LMM including directly measured V̇O 2peak from 4525 participants in the HUNT3 Fitness study and 14 million single-nucleotide polymorphisms (SNP). For validation, similar analyses were performed in the United Kingdom Biobank (UKB), where CRF was assessed through a submaximal bicycle test, including ~60,000 participants and ~60 million SNP. Functional mapping and annotation of the genome-wide association study results was conducted using FUMA.

RESULTS

In HUNT, two genome-wide significant SNP associated with V̇O 2peak were identified in the total population, two in males, and 35 in females. Two SNP in the female population showed nominally significant association in the UKB. One of the replicated SNP is located in PIK3R5 , shown to be of importance for cardiac function and CVD. Bioinformatic analyses of the total and male population revealed candidate SNP in PPP3CA , previously associated with CRF.

CONCLUSIONS

We identified 38 novel SNP associated with V̇O 2peak in HUNT. Two SNP were nominally replicated in UKB. Several interesting genes emerged from the functional analyses, among them one previously reported to be associated with CVD and another with CRF.

Identification and single-base gene-editing functional validation of a cis-EPO variant as a genetic predictor for EPO-increasing therapies

Hypoxia-inducible factor prolyl hydroxylase inhibitors (HIF-PHIs) are currently under clinical development for treating anemia in chronic kidney disease (CKD), but it is important to monitor their cardiovascular safety. Genetic variants can be used as predictors to help inform the potential risk of adverse effects associated with drug treatments. We therefore aimed to use human genetics to help assess the risk of adverse cardiovascular events associated with therapeutically altered EPO levels to help inform clinical trials studying the safety of HIF-PHIs. By performing a genome-wide association meta-analysis of EPO (n = 6,127), we identified a cis-EPO variant (rs1617640) lying in the EPO promoter region. We validated this variant as most likely causal in controlling EPO levels by using genetic and functional approaches, including single-base gene editing. Using this variant as a partial predictor for therapeutic modulation of EPO and large genome-wide association data in Mendelian randomization tests, we found no evidence (at p < 0.05) that genetically predicted long-term rises in endogenous EPO, equivalent to a 2.2-unit increase, increased risk of coronary artery disease (CAD, OR [95% CI] = 1.01 [0.93, 1.07]), myocardial infarction (MI, OR [95% CI] = 0.99 [0.87, 1.15]), or stroke (OR [95% CI] = 0.97 [0.87, 1.07]). We could exclude increased odds of 1.15 for cardiovascular disease for a 2.2-unit EPO increase. A combination of genetic and functional studies provides a powerful approach to investigate the potential therapeutic profile of EPO-increasing therapies for treating anemia in CKD.

Genetic evidence for a causal relationship between type 2 diabetes and peripheral artery disease in both Europeans and East Asians

BACKGROUND

Observational studies have revealed that type 2 diabetes (T2D) is associated with an increased risk of peripheral artery disease (PAD). However, whether the two diseases share a genetic basis and whether the relationship is causal remain unclear. It is also unclear as to whether these relationships differ between ethnic groups.

METHODS

By leveraging large-scale genome-wide association study (GWAS) summary statistics of T2D (European-based: N_case = 21,926, N_control = 342,747; East Asian-based: N_case = 36,614, N_control = 155,150) and PAD (European-based: N_case = 5673, N_control = 359,551; East Asian-based: N_case = 3593, N_control = 208,860), we explored the genetic correlation and putative causal relationship between T2D and PAD in both Europeans and East Asians using linkage disequilibrium score regression and seven Mendelian randomization (MR) models. We also performed multi-trait analysis of GWAS and two gene-based analyses to reveal candidate variants and risk genes involved in the shared genetic basis between T2D and PAD.

RESULTS

We observed a strong genetic correlation (r_g) between T2D and PAD in both Europeans (r_g = 0.51; p-value = 9.34 × 10^-15) and East Asians (r_g = 0.46; p-value = 1.67 × 10^-12). The MR analyses provided consistent evidence for a causal effect of T2D on PAD in both ethnicities (odds ratio [OR] = 1.05 to 1.28 for Europeans and 1.15 to 1.27 for East Asians) but not PAD on T2D. This putative causal effect was not influenced by total cholesterol, body mass index, systolic blood pressure, or smoking initiation according to multivariable MR analysis, and the genetic overlap between T2D and PAD was further explored employing an independent European sample through polygenic risk score regression. Multi-trait analysis of GWAS revealed two novel European-specific single nucleotide polymorphisms (rs927742 and rs1734409) associated with the shared genetic basis of T2D and PAD. Gene-based analyses consistently identified one gene ANKFY1 and gene-gene interactions (e.g., STARD10 [European-specific] to AP3S2 [East Asian-specific]; KCNJ11 [European-specific] to KCNQ1 [East Asian-specific]) associated with the trans-ethnic genetic overlap between T2D and PAD, reflecting a common genetic basis for the co-occurrence of T2D and PAD in both Europeans and East Asians.

CONCLUSIONS

Our study provides the first evidence for a genetically causal effect of T2D on PAD in both Europeans and East Asians. Several candidate variants and risk genes were identified as being associated with this genetic overlap. Our findings emphasize the importance of monitoring PAD status in T2D patients and suggest new genetic biomarkers for screening PAD risk among patients with T2D.

Polygenic Risk Scores for Cardiovascular Disease: A Scientific Statement From the American Heart Association

Cardiovascular disease is the leading contributor to years lost due to disability or premature death among adults. Current efforts focus on risk prediction and risk factor mitigation' which have been recognized for the past half-century. However, despite advances, risk prediction remains imprecise with persistently high rates of incident cardiovascular disease. Genetic characterization has been proposed as an approach to enable earlier and potentially tailored prevention. Rare mendelian pathogenic variants predisposing to cardiometabolic conditions have long been known to contribute to disease risk in some families. However, twin and familial aggregation studies imply that diverse cardiovascular conditions are heritable in the general population. Significant technological and methodological advances since the Human Genome Project are facilitating population-based comprehensive genetic profiling at decreasing costs. Genome-wide association studies from such endeavors continue to elucidate causal mechanisms for cardiovascular diseases. Systematic cataloging for cardiovascular risk alleles also enabled the development of polygenic risk scores. Genetic profiling is becoming widespread in large-scale research, including in health care-associated biobanks, randomized controlled trials, and direct-to-consumer profiling in tens of millions of people. Thus, individuals and their physicians are increasingly presented with polygenic risk scores for cardiovascular conditions in clinical encounters. In this scientific statement, we review the contemporary science, clinical considerations, and future challenges for polygenic risk scores for cardiovascular diseases. We selected 5 cardiometabolic diseases (coronary artery disease, hypercholesterolemia, type 2 diabetes, atrial fibrillation, and venous thromboembolic disease) and response to drug therapy and offer provisional guidance to health care professionals, researchers, policymakers, and patients.

The causes and consequences of Alzheimer's disease: phenome-wide evidence from Mendelian randomization

Alzheimer's disease (AD) has no proven causal and modifiable risk factors, or effective interventions. We report a phenome-wide association study (PheWAS) of genetic liability for AD in 334,968 participants of the UK Biobank study, stratified by age. We also examined the effects of AD genetic liability on previously implicated risk factors. We replicated these analyses in the HUNT study. PheWAS hits and previously implicated risk factors were followed up in a Mendelian randomization (MR) framework to identify the causal effect of each risk factor on AD risk. A higher genetic liability for AD was associated with medical history and cognitive, lifestyle, physical and blood-based measures as early as 39 years of age. These effects were largely driven by the APOE gene. The follow-up MR analyses were primarily null, implying that most of these associations are likely to be a consequence of prodromal disease or selection bias, rather than the risk factor causing the disease.

Genome-wide associations of aortic distensibility suggest causality for aortic aneurysms and brain white matter hyperintensities

Aortic dimensions and distensibility are key risk factors for aortic aneurysms and dissections, as well as for other cardiovascular and cerebrovascular diseases. We present genome-wide associations of ascending and descending aortic distensibility and area derived from cardiac magnetic resonance imaging (MRI) data of up to 32,590 Caucasian individuals in UK Biobank. We identify 102 loci (including 27 novel associations) tagging genes related to cardiovascular development, extracellular matrix production, smooth muscle cell contraction and heritable aortic diseases. Functional analyses highlight four signalling pathways associated with aortic distensibility (TGF-β, IGF, VEGF and PDGF). We identify distinct sex-specific associations with aortic traits. We develop co-expression networks associated with aortic traits and apply phenome-wide Mendelian randomization (MR-PheWAS), generating evidence for a causal role for aortic distensibility in development of aortic aneurysms. Multivariable MR suggests a causal relationship between aortic distensibility and cerebral white matter hyperintensities, mechanistically linking aortic traits and brain small vessel disease.

Improving GWAS discovery and genomic prediction accuracy in biobank data

Genetically informed, deep-phenotyped biobanks are an important research resource and it is imperative that the most powerful, versatile, and efficient analysis approaches are used. Here, we apply our recently developed Bayesian grouped mixture of regressions model (GMRM) in the UK and Estonian Biobanks and obtain the highest genomic prediction accuracy reported to date across 21 heritable traits. When compared to other approaches, GMRM accuracy was greater than annotation prediction models run in the LDAK or LDPred-funct software by 15% (SE 7%) and 14% (SE 2%), respectively, and was 18% (SE 3%) greater than a baseline BayesR model without single-nucleotide polymorphism (SNP) markers grouped into minor allele frequency-linkage disequilibrium (MAF-LD) annotation categories. For height, the prediction accuracy R² was 47% in a UK Biobank holdout sample, which was 76% of the estimated [Formula: see text]. We then extend our GMRM prediction model to provide mixed-linear model association (MLMA) SNP marker estimates for genome-wide association (GWAS) discovery, which increased the independent loci detected to 16,162 in unrelated UK Biobank individuals, compared to 10,550 from BoltLMM and 10,095 from Regenie, a 62 and 65% increase, respectively. The average [Formula: see text] value of the leading markers increased by 15.24 (SE 0.41) for every 1% increase in prediction accuracy gained over a baseline BayesR model across the traits. Thus, we show that modeling genetic associations accounting for MAF and LD differences among SNP markers, and incorporating prior knowledge of genomic function, is important for both genomic prediction and discovery in large-scale individual-level studies.

Genome-wide analyses of 200,453 individuals yield new insights into the causes and consequences of clonal hematopoiesis

Clonal hematopoiesis (CH), the clonal expansion of a blood stem cell and its progeny driven by somatic driver mutations, affects over a third of people, yet remains poorly understood. Here we analyze genetic data from 200,453 UK Biobank participants to map the landscape of inherited predisposition to CH, increasing the number of germline associations with CH in European-ancestry populations from 4 to 14. Genes at new loci implicate DNA damage repair (PARP1, ATM, CHEK2), hematopoietic stem cell migration/homing (CD164) and myeloid oncogenesis (SETBP1). Several associations were CH-subtype-specific including variants at TCL1A and CD164 that had opposite associations with DNMT3A- versus TET2-mutant CH, the two most common CH subtypes, proposing key roles for these two loci in CH development. Mendelian randomization analyses showed that smoking and longer leukocyte telomere length are causal risk factors for CH and that genetic predisposition to CH increases risks of myeloproliferative neoplasia, nonhematological malignancies, atrial fibrillation and blood epigenetic ageing.

CellRegMap: a statistical framework for mapping context-specific regulatory variants using scRNA-seq

Single‐cell RNA sequencing (scRNA‐seq) enables characterizing the cellular heterogeneity in human tissues. Recent technological advances have enabled the first population‐scale scRNA‐seq studies in hundreds of individuals, allowing to assay genetic effects with single‐cell resolution. However, existing strategies to analyze these data remain based on principles established for the genetic analysis of bulk RNA‐seq. In particular, current methods depend on a priori definitions of discrete cell types, and hence cannot assess allelic effects across subtle cell types and cell states. To address this, we propose the Cell Regulatory Map (CellRegMap), a statistical framework to test for and quantify genetic effects on gene expression in individual cells. CellRegMap provides a principled approach to identify and characterize genotype–context interactions of known eQTL variants using scRNA‐seq data. This model‐based approach resolves allelic effects across cellular contexts of different granularity, including genetic effects specific to cell subtypes and continuous cell transitions. We validate CellRegMap using simulated data and apply it to previously identified eQTL from two recent studies of differentiating iPSCs, where we uncover hundreds of eQTL displaying heterogeneity of genetic effects across cellular contexts. Finally, we identify fine‐grained genetic regulation in neuronal subtypes for eQTL that are colocalized with human disease variants.

HOTAIR interacts with PRC2 complex regulating the regional preadipocyte transcriptome and human fat distribution

Mechanisms governing regional human adipose tissue (AT) development remain undefined. Here, we show that the long non-coding RNA HOTAIR (HOX transcript antisense RNA) is exclusively expressed in gluteofemoral AT, where it is essential for adipocyte development. We find that HOTAIR interacts with polycomb repressive complex 2 (PRC2) and we identify core HOTAIR-PRC2 target genes involved in adipocyte lineage determination. Repression of target genes coincides with PRC2 promoter occupancy and H3K27 trimethylation. HOTAIR is also involved in modifying the gluteal adipocyte transcriptome through alternative splicing. Gluteal-specific expression of HOTAIR is maintained by defined regions of open chromatin across the HOTAIR promoter. HOTAIR expression levels can be modified by hormonal (estrogen, glucocorticoids) and genetic variation (rs1443512 is a HOTAIR eQTL associated with reduced gynoid fat mass). These data identify HOTAIR as a dynamic regulator of the gluteal adipocyte transcriptome and epigenome with functional importance for human regional AT development.

Transcriptome-wide association study in UK Biobank Europeans identifies associations with blood cell traits

Previous genome-wide association studies (GWAS) of hematological traits have identified over 10 000 distinct trait-specific risk loci. However, at these loci, the underlying causal mechanisms remain incompletely characterized. To elucidate novel biology and better understand causal mechanisms at known loci, we performed a transcriptome-wide association study (TWAS) of 29 hematological traits in 399 835 UK Biobank (UKB) participants of European ancestry using gene expression prediction models trained from whole blood RNA-seq data in 922 individuals. We discovered 557 gene-trait associations for hematological traits distinct from previously reported GWAS variants in European populations. Among the 557 associations, 301 were available for replication in a cohort of 141 286 participants of European ancestry from the Million Veteran Program. Of these 301 associations, 108 replicated at a strict Bonferroni adjusted threshold ($\alpha$= 0.05/301). Using our TWAS results, we systematically assigned 4261 out of 16 900 previously identified hematological trait GWAS variants to putative target genes. Compared to coloc, our TWAS results show reduced specificity and increased sensitivity in external datasets to assign variants to target genes.

Improving polygenic prediction with genetically inferred ancestry

Genome-wide association studies (GWASs) have demonstrated that most common diseases have a strong genetic component from many genetic variants each with a small effect size. GWAS summary statistics have allowed the construction of polygenic scores (PGSs) estimating part of the individual risk for common diseases. Here, we propose to improve PGS-based risk estimation by incorporating genetic ancestry derived from genome-wide genotyping data. Our method involves three cohorts: a base (or discovery) for association studies, a target for phenotype/risk prediction, and a map for ancestry mapping; successively, (1) it generates for each individual in the base and target cohorts a set of principal components based on the map cohort-called mapped PCs, (2) it associates in the base cohort the phenotype with the mapped-PCs, and (3) it uses the mapped PCs in the target cohort to generate a phenotypic predictor called the ancestry score. We evaluated the ancestry score by comparing a predictive model using a PGS with one combining a PGS and an ancestry score. First, we performed simulations and found that the ancestry score has a greater impact on traits that correlate with ancestry-specific variants. Second, we showed, using UK Biobank data, that the ancestry score improves genetic prediction for our nine phenotypes to very different degrees. Third, we performed simulations and found that the more heterogeneous the base and target cohorts, the more beneficial the ancestry score is. Finally, we validated our approach under realistic conditions with UK Biobank as the base cohort and Swiss individuals from the CoLaus|PsyCoLaus study as the target cohort.

SNP-to-gene linking strategies reveal contributions of enhancer-related and candidate master-regulator genes to autoimmune disease

We assess contributions to autoimmune disease of genes whose regulation is driven by enhancer regions (enhancer-related) and genes that regulate other genes in trans (candidate master-regulator). We link these genes to SNPs using several SNP-to-gene (S2G) strategies and apply heritability analyses to draw three conclusions about 11 autoimmune/blood-related diseases/traits. First, several characterizations of enhancer-related genes using functional genomics data are informative for autoimmune disease heritability after conditioning on a broad set of regulatory annotations. Second, candidate master-regulator genes defined using trans-eQTL in blood are also conditionally informative for autoimmune disease heritability. Third, integrating enhancer-related and master-regulator gene sets with protein-protein interaction (PPI) network information magnified their disease signal. The resulting PPI-enhancer gene score produced >2-fold stronger heritability signal and >2-fold stronger enrichment for drug targets, compared with the recently proposed enhancer domain score. In each case, functionally informed S2G strategies produced 4.1- to 13-fold stronger disease signals than conventional window-based strategies.

Incorporating family history of disease improves polygenic risk scores in diverse populations

Polygenic risk scores (PRSs) derived from genotype data and family history (FH) of disease provide valuable information for predicting disease risk, but PRSs perform poorly when applied to diverse populations. Here, we explore methods for combining both types of information (PRS-FH) in UK Biobank data. PRSs were trained using all British individuals (n = 409,000), and target samples consisted of unrelated non-British Europeans (n = 42,000), South Asians (n = 7,000), or Africans (n = 7,000). We evaluated PRS, FH, and PRS-FH using liability-scale R 2, primarily focusing on 3 well-powered diseases (type 2 diabetes, hypertension, and depression). PRS attained average prediction R 2s of 5.8%, 4.0%, and 0.53% in non-British Europeans, South Asians, and Africans, confirming poor cross-population transferability. In contrast, PRS-FH attained average prediction R 2s of 13%, 12%, and 10%, respectively, representing a large improvement in Europeans and an extremely large improvement in Africans. In conclusion, including family history improves the accuracy of polygenic risk scores, particularly in diverse populations.

Mendelian randomization for causal inference accounting for pleiotropy and sample structure using genome-wide summary statistics

Mendelian randomization (MR) is a valuable tool for inferring causal relationships among a wide range of traits using summary statistics from genome-wide association studies (GWASs). Existing summary-level MR methods often rely on strong assumptions, resulting in many false-positive findings. To relax MR assumptions, ongoing research has been primarily focused on accounting for confounding due to pleiotropy. Here, we show that sample structure is another major confounding factor, including population stratification, cryptic relatedness, and sample overlap. We propose a unified MR approach, MR-APSS, which 1) accounts for pleiotropy and sample structure simultaneously by leveraging genome-wide information; and 2) allows the inclusion of more genetic variants with moderate effects as instrument variables (IVs) to improve statistical power without inflating type I errors. We first evaluated MR-APSS using comprehensive simulations and negative controls and then applied MR-APSS to study the causal relationships among a collection of diverse complex traits. The results suggest that MR-APSS can better identify plausible causal relationships with high reliability. In particular, MR-APSS can perform well for highly polygenic traits, where the IV strengths tend to be relatively weak and existing summary-level MR methods for causal inference are vulnerable to confounding effects.

transferGWAS: GWAS of images using deep transfer learning

MOTIVATION

Medical images can provide rich information about diseases and their biology. However, investigating their association with genetic variation requires non-standard methods. We propose transferGWAS, a novel approach to perform genome-wide association studies directly on full medical images. First, we learn semantically meaningful representations of the images based on a transfer learning task, during which a deep neural network is trained on independent but similar data. Then, we perform genetic association tests with these representations.

RESULTS

We validate the type I error rates and power of transferGWAS in simulation studies of synthetic images. Then we apply transferGWAS in a genome-wide association study of retinal fundus images from the UK Biobank. This first-of-a-kind GWAS of full imaging data yielded 60 genomic regions associated with retinal fundus images, of which 7 are novel candidate loci for eye-related traits and diseases.

AVAILABILITY AND IMPLEMENTATION

Our method is implemented in Python and available at https://github.com/mkirchler/transferGWAS/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Leveraging the local genetic structure for trans-ancestry association mapping

Over the past two decades, genome-wide association studies (GWASs) have successfully advanced our understanding of the genetic basis of complex traits. Despite the fruitful discovery of GWASs, most GWAS samples are collected from European populations, and these GWASs are often criticized for their lack of ancestry diversity. Trans-ancestry association mapping (TRAM) offers an exciting opportunity to fill the gap of disparities in genetic studies between non-Europeans and Europeans. Here, we propose a statistical method, LOG-TRAM, to leverage the local genetic architecture for TRAM. By using biobank-scale datasets, we showed that LOG-TRAM can greatly improve the statistical power of identifying risk variants in under-represented populations while producing well-calibrated p values. We applied LOG-TRAM to the GWAS summary statistics of various complex traits/diseases from BioBank Japan, UK Biobank, and African populations. We obtained substantial gains in power and achieved effective correction of confounding biases in TRAM. Finally, we showed that LOG-TRAM can be successfully applied to identify ancestry-specific loci and the LOG-TRAM output can be further used for construction of more accurate polygenic risk scores in under-represented populations.

The sequences of 150,119 genomes in the UK Biobank

Detailed knowledge of how diversity in the sequence of the human genome affects phenotypic diversity depends on a comprehensive and reliable characterization of both sequences and phenotypic variation. Over the past decade, insights into this relationship have been obtained from whole-exome sequencing or whole-genome sequencing of large cohorts with rich phenotypic data^1,2. Here we describe the analysis of whole-genome sequencing of 150,119 individuals from the UK Biobank³. This constitutes a set of high-quality variants, including 585,040,410 single-nucleotide polymorphisms, representing 7.0% of all possible human single-nucleotide polymorphisms, and 58,707,036 indels. This large set of variants allows us to characterize selection based on sequence variation within a population through a depletion rank score of windows along the genome. Depletion rank analysis shows that coding exons represent a small fraction of regions in the genome subject to strong sequence conservation. We define three cohorts within the UK Biobank: a large British Irish cohort, a smaller African cohort and a South Asian cohort. A haplotype reference panel is provided that allows reliable imputation of most variants carried by three or more sequenced individuals. We identified 895,055 structural variants and 2,536,688 microsatellites, groups of variants typically excluded from large-scale whole-genome sequencing studies. Using this formidable new resource, we provide several examples of trait associations for rare variants with large effects not found previously through studies based on whole-exome sequencing and/or imputation.

Unbiased immune profiling reveals a natural killer cell-peripheral nerve axis in fibromyalgia

ABSTRACT

The pathophysiology of fibromyalgia syndrome (FMS) remains elusive, leading to a lack of objective diagnostic criteria and targeted treatment. We globally evaluated immune system changes in FMS by conducting multiparametric flow cytometry analyses of peripheral blood mononuclear cells and identified a natural killer (NK) cell decrease in patients with FMS. Circulating NK cells in FMS were exhausted yet activated, evidenced by lower surface expression of CD16, CD96, and CD226 and more CD107a and TIGIT. These NK cells were hyperresponsive, with increased CCL4 production and expression of CD107a when co-cultured with human leukocyte antigen null target cells. Genetic and transcriptomic pathway analyses identified significant enrichment of cell activation pathways in FMS driven by NK cells. Skin biopsies showed increased expression of NK activation ligand, unique long 16-binding protein, on subepidermal nerves of patients FMS and the presence of NK cells near peripheral nerves. Collectively, our results suggest that chronic activation and redistribution of circulating NK cells to the peripheral nerves contribute to the immunopathology associated with FMS.

Computational estimation of quality and clinical relevance of cancer cell lines

Immortal cancer cell lines (CCLs) are the most widely used system for investigating cancer biology and for the preclinical development of oncology therapies. Pharmacogenomic and genome‐wide editing screenings have facilitated the discovery of clinically relevant gene–drug interactions and novel therapeutic targets via large panels of extensively characterised CCLs. However, tailoring pharmacological strategies in a precision medicine context requires bridging the existing gaps between tumours and in vitro models. Indeed, intrinsic limitations of CCLs such as misidentification, the absence of tumour microenvironment and genetic drift have highlighted the need to identify the most faithful CCLs for each primary tumour while addressing their heterogeneity, with the development of new models where necessary. Here, we discuss the most significant limitations of CCLs in representing patient features, and we review computational methods aiming at systematically evaluating the suitability of CCLs as tumour proxies and identifying the best patient representative in vitro models. Additionally, we provide an overview of the applications of these methods to more complex models and discuss future machine‐learning‐based directions that could resolve some of the arising discrepancies.

Inherited basis of visceral, abdominal subcutaneous and gluteofemoral fat depots

For any given level of overall adiposity, individuals vary considerably in fat distribution. The inherited basis of fat distribution in the general population is not fully understood. Here, we study up to 38,965 UK Biobank participants with MRI-derived visceral (VAT), abdominal subcutaneous (ASAT), and gluteofemoral (GFAT) adipose tissue volumes. Because these fat depot volumes are highly correlated with BMI, we additionally study six local adiposity traits: VAT adjusted for BMI and height (VATadj), ASATadj, GFATadj, VAT/ASAT, VAT/GFAT, and ASAT/GFAT. We identify 250 independent common variants (39 newly-identified) associated with at least one trait, with many associations more pronounced in female participants. Rare variant association studies extend prior evidence for PDE3B as an important modulator of fat distribution. Local adiposity traits (1) highlight depot-specific genetic architecture and (2) enable construction of depot-specific polygenic scores that have divergent associations with type 2 diabetes and coronary artery disease. These results - using MRI-derived, BMI-independent measures of local adiposity - confirm fat distribution as a highly heritable trait with important implications for cardiometabolic health outcomes.

Impact of linkage disequilibrium heterogeneity along the genome on genomic prediction and heritability estimation

BACKGROUND

Compared to medium-density single nucleotide polymorphism (SNP) data, high-density SNP data contain abundant genetic variants and provide more information for the genetic evaluation of livestock, but it has been shown that they do not confer any advantage for genomic prediction and heritability estimation. One possible reason is the uneven distribution of the linkage disequilibrium (LD) along the genome, i.e., LD heterogeneity among regions. The aim of this study was to effectively use genome-wide SNP data for genomic prediction and heritability estimation by using models that control LD heterogeneity among regions.

METHODS

The LD-adjusted kinship (LDAK) and LD-stratified multicomponent (LDS) models were used to control LD heterogeneity among regions and were compared with the classical model that has no such control. Simulated and real traits of 2000 dairy cattle individuals with imputed high-density (770K) SNP data were used. Five types of phenotypes were simulated, which were controlled by very strongly, strongly, moderately, weakly and very weakly tagged causal variants, respectively. The performances of the models with high- and medium-density (50K) panels were compared to verify that the models that controlled LD heterogeneity among regions were more effective with high-density data.

RESULTS

Compared to the medium-density panel, the use of the high-density panel did not improve and even decreased prediction accuracies and heritability estimates from the classical model for both simulated and real traits. Compared to the classical model, LDS effectively improved the accuracy of genomic predictions and unbiasedness of heritability estimates, regardless of the genetic architecture of the trait. LDAK applies only to traits that are mainly controlled by weakly tagged causal variants, but is still less effective than LDS for this type of trait. Compared with the classical model, LDS improved prediction accuracy by about 13% for simulated phenotypes and by 0.3 to ~ 10.7% for real traits with the high-density panel, and by ~ 1% for simulated phenotypes and by - 0.1 to ~ 6.9% for real traits with the medium-density panel.

CONCLUSIONS

Grouping SNPs based on regional LD to construct the LD-stratified multicomponent model can effectively eliminate the adverse effects of LD heterogeneity among regions, and greatly improve the efficiency of high-density SNP data for genomic prediction and heritability estimation.

PolarMorphism enables discovery of shared genetic variants across multiple traits from GWAS summary statistics

MOTIVATION

Pleiotropic SNPs are associated with multiple traits. Such SNPs can help pinpoint biological processes with an effect on multiple traits or point to a shared etiology between traits. We present PolarMorphism, a new method for the identification of pleiotropic SNPs from genome-wide association studies (GWAS) summary statistics. PolarMorphism can be readily applied to more than two traits or whole trait domains. PolarMorphism makes use of the fact that trait-specific SNP effect sizes can be seen as Cartesian coordinates and can thus be converted to polar coordinates r (distance from the origin) and theta (angle with the Cartesian x-axis, in the case of two traits). r describes the overall effect of a SNP, while theta describes the extent to which a SNP is shared. r and theta are used to determine the significance of SNP sharedness, resulting in a P-value per SNP that can be used for further analysis.

RESULTS

We apply PolarMorphism to a large collection of publicly available GWAS summary statistics enabling the construction of a pleiotropy network that shows the extent to which traits share SNPs. We show how PolarMorphism can be used to gain insight into relationships between traits and trait domains and contrast it with genetic correlation. Furthermore, pathway analysis of the newly discovered pleiotropic SNPs demonstrates that analysis of more than two traits simultaneously yields more biologically relevant results than the combined results of pairwise analysis of the same traits. Finally, we show that PolarMorphism is more efficient and more powerful than previously published methods.

AVAILABILITY AND IMPLEMENTATION

code: https://github.com/UMCUGenetics/PolarMorphism, results: 10.5281/zenodo.5844193.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome-wide meta-analysis of iron status biomarkers and the effect of iron on all-cause mortality in HUNT

Iron is essential for many biological processes, but iron levels must be tightly regulated to avoid harmful effects of both iron deficiency and overload. Here, we perform genome-wide association studies on four iron-related biomarkers (serum iron, serum ferritin, transferrin saturation, total iron-binding capacity) in the Trøndelag Health Study (HUNT), the Michigan Genomics Initiative (MGI), and the SardiNIA study, followed by their meta-analysis with publicly available summary statistics, analyzing up to 257,953 individuals. We identify 123 genetic loci associated with iron traits. Among 19 novel protein-altering variants, we observe a rare missense variant (rs367731784) in HUNT, which suggests a role for DNAJC13 in transferrin recycling. We further validate recently published results using genetic risk scores for each biomarker in HUNT (6% variance in serum iron explained) and present linear and non-linear Mendelian randomization analyses of the traits on all-cause mortality. We find evidence of a harmful effect of increased serum iron and transferrin saturation in linear analyses that estimate population-averaged effects. However, there was weak evidence of a protective effect of increasing serum iron at the very low end of its distribution. Our findings contribute to our understanding of the genes affecting iron status and its consequences on human health.

A GWAS on the four major iron-related biomarkers is conducted using data from the Trøndelag Health Study (HUNT), Michigan Genomics Initiative (MGI), and SardiNIA study, to identify 123 loci associated with iron homeostasis and effect on all-cause mortality.