An integrative model for alternative polyadenylation, IntMAP, delineates mTOR-modulated endoplasmic reticulum stress response

Melatonin alleviates endoplasmic reticulum stress to improve ovarian function by regulating the mTOR pathway in aged laying hens

Granular cell apoptosis is a key factor leading to follicular atresia and decreased laying rate in aged laying hens. Endoplasmic reticulum stress (ERS) induced cell apoptosis is a new type of apoptosis pathway. Previous studies have shown that the ERS pathway is involved in the regulation of follicular development and atresia, and can be regulated by mTOR. Melatonin (MEL) can protect the normal development of follicles, but the precise mechanism by which MEL regulates follicular development is not yet clear. So, we investigated the potential relationship between MEL and ERS and mTOR signaling pathway in vivo through intraperitoneal injection of MEL in aged laying hens. The results show that the laying rate, ovarian follicle number, plasma MEL, E2, LH, FSH concentrations, as well as the mRNA expression of mTOR signaling-associated genes TSC1, TSC2, mTOR, 4E-BP1, and S6K in old later-period chicken control (Old-CN) group was significantly decreased (P < 0.01). In contrast, the ERS-related of plasma and granular cell layer mRNA expression of Grp78, CHOP, and Caspase-3 was significantly increased (P < 0.01). While both of the effects were reversed by MEL. Then, aging granulosa cells were treated with MEL in vitro, followed by RNA seq analysis, and it was found that 259 and 322 genes were upregulated and downregulated. After performing GO enrichment analysis, it was found that DEGs significantly contribute to the biological processes including cell growth and apoptosis. Using pathway enrichment analysis, we found significant overrepresentation of cellular processes related to mTOR signaling and endoplasmic reticulum (ER) stress, involving genes such as GRB10, SGK1, PRKCA, RPS6KA2, RAF1, PIK3R3, FOXO1, DERL3, HMOX1, TLR7, VAMP7 and INSIG2. The obtained results of RT-PCR showed consistency with the RNA-Seq data. In summary, the underlined results revealed that MEL has significantly contributed to follicular development via activating the mTOR signaling pathway-related genes and alleviating ERS-related genes in laying hens. The current study provides a theoretical background for enhancing the egg-laying capability of hens and also providing a basis for elucidating the molecular mechanism of follicular selection.

The emerging theme of 3'UTR mRNA isoform regulation in reprogramming of cell metabolism

The 3' untranslated region (3'UTR) of mRNA plays a key role in the post-transcriptional regulation of gene expression. Most eukaryotic protein-coding genes express 3'UTR isoforms owing to alternative cleavage and polyadenylation (APA). The 3'UTR isoform expression profile of a cell changes in cell proliferation, differentiation, and stress conditions. Here, we review the emerging theme of regulation of 3'UTR isoforms in cell metabolic reprogramming, focusing on cell growth and autophagy responses through the mTOR pathway. We discuss regulatory events that converge on the Cleavage Factor I complex, a master regulator of APA in 3'UTRs, and recent understandings of isoform-specific m6A modification and endomembrane association in determining differential metabolic fates of 3'UTR isoforms.

CEBPG suppresses ferroptosis through transcriptional control of SLC7A11 in ovarian cancer

BACKGROUND

Ovarian cancer (OC) has high mortality and poor prognosis for lacking of specific biomarkers and typical clinical symptoms in the early stage. CEBPG is an important regulator in tumor development, yet it is unclear exactly how it contributes to the progression of OC.

METHODS

TCGA and tissue microarrays with immunohistochemical staining (IHC) were used to examine CEBPG expression in OC. A variety of in vitro assays were conducted, including colony formation, proliferation, migration, and invasion. The orthotopic OC mouse model was established for in vivo studies. Ferroptosis was detected by observing mitochondrial changes with electron microscopy, detecting ROS expression, and detecting cell sensitivity to drugs by CCK8 assay. The interaction between CEBPG and SLC7A11 was confirmed by CUT&Tag and dual luciferase reporter assays.

RESULTS

A significantly higher expression level of CEBPG in OC when compared with benign tissues of ovary, and that high CEBPG expression level was also tightly associated with poor prognosis of patients diagnosed with OC, as determined by analysis of datasets and patient samples. Conversely, knockdown of CEBPG inhibited OC progression using experiments of OC cell lines and in vivo orthotopic OC-bearing mouse model. Importantly, CEBPG was identified as a new participator mediating ferroptosis evasion in OC cells using RNA-sequencing, which could contribute to OC progression. The CUT&Tag and dua luciferase reporter assays further revealed the inner mechanism that CEBPG regulated OC cell ferroptosis through transcriptional control of SLC7A11.

CONCLUSIONS

Our findings established CEBPG as a novel transcriptional regulator of OC ferroptosis, with potential value in predicting clinical outcomes and as a therapeutic candidate.

Regulation of the hypertonic stress response by the 3' mRNA cleavage and polyadenylation complex

Maintenance of osmotic homeostasis is one of the most aggressively defended homeostatic set points in physiology. One major mechanism of osmotic homeostasis involves the upregulation of proteins that catalyze the accumulation of solutes called organic osmolytes. To better understand how osmolyte accumulation proteins are regulated, we conducted a forward genetic screen in Caenorhabditis elegans for mutants with no induction of osmolyte biosynthesis gene expression (Nio mutants). The nio-3 mutant encoded a missense mutation in cpf-2/CstF64, while the nio-7 mutant encoded a missense mutation in symk-1/Symplekin. Both cpf-2 and symk-1 are nuclear components of the highly conserved 3' mRNA cleavage and polyadenylation complex. cpf-2 and symk-1 block the hypertonic induction of gpdh-1 and other osmotically induced mRNAs, suggesting they act at the transcriptional level. We generated a functional auxin-inducible degron (AID) allele for symk-1 and found that acute, post-developmental degradation in the intestine and hypodermis was sufficient to cause the Nio phenotype. symk-1 and cpf-2 exhibit genetic interactions that strongly suggest they function through alterations in 3' mRNA cleavage and/or alternative polyadenylation. Consistent with this hypothesis, we find that inhibition of several other components of the mRNA cleavage complex also cause a Nio phenotype. cpf-2 and symk-1 specifically affect the osmotic stress response since heat shock-induced upregulation of a hsp-16.2::GFP reporter is normal in these mutants. Our data suggest a model in which alternative polyadenylation of 1 or more mRNAs is essential to regulate the hypertonic stress response.

mTOR Contributes to the Proteome Diversity through Transcriptome-Wide Alternative Splicing

The mammalian target of rapamycin (mTOR) pathway is crucial in energy metabolism and cell proliferation. Previously, we reported transcriptome-wide 3′-untranslated region (UTR) shortening by alternative polyadenylation upon mTOR activation and its impact on the proteome. Here, we further interrogated the mTOR-activated transcriptome and found that hyperactivation of mTOR promotes transcriptome-wide exon skipping/exclusion, producing short isoform transcripts from genes. This widespread exon skipping confers multifarious regulations in the mTOR-controlled functional proteomics: AS in coding regions widely affects the protein length and functional domains. They also alter the half-life of proteins and affect the regulatory post-translational modifications. Among the RNA processing factors differentially regulated by mTOR signaling, we found that SRSF3 mechanistically facilitates exon skipping in the mTOR-activated transcriptome. This study reveals a role of mTOR in AS regulation and demonstrates that widespread AS is a multifaceted modulator of the mTOR-regulated functional proteome.

APA-Scan: detection and visualization of 3'-UTR alternative polyadenylation with RNA-seq and 3'-end-seq data

Background

The eukaryotic genome is capable of producing multiple isoforms from a gene by alternative polyadenylation (APA) during pre-mRNA processing. APA in the 3′-untranslated region (3′-UTR) of mRNA produces transcripts with shorter or longer 3′-UTR. Often, 3′-UTR serves as a binding platform for microRNAs and RNA-binding proteins, which affect the fate of the mRNA transcript. Thus, 3′-UTR APA is known to modulate translation and provides a mean to regulate gene expression at the post-transcriptional level. Current bioinformatics pipelines have limited capability in profiling 3′-UTR APA events due to incomplete annotations and a low-resolution analyzing power: widely available bioinformatics pipelines do not reference actionable polyadenylation (cleavage) sites but simulate 3′-UTR APA only using RNA-seq read coverage, causing false positive identifications. To overcome these limitations, we developed APA-Scan, a robust program that identifies 3′-UTR APA events and visualizes the RNA-seq short-read coverage with gene annotations.

Methods

APA-Scan utilizes either predicted or experimentally validated actionable polyadenylation signals as a reference for polyadenylation sites and calculates the quantity of long and short 3′-UTR transcripts in the RNA-seq data. APA-Scan works in three major steps: (i) calculate the read coverage of the 3′-UTR regions of genes; (ii) identify the potential APA sites and evaluate the significance of the events among two biological conditions; (iii) graphical representation of user specific event with 3′-UTR annotation and read coverage on the 3′-UTR regions. APA-Scan is implemented in Python3. Source code and a comprehensive user’s manual are freely available at https://github.com/compbiolabucf/APA-Scan.

Result

APA-Scan was applied to both simulated and real RNA-seq datasets and compared with two widely used baselines DaPars and APAtrap. In simulation APA-Scan significantly improved the accuracy of 3′-UTR APA identification compared to the other baselines. The performance of APA-Scan was also validated by 3′-end-seq data and qPCR on mouse embryonic fibroblast cells. The experiments confirm that APA-Scan can detect unannotated 3′-UTR APA events and improve genome annotation.

Conclusion

APA-Scan is a comprehensive computational pipeline to detect transcriptome-wide 3′-UTR APA events. The pipeline integrates both RNA-seq and 3′-end-seq data information and can efficiently identify the significant events with a high-resolution short reads coverage plots.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04939-w.

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.

In silico model for miRNA-mediated regulatory network in cancer

Deregulation of gene expression is associated with the pathogenesis of numerous human diseases including cancer. Current data analyses on gene expression are mostly focused on differential gene/transcript expression in big data-driven studies. However, a poor connection to the proteome changes is a widespread problem in current data analyses. This is partly due to the complexity of gene regulatory pathways at the post-transcriptional level. In this study, we overcome these limitations and introduce a graph-based learning model, PTNet, which simulates the microRNAs (miRNAs) that regulate gene expression post-transcriptionally in silico. Our model does not require large-scale proteomics studies to measure the protein expression and can successfully predict the protein levels by considering the miRNA-mRNA interaction network, the mRNA expression, and the miRNA expression. Large-scale experiments on simulations and real cancer high-throughput datasets using PTNet validated that (i) the miRNA-mediated interaction network affects the abundance of corresponding proteins and (ii) the predicted protein expression has a higher correlation with the proteomics data (ground-truth) than the mRNA expression data. The classification performance also shows that the predicted protein expression has an improved prediction power on cancer outcomes compared to the prediction done by the mRNA expression data only or considering both mRNA and miRNA. Availability: PTNet toolbox is available at http://github.com/CompbioLabUCF/PTNet.

Alternative polyadenylation: methods, mechanism, function, and role in cancer

Occurring in over 60% of human genes, alternative polyadenylation (APA) results in numerous transcripts with differing 3'ends, thus greatly expanding the diversity of mRNAs and of proteins derived from a single gene. As a key molecular mechanism, APA is involved in various gene regulation steps including mRNA maturation, mRNA stability, cellular RNA decay, and protein diversification. APA is frequently dysregulated in cancers leading to changes in oncogenes and tumor suppressor gene expressions. Recent studies have revealed various APA regulatory mechanisms that promote the development and progression of a number of human diseases, including cancer. Here, we provide an overview of four types of APA and their impacts on gene regulation. We focus particularly on the interaction of APA with microRNAs, RNA binding proteins and other related factors, the core pre-mRNA 3'end processing complex, and 3'UTR length change. We also describe next-generation sequencing methods and computational tools for use in poly(A) signal detection and APA repositories and databases. Finally, we summarize the current understanding of APA in cancer and provide our vision for future APA related research.

CEBPG promotes esophageal squamous cell carcinoma progression by enhancing PI3K-AKT signaling

CCAAT/enhancer binding proteins (CEBPs, including CEBPA, CEBPB, CEBPD, CEBPE, CEBPG, and CEBPZ) play critical roles in a variety of physiological and pathological processes. However, the molecular characteristics and biological significance of CEBPs in esophageal squamous cell carcinoma (ESCC) have rarely been reported. Here, we show that most of the CEBPs are upregulated and accompanied with copy number amplifications in ESCC. Of note, high CEBPG expression is regulated by the ESCC specific transcription factor TP63 and serves as a prognostic factor for poor survival in ESCC patients. Functionally, CEBPG significantly promotes the proliferation and migration of ESCC cells both in vitro and in vivo. Mechanistically, CEBPG activates the PI3K-AKT signaling pathway through directly binding to distal enhancers and/or promoters of genes involved in this pathway, including genes of CCND1, MYC, CDK2, etc. These findings provide new insights into CEBPs dysregulation in ESCC and elucidate a crucial role for CEBPG in the progression of ESCC, highlighting its potential therapeutic value for ESCC treatment.

Discovery of alternative polyadenylation dynamics from single cell types

Alternative polyadenylation (APA) occurs in the process of mRNA maturation by adding a poly(A) tail at different locations, resulting increased diversity of mRNA isoforms and contributing to the complexity of gene regulatory network. Benefit from the development of high-throughput sequencing technologies, we could now delineate APA profiles of transcriptomes at an unprecedented pace. Especially the single cell RNA sequencing (scRNA-seq) technologies provide us opportunities to interrogate biological details of diverse and rare cell types. Despite increasing evidence showing that APA is involved in the cell type-specific regulation and function, efficient and specific laboratory methods for capturing poly(A) sites at single cell resolution are underdeveloped to date. In this review, we summarize existing experimental and computational methods for the identification of APA dynamics from diverse single cell types. A future perspective is also provided.

Platform-integrated mRNA isoform quantification

MOTIVATION

Accurate estimation of transcript isoform abundance is critical for downstream transcriptome analyses and can lead to precise molecular mechanisms for understanding complex human diseases, like cancer. Simplex mRNA Sequencing (RNA-Seq) based isoform quantification approaches are facing the challenges of inherent sampling bias and unidentifiable read origins. A large-scale experiment shows that the consistency between RNA-Seq and other mRNA quantification platforms is relatively low at the isoform level compared to the gene level. In this project, we developed a platform-integrated model for transcript quantification (IntMTQ) to improve the performance of RNA-Seq on isoform expression estimation. IntMTQ, which benefits from the mRNA expressions reported by the other platforms, provides more precise RNA-Seq-based isoform quantification and leads to more accurate molecular signatures for disease phenotype prediction.

RESULTS

In the experiments to assess the quality of isoform expression estimated by IntMTQ, we designed three tasks for clustering and classification of 46 cancer cell lines with four different mRNA quantification platforms, including newly developed NanoString's nCounter technology. The results demonstrate that the isoform expressions learned by IntMTQ consistently provide more and better molecular features for downstream analyses compared with five baseline algorithms which consider RNA-Seq data only. An independent RT-qPCR experiment on seven genes in twelve cancer cell lines showed that the IntMTQ improved overall transcript quantification. The platform-integrated algorithms could be applied to large-scale cancer studies, such as The Cancer Genome Atlas (TCGA), with both RNA-Seq and array-based platforms available.

AVAILABILITY AND IMPLEMENTATION

Source code is available at: https://github.com/CompbioLabUcf/IntMTQ.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

mTOR-regulated U2af1 tandem exon splicing specifies transcriptome features for translational control

U2 auxiliary factor 1 (U2AF1) functions in 3′-splice site selection during pre-mRNA processing. Alternative usage of duplicated tandem exons in U2AF1 produces two isoforms, U2AF1a and U2AF1b, but their functional differences are unappreciated due to their homology. Through integrative approaches of genome editing, customized-transcriptome profiling and crosslinking-mediated interactome analyses, we discovered that the expression of U2AF1 isoforms is controlled by mTOR and they exhibit a distinctive molecular profile for the splice site and protein interactomes. Mechanistic dissection of mutually exclusive alternative splicing events revealed that U2AF1 isoforms’ inherent differential preferences of nucleotide sequences and their stoichiometry determine the 3′-splice site. Importantly, U2AF1a-driven transcriptomes feature alternative splicing events in the 5′-untranslated region (5′-UTR) that are favorable for translation. These findings unveil distinct roles of duplicated tandem exon-derived U2AF1 isoforms in the regulation of the transcriptome and suggest U2AF1a-driven 5′-UTR alternative splicing as a molecular mechanism of mTOR-regulated translational control.

Alternative polyadenylation of mRNA and its role in cancer

Alternative polyadenylation (APA) is a molecular process that generates diversity at the 3′ end of RNA polymerase II transcripts from over 60% of human genes. APA is derived from the existence of multiple polyadenylation signals (PAS) within the same transcript, and results in the differential inclusion of sequence information at the 3′ end. While APA can occur between two PASs allowing for generation of transcripts with distinct coding potential from a single gene, most APA occurs within the untranslated region (3′UTR) and changes the length and content of these non-coding sequences. APA within the 3′UTR can have tremendous impact on its regulatory potential of the mRNA through a variety of mechanisms, and indeed this layer of gene expression regulation has profound impact on processes vital to cell growth and development. Recent studies have particularly highlighted the importance of APA dysregulation in cancer onset and progression. Here, we review the current knowledge of APA and its impacts on mRNA stability, translation, localization and protein localization. We also discuss the implications of APA dysregulation in cancer research and therapy.

mTOR-coordinated Post-Transcriptional Gene Regulations: from Fundamental to Pathogenic Insights

Post-transcriptional regulations of mRNA transcripts such as alternative splicing and alternative polyadenylation can affect the expression of genes without changing the transcript levels. Recent studies have demonstrated that these post-transcriptional events can have significant physiological impacts on various biological systems and play important roles in the pathogenesis of a number of diseases, including cancers. Nevertheless, how cellular signaling pathways control these post-transcriptional processes in cells are not very well explored in the field yet. The mammalian target of rapamycin complex 1 (mTORC1) pathway plays a key role in sensing cellular nutrient and energy status and regulating the proliferation and growth of cells by controlling various anabolic and catabolic processes. Dysregulation of mTORC1 pathway can tip the metabolic balance of cells and is associated with a number of pathological conditions, including various types of cancers, diabetes, and cardiovascular diseases. Numerous reports have shown that mTORC1 controls its downstream pathways through translational and/or transcriptional regulation of the expression of key downstream effectors. And, recent studies have also shown that mTORC1 can control downstream pathways via post-transcriptional regulations. In this review, we will discuss the roles of post-transcriptional processes in gene expression regulations and how mTORC1-mediated post-transcriptional regulations contribute to cellular physiological changes. We highlight post-transcriptional regulation as an additional layer of gene expression control by mTORC1 to steer cellular biology. These emphasize the importance of studying post-transcriptional events in transcriptome datasets for gaining a fuller understanding of gene expression regulations in the biological systems of interest.

Alternative polyadenylation produces multiple 3' untranslated regions of odorant receptor mRNAs in mouse olfactory sensory neurons

BACKGROUND

Odorant receptor genes constitute the largest gene family in mammalian genomes and this family has been extensively studied in several species, but to date far less attention has been paid to the characterization of their mRNA 3' untranslated regions (3'UTRs). Given the increasing importance of UTRs in the understanding of RNA metabolism, and the growing interest in alternative polyadenylation especially in the nervous system, we aimed at identifying the alternative isoforms of odorant receptor mRNAs generated through 3'UTR variation.

RESULTS

We implemented a dedicated pipeline using IsoSCM instead of Cufflinks to analyze RNA-Seq data from whole olfactory mucosa of adult mice and obtained an extensive description of the 3'UTR isoforms of odorant receptor mRNAs. To validate our bioinformatics approach, we exhaustively analyzed the 3'UTR isoforms produced from 2 pilot genes, using molecular approaches including northern blot and RNA ligation mediated polyadenylation test. Comparison between datasets further validated the pipeline and confirmed the alternative polyadenylation patterns of odorant receptors. Qualitative and quantitative analyses of the annotated 3' regions demonstrate that 1) Odorant receptor 3'UTRs are longer than previously described in the literature; 2) More than 77% of odorant receptor mRNAs are subject to alternative polyadenylation, hence generating at least 2 detectable 3'UTR isoforms; 3) Splicing events in 3'UTRs are restricted to a limited subset of odorant receptor genes; and 4) Comparison between male and female data shows no sex-specific differences in odorant receptor 3'UTR isoforms.

CONCLUSIONS

We demonstrated for the first time that odorant receptor genes are extensively subject to alternative polyadenylation. This ground-breaking change to the landscape of 3'UTR isoforms of Olfr mRNAs opens new avenues for investigating their respective functions, especially during the differentiation of olfactory sensory neurons.

Alternative cleavage and polyadenylation in health and disease

Most human genes have multiple sites at which RNA 3' end cleavage and polyadenylation can occur, enabling the expression of distinct transcript isoforms under different conditions. Novel methods to sequence RNA 3' ends have generated comprehensive catalogues of polyadenylation (poly(A)) sites; their analysis using innovative computational methods has revealed how poly(A) site choice is regulated by core RNA 3' end processing factors, such as cleavage factor I and cleavage and polyadenylation specificity factor, as well as by other RNA-binding proteins, particularly splicing factors. Here, we review the experimental and computational methods that have enabled the global mapping of mRNA and of long non-coding RNA 3' ends, quantification of the resulting isoforms and the discovery of regulators of alternative cleavage and polyadenylation (APA). We highlight the different types of APA-derived isoforms and their functional differences, and illustrate how APA contributes to human diseases, including cancer and haematological, immunological and neurological diseases.

The histone variant H2A.Z in gene regulation

The histone variant H2A.Z is involved in several processes such as transcriptional control, DNA repair, regulation of centromeric heterochromatin and, not surprisingly, is implicated in diseases such as cancer. Here, we review the recent developments on H2A.Z focusing on its role in transcriptional activation and repression. H2A.Z, as a replication-independent histone, has been studied in several model organisms and inducible mammalian model systems. Its loading machinery and several modifying enzymes have been recently identified, and some of the long-standing discrepancies in transcriptional activation and/or repression are about to be resolved. The buffering functions of H2A.Z, as supported by genome-wide localization and analyzed in several dynamic systems, are an excellent example of transcriptional control. Posttranslational modifications such as acetylation and ubiquitination of H2A.Z, as well as its specific binding partners, are in our view central players in the control of gene expression. Understanding the key-mechanisms in either turnover or stabilization of H2A.Z-containing nucleosomes as well as defining the H2A.Z interactome will pave the way for therapeutic applications in the future.

DPAC: A Tool for Differential Poly(A)-Cluster Usage from Poly(A)-Targeted RNAseq Data

Poly(A)-tail targeted RNAseq approaches, such as 3′READS, PAS-Seq and Poly(A)-ClickSeq, are becoming popular alternatives to random-primed RNAseq to focus sequencing reads just to the 3′ ends of polyadenylated RNAs to identify poly(A)-sites and characterize changes in their usage. Additionally, we and others have demonstrated that these approaches perform similarly to other RNAseq strategies for differential gene expression analysis, while saving on the volume of sequencing data required and providing a simpler library synthesis strategy. Here, we present DPAC (Differential Poly(A)-Clustering); a streamlined pipeline for the preprocessing of poly(A)-tail targeted RNAseq data, mapping of poly(A)-sites, poly(A)-site clustering and annotation, and determination of differential poly(A)-cluster usage using DESeq2. Changes in poly(A)-cluster usage is simultaneously used to report differential gene expression, differential terminal exon usage and alternative polyadenylation (APA).