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

Nucleotide-resolution Mapping of RNA N6-Methyladenosine (m6A) modifications and comprehensive analysis of global polyadenylation events in mRNA 3' end processing in malaria pathogen Plasmodium falciparum

Plasmodium falciparum is an obligate human parasite of the phylum Apicomplexa and is the causative agent of the most lethal form of human malaria. Although N6-methyladenosine modification is thought to be one of the major post-transcriptional regulatory mechanisms for stage-specific gene expression in apicomplexan parasites, the precise base position of m6A in mRNAs or noncoding RNAs in these parasites remains unknown. Here, we report global nucleotide-resolution mapping of m6A residues in P. falciparum using DART-seq technology, which quantitatively displayed a stage-specific, dynamic distribution pattern with enrichment near mRNA 3' ends. In this process we identified 894, 788, and 1,762 m6A-modified genes in Ring, Trophozoite and Schizont stages respectively, with an average of 5-7 m6A sites per-transcript at the individual gene level. Notably, several genes involved in malaria pathophysiology, such as KAHRP, ETRAMPs, SERA and stress response genes, such as members of Heat Shock Protein (HSP) family are highly enriched in m6A and therefore could be regulated by this RNA modification. Since we observed preferential methylation at the 3' ends of P. falciparum transcripts and because malaria polyadenylation specificity factor PfCPSF30 harbors an m6A reader 'YTH' domain, we reasoned that m6A might play an important role in 3'-end processing of malaria mRNAs. To investigate this, we used two complementary high-throughput RNA 3'-end mapping approaches, which provided an initial framework to explore potential roles of m6A in the regulation of alternative polyadenylation (APA) during malaria development in human hosts.

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.

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.

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.

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.

Transcriptional dynamics in the protozoan parasite Sarcocystis neurona and mammalian host cells after treatment with a specific inhibitor of apicomplexan mRNA polyadenylation

In recent years, a class of chemical compounds (benzoxaboroles) that are active against a range of parasites has been shown to target mRNA polyadenylation by inhibiting the activity of CPSF73, the endonucleolytic core of the eukaryotic polyadenylation complex. One particular compound, termed AN3661, is active against several apicomplexan parasites that cause disease in humans. In this study, we report that AN3661 is active against an apicomplexan that causes disease in horses and marine mammals (Sarcocystis neurona), with an approximate IC50 value of 14.99 nM. Consistent with the reported mode of action of AN3661 against other apicomplexans, S. neurona mutants resistant to AN3661 had an alteration in CPSF73 that was identical to a mutation previously documented in AN3661-resistant Toxoplasma gondii and Plasmodium falciparum. AN3661 had a wide-ranging effect on poly(A) site choice in S. neurona, with more than half of all expressed genes showing some alteration in mRNA 3' ends. This was accompanied by changes in the relative expression of more than 25% of S. neurona genes and an overall 5-fold reduction of S. neurona transcripts in infected cells. In contrast, AN3661 had no discernible effect on poly(A) site choice or gene expression in the host cells. These transcriptomic studies indicate that AN3661 is exceedingly specific for the parasite CPSF73 protein, and has the potential to augment other therapies for the control of apicomplexan parasites in domestic animals.

QPAT-seq, a rapid and deduplicatable method for quantification of poly(A) site usages

Alternative polyadenylation (APA) is an essential regulatory mechanism for gene expression. The next generation sequencing provides ample opportunity to precisely delineate APA sites genome-wide. Various methods for profiling transcriptome-wide poly(A) sites were developed. By comparing available methods, the ways for adding sequencing adaptors to fit with the Illumina sequencing platform are different. These methods have identified more than 50% genes that undergo APA in eukaryotes. However, due to the unbalanced PCR during library preparation, accurate quantification of poly(A) sites is still a challenge. Here, we describe an updated poly(A) tag sequencing method that incorporates unique molecular identifier (UMI) into the adaptor for removing quantification bias induced by PCR duplicates. Hence, quantification of poly(A) site usages can be achieved by counting UMIs. This protocol, quantifying poly(A) tag sequencing (QPAT-seq), can be finished in 1 day with reduced cost, and is particularly useful for application with a large number of samples.

Genome-wide distribution and functions of the AAE complex in epigenetic regulation in Arabidopsis

Heterochromatin is widespread in eukaryotic genomes and has diverse impacts depending on its genomic context. Previous studies have shown that a protein complex, the ASI1-AIPP1-EDM2 (AAE) complex, participates in polyadenylation regulation of several intronic heterochromatin-containing genes. However, the genome-wide functions of AAE are still unknown. Here, we show that the ASI1 and EDM2 mostly target the common genomic regions on a genome-wide level and preferentially interacts with genetic heterochromatin. Polyadenylation (poly(A) sequencing reveals that AAE complex has a substantial influence on poly(A) site usage of heterochromatin-containing genes, including not only intronic heterochromatin-containing genes but also the genes showing overlap with heterochromatin. Intriguingly, AAE is also involved in the alternative splicing regulation of a number of heterochromatin-overlapping genes, such as the disease resistance gene RPP4. We provided evidence that genic heterochromatin is indispensable for the recruitment of AAE in polyadenylation and splicing regulation. In addition to conferring RNA processing regulation at genic heterochromatin-containing genes, AAE also targets some transposable elements (TEs) outside of genes (including TEs sandwiched by genes and island TEs) for epigenetic silencing. Our results reveal new functions of AAE in RNA processing and epigenetic silencing, and thus represent important advances in epigenetic regulation.

HDA6-dependent histone deacetylation regulates mRNA polyadenylation in Arabidopsis

Eukaryotic histone deacetylation, critical for maintaining nucleosome structure and regulating gene expression, is mediated by histone deacetylases (HDACs). Although nucleosomes have been reported to regulate mRNA polyadenylation in humans, the role of HDACs in regulating polyadenylation has not been uncovered. Taking advantage of phenotypic studies on Arabidopsis, HDA6 (one of HDACs) was found to be a critical part of many biological processes. Here, we report that HDA6 affects mRNA polyadenylation in Arabidopsis. Poly(A) sites of up-regulated transcripts are closer to the histone acetylation peaks in hda6 compared to the wild-type Col-0. HDA6 is required for the deacetylation of histones around DNA on nucleosomes, which solely coincides with up-regulated or uniquely presented poly(A) sites in hda6. Furthermore, defective HDA6 results in an overrepresentation of the canonical poly(A) signal (AAUAAA) usage. Chromatin loci for generating AAUAAA-type transcripts have a comparatively low H3K9K14ac around poly(A) sites when compared to other noncanonical poly(A) signal–containing transcripts. These results indicate that HDA6 regulates polyadenylation in a histone deacetylation–dependent manner in Arabidopsis.

Histone Deacetylase 11 Knockdown Blocks Larval Development and Metamorphosis in the Red Flour Beetle, Tribolium castaneum

Post-translational modifications (PTM) such as methylation, acetylation, phosphorylation, and ubiquitination of histones and other proteins regulate expression of genes. The acetylation levels of these proteins are determined by the balance of expression of histone acetyltransferase (HATs) and histone deacetylases (HDACs). We recently reported that class I HDACs (HDAC1 and HDAC3) play important roles in juvenile hormone (JH) suppression of metamorphosis in the red flour beetle, Tribolium castaneum. Here, we report on the function of a single class IV HDAC member, HDAC11. Injection of dsRNA targeting T. castaneum HDAC11 gene into newly molted last instar larvae induced knockdown of the target gene and arrested larval development and prevented metamorphosis into the pupal stage. Dark melanized areas were detected in larvae that showed developmental arrest and mortality. Developmental expression studies showed an increase in HDAC11 mRNA levels beginning at the end of the penultimate larval stage. These higher levels were maintained during the final instar larval and pupal stages. A JH analog, hydroprene, suppressed HDAC11 expression in the larvae. Sequencing of RNA isolated from control and dsHDAC11 injected larvae identified several differentially expressed genes, including those involved in JH action, ecdysone response, and melanization. The acetylation levels of core histones showed an increase in TcA cells exposed to dsHDAC11. Also, an increase in histone H3 acetylation, specifically H3K9, H3K18 and H3K27, were detected in HDAC11 knockdown larvae. These studies report the function of HDAC11 in insects other than Drosophila for the first time and show that HDAC11 influences the acetylation levels of histones and expression of multiple genes involved in T. castaneum larval development.

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.

Epigenetic modifications acetylation and deacetylation play important roles in juvenile hormone action

BACKGROUND

Epigenetic modifications including DNA methylation and post-translational modifications of histones are known to regulate gene expression. Antagonistic activities of histone acetyltransferases (HATs) and histone deacetylases (HDACs) mediate transcriptional reprogramming during insect development as shown in Drosophila melanogaster and other insects. Juvenile hormones (JH) play vital roles in the regulation of growth, development, metamorphosis, reproduction and other physiological processes. However, our current understanding of epigenetic regulation of JH action is still limited. Hence, we studied the role of CREB binding protein (CBP, contains HAT domain) and Trichostatin A (TSA, HDAC inhibitor) on JH action.

RESULTS

Exposure of Tribolium castaneum cells (TcA cells) to JH or TSA caused an increase in expression of Kr-h1 (a known JH-response gene) and 31 or 698 other genes respectively. Knockdown of the gene coding for CBP caused a decrease in the expression of 456 genes including Kr-h1. Interestingly, the expression of several genes coding for transcription factors, nuclear receptors, P450 and fatty acid synthase family members that are known to mediate JH action were affected by CBP knockdown or TSA treatment.

CONCLUSIONS

These data suggest that acetylation and deacetylation mediated by HATs and HDACs play an important role in JH action.

Alternative polyadenylation analysis in animals and plants: newly developed strategies for profiling, processing and validation

Alternative polyadenylation is an essential RNA processing event that contributes significantly to regulation of transcriptome diversity and functional dynamics in both animals and plants. Here we review newly developed next generation sequencing methods for genome-wide profiling of alternative polyadenylation (APA) sites, bioinformatics pipelines for data processing and both wet and dry laboratory approaches for APA validation. The library construction methods LITE-Seq (Low-Input 3'-Terminal sequencing) and PAC-seq (PolyA Click sequencing) tag polyA+ cDNA, while BAT-seq (BArcoded, three-prime specific sequencing) and PAPERCLIP (Poly(A) binding Protein-mediated mRNA 3'End Retrieval by CrossLinking ImmunoPrecipitation) enrich polyA+ RNA. Interestingly, only WTTS-seq (Whole Transcriptome Termini Site sequencing) targets both polyA+ RNA and polyA+ cDNA. Varieties of bioinformatics pipelines are well established to pursue read quality control, mapping, clustering, characterization and pathway analysis. The RHAPA (RNase H alternative polyadenylation assay) and 3'RACE-seq (3' rapid amplification of cDNA end sequencing) methods directly validate APA sites, while WTSS-seq (whole transcriptome start site sequencing), RNA-seq (RNA sequencing) and public APA databases can serve as indirect validation methods. We hope that these tools, pipelines and resources trigger huge waves of interest in the research community to investigate APA events underlying physiological, pathological and psychological changes and thus understand the information transfer events from genome to phenome relevant to economically important traits in both animals and plants.

Characterization of mRNA polyadenylation in the apicomplexa

Messenger RNA polyadenylation is a universal aspect of gene expression in eukaryotes. In well-established model organisms, this process is mediated by a conserved complex of 15–20 subunits. To better understand this process in apicomplexans, a group of unicellular parasites that causes serious disease in humans and livestock, a computational and high throughput sequencing study of the polyadenylation complex and poly(A) sites in several species was conducted. BLAST-based searches for orthologs of the human polyadenylation complex yielded clear matches to only two—poly(A) polymerase and CPSF73—of the 19 proteins used as queries in this analysis. As the human subunits that recognize the AAUAAA polyadenylation signal (PAS) were not immediately obvious, a computational analysis of sequences adjacent to experimentally-determined apicomplexan poly(A) sites was conducted. The results of this study showed that there exists in apicomplexans an A-rich region that corresponds in position to the AAUAAA PAS. The set of experimentally-determined sites in one species, Sarcocystis neurona, was further analyzed to evaluate the extent and significance of alternative poly(A) site choice in this organism. The results showed that almost 80% of S. neurona genes possess more than one poly(A) site, and that more than 780 sites showed differential usage in the two developmental stages–extracellular merozoites and intracellular schizonts–studied. These sites affected more than 450 genes, and included a disproportionate number of genes that encode membrane transporters and ribosomal proteins. Taken together, these results reveal that apicomplexan species seem to possess a poly(A) signal analogous to AAUAAA even though genes that may encode obvious counterparts of the AAUAAA-recognizing proteins are absent in these organisms. They also indicate that, as is the case in other eukaryotes, alternative polyadenylation is a widespread phenomenon in S. neurona that has the potential to impact growth and development.

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.

An integrative model for alternative polyadenylation, IntMAP, delineates mTOR-modulated endoplasmic reticulum stress response

3'-untranslated regions (UTRs) can vary through the use of alternative polyadenylation sites during pre-mRNA processing. Multiple publically available pipelines combining high profiling technologies and bioinformatics tools have been developed to catalog changes in 3'-UTR lengths. In our recent RNA-seq experiments using cells with hyper-activated mammalian target of rapamycin (mTOR), we found that cellular mTOR activation leads to transcriptome-wide alternative polyadenylation (APA), resulting in the activation of multiple cellular pathways. Here, we developed a novel bioinformatics algorithm, IntMAP, which integrates RNA-Seq and PolyA Site (PAS)-Seq data for a comprehensive characterization of APA events. By applying IntMAP to the datasets from cells with hyper-activated mTOR, we identified novel APA events that could otherwise not be identified by either profiling method alone. Several transcription factors including Cebpg (CCAAT/enhancer binding protein gamma) were among the newly discovered APA transcripts, indicating that diverse transcriptional networks may be regulated by mTOR-coordinated APA. The prevention of APA in Cebpg using the CRISPR/cas9-mediated genome editing tool showed that mTOR-driven 3'-UTR shortening in Cebpg is critical in protecting cells from endoplasmic reticulum (ER) stress. Taken together, we present IntMAP as a new bioinformatics algorithm for APA analysis by which we expand our understanding of the physiological role of mTOR-coordinated APA events to ER stress response. IntMAP toolbox is available at http://compbio.cs.umn.edu/IntMAP/.

Pipecolic acid confers systemic immunity by regulating free radicals

Pipecolic acid (Pip), a non-proteinaceous product of lysine catabolism, is an important regulator of immunity in plants and humans alike. In plants, Pip accumulates upon pathogen infection and has been associated with systemic acquired resistance (SAR). However, the molecular mechanisms underlying Pip-mediated signaling and its relationship to other known SAR inducers remain unknown. We show that in plants, Pip confers SAR by increasing levels of the free radicals, nitric oxide (NO), and reactive oxygen species (ROS), which act upstream of glycerol-3-phosphate (G3P). Plants defective in NO, ROS, G3P, or salicylic acid (SA) biosynthesis accumulate reduced Pip in their distal uninfected tissues although they contain wild-type-like levels of Pip in their infected leaves. These data indicate that de novo synthesis of Pip in distal tissues is dependent on both SA and G3P and that distal levels of SA and G3P play an important role in SAR. These results also suggest a unique scenario whereby metabolites in a signaling cascade can stimulate each other's biosynthesis depending on their relative levels and their site of action.

Multiple functions of CREB-binding protein during postembryonic development: identification of target genes

Background

Juvenile hormones (JH) and ecdysteroids control postembryonic development in insects. They serve as valuable targets for pest management. Hence, understanding the molecular mechanisms of their action is of crucial importance. CREB-binding protein (CBP) is a universal transcriptional co-regulator. It controls the expression of several genes including those from hormone signaling pathways through co-activation of many transcription factors. However, the role of CBP during postembryonic development in insects is not well understood. Therefore, we have studied the role of CBP in postembryonic development in Tribolium, a model coleopteran insect.

Results

CBP is ubiquitously expressed in the red flour beetle, Tribolium castaneum. RNA interference (RNAi) mediated knockdown of CBP resulted in a decrease in JH induction of Kr-h1 gene expression in Tribolium larvae and led to a block in their development. Moreover, the injection of CBP double-stranded RNA (dsRNA) showed lethal phenotypes within 8 days of injection. RNA-seq and subsequent differential gene expression analysis identified CBP target genes in Tribolium. Knockdown of CBP caused a decrease in the expression of 1306 genes coding for transcription factors and other proteins associated with growth and development. Depletion of CBP impaired the expression of several JH response genes (e.g., Kr-h1, Hairy, early trypsin) and ecdysone response genes (EcR, E74, E75, and broad complex). Further, GO enrichment analyses of the downregulated genes showed enrichment in different functions including developmental processes, pigmentation, anatomical structure development, regulation of biological and cellular processes, etc.

Conclusion

These data suggest diverse but crucial roles for CBP during postembryonic development in the coleopteran model insect, Tribolium. It can serve as a target for RNAi mediated pest management of this stored product pest.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4373-3) contains supplementary material, which is available to authorized users.

Alternative Polyadenylation: Methods, Findings, and Impacts

Alternative polyadenylation (APA), a phenomenon that RNA molecules with different 3' ends originate from distinct polyadenylation sites of a single gene, is emerging as a mechanism widely used to regulate gene expression. In the present review, we first summarized various methods prevalently adopted in APA study, mainly focused on the next-generation sequencing (NGS)-based techniques specially designed for APA identification, the related bioinformatics methods, and the strategies for APA study in single cells. Then we summarized the main findings and advances so far based on these methods, including the preferences of alternative polyA (pA) site, the biological processes involved, and the corresponding consequences. We especially categorized the APA changes discovered so far and discussed their potential functions under given conditions, along with the possible underlying molecular mechanisms. With more in-depth studies on extensive samples, more signatures and functions of APA will be revealed, and its diverse roles will gradually heave in sight.

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.

Comprehensive profiling of rhizome-associated alternative splicing and alternative polyadenylation in moso bamboo (Phyllostachys edulis)

Moso bamboo (Phyllostachys edulis) represents one of the fastest-spreading plants in the world, due in part to its well-developed rhizome system. However, the post-transcriptional mechanism for the development of the rhizome system in bamboo has not been comprehensively studied. We therefore used a combination of single-molecule long-read sequencing technology and polyadenylation site sequencing (PAS-seq) to re-annotate the bamboo genome, and identify genome-wide alternative splicing (AS) and alternative polyadenylation (APA) in the rhizome system. In total, 145 522 mapped full-length non-chimeric (FLNC) reads were analyzed, resulting in the correction of 2241 mis-annotated genes and the identification of 8091 previously unannotated loci. Notably, more than 42 280 distinct splicing isoforms were derived from 128 667 intron-containing full-length FLNC reads, including a large number of AS events associated with rhizome systems. In addition, we characterized 25 069 polyadenylation sites from 11 450 genes, 6311 of which have APA sites. Further analysis of intronic polyadenylation revealed that LTR/Gypsy and LTR/Copia were two major transposable elements within the intronic polyadenylation region. Furthermore, this study provided a quantitative atlas of poly(A) usage. Several hundred differential poly(A) sites in the rhizome-root system were identified. Taken together, these results suggest that post-transcriptional regulation may potentially have a vital role in the underground rhizome-root system.

Genomic and functional characterization of coleopteran insect-specific α-amylase inhibitor gene from Amaranthus species

The smallest 32 amino acid α-amylase inhibitor from Amaranthus hypochondriacus (AAI) is reported. The complete gene of pre-protein (AhAI) encoding a 26 amino acid (aa) signal peptide followed by the 43 aa region and the previously identified 32 aa peptide was cloned successfully. Three cysteine residues and one disulfide bond conserved within known α-amylase inhibitors were present in AhAI. Identical genomic and open reading frame was found to be present in close relatives of A. hypochondriacus namely Amaranthus paniculatus, Achyranthes aspera and Celosia argentea. Interestingly, the 3'UTR of AhAI varied in these species. The highest expression of AhAI was observed in A. hypochondriacus inflorescence; however, it was not detected in the seed. We hypothesized that the inhibitor expressed in leaves and inflorescence might be transported to the seeds. Sub-cellular localization studies clearly indicated the involvement of AhAI signal peptide in extracellular secretion. Full length rAhAI showed differential inhibition against α-amylases from human, insects, fungi and bacteria. Particularly, α-amylases from Helicoverpa armigera (Lepidoptera) were not inhibited by AhAI while Tribolium castaneum and Callosobruchus chinensis (Coleoptera) α-amylases were completely inhibited. Molecular docking of AhAI revealed tighter interactions with active site residues of T. castaneum α-amylase compared to C. chinensis α-amylase, which could be the rationale behind the disparity in their IC₅₀. Normal growth, development and adult emergence of C. chinensis were hampered after feeding on rAhAI. Altogether, the ability of AhAI to affect the growth of C. chinensis demonstrated its potential as an efficient bio-control agent, especially against stored grain pests.

Meeting report: processing, translation, decay - three ways to keep RNA sizzling

This meeting report highlights key trends that emerged from a conference entitled Post-Transcriptional Gene Regulation in Plants, which was held 14-15 July 2016, as a satellite meeting of the annual meeting of the American Society of Plant Biologists in Austin, Texas. The molecular biology of RNA is emerging as an integral part of the framework for plants' responses to environmental challenges such as drought and heat, hypoxia, nutrient deprivation, light and pathogens. Moreover, the conference illustrated how a multitude of customized and pioneering omics-related technologies are being applied, more and more often in combination, to describe and dissect the complexities of gene expression at the post-transcriptional level.

Accurate Profiling of Gene Expression and Alternative Polyadenylation with Whole Transcriptome Termini Site Sequencing (WTTS-Seq)

Construction of next-generation sequencing (NGS) libraries involves RNA manipulation, which often creates noisy, biased, and artifactual data that contribute to errors in transcriptome analysis. In this study, a total of 19 whole transcriptome termini site sequencing (WTTS-seq) and seven RNA sequencing (RNA-seq) libraries were prepared from Xenopus tropicalis adult and embryo samples to determine the most effective library preparation method to maximize transcriptomics investigation. We strongly suggest that appropriate primers/adaptors are designed to inhibit amplification detours and that PCR overamplification is minimized to maximize transcriptome coverage. Furthermore, genome annotation must be improved so that missing data can be recovered. In addition, a complete understanding of sequencing platforms is critical to limit the formation of false-positive results. Technically, the WTTS-seq method enriches both poly(A)+ RNA and complementary DNA, adds 5′- and 3′-adaptors in one step, pursues strand sequencing and mapping, and profiles both gene expression and alternative polyadenylation (APA). Although RNA-seq is cost prohibitive, tends to produce false-positive results, and fails to detect APA diversity and dynamics, its combination with WTTS-seq is necessary to validate transcriptome-wide APA.

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 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.

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.

Whole transcriptome analysis with sequencing: methods, challenges and potential solutions

Whole transcriptome analysis plays an essential role in deciphering genome structure and function, identifying genetic networks underlying cellular, physiological, biochemical and biological systems and establishing molecular biomarkers that respond to diseases, pathogens and environmental challenges. Here, we review transcriptome analysis methods and technologies that have been used to conduct whole transcriptome shotgun sequencing or whole transcriptome tag/target sequencing analyses. We focus on how adaptors/linkers are added to both 5' and 3' ends of mRNA molecules for cloning or PCR amplification before sequencing. Challenges and potential solutions are also discussed. In brief, next generation sequencing platforms have accelerated releases of the large amounts of gene expression data. It is now time for the genome research community to assemble whole transcriptomes of all species and collect signature targets for each gene/transcript, and thus use known genes/transcripts to determine known transcriptomes directly in the near future.

Alternative polyadenylation: less than meets the eye?

With the advances in deep-sequencing techniques over the last decade, the study of alternative cleavage and polyadenylation (APA) has shifted from individual gene to whole transcriptome analysis. Findings from such global studies have elevated APA to its currently accepted status as a major player in the regulation of eukaryotic gene expression. Although ~70% of human genes have been shown to contain multiple cleavage and polyadenylation sites, the extent of the consequences of APA and its role in regulating physiological processes are still largely unknown. The present review aims to summarize the experimental evidence that supports a physiological role of APA and highlights some of the shortcomings that need addressing to substantiate the widely proposed claim that APA is a key player in global gene regulation.

Poly(A)-tag deep sequencing data processing to extract poly(A) sites

Polyadenylation [poly(A)] is an essential posttranscriptional processing step in the maturation of eukaryotic mRNA. The advent of next-generation sequencing (NGS) technology has offered feasible means to generate large-scale data and new opportunities for intensive study of polyadenylation, particularly deep sequencing of the transcriptome targeting the junction of 3'-UTR and the poly(A) tail of the transcript. To take advantage of this unprecedented amount of data, we present an automated workflow to identify polyadenylation sites by integrating NGS data cleaning, processing, mapping, normalizing, and clustering. In this pipeline, a series of Perl scripts are seamlessly integrated to iteratively map the single- or paired-end sequences to the reference genome. After mapping, the poly(A) tags (PATs) at the same genome coordinate are grouped into one cleavage site, and the internal priming artifacts removed. Then the ambiguous region is introduced to parse the genome annotation for cleavage site clustering. Finally, cleavage sites within a close range of 24 nucleotides and from different samples can be clustered into poly(A) clusters. This procedure could be used to identify thousands of reliable poly(A) clusters from millions of NGS sequences in different tissues or treatments.

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.

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.