Signals for pre-mRNA cleavage and polyadenylation

Post-Transcriptional Regulation of Gene Expression and the Intricate Life of Eukaryotic mRNAs

In recent years, there has been a growing appreciation for how regulatory events that occur either co- or post-transcriptionally contribute to the control of gene expression. Messenger RNAs (mRNAs) are extensively regulated throughout their metabolism in a precise spatiotemporal manner that requires sophisticated molecular mechanisms for cell-type-specific gene expression, which dictates cell function. Moreover, dysfunction at any of these steps can result in a variety of human diseases, including cancers, muscular atrophies, and neurological diseases. This review summarizes the steps of the central dogma of molecular biology, focusing on the post-transcriptional regulation of gene expression.

A special focus on polyadenylation and alternative polyadenylation in neurodegenerative diseases: A systematic review

Neurodegenerative diseases (NDDs) are one of the prevailing conditions characterized by progressive neuronal loss. Polyadenylation (PA) and alternative polyadenylation (APA) are the two main post-transcriptional events that regulate neuronal gene expression and protein production. This systematic review analyzed the available literature on the role of PA and APA in NDDs, with an emphasis on their contributions to disease development. A comprehensive literature search was performed using the PubMed, Scopus, Cochrane, Google Scholar, Embase, Web of Science, and ProQuest databases. The search strategy was developed based on the framework introduced by Arksey and O'Malley and supplemented by the inclusion and exclusion criteria. The study selection was performed by two independent reviewers. Extraction and data organization were performed in accordance with the predefined variables. Subsequently, quantitative and qualitative analyses were performed. Forty-seven studies were included, related to a variety of NDDs, namely Alzheimer's disease, Parkinson's disease, Huntington's disease, and amyotrophic lateral sclerosis. Disease induction was performed using different models, including human tissues, animal models, and cultured cells. Most investigations were related to PA, although some were related to APA or both. Amyloid precursor protein (APP), Tau, SNCA, and STMN2 were the major genes identified; most of the altered PA patterns were related to mRNA stability and translation efficiency. This review particularly underscores the key roles of PA and APA in the pathogenesis of NDDs through their mechanisms that contribute to gene expression dysregulation, protein aggregation, and neuronal dysfunction. Insights into these mechanisms may lead to new therapeutic strategies focused on the modulation of PA and APA activities. Further research is required to investigate the translational potential of targeting these pathways for NDD treatment.

An antisense-long-noncoding-RNA modulates p75NTR expression levels during neuronal polarization

Proper polarization of newly generated neurons is a critical process for neural network formation and brain development. The pan-neurotrophin p75NTR receptor plays a key role in this process localizing asymmetrically in one of the differentiating neurites and specifying its axonal identity in response to neurotrophins. During axonal specification, p75NTR levels are transiently modulated, yet the molecular mechanisms underlying this process are not known. Here, we identified a previously uncharacterized natural antisense transcript, AS-p75, encoded within the p75NGFR mouse gene. Using an in vitro model of polarizing murine neurons, we found that AS-p75 and p75NTR display divergent expression profiles and that p75NTR expression levels increase upon competition or depletion of AS-p75, indicating that AS-p75 is a negative regulator of p75NTR expression. Depletion of AS-p75 also results in altered p75NTR subcellular distribution and affects the polarization process. Overall, our data uncovered AS-p75 as a modulator of p75NTR expression, offering new insights into the regulation of this neurotrophin receptor during in vitro neuronal polarization.

Single-cell landscape of alternative polyadenylation in human lymphoid hematopoiesis

Alternative polyadenylation (APA) is an essential post-transcriptional process that produces mature mRNA isoforms by regulating the usage of polyadenylation sites (PASs). APA is involved in lymphocyte activation; however, its role throughout the entire differentiation trajectory remains elusive. Here, we analyzed single-cell 3'-end transcriptome data from healthy subjects to construct a dynamic-APA landscape from hematopoietic stem and progenitor cells (HSPCs) to terminally differentiated lymphocytes. This analysis covered 19973 cells of 12 clusters from five lineages (B cells, CD4+ T cells, CD8+ T cells, natural killer cells, and plasmacytoid dendritic cells). A total of 2364 genes exhibited differential 3'-untranslated region (3'UTR) PAS usage, and 3021 genes displayed differential intronic cleavage during lymphoid differentiation. We observed a global trend of 3'UTR shortening during lymphoid differentiation. Nevertheless, specific events of both 3'UTR shortening and lengthening were also identified within each cluster. The APA patterns delineated three differentiation stages: HSPCs, precursor cells, and mature cells. Moreover, we demonstrated that the conversion of naïve T cells to memory T cells was accompanied by dynamic APA in transcription factor-encoding genes (TCF7 and NFATC2IP), immune function-related genes (BCL2, CD5, CD28, GOLT1B, and TMEM59), and protein ubiquitination-related genes (UBE2G1, YPEL5, and SUMO3). These findings expand our understanding of the underlying molecular mechanisms of APA and facilitate studies on the regulatory role of APA in lymphoid hematopoiesis.

InPAS: An R/Bioconductor Package for Identifying Novel Polyadenylation Sites and Alternative Polyadenylation from Bulk RNA-seq Data

BACKGROUND

Alternative cleavage and polyadenylation (APA) is a crucial post-transcriptional gene regulation mechanism that regulates gene expression in eukaryotes by increasing the diversity and complexity of both the transcriptome and proteome. Despite the development of more than a dozen experimental methods over the last decade to identify and quantify APA events, widespread adoption of these methods has been limited by technical, financial, and time constraints. Consequently, APA remains poorly understood in most eukaryotes. However, RNA sequencing (RNA-seq) technology has revolutionized transcriptome profiling and recent studies have shown that RNA-seq data can be leveraged to identify and quantify APA events.

RESULTS

To fully capitalize on the exponentially growing RNA-seq data, we developed InPAS (Identification of Novel alternative PolyAdenylation Sites), an R/Bioconductor package for accurate identification of novel and known cleavage and polyadenylation sites (CPSs), as well as quantification of APA from RNA-seq data of various experimental designs. Compared to other APA analysis tools, InPAS offers several important advantages, including the ability to detect both novel proximal and distal CPSs, to fine tune positions of CPSs using a naïve Bayes classifier based on flanking sequence features, and to identify APA events from RNA-seq data of complex experimental designs using linear models. We benchmarked the performance of InPAS and other leading tools using simulated and experimental RNA-seq data with matched 3'-end RNA-seq data. Our results reveal that InPAS frequently outperforms existing tools in terms of precision, sensitivity, and specificity. Furthermore, we demonstrate its scalability and versatility by applying it to large, diverse RNA-seq datasets.

CONCLUSIONS

InPAS is an efficient and robust tool for identifying and quantifying APA events using readily accessible conventional RNA-seq data. Its versatility opens doors to explore APA regulation across diverse eukaryotic systems with various experimental designs. We believe that InPAS will drive APA research forward, deepening our understanding of its role in regulating gene expression, and potentially leading to the discovery of biomarkers or therapeutics for diseases.

Involvement of oncomiRs miR-23, miR-24, and miR-27 in the regulation of alternative polyadenylation in glioblastoma via CFIm25 cleavage factor

Glioblastoma multiforme (GBM) is a highly aggressive brain tumor with a poor prognosis. The cleavage factor Im 25 (CFIm25), a crucial component of the CFIm complex, plays a key role in regulating the length of the mRNA 3'-UTR and has been implicated in various cancers, including GBM. This study sought to investigate the regulatory influence of specific microRNAs (miRNAs) on CFIm25 expression in GBM, a highly aggressive brain tumor. Bioinformatics analysis identified miRNA candidates targeting CFIm25 mRNA, and gene expression profiles from the NCBI database (GSE90603) were used for further analysis. Expression levels of CFIm25 and selected miRNAs were assessed using qRT-PCR in GBM clinical samples (n = 20) and non-malignant brain tissues (n = 5). Additionally, the MTT assay was performed to examine the effect of miRNA overexpression on U251 cell viability. Lentivectors expressing the identified miRNAs were employed to experimentally validate their regulatory role on CFIm25 in U251 cell lines, and Western blot analysis was conducted to determine CFIm25 protein levels. We observed significantly increased levels of miR-23, miR-24, and miR-27 expression, associated with a marked reduction in CFIm25 expression in GBM samples compared to non-malignant brain tissues. In particular, overexpression of miR-23, miR-24, and miR-27 in U251 cells resulted in CFIm25 downregulation at both the mRNA and protein levels, while their inhibition increased CFIm25 and reduced cell proliferation. These observations strongly implicate miR-23, miR-24, and miR-27 in regulating CFIm25 expression in GBM, emphasizing their potential as promising therapeutic targets for enhancing treatment responses in glioblastoma.

Defenders of the Transcriptome: Guard Protein-Mediated mRNA Quality Control in Saccharomyces cerevisiae

Cell survival depends on precise gene expression, which is controlled sequentially. The guard proteins surveil mRNAs from their synthesis in the nucleus to their translation in the cytoplasm. Although the proteins within this group share many similarities, they play distinct roles in controlling nuclear mRNA maturation and cytoplasmic translation by supporting the degradation of faulty transcripts. Notably, this group is continuously expanding, currently including the RNA-binding proteins Npl3, Gbp2, Hrb1, Hrp1, and Nab2 in Saccharomyces cerevisiae. Some of the human serine-arginine (SR) splicing factors (SRSFs) show remarkable similarities to the yeast guard proteins and may be considered as functional homologues. Here, we provide a comprehensive summary of their crucial mRNA surveillance functions and their implications for cellular health.

The separation between mRNA-ends is more variable than expected

Effective circularization of mRNA molecules is a key step for the efficient initiation of translation. Research has shown that the intrinsic separation of the ends of mRNA molecules is rather small, suggesting that intramolecular arrangements could provide this effective circularization. Considering that the innate proximity of RNA ends might have important unknown biological implications, we aimed to determine whether the close proximity of the ends of mRNA molecules is a conserved feature across organisms and gain further insights into the functional effects of the proximity of RNA ends. To do so, we studied the secondary structure of 274 full native mRNA molecules from 17 different organisms to calculate the contour length (CL) of the external loop as an index of their end-to-end separation. Our computational predictions show bigger variations (from 0.59 to 31.8 nm) than previously reported and also than those observed in random sequences. Our results suggest that separations larger than 18.5 nm are not favored, whereas short separations could be related to phenotypical stability. Overall, our work implies the existence of a biological mechanism responsible for the increase in the observed variability, suggesting that the CL features of the exterior loop could be relevant for the initiation of translation and that a short CL could contribute to the stability of phenotypes.

Delineating yeast cleavage and polyadenylation signals using deep learning

3'-end cleavage and polyadenylation is an essential process for eukaryotic mRNA maturation. In yeast species, the polyadenylation signals that recruit the processing machinery are degenerate and remain poorly characterized compared with the well-defined regulatory elements in mammals. Here we address this issue by developing deep learning models to deconvolute degenerate cis-regulatory elements and quantify their positional importance in mediating yeast poly(A) site formation, cleavage heterogeneity, and strength. In S. cerevisiae, cleavage heterogeneity is promoted by the depletion of U-rich elements around poly(A) sites as well as multiple occurrences of upstream UA-rich elements. Sites with high cleavage heterogeneity show overall lower strength. The site strength and tandem site distances modulate alternative polyadenylation (APA) under the diauxic stress. Finally, we develop a deep learning model to reveal the distinct motif configuration of S. pombe poly(A) sites, which show more precise cleavage than S. cerevisiae Altogether, our deep learning models provide unprecedented insights into poly(A) site formation of yeast species, and our results highlight divergent poly(A) signals across distantly related species.

A comprehensive analysis of 3'UTRs in Caenorhabditis elegans

3'Untranslated regions (3'UTRs) are essential portions of genes containing elements necessary for pre-mRNA 3'end processing and are involved in post-transcriptional gene regulation. Despite their importance, they remain poorly characterized in eukaryotes. Here, we have used a multi-pronged approach to extract and curate 3'UTR data from 11533 publicly available datasets, corresponding to the entire collection of Caenorhabditis elegans transcriptomes stored in the NCBI repository from 2009 to 2023. We have also performed high throughput cloning pipelines to identify and validate rare 3'UTR isoforms and incorporated and manually curated 3'UTR isoforms from previously published datasets. This updated C. elegans 3'UTRome (v3) is the most comprehensive resource in any metazoan to date, covering 97.4% of the 20362 experimentally validated protein-coding genes with refined and updated 3'UTR boundaries for 23489 3'UTR isoforms. We also used this novel dataset to identify and characterize sequence elements involved in pre-mRNA 3'end processing and update miRNA target predictions. This resource provides important insights into the 3'UTR formation, function, and regulation in eukaryotes.

Arabidopsis and maize terminator strength is determined by GC content, polyadenylation motifs and cleavage probability

The 3' end of a gene, often called a terminator, modulates mRNA stability, localization, translation, and polyadenylation. Here, we adapted Plant STARR-seq, a massively parallel reporter assay, to measure the activity of over 50,000 terminators from the plants Arabidopsis thaliana and Zea mays. We characterize thousands of plant terminators, including many that outperform bacterial terminators commonly used in plants. Terminator activity is species-specific, differing in tobacco leaf and maize protoplast assays. While recapitulating known biology, our results reveal the relative contributions of polyadenylation motifs to terminator strength. We built a computational model to predict terminator strength and used it to conduct in silico evolution that generated optimized synthetic terminators. Additionally, we discover alternative polyadenylation sites across tens of thousands of terminators; however, the strongest terminators tend to have a dominant cleavage site. Our results establish features of plant terminator function and identify strong naturally occurring and synthetic terminators.

FUS reads histone H3K36me3 to regulate alternative polyadenylation

Complex organisms generate differential gene expression through the same set of DNA sequences in distinct cells. The communication between chromatin and RNA regulates cellular behavior in tissues. However, little is known about how chromatin, especially histone modifications, regulates RNA polyadenylation. In this study, we found that FUS was recruited to chromatin by H3K36me3 at gene bodies. The H3K36me3 recognition of FUS was mediated by the proline residues in the ZNF domain. After these proline residues were mutated or H3K36me3 was abolished, FUS dissociated from chromatin and bound more to RNA, resulting in an increase in polyadenylation sites far from stop codons genome-wide. A proline mutation corresponding to a mutation in amyotrophic lateral sclerosis contributed to the hyperactivation of mitochondria and hyperdifferentiation in mouse embryonic stem cells. These findings reveal that FUS is an H3K36me3 reader protein that links chromatin-mediated alternative polyadenylation to human disease.

A Comparative Full-Length Transcriptome Analysis Using Oxford Nanopore Technologies (ONT) in Four Tissues of Bovine Origin

The transcriptome complexity and splicing patterns in male and female cattle are ambiguous, presenting a substantial obstacle to genomic selection programs that seek to improve productivity, disease resistance, and reproduction in cattle. A comparative transcriptomic analysis using Oxford Nanopore Technologies (ONT) was conducted in bovine testes (TESTs), ovaries (OVAs), muscles (MUSCs), and livers (LIVs). An average of 5,144,769 full-length reads were obtained from each sample. The TESTs were found to have the greatest number of alternative polyadenylation (APA) events involved in processes such as sperm flagellum development and fertilization in male reproduction. In total, 438 differentially expressed transcripts (DETs) were identified in the LIVs in a comparison of females vs. males, and 214 DETs were identified in the MUSCs between females and males. Additionally, 14,735, 36,347, and 33,885 DETs were detected in MUSC vs. LIV, MUSC vs. TEST, and OVA vs. TEST comparisons, respectively, revealing the complexity of the TEST. Gene Set Enrichment Analysis (GSEA) showed that these DETs were mainly involved in the "spermatogenesis", "flagellated sperm motility", "spermatid development", "reproduction", "reproductive process", and "microtubule-based movement" KEGG pathways. Additional studies are necessary to further characterize the transcriptome in different cell types, developmental stages, and physiological conditions in bovines and ascertain the functions of the novel transcripts.

RNA binding proteins in cardiovascular development and disease

Congenital heart disease (CHD) is the most common birth defect affecting>1.35 million newborn babies worldwide. CHD can lead to prenatal, neonatal, postnatal lethality or life-long cardiac complications. RNA binding protein (RBP) mutations or variants are emerging as contributors to CHDs. RBPs are wizards of gene regulation and are major contributors to mRNA and protein landscape. However, not much is known about RBPs in the developing heart and their contributions to CHD. In this chapter, we will discuss our current knowledge about specific RBPs implicated in CHDs. We are in an exciting era to study RBPs using the currently available and highly successful RNA-based therapies and methodologies. Understanding how RBPs shape the developing heart will unveil their contributions to CHD. Identifying their target RNAs in the embryonic heart will ultimately lead to RNA-based treatments for congenital heart disease.

Arabidopsis and Maize Terminator Strength is Determined by GC Content, Polyadenylation Motifs and Cleavage Probability

The 3' end of a gene, often called a terminator, modulates mRNA stability, localization, translation, and polyadenylation. Here, we adapted Plant STARR-seq, a massively parallel reporter assay, to measure the activity of over 50,000 terminators from the plants Arabidopsis thaliana and Zea mays. We characterize thousands of plant terminators, including many that outperform bacterial terminators commonly used in plants. Terminator activity is species-specific, differing in tobacco leaf and maize protoplast assays. While recapitulating known biology, our results reveal the relative contributions of polyadenylation motifs to terminator strength. We built a computational model to predict terminator strength and used it to conduct in silico evolution that generated optimized synthetic terminators. Additionally, we discover alternative polyadenylation sites across tens of thousands of terminators; however, the strongest terminators tend to have a dominant cleavage site. Our results establish features of plant terminator function and identify strong naturally occurring and synthetic terminators.

Quantifying negative selection in human 3' UTRs uncovers constrained targets of RNA-binding proteins

Many non-coding variants associated with phenotypes occur in 3' untranslated regions (3' UTRs), and may affect interactions with RNA-binding proteins (RBPs) to regulate gene expression post-transcriptionally. However, identifying functional 3' UTR variants has proven difficult. We use allele frequencies from the Genome Aggregation Database (gnomAD) to identify classes of 3' UTR variants under strong negative selection in humans. We develop intergenic mutability-adjusted proportion singleton (iMAPS), a generalized measure related to MAPS, to quantify negative selection in non-coding regions. This approach, in conjunction with in vitro and in vivo binding data, identifies precise RBP binding sites, miRNA target sites, and polyadenylation signals (PASs) under strong selection. For each class of sites, we identify thousands of gnomAD variants under selection comparable to missense coding variants, and find that sites in core 3' UTR regions upstream of the most-used PAS are under strongest selection. Together, this work improves our understanding of selection on human genes and validates approaches for interpreting genetic variants in human 3' UTRs.

Deep learning of human polyadenylation sites at nucleotide resolution reveals molecular determinants of site usage and relevance in disease

The genomic distribution of cleavage and polyadenylation (polyA) sites should be co-evolutionally optimized with the local gene structure. Otherwise, spurious polyadenylation can cause premature transcription termination and generate aberrant proteins. To obtain mechanistic insights into polyA site optimization across the human genome, we develop deep/machine learning models to identify genome-wide putative polyA sites at unprecedented nucleotide-level resolution and calculate their strength and usage in the genomic context. Our models quantitatively measure position-specific motif importance and their crosstalk in polyA site formation and cleavage heterogeneity. The intronic site expression is governed by the surrounding splicing landscape. The usage of alternative polyA sites in terminal exons is modulated by their relative locations and distance to downstream genes. Finally, we apply our models to reveal thousands of disease- and trait-associated genetic variants altering polyadenylation activity. Altogether, our models represent a valuable resource to dissect molecular mechanisms mediating genome-wide polyA site expression and characterize their functional roles in human diseases.

Delineating yeast cleavage and polyadenylation signals using deep learning

3'-end cleavage and polyadenylation is an essential process for eukaryotic mRNA maturation. In yeast species, the polyadenylation signals that recruit the processing machinery are degenerate and remain poorly characterized compared to well-defined regulatory elements in mammals. Especially, recent deep sequencing experiments showed extensive cleavage heterogeneity for some mRNAs in Saccharomyces cerevisiae and uncovered the polyA motif differences between S. cerevisiae vs. Schizosaccharomyces pombe . The findings raised the fundamental question of how polyadenylation signals are formed in yeast. Here we addressed this question by developing deep learning models to deconvolute degenerate cis -regulatory elements and quantify their positional importance in mediating yeast polyA site formation, cleavage heterogeneity, and strength. In S. cerevisiae , cleavage heterogeneity is promoted by the depletion of U-rich elements around polyA sites as well as multiple occurrences of upstream UA-rich elements. Sites with high cleavage heterogeneity show overall lower strength. The site strength and tandem site distances modulate alternative polyadenylation (APA) under the diauxic stress. Finally, we developed a deep learning model to reveal the distinct motif configuration of S. pombe polyA sites which show more precise cleavage than S. cerevisiae . Altogether, our deep learning models provide unprecedented insights into polyA site formation across yeast species.

3'-end mRNA processing within apicomplexan parasites, a patchwork of classic, and unexpected players

The 3'-end processing of mRNA is a co-transcriptional process that leads to the formation of a poly-adenosine tail on the mRNA and directly controls termination of the RNA polymerase II juggernaut. This process involves a megadalton complex composed of cleavage and polyadenylation specificity factors (CPSFs) that are able to recognize cis-sequence elements on nascent mRNA to then carry out cleavage and polyadenylation reactions. Recent structural and biochemical studies have defined the roles played by different subunits of the complex and provided a comprehensive mechanistic understanding of this machinery in yeast or metazoans. More recently, the discovery of small molecule inhibitors of CPSF function in Apicomplexa has stimulated interest in studying the specificities of this ancient eukaryotic machinery in these organisms. Although its function is conserved in Apicomplexa, the CPSF complex integrates a novel reader of the N6-methyladenosine (m6A). This feature, inherited from the plant kingdom, bridges m6A metabolism directly to 3'-end processing and by extension, to transcription termination. In this review, we will examine convergence and divergence of CPSF within the apicomplexan parasites and explore the potential of small molecule inhibition of this machinery within these organisms. This article is categorized under: RNA Processing > 3' End Processing RNA Processing > RNA Editing and Modification.

Musashi Exerts Control of Gonadotrope Target mRNA Translation During the Mouse Estrous Cycle

The anterior pituitary controls key biological processes, including growth, metabolism, reproduction, and stress responses through distinct cell types that each secrete specific hormones. The anterior pituitary cells show a remarkable level of cell type plasticity that mediates the shifts in hormone-producing cell populations that are required to meet organismal needs. The molecular mechanisms underlying pituitary cell plasticity are not well understood. Recent work has implicated the pituitary stem cell populations and specifically, the mRNA binding proteins of the Musashi family in control of pituitary cell type identity. In this study we have identified the target mRNAs that mediate Musashi function in the adult mouse pituitary and demonstrate the requirement for Musashi function in vivo. Using Musashi RNA immunoprecipitation, we identify a cohort of 1184 mRNAs that show specific Musashi binding. Identified Musashi targets include the Gnrhr mRNA, which encodes the gonadotropin-releasing hormone receptor (GnRHR), and the Fshb mRNA, encoding follicle-stimulating hormone (FSH). Reporter assays reveal that Musashi functions to exert repression of translation of the Fshb mRNA, in addition to the previously observed repression of the Gnrhr mRNA. Importantly, mice engineered to lack Musashi in gonadotropes demonstrate a failure to repress translation of the endogenous Gnrhr and Fshb mRNAs during the estrous cycle and display a significant heterogeneity in litter sizes. The range of identified target mRNAs suggests that, in addition to these key gonadotrope proteins, Musashi may exert broad regulatory control over the pituitary proteome in a cell type-specific manner.

Elevated pre-mRNA 3' end processing activity in cancer cells renders vulnerability to inhibition of cleavage and polyadenylation

Cleavage and polyadenylation (CPA) is responsible for 3' end processing of eukaryotic poly(A)+ RNAs and preludes transcriptional termination. JTE-607, which targets CPSF-73, is the first known CPA inhibitor (CPAi) in mammalian cells. Here we show that JTE-607 perturbs gene expression through both transcriptional readthrough and alternative polyadenylation (APA). Sensitive genes are associated with features similar to those previously identified for PCF11 knockdown, underscoring a unified transcriptomic signature of CPAi. The degree of inhibition of an APA site by JTE-607 correlates with its usage level and, consistently, cells with elevated CPA activities, such as those with induced overexpression of FIP1, display greater transcriptomic disturbances when treated with JTE-607. Moreover, JTE-607 causes S phase crisis and is hence synergistic with inhibitors of DNA damage repair pathways. Together, our data reveal CPA activity and proliferation rate as determinants of CPAi-mediated cell death, raising the possibility of using CPAi as an adjunct therapy to suppress certain cancers.

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.

3'-End Processing of Eukaryotic mRNA: Machinery, Regulation, and Impact on Gene Expression

Formation of the 3' end of a eukaryotic mRNA is a key step in the production of a mature transcript. This process is mediated by a number of protein factors that cleave the pre-mRNA, add a poly(A) tail, and regulate transcription by protein dephosphorylation. Cleavage and polyadenylation specificity factor (CPSF) in humans, or cleavage and polyadenylation factor (CPF) in yeast, coordinates these enzymatic activities with each other, with RNA recognition, and with transcription. The site of pre-mRNA cleavage can strongly influence the translation, stability, and localization of the mRNA. Hence, cleavage site selection is highly regulated. The length of the poly(A) tail is also controlled to ensure that every transcript has a similar tail when it is exported from the nucleus. In this review, we summarize new mechanistic insights into mRNA 3'-end processing obtained through structural studies and biochemical reconstitution and outline outstanding questions in the field.

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

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

Mechanism of STMN2 cryptic splice-polyadenylation and its correction for TDP-43 proteinopathies

Loss of nuclear TDP-43 is a hallmark of neurodegeneration in TDP-43 proteinopathies, including amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). TDP-43 mislocalization results in cryptic splicing and polyadenylation of pre-messenger RNAs (pre-mRNAs) encoding stathmin-2 (also known as SCG10), a protein that is required for axonal regeneration. We found that TDP-43 binding to a GU-rich region sterically blocked recognition of the cryptic 3' splice site in STMN2 pre-mRNA. Targeting dCasRx or antisense oligonucleotides (ASOs) suppressed cryptic splicing, which restored axonal regeneration and stathmin-2-dependent lysosome trafficking in TDP-43-deficient human motor neurons. In mice that were gene-edited to contain human STMN2 cryptic splice-polyadenylation sequences, ASO injection into cerebral spinal fluid successfully corrected Stmn2 pre-mRNA misprocessing and restored stathmin-2 expression levels independently of TDP-43 binding.

The Dynamic Poly(A) Tail Acts as a Signal Hub in mRNA Metabolism

In eukaryotes, mRNA metabolism requires a sophisticated signaling system. Recent studies have suggested that polyadenylate tail may play a vital role in such a system. The poly(A) tail used to be regarded as a common modification at the 3' end of mRNA, but it is now known to be more than just that. It appears to act as a platform or hub that can be understood in two ways. On the one hand, polyadenylation and deadenylation machinery constantly regulates its dynamic activity; on the other hand, it exhibits the ability to recruit RNA-binding proteins and then interact with diverse factors to send various signals to regulate mRNA metabolism. In this paper, we outline the main complexes that regulate the dynamic activities of poly(A) tails, explain how these complexes participate polyadenylation/deadenylation process and summarize the diverse signals this hub emit. We are trying to make a point that the poly(A) tail can metaphorically act as a "flagman" who is supervised by polyadenylation and deadenylation and sends out signals to regulate the orderly functioning of mRNA metabolism.

Slight Variations in the Sequence Downstream of the Polyadenylation Signal Significantly Increase Transgene Expression in HEK293T and CHO Cells

Compared to transcription initiation, much less is known about transcription termination. In particular, large-scale mutagenesis studies have, so far, primarily concentrated on promoter and enhancer, but not terminator sequences. Here, we used a massively parallel reporter assay (MPRA) to systematically analyze the influence of short (8 bp) sequence variants (mutations) located downstream of the polyadenylation signal (PAS) on the steady-state mRNA level of the upstream gene, employing an eGFP reporter and human HEK293T cells as a model system. In total, we evaluated 227,755 mutations located at different overlapping positions within +17..+56 bp downstream of the PAS for their ability to regulate the reporter gene expression. We found that the positions +17..+44 bp downstream of the PAS are more essential for gene upregulation than those located more distal to the PAS, and that the mutation sequences ensuring high levels of eGFP mRNA expression are extremely T-rich. Next, we validated the positive effect of a couple of mutations identified in the MPRA screening on the eGFP and luciferase protein expression. The most promising mutation increased the expression of the reporter proteins 13-fold and sevenfold on average in HEK293T and CHO cells, respectively. Overall, these findings might be useful for further improving the efficiency of production of therapeutic products, e.g., recombinant antibodies.

Human histone pre-mRNA assembles histone or canonical mRNA-processing complexes by overlapping 3'-end sequence elements

Human pre-mRNA processing relies on multi-subunit macromolecular complexes, which recognize specific RNA sequence elements essential for assembly and activity. Canonical pre-mRNA processing proceeds via the recognition of a polyadenylation signal (PAS) and a downstream sequence element (DSE), and produces polyadenylated mature mRNAs, while replication-dependent (RD) histone pre-mRNA processing requires association with a stem-loop (SL) motif and a histone downstream element (HDE), and produces cleaved but non-polyadenylated mature mRNAs. H2AC18 mRNA, a specific H2A RD histone pre-mRNA, can be processed to give either a non-polyadenylated mRNA, ending at the histone SL, or a polyadenylated mRNA. Here, we reveal how H2AC18 captures the two human pre-mRNA processing complexes in a mutually exclusive mode by overlapping a canonical PAS (AAUAAA) sequence element with a HDE. Disruption of the PAS sequence on H2AC18 pre-mRNA prevents recruitment of the canonical complex in vitro, without affecting the histone machinery. This shows how the relative position of cis-acting elements in histone pre-mRNAs allows the selective recruitment of distinct human pre-mRNA complexes, thereby expanding the capability to regulate 3' processing and polyadenylation.

The genetic and biochemical determinants of mRNA degradation rates in mammals

BACKGROUND

Degradation rate is a fundamental aspect of mRNA metabolism, and the factors governing it remain poorly characterized. Understanding the genetic and biochemical determinants of mRNA half-life would enable more precise identification of variants that perturb gene expression through post-transcriptional gene regulatory mechanisms.

RESULTS

We establish a compendium of 39 human and 27 mouse transcriptome-wide mRNA decay rate datasets. A meta-analysis of these data identified a prevalence of technical noise and measurement bias, induced partially by the underlying experimental strategy. Correcting for these biases allowed us to derive more precise, consensus measurements of half-life which exhibit enhanced consistency between species. We trained substantially improved statistical models based upon genetic and biochemical features to better predict half-life and characterize the factors molding it. Our state-of-the-art model, Saluki, is a hybrid convolutional and recurrent deep neural network which relies only upon an mRNA sequence annotated with coding frame and splice sites to predict half-life (r=0.77). The key novel principle learned by Saluki is that the spatial positioning of splice sites, codons, and RNA-binding motifs within an mRNA is strongly associated with mRNA half-life. Saluki predicts the impact of RNA sequences and genetic mutations therein on mRNA stability, in agreement with functional measurements derived from massively parallel reporter assays.

CONCLUSIONS

Our work produces a more robust ground truth for transcriptome-wide mRNA half-lives in mammalian cells. Using these revised measurements, we trained Saluki, a model that is over 50% more accurate in predicting half-life from sequence than existing models. Saluki succinctly captures many of the known determinants of mRNA half-life and can be rapidly deployed to predict the functional consequences of arbitrary mutations in the transcriptome.

Deciphering the impact of genetic variation on human polyadenylation using APARENT2

BACKGROUND

3'-end processing by cleavage and polyadenylation is an important and finely tuned regulatory process during mRNA maturation. Numerous genetic variants are known to cause or contribute to human disorders by disrupting the cis-regulatory code of polyadenylation signals. Yet, due to the complexity of this code, variant interpretation remains challenging.

RESULTS

We introduce a residual neural network model, APARENT2, that can infer 3'-cleavage and polyadenylation from DNA sequence more accurately than any previous model. This model generalizes to the case of alternative polyadenylation (APA) for a variable number of polyadenylation signals. We demonstrate APARENT2's performance on several variant datasets, including functional reporter data and human 3' aQTLs from GTEx. We apply neural network interpretation methods to gain insights into disrupted or protective higher-order features of polyadenylation. We fine-tune APARENT2 on human tissue-resolved transcriptomic data to elucidate tissue-specific variant effects. By combining APARENT2 with models of mRNA stability, we extend aQTL effect size predictions to the entire 3' untranslated region. Finally, we perform in silico saturation mutagenesis of all human polyadenylation signals and compare the predicted effects of [Formula: see text] million variants against gnomAD. While loss-of-function variants were generally selected against, we also find specific clinical conditions linked to gain-of-function mutations. For example, we detect an association between gain-of-function mutations in the 3'-end and autism spectrum disorder. To experimentally validate APARENT2's predictions, we assayed clinically relevant variants in multiple cell lines, including microglia-derived cells.

CONCLUSIONS

A sequence-to-function model based on deep residual learning enables accurate functional interpretation of genetic variants in polyadenylation signals and, when coupled with large human variation databases, elucidates the link between functional 3'-end mutations and human health.

The transcriptional terminator XRN2 and the RNA-binding protein Sam68 link alternative polyadenylation to cell cycle progression in prostate cancer

Alternative polyadenylation (APA) yields transcripts differing in their 3'-end, and its regulation is altered in cancer, including prostate cancer. Here we have uncovered a mechanism of APA regulation impinging on the interaction between the exonuclease XRN2 and the RNA-binding protein Sam68, whose increased expression in prostate cancer is promoted by the transcription factor MYC. Genome-wide transcriptome profiling revealed a widespread impact of the Sam68/XRN2 complex on APA. XRN2 promotes recruitment of Sam68 to its target transcripts, where it competes with the cleavage and polyadenylation specificity factor for binding to strong polyadenylation signals at distal ends of genes, thus promoting usage of suboptimal proximal polyadenylation signals. This mechanism leads to 3' untranslated region shortening and translation of transcripts encoding proteins involved in G1/S progression and proliferation. Thus, our findings indicate that the APA program driven by Sam68/XRN2 promotes cell cycle progression and may represent an actionable target for therapeutic intervention.

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.

Sequence features around cleavage sites are highly conserved among different species and a critical determinant for RNA cleavage position across eukaryotes

RNA degradation is one of the critical steps for control of gene expression, and endonucleolytic cleavage-dependent RNA degradation is conserved among eukaryotes. Some cleavage sites are secondarily capped in the cytoplasm and identified using the Cap analysis of gene expression (CAGE) method. Although uncapped cleavage sites are widespread in eukaryotes, comparatively little information has been obtained about these sites using CAGE-based degradome analysis. Previously, we developed the truncated RNA-end sequencing (TREseq) method in plant species and used it to acquire comprehensive information about uncapped cleavage sites; we observed G-rich sequences near cleavage sites. However, it remains unclear whether this finding is general to other eukaryotes. In this study, we conducted TREseq analyses in fruit flies (Drosophila melanogaster) and budding yeast (Saccharomyces cerevisiae). The results revealed specific sequence features related to RNA cleavage in D. melanogaster and S. cerevisiae that were similar to sequence patterns in Arabidopsis thaliana. Although previous studies suggest that ribosome movements are important for determining cleavage position, feature selection using a random forest classifier showed that sequences around cleavage sites were major determinant for cleaved or uncleaved sites. Together, our results suggest that sequence features around cleavage sites are critical for determining cleavage position, and that sequence-specific endonucleolytic cleavage-dependent RNA degradation is highly conserved across eukaryotes.

Mpe1 senses the binding of pre-mRNA and controls 3' end processing by CPF

Most eukaryotic messenger RNAs (mRNAs) are processed at their 3' end by the cleavage and polyadenylation specificity factor (CPF/CPSF). CPF mediates the endonucleolytic cleavage of the pre-mRNA and addition of a polyadenosine (poly(A)) tail, which together define the 3' end of the mature transcript. The activation of CPF is highly regulated to maintain the fidelity of RNA processing. Here, using cryo-EM of yeast CPF, we show that the Mpe1 subunit directly contacts the polyadenylation signal sequence in nascent pre-mRNA. The region of Mpe1 that contacts RNA also promotes the activation of CPF endonuclease activity and controls polyadenylation. The Cft2 subunit of CPF antagonizes the RNA-stabilized configuration of Mpe1. In vivo, the depletion or mutation of Mpe1 leads to widespread defects in transcription termination by RNA polymerase II, resulting in transcription interference on neighboring genes. Together, our data suggest that Mpe1 plays a major role in accurate 3' end processing, activating CPF, and ensuring timely transcription termination.

Mix and Match: Promoters and Terminators for Tuning Gene Expression in the Methylotrophic Yeast Ogataea polymorpha

The yeast Ogataea polymorpha is an upcoming host for bio-manufacturing due to its unique physiological properties, including its broad substrate spectrum, and particularly its ability to utilize methanol as the sole carbon and energy source. However, metabolic engineering tools for O. polymorpha are still rare. In this study we characterized the influence of 6 promoters and 15 terminators on gene expression throughout batch cultivations with glucose, glycerol, and methanol as carbon sources as well as mixes of these carbon sources. For this characterization, a short half-life Green Fluorescent Protein (GFP) variant was chosen, which allows a precise temporal resolution of gene expression. Our promoter studies revealed how different promoters do not only influence the expression strength but also the timepoint of maximal expression. For example, the expression strength of the catalase promoter (pCAT) and the methanol oxidase promoter (pMOX) are comparable on methanol, but the maximum expression level of the pCAT is reached more than 24 h earlier. By varying the terminators, a 6-fold difference in gene expression was achieved with the MOX terminator boosting gene expression on all carbon sources by around 50% compared to the second-strongest terminator. It was shown that this exceptional increase in gene expression is achieved by the MOX terminator stabilizing the mRNA, which results in an increased transcript level in the cells. We further found that different pairing of promoters and terminators or the expression of a different gene (β-galactosidase gene) did not influence the performance of the genetic parts. Consequently, it is possible to mix and match promoters and terminators as independent elements to tune gene expression in O. polymorpha.

Deep learning modeling m6A deposition reveals the importance of downstream cis-element sequences

The N⁶-methyladenosine (m⁶A) modification is deposited to nascent transcripts on chromatin, but its site-specificity mechanism is mostly unknown. Here we model the m⁶A deposition to pre-mRNA by iM6A (intelligent m⁶A), a deep learning method, demonstrating that the site-specific m⁶A methylation is primarily determined by the flanking nucleotide sequences. iM6A accurately models the m⁶A deposition (AUROC = 0.99) and uncovers surprisingly that the cis-elements regulating the m⁶A deposition preferentially reside within the 50 nt downstream of the m⁶A sites. The m⁶A enhancers mostly include part of the RRACH motif and the m⁶A silencers generally contain CG/GT/CT motifs. Our finding is supported by both independent experimental validations and evolutionary conservation. Moreover, our work provides evidences that mutations resulting in synonymous codons can affect the m⁶A deposition and the TGA stop codon favors m⁶A deposition nearby. Our iM6A deep learning modeling enables fast paced biological discovery which would be cost-prohibitive and unpractical with traditional experimental approaches, and uncovers a key cis-regulatory mechanism for m⁶A site-specific deposition.

Nuclear and cytoplasmic poly(A) binding proteins (PABPs) favor distinct transcripts and isoforms

The poly(A)-tail appended to the 3'-end of most eukaryotic transcripts plays a key role in their stability, nuclear transport, and translation. These roles are largely mediated by Poly(A) Binding Proteins (PABPs) that coat poly(A)-tails and interact with various proteins involved in the biogenesis and function of RNA. While it is well-established that the nuclear PABP (PABPN) binds newly synthesized poly(A)-tails and is replaced by the cytoplasmic PABP (PABPC) on transcripts exported to the cytoplasm, the distribution of transcripts for different genes or isoforms of the same gene on these PABPs has not been investigated on a genome-wide scale. Here, we analyzed the identity, splicing status, poly(A)-tail size, and translation status of RNAs co-immunoprecipitated with endogenous PABPN or PABPC in human cells. At steady state, many protein-coding and non-coding RNAs exhibit strong bias for association with PABPN or PABPC. While PABPN-enriched transcripts more often were incompletely spliced and harbored longer poly(A)-tails and PABPC-enriched RNAs had longer half-lives and higher translation efficiency, there are curious outliers. Overall, our study reveals the landscape of RNAs bound by PABPN and PABPC, providing new details that support and advance the current understanding of the roles these proteins play in poly(A)-tail synthesis, maintenance, and function.

Emerging value of the viroid model in molecular biology and beyond

Viroids are single-stranded circular noncoding RNAs that infect plants. Research in the past five decades has deciphered the viroid genome structures, viroid replication cycles, numerous host factors for viroid infection, viroid motifs for intracellular and intercellular trafficking, interactions with host defense machinery, etc. In this review, we mainly focus on some significant questions that remain to be tackled, centered around (1) how the RNA polymerase II machinery performs transcription on RNA templates of nuclear-replicating viroids, (2) how viroid RNAs coordinate multiple structural elements for diverse functions, and (3) how viroid RNAs activate plant immunity. Research on viroids has led to seminal discoveries in biology, and we expect the research directions outlined in this review to continue providing key knowledge inspiring other areas of biology.

Precise gene models using long-read sequencing reveal a unique poly(A) signal in Giardia lamblia

During pre-mRNA processing, the poly(A) signal is recognized by a protein complex that ensures precise cleavage and polyadenylation of the nascent transcript. The location of this cleavage event establishes the length and sequence of the 3' UTR of an mRNA, thus determining much of its post-transcriptional fate. Using long-read sequencing, we characterize the polyadenylation signal and related sequences surrounding Giardia lamblia cleavage sites for over 2600 genes. We find that G. lamblia uses an AGURAA poly(A) signal, which differs from the mammalian AAUAAA. We also describe how G. lamblia lacks common auxiliary elements found in other eukaryotes, along with the proteins that recognize them. Further, we identify 133 genes with evidence of alternative polyadenylation. These results suggest that despite pared-down cleavage and polyadenylation machinery, 3' end formation still appears to be an important regulatory step for gene expression in G. lamblia.

Novel Circoviruses from Birds Share Common Evolutionary Roots with Fish Origin Circoviruses

Circoviruses occur in a variety of animal species and are common pathogens of mammalian and avian hosts. In our study internal organ samples of wild birds were processed for screening of circoviral sequences. Two novel viruses were identified and characterized in specimens of a little bittern and a European bee-eater that suffered from wing injuries, were weakened, had liver or kidney failures, and finally succumbed at a rescue station. The 1935 nt and 1960 nt long viral DNA genomes exhibited a genomic structure typical for circoviruses and were predicted to encode replication-associated protein in the viral strand, and a capsid protein in the complementary strand of the replicative intermediate DNA form. The genome of the newly described viruses showed 37.6% pairwise identity with each other and ≤41.5% identity with circovirus sequences, and shared a common branch with fish, human and Weddel seal circoviruses in the phylogenetic tree, implying evolutionary relationship among the ancestors of these viruses. Based on the results the little bittern and European bee-eater circoviruses represent two distinct species of the Circovirus genus, Circoviridae family.

Differential Effects on the Translation of Immune-Related Alternatively Polyadenylated mRNAs in Melanoma and T Cells by eIF4A Inhibition

Targeting the translation initiation complex eIF4F, which binds the 5' cap of mRNAs, is a promising anti-cancer approach. Silvestrol, a small molecule inhibitor of eIF4A, the RNA helicase component of eIF4F, inhibits the translation of the mRNA encoding the signal transducer and activator of transcription 1 (STAT1) transcription factor, which, in turn, reduces the transcription of the gene encoding one of the major immune checkpoint proteins, i.e., programmed death ligand-1 (PD-L1) in melanoma cells. A large proportion of human genes produce multiple mRNAs differing in their 3'-ends through the use of alternative polyadenylation (APA) sites, which, when located in alternative last exons, can generate protein isoforms, as in the STAT1 gene. Here, we provide evidence that the STAT1α, but not STAT1β protein isoform generated by APA, is required for silvestrol-dependent inhibition of PD-L1 expression in interferon-γ-treated melanoma cells. Using polysome profiling in activated T cells we find that, beyond STAT1, eIF4A inhibition downregulates the translation of some important immune-related mRNAs, such as the ones encoding TIM-3, LAG-3, IDO1, CD27 or CD137, but with little effect on the ones for BTLA and ADAR-1 and no effect on the ones encoding CTLA-4, PD-1 and CD40-L. We next apply RT-qPCR and 3'-seq (RNA-seq focused on mRNA 3' ends) on polysomal RNAs to analyze in a high throughput manner the effect of eIF4A inhibition on the translation of APA isoforms. We identify about 150 genes, including TIM-3, LAG-3, AHNAK and SEMA4D, for which silvestrol differentially inhibits the translation of APA isoforms in T cells. It is therefore crucial to consider 3'-end mRNA heterogeneity in the understanding of the anti-tumor activities of eIF4A inhibitors.

RBBP6 activates the pre-mRNA 3' end processing machinery in humans

3' end processing of most human mRNAs is carried out by the cleavage and polyadenylation specificity factor (CPSF; CPF in yeast). Endonucleolytic cleavage of the nascent pre-mRNA defines the 3' end of the mature transcript, which is important for mRNA localization, translation, and stability. Cleavage must therefore be tightly regulated. Here, we reconstituted specific and efficient 3' endonuclease activity of human CPSF with purified proteins. This required the seven-subunit CPSF as well as three additional protein factors: cleavage stimulatory factor (CStF), cleavage factor IIm (CFIIm), and, importantly, the multidomain protein RBBP6. Unlike its yeast homolog Mpe1, which is a stable subunit of CPF, RBBP6 does not copurify with CPSF and is recruited in an RNA-dependent manner. Sequence and mutational analyses suggest that RBBP6 interacts with the WDR33 and CPSF73 subunits of CPSF. Thus, it is likely that the role of RBBP6 is conserved from yeast to humans. Overall, our data are consistent with CPSF endonuclease activation and site-specific pre-mRNA cleavage being highly controlled to maintain fidelity in mRNA processing.

Integrative profiling of Epstein-Barr virus transcriptome using a multiplatform approach

BACKGROUND

Epstein-Barr virus (EBV) is an important human pathogenic gammaherpesvirus with carcinogenic potential. The EBV transcriptome has previously been analyzed using both Illumina-based short read-sequencing and Pacific Biosciences RS II-based long-read sequencing technologies. Since the various sequencing methods have distinct strengths and limitations, the use of multiplatform approaches have proven to be valuable. The aim of this study is to provide a more complete picture on the transcriptomic architecture of EBV.

METHODS

In this work, we apply the Oxford Nanopore Technologies MinION (long-read sequencing) platform for the generation of novel transcriptomic data, and integrate these with other's data generated by another LRS approach, Pacific BioSciences RSII sequencing and Illumina CAGE-Seq and Poly(A)-Seq approaches. Both amplified and non-amplified cDNA sequencings were applied for the generation of sequencing reads, including both oligo-d(T) and random oligonucleotide-primed reverse transcription. EBV transcripts are identified and annotated using the LoRTIA software suite developed in our laboratory.

RESULTS

This study detected novel genes embedded into longer host genes containing 5'-truncated in-frame open reading frames, which potentially encode N-terminally truncated proteins. We also detected a number of novel non-coding RNAs and transcript length isoforms encoded by the same genes but differing in their start and/or end sites. This study also reports the discovery of novel splice isoforms, many of which may represent altered coding potential, and of novel replication-origin-associated transcripts. Additionally, novel mono- and multigenic transcripts were identified. An intricate meshwork of transcriptional overlaps was revealed.

CONCLUSIONS

An integrative approach applying multi-technique sequencing technologies is suitable for reliable identification of complex transcriptomes because each techniques has different advantages and limitations, and the they can be used for the validation of the results obtained by a particular approach.

Architectural and functional details of CF IA proteins involved in yeast 3'-end pre-mRNA processing and its significance for eukaryotes: A concise review

In eukaryotes, maturation of pre-mRNA relies on its precise 3'-end processing. This processing involves co-transcriptional steps regulated by sequence elements and other proteins. Although, it holds tremendous importance, defect in the processing machinery will result in erroneous pre-mRNA maturation leading to defective translation. Remarkably, more than 20 proteins in humans and yeast share homology and execute this processing. The defects in this processing are associated with various diseases in humans. We shed light on the CF IA subunit of yeast Saccharomyces cerevisiae that contains four proteins (Pcf11, Clp1, Rna14 and Rna15) involved in this processing. Structural details of various domains of CF IA and their roles during 3'-end processing, like cleavage and polyadenylation at 3'-UTR of pre-mRNA and other cellular events are explained. Further, the chronological development and important discoveries associated with 3'-end processing are summarized. Moreover, the mammalian homologues of yeast CF IA proteins, along with their key roles are described. This knowledge would be helpful for better comprehension of the mechanism associated with this marvel; thus opening up vast avenues in this area.

Benchmarking sequencing methods and tools that facilitate the study of alternative polyadenylation

BACKGROUND

Alternative cleavage and polyadenylation (APA), an RNA processing event, occurs in over 70% of human protein-coding genes. APA results in mRNA transcripts with distinct 3' ends. Most APA occurs within 3' UTRs, which harbor regulatory elements that can impact mRNA stability, translation, and localization.

RESULTS

APA can be profiled using a number of established computational tools that infer polyadenylation sites from standard, short-read RNA-seq datasets. Here, we benchmarked a number of such tools-TAPAS, QAPA, DaPars2, GETUTR, and APATrap- against 3'-Seq, a specialized RNA-seq protocol that enriches for reads at the 3' ends of genes, and Iso-Seq, a Pacific Biosciences (PacBio) single-molecule full-length RNA-seq method in their ability to identify polyadenylation sites and quantify polyadenylation site usage. We demonstrate that 3'-Seq and Iso-Seq are able to identify and quantify the usage of polyadenylation sites more reliably than computational tools that take short-read RNA-seq as input. However, we find that running one such tool, QAPA, with a set of polyadenylation site annotations derived from small quantities of 3'-Seq or Iso-Seq can reliably quantify variation in APA across conditions, such asacross genotypes, as demonstrated by the successful mapping of alternative polyadenylation quantitative trait loci (apaQTL).

CONCLUSIONS

We envisage that our analyses will shed light on the advantages of studying APA with more specialized sequencing protocols, such as 3'-Seq or Iso-Seq, and the limitations of studying APA with short-read RNA-seq. We provide a computational pipeline to aid in the identification of polyadenylation sites and quantification of polyadenylation site usages using Iso-Seq data as input.

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.

Emerging roles of alternative cleavage and polyadenylation (APA) in human disease

In the messenger RNA (mRNA) maturation process, the 3'-end of pre-mRNA is cleaved and a poly(A) sequence is added, this is an important determinant of mRNA stability and its cellular functions. More than 60%-70% of human genes have three or more polyadenylation (APA) sites and can be cleaved at different sites, generating mRNA transcripts of varying lengths. This phenomenon is termed as alternative cleavage and polyadenylation (APA) and it plays role in key biological processes like gene regulation, cell proliferation, senescence, and also in various human diseases. Loss of regulatory microRNA binding sites and interactions with RNA-binding proteins leading to APA are largely investigated in human diseases. However, the functions of the core APA machinery and related factors during disease conditions remain largely unknown. In this review, we discuss the roles of polyadenylation machinery in relation to brain disease, cardiac failure, pulmonary fibrosis, cancer, infectious conditions, and other human diseases. Collectively, we believe this review will be a useful avenue for understanding the emerging role of APA in the pathobiology of various human diseases.

Interplay between splicing and transcriptional pausing exerts genome-wide control over alternative polyadenylation

Recent studies have identified multiple polyadenylation sites in nearly all mammalian genes. Although these are interpreted as evidence for alternative polyadenylation, our knowledge of the underlying mechanisms is still limited. Most studies only consider the immediate surroundings of gene ends, even though in vitro experiments have uncovered the involvement of external factors such as splicing. Whereas in vivo splicing manipulation was impracticable until recently, we now used mutants in the Death Inducer Obliterator (DIDO) gene to study their impact on 3ʹ end processing. We observe multiple rounds of readthrough and gene fusions, suggesting that no arbitration between polyadenylation sites occurs. Instead, a window of opportunity seems to control end processing. Through the identification of T-rich sequence motifs, our data indicate that splicing and transcriptional pausing interact to regulate alternative polyadenylation. We propose that 3ʹ splice site activation comprises a variable timer, which determines how long transcription proceeds before polyadenylation signals are recognized. Thus, the role of core polyadenylation signals could be more passive than commonly believed. Our results provide new insights into the mechanisms of alternative polyadenylation and expand the catalog of related aberrations.

Abbreviations APA: alternative polyadenylation; bp: basepair; MEF: mouse embryonic fibroblasts; PA: polyadenylation; PAS: polyadenylation site; Pol II: (RNA) polymerase II ; RT-PCR:reverse-transcriptase PCR; SF:splicing factor; SFPQ:splicing factor rich in proline and glutamine; SS:splice site; TRSM:Thymidine rich sequence motif; UTR:untranslated terminal region

MAAPER: model-based analysis of alternative polyadenylation using 3' end-linked reads

Most eukaryotic genes express alternative polyadenylation (APA) isoforms. A growing number of RNA sequencing methods, especially those used for single-cell transcriptome analysis, generate reads close to the polyadenylation site (PAS), termed nearSite reads, hence inherently containing information about APA isoform abundance. Here, we present a probabilistic model-based method named MAAPER to utilize nearSite reads for APA analysis. MAAPER predicts PASs with high accuracy and sensitivity and examines different types of APA events with robust statistics. We show MAAPER's performance with both bulk and single-cell data and its applicability in unpaired or paired experimental designs.

CRISPRpas: programmable regulation of alternative polyadenylation by dCas9

Most human protein-coding genes produce alternative polyadenylation (APA) isoforms that differ in 3' UTR size or, when coupled with splicing, have variable coding sequences. APA is an important layer of gene expression program critical for defining cell identity. Here, by using a catalytically dead Cas9 and coupling its target site with polyadenylation site (PAS), we develop a method, named CRISPRpas, to alter APA isoform abundance. CRISPRpas functions by enhancing proximal PAS usage, whose efficiency is influenced by several factors, including targeting strand of DNA, distance between PAS and target sequence and strength of the PAS. For intronic polyadenylation (IPA), splicing features, such as strengths of 5' splice site and 3' splice site, also affect CRISPRpas efficiency. We show modulation of APA of multiple endogenous genes, including IPA of PCF11, a master regulator of APA and gene expression. In sum, CRISPRpas offers a programmable tool for APA regulation that impacts gene expression.