A transcriptome-wide association study of 229,000 women identifies new candidate susceptibility genes for breast cancer

Prognostic significance of neutrophil extracellular trap-related genes in childhood acute lymphoblastic leukemia: insights from multi-omics and in vitro experiment

BACKGROUND

This study aimed to develop a prognostic model based on extracellular trap-related genes (NETRGs) for patients with cALL.

METHODS

Data from the TARGET-ALL-P2 and TARGET-ALL-P3 cohorts in the Genomic Data Commons database, the transcriptome dataset GSE26713, the single-cell transcriptome dataset GSE130116 from the Gene Expression Omnibus database and 306 NETRGs identified were analysed. Differentially expressed genes (DEGs) were identified from GSE26713 and differentially expressed NETRGs (DE-NETRGs) were obtained by overlapping DEGs with NETRGs. Functional analyses were conducted. Key feature genes were identified through univariate and least absolute shrinkage and selection operator (LASSO) regression. Prognostic genes were determined via multivariate Cox regression analysis, followed by the construction and validation of a risk model and nomogram. Additional analyses included immune profiling, drug sensitivity, functional differences, cell-type-specific expression, enrichment analysis and RT-qPCR.

RESULTS

A total of 1,270 DEGs were identified in GSE26713, of which 74 overlapped with NETRGs. Seven prognostic genes were identified using univariate, LASSO and multivariate Cox regression analyses. Survival analysis revealed lower survival rates in the high-risk group. Independent prognostic analysis identified risk scores and primary diagnosis as independent predictors of prognosis. Immune cell profiling showed significant differences in cell populations such as aDCs, eosinophils and Th2 cells between risk groups. Six cell subtypes were annotated, with prognostic genes predominantly expressed in myeloid cells. RT-qPCR revealed that PTAFR, FCGR2A, RETN and CAT were significantly downregulated, while TLR2 and S100A12 were upregulated in cALL.

CONCLUSION

TLR2, PTAFR, FCGR2A, RETN, S100A12 and CAT may serve as potential therapeutic targets.

Transcriptome-wide association study of alternative polyadenylation identifies susceptibility genes in non-small cell lung cancer

Alternative polyadenylation (APA) plays a crucial role in cancer development and prognosis. However, the molecular characteristics of APA related to non-small cell lung cancer (NSCLC) susceptibility remain understudied, especially in East Asian populations. In this study, we constructed an atlas of APA-regulated 3' untranslated region (3'UTR) and profiled its genetic regulation in 747 lung tissue samples (including tumors and paired normal tissues) from 417 NSCLC Chinese patients. We verified a significant global shortening of 3'UTRs in tumor samples compared to normal samples and underscored the value of APA-regulation as a prognostic marker. The 3'UTR APA quantitative trait loci (3'aQTL) was identified by regressing the percentage of distal poly(A) site usage index (PDUI) value on genetic variants. We found that a significant proportion 3'aQTLs are independent of genetic regulation of expression and are specific in Chinese. We also conducted a 3'UTR APA transcriptome-wide association study (3'aTWAS) by integrating the APA regulation atlas with a genome-wide association study (GWAS) for NSCLC involving 7035 cases and 185,413 cancer-free controls. We identified NSCLC-associated genes, highlighting TUBB, TEAD3, and PPP1R10. Combining the consistent results from colocalization analysis, differential APA analysis, and survival analysis, we provide novel evidence for the role TUBB APA regulation in NSCLC and identified potential upstream regulators. Overall, our study profiled the APA regulation and highlighted the substantial role of APA in NSCLC carcinogenesis and prognosis in East Asian populations.

A framework for studying multi-omic risk factors and their interplay: application to coronary artery disease

Background

Transcriptome-wide and proteome-wide association studies (TWAS and PWAS) have identified risk genes across complex diseases. However, the contributions of proximal risk variants and cross-omic interplay remain to be understood.

Methods

We propose an integrative framework to characterize disease-associated transcripts and proteins, and apply it to coronary artery disease (CAD). We employed S-PrediXcan on a large-scale genome-wide association study (GWAS) for CAD, with prediction models from GTEx v8 whole blood tissue and the Atherosclerosis Risk in Communities (ARIC) plasma protein. Conditional analyses adjusting for nearby CAD risk variants were performed to retain significant associations. Genes identified by both TWAS and PWAS were subsequently examined for associations with CAD risk factors, and colocalization analyses were performed for expression quantitative trait loci (eQTLs) and protein QTLs (pQTLs).

Results

TWAS identified 294 genes, and PWAS identified 79 genes, with 10 genes common between two analyses: CHMP2B, CLIC4, IL6R, MIF, MXRA7, NME2, NUDT5, PCSK9, TAGLN2 , and WARS . The predicted transcripts and proteins of these genes exhibited consistent associations with CAD risk factors, and colocalization tests revealed shared signals between eQTLs and pQTLs.

Conclusion

Our summary-level data framework enables the construction of multi-omic risk profiles for CAD, advancing the understanding of its genetic etiology through integrated transcriptomic and proteomic analyses.

A Multi-omic study integrating plasma protein, multiple tissues, and single-cell identifies RNASET2 as a key gene for lung cancer

INTRODUCTION

Lung cancer (LC) has the highest cancer-related mortality rate. Even though genome-wide association studies (GWAS) have identified numerous loci linked to LC risk, the underlying causal genes and biological processes are still mostly unknown.

METHODS

The LC GWAS summary data comprised 29,863 cases and 55,586 controls of European ancestry. The weight file and related files of plasma protein, multi-tissue, and single-cell were obtained from Zhang's study, Mancuso lab, and Thompson's study, respectively. We conducted transcriptome association studies (TWAS) employing functional Summary-based Imputation (FUSION) from two levels, which were multiple tissues and single cell. We conducted proteome-wide association studies (PWAS) from plasma protein. Conditional and joint (COJO) analysis and multi-marker analysis of genomic annotation (MAGMA) analysis were used to further screen the PWAS/TWAS results. Summary-data-based Mendelian randomization (SMR) and colocalization analysis were utilized to explain the causal association between variables and results.

RESULTS

A total of 13, 251, and 16 genes were calculated from the three dimensions, which were plasma protein, multiple tissues, and single cell, respectively. RNASET2 and IREB2 were obtained through intersecting these three sets of genes. COJO analysis and MAGMA analysis were replicated the two genes successfully. Then, RNASET2 was replicated in both eQTL-SMR and mQTL-SMR and following colocalization analysis.

CONCLUSION

In summary, we conducted a multi-omic studies, which integrated three levels to investigate the novel targets for LC. Through a series of verifications, RNASET2 was identified as the key gene for LC in the current research.

Leveraging genome-wide association studies to better understand the etiology of cancers

Genome-wide association studies (GWAS) statistically assess the association between tens of millions of genetic variants in the whole genome and a phenotype of interest. Genome-wide association studies enable the elucidation of polygenic inheritance of cancer, in which myriad low-penetrance genetic variants collectively contribute to a substantial proportion of the heritable susceptibility. In addition to the robust genotype-phenotype associations provided by GWAS, combining GWAS data with functional genomic datasets or sophisticated statistical genetic methods unlocks deeper insights. Integrating genotype and molecular phenotyping data facilitates functional characterization of GWAS association signals through molecular quantitative trait loci mapping and transcriptome-wide association studies. Furthermore, aggregating genome-wide polygenic signals, including subthreshold associations, enables one to estimate genetic correlations across diverse phenotypes and helps in clinical risk predictions by evaluating polygenic risk scores. In this review, we begin by summarizing the rationale for GWAS of cancer, introduce recent methodological updates in the GWAS-derived downstream analyses, and demonstrate their applications to GWAS of cancers.

Enhancing disease risk gene discovery by integrating transcription factor-linked trans-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-variants for predicting gene expression remain largely unexplored. Here, we introduce transTF-TWAS, which incorporates transcription factor (TF)-linked trans-variants to enhance model building for TF downstream target genes. Using data from the Genotype-Tissue Expression project, we predict gene expression and alternative splicing and applied these prediction models to large GWAS datasets for breast, prostate, lung cancers and other diseases. We demonstrate that transTF-TWAS outperforms other existing TWAS approaches in both constructing gene expression prediction models and identifying disease-associated genes, as shown by simulations and real data analysis. Our transTF-TWAS approach significantly contributes to the discovery of disease risk genes. Findings from this study shed new light on several genetically driven key TF regulators and their associated TF-gene regulatory networks underlying disease susceptibility.

webTWAS 2.0: update platform for identifying complex disease susceptibility genes through transcriptome-wide association study

Transcriptome-wide association study (TWAS) has successfully identified numerous complex disease susceptibility genes in the post-genome-wide association study (GWAS) era. Over the past 3 years, the focus of TWAS algorithms has shifted from merely identifying associations to understanding how single nucleotide polymorphisms (SNPs) regulate gene expression, with a growing emphasis on incorporating fine-mapping techniques. Additionally, the rapid increase in GWAS summary statistics, driven largely by the UK Biobank and other consortia, has made it essential to update our webTWAS resource. To address these challenges and meet the growing needs of researchers, we developed webTWAS 2.0, an updated platform for identifying susceptibility genes for human complex diseases using TWAS. Additionally, webTWAS 2.0 provides an online TWAS analysis tool that simplifies conducting TWAS analyses. The updated resource includes 7247 GWAS summary statistics covering 1588 complex human diseases from 192 publications. It also incorporates multiple TWAS methods, such as sTF-TWAS, 3'aTWAS and GIFT, along with an updated interactive visualization tool that allows users to easily explore significant associations across different methods. Other upgrades include a personalized online analysis tool for user-submitted GWAS data and a refined search function that makes it easier to identify relevant associations and meet diverse user needs more efficiently. webTWAS 2.0 is freely accessible at http://www.webtwas.net.

Refining breast cancer genetic risk and biology through multi-ancestry fine-mapping analyses of 192 risk regions

Genome-wide association studies have identified approximately 200 genetic risk loci for breast cancer, but the causal variants and target genes are mostly unknown. We sought to fine-map all known breast cancer risk loci using genome-wide association study data from 172,737 female breast cancer cases and 242,009 controls of African, Asian and European ancestry. We identified 332 independent association signals for breast cancer risk, including 131 signals not reported previously, and for 50 of them, we narrowed the credible causal variants down to a single variant. Analyses integrating functional genomics data identified 195 putative susceptibility genes, enriched in PI3K/AKT, TNF/NF-κB, p53 and Wnt/β-catenin pathways. Single-cell RNA sequencing or in vitro experiment data provided additional functional evidence for 105 genes. Our study uncovered large numbers of association signals and candidate susceptibility genes for breast cancer, uncovered breast cancer genetics and biology, and supported the value of including multi-ancestry data in fine-mapping analyses.

Identification of proteins associated with type 2 diabetes risk in diverse racial and ethnic populations

AIMS/HYPOTHESIS

Several studies have reported associations between specific proteins and type 2 diabetes risk in European populations. To better understand the role played by proteins in type 2 diabetes aetiology across diverse populations, we conducted a large proteome-wide association study using genetic instruments across four racial and ethnic groups: African; Asian; Hispanic/Latino; and European.

METHODS

Genome and plasma proteome data from the Multi-Ethnic Study of Atherosclerosis (MESA) study involving 182 African, 69 Asian, 284 Hispanic/Latino and 409 European individuals residing in the USA were used to establish protein prediction models by using potentially associated cis- and trans-SNPs. The models were applied to genome-wide association study summary statistics of 250,127 type 2 diabetes cases and 1,222,941 controls from different racial and ethnic populations.

RESULTS

We identified three, 44 and one protein associated with type 2 diabetes risk in Asian, European and Hispanic/Latino populations, respectively. Meta-analysis identified 40 proteins associated with type 2 diabetes risk across the populations, including well-established as well as novel proteins not yet implicated in type 2 diabetes development.

CONCLUSIONS/INTERPRETATION

Our study improves our understanding of the aetiology of type 2 diabetes in diverse populations.

DATA AVAILABILITY

The summary statistics of multi-ethnic type 2 diabetes GWAS of MVP, DIAMANTE, Biobank Japan and other studies are available from The database of Genotypes and Phenotypes (dbGaP) under accession number phs001672.v3.p1. MESA genetic, proteome and covariate data can be accessed through dbGaP under phs000209.v13.p3. All code is available on GitHub ( https://github.com/Arthur1021/MESA-1K-PWAS ).

A proteome-wide association study identifies putative causal proteins for breast cancer risk

BACKGROUND

Genome-wide association studies (GWAS) have identified more than 200 breast cancer risk-associated genetic loci, yet the causal genes and biological mechanisms for most loci remain elusive. Proteins, as final gene products, are pivotal in cellular function. In this study, we conducted a proteome-wide association study (PWAS) to identify proteins in breast tissue related to breast cancer risk.

METHODS

We profiled the proteome in fresh frozen breast tissue samples from 120 cancer-free European-ancestry women from the Susan G. Komen Tissue Bank (KTB). Protein expression levels were log2-transformed then normalized via quantile and inverse-rank transformations. GWAS data were also generated for these 120 samples. These data were used to build statistical models to predict protein expression levels via cis-genetic variants using the elastic net method. The prediction models were then applied to the GWAS summary statistics data of 133,384 breast cancer cases and 113,789 controls to assess the associations of genetically predicted protein expression levels with breast cancer risk overall and its subtypes using the S-PrediXcan method.

RESULTS

A total of 6388 proteins were detected in the normal breast tissue samples from 120 women with a high detection false discovery rate (FDR) p value < 0.01. Among the 5820 proteins detected in more than 80% of participants, prediction models were successfully built for 2060 proteins with R > 0.1 and P < 0.05. Among these 2060 proteins, five proteins were significantly associated with overall breast cancer risk at an FDR p value < 0.1. Among these five proteins, the corresponding genes for proteins COPG1, DCTN3, and DDX6 were located at least 1 Megabase away from the GWAS-identified breast cancer risk variants. COPG1 was associated with an increased risk of breast cancer with a p value of 8.54 × 10-4. Both DCTN3 and DDX6 were associated with a decreased risk of breast cancer with p values of 1.01 × 10-3 and 3.25 × 10-4, respectively. The corresponding genes for the remaining two proteins, LSP1 and DNAJA3, were located in previously GWAS-identified breast cancer risk loci. After adjusting for GWAS-identified risk variants, the association for DNAJA3 was still significant (p value of 9.15 × 10-5 and adjusted p value of 1.94 × 10-4). However, the significance for LSP1 became weaker with a p value of 0.62. Stratification analyses by breast cancer subtypes identified three proteins, SMARCC1, LSP1, and NCKAP1L, associated with luminal A, luminal B, and ER-positive breast cancer. NCKAP1L was located at least 1Mb away from the GWAS-identified breast cancer risk variants. After adjusting for GWAS-identified breast cancer risk variants, the association for protein LSP1 was still significant (adjusted p value of 6.43 × 10-3 for luminal B subtype).

CONCLUSION

We conducted the first breast-tissue-based PWAS and identified seven proteins associated with breast cancer, including five proteins not previously implicated. These findings help improve our understanding of the underlying genetic mechanism of breast cancer development.

Age-Related Hearing Impairment: Genome and Blood Methylome Data Integration Reveals Candidate Epigenetic Biomarkers

Age-related hearing impairment (ARHI) is a major planetary health burden that is in need of precision medicine for prevention, diagnosis, and treatment. The present study was set out to identify candidate epigenetic markers for ARHI. Associations of genetically predicted DNA methylation levels with ARHI risk were evaluated using two sets of blood DNA methylation genetic prediction models in 147,997 cases and 575,269 controls of European descent. A total of 1314 CpG sites (CpGs) were significantly associated with ARHI risk at a false discovery rate (FDR) <0.05, including 12 putatively causal CpGs based on fine-mapping analysis. Measured methylation levels of 247 of the associated CpGs were significantly correlated with measured expression levels of 127 nearby genes in blood at an FDR <0.05. A total of 37 CpGs and their 18 nearby genes showed consistent association directions for the methylation-gene expression-ARHI risk pathway. Importantly, three genes (PEX6, TCF19, and SPTBN1) were enriched in auditory disease categories. Our results indicate that specific CpGs may modulate ARHI risk by regulating the expression of candidate ARHI target genes. Future precision medicine and biomarker development research on ARHI are called for.

An atlas of genetic effects on the monocyte methylome across European and African populations

Elucidating the genetic architecture of DNA methylation is crucial for decoding complex disease etiology. However, current epigenomic studies are often limited by incomplete methylation coverage and heterogeneous tissue samples. Here, we present the first comprehensive, multi-ancestry human methylome atlas of purified human monocytes, generated through integrated whole-genome bisulfite sequencing and whole-genome sequencing from 298 European Americans (EA) and 160 African Americans (AA). By analyzing over 25 million methylation sites, we identified 1,383,250 and 1,721,167 methylation quantitative trait loci (meQTLs) in cis- regions for EA and AA populations, respectively, revealing both shared (880,108 sites) and population-specific regulatory patterns. Furthermore, we developed population-specific DNAm imputation models, enabling methylome-wide association studies (MWAS) for 1,976,046 and 2,657,581 methylation sites in EA and AA, respectively. These models were validated through multi-ancestry analysis of 41 complex traits from the Million Veteran Program. The identified meQTLs, MWAS models, and data resources are freely available at www.gcbhub.org and https://osf.io/gct57/ .

Cross-population enhancement of PrediXcan predictions with a gnomAD-based east Asian reference framework

Over the past decade, genome-wide association studies have identified thousands of variants significantly associated with complex traits. For each locus, gene expression levels are needed to further explore its biological functions. To address this, the PrediXcan algorithm leverages large-scale reference data to impute the gene expression level from single nucleotide polymorphisms, and thus the gene-trait associations can be tested to identify the candidate causal genes. However, a challenge arises due to the fact that most reference data are from subjects of European ancestry, and the accuracy and robustness of predicted gene expression in subjects of East Asian (EAS) ancestry remains unclear. Here, we first simulated a variety of scenarios to explore the impact of the level of population diversity on gene expression. Population differentiated variants were estimated by using the allele frequency information from The Genome Aggregation Database. We found that the weights of a variants was the main factor that affected the gene expression predictions, and that ~70% of variants were significantly population differentiated based on proportion tests. To provide insights into this population effect on gene expression levels, we utilized the allele frequency information to develop a gene expression reference panel, Predict Asian-Population (PredictAP), for EAS ancestry. PredictAP can be viewed as an auxiliary tool for PrediXcan when using genotype data from EAS subjects.

A bootstrap model comparison test for identifying genes with context-specific patterns of genetic regulation

Understanding how genetic variation affects gene expression is essential for a complete picture of the functional pathways that give rise to complex traits. Although numerous studies have established that many genes are differentially expressed in distinct human tissues and cell types, no tools exist for identifying the genes whose expression is differentially regulated. Here we introduce DRAB (differential regulation analysis by bootstrapping), a gene-based method for testing whether patterns of genetic regulation are significantly different between tissues or other biological contexts. DRAB first leverages the elastic net to learn context-specific models of local genetic regulation and then applies a novel bootstrap-based model comparison test to check their equivalency. Unlike previous model comparison tests, our proposed approach can determine whether population-level models have equal predictive performance by accounting for the variability of feature selection and model training. We validated DRAB on mRNA expression data from a variety of human tissues in the Genotype-Tissue Expression (GTEx) Project. DRAB yielded biologically reasonable results and had sufficient power to detect genes with tissue-specific regulatory profiles while effectively controlling false positives. By providing a framework that facilitates the prioritization of differentially regulated genes, our study enables future discoveries on the genetic architecture of molecular phenotypes.

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.

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.

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.

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.

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.

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.

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.

Large-scale integration of omics and electronic health records to identify potential risk protein biomarkers and therapeutic drugs for cancer prevention and intervention

Identifying risk protein targets and their therapeutic drugs is crucial for effective cancer prevention. Here, we conduct integrative and fine-mapping analyses of large genome-wide association studies data for breast, colorectal, lung, ovarian, pancreatic, and prostate cancers, and characterize 710 lead variants independently associated with cancer risk. Through mapping protein quantitative trait loci (pQTL) for these variants using plasma proteomics data from over 75,000 participants, we identify 365 proteins associated with cancer risk. Subsequent colocalization analysis identifies 101 proteins, including 74 not reported in previous studies. We further characterize 36 potential druggable proteins for cancers or other disease indications. Analyzing >3.5 million electronic health records, we uncover five drugs (Haloperidol, Trazodone, Tranexamic Acid, Haloperidol, and Captopril) associated with increased cancer risk and two drugs (Caffeine and Acetazolamide) linked to reduced colorectal cancer risk. This study offers novel insights into therapeutic drugs targeting risk proteins for cancer prevention and intervention.

Multi-Omic Analysis Reveals Genetic Determinants and Therapeutic Targets of Chronic Kidney Disease and Kidney Function

Chronic kidney disease (CKD) presents a significant global health challenge, characterized by complex pathophysiology. This study utilized a multi-omic approach, integrating genomic data from the CKDGen consortium alongside transcriptomic, metabolomic, and proteomic data to elucidate the genetic underpinnings and identify therapeutic targets for CKD and kidney function. We employed a range of analytical methods including cross-tissue transcriptome-wide association studies (TWASs), Mendelian randomization (MR), summary-based MR (SMR), and molecular docking. These analyses collectively identified 146 cross-tissue genetic associations with CKD and kidney function. Key Golgi apparatus-related genes (GARGs) and 41 potential drug targets were highlighted, with MAP3K11 emerging as a significant gene from the TWAS and MR data, underscoring its potential as a therapeutic target. Capsaicin displayed promising drug-target interactions in molecular docking analyses. Additionally, metabolome- and proteome-wide MR (PWMR) analyses revealed 33 unique metabolites and critical inflammatory proteins such as FGF5 that are significantly linked to and colocalized with CKD and kidney function. These insights deepen our understanding of CKD pathogenesis and highlight novel targets for treatment and prevention.

SR-TWAS: Leveraging Multiple Reference Panels to Improve TWAS 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 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 showed that SR-TWAS improved power, due to increased effective 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 application studies identified 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 were found for these significant risk genes.

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.

DNA methylation marker identification and poly-methylation risk score in prediction of healthspan termination

Aim: To elucidate the epigenetic consequences of DNA methylation in healthspan termination (HST), considering the current limited understanding. Materials & methods: Genetically predicted DNA methylation models were established (n = 2478). These models were applied to genome-wide association study data on HST. Then, a poly-methylation risk score (PMRS) was established in 241,008 individuals from the UK Biobank. Results: Of the 63,046 CpGs from the prediction models, 13 novel CpGs were associated with HST. Furthermore, people with high PMRSs showed higher HST risk (hazard ratio: 1.18; 95% CI: 1.13-1.25). Conclusion: The study indicates that DNA methylation may influence HST by regulating the expression of genes (e.g., PRMT6, CTSK). PMRSs have a promising application in discriminating subpopulations to facilitate early prevention.

Expression- and splicing-based multi-tissue transcriptome-wide association studies identified multiple genes for breast cancer by estrogen-receptor status

BACKGROUND

Although several transcriptome-wide association studies (TWASs) have been performed to identify genes associated with overall breast cancer (BC) risk, only a few TWAS have explored the differences in estrogen receptor-positive (ER+) and estrogen receptor-negative (ER-) breast cancer. Additionally, these studies were based on gene expression prediction models trained primarily in breast tissue, and they did not account for alternative splicing of genes.

METHODS

In this study, we utilized two approaches to perform multi-tissue TWASs of breast cancer by ER subtype: (1) an expression-based TWAS that combined TWAS signals for each gene across multiple tissues and (2) a splicing-based TWAS that combined TWAS signals of all excised introns for each gene across tissues. To perform this TWAS, we utilized summary statistics for ER + BC from the Breast Cancer Association Consortium (BCAC) and for ER- BC from a meta-analysis of BCAC and the Consortium of Investigators of Modifiers of BRCA1 and BRCA2 (CIMBA).

RESULTS

In total, we identified 230 genes in 86 loci that were associated with ER + BC and 66 genes in 29 loci that were associated with ER- BC at a Bonferroni threshold of significance. Of these genes, 2 genes associated with ER + BC at the 1q21.1 locus were located at least 1 Mb from published GWAS hits. For several well-studied tumor suppressor genes such as TP53 and CHEK2 which have historically been thought to impact BC risk through rare, penetrant mutations, we discovered that common variants, which modulate gene expression, may additionally contribute to ER + or ER- etiology.

CONCLUSIONS

Our study comprehensively examined how differences in common variation contribute to molecular differences between ER + and ER- BC and introduces a novel, splicing-based framework that can be used in future TWAS studies.

Regulome-wide association study identifies genetically driven accessible regions associated with pancreatic cancer risk

Genome-wide association studies (GWAS) have identified two dozen genetic variants that are associated with the risk of pancreatic ductal adenocarcinoma (PDAC), a deadly malignancy. However, a majority of these variants are located in noncoding regions of the genome, which limits the translation of GWAS findings into clinical applications. The regulome-wide association study (RWAS) is a recently developed method for identifying TF binding-induced accessibility regions for diseases. However, their potential connection to PDAC has yet to be fully explored. We evaluated the associations between genetically predicted levels of chromatin accessibility and risk of PDAC by using pan-cancer chromatin accessibility genetic prediction models. Our analysis included 8275 cases and 6723 controls from the PanScan (I, II, and III) and PanC4 consortia. To further refine our results, we also integrated genes associated to allele-specific accessibility quantitative trait loci (as-aQTL) and TF motifs located in the as-aQTL. We found that 50 chromatin accessibility features were associated with PDAC risk at a false discovery rate (FDR) of less than 0.05. A total of 28 RWAS peaks were identified as conditionally significant. By integrating the results from as-aQTL, motif analysis, and RWAS, we identified candidate causal regulatory elements for two potential chromatin accessibility regions (THCA_89956 and ESCA_89167) that are associated with PDAC risk. Our study identified chromatin accessibility features in noncoding genomic regions that are associated with PDAC risk. We also predicted the associated genes and disrupt motifs. Our findings provide new insights into the regulatory mechanisms of noncoding regions for pancreatic tumorigenesis.

rvTWAS: identifying gene-trait association using sequences by utilizing transcriptome-directed feature selection

Toward the identification of genetic basis of complex traits, transcriptome-wide association study (TWAS) is successful in integrating transcriptome data. However, TWAS is only applicable for common variants, excluding rare variants in exome or whole-genome sequences. This is partly because of the inherent limitation of TWAS protocols that rely on predicting gene expressions. Our previous research has revealed the insight into TWAS: the 2 steps in TWAS, building and applying the expression prediction models, are essentially genetic feature selection and aggregations that do not have to involve predictions. Based on this insight disentangling TWAS, rare variants' inability of predicting expression traits is no longer an obstacle. Herein, we developed "rare variant TWAS," or rvTWAS, that first uses a Bayesian model to conduct expression-directed feature selection and then uses a kernel machine to carry out feature aggregation, forming a model leveraging expressions for association mapping including rare variants. We demonstrated the performance of rvTWAS by thorough simulations and real data analysis in 3 psychiatric disorders, namely schizophrenia, bipolar disorder, and autism spectrum disorder. We confirmed that rvTWAS outperforms existing TWAS protocols and revealed additional genes underlying psychiatric disorders. Particularly, we formed a hypothetical mechanism in which zinc finger genes impact all 3 disorders through transcriptional regulations. rvTWAS will open a door for sequence-based association mappings integrating gene expressions.

Big data and artificial intelligence in cancer research

The field of oncology has witnessed an extraordinary surge in the application of big data and artificial intelligence (AI). AI development has made multiscale and multimodal data fusion and analysis possible. A new era of extracting information from complex big data is rapidly evolving. However, challenges related to efficient data curation, in-depth analysis, and utilization remain. We provide a comprehensive overview of the current state of the art in big data and computational analysis, highlighting key applications, challenges, and future opportunities in cancer research. By sketching the current landscape, we seek to foster a deeper understanding and facilitate the advancement of big data utilization in oncology, call for interdisciplinary collaborations, ultimately contributing to improved patient outcomes and a profound understanding of cancer.

A large-scale microRNA transcriptome-wide association study identifies two susceptibility microRNAs, miR-1307-5p and miR-192-3p, for colorectal cancer risk

Transcriptome-wide association studies (TWAS) have identified many putative susceptibility genes for colorectal cancer (CRC) risk. However, susceptibility miRNAs, critical dysregulators of gene expression, remain unexplored. We genotyped DNA samples from 313 CRC East Asian patients and performed small RNA sequencing in their normal colon tissues distant from tumors to build genetic models for predicting miRNA expression. We applied these models and data from genome-wide association studies (GWAS) including 23 942 cases and 217 267 controls of East Asian ancestry to investigate associations of predicted miRNA expression with CRC risk. Perturbation experiments separately by promoting and inhibiting miRNAs expressions and further in vitro assays in both SW480 and HCT116 cells were conducted. At a Bonferroni-corrected threshold of P < 4.5 × 10-4, we identified two putative susceptibility miRNAs, miR-1307-5p and miR-192-3p, located in regions more than 500 kb away from any GWAS-identified risk variants in CRC. We observed that a high predicted expression of miR-1307-5p was associated with increased CRC risk, while a low predicted expression of miR-192-3p was associated with increased CRC risk. Our experimental results further provide strong evidence of their susceptible roles by showing that miR-1307-5p and miR-192-3p play a regulatory role, respectively, in promoting and inhibiting CRC cell proliferation, migration, and invasion, which was consistently observed in both SW480 and HCT116 cells. Our study provides additional insights into the biological mechanisms underlying CRC development.

Novel insights into genetic susceptibility for colorectal cancer from transcriptome-wide association and functional investigation

BACKGROUND

Transcriptome-wide association studies have been successful in identifying candidate susceptibility genes for colorectal cancer (CRC). To strengthen susceptibility gene discovery, we conducted a large transcriptome-wide association study and an alternative splicing transcriptome-wide association study in CRC using improved genetic prediction models and performed in-depth functional investigations.

METHODS

We analyzed RNA-sequencing data from normal colon tissues and genotype data from 423 European descendants to build genetic prediction models of gene expression and alternative splicing and evaluated model performance using independent RNA-sequencing data from normal colon tissues of the Genotype-Tissue Expression Project. We applied the verified models to genome-wide association studies (GWAS) summary statistics among 58 131 CRC cases and 67 347 controls of European ancestry to evaluate associations of genetically predicted gene expression and alternative splicing with CRC risk. We performed in vitro functional assays for 3 selected genes in multiple CRC cell lines.

RESULTS

We identified 57 putative CRC susceptibility genes, which included the 48 genes from transcriptome-wide association studies and 15 genes from splicing transcriptome-wide association studies, at a Bonferroni-corrected P value less than .05. Of these, 16 genes were not previously implicated in CRC susceptibility, including a gene PDE7B (6q23.3) at locus previously not reported by CRC GWAS. Gene knockdown experiments confirmed the oncogenic roles for 2 unreported genes, TRPS1 and METRNL, and a recently reported gene, C14orf166.

CONCLUSION

This study discovered new putative susceptibility genes of CRC and provided novel insights into the biological mechanisms underlying CRC development.

A statistical method for image-mediated association studies discovers genes and pathways associated with four brain disorders

Brain imaging and genomics are critical tools enabling characterization of the genetic basis of brain disorders. However, imaging large cohorts is expensive and may be unavailable for legacy datasets used for genome-wide association studies (GWASs). Using an integrated feature selection/aggregation model, we developed an image-mediated association study (IMAS), which utilizes borrowed imaging/genomics data to conduct association mapping in legacy GWAS cohorts. By leveraging the UK Biobank image-derived phenotypes (IDPs), the IMAS discovered genetic bases underlying four neuropsychiatric disorders and verified them by analyzing annotations, pathways, and expression quantitative trait loci (eQTLs). A cerebellar-mediated mechanism was identified to be common to the four disorders. Simulations show that, if the goal is identifying genetic risk, our IMAS is more powerful than a hypothetical protocol in which the imaging results were available in the GWAS dataset. This implies the feasibility of reanalyzing legacy GWAS datasets without conducting additional imaging, yielding cost savings for integrated analysis of genetics and imaging.

Splicing transcriptome-wide association study to identify splicing events for pancreatic cancer risk

A large proportion of the heritability of pancreatic cancer risk remains elusive, and the contribution of specific mRNA splicing events to pancreatic cancer susceptibility has not been systematically evaluated. In this study, we performed a large splicing transcriptome-wide association study (spTWAS) using three modeling strategies (Enet, LASSO and MCP) to develop alternative splicing genetic prediction models for identifying novel susceptibility loci and splicing introns for pancreatic cancer risk by assessing 8275 pancreatic cancer cases and 6723 controls of European ancestry. Data from 305 subjects of whom the majority are of European descent in the Genotype-Tissue Expression Project (GTEx) were used and both cis-acting and promoter-enhancer interaction regions were considered to build these models. We identified nine splicing events of seven genes (ABO, UQCRC1, STARD3, ETAA1, CELA3B, LGR4 and SFT2D1) that showed an association of genetically predicted expression with pancreatic cancer risk at a false discovery rate ≤0.05. Of these genes, UQCRC1 and LGR4 have not yet been reported to be associated with pancreatic cancer risk. Fine-mapping analyses supported likely causal associations corresponding to six splicing events of three genes (P4HTM, ABO and PGAP3). Our study identified novel genes and splicing events associated with pancreatic cancer risk, which can improve our understanding of the etiology of this deadly malignancy.

Identification of candidate DNA methylation biomarkers related to Alzheimer's disease risk by integrating genome and blood methylome data

Alzheimer disease (AD) is a common neurodegenerative disease with a late onset. It is critical to identify novel blood-based DNA methylation biomarkers to better understand the extent of the molecular pathways affected in AD. Two sets of blood DNA methylation genetic prediction models developed using different reference panels and modelling strategies were leveraged to evaluate associations of genetically predicted DNA methylation levels with AD risk in 111,326 (46,828 proxy) cases and 677,663 controls. A total of 1,168 cytosine-phosphate-guanine (CpG) sites showed a significant association with AD risk at a false discovery rate (FDR) < 0.05. Methylation levels of 196 CpG sites were correlated with expression levels of 130 adjacent genes in blood. Overall, 52 CpG sites of 32 genes showed consistent association directions for the methylation-gene expression-AD risk, including nine genes (CNIH4, THUMPD3, SERPINB9, MTUS1, CISD1, FRAT2, CCDC88B, FES, and SSH2) firstly reported as AD risk genes. Nine of 32 genes were enriched in dementia and AD disease categories (P values ranged from 1.85 × 10-4 to 7.46 × 10-6), and 19 genes in a neurological disease network (score = 54) were also observed. Our findings improve the understanding of genetics and etiology for AD.

RCCD1 promotes breast carcinogenesis through regulating hypoxia-associated mitochondrial homeostasis

Regulator of chromosome condensation domain-containing protein 1 (RCCD1), previously reported as a partner of histone H3K36 demethylase KDM8 involved in chromosome segregation, has been identified as a potential driver for breast cancer in a recent transcriptome-wide association study. We report here that, unexpectedly, RCCD1 is also localized in mitochondria. We show that RCCD1 resides in the mitochondrial matrix, where it interacts with the mitochondrial contact site/cristae organizing system (MICOS) and mitochondrial DNA (mtDNA) to regulate mtDNA transcription, oxidative phosphorylation, and the production of reactive oxygen species. Interestingly, RCCD1 is upregulated under hypoxic conditions, leading to decreased generation of reactive oxygen species and alleviated apoptosis favoring cancer cell survival. We show that RCCD1 promotes breast cancer cell proliferation in vitro and accelerates breast tumor growth in vivo. Indeed, RCCD1 is overexpressed in breast carcinomas, and its level of expression is associated with aggressive breast cancer phenotypes and poor patient survival. Our study reveals an additional dimension of RCCD1 functionality in regulating mitochondrial homeostasis, whose dysregulation inflicts pathologic states such as breast cancer.

Isoform-level transcriptome-wide association uncovers genetic risk mechanisms for neuropsychiatric disorders in the human brain

Methods integrating genetics with transcriptomic reference panels prioritize risk genes and mechanisms at only a fraction of trait-associated genetic loci, due in part to an overreliance on total gene expression as a molecular outcome measure. This challenge is particularly relevant for the brain, in which extensive splicing generates multiple distinct transcript-isoforms per gene. Due to complex correlation structures, isoform-level modeling from cis-window variants requires methodological innovation. Here we introduce isoTWAS, a multivariate, stepwise framework integrating genetics, isoform-level expression and phenotypic associations. Compared to gene-level methods, isoTWAS improves both isoform and gene expression prediction, yielding more testable genes, and increased power for discovery of trait associations within genome-wide association study loci across 15 neuropsychiatric traits. We illustrate multiple isoTWAS associations undetectable at the gene-level, prioritizing isoforms of AKT3, CUL3 and HSPD1 in schizophrenia and PCLO with multiple disorders. Results highlight the importance of incorporating isoform-level resolution within integrative approaches to increase discovery of trait associations, especially for brain-relevant traits.

Cis- and trans-eQTL TWAS of breast and ovarian cancer identify more than 100 risk associated genes in the BCAC and OCAC consortia

Transcriptome-wide association studies (TWAS) have investigated the role of genetically regulated transcriptional activity in the etiologies of breast and ovarian cancer. However, methods performed to date have only considered 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 TWAS 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 a first look into the role of trans-eQTLs in the complex molecular mechanisms underlying these diseases.

Large-scale alternative polyadenylation (APA)-wide association studies to identify putative susceptibility genes in human common cancers

Alternative polyadenylation (APA) modulates mRNA processing in the 3' untranslated regions (3'UTR), which affect mRNA stability and translation efficiency. Here, we build genetic models to predict APA levels in multiple tissues using sequencing data of 1,337 samples from the Genotype-Tissue Expression, and apply these models to assess associations between genetically predicted APA levels and cancer risk with data from large genome-wide association studies of six common cancers, including breast, ovary, prostate, colorectum, lung, and pancreas among European-ancestry populations. At a Bonferroni-corrected P □<□0.05, we identify 58 risk genes, including seven in newly identified loci. Using luciferase reporter assays, we demonstrate that risk alleles of 3'UTR variants, rs324015 ( STAT6 ), rs2280503 ( DIP2B ), rs1128450 ( FBXO38 ) and rs145220637 ( LDAH ), could significantly increase post-transcriptional activities of their target genes compared to reference alleles. Further gene knockdown experiments confirm their oncogenic roles. Our study provides additional insight into the genetic susceptibility of these common cancers.

A BOOTSTRAP MODEL COMPARISON TEST FOR IDENTIFYING GENES WITH CONTEXT-SPECIFIC PATTERNS OF GENETIC REGULATION

Understanding how genetic variation affects gene expression is essential for a complete picture of the functional pathways that give rise to complex traits. Although numerous studies have established that many genes are differentially expressed in distinct human tissues and cell types, no tools exist for identifying the genes whose expression is differentially regulated. Here we introduce DRAB (Differential Regulation Analysis by Bootstrapping), a gene-based method for testing whether patterns of genetic regulation are significantly different between tissues or other biological contexts. DRAB first leverages the elastic net to learn context-specific models of local genetic regulation and then applies a novel bootstrap-based model comparison test to check their equivalency. Unlike previous model comparison tests, our proposed approach can determine whether population-level models have equal predictive performance by accounting for the variability of feature selection and model training. We validated DRAB on mRNA expression data from a variety of human tissues in the Genotype-Tissue Expression (GTEx) Project. DRAB yielded biologically reasonable results and had sufficient power to detect genes with tissue-specific regulatory profiles while effectively controlling false positives. By providing a framework that facilitates the prioritization of differentially regulated genes, our study enables future discoveries on the genetic architecture of molecular phenotypes.

Systematic analyses of GWAS summary statistics from UK Biobank identified novel susceptibility loci and genes for upper gastrointestinal diseases

In recent decades, upper gastrointestinal (GI) diseases have been highly prevalent worldwide. Although genome-wide association studies (GWASs) have identified thousands of susceptibility loci, only a few of them were conducted for chronic upper GI disorders, and most of them were underpowered and with small sample sizes. Additionally, for the known loci, only a tiny fraction of heritability can be explained and the underlying mechanisms and related genes remain unclear. In this study, we conducted a multi-trait analysis by the MTAG software and a two-stage transcriptome-wide association study (TWAS) with UTMOST and FUSION for seven upper GI diseases (oesophagitis, gastro-oesophageal reflux disease, other diseases of oesophagus, gastric ulcer, duodenal ulcer, gastritis and duodenitis and other diseases of stomach and duodenum) based on summary GWAS statistics from UK Biobank. In the MTAG analysis, we identified 7 loci associated with these upper GI diseases, including 3 novel ones at 4p12 (rs10029980), 12q13.13 (rs4759317) and 18p11.32 (rs4797954). In the TWAS analysis, we revealed 5 susceptibility genes in known loci and identified 12 novel potential susceptibility genes, including HOXC9 at 12q13.13. Further functional annotations and colocalization analysis indicated that rs4759317 (A>G) driven the association for GWAS signals and expression quantitative trait loci (eQTL) simultaneously at 12q13.13. The identified variant acted by decreasing the expression of HOXC9 to affect the risk of gastro-oesophageal reflux disease. This study provided insights into the genetic nature of upper GI diseases.

Challenges and Opportunities for Data Science in Women's Health

The intersection of women's health and data science is a field of research that has historically trailed other fields, but more recently it has gained momentum. This growth is being driven not only by new investigators who are moving into this area but also by the significant opportunities that have emerged in new methodologies, resources, and technologies in data science. Here, we describe some of the resources and methods being used by women's health researchers today to meet challenges in biomedical data science. We also describe the opportunities and limitations of applying these approaches to advance women's health outcomes and the future of the field, with emphasis on repurposing existing methodologies for women's health.

A Splicing Transcriptome-Wide Association Study Identifies Candidate Altered Splicing for Prostate Cancer Risk

Prostate cancer (PCa) represents a huge public health burden among men. Many susceptibility genetic factors for PCa still remain unknown. In this study, we performed a large splicing transcriptome-wide association study (spTWAS) using three modeling strategies to develop alternative splicing genetic prediction models for identifying novel susceptibility loci and splicing introns for PCa risk by assessing 79,194 cases and 61,112 controls of European ancestry in the PRACTICAL, CRUK, CAPS, BPC3, and PEGASUS consortia. We identified 120 splicing introns of 97 genes showing an association with PCa risk at false discovery rate (FDR)-corrected threshold (FDR <0.05). Of them, 33 genes were enriched in PCa-related diseases and function categories. Fine-mapping analysis suggested that 21 splicing introns of 19 genes were likely causally associated with PCa risk. Thirty-five splicing introns of 34 novel genes were identified to be related to PCa susceptibility for the first time, and 11 of the genes were enriched in a cancer-related network. Our study identified novel loci and splicing introns associated with PCa risk, which can improve our understanding of the etiology of this common malignancy.

Clinical significance and potential application of cuproptosis-related genes in gastric cancer

BACKGROUND

Worldwide, gastric cancer (GC) is a common lethal solid malignancy with a poor prognosis. Cuproptosis is a novel type of cell death mediated by protein lipoylation and may be related to GC prognosis.

AIM

To offer new insights to predict GC prognosis and provide multiple therapeutic targets related to cuproptosis-related genes (CRGs) for future therapy.

METHODS

We collected data from several public data portals, systematically estimated the expression level and prognostic values of CRGs in GC samples, and investigated related mechanisms using public databases and bioinformatics.

RESULTS

Our results revealed that FDX1, LIAS, and MTF1 were differentially expressed in GC samples and exhibited important prognostic significance in The Cancer Genome Atlas (TCGA) cohort. We constructed a nomogram model for overall survival and disease-specific survival prediction and validated it via calibration plots. Mecha-nistically, immune cell infiltration and DNA methylation prominently affected the survival time of GC patients. Moreover, protein-protein interaction network, KEGG pathway and gene ontology enrichment analyses demonstrated that FDX1, LIAS, MTF1 and related proteins play key roles in the tricarboxylic acid cycle and cuproptosis. Gene Expression Omnibus database validation showed that the expression levels of FDX1, LIAS, and MTF1 were consistent with those in the TCGA cohort. Top 10 perturbagens has been filtered by Connectivity Map.

CONCLUSION

In conclusion, FDX1, LIAS, and MTF1 could serve as potential prognostic biomarkers for GC patients and provide novel targets for immunotarget therapy.

A splicing transcriptome-wide association study identifies novel altered splicing for Alzheimer's disease susceptibility

Alzheimer's disease (AD) is a common neurodegenerative disease in aging individuals. Alternative splicing is reported to be relevant to AD development while their roles in etiology of AD remain largely elusive. We performed a comprehensive splicing transcriptome-wide association study (spTWAS) using intronic excision expression genetic prediction models of 12 brain tissues developed through three modelling strategies, to identify candidate susceptibility splicing introns for AD risk. A total of 111,326 (46,828 proxy) cases and 677,663 controls of European ancestry were studied. We identified 343 associations of 233 splicing introns (143 genes) with AD risk after Bonferroni correction (0.05/136,884 = 3.65 × 10^-7). Fine-mapping analyses supported 155 likely causal associations corresponding to 83 splicing introns of 55 genes. Eighteen causal splicing introns of 15 novel genes (EIF2D, WDR33, SAP130, BYSL, EPHB6, MRPL43, VEGFB, PPP1R13B, TLN2, CLUHP3, LRRC37A4P, CRHR1, LINC02210, ZNF45-AS1, and XPNPEP3) were identified for the first time to be related to AD susceptibility. Our study identified novel genes and splicing introns associated with AD risk, which can improve our understanding of the etiology of AD.

Identification of target proteins for breast cancer genetic risk loci and blood risk biomarkers in a large study by integrating genomic and proteomic data

Genome-wide association studies (GWAS) have identified around 200 loci associated with breast cancer risk. However, protein targets for these loci remain largely unknown. Identifying protein targets and biomarkers can improve the understanding of cancer biology and etiology and identify high-risk individuals for cancer prevention. In this study, we investigated genetically predicted levels of 1142 circulating proteins with breast cancer risk in 133 384 cases and 113 789 controls of European ancestry included in the Breast Cancer Association Consortium (BCAC). We identified 22 blood protein biomarkers associated with the risk of overall breast cancer at a false discovery rate (FDR) <0.05, including nine proteins encoded by genes located at least 500 kb away from previously reported risk variants for breast cancer. Analyses focusing on 124 encoding genes located at GWAS-identified breast cancer risk loci found 20 proteins associated with overall breast cancer risk and one protein associated with triple-negative breast cancer risk at FDR <0.05. Adjustment for the GWAS-identified risk variants significantly attenuated the association for 13 of these proteins, suggesting that these proteins may be the targets of these GWAS-identified risk loci. The identified proteins are involved in various biological processes, including glutathione conjugation, STAT5 signaling and NF-κB signaling pathways. Our study identified novel protein targets and risk biomarkers for breast cancer risk.

Update Breast Cancer 2023 Part 1 - Early Stage Breast Cancer

With abemaciclib (monarchE study) and olaparib (OlympiA study) gaining approval in the adjuvant treatment setting, a significant change in the standard of care for patients with early stage breast cancer has been established for some time now. Accordingly, some diverse developments are slowly being transferred from the metastatic to the adjuvant treatment setting. Recently, there have also been positive reports of the NATALEE study. Other clinical studies are currently investigating substances that are already established in the metastatic setting. These include, for example, the DESTINY Breast05 study with trastuzumab deruxtecan and the SASCIA study with sacituzumab govitecan. In this review paper, we summarize and place in context the latest developments over the past months.

A joint transcriptome-wide association study across multiple tissues identifies candidate breast cancer susceptibility genes

Genome-wide association studies (GWASs) have identified more than 200 genomic loci for breast cancer risk, but specific causal genes in most of these loci have not been identified. In fact, transcriptome-wide association studies (TWASs) of breast cancer performed using gene expression prediction models trained in breast tissue have yet to clearly identify most target genes. To identify candidate genes, we performed a GWAS analysis in a breast cancer dataset from UK Biobank (UKB) and combined the results with the GWAS results of the Breast Cancer Association Consortium (BCAC) by a meta-analysis. Using the summary statistics from the meta-analysis, we performed a joint TWAS analysis that combined TWAS signals from multiple tissues. We used expression prediction models trained in 11 tissues that are potentially relevant to breast cancer from the Genotype-Tissue Expression (GTEx) data. In the GWAS analysis, we identified eight loci distinct from those reported previously. In the TWAS analysis, we identified 309 genes at 108 genomic loci to be significantly associated with breast cancer at the Bonferroni threshold. Of these, 17 genes were located in eight regions that were at least 1 Mb away from published GWAS hits. The remaining TWAS-significant genes were located in 100 known genomic loci from previous GWASs of breast cancer. We found that 21 genes located in known GWAS loci remained statistically significant after conditioning on previous GWAS index variants. Our study provides insights into breast cancer genetics through mapping candidate target genes in a large proportion of known GWAS loci and discovering multiple new loci.

Identifying Biomarkers Using Support Vector Machine to Understand the Racial Disparity in Triple-Negative Breast Cancer

With the properties of aggressive cancer and heterogeneous tumor biology, triple-negative breast cancer (TNBC) is a type of breast cancer known for its poor clinical outcome. The lack of estrogen, progesterone, and human epidermal growth factor receptor in the tumors of TNBC leads to fewer treatment options in clinics. The incidence of TNBC is higher in African American (AA) women compared with European American (EA) women with worse clinical outcomes. The significant factors responsible for the racial disparity in TNBC are socioeconomic lifestyle and tumor biology. The current study considered the open-source gene expression data of triple-negative breast cancer samples' racial information. We implemented a state-of-the-art classification Support Vector Machine (SVM) method with a recurrent feature elimination approach to the gene expression data to identify significant biomarkers deregulated in AA women and EA women. We also included Spearman's rho and Ward's linkage method in our feature selection workflow. Our proposed method generates 24 features/genes that can classify the AA and EA samples 98% accurately. We also performed the Kaplan-Meier analysis and log-rank test on the 24 features/genes. We only discussed the correlation between deregulated expression and cancer progression with a poor survival rate of 2 genes, KLK10 and LRRC37A2, out of 24 genes. We believe that further improvement of our method with a higher number of RNA-seq gene expression data will more accurately provide insight into racial disparity in TNBC.