Genome-wide control of polyadenylation site choice by CPSF30 in Arabidopsis

Experimental Verification and Evolutionary Origin of 5'-UTR Polyadenylation Sites in Arabidopsis thaliana

Messenger RNA (mRNA) polyadenylation is an indispensable step during post-transcriptional pre-mRNA processing for most genes in eukaryotes. The usage of one poly(A) site over another is known as alternative polyadenylation (APA). APA has been implicated in gene expression regulation through its role of selecting the ends of a transcript. Recent studies of polyadenylation profiles in the Arabidopsis database unexpectedly predicted that a portion of the poly(A) sites are located in the 5'-UTR, which remains to be experimentally verified. We selected 16 genes from a dataset of 744, based on criteria designed to minimize problems in interpretation. Here, we experimentally verify 5'-UTR-APA in Arabidopsis for 10 of the 16 selected genes, and show for the first time existence of independent polyadenylated 5'-UTR transcripts, arising due to alternative polyadenylation. We used 3'-RACE and sequencing to validate poly(A) sites and northern blot to show that the observed short upstream transcripts do not arise from the 3'-end of a previously unrecognized convergent gene. Evidence is reported showing that two of the independent upstream open reading frame (uORF) transcripts studied, one containing a complex dual uORF, very likely arose by exon shuffling following duplication of the 5'-end from the downstream major open reading frame (mORF). Finally, results are presented to show that the uORF in this gene may encode two short functional proteins, based on observation of amino acid sequence conservation encoded by the dual uORFs.

Transcript-level expression control of plant NLR genes

Plant NLR genes encode sensitive immune receptors that can mediate the specific recognition of pathogen avirulence effectors and activate a strong defence response, termed effector-triggered immunity. The expression of NLRs requires strict regulation, as their ability to trigger immunity is dependent on their dose, and overexpression of NLRs results in autoimmunity and massive fitness costs. An elaborate interplay of different mechanisms controlling NLR transcript levels allows plants to maximize their defence capacity, whilst limiting negative impact on their fitness. Global suppression of NLR transcripts may be a prerequisite for the fast evolution of new NLR variants and the expansion of this gene family. Here, we summarize recent progress made towards a comprehensive understanding of NLR transcript-level expression control. Multiple mechanistic steps, including transcription as well as co-/post-transcriptional processing and transcript turn-over, contribute to balanced base levels of NLR transcripts and allow for dynamic adjustments to defence situations.

UTR-Dependent Control of Gene Expression in Plants

Throughout their lives, plants sense many developmental and environmental stimuli, and activation of optimal responses against these stimuli requires extensive transcriptional reprogramming. To facilitate this activation, plant mRNA contains untranslated regions (UTRs) that significantly increase the coding capacity of the genome by producing multiple mRNA variants from the same gene. In this review we compare UTRs of arabidopsis (Arabidopsis thaliana) and rice (Oryza sativum) at the genome scale to highlight their complexity in crop plants. We discuss different modes of UTR-based regulation with emphasis on genes that regulate multiple plant processes, including flowering, stress responses, and nutrient homeostasis. We demonstrate functional specificity in genes with variable UTR length and propose future research directions.

Advances in Transcriptomics of Plants

The current global population of 7.3 billion is estimated to reach 9.7 billion in the year 2050. Rapid population growth is driving up global food demand. Additionally, global climate change, environmental degradation, drought, emerging diseases, and salty soils are the current threats to global food security. In order to mitigate the adverse effects of these diverse agricultural productivity constraints and enhance crop yield and stress-tolerance in plants, we need to go beyond traditional and molecular plant breeding. The powerful new tools for genome editing, Transcription Activator-Like Effector Nucleases (TALENs) and Clustered Regulatory Interspaced Short Palindromic Repeats (CRISPR)/Cas systems (CRISPR-Cas9), have been hailed as a quantum leap forward in the development of stress-resistant plants. Plant breeding techniques, however, have several drawbacks. Hence, identification of transcriptional regulatory elements and deciphering mechanisms underlying transcriptional regulation are crucial to avoiding unintended consequences in modified crop plants, which could ultimately have negative impacts on human health. RNA splicing as an essential regulated post-transcriptional process, alternative polyadenylation as an RNA-processing mechanism, along with non-coding RNAs (microRNAs, small interfering RNAs and long non-coding RNAs) have been identified as major players in gene regulation. In this chapter, we highlight new findings on the essential roles of alternative splicing and alternative polyadenylation in plant development and response to biotic and abiotic stresses. We also discuss biogenesis and the functions of microRNAs (miRNAs) and small interfering RNAs (siRNAs) in plants and recent advances in our knowledge of the roles of miRNAs and siRNAs in plant stress response. Graphical Abstract.

The Arabidopsis CPSF30-L gene plays an essential role in nitrate signaling and regulates the nitrate transceptor gene NRT1.1

Plants have evolved sophisticated mechanisms to adapt to fluctuating environmental nitrogen availability. However, more underlying genes regulating the response to nitrate have yet to be characterized. We report here the identification of a nitrate regulatory mutant whose mutation mapped to the Cleavage and Polyadenylation Specificity Factor 30 gene (CPSF30-L). In the mutant, induction of nitrate-responsive genes was inhibited independent of the ammonium conditions and was restored by expression of the wild-type 65 kDa encoded by CPSF30-L. Molecular and genetic evidence suggests that CPSF30-L works upstream of NRT1.1 and independently of NLP7 in response to nitrate. Analysis of the 3'-UTR of NRT1.1 showed that the pattern of polyadenylation sites was altered in the cpsf30 mutant. Transcriptome analysis revealed that four nitrogen-related clusters were enriched in the differentially expressed genes of the cpsf30 mutant. Nitrate uptake was decreased in the mutant along with reduced expression of the nitrate transporter/sensor gene NRT1.1, while nitrate reduction and amino acid content were enhanced in roots along with increased expression of several nitrate assimilatory genes. These findings indicate that the 65 kDa protein encoded by CPSF30-L mediates nitrate signaling in part by regulating NRT1.1 expression, thus adding an important component to the nitrate signaling network.

Poly(A)-ClickSeq: click-chemistry for next-generation 3΄-end sequencing without RNA enrichment or fragmentation

The recent emergence of alternative polyadenylation (APA) as an engine driving transcriptomic diversity has stimulated the development of sequencing methodologies designed to assess genome-wide polyadenylation events. The goal of these approaches is to enrich, partition, capture and ultimately sequence poly(A) site junctions. However, these methods often require poly(A) enrichment, 3΄ linker ligation steps, and RNA fragmentation, which can necessitate higher levels of starting RNA, increase experimental error and potentially introduce bias. We recently reported a click-chemistry based method for generating RNAseq libraries called ‘ClickSeq’. Here, we adapt this method to direct the cDNA synthesis specifically toward the 3΄UTR/poly(A) tail junction of cellular RNA. With this novel approach, we demonstrate sensitive and specific enrichment for poly(A) site junctions without the need for complex sample preparation, fragmentation or purification. Poly(A)-ClickSeq (PAC-seq) is therefore a simple procedure that generates high-quality RNA-seq poly(A) libraries. As a proof-of-principle, we utilized PAC-seq to explore the poly(A) landscape of both human and Drosophila cells in culture and observed outstanding overlap with existing poly(A) databases and also identified previously unannotated poly(A) sites. Moreover, we utilize PAC-seq to quantify and analyze APA events regulated by CFIm25 illustrating how this technology can be harnessed to identify alternatively polyadenylated RNA.

Role of cleavage and polyadenylation specificity factor 100: anchoring poly(A) sites and modulating transcription termination

CPSF100 is a core component of the cleavage and polyadenylation specificity factor (CPSF) complex for 3'-end formation of mRNA, but it still has no clear functional assignment. CPSF100 was reported to play a role in RNA silencing and promote flowering in Arabidopsis. However, the molecular mechanisms underlying these phenomena are not fully understood. Our genetics analyses indicate that plants with a hypomorphic mutant of CPSF100 (esp5) show defects in embryogenesis, reduced seed production or altered root morphology. To unravel this puzzle, we employed a poly(A) tag sequencing protocol and uncovered a different poly(A) profile in esp5. This transcriptome-wide analysis revealed alternative polyadenylation of thousands of genes, most of which result in transcriptional read-through in protein-coding genes. AtCPSF100 also affects poly(A) signal recognition on the far-upstream elements; in particular it prefers less U-rich sequences. Importantly, AtCPSF100 was found to exert its functions through the change of poly(A) sites on genes encoding binding proteins, such as nucleotide-binding, RNA-binding and poly(U)-binding proteins. In addition, through its interaction with RNA Polymerase II C-terminal domain (CTD) and affecting the expression level of CTD phosphatase-like 3 (CPL3), AtCPSF100 is shown to potentially ensure transcriptional termination by dephosphorylation of Ser2 on the CTD. These data suggest a key role for CPSF100 in locating poly(A) sites and affecting transcription termination.

PAcluster: Clustering polyadenylation site data using canonical correlation analysis

Alternative polyadenylation (APA) is a pervasive mechanism that contributes to gene regulation. Increasing sequenced poly(A) sites are placing new demands for the development of computational methods to investigate APA regulation. Cluster analysis is important to identify groups of co-expressed genes. However, clustering of poly(A) sites has not been extensively studied in APA, where most APA studies failed to consider the distribution, abundance, and variation of APA sites in each gene. Here we constructed a two-layer model based on canonical correlation analysis (CCA) to explore the underlying biological mechanisms in APA regulation. The first layer quantifies the general correlation of APA sites across various conditions between each gene and the second layer identifies genes with statistically significant correlation on their APA patterns to infer APA-specific gene clusters. Using hierarchical clustering, we comprehensively compared our method with four other widely used distance measures based on three performance indexes. Results showed that our method significantly enhanced the clustering performance for both synthetic and real poly(A) site data and could generate clusters with more biological meaning. We have implemented the CCA-based method as a publically available R package called PAcluster, which provides an efficient solution to the clustering of large APA-specific biological dataset.

Cleavage and polyadenylation: Ending the message expands gene regulation

Cleavage and polyadenylation (pA) is a fundamental step that is required for the maturation of primary protein encoding transcripts into functional mRNAs that can be exported from the nucleus and translated in the cytoplasm. 3'end processing is dependent on the assembly of a multiprotein processing complex on the pA signals that reside in the pre-mRNAs. Most eukaryotic genes have multiple pA signals, resulting in alternative cleavage and polyadenylation (APA), a widespread phenomenon that is important to establish cell state and cell type specific transcriptomes. Here, we review how pA sites are recognized and comprehensively summarize how APA is regulated and creates mRNA isoform profiles that are characteristic for cell types, tissues, cellular states and disease.

Noncanonical Alternative Polyadenylation Contributes to Gene Regulation in Response to Hypoxia

Stresses from various environmental challenges continually confront plants, and their responses are important for growth and survival. One molecular response to such challenges involves the alternative polyadenylation of mRNA. In plants, it is unclear how stress affects the production and fate of alternative mRNA isoforms. Using a genome-scale approach, we show that in Arabidopsis thaliana, hypoxia leads to increases in the number of mRNA isoforms with polyadenylated 3' ends that map to 5'-untranslated regions (UTRs), introns, and protein-coding regions. RNAs with 3' ends within protein-coding regions and introns were less stable than mRNAs that end at 3'-UTR poly(A) sites. Additionally, these RNA isoforms were underrepresented in polysomes isolated from control and hypoxic plants. By contrast, mRNA isoforms with 3' ends that lie within annotated 5'-UTRs were overrepresented in polysomes and were as stable as canonical mRNA isoforms. These results indicate that the generation of noncanonical mRNA isoforms is an important feature of the abiotic stress response. The finding that several noncanonical mRNA isoforms are relatively unstable suggests that the production of non-stop and intronic mRNA isoforms may represent a form of negative regulation in plants, providing a conceptual link with mechanisms that generate these isoforms (such as alternative polyadenylation) and RNA surveillance.

Dehydration stress extends mRNA 3' untranslated regions with noncoding RNA functions in Arabidopsis

The 3′ untranslated regions (3′ UTRs) of mRNAs play important roles in the regulation of mRNA localization, translation, and stability. Alternative cleavage and polyadenylation (APA) generates mRNAs with different 3′ UTRs, but the involvement of this process in stress response has not yet been clarified. Here, we report that a subset of stress-related genes exhibits 3′ UTR extensions of their mRNAs during dehydration stress. These extended 3′ UTRs have characteristics of long noncoding RNAs and likely do not interact with miRNAs. Functional studies using T-DNA insertion mutants reveal that they can act as antisense transcripts to repress expression levels of sense genes from the opposite strand or can activate the transcription or lead to read-through transcription of their downstream genes. Further analysis suggests that transcripts with 3′ UTR extensions have weaker poly(A) signals than those without 3′ UTR extensions. Finally, we show that their biogenesis is partially dependent on a trans-acting factor FPA. Taken together, we report that dehydration stress could induce transcript 3′ UTR extensions and elucidate a novel function for these stress-induced 3′ UTR extensions as long noncoding RNAs in the regulation of their neighboring genes.

Interplay between Alternative Splicing and Alternative Polyadenylation Defines the Expression Outcome of the Plant Unique OXIDATIVE TOLERANT-6 Gene

Pre-mRNA alternative splicing and alternative polyadenylation have been implicated to play important roles during eukaryotic gene expression. However, much remains unknown regarding the regulatory mechanisms and the interactions of these two processes in plants. Here we focus on an Arabidopsis gene OXT6 (Oxidative Tolerant-6) that has been demonstrated to encode two proteins through alternative splicing and alternative polyadenylation. Specifically, alternative polyadenylation at Intron-2 of OXT6 produces a transcript coding for AtCPSF30, an Arabidopsis ortholog of 30 kDa subunit of the Cleavage and Polyadenylation Specificity Factor. On the other hand, alternative splicing of Intron-2 generates a longer transcript encoding a protein named AtC30Y, a polypeptide including most part of AtCPSF30 and a YT521B domain. To investigate the expression outcome of OXT6 in plants, a set of mutations were constructed to alter the splicing and polyadenylation patterns of OXT6. Analysis of transgenic plants bearing these mutations by quantitative RT-PCR revealed a competition relationship between these two processes. Moreover, when both splice sites and poly(A) signals were mutated, polyadenylation became the preferred mode of OXT6 processing. These results demonstrate the interplay between alternative splicing and alternative polyadenylation, and it is their concerted actions that define a gene’s expression outcome.

The role of alternative polyadenylation in the antiviral innate immune response

Alternative polyadenylation (APA) is an important regulatory mechanism of gene functions in many biological processes. However, the extent of 3' UTR variation and the function of APA during the innate antiviral immune response are unclear. Here, we show genome-wide poly(A) sites switch and average 3' UTR length shortens gradually in response to vesicular stomatitis virus (VSV) infection in macrophages. Genes with APA and mRNA abundance change are enriched in immune-related categories such as the Toll-like receptor, RIG-I-like receptor, JAK-STAT and apoptosis-related signalling pathways. The expression of 3' processing factors is down-regulated upon VSV infection. When the core 3' processing factors are knocked down, viral replication is affected. Thus, our study reports the annotation of genes with APA in antiviral immunity and highlights the roles of 3' processing factors on 3' UTR variation upon viral infection.

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.

In Planta Determination of the mRNA-Binding Proteome of Arabidopsis Etiolated Seedlings

RNA binding proteins (RBPs) control the fate and expression of a transcriptome. Despite this fundamental importance, our understanding of plant RBPs is rudimentary, being mainly derived via bioinformatic extrapolation from other kingdoms. Here, we adapted the mRNA-protein interactome capture method to investigate the RNA binding proteome in planta. From Arabidopsis thaliana etiolated seedlings, we captured more than 700 proteins, including 300 with high confidence that we have defined as the At-RBP set. Approximately 75% of these At-RBPs are bioinformatically linked with RNA biology, containing a diversity of canonical RNA binding domains (RBDs). As no prior experimental RNA binding evidence exists for the majority of these proteins, their capture now authenticates them as RBPs. Moreover, we identified protein families harboring emerging and potentially novel RBDs, including WHIRLY, LIM, ALBA, DUF1296, and YTH domain-containing proteins, the latter being homologous to animal RNA methylation readers. Other At-RBP set proteins include major signaling proteins, cytoskeleton-associated proteins, membrane transporters, and enzymes, suggesting the scope and function of RNA-protein interactions within a plant cell is much broader than previously appreciated. Therefore, our foundation data set has provided an unbiased insight into the RNA binding proteome of plants, on which future investigations into plant RBPs can be based.

Deciphering the epitranscriptome: A green perspective

The advent of high-throughput sequencing technologies coupled with new detection methods of RNA modifications has enabled investigation of a new layer of gene regulation - the epitranscriptome. With over 100 known RNA modifications, understanding the repertoire of RNA modifications is a huge undertaking. This review summarizes what is known about RNA modifications with an emphasis on discoveries in plants. RNA ribose modifications, base methylations and pseudouridylation are required for normal development in Arabidopsis, as mutations in the enzymes modifying them have diverse effects on plant development and stress responses. These modifications can regulate RNA structure, turnover and translation. Transfer RNA and ribosomal RNA modifications have been mapped extensively and their functions investigated in many organisms, including plants. Recent work exploring the locations, functions and targeting of N⁶ -methyladenosine (m⁶ A), 5-methylcytosine (m⁵ C), pseudouridine (Ψ), and additional modifications in mRNAs and ncRNAs are highlighted, as well as those previously known on tRNAs and rRNAs. Many questions remain as to the exact mechanisms of targeting and functions of specific modified sites and whether these modifications have distinct functions in the different classes of RNAs.

Alternative polyadenylation of mRNA precursors

Alternative polyadenylation (APA) is an RNA-processing mechanism that generates distinct 3' termini on mRNAs and other RNA polymerase II transcripts. It is widespread across all eukaryotic species and is recognized as a major mechanism of gene regulation. APA exhibits tissue specificity and is important for cell proliferation and differentiation. In this Review, we discuss the roles of APA in diverse cellular processes, including mRNA metabolism, protein diversification and protein localization, and more generally in gene regulation. We also discuss the molecular mechanisms underlying APA, such as variation in the concentration of core processing factors and RNA-binding proteins, as well as transcription-based regulation.

PlantAPA: A Portal for Visualization and Analysis of Alternative Polyadenylation in Plants

Alternative polyadenylation (APA) is an important layer of gene regulation that produces mRNAs that have different 3' ends and/or encode diverse protein isoforms. Up to 70% of annotated genes in plants undergo APA. Increasing numbers of poly(A) sites collected in various plant species demand new methods and tools to access and mine these data. We have created an open-access web service called PlantAPA (http://bmi.xmu.edu.cn/plantapa) to visualize and analyze genome-wide poly(A) sites in plants. PlantAPA provides various interactive and dynamic graphics and seamlessly integrates a genome browser that can profile heterogeneous cleavage sites and quantify expression patterns of poly(A) sites across different conditions. Particularly, through PlantAPA, users can analyze poly(A) sites in extended 3' UTR regions, intergenic regions, and ambiguous regions owing to alternative transcription or RNA processing. In addition, it also provides tools for analyzing poly(A) site selections, 3' UTR lengthening or shortening, non-canonical APA site switching, and differential gene expression between conditions, making it more powerful for the study of APA-mediated gene expression regulation. More importantly, PlantAPA offers a bioinformatics pipeline that allows users to upload their own short reads or ESTs for poly(A) site extraction, enabling users to further explore poly(A) site selection using stored PlantAPA poly(A) sites together with their own poly(A) site datasets. To date, PlantAPA hosts the largest database of APA sites in plants, including Oryza sativa, Arabidopsis thaliana, Medicago truncatula, and Chlamydomonas reinhardtii. As a user-friendly web service, PlantAPA will be a valuable addition to the community of biologists studying APA mechanisms and gene expression regulation in plants.

A Genome-wide Study of "Non-3UTR" Polyadenylation Sites in Arabidopsis thaliana

Alternative polyadenylation has been recognized as a key contributor of gene expression regulation by generating different transcript isoforms with altered 3′ ends. Although polyadenylation is well known for marking the end of a 3′ UTR, an increasing number of studies have reported previously less-addressed polyadenylation events located in other parts of genes in many eukaryotic organisms. These other locations include 5′ UTRs, introns and coding sequences (termed herein as non-3UTR), as well as antisense and intergenic polyadenlation. Focusing on the non-3UTR polyadenylation sites (n3PASs), we detected and characterized more than 11000 n3PAS clusters in the Arabidopsis genome using poly(A)-tag sequencing data (PAT-Seq). Further analyses suggested that the occurrence of these n3PASs were positively correlated with certain characteristics of their respective host genes, including the presence of spliced, diminutive or diverse beginning of 5′ UTRs, number of introns and whether introns have extreme lengths. The interaction of the host genes with surrounding genetic elements, like a convergently overlapped gene and associated transposable element, may contribute to the generation of a n3PAS as well. Collectively, these results provide a better understanding of n3PASs, and offer some new insights of the underlying mechanisms for non-3UTR polyadenylation and its regulation in plants.

Dynamic Regulation of Tandem 3' Untranslated Regions in Zebrafish Spleen Cells during Immune Response

Alternative polyadenylation (APA) has been found to be involved in tumorigenesis, development, and cell differentiation, as well as in the activation of several subsets of immune cells in vitro. Whether APA takes place in immune responses in vivo is largely unknown. We profiled the variation in tandem 3' untranslated regions (UTRs) in pathogen-challenged zebrafish and identified hundreds of APA genes with ∼ 10% being immune response genes. The detected immune response APA genes were enriched in TLR signaling, apoptosis, and JAK-STAT signaling pathways. A greater number of microRNA target sites and AU-rich elements were found in the extended 3' UTRs than in the common 3' UTRs of these APA genes. Further analysis suggested that microRNA and AU-rich element-mediated posttranscriptional regulation plays an important role in modulating the expression of APA genes. These results indicate that APA is extensively involved in immune responses in vivo, and it may be a potential new paradigm for immune regulation.

Genome-wide determination of poly(A) site choice in plants

Second generation DNA sequencing technologies have been a great boon for the study of mRNA polyadenylation. The experimental determination of large numbers of polyadenylation sites using high-throughput sequencing strategies has provided the necessary platform for deeper understanding the regulation of gene expression in eukaryotes. For generating large sets of data to map poly-A sites, specialized sample preparations that target the junction of 3'-UTR and the poly(A) tail are usually employed. Here, we describe three different protocols that are effectively used for global determinations of poly(A) site choice in plants.

RADPRE: a computational program for identification of differential mRNA processing including alternative polyadenylation

Genome-wide studies revealed the prevalence of multiple transcripts resulting from alternative polyadenylation (APA) of a single given gene in higher eukaryotes. Several studies in the past few years attempted to address how those APA events are regulated and what the biological consequences of those regulations are. Common to these efforts is the comparison of unbiased transcriptome data, either derived from whole-genome tiling array or next generation sequencing, to identify the specific APA events in a given condition. RADPRE (Ratio-based Analysis of Differential mRNA Processing and Expression) is an R program, developed to serve such a purpose using data from the whole-genome tilling array. RADPRE took a set of tilling array data as input, performed a series of calculation including a correction of the probe affinity variation, a hierarchy of statistical tests and an estimation of the false discovery rate (FDR) of the differentially processed genes (DPG). The result was an output of a few tabular files including DPG and their corresponding FDR. This chapter is written for scientists with limited programming experiences.

Analysis of poly(A) site choice using a Java-based clustering algorithm

Modern high-throughput DNA sequencing has the potential to generate large volumes of data for analysis by investigators-including poly(A) site data. Here I describe a computational method to compare poly(A) site choice differences between two large data sets based on the relative abundance and position of tags within each reference sequence to which they are aligned. This method provides rapid quantification and visualization of differences and similarities in poly(A) site choice between the two datasets.

Genome-Wide Analysis of PAPS1-Dependent Polyadenylation Identifies Novel Roles for Functionally Specialized Poly(A) Polymerases in Arabidopsis thaliana

The poly(A) tail at 3’ ends of eukaryotic mRNAs promotes their nuclear export, stability and translational efficiency, and changes in its length can strongly impact gene expression. The Arabidopsis thaliana genome encodes three canonical nuclear poly(A) polymerases, PAPS1, PAPS2 and PAPS4. As shown by their different mutant phenotypes, these three isoforms are functionally specialized, with PAPS1 modifying organ growth and suppressing a constitutive immune response. However, the molecular basis of this specialization is largely unknown. Here, we have estimated poly(A)-tail lengths on a transcriptome-wide scale in wild-type and paps1 mutants. This identified categories of genes as particularly strongly affected in paps1 mutants, including genes encoding ribosomal proteins, cell-division factors and major carbohydrate-metabolic proteins. We experimentally verified two novel functions of PAPS1 in ribosome biogenesis and redox homoeostasis that were predicted based on the analysis of poly(A)-tail length changes in paps1 mutants. When overlaying the PAPS1-dependent effects observed here with coexpression analysis based on independent microarray data, the two clusters of transcripts that are most closely coexpressed with PAPS1 show the strongest change in poly(A)-tail length and transcript abundance in paps1 mutants in our analysis. This suggests that their coexpression reflects at least partly the preferential polyadenylation of these transcripts by PAPS1 versus the other two poly(A)-polymerase isoforms. Thus, transcriptome-wide analysis of poly(A)-tail lengths identifies novel biological functions and likely target transcripts for polyadenylation by PAPS1. Data integration with large-scale co-expression data suggests that changes in the relative activities of the isoforms are used as an endogenous mechanism to co-ordinately modulate plant gene expression.

The poly(A) tail of eukaryotic mRNAs promotes export from the nucleus, translation in the cytoplasm and stability of the mRNA, and changes in poly(A)-tail length can strongly impact on gene expression. The Arabidopsis thaliana genome encodes three nuclear canonical poly(A) polymerases (PAPS1, PAPS2, PAPS4) that fulfill different functions, presumably by preferentially polyadenylating certain subpopulations of pre-mRNAs. Here, we use a fractionation-based technique to assess the transcriptome-wide impact of reduced PAPS1 activity and identify functional classes of transcripts that are particularly sensitive to reduced PAPS1 activity. Analysis of these transcripts identifies two novel biological functions for PAPS1 in ribosome biogenesis and in redox homeostasis that we confirm experimentally. By overlaying our results with information about genome-wide co-expression, we demonstrate that genes co-expressed with PAPS1 are the most strongly affected in terms of poly(A)-tail length and total-abundance changes in the paps1 mutants. This provides strong evidence that the co-expression of these genes with PAPS1 that is seen across thousands of microarrays is at least partly caused by altered activity of the PAPS1 isoform, suggesting that the plant indeed uses modulation of the balance of isoform activities to coordinately regulate the expression of groups of genes.

Genome-wide characterization of intergenic polyadenylation sites redefines gene spaces in Arabidopsis thaliana

Background

Messenger RNA polyadenylation is an essential step for the maturation of most eukaryotic mRNAs. Accurate determination of poly(A) sites helps define the 3’-ends of genes, which is important for genome annotation and gene function research. Genomic studies have revealed the presence of poly(A) sites in intergenic regions, which may be attributed to 3’-UTR extensions and novel transcript units. However, there is no systematically evaluation of intergenic poly(A) sites in plants.

Results

Approximately 16,000 intergenic poly(A) site clusters (IPAC) in Arabidopsis thaliana were discovered and evaluated at the whole genome level. Based on the distributions of distance from IPACs to nearby sense and antisense genes, these IPACs were classified into three categories. About 70 % of them were from previously unannotated 3’-UTR extensions to known genes, which would extend 6985 transcripts of TAIR10 genome annotation beyond their 3’-ends, with a mean extension of 134 nucleotides. 1317 IPACs were originated from novel intergenic transcripts, 37 of which were likely to be associated with protein coding transcripts. 2957 IPACs corresponded to antisense transcripts for genes on the reverse strand, which might affect 2265 protein coding genes and 39 non-protein-coding genes, including long non-coding RNA genes. The rest of IPACs could be originated from transcriptional read-through or gene mis-annotations.

Conclusions

The identified IPACs corresponding to novel transcripts, 3’-UTR extensions, and antisense transcription should be incorporated into current Arabidopsis genome annotation. Comprehensive characterization of IPACs from this study provides insights of alternative polyadenylation and antisense transcription in plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1691-1) contains supplementary material, which is available to authorized users.

CPSF30 at the Interface of Alternative Polyadenylation and Cellular Signaling in Plants

Post-transcriptional processing, involving cleavage of precursor messenger RNA (pre mRNA), and further incorporation of poly(A) tail to the 3' end is a key step in the expression of genetic information. Alternative polyadenylation (APA) serves as an important check point for the regulation of gene expression. Recent studies have shown widespread prevalence of APA in diverse systems. A considerable amount of research has been done in characterizing different subunits of so-called Cleavage and Polyadenylation Specificity Factor (CPSF). In plants, CPSF30, an ortholog of the 30 kD subunit of mammalian CPSF is a key polyadenylation factor. CPSF30 in the model plant Arabidopsis thaliana was reported to possess unique biochemical properties. It was also demonstrated that poly(A) site choice in a vast majority of genes in Arabidopsis are CPSF30 dependent, suggesting a pivotal role of this gene in APA and subsequent regulation of gene expression. There are also indications of this gene being involved in oxidative stress and defense responses and in cellular signaling, suggesting a role of CPSF30 in connecting physiological processes and APA. This review will summarize the biochemical features of CPSF30, its role in regulating APA, and possible links with cellular signaling and stress response modules.

A 3' RACE protocol to confirm polyadenylation sites

A routine procedure in the study of polyadenylation involves the determinations of the junctions of the mRNA body and the poly(A) tail (the poly(A) site). This is typically accomplished by selectively amplifying the 3'-ends of cDNAs (3'-Rapid Amplification of cDNA Ends, or 3'-RACE), followed by sequencing of individual clones. We have developed a modification of the standard 3'-RACE protocol that couples high specificity with the ability to characterize thousands (or more) of individual sequences. This protocol may be used for numerous purposes, including the confirmation of results obtained using high throughput sequencing technologies.

A rapid, simple, and inexpensive method for the preparation of strand-specific RNA-Seq libraries

High-throughput sequencing of short cDNA tags, or RNA-Seq, has become a staple of genome-wide gene expression studies in plants. RNA-Seq libraries necessarily contain tags that correspond to the mRNA-poly(A) junction, or polyadenylation site, and thus may be mined for data that can help study alternative polyadenylation. This report presents a simple, rapid, and inexpensive method for preparing strand-specific RNA-Seq libraries from varying quantities of total RNA.

To die or not to die? Lessons from lesion mimic mutants

Programmed cell death (PCD) is a ubiquitous genetically regulated process consisting in an activation of finely controlled signaling pathways that lead to cellular suicide. Although some aspects of PCD control appear evolutionary conserved between plants, animals and fungi, the extent of conservation remains controversial. Over the last decades, identification and characterization of several lesion mimic mutants (LMM) has been a powerful tool in the quest to unravel PCD pathways in plants. Thanks to progress in molecular genetics, mutations causing the phenotype of a large number of LMM and their related suppressors were mapped, and the identification of the mutated genes shed light on major pathways in the onset of plant PCD such as (i) the involvements of chloroplasts and light energy, (ii) the roles of sphingolipids and fatty acids, (iii) a signal perception at the plasma membrane that requires efficient membrane trafficking, (iv) secondary messengers such as ion fluxes and ROS and (v) the control of gene expression as the last integrator of the signaling pathways.

Integration of developmental and environmental signals via a polyadenylation factor in Arabidopsis

The ability to integrate environmental and developmental signals with physiological responses is critical for plant survival. How this integration is done, particularly through posttranscriptional control of gene expression, is poorly understood. Previously, it was found that the 30 kD subunit of Arabidopsis cleavage and polyadenylation specificity factor (AtCPSF30) is a calmodulin-regulated RNA-binding protein. Here we demonstrated that mutant plants (oxt6) deficient in AtCPSF30 possess a novel range of phenotypes--reduced fertility, reduced lateral root formation, and altered sensitivities to oxidative stress and a number of plant hormones (auxin, cytokinin, gibberellic acid, and ACC). While the wild-type AtCPSF30 (C30G) was able to restore normal growth and responses, a mutant AtCPSF30 protein incapable of interacting with calmodulin (C30GM) could only restore wild-type fertility and responses to oxidative stress and ACC. Thus, the interaction with calmodulin is important for part of AtCPSF30 functions in the plant. Global poly(A) site analysis showed that the C30G and C30GM proteins can restore wild-type poly(A) site choice to the oxt6 mutant. Genes associated with hormone metabolism and auxin responses are also affected by the oxt6 mutation. Moreover, 19 genes that are linked with calmodulin-dependent CPSF30 functions, were identified through genome-wide expression analysis. These data, in conjunction with previous results from the analysis of the oxt6 mutant, indicate that the polyadenylation factor AtCPSF30 is a regulatory hub where different signaling cues are transduced, presumably via differential mRNA 3' end formation or alternative polyadenylation, into specified phenotypic outcomes. Our results suggest a novel function of a polyadenylation factor in environmental and developmental signal integration.

CPSF30 and Wdr33 directly bind to AAUAAA in mammalian mRNA 3' processing

AAUAAA is the most highly conserved motif in eukaryotic mRNA polyadenylation sites and, in mammals, is specifically recognized by the multisubunit CPSF (cleavage and polyadenylation specificity factor) complex. Despite its critical functions in mRNA 3' end formation, the molecular basis for CPSF-AAUAAA interaction remains poorly defined. The CPSF subunit CPSF160 has been implicated in AAUAAA recognition, but direct evidence has been lacking. Using in vitro and in vivo assays, we unexpectedly found that CPSF subunits CPSF30 and Wdr33 directly contact AAUAAA. Importantly, the CPSF30-RNA interaction is essential for mRNA 3' processing and is primarily mediated by its zinc fingers 2 and 3, which are specifically targeted by the influenza protein NS1A to suppress host mRNA 3' processing. Our data suggest that AAUAAA recognition in mammalian mRNA 3' processing is more complex than previously thought and involves multiple protein-RNA interactions.

CPSF30 and Wdr33 directly bind to AAUAAA in mammalian mRNA 3' processing

AAUAAA is the most highly conserved motif in eukaryotic mRNA polyadenylation sites and, in mammals, is specifically recognized by the multisubunit CPSF complex. Chan et al. found that CPSF subunits CPSF30 and Wdr33 directly contact AAUAAA. The CPSF30–RNA interaction is essential for mRNA 3′ processing and is primarily mediated by its zinc fingers 2 and 3, which are specifically targeted by the influenza protein NS1A to suppress host mRNA 3′ processing.

AAUAAA is the most highly conserved motif in eukaryotic mRNA polyadenylation sites and, in mammals, is specifically recognized by the multisubunit CPSF (cleavage and polyadenylation specificity factor) complex. Despite its critical functions in mRNA 3′ end formation, the molecular basis for CPSF–AAUAAA interaction remains poorly defined. The CPSF subunit CPSF160 has been implicated in AAUAAA recognition, but direct evidence has been lacking. Using in vitro and in vivo assays, we unexpectedly found that CPSF subunits CPSF30 and Wdr33 directly contact AAUAAA. Importantly, the CPSF30–RNA interaction is essential for mRNA 3′ processing and is primarily mediated by its zinc fingers 2 and 3, which are specifically targeted by the influenza protein NS1A to suppress host mRNA 3′ processing. Our data suggest that AAUAAA recognition in mammalian mRNA 3′ processing is more complex than previously thought and involves multiple protein–RNA interactions.

The Arabidopsis polyadenylation factor subunit CPSF30 as conceptual link between mRNA polyadenylation and cellular signaling

Alternative polyadenylation plays important roles in growth processes in plants. Although the scope and significance of the phenomenon have been described to considerable extent, the mechanisms that govern differential poly(A) site selection remain active areas of investigation. Of particular interest are the means by which the factors that control differential poly(A) site choice are themselves activated and inhibited. In this review, the case is made that one particular Arabidopsis polyadenylation factor subunit, termed AtCPSF30, stands out as a conceptual link between cellular signaling pathways and differential poly(A) site choice.

Genome-wide determination of poly(A) sites in Medicago truncatula: evolutionary conservation of alternative poly(A) site choice

Background

Alternative polyadenylation (APA) plays an important role in the post-transcriptional regulation of gene expression. Little is known about how APA sites may evolve in homologous genes in different plant species. To this end, comparative studies of APA sites in different organisms are needed. In this study, a collection of poly(A) sites in Medicago truncatula, a model system for legume plants, has been generated and compared with APA sites in Arabidopsis thaliana.

Results

The poly(A) tags from a deep-sequencing protocol were mapped to the annotated M. truncatula genome, and the identified poly(A) sites used to update the annotations of 14,203 genes. The results show that 64% of M. truncatula genes possess more than one poly(A) site, comparable to the percentages reported for Arabidopsis and rice. In addition, the poly(A) signals associated with M. truncatula genes were similar to those seen in Arabidopsis and other plants. The 3′-UTR lengths are correlated in pairs of orthologous genes between M. truncatula and Arabidopsis. Very little conservation of intronic poly(A) sites was found between Arabidopsis and M. truncatula, which suggests that such sites are likely to be species-specific in plants. In contrast, there is a greater conservation of CDS-localized poly(A) sites in these two species. A sizeable number of M. truncatula antisense poly(A) sites were found. A high percentage of the associated target genes possess Arabidopsis orthologs that are also associated with antisense sites. This is suggestive of important roles for antisense regulation of these target genes.

Conclusions

Our results reveal some distinct patterns of sense and antisense poly(A) sites in Arabidopsis and M. truncatula. In so doing, this study lends insight into general evolutionary trends of alternative polyadenylation in plants.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-615) contains supplementary material, which is available to authorized users.

CFIm25 links alternative polyadenylation to glioblastoma tumour suppression

The global shortening of messenger RNAs through alternative polyadenylation (APA) that occurs during enhanced cellular proliferation represents an important, yet poorly understood mechanism of regulated gene expression. The 3' untranslated region (UTR) truncation of growth-promoting mRNA transcripts that relieves intrinsic microRNA- and AU-rich-element-mediated repression has been observed to correlate with cellular transformation; however, the importance to tumorigenicity of RNA 3'-end-processing factors that potentially govern APA is unknown. Here we identify CFIm25 as a broad repressor of proximal poly(A) site usage that, when depleted, increases cell proliferation. Applying a regression model on standard RNA-sequencing data for novel APA events, we identified at least 1,450 genes with shortened 3' UTRs after CFIm25 knockdown, representing 11% of significantly expressed mRNAs in human cells. Marked increases in the expression of several known oncogenes, including cyclin D1, are observed as a consequence of CFIm25 depletion. Importantly, we identified a subset of CFIm25-regulated APA genes with shortened 3' UTRs in glioblastoma tumours that have reduced CFIm25 expression. Downregulation of CFIm25 expression in glioblastoma cells enhances their tumorigenic properties and increases tumour size, whereas CFIm25 overexpression reduces these properties and inhibits tumour growth. These findings identify a pivotal role of CFIm25 in governing APA and reveal a previously unknown connection between CFIm25 and glioblastoma tumorigenicity.

The Polyadenylation Factor Subunit CLEAVAGE AND POLYADENYLATION SPECIFICITY FACTOR30: A Key Factor of Programmed Cell Death and a Regulator of Immunity in Arabidopsis

Programmed cell death (PCD) is essential for several aspects of plant life, including development and stress responses. Indeed, incompatible plant-pathogen interactions are well known to induce the hypersensitive response, a localized cell death. Mutational analyses have identified several key PCD components, and we recently identified the mips1 mutant of Arabidopsis (Arabidopsis thaliana), which is deficient for the key enzyme catalyzing the limiting step of myoinositol synthesis. One of the most striking features of mips1 is the light-dependent formation of lesions on leaves due to salicylic acid (SA)-dependent PCD, revealing roles for myoinositol or inositol derivatives in the regulation of PCD. Here, we identified a regulator of plant PCD by screening for mutants that display transcriptomic profiles opposing that of the mips1 mutant. Our screen identified the oxt6 mutant, which has been described previously as being tolerant to oxidative stress. In the oxt6 mutant, a transfer DNA is inserted in the CLEAVAGE AND POLYADENYLATION SPECIFICITY FACTOR30 (CPSF30) gene, which encodes a polyadenylation factor subunit homolog. We show that CPSF30 is required for lesion formation in mips1 via SA-dependent signaling, that the prodeath function of CPSF30 is not mediated by changes in the glutathione status, and that CPSF30 activity is required for Pseudomonas syringae resistance. We also show that the oxt6 mutation suppresses cell death in other lesion-mimic mutants, including lesion-simulating disease1, mitogen-activated protein kinase4, constitutive expressor of pathogenesis-related genes5, and catalase2, suggesting that CPSF30 and, thus, the control of messenger RNA 3' end processing, through the regulation of SA production, is a key component of plant immune responses.

Structure and function of pre-mRNA 5'-end capping quality control and 3'-end processing

Messenger RNA precursors (pre-mRNAs) are produced as the nascent transcripts of RNA polymerase II (Pol II) in eukaryotes and must undergo extensive maturational processing, including 5′-end capping, splicing, and 3′-end cleavage and polyadenylation. This review will summarize the structural and functional information reported over the past few years on the large machinery required for the 3′-end processing of most pre-mRNAs, as well as the distinct machinery for the 3′-end processing of replication-dependent histone pre-mRNAs, which have provided great insights into the proteins and their subcomplexes in these machineries. Structural and biochemical studies have also led to the identification of a new class of enzymes (the DXO family enzymes) with activity toward intermediates of the 5′-end capping pathway. Functional studies demonstrate that these enzymes are part of a novel quality surveillance mechanism for pre-mRNA 5′-end capping. Incompletely capped pre-mRNAs are produced in yeast and human cells, in contrast to the general belief in the field that capping always proceeds to completion, and incomplete capping leads to defects in splicing and 3′-end cleavage in human cells. The DXO family enzymes are required for the detection and degradation of these defective RNAs.

Bioinformatics analysis of alternative polyadenylation in green alga Chlamydomonas reinhardtii using transcriptome sequences from three different sequencing platforms

Messenger RNA 3′-end formation is an essential posttranscriptional processing step for most eukaryotic genes. Different from plants and animals where AAUAAA and its variants routinely are found as the main poly(A) signal, Chlamydomonas reinhardtii uses UGUAA as the major poly(A) signal. The advance of sequencing technology provides an enormous amount of sequencing data for us to explore the variations of poly(A) signals, alternative polyadenylation (APA), and its relationship with splicing in this algal species. Through genome-wide analysis of poly(A) sites in C. reinhardtii, we identified a large number of poly(A) sites: 21,041 from Sanger expressed sequence tags, 88,184 from 454, and 195,266 from Illumina sequence reads. In comparison with previous collections, more new poly(A) sites are found in coding sequences and intron and intergenic regions by deep-sequencing. Interestingly, G-rich signals are particularly abundant in intron and intergenic regions. The prevalence of different poly(A) signals between coding sequences and a 3′-untranslated region implies potentially different polyadenylation mechanisms. Our data suggest that the APA occurs in about 68% of C. reinhardtii genes. Using Gene Ontolgy analysis, we found most of the APA genes are involved in RNA regulation and metabolic process, protein synthesis, hydrolase, and ligase activities. Moreover, intronic poly(A) sites are more abundant in constitutively spliced introns than retained introns, suggesting an interplay between polyadenylation and splicing. Our results support that APA, as in higher eukaryotes, may play significant roles in increasing transcriptome diversity and gene expression regulation in this algal species. Our datasets also provide useful information for accurate annotation of transcript ends in C. reinhardtii.

Genome-wide identification of members in the YTH domain-containing RNA-binding protein family in apple and expression analysis of their responsiveness to senescence and abiotic stresses

YT521-homology (YTH) domain-containing RNA-binding proteins (YTPs) are a small gene family involved in post-transcriptional regulation. We identified 26 putative YTP gene models in the apple genome. Although plant YTPs have been classified into three groups, those in multi-cellular organisms belong only to Groups A and B. The apple genome contains 22 YTP gene models in Group A and four in Group B. Duplication analysis showed that tandem and segmental duplications contributed only partially to an expansion in apple YTP numbers. YTH was the only recognizable domain in apple YTPs; its three-dimensional structure implied possible motifs for RNA-binding. After the assembly of expressed sequence tags (ESTs) and gene-cloning, we were able to identify 14 apple YTPs that were expressed in various tissues, especially senescing leaves. Expression analysis showed that these YTPs also responded to several abiotic stresses. Taken together, our genome-wide evaluation provides new insight for further research on the effects of those stresses.

The regulatory role of Pcf11-similar-4 (PCFS4) in Arabidopsis development by genome-wide physical interactions with target loci

Background

The yeast and human Pcf11 functions in both constitutive and regulated transcription and pre-mRNA processing. The constitutive roles of PCF11 are largely mediated by its direct interaction with RNA Polymerase II C-terminal domain and a polyadenylation factor, Clp1. However, little is known about the mechanism of the regulatory roles of Pcf11. Though similar to Pcf11 in multiple aspects, Arabidopsis Pcf11-similar-4 protein (PCFS4) plays only a regulatory role in Arabidopsis gene expression. Towards understanding how PCFS4 regulates the expression of its direct target genes in a genome level, ChIP-Seq approach was employed in this study to identify PCFS4 enrichment sites (ES) and the ES-linked genes within the Arabidopsis genome.

Results

A total of 892 PCFS4 ES sites linked to 839 genes were identified. Distribution analysis of the ES sites along the gene bodies suggested that PCFS4 is preferentially located on the coding sequences of the genes, consistent with its regulatory role in transcription and pre-mRNA processing. Gene ontology (GO) analysis revealed that the ES-linked genes were specifically enriched in a few GO terms, including those categories of known PCFS4 functions in Arabidopsis development. More interestingly, GO analysis suggested novel roles of PCFS4. An example is its role in circadian rhythm, which was experimentally verified herein. ES site sequences analysis identified some over-represented sequence motifs shared by subsets of ES sites. The motifs may explain the specificity of PCFS4 on its target genes and the PCFS4's functions in multiple aspects of Arabidopsis development and behavior.

Conclusions

Arabidopsis PCFS4 has been shown to specifically target on, and physically interact with, the subsets of genes. Its targeting specificity is likely mediated by cis-elements shared by the genes of each subset. The potential regulation on both transcription and mRNA processing levels of each subset of the genes may explain the functions of PCFS4 in multiple aspects of Arabidopsis development and behavior.

High throughput characterizations of poly(A) site choice in plants

The polyadenylation of mRNA in eukaryotes is an important biological process. In recent years, significant progress has been made in the field of mRNA polyadenylation owing to the advent of the next generation DNA sequencing technologies. The high-throughput sequencing capabilities have resulted in the direct experimental determinations of large numbers of polyadenylation sites, analysis of which has revealed a vast potential for the regulation of gene expression in eukaryotes. These collections have been generated using specialized sequencing methods that are targeted to the junction of 3'-UTR and the poly(A) tail. Here we present three variations of such a protocol that has been used for the analysis of alternative polyadenylation in plants. While all these methods use oligo-dT as an anchor to the 3'-end, they differ in the means of generating an anchor for the 5'-end in order to produce PCR products suitable for effective Illumina sequencing; the use of different methods to append 5' adapters expands the possible utility of these approaches. These methods are versatile, reproducible, and may be used for gene expression analysis as well as global determinations of poly(A) site choice.