Conditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traits

Hierarchical joint analysis of marginal summary statistics-Part I: Multipopulation fine mapping and credible set construction

Recent advancement in genome-wide association studies (GWAS) comes from not only increasingly larger sample sizes but also the shift in focus towards underrepresented populations. Multipopulation GWAS increase power to detect novel risk variants and improve fine-mapping resolution by leveraging evidence and differences in linkage disequilibrium (LD) from diverse populations. Here, we expand upon our previous approach for single-population fine-mapping through Joint Analysis of Marginal SNP Effects (JAM) to a multipopulation analysis (mJAM). Under the assumption that true causal variants are common across studies, we implement a hierarchical model framework that conditions on multiple SNPs while explicitly incorporating the different LD structures across populations. The mJAM framework can be used to first select index variants using the mJAM likelihood with different feature selection approaches. In addition, we present a novel approach leveraging the ideas of mediation to construct credible sets for these index variants. Construction of such credible sets can be performed given any existing index variants. We illustrate the implementation of the mJAM likelihood through two implementations: mJAM-SuSiE (a Bayesian approach) and mJAM-Forward selection. Through simulation studies based on realistic effect sizes and levels of LD, we demonstrated that mJAM performs well for constructing concise credible sets that include the underlying causal variants. In real data examples taken from the most recent multipopulation prostate cancer GWAS, we showed several practical advantages of mJAM over other existing multipopulation methods.

Genetic architecture of telomere length in 462,666 UK Biobank whole-genome sequences

Telomeres protect chromosome ends from damage and their length is linked with human disease and aging. We developed a joint telomere length metric, combining quantitative PCR and whole-genome sequencing measurements from 462,666 UK Biobank participants. This metric increased SNP heritability, suggesting that it better captures genetic regulation of telomere length. Exome-wide rare-variant and gene-level collapsing association studies identified 64 variants and 30 genes significantly associated with telomere length, including allelic series in ACD and RTEL1. Notably, 16% of these genes are known drivers of clonal hematopoiesis-an age-related somatic mosaicism associated with myeloid cancers and several nonmalignant diseases. Somatic variant analyses revealed gene-specific associations with telomere length, including lengthened telomeres in individuals with large SRSF2-mutant clones, compared with shortened telomeres in individuals with clonal expansions driven by other genes. Collectively, our findings demonstrate the impact of rare variants on telomere length, with larger effects observed among genes also associated with clonal hematopoiesis.

Multitrait Genetic Analysis Identifies Novel Pleiotropic Loci for Depression and Schizophrenia in East Asians

BACKGROUND AND HYPOTHESIS

While genetic correlations, pleiotropic loci, and shared genetic mechanisms of psychiatric disorders have been extensively studied in European populations, the investigation of these factors in East Asian populations has been relatively limited.

STUDY DESIGN

To identify novel pleiotropic risk loci for depression and schizophrenia (SCZ) in East Asians. We utilized the most comprehensive dataset available for East Asians and quantified the genetic overlap between depression, SCZ, and their related traits via a multitrait genome-wide association study. Global and local genetic correlations were estimated by LDSC and ρ-HESS. Pleiotropic loci were identified by the multitrait analysis of GWAS (MTAG).

STUDY RESULTS

Besides the significant correlation between depression and SCZ, our analysis revealed genetic correlations between depression and obesity-related traits, such as weight, BMI, T2D, and HDL. In SCZ, significant correlations were detected with HDL, heart diseases and use of various medications. Conventional meta-analysis of depression and SCZ identified a novel locus at 1q25.2 in East Asians. Further multitrait analysis of depression, SCZ and related traits identified ten novel pleiotropic loci for depression, and four for SCZ.

CONCLUSIONS

Our findings demonstrate shared genetic underpinnings between depression and SCZ in East Asians, as well as their associated traits, providing novel candidate genes for the identification and prioritization of therapeutic targets specific to this population.

Large-Scale Alternative Polyadenylation-Wide Association Studies to Identify Putative Cancer Susceptibility Genes

Alternative polyadenylation (APA) modulates mRNA processing in the 3'-untranslated regions (3' UTR), affecting mRNA stability and translation efficiency. Research into genetically regulated APA has the potential to provide insights into cancer risk. In this study, we conducted large APA-wide association studies to investigate associations between APA levels and cancer risk. Genetic models were built to predict APA levels in multiple tissues using genotype and RNA sequencing data from 1,337 samples from the Genotype-Tissue Expression project. Associations of genetically predicted APA levels with cancer risk were assessed by applying the prediction models to data from large genome-wide association studies of six common cancers among European ancestry populations: breast, ovarian, prostate, colorectal, lung, and pancreatic cancers. A total of 58 risk genes (corresponding to 76 APA sites) were associated with at least one type of cancer, including 25 genes previously not linked to cancer susceptibility. Of the identified risk APAs, 97.4% and 26.3% were supported by 3'-UTR APA quantitative trait loci and colocalization analyses, respectively. Luciferase reporter assays for four selected putative regulatory 3'-UTR variants demonstrated that the risk alleles of 3'-UTR variants, rs324015 (STAT6), rs2280503 (DIP2B), rs1128450 (FBXO38), and rs145220637 (LDHA), significantly increased the posttranscriptional activities of their target genes compared with reference alleles. Furthermore, knockdown of the target genes confirmed their ability to promote proliferation and migration. Overall, this study provides insights into the role of APA in the genetic susceptibility to common cancers. Significance: Systematic evaluation of associations of alternative polyadenylation with cancer risk reveals 58 putative susceptibility genes, highlighting the contribution of genetically regulated alternative polyadenylation of 3'UTRs to genetic susceptibility to cancer.

Genome-wide fine-mapping improves identification of causal variants

Fine-mapping refines genotype-phenotype association signals to identify causal variants underlying complex traits. However, current methods typically focus on individual genomic segments without considering the global genetic architecture. Here, we demonstrate the advantages of performing genome-wide fine-mapping (GWFM) and develop methods to facilitate GWFM. In simulations and real data analyses, GWFM outperforms current methods in error control, mapping power and precision, replication rate, and trans-ancestry phenotype prediction. For 48 well-powered traits in the UK Biobank, we identify causal variants that collectively explain 17% of the SNP-based heritability, and predict that fine-mapping 50% of that would require 2 million samples on average. We pinpoint a known causal variant, as proof-of-principle, at FTO for body mass index, unveil a hidden secondary variant with evolutionary conservation, and identify new missense causal variants for schizophrenia and Crohn's disease. Overall, we analyse 600 complex traits with 13 million SNPs, highlighting the efficacy of GWFM with functional annotations.

Genetic effects on the skin methylome in healthy older twins

Whole-skin DNA methylation variation has been implicated in several diseases, including melanoma, but its genetic basis has not yet been fully characterized. Using bulk skin tissue samples from 414 healthy female UK twins, we performed twin-based heritability and methylation quantitative trait loci (meQTL) analyses for >400,000 DNA methylation sites. We find that the human skin DNA methylome is on average less heritable than previously estimated in blood and other tissues (mean heritability: 10.02%). meQTL analysis identified local genetic effects influencing DNA methylation at 18.8% (76,442) of tested CpG sites, as well as 1,775 CpG sites associated with at least one distal genetic variant. As a functional follow-up, we performed skin expression QTL (eQTL) analyses in a partially overlapping sample of 604 female twins. Colocalization analysis identified over 3,500 shared genetic effects affecting thousands of CpG sites (10,067) and genes (4,475). Mediation analysis of putative colocalized gene-CpG pairs identified 114 genes with evidence for eQTL effects being mediated by DNA methylation in skin, including in genes implicating skin disease such as ALOX12 and CSPG4. We further explored the relevance of skin meQTLs to skin disease and found that skin meQTLs and CpGs under genetic influence were enriched for multiple skin-related genome-wide and epigenome-wide association signals, including for melanoma and psoriasis. Our findings give insights into the regulatory landscape of epigenomic variation in skin.

Genome-wide fine-mapping improves identification of causal variants

Fine-mapping refines genotype-phenotype association signals to identify causal variants underlying complex traits. However, current methods typically focus on individual genomic segments without considering the global genetic architecture. Here, we demonstrate the advantages of performing genome-wide fine-mapping (GWFM) and develop methods to facilitate GWFM. In simulations and real data analyses, GWFM outperforms current methods in error control, mapping power and precision, replication rate, and trans-ancestry phenotype prediction. For 48 well-powered traits in the UK Biobank, we identify causal variants that collectively explain 17% of the SNP-based heritability, and predict that fine-mapping 50% of that would require 2 million samples on average. We pinpoint a known causal variant, as proof-of-principle, at FTO for body mass index, unveil a hidden secondary variant with evolutionary conservation, and identify new missense causal variants for schizophrenia and Crohn's disease. Overall, we analyse 599 complex traits with 13 million SNPs, highlighting the efficacy of GWFM with functional annotations.

A novel framework with automated horizontal pleiotropy adjustment in mendelian randomization

The presence of horizontal pleiotropy in Mendelian randomization (MR) analysis has long been a concern due to its potential to induce substantial bias. In recent years, many robust MR methods have been proposed to address this by relaxing the "no horizontal pleiotropy" assumption. Here, we propose a novel two-stage framework called CMR, which integrates a conditional analysis of multiple genetic variants to remove pleiotropy induced by linkage disequilibrium, followed by the application of robust MR methods to model the conditional genetic effect estimates. We demonstrate how the conditional analysis can reduce horizontal pleiotropy and improve the performance of existing MR methods. Extensive simulation studies covering a wide range of scenarios of horizontal pleiotropy showcased the superior performance of the proposed CMR framework over the standard MR framework in which marginal genetic effects are modeled. Moreover, the application of CMR in a negative control outcome analysis and investigation into the causal role of body mass index across various diseases highlighted its potential to deliver more reliable results in real-world applications.

Case-Case Genome-Wide Analyses Identify Subtype-Informative Variants That Confer Risk for Breast Cancer

Breast cancer includes several subtypes with distinct characteristic biological, pathologic, and clinical features. Elucidating subtype-specific genetic etiology could provide insights into the heterogeneity of breast cancer to facilitate the development of improved prevention and treatment approaches. In this study, we conducted pairwise case-case comparisons among five breast cancer subtypes by applying a case-case genome-wide association study (CC-GWAS) approach to summary statistics data of the Breast Cancer Association Consortium. The approach identified 13 statistically significant loci and eight suggestive loci, the majority of which were identified from comparisons between triple-negative breast cancer (TNBC) and luminal A breast cancer. Associations of lead variants in 12 loci remained statistically significant after accounting for previously reported breast cancer susceptibility variants, among which, two were genome-wide significant. Fine mapping implicated putative functional/causal variants and risk genes at several loci, e.g., 3q26.31/TNFSF10, 8q22.3/NACAP1/GRHL2, and 8q23.3/LINC00536/TRPS1, for TNBC as compared with luminal cancer. Functional investigation further identified rs16867605 at 8q22.3 as a SNP that modulates the enhancer activity of GRHL2. Subtype-informative polygenic risk scores (PRS) were derived, and patients with a high subtype-informative PRS had an up to two-fold increased risk of being diagnosed with TNBC instead of luminal cancers. The CC-GWAS PRS remained statistically significant after adjusting for TNBC PRS derived from traditional case-control GWAS in The Cancer Genome Atlas and the African Ancestry Breast Cancer Genetic Consortium. The CC-GWAS PRS was also associated with overall survival and disease-specific survival among patients with breast cancer. Overall, these findings have advanced our understanding of the genetic etiology of breast cancer subtypes, particularly for TNBC. Significance: The discovery of subtype-informative genetic risk variants for breast cancer advances our understanding of the etiologic heterogeneity of breast cancer, which could accelerate the identification of targets and personalized strategies for prevention and treatment.

Genome-wide association meta-analysis identifies five loci associated with postpartum hemorrhage

Bleeding in early pregnancy and postpartum hemorrhage (PPH) bear substantial risks, with the former closely associated with pregnancy loss and the latter being the foremost cause of maternal death, underscoring the severe impact on maternal-fetal health. We identified five genetic loci linked to PPH in a meta-analysis. Functional annotation analysis indicated candidate genes HAND2, TBX3 and RAP2C/FRMD7 at three loci and showed that at each locus, associated variants were located within binding sites for progesterone receptors. There were strong genetic correlations with birth weight, gestational duration and uterine fibroids. Bleeding in early pregnancy yielded no genome-wide association signals but showed strong genetic correlation with various human traits, suggesting a potentially complex, polygenic etiology. Our results suggest that PPH is related to progesterone signaling dysregulation, whereas early bleeding is a complex trait associated with underlying health and possibly socioeconomic status and may include genetic factors that have not yet been identified.

Causal relationship between air pollution and infections: a two-sample Mendelian randomization study

Background

Traditional observational studies exploring the association between air pollution and infections have been limited by small sample sizes and potential confounding factors. To address these limitations, we applied Mendelian randomization (MR) to investigate the potential causal relationships between particulate matter (PM2.5, PM2.5-10, and PM10), nitrogen dioxide, and nitrogen oxide and the risks of infections.

Methods

Single nucleotide polymorphisms (SNPs) related to air pollution were selected from the genome-wide association study (GWAS) of the UK Biobank. Publicly available summary data for infections were obtained from the FinnGen Biobank and the COVID-19 Host Genetics Initiative. The inverse variance weighted (IVW) meta-analysis was used as the primary method for obtaining the Mendelian randomization (MR) estimates. Complementary analyses were performed using the weighted median method, MR-Egger method, and MR Pleiotropy Residual Sum and Outlier (MR-PRESSO) test.

Results

The fixed-effect IVW estimate showed that PM2.5, PM2.5-10 and Nitrogen oxides were suggestively associated with COVID-19 [for PM2.5: IVW (fe): OR 3.573(1.218,5.288), PIVW(fe) = 0.021; for PM2.5-10: IVW (fe): OR 2.940(1.385,6.239), PIVW(fe) = 0.005; for Nitrogen oxides, IVW (fe): OR 1.898(1.318,2.472), PIVW(fe) = 0.010]. PM2.5, PM2.5-10, PM10, and Nitrogen oxides were suggestively associated with bacterial pneumonia [for PM2.5: IVW(fe): OR 1.720 (1.007, 2.937), PIVW(fe) = 0.047; for PM2.5-10: IVW(fe): OR 1.752 (1.111, 2.767), P IVW(fe) = 0.016; for PM10: IVW(fe): OR 2.097 (1.045, 4.208), PIVW(fe) = 0.037; for Nitrogen oxides, IVW(fe): OR 3.907 (1.209, 5.987), PIVW(fe) = 0.023]. Furthermore, Nitrogen dioxide was suggestively associated with the risk of acute upper respiratory infections, while all air pollution were not associated with intestinal infections.

Conclusions

Our results support a role of related air pollution in the Corona Virus Disease 2019, bacterial pneumonia and acute upper respiratory infections. More work is need for policy formulation to reduce the air pollution and the emission of toxic and of harmful gas.

Effect modification by sex of genetic associations of vitamin C related metabolites in the Canadian Longitudinal study on aging

Introduction: Vitamin C is an essential nutrient. Sex differences in serum vitamin C concentrations have been observed but are not fully known. Investigation of levels of metabolites may help shed light on how dietary and other environmental exposures interact with molecular processes. O-methylascorbate and ascorbic acid 2-sulfate are two metabolites in the vitamin C metabolic pathway. Past research has found genetic factors that influence the levels of these two metabolites. Therefore, we investigated possible effect modification by sex of genetic variant-metabolite associations and characterized the biological function of these interactions. Methods: We included individuals of European descent from the Canadian Longitudinal Study on Aging with available genetic and metabolic data (n = 9004). We used linear mixed models to tests for genome-wide associations with O-methylascorbate and ascorbic acid 2-sulfate, with and without a sex interaction. We also investigated the biological function of the important genetic variant-sex interactions found for each metabolite. Results: Two genome-wide statistically significant (p value < 5 × 10-8) interaction effects and several suggestive (p value < 10-5) interaction effects were found. These suggestive interaction effects were mapped to several genes including HSD11B2, associated with sex hormones, and AGRP, associated with hunger drive. The genes mapped to O-methylascorbate were differently expressed in the testis tissues, and the genes mapped to ascorbic acid 2-sulfate were differently expressed in stomach tissues. Discussion: By understanding the genetic factors that impact metabolites associated with vitamin C, we can better understand its function in disease risk and the mechanisms behind sex differences in vitamin C concentrations.

The genetic effects of hormones modulated by the Pituitary-Thyroid/Adrenal/Gonadal axis on the risk of developing venous thromboembolism: a mendelian randomization study

BACKGROUND

The aim of this study was to explore the genetic effects of hormones modulated through the pituitary-thyroid/adrenal/gonadal axis on the risk of developing venous thromboembolism (VTE) and to investigate the potentially causal relationships between them.

METHODS

A two-sample Mendelian randomization (MR) design was used. The single-nucleotide polymorphisms (SNPs) used as instrumental variables for various hormones and hormone-mediated diseases were derived from published genome-wide association studies (GWASs). Summary statistics for the risk of developing VTE (including deep venous thrombosis [DVT] and pulmonary embolism [PE]) were obtained from the UK Biobank and the FinnGen consortium. Inverse-variance weighting (IVW) was applied as the primary method to analyse causal associations. Other MR methods were used for supplementary estimates and sensitivity analysis.

RESULTS

A genetic predisposition to greater free thyroxine (FT4) concentrations was associated with a greater risk of developing DVT (OR = 1.0007, 95%CI [1.0001-1.0013], p = 0.0174) and VTE (OR = 1.0008, 95%CI [1.0002-1.0013], p = 0.0123). Genetically predicted hyperthyroidism was significantly associated with an increased risk of developing DVT (OR = 1.0685, 95%CI [1.0139-1.1261], p = 0.0134) and VTE (OR = 1.0740, 95%CI [1.0165-1.1348], p = 0.0110). According to the initial MR analysis, testosterone concentrations were positively associated with the risk of developing VTE (OR = 1.0038, 95%CI [1.004-1.0072], p = 0.0285). After sex stratification, estradiol concentrations were positively associated with the risk of developing DVT (OR = 1.0143, 95%CI [1.0020-1.0267], p = 0.0226) and VTE (OR = 1.0156, 95%CI [1.0029-1.0285], p = 0.0158) in females, while the significant relationship between testosterone and VTE did not persist. SHBG rs858518 was identified as the only SNP that was associated with an increased risk of developing VTE, mediated by estradiol, in females.

CONCLUSIONS

Genetically predicted hyperthyroidism and increased FT4 concentrations were positively associated with the risk of developing VTE. The effects of genetically predicted sex hormones on the risk of developing VTE differed between males and females. Greater genetically predicted estradiol concentrations were associated with an increased risk of developing VTE in females, while the SHBG rs858518 variant may become a potential prevention and treatment target for female VTE.

Causal relevance of different blood pressure traits on risk of cardiovascular diseases: GWAS and Mendelian randomisation in 100,000 Chinese adults

Elevated blood pressure (BP) is major risk factor for cardiovascular diseases (CVD). Genome-wide association studies (GWAS) conducted predominantly in populations of European ancestry have identified >2,000 BP-associated loci, but other ancestries have been less well-studied. We conducted GWAS of systolic, diastolic, pulse, and mean arterial BP in 100,453 Chinese adults. We identified 128 non-overlapping loci associated with one or more BP traits, including 74 newly-reported associations. Despite strong genetic correlations between populations, we identified appreciably higher heritability and larger variant effect sizes in Chinese compared with European or Japanese ancestry populations. Using instruments derived from these GWAS, multivariable Mendelian randomisation demonstrated that BP traits contribute differently to the causal associations of BP with CVD. In particular, only pulse pressure was independently causally associated with carotid plaque. These findings reinforce the need for studies in diverse populations to understand the genetic determinants of BP traits and their roles in disease risk.

Gaining new insights into the etiology of ulcerative colitis through a cross-tissue transcriptome-wide association study

Background

Genome-wide association studies (GWASs) have identified 38 loci associated with ulcerative colitis (UC) susceptibility, but the risk genes and their biological mechanisms remained to be comprehensively elucidated.

Methods

Multi-marker analysis of genomic annotation (MAGMA) software was used to annotate genes on GWAS summary statistics of UC from FinnGen database. Genetic analysis was performed to identify risk genes. Cross-tissue transcriptome-wide association study (TWAS) using the unified test for molecular signatures (UTMOST) was performed to compare GWAS summary statistics with gene expression matrix (from Genotype-Tissue Expression Project) for data integration. Subsequently, we used FUSION software to select key genes from the individual tissues. Additionally, conditional and joint analysis was conducted to improve our understanding on UC. Fine-mapping of causal gene sets (FOCUS) software was employed to accurately locate risk genes. The results of the four genetic analyses (MAGMA, UTMOST, FUSION and FOCUS) were combined to obtain a set of UC risk genes. Finally, Mendelian randomization (MR) analysis and Bayesian colocalization analysis were conducted to determine the causal relationship between the risk genes and UC. To test the robustness of our findings, the same approaches were taken to verify the GWAS data of UC on IEU.

Results

Multiple correction tests screened PIM3 as a risk gene for UC. The results of Bayesian colocalization analysis showed that the posterior probability of hypothesis 4 was 0.997 and 0.954 in the validation dataset. MR was conducted using the inverse variance weighting method and two single nucleotide polymorphisms (SNPs, rs28645887 and rs62231924) were included in the analysis (p < 0.001, 95%CI: 1.45-1.89). In the validation dataset, MR result was p < 0.001, 95%CI: 1.19-1.72, indicating a clear causal relationship between PIM3 and UC.

Conclusion

Our study validated PIM3 as a key risk gene for UC and its expression level may be related to the risk of UC, providing a novel reference for further improving the current understanding on the genetic structure of UC.

Transcriptome-wide association analysis identifies candidate susceptibility genes for prostate-specific antigen levels in men without prostate cancer

Deciphering the genetic basis of prostate-specific antigen (PSA) levels may improve their utility for prostate cancer (PCa) screening. Using genome-wide association study (GWAS) summary statistics from 95,768 PCa-free men, we conducted a transcriptome-wide association study (TWAS) to examine impacts of genetically predicted gene expression on PSA. Analyses identified 41 statistically significant (p < 0.05/12,192 = 4.10 × 10-6) associations in whole blood and 39 statistically significant (p < 0.05/13,844 = 3.61 × 10-6) associations in prostate tissue, with 18 genes associated in both tissues. Cross-tissue analyses identified 155 statistically significantly (p < 0.05/22,249 = 2.25 × 10-6) genes. Out of 173 unique PSA-associated genes across analyses, we replicated 151 (87.3%) in a TWAS of 209,318 PCa-free individuals from the Million Veteran Program. Based on conditional analyses, we found 20 genes (11 single tissue, nine cross-tissue) that were associated with PSA levels in the discovery TWAS that were not attributable to a lead variant from a GWAS. Ten of these 20 genes replicated, and two of the replicated genes had colocalization probability of >0.5: CCNA2 and HIST1H2BN. Six of the 20 identified genes are not known to impact PCa risk. Fine-mapping based on whole blood and prostate tissue revealed five protein-coding genes with evidence of causal relationships with PSA levels. Of these five genes, four exhibited evidence of colocalization and one was conditionally independent of previous GWAS findings. These results yield hypotheses that should be further explored to improve understanding of genetic factors underlying PSA levels.

A robust cis-Mendelian randomization method with application to drug target discovery

Mendelian randomization (MR) uses genetic variants as instrumental variables (IVs) to investigate causal relationships between traits. Unlike conventional MR, cis-MR focuses on a single genomic region using only cis-SNPs. For example, using cis-pQTLs for a protein as exposure for a disease opens a cost-effective path for drug target discovery. However, few methods effectively handle pleiotropy and linkage disequilibrium (LD) of cis-SNPs. Here, we propose cisMR-cML, a method based on constrained maximum likelihood, robust to IV assumption violations with strong theoretical support. We further clarify the severe but largely neglected consequences of the current practice of modeling marginal, instead of conditional genetic effects, and only using exposure-associated SNPs in cis-MR analysis. Numerical studies demonstrated our method's superiority over other existing methods. In a drug-target analysis for coronary artery disease (CAD), including a proteome-wide application, we identified three potential drug targets, PCSK9, COLEC11 and FGFR1 for CAD.

Integrating muti-omics data to identify tissue-specific DNA methylation biomarkers for cancer risk

The relationship between tissue-specific DNA methylation and cancer risk remains inadequately elucidated. Leveraging resources from the Genotype-Tissue Expression consortium, here we develop genetic models to predict DNA methylation at CpG sites across the genome for seven tissues and apply these models to genome-wide association study data of corresponding cancers, namely breast, colorectal, renal cell, lung, ovarian, prostate, and testicular germ cell cancers. At Bonferroni-corrected P < 0.05, we identify 4248 CpGs that are significantly associated with cancer risk, of which 95.4% (4052) are specific to a particular cancer type. Notably, 92 CpGs within 55 putative novel loci retain significant associations with cancer risk after conditioning on proximal signals identified by genome-wide association studies. Integrative multi-omics analyses reveal 854 CpG-gene-cancer trios, suggesting that DNA methylation at 309 distinct CpGs might influence cancer risk through regulating the expression of 205 unique cis-genes. These findings substantially advance our understanding of the interplay between genetics, epigenetics, and gene expression in cancer etiology.

High incidence and geographic distribution of cleft palate cases in Finland are associated with a regulatory variant in IRF6

In Finland the frequency of isolated cleft palate (CP) is higher than that of isolated cleft lip with or without cleft palate (CL/P). This trend contrasts to that in other European countries but its genetic underpinnings are unknown. We performed a genome-wide association study for orofacial clefts, which include CL/P and CP, in the Finnish population. We identified rs570516915, a single nucleotide polymorphism that is highly enriched in Finns and Estonians, as being strongly associated with CP ( P = 5.25 × 10 -34 , OR = 8.65, 95% CI 6.11-12.25), but not with CL/P ( P = 7.2 × 10 -5 ), with genome-wide significance. The risk allele frequency of rs570516915 parallels the regional variation of CP prevalence in Finland, and the association was replicated in independent cohorts of CP cases from Finland ( P = 8.82 × 10 -28 ) and Estonia ( P = 1.25 × 10 -5 ). The risk allele of rs570516915 disrupts a conserved binding site for the transcription factor IRF6 within a previously characterized enhancer upstream of the IRF6 gene. Through reporter assay experiments we found that the risk allele of rs570516915 diminishes the enhancer activity. Oral epithelial cells derived from CRISPR-Cas9 edited induced pluripotent stem cells demonstrate that the CP-associated allele of rs570516915 concomitantly decreases the binding of IRF6 and the expression level of IRF6 , suggesting impaired IRF6 autoregulation as a molecular mechanism underlying the risk for CP.

Start codon variant in LAG3 is associated with decreased LAG-3 expression and increased risk of autoimmune thyroid disease

Autoimmune thyroid disease (AITD) is a common autoimmune disease. In a GWAS meta-analysis of 110,945 cases and 1,084,290 controls, 290 sequence variants at 225 loci are associated with AITD. Of these variants, 115 are previously unreported. Multiomics analysis yields 235 candidate genes outside the MHC-region and the findings highlight the importance of genes involved in T-cell regulation. A rare 5'-UTR variant (rs781745126-T, MAF = 0.13% in Iceland) in LAG3 has the largest effect (OR = 3.42, P = 2.2 × 10-16) and generates a novel start codon for an open reading frame upstream of the canonical protein translation initiation site. rs781745126-T reduces mRNA and surface expression of the inhibitory immune checkpoint LAG-3 co-receptor on activated lymphocyte subsets and halves LAG-3 levels in plasma among heterozygotes. All three homozygous carriers of rs781745126-T have AITD, of whom one also has two other T-cell mediated diseases, that is vitiligo and type 1 diabetes. rs781745126-T associates nominally with vitiligo (OR = 5.1, P = 6.5 × 10-3) but not with type 1 diabetes. Thus, the effect of rs781745126-T is akin to drugs that inhibit LAG-3, which unleash immune responses and can have thyroid dysfunction and vitiligo as adverse events. This illustrates how a multiomics approach can reveal potential drug targets and safety concerns.

Fine mapping of candidate effector genes for heart rate

An elevated resting heart rate (RHR) is associated with increased cardiovascular mortality. Genome-wide association studies (GWAS) have identified > 350 loci. Uniquely, in this study we applied genetic fine-mapping leveraging tissue specific chromatin segmentation and colocalization analyses to identify causal variants and candidate effector genes for RHR. We used RHR GWAS summary statistics from 388,237 individuals of European ancestry from UK Biobank and performed fine mapping using publicly available genomic annotation datasets. High-confidence causal variants (accounting for > 75% posterior probability) were identified, and we collated candidate effector genes using a multi-omics approach that combined evidence from colocalisation with molecular quantitative trait loci (QTLs), and long-range chromatin interaction analyses. Finally, we performed druggability analyses to investigate drug repurposing opportunities. The fine mapping pipeline indicated 442 distinct RHR signals. For 90 signals, a single variant was identified as a high-confidence causal variant, of which 22 were annotated as missense. In trait-relevant tissues, 39 signals colocalised with cis-expression QTLs (eQTLs), 3 with cis-protein QTLs (pQTLs), and 75 had promoter interactions via Hi-C. In total, 262 candidate genes were highlighted (79% had promoter interactions, 15% had a colocalised eQTL, 8% had a missense variant and 1% had a colocalised pQTL), and, for the first time, enrichment in nervous system pathways. Druggability analyses highlighted ACHE, CALCRL, MYT1 and TDP1 as potential targets. Our genetic fine-mapping pipeline prioritised 262 candidate genes for RHR that warrant further investigation in functional studies, and we provide potential therapeutic targets to reduce RHR and cardiovascular mortality.

Mendelian randomization analyses clarify the effects of height on cardiovascular diseases

An inverse correlation between stature and risk of coronary artery disease (CAD) has been observed in several epidemiologic studies, and recent Mendelian randomization (MR) experiments have suggested causal association. However, the extent to which the effect estimated by MR can be explained by cardiovascular, anthropometric, lung function, and lifestyle-related risk factors is unclear, with a recent report suggesting that lung function traits could fully explain the height-CAD effect. To clarify this relationship, we utilized a well-powered set of genetic instruments for human stature, comprising >1,800 genetic variants for height and CAD. In univariable analysis, we confirmed that a one standard deviation decrease in height (~6.5 cm) was associated with a 12.0% increase in the risk of CAD, consistent with previous reports. In multivariable analysis accounting for effects from up to 12 established risk factors, we observed a >3-fold attenuation in the causal effect of height on CAD susceptibility (3.7%, p = 0.02). However, multivariable analyses demonstrated independent effects of height on other cardiovascular traits beyond CAD, consistent with epidemiologic associations and univariable MR experiments. In contrast with published reports, we observed minimal effects of lung function traits on CAD risk in our analyses, indicating that these traits are unlikely to explain the residual association between height and CAD risk. In sum, these results suggest the impact of height on CAD risk beyond previously established cardiovascular risk factors is minimal and not explained by lung function measures.

Exome sequencing identifies novel genetic variants associated with varicose veins

BACKGROUND

Varicose veins (VV) are one of the common human diseases, but the role of genetics in its development is not fully understood.

METHODS

We conducted an exome-wide association study of VV using whole-exome sequencing data from the UK Biobank, and focused on common and rare variants using single-variant association analysis and gene-level collapsing analysis.

FINDINGS

A total of 13,823,269 autosomal genetic variants were obtained after quality control. We identified 36 VV-related independent common variants mapping to 34 genes by single-variant analysis and three rare variant genes (PIEZO1, ECE1, FBLN7) by collapsing analysis, and most associations between genes and VV were replicated in FinnGen. PIEZO1 was the closest gene associated with VV (P = 5.05 × 10-31), and it was found to reach exome-wide significance in both single-variant and collapsing analyses. Two novel rare variant genes (ECE1 and METTL21A) associated with VV were identified, of which METTL21A was associated only with females. The pleiotropic effects of VV-related genes suggested that body size, inflammation, and pulmonary function are strongly associated with the development of VV.

CONCLUSIONS

Our findings highlight the importance of causal genes for VV and provide new directions for treatment.

Anthropometrics and cancer prognosis: a multicenter cohort study

BACKGROUND

Anthropometric indicators have been shown to be associated with the prognosis of patients with cancer. However, any single anthropometric index has limitation in predicting the prognosis.

OBJECTIVES

This study aimed to observe the predictive role of 7 anthropometric indicators based on body size on the prognosis of patients with cancer.

METHODS

A principal component analysis (PCA) on 7 anthropometric measurements: height, weight, BMI, hand grip strength (HGS), triceps skinfold thickness (TSF), mid-upper arm circumference (MAC), and calf circumference (CAC) was conducted. Principal components (PCs) were derived from this analysis. Cox regression analysis was used to investigate the association between the prognosis of patients with cancer and the PCs. Subgroups and sensitivity analyses were also conducted.

RESULTS

Through PCA, 4 distinct PCs were identified, collectively explaining 88.3% of the variance. PC1, primarily characterized by general obesity, exhibited a significant inverse association with risk of cancer-related death (adjusted hazard ratio [HR]: 0.86; 95% confidence interval [CI]: 0.83, 0.88). PC2 (short stature with high TSF) was not significantly associated with cancer prognosis. PC3 (high BMI coupled with low HGS) demonstrated a significant increase with risk of cancer-related death (adjusted HR: 1.08; 95% CI: 1.05, 1.11). PC4 (tall stature with high TSF) exhibited a significant association with increased cancer risk (adjusted HR: 1.05; 95% CI: 1.02, 1.07). These associations varied across different cancer stages. The stability of the results was confirmed through sensitivity analyses.

CONCLUSIONS

Different body sizes are associated with distinct prognostic outcomes in patients with cancer. The impact of BMI on prognosis is influenced by both HGS and subcutaneous fat. This finding may influence the clinical care of cancer and improve the survival of cancer patients.

Multi-tissue transcriptome-wide association studies identified 235 genes for intrinsic subtypes of breast cancer

BACKGROUND

Although genome-wide association studies (GWAS) of breast cancer (BC) identified common variants which differ between intrinsic subtypes, genes through which these variants act to impact BC risk have not been fully established. Transcriptome-wide association studies (TWAS) have identified genes associated with overall BC risk, but subtype-specific differences are largely unknown.

METHODS

We performed two multi-tissue TWAS for each BC intrinsic subtype, including an expression-based approach that collated TWAS signals from expression quantitative trait loci (eQTLs) across multiple tissues and a novel splicing-based approach that collated signals from splicing QTLs (sQTLs) across intron clusters and subsequently across tissues. We used summary statistics for five intrinsic subtypes including Luminal A-like, Luminal B-like, Luminal B/HER2-negative-like, HER2-enriched-like, and triple-negative BC, generated from 106 278 BC cases and 91 477 controls in the Breast Cancer Association Consortium.

RESULTS

Overall, we identified 235 genes in 88 loci that were associated with at least one of the five intrinsic subtypes. Most genes were subtype-specific, and many have not been reported in previous TWAS. We discovered common variants that modulate expression of CHEK2 confer increased risk to Luminal A-like BC, in contrast to the viewpoint that CHEK2 primarily harbors rare, penetrant mutations. Additionally, our splicing-based TWAS provided population-level support for MDM4 splice variants that increased the risk of triple-negative BC.

CONCLUSION

Our comprehensive, multi-tissue TWAS corroborated previous GWAS loci for overall BC risk and intrinsic subtypes, while underscoring how common variation that impacts expression and splicing of genes in multiple tissue types can be used to further elucidate the etiology of BC.

Genetic drivers and cellular selection of female mosaic X chromosome loss

Mosaic loss of the X chromosome (mLOX) is the most common clonal somatic alteration in leukocytes of female individuals1,2, but little is known about its genetic determinants or phenotypic consequences. Here, to address this, we used data from 883,574 female participants across 8 biobanks; 12% of participants exhibited detectable mLOX in approximately 2% of leukocytes. Female participants with mLOX had an increased risk of myeloid and lymphoid leukaemias. Genetic analyses identified 56 common variants associated with mLOX, implicating genes with roles in chromosomal missegregation, cancer predisposition and autoimmune diseases. Exome-sequence analyses identified rare missense variants in FBXO10 that confer a twofold increased risk of mLOX. Only a small fraction of associations was shared with mosaic Y chromosome loss, suggesting that distinct biological processes drive formation and clonal expansion of sex chromosome missegregation. Allelic shift analyses identified X chromosome alleles that are preferentially retained in mLOX, demonstrating variation at many loci under cellular selection. A polygenic score including 44 allelic shift loci correctly inferred the retained X chromosomes in 80.7% of mLOX cases in the top decile. Our results support a model in which germline variants predispose female individuals to acquiring mLOX, with the allelic content of the X chromosome possibly shaping the magnitude of clonal expansion.

Proteome-wide association study using cis and trans variants and applied to blood cell and lipid-related traits in the Women's Health Initiative study

In most Proteome-Wide Association Studies (PWAS), variants near the protein-coding gene (±1 Mb), also known as cis single nucleotide polymorphisms (SNPs), are used to predict protein levels, which are then tested for association with phenotypes. However, proteins can be regulated through variants outside of the cis region. An intermediate GWAS step to identify protein quantitative trait loci (pQTL) allows for the inclusion of trans SNPs outside the cis region in protein-level prediction models. Here, we assess the prediction of 540 proteins in 1002 individuals from the Women's Health Initiative (WHI), split equally into a GWAS set, an elastic net training set, and a testing set. We compared the testing r2 between measured and predicted protein levels using this proposed approach, to the testing r2 using only cis SNPs. The two methods usually resulted in similar testing r2, but some proteins showed a significant increase in testing r2 with our method. For example, for cartilage acidic protein 1, the testing r2 increased from 0.101 to 0.351. We also demonstrate reproducible findings for predicted protein association with lipid and blood cell traits in WHI participants without proteomics data and in UK Biobank utilizing our PWAS weights.

Adaptive selection at G6PD and disparities in diabetes complications

Diabetes complications occur at higher rates in individuals of African ancestry. Glucose-6-phosphate dehydrogenase deficiency (G6PDdef), common in some African populations, confers malaria resistance, and reduces hemoglobin A1c (HbA1c) levels by shortening erythrocyte lifespan. In a combined-ancestry genome-wide association study of diabetic retinopathy, we identified nine loci including a G6PDdef causal variant, rs1050828 -T (Val98Met), which was also associated with increased risk of other diabetes complications. The effect of rs1050828 -T on retinopathy was fully mediated by glucose levels. In the years preceding diabetes diagnosis and insulin prescription, glucose levels were significantly higher and HbA1c significantly lower in those with versus without G6PDdef. In the Action to Control Cardiovascular Risk in Diabetes (ACCORD) trial, participants with G6PDdef had significantly higher hazards of incident retinopathy and neuropathy. At the same HbA1c levels, G6PDdef participants in both ACCORD and the Million Veteran Program had significantly increased risk of retinopathy. We estimate that 12% and 9% of diabetic retinopathy and neuropathy cases, respectively, in participants of African ancestry are due to this exposure. Across continentally defined ancestral populations, the differences in frequency of rs1050828 -T and other G6PDdef alleles contribute to disparities in diabetes complications. Diabetes management guided by glucose or potentially genotype-adjusted HbA1c levels could lead to more timely diagnoses and appropriate intensification of therapy, decreasing the risk of diabetes complications in patients with G6PDdef alleles.

Genome-wide study of half a million individuals with major depression identifies 697 independent associations, infers causal neuronal subtypes and biological targets for novel pharmacotherapies

In a genome-wide association study (GWAS) of 685,808 individuals with major depression (MD) and 4,364,225 controls from 29 countries and across diverse and admixed ancestries, we identify 697 independent associations at 636 genetic loci, 293 of which are novel. Using fine-mapping and functional genomic datasets, we find 308 high-confidence gene associations and enrichment of postsynaptic density and receptor clustering. Leveraging new single-cell gene expression data, we conducted a causal neural cell type enrichment analysis that implicated excitatory and inhibitory midbrain and forebrain neurons, peptidergic neurons, and medium spiny neurons in MD. Critically, our findings are enriched for the targets of antidepressants and provide potential antidepressant repurposing opportunities (e.g., pregabalin and modafinil). Polygenic scores (PGS) from European ancestries explained up to 5.7% of the variance in liability to MD in European samples and PGS trained using either European or multi-ancestry data significantly predicted case control status across all four diverse ancestries. These findings represent a major advance in our understanding of MD across global populations. We provide evidence that MD GWAS reveals known and novel biological targets that may be used to target and develop pharmacotherapies addressing the considerable unmet need for effective treatment.

Exploring the causal effects of sleep characteristics on TMD-related pain: a two-sample Mendelian randomization study

OBJECTIVES

The study was to explore the causal effects of sleep characteristics on temporomandibular disorder (TMD)-related pain using Mendelian randomization (MR) analysis.

MATERIALS AND METHODS

Five sleep characteristics (short sleep, insomnia, chronotype, snoring, sleep apnea) were designated as exposure factors. Data were obtained from previous publicized genome-wide association studies and single nucleotide polymorphisms (SNPs) strongly associated with them were utilized as instrumental variables (IVs). TMD-related pain was designed as outcome variable and sourced from the FinnGens database. MR analysis was employed to explore the causal effects of the five sleep characteristics on TMD-related pain. The causal effect was analyzed using the inverse variance-weighted (IVW), weighted median, and MR-Egger methods. Subsequently, sensitivity analyses were conducted using Cochran's Q tests, funnel plots, leave-one-out analyses, and MR-Egger intercept tests.

RESULTS

A causal effect of short sleep on TMD-related pain was revealed by IVW (OR: 1.60, 95% CI: 1.06-2.41, P = 0.026). No causal relationship was identified between other sleep characteristics (insomnia, chronotype, snoring, sleep apnea) and TMD-related pain.

CONCLUSIONS

Our study suggests that short sleep may increase the risk of TMD-related pain, while there was no causal relationship between other sleep characteristics and TMD-related pain. Further studies are warranted to deepen and definitively clarify their relationship.

CLINICAL RELEVANCE

These findings reveal that the short sleep may be a risk factor of TMD-related pain and highlight the potential therapeutical effect of extending sleep time on alleviating TMD-related pain.

Fast and flexible joint fine-mapping of multiple traits via the Sum of Single Effects model

We introduce mvSuSiE, a multi-trait fine-mapping method for identifying putative causal variants from genetic association data (individual-level or summary data). mvSuSiE learns patterns of shared genetic effects from data, and exploits these patterns to improve power to identify causal SNPs. Comparisons on simulated data show that mvSuSiE is competitive in speed, power and precision with existing multi-trait methods, and uniformly improves on single-trait fine-mapping (SuSiE) in each trait separately. We applied mvSuSiE to jointly fine-map 16 blood cell traits using data from the UK Biobank. By jointly analyzing the traits and modeling heterogeneous effect sharing patterns, we discovered a much larger number of causal SNPs (>3,000) compared with single-trait fine-mapping, and with narrower credible sets. mvSuSiE also more comprehensively characterized the ways in which the genetic variants affect one or more blood cell traits; 68% of causal SNPs showed significant effects in more than one blood cell type.

The causal relationship between autoimmune diseases and age-related macular degeneration: A two-sample mendelian randomization study

OBJECTIVE

The aim of this study is to investigate the potential causal relationship between autoimmune diseases, including systemic lupus erythematosus, rheumatoid arthritis, inflammatory bowel disease, multiple sclerosis, and Type 1 diabetes, and age-related macular degeneration (AMD). By utilizing the two-sample Mendelian Randomization (MR) approach, we endeavor to address this complex medical issue.

METHODS

Genome-wide association study (GWAS) data for autoimmune diseases and AMD were obtained from the IEU Open GWAS database and the FinnGen consortium. A series of stringent SNP filtering steps was applied to ensure the reliability of the genetic instruments. MR analyses were conducted using the TwoSampleMR and MR-PRESSO packages in R. The inverse-variance weighted (IVW) method served as the primary analysis, complemented by multiple supplementary analyses and sensitivity tests.

RESULTS

Within the discovery sample, only a statistically significant inverse causal relationship between multiple sclerosis (MS) and AMD was observed (OR = 0.92, 95% CI: 0.88-0.97, P = 0.003). This finding was confirmed in the replication sample (OR = 0.85, 95% CI: 0.80-0.89, P = 3.32×10-12). No statistically significant associations were detected between systemic lupus erythematosus, rheumatoid arthritis, inflammatory bowel disease, and Type 1 diabetes and AMD.

CONCLUSION

Strong evidence is provided by this study to support the existence of an inverse causal relationship between multiple sclerosis and age-related macular degeneration. However, no causal evidence was found linking other autoimmune diseases with AMD. These findings not only offer novel insights into the potential etiological mechanisms underlying AMD but also suggest possible directions for future clinical interventions.

Mapping and annotating genomic loci to prioritize genes and implicate distinct polygenic adaptations for skin color

Evidence for adaptation of human skin color to regional ultraviolet radiation suggests shared and distinct genetic variants across populations. However, skin color evolution and genetics in East Asians are understudied. We quantified skin color in 48,433 East Asians using image analysis and identified associated genetic variants and potential causal genes for skin color as well as their polygenic interplay with sun exposure. This genome-wide association study (GWAS) identified 12 known and 11 previously unreported loci and SNP-based heritability was 23-24%. Potential causal genes were determined through the identification of nonsynonymous variants, colocalization with gene expression in skin tissues, and expression levels in melanocytes. Genomic loci associated with pigmentation in East Asians substantially diverged from European populations, and we detected signatures of polygenic adaptation. This large GWAS for objectively quantified skin color in an East Asian population improves understanding of the genetic architecture and polygenic adaptation of skin color and prioritizes potential causal genes.

Genome-wide determinants of mortality and motor progression in Parkinson's disease

There are 90 independent genome-wide significant genetic risk variants for Parkinson's disease (PD) but currently only five nominated loci for PD progression. The biology of PD progression is likely to be of central importance in defining mechanisms that can be used to develop new treatments. We studied 6766 PD patients, over 15,340 visits with a mean follow-up of between 4.2 and 15.7 years and carried out genome-wide survival studies for time to a motor progression endpoint, defined by reaching Hoehn and Yahr stage 3 or greater, and death (mortality). There was a robust effect of the APOE ε4 allele on mortality in PD. We also identified a locus within the TBXAS1 gene encoding thromboxane A synthase 1 associated with mortality in PD. We also report 4 independent loci associated with motor progression in or near MORN1, ASNS, PDE5A, and XPO1. Only the non-Gaucher disease causing GBA1 PD risk variant E326K, of the known PD risk variants, was associated with mortality in PD. Further work is needed to understand the links between these genomic variants and the underlying disease biology. However, these may represent new candidates for disease modification in PD.

Genome-wide association study of glucocerebrosidase activity modifiers

One of the most common genetic risk factors for Parkinson's disease (PD) are variants in GBA1, which encodes the lysosomal enzyme glucocerebrosidase (GCase). GCase deficiency has been associated with an increased PD risk, but not all individuals with low GCase activity are carriers of GBA1 mutations, suggesting other factors may be acting as modifiers. We aimed to discover common variants associated with GCase activity, as well as replicate previously reported associations, by performing a genome-wide association study using two independent cohorts: a Columbia University cohort consisting of 697 PD cases and 347 controls and the Parkinson's Progression Markers Initiative (PPMI) cohort consisting of 357 PD cases and 163 controls. As expected, GBA1 variants have the strongest association with decreased activity, led by p.N370S (beta = -4.36, se = 0.32, p = 5.05e-43). We also identify a novel association in the GAA locus (encoding for acid alpha-glucosidase, beta = -0.96, se = 0.17, p = 5.23e-09) that may be the result of an interaction between GCase and acid alpha-glucosidase based on various interaction analyses. Lastly, we show that several PD-risk loci are potentially associated with GCase activity. Further research will be needed to replicate and validate our findings and to uncover the functional connection between acid alpha-glucosidase and GCase.

A multi-tissue, splicing-based joint transcriptome-wide association study identifies susceptibility genes for breast cancer

Splicing-based transcriptome-wide association studies (splicing-TWASs) of breast cancer have the potential to identify susceptibility genes. However, existing splicing-TWASs test the association of individual excised introns in breast tissue only and thus have limited power to detect susceptibility genes. In this study, we performed a multi-tissue joint splicing-TWAS that integrated splicing-TWAS signals of multiple excised introns in each gene across 11 tissues that are potentially relevant to breast cancer risk. We utilized summary statistics from a meta-analysis that combined genome-wide association study (GWAS) results of 424,650 women of European ancestry. Splicing-level prediction models were trained in GTEx (v.8) data. We identified 240 genes by the multi-tissue joint splicing-TWAS at the Bonferroni-corrected significance level; in the tissue-specific splicing-TWAS that combined TWAS signals of excised introns in genes in breast tissue only, we identified nine additional significant genes. Of these 249 genes, 88 genes in 62 loci have not been reported by previous TWASs, and 17 genes in seven loci are at least 1 Mb away from published GWAS index variants. By comparing the results of our splicing-TWASs with previous gene-expression-based TWASs that used the same summary statistics and expression prediction models trained in the same reference panel, we found that 110 genes in 70 loci that are identified only by the splicing-TWASs. Our results showed that for many genes, expression quantitative trait loci (eQTL) did not show a significant impact on breast cancer risk, whereas splicing quantitative trait loci (sQTL) showed a strong impact through intron excision events.

Meta-Analysis of Genome-Wide Association Studies Reveals Genetic Mechanisms of Supraventricular Arrhythmias

BACKGROUND

Substantial data support a heritable basis for supraventricular tachycardias, but the genetic determinants and molecular mechanisms of these arrhythmias are poorly understood. We sought to identify genetic loci associated with atrioventricular nodal reentrant tachycardia (AVNRT) and atrioventricular accessory pathways or atrioventricular reciprocating tachycardia (AVAPs/AVRT).

METHODS

We performed multiancestry meta-analyses of genome-wide association studies to identify genetic loci for AVNRT (4 studies) and AVAP/AVRT (7 studies). We assessed evidence supporting the potential causal effects of candidate genes by analyzing relations between associated variants and cardiac gene expression, performing transcriptome-wide analyses, and examining prior genome-wide association studies.

RESULTS

Analyses comprised 2384 AVNRT cases and 106 489 referents, and 2811 AVAP/AVRT cases and 1,483 093 referents. We identified 2 significant loci for AVNRT, which implicate NKX2-5 and TTN as disease susceptibility genes. A transcriptome-wide association analysis supported an association between reduced predicted cardiac expression of NKX2-5 and AVNRT. We identified 3 significant loci for AVAP/AVRT, which implicate SCN5A, SCN10A, and TTN/CCDC141. Variant associations at several loci have been previously reported for cardiac phenotypes, including atrial fibrillation, stroke, Brugada syndrome, and electrocardiographic intervals.

CONCLUSIONS

Our findings highlight gene regions associated with ion channel function (AVAP/AVRT), as well as cardiac development and the sarcomere (AVAP/AVRT and AVNRT) as important potential effectors of supraventricular tachycardia susceptibility.

Causal effects of endometriosis on cancer risk: A Mendelian randomization study

Endometriosis has been reported in epidemiological studies to be associated with certain types of cancer. However, the presence of reverse causality and residual confounding due to common risk factors introduces uncertainty regarding the extent to which endometriosis itself contributes to the development of cancer. We performed the Mendelian randomization (MR) to investigate the causal associations between endometriosis and 34 different types of cancers. The results of the inverse-variance-weighted (IVW) model suggested that genetic predisposition to endometriosis was causally associated with an increased risk for ovarian cancer (OR = 3.2913; p-value = .0320). The genetic liabilities to endometriosis had causal associations with the decreased risk for skin cancer (OR = 0.9973; p-value = .0219), hematological cancer (OR = 0.9953; p-value = .0175) and ER- breast cancer (OR = 0.6960; p-value = .0381). The causal association of the above combinations were robust by test of heterogeneity and pleiotropy. Together, our study suggests that endometriosis had causal effect on cancer risk.

Cross-ancestry genome-wide association studies of brain imaging phenotypes

Genome-wide association studies of brain imaging phenotypes are mainly performed in European populations, but other populations are severely under-represented. Here, we conducted Chinese-alone and cross-ancestry genome-wide association studies of 3,414 brain imaging phenotypes in 7,058 Chinese Han and 33,224 white British participants. We identified 38 new associations in Chinese-alone analyses and 486 additional new associations in cross-ancestry meta-analyses at P < 1.46 × 10-11 for discovery and P < 0.05 for replication. We pooled significant autosomal associations identified by single- or cross-ancestry analyses into 6,443 independent associations, which showed uneven distribution in the genome and the phenotype subgroups. We further divided them into 44 associations with different effect sizes and 3,557 associations with similar effect sizes between ancestries. Loci of these associations were shared with 15 brain-related non-imaging traits including cognition and neuropsychiatric disorders. Our results provide a valuable catalog of genetic associations for brain imaging phenotypes in more diverse populations.

Robust use of phenotypic heterogeneity at drug target genes for mechanistic insights: Application of cis-multivariable Mendelian randomization to GLP1R gene region

Phenotypic heterogeneity at genomic loci encoding drug targets can be exploited by multivariable Mendelian randomization to provide insight into the pathways by which pharmacological interventions may affect disease risk. However, statistical inference in such investigations may be poor if overdispersion heterogeneity in measured genetic associations is unaccounted for. In this work, we first develop conditional F statistics for dimension-reduced genetic associations that enable more accurate measurement of phenotypic heterogeneity. We then develop a novel extension for two-sample multivariable Mendelian randomization that accounts for overdispersion heterogeneity in dimension-reduced genetic associations. Our empirical focus is to use genetic variants in the GLP1R gene region to understand the mechanism by which GLP1R agonism affects coronary artery disease (CAD) risk. Colocalization analyses indicate that distinct variants in the GLP1R gene region are associated with body mass index and type 2 diabetes (T2D). Multivariable Mendelian randomization analyses that were corrected for overdispersion heterogeneity suggest that bodyweight lowering rather than T2D liability lowering effects of GLP1R agonism are more likely contributing to reduced CAD risk. Tissue-specific analyses prioritized brain tissue as the most likely to be relevant for CAD risk, of the tissues considered. We hope the multivariable Mendelian randomization approach illustrated here is widely applicable to better understand mechanisms linking drug targets to diseases outcomes, and hence to guide drug development efforts.

Dynamics of cognitive variability with age and its genetic underpinning in NIHR BioResource Genes and Cognition cohort participants

A leading explanation for translational failure in neurodegenerative disease is that new drugs are evaluated late in the disease course when clinical features have become irreversible. Here, to address this gap, we cognitively profiled 21,051 people aged 17-85 years as part of the Genes and Cognition cohort within the National Institute for Health and Care Research BioResource across England. We describe the cohort, present cognitive trajectories and show the potential utility. Surprisingly, when studied at scale, the APOE genotype had negligible impact on cognitive performance. Different cognitive domains had distinct genetic architectures, with one indicating brain region-specific activation of microglia and another with glycogen metabolism. Thus, the molecular and cellular mechanisms underpinning cognition are distinct from dementia risk loci, presenting different targets to slow down age-related cognitive decline. Participants can now be recalled stratified by genotype and cognitive phenotype for natural history and interventional studies of neurodegenerative and other disorders.

Genome-wide meta-analyses of restless legs syndrome yield insights into genetic architecture, disease biology and risk prediction

Restless legs syndrome (RLS) affects up to 10% of older adults. Their healthcare is impeded by delayed diagnosis and insufficient treatment. To advance disease prediction and find new entry points for therapy, we performed meta-analyses of genome-wide association studies in 116,647 individuals with RLS (cases) and 1,546,466 controls of European ancestry. The pooled analysis increased the number of risk loci eightfold to 164, including three on chromosome X. Sex-specific meta-analyses revealed largely overlapping genetic predispositions of the sexes (rg = 0.96). Locus annotation prioritized druggable genes such as glutamate receptors 1 and 4, and Mendelian randomization indicated RLS as a causal risk factor for diabetes. Machine learning approaches combining genetic and nongenetic information performed best in risk prediction (area under the curve (AUC) = 0.82-0.91). In summary, we identified targets for drug development and repurposing, prioritized potential causal relationships between RLS and relevant comorbidities and risk factors for follow-up and provided evidence that nonlinear interactions are likely relevant to RLS risk prediction.

Assessment of causal relationships between omega-3 and omega-6 polyunsaturated fatty acids in autoimmune rheumatic diseases: a brief research report from a Mendelian randomization study

Background

Currently, the association between the consumption of polyunsaturated fatty acids (PUFAs) and the susceptibility to autoimmune rheumatic diseases (ARDs) remains conflict and lacks substantial evidence in various clinical studies. To address this issue, we employed Mendelian randomization (MR) to establish causal links between six types of PUFAs and their connection to the risk of ARDs.

Methods

We retrieved summary-level data on six types of PUFAs, and five different types of ARDs from publicly accessible GWAS statistics. Causal relationships were determined using a two-sample MR analysis, with the IVW approach serving as the primary analysis method. To ensure the reliability of our research findings, we used four complementary approaches and conducted multivariable MR analysis (MVMR). Additionally, we investigated reverse causality through a reverse MR analysis.

Results

Our results indicate that a heightened genetic predisposition for elevated levels of EPA (ORIVW: 0.924, 95% CI: 0.666-1.283, P IVW = 0.025) was linked to a decreased susceptibility to psoriatic arthritis (PsA). Importantly, the genetically predicted higher levels of EPA remain significantly associated with an reduced risk of PsA, even after adjusting for multiple testing using the FDR method (P IVW-FDR-corrected = 0.033) and multivariable MR analysis (P MV-IVW < 0.05), indicating that EPA may be considered as the risk-protecting PUFAs for PsA. Additionally, high levels of LA showed a positive causal relationship with a higher risk of PsA (ORIVW: 1.248, 95% CI: 1.013-1.538, P IVW = 0.037). It is interesting to note, however, that the effects of these associations were weakened in our MVMR analyses, which incorporated adjustment for lipid profiles (P MV-IVW > 0.05) and multiple testing using the FDR method (P IVW-FDR-corrected = 0.062). Moreover, effects of total omega-3 PUFAs, DHA, EPA, and LA on PsA, were massively driven by SNP effects in the FADS gene region. Furthermore, no causal association was identified between the concentrations of other circulating PUFAs and the risk of other ARDs. Further analysis revealed no significant horizontal pleiotropy and heterogeneity or reverse causality.

Conclusion

Our comprehensive MR analysis indicated that EPA is a key omega-3 PUFA that may protect against PsA but not other ARDs. The FADS2 gene appears to play a central role in mediating the effects of omega-3 PUFAs on PsA risk. These findings suggest that EPA supplementation may be a promising strategy for preventing PsA onset. Further well-powered epidemiological studies and clinical trials are warranted to explore the potential mechanisms underlying the protective effects of EPA in PsA.

A genome-wide meta-analysis of palmoplantar pustulosis implicates TH2 responses and cigarette smoking in disease pathogenesis

BACKGROUND

Palmoplantar pustulosis (PPP) is an inflammatory skin disorder that mostly affects smokers and manifests with painful pustular eruptions on the palms and soles. Although the disease can present with concurrent plaque psoriasis, TNF and IL-17/IL-23 inhibitors show limited efficacy. There is therefore a pressing need to uncover PPP disease drivers and therapeutic targets.

OBJECTIVES

We sought to identify genetic determinants of PPP and investigate whether cigarette smoking contributes to disease pathogenesis.

METHODS

We performed a genome-wide association meta-analysis of 3 North-European cohorts (n = 1,456 PPP cases and 402,050 controls). We then used the scGWAS program to investigate the cell-type specificity of the association signals. We also undertook genetic correlation analyses to examine the similarities between PPP and other immune-mediated diseases. Finally, we applied Mendelian randomization to analyze the causal relationship between cigarette smoking and PPP.

RESULTS

We found that PPP is not associated with the main genetic determinants of plaque psoriasis. Conversely, we identified genome-wide significant associations with the FCGR3A/FCGR3B and CCHCR1 loci. We also observed 13 suggestive (P < 5 × 10-6) susceptibility regions, including the IL4/IL13 interval. Accordingly, we demonstrated a significant genetic correlation between PPP and TH2-mediated diseases such as atopic dermatitis and ulcerative colitis. We also found that genes mapping to PPP-associated intervals were preferentially expressed in dendritic cells and often implicated in T-cell activation pathways. Finally, we undertook a Mendelian randomization analysis, which supported a causal role of cigarette smoking in PPP.

CONCLUSIONS

The first genome-wide association study of PPP points to a pathogenic role for deregulated TH2 responses and cigarette smoking.

Brain cell-type shifts in Alzheimer's disease, autism, and schizophrenia interrogated using methylomics and genetics

Few neuropsychiatric disorders have replicable biomarkers, prompting high-resolution and large-scale molecular studies. However, we still lack consensus on a more foundational question: whether quantitative shifts in cell types-the functional unit of life-contribute to neuropsychiatric disorders. Leveraging advances in human brain single-cell methylomics, we deconvolve seven major cell types using bulk DNA methylation profiling across 1270 postmortem brains, including from individuals diagnosed with Alzheimer's disease, schizophrenia, and autism. We observe and replicate cell-type compositional shifts for Alzheimer's disease (endothelial cell loss), autism (increased microglia), and schizophrenia (decreased oligodendrocytes), and find age- and sex-related changes. Multiple layers of evidence indicate that endothelial cell loss contributes to Alzheimer's disease, with comparable effect size to APOE genotype among older people. Genome-wide association identified five genetic loci related to cell-type composition, involving plausible genes for the neurovascular unit (P2RX5 and TRPV3) and excitatory neurons (DPY30 and MEMO1). These results implicate specific cell-type shifts in the pathophysiology of neuropsychiatric disorders.

Identifying drug targets for schizophrenia through gene prioritization

Background

Schizophrenia genome-wide association studies (GWASes) have identified >250 significant loci and prioritized >100 disease-related genes. However, gene prioritization efforts have mostly been restricted to locus-based methods that ignore information from the rest of the genome.

Methods

To more accurately characterize genes involved in schizophrenia etiology, we applied a combination of highly-predictive tools to a published GWAS of 67,390 schizophrenia cases and 94,015 controls. We combined both locus-based methods (fine-mapped coding variants, distance to GWAS signals) and genome-wide methods (PoPS, MAGMA, ultra-rare coding variant burden tests). To validate our findings, we compared them with previous prioritization efforts, known neurodevelopmental genes, and results from the PsyOPS tool.

Results

We prioritized 62 schizophrenia genes, 41 of which were also highlighted by our validation methods. In addition to DRD2, the principal target of antipsychotics, we prioritized 9 genes that are targeted by approved or investigational drugs. These included drugs targeting glutamatergic receptors (GRIN2A and GRM3), calcium channels (CACNA1C and CACNB2), and GABAB receptor (GABBR2). These also included genes in loci that are shared with an addiction GWAS (e.g. PDE4B and VRK2).

Conclusions

We curated a high-quality list of 62 genes that likely play a role in the development of schizophrenia. Developing or repurposing drugs that target these genes may lead to a new generation of schizophrenia therapies. Rodent models of addiction more closely resemble the human disorder than rodent models of schizophrenia. As such, genes prioritized for both disorders could be explored in rodent addiction models, potentially facilitating drug development.

Genetic association and transferability for urinary albumin-creatinine ratio as a marker of kidney disease in four Sub-Saharan African populations and non-continental individuals of African ancestry

Background

Genome-wide association studies (GWAS) have predominantly focused on populations of European and Asian ancestry, limiting our understanding of genetic factors influencing kidney disease in Sub-Saharan African (SSA) populations. This study presents the largest GWAS for urinary albumin-to-creatinine ratio (UACR) in SSA individuals, including 8,970 participants living in different African regions and an additional 9,705 non-resident individuals of African ancestry from the UK Biobank and African American cohorts.

Methods

Urine biomarkers and genotype data were obtained from two SSA cohorts (AWI-Gen and ARK), and two non-resident African-ancestry studies (UK Biobank and CKD-Gen Consortium). Association testing and meta-analyses were conducted, with subsequent fine-mapping, conditional analyses, and replication studies. Polygenic scores (PGS) were assessed for transferability across populations.

Results

Two genome-wide significant (P < 5 × 10-8) UACR-associated loci were identified, one in the BMP6 region on chromosome 6, in the meta-analysis of resident African individuals, and another in the HBB region on chromosome 11 in the meta-analysis of non-resident SSA individuals, as well as the combined meta-analysis of all studies. Replication of previous significant results confirmed associations in known UACR-associated regions, including THB53, GATM, and ARL15. PGS estimated using previous studies from European ancestry, African ancestry, and multi-ancestry cohorts exhibited limited transferability of PGS across populations, with less than 1% of observed variance explained.

Conclusion

This study contributes novel insights into the genetic architecture of kidney disease in SSA populations, emphasizing the need for conducting genetic research in diverse cohorts. The identified loci provide a foundation for future investigations into the genetic susceptibility to chronic kidney disease in underrepresented African populations Additionally, there is a need to develop integrated scores using multi-omics data and risk factors specific to the African context to improve the accuracy of predicting disease outcomes.

The genetic dissection of fetal haemoglobin persistence in sickle cell disease in Nigeria

The clinical severity of sickle cell disease (SCD) is strongly influenced by the level of fetal haemoglobin (HbF) persistent in each patient. Three major HbF loci (BCL11A, HBS1L-MYB, and Xmn1-HBG2) have been reported, but a considerable hidden heritability remains. We conducted a genome-wide association study for HbF levels in 1006 Nigerian patients with SCD (HbSS/HbSβ0), followed by a replication and meta-analysis exercise in four independent SCD cohorts (3,582 patients). To dissect association signals at the major loci, we performed stepwise conditional and haplotype association analyses and included public functional annotation datasets. Association signals were detected for BCL11A (lead SNP rs6706648, β = -0.39, P = 4.96 × 10-34) and HBS1L-MYB (lead SNP rs61028892, β = 0.73, P = 1.18 × 10-9), whereas the variant allele for Xmn1-HBG2 was found to be very rare. In addition, we detected three putative new trait-associated regions. Genetically, dissecting the two major loci BCL11A and HBS1L-MYB, we defined trait-increasing haplotypes (P < 0.0001) containing so far unidentified causal variants. At BCL11A, in addition to a haplotype harbouring the putative functional variant rs1427407-'T', we identified a second haplotype, tagged by the rs7565301-'A' allele, where a yet-to-be-discovered causal DNA variant may reside. Similarly, at HBS1L-MYB, one HbF-increasing haplotype contains the likely functional small indel rs66650371, and a second tagged by rs61028892-'C' is likely to harbour a presently unknown functional allele. Together, variants at BCL11A and HBS1L-MYB SNPs explained 24.1% of the trait variance. Our findings provide a path for further investigation of the causes of variable fetal haemoglobin persistence in sickle cell disease.

Using genome and transcriptome data from African-ancestry female participants to identify putative breast cancer susceptibility genes

African-ancestry (AA) participants are underrepresented in genetics research. Here, we conducted a transcriptome-wide association study (TWAS) in AA female participants to identify putative breast cancer susceptibility genes. We built genetic models to predict levels of gene expression, exon junction, and 3' UTR alternative polyadenylation using genomic and transcriptomic data generated in normal breast tissues from 150 AA participants and then used these models to perform association analyses using genomic data from 18,034 cases and 22,104 controls. At Bonferroni-corrected P < 0.05, we identified six genes associated with breast cancer risk, including four genes not previously reported (CTD-3080P12.3, EN1, LINC01956 and NUP210L). Most of these genes showed a stronger association with risk of estrogen-receptor (ER) negative or triple-negative than ER-positive breast cancer. We also replicated the associations with 29 genes reported in previous TWAS at P < 0.05 (one-sided), providing further support for an association of these genes with breast cancer risk. Our study sheds new light on the genetic basis of breast cancer and highlights the value of conducting research in AA populations.

SharePro: an accurate and efficient genetic colocalization method accounting for multiple causal signals

MOTIVATION

Colocalization analysis is commonly used to assess whether two or more traits share the same genetic signals identified in genome-wide association studies (GWAS), and is important for prioritizing targets for functional follow-up of GWAS results. Existing colocalization methods can have suboptimal performance when there are multiple causal variants in one genomic locus.

RESULTS

We propose SharePro to extend the COLOC framework for colocalization analysis. SharePro integrates linkage disequilibrium (LD) modeling and colocalization assessment by grouping correlated variants into effect groups. With an efficient variational inference algorithm, posterior colocalization probabilities can be accurately estimated. In simulation studies, SharePro demonstrated increased power with a well-controlled false positive rate at a low computational cost. Compared to existing methods, SharePro provided stronger and more consistent colocalization evidence for known lipid-lowering drug target proteins and their corresponding lipid traits. Through an additional challenging case of the colocalization analysis of the circulating abundance of R-spondin 3 GWAS and estimated bone mineral density GWAS, we demonstrated the utility of SharePro in identifying biologically plausible colocalized signals.

AVAILABILITY AND IMPLEMENTATION

SharePro for colocalization analysis is written in Python and openly available at https://github.com/zhwm/SharePro_coloc.