Comprehensive polyadenylation site maps in yeast and human reveal pervasive alternative polyadenylation

The deubiquitylase USP9X controls ribosomal stalling

When a ribosome stalls during translation, it runs the risk of collision with a trailing ribosome. Such an encounter leads to the formation of a stable di-ribosome complex, which needs to be resolved by a dedicated machinery. The initial stalling and the subsequent resolution of di-ribosomal complexes requires activity of Makorin and ZNF598 ubiquitin E3 ligases, respectively, through ubiquitylation of the eS10 and uS10 subunits of the ribosome. We have developed a specific small-molecule inhibitor of the deubiquitylase USP9X. Proteomics analysis, following inhibitor treatment of HCT116 cells, confirms previous reports linking USP9X with centrosome-associated protein stability but also reveals a loss of Makorin 2 and ZNF598. We show that USP9X interacts with both these ubiquitin E3 ligases, regulating their abundance through the control of protein stability. In the absence of USP9X or following chemical inhibition of its catalytic activity, levels of Makorins and ZNF598 are diminished, and the ribosomal quality control pathway is impaired.

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

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

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

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

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

Non-coding Natural Antisense Transcripts: Analysis and Application

Non-coding natural antisense transcripts (ncNATs) are regulatory RNA sequences that are transcribed in the opposite direction to protein-coding or non-coding transcripts. These transcripts are implicated in a broad variety of biological and pathological processes, including tumorigenesis and oncogenic progression. With this complex field still in its infancy, annotations, expression profiling and functional characterisations of ncNATs are far less comprehensive than those for protein-coding genes, pointing out substantial gaps in the analysis and characterisation of these regulatory transcripts. In this review, we discuss ncNATs from an analysis perspective, in particular regarding the use of high-throughput sequencing strategies, such as RNA-sequencing, and summarize the unique challenges of investigating the antisense transcriptome. Finally, we elaborate on their potential as biomarkers and future targets for treatment, focusing on cancer.

Genetic and pharmacological evidence for kinetic competition between alternative poly(A) sites in yeast

Most eukaryotic mRNAs accommodate alternative sites of poly(A) addition in the 3’ untranslated region in order to regulate mRNA function. Here, we present a systematic analysis of 3’ end formation factors, which revealed 3’UTR lengthening in response to a loss of the core machinery, whereas a loss of the Sen1 helicase resulted in shorter 3’UTRs. We show that the anti-cancer drug cordycepin, 3’ deoxyadenosine, caused nucleotide accumulation and the usage of distal poly(A) sites. Mycophenolic acid, a drug which reduces GTP levels and impairs RNA polymerase II (RNAP II) transcription elongation, promoted the usage of proximal sites and reversed the effects of cordycepin on alternative polyadenylation. Moreover, cordycepin-mediated usage of distal sites was associated with a permissive chromatin template and was suppressed in the presence of an rpb1 mutation, which slows RNAP II elongation rate. We propose that alternative polyadenylation is governed by temporal coordination of RNAP II transcription and 3’ end processing and controlled by the availability of 3’ end factors, nucleotide levels and chromatin landscape.

LABRAT reveals association of alternative polyadenylation with transcript localization, RNA binding protein expression, transcription speed, and cancer survival

BACKGROUND

The sequence content of the 3' UTRs of many mRNA transcripts is regulated through alternative polyadenylation (APA). The study of this process using RNAseq data, though, has been historically challenging.

RESULTS

To combat this problem, we developed LABRAT, an APA isoform quantification method. LABRAT takes advantage of newly developed transcriptome quantification techniques to accurately determine relative APA site usage and how it varies across conditions. Using LABRAT, we found consistent relationships between gene-distal APA and subcellular RNA localization in multiple cell types. We also observed connections between transcription speed and APA site choice as well as tumor-specific transcriptome-wide shifts in APA isoform abundance in hundreds of patient-derived tumor samples that were associated with patient prognosis. We investigated the effects of APA on transcript expression and found a weak overall relationship, although many individual genes showed strong correlations between relative APA isoform abundance and overall gene expression. We interrogated the roles of 191 RNA-binding proteins in the regulation of APA isoforms, finding that dozens promote broad, directional shifts in relative APA isoform abundance both in vitro and in patient-derived samples. Finally, we find that APA site shifts in the two classes of APA, tandem UTRs and alternative last exons, are strongly correlated across many contexts, suggesting that they are coregulated.

CONCLUSIONS

We conclude that LABRAT has the ability to accurately quantify APA isoform ratios from RNAseq data across a variety of sample types. Further, LABRAT is able to derive biologically meaningful insights that connect APA isoform regulation to cellular and molecular phenotypes.

Set2-mediated H3K36 methylation states redundantly repress the production of antisense transcripts: role in transcription regulation

Methyltransferase Set2‐mediated methylation of histone H3 lysine 36 (H3K36), which involves the addition of up to three methyl groups at this site, has been demonstrated to function in many chromatin‐coupled events. The methylation of H3K36 is known to recruit different chromatin effector proteins, affecting transcription, mRNA splicing and DNA repair. In this study, we engineered two yeast set2 mutants that lack H3K36 mono/dimethylation (H3K36me1/2) and trimethylation (H3K36me3), respectively, and characterized their roles in the production of antisense transcripts under nutrient‐rich conditions. Using our new bioinformatics identification pipeline analysis, we are able to identify a larger number of antisense transcripts in set2∆ cells than has been published previously. We further show that H3K36me1/2 or H3K36me3 redundantly repressed the production of antisense transcripts. Moreover, gene ontology (GO) analysis implies that H3K36me3‐mediated antisense transcription might play a role in DNA replication and DNA damage repair, which is independent of regulation of the corresponding sense gene expression. Overall, our results validate a coregulatory mechanism of different H3K36 methylation states, particularly in the repression of antisense transcription.

TRENDseq-A highly multiplexed high throughput RNA 3' end sequencing for mapping alternative polyadenylation

Alternative polyadenylation (APA) is a widespread and highly dynamic mechanism of gene regulation. It affects more than 70% of all genes, resulting in transcript isoforms with distinct 3' end termini. APA thereby considerably expands the diversity of the transcriptome 3' end (TREND). This leads to mRNA isoforms with profoundly different physiological effects, by affecting protein output, production of distinct protein isoforms, or modulating protein localization. APA is globally regulated in various conditions, including developmental and adaptive programs. Since perturbations of APA can disrupt biological processes, ultimately resulting in most devastating disorders, querying the APA landscape is crucial to decipher underlying mechanisms, resulting consequences and potential diagnostic and therapeutic implications. Here we provide a detailed step-by-step protocol for TRENDseq, a method for transcriptome-wide high-throughput sequencing of polyadenylated RNA 3' ends in a highly multiplexed fashion. TRENDseq exploits linear amplification of the starting material to improve sensitivity while significantly reducing the amount of input material. It thereby represents a powerful tool to study APA in numerous experimental set-ups and/or limited human samples in a highly multiplexed and reproducible manner.

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

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

Detecting and Rescuing Stalled Ribosomes

Ribosomes that stall inappropriately during protein synthesis harbor proteotoxic components linked to cellular stress and neurodegenerative diseases. Molecular mechanisms that rescue stalled ribosomes must selectively detect rare aberrant translational complexes and process the heterogeneous components. Ribosome-associated quality control pathways eliminate problematic messenger RNAs and nascent proteins on stalled translational complexes. In addition, recent studies have uncovered general principles of stall recognition upstream of quality control pathways and fail-safe mechanisms that ensure nascent proteome integrity. Here, we discuss developments in our mechanistic understanding of the detection and rescue of stalled ribosomal complexes in eukaryotes.

HIV-1 Natural Antisense Transcription and Its Role in Viral Persistence

Natural antisense transcripts (NATs) represent a class of RNA molecules that are transcribed from the opposite strand of a protein-coding gene, and that have the ability to regulate the expression of their cognate protein-coding gene via multiple mechanisms. NATs have been described in many prokaryotic and eukaryotic systems, as well as in the viruses that infect them. The human immunodeficiency virus (HIV-1) is no exception, and produces one or more NAT from a promoter within the 3’ long terminal repeat. HIV-1 antisense transcripts have been the focus of several studies spanning over 30 years. However, a complete appreciation of the role that these transcripts play in the virus lifecycle is still lacking. In this review, we cover the current knowledge about HIV-1 NATs, discuss some of the questions that are still open and identify possible areas of future research.

Computational analysis of alternative polyadenylation from standard RNA-seq and single-cell RNA-seq data

Alternative polyadenylation (APA) is a major mechanism of post-transcriptional regulation in various cellular processes including cell proliferation and differentiation. Since conventional APA profiling methods have not been widely adopted, global APA studies are very limited. In this chapter, we summarize current computational methods for analyzing APA in standard RNA-seq and scRNA-seq data and describe two state-of-the-art bioinformatic algorithms DaPars and scDaPars in detail. The bioinformatic pipelines for both DaPars and scDaPars are presented and the application of both algorithms are highlighted.

Inter-dependent Centrosomal Co-localization of the cen and ik2 cis-Natural Antisense mRNAs in Drosophila

Overlapping genes are prevalent in most genomes, but the extent to which this organization influences regulatory events operating at the post-transcriptional level remains unclear. Studying the cen and ik2 genes of Drosophila melanogaster, which are convergently transcribed as cis-natural antisense transcripts (cis-NATs) with overlapping 3' UTRs, we found that their encoded mRNAs strikingly co-localize to centrosomes. These transcripts physically interact in a 3' UTR-dependent manner, and the targeting of ik2 requires its 3' UTR sequence and the presence of cen mRNA, which serves as the main driver of centrosomal co-localization. The cen transcript undergoes localized translation in proximity to centrosomes, and its localization is perturbed by polysome-disrupting drugs. By interrogating global fractionation-sequencing datasets generated from Drosophila and human cellular models, we find that RNAs expressed as cis-NATs tend to co-localize to specific subcellular fractions. This work suggests that post-transcriptional interactions between RNAs with complementary sequences can dictate their localization fate in the cytoplasm.

Investigation of the putative role of antisense transcripts as regulators of sense transcripts by correlation analysis of sense-antisense pairs in colorectal cancers

Antisense transcription occurs widely more expected than when it was first identified in bacteria in the 1980s. However, the functional relevance of antisense transcripts in transcription remains controversial. Here, we investigated the putative role of antisense transcripts in regulating their corresponding sense transcripts by analyzing changes in correlative relationships between sense-antisense pairs under tumor and normal conditions. A total of 3469 sense-antisense gene pairs (SAGPs) downloaded from BioMart mapped to a list of sense and antisense genes in RNA-seq data derived from 80 paired colorectal cancer (CRC) samples were analyzed. As a result, cancer-related genes were significantly enriched in the significantly correlated SAGPs (SCPs). Differentially expressed genes estimated between normal and tumor conditions were also significantly more enriched in SCPs than in non-SCPs. Interestingly, using differential correlation analysis, we found that tumor samples had a significantly larger density of genes with higher correlation coefficients than normal samples, as verified by various cancer transcriptomes from The Cancer Genome Atlas (TCGA). Moreover, we found that the magnitude of the correlation between SAGPs could distinguish poor prognostic CRCs from good prognostic CRCs, showing that correlation coefficients between the SAGPs of CRCs with a poor prognosis were significantly stronger than CRCs with a good prognosis. Consistent with this finding, the Cox proportion hazards model showed that the survival rates were significantly different between patients with high and low expression of genes in the SCPs. All these results strongly support the idea that antisense transcripts are important regulators of their corresponding sense transcripts.

Genome-wide identification and functional characterization of natural antisense transcripts in Salvia miltiorrhiza

Salvia miltiorrhiza is one of the most widely used traditional medicines. Natural antisense transcripts (NATs) are a class of long noncoding RNAs that can regulate gene expression. Here, we identified 812 NATs, including 168 cis-NATs and 644 trans-NATs from twelve root, flower, and leaf samples of S. miltiorrhiza using RNA-seq. The expression profiles for 41 of 50 NATs and their sense transcripts (STs) obtained from RNA-Seq were validated using qRT-PCR. The expression profiles of 17 NATs positively correlated with their STs. GO and KEGG pathway analyses mapped the STs for cis-NATs to pathways for biosynthesis of secondary metabolites. We characterized four NATs in detail, including NAT0001, NAT0002, NAT0004, and NAT00023. Their STs are kaurene synthase-like 1 and the homologs of UDP-glucose flavonoid 3-O-glucosyltransferase 6, UDP-glycosyltransferase 90A1, and beta-glucosidase 40, respectively. The first gene is involved in the biosynthesis of bioactive tanshinones, the next two are involved in anthocyanin biosynthesis, whereas the last is involved in phenylpropanoid biosynthesis. Besides, we found seven STs that are potential targets of miRNAs. And we found two miRNAs including miR156a and miR7208, might originate from NATs, NAT0112 and NAT0086. The results suggest that S. miltiorrhiza NATs might interact with STs, produce miRNAs, and be regulated by miRNAs. They potentially play significant regulatory roles in the biosynthesis of bioactive compounds.

TREND-DB-a transcriptome-wide atlas of the dynamic landscape of alternative polyadenylation

Alternative polyadenylation (APA) profoundly expands the transcriptome complexity. Perturbations of APA can disrupt biological processes, ultimately resulting in devastating disorders. A major challenge in identifying mechanisms and consequences of APA (and its perturbations) lies in the complexity of RNA 3′ end processing, involving poorly conserved RNA motifs and multi-component complexes consisting of far more than 50 proteins. This is further complicated in that RNA 3′ end maturation is closely linked to transcription, RNA processing and even epigenetic (histone/DNA/RNA) modifications. Here, we present TREND-DB (http://shiny.imbei.uni-mainz.de:3838/trend-db), a resource cataloging the dynamic landscape of APA after depletion of >170 proteins involved in various facets of transcriptional, co- and post-transcriptional gene regulation, epigenetic modifications and further processes. TREND-DB visualizes the dynamics of transcriptome 3′ end diversification (TREND) in a highly interactive manner; it provides a global APA network map and allows interrogating genes affected by specific APA-regulators and vice versa. It also permits condition-specific functional enrichment analyses of APA-affected genes, which suggest wide biological and clinical relevance across all RNAi conditions. The implementation of the UCSC Genome Browser provides additional customizable layers of gene regulation accounting for individual transcript isoforms (e.g. epigenetics, miRNA-binding sites and RNA-binding proteins). TREND-DB thereby fosters disentangling the role of APA for various biological programs, including potential disease mechanisms, and helps identify their diagnostic and therapeutic potential.

DMP-1 promoter-associated antisense strand non-coding RNA, panRNA-DMP-1, physically associates with EGFR to repress EGF-induced squamous cell carcinoma migration

Accumulating evidence suggests that specific non-coding RNAs exist in many types of malignant tissues, and are involved in cancer invasion and metastasis. However, little is known about the precise roles of non-coding RNAs in squamous cell carcinoma (SQCC) invasion and migration. Recently, the dentin matrix protein-1 (DMP-1) gene locus was identified as a transcriptionally active site in squamous cell carcinoma (SQCC) tissue and cells. However, it is unclear whether RNA associated with cell migration exist at the DMP-1 gene locus in SQCC cells. We identified a novel promoter-associated non-coding RNA in the antisense strand of DMP-1 gene locus, promoter-associated non-coding RNA (panRNA)-DMP-1, by the RACE method in SQCC cells and tissues, and characterized the functions of panRNA-DMP-1 in EGF-driven SQCC cell migration. The inhibition of endogenous panRNA-DMP-1 expression by specific siRNAs and exogenous over-expression of panRNA-DMP-1 resulted in increased and suppressed cellular migration toward EGF in SQCC cells, respectively, and nuclear expression of panRNA-DMP-1 was induced by EGF stimulation. Mechanistically, suppression of panRNA-DMP-1 expression increased EGFR nuclear localization upon EGF treatment and nuclear panRNA-DMP-1 physically interacted with EGFR, which was confirmed by RNA immunoprecipitation assay using a bacteriophage-delivered PP7 RNA labeling system. Furthermore, co-immunoprecipitation assay revealed that suppression of panRNA-DMP-1 stabilized EGFR interaction with STAT3, a known co-transcription factors of EGFR, to induce migratory properties in many cancer cells. Based on these findings, panRNA-DMP-1 is an EGFR-associating RNA that inhibits the EGF-induced migratory properties of SQCC possibly by regulating EGFR nuclear localization and EGFR binding to STAT3.

Functional roles of antisense enhancer RNA for promoting prostate cancer progression

Rationale: Enhancer RNA (eRNA) bi-directionally expresses from enhancer region and sense eRNA regulates adjacent mRNA in cis and in trans. However, it has remained unclear whether antisense eRNAs in different direction are functional or merely a reflection of enhancer activation. Methods: Strand-specific, ribosome-minus RNA sequencing (RNA-seq) were performed in AR positive prostate cancer cells. RNA-seq, GRO-seq, ChIP-seq, 4C-seq and DNA-methylation-seq that published in our and other labs were re-analyzed to define bi-directional enhancer RNA and DNA methylation regions. Molecular mechanisms were demonstrated by 3C, ChIP, ChIRP, CLIP, RT-PCR and western blot assays. The biological functions of antisense-eRNA were assessed using mice xenograft model and RT-PCR analysis in human tissues. Results: In this study, we identified that antisense eRNA was regulated by androgen receptor (AR) activity in prostate cancer cells. Antisense eRNA negatively regulated antisense ncRNA in AR-related target genes' loci, through recruiting DNMT1 on the antisense enhancer in the gene-ending regions and elevating DNA methylation. Importantly, the chromatin exhibited a double looping manner that facilitated sense-eRNA to promoter and antisense-eRNA to gene-ending region in cis. Depletion of antisense eRNA impaired its neighbor mRNA expression, cancer growth and invasion. The expressions of antisense eRNA were correlated with biochemical recurrence and clinical marker PSA's levels in patients' tissues. Conclusions: The findings indicated that antisense eRNA was a functional RNA and may be a novel target that when suppressed improved prostate cancer therapy and diagnosis. New chromatin interaction among enhancer, promoter and gene-ending region might provide new insight into the spatiotemporal mechanism of the gene transcription and acting of bi-directional eRNAs.

Ribosome-associated quality control of membrane proteins at the endoplasmic reticulum

Protein synthesis is an energetically costly, complex and risky process. Aberrant protein biogenesis can result in cellular toxicity and disease, with membrane-embedded proteins being particularly challenging for the cell. In order to protect the cell from consequences of defects in membrane proteins, quality control systems act to maintain protein homeostasis. The majority of these pathways act post-translationally; however, recent evidence reveals that membrane proteins are also subject to co-translational quality control during their synthesis in the endoplasmic reticulum (ER). This newly identified quality control pathway employs components of the cytosolic ribosome-associated quality control (RQC) machinery but differs from canonical RQC in that it responds to biogenesis state of the substrate rather than mRNA aberrations. This ER-associated RQC (ER-RQC) is sensitive to membrane protein misfolding and malfunctions in the ER insertion machinery. In this Review, we discuss the advantages of co-translational quality control of membrane proteins, as well as potential mechanisms of substrate recognition and degradation. Finally, we discuss some outstanding questions concerning future studies of ER-RQC of membrane proteins.

Transcriptome Analysis of the Liver and Muscle Tissues of Black Carp (Mylopharyngodon piceus) of Different Growth Rates

In this study, we used RNA-seq to analyze the muscle and liver tissues of black carps (Mylopharyngodon piceus) of different growth rates from the same batch to evaluate their growth traits. We have two groups; they are the black carp group with fast-growth rate and the slow-growth rate. A total of 23,132 genes were enriched in the Gene Ontology analysis, and 285 related pathways were found in the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. The KEGG pathway analysis showed significant differences in the expression of some genes involved in growth- and development-related metabolic pathways such as the FoxO signaling pathway, p53 signaling pathway, PI3K-Akt signaling pathway, apoptosis, TGF-β signaling pathway, and insulin signaling pathway. The numbers of differentially expressed genes in muscle and liver are 1913 and 1775. Nine of the differently expressed genes involved in the different growth traits and accuracy of the transcriptome data were validated using quantitative real-time PCR. We found that the expression levels of some growth-related genes were significantly higher in the fast-growth rate black carps than in the slow-growth rate black carps. The large number of transcriptome sequences obtained in this study has enriched the black carp gene resources, and the obtained differentially expressed genes and related pathway analysis provide valuable information for understanding the growth traits of the black carp.

Gene expression dynamics are a proxy for selective pressures on alternatively polyadenylated isoforms

Alternative polyadenylation (APA) produces isoforms with distinct 3′-ends, yet their functional differences remain largely unknown. Here, we introduce the APA-seq method to detect the expression levels of APA isoforms from 3′-end RNA-Seq data by exploiting both paired-end reads for gene isoform identification and quantification. We detected the expression levels of APA isoforms in individual Caenorhabditis elegans embryos at different stages throughout embryogenesis. Examining the correlation between the temporal profiles of isoforms led us to distinguish two classes of genes: those with highly correlated isoforms (HCI) and those with lowly correlated isoforms (LCI) across time. We hypothesized that variants with similar expression profiles may be the product of biological noise, while the LCI variants may be under tighter selection and consequently their distinct 3′ UTR isoforms are more likely to have functional consequences. Supporting this notion, we found that LCI genes have significantly more miRNA binding sites, more correlated expression profiles with those of their targeting miRNAs and a relative lack of correspondence between their transcription and protein abundances. Collectively, our results suggest that a lack of coherence among the regulation of 3′ UTR isoforms is a proxy for selective pressures acting upon APA usage and consequently for their functional relevance.

The transcriptional elongation rate regulates alternative polyadenylation in yeast

Yeast cells undergoing the diauxic response show a striking upstream shift in poly(A) site utilization, with increased use of ORF-proximal poly(A) sites resulting in shorter 3' mRNA isoforms for most genes. This altered poly(A) pattern is extremely similar to that observed in cells containing Pol II derivatives with slow elongation rates. Conversely, cells containing derivatives with fast elongation rates show a subtle downstream shift in poly(A) sites. Polyadenylation patterns of many genes are sensitive to both fast and slow elongation rates, and a global shift of poly(A) utilization is strongly linked to increased purine content of sequences flanking poly(A) sites. Pol II processivity is impaired in diauxic cells, but strains with reduced processivity and normal Pol II elongation rates have normal polyadenylation profiles. Thus, Pol II elongation speed is important for poly(A) site selection and for regulating poly(A) patterns in response to environmental conditions.

The nucleosome DNA entry-exit site is important for transcription termination and prevention of pervasive transcription

Compared to other stages in the RNA polymerase II transcription cycle, the role of chromatin in transcription termination is poorly understood. We performed a genetic screen in Saccharomyces cerevisiae to identify histone mutants that exhibit transcriptional readthrough of terminators. Amino acid substitutions identified by the screen map to the nucleosome DNA entry-exit site. The strongest H3 mutants revealed widespread genomic changes, including increased sense-strand transcription upstream and downstream of genes, increased antisense transcription overlapping gene bodies, and reduced nucleosome occupancy particularly at the 3' ends of genes. Replacement of the native sequence downstream of a gene with a sequence that increases nucleosome occupancy in vivo reduced readthrough transcription and suppressed the effect of a DNA entry-exit site substitution. Our results suggest that nucleosomes can facilitate termination by serving as a barrier to transcription and highlight the importance of the DNA entry-exit site in broadly maintaining the integrity of the transcriptome.

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

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

Inhibition of DUX4 expression with antisense LNA gapmers as a therapy for facioscapulohumeral muscular dystrophy

Facioscapulohumeral muscular dystrophy (FSHD), characterized by progressive muscle weakness and deterioration, is genetically linked to aberrant expression of DUX4 in muscle. DUX4, in its full-length form, is cytotoxic in nongermline tissues. Here, we designed locked nucleic acid (LNA) gapmer antisense oligonucleotides (AOs) to knock down DUX4 in immortalized FSHD myoblasts and the FLExDUX4 FSHD mouse model. Using a screening method capable of reliably evaluating the knockdown efficiency of LNA gapmers against endogenous DUX4 messenger RNA in vitro, we demonstrate that several designed LNA gapmers selectively and effectively reduced DUX4 expression with nearly complete knockdown. We also found potential functional benefits of AOs on muscle fusion and structure in vitro. Finally, we show that one of the LNA gapmers was taken up and induced effective silencing of DUX4 upon local treatment in vivo. The LNA gapmers developed here will help facilitate the development of FSHD therapies.

Ribosome collisions trigger cis-acting feedback inhibition of translation initiation

Translation of aberrant mRNAs can cause ribosomes to stall, leading to collisions with trailing ribosomes. Collided ribosomes are specifically recognised by ZNF598 to initiate protein and mRNA quality control pathways. Here we found using quantitative proteomics of collided ribosomes that EDF1 is a ZNF598-independent sensor of ribosome collisions. EDF1 stabilises GIGYF2 at collisions to inhibit translation initiation in cis via 4EHP. The GIGYF2 axis acts independently of the ZNF598 axis, but each pathway's output is more pronounced without the other. We propose that the widely conserved and highly abundant EDF1 monitors the transcriptome for excessive ribosome density, then triggers a GIGYF2-mediated response to locally and temporarily reduce ribosome loading. Only when collisions persist is translation abandoned to initiate ZNF598-dependent quality control. This tiered response to ribosome collisions would allow cells to dynamically tune translation rates while ensuring fidelity of the resulting protein products.

mRNA Editing, Processing and Quality Control in Caenorhabditis elegans

While DNA serves as the blueprint of life, the distinct functions of each cell are determined by the dynamic expression of genes from the static genome. The amount and specific sequences of RNAs expressed in a given cell involves a number of regulated processes including RNA synthesis (transcription), processing, splicing, modification, polyadenylation, stability, translation, and degradation. As errors during mRNA production can create gene products that are deleterious to the organism, quality control mechanisms exist to survey and remove errors in mRNA expression and processing. Here, we will provide an overview of mRNA processing and quality control mechanisms that occur in Caenorhabditis elegans, with a focus on those that occur on protein-coding genes after transcription initiation. In addition, we will describe the genetic and technical approaches that have allowed studies in C. elegans to reveal important mechanistic insight into these processes.

Detection and Degradation of Stalled Nascent Chains via Ribosome-Associated Quality Control

Stalled protein synthesis produces defective nascent chains that can harm cells. In response, cells degrade these nascent chains via a process called ribosome-associated quality control (RQC). Here, we review the irregularities in the translation process that cause ribosomes to stall as well as how cells use RQC to detect stalled ribosomes, ubiquitylate their tethered nascent chains, and deliver the ubiquitylated nascent chains to the proteasome. We additionally summarize how cells respond to RQC failure.

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

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

Wide-ranging transcriptome remodelling mediated by alternative polyadenylation in response to abiotic stresses in Sorghum

Alternative polyadenylation (APA) regulates diverse developmental and physiological processes through its effects on gene expression, mRNA stability, translatability, and transport. Sorghum is a major cereal crop in the world and, despite its importance, not much is known about the role of post-transcriptional regulation in mediating responses to abiotic stresses in Sorghum. A genome-wide APA analysis unveiled widespread occurrence of APA in Sorghum in response to drought, heat, and salt stress. Abiotic stress treatments incited changes in poly(A) site choice in a large number of genes. Interestingly, abiotic stresses led to the re-directing of transcriptional output into non-productive pathways defined by the class of poly(A) site utilized. This result revealed APA to be part of a larger global response of Sorghum to abiotic stresses that involves the re-direction of transcriptional output into non-productive transcriptional and translational pathways. Large numbers of stress-inducible poly(A) sites could not be linked with known, annotated genes, suggestive of the existence of numerous unidentified genes whose expression is strongly regulated by abiotic stresses. Furthermore, we uncovered a novel stress-specific cis-element in intronic poly(A) sites used in drought- and heat-stressed plants that might play an important role in non-canonical poly(A) site choice in response to abiotic stresses.

Transcriptome visualization and data availability at the Saccharomyces Genome Database

The Saccharomyces Genome Database (SGD; www.yeastgenome.org) maintains the official annotation of all genes in the Saccharomyces cerevisiae reference genome and aims to elucidate the function of these genes and their products by integrating manually curated experimental data. Technological advances have allowed researchers to profile RNA expression and identify transcripts at high resolution. These data can be configured in web-based genome browser applications for display to the general public. Accordingly, SGD has incorporated published transcript isoform data in our instance of JBrowse, a genome visualization platform. This resource will help clarify S. cerevisiae biological processes by furthering studies of transcriptional regulation, untranslated regions, genome engineering, and expression quantification in S. cerevisiae.

PolyASite 2.0: a consolidated atlas of polyadenylation sites from 3' end sequencing

Generated by 3′ end cleavage and polyadenylation at alternative polyadenylation (poly(A)) sites, alternative terminal exons account for much of the variation between human transcript isoforms. More than a dozen protocols have been developed so far for capturing and sequencing RNA 3′ ends from a variety of cell types and species. In previous studies, we have used these data to uncover novel regulatory signals and cell type-specific isoforms. Here we present an update of the PolyASite (https://polyasite.unibas.ch) resource of poly(A) sites, constructed from publicly available human, mouse and worm 3′ end sequencing datasets by enforcing uniform quality measures, including the flagging of putative internal priming sites. Through integrated processing of all data, we identified and clustered sites that are closely spaced and share polyadenylation signals, as these are likely the result of stochastic variations in processing. For each cluster, we identified the representative - most frequently processed - site and estimated the relative use in the transcriptome across all samples. We have established a modern web portal for efficient finding, exploration and export of data. Database generation is fully automated, greatly facilitating incorporation of new datasets and the updating of underlying genome resources.

Ipa1 Is an RNA Polymerase II Elongation Factor that Facilitates Termination by Maintaining Levels of the Poly(A) Site Endonuclease Ysh1

The yeast protein Ipa1 was recently discovered to interact with the Ysh1 endonuclease of the prem-RNA cleavage and polyadenylation (C/P) machinery, and Ipa1 mutation impairs 3′end processing. We report that Ipa1 globally promotes proper transcription termination and poly(A) site selection, but with variable effects on genes depending upon the specific configurations of polyadenylation signals. Our findings suggest that the role of Ipa1 in termination is mediated through interaction with Ysh1, since Ipa1 mutation leads to decrease in Ysh1 and poor recruitment of the C/P complex to a transcribed gene. The Ipa1 association with transcriptionally active chromatin resembles that of elongation factors, and the mutant shows defective Pol II elongation kinetics in vivo. Ysh1 overexpression in the Ipa1 mutant rescues the termination defect, but not the mutant’s sensitivity to 6-azauracil, an indicator of defective elongation. Our findings support a model in which an Ipa1/Ysh1 complex helps coordinate transcription elongation and 3′ end processing.

The essential, uncharacterized Ipa1 protein was recently discovered to interact with the Ysh1 endonuclease of the pre-mRNA cleavage and polyadenylation machinery. Pearson et al. propose that the Ipa1/Ysh1 interaction provides the cell with a means to coordinate and regulate transcription elongation with 3′ end processing in accordance with the cell’s needs.

Characterization of basal and estrogen-regulated antisense transcription in breast cancer cells: Role in regulating sense transcription

Estrogen-responsive breast cancer cells exhibit both basal and estrogen-regulated transcriptional programs, which lead to the transcription of many different transcription units (i.e., genes), including those that produce coding and non-coding sense (e.g., mRNA, lncRNA) and antisense (i.e., asRNA) transcripts. We have previously characterized the global basal and estrogen-regulated transcriptomes in estrogen receptor alpha (ERα)-positive MCF-7 breast cancer cells. Herein, we have mined genomic data to define three classes of antisense transcription in MCF-7 cells based on where their antisense transcription termination sites reside relative to their cognate sense mRNA and lncRNA genes. These three classes differ in their response to estrogen treatment, the enrichment of a number of genomic features associated with active promoters (H3K4me3, RNA polymerase II, open chromatin architecture), and the biological functions of their cognate sense genes as analyzed by DAVID gene ontology. We further characterized two estrogen-regulated antisense transcripts arising from the MYC gene in MCF-7 cells, showing that these antisense transcripts are 5'-capped, 3'-polyadenylated, and localized to different compartments of the cell. Together, our analyses have revealed distinct classes of antisense transcription correlated to different biological processes and response to estrogen stimulation, uncovering another layer of hormone-regulated gene regulation.

Re-annotation of the Theileria parva genome refines 53% of the proteome and uncovers essential components of N-glycosylation, a conserved pathway in many organisms

BACKGROUND

The apicomplexan parasite Theileria parva causes a livestock disease called East coast fever (ECF), with millions of animals at risk in sub-Saharan East and Southern Africa, the geographic distribution of T. parva. Over a million bovines die each year of ECF, with a tremendous economic burden to pastoralists in endemic countries. Comprehensive, accurate parasite genome annotation can facilitate the discovery of novel chemotherapeutic targets for disease treatment, as well as elucidate the biology of the parasite. However, genome annotation remains a significant challenge because of limitations in the quality and quantity of the data being used to inform the location and function of protein-coding genes and, when RNA data are used, the underlying biological complexity of the processes involved in gene expression. Here, we apply our recently published RNAseq dataset derived from the schizont life-cycle stage of T. parva to update structural and functional gene annotations across the entire nuclear genome.

RESULTS

The re-annotation effort lead to evidence-supported updates in over half of all protein-coding sequence (CDS) predictions, including exon changes, gene merges and gene splitting, an increase in average CDS length of approximately 50 base pairs, and the identification of 128 new genes. Among the new genes identified were those involved in N-glycosylation, a process previously thought not to exist in this organism and a potentially new chemotherapeutic target pathway for treating ECF. Alternatively-spliced genes were identified, and antisense and multi-gene family transcription were extensively characterized.

CONCLUSIONS

The process of re-annotation led to novel insights into the organization and expression profiles of protein-coding sequences in this parasite, and uncovered a minimal N-glycosylation pathway that changes our current understanding of the evolution of this post-translational modification in apicomplexan parasites.

RQT complex dissociates ribosomes collided on endogenous RQC substrate SDD1

Ribosome-associated quality control (RQC) represents a rescue pathway in eukaryotic cells that is triggered upon translational stalling. Collided ribosomes are recognized for subsequent dissociation followed by degradation of nascent peptides. However, endogenous RQC-inducing sequences and the mechanism underlying the ubiquitin-dependent ribosome dissociation remain poorly understood. Here, we identified SDD1 messenger RNA from Saccharomyces cerevisiae as an endogenous RQC substrate and reveal the mechanism of its mRNA-dependent and nascent peptide-dependent translational stalling. In vitro translation of SDD1 mRNA enabled the reconstitution of Hel2-dependent polyubiquitination of collided disomes and, preferentially, trisomes. The distinct trisome architecture, visualized using cryo-EM, provides the structural basis for the more-efficient recognition by Hel2 compared with that of disomes. Subsequently, the Slh1 helicase subunit of the RQC trigger (RQT) complex preferentially dissociates the first stalled polyubiquitinated ribosome in an ATP-dependent manner. Together, these findings provide fundamental mechanistic insights into RQC and its physiological role in maintaining cellular protein homeostasis.

Molecular mechanism of translational stalling by inhibitory codon combinations and poly(A) tracts

Inhibitory codon pairs and poly(A) tracts within the translated mRNA cause ribosome stalling and reduce protein output. The molecular mechanisms that drive these stalling events, however, are still unknown. Here, we use a combination of in vitro biochemistry, ribosome profiling, and cryo-EM to define molecular mechanisms that lead to these ribosome stalls. First, we use an in vitro reconstituted yeast translation system to demonstrate that inhibitory codon pairs slow elongation rates which are partially rescued by increased tRNA concentration or by an artificial tRNA not dependent on wobble base-pairing. Ribosome profiling data extend these observations by revealing that paused ribosomes with empty A sites are enriched on these sequences. Cryo-EM structures of stalled ribosomes provide a structural explanation for the observed effects by showing decoding-incompatible conformations of mRNA in the A sites of all studied stall- and collision-inducing sequences. Interestingly, in the case of poly(A) tracts, the inhibitory conformation of the mRNA in the A site involves a nucleotide stacking array. Together, these data demonstrate a novel mRNA-induced mechanisms of translational stalling in eukaryotic ribosomes.

Mechanism of ribosome stalling during translation of a poly(A) tail

Faulty or damaged mRNAs are detected by the cell when translating ribosomes stall during elongation and trigger pathways of mRNA decay, nascent protein degradation, and ribosome recycling. The most common mRNA defect in eukaryotes is probably inappropriate poly-adenylation at near-cognate sites within the coding region. How ribosomes stall selectively when they encounter poly(A) is unclear. Here, we use biochemical and structural approaches in mammalian systems to show that poly-lysine, encoded by poly(A), favors a peptidyl-tRNA conformation sub-optimal for peptide bond formation. This conformation partially slows elongation, permitting poly(A) mRNA in the ribosome’s decoding center to adopt an rRNA-stabilized single-stranded helix. The reconfigured decoding center clashes with incoming aminoacyl-tRNA, thereby precluding elongation. Thus, coincidence detection of poly-lysine in the exit tunnel and poly(A) in the decoding center allows ribosomes to detect aberrant mRNAs selectively, stall elongation, and trigger downstream quality control pathways essential for cellular homeostasis.

Regulation of mTOR signaling by long non-coding RNA

The mechanistic target of rapamycin (mTOR) is a major signaling hub that coordinates cellular and organismal responses, such as cell growth, proliferation, apoptosis, and metabolism. Dysregulation of mTOR signaling occurs in many human diseases, and there are significant ongoing efforts to pharmacologically target this pathway. Long noncoding RNAs (lncRNA), defined by a length > 200 nucleotides and absence of a long open-reading-frame, are a class of non-protein-coding RNAs. Mutations and dysregulations of lncRNAs are directly linked to the development and progression of many diseases, including cancer, diabetes, and neurologic disorders. Recent findings reveal diverse functions for lncRNA that include transcriptional regulation, organization of nuclear domains, and regulation of proteins or RNA molecules. Despite considerable development in our understanding of lncRNA over the past decade, only a fraction of annotated lncRNAs has been examined for biological function. In addition, lncRNAs have emerged as therapeutic targets due to their ability to modulate multiple pathways, including mTOR signaling. This review will provide an up-to-date summary of lncRNAs that are involved in regulating mTOR pathway.

Dynamic Transcriptional Landscape of the Early Chick Embryo

Defining the dynamic transcriptome of the early embryo at high resolution would assist greatly in understanding vertebrate development. Here, we describe the dynamic transcription landscape of early chick embryo development using advanced single-molecule long-read isoform sequencing (Iso-Seq) and RNA-Seq technology. Our transcriptomic profiling reflected the time course of chicken embryonic development from day 1 to day 8 of incubation, a period encompassing gastrulation, somitogenesis, and organogenesis. This analysis identified transcriptional isoforms, alternative splicing (AS) events, fusion transcripts, alternative polyadenylation (APA) sites, and novel genes. Our results showed that intron retention (IR) represented the most abundant AS type and displayed distinct features and dynamic modulation during development. Moreover, we constructed a high-resolution expression profile across embryonic development. Our combined expression dataset correlates distinct gene clusters with specific morphological changes, and provides the first framework for the molecular basis of early chicken embryogenesis. Analysis of gene expression in the developing chicken embryo highlighted the dynamic nature and complexity of the chicken transcriptome and demonstrated that dramatically increased IR events are associated with distinct gene sets.

Probing In Vivo Structure of Individual mRNA 3' Isoforms Using Dimethyl Sulfate

The DMS region extraction and deep sequencing (DREADS) procedure was designed to probe RNA structure in vivo and to link this structural information to specific 3' isoforms. Growing cells are treated with the alkylating agent dimethyl sulfate (DMS), which enters easily into cells and modifies RNA molecules at solvent-exposed A and C residues. RNA is isolated, and sequencing libraries are constructed in a manner that preserves the identities of individual mRNA isoforms arising from alternative cleavage/polyadenylation sites. During the cDNA synthesis step of library construction, the progress of reverse transcriptase (RT) is blocked when it encounters a DMS modification on the RNA, leading to disproportionate cDNA termination adjacent to DMS-modified positions. After paired-end deep sequencing, the downstream end of each sequenced fragment is mapped to a specific cleavage/poly(A) site representing an individual mRNA 3' isoform. The upstream mapped end of the sequenced fragment defines where the RT reaction stopped. Over the population of all sequenced fragments derived from a particular isoform, A and C positions that are overrepresented next to the upstream endpoints in the DMS sample (relative to a parallel untreated control) are inferred to have been DMS modified, and hence solvent exposed. This method thus allows in vivo structural information obtained using DMS to be linked to individual mRNA 3' isoforms. © 2019 by John Wiley & Sons, Inc.

Identification and characterization of a novel group of natural anti-sense transcripts from RNA1.2 gene locus of human cytomegalovirus

Background:

Natural anti-sense transcripts (NATs), which are transcribed from the complementary DNA strand of annotated genes, exert regulatory function of gene expression. Increasing studies recognized anti-sense transcription widespread throughout human cytomegalovirus (HCMV) genome, whereas the anti-sense transcription of RNA1.2 gene locus has never been investigated. In this study, the transcription of the RNA1.2 anti-sense strand was investigated in clinically isolated HCMV strain.

Methods:

Strand-specific high-through RNA-sequencing (RNA-seq) was performed to find possible anti-sense transcripts (ASTs). For analyzing and visualization of RNA-seq data sets, Integrative Genomics Viewer software was applied. To confirm these possibilities, Northern blotting and rapid amplification of cDNA ends (RACE) were used.

Results:

Transcription of the opposite strand of RNA1.2 gene locus was detected by RNA-sequencing using RNAs extracted from human embryonic lung fibroblasts infected with HCMV clinical isolate HAN. At least three HCMV NATs, named RNA1.2 AST 1, RNA1.2 AST2, and RNA1.2 AST3, were characterized by Northern blotting and RACE analyses. These RNA1.2 ASTs orientated from the complementary strand of RNA1.2 locus during the late phase of HCMV infection. The 5′- and 3′-termini of these transcripts were located within the opposite sequence of the predicted RNA1.2 gene.

Conclusion:

A cluster of novel NATs was transcribed from the opposite sequence of the HCMV RNA1.2 gene region.

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.

Alternative cleavage and polyadenylation in health and disease

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

The endonuclease Cue2 cleaves mRNAs at stalled ribosomes during No Go Decay

Translation of problematic sequences in mRNAs leads to ribosome collisions that trigger a series of quality control events including ribosome rescue, degradation of the stalled nascent polypeptide, and targeting of the mRNA for decay (No Go Decay or NGD). Using a reverse genetic screen in yeast, we identify Cue2 as the conserved endonuclease that is recruited to stalled ribosomes to promote NGD. Ribosome profiling and biochemistry provide strong evidence that Cue2 cleaves mRNA within the A site of the colliding ribosome. We demonstrate that NGD primarily proceeds via Xrn1-mediated exonucleolytic decay and Cue2-mediated endonucleolytic decay normally constitutes a secondary decay pathway. Finally, we show that the Cue2-dependent pathway becomes a major contributor to NGD in cells depleted of factors required for the resolution of stalled ribosome complexes. Together these results provide insights into how multiple decay processes converge to process problematic mRNAs in eukaryotic cells.​.

OTUD6B-AS1 Might Be a Novel Regulator of Apoptosis in Systemic Sclerosis

Antisense long non-coding RNAs (AS lncRNAs) have increasingly been recognized as important regulators of gene expression and they have been found to play key roles in several diseases. However, very little is known about the role of AS lncRNAs in fibrotic diseases such as systemic sclerosis (SSc). Our recent screening experiments by RNA sequencing showed that ovarian tumor domain containing 6B antisense RNA1 (OTUD6B-AS1) and its sense gene OTUD6B were significantly downregulated in SSc skin biopsies. Therefore, we aimed to identify key regulators of OTUD6B-AS1 and to analyze the functional relevance of OTUD6B-AS1 in SSc. OTUD6B-AS1 and OTUD6B expression in SSc and healthy control (HC) dermal fibroblasts (Fb) after stimulation with transforming growth factor-β (TGFβ), Interleukin (IL)-4, IL-13, and platelet-derived growth factor (PDGF) was analyzed by qPCR. To identify the functional role of OTUD6B-AS1, dermal Fb or human pulmonary artery smooth muscle cells (HPASMC) were transfected with a locked nucleic acid antisense oligonucleotide (ASO) targeting OTUD6B-AS1. Proliferation was measured by BrdU and real-time proliferation assay. Apoptosis was measured by Caspase 3/7 assay and Western blot for cleaved caspase 3. While no difference was recorded at the basal level between HC and SSc dermal Fb, the expression of OTUD6B-AS1 and OTUD6B was significantly downregulated in both SSc and HC dermal Fb after PDGF stimulation in a time-dependent manner. Only mild and inconsistent effects were observed with TGFβ, IL-4, and IL-13. OTUD6B-AS1 knockdown in Fb and HPASMC did not affect extracellular matrix or pro-fibrotic/proinflammatory cytokine production. However, OTUD6B-AS1 knockdown significantly increased Cyclin D1 expression at the mRNA and protein level. Moreover, silencing of OTUD6B-AS1 significantly reduced proliferation and suppressed apoptosis in both dermal Fb and HPASMC. OTUD6B-AS1 knockdown did not affect OTUD6B expression at the mRNA level and protein level. Our data suggest that OTUD6B-AS1 regulates proliferation and apoptosis via cyclin D1 expression in a sense gene independent manner. This is the first report investigating the function of OTUD6B-AS1. Our data shed light on a novel apoptosis resistance mechanism in Fb and vascular smooth muscle cells that might be relevant for pathogenesis of SSc.