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

The Arabidopsis U1 snRNP regulates mRNA 3'-end processing

The removal of introns by the spliceosome is a key gene regulatory mechanism in eukaryotes, with the U1 snRNP subunit playing a crucial role in the early stages of splicing. Studies in metazoans show that the U1 snRNP also conducts splicing-independent functions, but the lack of genetic tools and knowledge about U1 snRNP-associated proteins have limited the study of such splicing-independent functions in plants. Here we describe an RNA-centric approach that identified more than 200 proteins associated with the Arabidopsis U1 snRNP and revealed a tight link to mRNA cleavage and polyadenylation factors. Interestingly, we found that the U1 snRNP protects mRNAs against premature cleavage and polyadenylation within introns-a mechanism known as telescripting in metazoans-while also influencing alternative polyadenylation site selection in 3'-UTRs. Overall, our work provides a comprehensive view of U1 snRNP interactors and reveals novel functions in regulating mRNA 3'-end processing in Arabidopsis, laying the groundwork for understanding non-canonical functions of plant U1 snRNPs.

SIZ1-mediated SUMOylation of CPSF100 promotes plant thermomorphogenesis by controlling alternative polyadenylation

Under warm temperatures, plants adjust their morphologies for environmental adaption via precise gene expression regulation. However, the function and regulation of alternative polyadenylation (APA), an important fine-tuning of gene expression, remains unknown in plant thermomorphogenesis. In this study, we found that SUMOylation, a critical post-translational modification, is induced by a long-term treatment at warm temperatures via a SUMO ligase SIZ1 in Arabidopsis. Disruption of SIZ1 altered the global usage of polyadenylation signals and affected the APA dynamic of thermomorphogenesis-related genes. CPSF100, a key subunit of the CPSF complex for polyadenylation regulation, is SUMOylated by SIZ1. Importantly, we demonstrated that SUMOylation is essential for the function of CPSF100 in genome-wide polyadenylation site choice during thermomorphogenesis. Further analyses revealed that the SUMO conjugation on CPSF100 attenuates its interaction with two isoforms of its partner CPSF30, increasing the nuclear accumulation of CPSF100 for polyadenylation regulation. In summary, our study uncovers a regulatory mechanism of APA via SIZ1-mediated SUMOylation in plant thermomorphogenesis.

Modulation of diverse biological processes by CPSF, the master regulator of mRNA 3' ends

The cleavage and polyadenylation specificity factor (CPSF) complex plays a central role in the formation of mRNA 3' ends, being responsible for the recognition of the poly(A) signal sequence, the endonucleolytic cleavage step, and recruitment of poly(A) polymerase. CPSF has been extensively studied for over three decades, and its functions and those of its individual subunits are becoming increasingly well-defined, with much current research focusing on the impact of these proteins on the normal functioning or disease/stress states of cells. In this review, we provide an overview of the general functions of CPSF and its subunits, followed by a discussion of how they exert their functions in a surprisingly diverse variety of biological processes and cellular conditions. These include transcription termination, small RNA processing, and R-loop prevention/resolution, as well as more generally cancer, differentiation/development, and infection/immunity.

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.

Molecular details of the CPSF73-CPSF100 C-terminal heterodimer and interaction with Symplekin

Eukaryotic pre-mRNA is processed by a large multiprotein complex to accurately cleave the 3' end, and to catalyse the addition of the poly(A) tail. Within this cleavage and polyadenylation specificity factor (CPSF) machinery, the CPSF73/CPSF3 endonuclease subunit directly contacts both CPSF100/CPSF2 and the scaffold protein Symplekin to form a subcomplex known as the core cleavage complex or mammalian cleavage factor. Here we have taken advantage of a stable CPSF73-CPSF100 minimal heterodimer from Encephalitozoon cuniculi to determine the solution structure formed by the first and second C-terminal domain (CTD1 and CTD2) of both proteins. We find a large number of contacts between both proteins in the complex, and notably in the region between CTD1 and CTD2. A similarity is also observed between CTD2 and the TATA-box binding protein (TBP) domains. Separately, we have determined the structure of the terminal CTD3 domain of CPSF73, which also belongs to the TBP domain family and is connected by a flexible linker to the rest of CPSF73. Biochemical assays demonstrate a key role for the CTD3 of CPSF73 in binding Symplekin, and structural models of the trimeric complex from other species allow for comparative analysis and support an overall conserved architecture.

DNA-dependent RNA polymerases in plants

DNA-dependent RNA polymerases (Pols) transfer the genetic information stored in genomic DNA to RNA in all organisms. In eukaryotes, the typical products of nuclear Pol I, Pol II, and Pol III are ribosomal RNAs, mRNAs, and transfer RNAs, respectively. Intriguingly, plants possess two additional Pols, Pol IV and Pol V, which produce small RNAs and long noncoding RNAs, respectively, mainly for silencing transposable elements. The five plant Pols share some subunits, but their distinct functions stem from unique subunits that interact with specific regulatory factors in their transcription cycles. Here, we summarize recent advances in our understanding of plant nucleus-localized Pols, including their evolution, function, structures, and transcription cycles.

Inhibitor AN3661 reveals biological functions of Arabidopsis CLEAVAGE and POLYADENYLATION SPECIFICITY FACTOR 73

Cleavage and polyadenylation specificity factor (CPSF) is a protein complex that plays an essential biochemical role in mRNA 3'-end formation, including poly(A) signal recognition and cleavage at the poly(A) site. However, its biological functions at the organismal level are mostly unknown in multicellular eukaryotes. The study of plant CPSF73 has been hampered by the lethality of Arabidopsis (Arabidopsis thaliana) homozygous mutants of AtCPSF73-I and AtCPSF73-II. Here, we used poly(A) tag sequencing to investigate the roles of AtCPSF73-I and AtCPSF73-II in Arabidopsis treated with AN3661, an antimalarial drug with specificity for parasite CPSF73 that is homologous to plant CPSF73. Direct seed germination on an AN3661-containing medium was lethal; however, 7-d-old seedlings treated with AN3661 survived. AN3661 targeted AtCPSF73-I and AtCPSF73-II, inhibiting growth through coordinating gene expression and poly(A) site choice. Functional enrichment analysis revealed that the accumulation of ethylene and auxin jointly inhibited primary root growth. AN3661 affected poly(A) signal recognition, resulted in lower U-rich signal usage, caused transcriptional readthrough, and increased the distal poly(A) site usage. Many microRNA targets were found in the 3' untranslated region lengthened transcripts; these miRNAs may indirectly regulate the expression of these targets. Overall, this work demonstrates that AtCPSF73 plays important part in co-transcriptional regulation, affecting growth, and development in Arabidopsis.

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.

Coupling epigenetics and RNA polyadenylation: missing links

Precise regulation of gene expression is crucial for plant survival. As a cotranscriptional regulatory mechanism, pre-mRNA polyadenylation is essential for fine-tuning gene expression. Polyadenylation can be alternatively projected at various sites of a transcript, which contributes to transcriptome diversity. Epigenetic modification is another mechanism of transcriptional control. Recent studies have uncovered crosstalk between cotranscriptional polyadenylation processes and both epigenomic and epitranscriptomic markers. Genetic analyses have demonstrated that DNA methylation, histone modifications, and epitranscriptomic modification are involved in regulating polyadenylation in plants. Here we summarize current understanding of the links between epigenetics and polyadenylation and their novel biological efficacy for plant development and environmental responses. Unresolved issues and future directions are discussed to shed light on the field.

Review: Mechanisms underlying alternative polyadenylation in plants - looking in the right places

Recent years have seen an explosion of interest in the subject of alternative polyadenylation in plants. Connections between the polyadenylation complex and numerous developmental and stress responses are well-established. However, those that link stimuli with the functioning of the polyadenylation complex are less well understood. To this end, it is imperative to clearly delineate the roles of the polyadenylation complex in both plant growth AND alternative polyadenylation. It is also necessary to understand the ways by which other molecular processes may contribute to alternative polyadenylation. This review discusses these issues, with a focus on instances that reveal mechanisms by which mRNA polyadenylation may be regulated. Insights from from characterizations of mutants affected in the polyadenylation complex are discussed, as are the limitations of such characterizations when it comes to teasing out cause and effect. These limitations encourage explorations to other processes that are beyond the core polyadenylation complex. Two such processes that sculpt the plant transcriptome - transcription termination and the epigenetic control of transposon activity - also contribute to regulated poly(A) site choice. These subjects define "the right places" - molecular mechanisms that contribute to the wide-ranging control of gene expression via mRNA polyadenylation.

Genome-wide characterization of nascent RNA processing in plants

Following transcription initiation, RNA polymerase II (Pol II) elongates through the genic region and terminates after the polyadenylation signal. This process is accompanied by splicing, 3' cleavage, and polyadenylation, to eventually form a mature mRNA. Recent advances in short-read and long-read high-throughput sequencing methods have shed light on the global landscape of these co-transcriptional events at nucleotide resolution. In this mini review, we summarize recent developments in genome-wide approaches that broadened our understanding of nascent RNA processing in plants.

Signatures of mRNA Alternative Polyadenylation in Arabidopsis Leaf Development

Alternative polyadenylation (APA) of pre-mRNA is an important co-transcriptional mechanism that modulates gene expression, leading to transcriptomic and functional diversities. The role of APA in Arabidopsis leaf development, however, remains elusive. We applied a poly(A)-tag sequencing (PAT-seq) technique to characterize APA-mediated regulation events in cotyledon and in five stages of true leaf development. Over 60% APA was identified in genes expressed in leaves, consistent with the results in previous publications. However, a reduced APA level was detected in younger leaves, reaching 44% in the 18th true leaf. Importantly, we also found that >70% of the poly(A) site usages were altered in the second true leaf relative to the cotyledon. Compared with the cotyledon, more genes in the second true leaf tended to use the distal site of 3′UTR, but this was not found in pairwise comparison among other true leaves. In addition, a significant APA gene was found to be decreased in a pairwise comparison among true leaves, including differentially expressed genes. The APA genes identified herein were associated with specific biological processes, including metabolic and cellular processes and response to stimuli and hormones. These results provide a new insight into the regulation of Arabidopsis leaf development through APA.

Full-length transcriptome sequencing reveals the molecular mechanism of potato seedlings responding to low-temperature

BACKGROUND

Potato (Solanum tuberosum L.) is one of the world's most important crops, the cultivated potato is frost-sensitive, and low-temperature severely influences potato production. However, the mechanism by which potato responds to low-temperature stress is unclear. In this research, we apply a combination of second-generation sequencing and third-generation sequencing technologies to sequence full-length transcriptomes in low-temperature-sensitive cultivars to identify the important genes and main pathways related to low-temperature resistance.

RESULTS

In this study, we obtained 41,016 high-quality transcripts, which included 15,189 putative new transcripts. Amongst them, we identified 11,665 open reading frames, 6085 simple sequence repeats out of the potato dataset. We used public available genomic contigs to analyze the gene features, simple sequence repeat, and alternative splicing event of 24,658 non-redundant transcript sequences, predicted the coding sequence and identified the alternative polyadenylation. We performed cluster analysis, GO, and KEGG functional analysis of 4518 genes that were differentially expressed between the different low-temperature treatments. We examined 36 transcription factor families and identified 542 transcription factors in the differentially expressed genes, and 64 transcription factors were found in the AP2 transcription factor family which was the most. We measured the malondialdehyde, soluble sugar, and proline contents and the expression genes changed associated with low temperature resistance in the low-temperature treated leaves. We also tentatively speculate that StLPIN10369.5 and StCDPK16 may play a central coordinating role in the response of potatoes to low temperature stress.

CONCLUSIONS

Overall, this study provided the first large-scale full-length transcriptome sequencing of potato and will facilitate structure-function genetic and comparative genomics studies of this important crop.

Divergence in the Regulation of the Salt Tolerant Response Between Arabidopsis thaliana and Its Halophytic Relative Eutrema salsugineum by mRNA Alternative Polyadenylation

Salt tolerance is an important mechanism by which plants can adapt to a saline environment. To understand the process of salt tolerance, we performed global analyses of mRNA alternative polyadenylation (APA), an important regulatory mechanism during eukaryotic gene expression, in Arabidopsis thaliana and its halophytic relative Eutrema salsugineum with regard to their responses to salt stress. Analyses showed that while APA occurs commonly in both Arabidopsis and Eutrema, Eutrema possesses fewer APA genes than Arabidopsis (47% vs. 54%). However, the proportion of APA genes was significantly increased in Arabidopsis under salt stress but not in Eutrema. This indicated that Arabidopsis is more sensitive to salt stress and that Eutrema exhibits an innate response to such conditions. Both species utilized distal poly(A) sites under salt stress; however, only eight genes were found to overlap when their 3' untranslated region (UTR) lengthen genes were compared, thus revealing their distinct responses to salt stress. In Arabidopsis, genes that use distal poly(A) sites were enriched in response to salt stress. However, in Eutrema, the use of poly(A) sites was less affected and fewer genes were enriched. The transcripts with upregulated poly(A) sites in Arabidopsis showed enriched pathways in plant hormone signal transduction, starch and sucrose metabolism, and fatty acid elongation; in Eutrema, biosynthetic pathways (stilbenoid, diarylheptanoid, and gingerol) and metabolic pathways (arginine and proline) showed enrichment. APA was associated with 42% and 29% of the differentially expressed genes (DE genes) in Arabidopsis and Eutrema experiencing salt stress, respectively. Salt specific poly(A) sites and salt-inducible APA events were identified in both species; notably, some salt tolerance-related genes and transcription factor genes exhibited differential APA patterns, such as CIPK21 and LEA4-5. Our results suggest that adapted species exhibit more orderly response at the RNA maturation step under salt stress, while more salt-specific poly(A) sites were activated in Arabidopsis to cope with salinity conditions. Collectively, our findings not only highlight the importance of APA in the regulation of gene expression in response to salt stress, but also provide a new perspective on how salt-sensitive and salt-tolerant species perform differently under stress conditions through transcriptome diversity.

Alternative polyadenylated mRNAs behave as asynchronous rhythmic transcription in Arabidopsis

Alternative polyadenylation (APA) is a widespread post-transcriptional modification method that changes the 3' ends of transcripts by altering poly(A) site usage. However, the longitudinal transcriptomic 3' end profile and its mechanism of action are poorly understood. We applied diurnal time-course poly(A) tag sequencing (PAT-seq) for Arabidopsis and identified 3284 genes that generated both rhythmic and arrhythmic transcripts. These two classes of transcripts appear to exhibit dramatic differences in expression and translation activisty. The asynchronized transcripts derived by APA are embedded with different poly(A) signals, especially for rhythmic transcripts, which contain higher AAUAAA and UGUA signal proportions. The Pol II occupancy maximum is reached upstream of rhythmic poly(A) sites, while it is present directly at arrhythmic poly(A) sites. Integrating H3K9ac and H3K4me3 time-course data analyses revealed that transcriptional activation of histone markers may be involved in the differentiation of rhythmic and arrhythmic APA transcripts. These results implicate an interplay between histone modification and RNA 3'-end processing, shedding light on the mechanism of transcription rhythm and alternative polyadenylation.

Global gene expression signatures in response to citrate-coated silver nanoparticles exposure

Silver nanoparticles (AgNPs) are widely used in medical and commercial products for their unique antibacterial functions. However, the impact of AgNPs on human neural development is not well understood. To investigate the effect of AgNPs on human neural development, various doses of 20 nm citrate-coated AgNP (AgSC) were administered to human embryonic stem cell derived neural progenitors during the neuronal differentiation. Immunofluorescence staining with neuronal progenitor markers SOX2 (sex determining region Y-box 2) and Nestin (VI intermediate filament protein) showed that AgSC inhibited rosette formation, neuronal progenitor proliferation, and neurite outgrowth. Furthermore, AgSC promoted astrocyte activation and neuronal apoptosis. These adverse effects can be partially recovered with ascorbic acid. A genome-wide transcriptome analysis of both AgSC treated and untreated samples indicated that the most up-graduated genes were a group of Metallothionein (1F, 1E, 2A) proteins, a metal-binding protein that plays an essential role in metal homeostasis, heavy metal detoxification, and cellular anti-oxidative defence. The most significantly down-regulated genes were neuronal differentiation 6 (NEUROD6) and fork head box G1 (FOXG1). GO analyse indicated that the regulation of cholesterol biosynthetic process, neuron differentiation, synapse organization and pattern specification, oliogenesis, and neuronal apoptosis were the most impacted biological processes. KEGG pathway analyse showed that the most significantly impacted pathways were C5 isoprenoid, axon guidance, Notch, WNT, RAS-MAPK signalling pathways, lysosome, and apoptosis. Our data suggests that AgSCs interfered with metal homeostasis and cholesterol biosynthesis which induced oxidative stress, inhibited neurogenesis, axon guidance, and promoted apoptosis. Supplementation with ascorbic acid could act as an antioxidant to prevent AgSC-mediated neurotoxicity.

Arabidopsis N6-methyladenosine reader CPSF30-L recognizes FUE signals to control polyadenylation site choice in liquid-like nuclear bodies

The biological functions of the epitranscriptomic modification N⁶-methyladenosine (m⁶A) in plants are not fully understood. CPSF30-L is a predominant isoform of the polyadenylation factor CPSF30 and consists of CPSF30-S and an m⁶A-binding YTH domain. Little is known about the biological roles of CPSF30-L and the molecular mechanism underlying its m⁶A-binding function in alternative polyadenylation. Here, we characterized CPSF30-L as an Arabidopsis m⁶A reader whose m⁶A-binding function is required for the floral transition and abscisic acid (ABA) response. We found that the m⁶A-binding activity of CPSF30-L enhances the formation of liquid-like nuclear bodies, where CPSF30-L mainly recognizes m⁶A-modified far-upstream elements to control polyadenylation site choice. Deficiency of CPSF30-L lengthens the 3' untranslated region of three phenotypes-related transcripts, thereby accelerating their mRNA degradation and leading to late flowering and ABA hypersensitivity. Collectively, this study uncovers a new molecular mechanism for m⁶A-driven phase separation and polyadenylation in plants.

Alternative polyadenylation: methods, mechanism, function, and role in cancer

Occurring in over 60% of human genes, alternative polyadenylation (APA) results in numerous transcripts with differing 3'ends, thus greatly expanding the diversity of mRNAs and of proteins derived from a single gene. As a key molecular mechanism, APA is involved in various gene regulation steps including mRNA maturation, mRNA stability, cellular RNA decay, and protein diversification. APA is frequently dysregulated in cancers leading to changes in oncogenes and tumor suppressor gene expressions. Recent studies have revealed various APA regulatory mechanisms that promote the development and progression of a number of human diseases, including cancer. Here, we provide an overview of four types of APA and their impacts on gene regulation. We focus particularly on the interaction of APA with microRNAs, RNA binding proteins and other related factors, the core pre-mRNA 3'end processing complex, and 3'UTR length change. We also describe next-generation sequencing methods and computational tools for use in poly(A) signal detection and APA repositories and databases. Finally, we summarize the current understanding of APA in cancer and provide our vision for future APA related research.

Profiling Alternative 3' Untranslated Regions in Sorghum using RNA-seq Data

Sorghum is an important crop widely used for food, feed, and fuel. Transcriptome-wide studies of 3′ untranslated regions (3′UTR) using regular RNA-seq remain scarce in sorghum, while transcriptomes have been characterized extensively using Illumina short-read sequencing platforms for many sorghum varieties under various conditions or developmental contexts. 3′UTR is a critical regulatory component of genes, controlling the translation, transport, and stability of messenger RNAs. In the present study, we profiled the alternative 3′UTRs at the transcriptome level in three genetically related but phenotypically contrasting lines of sorghum: Rio, BTx406, and R9188. A total of 1,197 transcripts with alternative 3′UTRs were detected using RNA-seq data. Their categorization identified 612 high-confidence alternative 3′UTRs. Importantly, the high-confidence alternative 3′UTR genes significantly overlapped with the genesets that are associated with RNA N⁶-methyladenosine (m⁶A) modification, suggesting a clear indication between alternative 3′UTR and m⁶A methylation in sorghum. Moreover, taking advantage of sorghum genetics, we provided evidence of genotype specificity of alternative 3′UTR usage. In summary, our work exemplifies a transcriptome-wide profiling of alternative 3′UTRs using regular RNA-seq data in non-model crops and gains insights into alternative 3′UTRs and their genotype specificity.

Transcript isoform sequencing reveals widespread promoter-proximal transcriptional termination in Arabidopsis

RNA polymerase II (RNAPII) transcription converts the DNA sequence of a single gene into multiple transcript isoforms that may carry alternative functions. Gene isoforms result from variable transcription start sites (TSSs) at the beginning and polyadenylation sites (PASs) at the end of transcripts. How alternative TSSs relate to variable PASs is poorly understood. Here, we identify both ends of RNA molecules in Arabidopsis thaliana by transcription isoform sequencing (TIF-seq) and report four transcript isoforms per expressed gene. While intragenic initiation represents a large source of regulated isoform diversity, we observe that ~14% of expressed genes generate relatively unstable short promoter-proximal RNAs (sppRNAs) from nascent transcript cleavage and polyadenylation shortly after initiation. The location of sppRNAs correlates with the position of promoter-proximal RNAPII stalling, indicating that large pools of promoter-stalled RNAPII may engage in transcriptional termination. We propose that promoter-proximal RNAPII stalling-linked to premature transcriptional termination may represent a checkpoint that governs plant gene expression.

The circadian clock shapes the Arabidopsis transcriptome by regulating alternative splicing and alternative polyadenylation

The circadian clock in plants temporally coordinates biological processes throughout the day, synchronizing gene expression with diurnal environmental changes. Circadian oscillator proteins are known to regulate the expression of clock-controlled plant genes by controlling their transcription. Here, using a high-throughput RNA-Seq approach, we examined genome-wide circadian and diurnal control of the Arabidopsis transcriptome, finding that the oscillation patterns of different transcripts of multitranscript genes can exhibit substantial differences and demonstrating that the circadian clock affects posttranscriptional regulation. In parallel, we found that two major posttranscriptional mechanisms, alternative splicing (AS; especially intron retention) and alternative polyadenylation (APA), display circadian rhythmicity resulting from oscillation in the genes involved in AS and APA. Moreover, AS-related genes exhibited rhythmic AS and APA regulation, adding another layer of complexity to circadian regulation of gene expression. We conclude that the Arabidopsis circadian clock not only controls transcription of genes but also affects their posttranscriptional regulation by influencing alternative splicing and alternative polyadenylation.

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

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

Differential alternative polyadenylation contributes to the developmental divergence between two rice subspecies, japonica and indica

Alternative polyadenylation (APA) is a widespread post-transcriptional mechanism that regulates gene expression through mRNA metabolism, playing a pivotal role in modulating phenotypic traits in rice (Oryza sativa L.). However, little is known about the APA-mediated regulation underlying the distinct characteristics between two major rice subspecies, indica and japonica. Using a poly(A)-tag sequencing approach, polyadenylation (poly(A)) site profiles were investigated and compared pairwise from germination to the mature stage between indica and japonica, and extensive differentiation in APA profiles was detected genome-wide. Genes with subspecies-specific poly(A) sites were found to contribute to subspecies characteristics, particularly in disease resistance of indica and cold-stress tolerance of japonica. In most tissues, differential usage of APA sites exhibited an apparent impact on the gene expression profiles between subspecies, and genes with those APA sites were significantly enriched in quantitative trait loci (QTL) related to yield traits, such as spikelet number and 1000-seed weight. In leaves of the booting stage, APA site-switching genes displayed global shortening of 3' untranslated regions with increased expression in indica compared with japonica, and they were overrepresented in the porphyrin and chlorophyll metabolism pathways. This phenomenon may lead to a higher chlorophyll content and photosynthesis in indica than in japonica, being associated with their differential growth rates and yield potentials. We further constructed an online resource for querying and visualizing the poly(A) atlas in these two rice subspecies. Our results suggest that APA may be largely involved in developmental differentiations between two rice subspecies, especially in leaf characteristics and the stress response, broadening our knowledge of the post-transcriptional genetic basis underlying the divergence of rice traits.

RNA polymerase II activity revealed by GRO-seq and pNET-seq in Arabidopsis

RNA polymerase II (Pol II) plays an essential role in gene expression. We used plant native elongating transcript sequencing and global run-on sequencing to profile nascent RNAs genome wide in Arabidopsis. We found that Pol II tends to accumulate downstream of the transcription start site (TSS). Moreover, Pol II with an unphosphorylated carboxyl-terminal domain (CTD) mainly accumulates downstream of the TSS, while Pol II with a Ser 5P CTD associates with spliceosomes, and Pol II with a Ser 2P CTD presents a sharp peak within 250 base pairs downstream of the polyadenylation site (PAS). Pol II pausing both at promoter-proximal regions and after PAS affects the transcription rate. Interestingly, active genes can be classified into three clusters based on the different modes of transcription. We demonstrate that these two methods are suitable to study Pol II dynamics in planta. Although transcription is conserved overall within eukaryotes, there is plant-specific regulation.

Genome-wide atlas of alternative polyadenylation in the forage legume red clover

Studies on prevalence and significance of alternative polyadenylation (APA) in plants have been so far limited mostly to the model plants. Here, a genome-wide analysis of APA was carried out in different tissue types in the non-model forage legume red clover (Trifolium pratense L). A profile of poly(A) sites in different tissue types was generated using so-called 'poly(A)-tag sequencing' (PATseq) approach. Our analysis revealed tissue-wise dynamics of usage of poly(A) sites located at different genomic locations. We also identified poly(A) sites and underlying genes displaying APA in different tissues. Functional categories enriched in groups of genes manifesting APA between tissue types were determined. Analysis of spatial expression of genes encoding different poly(A) factors showed significant differential expression of genes encoding orthologs of FIP1(V) and PCFS4, suggesting that these two factors may play a role in regulating spatial APA in red clover. Our analysis also revealed a high degree of conservation in diverse plant species of APA events in mRNAs encoding two key polyadenylation factors, CPSF30 and FIP1(V). Together with our previously reported study of spatial gene expression in red clover, this study will provide a comprehensive account of transcriptome dynamics in this non-model forage legume.