A gene-based association method for mapping traits using reference transcriptome data

Optimal variable identification for accurate detection of causal expression Quantitative Trait Loci with applications in heart-related diseases

Gene expression plays a pivotal role in various diseases, contributing significantly to their mechanisms. Most GWAS risk loci are in non-coding regions, potentially affecting disease risk by altering gene expression in specific tissues. This expression is notably tissue-specific, with genetic variants substantially influencing it. However, accurately detecting the expression Quantitative Trait Loci (eQTL) is challenging due to limited heritability in gene expression, extensive linkage disequilibrium (LD), and multiple causal variants. The single variant association approach in eQTL analysis is limited by its susceptibility to capture the combined effects of multiple variants, and a bias towards common variants, underscoring the need for a more robust method to accurately identify causal eQTL variants. To address this, we developed an algorithm, CausalEQTL, which integrates L 0 +L 1 penalized regression with an ensemble approach to localize eQTL, thereby enhancing prediction performance precisely. Our results demonstrate that CausalEQTL outperforms traditional models, including LASSO, Elastic Net, Ridge, in terms of power and overall performance. Furthermore, analysis of heart tissue data from the GTEx project revealed that eQTL sites identified by our algorithm provide deeper insights into heart-related tissue eQTL detection. This advancement in eQTL mapping promises to improve our understanding of the genetic basis of tissue-specific gene expression and its implications in disease. The source code and identified causal eQTLs for CausalEQTL are available on GitHub: https://github.com/zhc-moushang/CausalEQTL.

MIMOSA: a resource consisting of improved methylome prediction models increases power to identify DNA methylation-phenotype associations

Although DNA methylation (DNAm) has been implicated in the pathogenesis of numerous complex diseases, from cancer to cardiovascular disease to autoimmune disease, the exact methylation sites that play key roles in these processes remain elusive. One strategy to identify putative causal CpG sites and enhance disease etiology understanding is to conduct methylome-wide association studies (MWASs), in which predicted DNA methylation that is associated with complex diseases can be identified. However, current MWAS models are primarily trained using the data from single studies, thereby limiting the methylation prediction accuracy and the power of subsequent association studies. Here, we introduce a new resource, MWAS Imputing Methylome Obliging Summary-level mQTLs and Associated LD matrices (MIMOSA), a set of models that substantially improve the prediction accuracy of DNA methylation and subsequent MWAS power through the use of a large summary-level mQTL dataset provided by the Genetics of DNA Methylation Consortium (GoDMC). Through the analyses of GWAS (genome-wide association study) summary statistics for 28 complex traits and diseases, we demonstrate that MIMOSA considerably increases the accuracy of DNA methylation prediction in whole blood, crafts fruitful prediction models for low heritability CpG sites, and determines markedly more CpG site-phenotype associations than preceding methods. Finally, we use MIMOSA to conduct a case study on high cholesterol, pinpointing 146 putatively causal CpG sites.

Differential interactions between gene expressions and stressors across the lifespan in major depressive disorder

BACKGROUND

Both genetic predispositions and exposures to stressors have collectively contributed to the development of major depressive disorder (MDD). To deep dive into their roles in MDD, our study aimed to examine which susceptible gene expression interacts with various dimensions of stressors in the MDD risk among a large population cohort.

METHODS

Data analyzed were from a longitudinal community-based cohort from Southwest Montreal, Canada (N = 1083). Latent profile models were used to identify distinct patterns of stressors for the study cohort. A transcriptome-wide association study (TWAS) method was performed to examine the interactive effects of three dimensions of stressors (threat, deprivation, and cumulative lifetime stress) and gene expression on the MDD risk in a total of 48 tissues from GTEx. Additional analyses were also conducted to further explore and specify these associations including colocalization, and fine-mapping analyses, in addition to enrichment analysis investigations based on TWAS.

RESULTS

We identified 3321 genes linked to MDD at the nominal p-value <0.05 and found that different patterns of stressors can amplify the genetic susceptibility to MDD. We also observed specific genes and pathways that interacted with deprivation and cumulative lifetime stressors, particularly in specific brain tissues including basal ganglia, prefrontal cortex, brain amygdala, brain cerebellum, brain cortex, and the whole blood. Colocalization analysis also identified these genes as having a high probability of sharing MDD causal variants.

LIMITATIONS

The study cohort was composed exclusively of individuals of Caucasians, which restricts the generalizability of the findings to other ethnic population groups.

CONCLUSIONS

The findings of the study unveiled significant interactions between potential tissue-specific gene expression × stressors in the MDD risk and shed light on the intricate etiological attributes of gene expression and specific stressors across the lifespan in MDD. These genetic and environmental attributes in MDD corroborate the vulnerability-stress theory and direct future stress research to have a closer examination of genetic predisposition and potential involvements of omics studies to specify the intricate relationships between genes and stressful environments.

Pleiotropy, epistasis and the genetic architecture of quantitative traits

Pleiotropy (whereby one genetic polymorphism affects multiple traits) and epistasis (whereby non-linear interactions between genetic polymorphisms affect the same trait) are fundamental aspects of the genetic architecture of quantitative traits. Recent advances in the ability to characterize the effects of polymorphic variants on molecular and organismal phenotypes in human and model organism populations have revealed the prevalence of pleiotropy and unexpected shared molecular genetic bases among quantitative traits, including diseases. By contrast, epistasis is common between polymorphic loci associated with quantitative traits in model organisms, such that alleles at one locus have different effects in different genetic backgrounds, but is rarely observed for human quantitative traits and common diseases. Here, we review the concepts and recent inferences about pleiotropy and epistasis, and discuss factors that contribute to similarities and differences between the genetic architecture of quantitative traits in model organisms and humans.

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 study identifies novel genes associated with bone toxicities in children with acute lymphoblastic leukaemia

Bone toxicities are common among paediatric patients treated for acute lymphoblastic leukaemia (ALL) with potentially major negative impact on patients' quality of life. To identify the underlying genetic contributors, we conducted a genome-wide association study (GWAS) and a transcriptome-wide association study (TWAS) in 260 patients of European-descent from the DFCI 05-001 ALL trial, with validation in 101 patients of European-descent from the DFCI 11-001 ALL trial. We identified a significant association between rs844882 on chromosome 20 and bone toxicities in the DFCI 05-001 trial (p = 1.7 × 10-8). In DFCI 11-001 trial, we observed a consistent trend of this variant with fracture. The variant was an eQTL for two nearby genes, CD93 and THBD. In TWAS, genetically predicted ACAD9 expression was associated with an increased risk of bone toxicities, which was confirmed by meta-analysis of the two cohorts (meta-p = 2.4 × 10-6). In addition, a polygenic risk score of heel quantitative ultrasound speed of sound was associated with fracture risk in both cohorts (meta-p = 2.3 × 10-3). Our findings highlight the genetic influence on treatment-related bone toxicities in this patient population. The genes we identified in our study provide new biological insights into the development of bone adverse events related to ALL treatment.

Transcriptome- and proteome-wide association studies identify genes associated with renal cell carcinoma

We performed a series of integrative analyses including transcriptome-wide association studies (TWASs) and proteome-wide association studies (PWASs) of renal cell carcinoma (RCC) to nominate and prioritize molecular targets for laboratory investigation. On the basis of a genome-wide association study (GWAS) of 29,020 affected individuals and 835,670 control individuals and prediction models trained in transcriptomic reference models, our TWAS across four kidney transcriptomes (GTEx kidney cortex, kidney tubules, TCGA-KIRC [The Cancer Genome Atlas kidney renal clear-cell carcinoma], and TCGA-KIRP [TCGA kidney renal papillary cell carcinoma]) identified 38 gene associations (false-discovery rate <5%) in at least two of four transcriptomic panels and identified 12 genes that were independent of GWAS susceptibility regions. Analyses combining TWAS associations across 48 tissues from GTEx identified associations that were replicable in tumor transcriptomes for 23 additional genes. Analyses by the two major histologic types (clear-cell RCC and papillary RCC) revealed subtype-specific associations, although at least three gene associations were common to both subtypes. PWAS identified 13 associated proteins, all mapping to GWAS-significant loci. TWAS-identified genes were enriched for active enhancer or promoter regions in RCC tumors and hypoxia-inducible factor binding sites in relevant cell lines. Using gene expression correlation, common cancers (breast and prostate) and RCC risk factors (e.g., hypertension and BMI) display genetic contributions shared with RCC. Our work identifies potential molecular targets for RCC susceptibility for downstream functional investigation.

High-resolution global diversity copy number variation maps and association with ctyper

Genetic analysis of copy number variations (CNVs), especially in complex regions, is challenging due to reference bias and ambiguous alignment of Next-Generation Sequencing (NGS) reads to repetitive DNA. Consequently, aggregate copy numbers are typically analyzed, overlooking variation between gene copies. Pangenomes contain diverse sequences of gene copies and enable the study of sequence-resolved CNVs. We developed a method, ctyper, to discover sequence-resolved CNVs in NGS data by leveraging CNV genes from pangenomes. From 118 public assemblies, we constructed a database of 3,351 CNV genes, distinguishing each gene copy as a resolved allele. We used phylogenetic trees to organize alleles into highly similar allele-types that revealed events of linked small variants due to stratification, structural variation, conversion, and duplication. Saturation analysis showed that new samples share an average of 97.8% CNV alleles with the database. The ctyper method traces individual gene copies in NGS data to their nearest alleles in the database and identifies allele-specific copy numbers using multivariate linear regression on k-mer counts and phylogenetic clustering. Applying ctyper to 1000 Genomes Project (1kgp) samples showed Hardy-Weinberg Equilibrium on 99.3% of alleles and a 97.6% F1 score on genotypes based on 641 1kgp trios. Leave-one-out analysis on 39 assemblies matched to 1kgp samples showed that 96.5% of variants in query sequences match the genotyped allele. Genotyping 1kgp data revealed 226 population-specific CNVs, including a conversion on SMN2 to SMN1, potentially impacting Spinal Muscular Atrophy diagnosis in Africans. Our results revealed two models of CNV: recent CNVs due to ongoing duplications and polymorphic CNVs from ancient paralogs missing from the reference. To measure the functional impact of CNVs, after merging allele-types, we conducted genome-wide Quantitative Trait Locus analysis on 451 1kgp samples with Geuvadis rRNA-seqs. Using a linear mixed model, our genotyping enables the inference of relative expression levels of paralogs within a gene family. In a global evolutionary context, 150 out of 1,890 paralogs (7.94%) and 546 out of 16,628 orthologs (3.28%) had significantly different expression levels, suggesting divergent expression from original genes. Specific examples include lower expression on the converted SMN and increased expression on translocated AMY2B (GTEx pancreas data). Our method enables large cohort studies on complex CNVs to uncover hidden health impacts and overcome reference bias.

Splicing-specific transcriptome-wide association uncovers genetic mechanisms for schizophrenia

Recent studies have highlighted the essential role of RNA splicing, a key mechanism of alternative RNA processing, in establishing connections between genetic variations and disease. Genetic loci influencing RNA splicing variations show considerable influence on complex traits, possibly surpassing those affecting total gene expression. Dysregulated RNA splicing has emerged as a major potential contributor to neurological and psychiatric disorders, likely due to the exceptionally high prevalence of alternatively spliced genes in the human brain. Nevertheless, establishing direct associations between genetically altered splicing and complex traits has remained an enduring challenge. We introduce Spliced-Transcriptome-Wide Associations (SpliTWAS) to integrate alternative splicing information with genome-wide association studies to pinpoint genes linked to traits through exon splicing events. We applied SpliTWAS to two schizophrenia (SCZ) RNA-sequencing datasets, BrainGVEX and CommonMind, revealing 137 and 88 trait-associated exons (in 84 and 67 genes), respectively. Enriched biological functions in the associated gene sets converged on neuronal function and development, immune cell activation, and cellular transport, which are highly relevant to SCZ. SpliTWAS variants impacted RNA-binding protein binding sites, revealing potential disruption of RNA-protein interactions affecting splicing. We extended the probabilistic fine-mapping method FOCUS to the exon level, identifying 36 genes and 48 exons as putatively causal for SCZ. We highlight VPS45 and APOPT1, where splicing of specific exons was associated with disease risk, eluding detection by conventional gene expression analysis. Collectively, this study supports the substantial role of alternative splicing in shaping the genetic basis of SCZ, providing a valuable approach for future investigations in this area.

A transcriptomic atlas of the human brain reveals genetically determined aspects of neuropsychiatric health

Regulation of gene expression is a vital component of neurological homeostasis. Cataloging the consequences of endogenous gene expression on the physical structure and connectivity of the brain offers a means of unifying trait-associated genetic variation with trait-associated neurological features. We perform tissue-specific transcriptome-wide association studies (TWASs) on over 3,400 neuroimaging phenotypes in the UK Biobank (N = 33,224) using our joint-tissue imputation (JTI)-TWAS method. We identify highly significant associations between predicted expression for 7,192 genes and a wide variety of measures of the brain derived from magnetic resonance imaging (MRI). Our approach generates reproducible results in internal and external replication datasets. Genetically determined expression alone is sufficient for high-fidelity reconstruction of brain structure and organization. We demonstrate complementary benefits of cross-tissue and single-tissue analyses toward an integrated neurobiology and provide evidence that gene expression outside the central nervous system provides unique insights into brain health. As an application, we provide evidence suggesting that the genetically regulated expression of schizophrenia risk genes causally affects over 73% of neurological phenotypes that are altered in individuals with schizophrenia (as identified by neuroimaging studies). Imaging features associated with neuropsychiatric traits can provide valuable insights into underlying pathophysiology. By linking neuroimaging-derived phenotypes with expression levels of specific genes, this resource represents a powerful gene prioritization schema that can improve our understanding of brain function, development, and disease. The use of multiple different cortical and subcortical atlases in the resource facilitates direct integration of these data with findings from a diverse range of clinical neuroimaging studies.

Inferring causal direction between two traits using R2 with application to transcriptome-wide association studies

In Mendelian randomization, two single SNP-trait correlation-based methods have been developed to infer the causal direction between an exposure (e.g., a gene) and an outcome (e.g., a trait), called MR Steiger's method and its recent extension called Causal Direction-Ratio (CD-Ratio). Here we propose an approach based on R2, the coefficient of determination, to combine information from multiple (possibly correlated) SNPs to simultaneously infer the presence and direction of a causal relationship between an exposure and an outcome. Our proposed method generalizes Steiger's method from using a single SNP to multiple SNPs as IVs. It is especially useful in transcriptome-wide association studies (TWASs) (and similar applications) with typically small sample sizes for gene expression (or another molecular trait) data, providing a more flexible and powerful approach to inferring causal directions. It can be applied to GWAS summary data with a reference panel. We also discuss the influence of invalid IVs and introduce a new approach called R2S to select and remove invalid IVs (if any) to enhance the robustness. We compared the performance of the proposed method with existing methods in simulations to demonstrate its advantages. We applied the methods to identify causal genes for high/low-density lipoprotein cholesterol (HDL/LDL) using the individual-level GTEx gene expression data and UK Biobank GWAS data. The proposed method was able to confirm some well-known causal genes while identifying some novel ones. Additionally, we illustrated an application of the proposed method to GWAS summary to infer causal relationships between HDL/LDL and stroke/coronary artery disease (CAD).

Leveraging Random Effects in Cistrome-Wide Association Studies for Decoding the Genetic Determinants of Prostate Cancer

Cistrome-wide association studies (CWAS) are pivotal for identifying genetic determinants of diseases by correlating genetically regulated cistrome states with phenotypes. Traditional CWAS typically develops a model based on cistrome and genotype data to associate predicted cistrome states with phenotypes. The random effect cistrome-wide association study (RECWAS), reevaluates the necessity of cistrome state prediction in CWAS. RECWAS utilizes either a linear model or marginal effect for initial feature selection, followed by kernel-based feature aggregation for association testing is introduced. Through simulations and analysis of prostate cancer data, a thorough evaluation of CWAS and RECWAS is conducted. The results suggest that RECWAS offers improved power compared to traditional CWAS, identifying additional genomic regions associated with prostate cancer. CWAS identified 102 significant regions, while RECWAS found 50 additional significant regions compared to CWAS, many of which are validated. Validation encompassed a range of biological evidence, including risk signals from the GWAS catalog, susceptibility genes from the DisGeNET database, and enhancer-domain scores. RECWAS consistently demonstrated improved performance over traditional CWAS in identifying genomic regions associated with prostate cancer. These findings demonstrate the benefits of incorporating kernel methods into CWAS and provide new insights for genetic discovery in complex diseases.

Deciphering the genetics and mechanisms of predisposition to multiple myeloma

Multiple myeloma (MM) is an incurable malignancy of plasma cells. Epidemiological studies indicate a substantial heritable component, but the underlying mechanisms remain unclear. Here, in a genome-wide association study totaling 10,906 cases and 366,221 controls, we identify 35 MM risk loci, 12 of which are novel. Through functional fine-mapping and Mendelian randomization, we uncover two causal mechanisms for inherited MM risk: longer telomeres; and elevated levels of B-cell maturation antigen (BCMA) and interleukin-5 receptor alpha (IL5RA) in plasma. The largest increase in BCMA and IL5RA levels is mediated by the risk variant rs34562254-A at TNFRSF13B. While individuals with loss-of-function variants in TNFRSF13B develop B-cell immunodeficiency, rs34562254-A exerts a gain-of-function effect, increasing MM risk through amplified B-cell responses. Our results represent an analysis of genetic MM predisposition, highlighting causal mechanisms contributing to MM development.

SR-TWAS: leveraging multiple reference panels to improve transcriptome-wide association study power by ensemble machine learning

Multiple reference panels of a given tissue or multiple tissues often exist, and multiple regression methods could be used for training gene expression imputation models for transcriptome-wide association studies (TWAS). To leverage expression imputation models (i.e., base models) trained with multiple reference panels, regression methods, and tissues, we develop a Stacked Regression based TWAS (SR-TWAS) tool which can obtain optimal linear combinations of base models for a given validation transcriptomic dataset. Both simulation and real studies show that SR-TWAS improves power, due to increased training sample sizes and borrowed strength across multiple regression methods and tissues. Leveraging base models across multiple reference panels, tissues, and regression methods, our real studies identify 6 independent significant risk genes for Alzheimer's disease (AD) dementia for supplementary motor area tissue and 9 independent significant risk genes for Parkinson's disease (PD) for substantia nigra tissue. Relevant biological interpretations are found for these significant risk genes.

From genetic associations to genes: methods, applications, and challenges

Genome-wide association studies (GWASs) have identified numerous genetic loci associated with human traits and diseases. However, pinpointing the causal genes remains a challenge, which impedes the translation of GWAS findings into biological insights and medical applications. In this review, we provide an in-depth overview of the methods and technologies used for prioritizing genes from GWAS loci, including gene-based association tests, integrative analysis of GWAS and molecular quantitative trait loci (xQTL) data, linking GWAS variants to target genes through enhancer-gene connection maps, and network-based prioritization. We also outline strategies for generating context-dependent xQTL data and their applications in gene prioritization. We further highlight the potential of gene prioritization in drug repurposing. Lastly, we discuss future challenges and opportunities in this field.

Transcriptome-wide association study and Mendelian randomization in pancreatic cancer identifies susceptibility genes and causal relationships with type 2 diabetes and venous thromboembolism

BACKGROUND

Two important questions regarding the genetics of pancreatic adenocarcinoma (PDAC) are 1. Which germline genetic variants influence the incidence of this cancer; and 2. Whether PDAC causally predisposes to associated non-malignant phenotypes, such as type 2 diabetes (T2D) and venous thromboembolism (VTE).

METHODS

In this study of 8803 patients with PDAC and 67,523 controls, we first performed a large-scale transcriptome-wide association study to investigate the association between genetically determined gene expression in normal pancreas tissue and PDAC risk. Secondly, we used Mendelian Randomization (MR) to analyse the causal relationships among PDAC, T2D (74,124 cases and 824,006 controls) and VTE (30,234 cases and 172,122 controls).

FINDINGS

Sixteen genes showed an association with PDAC risk (FDR <0.10), including six genes not yet reported for PDAC risk (PPIP5K2, TFR2, HNF4G, LRRC10B, PRC1 and FBXL20) and ten previously reported genes (INHBA, SMC2, ABO, PDX1, MTMR6, ACOT2, PGAP3, STARD3, GSDMB, ADAM33). MR provided support for a causal effect of PDAC on T2D using genetic instruments in the HNF4G and PDX1 loci, and unidirectional causality of VTE on PDAC involving the ABO locus (OR 2.12, P < 1e-7). No evidence of a causal effect of PDAC on VTE was found.

INTERPRETATION

These analyses identified candidate susceptibility genes and disease relationships for PDAC that warrant further investigation. HNF4G and PDX1 may induce PDAC-associated diabetes, whereas ABO may induce the causative effect of VTE on PDAC.

FUNDING

National Institutes of Health (USA).

Genomic Analysis Identifies Risk Factors in Restless Legs Syndrome

OBJECTIVE

Restless legs syndrome (RLS) is a neurological condition that causes uncomfortable sensations in the legs and an irresistible urge to move them, typically during periods of rest. The genetic basis and pathophysiology of RLS are incompletely understood. We sought to identify additional novel genetic risk factors associated with RLS susceptibility.

METHODS

We performed a whole-genome sequencing and genome-wide association meta-analysis of RLS cases (n = 9,851) and controls (n = 38,957) in 3 population-based biobanks (All of Us, Canadian Longitudinal Study on Aging, and CARTaGENE).

RESULTS

Genome-wide association analysis identified 9 independent risk loci, of which 8 had been previously reported, and 1 was a novel risk locus (LMX1B, rs35196838, OR 1.14, 95% CI 1.09-1.19, p value = 2.2 × 10-9). Furthermore, a transcriptome-wide association study also identified GLO1 and a previously unreported gene, ELFN1. A genetic correlation analysis revealed significant common variant overlaps between RLS and neuroticism (rg = 0.40, se = 0.08, p value = 5.4 × 10-7), depression (rg = 0.35, se = 0.06, p value = 2.17 × 10-8), and intelligence (rg = -0.20, se = 0.06, p value = 4.0 × 10-4).

INTERPRETATION

Our study expands the understanding of the genetic architecture of RLS, and highlights the contributions of common variants to this prevalent neurological disorder. ANN NEUROL 2024.

European and African-specific plasma protein-QTL and metabolite-QTL analyses identify ancestry-specific T2D effector proteins and metabolites

Initially focused on the European population, multiple genome-wide association studies (GWAS) of complex diseases, such as type-2 diabetes (T2D), have now extended to other populations. However, to date, few ancestry-matched omics datasets have been generated or further integrated with the disease GWAS to nominate the key genes and/or molecular traits underlying the disease risk loci. In this study, we generated and integrated plasma proteomics and metabolomics with array-based genotype datasets of European (EUR) and African (AFR) ancestries to identify ancestry-specific muti-omics quantitative trait loci (QTLs). We further applied these QTLs to ancestry-stratified T2D risk to pinpoint key proteins and metabolites underlying the disease-associated genetic loci. We nominated five proteins and four metabolites in the European group and one protein and one metabolite in the African group to be part of the molecular pathways of T2D risk in an ancestry-stratified manner. Our study demonstrates the integration of genetic and omic studies of different ancestries can be used to identify distinct effector molecular traits underlying the same disease across diverse populations. Specifically, in the AFR proteomic findings on T2D, we prioritized the protein QSOX2; while in the AFR metabolomic findings, we pinpointed the metabolite GlcNAc sulfate conjugate of C21H34O2 steroid. Neither of these findings overlapped with the corresponding EUR results.

Instrumental variable and colocalization analyses identify endotrophin and HTRA1 as potential therapeutic targets for coronary artery disease

Coronary artery disease (CAD) remains a leading cause of disease burden globally, and there is a persistent need for new therapeutic targets. Instrumental variable (IV) and genetic colocalization analyses can help identify novel therapeutic targets for human disease by nominating causal genes in genome-wide association study (GWAS) loci. We conducted cis-IV analyses for 20,125 genes and 1,746 plasma proteins with CAD using molecular trait quantitative trait loci variant (QTLs) data from three different studies. 19 proteins and 119 genes were significantly associated with CAD risk by IV analyses and demonstrated evidence of genetic colocalization. Notably, our analyses validated well-established targets such as PCSK9 and ANGPTL4 while also identifying HTRA1 and endotrophin (a cleavage product of COL6A3) as proteins whose levels are causally associated with CAD risk. Further experimental studies are needed to confirm the causal role of the genes and proteins identified through our multiomic cis-IV analyses on human disease.

Human brain proteome-wide association study provides insights into the genetic components of protein abundance in obesity

BACKGROUNDS

Genome-wide association studies have identified multiple genetic variants associated with obesity. However, most obesity-associated loci were waiting to be translated into new biological insights. Given the critical role of brain in obesity development, we sought to explore whether obesity-associated genetic variants could be mapped to brain protein abundances.

METHODS

We performed proteome-wide association studies (PWAS) and colocalization analyses to identify genes whose cis-regulated brain protein abundances were associated with obesity-related traits, including body fat percentage, trunk fat percentage, body mass index, visceral adipose tissue, waist circumference, and waist-to-hip ratio. We then assessed the druggability of the identified genes and conducted pathway enrichment analysis to explore their functional relevance. Finally, we evaluated the effects of the significant PWAS genes at the brain transcriptional level.

RESULTS

By integrating human brain proteomes from discovery (ROSMAP, N = 376) and validation datasets (BANNER, N = 198) with genome-wide summary statistics of obesity-related phenotypes (N ranged from 325,153 to 806,834), we identified 51 genes whose cis-regulated brain protein abundance was associated with obesity. These 51 genes were enriched in 11 metabolic processes, e.g., small molecule metabolic process and metabolic pathways. Fourteen of the 51 genes had high drug repurposing value. Ten of the 51 genes were also associated with obesity at the transcriptome level, suggesting that genetic variants likely confer risk of obesity by regulating mRNA expression and protein abundance of these genes.

CONCLUSIONS

Our study provides new insights into the genetic component of human brain protein abundance in obesity. The identified proteins represent promising therapeutic targets for future drug development.

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.

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.

Application of GWAS summary data and drug-induced gene expression profiles of neural progenitor cells in psychiatric drug prioritization analysis

Common psychiatric disorders constitute one of the most substantial healthcare burdens worldwide. However, drug development in psychiatry remains hampered partially due to the lack of approaches to estimating drugs that can simultaneously modulate the expression of a nontrivial fraction of disease susceptibility genes. We proposed a new drug prioritization strategy under the framework of our previously proposed phenotype-associated tissues estimation approach (DESE) by investigating the drugs' selective perturbation effect on disease susceptibility genes. Based on the genome-wide association study summary data and drug-induced gene expression profiles of neural progenitor cells, we applied this strategy to prioritize candidate drugs for schizophrenia, depression and bipolar I disorder and identified several known therapeutic drugs among the top-ranked drug candidates. Also, our results revealed that the disease susceptibility genes involved in the selective gene perturbation analysis were enriched with many biologically sensible function terms and interacted with known therapeutic drugs. Our results suggested that selective gene perturbation analysis could be a promising starting point to prioritize biologically sensible drug candidates under the "one drug, multiple targets" paradigm for the drug development of common psychiatric disorders.

Transcriptome-Wide Association Study Reveals New Molecular Interactions Associated with Melanoma Pathogenesis

A transcriptome-wide association study (TWAS) was conducted on genome-wide association study (GWAS) summary statistics of malignant melanoma of skin (UK Biobank dataset) and The Cancer Genome Atlas-Skin Cutaneous Melanoma (TCGA-SKCM) gene expression weights to identify melanoma susceptibility genes. The GWAS included 2465 cases and 449,799 controls, while the gene expression testing was conducted on 103 cases. Afterward, a gene enrichment analysis was applied to identify significant TWAS associations. The melanoma's gene-microRNA (miRNA) regulatory network was constructed from the TWAS genes and their corresponding miRNAs. At last, a disease enrichment analysis was conducted on the corresponding miRNAs. The TWAS detected 27 genes associated with melanoma with p-values less than 0.05 (the top three genes are LOC389458 (RBAK), C16orf73 (MEIOB), and EIF3CL). After the joint/conditional test, one gene (AMIGO1) was dropped, resulting in 26 significant genes. The Gene Ontology (GO) biological process associated the extended gene set (76 genes) with protein K11-linked ubiquitination and regulation of cell cycle phase transition. K11-linked ubiquitin chains regulate cell division. Interestingly, the extended gene set was related to different skin cancer subtypes. Moreover, the enriched pathways were nsp1 from SARS-CoV-2 that inhibit translation initiation in the host cell, cell cycle, translation factors, and DNA repair pathways full network. The gene-miRNA regulatory network identified 10 hotspot genes with the top three: TP53, BRCA1, and MDM2; and four hotspot miRNAs: mir-16, mir-15a, mir-125b, and mir-146a. Melanoma was among the top ten diseases associated with the corresponding (106) miRNAs. Our results shed light on melanoma pathogenesis and biologically significant molecular interactions.

Partitioning and aggregating cross-tissue and tissue-specific genetic effects to identify gene-trait associations

TWAS have shown great promise in extending GWAS loci to a functional understanding of disease mechanisms. In an effort to fully unleash the TWAS and GWAS information, we propose MTWAS, a statistical framework that partitions and aggregates cross-tissue and tissue-specific genetic effects in identifying gene-trait associations. We introduce a non-parametric imputation strategy to augment the inaccessible tissues, accommodating complex interactions and non-linear expression data structures across various tissues. We further classify eQTLs into cross-tissue eQTLs and tissue-specific eQTLs via a stepwise procedure based on the extended Bayesian information criterion, which is consistent under high-dimensional settings. We show that MTWAS significantly improves the prediction accuracy across all 47 tissues of the GTEx dataset, compared with other single-tissue and multi-tissue methods, such as PrediXcan, TIGAR, and UTMOST. Applying MTWAS to the DICE and OneK1K datasets with bulk and single-cell RNA sequencing data on immune cell types showcases consistent improvements in prediction accuracy. MTWAS also identifies more predictable genes, and the improvement can be replicated with independent studies. We apply MTWAS to 84 UK Biobank GWAS studies, which provides insights into disease etiology.

Dynamic stress- and inflammatory-based regulation of psychiatric risk loci in human neurons

The prenatal environment can alter neurodevelopmental and clinical trajectories, markedly increasing risk for psychiatric disorders in childhood and adolescence. To understand if and how fetal exposures to stress and inflammation exacerbate manifestation of genetic risk for complex brain disorders, we report a large-scale context-dependent massively parallel reporter assay (MPRA) in human neurons designed to catalogue genotype x environment (GxE) interactions. Across 240 genome-wide association study (GWAS) loci linked to ten brain traits/disorders, the impact of hydrocortisone, interleukin 6, and interferon alpha on transcriptional activity is empirically evaluated in human induced pluripotent stem cell (hiPSC)-derived glutamatergic neurons. Of ~3,500 candidate regulatory risk elements (CREs), 11% of variants are active at baseline, whereas cue-specific CRE regulatory activity range from a high of 23% (hydrocortisone) to a low of 6% (IL-6). Cue-specific regulatory activity is driven, at least in part, by differences in transcription factor binding activity, the gene targets of which show unique enrichments for brain disorders as well as co-morbid metabolic and immune syndromes. The dynamic nature of genetic regulation informs the influence of environmental factors, reveals a mechanism underlying pleiotropy and variable penetrance, and identifies specific risk variants that confer greater disorder susceptibility after exposure to stress or inflammation. Understanding neurodevelopmental GxE interactions will inform mental health trajectories and uncover novel targets for therapeutic intervention.

Metabolic gene function discovery platform GeneMAP identifies SLC25A48 as necessary for mitochondrial choline import

Organisms maintain metabolic homeostasis through the combined functions of small-molecule transporters and enzymes. While many metabolic components have been well established, a substantial number remains without identified physiological substrates. To bridge this gap, we have leveraged large-scale plasma metabolome genome-wide association studies (GWAS) to develop a multiomic Gene-Metabolite Association Prediction (GeneMAP) discovery platform. GeneMAP can generate accurate predictions and even pinpoint genes that are distant from the variants implicated by GWAS. In particular, our analysis identified solute carrier family 25 member 48 (SLC25A48) as a genetic determinant of plasma choline levels. Mechanistically, SLC25A48 loss strongly impairs mitochondrial choline import and synthesis of its downstream metabolite betaine. Integrative rare variant and polygenic score analyses in UK Biobank provide strong evidence that the SLC25A48 causal effects on human disease may in part be mediated by the effects of choline. Altogether, our study provides a discovery platform for metabolic gene function and proposes SLC25A48 as a mitochondrial choline transporter.

Machine Learning of Three-Dimensional Protein Structures to Predict the Functional Impacts of Genome Variation

Research in the human genome sciences generates a substantial amount of genetic data for hundreds of thousands of individuals, which concomitantly increases the number of variants of unknown significance (VUS). Bioinformatic analyses can successfully reveal rare variants and variants with clear associations with disease-related phenotypes. These studies have had a significant impact on how clinical genetic screens are interpreted and how patients are stratified for treatment. There are few, if any, computational methods for variants comparable to biological activity predictions. To address this gap, we developed a machine learning method that uses protein three-dimensional structures from AlphaFold to predict how a variant will influence changes to a gene's downstream biological pathways. We trained state-of-the-art machine learning classifiers to predict which protein regions will most likely impact transcriptional activities of two proto-oncogenes, nuclear factor erythroid 2 (NFE2L2)-related factor 2 (NRF2) and c-Myc. We have identified classifiers that attain accuracies higher than 80%, which have allowed us to identify a set of key protein regions that lead to significant perturbations in c-Myc or NRF2 transcriptional pathway activities.

Enhancing Disease Risk Gene Discovery by Integrating Transcription Factor-Linked Trans-located Variants into Transcriptome-Wide Association Analyses

Transcriptome-wide association studies (TWAS) have been successful in identifying disease susceptibility genes by integrating cis-variants predicted gene expression with genome-wide association studies (GWAS) data. However, trans-located variants for predicting gene expression remain largely unexplored. Here, we introduce transTF-TWAS, which incorporates transcription factor (TF)-linked trans-located variants to enhance model building. Using data from the Genotype-Tissue Expression project, we predict gene expression and alternative splicing and applied these models to large GWAS datasets for breast, prostate, and lung cancers. We demonstrate that transTF-TWAS outperforms other existing TWAS approaches in both constructing gene prediction models and identifying disease-associated genes, as evidenced by simulations and real data analysis. Our transTF-TWAS approach significantly contributes to the discovery of disease risk genes. Findings from this study have shed new light on several genetically driven key regulators and their associated regulatory networks underlying disease susceptibility.

Abdominal aortic aneurysm and cardiometabolic traits share strong genetic susceptibility to lipid metabolism and inflammation

Abdominal aortic aneurysm has a high heritability and often co-occurs with other cardiometabolic disorders, suggesting shared genetic susceptibility. We investigate this commonality leveraging recent GWAS studies of abdominal aortic aneurysm and 32 cardiometabolic traits. We find significant genetic correlations between abdominal aortic aneurysm and 21 of the cardiometabolic traits investigated, including causal relationships with coronary artery disease, hypertension, lipid traits, and blood pressure. For each trait pair, we identify shared causal variants, genes, and pathways, revealing that cholesterol metabolism and inflammation are shared most prominently. Additionally, we show the tissue and cell type specificity in the shared signals, with strong enrichment across traits in the liver, arteries, adipose tissues, macrophages, adipocytes, and fibroblasts. Finally, we leverage drug-gene databases to identify several lipid-lowering drugs and antioxidants with high potential to treat abdominal aortic aneurysm with comorbidities. Our study provides insight into the shared genetic mechanism between abdominal aortic aneurysm and cardiometabolic traits, and identifies potential targets for pharmacological intervention.

GRPa-PRS: A risk stratification method to identify genetically-regulated pathways in polygenic diseases

Background

Polygenic risk scores (PRS) are tools used to evaluate an individual's susceptibility to polygenic diseases based on their genetic profile. A considerable proportion of people carry a high genetic risk but evade the disease. On the other hand, some individuals with a low risk of eventually developing the disease. We hypothesized that unknown counterfactors might be involved in reversing the PRS prediction, which might provide new insights into the pathogenesis, prevention, and early intervention of diseases.

Methods

We built a novel computational framework to identify genetically-regulated pathways (GRPas) using PRS-based stratification for each cohort. We curated two AD cohorts with genotyping data; the discovery (disc) and the replication (rep) datasets include 2722 and 2854 individuals, respectively. First, we calculated the optimized PRS model based on the three recent AD GWAS summary statistics for each cohort. Then, we stratified the individuals by their PRS and clinical diagnosis into six biologically meaningful PRS strata, such as AD cases with low/high risk and cognitively normal (CN) with low/high risk. Lastly, we imputed individual genetically-regulated expression (GReX) and identified differential GReX and GRPas between risk strata using gene-set enrichment and variational analyses in two models, with and without APOE effects. An orthogonality test was further conducted to verify those GRPas are independent of PRS risk. To verify the generalizability of other polygenic diseases, we further applied a default model of GRPa-PRS for schizophrenia (SCZ).

Results

For each stratum, we conducted the same procedures in both the disc and rep datasets for comparison. In AD, we identified several well-known AD-related pathways, including amyloid-beta clearance, tau protein binding, and astrocyte response to oxidative stress. Additionally, we discovered resilience-related GRPs that are orthogonal to AD PRS, such as the calcium signaling pathway and divalent inorganic cation homeostasis. In SCZ, pathways related to mitochondrial function and muscle development were highlighted. Finally, our GRPa-PRS method identified more consistent differential pathways compared to another variant-based pathway PRS method.

Conclusions

We developed a framework, GRPa-PRS, to systematically explore the differential GReX and GRPas among individuals stratified by their estimated PRS. The GReX-level comparison among those strata unveiled new insights into the pathways associated with disease risk and resilience. Our framework is extendable to other polygenic complex diseases.

Analysis of Evolutionary Conservation, Expression Level, and Genetic Association at a Genome-wide Scale Reveals Heterogeneity Across Polygenic Phenotypes

Understanding the expression level and evolutionary rate of associated genes with human polygenic diseases provides crucial insights into their disease-contributing roles. In this work, we leveraged genome-wide association studies (GWASs) to investigate the relationship between the genetic association and both the evolutionary rate (dN/dS) and expression level of human genes associated with the two polygenic diseases of schizophrenia and coronary artery disease. Our findings highlight a distinct variation in these relationships between the two diseases. Genes associated with both diseases exhibit a significantly greater variance in evolutionary rate compared to those implicated in monogenic diseases. Expanding our analyses to 4,756 complex traits in the GWAS atlas database, we unraveled distinct trait categories with a unique interplay among the evolutionary rate, expression level, and genetic association of human genes. In most polygenic traits, highly expressed genes were more associated with the polygenic phenotypes compared to lowly expressed genes. About 69% of polygenic traits displayed a negative correlation between genetic association and evolutionary rate, while approximately 30% of these traits showed a positive correlation between genetic association and evolutionary rate. Our results demonstrate the presence of a spectrum among complex traits, shaped by natural selection. Notably, at opposite ends of this spectrum, we find metabolic traits being more likely influenced by purifying selection, and immunological traits that are more likely shaped by positive selection. We further established the polygenic evolution portal (evopolygen.de) as a resource for investigating relationships and generating hypotheses in the field of human polygenic trait evolution.

MICOS Complex Loss Governs Age-Associated Murine Mitochondrial Architecture and Metabolism in the Liver, While Sam50 Dictates Diet Changes

The liver, the largest internal organ and a metabolic hub, undergoes significant declines due to aging, affecting mitochondrial function and increasing the risk of systemic liver diseases. How the mitochondrial three-dimensional (3D) structure changes in the liver across aging, and the biological mechanisms regulating such changes confers remain unclear. In this study, we employed Serial Block Face-Scanning Electron Microscopy (SBF-SEM) to achieve high-resolution 3D reconstructions of murine liver mitochondria to observe diverse phenotypes and structural alterations that occur with age, marked by a reduction in size and complexity. We also show concomitant metabolomic and lipidomic changes in aged samples. Aged human samples reflected altered disease risk. To find potential regulators of this change, we examined the Mitochondrial Contact Site and Cristae Organizing System (MICOS) complex, which plays a crucial role in maintaining mitochondrial architecture. We observe that the MICOS complex is lost during aging, but not Sam50. Sam50 is a component of the sorting and assembly machinery (SAM) complex that acts in tandem with the MICOS complex to modulate cristae morphology. In murine models subjected to a high-fat diet, there is a marked depletion of the mitochondrial protein SAM50. This reduction in Sam50 expression may heighten the susceptibility to liver disease, as our human biobank studies corroborate that Sam50 plays a genetically regulated role in the predisposition to multiple liver diseases. We further show that changes in mitochondrial calcium dysregulation and oxidative stress accompany the disruption of the MICOS complex. Together, we establish that a decrease in mitochondrial complexity and dysregulated metabolism occur with murine liver aging. While these changes are partially be regulated by age-related loss of the MICOS complex, the confluence of a murine high-fat diet can also cause loss of Sam50, which contributes to liver diseases. In summary, our study reveals potential regulators that affect age-related changes in mitochondrial structure and metabolism, which can be targeted in future therapeutic techniques.

Identification of the Molecular Components of Enhancer-Mediated Gene Expression Variation in Multiple Tissues Regulating Blood Pressure

BACKGROUND

Inter-individual variation in blood pressure (BP) arises in part from sequence variants within enhancers modulating the expression of causal genes. We propose that these genes, active in tissues relevant to BP physiology, can be identified from tissue-level epigenomic data and genotypes of BP-phenotyped individuals.

METHODS

We used chromatin accessibility data from the heart, adrenal, kidney, and artery to identify cis-regulatory elements (CREs) in these tissues and estimate the impact of common human single-nucleotide variants within these CREs on gene expression, using machine learning methods. To identify causal genes, we performed a gene-wise association test. We conducted analyses in 2 separate large-scale cohorts: 77 822 individuals from the Genetic Epidemiology Research on Adult Health and Aging and 315 270 individuals from the UK Biobank.

RESULTS

We identified 309, 259, 331, and 367 genes (false discovery rate <0.05) for diastolic BP and 191, 184, 204, and 204 genes for systolic BP in the artery, kidney, heart, and adrenal, respectively, in Genetic Epidemiology Research on Adult Health and Aging; 50% to 70% of these genes were replicated in the UK Biobank, significantly higher than the 12% to 15% expected by chance (P<0.0001). These results enabled tissue expression prediction of these 988 to 2875 putative BP genes in individuals of both cohorts to construct an expression polygenic score. This score explained ≈27% of the reported single-nucleotide variant heritability, substantially higher than expected from prior studies.

CONCLUSIONS

Our work demonstrates the power of tissue-restricted comprehensive CRE analysis, followed by CRE-based expression prediction, for understanding BP regulation in relevant tissues and provides dual-modality supporting evidence, CRE and expression, for the causality genes.

Distinct genetic liability profiles define clinically relevant patient strata across common diseases

Stratified medicine holds great promise to tailor treatment to the needs of individual patients. While genetics holds great potential to aid patient stratification, it remains a major challenge to operationalize complex genetic risk factor profiles to deconstruct clinical heterogeneity. Contemporary approaches to this problem rely on polygenic risk scores (PRS), which provide only limited clinical utility and lack a clear biological foundation. To overcome these limitations, we develop the CASTom-iGEx approach to stratify individuals based on the aggregated impact of their genetic risk factor profiles on tissue specific gene expression levels. The paradigmatic application of this approach to coronary artery disease or schizophrenia patient cohorts identified diverse strata or biotypes. These biotypes are characterized by distinct endophenotype profiles as well as clinical parameters and are fundamentally distinct from PRS based groupings. In stark contrast to the latter, the CASTom-iGEx strategy discovers biologically meaningful and clinically actionable patient subgroups, where complex genetic liabilities are not randomly distributed across individuals but rather converge onto distinct disease relevant biological processes. These results support the notion of different patient biotypes characterized by partially distinct pathomechanisms. Thus, the universally applicable approach presented here has the potential to constitute an important component of future personalized medicine paradigms.

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.

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.

TWAS facilitates gene-scale trait genetic dissection through gene expression, structural variations, and alternative splicing in soybean

A genome-wide association study (GWAS) identifies trait-associated loci, but identifying the causal genes can be a bottleneck, due in part to slow decay of linkage disequilibrium (LD). A transcriptome-wide association study (TWAS) addresses this issue by identifying gene expression-phenotype associations or integrating gene expression quantitative trait loci with GWAS results. Here, we used self-pollinated soybean (Glycine max [L.] Merr.) as a model to evaluate the application of TWAS to the genetic dissection of traits in plant species with slow LD decay. We generated RNA sequencing data for a soybean diversity panel and identified the genetic expression regulation of 29 286 soybean genes. Different TWAS solutions were less affected by LD and were robust to the source of expression, identifing known genes related to traits from different tissues and developmental stages. The novel pod-color gene L2 was identified via TWAS and functionally validated by genome editing. By introducing a new exon proportion feature, we significantly improved the detection of expression variations that resulted from structural variations and alternative splicing. As a result, the genes identified through our TWAS approach exhibited a diverse range of causal variations, including SNPs, insertions or deletions, gene fusion, copy number variations, and alternative splicing. Using this approach, we identified genes associated with flowering time, including both previously known genes and novel genes that had not previously been linked to this trait, providing insights complementary to those from GWAS. In summary, this study supports the application of TWAS for candidate gene identification in species with low rates of LD decay.

Susceptibility gene identification and risk evaluation model construction by transcriptome-wide association analysis for salt sensitivity of blood pressure

BACKGROUND

Salt sensitivity of blood pressure (SSBP) is an intermediate phenotype of hypertension and is a predictor of long-term cardiovascular events and death. However, the genetic structures of SSBP are uncertain, and it is difficult to precisely diagnose SSBP in population. So, we aimed to identify genes related to susceptibility to the SSBP, construct a risk evaluation model, and explore the potential functions of these genes.

METHODS AND RESULTS

A genome-wide association study of the systemic epidemiology of salt sensitivity (EpiSS) cohort was performed to obtain summary statistics for SSBP. Then, we conducted a transcriptome-wide association study (TWAS) of 12 tissues using FUSION software to predict the genes associated with SSBP and verified the genes with an mRNA microarray. The potential roles of the genes were explored. Risk evaluation models of SSBP were constructed based on the serial P value thresholds of polygenetic risk scores (PRSs), polygenic transcriptome risk scores (PTRSs) and their combinations of the identified genes and genetic variants from the TWAS. The TWAS revealed that 2605 genes were significantly associated with SSBP. Among these genes, 69 were differentially expressed according to the microarray analysis. The functional analysis showed that the genes identified in the TWAS were enriched in metabolic process pathways. The PRSs were correlated with PTRSs in the heart atrial appendage, adrenal gland, EBV-transformed lymphocytes, pituitary, artery coronary, artery tibial and whole blood. Multiple logistic regression models revealed that a PRS of P < 0.05 had the best predictive ability compared with other PRSs and PTRSs. The combinations of PRSs and PTRSs did not significantly increase the prediction accuracy of SSBP in the training and validation datasets.

CONCLUSIONS

Several known and novel susceptibility genes for SSBP were identified via multitissue TWAS analysis. The risk evaluation model constructed with the PRS of susceptibility genes showed better diagnostic performance than the transcript levels, which could be applied to screen for SSBP high-risk individuals.

Fine-mapping causal tissues and genes at disease-associated loci

Heritable diseases often manifest in a highly tissue-specific manner, with different disease loci mediated by genes in distinct tissues or cell types. We propose Tissue-Gene Fine-Mapping (TGFM), a fine-mapping method that infers the posterior probability (PIP) for each gene-tissue pair to mediate a disease locus by analyzing GWAS summary statistics (and in-sample LD) and leveraging eQTL data from diverse tissues to build cis-predicted expression models; TGFM also assigns PIPs to causal variants that are not mediated by gene expression in assayed genes and tissues. TGFM accounts for both co-regulation across genes and tissues and LD between SNPs (generalizing existing fine-mapping methods), and incorporates genome-wide estimates of each tissue's contribution to disease as tissue-level priors. TGFM was well-calibrated and moderately well-powered in simulations; unlike previous methods, TGFM was able to attain correct calibration by modeling uncertainty in cis-predicted expression models. We applied TGFM to 45 UK Biobank diseases/traits (average N = 316K) using eQTL data from 38 GTEx tissues. TGFM identified an average of 147 PIP > 0.5 causal genetic elements per disease/trait, of which 11% were gene-tissue pairs. Implicated gene-tissue pairs were concentrated in known disease-critical tissues, and causal genes were strongly enriched in disease-relevant gene sets. Causal gene-tissue pairs identified by TGFM recapitulated known biology (e.g., TPO-thyroid for Hypothyroidism), but also included biologically plausible novel findings (e.g., SLC20A2-artery aorta for Diastolic blood pressure). Further application of TGFM to single-cell eQTL data from 9 cell types in peripheral blood mononuclear cells (PBMC), analyzed jointly with GTEx tissues, identified 30 additional causal gene-PBMC cell type pairs at PIP > 0.5-primarily for autoimmune disease and blood cell traits, including the biologically plausible example of CD52 in classical monocyte cells for Monocyte count. In conclusion, TGFM is a robust and powerful method for fine-mapping causal tissues and genes at disease-associated loci.

Hierarchical joint analysis of marginal summary statistics-Part II: High-dimensional instrumental analysis of omics data

Instrumental variable (IV) analysis has been widely applied in epidemiology to infer causal relationships using observational data. Genetic variants can also be viewed as valid IVs in Mendelian randomization and transcriptome-wide association studies. However, most multivariate IV approaches cannot scale to high-throughput experimental data. Here, we leverage the flexibility of our previous work, a hierarchical model that jointly analyzes marginal summary statistics (hJAM), to a scalable framework (SHA-JAM) that can be applied to a large number of intermediates and a large number of correlated genetic variants-situations often encountered in modern experiments leveraging omic technologies. SHA-JAM aims to estimate the conditional effect for high-dimensional risk factors on an outcome by incorporating estimates from association analyses of single-nucleotide polymorphism (SNP)-intermediate or SNP-gene expression as prior information in a hierarchical model. Results from extensive simulation studies demonstrate that SHA-JAM yields a higher area under the receiver operating characteristics curve (AUC), a lower mean-squared error of the estimates, and a much faster computation speed, compared to an existing approach for similar analyses. In two applied examples for prostate cancer, we investigated metabolite and transcriptome associations, respectively, using summary statistics from a GWAS for prostate cancer with more than 140,000 men and high dimensional publicly available summary data for metabolites and transcriptomes.

Hippocampal transcriptome-wide association study and pathway analysis of mitochondrial solute carriers in Alzheimer's disease

The etiopathogenesis of late-onset Alzheimer's disease (AD) is increasingly recognized as the result of the combination of the aging process, toxic proteins, brain dysmetabolism, and genetic risks. Although the role of mitochondrial dysfunction in the pathogenesis of AD has been well-appreciated, the interaction between mitochondrial function and genetic variability in promoting dementia is still poorly understood. In this study, by tissue-specific transcriptome-wide association study (TWAS) and further meta-analysis, we examined the genetic association between mitochondrial solute carrier family (SLC25) genes and AD in three independent cohorts and identified three AD-susceptibility genes, including SLC25A10, SLC25A17, and SLC25A22. Integrative analysis using neuroimaging data and hippocampal TWAS-predicted gene expression of the three susceptibility genes showed an inverse correlation of SLC25A22 with hippocampal atrophy rate in AD patients, which outweighed the impacts of sex, age, and apolipoprotein E4 (ApoE4). Furthermore, SLC25A22 downregulation demonstrated an association with AD onset, as compared with the other two transcriptome-wide significant genes. Pathway and network analysis related hippocampal SLC25A22 downregulation to defects in neuronal function and development, echoing the enrichment of SLC25A22 expression in human glutamatergic neurons. The most parsimonious interpretation of the results is that we have identified AD-susceptibility genes in the SLC25 family through the prediction of hippocampal gene expression. Moreover, our findings mechanistically yield insight into the mitochondrial cascade hypothesis of AD and pave the way for the future development of diagnostic tools for the early prevention of AD from a perspective of precision medicine by targeting the mitochondria-related genes.

An integrative framework to prioritize genes in more than 500 loci associated with body mass index

Obesity is a major risk factor for a myriad of diseases, affecting >600 million people worldwide. Genome-wide association studies (GWASs) have identified hundreds of genetic variants that influence body mass index (BMI), a commonly used metric to assess obesity risk. Most variants are non-coding and likely act through regulating genes nearby. Here, we apply multiple computational methods to prioritize the likely causal gene(s) within each of the 536 previously reported GWAS-identified BMI-associated loci. We performed summary-data-based Mendelian randomization (SMR), FINEMAP, DEPICT, MAGMA, transcriptome-wide association studies (TWASs), mutation significance cutoff (MSC), polygenic priority score (PoPS), and the nearest gene strategy. Results of each method were weighted based on their success in identifying genes known to be implicated in obesity, ranking all prioritized genes according to a confidence score (minimum: 0; max: 28). We identified 292 high-scoring genes (≥11) in 264 loci, including genes known to play a role in body weight regulation (e.g., DGKI, ANKRD26, MC4R, LEPR, BDNF, GIPR, AKT3, KAT8, MTOR) and genes related to comorbidities (e.g., FGFR1, ISL1, TFAP2B, PARK2, TCF7L2, GSK3B). For most of the high-scoring genes, however, we found limited or no evidence for a role in obesity, including the top-scoring gene BPTF. Many of the top-scoring genes seem to act through a neuronal regulation of body weight, whereas others affect peripheral pathways, including circadian rhythm, insulin secretion, and glucose and carbohydrate homeostasis. The characterization of these likely causal genes can increase our understanding of the underlying biology and offer avenues to develop therapeutics for weight loss.

Cis- and trans-eQTL TWASs of breast and ovarian cancer identify more than 100 susceptibility genes in the BCAC and OCAC consortia

Transcriptome-wide association studies (TWASs) have investigated the role of genetically regulated transcriptional activity in the etiologies of breast and ovarian cancer. However, methods performed to date have focused on the regulatory effects of risk-associated SNPs thought to act in cis on a nearby target gene. With growing evidence for distal (trans) regulatory effects of variants on gene expression, we performed TWASs of breast and ovarian cancer using a Bayesian genome-wide TWAS method (BGW-TWAS) that considers effects of both cis- and trans-expression quantitative trait loci (eQTLs). We applied BGW-TWAS to whole-genome and RNA sequencing data in breast and ovarian tissues from the Genotype-Tissue Expression project to train expression imputation models. We applied these models to large-scale GWAS summary statistic data from the Breast Cancer and Ovarian Cancer Association Consortia to identify genes associated with risk of overall breast cancer, non-mucinous epithelial ovarian cancer, and 10 cancer subtypes. We identified 101 genes significantly associated with risk with breast cancer phenotypes and 8 with ovarian phenotypes. These loci include established risk genes and several novel candidate risk loci, such as ACAP3, whose associations are predominantly driven by trans-eQTLs. We replicated several associations using summary statistics from an independent GWAS of these cancer phenotypes. We further used genotype and expression data in normal and tumor breast tissue from the Cancer Genome Atlas to examine the performance of our trained expression imputation models. This work represents an in-depth look into the role of trans eQTLs in the complex molecular mechanisms underlying these diseases.

A cross-tissue transcriptome-wide association study reveals novel susceptibility genes for migraine

BACKGROUND

Migraine is a common neurological disorder with a strong genetic component. Despite the identification of over 100 loci associated with migraine susceptibility through genome-wide association studies (GWAS), the underlying causative genes and biological mechanisms remain predominantly elusive.

METHODS

The FinnGen R10 dataset, consisting of 333,711 subjects (20,908 cases and 312,803 controls), was utilized in conjunction with the Genotype-Tissue Expression Project (GTEx) v8 EQTls files to conduct cross-tissue transcriptome association studies (TWAS). Functional Summary-based Imputation (FUSION) was employed to validate these findings in single tissues. Additionally, candidate susceptibility genes were screened using Gene Analysis combined with Multi-marker Analysis of Genomic Annotation (MAGMA). Subsequent Mendelian randomization (MR) and colocalization analyses were conducted. Furthermore, GeneMANIA analysis was employed to enhance our understanding of the functional implications of these susceptibility genes.

RESULTS

We identified a total of 19 susceptibility genes associated with migraine in the cross-tissue TWAS analysis. Two novel susceptibility genes, REV1 and SREBF2, were validated through both single tissue TWAS and MAGMA analysis. Mendelian randomization and colocalization analyses further confirmed these findings. REV1 may reduce the migraine risk by regulating DNA damage repair, while SREBF2 may increase the risk of migraine by regulating cholesterol metabolism.

CONCLUSION

Our study identified two novel genes whose predicted expression was associated with the risk of migraine, providing new insights into the genetic framework of migraine.

Exploring the Interplay between the Hologenome and Complex Traits in Bovine and Porcine Animals Using Genome-Wide Association Analysis

Genome-wide association studies (GWAS) significantly enhance our ability to identify trait-associated genomic variants by considering the host genome. Moreover, the hologenome refers to the host organism's collective genetic material and its associated microbiome. In this study, we utilized the hologenome framework, called Hologenome-wide association studies (HWAS), to dissect the architecture of complex traits, including milk yield, methane emissions, rumen physiology in cattle, and gut microbial composition in pigs. We employed four statistical models: (1) GWAS, (2) Microbial GWAS (M-GWAS), (3) HWAS-CG (hologenome interaction estimated using COvariance between Random Effects Genome-based restricted maximum likelihood (CORE-GREML)), and (4) HWAS-H (hologenome interaction estimated using the Hadamard product method). We applied Bonferroni correction to interpret the significant associations in the complex traits. The GWAS and M-GWAS detected one and sixteen significant SNPs for milk yield traits, respectively, whereas the HWAS-CG and HWAS-H each identified eight SNPs. Moreover, HWAS-CG revealed four, and the remaining models identified three SNPs each for methane emissions traits. The GWAS and HWAS-CG detected one and three SNPs for rumen physiology traits, respectively. For the pigs' gut microbial composition traits, the GWAS, M-GWAS, HWAS-CG, and HWAS-H identified 14, 16, 13, and 12 SNPs, respectively. We further explored these associations through SNP annotation and by analyzing biological processes and functional pathways. Additionally, we integrated our GWA results with expression quantitative trait locus (eQTL) data using transcriptome-wide association studies (TWAS) and summary-based Mendelian randomization (SMR) methods for a more comprehensive understanding of SNP-trait associations. Our study revealed hologenomic variability in agriculturally important traits, enhancing our understanding of host-microbiome interactions.

Development of a Longitudinal Prostate Cancer Transcriptomic and Clinical Data Linkage

Importance

Although tissue-based gene expression testing has become widely used for prostate cancer risk stratification, its prognostic performance in the setting of clinical care is not well understood.

Objective

To develop a linkage between a prostate genomic classifier (GC) and clinical data across payers and sites of care in the US.

Design, Setting, and Participants

In this cohort study, clinical and transcriptomic data from clinical use of a prostate GC between 2016 and 2022 were linked with data aggregated from insurance claims, pharmacy records, and electronic health record (EHR) data. Participants were anonymously linked between datasets by deterministic methods through a deidentification engine using encrypted tokens. Algorithms were developed and refined for identifying prostate cancer diagnoses, treatment timing, and clinical outcomes using diagnosis codes, Common Procedural Terminology codes, pharmacy codes, Systematized Medical Nomenclature for Medicine clinical terms, and unstructured text in the EHR. Data analysis was performed from January 2023 to January 2024.

Exposure

Diagnosis of prostate cancer.

Main Outcomes and Measures

The primary outcomes were biochemical recurrence and development of prostate cancer metastases after diagnosis or radical prostatectomy (RP). The sensitivity of the linkage and identification algorithms for clinical and administrative data were calculated relative to clinical and pathological information obtained during the GC testing process as the reference standard.

Results

A total of 92 976 of 95 578 (97.2%) participants who underwent prostate GC testing were successfully linked to administrative and clinical data, including 53 871 who underwent biopsy testing and 39 105 who underwent RP testing. The median (IQR) age at GC testing was 66.4 (61.0-71.0) years. The sensitivity of the EHR linkage data for prostate cancer diagnoses was 85.0% (95% CI, 84.7%-85.2%), including 80.8% (95% CI, 80.4%-81.1%) for biopsy-tested participants and 90.8% (95% CI, 90.5%-91.0%) for RP-tested participants. Year of treatment was concordant in 97.9% (95% CI, 97.7%-98.1%) of those undergoing GC testing at RP, and 86.0% (95% CI, 85.6%-86.4%) among participants undergoing biopsy testing. The sensitivity of the linkage was 48.6% (95% CI, 48.1%-49.1%) for identifying RP and 50.1% (95% CI, 49.7%-50.5%) for identifying prostate biopsy.

Conclusions and Relevance

This study established a national-scale linkage of transcriptomic and longitudinal clinical data yielding high accuracy for identifying key clinical junctures, including diagnosis, treatment, and early cancer outcome. This resource can be leveraged to enhance understandings of disease biology, patterns of care, and treatment effectiveness.

Bayesian genome-wide TWAS with reference transcriptomic data of brain and blood tissues identified 141 risk genes for Alzheimer's disease dementia

BACKGROUND

Transcriptome-wide association study (TWAS) is an influential tool for identifying genes associated with complex diseases whose genetic effects are likely mediated through transcriptome. TWAS utilizes reference genetic and transcriptomic data to estimate effect sizes of genetic variants on gene expression (i.e., effect sizes of a broad sense of expression quantitative trait loci, eQTL). These estimated effect sizes are employed as variant weights in gene-based association tests, facilitating the mapping of risk genes with genome-wide association study (GWAS) data. However, most existing TWAS of Alzheimer's disease (AD) dementia are limited to studying only cis-eQTL proximal to the test gene. To overcome this limitation, we applied the Bayesian Genome-wide TWAS (BGW-TWAS) method to leveraging both cis- and trans- eQTL of brain and blood tissues, in order to enhance mapping risk genes for AD dementia.

METHODS

We first applied BGW-TWAS to the Genotype-Tissue Expression (GTEx) V8 dataset to estimate cis- and trans- eQTL effect sizes of the prefrontal cortex, cortex, and whole blood tissues. Estimated eQTL effect sizes were integrated with the summary data of the most recent GWAS of AD dementia to obtain BGW-TWAS (i.e., gene-based association test) p-values of AD dementia per gene per tissue type. Then we used the aggregated Cauchy association test to combine TWAS p-values across three tissues to obtain omnibus TWAS p-values per gene.

RESULTS

We identified 85 significant genes in prefrontal cortex, 82 in cortex, and 76 in whole blood that were significantly associated with AD dementia. By combining BGW-TWAS p-values across these three tissues, we obtained 141 significant risk genes including 34 genes primarily due to trans-eQTL and 35 mapped risk genes in GWAS Catalog. With these 141 significant risk genes, we detected functional clusters comprised of both known mapped GWAS risk genes of AD in GWAS Catalog and our identified TWAS risk genes by protein-protein interaction network analysis, as well as several enriched phenotypes related to AD.

CONCLUSION

We applied BGW-TWAS and aggregated Cauchy test methods to integrate both cis- and trans- eQTL data of brain and blood tissues with GWAS summary data, identifying 141 TWAS risk genes of AD dementia. These identified risk genes provide novel insights into the underlying biological mechanisms of AD dementia and potential gene targets for therapeutics development.

Mendelian randomization: causal inference leveraging genetic data

Mendelian randomization (MR) leverages genetic information to examine the causal relationship between phenotypes allowing for the presence of unmeasured confounders. MR has been widely applied to unresolved questions in epidemiology, making use of summary statistics from genome-wide association studies on an increasing number of human traits. However, an understanding of essential concepts is necessary for the appropriate application and interpretation of MR. This review aims to provide a non-technical overview of MR and demonstrate its relevance to psychiatric research. We begin with the origins of MR and the reasons for its recent expansion, followed by an overview of its statistical methodology. We then describe the limitations of MR, and how these are being addressed by recent methodological advances. We showcase the practical use of MR in psychiatry through three illustrative examples - the connection between cannabis use and psychosis, the link between intelligence and schizophrenia, and the search for modifiable risk factors for depression. The review concludes with a discussion of the prospects of MR, focusing on the integration of multi-omics data and its extension to delineating complex causal networks.