Mechanisms and consequences of alternative polyadenylation

Alternative polyadenylation: An enigma of transcript length variation in health and disease

Alternative polyadenylation (APA) is a molecular mechanism during a pre-mRNA processing that involves usage of more than one polyadenylation site (PA-site) generating transcripts of varying length from a single gene. The location of a PA-site affects transcript length and coding potential of an mRNA contributing to both mRNA and protein diversification. This variation in the transcript length affects mRNA stability and translation, mRNA subcellular and tissue localization, and protein function. APA is now considered as an important regulatory mechanism in the pathophysiology of human diseases. An important consequence of the changes in the length of 3'-untranslated region (UTR) from disease-induced APA is altered protein expression. Yet, the relationship between 3'-UTR length and protein expression remains a paradox in a majority of diseases. Here, we review occurrence of APA, mechanism of PA-site selection, and consequences of transcript length variation in different diseases. Emerging evidence reveals coordinated involvement of core RNA processing factors including poly(A) polymerases in the PA-site selection in diseases-associated APAs. Targeting such APA regulators will be therapeutically significant in combating drug resistance in cancer and other complex diseases. This article is categorized under: RNA Processing > 3' End Processing RNA in Disease and Development > RNA in Disease Translation > Regulation.

Seeking a Role for Translational Control by Alternative Polyadenylation in Saccharomyces cerevisiae

Alternative polyadenylation (APA) represents an important mechanism for regulating isoform-specific translation efficiency, stability, and localisation. Though some progress has been made in understanding its consequences in metazoans, the role of APA in the model organism Saccharomyces cerevisiae remains a relative mystery because, despite abundant studies on the translational state of mRNA, none differentiate mRNA isoforms’ alternative 3′-end. This review discusses the implications of alternative polyadenylation in S. cerevisiae using other organisms to draw inferences. Given the foundational role that research in this yeast has played in the discovery of the mechanisms of cleavage and polyadenylation and in the drivers of APA, it is surprising that such an inference is required. However, because advances in ribosome profiling are insensitive to APA, how it impacts translation is still unclear. To bridge the gap between widespread observed APA and the discovery of any functional consequence, we also provide a review of the experimental techniques used to uncover the functional importance of 3′ UTR isoforms on translation.

Regulation of closely juxtaposed proto-oncogene c-fms and HMGXB3 gene expression by mRNA 3' end polymorphism in breast cancer cells

Sense-antisense mRNA pairs generated by convergent transcription is a way of gene regulation. c-fms gene is closely juxtaposed to the HMGXB3 gene in the opposite orientation, in chromosome 5. The intergenic region (IR) between c-fms and HMGXB3 genes is 162 bp. We found that a small portion (∼4.18%) of HMGXB3 mRNA is transcribed further downstream, including the end of the c-fms gene generating antisense mRNA against c-fms mRNA. Similarly, a small portion (∼1.1%) of c-fms mRNA is transcribed further downstream, including the end of the HMGXB3 gene generating antisense mRNA against the HMGXB3 mRNA. Insertion of the strong poly(A) signal sequence in the IR results in decreased c-fms and HMGXB3 antisense mRNAs, resulting in up-regulation of both c-fms and HMGXB3 mRNA expression. miR-324-5p targets HMGXB3 mRNA 3' UTR, and as a result, regulates c-fms mRNA expression. HuR stabilizes c-fms mRNA, and as a result, down-regulates HMGXB3 mRNA expression. UALCAN analysis indicates that the expression pattern between c-fms and HMGXB3 proteins are opposite in vivo in breast cancer tissues. Together, our results indicate that the mRNA encoded by the HMGXB3 gene can influence the expression of adjacent c-fms mRNA, or vice versa.

The landscape of alternative polyadenylation in single cells of the developing mouse embryo

3′ untranslated regions (3′ UTRs) post-transcriptionally regulate mRNA stability, localization, and translation rate. While 3′-UTR isoforms have been globally quantified in limited cell types using bulk measurements, their differential usage among cell types during mammalian development remains poorly characterized. In this study, we examine a dataset comprising ~2 million nuclei spanning E9.5–E13.5 of mouse embryonic development to quantify transcriptome-wide changes in alternative polyadenylation (APA). We observe a global lengthening of 3′ UTRs across embryonic stages in all cell types, although we detect shorter 3′ UTRs in hematopoietic lineages and longer 3′ UTRs in neuronal cell types within each stage. An analysis of RNA-binding protein (RBP) dynamics identifies ELAV-like family members, which are concomitantly induced in neuronal lineages and developmental stages experiencing 3′-UTR lengthening, as putative regulators of APA. By measuring 3′-UTR isoforms in an expansive single cell dataset, our work provides a transcriptome-wide and organism-wide map of the dynamic landscape of alternative polyadenylation during mammalian organogenesis.

Alternative polyadenylation regulates localization, half-life and translation of mRNA isoforms. Here the authors investigate alternative polyadenylation using single cell RNA sequencing data from mouse embryos and identify 3’-UTR isoforms that are regulated across cell types and developmental time.

The endogenous mex-3 3´UTR is required for germline repression and contributes to optimal fecundity in C. elegans

RNA regulation is essential to successful reproduction. Messenger RNAs delivered from parent to progeny govern early embryonic development. RNA-binding proteins (RBPs) are the key effectors of this process, regulating the translation and stability of parental transcripts to control cell fate specification events prior to zygotic gene activation. The KH-domain RBP MEX-3 is conserved from nematode to human. It was first discovered in Caenorhabditis elegans, where it is essential for anterior cell fate and embryo viability. Here, we show that loss of the endogenous mex-3 3´UTR disrupts its germline expression pattern. An allelic series of 3´UTR deletion variants identify repressing regions of the UTR and demonstrate that repression is not precisely coupled to reproductive success. We also show that several RBPs regulate mex-3 mRNA through its 3´UTR to define its unique germline spatiotemporal expression pattern. Additionally, we find that both poly(A) tail length control and the translation initiation factor IFE-3 contribute to its expression pattern. Together, our results establish the importance of the mex-3 3´UTR to reproductive health and its expression in the germline. Our results suggest that additional mechanisms control MEX-3 function when 3´UTR regulation is compromised.

In sexually reproducing organisms, germ cells undergo meiosis and differentiate to form oocytes or sperm. Coordination of this process requires a gene regulatory program that acts while the genome is undergoing chromatin condensation. As such, RNA regulatory pathways are an important contributor. The germline of the nematode Caenorhabditis elegans is a suitable model system to study germ cell differentiation. Several RNA-binding proteins (RBPs) coordinate each transition in the germline such as the transition from mitosis to meiosis. MEX-3 is a conserved RNA-binding protein found in most animals including humans. In C. elegans, MEX-3 displays a highly restricted pattern of expression. Here, we define the importance of the 3´UTR in regulating MEX-3 expression pattern in vivo and characterize the RNA-binding proteins involved in this regulation. Our results show that deleting various mex-3 3´UTR regions alter the pattern of expression in the germline in various ways. These mutations also reduced—but did not eliminate—reproductive capacity. Finally, we demonstrate that multiple post-transcriptional mechanisms control MEX-3 levels in different domains of the germline. Our data suggest that coordination of MEX-3 activity requires multiple layers of regulation to ensure reproductive robustness.

CTCF as a regulator of alternative splicing: new tricks for an old player

Three decades of research have established the CCCTC-binding factor (CTCF) as a ubiquitously expressed chromatin organizing factor and master regulator of gene expression. A new role for CTCF as a regulator of alternative splicing (AS) has now emerged. CTCF has been directly and indirectly linked to the modulation of AS at the individual transcript and at the transcriptome-wide level. The emerging role of CTCF-mediated regulation of AS involves diverse mechanisms; including transcriptional elongation, DNA methylation, chromatin architecture, histone modifications, and regulation of splicing factor expression and assembly. CTCF thereby appears to not only co-ordinate gene expression regulation but contributes to the modulation of transcriptomic complexity. In this review, we highlight previous discoveries regarding the role of CTCF in AS. In addition, we summarize detailed mechanisms by which CTCF mediates AS regulation. We propose opportunities for further research designed to examine the possible fate of CTCF-mediated alternatively spliced genes and associated biological consequences. CTCF has been widely acknowledged as the 'master weaver of the genome'. Given its multiple connections, further characterization of CTCF's emerging role in splicing regulation might extend its functional repertoire towards a 'conductor of the splicing orchestra'.

The Detection and Bioinformatic Analysis of Alternative 3' UTR Isoforms as Potential Cancer Biomarkers

Alternative transcript cleavage and polyadenylation is linked to cancer cell transformation, proliferation and outcome. This has led researchers to develop methods to detect and bioinformatically analyse alternative polyadenylation as potential cancer biomarkers. If incorporated into standard prognostic measures such as gene expression and clinical parameters, these could advance cancer prognostic testing and possibly guide therapy. In this review, we focus on the existing methodologies, both experimental and computational, that have been applied to support the use of alternative polyadenylation as cancer biomarkers.

Single cell RNA sequencing for breast cancer: present and future

Breast cancer is one of the most common malignant tumors in women. It is a heterogeneous disease related to genetic and environmental factors. Presently, the treatment of breast cancer still faces challenges due to recurrence and metastasis. The emergence of single-cell RNA sequencing (scRNA-seq) technology has brought new strategies to deeply understand the biological behaviors of breast cancer. By analyzing cell phenotypes and transcriptome differences at the single-cell level, scRNA-seq reveals the heterogeneity, dynamic growth and differentiation process of cells. This review summarizes the application of scRNA-seq technology in breast cancer research, such as in studies on cell heterogeneity, cancer cell metastasis, drug resistance, and prognosis. scRNA-seq technology is of great significance to deeply analyze the mechanism of breast cancer occurrence and development, identify new therapeutic targets and develop new therapeutic approaches for breast cancer.

Alternative Polyadenylation in Stem Cell Self-Renewal and Differentiation

Cellular function is shaped by transcriptional and post-transcriptional mechanisms, including alternative polyadenylation (APA). By directly controlling 3'- untranslated region (UTR) length and the selection of the last exon, APA regulates up to 70% of all cellular transcripts influencing RNA stability, output, and protein isoform expression. Cell-state-dependent 3'-UTR shortening has been identified as a hallmark of cellular proliferation. Hence, quiescent/dormant stem cells are characterized by long 3'-UTRs, whereas proliferative stem/progenitor cells exhibit 3'-UTR shortening. Here, the latest studies analyzing the role of APA in regulating stem cell state, self-renewal, differentiation, and metabolism are reviewed. The new role of APA in controlling stem cell fate opens novel potential therapeutic avenues in the field of regenerative medicine.

Influence of Age on Skeletal Muscle Hypertrophy and Atrophy Signaling: Established Paradigms and Unexpected Links

Skeletal muscle atrophy in an inevitable occurrence with advancing age, and a consequence of disease including cancer. Muscle atrophy in the elderly is managed by a regimen of resistance exercise and increased protein intake. Understanding the signaling that regulates muscle mass may identify potential therapeutic targets for the prevention and reversal of muscle atrophy in metabolic and neuromuscular diseases. This review covers the major anabolic and catabolic pathways that regulate skeletal muscle mass, with a focus on recent progress and potential new players.

Systematic refinement of gene annotations by parsing mRNA 3' end sequencing datasets

Alternative cleavage and polyadenylation generates mRNA 3' isoforms in a cell type-specific manner. Due to finite available RNA sequencing data of organisms with vast cell type complexity, currently available gene annotation resources are incomplete, which poses significant challenges to the comprehensive interpretation and quantification of transcriptomes. In this chapter, we introduce 3'GAmES, a stand-alone computational pipeline for the identification and quantification of novel mRNA 3'end isoforms from 3'mRNA sequencing data. 3'GAmES expands available repositories and improves comprehensive gene-tag counting by cost-effective 3' mRNA sequencing, faithfully mirroring whole-transcriptome RNAseq measurements. By employing R and bash shell scripts (assembled in a Singularity container) 3'GAmES systematically augments cell type-specific 3' ends of RNA polymerase II transcripts and increases the sensitivity of quantitative gene expression profiling by 3' mRNA sequencing. Public access: https://github.com/AmeresLab/3-GAmES.git.

PAPOLA contributes to cyclin D1 mRNA alternative polyadenylation and promotes breast cancer cell proliferation

Poly(A) polymerases add the poly(A) tail at the 3' end of nearly all eukaryotic mRNA, and are associated with proliferation and cancer. To elucidate the role of the most-studied mammalian poly(A) polymerase, poly(A) polymerase α (PAPOLA), in cancer, we assessed its expression in 221 breast cancer samples and found it to correlate strongly with the aggressive triple-negative subtype. Silencing PAPOLA in MCF-7 and MDA-MB-231 breast cancer cells reduced proliferation and anchorage-independent growth by decreasing steady-state cyclin D1 (CCND1) mRNA and protein levels. Whereas the length of the CCND1 mRNA poly(A) tail was not affected, its 3' untranslated region (3'UTR) lengthened. Overexpressing PAPOLA caused CCND1 mRNA 3'UTR shortening with a concomitant increase in the amount of corresponding transcript and protein, resulting in growth arrest in MCF-7 cells and DNA damage in HEK-293 cells. Such overexpression of PAPOLA promoted proliferation in the p53 mutant MDA-MB-231 cells. Our data suggest that PAPOLA is a possible candidate target for the control of tumor growth that is mostly relevant to triple-negative tumors, a group characterized by PAPOLA overexpression and lack of alternative targeted therapies.

On the function and relevance of alternative 3'-UTRs in gene expression regulation

Messanger RNA (mRNA) isoforms with alternative 3'-untranslated regions (3'-UTRs) are produced by alternative polyadenylation (APA), which occurs during transcription in most eukaryotic genes. APA fine-tunes gene expression in a cell-type- and cellular state-dependent manner. Selection of an APA site entails the binding of core cleavage and polyadenylation factors to a particular polyadenylation site localized in the pre-mRNA and is controlled by multiple regulatory determinants, including transcription, pre-mRNA cis-regulatory sequences, and protein factors. Alternative 3'-UTRs serve as platforms for specific RNA binding proteins and microRNAs, which regulate gene expression in a coordinated manner by controlling mRNA fate and function in the cell. Genome-wide studies illustrated the full extent of APA prevalence and revealed that specific 3'-UTR profiles are associated with particular cellular states and diseases. Generally, short 3'-UTRs are associated with proliferative and cancer cells, and long 3'-UTRs are mostly found in polarized and differentiated cells. Fundamental new insights on the physiological consequences of this widespread event and the molecular mechanisms involved have been revealed through single-cell studies. Publicly available comprehensive databases that cover all APA mRNA isoforms identified in many cellular states and diseases reveal specific APA signatures. Therapies tackling APA mRNA isoforms or APA regulators may be regarded as innovative and attractive tools for diagnostics or treatment of several pathologies. We highlight the function of APA and alternative 3'-UTRs in gene expression regulation, the control of these mechanisms, their physiological consequences, and their potential use as new biomarkers and therapeutic tools. This article is categorized under: RNA Processing > 3' End Processing RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA in Disease and Development > RNA in Disease.

Revisiting poly(A)-binding proteins: Multifaceted regulators during gametogenesis and early embryogenesis

Post-transcriptional regulation faces a distinctive challenge in gametes. Transcription is limited when the germ cells enter the division phase due to condensed chromatin, while gene expression during gamete maturation, fertilization, and early cleavage depends on existing mRNA post-transcriptional coordination. The dynamics of the 3'-poly(A) tail play crucial roles in defining mRNA fate. The 3'-poly(A) tail is covered with poly(A)-binding proteins (PABPs) that help to mediate mRNA metabolism and recent work has shed light on the number and function of germ cell-specific expressed PABPs. There are two structurally different PABP groups distinguished by their cytoplasmic and nuclear localization. Both lack catalytic activity but are coupled with various roles through their interaction with multifunctional partners during mRNA metabolism. Here, we present a synopsis of PABP function during gametogenesis and early embryogenesis and describe both conventional and current models of the functions and regulation of PABPs, with an emphasis on the physiological significance of how germ cell-specific PABPs potentially affect human fertility.

De novo full length transcriptome analysis of Arachis glabrata provides insights into gene expression dynamics in response to biotic and abiotic stresses

The perennial ornamental peanut Arachis glabrata represents one of the most adaptable wild Arachis species. This study used PacBio combined with BGISEQ-500 RNA-seq technology to study the transcriptome and gene expression dynamics of A. glabrata. Of the total 109,747 unique transcripts obtained, >90,566 transcripts showed significant homology to known proteins and contained the complete coding sequence (CDS). RNA-seq revealed that 1229, 1039, 1671, 3923, 1521 and 1799 transcripts expressed specifically in the root, stem, leaf, flower, peg and pod, respectively. We also identified thousands of differentially expressed transcripts in response to drought, salt, cold and leaf spot disease. Furthermore, we identified 30 polyphenol oxidase encoding genes associated with the quality of forage, making A. glabrata suitable as a forage crop. Our findings presented the first transcriptome study of A. glabrata which will facilitate genetic and genomics studies and lays the groundwork for a deeper understanding of the A. glabrata genome.

POINT technology illuminates the processing of polymerase-associated intact nascent transcripts

Mammalian chromatin is the site of both RNA polymerase II (Pol II) transcription and coupled RNA processing. However, molecular details of such co-transcriptional mechanisms remain obscure, partly because of technical limitations in purifying authentic nascent transcripts. We present a new approach to characterize nascent RNA, called polymerase intact nascent transcript (POINT) technology. This three-pronged methodology maps nascent RNA 5′ ends (POINT-5), establishes the kinetics of co-transcriptional splicing patterns (POINT-nano), and profiles whole transcription units (POINT-seq). In particular, we show by depletion of the nuclear exonuclease Xrn2 that this activity acts selectively on cleaved 5′ P-RNA at polyadenylation sites. Furthermore, POINT-nano reveals that co-transcriptional splicing either occurs immediately after splice site transcription or is delayed until Pol II transcribes downstream sequences. Finally, we connect RNA cleavage and splicing with either premature or full-length transcript termination. We anticipate that POINT technology will afford full dissection of the complexity of co-transcriptional RNA processing.

•

POINT methodology dissects intact nascent RNA processing

•

Specificity of Xrn2 exonuclease in co-transcriptional RNA degradation

•

Splicing suppresses Xrn2-dependent premature termination

•

Different kinetic classes of co-transcriptional splicing in human genes

Aptardi predicts polyadenylation sites in sample-specific transcriptomes using high-throughput RNA sequencing and DNA sequence

Annotation of polyadenylation sites from short-read RNA sequencing alone is a challenging computational task. Other algorithms rooted in DNA sequence predict potential polyadenylation sites; however, in vivo expression of a particular site varies based on a myriad of conditions. Here, we introduce aptardi (alternative polyadenylation transcriptome analysis from RNA-Seq data and DNA sequence information), which leverages both DNA sequence and RNA sequencing in a machine learning paradigm to predict expressed polyadenylation sites. Specifically, as input aptardi takes DNA nucleotide sequence, genome-aligned RNA-Seq data, and an initial transcriptome. The program evaluates these initial transcripts to identify expressed polyadenylation sites in the biological sample and refines transcript 3'-ends accordingly. The average precision of the aptardi model is twice that of a standard transcriptome assembler. In particular, the recall of the aptardi model (the proportion of true polyadenylation sites detected by the algorithm) is improved by over three-fold. Also, the model-trained using the Human Brain Reference RNA commercial standard-performs well when applied to RNA-sequencing samples from different tissues and different mammalian species. Finally, aptardi's input is simple to compile and its output is easily amenable to downstream analyses such as quantitation and differential expression.

Cleavage and Polyadenylation Specific Factor 1 Promotes Tumor Progression via Alternative Polyadenylation and Splicing in Hepatocellular Carcinoma

Alternative polyadenylation (APA) is an important post-transcriptional regulatory mechanism required for cleavage and polyadenylation (CPA) of the 3′ untranslated region (3′ UTR) of mRNAs. Several aberrant APA events have been reported in hepatocellular carcinoma (HCC). However, the regulatory mechanisms underlying APA remain unclear. In this study, we found that the expression of cleavage and polyadenylation specific factor 1 (CPSF1), a major component of the CPA complex, was significantly increased in HCC tissues and correlated with unfavorable survival outcomes. Knockdown of CPSF1 inhibited HCC cell proliferation and migration, whereas overexpression of CPSF1 caused the opposite effect. Based on integrative analysis of Iso-Seq and RNA-seq data from HepG2.2.15 cells, we identified a series of transcripts with differential 3′ UTR lengths following the knockdown of CPSF1. These transcripts were related to the biological functions of gene transcription, cytoskeleton maintenance, and endomembrane system transportation. Moreover, knockdown of CPSF1 induced an increase in alternative splicing (AS) events in addition to APA. Taken together, this study provides new insights into our understanding of the post-transcriptional regulatory mechanisms in HCC and implies that CPSF1 may be a potential prognostic biomarker and therapeutic target for HCC.

Mechanism and consequences of herpes simplex virus 1-mediated regulation of host mRNA alternative polyadenylation

Eukaryotic gene expression is extensively regulated by cellular stress and pathogen infections. We have previously shown that herpes simplex virus 1 (HSV-1) and several cellular stresses cause widespread disruption of transcription termination (DoTT) of RNA polymerase II (RNAPII) in host genes and that the viral immediate early factor ICP27 plays an important role in HSV-1-induced DoTT. Here, we show that HSV-1 infection also leads to widespread changes in alternative polyadenylation (APA) of host mRNAs. In the majority of cases, polyadenylation shifts to upstream poly(A) sites (PAS), including many intronic PAS. Mechanistically, ICP27 contributes to HSV-1-mediated APA regulation. HSV-1- and ICP27-induced activation of intronic PAS is sequence-dependent and does not involve general inhibition of U1 snRNP. HSV1-induced intronic polyadenylation is accompanied by early termination of RNAPII. HSV-1-induced mRNAs polyadenylated at intronic PAS (IPA) are exported into the cytoplasm while APA isoforms with extended 3' UTRs are sequestered in the nuclei, both preventing the expression of the full-length gene products. Finally we provide evidence that HSV-induced IPA isoforms are translated. Together with other recent studies, our results suggest that viral infection and cellular stresses induce a multi-faceted host response that includes DoTT and changes in APA profiles.

Divergence in alternative polyadenylation contributes to gene regulatory differences between humans and chimpanzees

While comparative functional genomic studies have shown that inter-species differences in gene expression can be explained by corresponding inter-species differences in genetic and epigenetic regulatory mechanisms, co-transcriptional mechanisms, such as alternative polyadenylation (APA), have received little attention. We characterized APA in lymphoblastoid cell lines from six humans and six chimpanzees by identifying and estimating the usage for 44,432 polyadenylation sites (PAS) in 9518 genes. Although APA is largely conserved, 1705 genes showed significantly different PAS usage (FDR 0.05) between species. Genes with divergent APA also tend to be differentially expressed, are enriched among genes showing differences in protein translation, and can explain a subset of observed inter-species protein expression differences that do not differ at the transcript level. Finally, we found that genes with a dominant PAS, which is used more often than other PAS, are particularly enriched for differentially expressed genes.

Cross-talk of four types of RNA modification writers defines tumor microenvironment and pharmacogenomic landscape in colorectal cancer

BACKGROUND

The four major RNA adenosine modifications, i.e., m6A, m1A, alternative polyadenylation, and adenosine-to-inosine RNA editing, are mediated mostly by the "writer" enzymes and constitute critical mechanisms of epigenetic regulation in immune response and tumorigenesis. However, the cross-talk and potential roles of these "writers" in the tumor microenvironment (TME), drug sensitivity, and immunotherapy remain unknown.

METHODS

We systematically characterized mRNA expression and genetic alterations of 26 RNA modification "writers" in colorectal cancer (CRC), and evaluated their expression pattern in 1697 CRC samples from 8 datasets. We used an unsupervised clustering method to assign the samples into two patterns of expression of RNA modification "writers". Subsequently, we constructed the RNA modification "writer" Score (WM_Score) model based on differentially expressed genes (DEGs) responsible for the RNA modification patterns to quantify the RNA modification-related subtypes of individual tumors. Furthermore, we performed association analysis for WM_Score and characteristics of TME, consensus molecular subtypes (CMSs), clinical features, transcriptional and post-transcriptional regulation, drug response, and the efficacy of immunotherapy.

RESULTS

We demonstrated that multi-layer alterations of RNA modification "writer" are associated with patient survival and TME cell-infiltrating characteristics. We identified two distinct RNA modification patterns, characterized by a high and a low WM_Score. The WM_Score-high group was associated with worse patient overall survival and with the infiltration of inhibitory immune cells, such as M2 macrophages, EMT activation, and metastasis, while the WM_Score-low group was associated with a survival advantage, apoptosis, and cell cycle signaling pathways. WM_Score correlated highly with the regulation of transcription and post-transcriptional events contributing to the development of CRC. In response to anti-cancer drugs, WM_Score highly negatively correlated (drug sensitive) with drugs which targeted oncogenic related pathways, such as MAPK, EGFR, and mTOR signaling pathways, positively correlated (drug resistance) with drugs which targeted in apoptosis and cell cycle. Importantly, the WM_Score was associated with the therapeutic efficacy of PD-L1 blockade, suggesting that the development of potential drugs targeting these "writers" to aid the clinical benefits of immunotherapy.

CONCLUSIONS

Our study is the first to provide a comprehensive analysis of four RNA modifications in CRC. We revealed the potential function of these writers in TME, transcriptional and post-transcriptional events, and identified their therapeutic liability in targeted therapy and immunotherapy. This work highlights the cross-talk and potential clinical utility of RNA modification "writers" in cancer therapy.

LARP1 and LARP4: up close with PABP for mRNA 3' poly(A) protection and stabilization

La-related proteins (LARPs) share a La motif (LaM) followed by an RNA recognition motif (RRM). Together these are termed the La-module that, in the prototypical nuclear La protein and LARP7, mediates binding to the UUU-3'OH termination motif of nascent RNA polymerase III transcripts. We briefly review La and LARP7 activities for RNA 3' end binding and protection from exonucleases before moving to the more recently uncovered poly(A)-related activities of LARP1 and LARP4. Two features shared by LARP1 and LARP4 are direct binding to poly(A) and to the cytoplasmic poly(A)-binding protein (PABP, also known as PABPC1). LARP1, LARP4 and other proteins involved in mRNA translation, deadenylation, and decay, contain PAM2 motifs with variable affinities for the MLLE domain of PABP. We discuss a model in which these PABP-interacting activities contribute to poly(A) pruning of active mRNPs. Evidence that the SARS-CoV-2 RNA virus targets PABP, LARP1, LARP 4 and LARP 4B to control mRNP activity is also briefly reviewed. Recent data suggests that LARP4 opposes deadenylation by stabilizing PABP on mRNA poly(A) tails. Other data suggest that LARP1 can protect mRNA from deadenylation. This is dependent on a PAM2 motif with unique characteristics present in its La-module. Thus, while nuclear La and LARP7 stabilize small RNAs with 3' oligo(U) from decay, LARP1 and LARP4 bind and protect mRNA 3' poly(A) tails from deadenylases through close contact with PABP.Abbreviations: 5'TOP: 5' terminal oligopyrimidine, LaM: La motif, LARP: La-related protein, LARP1: La-related protein 1, MLLE: mademoiselle, NTR: N-terminal region, PABP: cytoplasmic poly(A)-binding protein (PABPC1), Pol III: RNA polymerase III, PAM2: PABP-interacting motif 2, PB: processing body, RRM: RNA recognition motif, SG: stress granule.

Fine gene expression regulation by minor sequence variations downstream of the polyadenylation signal

The termination of transcription is a complex process that substantially contributes to gene regulation in eukaryotes. Previously, it was noted that a single cytosine deletion at the position + 32 bp relative to the single polyadenylation signal AAUAAA (hereafter the dC mutation) causes a 2-fold increase in the transcription level of the upstream eGFP reporter in mouse embryonic stem cells. Here, we analyzed the conservation of this phenomenon in immortalized mouse, human and drosophila cell lines and the influence of the dC mutation on the choice of the pre-mRNA cleavage sites. We have constructed dual-reporter plasmids to accurately measure the effect of the dC and other nearby located mutations on eGFP mRNA level by RT-qPCR. In this way, we found that the dC mutation leads to a 2-fold increase in the expression level of the upstream eGFP reporter gene in cultured mouse and human, but not in drosophila cells. In addition, 3' RACE analysis demonstrated that eGFP pre-mRNAs are cut at multiple positions between + 14 to + 31, and that the most proximal cleavage site becomes almost exclusively utilized in the presence of the dC mutation. We also identified new short sequence variations located within positions + 25.. + 40 and + 33.. + 48 that increase eGFP expression up to ~2-4-fold. Altogether, the positive effect of the dC mutation seems to be conserved in mouse embryonic stem cells, mouse embryonic 3T3 fibroblasts and human HEK293T cells. In the latter cells, the dC mutation appears to be involved in regulating pre-mRNA cleavage site selection. Finally, a multiplexed approach is proposed to identify motifs located downstream of cleavage site(s) that are essential for transcription termination.

Animal-APAdb: a comprehensive animal alternative polyadenylation database

Alternative polyadenylation (APA) is an important post-transcriptional regulatory mechanism that recognizes different polyadenylation signals on transcripts, resulting in transcripts with different lengths of 3′ untranslated regions and thereby influencing a series of biological processes. Recent studies have highlighted the important roles of APA in human. However, APA profiles in other animals have not been fully recognized, and there is no database that provides comprehensive APA information for other animals except human. Here, by using the RNA sequencing data collected from public databases, we systematically characterized the APA profiles in 9244 samples of 18 species. In total, we identified 342 952 APA events with a median of 17 020 per species using the DaPars2 algorithm, and 315 691 APA events with a median of 17 953 per species using the QAPA algorithm in these 18 species, respectively. In addition, we predicted the polyadenylation sites (PAS) and motifs near PAS of these species. We further developed Animal-APAdb, a user-friendly database (http://gong_lab.hzau.edu.cn/Animal-APAdb/) for data searching, browsing and downloading. With comprehensive information of APA events in different tissues of different species, Animal-APAdb may greatly facilitate the exploration of animal APA patterns and novel mechanisms, gene expression regulation and APA evolution across tissues and species.

NUDT21 Links Mitochondrial IPS-1 to RLR-Containing Stress Granules and Activates Host Antiviral Defense

Viral RNA in the cytoplasm of mammalian host cells is recognized by retinoic acid-inducible protein-I-like receptors (RLRs), which localize to cytoplasmic stress granules (SGs). Activated RLRs associate with the mitochondrial adaptor protein IPS-1, which activates antiviral host defense mechanisms, including type I IFN induction. It has remained unclear, however, how RLRs in SGs and IPS-1 in the mitochondrial outer membrane associate physically and engage in information transfer. In this study, we show that NUDT21, an RNA-binding protein that regulates alternative transcript polyadenylation, physically associates with IPS-1 and mediates its localization to SGs in response to transfection with polyinosinic-polycytidylic acid [poly(I:C)], a mimic of viral dsRNA. We found that despite its well-established function in the nucleus, a fraction of NUDT21 localizes to mitochondria in resting cells and becomes localized to SGs in response to poly(I:C) transfection. NUDT21 was also found to be required for efficient type I IFN induction in response to viral infection in both human HeLa cells and mouse macrophage cell line RAW264.7 cells. Our results together indicate that NUDT21 links RLRs in SGs to mitochondrial IPS-1 and thereby activates host defense responses to viral infection.

The Role of cis- and trans-Acting RNA Regulatory Elements in Leukemia

RNA molecules are a source of phenotypic diversity and an operating system that connects multiple genetic and metabolic processes in the cell. A dysregulated RNA network is a common feature of cancer. Aberrant expression of long non-coding RNA (lncRNA), micro RNA (miRNA), and circular RNA (circRNA) in tumors compared to their normal counterparts, as well as the recurrent mutations in functional regulatory cis-acting RNA motifs have emerged as biomarkers of disease development and progression, opening avenues for the design of novel therapeutic approaches. This review looks at the progress, challenges and future prospects of targeting cis-acting and trans-acting RNA elements for leukemia diagnosis and treatment.

Profiling Novel Alternative Splicing within Multiple Tissues Provides Useful Insights into Porcine Genome Annotation

Alternative splicing (AS) is a process during gene expression that results in a single gene coding for different protein variants. AS contributes to transcriptome and proteome diversity. In order to characterize AS in pigs, genome-wide transcripts and AS events were detected using RNA sequencing of 34 different tissues in Duroc pigs. In total, 138,403 AS events and 29,270 expressed genes were identified. An alternative donor site was the most common AS form and accounted for 44% of the total AS events. The percentage of the other three AS forms (exon skipping, alternative acceptor site, and intron retention) was approximately 19%. The results showed that the most common AS events involving alternative donor sites could produce different transcripts or proteins that affect the biological processes. The expression of genes with tissue-specific AS events showed that gene functions were consistent with tissue functions. AS increased proteome diversity and resulted in novel proteins that gained or lost important functional domains. In summary, these findings extend porcine genome annotation and highlight roles that AS could play in determining tissue identity.

Alternative Polyadenylation: a new frontier in post transcriptional regulation

Polyadenylation of pre-messenger RNA (pre-mRNA) specific sites and termination of their downstream transcriptions are signaled by unique sequence motif structures such as AAUAAA and its auxiliary elements. Alternative polyadenylation (APA) is an important post-transcriptional regulatory mechanism that processes RNA products depending on its 3'-untranslated region (3'-UTR) specific sequence signal. APA processing can generate several mRNA isoforms from a single gene, which may have different biological functions on their target gene. As a result, cellular genomic stability, proliferation capability, and transformation feasibility could all be affected. Furthermore, APA modulation regulates disease initiation and progression. APA status could potentially act as a biomarker for disease diagnosis, severity stratification, and prognosis forecast. While the advance of modern throughout technologies, such as next generation-sequencing (NGS) and single-cell sequencing techniques, have enriched our knowledge about APA, much of APA biological process is unknown and pending for further investigation. Herein, we review the current knowledge on APA and how its regulatory complex factors (CFI/IIm, CPSF, CSTF, and RBPs) work together to determine RNA splicing location, cell cycle velocity, microRNA processing, and oncogenesis regulation. We also discuss various APA experiment strategies and the future direction of APA research.

Aptamer-Mediated Control of Polyadenylation for Gene Expression Regulation in Mammalian Cells

Small aptamer-based regulatory devices can be designed to control a range of RNA-dependent cellular processes and emerged as promising tools for fine-tuning gene expression in synthetic biology. Here, we design a conceptually new riboswitch device that allows for the conditional regulation of polyadenylation. By making use of ligand-induced sequence occlusion, the system efficiently controls the accessibility of the eukaryotic polyadenylation signal. Undesirable 3'-extended read-through products are counteracted by the downstream insertion of a microRNA target site. We demonstrate the modularity of the system with regard to sensor aptamers and polyadenylation signals used and combine the newly designed riboswitch with well-known aptazymes to yield superior composite systems. In addition, we show that the switches can be used to control alternative polyadenylation. The presented genetic switches require very little coding space and can be easily optimized by rational adjustments of the thermodynamic stability. The polyadenylation riboswitch extends the repertoire of RNA-based regulators and opens new possibilities for the generation of complex synthetic circuits.

Poly(A)-DG: A deep-learning-based domain generalization method to identify cross-species Poly(A) signal without prior knowledge from target species

In eukaryotes, polyadenylation (poly(A)) is an essential process during mRNA maturation. Identifying the cis-determinants of poly(A) signal (PAS) on the DNA sequence is the key to understand the mechanism of translation regulation and mRNA metabolism. Although machine learning methods were widely used in computationally identifying PAS, the need for tremendous amounts of annotation data hinder applications of existing methods in species without experimental data on PAS. Therefore, cross-species PAS identification, which enables the possibility to predict PAS from untrained species, naturally becomes a promising direction. In our works, we propose a novel deep learning method named Poly(A)-DG for cross-species PAS identification. Poly(A)-DG consists of a Convolution Neural Network-Multilayer Perceptron (CNN-MLP) network and a domain generalization technique. It learns PAS patterns from the training species and identifies PAS in target species without re-training. To test our method, we use four species and build cross-species training sets with two of them and evaluate the performance of the remaining ones. Moreover, we test our method against insufficient data and imbalanced data issues and demonstrate that Poly(A)-DG not only outperforms state-of-the-art methods but also maintains relatively high accuracy when it comes to a smaller or imbalanced training set.

Full length transcriptome profiling reveals novel immune-related genes in black rockfish (Sebastes schlegelii)

Lacking full-length transcriptome for black rockfish (Sebastes schlegelii) limits novel gene discoveries and gene structures analysis. Therefore, we constructed the full-length transcriptome of black rockfish using Single-Molecule Real-Time Sequencing technology. Totally, we produced 21.73 Gb raw reads containing 298,904 circular consensus sequence (CCS) reads. Full-length (FL) and Non-full-length (NFL) isoforms were obtained based on the presence of 5' and 3' primers as well as poly (A) tails. The results showed 70.71% reads were identified as FL isoforms. Moreover, the average length of these PacBio isoforms is 2,632 bp, which is much longer than the length of the unigenes with the average length of 589 bp which generated from Illumina platform. Meanwhile, we identified 43,068 non-redundant transcripts, 12,485 alternative splicing (AS), 6,320 polyadenylation (APA) and 499 gene fusions as well as numerous long non-coding RNAs based on mapped FL isoforms. In addition, we identified 147 and 528 immune-related genes from novel genes and unmapped transcripts. The provided dataset can be utilized to discover novel genes and construct a comprehensive transcript dataset for black rockfish.

Single-cell RNA cap and tail sequencing (scRCAT-seq) reveals subtype-specific isoforms differing in transcript demarcation

The differences in transcription start sites (TSS) and transcription end sites (TES) among gene isoforms can affect the stability, localization, and translation efficiency of mRNA. Gene isoforms allow a single gene diverse functions across different cell types, and isoform dynamics allow different functions over time. However, methods to efficiently identify and quantify RNA isoforms genome-wide in single cells are still lacking. Here, we introduce single cell RNA Cap And Tail sequencing (scRCAT-seq), a method to demarcate the boundaries of isoforms based on short-read sequencing, with higher efficiency and lower cost than existing long-read sequencing methods. In conjunction with machine learning algorithms, scRCAT-seq demarcates RNA transcripts with unprecedented accuracy. We identified hundreds of previously uncharacterized transcripts and thousands of alternative transcripts for known genes, revealed cell-type specific isoforms for various cell types across different species, and generated a cell atlas of isoform dynamics during the development of retinal cones.

Nudt21-mediated alternative polyadenylation of HMGA2 3'-UTR impairs stemness of human tendon stem cell

Tendon-derived stem cells (TSCs) play a primary role in tendon physiology, pathology, as well as tendon repair and regeneration after injury. TSCs are often exposed to mechanical loading-related cellular stresses such as oxidative stress, resulting in loss of stemness and multipotent differentiation potential. Cytoprotective autophagy has previously been identified as an important mechanism to protect human TSCs (hTSCs) from oxidative stress induced impairments. In this study, we found that high-mobility AT-hook 2 (HMGA2) overexpression protects hTSCs against H₂O₂-induced loss of stemness through autophagy activation. Evidentially, H₂O₂ treatment increases the expression of Nudt21, a protein critical to polyadenylation site selection in alternative polyadenylation (APA) of mRNA transcripts. This leads to increased cleavage and polyadenylation of HMGA2 3'-UTR at the distal site, resulting in increased HMGA2 silencing by the microRNA let-7 and reduced HMGA2 expression. In conclusion, Nudt21-regulated APA of HMGA2 3'-UTR and subsequent HMGA2 downregulation mediates oxidative stress induced hTSC impairments.

Downregulation of CFIm25 amplifies dermal fibrosis through alternative polyadenylation

This study implicates the key regulator of alternative polyadenylation, CFIm25 in dermal fibrosis and in systemic sclerosis (scleroderma) pathogenesis. CFIm25 downregulation promotes the expression of profibrotic factors, exaggerates bleomycin-induced skin fibrosis, while CFIm25 restoration attenuates skin fibrosis.

Systemic sclerosis (SSc; scleroderma) is a multisystem fibrotic disease. The mammalian cleavage factor I 25-kD subunit (CFIm25; encoded by NUDT21) is a key regulator of alternative polyadenylation, and its depletion causes predominantly 3′UTR shortening through loss of stimulation of distal polyadenylation sites. A shortened 3′UTR will often lack microRNA target sites, resulting in increased mRNA translation due to evasion of microRNA-mediated repression. Herein, we report that CFlm25 is downregulated in SSc skin, primary dermal fibroblasts, and two murine models of dermal fibrosis. Knockdown of CFIm25 in normal skin fibroblasts is sufficient to promote the 3′UTR shortening of key TGFβ-regulated fibrotic genes and enhance their protein expression. Moreover, several of these fibrotic transcripts show 3′UTR shortening in SSc skin. Finally, mice with CFIm25 deletion in fibroblasts show exaggerated skin fibrosis upon bleomycin treatment, and CFIm25 restoration attenuates bleomycin-induced skin fibrosis. Overall, our data link this novel RNA-processing mechanism to dermal fibrosis and SSc pathogenesis.

Role of Arginine Methylation in Alternative Polyadenylation of VEGFR-1 (Flt-1) pre-mRNA

Mature mRNA is generated by the 3ʹ end cleavage and polyadenylation of its precursor pre-mRNA. Eukaryotic genes frequently have multiple polyadenylation sites, resulting in mRNA isoforms with different 3ʹ-UTR lengths that often encode different C-terminal amino acid sequences. It is well-known that this form of post-transcriptional modification, termed alternative polyadenylation, can affect mRNA stability, localization, translation, and nuclear export. We focus on the alternative polyadenylation of pre-mRNA for vascular endothelial growth factor receptor-1 (VEGFR-1), the receptor for VEGF. VEGFR-1 is a transmembrane protein with a tyrosine kinase in the intracellular region. Secreted forms of VEGFR-1 (sVEGFR-1) are also produced from the same gene by alternative polyadenylation, and sVEGFR-1 has a function opposite to that of VEGFR-1 because it acts as a decoy receptor for VEGF. However, the mechanism that regulates the production of sVEGFR-1 by alternative polyadenylation remains poorly understood. In this review, we introduce and discuss the mechanism of alternative polyadenylation of VEGFR-1 mediated by protein arginine methylation.

So close, no matter how far: multiple paths connecting transcription to mRNA translation in eukaryotes

Transcription of DNA into mRNA and translation of mRNA into proteins are two major processes underlying gene expression. Due to the distinct molecular mechanisms, timings, and locales of action, these processes are mainly considered to be independent. During the last two decades, however, multiple factors and elements were shown to coordinate transcription and translation, suggesting an intricate level of synchronization. This review discusses the molecular mechanisms that impact both processes in eukaryotic cells of different origins. The emerging global picture suggests evolutionarily conserved regulation and coordination between transcription and mRNA translation, indicating the importance of this phenomenon for the fine-tuning of gene expression and the adjustment to constantly changing conditions.

Dissecting the heterogeneity of the alternative polyadenylation profiles in triple-negative breast cancers

Background: Triple-negative breast cancer (TNBC) is an aggressive malignancy with high heterogeneity. However, the alternative polyadenylation (APA) profiles of TNBC remain unknown. Here, we aimed to define the characteristics of the APA events at post-transcription level among TNBCs. Methods: Using transcriptome microarray data, we analyzed APA profiles of 165 TNBC samples and 33 paired normal tissues. A pooled short hairpin RNA screen targeting 23 core cleavage and polyadenylation (C/P) genes was used to identify key C/P factors. Results: We established an unconventional APA subtyping system composed of four stable subtypes: 1) luminal androgen receptor (LAR), 2) mesenchymal-like immune-activated (MLIA), 3) basal-like (BL), 4) suppressed (S) subtypes. Patients in the S subtype had the worst disease-free survival comparing to other patients (log-rank p = 0.021). Enriched clinically actionable pathways and putative therapeutic APA events were analyzed among each APA subtype. Furthermore, CPSF1 and PABPN1 were identified as the master C/P factors in regulating APA events and TNBC proliferation. The depletion of CPSF1 or PABPN1 weakened cell proliferation, enhanced apoptosis, resulted in cell cycle redistribution and a reversion of APA events of genes associated with tumorigenesis, proliferation, metastasis and chemosensitivity in breast cancer. Conclusions: Our findings advance the understanding of tumor heterogeneity regulation in APA and yield new insights into therapeutic target identification in TNBC.

A genome-wide scan for candidate lethal variants in Thoroughbred horses

Domestic animal populations are often characterised by high rates of inbreeding and low effective population sizes due to selective breeding practices. These practices can result in otherwise rare recessive deleterious alleles drifting to high frequencies, resulting in reduced fertility rates. This study aimed to identify potential recessive lethal haplotypes in the Thoroughbred horse breed, a closed population that has been selectively bred for racing performance. In this study, we identified a haplotype in the LY49B gene that shows strong evidence of being homozygous lethal, despite having high frequencies of heterozygotes in Thoroughbreds and other domestic horse breeds. Variant analysis of whole-genome sequence data identified two SNPs in the 3'UTR of the LY49B gene that may result in loss of function. Analysis of transcriptomic data from equine embryonic tissue revealed that LY49B is expressed in the trophoblast during placentation stage of development. These findings suggest that LY49B may have an essential, but as yet unknown function in the implantation stage of equine development. Further investigation of this region may allow for the development of a genetic test to improve fertility rates in horse populations. Identification of other lethal variants could assist in improving natural levels of fertility in horse populations.

HNRNPA1-mediated 3' UTR length changes of HN1 contributes to cancer- and senescence-associated phenotypes

Cellular senescence has been regarded as a mechanism of tumor suppression. Studying the regulation of gene expression at various levels in cell senescence will shed light on cancer therapy. Alternative polyadenylation (APA) regulates gene expression by altering 3′ untranslated regions (3′ UTR) and plays important roles in diverse biological processes. However, whether APA of a specific gene functions in both cancer and senescence remains unclear. Here, we discovered that 3′ UTR of HN1 (or JPT1) showed shortening in cancers and lengthening in senescence, correlated well with its high expression in cancer cells and low expression in senescent cells, respectively. HN1 transcripts with longer 3′ UTR were less stable and produced less protein. Down-regulation of HN1 induced senescence-associated phenotypes in both normal and cancer cells. Patients with higher HN1 expression had lower survival rates in various carcinomas. Interestingly, down-regulating the splicing factor HNRNPA1 induced 3′ UTR lengthening of HN1 and senescence-associated phenotypes, which could be partially reversed by overexpressing HN1. Together, we revealed for the first time that HNRNPA1-mediated APA of HN1 contributed to cancer- and senescence-related phenotypes. Given senescence is a cancer prevention mechanism, our discovery indicates the HNRNPA1-HN1 axis as a potential target for cancer treatment.

A Generative Neural Network for Maximizing Fitness and Diversity of Synthetic DNA and Protein Sequences

Engineering gene and protein sequences with defined functional properties is a major goal of synthetic biology. Deep neural network models, together with gradient ascent-style optimization, show promise for sequence design. The generated sequences can however get stuck in local minima and often have low diversity. Here, we develop deep exploration networks (DENs), a class of activation-maximizing generative models, which minimize the cost of a neural network fitness predictor by gradient descent. By penalizing any two generated patterns on the basis of a similarity metric, DENs explicitly maximize sequence diversity. To avoid drifting into low-confidence regions of the predictor, we incorporate variational autoencoders to maintain the likelihood ratio of generated sequences. Using DENs, we engineered polyadenylation signals with more than 10-fold higher selection odds than the best gradient ascent-generated patterns, identified splice regulatory sequences predicted to result in highly differential splicing between cell lines, and improved on state-of-the-art results for protein design tasks.

Involved microRNAs in alternative polyadenylation intervene in breast cancer via regulation of cleavage factor "CFIm25"

Cleavage factor “CFIm25”, as a key repressor at proximal poly (A) site, negatively correlates to cell proliferation and tumorigenicity in various cancers. Hence, understanding CFIm25 mechanism of action in breast cancer would be a great benefit. To this aim four steps were designed. First, potential miRNAs that target 3′-UTR of CFIm25 mRNA, retrieved from Targetscan web server. Second, screened miRNAs were profiled in 100 breast cancer and 100 normal adjacent samples. Third, miRNAs that their expression was inversely correlated to the CFIm25, overexpressed in MDA-MB-231 cell line, and their effect on proliferation and migration monitored via MTT and wound healing assays, respectively. Fourth, interaction of miRNAs of interest with 3′-UTR of CFIm25 confirmed via luciferase assay and western blot. Our results indicate that CFIm25 considerably down-regulates in human breast cancer tissue. qRT-PCR assay, luciferase test, and western blotting confirm that CFIm25 itself could be directly regulated by oncomiRs such as miR-23, -24, -27, -135, -182 and -374. Besides, according to MTT and wound healing assays of cell lines, CFIm25 knockdown intensifies cell growth, proliferation and migration. Our results also confirm indirect impact of CFIm25 on regulation of mRNA’s 3′–UTR length, which then control corresponding miRNAs’ action. miRNAs directly control CFIm25 expression level, which then tunes expression of the oncogenes and tumor proliferation. Therefore, regulation of CFIm25 expression level via miRNAs is expected to improve treatment responses in breast cancer.

DeepPASTA: deep neural network based polyadenylation site analysis

MOTIVATION

Alternative polyadenylation (polyA) sites near the 3' end of a pre-mRNA create multiple mRNA transcripts with different 3' untranslated regions (3' UTRs). The sequence elements of a 3' UTR are essential for many biological activities such as mRNA stability, sub-cellular localization, protein translation, protein binding and translation efficiency. Moreover, numerous studies in the literature have reported the correlation between diseases and the shortening (or lengthening) of 3' UTRs. As alternative polyA sites are common in mammalian genes, several machine learning tools have been published for predicting polyA sites from sequence data. These tools either consider limited sequence features or use relatively old algorithms for polyA site prediction. Moreover, none of the previous tools consider RNA secondary structures as a feature to predict polyA sites.

RESULTS

In this paper, we propose a new deep learning model, called DeepPASTA, for predicting polyA sites from both sequence and RNA secondary structure data. The model is then extended to predict tissue-specific polyA sites. Moreover, the tool can predict the most dominant (i.e. frequently used) polyA site of a gene in a specific tissue and relative dominance when two polyA sites of the same gene are given. Our extensive experiments demonstrate that DeepPASTA signisficantly outperforms the existing tools for polyA site prediction and tissue-specific relative and absolute dominant polyA site prediction.

AVAILABILITY AND IMPLEMENTATION

https://github.com/arefeen/DeepPASTA.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Extensive Involvement of Alternative Polyadenylation in Single-Nucleus Neurons

Cleavage and polyadenylation are essential processes that can impact many aspects of mRNA fate. Most eukaryotic genes have alternative polyadenylation (APA) events. While the heterogeneity of mRNA polyadenylation isoform choice has been studied in specific tissues, less attention has been paid to the neuronal heterogeneity of APA selection at single-nucleus resolution. APA is highly controlled during development and neuronal activation, however, to what extent APA events vary in a specific neuronal cell population and the regulatory mechanisms are still unclear. In this paper, we investigated dynamic APA usage in different cell types using snRNA-seq data of 1424 human brain cells generated by single-cell 3' RNA sequencing. We found that distal APA sites are not only favored by global neuronal cells, but that their usage also varies between the principal types of neuronal cell populations (excitatory neurons and inhibitory neurons). A motif analysis and a gene functional analysis indicated the enrichment of RNA-binding protein (RBP) binding sites and neuronal functions for the set of genes with neuron-enhanced distal PAS usage. Our results revealed the extensive involvement of APA regulation in neuronal populations at the single-nucleus level, providing new insights into roles for APA in specific neuronal cell populations, as well as utility in future functional studies.

Emerging Roles of RNA 3'-end Cleavage and Polyadenylation in Pathogenesis, Diagnosis and Therapy of Human Disorders

A crucial feature of gene expression involves RNA processing to produce 3′ ends through a process termed 3′ end cleavage and polyadenylation (CPA). This ensures the nascent RNA molecule can exit the nucleus and be translated to ultimately give rise to a protein which can execute a function. Further, alternative polyadenylation (APA) can produce distinct transcript isoforms, profoundly expanding the complexity of the transcriptome. CPA is carried out by multi-component protein complexes interacting with multiple RNA motifs and is tightly coupled to transcription, other steps of RNA processing, and even epigenetic modifications. CPA and APA contribute to the maintenance of a multitude of diverse physiological processes. It is therefore not surprising that disruptions of CPA and APA can lead to devastating disorders. Here, we review potential CPA and APA mechanisms involving both loss and gain of function that can have tremendous impacts on health and disease. Ultimately we highlight the emerging diagnostic and therapeutic potential CPA and APA offer.

Aminoacyl-tRNA synthetases in cell signaling

Aminoacyl-tRNA synthetases (ARSs) are a family of essential "housekeeping" enzymes ubiquitous in the three major domains of life. ARSs uniquely connect the essential minimal units of both major oligomer classes-the 3-nucleotide codons of oligonucleotides and the amino acids of proteins. They catalyze the esterification of amino acids to the 3'-end of cognate transfer RNAs (tRNAs) bearing the correct anticodon triplet to ensure accurate transfer of information from mRNA to protein according to the genetic code. As an essential translation factor responsible for the first biochemical reaction in protein biosynthesis, ARSs control protein production by catalyzing aminoacylation, and by editing of mischarged aminoacyl-tRNAs to maintain translational fidelity. In addition to their primary enzymatic activities, many ARSs have noncanonical functions unrelated to their catalytic activity in protein synthesis. Among the ARSs with "moonlighting" activities, several, including GluProRS (or EPRS), LeuRS, LysRS, SerRS, TyrRS, and TrpRS, exhibit cell signaling-related activities that sense environmental signals, regulate gene expression, and modulate cellular functions. ARS signaling functions generally depend on catalytically-inactive, appended domains not present in ancient enzyme forms, and are activated by stimulus-dependent post-translational modification. Activation often results in cellular re-localization and gain of new interacting partners. The newly formed ARS-bearing complexes conduct a host of signal transduction functions, including immune response, mTORC1 pathway signaling, and fibrogenic and angiogenic signaling, among others. Because noncanonical functions of ARSs in signal transduction are uncoupled from canonical aminoacylation functions, function-specific inhibitors can be developed, thus providing promising opportunities and therapeutic targets for treatment of human disease.

Recent molecular insights into canonical pre-mRNA 3'-end processing

The majority of eukaryotic messenger RNA precursors (pre-mRNAs) undergo cleavage and polyadenylation at their 3' end. This canonical 3'-end processing depends on sequence elements in the pre-mRNA as well as a mega-dalton protein machinery. The cleavage site in mammalian pre-mRNAs is located between an upstream poly(A) signal, most frequently an AAUAAA hexamer, and a GU-rich downstream sequence element. This review will summarize recent advances from the studies on this canonical 3'-end processing machinery. They have revealed the molecular mechanism for the recognition of the poly(A) signal and provided the first glimpse into the overall architecture of the machinery. The studies also show that the machinery is highly dynamic conformationally, and extensive re-arrangements are necessary for its activation. Inhibitors targeting the active site of the CPSF73 nuclease of this machinery have anti-cancer, anti-inflammatory and anti-protozoal effects, indicating that CPSF73 and pre-mRNA 3'-end processing in general are attractive targets for drug discovery.

ABBREVIATIONS

APA: alternative polyadenylation; β-CASP: metallo-β-lactamase-associated CPSF Artemis SNM1/PSO2; CTD: C-terminal domain; CF: cleavage factor; CPF: cleavage and polyadenylation factor; CPSF: cleavage and polyadenylation specificity factor; CstF: cleavage stimulation factor; DSE: downstream element; HAT: half a TPR; HCC: histone pre-mRNA cleavage complex; mCF: mammalian cleavage factor; mPSF: mammalian polyadenylation specificity factor; mRNA: messenger RNA; nt: nucleotide; NTD: N-terminal domain; PAP: polyadenylate polymerase; PAS: polyadenylation signal; PIM: mPSF interaction motif; Poly(A): polyadenylation, polyadenylate; Pol II: RNA polymerase II; pre-mRNA: messenger RNA precursor; RRM: RNA recognition module, RNA recognition motif; snRNP: small nuclear ribonucleoprotein; TPR: tetratricopeptide repeat; UTR: untranslated region; ZF: zinc finger.

U1 snRNP Telescripting Roles in Transcription and Its Mechanism

Telescripting is a fundamental cotranscriptional gene regulation process that relies on U1 snRNP (U1) to suppress premature 3'-end cleavage and polyadenylation (PCPA) in RNA polymerase II (Pol II) transcripts, which is necessary for full-length transcription of thousands of protein-coding (pre-mRNAs) and long noncoding (lncRNA) genes. Like U1 role in splicing, telescripting requires U1 snRNA base-pairing with nascent transcripts. Inhibition of U1 base-pairing with U1 snRNA antisense morpholino oligonucleotide (U1 AMO) mimics widespread PCPA from cryptic polyadenylation signals (PASs) in human tissues, including PCPA in introns and last exons' 3'-untranslated regions (3' UTRs). U1 telescripting-PCPA balance changes generate diverse RNAs depending on where in a gene it occurs. Long genes are highly U1-telescripting-dependent because of PASs in introns compared to short genes. Enrichment of cell cycle control, differentiation, and developmental functions in long genes, compared to housekeeping and acute cell stress response genes in short genes, reveals a gene size-function relationship in mammalian genomes. This polarization increased in metazoan evolution by previously unexplained intron expansion, suggesting that U1 telescripting could shift global gene expression priorities. We show that that modulating U1 availability can profoundly alter cell phenotype, such as cancer cell migration and invasion, underscoring the critical role of U1 homeostasis and suggesting it as a potential target for therapies. We describe a complex of U1 with cleavage and polyadenylation factors that silences PASs in introns and 3' UTR, which gives insights into U1 telescripting mechanism and transcription elongation regulation.

Alternative polyadenylation dependent function of splicing factor SRSF3 contributes to cellular senescence

Down-regulated splicing factor SRSF3 is known to promote cellular senescence, an important biological process in preventing cancer and contributing to individual aging, via its alternative splicing dependent function in human cells. Here we discovered alternative polyadenylation (APA) dependent function of SRSF3 as a novel mechanism explaining SRSF3 downregulation induced cellular senescence. Knockdown of SRSF3 resulted in preference usage of proximal poly(A) sites and thus global shortening of 3′ untranslated regions (3′ UTRs) of mRNAs. SRSF3-depletion also induced senescence-related phenotypes in both human and mouse cells. These 3′ UTR shortened genes were enriched in senescence-associated pathways. Shortened 3′ UTRs tended to produce more proteins than the longer ones. Simulating the effects of 3′ UTR shortening by overexpression of three candidate genes (PTEN, PIAS1 and DNMT3A) all led to senescence-associated phenotypes. Mechanistically, SRSF3 has higher binding density near proximal poly(A) site than distal one in 3′ UTR shortened genes. Further, upregulation of PTEN by either ectopic overexpression or SRSF3-knockdown induction both led to reduced phosphorylation of AKT and ultimately senescence-associated phenotypes. We revealed for the first time that reduced SRSF3 expression could promote cellular senescence through its APA-dependent function, largely extending our mechanistic understanding in splicing factor regulated cellular senescence.

SNP2APA: a database for evaluating effects of genetic variants on alternative polyadenylation in human cancers

Alternative polyadenylation (APA) is an important post-transcriptional regulation that recognizes different polyadenylation signals (PASs), resulting in transcripts with different 3' untranslated regions, thereby influencing a series of biological processes and functions. Recent studies have revealed that some single nucleotide polymorphisms (SNPs) could contribute to tumorigenesis and development through dysregulating APA. However, the associations between SNPs and APA in human cancers remain largely unknown. Here, using genotype and APA data of 9082 samples from The Cancer Genome Atlas (TCGA) and The Cancer 3'UTR Altas (TC3A), we systematically identified SNPs affecting APA events across 32 cancer types and defined them as APA quantitative trait loci (apaQTLs). As a result, a total of 467 942 cis-apaQTLs and 30 721 trans-apaQTLs were identified. By integrating apaQTLs with survival and genome-wide association studies (GWAS) data, we further identified 2154 apaQTLs associated with patient survival time and 151 342 apaQTLs located in GWAS loci. In addition, we designed an online tool to predict the effects of SNPs on PASs by utilizing PAS motif prediction tool. Finally, we developed SNP2APA, a user-friendly and intuitive database (http://gong_lab.hzau.edu.cn/SNP2APA/) for data browsing, searching, and downloading. SNP2APA will significantly improve our understanding of genetic variants and APA in human cancers.