Formation of mRNA 3' ends in eukaryotes: mechanism, regulation, and interrelationships with other steps in mRNA synthesis

A Developmental Mechanism to Regulate Alternative Polyadenylation in an Adult Stem Cell Lineage

Alternative Cleavage and Polyadenylation (APA) often results in production of mRNA isoforms with either longer or shorter 3'UTRs from the same genetic locus, potentially impacting mRNA translation, localization and stability. Developmentally regulated APA can thus make major contributions to cell-type-specific gene expression programs as cells differentiate. During Drosophila spermatogenesis, approximately 500 genes undergo APA when proliferating spermatogonia differentiate into spermatocytes, producing transcripts with shortened 3' UTRs, leading to profound stage-specific changes in the proteins expressed. The molecular mechanisms that specify usage of upstream polyadenylation sites in spermatocytes are thus key to understanding the changes in cell state. Here, we show that upregulation of PCF11 and Cbc, the two components of Cleavage Factor II (CFII), orchestrates APA during Drosophila spermatogenesis. Knock down of PCF11 or cbc in spermatocytes caused dysregulation of APA, with many transcripts normally cleaved at a proximal site in spermatocytes now cleaved at their distal site, as in spermatogonia. Forced overexpression of CFII components in spermatogonia switched cleavage of some transcripts to the proximal site normally used in spermatocytes. Our findings reveal a developmental mechanism where changes in expression of specific cleavage factors can direct cell-type-specific APA at selected genes.

Deciphering the cis-regulatory landscape of natural yeast Transcript Leaders

Protein synthesis is a vital process that is highly regulated at the initiation step of translation. Eukaryotic 5' transcript leaders (TLs) contain a variety of cis-regulatory features that influence translation and mRNA stability. However, the relative influences of these features in natural TLs are poorly characterized. To address this, we used massively parallel reporter assays (MPRAs) to quantify RNA levels, ribosome loading, and protein levels from 11,027 natural yeast TLs in vivo and systematically compared the relative impacts of their sequence features on gene expression. We found that yeast TLs influence gene expression over two orders of magnitude. While a leaky scanning model using Kozak contexts and uAUGs explained half of the variance in expression across transcript leaders, the addition of other features explained ~70% of gene expression variation. Our analyses detected key cis-acting sequence features, quantified their effects in vivo, and compared their roles to motifs reported from an in vitro study of ribosome recruitment. In addition, our work quantitated the effects of alternative transcription start site usage on gene expression in yeast. Thus, our study provides new quantitative insights into the roles of TL cis-acting sequences in regulating gene expression.

Chimera RNA transcribed from integrated HPV18 genome with adjacent host genomic region promotes oncogenic gene expression through condensate formation

Most cervical cancers are caused by human papillomavirus (HPV) infection. In HeLa cells, the HPV18 viral genome is integrated at chromosome 8q24.21 and activates transcription of the proto-oncogene c-Myc. However, the mechanism of how the integrated HPV genome and its transcribed RNAs exhibit transcription activation function has not been fully elucidated. In this study, we found that HPV18 transcripts contain an enhancer RNA-like function to activate proximal genes including CCAT1-5L and c-Myc. We showed that the human genome-integrated HPV18 genes are activated by transcription coregulators including BRD4 and Mediator. The transcribed HPV18 RNAs form a liquid-like condensate at chromosome 8q24.21 locus, which in turn accumulates RNA polymerase II. Moreover, we focused on a relatively uncharacterized transcript from the upstream region of CCAT1, named URC. The URC RNA is transcribed as a chimera RNA with HPV18 and is composed of the 3'-untranslated region of the HPV18 transcript. We experimentally showed that the URC contributes to stabilization of HPV18 RNAs by supplying a polyadenylation site for the HPV18 transcript. Our findings suggest that integrated HPV18 at 8q24.21 locus produces HPV18-URC chimera RNA and promotes tumorigenesis through RNA-based condensate formation.

Helicases DDX5 and DDX17 promote heterogeneity in HBV transcription termination in infected human hepatocytes

BACKGROUND & AIMS

Transcription termination fine-tunes gene expression and contributes to the specification of RNA function in eukaryotic cells. Transcription termination of HBV is subject to the recognition of the canonical polyadenylation signal (cPAS) common to all viral transcripts. However, the regulation of this cPAS and its impact on viral gene expression and replication is currently unknown.

METHODS

To unravel the regulation of HBV transcript termination, we implemented a 3' RACE (rapid amplification of cDNA ends)-PCR assay coupled to single molecule sequencing both in in vitro-infected hepatocytes and in chronically infected patients.

RESULTS

The detection of a previously unidentified transcriptional readthrough indicated that the cPAS was not systematically recognized during HBV replication in vitro and in vivo. Gene expression downregulation experiments demonstrated a role for the RNA helicases DDX5 and DDX17 in promoting viral transcriptional readthrough, which was, in turn, associated with HBV RNA destabilization and decreased HBx protein expression. RNA and chromatin immunoprecipitation, together with mutation of the cPAS sequence, suggested a direct role of DDX5 and DDX17 in functionally linking cPAS recognition to transcriptional readthrough, HBV RNA stability and replication.

CONCLUSIONS

Our findings identify DDX5 and DDX17 as crucial determinants of HBV transcriptional fidelity and as host restriction factors for HBV replication.

IMPACT AND IMPLICATIONS

HBV covalently closed circular (ccc)DNA degradation or functional inactivation remains the holy grail for the achievement of HBV cure. Transcriptional fidelity is a cornerstone in the regulation of gene expression. Here, we demonstrate that two helicases, DDX5 and DDX17, inhibit recognition of the HBV polyadenylation signal and thereby transcriptional termination, thus decreasing HBV RNA stability and acting as restriction factors for efficient cccDNA transcription and viral replication. The observation that DDX5 and DDX17 are downregulated in patients chronically infected with HBV suggests a role for these helicases in HBV persistence in vivo. These results open new perspectives for researchers aiming at identifying new targets to neutralise cccDNA transcription.

Identifying human pre-mRNA cleavage and polyadenylation factors by genome-wide CRISPR screens using a dual fluorescence readthrough reporter

Messenger RNA precursors (pre-mRNA) generally undergo 3' end processing by cleavage and polyadenylation (CPA), which is specified by a polyadenylation site (PAS) and adjacent RNA sequences and regulated by a large variety of core and auxiliary CPA factors. To date, most of the human CPA factors have been discovered through biochemical and proteomic studies. However, genetic identification of the human CPA factors has been hampered by the lack of a reliable genome-wide screening method. We describe here a dual fluorescence readthrough reporter system with a PAS inserted between two fluorescent reporters. This system enables measurement of the efficiency of 3' end processing in living cells. Using this system in combination with a human genome-wide CRISPR/Cas9 library, we conducted a screen for CPA factors. The screens identified most components of the known core CPA complexes and other known CPA factors. The screens also identified CCNK/CDK12 as a potential core CPA factor, and RPRD1B as a CPA factor that binds RNA and regulates the release of RNA polymerase II at the 3' ends of genes. Thus, this dual fluorescence reporter coupled with CRISPR/Cas9 screens reliably identifies bona fide CPA factors and provides a platform for investigating the requirements for CPA in various contexts.

How is polyadenylation restricted to 3'-untranslated regions?

Polyadenylation occurs at numerous sites within 3'-untranslated regions (3'-UTRs) but rarely within coding regions. How does Pol II travel through long coding regions without generating poly(A) sites, yet then permits promiscuous polyadenylation once it reaches the 3'-UTR? The cleavage/polyadenylation (CpA) machinery preferentially associates with 3'-UTRs, but it is unknown how its recruitment is restricted to 3'-UTRs during Pol II elongation. Unlike coding regions, 3'-UTRs have long AT-rich stretches of DNA that may be important for restricting polyadenylation to 3'-UTRs. Recognition of the 3'-UTR could occur at the DNA (AT-rich), RNA (AU-rich), or RNA:DNA hybrid (rU:dA- and/or rA:dT-rich) level. Based on the nucleic acid critical for 3'-UTR recognition, there are three classes of models, not mutually exclusive, for how the CpA machinery is selectively recruited to 3'-UTRs, thereby restricting where polyadenylation occurs: (1) RNA-based models suggest that the CpA complex directly (or indirectly through one or more intermediary proteins) binds long AU-rich stretches that are exposed after Pol II passes through these regions. (2) DNA-based models suggest that the AT-rich sequence affects nucleosome depletion or the elongating Pol II machinery, resulting in dissociation of some elongation factors and subsequent recruitment of the CpA machinery. (3) RNA:DNA hybrid models suggest that preferential destabilization of the Pol II elongation complex at rU:dA- and/or rA:dT-rich duplexes bridging the nucleotide addition and RNA exit sites permits preferential association of the CpA machinery with 3'-UTRs. Experiments to provide evidence for one or more of these models are suggested.

Transcription termination and readthrough in African swine fever virus

Introduction

African swine fever virus (ASFV) is a nucleocytoplasmic large DNA virus (NCLDV) that encodes its own host-like RNA polymerase (RNAP) and factors required to produce mature mRNA. The formation of accurate mRNA 3' ends by ASFV RNAP depends on transcription termination, likely enabled by a combination of sequence motifs and transcription factors, although these are poorly understood. The termination of any RNAP is rarely 100% efficient, and the transcriptional "readthrough" at terminators can generate long mRNAs which may interfere with the expression of downstream genes. ASFV transcriptome analyses reveal a landscape of heterogeneous mRNA 3' termini, likely a combination of bona fide termination sites and the result of mRNA degradation and processing. While short-read sequencing (SRS) like 3' RNA-seq indicates an accumulation of mRNA 3' ends at specific sites, it cannot inform about which promoters and transcription start sites (TSSs) directed their synthesis, i.e., information about the complete and unprocessed mRNAs at nucleotide resolution.

Methods

Here, we report a rigorous analysis of full-length ASFV transcripts using long-read sequencing (LRS). We systematically compared transcription termination sites predicted from SRS 3' RNA-seq with 3' ends mapped by LRS during early and late infection.

Results

Using in-vitro transcription assays, we show that recombinant ASFV RNAP terminates transcription at polyT stretches in the non-template strand, similar to the archaeal RNAP or eukaryotic RNAPIII, unaided by secondary RNA structures or predicted viral termination factors. Our results cement this T-rich motif (U-rich in the RNA) as a universal transcription termination signal in ASFV. Many genes share the usage of the same terminators, while genes can also use a range of terminators to generate transcript isoforms varying enormously in length. A key factor in the latter phenomenon is the highly abundant terminator readthrough we observed, which is more prevalent during late compared with early infection.

Discussion

This indicates that ASFV mRNAs under the control of late gene promoters utilize different termination mechanisms and factors to early promoters and/or that cellular factors influence the viral transcriptome landscape differently during the late stages of infection.

Application of the iPLUS non-coding sequence in improving biopharmaceuticals production

The biotechnological landscape has witnessed significant growth in biological therapeutics particularly in the field of recombinant protein production. Here we investigate the function of 3'UTR cis-regulatory elements in increasing mRNA and protein levels in different biological therapeutics and model systems, spanning from monoclonal antibodies to mRNA vaccines. We explore the regulatory function of iPLUS - a universal sequence capable of consistently augmenting recombinant protein levels. By incorporating iPLUS in a vector to express a monoclonal antibody used in immunotherapy, in a mammalian cell line used by the industry (ExpiCHO), trastuzumab production increases by 2-fold. As yeast Pichia pastoris is widely used in the manufacture of industrial enzymes and pharmaceuticals, we then used iPLUS in tandem (3x) and iPLUSv2 (a variant of iPLUS) to provide proof-of-concept data that it increases the production of a reporter protein more than 100-fold. As iPLUS functions by also increasing mRNA levels, we hypothesize that these sequences could be used as an asset in the mRNA vaccine industry. In fact, by including iPLUSv2 downstream of Spike we were able to double its production. Moreover, the same effect was observed when we introduced iPLUSv2 downstream of MAGEC2, a tumor-specific antigen tested for cancer mRNA vaccines. Taken together, our study provides data (TLR4) showing that iPLUS may be used as a valuable asset in a variety of systems used by the biotech and biopharmaceutical industry. Our results underscore the critical role of non-coding sequences in controlling gene expression, offering a promising avenue to accelerate, enhance, and cost-effectively optimize biopharmaceutical production processes.

Emerging Drug Targets for Antimalarial Drug Discovery: Validation and Insights into Molecular Mechanisms of Function

Approximately 619,000 malaria deaths were reported in 2021, and resistance to recommended drugs, including artemisinin-combination therapies (ACTs), threatens malaria control. Treatment failure with ACTs has been found to be as high as 93% in northeastern Thailand, and parasite mutations responsible for artemisinin resistance have already been reported in some African countries. Therefore, there is an urgent need to identify alternative treatments with novel targets. In this Perspective, we discuss some promising antimalarial drug targets, including enzymes involved in proteolysis, DNA and RNA metabolism, protein synthesis, and isoprenoid metabolism. Other targets discussed are transporters, Plasmodium falciparum acetyl-coenzyme A synthetase, N-myristoyltransferase, and the cyclic guanosine monophosphate-dependent protein kinase G. We have outlined mechanistic details, where these are understood, underpinning the biological roles and hence druggability of such targets. We believe that having a clear understanding of the underlying chemical interactions is valuable to medicinal chemists in their quest to design appropriate inhibitors.

Natural variation in the plant polyadenylation complex

Messenger RNA polyadenylation, the process wherein the primary RNA polymerase II transcript is cleaved and a poly(A) tract added, is a key step in the expression of genes in plants. Moreover, it is a point at which gene expression may be regulated by determining the functionality of the mature mRNA. Polyadenylation is mediated by a complex (the polyadenylation complex, or PAC) that consists of between 15 and 20 subunits. While the general functioning of these subunits may be inferred by extending paradigms established in well-developed eukaryotic models, much remains to be learned about the roles of individual subunits in the regulation of polyadenylation in plants. To gain further insight into this, we conducted a survey of variability in the plant PAC. For this, we drew upon a database of naturally-occurring variation in numerous geographic isolates of Arabidopsis thaliana. For a subset of genes encoding PAC subunits, the patterns of variability included the occurrence of premature stop codons in some Arabidopsis accessions. These and other observations lead us to conclude that some genes purported to encode PAC subunits in Arabidopsis are actually pseudogenes, and that others may encode proteins with dispensable functions in the plant. Many subunits of the PAC showed patterns of variability that were consistent with their roles as essential proteins in the cell. Several other PAC subunits exhibit patterns of variability consistent with selection for new or altered function. We propose that these latter subunits participate in regulatory interactions important for differential usage of poly(A) sites.

RNA Folding, Mutation, and Detection

The structure of RNA molecules is absolutely critical to their functions in a biological system. RNA structure is dynamic and changes in response to cellular needs. Within the last few decades, there has been an increased interest in studying the structure of RNA molecules and how they change to support the needs of the cell in different conditions. Selective 2'-hydroxyl acylation-based mutational profiling using high-throughput sequencing is a powerful method to predict the secondary structure of RNA molecules both in vivo and in immunopurified samples. Selective 2'-hydroxyl acylation-based mutational profiling using high-throughput sequencing works by adding bulky groups onto accessible "flexible" bases in an RNA molecule that are not involved in any base-pairing or RNA-protein interactions. When the RNA is reverse transcribed into cDNA, the bulky groups are incorporated as base mutations, which can be compared to an unmodified control to identify the locations of flexible bases. The comparison of sequence data between modified and unmodified samples allows the computer software program (developed to generate reactivity profiles) to generate RNA secondary structure models. These models can be compared in a variety of conditions to determine how specific stimuli influence RNA secondary structures.

A direct interaction between CPF and RNA Pol II links RNA 3' end processing to transcription

Transcription termination by RNA polymerase II (RNA Pol II) is linked to RNA 3' end processing by the cleavage and polyadenylation factor (CPF or CPSF). CPF contains endonuclease, poly(A) polymerase, and protein phosphatase activities, which cleave and polyadenylate pre-mRNAs and dephosphorylate RNA Pol II to control transcription. Exactly how the RNA 3' end processing machinery is coupled to transcription remains unclear. Here, we combine in vitro reconstitution, structural studies, and genome-wide analyses to show that yeast CPF physically and functionally interacts with RNA Pol II. Surprisingly, CPF-mediated dephosphorylation promotes the formation of an RNA Pol II stalk-to-stalk homodimer in vitro. This dimer is compatible with transcription but not with the binding of transcription elongation factors. Disruption of the dimerization interface in cells causes transcription defects, including altered RNA Pol II abundance on protein-coding genes, tRNA genes, and intergenic regions. We hypothesize that RNA Pol II dimerization may provide a mechanistic basis for the allosteric model of transcription termination.

The polyA tail facilitates splicing of last introns with weak 3' splice sites via PABPN1

The polyA tail of mRNAs is important for many aspects of RNA metabolism. However, whether and how it regulates pre-mRNA splicing is still unknown. Here, we report that the polyA tail acts as a splicing enhancer for the last intron via the nuclear polyA binding protein PABPN1 in HeLa cells. PABPN1-depletion induces the retention of a group of introns with a weaker 3' splice site, and they show a strong 3'-end bias and mainly locate in nuclear speckles. The polyA tail is essential for PABPN1-enhanced last intron splicing and functions in a length-dependent manner. Tethering PABPN1 to nonpolyadenylated transcripts also promotes splicing, suggesting a direct role for PABPN1 in splicing regulation. Using TurboID-MS, we construct the PABPN1 interactome, including many spliceosomal and RNA-binding proteins. Specifically, PABPN1 can recruit RBM26&27 to promote splicing by interacting with the coiled-coil and RRM domain of RBM27. PABPN1-regulated terminal intron splicing is conserved in mice. Together, our study establishes a novel mode of post-transcriptional splicing regulation via the polyA tail and PABPN1.

Advancements of in vitro transcribed mRNA (IVT mRNA) to enable translation into the clinics

The accelerated progress and approval of two mRNA-based vaccines to address the SARS-CoV-2 virus were unprecedented. This record-setting feat was made possible through the solid foundation of research on in vitro transcribed mRNA (IVT mRNA) which could be utilized as a therapeutic modality. Through decades of thorough research to overcome barriers to implementation, mRNA-based vaccines or therapeutics offer many advantages to rapidly address a broad range of applications including infectious diseases, cancers, and gene editing. Here, we describe the advances that have supported the adoption of IVT mRNA in the clinics, including optimization of the IVT mRNA structural components, synthesis, and lastly concluding with different classes of IVT RNA. Continuing interest in driving IVT mRNA technology will enable a safer and more efficacious therapeutic modality to address emerging and existing diseases.

Plant Promoters and Terminators for High-Precision Bioengineering

High-precision bioengineering and synthetic biology require fine-tuning gene expression at both transcriptional and posttranscriptional levels. Gene transcription is tightly regulated by promoters and terminators. Promoters determine the timing, tissues and cells, and levels of the expression of genes. Terminators mediate transcription termination of genes and affect mRNA levels posttranscriptionally, e.g., the 3'-end processing, stability, translation efficiency, and nuclear to cytoplasmic export of mRNAs. The promoter and terminator combination affects gene expression. In the present article, we review the function and features of plant core promoters, proximal and distal promoters, and terminators, and their effects on and benchmarking strategies for regulating gene expression.

Alternative polyadenylation regulation in cardiac development and cardiovascular disease

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

Differential alternative polyadenylation of homoeologous genes of allohexaploid wheat ABD subgenomes during drought stress response

Because allohexaploid wheat genome contains ABD subgenomes, how the expression of homoeologous genes is coordinated remains largely unknown, particularly at the co-transcriptional level. Alternative polyadenylation (APA) is an important part of co-transcriptional regulation, which is crucial in developmental processes and stress responses. Drought stress is a major threat to the stable yield of wheat. Focusing on APA, we used poly(A) tag sequencing to track poly(A) site dynamics in wheat under drought stress. The results showed that drought stress led to extensive APA involving 37-47% of differentially expressed genes in wheat. Significant poly(A) site switching was found in stress-responsive genes. Interestingly, homoeologous genes exhibit unequal numbers of poly(A) sites, divergent APA patterns with tissue specificity and time-course dynamics, and distinct 3'-UTR length changes. Moreover, differentially expressed transcripts in leaves and roots used different poly(A) signals, the up- and downregulated isoforms had distinct preferences for non-canonical poly(A) sites. Genes that encode key polyadenylation factors showed differential expression patterns under drought stress. In summary, poly(A) signals and the changes in core poly(A) factors may widely affect the selection of poly(A) sites and gene expression levels during the response to drought stress, and divergent APA patterns among homoeologous genes add extensive plasticity to this responsive network. These results not only reveal the significant role of APA in drought stress response, but also provide a fresh perspective on how homoeologous genes contribute to adaptability through transcriptome diversity. In addition, this work provides information about the ends of transcripts for a better annotation of the wheat genome.

The emerging roles of CFIm25 (NUDT21/CPSF5) in human biology and disease

The mammalian cleavage factor I subunit CFIm25 (NUDT21) binds to the UGUA sequences of precursor RNAs. Traditionally, CFIm25 is known to facilitate 3' end formation of pre-mRNAs resulting in the formation of polyadenylated transcripts. Recent studies suggest that CFIm25 may be involved in the cyclization and hence generation of circular RNAs (circRNAs) that contain UGUA motifs. These circRNAs act as competing endogenous RNAs (ceRNAs) that disrupt the ceRNA-miRNA-mRNA axis. Other emerging roles of CFIm25 include regulating both alternative splicing and alternative polyadenylation (APA). APA generates different sized transcripts that may code for different proteins, or more commonly transcripts that code for the same protein but differ in the length and sequence content of their 3' UTRs (3' UTR-APA). CFIm25 mediated global changes in 3' UTR-APA affect human physiology including spermatogenesis and the determination of cell fate. Deregulation of CFIm25 and changes in 3' UTR-APA have been implicated in several human diseases including cancer. In many cancers, CFIm25 acts as a tumor suppressor. However, there are some cancers where CFIm25 has the opposite effect. Alterations in CFIm25-driven 3' UTR-APA may also play a role in neural dysfunction and fibrosis. CFIm25 mediated 3' UTR-APA changes can be used to generate specific signatures that can be used as potential biomarkers in development and disease. Due to the emerging role of CFIm25 as a regulator of the aforementioned RNA processing events, modulation of CFIm25 levels may be a novel viable therapeutic approach. This article is categorized under: RNA Processing > 3' End Processing RNA in Disease and Development > RNA in Disease.

An examination of the metal ion content in the active sites of human endonucleases CPSF73 and INTS11

Human cleavage and polyadenylation specificity factor (CPSF)73 (also known as CPSF3) is the endoribonuclease that catalyzes the cleavage reaction for the 3'-end processing of pre-mRNAs. The active site of CPSF73 is located at the interface between a metallo-β-lactamase domain and a β-CASP domain. Two metal ions are coordinated by conserved residues, five His and two Asp, in the active site, and they are critical for the nuclease reaction. The metal ions have long been thought to be zinc ions, but their exact identity has not been examined. Here we present evidence from inductively coupled plasma mass spectrometry and X-ray diffraction analyses that a mixture of metal ions, including Fe, Zn, and Mn, is present in the active site of CPSF73. The abundance of the various metal ions is different in samples prepared from different expression hosts. Zinc is present at less than 20% abundance in a sample expressed in insect cells, but the sample is active in cleaving a pre-mRNA substrate in a reconstituted canonical 3'-end processing machinery. Zinc is present at 75% abundance in a sample expressed in human cells, which has comparable endonuclease activity. We also observe a mixture of metal ions in the active site of the CPSF73 homolog INTS11, the endonuclease for Integrator. Taken together, our results provide further insights into the role of metal ions in the activity of CPSF73 and INTS11 for RNA 3'-end processing.

Mutation of the polyadenylation complex subunit CstF77 reveals that mRNA 3' end formation and HSP101 levels are critical for a robust heat stress response

Heat shock protein 101 (HSP101) in plants, and bacterial and yeast orthologs, is essential for thermotolerance. To investigate thermotolerance mechanisms involving HSP101, we performed a suppressor screen in Arabidopsis thaliana of a missense HSP101 allele (hot1-4). hot1-4 plants are sensitive to acclimation heat treatments that are otherwise permissive for HSP101 null mutants, indicating that the hot1-4 protein is toxic. We report one suppressor (shot2, suppressor of hot1-4 2) has a missense mutation of a conserved residue in CLEAVAGE STIMULATION FACTOR77 (CstF77), a subunit of the polyadenylation complex critical for mRNA 3' end maturation. We performed ribosomal RNA depletion RNA-Seq and captured transcriptional readthrough with a custom bioinformatics pipeline. Acclimation heat treatment caused transcriptional readthrough in hot1-4 shot2, with more readthrough in heat-induced genes, reducing the levels of toxic hot1-4 protein and suppressing hot1-4 heat sensitivity. Although shot2 mutants develop like the wild type in the absence of stress and survive mild heat stress, reduction of heat-induced genes and decreased HSP accumulation makes shot2 in HSP101 null and wild-type backgrounds sensitive to severe heat stress. Our study reveals the critical function of CstF77 for 3' end formation of mRNA and the dominant role of HSP101 in dictating the outcome of severe heat stress.

Knowing when to stop: Transcription termination on protein-coding genes by eukaryotic RNAPII

Gene expression is controlled in a dynamic and regulated manner to allow for the consistent and steady expression of some proteins as well as the rapidly changing production of other proteins. Transcription initiation has been a major focus of study because it is highly regulated. However, termination of transcription also plays an important role in controlling gene expression. Transcription termination on protein-coding genes is intimately linked with 3' end cleavage and polyadenylation of transcripts, and it generally results in the production of a mature mRNA that is exported from the nucleus. Termination on many non-coding genes can also result in the production of a mature transcript. Termination is dynamically regulated-premature termination and transcription readthrough occur in response to a number of cellular signals, and these can have varied consequences on gene expression. Here, we review eukaryotic transcription termination by RNA polymerase II (RNAPII), focusing on protein-coding genes.

Relationship between NUDT21 mediated alternative polyadenylation process and tumor

Alternative polyadenylation (APA) is a molecular process that generates diversity at the 3' end of RNA polymerase II transcripts from over 60% of human genes. APA and microRNA regulation are both mechanisms of post-transcriptional regulation of gene expression. As a key molecular mechanism, Alternative polyadenylation often results in mRNA isoforms with the same coding sequence but different lengths of 3' UTRs, while microRNAs regulate gene expression by binding to specific mRNA 3' UTRs. Nudix Hydrolase 21 (NUDT21) is a crucial mediator involved in alternative polyadenylation (APA). Different studies have reported a dual role of NUDT21 in cancer (both oncogenic and tumor suppressor). The present review focuses on the functions of APA, miRNA and their interaction and roles in development of different types of tumors.NUDT21 mediated 3' UTR-APA changes can be used to generate specific signatures that can be used as potential biomarkers in development and disease. Due to the emerging role of NUDT21 as a regulator of the aforementioned RNA processing events, modulation of NUDT21 levels may be a novel viable therapeutic approach.

Identifying transcript-level differential expression in primary human immune cells

BACKGROUND

Multipotential hematopoietic stem cells differentiate into a wide variety of immune cells with a diversity of functions, including the ability to respond to a variety of stimuli. Importantly, numerous studies have demonstrated the importance of gene transcription in defining cell identity and functions. While these studies have primarily been performed at the level of the gene, it is known that key immune genes such as CD44 and CD45 generate multiple different transcripts that are differentially expressed across different immune cells, and that encode proteins with different sequences and functions. Prior genomic surveys have shown that the mechanisms for generating diversity in expressed transcripts (alternate splicing, alternate transcription start sites, etc.) are very active in immune cells, but have been lacking in terms of identifying genes with multiple transcripts, that are differentially expressed, and likely to affect cell functions.

METHODS

We first identified the set of genes that had at least two transcripts expressed in our RNA sequencing dataset generated from purified populations of neutrophils, monocytes and five lymphocyte populations (B, NK, γδ T, CD4 + T and CD8 + T) from twelve healthy donors. Next, we developed a heuristic approach to identify genes where two or more transcripts have distinct expression patterns across lymphoid and/or myeloid populations. We then focused our annotation and interpretation on differentially expressed transcripts that affect the coding sequence. This process was repeated to identify transcripts that were differentially expressed between monocytes and populations of macrophages and LPS-stimulated macrophages derived from these monocytes in vitro.

RESULTS

We found that over 55 % of genes had two or more expressed transcripts, with an average ∼3 transcripts per gene, and that 70 % of these had at least two of the transcripts that encoded proteins with different sequences. As expected, we identified a complex pattern of differential expression for multiple transcripts encoding the CD45 transmembrane protein, but we also found similar evidence for ten other genes (CD300A, FYB1, GPI, LITAF, PSMA1, PTMA, RPL32, SEPTIN9, SH3BP2, SH3KBP1) when comparing the expression patterns of transcripts within myeloid and lymphoid cells. We also identified five genes with differentially expressed transcripts associated with the transition from monocytes to macrophages (FNBP1, KLF6, and SEPTIN9) or between macrophages and LPS-stimulated macrophages (CD44, OAZ2, and SEPTIN9). For the most part, we found that the different transcripts of these genes are expected to impact specific biological functions, for example the different transcripts of SEPTIN9 likely regulate the cytoskeleton in immune cells via their interactions with actins filaments and microtubules.

CONCLUSIONS

This analytic approach successfully identified multi-transcript genes that are differentially expressed across immune cells and could be applied to other transcriptomic data.

DATA AVAILABILITY STATEMENT

Researchers can request access to the individual-level data from the current study by contacting the Montreal Heart Institute ethics committee at the following institutional email address: cer.icm@icm-mhi.org.

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

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

Allosteric links between the hydrophilic N-terminus and transmembrane core of human Na+ /H+ antiporter NHA2

The human Na+ /H+ antiporter NHA2 (SLC9B2) transports Na+ or Li+ across the plasma membrane in exchange for protons, and is implicated in various pathologies. It is a 537 amino acids protein with an 82 residues long hydrophilic cytoplasmic N-terminus followed by a transmembrane part comprising 14 transmembrane helices. We optimized the functional expression of HsNHA2 in the plasma membrane of a salt-sensitive Saccharomyces cerevisiae strain and characterized in vivo a set of mutated or truncated versions of HsNHA2 in terms of their substrate specificity, transport activity, localization, and protein stability. We identified a highly conserved proline 246, located in the core of the protein, as being crucial for ion selectivity. The replacement of P246 with serine or threonine resulted in antiporters with altered substrate specificity that were not only highly active at acidic pH 4.0 (like the native antiporter), but also at neutral pH. P246T/S versions also exhibited increased resistance to the HsNHA2-specific inhibitor phloretin. We experimentally proved that a putative salt bridge between E215 and R432 is important for antiporter function, but also structural integrity. Truncations of the first 50-70 residues of the N-terminus doubled the transport activity of HsNHA2, while changes in the charge at positions E47, E56, K57, or K58 decreased the antiporter's transport activity. Thus, the hydrophilic N-terminal part of the protein appears to allosterically auto-inhibit cation transport of HsNHA2. Our data also show this in vivo approach to be useful for a rapid screening of SNP's effect on HsNHA2 activity.

Fip1 is a multivalent interaction scaffold for processing factors in human mRNA 3' end biogenesis

3′ end formation of most eukaryotic mRNAs is dependent on the assembly of a ~1.5 MDa multiprotein complex, that catalyzes the coupled reaction of pre-mRNA cleavage and polyadenylation. In mammals, the cleavage and polyadenylation specificity factor (CPSF) constitutes the core of the 3′ end processing machinery onto which the remaining factors, including cleavage stimulation factor (CstF) and poly(A) polymerase (PAP), assemble. These interactions are mediated by Fip1, a CPSF subunit characterized by high degree of intrinsic disorder. Here, we report two crystal structures revealing the interactions of human Fip1 (hFip1) with CPSF30 and CstF77. We demonstrate that CPSF contains two copies of hFip1, each binding to the zinc finger (ZF) domains 4 and 5 of CPSF30. Using polyadenylation assays we show that the two hFip1 copies are functionally redundant in recruiting one copy of PAP, thereby increasing the processivity of RNA polyadenylation. We further show that the interaction between hFip1 and CstF77 is mediated via a short motif in the N-terminal ‘acidic’ region of hFip1. In turn, CstF77 competitively inhibits CPSF-dependent PAP recruitment and 3′ polyadenylation. Taken together, these results provide a structural basis for the multivalent scaffolding and regulatory functions of hFip1 in 3′ end processing.

Controlling gene expression with deep generative design of regulatory DNA

Design of de novo synthetic regulatory DNA is a promising avenue to control gene expression in biotechnology and medicine. Using mutagenesis typically requires screening sizable random DNA libraries, which limits the designs to span merely a short section of the promoter and restricts their control of gene expression. Here, we prototype a deep learning strategy based on generative adversarial networks (GAN) by learning directly from genomic and transcriptomic data. Our ExpressionGAN can traverse the entire regulatory sequence-expression landscape in a gene-specific manner, generating regulatory DNA with prespecified target mRNA levels spanning the whole gene regulatory structure including coding and adjacent non-coding regions. Despite high sequence divergence from natural DNA, in vivo measurements show that 57% of the highly-expressed synthetic sequences surpass the expression levels of highly-expressed natural controls. This demonstrates the applicability and relevance of deep generative design to expand our knowledge and control of gene expression regulation in any desired organism, condition or tissue.

Design of de novo synthetic regulatory DNA is a promising avenue to control gene expression in biotechnology and medicine. Here the authors present EspressionGAN, a generative adversarial network that uses genomic and transcriptomic data to generate regulatory sequences.

β-Thalassemia Intermedia Caused by the β-Globin Gene 3' Untranslated Region: Another Case Report

The 3' untranslated region (3'UTR) is associated with mRNA stability because of its involvement in 3' end processing, polyadenylation, and mRNA capping. Mutations located in this area can cause a phenotype compatible with β⁺-thalassemia (β⁺-thal). We report a Chinese subject with β-thal intermedia (β-TI) who developed transfusion-dependent anemia. Molecular studies revealed that the patient was a compound heterozygote for two β-thal alleles: codons 41/42 (-TTCT) (HBB: c.126_129delCTTT) and term codon +32 (A>C) (HBB: c.*32A>C).

NBBt-test: a versatile method for differential analysis of multiple types of RNA-seq data

Rapid development of transcriptome sequencing technologies has resulted in a data revolution and emergence of new approaches to study transcriptomic regulation such as alternative splicing, alternative polyadenylation, CRISPR knockout screening in addition to the regular gene expression. A full characterization of the transcriptional landscape of different groups of cells or tissues holds enormous potential for both basic science as well as clinical applications. Although many methods have been developed in the realm of differential gene expression analysis, they all geared towards a particular type of sequencing data and failed to perform well when applied in different types of transcriptomic data. To fill this gap, we offer a negative beta binomial t-test (NBBt-test). NBBt-test provides multiple functions to perform differential analyses of alternative splicing, polyadenylation, CRISPR knockout screening, and gene expression datasets. Both real and large-scale simulation data show superior performance of NBBt-test with higher efficiency, and lower type I error rate and FDR to identify differential isoforms and differentially expressed genes and differential CRISPR knockout screening genes with different sample sizes when compared against the current very popular statistical methods. An R-package implementing NBBt-test is available for downloading from CRAN (https://CRAN.R-project.org/package=NBBttest).

An in vitro Assay of mRNA 3' end Using the E. coli Cell-free Expression System

At the end of about 80% of the operon in Escherichia coli, translation termination decouples transcription, leading to Rho-dependent transcription termination (RDT). However, no in vitro or in vivo assay system has proven to be good enough to see the 3' end of the mRNA generated by RDT. Here, we present a cell-free assay system that could provide detailed information on the 3' end of a transcript RNA generated by RDT. Our protocol shows how to extract transcript RNA generated by transcription reactions from a cell-free extract, followed by an RNA oligomer ligation to the 3' end of a transcript RNA of interest. The 3' end of the RNA is amplified using RT-PCR. Its genetic location can be determined using a gene-specific primer extension reaction. The 3' ends of mRNA can be visualized and quantified by polyacrylamide gel electrophoresis. One significant advantage of a cell-free assay system is that factors involved in the generation of the 3' end, such as proteins and sRNA, can be directly assayed by exogenously adding factor(s) to the reaction. Graphic abstract: An illustration of the experimental methodology.

Electrostatic Interactions between CSTF2 and pre-mRNA Drive Cleavage and Polyadenylation

Nascent pre-mRNA 3'-end cleavage and polyadenylation (C/P) involves numerous proteins that recognize multiple RNA elements. Human CSTF2 binds to a downstream U- or G/U-rich sequence through its RNA recognition motif (RRM) regulating C/P. We previously reported the only known disease-related CSTF2 RRM mutant (CSTF2^D50A) and showed that it changed the on-rate of RNA binding, leading to alternative polyadenylation in brains of mice carrying the same mutation. In this study, we further investigated the role of electrostatic interactions in the thermodynamics and kinetics of RNA binding for the CSTF2 RRM and the downstream consequences for regulation of C/P. By combining mutagenesis with NMR spectroscopy and biophysical assays, we confirmed that electrostatic attraction is the dominant factor in RRM binding to a naturally occurring U-rich RNA sequence. Moreover, we demonstrate that RNA binding is accompanied by an enthalpy-entropy compensation mechanism that is supported by changes in pico-to-nanosecond timescale RRM protein dynamics. We suggest that the dynamic binding of the RRM to U-rich RNA supports the diversity of sequences it encounters in the nucleus. Lastly, in vivo C/P assays demonstrate a competition between fast, high affinity RNA binding and efficient, correct C/P. These results highlight the importance of the surface charge of the RRM in RNA binding and the balance between nascent mRNA binding and C/P in vivo.

Implications of Poly(A) Tail Processing in Repeat Expansion Diseases

Repeat expansion diseases are a group of more than 40 disorders that affect mainly the nervous and/or muscular system and include myotonic dystrophies, Huntington’s disease, and fragile X syndrome. The mutation-driven expanded repeat tract occurs in specific genes and is composed of tri- to dodeca-nucleotide-long units. Mutant mRNA is a pathogenic factor or important contributor to the disease and has great potential as a therapeutic target. Although repeat expansion diseases are quite well known, there are limited studies concerning polyadenylation events for implicated transcripts that could have profound effects on transcript stability, localization, and translation efficiency. In this review, we briefly present polyadenylation and alternative polyadenylation (APA) mechanisms and discuss their role in the pathogenesis of selected diseases. We also discuss several methods for poly(A) tail measurement (both transcript-specific and transcriptome-wide analyses) and APA site identification—the further development and use of which may contribute to a better understanding of the correlation between APA events and repeat expansion diseases. Finally, we point out some future perspectives on the research into repeat expansion diseases, as well as APA studies.

Targeted RNA editing: novel tools to study post-transcriptional regulation

RNA binding proteins (RBPs) regulate nearly all post-transcriptional processes within cells. To fully understand RBP function, it is essential to identify their in vivo targets. Standard techniques for profiling RBP targets, such as crosslinking immunoprecipitation (CLIP) and its variants, are limited or suboptimal in some situations, e.g. when compatible antibodies are not available and when dealing with small cell populations such as neuronal subtypes and primary stem cells. This review summarizes and compares several genetic approaches recently designed to identify RBP targets in such circumstances. TRIBE (targets of RNA binding proteins identified by editing), RNA tagging, and STAMP (surveying targets by APOBEC-mediated profiling) are new genetic tools useful for the study of post-transcriptional regulation and RBP identification. We describe the underlying RNA base editing technology, recent applications, and therapeutic implications.

LncRNAs at the heart of development and disease

Long noncoding RNAs (LncRNAs) have emerged as a diverse class of functional molecules that contribute to nearly every facet of mammalian cardiac development and disease. Recent examples show that lncRNAs can be important co-regulators of cardiac patterning and morphogenesis and modulators of the pathogenic signaling that drives heart disease. The flexibility and chemical nature of RNA allows lncRNAs to utilize diverse mechanisms, mediating their effects through their sequence, structure, and molecular interactions with DNA, protein, and other RNAs. In vivo, i.e., animal, studies of individual lncRNAs highlight their ability to balance conserved cardiac gene expression networks, serve as specific and early biomarkers, and indicate their promise as useful therapeutic targets to treat human heart disease. Here, we review recent functionally characterized lncRNAs in cardiac biology and pathology and provide a perspective on emerging approaches to decipher the role of lncRNAs in the heart.

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

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

Diverse functional elements in RNA predicted transcriptome-wide by orthogonal RNA structure probing

RNA molecules can fold into complex structures and interact with trans-acting factors to control their biology. Recent methods have been focused on developing novel tools to measure RNA structure transcriptome-wide, but their utility to study and predict RNA-protein interactions or RNA processing has been limited thus far. Here, we extend these studies with the first transcriptome-wide mapping method for cataloging RNA solvent accessibility, icLASER. By combining solvent accessibility (icLASER) with RNA flexibility (icSHAPE) data, we efficiently predict RNA-protein interactions transcriptome-wide and catalog RNA polyadenylation sites by RNA structure alone. These studies showcase the power of designing novel chemical approaches to studying RNA biology. Further, our study exemplifies merging complementary methods to measure RNA structure inside cells and its utility for predicting transcriptome-wide interactions that are critical for control of and regulation by RNA structure. We envision such approaches can be applied to studying different cell types or cells under varying conditions, using RNA structure and footprinting to characterize cellular interactions and processing involving RNA.

PABPN1L assemble into "ring-like" aggregates in the cytoplasm of MII oocytes and is associated with female infertility†

Infertility affects 10% - 15% of families worldwide. However, the pathogenesis of female infertility caused by abnormal early embryonic development is not clear. A resent study showed that PABPN1L recruited BTG4 to mRNA 3'-poly(A) tails and was essential for maternal mRNA degradation. Here, we generated an PABPN1L-antibody and found "ring-like" PABPN1L aggregates in the cytoplasm of MII oocytes. PABPN1L-EGFP proteins spontaneously formed"ring-like" aggregates in vitro. This phenomenon is similar with CCR4-NOT catalytic subunit, CNOT7, when it starts deadenylation process in vitro. We constructed two mouse model (Pabpn1l  -/- and Pabpn1l  tm1a/tm1a) simulating the intron1-exon2 abnormality of human PABPN1L and found that the female was sterile and the male was fertile. Using RNA-Seq, we observed a large-scale up-regulation of RNA in zygotes derived from Pabpn1l-/- MII oocytes. We found that 9222 genes were up-regulated instead of being degraded in the Pabpn1l-♀/+♂zygote. Both the Btg4 and Cnot61 genes are necessary for the deadenylation process and Pabpn1l  -/- resembled both the Btg4 and Cnot6l knockouts, where 71.2% genes stabilized in the Btg4-♀/+♂ zygote and 84.2% genes stabilized in the Cnot6l-♀/+♂zygote were also stabilized in Pabpn1l-♀/+♂ zygote. BTG4/CNOT7/CNOT6L was partially co-located with PABPN1L in MII oocytes. The above results suggest that PABPN1L is widely associated with CCR4-NOT-mediated maternal mRNA degradation and PABPN1L variants on intron1-exon2 could be a genetic marker of female infertility. Summary sentence. "Ring-like" PABPN1L aggregates was found in the cytoplasm of MII oocytes and in vitro; intron1-exon2 abnormality of Pabpn1l leads female sterile in mice.

Dynamics in Fip1 regulate eukaryotic mRNA 3' end processing

In this study, Kumar et al. characterized the structure–function relationship of the essential poly(A) factor Fip1. Using in vitro reconstitution and structural studies, the authors report that Fip1 dynamics within the 3′ end processing machinery are required to coordinate cleavage and polyadenylation.

Cleavage and polyadenylation factor (CPF/CPSF) is a multiprotein complex essential for mRNA 3′ end processing in eukaryotes. It contains an endonuclease that cleaves pre-mRNAs, and a polymerase that adds a poly(A) tail onto the cleaved 3′ end. Several CPF subunits, including Fip1, contain intrinsically disordered regions (IDRs). IDRs within multiprotein complexes can be flexible, or can become ordered upon interaction with binding partners. Here, we show that yeast Fip1 anchors the poly(A) polymerase Pap1 onto CPF via an interaction with zinc finger 4 of another CPF subunit, Yth1. We also reconstitute a fully recombinant 850-kDa CPF. By incorporating selectively labeled Fip1 into recombinant CPF, we could study the dynamics of Fip1 within the megadalton complex using nuclear magnetic resonance (NMR) spectroscopy. This reveals that a Fip1 IDR that connects the Yth1- and Pap1-binding sites remains highly dynamic within CPF. Together, our data suggest that Fip1 dynamics within the 3′ end processing machinery are required to coordinate cleavage and polyadenylation.

A real-time fluorescence assay for CPSF73, the nuclease for pre-mRNA 3'-end processing

CPSF73 is the endonuclease that catalyzes the cleavage reaction for 3'-end processing of mRNA precursors (pre-mRNAs) in two distinct machineries, a canonical machinery for the majority of pre-mRNAs and a U7 snRNP (U7 machinery) for replication-dependent histone pre-mRNAs in animal cells. CPSF73 also possesses 5'-3' exonuclease activity in the U7 machinery, degrading the downstream cleavage product after the endonucleolytic cleavage. Recent studies show that CPSF73 is a potential target for developing anticancer, antimalarial, and antiprotozoal drugs, spurring interest in identifying new small-molecule inhibitors against this enzyme. CPSF73 nuclease activity has so far been demonstrated using a gel-based end-point assay, using radiolabeled or fluorescently labeled RNA substrates. By taking advantage of unique properties of the U7 machinery, we have developed a novel, real-time fluorescence assay for the nuclease activity of CPSF73. This assay is facile and high-throughput, and should also be helpful for the discovery of new CPSF73 inhibitors.

RiboRid: A low cost, advanced, and ultra-efficient method to remove ribosomal RNA for bacterial transcriptomics

RNA sequencing techniques have enabled the systematic elucidation of gene expression (RNA-Seq), transcription start sites (differential RNA-Seq), transcript 3′ ends (Term-Seq), and post-transcriptional processes (ribosome profiling). The main challenge of transcriptomic studies is to remove ribosomal RNAs (rRNAs), which comprise more than 90% of the total RNA in a cell. Here, we report a low-cost and robust bacterial rRNA depletion method, RiboRid, based on the enzymatic degradation of rRNA by thermostable RNase H. This method implemented experimental considerations to minimize nonspecific degradation of mRNA and is capable of depleting pre-rRNAs that often comprise a large portion of RNA, even after rRNA depletion. We demonstrated the highly efficient removal of rRNA up to a removal efficiency of 99.99% for various transcriptome studies, including RNA-Seq, Term-Seq, and ribosome profiling, with a cost of approximately $10 per sample. This method is expected to be a robust method for large-scale high-throughput bacterial transcriptomic studies.

Removal of ribosomal RNAs, a major constituent (over 90%) of cellular RNA is a critical experimental step for transcriptomic studies that deal with messenger RNAs. In this manuscript, we describe a robust method to subtract ribosomal RNA from various RNA samples. The method is based on the enzymatic degradation of target RNA by short complementary DNA and RNA:DNA duplex specific nuclease. The method comprises carefully designed experimental procedures to minimize experimental bias and unwanted removal of messenger RNAs. We validate the method on various types of transcriptomic studies for seven diverse bacterial species. This method successfully removed ribosomal RNA with over 99% of efficiency and it was comparable to commercial systems even for degraded RNA samples at a fraction of a cost.

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

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

A comprehensive map of alternative polyadenylation in African American and European American lung cancer patients

Deciphering the post-transcriptional mechanisms (PTM) regulating gene expression is critical to understand the dynamics underlying transcriptomic regulation in cancer. Alternative polyadenylation (APA)-regulation of mRNA 3'UTR length by alternating poly(A) site usage-is a key PTM mechanism whose comprehensive analysis in cancer remains an important open challenge. Here we use a method and analysis pipeline that sequences 3'end-enriched RNA directly to overcome the saturation limitation of traditional 5'-3' based sequencing. We comprehensively map the APA landscape in lung cancer in a cohort of 98 tumor/non-involved tissues derived from European American and African American patients. We identify a global shortening of 3'UTR transcripts in lung cancer, with notable functional implications on the expression of both coding and noncoding genes. We find that APA of non-coding RNA transcripts (long non-coding RNAs and microRNAs) is a recurrent event in lung cancer and discover that the selection of alternative polyA sites is a form of non-coding RNA expression control. Our results indicate that mRNA transcripts from EAs are two times more likely than AAs to undergo APA in lung cancer. Taken together, our findings comprehensively map and identify the important functional role of alternative polyadenylation in determining transcriptomic heterogeneity in lung cancer.

Star-PAP RNA Binding Landscape Reveals Novel Role of Star-PAP in mRNA Metabolism That Requires RBM10-RNA Association

Star-PAP is a non-canonical poly(A) polymerase that selects mRNA targets for polyadenylation. Yet, genome-wide direct Star-PAP targets or the mechanism of specific mRNA recognition is still vague. Here, we employ HITS-CLIP to map the cellular Star-PAP binding landscape and the mechanism of global Star-PAP mRNA association. We show a transcriptome-wide association of Star-PAP that is diminished on Star-PAP depletion. Consistent with its role in the 3'-UTR processing, we observed a high association of Star-PAP at the 3'-UTR region. Strikingly, there is an enrichment of Star-PAP at the coding region exons (CDS) in 42% of target mRNAs. We demonstrate that Star-PAP binding de-stabilises these mRNAs indicating a new role of Star-PAP in mRNA metabolism. Comparison with earlier microarray data reveals that while UTR-associated transcripts are down-regulated, CDS-associated mRNAs are largely up-regulated on Star-PAP depletion. Strikingly, the knockdown of a Star-PAP coregulator RBM10 resulted in a global loss of Star-PAP association on target mRNAs. Consistently, RBM10 depletion compromises 3'-end processing of a set of Star-PAP target mRNAs, while regulating stability/turnover of a different set of mRNAs. Our results establish a global profile of Star-PAP mRNA association and a novel role of Star-PAP in the mRNA metabolism that requires RBM10-mRNA association in the cell.

Subcellular Localization of Fad1p in Saccharomyces cerevisiae: A Choice at Post-Transcriptional Level?

FAD synthase is the last enzyme in the pathway that converts riboflavin into FAD. In Saccharomyces cerevisiae, the gene encoding for FAD synthase is FAD1, from which a sole protein product (Fad1p) is expected to be generated. In this work, we showed that a natural Fad1p exists in yeast mitochondria and that, in its recombinant form, the protein is able, per se, to both enter mitochondria and to be destined to cytosol. Thus, we propose that FAD1 generates two echoforms-that is, two identical proteins addressed to different subcellular compartments. To shed light on the mechanism underlying the subcellular destination of Fad1p, the 3' region of FAD1 mRNA was analyzed by 3'RACE experiments, which revealed the existence of (at least) two FAD1 transcripts with different 3'UTRs, the short one being 128 bp and the long one being 759 bp. Bioinformatic analysis on these 3'UTRs allowed us to predict the existence of a cis-acting mitochondrial localization motif, present in both the transcripts and, presumably, involved in protein targeting based on the 3'UTR context. Here, we propose that the long FAD1 transcript might be responsible for the generation of mitochondrial Fad1p echoform.

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

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.

Oocyte meiosis-coupled poly(A) polymerase α phosphorylation and activation trigger maternal mRNA translation in mice

Mammalian oocyte maturation is driven by strictly regulated polyadenylation and translational activation of maternal mRNA stored in the cytoplasm. However, the poly(A) polymerase (PAP) that directly mediates cytoplasmic polyadenylation in mammalian oocytes has not been determined. In this study, we identified PAPα as the elusive enzyme that catalyzes cytoplasmic mRNA polyadenylation implicated in mouse oocyte maturation. PAPα was mainly localized in the germinal vesicle (GV) of fully grown oocytes but was distributed to the ooplasm after GV breakdown. Inhibition of PAPα activity impaired cytoplasmic polyadenylation and translation of maternal transcripts, thus blocking meiotic cell cycle progression. Once an oocyte resumes meiosis, activated CDK1 and ERK1/2 cooperatively mediate the phosphorylation of three serine residues of PAPα, 537, 545 and 558, thereby leading to increased activity. This mechanism is responsible for translational activation of transcripts lacking cytoplasmic polyadenylation elements in their 3′-untranslated region (3′-UTR). In turn, activated PAPα stimulated polyadenylation and translation of the mRNA encoding its own (Papola) through a positive feedback circuit. ERK1/2 promoted Papola mRNA translation in a 3′-UTR polyadenylation signal-dependent manner. Through these mechanisms, PAPα activity and levels were significantly amplified, improving the levels of global mRNA polyadenylation and translation, thus, benefiting meiotic cell cycle progression.

Learning the Regulatory Code of Gene Expression

Data-driven machine learning is the method of choice for predicting molecular phenotypes from nucleotide sequence, modeling gene expression events including protein-DNA binding, chromatin states as well as mRNA and protein levels. Deep neural networks automatically learn informative sequence representations and interpreting them enables us to improve our understanding of the regulatory code governing gene expression. Here, we review the latest developments that apply shallow or deep learning to quantify molecular phenotypes and decode the cis-regulatory grammar from prokaryotic and eukaryotic sequencing data. Our approach is to build from the ground up, first focusing on the initiating protein-DNA interactions, then specific coding and non-coding regions, and finally on advances that combine multiple parts of the gene and mRNA regulatory structures, achieving unprecedented performance. We thus provide a quantitative view of gene expression regulation from nucleotide sequence, concluding with an information-centric overview of the central dogma of molecular biology.

Known sequence features can explain half of all human gene ends

Cleavage and polyadenylation (CPA) sites define eukaryotic gene ends. CPA sites are associated with five key sequence recognition elements: the upstream UGUA, the polyadenylation signal (PAS), and U-rich sequences; the CA/UA dinucleotide where cleavage occurs; and GU-rich downstream elements (DSEs). Currently, it is not clear whether these sequences are sufficient to delineate CPA sites. Additionally, numerous other sequences and factors have been described, often in the context of promoting alternative CPA sites and preventing cryptic CPA site usage. Here, we dissect the contributions of individual sequence features to CPA using standard discriminative models. We show that models comprised only of the five primary CPA sequence features give highest probability scores to constitutive CPA sites at the ends of coding genes, relative to the entire pre-mRNA sequence, for 41% of all human genes. U1-hybridizing sequences provide a small boost in performance. The addition of all known RBP RNA binding motifs to the model, however, increases this figure to 49%, and suggests an involvement of both known and suspected CPA regulators as well as potential new factors in delineating constitutive CPA sites. To our knowledge, this high effectiveness of established features to predict human gene ends has not previously been documented.

Reconstitution and biochemical assays of an active human histone pre-mRNA 3'-end processing machinery

In animal cells, replication-dependent histone pre-mRNAs are processed at the 3'-end by an endonucleolytic cleavage carried out by the U7 snRNP, a machinery that contains the U7 snRNA and many protein subunits. Studies on the composition of this machinery and understanding of its role in 3'-end processing were greatly facilitated by the development of an in vitro system utilizing nuclear extracts from mammalian cells 35 years ago and later from Drosophila cells. Most recently, recombinant expression and purification of the components of the machinery have enabled the full reconstitution of an active machinery and its complex with a model pre-mRNA substrate, using 13 proteins and 2 RNAs, and the determination of the structure of this active machinery. This chapter presents protocols for preparing nuclear extracts containing endogenous processing machinery, for assembling semi-recombinant and fully reconstituted machineries, and for histone pre-mRNA 3'-end processing assays with these samples.