Eliciting priors and relaxing the single causal variant assumption in colocalisation analyses

Integrative common and rare variant analyses provide insights into the genetic architecture of liver cirrhosis

We report a multi-ancestry genome-wide association study on liver cirrhosis and its associated endophenotypes, alanine aminotransferase (ALT) and γ-glutamyl transferase. Using data from 12 cohorts, including 18,265 cases with cirrhosis, 1,782,047 controls, up to 1 million individuals with liver function tests and a validation cohort of 21,689 cases and 617,729 controls, we identify and validate 14 risk associations for cirrhosis. Many variants are located near genes involved in hepatic lipid metabolism. One of these, PNPLA3 p.Ile148Met, interacts with alcohol intake, obesity and diabetes on the risk of cirrhosis and hepatocellular carcinoma (HCC). We develop a polygenic risk score that associates with the progression from cirrhosis to HCC. By focusing on prioritized genes from common variant analyses, we find that rare coding variants in GPAM associate with lower ALT, supporting GPAM as a potential target for therapeutic inhibition. In conclusion, this study provides insights into the genetic underpinnings of cirrhosis.

Large-scale cross-ancestry genome-wide meta-analysis of serum urate

Hyperuricemia is an essential causal risk factor for gout and is associated with cardiometabolic diseases. Given the limited contribution of East Asian ancestry to genome-wide association studies of serum urate, the genetic architecture of serum urate requires exploration. A large-scale cross-ancestry genome-wide association meta-analysis of 1,029,323 individuals and ancestry-specific meta-analysis identifies a total of 351 loci, including 17 previously unreported loci. The genetic architecture of serum urate control is similar between European and East Asian populations. A transcriptome-wide association study, enrichment analysis, and colocalization analysis in relevant tissues identify candidate serum urate-associated genes, including CTBP1, SKIV2L, and WWP2. A phenome-wide association study using polygenic risk scores identifies serum urate-correlated diseases including heart failure and hypertension. Mendelian randomization and mediation analyses show that serum urate-associated genes might have a causal relationship with serum urate-correlated diseases via mediation effects. This study elucidates our understanding of the genetic architecture of serum urate control.

The Impact of Inherited Genetic Variation on DNA Methylation in Prostate Cancer and Benign Tissues of African American and European American Men

BACKGROUND

American men of African ancestry (AA) have higher prostate cancer incidence and mortality rates compared with American men of European ancestry (EA). Differences in genetic susceptibility mechanisms may contribute to this disparity.

METHODS

To gain insights into the regulatory mechanisms of prostate cancer susceptibility variants, we tested the association between SNPs and DNA methylation (DNAm) at nearby CpG sites across the genome in benign and cancer prostate tissue from 74 AA and 74 EA men. Genome-wide SNP data (from benign tissue) and DNAm were generated using Illumina arrays.

RESULTS

Among AA men, we identified 6,298 and 2,641 cis-methylation QTLs (meQTL; FDR of 0.05) in benign and tumor tissue, respectively, with 6,960 and 1,700 detected in EA men. We leveraged genome-wide association study (GWAS) summary statistics to identify previously reported prostate cancer GWAS signals likely to share a common causal variant with a detected meQTL. We identified nine GWAS-meQTL pairs with strong evidence of colocalization (four in EA benign, three in EA tumor, two in AA benign, and three in AA tumor). Among these colocalized GWAS-meQTL pairs, we identified colocalizing expression quantitative trait loci (eQTL) impacting four eGenes with known roles in tumorigenesis.

CONCLUSIONS

These findings highlight epigenetic regulatory mechanisms by which prostate cancer-risk SNPs can modify local DNAm and/or gene expression in prostate tissue.

IMPACT

Overall, our findings showed general consistency in the meQTL landscape of AA and EA men, but meQTLs often differ by tissue type (normal vs. cancer). Ancestry-based linkage disequilibrium differences and lack of AA representation in GWAS decrease statistical power to detect colocalization for some regions.

Genetic determinants of plasma protein levels in the Estonian population

The proteome holds great potential as an intermediate layer between the genome and phenome. Previous protein quantitative trait locus studies have focused mainly on describing the effects of common genetic variations on the proteome. Here, we assessed the impact of the common and rare genetic variations as well as the copy number variants (CNVs) on 326 plasma proteins measured in up to 500 individuals. We identified 184 cis and 94 trans signals for 157 protein traits, which were further fine-mapped to credible sets for 101 cis and 87 trans signals for 151 proteins. Rare genetic variation contributed to the levels of 7 proteins, with 5 cis and 14 trans associations. CNVs were associated with the levels of 11 proteins (7 cis and 5 trans), examples including a 3q12.1 deletion acting as a hub for multiple trans associations; and a CNV overlapping NAIP, a sensor component of the NAIP-NLRC4 inflammasome which is affecting pro-inflammatory cytokine interleukin 18 levels. In summary, this work presents a comprehensive resource of genetic variation affecting the plasma protein levels and provides the interpretation of identified effects.

Secreted Protein Profiling of Human Aortic Smooth Muscle Cells Identifies Vascular Disease Associations

BACKGROUND

Smooth muscle cells (SMCs), which make up the medial layer of arteries, are key cell types involved in cardiovascular disease, the leading cause of mortality and morbidity worldwide. In response to microenvironment alterations, SMCs dedifferentiate from a contractile to a synthetic phenotype characterized by an increased proliferation, migration, production of ECM (extracellular matrix) components, and decreased expression of SMC-specific contractile markers. These phenotypic changes result in vascular remodeling and contribute to the pathogenesis of cardiovascular disease, including coronary artery disease, stroke, hypertension, and aortic aneurysms. Here, we aim to identify the genetic variants that regulate ECM secretion in SMCs and predict the causal proteins associated with vascular disease-related loci identified in genome-wide association studies.

METHODS

Using human aortic SMCs from 123 multiancestry healthy heart transplant donors, we collected the serum-free media in which the cells were cultured for 24 hours and conducted liquid chromatography-tandem mass spectrometry-based proteomic analysis of the conditioned media.

RESULTS

We measured the abundance of 270 ECM and related proteins. Next, we performed protein quantitative trait locus mapping and identified 20 loci associated with secreted protein abundance in SMCs. We functionally annotated these loci using a colocalization approach. This approach prioritized the genetic variant rs6739323-A at the 2p22.3 locus, which is associated with lower expression of LTBP1 (latent-transforming growth factor beta-binding protein 1) in SMCs and atherosclerosis-prone areas of the aorta, and increased risk for SMC calcification. We found that LTBP1 expression is abundant in SMCs, and its expression at mRNA and protein levels was reduced in unstable and advanced atherosclerotic plaque lesions.

CONCLUSIONS

Our results unravel the SMC proteome signature associated with vascular disorders, which may help identify potential therapeutic targets to accelerate the pathway to translation.

Integrative multi-omics analysis identifies genetically supported druggable targets and immune cell specificity for myasthenia gravis

BACKGROUND

Myasthenia gravis (MG) is a chronic autoimmune disorder characterized by fluctuating muscle weakness. Despite the availability of established therapies, the management of MG symptoms remains suboptimal, partially attributed to lack of efficacy or intolerable side-effects. Therefore, new effective drugs are warranted for treatment of MG.

METHODS

By employing an analytical framework that combines Mendelian randomization (MR) and colocalization analysis, we estimate the causal effects of blood druggable expression quantitative trait loci (eQTLs) and protein quantitative trait loci (pQTLs) on the susceptibility of MG. We subsequently investigated whether potential genetic effects exhibit cell-type specificity by utilizing genetic colocalization analysis to assess the interplay between immune-cell-specific eQTLs and MG risk.

RESULTS

We identified significant MR results for four genes (CDC42BPB, CD226, PRSS36, and TNFSF12) using cis-eQTL genetic instruments and three proteins (CTSH, PRSS8, and CPN2) using cis-pQTL genetic instruments. Six of these loci demonstrated evidence of colocalization with MG susceptibility (posterior probability > 0.80). We next undertook genetic colocalization to investigate cell-type-specific effects at these loci. Notably, we identified robust evidence of colocalization, with a posterior probability of 0.854, linking CTSH expression in TH2 cells and MG risk.

CONCLUSIONS

This study provides crucial insights into the genetic and molecular factors associated with MG susceptibility, singling out CTSH as a potential candidate for in-depth investigation and clinical consideration. It additionally sheds light on the immune-cell regulatory mechanisms related to the disease. However, further research is imperative to validate these targets and evaluate their feasibility for drug development.

Genome-wide analyses identify 21 infertility loci and over 400 reproductive hormone loci across the allele frequency spectrum

Genome-wide association studies (GWASs) may help inform treatments for infertility, whose causes remain unknown in many cases. Here we present GWAS meta-analyses across six cohorts for male and female infertility in up to 41,200 cases and 687,005 controls. We identified 21 genetic risk loci for infertility (P≤5E-08), of which 12 have not been reported for any reproductive condition. We found positive genetic correlations between endometriosis and all-cause female infertility (r g=0.585, P=8.98E-14), and between polycystic ovary syndrome and anovulatory infertility (r g=0.403, P=2.16E-03). The evolutionary persistence of female infertility-risk alleles in EBAG9 may be explained by recent directional selection. We additionally identified up to 269 genetic loci associated with follicle-stimulating hormone (FSH), luteinising hormone, oestradiol, and testosterone through sex-specific GWAS meta-analyses (N=6,095-246,862). While hormone-associated variants near FSHB and ARL14EP colocalised with signals for anovulatory infertility, we found no r g between female infertility and reproductive hormones (P>0.05). Exome sequencing analyses in the UK Biobank (N=197,340) revealed that women carrying testosterone-lowering rare variants in GPC2 were at higher risk of infertility (OR=2.63, P=1.25E-03). Taken together, our results suggest that while individual genes associated with hormone regulation may be relevant for fertility, there is limited genetic evidence for correlation between reproductive hormones and infertility at the population level. We provide the first comprehensive view of the genetic architecture of infertility across multiple diagnostic criteria in men and women, and characterise its relationship to other health conditions.

Genetic overlap and causality between COVID-19 and multi-site chronic pain: the importance of immunity

Background

The existence of chronic pain increases susceptibility to virus and is now widely acknowledged as a prominent feature recognized as a major manifestation of long-term coronavirus disease 2019 (COVID-19) infection. Given the ongoing COVID-19 pandemic, it is imperative to explore the genetic associations between chronic pain and predisposition to COVID-19.

Methods

We conducted genetic analysis at the single nucleotide polymorphism (SNP), gene, and molecular levels using summary statistics of genome-wide association study (GWAS) and analyzed the drug targets by summary data-based Mendelian randomization analysis (SMR) to alleviate the multi-site chronic pain in COVID-19. Additionally, we performed a latent causal variable (LCV) method to investigate the causal relationship between chronic pain and susceptibility to COVID-19.

Results

The cross-trait meta-analysis identified 19 significant SNPs shared between COVID-19 and chronic pain. Coloc analysis indicated that the posterior probability of association (PPH4) for three loci was above 70% in both critical COVID-19 and COVID-19, with the corresponding top three SNPs being rs13135092, rs7588831, and rs13135092. A total of 482 significant overlapped genes were detected from MAGMA and CPASSOC results. Additionally, the gene ANAPC4 was identified as a potential drug target for treating chronic pain (P=7.66E-05) in COVID-19 (P=8.23E-03). Tissue enrichment analysis highlighted that the amygdala (P=7.81E-04) and prefrontal cortex (P=8.19E-05) as pivotal in regulating chronic pain of critical COVID-19. KEGG pathway enrichment further revealed the enrichment of pleiotropic genes in both COVID-19 (P=3.20E-03,Padjust=4.77E-02,hsa05171) and neurotrophic pathways (P=9.03E-04,Padjust =2.55E-02,hsa04621). Finally, the latent causal variable (LCV) model was applied to find the genetic component of critical COVID-19 was causal for multi-site chronic pain (P=0.015), with a genetic causality proportion (GCP) of was 0.60.

Conclusions

In this study, we identified several functional genes and underscored the pivotal role of the inflammatory system in the correlation between the paired traits. Notably, heat shock proteins emerged as potential objective biomarkers for chronic pain symptoms in individuals with COVID-19. Additionally, the ubiquitin system might play a role in mediating the impact of COVID-19 on chronic pain. These findings contribute to a more comprehensive understanding of the pleiotropy between COVID-19 and chronic pain, offering insights for therapeutic trials.

Genome-wide association study identifies 30 obsessive-compulsive disorder associated loci

Obsessive-compulsive disorder (OCD) affects ~1% of the population and exhibits a high SNP-heritability, yet previous genome-wide association studies (GWAS) have provided limited information on the genetic etiology and underlying biological mechanisms of the disorder. We conducted a GWAS meta-analysis combining 53,660 OCD cases and 2,044,417 controls from 28 European-ancestry cohorts revealing 30 independent genome-wide significant SNPs and a SNP-based heritability of 6.7%. Separate GWAS for clinical, biobank, comorbid, and self-report sub-groups found no evidence of sample ascertainment impacting our results. Functional and positional QTL gene-based approaches identified 249 significant candidate risk genes for OCD, of which 25 were identified as putatively causal, highlighting WDR6, DALRD3, CTNND1 and genes in the MHC region. Tissue and single-cell enrichment analyses highlighted hippocampal and cortical excitatory neurons, along with D1- and D2-type dopamine receptor-containing medium spiny neurons, as playing a role in OCD risk. OCD displayed significant genetic correlations with 65 out of 112 examined phenotypes. Notably, it showed positive genetic correlations with all included psychiatric phenotypes, in particular anxiety, depression, anorexia nervosa, and Tourette syndrome, and negative correlations with a subset of the included autoimmune disorders, educational attainment, and body mass index.. This study marks a significant step toward unraveling its genetic landscape and advances understanding of OCD genetics, providing a foundation for future interventions to address this debilitating disorder.

Multivariate canonical correlation analysis identifies additional genetic variants for chronic kidney disease

Chronic kidney diseases (CKD) have genetic associations with kidney function. Univariate genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) associated with estimated glomerular filtration rate (eGFR) and blood urea nitrogen (BUN), two complementary kidney function markers. However, it is unknown whether additional SNPs for kidney function can be identified by multivariate statistical analysis. To address this, we applied canonical correlation analysis (CCA), a multivariate method, to two individual-level CKD genotype datasets, and metaCCA to two published GWAS summary statistics datasets. We identified SNPs previously associated with kidney function by published univariate GWASs with high replication rates, validating the metaCCA method. We then extended discovery and identified previously unreported lead SNPs for both kidney function markers, jointly. These showed expression quantitative trait loci (eQTL) colocalisation with genes having significant differential expression between CKD and healthy individuals. Several of these identified lead missense SNPs were predicted to have a functional impact, including in SLC14A2. We also identified previously unreported lead SNPs that showed significant correlation with both kidney function markers, jointly, in the European ancestry CKDGen, National Unified Renal Translational Research Enterprise (NURTuRE)-CKD and Salford Kidney Study (SKS) datasets. Of these, rs3094060 colocalised with FLOT1 gene expression and was significantly more common in CKD cases in both NURTURE-CKD and SKS, than in the general population. Overall, by using multivariate analysis by CCA, we identified additional SNPs and genes for both kidney function and CKD, that can be prioritised for further CKD analyses.

Epigenetic variation impacts individual differences in the transcriptional response to influenza infection

Humans display remarkable interindividual variation in their immune response to identical challenges. Yet, our understanding of the genetic and epigenetic factors contributing to such variation remains limited. Here we performed in-depth genetic, epigenetic and transcriptional profiling on primary macrophages derived from individuals of European and African ancestry before and after infection with influenza A virus. We show that baseline epigenetic profiles are strongly predictive of the transcriptional response to influenza A virus across individuals. Quantitative trait locus (QTL) mapping revealed highly coordinated genetic effects on gene regulation, with many cis-acting genetic variants impacting concomitantly gene expression and multiple epigenetic marks. These data reveal that ancestry-associated differences in the epigenetic landscape can be genetically controlled, even more than gene expression. Lastly, among QTL variants that colocalized with immune-disease loci, only 7% were gene expression QTL, while the remaining genetic variants impact epigenetic marks, stressing the importance of considering molecular phenotypes beyond gene expression in disease-focused studies.

A distinct class of pan-cancer susceptibility genes revealed by an alternative polyadenylation transcriptome-wide association study

Alternative polyadenylation plays an important role in cancer initiation and progression; however, current transcriptome-wide association studies mostly ignore alternative polyadenylation when identifying putative cancer susceptibility genes. Here, we perform a pan-cancer 3' untranslated region alternative polyadenylation transcriptome-wide association analysis by integrating 55 well-powered (n > 50,000) genome-wide association studies datasets across 22 major cancer types with alternative polyadenylation quantification from 23,955 RNA sequencing samples across 7,574 individuals. We find that genetic variants associated with alternative polyadenylation are co-localized with 28.57% of cancer loci and contribute a significant portion of cancer heritability. We further identify 642 significant cancer susceptibility genes predicted to modulate cancer risk via alternative polyadenylation, 62.46% of which have been overlooked by traditional expression- and splicing- studies. As proof of principle validation, we show that alternative alleles facilitate 3' untranslated region lengthening of CRLS1 gene leading to increased protein abundance and promoted proliferation of breast cancer cells. Together, our study highlights the significant role of alternative polyadenylation in discovering new cancer susceptibility genes and provides a strong foundational framework for enhancing our understanding of the etiology underlying human cancers.

Novel insights into causal effects of serum lipids and lipid-modifying targets on cholelithiasis

OBJECTIVE

Different serum lipids and lipid-modifying targets should affect the risk of cholelithiasis differently, however, whether such effects are causal is still controversial and we aimed to answer this question.

DESIGN

We prospectively estimated the associations of four serum lipids with cholelithiasis in UK Biobank using the Cox proportional hazard model, including total cholesterol (TC), low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C) and triglycerides (TG). Furthermore, we estimated the causal associations of the genetically predicted serum lipids with cholelithiasis in Europeans using the Mendelian randomisation (MR) design. Finally, both drug-target MR and colocalisation analyses were performed to estimate the lipid-modifying targets' effects on cholelithiasis, including HMGCR, NPC1L1, PCSK9, APOB, LDLR, ACLY, ANGPTL3, MTTP, PPARA, PPARD and PPARG.

RESULTS

We found that serum levels of LDL-C and HDL-C were inversely associated with cholelithiasis risk and such associations were linear. However, the serum level of TC was non-linearly associated with cholelithiasis risk where lower TC was associated with higher risk of cholelithiasis, and the serum TG should be in an inverted 'U-shaped' relationship with it. The MR analyses supported that lower TC and higher TG levels were two independent causal risk factors. The drug-target MR analysis suggested that HMGCR inhibition should reduce the risk of cholelithiasis, which was corroborated by colocalisation analysis.

CONCLUSION

Lower serum TC can causally increase the risk of cholelithiasis. The cholelithiasis risk would increase with the elevation of serum TG but would decrease when exceeding 2.57 mmol/L. The use of HMGCR inhibitors should prevent its risk.

Identification of the shared genes in type 2 diabetes mellitus and osteoarthritis and the role of quercetin

This study investigated the underlying comorbidity mechanism between type 2 diabetes mellitus (T2DM) and osteoarthritis (OA), while also assessing the therapeutic potential of quercetin for early intervention and treatment of these two diseases. The shared genes were obtained through GEO2R, limma and weighted gene co-expression network analysis (WGCNA), and validated using clinical databases and the area under the curves (ROC). Functional enrichment analysis was conducted to elucidate the underlying mechanisms of comorbidity between T2DM and OA. The infiltration of immune cells was analysed using the CIBERSORT algorithm in conjunction with ESTIMATE algorithm. Subsequently, transcriptional regulation analysis, potential chemical prediction, gene-disease association, relationships between the shared genes and ferroptosis as well as immunity-related genes were investigated along with molecular docking. We identified the 12 shared genes (EPHA3, RASIP1, PENK, LRRC17, CEBPB, EFEMP2, UBAP1, PPP1R15A, SPEN, MAFF, GADD45B and KLF4) across the four datasets. Our predictions suggested that targeting these shared genes could potentially serve as therapeutic interventions for both T2DM and OA. Specifically, they are involved in key signalling pathways such as p53, IL-17, NF-kB and MAPK signalling pathways. Furthermore, the regulation of ferroptosis and immunity appears to be interconnected in both diseases. Notably, in this context quercetin emerges as a promising drug candidate for treating T2DM and OA by specifically targeting the shared genes. We conducted a bioinformatics analysis to identify potential therapeutic targets, mechanisms and drugs for T2DM and OA, thereby offering novel insights into molecular therapy for these two diseases.

Investigating the shared genetic architecture between COVID-19 and obesity: a large-scale genome wide cross-trait analysis

Observational studies have reported high comorbidity between obesity and severe COVID-19. The aim of this study is to explore whether genetic factors are involved in the co-occurrence of the two traits. Based on the available genome-wide association studies (GWAS) summary statistics, we explored the genetic correlation and performed cross-trait meta-analysis (CPASSOC) and colocalization analysis (COLOC) to detect pleiotropic single nucleotide polymorphisms (SNPs). At the genetic level, we obtained genes detected by Functional mapping and annotation (FUMA) and the Multi-marker Analysis of GenoMic Annotation (MAGMA). Potential functional genes were further investigated by summary-data-based Mendelian randomization (SMR). Finally, the casualty was identiied using the latent causal variable model (LCV). A significant positive genetic correlation was revealed between obesity and COVID-19. We found 331 shared genetic SNPs by CPASSOC and 13 shared risk loci by COLOC. At the genetic level, We obtained 3546 pleiotropic genes, among which 107 genes were found to be significantly expressed by SMR. Lastly, we observed these genes were mainly enriched in immune pathways and signaling transduction. These indings could provide new insights into the etiology of comorbidity and have implications for future therapeutic trial.

Integrating genetic regulation and single-cell expression with GWAS prioritizes causal genes and cell types for glaucoma

Primary open-angle glaucoma (POAG), characterized by retinal ganglion cell death, is a leading cause of irreversible blindness worldwide. However, its molecular and cellular causes are not well understood. Elevated intraocular pressure (IOP) is a major risk factor, but many patients have normal IOP. Colocalization and Mendelian randomization analysis of >240 POAG and IOP genome-wide association study (GWAS) loci and overlapping expression and splicing quantitative trait loci (e/sQTLs) in 49 GTEx tissues and retina prioritizes causal genes for 60% of loci. These genes are enriched in pathways implicated in extracellular matrix organization, cell adhesion, and vascular development. Analysis of single-nucleus RNA-seq of glaucoma-relevant eye tissues reveals that the POAG and IOP colocalizing genes and genome-wide associations are enriched in specific cell types in the aqueous outflow pathways, retina, optic nerve head, peripapillary sclera, and choroid. This study nominates IOP-dependent and independent regulatory mechanisms, genes, and cell types that may contribute to POAG pathogenesis.

Cerebrospinal Fluid C1-Esterase Inhibitor and Tie-1 Levels Affect Cognitive Performance: Evidence from Proteome-Wide Mendelian Randomization

OBJECTIVE

The association of cerebrospinal fluid (CSF) protein levels with cognitive function in the general population remains largely unexplored. We performed Mendelian randomization (MR) analyses to query which CSF proteins may have potential causal effects on cognitive performance.

METHODS AND ANALYSIS

Genetic associations with CSF proteins were obtained from a genome-wide association study conducted in up to 835 European-ancestry individuals and for cognitive performance from a meta-analysis of GWAS including 257,841 European-ancestry individuals. We performed Mendelian randomization (MR) analyses to test the effect of randomly allocated variation in 154 genetically predicted CSF protein levels on cognitive performance. Findings were validated by performing colocalization analyses and considering cognition-related phenotypes.

RESULTS

Genetically predicted C1-esterase inhibitor levels in the CSF were associated with a better cognitive performance (SD units of cognitive performance per 1 log-relative fluorescence unit (RFU): 0.23, 95% confidence interval: 0.12 to 0.35, p = 7.91 × 10-5), while tyrosine-protein kinase receptor Tie-1 (sTie-1) levels were associated with a worse cognitive performance (-0.43, -0.62 to -0.23, p = 2.08 × 10-5). These findings were supported by colocalization analyses and by concordant effects on distinct cognition-related and brain-volume measures.

CONCLUSIONS

Human genetics supports a role for the C1-esterase inhibitor and sTie-1 in cognitive performance.

Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects

Bulk-tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, and context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from the blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell-type proportions, we demonstrate that cell-type iQTLs could be considered as proxies for cell-type-specific QTL effects, particularly for the most abundant cell type in the tissue. The interpretation of age iQTLs, however, warrants caution because the moderation effect of age on the genotype and molecular phenotype association could be mediated by changes in cell-type composition. Finally, we show that cell-type iQTLs contribute to cell-type-specific enrichment of diseases that, in combination with additional functional data, could guide future functional studies. Overall, this study highlights the use of iQTLs to gain insights into the context specificity of regulatory effects.

DNA methylation at birth and lateral ventricular volume in childhood: a neuroimaging epigenetics study

BACKGROUND

Lateral ventricular volume (LVV) enlargement has been repeatedly linked to schizophrenia; yet, what biological factors shape LVV during early development remain unclear. DNA methylation (DNAm), an essential process for neurodevelopment that is altered in schizophrenia, is a key molecular system of interest.

METHODS

In this study, we conducted the first epigenome-wide association study of neonatal DNAm in cord blood with LVV in childhood (measured using T1-weighted brain scans at 10 years), based on data from a large population-based birth cohort, the Generation R Study (N = 840). Employing both probe-level and methylation profile score (MPS) approaches, we further examined whether epigenetic modifications identified at birth in cord blood are: (a) also observed cross-sectionally in childhood using peripheral blood DNAm at age of 10 years (Generation R, N = 370) and (b) prospectively associated with LVV measured in young adulthood in an all-male sample from the Avon Longitudinal Study of Parents and Children (ALSPAC, N = 114).

RESULTS

At birth, DNAm levels at four CpGs (annotated to potassium channel tetramerization domain containing 3, KCTD3; SHH signaling and ciliogenesis regulator, SDCCAG8; glutaredoxin, GLRX) prospectively associated with childhood LVV after genome-wide correction; these genes have been implicated in brain development and psychiatric traits including schizophrenia. An MPS capturing a broader epigenetic profile of LVV - but not individual top hits - showed significant cross-sectional associations with LVV in childhood in Generation R and prospectively associated with LVV in early adulthood within ALSPAC.

CONCLUSIONS

This study finds suggestive evidence that DNAm at birth prospectively associates with LVV at different life stages, albeit with small effect sizes. The prediction of MPS on LVV in a childhood sample and an independent male adult sample further underscores the stability and reproducibility of DNAm as a potential marker for LVV. Future studies with larger samples and comparable time points across development are needed to further elucidate how DNAm associates with this clinically relevant brain structure and risk for neuropsychiatric disorders, and what factors explain the identified DNAm profile of LVV at birth.

Open Science Practices in Psychiatric Genetics: A Primer

Open science ensures that research is transparently reported and freely accessible for all to assess and collaboratively build on. Psychiatric genetics has led among the health sciences in implementing some open science practices in common study designs, such as replication as part of genome-wide association studies. However, thorough open science implementation guidelines are limited and largely not specific to data, privacy, and research conduct challenges in psychiatric genetics. Here, we present a primer of open science practices, including selection of a research topic with patients/nonacademic collaborators, equitable authorship and citation practices, design of replicable, reproducible studies, preregistrations, open data, and privacy issues. We provide tips for informative figures and inclusive, precise reporting. We discuss considerations in working with nonacademic collaborators and distributing research through preprints, blogs, social media, and accessible lecture materials. Finally, we provide extra resources to support every step of the research process.

Endothelial lipid droplets suppress eNOS to link high fat consumption to blood pressure elevation

Metabolic syndrome, today affecting more than 20% of the US population, is a group of 5 conditions that often coexist and that strongly predispose to cardiovascular disease. How these conditions are linked mechanistically remains unclear, especially two of these: obesity and elevated blood pressure. Here, we show that high fat consumption in mice leads to the accumulation of lipid droplets in endothelial cells throughout the organism and that lipid droplet accumulation in endothelium suppresses endothelial nitric oxide synthase (eNOS), reduces NO production, elevates blood pressure, and accelerates atherosclerosis. Mechanistically, the accumulation of lipid droplets destabilizes eNOS mRNA and activates an endothelial inflammatory signaling cascade that suppresses eNOS and NO production. Pharmacological prevention of lipid droplet formation reverses the suppression of NO production in cell culture and in vivo and blunts blood pressure elevation in response to a high-fat diet. These results highlight lipid droplets as a critical and unappreciated component of endothelial cell biology, explain how lipids increase blood pressure acutely, and provide a mechanistic account for the epidemiological link between obesity and elevated blood pressure.

TIE1 and TEK signalling, intraocular pressure, and primary open-angle glaucoma: a Mendelian randomization study

BACKGROUND

In primary open-angle glaucoma (POAG), lowering intraocular pressure (IOP) is the only proven way of slowing vision loss. Schlemm's canal (SC) is a hybrid vascular and lymphatic vessel that mediates aqueous humour drainage from the anterior ocular chamber. Animal studies support the importance of SC endothelial angiopoietin-TEK signalling, and more recently TIE1 signalling, in maintaining normal IOP. However, human genetic support for a causal role of TIE1 and TEK signalling in lowering IOP is currently lacking.

METHODS

GWAS summary statistics were obtained for plasma soluble TIE1 (sTIE1) protein levels (N = 35,559), soluble TEK (sTEK) protein levels (N = 35,559), IOP (N = 139,555) and POAG (Ncases = 16,677, Ncontrols = 199,580). Mendelian randomization (MR) was performed to estimate the association of genetically proxied TIE1 and TEK protein levels with IOP and POAG liability. Where significant MR estimates were obtained, genetic colocalization was performed to assess the probability of a shared causal variant (PPshared) versus distinct (PPdistinct) causal variants underlying TIE1/TEK signalling and the outcome. Publicly available single-nucleus RNA-sequencing data were leveraged to investigate differential expression of TIE1 and TEK in the human ocular anterior segment.

RESULTS

Increased genetically proxied TIE1 signalling and TEK signalling associated with a reduction in IOP (- 0.21 mmHg per SD increase in sTIE1, 95% CI = - 0.09 to - 0.33 mmHg, P = 6.57 × 10-4, and - 0.14 mmHg per SD decrease in sTEK, 95% CI = - 0.03 to - 0.25 mmHg, P = 0.011), but not with POAG liability. Colocalization analysis found that the probability of a shared causal variant was greater for TIE1 and IOP than for TEK and IOP (PPshared/(PPdistinct + PPshared) = 0.98 for TIE1 and 0.30 for TEK). In the anterior segment, TIE1 and TEK were preferentially expressed in SC, lymphatic, and vascular endothelium.

CONCLUSIONS

This study provides novel human genetic support for a causal role of both TIE1 and TEK signalling in regulating IOP. Here, combined evidence from cis-MR and colocalization analyses provide stronger support for TIE1 than TEK as a potential IOP-lowering therapeutic target.

Secreted protein profiling of human aortic smooth muscle cells identifies vascular disease associations

Background

Smooth muscle cells (SMCs), which make up the medial layer of arteries, are key cell types involved in cardiovascular diseases (CVD), the leading cause of mortality and morbidity worldwide. In response to microenvironment alterations, SMCs dedifferentiate from a "contractile" to a "synthetic" phenotype characterized by an increased proliferation, migration, production of extracellular matrix (ECM) components, and decreased expression of SMC-specific contractile markers. These phenotypic changes result in vascular remodeling and contribute to the pathogenesis of CVD, including coronary artery disease (CAD), stroke, hypertension, and aortic aneurysms. Here, we aim to identify the genetic variants that regulate ECM secretion in SMCs and predict the causal proteins associated with vascular disease-related loci identified in genome-wide association studies (GWAS).

Methods

Using human aortic SMCs from 123 multi-ancestry healthy heart transplant donors, we collected the serum-free media in which the cells were cultured for 24 hours and conducted Liquid chromatography-tandem mass spectrometry (LC-MS/MS)-based proteomic analysis of the conditioned media.

Results

We measured the abundance of 270 ECM and related proteins. Next, we performed protein quantitative trait locus mapping (pQTL) and identified 20 loci associated with secreted protein abundance in SMCs. We functionally annotated these loci using a colocalization approach. This approach prioritized the genetic variant rs6739323-A at the 2p22.3 locus, which is associated with lower expression of LTBP1 in SMCs and atherosclerosis-prone areas of the aorta, and increased risk for SMC calcification. We found that LTBP1 expression is abundant in SMCs, and its expression at mRNA and protein levels was reduced in unstable and advanced atherosclerotic plaque lesions.

Conclusions

Our results unravel the SMC proteome signature associated with vascular disorders, which may help identify potential therapeutic targets to accelerate the pathway to translation.

Genome-Wide Association Study of Pericardial Fat Area in 28 161 UK Biobank Participants

BACKGROUND Pericardial adipose tissue (PAT) is the visceral adipose tissue compartment surrounding the heart. Experimental and observational research has suggested that greater PAT deposition might mediate cardiovascular disease, independent of general or subcutaneous adiposity. We characterize the genetic architecture of adiposity-adjusted PAT and identify causal associations between PAT and adverse cardiac magnetic resonance imaging measures of cardiac structure and function in 28 161 UK Biobank participants. METHODS AND RESULTS The PAT phenotype was extracted from cardiac magnetic resonance images using an automated image analysis tool previously developed and validated in this cohort. A genome-wide association study was performed with PAT area set as the phenotype, adjusting for age, sex, and other measures of obesity. Functional mapping and Bayesian colocalization were used to understand the biologic role of identified variants. Mendelian randomization analysis was used to examine potential causal links between genetically determined PAT and cardiac magnetic resonance-derived measures of left ventricular structure and function. We discovered 12 genome-wide significant variants, with 2 independent sentinel variants (rs6428792, P=4.20×10-9 and rs11992444, P=1.30×10-12) at 2 distinct genomic loci, that were mapped to 3 potentially causal genes: T-box transcription factor 15 (TBX15), tryptophanyl tRNA synthetase 2, mitochondrial (WARS2) and early B-cell factor-2 (EBF2) through functional annotation. Bayesian colocalization additionally suggested a role of RP4-712E4.1. Genetically predicted differences in adiposity-adjusted PAT were causally associated with adverse left ventricular remodeling. CONCLUSIONS This study provides insights into the genetic architecture determining differential PAT deposition, identifies causal links with left structural and functional parameters, and provides novel data about the pathophysiological importance of adiposity distribution.

Pulmonary emphysema subtypes defined by unsupervised machine learning on CT scans

BACKGROUND

Treatment and preventative advances for chronic obstructive pulmonary disease (COPD) have been slow due, in part, to limited subphenotypes. We tested if unsupervised machine learning on CT images would discover CT emphysema subtypes with distinct characteristics, prognoses and genetic associations.

METHODS

New CT emphysema subtypes were identified by unsupervised machine learning on only the texture and location of emphysematous regions on CT scans from 2853 participants in the Subpopulations and Intermediate Outcome Measures in COPD Study (SPIROMICS), a COPD case-control study, followed by data reduction. Subtypes were compared with symptoms and physiology among 2949 participants in the population-based Multi-Ethnic Study of Atherosclerosis (MESA) Lung Study and with prognosis among 6658 MESA participants. Associations with genome-wide single-nucleotide-polymorphisms were examined.

RESULTS

The algorithm discovered six reproducible (interlearner intraclass correlation coefficient, 0.91-1.00) CT emphysema subtypes. The most common subtype in SPIROMICS, the combined bronchitis-apical subtype, was associated with chronic bronchitis, accelerated lung function decline, hospitalisations, deaths, incident airflow limitation and a gene variant near DRD1, which is implicated in mucin hypersecretion (p=1.1 ×10-8). The second, the diffuse subtype was associated with lower weight, respiratory hospitalisations and deaths, and incident airflow limitation. The third was associated with age only. The fourth and fifth visually resembled combined pulmonary fibrosis emphysema and had distinct symptoms, physiology, prognosis and genetic associations. The sixth visually resembled vanishing lung syndrome.

CONCLUSION

Large-scale unsupervised machine learning on CT scans defined six reproducible, familiar CT emphysema subtypes that suggest paths to specific diagnosis and personalised therapies in COPD and pre-COPD.

Adipose tissue eQTL meta-analysis reveals the contribution of allelic heterogeneity to gene expression regulation and cardiometabolic traits

Complete characterization of the genetic effects on gene expression is needed to elucidate tissue biology and the etiology of complex traits. Here, we analyzed 2,344 subcutaneous adipose tissue samples and identified 34K conditionally distinct expression quantitative trait locus (eQTL) signals in 18K genes. Over half of eQTL genes exhibited at least two eQTL signals. Compared to primary signals, non-primary signals had lower effect sizes, lower minor allele frequencies, and less promoter enrichment; they corresponded to genes with higher heritability and higher tolerance for loss of function. Colocalization of eQTL with conditionally distinct genome-wide association study signals for 28 cardiometabolic traits identified 3,605 eQTL signals for 1,861 genes. Inclusion of non-primary eQTL signals increased colocalized signals by 46%. Among 30 genes with ≥2 pairs of colocalized signals, 21 showed a mediating gene dosage effect on the trait. Thus, expanded eQTL identification reveals more mechanisms underlying complex traits and improves understanding of the complexity of gene expression regulation.

Brain-wide genome-wide colocalization study for integrating genetics, transcriptomics and brain morphometry in Alzheimer's disease

Alzheimer's disease (AD) is one of the most common neurodegenerative diseases. However, the AD mechanism has not yet been fully elucidated to date, hindering the development of effective therapies. In our work, we perform a brain imaging genomics study to link genetics, single-cell gene expression data, tissue-specific gene expression data, brain imaging-derived volumetric endophenotypes, and disease diagnosis to discover potential underlying neurobiological pathways for AD. To do so, we perform brain-wide genome-wide colocalization analyses to integrate multidimensional imaging genomic biobank data. Specifically, we use (1) the individual-level imputed genotyping data and magnetic resonance imaging (MRI) data from the UK Biobank, (2) the summary statistics of the genome-wide association study (GWAS) from multiple European ancestry cohorts, and (3) the tissue-specific cis-expression quantitative trait loci (cis-eQTL) summary statistics from the GTEx project. We apply a Bayes factor colocalization framework and mediation analysis to these multi-modal imaging genomic data. As a result, we derive the brain regional level GWAS summary statistics for 145 brain regions with 482,831 single nucleotide polymorphisms (SNPs) followed by posthoc functional annotations. Our analysis yields the discovery of a potential AD causal pathway from a systems biology perspective: the SNP chr10:124165615:G>A (rs6585827) mutation upregulates the expression of BTBD16 gene in oligodendrocytes, a specialized glial cells, in the brain cortex, leading to a reduced risk of volumetric loss in the entorhinal cortex, resulting in the protective effect on AD. We substantiate our findings with multiple evidence from existing imaging, genetic and genomic studies in AD literature. Our study connects genetics, molecular and cellular signatures, regional brain morphologic endophenotypes, and AD diagnosis, providing new insights into the mechanistic understanding of the disease. Our findings can provide valuable guidance for subsequent therapeutic target identification and drug discovery in AD.

Integration of multiple-omics data to reveal the shared genetic architecture of educational attainment, intelligence, cognitive performance, and Alzheimer's disease

Growing evidence suggests the effect of educational attainment (EA) on Alzheimer's disease (AD), but less is known about the shared genetic architecture between them. Here, leveraging genome-wide association studies (GWAS) for AD (N = 21,982/41,944), EA (N = 1,131,881), cognitive performance (N = 257,828), and intelligence (N = 78,308), we investigated their causal association with the linkage disequilibrium score (LDSC) and Mendelian randomization and their shared loci with the conjunctional false discovery rate (conjFDR), transcriptome-wide association studies (TWAS), and colocalization. We observed significant genetic correlations of EA (rg = -0.22, p = 5.07E-05), cognitive performance (rg = -0.27, p = 2.44E-05), and intelligence (rg = -0.30, p = 3.00E-04) with AD, and a causal relationship between EA and AD (OR = 0.74, 95% CI: 0.58-0.94, p = 0.013). We identified 13 shared loci at conjFDR <0.01, of which five were novel, and prioritized three causal genes. These findings inform early prevention strategies for AD.

ADRA2A and IRX1 are putative risk genes for Raynaud's phenomenon

Raynaud's phenomenon (RP) is a common vasospastic disorder that causes severe pain and ulcers, but despite its high reported heritability, no causal genes have been robustly identified. We conducted a genome-wide association study including 5,147 RP cases and 439,294 controls, based on diagnoses from electronic health records, and identified three unreported genomic regions associated with the risk of RP (p < 5 × 10-8). We prioritized ADRA2A (rs7090046, odds ratio (OR) per allele: 1.26; 95%-CI: 1.20-1.31; p < 9.6 × 10-27) and IRX1 (rs12653958, OR: 1.17; 95%-CI: 1.12-1.22, p < 4.8 × 10-13) as candidate causal genes through integration of gene expression in disease relevant tissues. We further identified a likely causal detrimental effect of low fasting glucose levels on RP risk (rG = -0.21; p-value = 2.3 × 10-3), and systematically highlighted drug repurposing opportunities, like the antidepressant mirtazapine. Our results provide the first robust evidence for a strong genetic contribution to RP and highlight a so far underrated role of α2A-adrenoreceptor signalling, encoded at ADRA2A, as a possible mechanism for hypersensitivity to catecholamine-induced vasospasms.

Integrated Analyses of Single-Cell Transcriptome and Mendelian Randomization Reveal the Protective Role of Resistin in Sepsis Survival in Intensive Care Unit

The high morbidity and mortality rates associated with sepsis highlight the challenges of finding specific remedies for this condition in the intensive care unit (ICU). This study aimed to explore the differentially expressed genes (DEGs) specific to cell types in sepsis and investigate the role of resistin in the survival of sepsis patients through Mendelian randomization (MR) analyses. We used single-cell and bulk transcriptome data to identify cell type-specific DEGs between sepsis and healthy controls. MR analyses were then conducted to investigate the causal relationships between resistin (one of the identified DEGs) levels and the survival of sepsis patients. Additionally, we utilized meQTL (methylation quantitative trait loci) to identify cytosine-phosphate-guanine (CpG) sites that may directly affect sepsis. We identified 560 cell type-specific DEGs between sepsis and healthy controls. Notably, we observed the upregulation of resistin levels in macrophages during sepsis. In bulk transcriptome, RETN is also upregulated in sepsis samples compared with healthy controls. MR analyses revealed a negative association existed between the expression of resistin, at both gene and protein levels, and the mortality or severity of sepsis patients in ICU. Moreover, there were no associations observed between resistin levels and death or organ failure due to other causes. We also identified three methylation CpG sites, located in RETN or its promoter region-cg06633066, cg22322184, and cg02346997-that directly affected both resistin protein levels and sepsis death in the ICU. Our findings suggest that resistin may provide feasible protection for sepsis patients, particularly those with severe cases, without serious side effects. Therefore, resistin could be a potential drug candidate for sepsis treatment. Additionally, we identified two CpG sites, cg06633066 and cg22322184, that were associated with RETN protein levels and sepsis death, providing novel insights into the underlying mechanisms of sepsis.

Functional genomics identify causal variant underlying the protective CTSH locus for Alzheimer's disease

Alzheimer's disease (AD) is the most prevalent age-related neurodegenerative disease, which has a high heritability of up to 79%. Exploring the genetic basis is essential for understanding the pathogenic mechanisms underlying AD development. Recent genome-wide association studies (GWASs) reported an AD-associated signal in the Cathepsin H (CTSH) gene in European populations. However, the exact functional/causal variant(s), and the genetic regulating mechanism of CTSH in AD remain to be determined. In this study, we carried out a comprehensive study to characterize the role of CTSH variants in the pathogenesis of AD. We identified rs2289702 in CTSH as the most significant functional variant that is associated with a protective effect against AD. The genetic association between rs2289702 and AD was validated in independent cohorts of the Han Chinese population. The CTSH mRNA expression level was significantly increased in AD patients and AD animal models, and the protective allele T of rs2289702 was associated with a decreased expression level of CTSH through the disruption of the binding affinity of transcription factors. Human microglia cells with CTSH knockout showed a significantly increased phagocytosis of Aβ peptides. Our study identified CTSH as being involved in AD genetic susceptibility and uncovered the genetic regulating mechanism of CTSH in pathogenesis of AD.

Identifying proteomic risk factors for overall, aggressive and early onset prostate cancer using Mendelian randomization and tumor spatial transcriptomics

Background

Understanding the role of circulating proteins in prostate cancer risk can reveal key biological pathways and identify novel targets for cancer prevention.

Methods

We investigated the association of 2,002 genetically predicted circulating protein levels with risk of prostate cancer overall, and of aggressive and early onset disease, using cis-pQTL Mendelian randomization (MR) and colocalization. Findings for proteins with support from both MR, after correction for multiple-testing, and colocalization were replicated using two independent cancer GWAS, one of European and one of African ancestry. Proteins with evidence of prostate-specific tissue expression were additionally investigated using spatial transcriptomic data in prostate tumor tissue to assess their role in tumor aggressiveness. Finally, we mapped risk proteins to drug and ongoing clinical trials targets.

Results

We identified 20 proteins genetically linked to prostate cancer risk (14 for overall [8 specific], 7 for aggressive [3 specific], and 8 for early onset disease [2 specific]), of which a majority were novel and replicated. Among these were proteins associated with aggressive disease, such as PPA2 [Odds Ratio (OR) per 1 SD increment = 2.13, 95% CI: 1.54-2.93], PYY [OR = 1.87, 95% CI: 1.43-2.44] and PRSS3 [OR = 0.80, 95% CI: 0.73-0.89], and those associated with early onset disease, including EHPB1 [OR = 2.89, 95% CI: 1.99-4.21], POGLUT3 [OR = 0.76, 95% CI: 0.67-0.86] and TPM3 [OR = 0.47, 95% CI: 0.34-0.64]. We confirm an inverse association of MSMB with prostate cancer overall [OR = 0.81, 95% CI: 0.80-0.82], and also find an inverse association with both aggressive [OR = 0.84, 95% CI: 0.82-0.86] and early onset disease [OR = 0.71, 95% CI: 0.68-0.74]. Using spatial transcriptomics data, we identified MSMB as the genome-wide top-most predictive gene to distinguish benign regions from high grade cancer regions that had five-fold lower MSMB expression. Additionally, ten proteins that were associated with prostate cancer risk mapped to existing therapeutic interventions.

Conclusion

Our findings emphasize the importance of proteomics for improving our understanding of prostate cancer etiology and of opportunities for novel therapeutic interventions. Additionally, we demonstrate the added benefit of in-depth functional analyses to triangulate the role of risk proteins in the clinical aggressiveness of prostate tumors. Using these integrated methods, we identify a subset of risk proteins associated with aggressive and early onset disease as priorities for investigation for the future prevention and treatment of prostate cancer.

Potential drug targets for asthma identified in the plasma and brain through Mendelian randomization analysis

Background

Asthma is a heterogeneous disease, and the involvement of neurogenic inflammation is crucial in its development. The standardized treatments focus on alleviating symptoms. Despite the availability of medications for asthma, they have proven to be inadequate in controlling relapses and halting the progression of the disease. Therefore, there is a need for novel drug targets to prevent asthma.

Methods

We utilized Mendelian randomization to investigate potential drug targets for asthma. We analyzed summary statistics from the UK Biobank and then replicated our findings in GWAS data by Demenais et al. and the FinnGen cohort. We obtained genetic instruments for 734 plasma and 73 brain proteins from recently reported GWAS. Next, we utilized reverse causal relationship analysis, Bayesian co-localization, and phenotype scanning as part of our sensitivity analysis. Furthermore, we performed a comparison and protein-protein interaction analysis to identify causal proteins. We also analyzed the possible consequences of our discoveries by the given existing asthma drugs and their targets.

Results

Using Mendelian randomization analysis, we identified five protein-asthma pairs that were significant at the Bonferroni level (P < 6.35 × 10-5). Specifically, in plasma, we found that an increase of one standard deviation in IL1R1 and ECM1 was associated with an increased risk of asthma, while an increase in ADAM19 was found to be protective. The corresponding odds ratios were 1.03 (95% CI, 1.02-1.04), 1.00 (95% CI, 1.00-1.01), and 0.99 (95% CI, 0.98-0.99), respectively. In the brain, per 10-fold increase in ECM1 (OR, 1.05; 95% CI, 1.03-1.08) and PDLIM4 (OR, 1.05; 95% CI, 1.04-1.07) increased the risk of asthma. Bayesian co-localization found that ECM1 in the plasma (coloc.abf-PPH4 = 0.965) and in the brain (coloc.abf-PPH4 = 0.931) shared the same mutation with asthma. The target proteins of current asthma medications were found to interact with IL1R1. IL1R1 and PDLIM4 were validated in two replication cohorts.

Conclusion

Our integrative analysis revealed that asthma risk is causally affected by the levels of IL1R1, ECM1, and PDLIM4. The results suggest that these three proteins have the potential to be used as drug targets for asthma, and further investigation through clinical trials is needed.

Alcohol intake and endogenous sex hormones in women: meta-analysis of cohort studies and Mendelian randomization

Background

The mechanisms underlying alcohol-induced breast carcinogenesis are not fully understood but may involve hormonal changes.

Methods

We investigated cross-sectional associations between self-reported alcohol intake and serum or plasma concentrations of oestradiol, oestrone, progesterone (in pre-menopausal women only), testosterone, androstenedione, DHEAS (dehydroepiandrosterone sulphate) and SHBG (sex hormone binding globulin) in 45 431 pre-menopausal and 173 476 post-menopausal women. We performed multivariable linear regression separately for UK Biobank, EPIC (European Prospective Investigation into Cancer and Nutrition) and EHBCCG (Endogenous Hormones and Breast Cancer Collaborative Group), and meta-analysed the results. For testosterone and SHBG, we also conducted two-sample Mendelian Randomization (MR) and colocalisation using the ADH1B (Alcohol Dehydrogenase 1B) variant (rs1229984).

Results

Alcohol intake was positively, though weakly, associated with all hormones (except progesterone in pre-menopausal women), with increments in concentrations per 10 g/day increment in alcohol intake ranging from 1.7% for luteal oestradiol to 6.6% for post-menopausal DHEAS. There was an inverse association of alcohol with SHBG in post-menopausal women but a small positive association in pre-menopausal women. MR identified positive associations of alcohol intake with total testosterone (difference per 10 g/day increment: 4.1%; 95% CI: 0.6%, 7.6%) and free testosterone (7.8%; 4.1%, 11.5%), and an inverse association with SHBG (-8.1%; -11.3%, -4.9%). Colocalisation suggested a shared causal locus at ADH1B between alcohol intake and higher free testosterone and lower SHBG (PP4: 0.81 and 0.97 respectively).

Conclusions

Alcohol intake was associated with small increases in sex hormone concentrations, including bioavailable fractions, which may contribute to its effect on breast cancer risk.

Time-to-Event Genome-Wide Association Study for Incident Cardiovascular Disease in People with Type 2 Diabetes Mellitus

BACKGROUND

Type 2 diabetes mellitus (T2D) confers a two- to three-fold increased risk of cardiovascular disease (CVD). However, the mechanisms underlying increased CVD risk among people with T2D are only partially understood. We hypothesized that a genetic association study among people with T2D at risk for developing incident cardiovascular complications could provide insights into molecular genetic aspects underlying CVD.

METHODS

From 16 studies of the Cohorts for Heart & Aging Research in Genomic Epidemiology (CHARGE) Consortium, we conducted a multi-ancestry time-to-event genome-wide association study (GWAS) for incident CVD among people with T2D using Cox proportional hazards models. Incident CVD was defined based on a composite of coronary artery disease (CAD), stroke, and cardiovascular death that occurred at least one year after the diagnosis of T2D. Cohort-level estimated effect sizes were combined using inverse variance weighted fixed effects meta-analysis. We also tested 204 known CAD variants for association with incident CVD among patients with T2D.

RESULTS

A total of 49,230 participants with T2D were included in the analyses (31,118 European ancestries and 18,112 non-European ancestries) which consisted of 8,956 incident CVD cases over a range of mean follow-up duration between 3.2 and 33.7 years (event rate 18.2%). We identified three novel, distinct genetic loci for incident CVD among individuals with T2D that reached the threshold for genome-wide significance (P<5.0×10-8): rs147138607 (intergenic variant between CACNA1E and ZNF648) with a hazard ratio (HR) 1.23, 95% confidence interval (CI) 1.15 - 1.32, P=3.6×10-9, rs11444867 (intergenic variant near HS3ST1) with HR 1.89, 95% CI 1.52 - 2.35, P=9.9×10-9, and rs335407 (intergenic variant between TFB1M and NOX3) HR 1.25, 95% CI 1.16 - 1.35, P=1.5×10-8. Among 204 known CAD loci, 32 were associated with incident CVD in people with T2D with P<0.05, and 5 were significant after Bonferroni correction (P<0.00024, 0.05/204). A polygenic score of these 204 variants was significantly associated with incident CVD with HR 1.14 (95% CI 1.12 - 1.16) per 1 standard deviation increase (P=1.0×10-16).

CONCLUSIONS

The data point to novel and known genomic regions associated with incident CVD among individuals with T2D.

Genetically mimicked effects of ASGR1 inhibitors on all-cause mortality and health outcomes: a drug-target Mendelian randomization study and a phenome-wide association study

BACKGROUND

Asialoglycoprotein receptor 1 (ASGR1) is emerging as a potential drug target to reduce low-density lipoprotein (LDL)-cholesterol and coronary artery disease (CAD) risk. Here, we investigated genetically mimicked ASGR1 inhibitors on all-cause mortality and any possible adverse effects.

METHODS

We conducted a drug-target Mendelian randomization study to assess genetically mimicked effects of ASGR1 inhibitors on all-cause mortality and 25 a priori outcomes relevant to lipid traits, CAD, and possible adverse effects, i.e. liver function, cholelithiasis, adiposity and type 2 diabetes. We also performed a phenome-wide association study of 1951 health-related phenotypes to identify any novel effects. Associations found were compared with those for currently used lipid modifiers, assessed using colocalization, and replicated where possible.

RESULTS

Genetically mimicked ASGR1 inhibitors were associated with a longer lifespan (3.31 years per standard deviation reduction in LDL-cholesterol, 95% confidence interval 1.01 to 5.62). Genetically mimicked ASGR1 inhibitors were inversely associated with apolipoprotein B (apoB), triglycerides (TG) and CAD risk. Genetically mimicked ASGR1 inhibitors were positively associated with alkaline phosphatase, gamma glutamyltransferase, erythrocyte traits, insulin-like growth factor 1 (IGF-1) and C-reactive protein (CRP), but were inversely associated with albumin and calcium. Genetically mimicked ASGR1 inhibitors were not associated with cholelithiasis, adiposity or type 2 diabetes. Associations with apoB and TG were stronger for ASGR1 inhibitors compared with currently used lipid modifiers, and most non-lipid effects were specific to ASGR1 inhibitors. The probabilities for colocalization were > 0.80 for most of these associations, but were 0.42 for lifespan and 0.30 for CAD. These associations were replicated using alternative genetic instruments and other publicly available genetic summary statistics.

CONCLUSIONS

Genetically mimicked ASGR1 inhibitors reduced all-cause mortality. Beyond lipid-lowering, genetically mimicked ASGR1 inhibitors increased liver enzymes, erythrocyte traits, IGF-1 and CRP, but decreased albumin and calcium.

Shared genetic architecture of blood eosinophil counts and asthma in UK Biobank

Rationale

Asthma is a complex, heterogeneous disease strongly associated with type 2 inflammation, and blood eosinophil counts guide therapeutic interventions in moderate and severe asthma. Eosinophils are leukocytes involved in type 2 immune responses. Despite these critical associations between asthma and blood eosinophil counts, the shared genetic architecture of these two traits remains unknown. The objective of the present study was to characterise the genetic architecture of blood eosinophil counts and asthma in the UK Biobank.

Methods

We performed genome-wide association studies (GWAS) of doctor-diagnosed asthma, blood eosinophil, neutrophil, lymphocyte and monocyte counts in the UK Biobank. Genetic correlation analysis was performed on GWAS results and validated in the Trans-National Asthma Genetic Consortium (TAGC) study of asthma.

Results

GWAS of doctor-diagnosed asthma and blood eosinophil counts in the UK Biobank identified 585 and 3429 significant variants, respectively. STAT6, a transcription factor involved in interleukin-4 signalling, was a key shared pathway between asthma and blood eosinophil counts. Genetic correlation analysis demonstrated a positive correlation between doctor-diagnosed asthma and blood eosinophil counts (r=0.38±0.10, correlation±se; p=4.7×10-11). As a validation of this association, we found a similar correlation between TAGC and blood eosinophil counts in the UK Biobank (0.37±0.08, correlation±se; p=1.2×10-6).

Conclusions

These findings define the shared genetic architecture between blood eosinophil counts and asthma risk in subjects of European ancestry and point to a genetic link to the STAT6 signalling pathway in these two traits.

Interaction molecular QTL mapping discovers cellular and environmental modifiers of genetic regulatory effects

Bulk tissue molecular quantitative trait loci (QTLs) have been the starting point for interpreting disease-associated variants, while context-specific QTLs show particular relevance for disease. Here, we present the results of mapping interaction QTLs (iQTLs) for cell type, age, and other phenotypic variables in multi-omic, longitudinal data from blood of individuals of diverse ancestries. By modeling the interaction between genotype and estimated cell type proportions, we demonstrate that cell type iQTLs could be considered as proxies for cell type-specific QTL effects. The interpretation of age iQTLs, however, warrants caution as the moderation effect of age on the genotype and molecular phenotype association may be mediated by changes in cell type composition. Finally, we show that cell type iQTLs contribute to cell type-specific enrichment of diseases that, in combination with additional functional data, may guide future functional studies. Overall, this study highlights iQTLs to gain insights into the context-specificity of regulatory effects.

Predicting mechanisms of action at genetic loci associated with discordant effects on type 2 diabetes and abdominal fat accumulation

Obesity is a major risk factor for cardiovascular disease, stroke, and type 2 diabetes (T2D). Excessive accumulation of fat in the abdomen further increases T2D risk. Abdominal obesity is measured by calculating the ratio of waist-to-hip circumference adjusted for the body-mass index (WHRadjBMI), a trait with a significant genetic inheritance. Genetic loci associated with WHRadjBMI identified in genome-wide association studies are predicted to act through adipose tissues, but many of the exact molecular mechanisms underlying fat distribution and its consequences for T2D risk are poorly understood. Further, mechanisms that uncouple the genetic inheritance of abdominal obesity from T2D risk have not yet been described. Here we utilize multi-omic data to predict mechanisms of action at loci associated with discordant effects on abdominal obesity and T2D risk. We find six genetic signals in five loci associated with protection from T2D but also with increased abdominal obesity. We predict the tissues of action at these discordant loci and the likely effector Genes (eGenes) at three discordant loci, from which we predict significant involvement of adipose biology. We then evaluate the relationship between adipose gene expression of eGenes with adipogenesis, obesity, and diabetic physiological phenotypes. By integrating these analyses with prior literature, we propose models that resolve the discordant associations at two of the five loci. While experimental validation is required to validate predictions, these hypotheses provide potential mechanisms underlying T2D risk stratification within abdominal obesity.

Gene expression in African Americans, Puerto Ricans and Mexican Americans reveals ancestry-specific patterns of genetic architecture

We explored ancestry-related differences in the genetic architecture of whole-blood gene expression using whole-genome and RNA sequencing data from 2,733 African Americans, Puerto Ricans and Mexican Americans. We found that heritability of gene expression significantly increased with greater proportions of African genetic ancestry and decreased with higher proportions of Indigenous American ancestry, reflecting the relationship between heterozygosity and genetic variance. Among heritable protein-coding genes, the prevalence of ancestry-specific expression quantitative trait loci (anc-eQTLs) was 30% in African ancestry and 8% for Indigenous American ancestry segments. Most anc-eQTLs (89%) were driven by population differences in allele frequency. Transcriptome-wide association analyses of multi-ancestry summary statistics for 28 traits identified 79% more gene-trait associations using transcriptome prediction models trained in our admixed population than models trained using data from the Genotype-Tissue Expression project. Our study highlights the importance of measuring gene expression across large and ancestrally diverse populations for enabling new discoveries and reducing disparities.

Biomarker and genomic analyses reveal molecular signatures of non-cardioembolic ischemic stroke

Acute ischemic stroke (AIS) is a major cause of disability and mortality worldwide. Non-cardioembolic ischemic stroke (NCIS), which constitutes the majority of AIS cases, is highly heterogeneous, thus requiring precision medicine treatments. This study aimed to investigate the molecular mechanisms underlying NCIS heterogeneity. We integrated data from the Third China National Stroke Registry, including clinical phenotypes, biomarkers, and whole-genome sequencing data for 7695 patients with NCIS. We identified 30 molecular clusters based on 63 biomarkers and explored the comprehensive landscape of biological heterogeneity and subpopulations in NCIS. Dimensionality reduction revealed fine-scale subpopulation structures associated with specific biomarkers. The subpopulations with biomarkers for inflammation, abnormal liver and kidney function, homocysteine metabolism, lipid metabolism, and gut microbiota metabolism were associated with a high risk of unfavorable clinical outcomes, including stroke recurrence, disability, and mortality. Several genes encoding potential drug targets were identified as putative causal genes that drive the clusters, such as CDK10, ERCC3, and CHEK2. We comprehensively characterized the genetic architecture of these subpopulations, identified their molecular signatures, and revealed the potential of the polybiomarkers and polygenic prediction for assessing clinical outcomes. Our study demonstrates the power of large-scale molecular biomarkers and genomics to understand the underlying biological mechanisms of and advance precision medicine for NCIS.

Body Fat Distribution, Fasting Insulin Levels, and Insulin Secretion: A Bidirectional Mendelian Randomization Study

CONTEXT

Hyperinsulinemia and adiposity are associated with one another, but the directionality of this relation is debated.

OBJECTIVE

Here, we tested the direction of the causal effects of fasting insulin (FI) levels and body fat accumulation/distribution using 2-sample bidirectional Mendelian randomization (MR).

METHODS

We included summary statistics from large-scale genome-wide association studies for body mass index (BMI, n = 806 834), waist to hip ratio adjusted for BMI (WHRadjBMI, n = 694 649), abdominal subcutaneous, visceral and gluteofemoral adipose tissue (n = 38 965), FI levels (n = 98 210), pancreatic islets gene expression (n = 420), and hypothalamus gene expression (n = 155). We used inverse variance-weighted and robust MR methods that relied on statistically and biologically driven genetic instruments.

RESULTS

Both BMI and WHRadjBMI were positively associated with FI. Results were consistent across all robust MR methods and when variants mapped to the hypothalamus (presumably associated with food behavior) were included. In multivariable MR analyses, when waist circumference and BMI were mutually adjusted, the direct effect of waist circumference on FI was 2.43 times larger than the effect of BMI on FI. FI was not associated with adiposity. By contrast, using genetic instruments mapped to gene expression in pancreatic islets (presumably more specific to insulin secretion), insulin was positively associated with BMI and abdominal subcutaneous and gluteofemoral adipose tissue, but not with visceral adipose tissue.

CONCLUSION

Although these results will need to be supported by experimental investigations, results of this MR study suggest that abdominal adiposity may be a key determinant of circulating insulin levels. Alternatively, insulin secretion may promote peripheral adipose tissue accumulation.

FHL5 Controls Vascular Disease-Associated Gene Programs in Smooth Muscle Cells

BACKGROUND

Genome-wide association studies have identified hundreds of loci associated with common vascular diseases, such as coronary artery disease, myocardial infarction, and hypertension. However, the lack of mechanistic insights for many GWAS loci limits their translation into the clinic. Among these loci with unknown functions is UFL1-four-and-a-half LIM (LIN-11, Isl-1, MEC-3) domain 5 (FHL5; chr6q16.1), which reached genome-wide significance in a recent coronary artery disease/ myocardial infarction GWAS meta-analysis. UFL1-FHL5 is also associated with several vascular diseases, consistent with the widespread pleiotropy observed for GWAS loci.

METHODS

We apply a multimodal approach leveraging statistical fine-mapping, epigenomic profiling, and ex vivo analysis of human coronary artery tissues to implicate FHL5 as the top candidate causal gene. We unravel the molecular mechanisms of the cross-phenotype genetic associations through in vitro functional analyses and epigenomic profiling experiments in coronary artery smooth muscle cells.

RESULTS

We prioritized FHL5 as the top candidate causal gene at the UFL1-FHL5 locus through expression quantitative trait locus colocalization methods. FHL5 gene expression was enriched in the smooth muscle cells and pericyte population in human artery tissues with coexpression network analyses supporting a functional role in regulating smooth muscle cell contraction. Unexpectedly, under procalcifying conditions, FHL5 overexpression promoted vascular calcification and dysregulated processes related to extracellular matrix organization and calcium handling. Lastly, by mapping FHL5 binding sites and inferring FHL5 target gene function using artery tissue gene regulatory network analyses, we highlight regulatory interactions between FHL5 and downstream coronary artery disease/myocardial infarction loci, such as FOXL1 and FN1 that have roles in vascular remodeling.

CONCLUSIONS

Taken together, these studies provide mechanistic insights into the pleiotropic genetic associations of UFL1-FHL5. We show that FHL5 mediates vascular disease risk through transcriptional regulation of downstream vascular remodeling gene programs. These transacting mechanisms may explain a portion of the heritable risk for complex vascular diseases.

Common genetic variations in telomere length genes and lung cancer: a Mendelian randomisation study and its novel application in lung tumour transcriptome

Background

Genome-wide association studies (GWASs) have identified genetic susceptibility variants for both leukocyte telomere length (LTL) and lung cancer susceptibility. Our study aims to explore the shared genetic basis between these traits and investigate their impact on somatic environment of lung tumours.

Methods

We performed genetic correlation, Mendelian randomisation (MR), and colocalisation analyses using the largest available GWASs summary statistics of LTL (N=464,716) and lung cancer (N=29,239 cases and 56,450 controls). Principal components analysis based on RNA-sequencing data was used to summarise gene expression profile in lung adenocarcinoma cases from TCGA (N=343).

Results

Although there was no genome-wide genetic correlation between LTL and lung cancer risk, longer LTL conferred an increased risk of lung cancer regardless of smoking status in the MR analyses, particularly for lung adenocarcinoma. Of the 144 LTL genetic instruments, 12 colocalised with lung adenocarcinoma risk and revealed novel susceptibility loci, including MPHOSPH6, PRPF6, and POLI. The polygenic risk score for LTL was associated with a specific gene expression profile (PC2) in lung adenocarcinoma tumours. The aspect of PC2 associated with longer LTL was also associated with being female, never smokers, and earlier tumour stages. PC2 was strongly associated with cell proliferation score and genomic features related to genome stability, including copy number changes and telomerase activity.

Conclusions

This study identified an association between longer genetically predicted LTL and lung cancer and sheds light on the potential molecular mechanisms related to LTL in lung adenocarcinomas.

Funding

Institut National du Cancer (GeniLuc2017-1-TABAC-03-CIRC-1-TABAC17-022), INTEGRAL/NIH (5U19CA203654-03), CRUK (C18281/A29019), and Agence Nationale pour la Recherche (ANR-10-INBS-09).

Multi-omic underpinnings of epigenetic aging and human longevity

Biological aging is accompanied by increasing morbidity, mortality, and healthcare costs; however, its molecular mechanisms are poorly understood. Here, we use multi-omic methods to integrate genomic, transcriptomic, and metabolomic data and identify biological associations with four measures of epigenetic age acceleration and a human longevity phenotype comprising healthspan, lifespan, and exceptional longevity (multivariate longevity). Using transcriptomic imputation, fine-mapping, and conditional analysis, we identify 22 high confidence associations with epigenetic age acceleration and seven with multivariate longevity. FLOT1, KPNA4, and TMX2 are novel, high confidence genes associated with epigenetic age acceleration. In parallel, cis-instrument Mendelian randomization of the druggable genome associates TPMT and NHLRC1 with epigenetic aging, supporting transcriptomic imputation findings. Metabolomics Mendelian randomization identifies a negative effect of non-high-density lipoprotein cholesterol and associated lipoproteins on multivariate longevity, but not epigenetic age acceleration. Finally, cell-type enrichment analysis implicates immune cells and precursors in epigenetic age acceleration and, more modestly, multivariate longevity. Follow-up Mendelian randomization of immune cell traits suggests lymphocyte subpopulations and lymphocytic surface molecules affect multivariate longevity and epigenetic age acceleration. Our results highlight druggable targets and biological pathways involved in aging and facilitate multi-omic comparisons of epigenetic clocks and human longevity.

Genetic effects on the timing of parturition and links to fetal birth weight

The timing of parturition is crucial for neonatal survival and infant health. Yet, its genetic basis remains largely unresolved. We present a maternal genome-wide meta-analysis of gestational duration (n = 195,555), identifying 22 associated loci (24 independent variants) and an enrichment in genes differentially expressed during labor. A meta-analysis of preterm delivery (18,797 cases, 260,246 controls) revealed six associated loci and large genetic similarities with gestational duration. Analysis of the parental transmitted and nontransmitted alleles (n = 136,833) shows that 15 of the gestational duration genetic variants act through the maternal genome, whereas 7 act both through the maternal and fetal genomes and 2 act only via the fetal genome. Finally, the maternal effects on gestational duration show signs of antagonistic pleiotropy with the fetal effects on birth weight: maternal alleles that increase gestational duration have negative fetal effects on birth weight. The present study provides insights into the genetic effects on the timing of parturition and the complex maternal-fetal relationship between gestational duration and birth weight.

Insights into the pleiotropic relationships between chronic back pain and inflammation-related musculoskeletal conditions: rheumatoid arthritis and osteoporotic abnormalities

ABSTRACT

The ageing process includes the development of debilitating musculoskeletal (MSK) conditions, including chronic back pain (CBP), rheumatoid arthritis (RA), and osteoporosis (OP). The mechanisms involved in the genetic-epidemiological relationships between these MSK phenotypes are controversial and limited and thus require clarification, in particular, between CBP and the other MSK phenotypes. A cross-sectional statistical analysis was conducted using Europeans from the UK Biobank data collection, including 73,794 CBP, 4883 RA, and 7153 OP cases as well as 242,216 calcaneus bone mineral density scores. C-reactive protein (CRP) was measured for 402,165 subjects in this sample. Genetic correlations were assessed to evaluate shared genetic background between traits. Mendelian randomization was performed to assess a causal relationship between CBP and RA and OP along with other risk factors, such as CRP. Colocalization analysis was conducted to identify shared pleiotropic regions between the examined traits. Bayesian modelling was performed to determine a potential pathway that may explain the interrelationships among these traits. Mendelian randomization analyses revealed that CRP causally predicts CBP only (β = 0.183, 95% CI = 0.077-0.290, P -value = 0.001). Horizontally pleiotropy appeared to explain the relationship between CBP and RA and OP. Through colocalization analysis, several genomic regions emerged describing common genetic influences between CBP and its proposed risk factors, including HLA-DQA1/HLA-DQB1, APOE , SOX5, and MYH7B as well as Histone 1 genes. We speculate that among other factors, CBP and its MSK comorbidities may arise from common inflammatory mechanisms. Colocalized identified genes may aid in advancing or improving the mode of treatment in patients with CBP.

Brain expression quantitative trait locus and network analyses reveal downstream effects and putative drivers for brain-related diseases

Identification of therapeutic targets from genome-wide association studies (GWAS) requires insights into downstream functional consequences. We harmonized 8,613 RNA-sequencing samples from 14 brain datasets to create the MetaBrain resource and performed cis- and trans-expression quantitative trait locus (eQTL) meta-analyses in multiple brain region- and ancestry-specific datasets (n ≤ 2,759). Many of the 16,169 cortex cis-eQTLs were tissue-dependent when compared with blood cis-eQTLs. We inferred brain cell types for 3,549 cis-eQTLs by interaction analysis. We prioritized 186 cis-eQTLs for 31 brain-related traits using Mendelian randomization and co-localization including 40 cis-eQTLs with an inferred cell type, such as a neuron-specific cis-eQTL (CYP24A1) for multiple sclerosis. We further describe 737 trans-eQTLs for 526 unique variants and 108 unique genes. We used brain-specific gene-co-regulation networks to link GWAS loci and prioritize additional genes for five central nervous system diseases. This study represents a valuable resource for post-GWAS research on central nervous system diseases.

Genome-wide meta-analysis identifies 93 risk loci and enables risk prediction equivalent to monogenic forms of venous thromboembolism

We report a genome-wide association study of venous thromboembolism (VTE) incorporating 81,190 cases and 1,419,671 controls sampled from six cohorts. We identify 93 risk loci, of which 62 are previously unreported. Many of the identified risk loci are at genes encoding proteins with functions converging on the coagulation cascade or platelet function. A VTE polygenic risk score (PRS) enabled effective identification of both high- and low-risk individuals. Individuals within the top 0.1% of PRS distribution had a VTE risk similar to homozygous or compound heterozygous carriers of the variants G20210A (c.*97 G > A) in F2 and p.R534Q in F5. We also document that F2 and F5 mutation carriers in the bottom 10% of the PRS distribution had a risk similar to that of the general population. We further show that PRS improved individual risk prediction beyond that of genetic and clinical risk factors. We investigated the extent to which venous and arterial thrombosis share clinical risk factors using Mendelian randomization, finding that some risk factors for arterial thrombosis were directionally concordant with VTE risk (for example, body mass index and smoking) whereas others were discordant (for example, systolic blood pressure and triglyceride levels).

ANKLE1 cleaves mitochondrial DNA and contributes to cancer risk by promoting apoptosis resistance and metabolic dysregulation

Alleles within the chr19p13.1 locus are associated with increased risk of both ovarian and breast cancer and increased expression of the ANKLE1 gene. ANKLE1 is molecularly characterized as an endonuclease that efficiently cuts branched DNA and shuttles between the nucleus and cytoplasm. However, the role of ANKLE1 in mammalian development and homeostasis remains unknown. In normal development ANKLE1 expression is limited to the erythroblast lineage and we found that ANKLE1's role is to cleave the mitochondrial genome during erythropoiesis. We show that ectopic expression of ANKLE1 in breast epithelial-derived cells leads to genome instability and mitochondrial DNA (mtDNA) cleavage. mtDNA degradation then leads to mitophagy and causes a shift from oxidative phosphorylation to glycolysis (Warburg effect). Moreover, mtDNA degradation activates STAT1 and expression of epithelial-mesenchymal transition (EMT) genes. Reduction in mitochondrial content contributes to apoptosis resistance, which may allow precancerous cells to avoid apoptotic checkpoints and proliferate. These findings provide evidence that ANKLE1 is the causal cancer susceptibility gene in the chr19p13.1 locus and describe mechanisms by which higher ANKLE1 expression promotes cancer risk.