Identification of high-confidence human poly(A) RNA isoform scaffolds using nanopore sequencing.

Nanopore sequencing devices read individual RNA strands directly. This facilitates identification of exon linkages and nucleotide modifications; however, using conventional direct RNA nanopore sequencing, the 5' and 3' ends of poly(A) RNA cannot be identified unambiguously. This is due in part to RNA degradation in vivo and in vitro that can obscure transcription start and end sites. In this study, we aimed to identify individual full-length human RNA isoforms among ∼4 million nanopore poly(A)-selected RNA reads. First, to identify RNA strands bearing 5' m7G caps, we exchanged the biological cap for a modified cap attached to a 45-nt oligomer. This oligomer adaptation method improved 5' end sequencing and ensured correct identification of the 5' m7G capped ends. Second, among these 5'-capped nanopore reads, we screened for features consistent with a 3' polyadenylation site. Combining these two steps, we identified 294,107 individual high-confidence full-length RNA scaffolds from human GM12878 cells, most of which (257,721) aligned to protein-coding genes. Of these, 4876 scaffolds indicated unannotated isoforms that were often internal to longer, previously identified RNA isoforms. Orthogonal data for m7G caps and open chromatin, such as CAGE and DNase-HS seq, confirmed the validity of these high-confidence RNA scaffolds.

3'aQTL-atlas: an atlas of 3'UTR alternative polyadenylation quantitative trait loci across human normal tissues.

Genome-wide association studies (GWAS) have identified thousands of non-coding single-nucleotide polymorphisms (SNPs) associated with human traits and diseases. However, functional interpretation of these SNPs remains a significant challenge. Our recent study established the concept of 3' untranslated region (3'UTR) alternative polyadenylation (APA) quantitative trait loci (3'aQTLs), which can be used to interpret ∼16.1% of GWAS SNPs and are distinct from gene expression QTLs and splicing QTLs. Despite the growing interest in 3'aQTLs, there is no comprehensive database for users to search and visualize them across human normal tissues. In the 3'aQTL-atlas (https://wlcb.oit.uci.edu/3aQTLatlas), we provide a comprehensive list of 3'aQTLs containing ∼1.49 million SNPs associated with APA of target genes, based on 15,201 RNA-seq samples across 49 human Genotype-Tissue Expression (GTEx v8) tissues isolated from 838 individuals. The 3'aQTL-atlas provides a ∼2-fold increase in sample size compared with our published study. It also includes 3'aQTL searches by Gene/SNP across tissues, a 3'aQTL genome browser, 3'aQTL boxplots, and GWAS-3'aQTL colocalization event visualization. The 3'aQTL-atlas aims to establish APA as an emerging molecular phenotype to explain a large fraction of GWAS risk SNPs, leading to significant novel insights into the genetic basis of APA and APA-linked susceptibility genes in human traits and diseases.

scAPAtrap: identification and quantification of alternative polyadenylation sites from single-cell RNA-seq data.

Alternative polyadenylation (APA) generates diverse mRNA isoforms, which contributes to transcriptome diversity and gene expression regulation by affecting mRNA stability, translation and localization in cells. The rapid development of 3' tag-based single-cell RNA-sequencing (scRNA-seq) technologies, such as CEL-seq and 10x Genomics, has led to the emergence of computational methods for identifying APA sites and profiling APA dynamics at single-cell resolution. However, existing methods fail to detect the precise location of poly(A) sites or sites with low read coverage. Moreover, they rely on priori genome annotation and can only detect poly(A) sites located within or near annotated genes. Here we proposed a tool called scAPAtrap for detecting poly(A) sites at the whole genome level in individual cells from 3' tag-based scRNA-seq data. scAPAtrap incorporates peak identification and poly(A) read anchoring, enabling the identification of the precise location of poly(A) sites, even for sites with low read coverage. Moreover, scAPAtrap can identify poly(A) sites without using priori genome annotation, which helps locate novel poly(A) sites in previously overlooked regions and improve genome annotation. We compared scAPAtrap with two latest methods, scAPA and Sierra, using scRNA-seq data from different experimental technologies and species. Results show that scAPAtrap identified poly(A) sites with higher accuracy and sensitivity than competing methods and could be used to explore APA dynamics among cell types or the heterogeneous APA isoform expression in individual cells. scAPAtrap is available at https://github.com/BMILAB/scAPAtrap.

Breaking the code of polyadenylation-induced translation.

The translation of many maternal mRNAs is regulated by dynamic changes in poly(A) tail length. During maturation of Xenopus oocytes, polyadenylation is mediated by three different cis elements in the 3' untranslated region (UTR) of maternal mRNAs. In this issue, Piqué et al. (2008) explore the interplay of these elements to elucidate a combinatorial code that predicts the timing of polyadenylation and translation of maternal mRNAs.

A combinatorial code for CPE-mediated translational control.

Cytoplasmic polyadenylation plays a key role in the translational control of mRNAs driving biological processes such as gametogenesis, cell-cycle progression, and synaptic plasticity. What determines the distinct time of polyadenylation and extent of translational control of a given mRNA, however, is poorly understood. The polyadenylation-regulated translation is controlled by the cytoplasmic polyadenylation element (CPE) and its binding protein, CPEB, which can assemble both translational repression or activation complexes. Using a combination of mutagenesis and experimental validation of genome-wide computational predictions, we show that the number and relative position of two elements, the CPE and the Pumilio-binding element, with respect to the polyadenylation signal define a combinatorial code that determines whether an mRNA will be translationally repressed by CPEB, as well as the extent and time of cytoplasmic polyadenylation-dependent translational activation.

Systematic identification of functional SNPs interrupting 3'UTR polyadenylation signals.

Alternative polyadenylation (APA) is emerging as a widespread regulatory layer since the majority of human protein-coding genes contain several polyadenylation (p(A)) sites in their 3'UTRs. By generating isoforms with different 3'UTR length, APA potentially affects mRNA stability, translation efficiency, nuclear export, and cellular localization. Polyadenylation sites are regulated by adjacent RNA cis-regulatory elements, the principals among them are the polyadenylation signal (PAS) AAUAAA and its main variant AUUAAA, typically located ~20-nt upstream of the p(A) site. Mutations in PAS and other auxiliary poly(A) cis-elements in the 3'UTR of several genes have been shown to cause human Mendelian diseases, and to date, only a few common SNPs that regulate APA were associated with complex diseases. Here, we systematically searched for SNPs that affect gene expression and human traits by modulation of 3'UTR APA. First, focusing on the variants most likely to exert the strongest effect, we identified 2,305 SNPs that interrupt the canonical PAS or its main variant. Implementing pA-QTL tests using GTEx RNA-seq data, we identified 330 PAS SNPs (called PAS pA-QTLs) that were significantly associated with the usage of their p(A) site. As expected, PAS-interrupting alleles were mostly linked with decreased cleavage at their p(A) site and the consequential 3'UTR lengthening. However, interestingly, in ~10% of the cases, the PAS-interrupting allele was associated with increased usage of an upstream p(A) site and 3'UTR shortening. As an indication of the functional effects of these PAS pA-QTLs on gene expression and complex human traits, we observed for few dozens of them marked colocalization with eQTL and/or GWAS signals. The PAS-interrupting alleles linked with 3'UTR lengthening were also strongly associated with decreased gene expression, indicating that shorter isoforms generated by APA are generally more stable than longer ones. Last, we carried out an extended, genome-wide analysis of 3'UTR variants and detected thousands of additional pA-QTLs having weaker effects compared to the PAS pA-QTLs.

Forced isoform switching of Neat1_1 to Neat1_2 leads to the loss of Neat1_1 and the hyperformation of paraspeckles but does not affect the development and growth of mice.

Neat1 is a long noncoding RNA (lncRNA) that serves as an architectural component of the nuclear bodies known as paraspeckles. Two isoforms of Neat1, the short isoform Neat1_1 and the long isoform Neat1_2, are generated from the same gene locus by alternative 3' processing. Neat1_1 is the most abundant and the best conserved isoform expressed in various cell types, whereas Neat1_2 is expressed in a small population of particular cell types, including the tip cells of the intestinal epithelium. To investigate the physiological significance of isoform switching, we created mutant mice that solely expressed Neat1_2 by deleting the upstream polyadenylation (poly-A) signal (PAS) required for the production of Neat1_1. We observed the loss of Neat1_1 and strong up-regulation of Neat1_2 in various tissues and cells and the subsequent hyperformation of paraspeckles, especially in cells that normally express Neat1_2. However, the mutant mice were born at the expected Mendelian ratios and did not exhibit obvious external and histological abnormalities. These observations suggested that the hyperformation of paraspeckles does not interfere with the development and growth of these animals under normal laboratory conditions.

Transcriptome Analyses of FY Mutants Reveal Its Role in mRNA Alternative Polyadenylation.

A crucial step for mRNA polyadenylation is poly(A) signal recognition by trans-acting factors. The mammalian cleavage and polyadenylation specificity factor (CPSF) complex components CPSF30 and WD repeat-containing protein33 (WDR33) recognize the canonical AAUAAA for polyadenylation. In Arabidopsis (Arabidopsis thaliana), the flowering time regulator FY is the homolog of WDR33. However, its role in mRNA polyadenylation is poorly understood. Using poly(A) tag sequencing, we found that >50% of alternative polyadenylation (APA) events are altered in fy single mutants or double mutants with oxt6 (a null mutant of AtCPSF30), but mutation of the FY WD40-repeat has a stronger effect than deletion of the plant-unique Pro-Pro-Leu-Pro-Pro (PPLPP) domain. fy mutations disrupt AAUAAA or AAUAAA-like poly(A) signal recognition. Notably, A-rich signal usage is suppressed in the WD40-repeat mutation but promoted in PPLPP-domain deficiency. However, fy mutations do not aggravate the altered signal usage in oxt6 Furthermore, the WD40-repeat mutation shows a preference for 3' untranslated region shortening, but the PPLPP-domain deficiency shows a preference for lengthening. Interestingly, the WD40-repeat mutant exhibits shortened primary roots and late flowering with alteration of APA of related genes. Importantly, the long transcripts of two APA genes affected in fy are related to abiotic stress responses. These results reveal a conserved and specific role of FY in mRNA polyadenylation.

Determinants of heterochromatic siRNA biogenesis and function.

Endogenous small interfering RNAs (siRNAs) and other classes of small RNA provide the specificity signals for silencing of transposons and repeated DNA elements at the posttranscriptional and transcriptional levels. However, the determinants that define an siRNA-producing region or control the silencing function of siRNAs are poorly understood. Here we show that convergent antisense transcription and availability of the Dicer ribonuclease are the key determinants for primary siRNA generation. Surprisingly, Dicer makes dual contributions to heterochromatin formation, promoting histone H3 lysine 9 methylation independently of its catalytic activity, in addition to its well-known role in catalyzing siRNA generation. Furthermore, sequences in the 3' UTR of an mRNA-coding gene inhibit the ability of siRNAs to promote heterochromatin formation, providing another layer of control that prevents the silencing of protein-coding RNAs. Our results reveal distinct mechanisms that limit siRNA generation to centromeric DNA repeats and prevent spurious siRNA-mediated silencing at euchromatic loci.

Crosstalk between alternative polyadenylation and miRNAs in the regulation of protein translational efficiency.

3' UTRs play important roles in the gene regulation network via their influence on mRNA stability, translational efficiency, and subcellular localization. For a given gene, 3' UTRs of different lengths generated by alternative polyadenylation (APA) may result in functional differences in regulation. The mechanistic details of how length changes of 3' UTRs alter gene function remain unclear. By combining APA sequencing and polysome profiling, we observed that mRNA isoforms with shorter 3' UTRs were bound with more polysomes in six cell lines but not in NIH3T3 cells, suggesting that changing 3' UTRs to shorter isoforms may lead to a higher gene translational efficiency. By interfering with the expression of TNRC6A and analyzing AGO2-PAR-CLIP data, we revealed that the APA effect on translational efficiency was mainly regulated by miRNAs, and this regulation was cell cycle dependent. The discrepancy between NIH3T3 and other cell lines was due to contact inhibition of NIH3T3. Thus, the crosstalk between APA and miRNAs may be needed for the regulation of protein translational efficiency.

Analysis of SINE Families B2, Dip, and Ves with Special Reference to Polyadenylation Signals and Transcription Terminators.

Short Interspersed Elements (SINEs) are eukaryotic non-autonomous retrotransposons transcribed by RNA polymerase III (pol III). The 3'-terminus of many mammalian SINEs has a polyadenylation signal (AATAAA), pol III transcription terminator, and A-rich tail. The RNAs of such SINEs can be polyadenylated, which is unique for pol III transcripts. Here, B2 (mice and related rodents), Dip (jerboas), and Ves (vespertilionid bats) SINE families were thoroughly studied. They were divided into subfamilies reliably distinguished by relatively long indels. The age of SINE subfamilies can be estimated, which allows us to reconstruct their evolution. The youngest and most active variants of SINE subfamilies were given special attention. The shortest pol III transcription terminators are TCTTT (B2), TATTT (Ves and Dip), and the rarer TTTT. The last nucleotide of the terminator is often not transcribed; accordingly, the truncated terminator of its descendant becomes nonfunctional. The incidence of complete transcription of the TCTTT terminator is twice higher compared to TTTT and thus functional terminators are more likely preserved in daughter SINE copies. Young copies have long poly(A) tails; however, they gradually shorten in host generations. Unexpectedly, the tail shortening below A10 increases the incidence of terminator elongation by Ts thus restoring its efficiency. This process can be critical for the maintenance of SINE activity in the genome.

Dicer-2 promotes mRNA activation through cytoplasmic polyadenylation.

Cytoplasmic polyadenylation is a widespread mechanism to regulate mRNA translation. In vertebrates, this process requires two sequence elements in target 3' UTRs: the U-rich cytoplasmic polyadenylation element and the AAUAAA hexanucleotide. In Drosophila melanogaster, cytoplasmic polyadenylation of Toll mRNA occurs independently of these canonical elements and requires a machinery that remains to be characterized. Here we identify Dicer-2 as a component of this machinery. Dicer-2, a factor previously involved in RNA interference (RNAi), interacts with the cytoplasmic poly(A) polymerase Wispy. Depletion of Dicer-2 from polyadenylation-competent embryo extracts and analysis of wispy mutants indicate that both factors are necessary for polyadenylation and translation of Toll mRNA. We further identify r2d2 mRNA, encoding a Dicer-2 partner in RNAi, as a Dicer-2 polyadenylation target. Our results uncover a novel function of Dicer-2 in activation of mRNA translation through cytoplasmic polyadenylation.

Pumilio directs deadenylation-associated translational repression of the cyclin-dependent kinase 1 activator RGC-32.

Response gene to complement-32 (RGC-32) activates cyclin-dependent kinase 1, regulates the cell cycle and is deregulated in many human tumours. We previously showed that RGC-32 expression is upregulated by the cancer-associated Epstein-Barr virus (EBV) in latently infected B cells through the relief of translational repression. We now show that EBV infection of naïve primary B cells also induces RGC-32 protein translation. In EBV-immortalised cell lines, we found that RGC-32 depletion resulted in cell death, indicating a key role in B cell survival. Studying RGC-32 translational control in EBV-infected cells, we found that the RGC-32 3'untranslated region (3'UTR) mediates translational repression. Repression was dependent on a single Pumilio binding element (PBE) adjacent to the polyadenylation signal. Mutation of this PBE did not affect mRNA cleavage, but resulted in increased polyA tail length. Consistent with Pumilio-dependent recruitment of deadenylases, we found that depletion of Pumilio in EBV-infected cells increased RGC-32 protein expression and polyA tail length. The extent of Pumilio binding to the endogenous RGC-32 mRNA in EBV-infected cell lines also correlated with RGC-32 protein expression. Our data demonstrate the importance of RGC-32 for the survival of EBV-immortalised B cells and identify Pumilio as a key regulator of RGC-32 translation.

An allelic series at the KARRIKIN INSENSITIVE 2 locus of Arabidopsis thaliana decouples ligand hydrolysis and receptor degradation from downstream signalling.

Karrikins are butenolide compounds present in post-fire environments that can stimulate seed germination in many species, including Arabidopsis thaliana. Plants also produce endogenous butenolide compounds that serve as hormones, namely strigolactones (SLs). The receptor for karrikins (KARRIKIN INSENSITIVE 2; KAI2) and the receptor for SLs (DWARF14; D14) are homologous proteins that share many similarities. The mode of action of D14 as a dual enzyme receptor protein is well established, but the nature of KAI2-dependent signalling and its function as a receptor are not fully understood. To expand our knowledge of how KAI2 operates, we screened ethyl methanesulphonate (EMS)-mutagenized populations of A. thaliana for mutants with kai2-like phenotypes and isolated 13 new kai2 alleles. Among these alleles, kai2-10 encoded a D184N protein variant that was stable in planta. Differential scanning fluorimetry assays indicated that the KAI2 D184N protein could interact normally with bioactive ligands. We developed a KAI2-active version of the fluorescent strigolactone analogue Yoshimulactone Green to show that KAI2 D184N exhibits normal rates of ligand hydrolysis. KAI2 D184N degraded in response to treatment with exogenous ligands, suggesting that receptor degradation is a consequence of ligand binding and hydrolysis, but is insufficient for signalling activity. Remarkably, KAI2 D184N degradation was hypersensitive to karrikins, but showed a normal response to strigolactone analogues, implying that these butenolides may interact differently with KAI2. These results demonstrate that the enzymatic and signalling functions of KAI2 can be decoupled, and provide important insights into the mechanistic events that underpin butenolide signalling in plants.

Genome-wide dynamics of alternative polyadenylation in rice.

Alternative polyadenylation (APA), in which a transcript uses one of the poly(A) sites to define its 3'-end, is a common regulatory mechanism in eukaryotic gene expression. However, the potential of APA in determining crop agronomic traits remains elusive. This study systematically tallied poly(A) sites of 14 different rice tissues and developmental stages using the poly(A) tag sequencing (PAT-seq) approach. The results indicate significant involvement of APA in developmental and quantitative trait loci (QTL) gene expression. About 48% of all expressed genes use APA to generate transcriptomic and proteomic diversity. Some genes switch APA sites, allowing differentially expressed genes to use alternate 3' UTRs. Interestingly, APA in mature pollen is distinct where differential expression levels of a set of poly(A) factors and different distributions of APA sites are found, indicating a unique mRNA 3'-end formation regulation during gametophyte development. Equally interesting, statistical analyses showed that QTL tends to use APA for regulation of gene expression of many agronomic traits, suggesting a potential important role of APA in rice production. These results provide thus far the most comprehensive and high-resolution resource for advanced analysis of APA in crops and shed light on how APA is associated with trait formation in eukaryotes.

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

BACKGROUND: Polyadenylation influences many aspects of mRNA as well as viral RNA. variable polyadenylation at the 3' end have been reported in RNA viruses. It is interesting to identify the characteristic and potential role of 3' polyadenylation of Wheat yellow mosaic virus (WYMV), which has been reported to contain two genomic RNAs with 3' poly(A) tails and caused severe disease on wheat in East Asia region. METHODS: 3' RACE was used to identify sequences of the 3' end in WYMV RNAs from naturally infected wheat by WYMV. In vitro translation assay was performed to analyze effect of UTRs of WYMV with or without 3'polyadenylation on translation. In vitro replication mediated by WYMV NIb protein were performed to evaluate effect of variable polyadenylation on replication. RESULTS: Variable polyadenylation in WYMV RNAs was identified via 3' RACE. WYMV RNAs in naturally infected wheat in China simultaneously present with regions of long, short, or no adenylation at the 3' ends. The effects of variable polyadenylation on translation and replication of WYMV RNAs were evaluated. 5'UTR and 3'UTR of WYMV RNA1 or RNA2 synergistically enhanced the translation of the firefly luciferase (Fluc) gene in in vitro WGE system, whereas additional adenylates had an oppositive effect on this enhancement on translation mediated by UTRs of WYMV. Additional adenylates remarkably inhibited the synthesis of complementary strand from viral genome RNA during the in vitro replication mediated by WYMV NIb protein. CONCLUSIONS: 3' end of WYMV RNAs present variable polyadenylation even no polyadenylation. 3' polyadenylation have opposite effect on translation mediated by UTRs of WYMV RNA1 or RNA2. 3' polyadenylation have negative effect on minus-strand synthesis of WYMV RNA in vitro. Variable polyadenylation of WYMV RNAs may provide sufficient selection on the template for translation and replication.

Therapeutic Genome Editing for Myotonic Dystrophy Type 1 Using CRISPR/Cas9.

Myotonic dystrophy type 1 (DM1) is caused by a CTG nucleotide repeat expansion within the 3' UTR of the Dystrophia Myotonica protein kinase gene. In this study, we explored therapeutic genome editing using CRISPR/Cas9 via targeted deletion of expanded CTG repeats and targeted insertion of polyadenylation signals in the 3' UTR upstream of the CTG repeats to eliminate toxic RNA CUG repeats. We found paired SpCas9 or SaCas9 guide RNA induced deletion of expanded CTG repeats. However, this approach incurred frequent inversion in both the mutant and normal alleles. In contrast, the insertion of polyadenylation signals in the 3' UTR upstream of the CTG repeats eliminated toxic RNA CUG repeats, which led to phenotype reversal in differentiated neural stem cells, forebrain neurons, cardiomyocytes, and skeletal muscle myofibers. We concluded that targeted insertion of polyadenylation signals in the 3' UTR is a viable approach to develop therapeutic genome editing for DM1.

Cleavage and polyadenylation specificity factor 30: An RNA-binding zinc-finger protein with an unexpected 2Fe-2S cluster.

Cleavage and polyadenylation specificity factor 30 (CPSF30) is a key protein involved in pre-mRNA processing. CPSF30 contains five Cys3His domains (annotated as "zinc-finger" domains). Using inductively coupled plasma mass spectrometry, X-ray absorption spectroscopy, and UV-visible spectroscopy, we report that CPSF30 is isolated with iron, in addition to zinc. Iron is present in CPSF30 as a 2Fe-2S cluster and uses one of the Cys3His domains; 2Fe-2S clusters with a Cys3His ligand set are rare and notably have also been identified in MitoNEET, a protein that was also annotated as a zinc finger. These findings support a role for iron in some zinc-finger proteins. Using electrophoretic mobility shift assays and fluorescence anisotropy, we report that CPSF30 selectively recognizes the AU-rich hexamer (AAUAAA) sequence present in pre-mRNA, providing the first molecular-based evidence to our knowledge for CPSF30/RNA binding. Removal of zinc, or both zinc and iron, abrogates binding, whereas removal of just iron significantly lessens binding. From these data we propose a model for RNA recognition that involves a metal-dependent cooperative binding mechanism.

Ending the message: poly(A) signals then and now.

Polyadenylation [poly(A)] signals (PAS) are a defining feature of eukaryotic protein-coding genes. The central sequence motif AAUAAA was identified in the mid-1970s and subsequently shown to require flanking, auxiliary elements for both 3'-end cleavage and polyadenylation of premessenger RNA (pre-mRNA) as well as to promote downstream transcriptional termination. More recent genomic analysis has established the generality of the PAS for eukaryotic mRNA. Evidence for the mechanism of mRNA 3'-end formation is outlined, as is the way this RNA processing reaction communicates with RNA polymerase II to terminate transcription. The widespread phenomenon of alternative poly(A) site usage and how this interrelates with pre-mRNA splicing is then reviewed. This shows that gene expression can be drastically affected by how the message is ended. A central theme of this review is that while genomic analysis provides generality for the importance of PAS selection, detailed mechanistic understanding still requires the direct analysis of specific genes by genetic and biochemical approaches.

Inflammatory facioscapulohumeral muscular dystrophy type 2 in 18p deletion syndrome.

Facioscapulohumeral muscular dystrophy (FSHD) has been shown to be related to genetic and epigenetic derepression of DUX4 (mapping to chromosome 4), a gene located within a repeat array of D4Z4 sequences of polymorphic length. FSHD type 1 (FSHD1) is associated with pathogenic D4Z4 repeat array contraction, while FSHD type 2 (FSHD2) is associated with SMCHD1 variants (a chromatin modifier gene that maps to the short arm of chromosome 18). Both FSHD types require permissive polyadenylation signal (4qA) downstream of the D4Z4 array.