Mechanisms and consequences of alternative polyadenylation

DNA Methylation Regulates Alternative Polyadenylation via CTCF and the Cohesin Complex

Dysregulation of DNA methylation and mRNA alternative cleavage and polyadenylation (APA) are both prevalent in cancer and have been studied as independent processes. We discovered a DNA methylation-regulated APA mechanism when we compared genome-wide DNA methylation and polyadenylation site usage between DNA methylation-competent and DNA methylation-deficient cells. Here, we show that removal of DNA methylation enables CTCF binding and recruitment of the cohesin complex, which, in turn, form chromatin loops that promote proximal polyadenylation site usage. In this DNA demethylated context, either deletion of the CTCF binding site or depletion of RAD21 cohesin complex protein can recover distal polyadenylation site usage. Using data from The Cancer Genome Atlas, we authenticated the relationship between DNA methylation and mRNA polyadenylation isoform expression in vivo. This DNA methylation-regulated APA mechanism demonstrates how aberrant DNA methylation impacts transcriptome diversity and highlights the potential sequelae of global DNA methylation inhibition as a cancer treatment.

Star-PAP controlled alternative polyadenylation coupled poly(A) tail length regulates protein expression in hypertrophic heart

Alternative polyadenylation (APA)-mediated 3′-untranslated region (UTR) shortening is known to increase protein expression due to the loss of miRNA regulatory sites. Yet, mRNAs with longer 3′-UTR also show enhanced protein expression. Here, we identify a mechanism by which longer transcripts generated by the distal-most APA site leads to increased protein expression compared to the shorter transcripts and the longer transcripts are positioned to regulate heart failure (HF). A Star-PAP target gene, NQO1 has three poly(A) sites (PA-sites) at the terminal exon on the pre-mRNA. Star-PAP selects the distal-most site that results in the expression of the longest isoform. We show that the NQO1 distal-specific mRNA isoform accounts for the majority of cellular NQO1 protein. Star-PAP control of the distal-specific isoform is stimulated by oxidative stress and the toxin dioxin. The longest NQO1 transcript has increased poly(A) tail (PA-tail) length that accounts for the difference in translation potentials of the three NQO1 isoforms. This mechanism is involved in the regulation of cardiac hypertrophy (CH), an antecedent condition to HF where NQO1 downregulation stems from the loss of the distal-specific transcript. The loss of NQO1 during hypertrophy was rescued by ectopic expression of the distal- but not the proximal- or middle-specific NQO1 mRNA isoforms in the presence of Star-PAP expression, and reverses molecular events of hypertrophy in cardiomyocytes.

RNA modifications in brain tumorigenesis

RNA modifications are emerging as critical regulators in cancer biology, thanks to their ability to influence gene expression and the predominant protein isoforms expressed during cell proliferation, migration, and other pro-oncogenic properties. The reversibility and dynamic nature of post-transcriptional RNA modifications allow cells to quickly adapt to microenvironmental changes. Recent literature has revealed that the deregulation of RNA modifications can promote a plethora of developmental diseases, including tumorigenesis. In this review, we will focus on four key post-transcriptional RNA modifications which have been identified as contributors to the pathogenesis of brain tumors: m6A, alternative polyadenylation, alternative splicing and adenosine to inosine modifications. In addition to the role of RNA modifications in brain tumor progression, we will also discuss potential opportunities to target these processes to improve the dismal prognosis for brain tumors.

Characterization and Functional Analysis of Polyadenylation Sites in Fast and Slow Muscles

Many increasing documents have proved that alternative polyadenylation (APA) events with different polyadenylation sites (PAS) contribute to posttranscriptional regulation. However, little is known about the detailed molecular features of PASs and its role in porcine fast and slow skeletal muscles through microRNAs (miRNAs) and RNA binding proteins (RBPs). In this study, we combined single-molecule real-time sequencing and Illumina RNA-seq datasets to comprehensively analyze polyadenylation in pigs. We identified a total of 10,334 PASs, of which 8734 were characterized by reference genome annotation. 32.86% of PAS-associated genes were determined to have more than one PAS. Further analysis demonstrated that tissue-specific PASs between fast and slow muscles were enriched in skeletal muscle development pathways. In addition, we obtained 1407 target genes regulated by APA events through potential binding 69 miRNAs and 28 RBPs in variable 3′ UTR regions and some are involved in myofiber transformation. Furthermore, the de novo motif search confirmed that the most common usage of canonical motif AAUAAA and three types of PASs may be related to the strength of motifs. In summary, our results provide a useful annotation of PASs for pig transcriptome and suggest that APA may serve as a role in fast and slow muscle development under the regulation of miRNAs and RBPs.

Differential Alternative Polyadenylation Landscapes Mediate Hematopoietic Stem Cell Activation and Regulate Glutamine Metabolism

Alternative polyadenylation (APA) is emerging as an important regulatory mechanism of RNA and protein isoform expression by controlling 3' untranslated region (3'-UTR) composition. The relevance of APA in stem cell hierarchies remains elusive. Here, we first demonstrate the requirement of the APA regulator Pabpn1 for hematopoietic stem cell (HSC) function. We then determine the genome-wide APA landscape (APAome) of HSCs and progenitors by performing low-input 3' sequencing paired with bioinformatic pipelines. This reveals transcriptome-wide dynamic APA patterns and an overall shortening of 3'-UTRs during differentiation and upon homeostatic or stress-induced transition from quiescence to proliferation. Specifically, we show that APA regulates activation-induced Glutaminase (Gls) isoform switching by Nudt21. This adaptation of the glutamine metabolism by increasing the GAC:KGA isoform ratio fuels versatile metabolic pathways necessary for HSC self-renewal and proper stress response. Our study establishes APA as a critical regulatory layer orchestrating HSC self-renewal, behavior, and commitment.

hPCL3S promotes proliferation and migration of androgen-independent prostate cancer cells

Polycomb repressive complex 2 (PRC2) allows the deposition of H3K27me3. PRC2 facultative subunits modulate its activity and recruitment such as hPCL3/PHF19, a human ortholog of Drosophila Polycomb-like protein (PCL). These proteins contain a TUDOR domain binding H3K36me3, two PHD domains and a “Winged-helix” domain involved in GC-rich DNA binding. The human PCL3 locus encodes the full-length hPCL3L protein and a shorter isoform, hPCL3S containing the TUDOR and PHD1 domains only.

In this study, we demonstrated by RT-qPCR analyses of 25 prostate tumors that hPCL3S is frequently up-regulated. In addition, hPCL3S is overexpressed in the androgen-independent DU145 and PC3 cells, but not in the androgen-dependent LNCaP cells. hPCL3S knockdown decreased the proliferation and migration of DU145 and PC3 whereas its forced expression into LNCaP increased these properties. A mutant hPCL3S unable to bind H3K36me3 (TUDOR-W50A) increased proliferation and migration of LNCaP similarly to wt hPCL3S whereas inactivation of its PHD1 domain decreased proliferation. These effects partially relied on the up-regulation of genes known to be important for the proliferation and/or migration of prostate cancer cells such as S100A16, PlexinA2, and Spondin1.

Collectively, our results suggest hPCL3S as a new potential therapeutic target in castration resistant prostate cancers.

Alternative Polyadenylation and Salicylic Acid Modulate Root Responses to Low Nitrogen Availability

Nitrogen (N) is probably the most important macronutrient and its scarcity limits plant growth, development and fitness. N starvation response has been largely studied by transcriptomic analyses, but little is known about the role of alternative polyadenylation (APA) in such response. In this work, we show that N starvation modifies poly(A) usage in a large number of transcripts, some of them mediated by FIP1, a component of the polyadenylation machinery. Interestingly, the number of mRNAs isoforms with poly(A) tags located in protein-coding regions or 5'-UTRs significantly increases in response to N starvation. The set of genes affected by APA in response to N deficiency is enriched in N-metabolism, oxidation-reduction processes, response to stresses, and hormone responses, among others. A hormone profile analysis shows that the levels of salicylic acid (SA), a phytohormone that reduces nitrate accumulation and root growth, increase significantly upon N starvation. Meta-analyses of APA-affected and fip1-2-deregulated genes indicate a connection between the nitrogen starvation response and salicylic acid (SA) signaling. Genetic analyses show that SA may be important for preventing the overgrowth of the root system in low N environments. This work provides new insights on how plants interconnect different pathways, such as defense-related hormonal signaling and the regulation of genomic information by APA, to fine-tune the response to low N availability.

Stiff matrix instigates type I collagen biogenesis by mammalian cleavage factor I complex-mediated alternative polyadenylation

Alternative polyadenylation (APA) is a widespread and important mechanism in regulation of gene expression. Dysregulation of the 3' UTR cleavage and polyadenylation represents a common characteristic among many disease states, including lung fibrosis. In this study, we investigated the role of mammalian cleavage factor I-mediated (CFIm-mediated) APA in regulating extracellular matrix production in response to mechanical stimuli from stiffened matrix simulating the fibrotic lungs. We found that stiff matrix downregulated expression of CFIm68, CFIm59 and CFIm25 subunits and promoted APA in favor of the proximal poly(A) site usage in the 3' UTRs of type I collagen (COL1A1) and fibronectin (FN1) in primary human lung fibroblasts. Knockdown and overexpression of each individual CFIm subunit demonstrated that CFIm68 and CFIm25 are indispensable attributes of stiff matrix-induced APA and overproduction of COL1A1, whereas CFIm did not appear to mediate stiffness-regulated FN1 APA. Furthermore, expression of the CFIm subunits was associated with matrix stiffness in vivo in a bleomycin-induced mouse model of pulmonary fibrosis. These data suggest that stiff matrix instigates type I collagen biogenesis by selectively targeting mRNA transcripts for 3' UTR shortening. The current study uncovered a potential mechanism for regulation of the CFIm complex by mechanical cues under fibrotic conditions.

Alternative polyadenylation of the stacyose synthase gene mediates source-sink regulation in cucumber

Alternative polyadenylation (APA) is a pervasive mechanism for gene regulation in eukaryotes. Stachyose is the main assimilate translocated in the cucumber phloem. Stachyose synthase (CsSTS) catalyzes the last step of stachyose biosynthesis in cucumber leaves and plays a key role in the regulation of assimilate partitioning between source and sink. In this study, three CsSTS mRNAs with the same open reading frame and the 5`untranslated region (UTR), but differing in their 3`UTRs, named CsSTS1 (short), CsSTS2 (medium), and CsSTS3 (long), were identified. Southern blot and sequence analysis of the cucumber genome confirmed that these transcripts are regulated through APA from a single gene. No significant difference of in vitro translation efficiency was found among three mRNAs. However, the relative stabilities of three transcripts varied among different tissues and different leaf development stages of cucumber. CsSTS1 expression in cucumber calli was up-regulated by the raffinose (substrate of CsSTS) and down-regulated by stachyose (product of CsSTS), respectively. In cucumber plants, all three isoforms have considerable expression in non-fruit node leaves. However, in fruit-carrying node leaves, the expression of CsSTS2 and CsSTS3 was severely inhibited and only CsSTS1 was highly expressed, indicating fruit setting has a remarkable effect on the relative expression level of three transcripts. This "fruit setting" effect could be observed until at least 36 h after the fruit was removed from the node. Our results suggest that abundant expression of CsSTS1 is beneficial for stachyose loading in source leaves, and APA is a delicate mechanism for CsSTS to regulate cucumber source-sink balance.

U1 snRNP regulates cancer cell migration and invasion in vitro

Stimulated cells and cancer cells have widespread shortening of mRNA 3'-untranslated regions (3'UTRs) and switches to shorter mRNA isoforms due to usage of more proximal polyadenylation signals (PASs) in introns and last exons. U1 snRNP (U1), vertebrates' most abundant non-coding (spliceosomal) small nuclear RNA, silences proximal PASs and its inhibition with antisense morpholino oligonucleotides (U1 AMO) triggers widespread premature transcription termination and mRNA shortening. Here we show that low U1 AMO doses increase cancer cells' migration and invasion in vitro by up to 500%, whereas U1 over-expression has the opposite effect. In addition to 3'UTR length, numerous transcriptome changes that could contribute to this phenotype are observed, including alternative splicing, and mRNA expression levels of proto-oncogenes and tumor suppressors. These findings reveal an unexpected role for U1 homeostasis (available U1 relative to transcription) in oncogenic and activated cell states, and suggest U1 as a potential target for their modulation.

Transcription Terminator-Mediated Enhancement in Transgene Expression in Maize: Preponderance of the AUGAAU Motif Overlapping With Poly(A) Signals

The selection of transcription terminators (TTs) for pairing with high expressing constitutive promoters in chimeric constructs is crucial to deliver optimal transgene expression in plants. In this study, the use of the native combinations of four polyubiquitin gene promoters and corresponding TTs resulted in up to >3-fold increase in transgene expression in maize. Of the eight polyubiquitin promoter and TT regulatory elements utilized, seven were novel and identified from the polyubiquitin genes of Brachypodium distachyon, Setaria italica, and Zea mays. Furthermore, gene expression driven by the Cassava mosaic virus promoter was studied by pairing the promoter with distinct TTs derived from the high expressing genes of Arabidopsis. Of the three TTs studied, the polyubiquitin10 gene TT produced the highest transgene expression in maize. Polyadenylation patterns and mRNA abundance from eight distinct TTs were analyzed using 3'-RACE and next-generation sequencing. The results exhibited one to three unique polyadenylation sites in the TTs. The poly(A) site patterns for the StPinII TT were consistent when the same TT was deployed in chimeric constructs irrespective of the reporter gene and promoter used. Distal to the poly(A) sites, putative polyadenylation signals were identified in the near-upstream regions of the TTs based on previously reported mutagenesis and bioinformatics studies in rice and Arabidopsis. The putative polyadenylation signals were 9 to 11 nucleotides in length. Six of the eight TTs contained the putative polyadenylation signals that were overlaps of either canonical AAUAAA or AAUAAA-like polyadenylation signals and AUGAAU, a top-ranking-hexamer of rice and Arabidopsis gene near-upstream regions. Three of the polyubiquitin gene TTs contained the identical 9-nucleotide overlap, AUGAAUAAG, underscoring the functional significance of such overlaps in mRNA 3' end processing. In addition to identifying new combinations of regulatory elements for high constitutive trait gene expression in maize, this study demonstrated the importance of TTs for optimizing gene expression in plants. Learning from this study could be applied to other dicotyledonous and monocotyledonous plant species for transgene expression. Research on TTs is not limited to transgene expression but could be extended to the introduction of appropriate mutations into TTs via genome editing, paving the way for expression modulation of endogenous genes.

SAM68-Specific Splicing Is Required for Proper Selection of Alternative 3' UTR Isoforms in the Nervous System

Neuronal alternative splicing is a core mechanism for functional diversification. We previously found that STAR family proteins (SAM68, SLM1, SLM2) regulate spatiotemporal alternative splicing in the nervous system. However, the whole aspect of alternative splicing programs by STARs remains unclear. Here, we performed a transcriptomic analysis using SAM68 knockout and SAM68/SLM1 double-knockout midbrains. We revealed different alternative splicing activity between SAM68 and SLM1; SAM68 preferentially targets alternative 3′ UTR exons. SAM68 knockout causes a long-to-short isoform switch of a number of neuronal targets through the alteration in alternative last exon (ALE) selection or alternative polyadenylation. The altered ALE usage of a novel target, interleukin 1 receptor accessory protein (Il1rap), results in remarkable conversion from a membrane-bound type to a secreted type in Sam68^KO brains. Proper ALE selection is necessary for IL1RAP neuronal function. Thus the SAM68-specific splicing program provides a mechanism for neuronal selection of alternative 3′ UTR isoforms.

•

SAM68 and the related protein SLM1 exhibit distinct alternative splicing activity

•

SAM68 specifically controls 3′ UTR selection of multiple neuronal genes

•

Proper 3′ UTR selection is necessary for IL1RAP neuronal function

•

Neuronal expression of SAM68 requires proper 3′ UTR selection in the nervous system

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

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

Inference of the human polyadenylation code

Motivation

Processing of transcripts at the 3′-end involves cleavage at a polyadenylation site followed by the addition of a poly(A)-tail. By selecting which site is cleaved, the process of alternative polyadenylation enables genes to produce transcript isoforms with different 3′-ends. To facilitate the identification and treatment of disease-causing mutations that affect polyadenylation and to understand the sequence determinants underlying this regulatory process, a computational model that can accurately predict polyadenylation patterns from genomic features is desirable.

Results

Previous works have focused on identifying candidate polyadenylation sites and classifying tissue-specific sites. By training on how multiple sites in genes are competitively selected for polyadenylation from 3′-end sequencing data, we developed a deep learning model that can predict the tissue-specific strength of a polyadenylation site in the 3′ untranslated region of the human genome given only its genomic sequence. We demonstrate the model’s broad utility on multiple tasks, without any application-specific training. The model can be used to predict which polyadenylation site is more likely to be selected in genes with multiple sites. It can be used to scan the 3′ untranslated region to find candidate polyadenylation sites. It can be used to classify the pathogenicity of variants near annotated polyadenylation sites in ClinVar. It can also be used to anticipate the effect of antisense oligonucleotide experiments to redirect polyadenylation. We provide analysis on how different features affect the model’s predictive performance and a method to identify sensitive regions of the genome at the single-based resolution that can affect polyadenylation regulation.

Supplementary information

Supplementary data are available at Bioinformatics online.

NUDT21 inhibits bladder cancer progression through ANXA2 and LIMK2 by alternative polyadenylation

Purpose: Nudix Hydrolase 21 (NUDT21) is a crucial mediator involved in alternative polyadenylation (APA), and this molecule has been reported to be a tumor suppressor in human cancers. However, neither the role NUDT21 plays in bladder cancer (BC) nor the mechanisms which are involved have been investigated.

Methods: Expression levels of NUDT21 in BC were evaluated with real-time PCR, western blotting, and immunohistochemistry (IHC). In vitro and in vivo assays were performed to investigate the function of NUDT21 in tumorigenesis in bladder cancer cells. The TOP/FOP flash reporter assay, western blot, and global APA site profiling analysis were used to identify the pathway which mediates the biologic roles of NUDT21 in BC.

Results: NUDT21 expression is reduced in BC tissue and cells, and BC patients with lower NUDT21 expression have shorter overall and recurrent-free survival than patients with higher NUDT21 expression. NUDT21 ectopic expression or knockdown respectively profoundly inhibited or promoted the capacity of BC cells for proliferation, migration and invasion. We also identified a number of genes with shortened 3'UTRs through modulation of NUDT21 expression, and further characterized the NUDT21-regulated genes ANXA2 and LIMK2. We found NUDT21 modulates the expression of ANXA2 and LIMK2 in the Wnt/β-catenin and NF-κB signaling pathways.

Conclusions: These findings show NUDT21 plays a crucial role in BC progression, at least in part through ANXA2 and LIMK2 which act by alternative polyadenylation. NUDT21 may thus have potential as a diagnostic and therapeutic target in treatment of BC.

UPF1/SMG7-dependent microRNA-mediated gene regulation

The stability and quality of metazoan mRNAs are under microRNA (miRNA)-mediated and nonsense-mediated control. Although UPF1, a core mediator of nonsense-mediated mRNA decay (NMD), mediates the decay of target mRNA in a 3′UTR-length-dependent manner, the detailed mechanism remains unclear. Here, we suggest that 3′UTR-length-dependent mRNA decay is not mediated by nonsense mRNAs but rather by miRNAs that downregulate target mRNAs via Ago-associated UPF1/SMG7. Global analyses of mRNAs in response to UPF1 RNA interference in miRNA-deficient cells reveal that 3′UTR-length-dependent mRNA decay by UPF1 requires canonical miRNA targeting. The destabilization of miRNA targets is accomplished by the combination of Ago2 and UPF1/SMG7, which may recruit the CCR4-NOT deadenylase complex. Indeed, loss of the SMG7-deadenylase complex interaction increases the levels of transcripts regulated by UPF1-SMG7. This UPF1/SMG7-dependent miRNA-mediated mRNA decay pathway may enable miRNA targeting to become more predictable and expand the miRNA-mRNA regulatory network.

UPF1 mediates the decay of target mRNA in a 3′ untranslated region (UTR)-length-dependent manner. Here the authors reveal that the 3′UTR-length-dependent regulation of UPF1-dependent mRNA decay occurs through EJC-independent but miRNA-dependent regulation.

Transforming growth factor β1 alters the 3'-UTR of mRNA to promote lung fibrosis

Idiopathic pulmonary fibrosis (IPF) is a chronic disease characterized by the pathological remodeling of air sacs as a result of excessive accumulation of extracellular matrix (ECM) proteins, but the mechanism governing the robust protein expression is poorly understood. Our recent findings demonstrate that alternative polyadenylation (APA) caused by NUDT21 reduction is important for the increased expression of fibrotic mediators and ECM proteins in lung fibroblasts by shortening the 3'-untranslated regions (3'-UTRs) of mRNAs and stabilizing their transcripts, therefore activating pathological signaling pathways. Despite the importance of NUDT21 reduction in the regulation of fibrosis, the underlying mechanisms for the depletion are unknown. We demonstrate here that NUDT21 is depleted by TGFβ1. We found that miR203, which is increased in IPF, was induced by TGFβ1 to target the NUDT21 3'-UTR, thus depleting NUDT21 in human and mouse lung fibroblasts. TGFβ1-mediated NUDT21 reduction was attenuated by the miR203 inhibitor antagomiR203 in fibroblasts. TGFβ1 transgenic mice revealed that TGFβ1 down-regulates NUDT21 in fibroblasts in vivo Furthermore, TGFβ1 promoted differential APA of fibrotic genes, including FGF14, RICTOR, TMOD2, and UCP5, in association with increased protein expression. This unique differential APA signature was also observed in IPF fibroblasts. Altogether, our results identified TGFβ1 as an APA regulator through NUDT21 depletion amplifying pulmonary fibrosis.

The polyadenylation factor FIP1 is important for plant development and root responses to abiotic stresses

Root development and its response to environmental changes is crucial for whole plant adaptation. These responses include changes in transcript levels. Here, we show that the alternative polyadenylation (APA) of mRNA is important for root development and responses. Mutations in FIP1, a component of polyadenylation machinery, affects plant development, cell division and elongation, and response to different abiotic stresses. Salt treatment increases the amount of poly(A) site usage within the coding region and 5' untranslated regions (5'-UTRs), and the lack of FIP1 activity reduces the poly(A) site usage within these non-canonical sites. Gene ontology analyses of transcripts displaying APA in response to salt show an enrichment in ABA signaling, and in the response to stresses such as salt or cadmium (Cd), among others. Root growth assays show that fip1-2 is more tolerant to salt but is hypersensitive to ABA or Cd. Our data indicate that FIP1-mediated alternative polyadenylation is important for plant development and stress responses.

Genome-Wide Profiling of Polyadenylation Events in Maize Using High-Throughput Transcriptomic Sequences

Polyadenylation is an essential post-transcriptional modification of eukaryotic transcripts that plays critical role in transcript stability, localization, transport, and translational efficiency. About 70% genes in plants contain alternative polyadenylation (APA) sites. Despite availability of vast amount of sequencing data, to date, a comprehensive map of the polyadenylation events in maize is not available. Here, 9.48 billion RNA-Seq reads were analyzed to characterize 95,345 Poly(A) Clusters (PAC) in 23,705 (51%) maize genes. Of these, 76% were APA genes. However, most APA genes (55%) expressed a dominant PAC rather than favoring multiple PACs equally. The lincRNA genes with PACs were significantly longer in length than the genes without any PAC and about 48% genes had APA sites. Heterogeneity was observed in 52% of the PACs supporting the imprecise nature of the polyadenylation process. Genomic distribution revealed that the majority of the PACs (78%) were located in the genic regions. Unlike previous studies, large number of PACs were observed in the intergenic (n = 21,264), 5′-UTR (735), CDS (2,542), and the intronic regions (12,841). The CDS and introns with PACs were longer in length than without PACs, whereas intergenic PACs were more often associated with transcripts that lacked annotated 3′-UTRs. Nucleotide composition around PACs demonstrated AT-richness and the common upstream motif was AAUAAA, which is consistent with other plants. According to this study, only 2,830 genes still maintained the use of AAUAAA motif. This large-scale data provides useful insights about the gene expression regulation and could be utilized as evidence to validate the annotation of transcript ends.

TC3A: The Cancer 3' UTR Atlas

Widespread alternative polyadenylation (APA) occurs during enhanced cellular proliferation and transformation. Recently, we demonstrated that CFIm25-mediated 3′ UTR shortening through APA promotes glioblastoma tumor growth in vitro and in vivo, further underscoring its significance to tumorigenesis. Here, we report The Cancer 3′ UTR Atlas (TC3A), a comprehensive resource of APA usage for 10,537 tumors across 32 cancer types. These APA events represent potentially novel prognostic biomarkers and may uncover novel mechanisms for the regulation of cancer driver genes. TC3A is built on top of the now de facto standard cBioPortal. Therefore, the large community of existing cBioPortal users and clinical researchers will find TC3A familiar and immediately usable. TC3A is currently fully functional and freely available at http://tc3a.org.

A-to-I RNA editing - immune protector and transcriptome diversifier

Modifications of RNA affect its function and stability. RNA editing is unique among these modifications because it not only alters the cellular fate of RNA molecules but also alters their sequence relative to the genome. The most common type of RNA editing is A-to-I editing by double-stranded RNA-specific adenosine deaminase (ADAR) enzymes. Recent transcriptomic studies have identified a number of 'recoding' sites at which A-to-I editing results in non-synonymous substitutions in protein-coding sequences. Many of these recoding sites are conserved within (but not usually across) lineages, are under positive selection and have functional and evolutionary importance. However, systematic mapping of the editome across the animal kingdom has revealed that most A-to-I editing sites are located within mobile elements in non-coding parts of the genome. Editing of these non-coding sites is thought to have a critical role in protecting against activation of innate immunity by self-transcripts. Both recoding and non-coding events have implications for genome evolution and, when deregulated, may lead to disease. Finally, ADARs are now being adapted for RNA engineering purposes.

The RNA encoding the microtubule-associated protein tau has extensive structure that affects its biology

Tauopathies are neurodegenerative diseases that affect millions of people worldwide including those with Alzheimer’s disease. While many efforts have focused on understanding the role of tau protein in neurodegeneration, there has been little done to systematically analyze and study the structures within tau’s encoding RNA and their connection to disease pathology. Knowledge of RNA structure can provide insights into disease mechanisms and how to affect protein production for therapeutic benefit. Using computational methods based on thermodynamic stability and evolutionary conservation, we identified structures throughout the tau pre-mRNA, especially at exon-intron junctions and within the 5′ and 3′ untranslated regions (UTRs). In particular, structures were identified at twenty exon-intron junctions. The 5′ UTR contains one structured region, which lies within a known internal ribosome entry site. The 3′ UTR contains eight structured regions, including one that contains a polyadenylation signal. A series of functional experiments were carried out to assess the effects of mutations associated with mis-regulation of alternative splicing of exon 10 and to identify regions of the 3′ UTR that contain cis-regulatory elements. These studies defined novel structural regions within the mRNA that affect stability and pre-mRNA splicing and may lead to new therapeutic targets for treating tau-associated diseases.

Suboptimal RNA-RNA interaction limits U1 snRNP inhibition of canonical mRNA 3' processing

It is increasingly appreciated that U1 snRNP transcriptomically suppresses the usage of intronic polyadenylation site (PAS) of mRNAs, an outstanding question is why frequently used PASs are not suppressed. Here we found that U1 snRNP could be transiently associated with sequences upstream of actionable PASs in human cells, and RNA-RNA interaction might contribute to the association. By focusing on individual PAS, we showed that the stable assembly of U1 snRNP near PAS might be generally required for U1 inhibition of mRNA 3' processing. Therefore, actionable PASs that often lack optimal U1 snRNP docking site nearby is free from U1 inhibitory effect. Consistently, natural 5' splicing site (5'-SS) is moderately enriched ~250 nt upstream of intronic PASs whose usage is sensitive to functional knockdown of U1 snRNA. Collectively, our results provided an insight into how U1 snRNP selectively inhibits the usage of PASs in a cellular context, and supported a prevailing model that U1 snRNP scans pre-mRNA through RNA-RNA interaction to find a stable interaction site to exercise its function in pre-mRNA processing, including repressing the usage of cryptic PASs.

Modeling of Genome-Wide Polyadenylation Signals in Xenopus tropicalis

Alternative polyadenylation (APA) is an important post-transcriptional modification event to process messenger RNA (mRNA) for transcriptional termination, transport, and translation. In the present study, we characterized poly(A) signals in Xenopus tropicalis using 70,918 highly confident poly(A) sites derived from 16,511 protein-coding genes to understand their roles in the regulation of embryo development and gender difference. We examined potential factors, including the gene length, the number of introns in a gene, and the intron length, that may affect the prevalence of APA. We observed 12 prominent poly(A) signal patterns, which accounted for approximately 92% of total APA sites in Xenopus tropicalis. Among them, three patterns are specific to X. tropicalis, so they are absent in other animals such as humans or mice. We catalogued APA sites based on their genomic regions and developed a bioinformatics pipeline to identify over-represented signal patterns for each class. Then the schema of cis elements for APA sites in each genomic region was proposed. More importantly, APA usage is dramatically dynamic in embryos along five developmental stages and well-coordinated with the maternal-to-zygotic transition event. We used an entropy-based method to identify developmental stage-specific APA sites and identified significant signal patterns around specific sites and constitutive sites. We found that the APA frequency in different genomic regions varies with developmental stages and that those sites located in intron or coding sequence regions contribute most to the dynamics of gene expression during developmental stages. This study deciphers the characteristics and poly(A) signal patterns for both canonical APA sites and non-canonical APA sites across different developmental stages and gender dimorphisms in X. tropicalis, providing new insights into the dynamic regulation of distal and proximal APA.

A survey on identification and quantification of alternative polyadenylation sites from RNA-seq data

Alternative polyadenylation (APA) has been implicated to play an important role in post-transcriptional regulation by regulating mRNA abundance, stability, localization and translation, which contributes considerably to transcriptome diversity and gene expression regulation. RNA-seq has become a routine approach for transcriptome profiling, generating unprecedented data that could be used to identify and quantify APA site usage. A number of computational approaches for identifying APA sites and/or dynamic APA events from RNA-seq data have emerged in the literature, which provide valuable yet preliminary results that should be refined to yield credible guidelines for the scientific community. In this review, we provided a comprehensive overview of the status of currently available computational approaches. We also conducted objective benchmarking analysis using RNA-seq data sets from different species (human, mouse and Arabidopsis) and simulated data sets to present a systematic evaluation of 11 representative methods. Our benchmarking study showed that the overall performance of all tools investigated is moderate, reflecting that there is still lot of scope to improve the prediction of APA site or dynamic APA events from RNA-seq data. Particularly, prediction results from individual tools differ considerably, and only a limited number of predicted APA sites or genes are common among different tools. Accordingly, we attempted to give some advice on how to assess the reliability of the obtained results. We also proposed practical recommendations on the appropriate method applicable to diverse scenarios and discussed implications and future directions relevant to profiling APA from RNA-seq data.

Tissue-specific mechanisms of alternative polyadenylation: Testis, brain, and beyond (2018 update)

Alternative polyadenylation (APA) is how genes choose different sites for 3' end formation for mRNAs during transcription. APA often occurs in a tissue- or developmental stage-specific manner that can significantly affect gene activity by changing the protein product generated, the stability of the transcript, its localization within the cell, or its translatability. Despite the important regulatory effects that APA has on tissue-specific gene expression, only a few examples have been characterized mechanistically. In this 2018 update to our 2010 review, we examine mechanisms for the control of APA and update our understanding of the older mechanisms since 2010. We once postulated the existence of tissue-specific factors in APA. However, while a few tissue-specific polyadenylation factors are known, the emerging conclusion is that the majority of APA is accomplished by altering levels of core polyadenylation proteins. Examples of those core proteins include CSTF2, CPSF1, and subunits of mammalian cleavage factor I. But despite support for these mechanisms, no one has yet documented any of these proteins changing in either a tissue-specific or developmental manner. Given the profound effect that APA can have on gene expression and human health, improved understanding of tissue-specific APA could lead to numerous advances in gene activity control. This article is categorized under: RNA Processing > 3' End Processing RNA in Disease and Development > RNA in Development.

A Deep Neural Network for Predicting and Engineering Alternative Polyadenylation

Alternative polyadenylation (APA) is a major driver of transcriptome diversity in human cells. Here, we use deep learning to predict APA from DNA sequence alone. We trained our model (APARENT, APA REgression NeT) on isoform expression data from over 3 million APA reporters. APARENT's predictions are highly accurate when tasked with inferring APA in synthetic and human 3'UTRs. Visualizing features learned across all network layers reveals that APARENT recognizes sequence motifs known to recruit APA regulators, discovers previously unknown sequence determinants of 3' end processing, and integrates these features into a comprehensive, interpretable, cis-regulatory code. We apply APARENT to forward engineer functional polyadenylation signals with precisely defined cleavage position and isoform usage and validate predictions experimentally. Finally, we use APARENT to quantify the impact of genetic variants on APA. Our approach detects pathogenic variants in a wide range of disease contexts, expanding our understanding of the genetic origins of disease.

APAtrap: identification and quantification of alternative polyadenylation sites from RNA-seq data

Motivation

Alternative polyadenylation (APA) has been increasingly recognized as a crucial mechanism that contributes to transcriptome diversity and gene expression regulation. As RNA-seq has become a routine protocol for transcriptome analysis, it is of great interest to leverage such unprecedented collection of RNA-seq data by new computational methods to extract and quantify APA dynamics in these transcriptomes. However, research progress in this area has been relatively limited. Conventional methods rely on either transcript assembly to determine transcript 3' ends or annotated poly(A) sites. Moreover, they can neither identify more than two poly(A) sites in a gene nor detect dynamic APA site usage considering more than two poly(A) sites.

Results

We developed an approach called APAtrap based on the mean squared error model to identify and quantify APA sites from RNA-seq data. APAtrap is capable of identifying novel 3' UTRs and 3' UTR extensions, which contributes to locating potential poly(A) sites in previously overlooked regions and improving genome annotations. APAtrap also aims to tally all potential poly(A) sites and detect genes with differential APA site usages between conditions. Extensive comparisons of APAtrap with two other latest methods, ChangePoint and DaPars, using various RNA-seq datasets from simulation studies, human and Arabidopsis demonstrate the efficacy and flexibility of APAtrap for any organisms with an annotated genome.

Availability and implementation

Freely available for download at https://apatrap.sourceforge.io.

Contact

liqq@xmu.edu.cn or xhuister@xmu.edu.cn.

Supplementary information

Supplementary data are available at Bioinformatics online.

CFIm25 and alternative polyadenylation: Conflicting roles in cancer

Alternative polyadenylation (APA) is now widely recognized to regulate gene expression. APA is an RNA-processing mechanism that generates distinct 3' termini on mRNAs, producing mRNA isoforms. Different factors influence the initiation and development of this process. CFIm25 (among others) is a cleavage and polyadenylation factor that plays a key role in the regulation of APA. Shortening of the 3'UTRs on mRNAs leads to enhanced cellular proliferation and tumorigenicity. One reason may be the up-regulation of growth promoting factors, such as Cyclin D1. Different studies have reported a dual role of CFIm25 in cancer (both oncogenic and tumor suppressor). microRNAs (miRNAs) may be involved in CFIm25 function as well as competing endogenous RNAs (ceRNAs). The present review focuses on the role of CFIm25 in cancer, cancer treatment, and possible involvement in other human diseases. We highlight the involvement of miRNAs and ceRNAs in the function of CFIm25 to affect gene expression. The lack of understanding of the mechanisms and regulation of CFIm25 and APA has underscored the need for further research regarding their role in cancer and other diseases.

Cross-Regulation between TDP-43 and Paraspeckles Promotes Pluripotency-Differentiation Transition

RNA-binding proteins (RBPs) and long non-coding RNAs (lncRNAs) are key regulators of gene expression, but their joint functions in coordinating cell fate decisions are poorly understood. Here we show that the expression and activity of the RBP TDP-43 and the long isoform of the lncRNA Neat1, the scaffold of the nuclear compartment "paraspeckles," are reciprocal in pluripotent and differentiated cells because of their cross-regulation. In pluripotent cells, TDP-43 represses the formation of paraspeckles by enhancing the polyadenylated short isoform of Neat1. TDP-43 also promotes pluripotency by regulating alternative polyadenylation of transcripts encoding pluripotency factors, including Sox2, which partially protects its 3' UTR from miR-21-mediated degradation. Conversely, paraspeckles sequester TDP-43 and other RBPs from mRNAs and promote exit from pluripotency and embryonic patterning in the mouse. We demonstrate that cross-regulation between TDP-43 and Neat1 is essential for their efficient regulation of a broad network of genes and, therefore, of pluripotency and differentiation.

Nuclei of HeLa cells interactomes unravel a network of ghost proteins involved in proteins translation

Ghost proteins are issued from alternative Open Reading Frames (ORFs) and are missing a genome annotation. Indeed, historical filters applied for the detection of putative translated ORFs led to a wrong classification of transcripts considered as non-coding although translated proteins can be detected by proteomics. This Ghost (also called Alternative) proteome was neglected, and one major issue is to identify the implication of the Ghost proteins in the biological processes. In this context, we aimed to identify the protein-protein interactions (PPIs) of the Ghost proteins. For that, we re-explored a cross-link MS study performed on nuclei of HeLa cells using cross-linking mass spectrometry (XL-MS) associated with the HaltOrf database. Among 1679 cross-link interactions identified, 292 are involving Ghost Proteins. Forty-Four of these Ghost proteins are found to interact with 7 Reference proteins related to ribonucleoproteins, ribosome subunits and zinc finger proteins network. We, thus, have focused our attention on the heterotrimer between the RE/poly(U)-binding/degradation factor 1 (AUF1), the Ribosomal protein 10 (RPL10) and AltATAD2. Using I-Tasser software we performed docking models from which we could suggest the attachment of AUF1 on the external part of RPL10 and the interaction of AltATAD2 on the RPL10 region interacting with 5S ribosomal RNA as a mechanism of regulation of the ribosome. Taken together, these results reveal the importance of Ghost Proteins within known protein interaction networks.

Alternative polyadenylation of single cells delineates cell types and serves as a prognostic marker in early stage breast cancer

Alternative polyadenylation (APA) in 3’ untranslated regions (3’ UTR) plays an important role in regulating transcript abundance, localization, and interaction with microRNAs. Length-variation of 3’UTRs by APA contributes to efficient proliferation of cancer cells. In this study, we investigated APA in single cancer cells and tumor microenvironment cells to understand the physiological implication of APA in different cell types. We analyzed APA patterns and the expression level of genes from the 515 single-cell RNA sequencing (scRNA-seq) dataset from 11 breast cancer patients. Although the overall 3’UTR length of individual genes was distributed equally in tumor and non-tumor cells, we found a differential pattern of polyadenylation in gene sets between tumor and non-tumor cells. In addition, we found a differential pattern of APA across tumor types using scRNA-seq data from 3 glioblastoma patients and 1 renal cell carcinoma patients. In detail, 1,176 gene sets and 53 genes showed the distinct pattern of 3’UTR shortening and over-expression as signatures for five cell types including B lymphocytes, T lymphocytes, myeloid cells, stromal cells, and breast cancer cells. Functional categories of gene sets for cellular proliferation demonstrated concordant regulation of APA and gene expression specific to cell types. The expression of APA genes in breast cancer was significantly correlated with the clinical outcome of earlier stage breast cancer patients. We identified cell type-specific APA in single cells, which allows the identification of cell types based on 3’UTR length variation in combination with gene expression. Specifically, an immune-specific APA signature in breast cancer could be utilized as a prognostic marker of early stage breast cancer.

A PolH Transcript with a Short 3'UTR Enhances PolH Expression and Mediates Cisplatin Resistance

Platinum-based anticancer drugs are widely used as a first-line drug for cancers, such as non-small cell lung carcinoma (NSCLC) and bladder cancer. However, the efficacy is limited due to intrinsic or acquired resistance to these drugs. DNA polymerase eta (PolH, Polη) belongs to the Y-family of DNA polymerases and mediates DNA translesion synthesis, a major mechanism for DNA damage tolerance. Here, we showed that a high level of PolH is associated with cisplatin resistance in lung and bladder cancer. Consistent with this, loss of PolH markedly attenuates cisplatin resistance in both cisplatin-sensitive and cisplatin-resistant lung cancer cells. Interestingly, we found that due to the presence of multiple polyadenylation sites, alternative polyadenylation (APA) produces three major PolH transcripts with various lengths of 3'untranslated region (3'UTR; 427-/2516-/6245-nt). We showed that the short PolH transcript with 427-nt 3'UTR is responsible for high expression of PolH in various cisplatin-resistant lung and bladder cancer cell lines. Importantly, loss of the short PolH transcript significantly sensitizes cancer cells to cisplatin treatment. Moreover, we found that miR-619 selectively inhibits the ability of the long PolH transcript with 6245-nt 3'UTR to produce PolH protein and, subsequently, PolH-dependent cell growth. Together, our data suggest that PolH expression is controlled by APA and that the short PolH transcript produced by APA can escape miR-619-mediated repression and, subsequently, confers PolH-mediated cisplatin resistance. SIGNIFICANCE: A short PolH transcript produced by alternative polyadenylation escapes repression by miR-619 and confers resistance to cisplatin.

CEP135 isoform dysregulation promotes centrosome amplification in breast cancer cells

The centrosome, composed of two centrioles surrounded by pericentriolar material, is the cell's central microtubule-organizing center. Centrosome duplication is coupled with the cell cycle such that centrosomes duplicate once in S phase. Loss of such coupling produces supernumerary centrosomes, a condition called centrosome amplification (CA). CA promotes cell invasion and chromosome instability, two hallmarks of cancer. We examined the contribution of centriole overduplication to CA and the consequences for genomic stability in breast cancer cells. CEP135, a centriole assembly protein, is dysregulated in some breast cancers. We previously identified a short isoform of CEP135, CEP135mini, that represses centriole duplication. Here, we show that the relative level of full-length CEP135 (CEP135full) to CEP135mini (the CEP135full:mini ratio) is increased in breast cancer cell lines with high CA. Inducing expression of CEP135full in breast cancer cells increases the frequency of CA, multipolar spindles, anaphase-lagging chromosomes, and micronuclei. Conversely, inducing expression of CEP135mini reduces centrosome number. The differential expression of the CEP135 isoforms in vivo is generated by alternative polyadenylation. Directed genetic mutations near the CEP135mini alternative polyadenylation signal reduces the CEP135full:mini ratio and decreases CA. We conclude that dysregulation of CEP135 isoforms promotes centriole overduplication and contributes to chromosome segregation errors in breast cancer cells.

Polyadenylation sites and their characteristics in the genome of channel catfish (Ictalurus punctatus) as revealed by using RNA-Seq data

Polyadenylation plays important roles in gene expression regulation in eukaryotes, which typically involves cleavage and poly(A) tail addition at the polyadenylation site (PAS) of the pre-mature mRNA. Many eukaryotic genes contain more than one PASs, termed as alternative polyadenylation (APA). As a crucial post-transcriptional regulation, polyadenylation affects various aspects of RNA metabolism such as mRNA stability, translocation, and translation. However, polyadenylation has been rarely studied in teleosts. Here we conducted polyadenylation analysis in channel catfish, a commercially important aquaculture species around the world. Using RNA-Seq data, we identified 20,320 PASs which were classified into 14,500 clusters by merging adjacent PASs. Most of the PASs were found in 3' UTRs, followed by intron regions based on the annotation of channel catfish reference genome. No apparent difference in PAS distribution was observed between the sense and antisense strand of the channel catfish genome. The sequence analysis of nucleotide composition and motif around PASs yielded a highly similar profile among various organisms, suggesting the conservation and importance of polyadenylation in evolution. Using APA genes with more than two PASs, gene ontology enrichment revealed genes particularly involved in RNA binding. Reactome pathway analysis showed the enrichment of the innate immune system, especially neutrophil degranulation.

Transposable elements drive rapid phenotypic variation in Capsella rubella

The mechanisms underlying rapid adaptation to changing environments in species with reduced genetic variation, referred to as the “genetic paradox of invasion,” are unknown. We report that transposable elements (TEs) are highly enriched in the gene promoter regions of Capsella rubella compared with its outcrossing sister species Capsella grandiflora. We also show that a number of polymorphic TEs in C. rubella are associated with changes in gene expression. Frequent TE insertions at FLOWERING LOCUS C of C. rubella affect flowering-time variation, an important life history trait correlated with fitness. These results indicate that TE insertions drive rapid phenotypic variation, which could potentially help adapting to novel environments in species with limited genetic variation.

Rapid phenotypic changes in traits of adaptive significance are crucial for organisms to thrive in changing environments. How such phenotypic variation is achieved rapidly, despite limited genetic variation in species that experience a genetic bottleneck is unknown. Capsella rubella, an annual and inbreeding forb (Brassicaceae), is a great system for studying this basic question. Its distribution is wider than those of its congeneric species, despite an extreme genetic bottleneck event that severely diminished its genetic variation. Here, we demonstrate that transposable elements (TEs) are an important source of genetic variation that could account for its high phenotypic diversity. TEs are (i) highly enriched in C. rubella compared with its outcrossing sister species Capsella grandiflora, and (ii) 4.2% of polymorphic TEs in C. rubella are associated with variation in the expression levels of their adjacent genes. Furthermore, we show that frequent TE insertions at FLOWERING LOCUS C (FLC) in natural populations of C. rubella could explain 12.5% of the natural variation in flowering time, a key life history trait correlated with fitness and adaptation. In particular, we show that a recent TE insertion at the 3′ UTR of FLC affects mRNA stability, which results in reducing its steady-state expression levels, to promote the onset of flowering. Our results highlight that TE insertions can drive rapid phenotypic variation, which could potentially help with adaptation to changing environments in a species with limited standing genetic variation.

Detection of Differentially Expressed Cleavage Site Intervals Within 3' Untranslated Regions Using CSI-UTR Reveals Regulated Interaction Motifs

The length of untranslated regions at the 3' end of transcripts (3'UTRs) is regulated by alternate polyadenylation (APA). 3'UTRs contain regions that harbor binding motifs for regulatory molecules. However, the mechanisms that coordinate the 3'UTR length of specific groups of transcripts are not well-understood. We therefore developed a method, CSI-UTR, that models 3'UTR structure as tandem segments between functional alternative-polyadenylation sites (termed cleavage site intervals-CSIs). This approach facilitated (1) profiling of 3'UTR isoform expression changes and (2) statistical enrichment of putative regulatory motifs. CSI-UTR analysis is UTR-annotation independent and can interrogate legacy data generated from standard RNA-Seq libraries. CSI-UTR identified a set of CSIs in human and rodent transcriptomes. Analysis of RNA-Seq datasets from neural tissue identified differential expression events within 3'UTRs not detected by standard gene-based differential expression analyses. Further, in many instances 3'UTR and CDS from the same gene were regulated differently. This modulation of motifs for RNA-interacting molecules with potential condition-dependent and tissue-specific RNA binding partners near the polyA signal and CSI junction may play a mechanistic role in the specificity of alternative polyadenylation. Source code, CSI BED files and example datasets are available at: https://github.com/UofLBioinformatics/CSI-UTR.

Writing a wrong: Coupled RNA polymerase II transcription and RNA quality control

Processing and maturation of precursor RNA species is coupled to RNA polymerase II transcription. Co-transcriptional RNA processing helps to ensure efficient and proper capping, splicing, and 3' end processing of different RNA species to help ensure quality control of the transcriptome. Many improperly processed transcripts are not exported from the nucleus, are restricted to the site of transcription, and are in some cases degraded, which helps to limit any possibility of aberrant RNA causing harm to cellular health. These critical quality control pathways are regulated by the highly dynamic protein-protein interaction network at the site of transcription. Recent work has further revealed the extent to which the processes of transcription and RNA processing and quality control are integrated, and how critically their coupling relies upon the dynamic protein interactions that take place co-transcriptionally. This review focuses specifically on the intricate balance between 3' end processing and RNA decay during transcription termination. This article is categorized under: RNA Turnover and Surveillance > Turnover/Surveillance Mechanisms RNA Processing > 3' End Processing RNA Processing > Splicing Mechanisms RNA Processing > Capping and 5' End Modifications.

Comprehensive Characterization of Alternative Polyadenylation in Human Cancer

Background

Alternative polyadenylation (APA) is emerging as a major post-transcriptional mechanism for gene regulation, and dysregulation of APA contributes to several human diseases. However, the functional consequences of APA in human cancer are not fully understood. Particularly, there is no large-scale analysis in cancer cell lines.

Methods

We characterized the global APA profiles of 6398 patient samples across 17 cancer types from The Cancer Genome Atlas and 739 cancer cell lines from the Cancer Cell Line Encyclopedia. We built a linear regression model to explore the correlation between APA factors and APA events across different cancer types. We used Spearman correlation to assess the effects of APA events on drug sensitivity and the Wilcoxon rank-sum test or Cox proportional hazards model to identify clinically relevant APA events.

Results

We revealed a striking global 3'UTR shortening in cancer cell lines compared with tumor samples. Our analysis further suggested PABPN1 as the master regulator in regulating APA profile across different cancer types. Furthermore, we showed that APA events could affect drug sensitivity, especially of drugs targeting chromatin modifiers. Finally, we identified 1971 clinically relevant APA events, as well as alterations of APA in clinically actionable genes, suggesting that analysis of the complexity of APA profiles could have clinical utility.

Conclusions

Our study highlights important roles for APA in human cancer, including reshaping cellular pathways and regulating specific gene expression, exemplifying the complex interplay between APA and other biological processes and yielding new insights into the action mechanism of cancer drugs.

Cleavage factor 25 deregulation contributes to pulmonary fibrosis through alternative polyadenylation

Idiopathic pulmonary fibrosis (IPF) is a chronic and deadly disease with a poor prognosis and few treatment options. Pathological remodeling of the extracellular matrix (ECM) by myofibroblasts is a key factor that drives disease pathogenesis, although the underlying mechanisms remain unknown. Alternative polyadenylation (APA) has recently been shown to play a major role in cellular responses to stress by driving the expression of fibrotic factors and ECMs through altering microRNA sensitivity, but a connection to IPF has not been established. Here, we demonstrate that CFIm25, a global regulator of APA, is down-regulated in the lungs of patients with IPF and mice with pulmonary fibrosis, with its expression selectively reduced in alpha-smooth muscle actin (α-SMA) positive fibroblasts. Following the knockdown of CFIm25 in normal human lung fibroblasts, we identified 808 genes with shortened 3'UTRs, including those involved in the transforming growth factor-β signaling pathway, the Wnt signaling pathway, and cancer pathways. The expression of key pro-fibrotic factors can be suppressed by CFIm25 overexpression in IPF fibroblasts. Finally, we demonstrate that deletion of CFIm25 in fibroblasts or myofibroblast precursors using either the Col1a1 or the Foxd1 promoter enhances pulmonary fibrosis after bleomycin exposure in mice. Taken together, our results identified CFIm25 down-regulation as a novel mechanism to elevate pro-fibrotic gene expression in pulmonary fibrosis.

Variable 3'polyadenylation of Wheat yellow mosaic virus and its novel effects on translation and replication

BACKGROUND

Polyadenylation influences many aspects of mRNA as well as viral RNA. variable polyadenylation at the 3' end have been reported in RNA viruses. It is interesting to identify the characteristic and potential role of 3' polyadenylation of Wheat yellow mosaic virus (WYMV), which has been reported to contain two genomic RNAs with 3' poly(A) tails and caused severe disease on wheat in East Asia region.

METHODS

3' RACE was used to identify sequences of the 3' end in WYMV RNAs from naturally infected wheat by WYMV. In vitro translation assay was performed to analyze effect of UTRs of WYMV with or without 3'polyadenylation on translation. In vitro replication mediated by WYMV NIb protein were performed to evaluate effect of variable polyadenylation on replication.

RESULTS

Variable polyadenylation in WYMV RNAs was identified via 3' RACE. WYMV RNAs in naturally infected wheat in China simultaneously present with regions of long, short, or no adenylation at the 3' ends. The effects of variable polyadenylation on translation and replication of WYMV RNAs were evaluated. 5'UTR and 3'UTR of WYMV RNA1 or RNA2 synergistically enhanced the translation of the firefly luciferase (Fluc) gene in in vitro WGE system, whereas additional adenylates had an oppositive effect on this enhancement on translation mediated by UTRs of WYMV. Additional adenylates remarkably inhibited the synthesis of complementary strand from viral genome RNA during the in vitro replication mediated by WYMV NIb protein.

CONCLUSIONS

3' end of WYMV RNAs present variable polyadenylation even no polyadenylation. 3' polyadenylation have opposite effect on translation mediated by UTRs of WYMV RNA1 or RNA2. 3' polyadenylation have negative effect on minus-strand synthesis of WYMV RNA in vitro. Variable polyadenylation of WYMV RNAs may provide sufficient selection on the template for translation and replication.

It's all about the T: transcription termination in archaea

One of the most fundamental biological processes driving all life on earth is transcription. The, at first glance, relatively simple cycle is divided into three stages: initiation at the promoter site, elongation throughout the open reading frame, and finally termination and product release at the terminator. In all three processes, motifs of the template DNA and protein factors of the transcription machinery including the multisubunit polymerase itself as well as a broad range of associated transcription factors work together and mutually influence each other. Despite several decades of research, this interplay holds delicate mechanistic and structural details as well as interconnections yet to be explored. One of the surprising characteristics of archaeal biology is the use of eukaryotic-like information processing systems against a backdrop of a bacterial-like genome. Archaeal genomes usually comprise main chromosomes alongside chromosomal plasmids, and the genetic information is encoded in single transcriptional units as well as in multicistronic operons alike their bacterial counterparts. Moreover, archaeal genomes are densely packed and this necessitates a tight regulation of transcription and especially assured termination events in order to prevent read-through into downstream coding regions and the accumulation of antisense transcripts.

Alternative Splicing of MXD3 and Its Regulation of MXD3 Levels in Glioblastoma

The transcription factor MXD3 is an atypical member of the MYC/MAX/MXD transcriptional network and has been previously shown to be an important regulator of cell proliferation. MXD3 has been shown to be overexpressed and to be required for medulloblastoma and acute lymphoblastic leukemia cell proliferation. In this study we leveraged datasets from The Cancer Genome Atlas to examine MXD3 across several cancers. We find that MXD3 transcripts are significantly overexpressed in ~72% of the available datasets. The gene itself is not frequently altered, while the promoter appears to be hypomethylated. We examine the possibility that aberrant regulation of the MXD3 message is the cause of abnormal MXD3 expression across cancers. Specifically, we looked at MXD3 alternative splicing in glioblastoma multiforme (GBM) and find notable functional differences between the splice variants. The 3′UTR confers differential message stability. Furthermore, the different coding sequences lead to different protein stabilities and localizations. Altogether, these data extend our knowledge of MXD3 in the context of human cancers while characterizing a previously unstudied splice variant of MXD3.

Arabidopsis CRL4 Complexes: Surveying Chromatin States and Gene Expression

CULLIN4 (CUL4) RING ligase (CRL4) complexes contain a CUL4 scaffold protein, associated to RBX1 and to DDB1 proteins and have traditionally been associated to protein degradation events. Through DDB1, these complexes can associate with numerous DCAF proteins, which directly interact with specific targets promoting their ubiquitination and subsequent degradation by the proteasome. A characteristic feature of the majority of DCAF proteins that associate with DDB1 is the presence of the DWD motif. DWD-containing proteins sum up to 85 in the plant model species Arabidopsis. In the last decade, numerous Arabidopsis DWD proteins have been studied and their molecular functions uncovered. Independently of whether their association with CRL4 has been confirmed or not, DWD proteins are often found as components of additional multimeric protein complexes that play key roles in essential nuclear events. For most of them, the significance of their complex partnership is still unexplored. Here, we summarize recent findings involving both confirmed and putative CRL4-associated DCAF proteins in regulating nuclei architecture remodelling, DNA damage repair, histone post-translational modification, mRNA processing and export, and ribosome biogenesis, that definitely have an impact in gene expression and de novo protein synthesis. We hypothesized that, by maintaining accurate levels of regulatory proteins through targeted degradation and transcriptional control, CRL4 complexes help to surveil nuclear processes essential for plant development and survival.

The histone demethylase Jarid1b mediates angiotensin II-induced endothelial dysfunction by controlling the 3'UTR of soluble epoxide hydrolase

AIM

The histone demethylase Jarid1b limits gene expression by removing the active methyl mark from histone3 lysine4 at gene promoter regions. A vascular function of Jarid1b is unknown, but a vasoprotective function to inflammatory and hypertrophic stimuli, like angiotensin II (AngII) could be inferred. This hypothesis was tested using Jarid1b knockout mice and the inhibitor PBIT.

METHODS

Mice or aortic segments were treated with AngII to induce endothelial dysfunction. Aortae from WT and Jarid1b knockout were studied in organ chambers and endothelium-dependent dilator responses to acetylcholine and endothelium-independent responses to DetaNONOate were recorded after pre-constriction with phenylephrine in the presence or absence of the NO-synthase inhibitor nitro-L-arginine. Molecular mechanisms were investigated with chromatin immunoprecipitation, RNA-Seq, RNA-3'-adaptor-ligation, actinomycin D and RNA-immunoprecipitation.

RESULTS

Knockout or inhibition of Jarid1b prevented the development of endothelial dysfunction in response to AngII. This effect was not a consequence of altered nitrite oxide availability but accompanied by a loss of the inflammatory response to AngII. As Jarid1b mainly inhibits gene expression, an indirect effect should account for this observation. AngII induced the soluble epoxide hydrolase (sEH), which degrades anti-inflammatory lipids, and thus promotes inflammation. Knockout or inhibition of Jarid1b prevented the AngII-mediated sEH induction. Mechanistically, Jarid1b maintained the length of the 3'untranslated region of the sEH mRNA, thereby increasing its stability and thus sEH protein expression. Loss of Jarid1b activity therefore resulted in sEH mRNA destabilization.

CONCLUSION

Jarid1b contributes to the pro-inflammatory effects of AngII by stabilizing sEH expression. Jarid1b inhibition might be an option for future therapeutics against cardiovascular dysfunction.