Context-specific regulation and function of mRNA alternative polyadenylation

Alternative cleavage and polyadenylation of the Ccnb1 mRNA defines accumulation of cyclin protein during the meiotic cell cycle

Progression through the mitotic and meiotic cell cycle is driven by fluctuations in the levels of cyclins, the regulatory subunits controlling the localization and activity of CDK1 kinases. Cyclin levels are regulated through a precise balance of synthesis and degradation. Here we demonstrate that the synthesis of Cyclin B1 during the oocyte meiotic cell cycle is defined by the selective translation of mRNA variants generated through alternative cleavage and polyadenylation (APA). Using gene editing in mice, we introduced mutations into the proximal and distal polyadenylation elements of the 3' untranslated region (UTR) of the Ccnb1 mRNA. Through in vivo loss-of-function experiments, we demonstrate that the translation of mRNA with a short 3' UTR specifies Cyclin B1 protein levels that set the timing of meiotic re-entry. In contrast, translation directed by a long 3' UTR is necessary to direct Cyclin B1 protein accumulation during the MI/MII transition. These findings establish that the progression through the cell cycle is dependent on the selective translation of multiple mRNA variants generated by APA.

Alternative polyadenylation and dynamic 3' UTR length is associated with polysome recruitment throughout the cardiomyogenic differentiation of hESCs

Alternative polyadenylation (APA) increases transcript diversity through the generation of isoforms with varying 3' untranslated region (3' UTR) lengths. As the 3' UTR harbors regulatory element target sites, such as miRNAs or RNA-binding proteins, changes in this region can impact post-transcriptional regulation and translation. Moreover, the APA landscape can change based on the cell type, cell state, or condition. Given that APA events can impact protein expression, investigating translational control is crucial for comprehending the overall cellular regulation process. Revisiting data from polysome profiling followed by RNA sequencing, we investigated the cardiomyogenic differentiation of pluripotent stem cells by identifying the transcripts that show dynamic 3' UTR lengthening or shortening, which are being actively recruited to ribosome complexes. Our findings indicate that dynamic 3' UTR lengthening is not exclusively associated with differential expression during cardiomyogenesis but rather with recruitment to polysomes. We confirm that the differentiated state of cardiomyocytes shows a preference for shorter 3' UTR in comparison to the pluripotent stage although preferences vary during the days of the differentiation process. The most distinct regulatory changes are seen in day 4 of differentiation, which is the mesoderm commitment time point of cardiomyogenesis. After identifying the miRNAs that would target specifically the alternative 3' UTR region of the isoforms, we constructed a gene regulatory network for the cardiomyogenesis process, in which genes related to the cell cycle were identified. Altogether, our work sheds light on the regulation and dynamic 3' UTR changes of polysome-recruited transcripts that take place during the cardiomyogenic differentiation of pluripotent stem cells.

Comparative alternative polyadenylation profiles in differentiated adipocytes of subcutaneous and intramuscular fat tissue in cattle

Alternative polyadenylation (APA) is a key molecular mechanism involved in the post-transcriptional regulation of gene expression, which has been proven to play a critical role in cell differentiation. In the present study, we performed IVT-SAPAS sequencing to profile the dynamic changes of APA sites in bovine subcutaneous preadipocytes and intramuscular preadipocytes during adipogenesis. A total of 52621 high quality APA sites were identified in preadipocytes and adipocytes. Compared with preadipocytes, the increased usage of canonical AATAAA was observed in the cell-biased APA sites of adipocytes. Furthermore, 1933 and 2140 differentially expressed APA (DE-APA) sites, as well as 341 and 337 untranslated region-APA (UTR-APA) switching genes were identified in subcutaneous preadipocytes and intramuscular preadipocytes during adipogenesis, respectively. The UTR-APA switching genes showed divergent trends in preadipocytes, among which UTR-APA switching genes in intramuscular preadipocytes tended to use shorter 3'UTR for differentiation into mature adipocytes. APA events mediated by UTR-APA switching in intramuscular adipocytes were enriched in lipid synthesis and adipocyte differentiation. TRIB3, WWTR1, and INSIG1 played important roles in the differentiation of intramuscular preadipocytes. Briefly, our results provided new insights into understanding the mechanisms of bovine adipocyte differentiation.

Downregulation of CPSF6 leads to global mRNA 3' UTR shortening and enhanced antiviral immune responses

Alternative polyadenylation (APA) is a widespread mechanism of gene regulation that generates mRNA isoforms with alternative 3' untranslated regions (3' UTRs). Our previous study has revealed the global 3' UTR shortening of host mRNAs through APA upon viral infection. However, how the dynamic changes in the APA landscape occur upon viral infection remains largely unknown. Here we further found that, the reduced protein abundance of CPSF6, one of the core 3' processing factors, promotes the usage of proximal poly(A) sites (pPASs) of many immune related genes in macrophages and fibroblasts upon viral infection. Shortening of the 3' UTR of these transcripts may improve their mRNA stability and translation efficiency, leading to the promotion of type I IFN (IFN-I) signalling-based antiviral immune responses. In addition, dysregulated expression of CPSF6 is also observed in many immune related physiological and pathological conditions, especially in various infections and cancers. Thus, the global APA dynamics of immune genes regulated by CPSF6, can fine-tune the antiviral response as well as the responses to other cellular stresses to maintain the tissue homeostasis, which may represent a novel regulatory mechanism for antiviral immunity.

TDP-43 loss induces extensive cryptic polyadenylation in ALS/FTD

Nuclear depletion and cytoplasmic aggregation of the RNA-binding protein TDP-43 is the hallmark of ALS, occurring in over 97% of cases. A key consequence of TDP-43 nuclear loss is the de-repression of cryptic exons. Whilst TDP-43 regulated cryptic splicing is increasingly well catalogued, cryptic alternative polyadenylation (APA) events, which define the 3' end of last exons, have been largely overlooked, especially when not associated with novel upstream splice junctions. We developed a novel bioinformatic approach to reliably identify distinct APA event types: alternative last exons (ALE), 3'UTR extensions (3'Ext) and intronic polyadenylation (IPA) events. We identified novel neuronal cryptic APA sites induced by TDP-43 loss of function by systematically applying our pipeline to a compendium of publicly available and in house datasets. We find that TDP-43 binding sites and target motifs are enriched at these cryptic events and that TDP-43 can have both repressive and enhancing action on APA. Importantly, all categories of cryptic APA can also be identified in ALS and FTD post mortem brain regions with TDP-43 proteinopathy underlining their potential disease relevance. RNA-seq and Ribo-seq analyses indicate that distinct cryptic APA categories have different downstream effects on transcript and translation. Intriguingly, cryptic 3'Exts occur in multiple transcription factors, such as ELK1, SIX3, and TLX1, and lead to an increase in wild-type protein levels and function. Finally, we show that an increase in RNA stability leading to a higher cytoplasmic localisation underlies these observations. In summary, we demonstrate that TDP-43 nuclear depletion induces a novel category of cryptic RNA processing events and we expand the palette of TDP-43 loss consequences by showing this can also lead to an increase in normal protein translation.

TDP-43 nuclear loss in FTD/ALS causes widespread alternative polyadenylation changes

In frontotemporal dementia and amyotrophic lateral sclerosis, the RNA-binding protein TDP-43 is depleted from the nucleus. TDP-43 loss leads to cryptic exon inclusion but a role in other RNA processing events remains unresolved. Here, we show that loss of TDP-43 causes widespread changes in alternative polyadenylation, impacting expression of disease-relevant genes (e.g., ELP1, NEFL, and TMEM106B) and providing evidence that alternative polyadenylation is a new facet of TDP-43 pathology.

Identification of sleep and circadian alternative polyadenylation sites associated with APA-linked human brain disorders

Sleep and circadian rhythm disruptions are comorbid features of many pathologies and can negatively influence numerous health conditions, including degenerative diseases, metabolic illnesses, cancer, and various neurological disorders. Genetic association studies linking sleep and circadian disturbances with disease susceptibility have mainly focused on changes in gene expression due to mutations, such as single-nucleotide polymorphisms. Thus, associations between sleep and/or circadian rhythm and alternative polyadenylation (APA), particularly in the context of other health challenges, are largely undescribed. APA is a process that generates various transcript isoforms from the same gene, resulting in effects on mRNA translation, stability, localization, and subsequent function. Here, we have identified unique APAs in rat brain that exhibit time-of-day-dependent oscillations in expression as well as APAs that are altered by sleep deprivation and the subsequent recovery period. Genes affected by APA usage include Mapt/Tau, Ntrk2, Homer1A, Sin3band Sorl. Sorl1 has two APAs which cycle with a 24 h period, one additional APA cycles with a 12 h period and one more that is reduced during recovery sleep. Finally, we compared sleep- or circadian-associated APAs with recently described APA-linked brain disorder susceptibility genes and found 46 genes in common.

WDR33 alternative polyadenylation is dependent on stochastic poly(a) site usage and splicing efficiencies

Transcripts from the human WDR33 gene, which encodes a central component of the mRNA polyadenylation (PA) machinery, are subject to alternative polyadenylation (APA) within promoter-proximal introns/exons. This APA, which itself involves usage of multiple PA sites, results in the production of two non-canonical protein isoforms, V2 and V3, that are functionally completely unrelated to the full-length protein, with roles in innate immunity. The mechanism and regulation of WDR33 APA are unclear. Here, we report that levels of the PA factor CFIm25 modulate V2 and V3 expression, and that PA site usage of both V2 and V3 varies in distinct immune responses. Using newly developed assays to measure splicing and PA site strength, we show that splicing of V2-associated intron 6 is inefficient, allowing V2 to be produced using weak PA sites. Usage of V3's strong PA sites, on the other hand, is relatively low, reflecting the high efficiency of intron 7 splicing coupled with dependency on usage of an alternative 3' splice site within the intron. Overall, our findings demonstrate that usage of WDR33 alternative PA sites is stochastic, dependent on a complex interplay between splicing and PA, and thus provide new insights into mechanisms underlying APA.

The landscape of alternative polyadenylation during EMT and its regulation by the RNA-binding protein Quaking

Epithelial-mesenchymal transition (EMT) plays important roles in tumour progression and is orchestrated by dynamic changes in gene expression. While it is well established that post-transcriptional regulation plays a significant role in EMT, the extent of alternative polyadenylation (APA) during EMT has not yet been explored. Using 3' end anchored RNA sequencing, we mapped the alternative polyadenylation (APA) landscape following Transforming Growth Factor (TGF)-β-mediated induction of EMT in human mammary epithelial cells and found APA generally causes 3'UTR lengthening during this cell state transition. Investigation of potential mediators of APA indicated the RNA-binding protein Quaking (QKI), a splicing factor induced during EMT, regulates a subset of events including the length of its own transcript. Analysis of QKI crosslinked immunoprecipitation (CLIP)-sequencing data identified the binding of QKI within 3' untranslated regions (UTRs) was enriched near cleavage and polyadenylation sites. Following QKI knockdown, APA of many transcripts is altered to produce predominantly shorter 3'UTRs associated with reduced gene expression. These findings reveal the changes in APA that occur during EMT and identify a potential role for QKI in this process.

The anticancer compound JTE-607 reveals hidden sequence specificity of the mRNA 3' processing machinery

JTE-607 is an anticancer and anti-inflammatory compound and its active form, compound 2, directly binds to and inhibits CPSF73, the endonuclease for the cleavage step in pre-messenger RNA (pre-mRNA) 3' processing. Surprisingly, compound 2-mediated inhibition of pre-mRNA cleavage is sequence specific and the drug sensitivity is predominantly determined by sequences flanking the cleavage site (CS). Using massively parallel in vitro assays, we identified key sequence features that determine drug sensitivity. We trained a machine learning model that can predict poly(A) site (PAS) relative sensitivity to compound 2 and provide the molecular basis for understanding the impact of JTE-607 on PAS selection and transcription termination genome wide. We propose that CPSF73 and associated factors bind to the CS region in a sequence-dependent manner and the interaction affinity determines compound 2 sensitivity. These results have not only elucidated the mechanism of action of JTE-607, but also unveiled an evolutionarily conserved sequence specificity of the mRNA 3' processing machinery.

Regulation of alternative splicing and polyadenylation in neurons

Cell-type-specific gene expression is a fundamental feature of multicellular organisms and is achieved by combinations of regulatory strategies. Although cell-restricted transcription is perhaps the most widely studied mechanism, co-transcriptional and post-transcriptional processes are also central to the spatiotemporal control of gene functions. One general category of expression control involves the generation of multiple transcript isoforms from an individual gene, whose balance and cell specificity are frequently tightly regulated via diverse strategies. The nervous system makes particularly extensive use of cell-specific isoforms, specializing the neural function of genes that are expressed more broadly. Here, we review regulatory strategies and RNA-binding proteins that direct neural-specific isoform processing. These include various classes of alternative splicing and alternative polyadenylation events, both of which broadly diversify the neural transcriptome. Importantly, global alterations of splicing and alternative polyadenylation are characteristic of many neural pathologies, and recent genetic studies demonstrate how misregulation of individual neural isoforms can directly cause mutant phenotypes.

Long non-coding RNA generated from CDKN1A gene by alternative polyadenylation regulates p21 expression during DNA damage response

Alternative Polyadenylation (APA) is an emerging mechanism for dynamic changes in gene expression. Previously, we described widespread APA occurrence in introns during the DNA damage response (DDR). Here, we show that a DDR-activated APA event occurs in the first intron of CDKN1A, inducing an alternate last exon-containing lncRNA. We named this lncRNA SPUD (Selective Polyadenylation Upon DNA Damage). SPUD localizes to polysomes in the cytoplasm and is detectable as multiple isoforms in available high-throughput studies. SPUD has low abundance compared to the CDKN1A full-length isoform under non-stress conditions, and SPUD is induced in cancer and normal cells under a variety of DNA damaging conditions in part through p53. The RNA binding protein HuR binds to and promotes the stability of SPUD precursor RNA. SPUD induction increases p21 protein, but not mRNA levels, affecting p21 functions in cell-cycle, CDK2 expression and cell growth. Like CDKN1A full-length isoform, SPUD can bind two competitive p21 translational regulators, the inhibitor calreticulin and the activator CUGBP1; SPUD alters their association with CDKN1A full-length in a DDR-dependent manner, promoting CDKN1A translation. Together, these results show a new regulatory mechanism by which a lncRNA controls p21 expression post-transcriptionally, highlighting lncRNA relevance in DDR progression and cell-cycle.

A 3'UTR Insertion Is a Candidate Causal Variant at the TMEM106B Locus Associated with Increased Risk for FTLD-TDP

Background and Objectives

Single nucleotide variants near TMEM106B associate with risk of frontotemporal lobar dementia with TDP-43 inclusions (FTLD-TDP) and Alzheimer's disease (AD) in genome-wide association studies (GWAS), but the causal variant at this locus remains unclear. Here we asked whether a novel structural variant on TMEM106B is the causal variant.

Methods

An exploratory analysis identified structural variants on neurodegeneration-related genes. Subsequent analyses focused on an Alu element insertion on the 3'UTR of TMEM106B. This study included data from longitudinal aging and neurogenerative disease cohorts at Stanford University, case-control cohorts in the Alzheimer's Disease Sequencing Project (ADSP), and expression and proteomics data from Washington University in St. Louis (WUSTL). 432 individuals from two Stanford aging cohorts were whole-genome long-read and short-read sequenced. 16,906 samples from ADSP were short-read sequenced. Genotypes, transcriptomics, and proteomics data were available in 1,979 participants from an aging and dementia cohort at WUSTL. Selection criteria were specific to each cohort. In primary analyses, the linkage disequilibrium between the TMEM106B locus variants in the FTLD-TDP GWAS and the 3'UTR insertion was estimated. We then estimated linkage by ancestry in the ADSP and evaluated the effect of the TMEM106B lead variant on mRNA and protein levels.

Results

The primary analysis included 432 participants (52.5% females, age range 45-92 years old). We identified a 316 bp Alu insertion overlapping the TMEM106B 3'UTR tightly linked with top GWAS variants rs3173615(C) and rs1990622(A). In ADSP European-ancestry participants, this insertion is in equivalent linkage with rs1990622(A) (R2=0.962, D'=0.998) and rs3173615(C) (R2=0.960, D'=0.996). In African-ancestry participants, the insertion is in stronger linkage with rs1990622(A) (R2=0.992, D'=0.998) than with rs3173615(C) (R2=0.811, D'=0.994). In public datasets, rs1990622 was consistently associated with TMEM106B protein levels but not with mRNA expression. In the WUSTL dataset, rs1990622 is associated with TMEM106B protein levels in plasma and cerebrospinal fluid, but not with TMEM106B mRNA expression.

Discussion

We identified a novel Alu element insertion in the 3'UTR of TMEM106B in tight linkage with the lead FTLD-TDP risk variant. The lead variant is associated with TMEM106B protein levels, but not expression. The 3'UTR insertion is a lead candidate for the causal variant at this complex locus, pending confirmation with functional studies.

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.

CPSF6 regulates alternative polyadenylation and proliferation of cancer cells through phase separation

Cancer cells usually exhibit shortened 3' untranslated regions (UTRs) due to alternative polyadenylation (APA) to promote cell proliferation and migration. Upregulated CPSF6 leads to a systematic prolongation of 3' UTRs, but CPSF6 expression in tumors is typically higher than that in healthy tissues. This contradictory observation suggests that it is necessary to investigate the underlying mechanism by which CPSF6 regulates APA switching in cancer. Here, we find that CPSF6 can undergo liquid-liquid phase separation (LLPS), and elevated LLPS is associated with the preferential usage of the distal poly(A) sites. CLK2, a kinase upregulated in cancer cells, destructs CPSF6 LLPS by phosphorylating its arginine/serine-like domain. The reduction of CPSF6 LLPS can lead to a shortened 3' UTR of cell-cycle-related genes and accelerate cell proliferation. These results suggest that CPSF6 LLPS, rather than its expression level, may be responsible for APA regulation in cancer cells.

Downregulation of apoptotic repressor AVEN exacerbates cardiac injury after myocardial infarction

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

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.

Effect of mRNA-LNP components of two globally-marketed COVID-19 vaccines on efficacy and stability

During the COVID-19 pandemic, Pfizer-BioNTech and Moderna successfully developed nucleoside-modified mRNA lipid nanoparticle (LNP) vaccines. SARS-CoV-2 spike protein expressed by those vaccines are identical in amino acid sequence, but several key components are distinct. Here, we compared the effect of ionizable lipids, untranslated regions (UTRs), and nucleotide composition of the two vaccines, focusing on mRNA delivery, antibody generation, and long-term stability. We found that the ionizable lipid, SM-102, in Moderna's vaccine performs better than ALC-0315 in Pfizer-BioNTech's vaccine for intramuscular delivery of mRNA and antibody production in mice and long-term stability at 4 °C. Moreover, Pfizer-BioNTech's 5' UTR and Moderna's 3' UTR outperform their counterparts in their contribution to transgene expression in mice. We further found that varying N1-methylpseudouridine content at the wobble position of mRNA has little effect on vaccine efficacy. These findings may contribute to the further improvement of nucleoside-modified mRNA-LNP vaccines and therapeutics.

Infernape uncovers cell type-specific and spatially resolved alternative polyadenylation in the brain

Differential polyadenylation sites (PAs) critically regulate gene expression, but their cell type-specific usage and spatial distribution in the brain have not been systematically characterized. Here, we present Infernape, which infers and quantifies PA usage from single-cell and spatial transcriptomic data and show its application in the mouse brain. Infernape uncovers alternative intronic PAs and 3'-UTR lengthening during cortical neurogenesis. Progenitor-neuron comparisons in the excitatory and inhibitory neuron lineages show overlapping PA changes in embryonic brains, suggesting that the neural proliferation-differentiation axis plays a prominent role. In the adult mouse brain, we uncover cell type-specific PAs and visualize such events using spatial transcriptomic data. Over two dozen neurodevelopmental disorder-associated genes such as Csnk2a1 and Mecp2 show differential PAs during brain development. This study presents Infernape to identify PAs from scRNA-seq and spatial data, and highlights the role of alternative PAs in neuronal gene regulation.

Versatile RNA: overlooked gems of the transcriptome

The critical role of RNA, its use and targetability concerning different aspects of human health are gaining more attention because our understanding of the versatility of RNA has dramatically evolved over the last decades. We now appreciate that RNA is far more critical than a messenger molecule and possesses many complicated functions. As a multifunctional molecule with its sequence, flexible structures and enzymatic abilities, RNA is genuinely powerful. Mammalian transcriptomes consist of a dynamically regulated plethora of coding and noncoding RNA types. However, some aspects of RNA metabolism remain to be explored. In this Viewpoint, we focus on the transcriptome's unconventional and possibly overlooked aspects to emphasize the importance of RNA in mammalian systems.

Near chromosome-level genome assembly of the microsporidium Hamiltosporidium tvaerminnensis

Microsporidia are intracellular parasitic fungi whose genomes rank among the smallest of all known eukaryotes. A number of outstanding questions remain concerning the evolution of their large-scale variation in genome architecture, responsible for genome size variation of more than an order of magnitude. This genome report presents the first near-chromosomal assembly of a large-genome microsporidium, Hamiltosporidium tvaerminnensis. Combined Oxford Nanopore, Pacific Biosciences (PacBio), and Illumina sequencing led to a genome assembly of 17 contigs, 11 of which represent complete chromosomes. Our assembly is 21.64 Mb in length, has an N50 of 1.44 Mb, and consists of 39.56% interspersed repeats. We introduce a novel approach in microsporidia, PacBio Iso-Seq, as part of a larger annotation pipeline for obtaining high-quality annotations of 3,573 protein-coding genes. Based on direct evidence from the full-length Iso-Seq transcripts, we present evidence for alternative polyadenylation and variation in splicing efficiency, which are potential regulation mechanisms for gene expression in microsporidia. The generated high-quality genome assembly is a necessary resource for comparative genomics that will help elucidate the evolution of genome architecture in response to intracellular parasitism.

Stress responses of plants through transcriptome plasticity by mRNA alternative polyadenylation

The sessile nature of plants confines their responsiveness to changing environmental conditions. Gene expression regulation becomes a paramount mechanism for plants to adjust their physiological and morphological behaviors. Alternative polyadenylation (APA) is known for its capacity to augment transcriptome diversity and plasticity, thereby furnishing an additional set of tools for modulating gene expression. APA has also been demonstrated to exhibit intimate associations with plant stress responses. In this study, we review APA dynamic features and consequences in plants subjected to both biotic and abiotic stresses. These stresses include adverse environmental stresses, and pathogenic attacks, such as cadmium toxicity, high salt, hypoxia, oxidative stress, cold, heat shock, along with bacterial, fungal, and viral infections. We analyzed the overarching research framework employed to elucidate plant APA response and the alignment of polyadenylation site transitions with the modulation of gene expression levels within the ambit of each stress condition. We also proposed a general APA model where transacting factors, including poly(A) factors, epigenetic regulators, RNA m6A modification factors, and phase separation proteins, assume pivotal roles in APA related transcriptome plasticity during stress response in plants.

A twin UGUA motif directs the balance between gene isoforms through CFIm and the mTORC1 signaling pathway

Alternative polyadenylation (APA) generates mRNA isoforms and diversifies gene expression. Here we report the discovery that the mTORC1 signaling pathway balances the expression of two Trim9/TRIM9 isoforms through APA regulation in human and mouse. We showed that CFIm components, CPSF6 and NUDT21, promote the short Trim9/TRIM9 isoform (Trim9-S/TRIM9-S) expression. In addition, we identified an evolutionarily conserved twin UGUA motif, UGUAYUGUA, in TRIM9-S polyadenylation site (PAS) that is critical for its regulation by CPSF6. We found additional CPSF6-regulated PASs with similar twin UGUA motifs in human and experimentally validated the twin UGUA motif functionality in BMPR1B, MOB4, and BRD4-L. Importantly, we showed that inserting a twin UGUA motif into a heterologous PAS was sufficient to confer regulation by CPSF6 and mTORC1. Our study reveals an evolutionarily conserved mechanism to regulate gene isoform expression by mTORC1 and implicates possible gene isoform imbalance in cancer and neurological disorders with mTORC1 pathway dysregulation.

Identification of experimentally-supported poly(A) sites in single-cell RNA-seq data with SCINPAS

Alternative polyadenylation is a main driver of transcriptome diversity in mammals, generating transcript isoforms with different 3' ends via cleavage and polyadenylation at distinct polyadenylation (poly(A)) sites. The regulation of cell type-specific poly(A) site choice is not completely resolved, and requires quantitative poly(A) site usage data across cell types. 3' end-based single-cell RNA-seq can now be broadly used to obtain such data, enabling the identification and quantification of poly(A) sites with direct experimental support. We propose SCINPAS, a computational method to identify poly(A) sites from scRNA-seq datasets. SCINPAS modifies the read deduplication step to favor the selection of distal reads and extract those with non-templated poly(A) tails. This approach improves the resolution of poly(A) site recovery relative to standard software. SCINPAS identifies poly(A) sites in genic and non-genic regions, providing complementary information relative to other tools. The workflow is modular, and the key read deduplication step is general, enabling the use of SCINPAS in other typical analyses of single cell gene expression. Taken together, we show that SCINPAS is able to identify experimentally-supported, known and novel poly(A) sites from 3' end-based single-cell RNA sequencing data.

The predicted RNA-binding protein regulome of axonal mRNAs

Neurons are morphologically complex cells that rely on the compartmentalization of protein expression to develop and maintain their cytoarchitecture. The targeting of RNA transcripts to axons is one of the mechanisms that allows rapid local translation of proteins in response to extracellular signals. 3' Untranslated regions (UTRs) of mRNA are noncoding sequences that play a critical role in determining transcript localization and translation by interacting with specific RNA-binding proteins (RBPs). However, how 3' UTRs contribute to mRNA metabolism and the nature of RBP complexes responsible for these functions remains elusive. We performed 3' end sequencing of RNA isolated from cell bodies and axons of sympathetic neurons exposed to either nerve growth factor (NGF) or neurotrophin 3 (NTF3, also known as NT-3). NGF and NTF3 are growth factors essential for sympathetic neuron development through distinct signaling mechanisms. Whereas NTF3 acts mostly locally, NGF signal is retrogradely transported from axons to cell bodies. We discovered that both NGF and NTF3 affect transcription and alternative polyadenylation in the nucleus and induce the localization of specific 3' UTR isoforms to axons, including short 3' UTR isoforms found exclusively in axons. The integration of our data with CLIP sequencing data supports a model whereby long 3' UTR isoforms associate with RBP complexes in the nucleus and, upon reaching the axons, are remodeled locally into shorter isoforms. Our findings shed new light into the complex relationship between nuclear polyadenylation, mRNA localization, and local 3' UTR remodeling in developing neurons.

Mycobacterium tuberculosis-dependent Monocyte Expression Quantitative Trait Loci and Tuberculosis Pathogenesis

The heterogeneity of outcomes after Mycobacterium tuberculosis (Mtb) exposure is a conundrum associated with millennia of host-pathogen co-evolution. We hypothesized that human myeloid cells contain genetically encoded, Mtb-specific responses that regulate critical steps in tuberculosis (TB) pathogenesis. We mapped genome-wide expression quantitative trait loci (eQTLs) in Mtb-infected monocytes with RNAseq from 80 Ugandan household contacts of pulmonary TB cases to identify monocyte-specific, Mtb-dependent eQTLs and their association with cytokine expression and clinical resistance to tuberculin skin test (TST) and interferon-γ release assay (IGRA) conversion. cis-eQTLs (n=1,567) were identified in Mtb-infected monocytes (FDR<0.01), including 29 eQTLs in 16 genes which were Mtb-dependent (significant for Mtb:genotype interaction [FDR<0.1], but not classified as eQTL in media condition [FDR≥0.01]). A subset of eQTLs were associated with Mtb-induced cytokine expression (n=8) and/or clinical resistance to TST/IGRA conversion (n=1). Expression of BMP6, an Mtb-dependent eQTL gene, was associated with IFNB1 induction in Mtb-infected and DNA ligand-induced cells. Network and enrichment analyses identified fatty acid metabolism as a pathway associated with eQTL genes. These findings suggest that monocyte genes contain Mtb-dependent eQTLs, including a subset associated with cytokine expression and/or clinical resistance to TST/IGRA conversion, providing insight into immunogenetic pathways regulating susceptibility to Mtb infection and TB pathogenesis.

mRNA 3'UTR lengthening by alternative polyadenylation attenuates inflammatory responses and correlates with virulence of Influenza A virus

Changes of mRNA 3'UTRs by alternative polyadenylation (APA) have been associated to numerous pathologies, but the mechanisms and consequences often remain enigmatic. By combining transcriptomics, proteomics and recombinant viruses we show that all tested strains of IAV, including A/PR/8/34(H1N1) (PR8) and A/Cal/07/2009 (H1N1) (Cal09), cause APA. We mapped the effect to the highly conserved glycine residue at position 184 (G184) of the viral non-structural protein 1 (NS1). Unbiased mass spectrometry-based analyses indicate that NS1 causes APA by perturbing the function of CPSF4 and that this function is unrelated to virus-induced transcriptional shutoff. Accordingly, IAV strain PR8, expressing an NS1 variant with weak CPSF binding, does not induce host shutoff but only APA. However, recombinant IAV (PR8) expressing NS1(G184R) lacks binding to CPSF4 and thereby also the ability to cause APA. Functionally, the impaired ability to induce APA leads to an increased inflammatory cytokine production and an attenuated phenotype in a mouse infection model. Investigating diverse viral infection models showed that APA induction is a frequent ability of many pathogens. Collectively, we propose that targeting of the CPSF complex, leading to widespread alternative polyadenylation of host transcripts, constitutes a general immunevasion mechanism employed by a variety of pathogenic viruses.

The promoter as a trip navigator: Guiding alternative polyadenylation site destinations

Alfonso-Gonzalez et al.1 present an innovative combination of long-read-sequencing approaches that reveal coupling of alternative transcription start sites and alternative polyadenylation site usage on a global level.

Birth of a poly(A) tail: mechanisms and control of mRNA polyadenylation

During their synthesis in the cell nucleus, most eukaryotic mRNAs undergo a two-step 3'-end processing reaction in which the pre-mRNA is cleaved and released from the transcribing RNA polymerase II and a polyadenosine (poly(A)) tail is added to the newly formed 3'-end. These biochemical reactions might appear simple at first sight (endonucleolytic RNA cleavage and synthesis of a homopolymeric tail), but their catalysis requires a multi-faceted enzymatic machinery, the cleavage and polyadenylation complex (CPAC), which is composed of more than 20 individual protein subunits. The activity of CPAC is further orchestrated by Poly(A) Binding Proteins (PABPs), which decorate the poly(A) tail during its synthesis and guide the mRNA through subsequent gene expression steps. Here, we review the structure, molecular mechanism, and regulation of eukaryotic mRNA 3'-end processing machineries with a focus on the polyadenylation step. We concentrate on the CPAC and PABPs from mammals and the budding yeast, Saccharomyces cerevisiae, because these systems are the best-characterized at present. Comparison of their functions provides valuable insights into the principles of mRNA 3'-end processing.

Differential translation of mRNA isoforms underlies oncogenic activation of cell cycle kinase Aurora A

Aurora Kinase A (AURKA) is an oncogenic kinase with major roles in mitosis, but also exerts cell cycle- and kinase-independent functions linked to cancer. Therefore, control of its expression, as well as its activity, is crucial. A short and a long 3'UTR isoform exist for AURKA mRNA, resulting from alternative polyadenylation (APA). We initially observed that in triple-negative breast cancer, where AURKA is typically overexpressed, the short isoform is predominant and this correlates with faster relapse times of patients. The short isoform is characterized by higher translational efficiency since translation and decay rate of the long isoform are targeted by hsa-let-7a tumor-suppressor miRNA. Additionally, hsa-let-7a regulates the cell cycle periodicity of translation of the long isoform, whereas the short isoform is translated highly and constantly throughout interphase. Finally, disrupted production of the long isoform led to an increase in proliferation and migration rates of cells. In summary, we uncovered a new mechanism dependent on the cooperation between APA and miRNA targeting likely to be a route of oncogenic activation of human AURKA.

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.

mRNA Regulation by RNA Modifications

Chemical modifications on mRNA represent a critical layer of gene expression regulation. Research in this area has continued to accelerate over the last decade, as more modifications are being characterized with increasing depth and breadth. mRNA modifications have been demonstrated to influence nearly every step from the early phases of transcript synthesis in the nucleus through to their decay in the cytoplasm, but in many cases, the molecular mechanisms involved in these processes remain mysterious. Here, we highlight recent work that has elucidated the roles of mRNA modifications throughout the mRNA life cycle, describe gaps in our understanding and remaining open questions, and offer some forward-looking perspective on future directions in the field.

Alternative polyadenylation regulation in cardiac development and cardiovascular disease

Cleavage and polyadenylation of pre-mRNAs is a necessary step for gene expression and function. Majority of human genes exhibit multiple polyadenylation sites, which can be alternatively used to generate different mRNA isoforms from a single gene. Alternative polyadenylation (APA) of pre-mRNAs is important for the proteome and transcriptome landscape. APA is tightly regulated during development and contributes to tissue-specific gene regulation. Mis-regulation of APA is linked to a wide range of pathological conditions. APA-mediated gene regulation in the heart is emerging as a new area of research. Here, we will discuss the impact of APA on gene regulation during heart development and in cardiovascular diseases. First, we will briefly review how APA impacts gene regulation and discuss molecular mechanisms that control APA. Then, we will address APA regulation during heart development and its dysregulation in cardiovascular diseases. Finally, we will discuss pre-mRNA targeting strategies to correct aberrant APA patterns of essential genes for the treatment or prevention of cardiovascular diseases. The RNA field is blooming due to advancements in RNA-based technologies. RNA-based vaccines and therapies are becoming the new line of effective and safe approaches for the treatment and prevention of human diseases. Overall, this review will be influential for understanding gene regulation at the RNA level via APA in the heart and will help design RNA-based tools for the treatment of cardiovascular diseases in the future.

Sites of transcription initiation drive mRNA isoform selection

The generation of distinct messenger RNA isoforms through alternative RNA processing modulates the expression and function of genes, often in a cell-type-specific manner. Here, we assess the regulatory relationships between transcription initiation, alternative splicing, and 3' end site selection. Applying long-read sequencing to accurately represent even the longest transcripts from end to end, we quantify mRNA isoforms in Drosophila tissues, including the transcriptionally complex nervous system. We find that in Drosophila heads, as well as in human cerebral organoids, 3' end site choice is globally influenced by the site of transcription initiation (TSS). "Dominant promoters," characterized by specific epigenetic signatures including p300/CBP binding, impose a transcriptional constraint to define splice and polyadenylation variants. In vivo deletion or overexpression of dominant promoters as well as p300/CBP loss disrupted the 3' end expression landscape. Our study demonstrates the crucial impact of TSS choice on the regulation of transcript diversity and tissue identity.

The anti-cancer compound JTE-607 reveals hidden sequence specificity of the mRNA 3' processing machinery

JTE-607 is a small molecule compound with anti-inflammation and anti-cancer activities. Upon entering the cell, it is hydrolyzed to Compound 2, which directly binds to and inhibits CPSF73, the endonuclease for the cleavage step in pre-mRNA 3' processing. Although CPSF73 is universally required for mRNA 3' end formation, we have unexpectedly found that Compound 2- mediated inhibition of pre-mRNA 3' processing is sequence-specific and that the sequences flanking the cleavage site (CS) are a major determinant for drug sensitivity. By using massively parallel in vitro assays, we have measured the Compound 2 sensitivities of over 260,000 sequence variants and identified key sequence features that determine drug sensitivity. A machine learning model trained on these data can predict the impact of JTE-607 on poly(A) site (PAS) selection and transcription termination genome-wide. We propose a biochemical model in which CPSF73 and other mRNA 3' processing factors bind to RNA of the CS region in a sequence-specific manner and the affinity of such interaction determines the Compound 2 sensitivity of a PAS. As the Compound 2-resistant CS sequences, characterized by U/A-rich motifs, are prevalent in PASs from yeast to human, the CS region sequence may have more fundamental functions beyond determining drug resistance. Together, our study not only characterized the mechanism of action of a compound with clinical implications, but also revealed a previously unknown and evolutionarily conserved sequence-specificity of the mRNA 3' processing machinery.

Iron restriction increases the expression of a cytotoxic cysteine proteinase TvCP2 by a novel mechanism of tvcp2 mRNA alternative polyadenylation in Trichomonas vaginalis

Trichomonas vaginalis TvCP2 (TVAG_057000) is a cytotoxic cysteine proteinase (CP) expressed under iron-limited conditions. This work aimed to identify one of the mechanisms of tvcp2 gene expression regulation by iron at the posttranscriptional level. We checked tvcp2 mRNA stability under both iron-restricted (IR) and high iron (HI) conditions in the presence of actinomycin D. Greater stability of the tvcp2 mRNA under the IR than in HI conditions was observed, as expected. In silico analysis of the 3' regulatory region showed the presence of two putative polyadenylation signals in the tvcp2 transcript. By 3'-RACE assays, we demonstrated the existence of two isoforms of the tvcp2 mRNA with different 3'-UTR that resulted in more TvCP2 protein under IR than in HI conditions detected by WB assays. Additionally, we searched for homologs of the trichomonad polyadenylation machinery by an in silico analysis in the genome database, TrichDB. 16 genes that encode proteins that could be part of the trichomonad polyadenylation machinery were found. qRT-PCR assays showed that most of these genes were positively regulated by iron. Thus, our results show the presence of alternative polyadenylation as a novel iron posttranscriptional regulatory mechanism in T. vaginalis for the tvcp2 gene expression.

Revisiting and corrections to the annotated SRSF3 (SRp20) gene structure and RefSeq sequences from the human and mouse genomes

SRSF3 (SRp20) is the smallest member of the serine/arginine (SR)-rich protein family. We found the annotated human SRSF3 and mouse Srsf3 RefSeq sequences are much larger than the detected SRSF3/Srsf3 RNA size by Northern blot. Mapping of RNA-seq reads from various human and mouse cell lines to the annotated SRSF3/Srsf3 gene illustrated only a partial coverage of its terminal exon 7. By 5' RACE and 3' RACE, we determined that SRSF3 gene spanning over 8422 bases and Srsf3 gene spanning over 9423 bases. SRSF3/Srsf3 gene has seven exons with exon 7 bearing two alternative polyadenylation signals (PAS). Through alternative PAS selection and exon 4 exclusion/inclusion by alternative RNA splicing, SRSF3/Srsf3 gene expresses four RNA isoforms. The major SRSF3 mRNA isoform with exon 4 exclusion by using a favorable distal PAS to encode a full-length protein is 1411 nt long (not annotated 4228 nt) and the same major mouse Srsf3 mRNA isoform is only 1295 nt (not annotated 2585 nt). The difference from the redefined RNA size of SRSF3/Srsf3 to the corresponding RefSeq sequence is at the 3' UTR region. Collectively, the redefined SRSF3/Srsf3 gene structure and expression will allow better understanding of SRSF3 functions and its regulations in health and diseases.

The CRISPR/Cas9 System Delivered by Extracellular Vesicles

Clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated protein (Cas) systems can precisely manipulate DNA sequences to change the characteristics of cells and organs, which has potential in the mechanistic research on genes and the treatment of diseases. However, clinical applications are restricted by the lack of safe, targeted and effective delivery vectors. Extracellular vesicles (EVs) are an attractive delivery platform for CRISPR/Cas9. Compared with viral and other vectors, EVs present several advantages, including safety, protection, capacity, penetrating ability, targeting ability and potential for modification. Consequently, EVs are profitably used to deliver the CRISPR/Cas9 in vivo. In this review, the advantages and disadvantages of the delivery form and vectors of the CRISPR/Cas9 are concluded. The favorable traits of EVs as vectors, such as the innate characteristics, physiological and pathological functions, safety and targeting ability of EVs, are summarized. Furthermore, in terms of the delivery of the CRISPR/Cas9 by EVs, EV sources and isolation strategies, the delivery form and loading methods of the CRISPR/Cas9 and applications have been concluded and discussed. Finally, this review provides future directions of EVs as vectors of the CRISPR/Cas9 system in clinical applications, such as the safety, capacity, consistent quality, yield and targeting ability of EVs.

Novel roles of RNA-binding proteins in drug resistance of breast cancer: from molecular biology to targeting therapeutics

Therapy resistance remains a huge challenge for current breast cancer treatments. Exploring molecular mechanisms of therapy resistance might provide therapeutic targets for patients with advanced breast cancer and improve their prognosis. RNA-binding proteins (RBPs) play an important role in regulating therapy resistance. Here we summarize the functions of RBPs, highlight their tremendously important roles in regulating therapy sensitivity and resistance and we also reveal current therapeutic approaches reversing abnormal functions of RBPs in breast cancer.

Sequential Polyadenylation to Enable Alternative mRNA 3' End Formation

In eukaryotic cells, a key RNA processing step to generate mature mRNA is the coupled reaction for cleavage and polyadenylation (CPA) at the 3' end of individual transcripts. Many transcripts are alternatively polyadenylated (APA) to produce mRNAs with different 3' ends that may either alter protein coding sequence (CDS-APA) or create different lengths of 3'UTR (tandem-APA). As the CPA reaction is intimately associated with transcriptional termination, it has been widely assumed that APA is regulated cotranscriptionally. Isoforms terminated at different regions may have distinct RNA stability under different conditions, thus altering the ratio of APA isoforms. Such differential impacts on different isoforms have been considered as post-transcriptional APA, but strictly speaking, this can only be considered "apparent" APA, as the choice is not made during the CPA reaction. Interestingly, a recent study reveals sequential APA as a new mechanism for post-transcriptional APA. This minireview will focus on this new mechanism to provide insights into various documented regulatory paradigms.

3'UTR Diversity: Expanding Repertoire of RNA Alterations in Human mRNAs

Genomic information stored in the DNA is transcribed to the mRNA and translated to proteins. The 3' untranslated regions (3'UTRs) of the mRNA serve pivotal roles in posttranscriptional gene expression, regulating mRNA stability, translation, and localization. Similar to DNA mutations producing aberrant proteins, RNA alterations expand the transcriptome landscape and change the cellular proteome. Recent global analyses reveal that many genes express various forms of altered RNAs, including 3'UTR length variants. Alternative polyadenylation and alternative splicing are involved in diversifying 3'UTRs, which could act as a hidden layer of eukaryotic gene expression control. In this review, we summarize the functions and regulations of 3'UTRs and elaborate on the generation and functional consequences of 3'UTR diversity. Given that dynamic 3'UTR length control contributes to phenotypic complexity, dysregulated 3'UTR diversity might be relevant to disease development, including cancers. Thus, 3'UTR diversity in cancer could open exciting new research areas and provide avenues for novel cancer theragnostics.

The internal ribosome entry site determines the neurotropic potential of enterovirus A71

The mechanisms underlying tissue-specific replication of enteroviruses are currently unclear. Although enterovirus A71 (EV-A71) and coxsackievirus A16 (CV-A16) are both common pathogens that cause hand-foot-mouth disease, they display quite different neurotropic properties. Herein, we characterized the role of the internal ribosomal entry site (IRES) in determining neurovirulence using an oral inoculation model of EV-A71. The receptor transgenic (hSCARB2-Tg) mice developed neurological symptoms after being infected with a mouse-adapted EV-A71 strain (MP4) via different administrative routes. Intragastric administration of the MP4 strain caused pathological changes such as neuronal loss and neuropil vacuolation, whereas replacing EV-A71 IRES with CV-A16 abolished the neuropathological phenotypes. The attenuated neurotropic potential of IRES-swapped EV-A71 was observed even in mice that received intraperitoneal and intracerebral inoculations. Fewer chimeric MP4 viral RNAs and proteins were detected in the mouse tissues, regardless of the inoculation route. Tissue-specific replication can be reflected in cell-based characterization. While chimeric MP4 virus replicated poorly in human intestinal C2BBe1 and neuroblastoma SH-SY5Y cells, its replication in susceptible rhabdomyosarcoma cells was not affected. Overall, our results demonstrated that the IRES determined the neurotropic potential of EV-A71 and CV-A16, emphasizing the importance of the IRES in tissue tropism, along with the host receptors.

PolyAMiner-Bulk: A Machine Learning Based Bioinformatics Algorithm to Infer and Decode Alternative Polyadenylation Dynamics from bulk RNA-seq data

More than half of human genes exercise alternative polyadenylation (APA) and generate mRNA transcripts with varying 3' untranslated regions (UTR). However, current computational approaches for identifying cleavage and polyadenylation sites (C/PASs) and quantifying 3'UTR length changes from bulk RNA-seq data fail to unravel tissue- and disease-specific APA dynamics. Here, we developed a next-generation bioinformatics algorithm and application, PolyAMiner-Bulk, that utilizes an attention-based machine learning architecture and an improved vector projection-based engine to infer differential APA dynamics accurately. When applied to earlier studies, PolyAMiner-Bulk accurately identified more than twice the number of APA changes in an RBM17 knockdown bulk RNA-seq dataset compared to current generation tools. Moreover, on a separate dataset, PolyAMiner-Bulk revealed novel APA dynamics and pathways in scleroderma pathology and identified differential APA in a gene that was identified as being involved in scleroderma pathogenesis in an independent study. Lastly, we used PolyAMiner-Bulk to analyze the RNA-seq data of post-mortem prefrontal cortexes from the ROSMAP data consortium and unraveled novel APA dynamics in Alzheimer's Disease. Our method, PolyAMiner-Bulk, creates a paradigm for future alternative polyadenylation analysis from bulk RNA-seq data.

Long Non-Coding RNA Generated from CDKN1A Gene by Alternative Polyadenylation Regulates p21 Expression during DNA Damage Response

Alternative Polyadenylation (APA) is an emerging mechanism for dynamic changes in gene expression. Previously, we described widespread APA occurrence in introns during the DNA damage response (DDR). Here, we show that a DNA damage activated APA event occurs in the first intron of CDKN1A , inducing an alternate last exon (ALE)-containing lncRNA. We named this lncRNA SPUD (Selective Polyadenylation Upon Damage). SPUD localizes to polysomes in the cytoplasm and is detectable as multiple isoforms in available high throughput studies. SPUD has low abundance compared to the CDKN1A full-length isoform and is induced in cancer and normal cells under a variety of DNA damaging conditions in part through p53 transcriptional activation. RNA binding protein (RBP) HuR and the transcriptional repressor CTCF regulate SPUD levels. SPUD induction increases p21 protein, but not CDKN1A full-length levels, affecting p21 functions in cell-cycle, CDK2 expression, and cell viability. Like CDKN1A full-length isoform, SPUD can bind two competitive p21 translational regulators, the inhibitor calreticulin and the activator CUGBP1; SPUD can change their association with CDKN1A full-length in a DDR-dependent manner. Together, these results show a new regulatory mechanism by which a lncRNA controls p21 expression post-transcriptionally, highlighting lncRNA relevance in DDR progression and cellcycle.

Regulation and function of alternative polyadenylation in development and differentiation

Alternative processing of nascent mRNAs is widespread in eukaryotic organisms and greatly impacts the output of gene expression. Specifically, alternative cleavage and polyadenylation (APA) is a co-transcriptional molecular process that switches the polyadenylation site (PAS) at which a nascent mRNA is cleaved, resulting in mRNA isoforms with different 3'UTR length and content. APA can potentially affect mRNA translation efficiency, localization, stability, and mRNA seeded protein-protein interactions. APA naturally occurs during development and cellular differentiation, with around 70% of human genes displaying APA in particular tissues and cell types. For example, neurons tend to express mRNAs with long 3'UTRs due to preferential processing at PASs more distal than other PASs used in other cell types. In addition, changes in APA mark a variety of pathological states, including many types of cancer, in which mRNAs are preferentially cleaved at more proximal PASs, causing expression of mRNA isoforms with short 3'UTRs. Although APA has been widely reported, both the function of APA in development and the mechanisms that regulate the choice of 3'end cut sites in normal and pathogenic conditions are still poorly understood. In this review, we summarize current understanding of how APA is regulated during development and cellular differentiation and how the resulting change in 3'UTR content affects multiple aspects of gene expression. With APA being a widespread phenomenon, the advent of cutting-edge scientific techniques and the pressing need for in-vivo studies, there has never been a better time to delve into the intricate mechanisms of alternative cleavage and polyadenylation.

Interplay of four types of RNA modification writers revealed distinct tumor microenvironment and biological characteristics in pancreatic cancer

Background

Pancreatic cancer (PC) is one of the most lethal malignancies and carries a dismal mortality and morbidity. Four types of RNA modification (namely m6A, m1A, APA and A-to-I) could be catalyzed by distinct enzymatic compounds ("writers"), mediating numerous epigenetic events in carcinogenesis and immunomodulation. We aim to investigate the interplay mechanism of these writers in immunogenomic features and molecular biological characteristics in PC.

Methods

We first accessed the specific expression pattern and transcriptional variation of 26 RNA modification writers in The Cancer Genome Atlas (TCGA) dataset. Unsupervised consensus clustering was performed to divide patients into two RNA modification clusters. Then, based on the differentially expressed genes (DEGs) among two clusters, RNA modification score (WM_Score) model was established to determine RNA modification-based subtypes and was validated in International Cancer Genome Consortium (ICGC) dataset. What's more, we manifested the unique status of WM_Score in transcriptional and post-transcriptional regulation, molecular biological characteristics, targeted therapies and immunogenomic patterns.

Results

We documented the tight-knit correlations between transcriptional expression and variation of RNA modification writers. We classified patients into two distinct RNA modification patterns (WM_Score_high and _low), The WM_Score_high subgroup was correlated with worse prognosis, Th2/Th17 cell polarization and oncogenic pathways (e.g. EMT, TGF-β, and mTORC1 signaling pathways), whereas the WM_Score_low subgroup associated with favorable survival rate and Th1 cell trend. WM_Score model also proved robust predictive power in interpreting transcriptional and post-transcriptional events. Additionally, the potential targeted compounds with related pathways for the WM_Score model were further identified.

Conclusions

Our research unfolds a novel horizon on the interplay network of four RNA modifications in PC. This WM_Score model demonstrated powerful predictive capacity in epigenetic, immunological and biological landscape, providing a theoretical basis for future clinical judgments of PC.

RNA-Binding Proteins as Epigenetic Regulators of Brain Functions and Their Involvement in Neurodegeneration

A central aspect of nervous system development and function is the post-transcriptional regulation of mRNA fate, which implies time- and site-dependent translation, in response to cues originating from cell-to-cell crosstalk. Such events are fundamental for the establishment of brain cell asymmetry, as well as of long-lasting modifications of synapses (long-term potentiation: LTP), responsible for learning, memory, and higher cognitive functions. Post-transcriptional regulation is in turn dependent on RNA-binding proteins that, by recognizing and binding brief RNA sequences, base modifications, or secondary/tertiary structures, are able to control maturation, localization, stability, and translation of the transcripts. Notably, most RBPs contain intrinsically disordered regions (IDRs) that are thought to be involved in the formation of membrane-less structures, probably due to liquid-liquid phase separation (LLPS). Such structures are evidenced as a variety of granules that contain proteins and different classes of RNAs. The other side of the peculiar properties of IDRs is, however, that, under altered cellular conditions, they are also prone to form aggregates, as observed in neurodegeneration. Interestingly, RBPs, as part of both normal and aggregated complexes, are also able to enter extracellular vesicles (EVs), and in doing so, they can also reach cells other than those that produced them.

UTRdb 2.0: a comprehensive, expert curated catalog of eukaryotic mRNAs untranslated regions

The 5' and 3' untranslated regions of eukaryotic mRNAs (UTRs) play crucial roles in the post-transcriptional regulation of gene expression through the modulation of nucleo-cytoplasmic mRNA transport, translation efficiency, subcellular localization, and message stability. Since 1996, we have developed and maintained UTRdb, a specialized database of UTR sequences. Here we present UTRdb 2.0, a major update of UTRdb featuring an extensive collection of eukaryotic 5' and 3' UTR sequences, including over 26 million entries from over 6 million genes and 573 species, enriched with a curated set of functional annotations. Annotations include CAGE tags and polyA signals to label the completeness of 5' and 3'UTRs, respectively. In addition, uORFs and IRES are annotated in 5'UTRs as well as experimentally validated miRNA targets in 3'UTRs. Further annotations include evolutionarily conserved blocks, Rfam motifs, ADAR-mediated RNA editing events, and m6A modifications. A web interface allowing a flexible selection and retrieval of specific subsets of UTRs, selected according to a combination of criteria, has been implemented which also provides comprehensive download facilities. UTRdb 2.0 is accessible at http://utrdb.cloud.ba.infn.it/utrdb/.