PancanQTL: systematic identification of cis-eQTLs and trans-eQTLs in 33 cancer types

Integrating single-cell with transcriptome-proteome Mendelian randomization reveals colorectal cancer targets

BACKGROUND

Colorectal carcinogenesis involves dynamic interactions between genetic susceptibility and cellular heterogeneity, yet current studies rarely disentangle causal genes from passive associations. While GWAS have mapped numerous risk loci, only a minority colocalize with eQTL/pQTL. A multi-omics framework combining single-cell transcriptomics, transcriptomics, proteomics, and MR is urgently needed to resolve cell-type-specific drivers of colorectal cancer pathogenesis.

METHODS

We integrated GWAS data, eQTL data, pQTL data, and single-cell RNA sequencing differential gene expression profiles from public databases. Subsequent batch Two-sample Mendelian randomization and further SMR analysis aimed to identify key genes in the pathogenesis of colorectal cancer.

RESULTS

Cluster analysis identified 4909 DEGs across various cell types. We discovered that 428 DEGs had a causal association with colorectal cancer through eQTL, of which 38 genes met the FDR statistical standards, and four of these genes (CTSF, PCSK7, LYZ, LMAN2L) also had causal associations through pQTL. SMR analysis confirmed the reliability of PCSK7 as a disease target.

CONCLUSION

By integrating single-cell data, transcriptomic data, proteomic data and GWAS data for MR analysis, we identified CTSF, PCSK7, LYZ, LMAN2L as potential targets for colorectal cancer.

Multicancer analyses of short tandem repeat variations reveal shared gene regulatory mechanisms

Short tandem repeats (STRs) have been reported to influence gene expression across various human tissues. While STR variations are enriched in colorectal, stomach, and endometrial cancers, particularly in microsatellite instable tumors, their functional effects and regulatory mechanisms on gene expression remain poorly understood across these cancer types. Here, we leverage whole-exome sequencing and gene expression data to identify STRs for which repeat lengths are associated with the expression of nearby genes (eSTRs) in colorectal, stomach, and endometrial tumors. While most eSTRs are cancer-specific, shared eSTRs across multiple cancers exhibit consistent effects on gene expression. Notably, coding-region eSTRs identified in all three cancer types show positive correlations with nearby gene expression. We further validate the functional effects of eSTRs by demonstrating associations between somatic eSTR mutations and gene expression changes during the transition from normal to tumor tissues, suggesting their potential roles in tumorigenesis. Combined with DNA methylation data, we perform the first quantitative analysis of the interplay between STR variations and DNA methylation in tumors. We identify eSTRs where repeat lengths are associated with methylation levels of nearby CpG sites (meSTRs) and show that >70% of eSTRs are significantly linked to local DNA methylation. Importantly, the effects of meSTRs on DNA methylation remain consistent across cancer types. Overall, our findings enhance the understanding of how functional STR variations influence gene expression and DNA methylation. Our study highlights shared regulatory mechanisms of STRs across multiple cancers, offering a foundation for future research into their broader implications in tumor biology.

Genetic Regulation of Alternative Polyadenylation Provides Novel Insights into Molecular Mechanisms Underlying Non-small Cell Lung Cancer

Emerging evidence emphasizes the critical role of alternative polyadenylation (APA) in posttranscriptional regulation of genes, and APA-associated genetic variants (apaQTLs) show particular relevance for multiple disease. However, genetic regulation of APA and its role in non-small cell lung cancer (NSCLC) risk have not been thoroughly studied. Here, by leveraging genotype and APA data from The Cancer Genome Atlas, the association between genetic variation and APA is determined in NSCLC samples. The identified apaQTLs are distinct from eQTLs and are preferentially enriched in functionally relevant characteristics, including poly(A) motifs, APA-associated RBP binding sites, functional elements, and known NSCLC risk loci. Moreover, genes associated with apaQTLs are broadly involved in cancer-related biological process. Of note, integration of apaQTL variants with traditional GWAS-derived PRS is proved as a potential screening tool for NSCLC. By integrating large-scale population and biological experiments, a functional apaQTL variant rs9606 in LYRM4 is identified. Mechanistically, rs9606 induces aberrant APA process of LYRM4 via allele-specific interacting with NUDT21, which lead to increased expression of oncogene LYRM4 and thus contribute to NSCLC risk. This study demonstrates the distinct contribution of APA-associated genetic variants in NSCLC risk, providing critical clues and potential targets for NSCLC etiology and clinical intervention.

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.

Deciphering genetic regulation at single-cell resolution in gastric cancer

Understanding cell-type-specific genetic regulation in gastric cancer is essential for uncovering disease susceptibility. By performing single-cell eQTL mapping in gastric tissues, Bian et al.1 identified previously uncharacterized regulatory genetic mechanisms, risk genes, and co-localization signals associated with gastric cancer susceptibility, providing insights into its pathogenesis and potential therapeutic approaches.

Enhancer RNA Transcriptome-Wide Association Study Reveals a Distinctive Class of Pan-Cancer Susceptibility eRNAs

Many cancer risk variants are located within enhancer regions and lack sufficient molecular interpretation. Here, we constructed the first comprehensive atlas of enhancer RNA (eRNA)-mediated genetic effects from 28 033 RNA sequencing samples across 11 606 individuals, identifying 21 073 eRNA quantitative trait loci (eRNA-QTLs) significantly associated with eRNA expression. Mechanistically, eRNA-QTLs frequently altered binding motifs of transcription factors. In addition, 28.48% of cancer risk variants are strongly colocalized with eRNA-QTLs. A pan-cancer eRNA-based transcriptome-wide association study is conducted across 23 major cancer types, identifying 626 significant cancer susceptibility eRNAs predicted to modulate cancer risk via eRNA, from which 54.90% of the eRNA target genes are overlooked by traditional gene expression studies, and most are essential for cancer cell proliferation. As proof of principle validation, the enhancer functionality of two newly identified susceptibility eRNAs, CCND1e and SNAPC1e, is confirmed through CRISPR inhibition and shRNA-mediated knockdown, resulting in a marked decrease in the expression of their respective target genes, consequently suppressing the proliferation of prostate cancer cells. The study underscores the essential role of eRNA in unveiling new cancer susceptibility genes and establishes a strong framework for enhancing our understanding of human cancer etiology.

Lessons learned from a candidate gene study investigating aromatase inhibitor treatment outcome in breast cancer

The role of germline genetics in adjuvant aromatase inhibitor (AI) treatment efficacy in ER-positive breast cancer is poorly understood. We employed a two-stage candidate gene approach to examine associations between survival endpoints and common germline variants in 753 endocrine resistance-related genes. For a discovery cohort, we screened the Breast Cancer Association Consortium database (n ≥ 90,000 cases) and retrieved 2789 AI-treated patients. Cox model-based analysis revealed 125 variants associated with overall, distant relapse-free, and relapse-free survival (p-value ≤ 1E-04). In validation analysis using five independent cohorts (n = 8857), none of the six selected candidates representing major linkage blocks at CELA2B/CASP9, NR1I2/GSK3B, LRP1B, and MIR143HG (CARMN) were validated. We discuss potential reasons for the failed validation and replication of published findings, including study/treatment heterogeneity and other limitations inherent to genomic treatment outcome studies. For the future, we envision prospective longitudinal studies with sufficiently long follow-up and endpoints that reflect the dynamic nature of endocrine resistance.

Overcoming collaboration barriers in quantitative trait loci analysis

In this issue of Cell Genomics, Choi et al.1 report a novel approach, privateQTL, which leverages secure multiparty computation (MPC) to enable federated expression quantitative trait loci (eQTL) mapping across institutions without compromising data privacy. Zhang et al. preview their approach and discuss its application in future genetic analyses.

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.

sc2GWAS: a comprehensive platform linking single cell and GWAS traits of human

Identifying cell populations associated with risk variants is essential for uncovering cell-specific mechanisms that drive disease development and progression. Integrating genome-wide association studies (GWAS) with single-cell RNA sequencing (scRNA-seq) has become an effective strategy for detecting trait-cell relationships. The accumulation of trait-related single cell data has led to an urgent need for its comprehensively processing. To address this, we developed sc2GWAS (https://bio.liclab.net/sc2GWAS/), which aims to document large-scale GWAS trait-cell regulatory pairs at single-cell resolution and provide comprehensive annotations and enrichment analyses for these related pairs. The current version of sc2GWAS curates a total of 15 078 310 candidate trait-cell pairs from > 6 300 000 individual cells, offering a valuable resource for exploring complex regulatory relationships between traits and cells. We applied strict quality control measures on both scRNA-seq data and GWAS data, ensuring the reliability and accuracy of the datasets for the identification of trait-relevant cells and genes. In addition, sc2GWAS provides ranked lists of trait-relevant genes and extensive (epi) genetic annotations, making it a valuable resource for downstream analyses. We demonstrate the utility of the platform by investigating Alzheimer's disease, where we identified significant associations between the disease and microglial cells, with the APOE gene emerging as particularly significant. This platform facilitates detailed research into complex trait-cell and trait-gene interactions, we anticipate that sc2GWAS will become a comprehensive and valuable platform for exploring GWAS trait-cell regulatory mechanisms.

miRNASNP-v4: a comprehensive database for miRNA-related SNPs across 17 species

Single nucleotide polymorphisms (SNPs) within microRNAs (miRNAs) and their target binding sites can influence miRNA biogenesis and target regulation, thereby participating in a variety of diseases and biological processes. Current miRNA-related SNP databases are often species-limited or based on outdated data. Therefore, we updated our miRNASNP database to version 4 by updating data, expanding the species from Homo sapiens to 17 species, and introducing several new features. In miRNASNP-v4, 82 580 SNPs in miRNAs and 24 836 179 SNPs in 3'UTRs of genes across 17 species were identified and their potential effects on miRNA secondary structure and target binding were characterized. In addition, compared to the last release, miRNASNP-v4 includes the following improvements: (i) gene enrichment analysis for gained or lost miRNA target genes; (ii) identification of miRNA-related SNPs associated with drug response and immune infiltration in human cancers; (iii) inclusion of experimentally supported immune-related miRNAs and (iv) online prediction tools for 17 animal species. With the extensive data and user-friendly web interface, miRNASNP-v4 will serve as an invaluable resource for functional studies of SNPs and miRNAs in multiple species. The database is freely accessible at http://gong_lab.hzau.edu.cn/miRNASNP/.

Comprehensive characterization of the transcriptional landscape in Alzheimer's disease (AD) brains

Alzheimer's disease (AD) is the leading dementia among the elderly with complex origins. Despite extensive investigation into the AD-associated protein-coding genes, the involvement of noncoding RNAs (ncRNAs) and posttranscriptional modification (PTM) in AD pathogenesis remains unclear. Here, we comprehensively characterized the landscape of ncRNAs and PTM events in 1460 samples across six brain regions sourced from the Mount Sinai/JJ Peters VA Medical Center Brain Bank Study and Mayo cohorts, encompassing 33,321 long ncRNAs, 92,897 enhancer RNAs, 53,763 alternative polyadenylation events, and 900,221 A-to-I RNA editing events. We additionally identified 25,351 aberrantly expressed ncRNAs and altered PTM events associated with AD traits and further identified the corresponding protein-coding genes to construct regulatory networks. Furthermore, we developed a user-friendly data portal, ADatlas, facilitating users in exploring our results. Our study aims to establish a comprehensive data platform for ncRNAs and PTMs in AD to advance related research.

The three YTHDF paralogs and VIRMA are strong cross-histotype tumor driver candidates among m6A core genes

N6-Methyladenosine (m6A) is the most abundant internal modification in mRNAs. Despite accumulating evidence for the profound impact of m6A on cancer biology, there are conflicting reports that alterations in genes encoding the m6A machinery proteins can either promote or suppress cancer, even in the same tumor type. Using data from The Cancer Genome Atlas, we performed a pan-cancer investigation of 15 m6A core factors in nearly 10000 samples from 31 tumor types to reveal underlying cross-tumor patterns. Altered expression, largely driven by copy number variations at the chromosome arm level, results in the most common mode of dysregulation of these factors. YTHDF1, YTHDF2, YTHDF3 and VIRMA are the most frequently altered factors and the only ones to be uniquely altered when tumors are grouped according to the expression pattern of the m6A factors. These genes are also the only ones with coherent, pan-cancer predictive power for progression-free survival. On the contrary, METTL3, the most intensively studied m6A factor as a cancer target, shows much lower levels of alteration and no predictive power for patient survival. Therefore, we propose the non-enzymatic YTHDF and VIRMA genes as preferred subjects to dissect the role of m6A in cancer and as priority cancer targets.

Machine learning model identifies genetic predictors of cisplatin-induced ototoxicity in CERS6 and TLR4

BACKGROUND

Cisplatin-induced ototoxicity remains a significant concern in pediatric cancer treatment due to its permanent impact on quality of life. Previously, genetic association analyses have been performed to detect genetic variants associated with this adverse reaction.

METHODS

In this study, a combination of interpretable neural networks and Generative Adversarial Networks (GANs) was employed to identify genetic markers associated with cisplatin-induced ototoxicity. The applied method, BRI-Net, incorporates biological domain knowledge to define the network structure and employs adversarial training to learn an unbiased representation of the data, which is robust to known confounders. Leveraging genomic data from a cohort of 362 cisplatin-treated pediatric cancer patients recruited by the CPNDS (Canadian Pharmacogenomics Network for Drug Safety), this model revealed two statistically significant single nucleotide polymorphisms to be associated with cisplatin-induced ototoxicity.

RESULTS

Two markers within the CERS6 (rs13022792, p-value: 3 × 10-4) and TLR4 (rs10759932, p-value: 7 × 10-4) genes were associated with this cisplatin-induced adverse reaction. CERS6, a ceramide synthase, contributes to elevated ceramide levels, a known initiator of apoptotic signals in mouse models of inner ear hair cells. TLR4, a pattern-recognition protein, initiates inflammation in response to cisplatin, and reduced TLR4 expression has been shown in murine hair cells to confer protection from ototoxicity.

CONCLUSION

Overall, these findings provide a foundation for understanding the genetic landscape of cisplatin-induced ototoxicity, with implications for improving patient care and treatment outcomes.

VCAT: an integrated variant function annotation tools

The development of sequencing technology has promoted discovery of variants in the human genome. Identifying functions of these variants is important for us to link genotype to phenotype, and to diagnose diseases. However, it usually requires researchers to visit multiple databases. Here, we presented a one-stop webserver for variant function annotation tools (VCAT, https://biomed.nscc-gz.cn/zhaolab/VCAT/ ) that is the first one connecting variant to functions via the epigenome, protein, drug and RNA. VCAT is also the first one to make all annotations visualized in interactive charts or molecular structures. VCAT allows users to upload data in VCF format, and download results via a URL. Moreover, VCAT has annotated a huge number (1,262,041,068) of variants collected from dbSNP, 1000 Genomes projects, gnomAD, ICGC, TCGA, and HPRC Pangenome project. For these variants, users are able to searcher their functions, related diseases and drugs from VCAT. In summary, VCAT provides a one-stop webserver to explore the potential functions of human genomic variants including their relationship with diseases and drugs.

Towards functional maps of non-coding variants in cancer

Large scale cancer genomic studies in patients have unveiled millions of non-coding variants. While a handful have been shown to drive cancer development, the vast majority have unknown function. This review describes the challenges of functionally annotating non-coding cancer variants and understanding how they contribute to cancer. We summarize recently developed high-throughput technologies to address these challenges. Finally, we outline future prospects for non-coding cancer genetics to help catalyze personalized cancer therapy.

Integrated systematic functional screen and fine-mapping decipher the role and genetic regulation of RPS19 in colorectal cancer development

Despite genome-wide association studies (GWAS) have identified more than 200 risk loci associated with colorectal cancer (CRC), the causal genes or risk variants within these loci and their biological functions remain not fully revealed. Recently, the genomic locus 19q13.2, with the lead SNP rs1800469 was identified as a crucial CRC risk locus in Asian populations. However, the functional mechanism of this region has not been fully elucidated. Here we employed an RNA interfering-based on-chip approach to screen for the genes essential for cell proliferation in the CRC risk locus 19q13.2. Notably, we found that RPS19 exhibited the most significant effect among the identified genes and acted as a critical oncogene facilitating CRC cell proliferation. Subsequently, combining integrative fine-mapping analysis and a large-scale population study consisting of 6027 cases and 6099 controls, we prioritized rs1025497 as a potential causal candidate for CRC risk, demonstrating that rs1025497[A] allele significantly reduced the risk of CRC (OR 0.70, 95% confidence interval = 0.56-0.83, P = 1.12 × 10-6), which was further validated in UK Biobank cohort comprising 5,313 cases and 21,252 controls. Mechanistically, we experimentally elucidated that variant rs1025497 might acted as an allele-specific silencer, inhibiting the expression level of oncogene RPS19 mediated by the transcription suppressive factor HBP1. Taken together, our sturdy unveils the significant role of RPS19 during CRC pathogenesis and delineates its distal regulatory mechanism mediated by rs1025497, advancing our understanding of the etiology of CRC and provided new insights into the personalized medicine of human cancer.

Systematic prioritization of functional variants and effector genes underlying colorectal cancer risk

Genome-wide association studies of colorectal cancer (CRC) have identified 170 autosomal risk loci. However, for most of these, the functional variants and their target genes are unknown. Here, we perform statistical fine-mapping incorporating tissue-specific epigenetic annotations and massively parallel reporter assays to systematically prioritize functional variants for each CRC risk locus. We identify plausible causal variants for the 170 risk loci, with a single variant for 40. We link these variants to 208 target genes by analyzing colon-specific quantitative trait loci and implementing the activity-by-contact model, which integrates epigenomic features and Micro-C data, to predict enhancer-gene connections. By deciphering CRC risk loci, we identify direct links between risk variants and target genes, providing further insight into the molecular basis of CRC susceptibility and highlighting potential pharmaceutical targets for prevention and treatment.

Bioinformatics tools and resources for cancer and application

ABSTRACT

Tumor bioinformatics plays an important role in cancer research and precision medicine. The primary focus of traditional cancer research has been molecular and clinical studies of a number of fundamental pathways and genes. In recent years, driven by breakthroughs in high-throughput technologies, large-scale cancer omics data have accumulated rapidly. How to effectively utilize and share these data is particularly important. To address this crucial task, many computational tools and databases have been developed over the past few years. To help researchers quickly learn and understand the functions of these tools, in this review, we summarize publicly available bioinformatics tools and resources for pan-cancer multi-omics analysis, regulatory analysis of tumorigenesis, tumor treatment and prognosis, immune infiltration analysis, immune repertoire analysis, cancer driver gene and driver mutation analysis, and cancer single-cell analysis, which may further help researchers find more suitable tools for their research.

eQTL-Detect: nextflow-based pipeline for eQTL detection in modular format with sharable and parallelizable scripts

Bioinformatic pipelines are becoming increasingly complex with the ever-accumulating amount of Next-generation sequencing (NGS) data. Their orchestration is difficult with a simple Bash script, but bioinformatics workflow managers such as Nextflow provide a framework to overcome respective problems. This study used Nextflow to develop a bioinformatic pipeline for detecting expression quantitative trait loci (eQTL) using a DSL2 Nextflow modular syntax, to enable sharing the huge demand for computing power as well as data access limitation across different partners often associated with eQTL studies. Based on the results from a test run with pilot data by measuring the required runtime and computational resources, the new pipeline should be suitable for eQTL studies in large scale analyses.

Multiomic analysis of uterine leiomyomas in self-described Black and White women: molecular insights into health disparities

BACKGROUND

Black women are at an increased risk of developing uterine leiomyomas and experiencing worse disease prognosis than White women. Epidemiologic and molecular factors have been identified as underlying these disparities, but there remains a paucity of deep, multiomic analysis investigating molecular differences in uterine leiomyomas from Black and White patients.

OBJECTIVE

To identify molecular alterations within uterine leiomyoma tissues correlating with patient race by multiomic analyses of uterine leiomyomas collected from cohorts of Black and White women.

STUDY DESIGN

We performed multiomic analysis of uterine leiomyomas from Black (42) and White (47) women undergoing hysterectomy for symptomatic uterine leiomyomata. In addition, our analysis included the application of orthogonal methods to evaluate fibroid biomechanical properties, such as second harmonic generation microscopy, uniaxial compression testing, and shear-wave ultrasonography analyses.

RESULTS

We found a greater proportion of MED12 mutant uterine leiomyomas from Black women (>35% increase; Mann-Whitney U, P<.001). MED12 mutant tumors exhibited an elevated abundance of extracellular matrix proteins, including several collagen isoforms, involved in the regulation of the core matrisome. Histologic analysis of tissue fibrosis using trichrome staining and secondary harmonic generation microscopy confirmed that MED12 mutant tumors are more fibrotic than MED12 wild-type tumors. Using shear-wave ultrasonography in a prospectively collected cohort, Black patients had fibroids that were firmer than White patients, even when similar in size. In addition, these analyses uncovered ancestry-linked expression quantitative trait loci with altered allele frequencies in African and European populations correlating with differential abundance of several proteins in uterine leiomyomas independently of MED12 mutation status, including tetratricopeptide repeat protein 38.

CONCLUSION

Our study shows that Black women have a higher prevalence of uterine leiomyomas harboring mutations in MED12 and that this mutational status correlates with increased tissue fibrosis compared with wild-type uterine leiomyomas. Our study provides insights into molecular alterations correlating with racial disparities in uterine leiomyomas and improves our understanding of the molecular etiology underlying uterine leiomyoma development within these populations.

Uncovering methylation-dependent genetic effects on regulatory element function in diverse genomes

A major goal in evolutionary biology and biomedicine is to understand the complex interactions between genetic variants, the epigenome, and gene expression. However, the causal relationships between these factors remain poorly understood. mSTARR-seq, a methylation-sensitive massively parallel reporter assay, is capable of identifying methylation-dependent regulatory activity at many thousands of genomic regions simultaneously, and allows for the testing of causal relationships between DNA methylation and gene expression on a region-by-region basis. Here, we developed a multiplexed mSTARR-seq protocol to assay naturally occurring human genetic variation from 25 individuals sampled from 10 localities in Europe and Africa. We identified 6,957 regulatory elements in either the unmethylated or methylated state, and this set was enriched for enhancer and promoter annotations, as expected. The expression of 58% of these regulatory elements was modulated by methylation, which was generally associated with decreased RNA expression. Within our set of regulatory elements, we used allele-specific expression analyses to identify 8,020 sites with genetic effects on gene regulation; further, we found that 42.3% of these genetic effects varied between methylated and unmethylated states. Sites exhibiting methylation-dependent genetic effects were enriched for GWAS and EWAS annotations, implicating them in human disease. Compared to datasets that assay DNA from a single European individual, our multiplexed assay uncovers dramatically more genetic effects and methylation-dependent genetic effects, highlighting the importance of including diverse individuals in assays which aim to understand gene regulatory processes.

Mechanistic analysis of enhancer sequences in the estrogen receptor transcriptional program

Estrogen Receptor α (ERα) is a major lineage determining transcription factor (TF) in mammary gland development. Dysregulation of ERα-mediated transcriptional program results in cancer. Transcriptomic and epigenomic profiling of breast cancer cell lines has revealed large numbers of enhancers involved in this regulatory program, but how these enhancers encode function in their sequence remains poorly understood. A subset of ERα-bound enhancers are transcribed into short bidirectional RNA (enhancer RNA or eRNA), and this property is believed to be a reliable marker of active enhancers. We therefore analyze thousands of ERα-bound enhancers and build quantitative, mechanism-aware models to discriminate eRNAs from non-transcribing enhancers based on their sequence. Our thermodynamics-based models provide insights into the roles of specific TFs in ERα-mediated transcriptional program, many of which are supported by the literature. We use in silico perturbations to predict TF-enhancer regulatory relationships and integrate these findings with experimentally determined enhancer-promoter interactions to construct a gene regulatory network. We also demonstrate that the model can prioritize breast cancer-related sequence variants while providing mechanistic explanations for their function. Finally, we experimentally validate the model-proposed mechanisms underlying three such variants.

Decoding the interplay between genetic and non-genetic drivers of metastasis

Metastasis is a multistep process by which cancer cells break away from their original location and spread to distant organs, and is responsible for the vast majority of cancer-related deaths. Preventing early metastatic dissemination would revolutionize the ability to fight cancer. Unfortunately, the relatively poor understanding of the molecular underpinnings of metastasis has hampered the development of effective anti-metastatic drugs. Although it is now accepted that disseminating tumour cells need to acquire multiple competencies to face the many obstacles they encounter before reaching their metastatic site(s), whether these competencies are acquired through an accumulation of metastasis-specific genetic alterations and/or non-genetic events is often debated. Here we review a growing body of literature highlighting the importance of both genetic and non-genetic reprogramming events during the metastatic cascade, and discuss how genetic and non-genetic processes act in concert to confer metastatic competencies. We also describe how recent technological advances, and in particular the advent of single-cell multi-omics and barcoding approaches, will help to better elucidate the cross-talk between genetic and non-genetic mechanisms of metastasis and ultimately inform innovative paths for the early detection and interception of this lethal process.

Causal effect of gut microbiota on pancreatic cancer: A Mendelian randomization and colocalization study

The causal relationship between gut microbiota (GM) and pancreatic cancer (PC) remains unclear. This study aimed to investigate the potential genes underlying this mechanism. GM Genome-wide association study (GWAS) summary data were from the MiBioGen consortium. PC GWAS data were from the National Human Genome Research Institute-European Bioinformatics Institute (NHGRI-EBI) GWAS Catalogue. To detect the causal relationship between GM and PC, we implemented three complementary Mendelian randomization (MR) methods: Inverse Variance Weighting (IVW), MR-Egger and Weighted Median, followed by sensitivity analyses. Furthermore, we integrated GM GWAS data with blood cis-expression quantitative trait loci (eQTLs) and blood cis-DNA methylation QTL (mQTLs) using Summary data-based Mendelian Randomization (SMR) methods. This integration aimed to prioritize potential GM-affecting genes through SMR analysis of two molecular traits. PC cis-eQTLs and cis-mQTLs were summarized from The Cancer Genome Atlas (TCGA) data. Through colocalization analysis of GM cis-QTLs and PC cis-QTLs data, we identified common genes that influence both GM and PC. Our study found a causal association between GM and PC, including four protective and five risk-associated GM [Inverse Variance Weighted (IVW), p < 0.05]. No significant heterogeneity of instrumental variables (IVs) or horizontal pleiotropy was found. The gene SVBP was identified as a GM-affecting gene using SMR analysis of two molecular traits (FDR<0.05, P_HEIDI>0.05). Additionally, two genes, MCM6 and RPS26, were implicated in the interaction between GM and PC based on colocalization analysis (PPH4>0.5). In summary, this study provides evidence for future research aimed at developing suitable therapeutic interventions and disease prevention.

New role of fat-free mass in cancer risk linked with genetic predisposition

Cancer risk is associated with the widely debated measure body mass index (BMI). Fat mass and fat-free mass measurements from bioelectrical impedance may further clarify this association. The UK Biobank is a rare resource in which bioelectrical impedance and BMI data was collected on ~ 500,000 individuals. Using this dataset, a comprehensive analysis using regression, principal component and genome-wide genetic association, provided multiple levels of evidence that increasing whole body fat (WBFM) and fat-free mass (WBFFM) are both associated with increased post-menopausal breast cancer risk, and colorectal cancer risk in men. WBFM was inversely associated with prostate cancer. We also identified rs615029[T] and rs1485995[G] as associated in independent analyses with both PMBC (p = 1.56E-17 and 1.78E-11) and WBFFM (p = 2.88E-08 and 8.24E-12), highlighting splice variants of the intriguing long non-coding RNA CUPID1 (LINC01488) as a potential link between PMBC risk and fat-free mass.

Incorporating genetic similarity of auxiliary samples into eGene identification under the transfer learning framework

BACKGROUND

The term eGene has been applied to define a gene whose expression level is affected by at least one independent expression quantitative trait locus (eQTL). It is both theoretically and empirically important to identify eQTLs and eGenes in genomic studies. However, standard eGene detection methods generally focus on individual cis-variants and cannot efficiently leverage useful knowledge acquired from auxiliary samples into target studies.

METHODS

We propose a multilocus-based eGene identification method called TLegene by integrating shared genetic similarity information available from auxiliary studies under the statistical framework of transfer learning. We apply TLegene to eGene identification in ten TCGA cancers which have an explicit relevant tissue in the GTEx project, and learn genetic effect of variant in TCGA from GTEx. We also adopt TLegene to the Geuvadis project to evaluate its usefulness in non-cancer studies.

RESULTS

We observed substantial genetic effect correlation of cis-variants between TCGA and GTEx for a larger number of genes. Furthermore, consistent with the results of our simulations, we found that TLegene was more powerful than existing methods and thus identified 169 distinct candidate eGenes, which was much larger than the approach that did not consider knowledge transfer across target and auxiliary studies. Previous studies and functional enrichment analyses provided empirical evidence supporting the associations of discovered eGenes, and it also showed evidence of allelic heterogeneity of gene expression. Furthermore, TLegene identified more eGenes in Geuvadis and revealed that these eGenes were mainly enriched in cells EBV transformed lymphocytes tissue.

CONCLUSION

Overall, TLegene represents a flexible and powerful statistical method for eGene identification through transfer learning of genetic similarity shared across auxiliary and target studies.

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

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

A Bayesian fine-mapping model using a continuous global-local shrinkage prior with applications in prostate cancer analysis

The aim of fine mapping is to identify genetic variants causally contributing to complex traits or diseases. Existing fine-mapping methods employ Bayesian discrete mixture priors and depend on a pre-specified maximum number of causal variants, which may lead to sub-optimal solutions. In this work, we propose a Bayesian fine-mapping method called h2-D2, utilizing a continuous global-local shrinkage prior. We also present an approach to define credible sets of causal variants in continuous prior settings. Simulation studies demonstrate that h2-D2 outperforms current state-of-the-art fine-mapping methods such as SuSiE and FINEMAP in accurately identifying causal variants and estimating their effect sizes. We further applied h2-D2 to prostate cancer analysis and discovered some previously unknown causal variants. In addition, we inferred 369 target genes associated with the detected causal variants and several pathways that were significantly over-represented by these genes, shedding light on their potential roles in prostate cancer development and progression.

NCAD v1.0: a database for non-coding variant annotation and interpretation

The application of whole genome sequencing is expanding in clinical diagnostics across various genetic disorders, and the significance of non-coding variants in penetrant diseases is increasingly being demonstrated. Therefore, it is urgent to improve the diagnostic yield by exploring the pathogenic mechanisms of variants in non-coding regions. However, the interpretation of non-coding variants remains a significant challenge, due to the complex functional regulatory mechanisms of non-coding regions and the current limitations of available databases and tools. Hence, we develop the non-coding variant annotation database (NCAD, http://www.ncawdb.net/), encompassing comprehensive insights into 665,679,194 variants, regulatory elements, and element interaction details. Integrating data from 96 sources, spanning both GRCh37 and GRCh38 versions, NCAD v1.0 provides vital information to support the genetic diagnosis of non-coding variants, including allele frequencies of 12 diverse populations, with a particular focus on the population frequency information for 230,235,698 variants in 20,964 Chinese individuals. Moreover, it offers prediction scores for variant functionality, five categories of regulatory elements, and four types of non-coding RNAs. With its rich data and comprehensive coverage, NCAD serves as a valuable platform, empowering researchers and clinicians with profound insights into non-coding regulatory mechanisms while facilitating the interpretation of non-coding variants.

A Multiethnic Germline-Somatic Association Database Deciphers Multilayered and Interconnected Genetic Mutations in Cancer

Inherited germline and acquired somatic alterations can both promote human tumor development. Elucidating the cooperation between somatic and germline genetic alterations that drive tumorigenesis could help inform precision cancer prevention and treatment strategies. Here, leveraging genomic genotyping and sequencing data from 9,029 patients with cancer with European, East Asian, and African ancestry, we performed a pan-cancer analysis to evaluate the associations between germline SNPs and somatic alterations, including single-nucleotide variant and small insertion/deletion mutations, copy-number variation, tumor mutational burden, and mutational signatures. Genome-wide significant germline-somatic pairs were abundant, and most of the associations were observed in one cancer type and one ancestry group. A user-friendly interactive Multiethnic Germline-Somatic Association (MGSA) database (http://wanglab-hkust.cn:3838/MGSA/) was developed, which can be used to query, browse, and download the results of the association analyses. Moreover, the MGSA database offers additional survival analysis and functional annotation. Together, this work provides a resource for uncovering the clinical and biological roles of associations between germline variants and somatic alterations in human cancer.

SIGNIFICANCE

Comprehensive analysis of connections between germline variants and somatic events in cancer offers a resource for investigating the functional significance of genetic mutations and exploring genetic factors contributing to racial disparities.

KnockTF 2.0: a comprehensive gene expression profile database with knockdown/knockout of transcription (co-)factors in multiple species

Transcription factors (TFs), transcription co-factors (TcoFs) and their target genes perform essential functions in diseases and biological processes. KnockTF 2.0 (http://www.licpathway.net/KnockTF/index.html) aims to provide comprehensive gene expression profile datasets before/after T(co)F knockdown/knockout across multiple tissue/cell types of different species. Compared with KnockTF 1.0, KnockTF 2.0 has the following improvements: (i) Newly added T(co)F knockdown/knockout datasets in mice, Arabidopsis thaliana and Zea mays and also an expanded scale of datasets in humans. Currently, KnockTF 2.0 stores 1468 manually curated RNA-seq and microarray datasets associated with 612 TFs and 172 TcoFs disrupted by different knockdown/knockout techniques, which are 2.5 times larger than those of KnockTF 1.0. (ii) Newly added (epi)genetic annotations for T(co)F target genes in humans and mice, such as super-enhancers, common SNPs, methylation sites and chromatin interactions. (iii) Newly embedded and updated search and analysis tools, including T(co)F Enrichment (GSEA), Pathway Downstream Analysis and Search by Target Gene (BLAST). KnockTF 2.0 is a comprehensive update of KnockTF 1.0, which provides more T(co)F knockdown/knockout datasets and (epi)genetic annotations across multiple species than KnockTF 1.0. KnockTF 2.0 facilitates not only the identification of functional T(co)Fs and target genes but also the investigation of their roles in the physiological and pathological processes.

eRNAbase: a comprehensive database for decoding the regulatory eRNAs in human and mouse

Enhancer RNAs (eRNAs) transcribed from distal active enhancers serve as key regulators in gene transcriptional regulation. The accumulation of eRNAs from multiple sequencing assays has led to an urgent need to comprehensively collect and process these data to illustrate the regulatory landscape of eRNAs. To address this need, we developed the eRNAbase (http://bio.liclab.net/eRNAbase/index.php) to store the massive available resources of human and mouse eRNAs and provide comprehensive annotation and analyses for eRNAs. The current version of eRNAbase cataloged 10 399 928 eRNAs from 1012 samples, including 858 human samples and 154 mouse samples. These eRNAs were first identified and uniformly processed from 14 eRNA-related experiment types manually collected from GEO/SRA and ENCODE. Importantly, the eRNAbase provides detailed and abundant (epi)genetic annotations in eRNA regions, such as super enhancers, enhancers, common single nucleotide polymorphisms, expression quantitative trait loci, transcription factor binding sites, CRISPR/Cas9 target sites, DNase I hypersensitivity sites, chromatin accessibility regions, methylation sites, chromatin interactions regions, topologically associating domains and RNA spatial interactions. Furthermore, the eRNAbase provides users with three novel analyses including eRNA-mediated pathway regulatory analysis, eRNA-based variation interpretation analysis and eRNA-mediated TF-target gene analysis. Hence, eRNAbase is a powerful platform to query, browse and visualize regulatory cues associated with eRNAs.

EndoQuad: a comprehensive genome-wide experimentally validated endogenous G-quadruplex database

G-quadruplexes (G4s) are non-canonical four-stranded structures and are emerging as novel genetic regulatory elements. However, a comprehensive genomic annotation of endogenous G4s (eG4s) and systematic characterization of their regulatory network are still lacking, posing major challenges for eG4 research. Here, we present EndoQuad (https://EndoQuad.chenzxlab.cn/) to address these pressing issues by integrating high-throughput experimental data. First, based on high-quality genome-wide eG4s mapping datasets (human: 1181; mouse: 24; chicken: 2) generated by G4 ChIP-seq/CUT&Tag, we generate a reference set of genome-wide eG4s. Our multi-omics analyses show that most eG4s are identified in one or a few cell types. The eG4s with higher occurrences across samples are more structurally stable, evolutionarily conserved, enriched in promoter regions, mark highly expressed genes and associate with complex regulatory programs, demonstrating higher confidence level for further experiments. Finally, we integrate millions of functional genomic variants and prioritize eG4s with regulatory functions in disease and cancer contexts. These efforts have culminated in the comprehensive and interactive database of experimentally validated DNA eG4s. As such, EndoQuad enables users to easily access, download and repurpose these data for their own research. EndoQuad will become a one-stop resource for eG4 research and lay the foundation for future functional studies.

PancanQTLv2.0: a comprehensive resource for expression quantitative trait loci across human cancers

Expression quantitative trait locus (eQTL) analysis is a powerful tool used to investigate genetic variations in complex diseases, including cancer. We previously developed a comprehensive database, PancanQTL, to characterize cancer eQTLs using The Cancer Genome Atlas (TCGA) dataset, and linked eQTLs with patient survival and GWAS risk variants. Here, we present an updated version, PancanQTLv2.0 (https://hanlaboratory.com/PancanQTLv2/), with advancements in fine-mapping causal variants for eQTLs, updating eQTLs overlapping with GWAS linkage disequilibrium regions and identifying eQTLs associated with drug response and immune infiltration. Through fine-mapping analysis, we identified 58 747 fine-mapped eQTLs credible sets, providing mechanic insights of gene regulation in cancer. We further integrated the latest GWAS Catalog and identified a total of 84 592 135 linkage associations between eQTLs and the existing GWAS loci, which represents a remarkable ∼50-fold increase compared to the previous version. Additionally, PancanQTLv2.0 uncovered 659516 associations between eQTLs and drug response and identified 146948 associations between eQTLs and immune cell abundance, providing potentially clinical utility of eQTLs in cancer therapy. PancanQTLv2.0 expanded the resources available for investigating gene expression regulation in human cancers, leading to advancements in cancer research and precision oncology.

Pan-cancer analysis of promoter activity quantitative trait loci

Altered promoter activity has been generally observed in diverse biological processes, including tumorigenesis. Accumulating evidence suggests that employing a quantitative trait locus mapping approach is effective in comprehending the genetic basis of promoter activity. By utilizing genotype data from The Cancer Genome Atlas and calculating corresponding promoter activity values using proActiv, we systematically evaluated the impact of genetic variants on promoter activity and identified >1.0 million promoter activity quantitative trait loci (paQTLs) as both cis- and trans-acting. Additionally, leveraging data from the genome-wide association study (GWAS) catalog, we discovered >1.3 million paQTLs that overlap with known GWAS linkage disequilibrium regions. Remarkably, ∼9324 paQTLs exhibited significant associations with patient prognosis. Moreover, investigating the impact of promoter activity on >1000 imputed antitumor therapy responses among pan-cancer patients revealed >43 000 million significant associations. Furthermore, ∼25 000 significant associations were identified between promoter activity and immune cell abundance. Finally, a user-friendly data portal, Pancan-paQTL (https://www.hbpding.com/PancanPaQTL/), was constructed for users to browse, search and download data of interest. Pancan-paQTL serves as a comprehensive multidimensional database, enabling functional and clinical investigations into genetic variants associated with promoter activity, drug responses and immune infiltration across multiple cancer types.

Characterizing prostate cancer risk through multi-ancestry genome-wide discovery of 187 novel risk variants

The transferability and clinical value of genetic risk scores (GRSs) across populations remain limited due to an imbalance in genetic studies across ancestrally diverse populations. Here we conducted a multi-ancestry genome-wide association study of 156,319 prostate cancer cases and 788,443 controls of European, African, Asian and Hispanic men, reflecting a 57% increase in the number of non-European cases over previous prostate cancer genome-wide association studies. We identified 187 novel risk variants for prostate cancer, increasing the total number of risk variants to 451. An externally replicated multi-ancestry GRS was associated with risk that ranged from 1.8 (per standard deviation) in African ancestry men to 2.2 in European ancestry men. The GRS was associated with a greater risk of aggressive versus non-aggressive disease in men of African ancestry (P = 0.03). Our study presents novel prostate cancer susceptibility loci and a GRS with effective risk stratification across ancestry groups.

Finding Needles in the Haystack: Strategies for Uncovering Noncoding Regulatory Variants

Despite accumulating evidence implicating noncoding variants in human diseases, unraveling their functionality remains a significant challenge. Systematic annotations of the regulatory landscape and the growth of sequence variant data sets have fueled the development of tools and methods to identify causal noncoding variants and evaluate their regulatory effects. Here, we review the latest advances in the field and discuss potential future research avenues to gain a more in-depth understanding of noncoding regulatory variants.

Genetic Variants That Impact Alternative Polyadenylation in Cancer Represent Candidate Causal Risk Loci

Alternative polyadenylation (APA) is emerging as a major mechanism of posttranscriptional regulation. APA can impact the development and progression of cancer, suggesting that the genetic determinants of APA might play an important role in regulating cancer risk. Here, we depicted a pan-cancer atlas of human APA quantitative trait loci (apaQTL), containing approximately 0.7 million apaQTLs across 32 cancer types. Systematic multiomics analyses indicated that cancer apaQTLs could contribute to APA regulation by altering poly(A) motifs, RNA-binding proteins (RBP), and chromatin regulatory elements and were preferentially enriched in genome-wide association studies (GWAS)-identified cancer susceptibility loci. Moreover, apaQTL-related genes (aGene) were broadly related to cancer signaling pathways, high mutational burden, immune infiltration, and drug response, implicating their potential as therapeutic targets. Furthermore, apaQTLs were mapped in Chinese colorectal cancer tumor tissues and then screened for functional apaQTLs associated with colorectal cancer risk in 17,789 cases and 19,951 controls using GWAS-ChIP data, with independent validation in a large-scale population consisting of 6,024 cases and 10,022 controls. A multi-ancestry-associated apaQTL variant rs1020670 with a C>G change in DNM1L was identified, and the G allele contributed to an increased risk of colorectal cancer. Mechanistically, the risk variant promoted aberrant APA and facilitated higher usage of DNM1L proximal poly(A) sites mediated by the RBP CSTF2T, which led to higher expression of DNM1L with a short 3'UTR. This stabilized DNM1L to upregulate its expression, provoking colorectal cancer cell proliferation. Collectively, these findings generate a resource for understanding APA regulation and the genetic basis of human cancers, providing insights into cancer etiology.

SIGNIFICANCE

Cancer risk is mediated by alternative polyadenylation quantitative trait loci, including the rs1020670-G variant that promotes alternative polyadenylation of DNM1L and increases colorectal cancer risk.

Leveraging molecular quantitative trait loci to comprehend complex diseases/traits from the omics perspective

Comprehending the molecular basis of quantitative genetic variation is a principal goal for complex diseases or traits. Molecular quantitative trait loci (molQTLs) have made it possible to investigate the effects of genetic variants hiding behind large-scale omics data. A deeper understanding of molQTL is urgently required in light of the multi-dimensionalization of omics data to more fully elucidate the pertinent biological mechanisms. Herein, we reviewed molQTLs with the corresponding resource from the omics perspective and further discussed the integrative strategy of GWAS-molQTL to infer their causal effects. Subsequently, we described the opportunities and challenges encountered by molQTL. The case studies showed that molQTL is essential for complex diseases and traits, whether single- or multi-omics QTLs. Overall, we highlighted the functional significance of genetic variants to employ the discovery of molQTL in complex diseases and traits.

AgingReG: a curated database of aging regulatory relationships in humans

Aging and cellular senescence are characterized by a progressive loss of physiological integrity, which could be triggered by aging factors such as physiological, pathological and external factors. Numerous studies have shown that gene regulatory events play crucial roles in aging, increasing the need for a comprehensive repository of regulatory relationships during aging. Here, we established a manually curated database of aging factors (AgingReG, https://bio.liclab.net/Aging-ReG/), focusing on the regulatory relationships during aging with experimental evidence in humans. By curating thousands of published literature, 2157 aging factor entries (1345 aging gene entries, 804 external factor entries and eight aging-related pathway entries) and related regulatory information were manually curated. The regulatory relationships were classified into four types according to their functions: (i) upregulation, which indicates that aging factors upregulate the expression of target genes during aging; (ii) downregulation, which indicates that aging factors downregulate the expression of target genes during aging; (iii) activation, which indicates that aging factors influence the activity of target genes during aging and (iv) inhibition, which indicates that aging factors inhibit the activation of target molecule activity, leading to declined or lost target activity. AgingReG involves 651 upregulating pairs, 632 downregulating pairs, 330 activation-regulating pairs and 34 inhibition-regulating pairs, covering 195 disease types and more than 800 kinds of cells and tissues from 1784 published literature studies. AgingReG provides a user-friendly interface to query, browse and visualize detailed information about the regulatory relationships during aging. We believe that AgingReG will serve as a valuable resource database in the field of aging research. Database URL: https://bio.liclab.net/Aging-ReG/.

Genome-wide enhancer-gene regulatory maps link causal variants to target genes underlying human cancer risk

Genome-wide association studies have identified numerous variants associated with human complex traits, most of which reside in the non-coding regions, but biological mechanisms remain unclear. However, assigning function to the non-coding elements is still challenging. Here we apply Activity-by-Contact (ABC) model to evaluate enhancer-gene regulation effect by integrating multi-omics data and identified 544,849 connections across 20 cancer types. ABC model outperforms previous approaches in linking regulatory variants to target genes. Furthermore, we identify over 30,000 enhancer-gene connections in colorectal cancer (CRC) tissues. By integrating large-scale population cohorts (23,813 cases and 29,973 controls) and multipronged functional assays, we demonstrate an ABC regulatory variant rs4810856 associated with CRC risk (Odds Ratio = 1.11, 95%CI = 1.05-1.16, P = 4.02 × 10-5) by acting as an allele-specific enhancer to distally facilitate PREX1, CSE1L and STAU1 expression, which synergistically activate p-AKT signaling. Our study provides comprehensive regulation maps and illuminates a single variant regulating multiple genes, providing insights into cancer etiology.

Cis-Cardio: A comprehensive analysis platform for cardiovascular-relavant cis-regulation in human and mouse

Cis-regulatory elements are important molecular switches in controlling gene expression and are regarded as determinant hubs in the transcriptional regulatory network. Collection and processing of large-scale cis-regulatory data are urgent to decipher the potential mechanisms of cardiovascular diseases from a cis-regulatory element aspect. Here, we developed a novel web server, Cis-Cardio, which aims to document a large number of available cardiovascular-related cis-regulatory data and to provide analysis for unveiling the comprehensive mechanisms at a cis-regulation level. The current version of Cis-Cardio catalogs a total of 45,382,361 genomic regions from 1,013 human and mouse epigenetic datasets, including ATAC-seq, DNase-seq, Histone ChIP-seq, TF/TcoF ChIP-seq, RNA polymerase ChIP-seq, and Cohesin ChIP-seq. Importantly, Cis-Cardio provides six analysis tools, including region overlap analysis, element upstream/downstream analysis, transcription regulator enrichment analysis, variant interpretation, and protein-protein interaction-based co-regulatory analysis. Additionally, Cis-Cardio provides detailed and abundant (epi-) genetic annotations in cis-regulatory regions, such as super-enhancers, enhancers, transcription factor binding sites (TFBSs), methylation sites, common SNPs, risk SNPs, expression quantitative trait loci (eQTLs), motifs, DNase I hypersensitive sites (DHSs), and 3D chromatin interactions. In summary, Cis-Cardio is a valuable resource for elucidating and analyzing regulatory cues of cardiovascular-specific cis-regulatory elements. The platform is freely available at http://www.licpathway.net/Cis-Cardio/index.html.

Integration of eQTL and GWAS analysis uncovers a genetic regulation of natural ionomic variation in Arabidopsis

KEY MESSAGE

This study provided important insights into the genetic architecture of variations in A. thaliana leaf ionome in a cell-type-specific manner. The functional interpretation of traits associated variants by expression quantitative trait loci (eQTL) analysis is usually performed in bulk tissue samples. While the regulation of gene expression is context-dependent, such as cell-type-specific manner. In this study, we estimated cell-type abundances from 728 bulk tissue samples using single-cell RNA-sequencing dataset, and performed cis-eQTL mapping to identify cell-type-interaction eQTL (cis-eQTLs(ci)) in A. thaliana. Also, we performed Genome-wide association studies (GWAS) analyses for 999 accessions to identify the genetic basis of variations in A. thaliana leaf ionome. As a result, a total of 5,664 unique eQTL genes and 15,038 unique cis-eQTLs(ci) were significant. The majority (62.83%) of cis-eQTLs(ci) were cell-type-specific eQTLs. Using colocalization, we uncovered one interested gene AT2G25590 in Phloem cell, encoding a kind of plant Tudor-like protein with possible chromatin-associated functions, which colocalized with the most significant cis-eQTL(ci) of a Mo-related locus (Chr2:10,908,806:A:C; P = 3.27 × 10-27). Furthermore, we prioritized eight target genes associated with AT2G25590, which were previously reported in regulating the concentration of Mo element in A. thaliana. This study revealed the genetic regulation of ionomic variations and provided a foundation for further studies on molecular mechanisms of genetic variants controlling the A. thaliana ionome.

GATA2 co-opts TGFβ1/SMAD4 oncogenic signaling and inherited variants at 6q22 to modulate prostate cancer progression

BACKGROUND

Aberrant somatic genomic alteration including copy number amplification is a hallmark of cancer genomes. We previously profiled genomic landscapes of prostate cancer (PCa), yet the underlying causal genes with prognostic potential has not been defined. It remains unclear how a somatic genomic event cooperates with inherited germline variants contribute to cancer predisposition and progression.

METHODS

We applied integrated genomic and clinical data, experimental models and bioinformatic analysis to identify GATA2 as a highly prevalent metastasis-associated genomic amplification in PCa. Biological roles of GATA2 in PCa metastasis was determined in vitro and in vivo. Global chromatin co-occupancy and co-regulation of GATA2 and SMAD4 was investigated by coimmunoprecipitation, ChIP-seq and RNA-seq assays. Tumor cellular assays, qRT-PCR, western blot, ChIP, luciferase assays and CRISPR-Cas9 editing methods were performed to mechanistically understand the cooperation of GATA2 with SMAD4 in promoting TGFβ1 and AR signaling and mediating inherited PCa risk and progression.

RESULTS

In this study, by integrated genomics and experimental analysis, we identified GATA2 as a prevalent metastasis-associated genomic amplification to transcriptionally augment its own expression in PCa. Functional experiments demonstrated that GATA2 physically interacted and cooperated with SMAD4 for genome-wide chromatin co-occupancy and co-regulation of PCa genes and metastasis pathways like TGFβ signaling. Mechanistically, GATA2 was cooperative with SMAD4 to enhance TGFβ and AR signaling pathways, and activated the expression of TGFβ1 via directly binding to a distal enhancer of TGFβ1. Strinkingly, GATA2 and SMAD4 globally mediated inherited PCa risk and formed a transcriptional complex with HOXB13 at the PCa risk-associated rs339331/6q22 enhancer, leading to increased expression of the PCa susceptibility gene RFX6.

CONCLUSIONS

Our study prioritizes causal genomic amplification genes with prognostic values in PCa and reveals the pivotal roles of GATA2 in transcriptionally activating the expression of its own and TGFβ1, thereby co-opting to TGFβ1/SMAD4 signaling and RFX6 at 6q22 to modulate PCa predisposition and progression.

Combined CRISPRi and proteomics screening reveal a cohesin-CTCF-bound allele contributing to increased expression of RUVBL1 and prostate cancer progression

Genome-wide association studies along with expression quantitative trait locus (eQTL) mapping have identified hundreds of single-nucleotide polymorphisms (SNPs) and their target genes in prostate cancer (PCa), yet functional characterization of these risk loci remains challenging. To screen for potential regulatory SNPs, we designed a CRISPRi library containing 9,133 guide RNAs (gRNAs) to cover 2,166 candidate SNP loci implicated in PCa and identified 117 SNPs that could regulate 90 genes for PCa cell growth advantage. Among these, rs60464856 was covered by multiple gRNAs significantly depleted in screening (FDR < 0.05). Pooled SNP association analysis in the PRACTICAL and FinnGen cohorts showed significantly higher PCa risk for the rs60464856 G allele (p value = 1.2 × 10-16 and 3.2 × 10-7, respectively). Subsequent eQTL analysis revealed that the G allele is associated with increased RUVBL1 expression in multiple datasets. Further CRISPRi and xCas9 base editing confirmed that the rs60464856 G allele leads to elevated RUVBL1 expression. Furthermore, SILAC-based proteomic analysis demonstrated allelic binding of cohesin subunits at the rs60464856 region, where the HiC dataset showed consistent chromatin interactions in prostate cell lines. RUVBL1 depletion inhibited PCa cell proliferation and tumor growth in a xenograft mouse model. Gene-set enrichment analysis suggested an association of RUVBL1 expression with cell-cycle-related pathways. Increased expression of RUVBL1 and activation of cell-cycle pathways were correlated with poor PCa survival in TCGA datasets. Our CRISPRi screening prioritized about one hundred regulatory SNPs essential for prostate cell proliferation. In combination with proteomics and functional studies, we characterized the mechanistic role of rs60464856 and RUVBL1 in PCa progression.

Multi-omics data analysis reveals the biological implications of alternative splicing events in lung adenocarcinoma

Cancer is characterized by the dysregulation of alternative splicing (AS). However, the comprehensive regulatory mechanisms of AS in lung adenocarcinoma (LUAD) are poorly understood. Here, we displayed the AS landscape in LUAD based on the integrated analyses of LUAD's multi-omics data. We identified 13,995 AS events in 6309 genes as differentially expressed alternative splicing events (DEASEs) mainly covering protein-coding genes. These DEASEs were strongly linked to "cancer hallmarks", such as apoptosis, DNA repair, cell cycle, cell proliferation, angiogenesis, immune response, generation of precursor metabolites and energy, p53 signaling pathway and PI3K-AKT signaling pathway. We further built a regulatory network connecting splicing factors (SFs) and DEASEs. In addition, RNA-binding protein (RBP) mutations that can affect DEASEs were investigated to find some potential cancer drivers. Further association analysis demonstrated that DNA methylation levels were highly correlated with DEASEs. In summary, our results can bring new insight into understanding the mechanism of AS and provide novel biomarkers for personalized medicine of LUAD.

CenhANCER: a comprehensive cancer enhancer database for primary tissues and cell lines

Enhancers, which are key tumorigenic factors with wide applications for subtyping, diagnosis and treatment of cancer, are attracting increasing attention in the cancer research. However, systematic analysis of cancer enhancers poses a challenge due to the lack of integrative data resources, especially those from tumor primary tissues. To provide a comprehensive enhancer profile across cancer types, we developed a cancer enhancer database CenhANCER by curating public resources including all the public H3K27ac ChIP-Seq data from 805 primary tissue samples and 671 cell line samples across 41 cancer types. In total, 57 029 408 typical enhancers, 978 411 super-enhancers and 226 726 enriched transcription factors were identified. We annotated the super-enhancers with chromatin accessibility regions, cancer expression quantitative trait loci (eQTLs), genotype-tissue expression eQTLs and genome-wide association study risk single nucleotide polymorphisms (SNPs) for further functional analysis. The identified enhancers were highly consistent with accessible chromatin regions in the corresponding cancer types, and all the 10 super-enhancer regions identified from one colorectal cancer study were recapitulated in our CenhANCER, both of which testified the high quality of our data. CenhANCER with high-quality cancer enhancer candidates and transcription factors that are potential therapeutic targets across multiple cancer types provides a credible resource for single cancer analysis and for comparative studies of various cancer types. Database URL http://cenhancer.chenzxlab.cn/.

Refined expression quantitative trait locus analysis on adenocarcinoma at the gastroesophageal junction reveals susceptibility and prognostic markers

Objectives: This study aimed to explore cell type level expression quantitative trait loci (eQTL) in adenocarcinoma at the gastroesophageal junction (ACGEJ) and identify susceptibility and prognosis markers. Methods: Whole-genome sequencing (WGS) was performed on 120 paired samples from Chinese ACGEJ patients. Germline mutations were detected by GATK tools. RNA sequencing (RNA-seq) data on ACGEJ samples were taken from our previous studies. Public single-cell RNA sequencing (scRNA-seq) data were used to produce the proportion of epithelial cells. Matrix eQTL and a linear mixed model were used to identify condition-specific cis-eQTLs. The R package coloc was used to perform co-localization analysis with the public data of genome-wide association studies (GWASs). Log-rank and Cox regression tests were used to identify survival-associated eQTL and genes. Functions of candidate risk loci were explored by experimental validation. Results: Refined eQTL analyses of paired ACGEJ samples were performed and 2,036 potential ACGEJ-specific eQTLs with East Asian specificity were identified in total. ACGEJ-gain eQTLs were enriched at promoter regions more than ACGEJ-loss eQTLs. rs658524 was identified as the top eQTL close to the transcription start site of its paired gene (CTSW). rs2240191-RASAL1, rs4236599-FOXP2, rs4947311-PSORS1C1, rs13134812-LOC391674, and rs17508585-CDK13-DT were identified as ACGEJ-specific susceptibility eQTLs. rs309483-LINC01355 was associated with the overall survival of ACGEJ patients. We explored functions of candidate eQTLs such as rs658524, rs309483, rs2240191, and rs4947311 by experimental validation. Conclusion: This study provides new risk loci for ACGEJ susceptibility and effective disease prognosis biomarkers.

Combined CRISPRi and proteomics screening reveal a cohesin-CTCF-bound allele contributing to increased expression of RUVBL1 and prostate cancer progression

Genome-wide association studies along with expression quantitative trait loci (eQTL) mapping have identified hundreds of single nucleotide polymorphisms (SNPs) and their target genes in prostate cancer (PCa), yet functional characterization of these risk loci remains challenging. To screen for potential regulatory SNPs, we designed a CRISPRi library containing 9133 guide RNAs (gRNAs) to target 2,166 candidate SNP sites implicated in PCa and identified 117 SNPs that could regulate 90 genes for PCa cell growth advantage. Among these, rs60464856 was covered by multiple gRNAs significantly depleted in the screening (FDR<0.05). Pooled SNP association analysis in the PRACTICAL and FinnGen cohorts showed significantly higher PCa risk for the rs60464856 G allele (pvalue=1.2E-16 and 3.2E-7). Subsequent eQTL analysis revealed that the G allele is associated with increased RUVBL1 expression in multiple datasets. Further CRISPRi and xCas9 base editing proved the rs60464856 G allele leading to an elevated RUVBL1 expression. Furthermore, SILAC-based proteomic analysis demonstrated allelic binding of cohesin subunits at the rs60464856 region, where HiC dataset showed consistent chromatin interactions in prostate cell lines. RUVBL1 depletion inhibited PCa cell proliferation and tumor growth in xenograft mouse model. Gene set enrichment analysis suggested an association of RUVBL1 expression with cell-cycle-related pathways. An increased expression of RUVBL1 and activations of cell-cycle pathways were correlated with poor PCa survival in TCGA datasets. Together, our CRISPRi screening prioritized about one hundred regulatory SNPs essential for prostate cell proliferation. In combination with proteomics and functional studies, we characterized the mechanistic role of rs60464856 and RUVBL1 in PCa progression.