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

Insights into the role of N6-methyladenosine (m6A) in plant-virus interactions

N6-methyladenosine (m6A) is a common and dynamic epitranscriptomic modification in eukaryotic RNAs, affecting stability, splicing, translation, and degradation. Recent technological advancements have revealed the complex nature of m6A modifications, highlighting their importance in plant and animal species. The m6A modification is a reversible process, with "writers" depositing methylation, "erasers" demethylating it, and "reader" proteins recognizing m6A and executing various biological functions. Studying the relationship between m6A methylation and viral infection is crucial. Animal viruses, including retroviruses, RNA viruses, and DNA viruses, often employ the host's m6A machinery to replicate or avoid immune responses. In plant viruses, host methyltransferases or demethylases can stabilize or degrade viral RNA, depending on the virus-host interaction. Additionally, viral infections can modify the host's m6A machinery, impacting the viral life cycle. This review examines the role of m6A modifications in plant viral pathogenesis, focussing on RNA viruses infecting crops like alfalfa, turnip, wheat, rice, and potato. Understanding the role of m6A in virus-host interactions can aid in studying plant viral disease development and discovering novel antiviral targets for crop protection. In this review, we summarize current information on m6A in RNA biology, focussing on its function in viral infections and plant-virus interactions.

Characterizing the impact of CPSF30 gene disruption on TuMV infection in Arabidopsis thaliana

CPSF30, a key polyadenylation factor, also serves as an m6A reader, playing a crucial role in determining RNA fate post-transcription. While its homologs mammals are known to be vital for viral replication and immune evasion, the full scope of CPSF30 in plant, particular in viral regulation, remains less explored. Our study demonstrates that CPSF30 significantly facilitates the infection of turnip mosaic virus (TuMV) in Arabidopsis thaliana, as evidenced by infection experiments on the engineered cpsf30 mutant. Among the two isoforms, CPSF30-L, which were characterized with m6A binding activity, emerged as the primary isoform responding to TuMV infection. Analysis of m6A components revealed potential involvement of the m6A machinery in regulating TuMV infection. In contrast, CPSF30-S exhibited distinct subcellular localization, coalescing with P-body markers (AtDCP1 and AtDCP2) in cytoplasmic granules, suggesting divergent regulatory mechanisms between the isoforms. Furthermore, comprehensive mRNA-Seq and miRNA-Seq analysis of Col-0 and cpsf30 mutants revealed global transcriptional reprogramming, highlighting CPSF30's role in selectively modulating gene expression during TuMV infection. In conclusion, this research underscores CPSF30's critical role in the TuMV lifecycle and sets the stage for further exploration of its function in plant viral regulation.

Functional Analysis of CPSF30 in Nilaparvata lugens Using RNA Interference Reveals Its Essential Role in Development and Survival

The brown planthopper (Nilaparvata lugens) is a major pest threatening global rice production, significantly reducing yields annually. As N. lugens increasingly develops resistance to conventional control methods, such as chemical pesticides, there is an urgent need for innovative and sustainable pest management strategies. Cleavage and Polyadenylation Specificity Factor 30 (CPSF30) is a key protein involved in mRNA 3' end processing, yet its function in N. lugens remains poorly understood. This study aims to elucidate the role of CPSF30 in the growth and development of N. lugens and evaluate its potential as a target for RNA interference (RNAi)-based pest control strategies. We cloned and characterized the cDNA sequence of NlCPSF30, which encodes a protein of 341 amino acids containing five CCCH zinc-finger domains and two CCHC zinc-knuckle domains. Sequence alignment revealed that NlCPSF30 is highly conserved among insect species, particularly in the zinc-finger domains essential for RNA binding and processing. Phylogenetic analysis showed that NlCPSF30 is closely related to CPSF30 proteins from other hemipteran species. Expression analysis indicated that NlCPSF30 is most highly expressed in the fat body and during the adult stage, with significantly higher expression in females than in males. RNAi-mediated silencing of NlCPSF30 in third-instar nymphs resulted in severe phenotypic abnormalities, including disrupted molting and increased mortality following injection of double-stranded RNA (dsRNA) targeting NlCPSF30. Moreover, it influenced the expression of genes associated with hormone regulation, namely NlHry, NlE93, and NlKr-h1. These results suggest that NlCPSF30 is integral to critical physiological processes, with its disruption leading to increased mortality. Our findings identify NlCPSF30 as an essential gene for N. lugens' survival and a promising target for RNAi-based pest management strategies. This study provides a valuable molecular target and theoretical insights for developing RNAi-based control methods against N. lugens.

Dynamics of alternative polyadenylation in single root cells of Arabidopsis thaliana

Introduction

Single-cell RNA-seq (scRNA-seq) technologies have been widely used to reveal the diversity and complexity of cells, and pioneering studies on scRNA-seq in plants began to emerge since 2019. However, existing studies on plants utilized scRNA-seq focused only on the gene expression regulation. As an essential post-transcriptional mechanism for regulating gene expression, alternative polyadenylation (APA) generates diverse mRNA isoforms with distinct 3' ends through the selective use of different polyadenylation sites in a gene. APA plays important roles in regulating multiple developmental processes in plants, such as flowering time and stress response.

Methods

In this study, we developed a pipeline to identify and integrate APA sites from different scRNA-seq data and analyze APA dynamics in single cells. First, high-confidence poly(A) sites in single root cells were identified and quantified. Second, three kinds of APA markers were identified for exploring APA dynamics in single cells, including differentially expressed poly(A) sites based on APA site expression, APA markers based on APA usages, and APA switching genes based on 3' UTR (untranslated region) length change. Moreover, cell type annotations of single root cells were refined by integrating both the APA information and the gene expression profile.

Results

We comprehensively compiled a single-cell APA atlas from five scRNA-seq studies, covering over 150,000 cells spanning four major tissue branches, twelve cell types, and three developmental stages. Moreover, we quantified the dynamic APA usages in single cells and identified APA markers across tissues and cell types. Further, we integrated complementary information of gene expression and APA profiles to annotate cell types and reveal subtle differences between cell types.

Discussion

This study reveals that APA provides an additional layer of information for determining cell identity and provides a landscape of APA dynamics during Arabidopsis root development.

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.

Global analysis of key post-transcriptional regulation in early leaf development of Limonium bicolor identifies a long non-coding RNA that promotes salt gland development and salt resistance

Limonium bicolor, known horticulturally as sea lavender, is a typical recretohalophyte with salt glands in its leaf epidermis that secrete excess Na+ out of the plant. Although many genes have been proposed to contribute to salt gland initiation and development, a detailed analysis of alternative splicing, alternative polyadenylation patterns, and long non-coding RNAs (lncRNAs) has been lacking. Here, we applied single-molecule long-read mRNA isoform sequencing (Iso-seq) to explore the complexity of the L. bicolor transcriptome in leaves during salt gland initiation (stage A) and salt gland differentiation (stage B) based on the reference genome. We identified alternative splicing events and the use of alternative poly(A) sites unique to stage A or stage B, leading to the hypothesis that they might contribute to the differentiation of salt glands. Based on the Iso-seq data and RNA in situ hybridization of candidate genes, we selected the lncRNA Btranscript_153392 for gene editing and virus-induced gene silencing to dissect its function. In the absence of this transcript, we observed fewer salt glands on the leaf epidermis, leading to diminished salt secretion and salt tolerance. Our data provide transcriptome resources for unraveling the mechanisms behind salt gland development and furthering crop transformation efforts towards enhanced survivability in saline soils.

Light regulates widespread plant alternative polyadenylation through the chloroplast

Transcription of eukaryotic protein-coding genes generates immature mRNAs that are subjected to a series of processing events, including capping, splicing, cleavage, and polyadenylation (CPA), and chemical modifications of bases. Alternative polyadenylation (APA) greatly contributes to mRNA diversity in the cell. By determining the length of the 3' untranslated region, APA generates transcripts with different regulatory elements, such as miRNA and RBP binding sites, which can influence mRNA stability, turnover, and translation. In the model plant Arabidopsis thaliana, APA is involved in the control of seed dormancy and flowering. In view of the physiological importance of APA in plants, we decided to investigate the effects of light/dark conditions and compare the underlying mechanisms to those elucidated for alternative splicing (AS). We found that light controls APA in approximately 30% of Arabidopsis genes. Similar to AS, the effect of light on APA requires functional chloroplasts, is not affected in mutants of the phytochrome and cryptochrome photoreceptor pathways, and is observed in roots only when the communication with the photosynthetic tissues is not interrupted. Furthermore, mitochondrial and TOR kinase activities are necessary for the effect of light. However, unlike AS, coupling with transcriptional elongation does not seem to be involved since light-dependent APA regulation is neither abolished in mutants of the TFIIS transcript elongation factor nor universally affected by chromatin relaxation caused by histone deacetylase inhibition. Instead, regulation seems to correlate with changes in the abundance of constitutive CPA factors, also mediated by the chloroplast.

Arabidopsis and maize terminator strength is determined by GC content, polyadenylation motifs and cleavage probability

The 3' end of a gene, often called a terminator, modulates mRNA stability, localization, translation, and polyadenylation. Here, we adapted Plant STARR-seq, a massively parallel reporter assay, to measure the activity of over 50,000 terminators from the plants Arabidopsis thaliana and Zea mays. We characterize thousands of plant terminators, including many that outperform bacterial terminators commonly used in plants. Terminator activity is species-specific, differing in tobacco leaf and maize protoplast assays. While recapitulating known biology, our results reveal the relative contributions of polyadenylation motifs to terminator strength. We built a computational model to predict terminator strength and used it to conduct in silico evolution that generated optimized synthetic terminators. Additionally, we discover alternative polyadenylation sites across tens of thousands of terminators; however, the strongest terminators tend to have a dominant cleavage site. Our results establish features of plant terminator function and identify strong naturally occurring and synthetic terminators.

Identification of m6A RNA methylation genes in Oryza sativa and expression profiling in response to different developmental and environmental stimuli

Eukaryotic messenger RNAs (mRNAs) transcend their predominant function of protein encoding by incorporating auxiliary components that ultimately contribute to their processing, transportation, translation, and decay. In doing so, additional layers of modifications are incorporated in mRNAs at post-transcriptional stage. Among them, N6-methyladenosine (m6A) is the most frequently found mRNA modification that plays crucial roles in plant development and stress response. In the overall mechanism of m6A methylation, key proteins classified based on their functions such as writers, readers, and erasers dynamically add, read, and subtract methyl groups respectively to deliver relevant functions in response to external stimuli. In this study, we identified 30 m6A regulatory genes (9 writers, 5 erasers, and 16 readers) in rice that encode 53 proteins (13 writers, 7 erasers, and 33 readers) where segmental duplication was found in one writer and four reader gene pairs. Reproductive cells such as sperm, anther and panicle showed high levels of expression for most of the m6A regulatory genes. Notably, writers like OsMTA, OsMTD, and OsMTC showed varied responses in different stress and infection contexts, with initial upregulation in response to early exposure followed by downregulation later. OsALKBH9A, a noteworthy eraser, displayed varied expression in response to different stresses at different time intervals, but upregulation in certain infections. Reader genes like OsECT5, OsCPSF30-L3, and OsECT8 showed continuous upregulation in exertion of all kinds of stress relevant here. Conversely, other reader genes along with OsECT11 and OsCPSF30-L2 were observed to be consistently downregulated. The apparent correlation between the expression patterns of m6A regulatory genes and stress modulation pathways in this study underscores the need for additional research to unravel their intricate regulatory mechanisms that could ultimately contribute to the substantial development of enhanced stress tolerance in rice through mRNA modification.

Alternative polyadenylation profiles of susceptible and resistant rice (Oryza sativa L.) in response to bacterial leaf blight using RNA-seq

BACKGROUND

Alternative polyadenylation (APA) is an important pattern of post-transcriptional regulation of genes widely existing in eukaryotes, involving plant physiological and pathological processes. However, there is a dearth of studies investigating the role of APA profile in rice leaf blight.

RESULTS

In this study, we compared the APA profile of leaf blight-susceptible varieties (CT 9737-613P-M) and resistant varieties (NSIC RC154) following bacterial blight infection. Through gene enrichment analysis, we found that the genes of two varieties typically exhibited distal poly(A) (PA) sites that play different roles in two kinds of rice, indicating differential APA regulatory mechanisms. In this process, many disease-resistance genes displayed multiple transcripts via APA. Moreover, we also found five polyadenylation factors of similar expression patterns of rice, highlighting the critical roles of these five factors in rice response to leaf blight about PA locus diversity.

CONCLUSION

Notably, the present study provides the first dynamic changes of APA in rice in early response to biotic stresses and proposes a possible functional conjecture of APA in plant immune response, which lays the theoretical foundation for in-depth determination of the role of APA events in plant stress response and other life processes.

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.

Arabidopsis and Maize Terminator Strength is Determined by GC Content, Polyadenylation Motifs and Cleavage Probability

The 3' end of a gene, often called a terminator, modulates mRNA stability, localization, translation, and polyadenylation. Here, we adapted Plant STARR-seq, a massively parallel reporter assay, to measure the activity of over 50,000 terminators from the plants Arabidopsis thaliana and Zea mays. We characterize thousands of plant terminators, including many that outperform bacterial terminators commonly used in plants. Terminator activity is species-specific, differing in tobacco leaf and maize protoplast assays. While recapitulating known biology, our results reveal the relative contributions of polyadenylation motifs to terminator strength. We built a computational model to predict terminator strength and used it to conduct in silico evolution that generated optimized synthetic terminators. Additionally, we discover alternative polyadenylation sites across tens of thousands of terminators; however, the strongest terminators tend to have a dominant cleavage site. Our results establish features of plant terminator function and identify strong naturally occurring and synthetic terminators.

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.

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.

Coupling of co-transcriptional splicing and 3' end Pol II pausing during termination in Arabidopsis

BACKGROUND

In Arabidopsis, RNA Polymerase II (Pol II) often pauses within a few hundred base pairs downstream of the polyadenylation site, reflecting efficient transcriptional termination, but how such pausing is regulated remains largely elusive.

RESULT

Here, we analyze Pol II dynamics at 3' ends by combining comprehensive experiments with mathematical modelling. We generate high-resolution serine 2 phosphorylated (Ser2P) Pol II positioning data specifically enriched at 3' ends and define a 3' end pause index (3'PI). The position but not the extent of the 3' end pause correlates with the termination window size. The 3'PI is not decreased but even mildly increased in the termination deficient mutant xrn3, indicating 3' end pause is a regulatory step early during the termination and before XRN3-mediated RNA decay that releases Pol II. Unexpectedly, 3'PI is closely associated with gene exon numbers and co-transcriptional splicing efficiency. Multiple exons genes often display stronger 3' end pauses and more efficient on-chromatin splicing than genes with fewer exons. Chemical inhibition of splicing strongly reduces the 3'PI and disrupts its correlation with exon numbers but does not globally impact 3' end readthrough levels. These results are further confirmed by fitting Pol II positioning data with a mathematical model, which enables the estimation of parameters that define Pol II dynamics.

CONCLUSION

Our work highlights that the number of exons via co-transcriptional splicing is a major determinant of Pol II pausing levels at the 3' end of genes in plants.

Inhibition of DNA and RNA methylation disturbs root development of moso bamboo

DNA methylation (5mC) and N6-methyladenosine (m6A) are two important epigenetics regulators, which have a profound impact on plant growth development. Phyllostachys edulis (P. edulis) is one of the fastest spreading plants due to its well-developed root system. However, the association between 5mC and m6A has seldom been reported in P. edulis. In particular, the connection between m6A and several post-transcriptional regulators remains uncharacterized in P. edulis. Here, our morphological and electron microscope observations showed the phenotype of increased lateral root under RNA methylation inhibitor (DZnepA) and DNA methylation inhibitor (5-azaC) treatment. RNA epitranscriptome based on Nanopore direct RNA sequencing revealed that DZnepA treatment exhibits significantly decreased m6A level in the 3'-untranslated region (3'-UTR), which was accompanied by increased gene expression, full-length ratio, higher proximal poly(A) site usage and shorter poly(A) tail length. DNA methylation levels of CG and CHG were reduced in both coding sequencing and transposable element upon 5-azaC treatment. Cell wall synthesis was impaired under methylation inhibition. In particular, differentially expressed genes showed a high percentage of overlap between DZnepA and 5-azaC treatment, which suggested a potential correlation between two methylations. This study provides preliminary information for a better understanding of the link between m6A and 5mC in root development of moso bamboo.

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 3'-untranslated regions regulate high-salt tolerance of Spartina alterniflora

High-salt stress continues to challenge the growth and survival of many plants. Alternative polyadenylation (APA) produces mRNAs with different 3'-untranslated regions (3' UTRs) to regulate gene expression at the post-transcriptional level. However, the roles of alternative 3' UTRs in response to salt stress remain elusive. Here, we report the function of alternative 3' UTRs in response to high-salt stress in S. alterniflora (Spartina alterniflora), a monocotyledonous halophyte tolerant of high-salt environments. We found that high-salt stress induced global APA dynamics, and ∼42% of APA genes responded to salt stress. High-salt stress led to 3' UTR lengthening of 207 transcripts through increasing the usage of distal poly(A) sites. Transcripts with alternative 3' UTRs were mainly enriched in salt stress-related ion transporters. Alternative 3' UTRs of HIGH-AFFINITY K+ TRANSPORTER 1 (SaHKT1) increased RNA stability and protein synthesis in vivo. Regulatory AU-rich elements were identified in alternative 3' UTRs, boosting the protein level of SaHKT1. RNAi-knock-down experiments revealed that the biogenesis of 3' UTR lengthening in SaHKT1 was controlled by the poly(A) factor CLEAVAGE AND POLYADENYLATION SPECIFICITY FACTOR 30 (SaCPSF30). Over-expression of SaHKT1 with an alternative 3' UTR in rice (Oryza sativa) protoplasts increased mRNA accumulation of salt-tolerance genes in an AU-rich element-dependent manner. These results suggest that mRNA 3' UTR lengthening is a potential mechanism in response to high-salt stress. These results also reveal complex regulatory roles of alternative 3' UTRs coupling APA and regulatory elements at the post-transcriptional level in plants.

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.

Drought induces epitranscriptome and proteome changes in stem-differentiating xylem of Populus trichocarpa

Understanding gene expression and regulation requires insights into RNA transcription, processing, modification, and translation. However, the relationship between the epitranscriptome and the proteome under drought stress remains undetermined in poplar (Populus trichocarpa). In this study, we used Nanopore direct RNA sequencing and tandem mass tag-based proteomic analysis to examine epitranscriptomic and proteomic regulation induced by drought treatment in stem-differentiating xylem (SDX). Our results revealed a decreased full-length read ratio under drought treatment and, especially, a decreased association between transcriptome and proteome changes in response to drought. Epitranscriptome analysis of cellulose- and lignin-related genes revealed an increased N6-Methyladenosine (m6A) ratio, which was accompanied by decreased RNA abundance and translation, under drought stress. Interestingly, usage of the distal poly(A) site increased during drought stress. Finally, we found that transcripts of highly expressed genes tend to have shorter poly(A) tail length (PAL), and drought stress increased the percentage of transcripts with long PAL. These findings provide insights into the interplay among m6A, polyadenylation, PAL, and translation under drought stress in P. trichocarpa SDX.

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.

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.

Genome-Wide Profiling Reveals Alternative Polyadenylation of Innate Immune-Related mRNA in Patients With COVID-19

The coronavirus disease 2019 (COVID-19) pandemic has caused many deaths worldwide. To date, the mechanism of viral immune escape remains unclear, which is a great obstacle to developing effective clinical treatment. RNA processing mechanisms, including alternative polyadenylation (APA) and alternative splicing (AS), are crucial in the regulation of most human genes in many types of infectious diseases. Because the role of APA and AS in response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection remains unknown, we performed de novo identification of dynamic APA sites using a public dataset of human peripheral blood mononuclear cell (PBMC) RNA-Seq data in COVID-19 patients. We found that genes with APA were enriched in innate immunity -related gene ontology categories such as neutrophil activation, regulation of the MAPK cascade and cytokine production, response to interferon-gamma and the innate immune response. We also reported genome-wide AS events and enriched viral transcription-related categories upon SARS-CoV-2 infection. Interestingly, we found that APA events may give better predictions than AS in COVID-19 patients, suggesting that APA could act as a potential therapeutic target and novel biomarker in those patients. Our study is the first to annotate genes with APA and AS in COVID-19 patients and highlights the roles of APA variation in SARS-CoV-2 infection.

Phase separation in plants: New insights into cellular compartmentalization

A fundamental challenge for cells is how to coordinate various biochemical reactions in space and time. To achieve spatiotemporal control, cells have developed organelles that are surrounded by lipid bilayer membranes. Further, membraneless compartmentalization, a process induced by dynamic physical association of biomolecules through phase transition offers another efficient mechanism for intracellular organization. While our understanding of phase separation was predominantly dependent on yeast and animal models, recent findings have provided compelling evidence for emerging roles of phase separation in plants. In this review, we first provide an overview of the current knowledge of phase separation, including its definition, biophysical principles, molecular features and regulatory mechanisms. Then we summarize plant-specific phase separation phenomena and describe their functions in plant biological processes in great detail. Moreover, we propose that phase separation is an evolutionarily conserved and efficient mechanism for cellular compartmentalization which allows for distinct metabolic processes and signaling pathways, and is especially beneficial for the sessile lifestyle of plants to quickly and efficiently respond to the changing environment.

m6A RNA methylation facilitates pre-mRNA 3'-end formation and is essential for viability of Toxoplasma gondii

Toxoplasma gondii is an obligate intracellular parasite that can cause serious opportunistic disease in the immunocompromised or through congenital infection. To progress through its life cycle, Toxoplasma relies on multiple layers of gene regulation that includes an array of transcription and epigenetic factors. Over the last decade, the modification of mRNA has emerged as another important layer of gene regulation called epitranscriptomics. Here, we report that epitranscriptomics machinery exists in Toxoplasma, namely the methylation of adenosines (m6A) in mRNA transcripts. We identified novel components of the m6A methyltransferase complex and determined the distribution of m6A marks within the parasite transcriptome. m6A mapping revealed the modification to be preferentially located near the 3’-boundary of mRNAs. Knockdown of the m6A writer components METTL3 and WTAP resulted in diminished m6A marks and a complete arrest of parasite replication. Furthermore, we examined the two proteins in Toxoplasma that possess YTH domains, which bind m6A marks, and showed them to be integral members of the cleavage and polyadenylation machinery that catalyzes the 3’-end processing of pre-mRNAs. Loss of METTL3, WTAP, or YTH1 led to a defect in transcript 3’-end formation. Together, these findings establish that the m6A epitranscriptome is essential for parasite viability by contributing to the processing of mRNA 3’-ends.

Toxoplasma gondii is a parasite of medical importance that causes disease upon immuno-suppression. Uncovering essential pathways that the parasite uses for its basic biological processes may reveal opportunities for new anti-parasitic drug therapies. Here, we describe the machinery that Toxoplasma uses to modify specific adenosine residues within its messenger RNAs (mRNA) by N6-adenosine methylation (m6A). We discovered that m6A mRNA methylation is prevalent in multiple stages of the parasite life cycle and is required for parasite replication. We also establish that m6A plays a major role in the proper maturation of mRNA. Two proteins that bind m6A modifications on mRNA associate with factors responsible for the cleavage and final processing steps of mRNA maturation. Since all of the machinery is conserved from plants to Toxoplasma and other related parasites, we propose that this system operates similarly in these organisms.

QuantifyPoly(A): reshaping alternative polyadenylation landscapes of eukaryotes with weighted density peak clustering

The dynamic choice of different polyadenylation sites in a gene is referred to as alternative polyadenylation, which functions in many important biological processes. Large-scale messenger RNA 3' end sequencing has revealed that cleavage sites for polyadenylation are presented with microheterogeneity. To date, the conventional determination of polyadenylation site clusters is subjective and arbitrary, leading to inaccurate annotations. Here, we present a weighted density peak clustering method, QuantifyPoly(A), to accurately quantify genome-wide polyadenylation choices. Applying QuantifyPoly(A) on published 3' end sequencing datasets from both animals and plants, their polyadenylation profiles are reshaped into myriads of novel polyadenylation site clusters. Most of these novel polyadenylation site clusters show significantly dynamic usage across different biological samples or associate with binding sites of trans-acting factors. Upstream sequences of these clusters are enriched with polyadenylation signals UGUA, UAAA and/or AAUAAA in a species-dependent manner. Polyadenylation site clusters also exhibit species specificity, while plants ones generally show higher microheterogeneity than that of animals. QuantifyPoly(A) is broadly applicable to any types of 3' end sequencing data and species for accurate quantification and construction of the complex and dynamic polyadenylation landscape and enables us to decode alternative polyadenylation events invisible to conventional methods at a much higher resolution.

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.

CPSF30-L-mediated recognition of mRNA m6A modification controls alternative polyadenylation of nitrate signaling-related gene transcripts in Arabidopsis

N⁶-methyladenosine (m⁶A), a ubiquitous internal modification of eukaryotic mRNAs, plays a vital role in almost every aspect of mRNA metabolism. However, there is little evidence documenting the role of m⁶A in regulating alternative polyadenylation (APA) in plants. APA is controlled by a large protein-RNA complex with many components, including CLEAVAGE AND POLYADENYLATION SPECIFICITY FACTOR30 (CPSF30). In Arabidopsis, CPSF30 has two isoforms and the longer isoform (CPSF30-L) contains a YT512-B Homology (YTH) domain, which is unique to plants. In this study, we showed that CPSF30-L YTH domain binds to m⁶A in vitro. In the cpsf30-2 mutant, the transcripts of many genes including several important nitrate signaling-related genes had shifts in polyadenylation sites that were correlated with m⁶A peaks, indicating that these gene transcripts carrying m⁶A tend to be regulated by APA. Wild-type CPSF30-L could rescue the defects in APA and nitrate metabolism in cpsf30-2, but m⁶A-binding-defective mutants of CPSF30-L could not. Taken together, our results demonstrated that m⁶A modification regulates APA in Arabidopsis and revealed that the m⁶A reader CPSF30-L affects nitrate signaling by controlling APA, shedding new light on the roles of the m⁶A modification during RNA 3'-end processing in nitrate metabolism.

Co-Transcriptional RNA Processing in Plants: Exploring from the Perspective of Polyadenylation

Most protein-coding genes in eukaryotes possess at least two poly(A) sites, and alternative polyadenylation is considered a contributing factor to transcriptomic and proteomic diversity. Following transcription, a nascent RNA usually undergoes capping, splicing, cleavage, and polyadenylation, resulting in a mature messenger RNA (mRNA); however, increasing evidence suggests that transcription and RNA processing are coupled. Plants, which must produce rapid responses to environmental changes because of their limited mobility, exhibit such coupling. In this review, we summarize recent advances in our understanding of the coupling of transcription with RNA processing in plants, and we describe the possible spatial environment and important proteins involved. Moreover, we describe how liquid–liquid phase separation, mediated by the C-terminal domain of RNA polymerase II and RNA processing factors with intrinsically disordered regions, enables efficient co-transcriptional mRNA processing in plants.

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.

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.

Marking RNA: m6A writers, readers, and functions in Arabidopsis

N6-methyladenosine (m6A) emerges as an important modification in eukaryotic mRNAs. m6A has first been reported in 1974, and its functional significance in mammalian gene regulation and importance for proper development have been well established. An arsenal of writer, eraser, and reader proteins accomplish deposition, removal, and interpretation of the m6A mark, resulting in dynamic function. This led to the concept of an epitranscriptome, the compendium of RNA species with chemical modification of the nucleobases in the cell, in analogy to the epigenome. While m6A has long been known to also exist in plant mRNAs, proteins involved in m6A metabolism have only recently been detected by mutant analysis, homology search, and mRNA interactome capture in the reference plant Arabidopsis thaliana. Dysregulation of the m6A modification causes severe developmental abnormalities of leaves and roots and altered timing of reproductive development. Furthermore, m6A modification affects viral infection. Here, we discuss recent progress in identifying m6A sites transcriptome-wide, in identifying the molecular players involved in writing, removing, and reading the mark, and in assigning functions to this RNA modification in A. thaliana. We highlight similarities and differences to m6A modification in mammals and provide an outlook on important questions that remain to be addressed.

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.

Occurrence and Functions of m6A and Other Covalent Modifications in Plant mRNA

Posttranscriptional control of gene expression is indispensable for the execution of developmental programs and environmental adaptation. Among the many cellular mechanisms that regulate mRNA fate, covalent nucleotide modification has emerged as a major way of controlling the processing, localization, stability, and translatability of mRNAs. This powerful mechanism is conserved across eukaryotes and controls the cellular events that lead to development and growth. As in other eukaryotes, N 6-methylation of adenosine is the most abundant and best studied mRNA modification in flowering plants. It is essential for embryonic and postembryonic plant development and it affects growth rate and stress responses, including susceptibility to plant RNA viruses. Although the mRNA modification field is young, the intense interest triggered by its involvement in stem cell differentiation and cancer has led to rapid advances in understanding how mRNA modifications control gene expression in mammalian systems. An equivalent effort from plant molecular biologists has been lagging behind, but recent work in Arabidopsis (Arabidopsis thaliana) and other plant species is starting to give insights into how this essential layer of posttranscriptional regulation works in plants, and both similarities and differences with other eukaryotes are emerging. In this Update, we summarize, connect, and evaluate the experimental work that supports our current knowledge of the biochemistry, molecular mechanisms, and biological functions of mRNA modifications in plants. We devote particular attention to N 6-methylation of adenosine and attempt to place the knowledge gained from plant studies within the context of a more general framework derived from studies in other eukaryotes.

Genome-wide alternative polyadenylation dynamics in response to biotic and abiotic stresses in rice

Alternative polyadenylation (APA) is an important way to regulate gene expression at the post-transcriptional level, and is extensively involved in plant stress responses. However, the systematic roles of APA regulation in response to abiotic and biotic stresses in rice at the genome scale remain unknown. To take advantage of available RNA-seq datasets, using a novel tool APAtrap, we identified thousands of genes with significantly differential usage of polyadenylation [poly(A)] sites in response to the abiotic stress (drought, heat shock, and cadmium) and biotic stress [bacterial blight (BB), rice blast, and rice stripe virus (RSV)]. Genes with stress-responsive APA dynamics commonly exhibited higher expression levels when their isoforms with short 3' untranslated region (3' UTR) were more abundant. The stress-responsive APA events were widely involved in crucial stress-responsive genes and pathways: e.g. APA acted as a negative regulator in heat stress tolerance; APA events were involved in DNA repair and cell wall formation under Cd stress; APA regulated chlorophyll metabolism, being associated with the pathogenesis of leaf diseases under RSV and BB challenges. Furthermore, APA events were found to be involved in glutathione metabolism and MAPK signaling pathways, mediating a crosstalk among the abiotic and biotic stress-responsive regulatory networks in rice. Analysis of large-scale datasets revealed that APA may regulate abiotic and biotic stress-responsive processes in rice. Such post-transcriptome diversities contribute to rice adaption to various environmental challenges. Our study would supply useful resource for further molecular assisted breeding of multiple stress-tolerant cultivars for rice.

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

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

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.

The m6A pathway protects the transcriptome integrity by restricting RNA chimera formation in plants

This study reveals that an m⁶A-assisted polyadenylation pathway comprising conserved m⁶A writer proteins and a plant-specific m⁶A reader contributes to transcriptome integrity in Arabidopsis thaliana by restricting RNA chimera formation at rearranged loci.

Global, segmental, and gene duplication–related processes are driving genome size and complexity in plants. Despite their evolutionary potentials, those processes can also have adverse effects on genome regulation, thus implying the existence of specialized corrective mechanisms. Here, we report that an N6-methyladenosine (m⁶A)–assisted polyadenylation (m-ASP) pathway ensures transcriptome integrity in Arabidopsis thaliana. Efficient m-ASP pathway activity requires the m⁶A methyltransferase-associated factor FIP37 and CPSF30L, an m⁶A reader corresponding to an YT512-B Homology Domain-containing protein (YTHDC)-type domain containing isoform of the 30-kD subunit of cleavage and polyadenylation specificity factor. Targets of the m-ASP pathway are enriched in recently rearranged gene pairs, displayed an atypical chromatin signature, and showed transcriptional readthrough and mRNA chimera formation in FIP37- and CPSF30L-deficient plants. Furthermore, we showed that the m-ASP pathway can also restrict the formation of chimeric gene/transposable-element transcript, suggesting a possible implication of this pathway in the control of transposable elements at specific locus. Taken together, our results point to selective recognition of 3′-UTR m⁶A as a safeguard mechanism ensuring transcriptome integrity at rearranged genomic loci in plants.

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.

Cluster analysis of replicated alternative polyadenylation data using canonical correlation analysis

BACKGROUND

Alternative polyadenylation (APA) has emerged as a pervasive mechanism that contributes to the transcriptome complexity and dynamics of gene regulation. The current tsunami of whole genome poly(A) site data from various conditions generated by 3' end sequencing provides a valuable data source for the study of APA-related gene expression. Cluster analysis is a powerful technique for investigating the association structure among genes, however, conventional gene clustering methods are not suitable for APA-related data as they fail to consider the information of poly(A) sites (e.g., location, abundance, number, etc.) within each gene or measure the association among poly(A) sites between two genes.

RESULTS

Here we proposed a computational framework, named PASCCA, for clustering genes from replicated or unreplicated poly(A) site data using canonical correlation analysis (CCA). PASCCA incorporates multiple layers of gene expression data from both the poly(A) site level and gene level and takes into account the number of replicates and the variability within each experimental group. Moreover, PASCCA characterizes poly(A) sites in various ways including the abundance and relative usage, which can exploit the advantages of 3' end deep sequencing in quantifying APA sites. Using both real and synthetic poly(A) site data sets, the cluster analysis demonstrates that PASCCA outperforms other widely-used distance measures under five performance metrics including connectivity, the Dunn index, average distance, average distance between means, and the biological homogeneity index. We also used PASCCA to infer APA-specific gene modules from recently published poly(A) site data of rice and discovered some distinct functional gene modules. We have made PASCCA an easy-to-use R package for APA-related gene expression analyses, including the characterization of poly(A) sites, quantification of association between genes, and clustering of genes.

CONCLUSIONS

By providing a better treatment of the noise inherent in repeated measurements and taking into account multiple layers of poly(A) site data, PASCCA could be a general tool for clustering and analyzing APA-specific gene expression data. PASCCA could be used to elucidate the dynamic interplay of genes and their APA sites among various biological conditions from emerging 3' end sequencing data to address the complex biological phenomenon.

Readers of the m6A epitranscriptomic code

N⁶-methyl adenosine (m⁶A) is the most prevalent and evolutionarily conserved, modification of polymerase II transcribed RNAs. By post-transcriptionally controlling patterns of gene expression, m⁶A deposition is crucial for organism reproduction, development and likely stress responses. m⁶A mostly mediates its effect by recruiting reader proteins that either directly accommodate the modified residue in a hydrophobic pocket formed by their YTH domain, or otherwise have their affinity positively influenced by the presence of m⁶A. We firstly describe here the evolutionary history, and review known molecular and physiological roles of eukaryote YTH readers. In the second part, we present non YTH-proteins whose roles as m⁶A readers largely remain to be explored. The diversity and multiplicity of m⁶A readers together with the possibility to regulate their expression and function in response to various cues, offers a multitude of possible combinations to rapidly and finely tune gene expression patterns and hence cellular plasticity. This article is part of a Special Issue entitled: mRNA modifications in gene expression control edited by Dr. Soller Matthias and Dr. Fray Rupert.

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.