Cancer-associated dynamics and potential regulators of intronic polyadenylation revealed by IPAFinder using standard RNA-seq data

DIPAN: Detecting personalized intronic polyadenylation derived neoantigens from RNA sequencing data

Intronic polyadenylation (IPA) refers to a particular type of alternative polyadenylation where a gene makes use of a polyadenylation site located within its introns. Aberrant IPA events have been observed in various types of cancer. IPA can produce noncoding transcripts or truncated protein-coding transcripts with altered coding sequences in the resulting protein product. Therefore, IPA events hold the potential to act as a reservoir of tumor neoantigens. Here, we developed a computational method termed DIPAN, which incorporates IPA detection, protein fragmentation, and MHC binding prediction to predict IPA-derived neoantigens. Utilizing RNA-seq from breast cancer cell lines and ovarian cancer clinical samples, we demonstrated the significant contribution of IPA events to the neoantigen repertoire. Through mass spectrometry immunopeptidome analysis, we further illustrated the processing and presentation of IPA-derived neoantigens on the surface of cancer cells. While most IPA-derived neoantigens are sample-specific, shared neoantigens were identified in both cancer cell lines and clinical samples. Furthermore, we demonstrated an association between IPA-derived neoantigen burden and overall survival in cancer patients.

The landscape and clinical relevance of intronic polyadenylation in human cancers

Intronic polyadenylation (IPA) is an RNA 3' end processing event which has been reported to play important roles in cancer development. However, the comprehensive landscape of IPA events across various cancer types is lacking. Here, we apply IPAFinder to identify and quantify IPA events in 10,383 samples covering all 33 cancer types from The Cancer Genome Atlas (TCGA) project. We identify a total of 21,835 IPA events, almost half of which are ubiquitously expressed. We identify 2761 unique dynamically changed IPA events across cancer types. Furthermore, we observe 8855 non-redundant clinically relevant IPA events, which could potentially be used as prognostic indicators. Our analysis also reveals that dynamic IPA usage within cancer signaling pathways may affect drug response. Finally, we develop a user-friendly data portal, IPACancer Atlas (http://www.tingni-lab.com/Pancan_IPA/), to search and explore IPAs in cancer.

Pan-cancer transcriptome analysis reveals widespread regulation through alternative tandem transcription initiation

Abnormal transcription initiation from alternative first exon has been reported to promote tumorigenesis. However, the prevalence and impact of gene expression regulation mediated by alternative tandem transcription initiation were mostly unknown in cancer. Here, we developed a robust computational method to analyze alternative tandem transcription start site (TSS) usage from standard RNA sequencing data. Applying this method to pan-cancer RNA sequencing datasets, we observed widespread dysregulation of tandem TSS usage in tumors, many of which were independent of changes in overall expression level or alternative first exon usage. We showed that the dynamics of tandem TSS usage was associated with epigenomic modulation. We found that significant 5' untranslated region shortening of gene TIMM13 contributed to increased protein production, and up-regulation of TIMM13 by CRISPR-mediated transcriptional activation promoted proliferation and migration of lung cancer cells. Our findings suggest that dysregulated tandem TSS usage represents an addtional layer of cancer-associated transcriptome alterations.

Alternative polyadenylation quantitative trait loci contribute to acute myeloid leukemia risk genes regulation

Acute myeloid leukemia (AML) is a hematopoietic malignancy with a high relapse rate and progressive drug resistance. Alternative polyadenylation (APA) contributes to post-transcriptional dysregulation, but little is known about the association between APA and AML. The APA quantitative trait locus (apaQTL) is a powerful method to investigate the relationship between APA and single nucleotide polymorphisms (SNPs). We quantified APA usage in 195 Chinese AML patients and identified 4922 cis-apaQTLs related to 1875 genes, most of which were newly reported. Cis-apaQTLs may modulate the APA selection of 115 genes through poly(A) signals. Colocalization analysis revealed that cis-apaQTLs colocalized with cis-eQTLs may regulate gene expression by affecting miRNA binding sites or RNA secondary structures. We discovered 207 cis-apaQTLs related to AML risk by comparing genotype frequency with the East Asian healthy controls from the 1000 Genomes Project. Genes with cis-apaQTLs were associated with hematological phenotypes and tumor incidence according to the PHARMGKB and MGI databases. Collectively, we profiled an atlas of cis-apaQTLs in Asian AML patients and found their association with APA selection, gene expression, AML risk, and complex traits. Cis-apaQTLs may provide insights into the regulatory mechanisms related to APA in AML occurrence, progression, and prognosis.

InPACT: a computational method for accurate characterization of intronic polyadenylation from RNA sequencing data

Alternative polyadenylation can occur in introns, termed intronic polyadenylation (IPA), has been implicated in diverse biological processes and diseases, as it can produce noncoding transcripts or transcripts with truncated coding regions. However, a reliable method is required to accurately characterize IPA. Here, we propose a computational method called InPACT, which allows for the precise characterization of IPA from conventional RNA-seq data. InPACT successfully identifies numerous previously unannotated IPA transcripts in human cells, many of which are translated, as evidenced by ribosome profiling data. We have demonstrated that InPACT outperforms other methods in terms of IPA identification and quantification. Moreover, InPACT applied to monocyte activation reveals temporally coordinated IPA events. Further application on single-cell RNA-seq data of human fetal bone marrow reveals the expression of several IPA isoforms in a context-specific manner. Therefore, InPACT represents a powerful tool for the accurate characterization of IPA from RNA-seq data.

FLIBase: a comprehensive repository of full-length isoforms across human cancers and tissues

Regulatory processes at the RNA transcript level play a crucial role in generating transcriptome diversity and proteome composition in human cells, impacting both physiological and pathological states. This study introduces FLIBase (www.FLIBase.org), a specialized database that focuses on annotating full-length isoforms using long-read sequencing techniques. We collected and integrated long-read (351 samples) and short-read (12 469 samples) RNA sequencing data from diverse normal and cancerous human tissues and cells. The current version of FLIBase comprises a total of 983 789 full-length spliced isoforms, identified through long-read sequences and verified using short-read exon-exon splice junctions. Of these, 188 248 isoforms have been annotated, while 795 541 isoforms remain unannotated. By overcoming the limitations of short-read RNA sequencing methods, FLIBase provides an accurate and comprehensive representation of full-length transcripts. These comprehensive annotations empower researchers to undertake various downstream analyses and investigations. Importantly, FLIBase exhibits a significant advantage in identifying a substantial number of previously unannotated isoforms and tumor-specific RNA transcripts. These tumor-specific RNA transcripts have the potential to serve as a source of immunogenic recurrent neoantigens. This remarkable discovery holds tremendous promise for advancing the development of tailored RNA-based diagnostic and therapeutic strategies for various types of human cancer.

Downregulation of apoptotic repressor AVEN exacerbates cardiac injury after myocardial infarction

Myocardial infarction (MI) is a leading cause of heart failure (HF), associated with morbidity and mortality worldwide. As an essential part of gene expression regulation, the role of alternative polyadenylation (APA) in post-MI HF remains elusive. Here, we revealed a global, APA-mediated, 3' untranslated region (3' UTR)-lengthening pattern in both human and murine post-MI HF samples. Furthermore, the 3' UTR of apoptotic repressor gene, AVEN, is lengthened after MI, contributing to its downregulation. AVEN knockdown increased cardiomyocyte apoptosis, whereas restoration of AVEN expression substantially improved cardiac function. Mechanistically, AVEN 3' UTR lengthening provides additional binding sites for miR-30b-5p and miR-30c-5p, thus reducing AVEN expression. Additionally, PABPN1 (poly(A)-binding protein 1) was identified as a potential regulator of AVEN 3' UTR lengthening after MI. Altogether, our findings revealed APA as a unique mechanism regulating cardiac injury in response to MI and also indicated that the APA-regulated gene, AVEN, holds great potential as a critical therapeutic target for treating post-MI HF.

Transcriptome sequencing suggests that pre-mRNA splicing counteracts widespread intronic cleavage and polyadenylation

Alternative splicing (AS) and alternative polyadenylation (APA) are two crucial steps in the post-transcriptional regulation of eukaryotic gene expression. Protocols capturing and sequencing RNA 3'-ends have uncovered widespread intronic polyadenylation (IPA) in normal and disease conditions, where it is currently attributed to stochastic variations in the pre-mRNA processing. Here, we took advantage of the massive amount of RNA-seq data generated by the Genotype Tissue Expression project (GTEx) to simultaneously identify and match tissue-specific expression of intronic polyadenylation sites with tissue-specific splicing. A combination of computational methods including the analysis of short reads with non-templated adenines revealed that APA events are more abundant in introns than in exons. While the rate of IPA in composite terminal exons and skipped terminal exons expectedly correlates with splicing, we observed a considerable fraction of IPA events that lack AS support and attributed them to spliced polyadenylated introns (SPI). We hypothesize that SPIs represent transient byproducts of a dynamic coupling between APA and AS, in which the spliceosome removes the intron while it is being cleaved and polyadenylated. These findings indicate that cotranscriptional pre-mRNA splicing could serve as a rescue mechanism to suppress premature transcription termination at intronic polyadenylation sites.

The peptide woods are lovely, dark and deep: Hunting for novel cancer antigens

Harnessing the patient's immune system to control a tumor is a proven avenue for cancer therapy. T cell therapies as well as therapeutic vaccines, which target specific antigens of interest, are being explored as treatments in conjunction with immune checkpoint blockade. For these therapies, selecting the best suited antigens is crucial. Most of the focus has thus far been on neoantigens that arise from tumor-specific somatic mutations. Although there is clear evidence that T-cell responses against mutated neoantigens are protective, the large majority of these mutations are not immunogenic. In addition, most somatic mutations are unique to each individual patient and their targeting requires the development of individualized approaches. Therefore, novel antigen types are needed to broaden the scope of such treatments. We review high throughput approaches for discovering novel tumor antigens and some of the key challenges associated with their detection, and discuss considerations when selecting tumor antigens to target in the clinic.

The landscape of cryptic antisense transcription in human cancers reveals an oncogenic noncoding RNA in lung cancer

Cryptic transcription initiation has been previously linked to activation of oncogenic transcripts. However, the prevalence and impact of cryptic antisense transcription from the opposite strand of protein-coding genes were mostly unknown in cancer. Applying a robust computational pipeline to publicly available transcriptome and epigenome datasets, we identified hundreds of previously unannotated cryptic antisense polyadenylated transcripts (CAPTs) that were enriched in tumor samples. We showed that the activation of cryptic antisense transcription was associated with increased chromatin accessibility and active histone marks. Accordingly, we found that many of the antisense transcripts were inducible by treatment of epigenetic drugs. Moreover, CRISPR-mediated epigenetic editing assays revealed that transcription of a noncoding RNA LRRK1-CAPT promoted LUSC cell proliferation, suggesting its oncogenic role. Our findings largely expand our understanding of cancer-associated transcription events, which may facilitate the development of novel strategies for cancer diagnosis and treatment.

Alternative Polyadenylation Is a Novel Strategy for the Regulation of Gene Expression in Response to Stresses in Plants

Precursor message RNA requires processing to generate mature RNA. Cleavage and polyadenylation at the 3'-end in the maturation of mRNA is one of key processing steps in eukaryotes. The polyadenylation (poly(A)) tail of mRNA is an essential feature that is required to mediate its nuclear export, stability, translation efficiency, and subcellular localization. Most genes have at least two mRNA isoforms via alternative splicing (AS) or alternative polyadenylation (APA), which increases the diversity of transcriptome and proteome. However, most previous studies have focused on the role of alternative splicing on the regulation of gene expression. In this review, we summarize the recent advances concerning APA in the regulation of gene expression and in response to stresses in plants. We also discuss the mechanisms for the regulation of APA for plants in the adaptation to stress responses, and suggest that APA is a novel strategy for the adaptation to environmental changes and response to stresses in plants.

Intergenic splicing-stimulated transcriptional readthrough is suppressed by nonsense-mediated mRNA decay in Arabidopsis

Recent emerging evidence has shown that readthrough transcripts (RTs), including polycistronic mRNAs, are also transcribed in eukaryotes. However, the post-transcriptional regulation for these remains to be elucidated. Here, we identify 271 polycistronic RT-producing loci in Arabidopsis. Increased accumulation of RTs is detected in the nonsense-mediated mRNA decay (NMD)-deficient mutants compared with wild type, and the second open reading frames (ORFs) of bicistronic mRNAs are rarely translated in contrast to the first ORFs. Intergenic splicing (IS) events which occur between first and second genes are seen in 158 RTs. Splicing inhibition assays suggest that IS eliminates the chance of transcription termination at the polyadenylation sites of the first gene and promotes accumulation of RTs. These results indicate that RTs arise from genes whose transcription termination is relatively weak or attenuated by IS, but NMD selectively degrades them. Ultimately, this report presents a eukaryotic strategy for RNA metabolism.

A Survey on Methods for Predicting Polyadenylation Sites from DNA Sequences, Bulk RNA-seq, and Single-cell RNA-seq

Alternative polyadenylation (APA) plays important roles in modulating mRNA stability, translation, and subcellular localization, and contributes extensively to shaping eukaryotic transcriptome complexity and proteome diversity. Identification of poly(A) sites (pAs) on a genome-wide scale is a critical step toward understanding the underlying mechanism of APA-mediated gene regulation. A number of established computational tools have been proposed to predict pAs from diverse genomic data. Here we provided an exhaustive overview of computational approaches for predicting pAs from DNA sequences, bulk RNA sequencing (RNA-seq) data, and single-cell RNA sequencing (scRNA-seq) data. Particularly, we examined several representative tools using bulk RNA-seq and scRNA-seq data from peripheral blood mononuclear cells and put forward operable suggestions on how to assess the reliability of pAs predicted by different tools. We also proposed practical guidelines on choosing appropriate methods applicable to diverse scenarios. Moreover, we discussed in depth the challenges in improving the performance of pA prediction and benchmarking different methods. Additionally, we highlighted outstanding challenges and opportunities using new machine learning and integrative multi-omics techniques, and provided our perspective on how computational methodologies might evolve in the future for non-3' untranslated region, tissue-specific, cross-species, and single-cell pA prediction.

ipaQTL-atlas: an atlas of intronic polyadenylation quantitative trait loci across human tissues

Functional interpretation of disease-associated non-coding variants remains a significant challenge in the post-GWAS era. Our recent study has identified 3'UTR alternative polyadenylation (APA) quantitative trait loci (3'aQTLs) and connects APA events with QTLs as a major driver of human traits and diseases. Besides 3'UTR, APA events can also occur in intron regions, and increasing evidence has connected intronic polyadenylation with disease risk. However, systematic investigation of the roles of intronic polyadenylation in human diseases remained challenging due to the lack of a comprehensive database across a variety of human tissues. Here, we developed ipaQTL-atlas (http://bioinfo.szbl.ac.cn/ipaQTL) as the first comprehensive portal for intronic polyadenylation. The ipaQTL-atlas is based on the analysis of 15 170 RNA-seq data from 838 individuals across 49 Genotype-Tissue Expression (GTEx v8) tissues and contains ∼0.98 million SNPs associated with intronic APA events. It provides an interface for ipaQTLs search, genome browser, boxplots, and data download, as well as the visualization of GWAS and ipaQTL colocalization results. ipaQTL-atlas provides a one-stop portal to access intronic polyadenylation information and could significantly advance the discovery of APA-associated disease susceptibility genes.

Dysregulated RNA polyadenylation contributes to metabolic impairment in non-alcoholic fatty liver disease

Pre-mRNA processing is an essential mechanism for the generation of mature mRNA and the regulation of gene expression in eukaryotic cells. While defects in pre-mRNA processing have been implicated in a number of diseases their involvement in metabolic pathologies is still unclear. Here, we show that both alternative splicing and alternative polyadenylation, two major steps in pre-mRNA processing, are significantly altered in non-alcoholic fatty liver disease (NAFLD). Moreover, we find that Serine and Arginine Rich Splicing Factor 10 (SRSF10) binding is enriched adjacent to consensus polyadenylation motifs and its expression is significantly decreased in NAFLD, suggesting a role mediating pre-mRNA dysregulation in this condition. Consistently, inactivation of SRSF10 in mouse and human hepatocytes in vitro, and in mouse liver in vivo, was found to dysregulate polyadenylation of key metabolic genes such as peroxisome proliferator-activated receptor alpha (PPARA) and exacerbate diet-induced metabolic dysfunction. Collectively our work implicates dysregulated pre-mRNA polyadenylation in obesity-induced liver disease and uncovers a novel role for SRSF10 in this process.

Intron mutations and early transcription termination in Duchenne and Becker muscular dystrophy

DMD pathogenic variants for Duchenne and Becker muscular dystrophy are detectable with high sensitivity by standard clinical exome analyses of genomic DNA. However, up to 7% of DMD mutations are deep intronic and analysis of muscle-derived RNA is an important diagnostic step for patients who have negative genomic testing but abnormal dystrophin expression in muscle. In this study, muscle biopsies were evaluated from 19 patients with clinical features of a dystrophinopathy, but negative clinical DMD mutation analysis. Reverse transcription-polymerase chain reaction or high-throughput RNA sequencing methods identified 19 mutations with one of three pathogenic pseudoexon types: deep intronic point mutations, deletions or insertions, and translocations. In association with point mutations creating intronic splice acceptor sites, we observed the first examples of DMD pseudo 3'-terminal exon mutations causing high efficiency transcription termination within introns. This connection between splicing and premature transcription termination is reminiscent of U1 snRNP-mediating telescripting in sustaining RNA polymerase II elongation across large genes, such as DMD. We propose a novel classification of three distinct types of mutations identifiable by muscle RNA analysis, each of which differ in potential treatment approaches. Recognition and appropriate characterization may lead to therapies directed toward full-length dystrophin expression for some patients.

Heterogeneity and Molecular Markers for CNS Glial Cells Revealed by Single-Cell Transcriptomics

Glial cells, including astrocytes, oligodendrocytes, and microglia, are the major components in the central nervous system (CNS). Studies have revealed the heterogeneity of each glial cell type and that they each may play distinct roles in physiological processes and/or neurological diseases. Single-cell sequencing (scRNA-seq) technology developed in recent years has extended our understanding of glial cell heterogeneity from the perspective of transcriptome profiling. This review summarizes the marker genes of major glial cells in the CNS and reveals their heterogeneity in different species, CNS regions, developmental stages, and pathological states (Alzheimer's disease and spinal cord injury), expanding our knowledge of glial cell heterogeneity on both molecular and functional levels.

Comprehensive characterization of somatic variants associated with intronic polyadenylation in human cancers

Somatic single nucleotide variants (SNVs) in cancer genome affect gene expression through various mechanisms depending on their genomic location. While somatic SNVs near canonical splice sites have been reported to cause abnormal splicing of cancer-related genes, whether these SNVs can affect gene expression through other mechanisms remains an open question. Here, we analyzed RNA sequencing and exome data from 4,998 cancer patients covering ten cancer types and identified 152 somatic SNVs near splice sites that were associated with abnormal intronic polyadenylation (IPA). IPA-associated somatic variants favored the localization near the donor splice sites compared to the acceptor splice sites. A proportion of SNV-associated IPA events overlapped with premature cleavage and polyadenylation events triggered by U1 small nuclear ribonucleoproteins (snRNP) inhibition. GC content, intron length and polyadenylation signal were three genomic features that differentiated between SNV-associated IPA and intron retention. Notably, IPA-associated SNVs were enriched in tumor suppressor genes (TSGs), including the well-known TSGs such as PTEN and CDH1 with recurrent SNV-associated IPA events. Minigene assay confirmed that SNVs from PTEN, CDH1, VEGFA, GRHL2, CUL3 and WWC2 could lead to IPA. This work reveals that IPA acts as a novel mechanism explaining the functional consequence of somatic SNVs in human cancer.