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

Stress responses of plants through transcriptome plasticity by mRNA alternative polyadenylation

The sessile nature of plants confines their responsiveness to changing environmental conditions. Gene expression regulation becomes a paramount mechanism for plants to adjust their physiological and morphological behaviors. Alternative polyadenylation (APA) is known for its capacity to augment transcriptome diversity and plasticity, thereby furnishing an additional set of tools for modulating gene expression. APA has also been demonstrated to exhibit intimate associations with plant stress responses. In this study, we review APA dynamic features and consequences in plants subjected to both biotic and abiotic stresses. These stresses include adverse environmental stresses, and pathogenic attacks, such as cadmium toxicity, high salt, hypoxia, oxidative stress, cold, heat shock, along with bacterial, fungal, and viral infections. We analyzed the overarching research framework employed to elucidate plant APA response and the alignment of polyadenylation site transitions with the modulation of gene expression levels within the ambit of each stress condition. We also proposed a general APA model where transacting factors, including poly(A) factors, epigenetic regulators, RNA m6A modification factors, and phase separation proteins, assume pivotal roles in APA related transcriptome plasticity during stress response in plants.

stAPAminer: Mining Spatial Patterns of Alternative Polyadenylation for Spatially Resolved Transcriptomic Studies

Alternative polyadenylation (APA) contributes to transcriptome complexity and gene expression regulation and has been implicated in various cellular processes and diseases. Single-cell RNA sequencing (scRNA-seq) has enabled the profiling of APA at the single-cell level; however, the spatial information of cells is not preserved in scRNA-seq. Alternatively, spatial transcriptomics (ST) technologies provide opportunities to decipher the spatial context of the transcriptomic landscape. Pioneering studies have revealed potential spatially variable genes and/or splice isoforms; however, the pattern of APA usage in spatial contexts remains unappreciated. In this study, we developed a toolkit called stAPAminer for mining spatial patterns of APA from spatially barcoded ST data. APA sites were identified and quantified from the ST data. In particular, an imputation model based on the k-nearest neighbors algorithm was designed to recover APA signals, and then APA genes with spatial patterns of APA usage variation were identified. By analyzing well-established ST data of the mouse olfactory bulb (MOB), we presented a detailed view of spatial APA usage across morphological layers of the MOB. We compiled a comprehensive list of genes with spatial APA dynamics and obtained several major spatial expression patterns that represent spatial APA dynamics in different morphological layers. By extending this analysis to two additional replicates of the MOB ST data, we observed that the spatial APA patterns of several genes were reproducible among replicates. stAPAminer employs the power of ST to explore the transcriptional atlas of spatial APA patterns with spatial resolution. This toolkit is available at https://github.com/BMILAB/stAPAminer and https://ngdc.cncb.ac.cn/biocode/tools/BT007320.

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.

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

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

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.

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

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

Plant 3' Regulatory Regions From mRNA-Encoding Genes and Their Uses to Modulate Expression

Molecular biotechnology has made it possible to explore the potential of plants for different purposes. The 3' regulatory regions have a great diversity of cis-regulatory elements directly involved in polyadenylation, stability, transport and mRNA translation, essential to achieve the desired levels of gene expression. A complex interaction between the cleavage and polyadenylation molecular complex and cis-elements determine the polyadenylation site, which may result in the choice of non-canonical sites, resulting in alternative polyadenylation events, involved in the regulation of more than 80% of the genes expressed in plants. In addition, after transcription, a wide array of RNA-binding proteins interacts with cis-acting elements located mainly in the 3' untranslated region, determining the fate of mRNAs in eukaryotic cells. Although a small number of 3' regulatory regions have been identified and validated so far, many studies have shown that plant 3' regulatory regions have a higher potential to regulate gene expression in plants compared to widely used 3' regulatory regions, such as NOS and OCS from Agrobacterium tumefaciens and 35S from cauliflower mosaic virus. In this review, we discuss the role of 3' regulatory regions in gene expression, and the superior potential that plant 3' regulatory regions have compared to NOS, OCS and 35S 3' regulatory regions.

PlantAPAdb: A Comprehensive Database for Alternative Polyadenylation Sites in Plants

Alternative cleavage and polyadenylation (APA) is increasingly recognized as an important regulatory mechanism in eukaryotic gene expression and is dynamically modulated in a developmental, tissue-specific, or environmentally responsive manner. Given the functional importance of APA and the rapid accumulation of APA sites in plants, a comprehensive and easily accessible APA site database is necessary for improved understanding of APA-mediated gene expression regulation. We present a database called PlantAPAdb that catalogs the most comprehensive APA site data derived from sequences from diverse 3' sequencing protocols and biological samples in plants. Currently, PlantAPAdb contains APA sites in six species, Oryza sativa (japonica and indica), Arabidopsis (Arabidopsis thaliana), Medicago truncatula, Trifolium pratense, Phyllostachys edulis, and Chlamydomonas reinhardtii APA sites in PlantAPAdb are available for bulk download and can be queried in a Google-like manner. PlantAPAdb provides rich information of the whole-genome APA sites, including genomic locations, heterogeneous cleavage sites, expression levels, and sample information. It also provides comprehensive poly(A) signals for APA sites in different genomic regions according to distinct profiles of cis-elements in plants. In addition, PlantAPAdb contains events of 3' untranslated region shortening/lengthening resulting from APA, which helps to understand the mechanisms underlying systematic changes in 3' untranslated region lengths. Additional information about conservation of APA sites in plants is also available, providing insights into the evolutionary polyadenylation configuration across species. As a user-friendly database, PlantAPAdb is a large and extendable resource for elucidating APA mechanisms, APA conservation, and gene expression regulation.

The polyadenylation factor FIP1 is important for plant development and root responses to abiotic stresses

Root development and its response to environmental changes is crucial for whole plant adaptation. These responses include changes in transcript levels. Here, we show that the alternative polyadenylation (APA) of mRNA is important for root development and responses. Mutations in FIP1, a component of polyadenylation machinery, affects plant development, cell division and elongation, and response to different abiotic stresses. Salt treatment increases the amount of poly(A) site usage within the coding region and 5' untranslated regions (5'-UTRs), and the lack of FIP1 activity reduces the poly(A) site usage within these non-canonical sites. Gene ontology analyses of transcripts displaying APA in response to salt show an enrichment in ABA signaling, and in the response to stresses such as salt or cadmium (Cd), among others. Root growth assays show that fip1-2 is more tolerant to salt but is hypersensitive to ABA or Cd. Our data indicate that FIP1-mediated alternative polyadenylation is important for plant development and stress responses.

A GmSIN1/GmNCED3s/GmRbohBs Feed-Forward Loop Acts as a Signal Amplifier That Regulates Root Growth in Soybean Exposed to Salt Stress

Abscisic acid (ABA) and reactive oxygen species (ROS) act as key signaling molecules in the plant response to salt stress; however, how these signals are transduced and amplified remains unclear. Here, a soybean (Glycine max) salinity-induced NAM/ATAF1/2/CUC2 (NAC) transcription factor encoded by SALT INDUCED NAC1 (GmSIN1) was shown to be a key component of this process. Overexpression of GmSIN1 in soybean promoted root growth and salt tolerance and increased yield under salt stress; RNA interference-mediated knockdown of GmSIN1 had the opposite effect. The rapid induction of GmSIN1 in response to salinity required ABA and ROS, and the effect of GmSIN1 on root elongation and salt tolerance was achieved by boosting cellular ABA and ROS contents. GmSIN1 upregulated 9-cis-epoxycarotenoid dioxygenase coding genes in soybean (GmNCED3s, associated with ABA synthesis) and Respiratory burst oxidase homolog B genes in soybean (GmRbohBs, associated with ROS generation) by binding to their promoters at a site that has not been described to date. Together, GmSIN1, GmNCED3s, and GmRbohBs constitute a positive feed-forward system that enables the rapid accumulation of ABA and ROS, effectively amplifying the initial salt stress signal. These findings suggest that the combined modulation of ABA and ROS contents enhances soybean salt tolerance.

Genome-Wide Profiling of Polyadenylation Events in Maize Using High-Throughput Transcriptomic Sequences

Polyadenylation is an essential post-transcriptional modification of eukaryotic transcripts that plays critical role in transcript stability, localization, transport, and translational efficiency. About 70% genes in plants contain alternative polyadenylation (APA) sites. Despite availability of vast amount of sequencing data, to date, a comprehensive map of the polyadenylation events in maize is not available. Here, 9.48 billion RNA-Seq reads were analyzed to characterize 95,345 Poly(A) Clusters (PAC) in 23,705 (51%) maize genes. Of these, 76% were APA genes. However, most APA genes (55%) expressed a dominant PAC rather than favoring multiple PACs equally. The lincRNA genes with PACs were significantly longer in length than the genes without any PAC and about 48% genes had APA sites. Heterogeneity was observed in 52% of the PACs supporting the imprecise nature of the polyadenylation process. Genomic distribution revealed that the majority of the PACs (78%) were located in the genic regions. Unlike previous studies, large number of PACs were observed in the intergenic (n = 21,264), 5′-UTR (735), CDS (2,542), and the intronic regions (12,841). The CDS and introns with PACs were longer in length than without PACs, whereas intergenic PACs were more often associated with transcripts that lacked annotated 3′-UTRs. Nucleotide composition around PACs demonstrated AT-richness and the common upstream motif was AAUAAA, which is consistent with other plants. According to this study, only 2,830 genes still maintained the use of AAUAAA motif. This large-scale data provides useful insights about the gene expression regulation and could be utilized as evidence to validate the annotation of transcript ends.

APAtrap: identification and quantification of alternative polyadenylation sites from RNA-seq data

Motivation

Alternative polyadenylation (APA) has been increasingly recognized as a crucial mechanism that contributes to transcriptome diversity and gene expression regulation. As RNA-seq has become a routine protocol for transcriptome analysis, it is of great interest to leverage such unprecedented collection of RNA-seq data by new computational methods to extract and quantify APA dynamics in these transcriptomes. However, research progress in this area has been relatively limited. Conventional methods rely on either transcript assembly to determine transcript 3' ends or annotated poly(A) sites. Moreover, they can neither identify more than two poly(A) sites in a gene nor detect dynamic APA site usage considering more than two poly(A) sites.

Results

We developed an approach called APAtrap based on the mean squared error model to identify and quantify APA sites from RNA-seq data. APAtrap is capable of identifying novel 3' UTRs and 3' UTR extensions, which contributes to locating potential poly(A) sites in previously overlooked regions and improving genome annotations. APAtrap also aims to tally all potential poly(A) sites and detect genes with differential APA site usages between conditions. Extensive comparisons of APAtrap with two other latest methods, ChangePoint and DaPars, using various RNA-seq datasets from simulation studies, human and Arabidopsis demonstrate the efficacy and flexibility of APAtrap for any organisms with an annotated genome.

Availability and implementation

Freely available for download at https://apatrap.sourceforge.io.

Contact

liqq@xmu.edu.cn or xhuister@xmu.edu.cn.

Supplementary information

Supplementary data are available at Bioinformatics online.

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.

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.

Capturing the Alternative Cleavage and Polyadenylation Sites of 14 NAC Genes in Populus Using a Combination of 3'-RACE and High-Throughput Sequencing

Detection of complex splice sites (SSs) and polyadenylation sites (PASs) of eukaryotic genes is essential for the elucidation of gene regulatory mechanisms. Transcriptome-wide studies using high-throughput sequencing (HTS) have revealed prevalent alternative splicing (AS) and alternative polyadenylation (APA) in plants. However, small-scale and high-depth HTS aimed at detecting genes or gene families are very few and limited. We explored a convenient and flexible method for profiling SSs and PASs, which combines rapid amplification of 3′-cDNA ends (3′-RACE) and HTS. Fourteen NAC (NAM, ATAF1/2, CUC2) transcription factor genes of Populus trichocarpa were analyzed by 3′-RACE-seq. Based on experimental reproducibility, boundary sequence analysis and reverse transcription PCR (RT-PCR) verification, only canonical SSs were considered to be authentic. Based on stringent criteria, candidate PASs without any internal priming features were chosen as authentic PASs and assumed to be PAS-rich markers. Thirty-four novel canonical SSs, six intronic/internal exons and thirty 3′-UTR PAS-rich markers were revealed by 3′-RACE-seq. Using 3′-RACE and real-time PCR, we confirmed that three APA transcripts ending in/around PAS-rich markers were differentially regulated in response to plant hormones. Our results indicate that 3′-RACE-seq is a robust and cost-effective method to discover SSs and label active regions subjected to APA for genes or gene families. The method is suitable for small-scale AS and APA research in the initial stage.

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.

Differential expression of two galactinol synthase isoforms LcGolS1 and LcGolS2 in developing lentil (Lens culinaris Medik. cv CDC Redberry) seeds

Galactinol synthase (GS, EC 2.4.1.123) catalyzes the transfer of a galactosyl residue from UDP-galactose to myo-inositol to synthesize galactinol, a precursor for raffinose family oligosaccharides (RFO) biosynthesis. Screening, a cDNA library constructed with RNA isolated from developing lentil seeds, with partial GS genes resulted in identification of cDNA clones for two isoforms of GS, LcGolS1 (1336 bp, ORF-1002 bp, 334 amino acids) and LcGolS2 (1324bp, ORF-975bp, 325 amino acids) with predicted molecular weights of 38.7 kDa and 37.6 kDa, respectively. During lentil seed development, LcGolS1 transcripts showed higher accumulation during 26-32 days after flowering (DAF) corresponding to seed desiccation, while LcGolS2 showed maximum accumulation at 24 DAF, prior to increase in LcGolS1 transcripts. GS enzyme activity was maximum at 26 and 28 DAF and corresponded to galactinol accumulation, which also increased rapidly at 22 DAF with maximum accumulation at 26 DAF. Substrates for GS activity, myo-inositol and glucose/galactose were present in high concentrations during early stages of seed development but gradually decreased from 20 DAF to 32 DAF when galactinol concentration increased coinciding with increased GS enzyme activity.

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.

Experimental Genome-Wide Determination of RNA Polyadenylation in Chlamydomonas reinhardtii

The polyadenylation of RNA is a near-universal feature of RNA metabolism in eukaryotes. This process has been studied in the model alga Chlamydomonas reinhardtii using low-throughput (gene-by-gene) and high-throughput (transcriptome sequencing) approaches that recovered poly(A)-containing sequence tags which revealed interesting features of this critical process in Chlamydomonas. In this study, RNA polyadenylation has been studied using the so-called Poly(A) Tag Sequencing (PAT-Seq) approach. Specifically, PAT-Seq was used to study poly(A) site choice in cultures grown in four different media types—Tris-Phosphate (TP), Tris-Phosphate-Acetate (TAP), High-Salt (HS), and High-Salt-Acetate (HAS). The results indicate that: 1. As reported before, the motif UGUAA is the primary, and perhaps sole, cis-element that guides mRNA polyadenylation in the nucleus; 2. The scope of alternative polyadenylation events with the potential to change the coding sequences of mRNAs is limited; 3. Changes in poly(A) site choice in cultures grown in the different media types are very few in number and do not affect protein-coding potential; 4. Organellar polyadenylation is considerable and affects primarily ribosomal RNAs in the chloroplast and mitochondria; and 5. Organellar RNA polyadenylation is a dynamic process that is affected by the different media types used for cell growth.

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.

PASPA: a web server for mRNA poly(A) site predictions in plants and algae

MOTIVATION

Polyadenylation is an essential process during eukaryotic gene expression. Prediction of poly(A) sites helps to define the 3' end of genes, which is important for gene annotation and elucidating gene regulation mechanisms. However, due to limited knowledge of poly(A) signals, it is still challenging to predict poly(A) sites in plants and algae. PASPA is a web server for P: oly( A: ) S: ite prediction in P: lants and A: lgae, which integrates many in-house tools as add-ons to facilitate poly(A) site prediction, visualization and mining. This server can predict poly(A) sites for ten species, including seven previously poly(A) signal non-characterized species, with sensitivity and specificity in a range between 0.80 and 0.95.