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

Untranslated yet indispensable-UTRs act as key regulators in the environmental control of gene expression

To survive and thrive in a dynamic environment, plants must continuously monitor their surroundings and adjust their development and physiology accordingly. Changes in gene expression underlie these developmental and physiological adjustments, and are traditionally attributed to widespread transcriptional reprogramming. Growing evidence, however, suggests that post-transcriptional mechanisms also play a vital role in tailoring gene expression to a plant's environment. Untranslated regions (UTRs) act as regulatory hubs for post-transcriptional control, harbouring cis-elements that affect an mRNA's processing, localization, translation, and stability, and thereby tune the abundance of the encoded protein. Here, we review recent advances made in understanding the critical function UTRs exert in the post-transcriptional control of gene expression in the context of a plant's abiotic environment. We summarize the molecular mechanisms at play, present examples of UTR-controlled signalling cascades, and discuss the potential that resides within UTRs to render plants more resilient to a changing climate.

Identification of Eleutherococcus senticosus NAC transcription factors and their mechanisms in mediating DNA methylation of EsFPS, EsSS, and EsSE promoters to regulate saponin synthesis

BACKGROUND

The formation of pharmacologically active components in medicinal plants is significantly impacted by DNA methylation. However, the exact mechanisms through which DNA methylation regulates secondary metabolism remain incompletely understood. Research in model species has demonstrated that DNA methylation at the transcription factor binding site within functional gene promoters can impact the binding of transcription factors to target DNA, subsequently influencing gene expression. These findings suggest that the interaction between transcription factors and target DNA could be a significant mechanism through which DNA methylation regulates secondary metabolism in medicinal plants.

RESULTS

This research conducted a comprehensive analysis of the NAC family in E. senticosus, encompassing genome-wide characterization and functional analysis. A total of 117 EsNAC genes were identified and phylogenetically divided into 15 subfamilies. Tandem duplications and chromosome segment duplications were found to be the primary replication modes of these genes. Motif 2 was identified as the core conserved motif of the genes, and the cis-acting elements, gene structures, and expression patterns of each EsNAC gene were different. EsJUB1, EsNAC047, EsNAC098, and EsNAC005 were significantly associated with the DNA methylation ratio in E. senticosus. These four genes were located in the nucleus or cytoplasm and exhibited transcriptional self-activation activity. DNA methylation in EsFPS, EsSS, and EsSE promoters significantly reduced their activity. The methyl groups added to cytosine directly hindered the binding of the promoters to EsJUB1, EsNAC047, EsNAC098, and EsNAC005 and altered the expression of EsFPS, EsSS, and EsSE genes, eventually leading to changes in saponin synthesis in E. senticosus.

CONCLUSIONS

NAC transcription factors that are hindered from binding by methylated DNA are found in E. senticosus. The incapacity of these NACs to bind to the promoter of the methylated saponin synthase gene leads to subsequent alterations in gene expression and saponin synthesis. This research is the initial evidence showcasing the involvement of EsNAC in governing the impact of DNA methylation on saponin production in E. senticosus.

Epitranscriptome insights into Riccia fluitans L. (Marchantiophyta) aquatic transition using nanopore direct RNA sequencing

BACKGROUND

Riccia fluitans, an amphibious liverwort, exhibits a fascinating adaptation mechanism to transition between terrestrial and aquatic environments. Utilizing nanopore direct RNA sequencing, we try to capture the complex epitranscriptomic changes undergone in response to land-water transition.

RESULTS

A significant finding is the identification of 45 differentially expressed genes (DEGs), with a split of 33 downregulated in terrestrial forms and 12 upregulated in aquatic forms, indicating a robust transcriptional response to environmental changes. Analysis of N6-methyladenosine (m6A) modifications revealed 173 m6A sites in aquatic and only 27 sites in the terrestrial forms, indicating a significant increase in methylation in the former, which could facilitate rapid adaptation to changing environments. The aquatic form showed a global elongation bias in poly(A) tails, which is associated with increased mRNA stability and efficient translation, enhancing the plant's resilience to water stress. Significant differences in polyadenylation signals were observed between the two forms, with nine transcripts showing notable changes in tail length, suggesting an adaptive mechanism to modulate mRNA stability and translational efficiency in response to environmental conditions. This differential methylation and polyadenylation underline a sophisticated layer of post-transcriptional regulation, enabling Riccia fluitans to fine-tune gene expression in response to its living conditions.

CONCLUSIONS

These insights into transcriptome dynamics offer a deeper understanding of plant adaptation strategies at the molecular level, contributing to the broader knowledge of plant biology and evolution. These findings underscore the sophisticated post-transcriptional regulatory strategies Riccia fluitans employs to navigate the challenges of aquatic versus terrestrial living, highlighting the plant's dynamic adaptation to environmental stresses and its utility as a model for studying adaptation mechanisms in amphibious plants.

The cAMP-dependent phosphorylation footprint in response to heat stress

KEY MESSAGE

cAMP modulates the phosphorylation status of highly conserved phosphosites in RNA-binding proteins crucial for mRNA metabolism and reprogramming in response to heat stress. In plants, 3',5'-cyclic adenosine monophosphate (3',5'-cAMP) is a second messenger that modulates multiple cellular targets, thereby participating in plant developmental and adaptive processes. Although its role in ameliorating heat-related damage has been demonstrated, mechanisms that govern cAMP-dependent responses to heat have remained elusive. Here we analyze the role cAMP-dependent phosphorylation during prolonged heat stress (HS) with a view to gain insight into processes that govern plant responses to HS. To do so, we performed quantitative phosphoproteomic analyses in Nicotiana tabacum Bright Yellow-2 cells grown at 27 °C or 35 °C for 3 days overexpressing a molecular "sponge" that reduces free intracellular cAMP levels. Our phosphorylation data and analyses reveal that the presence of cAMP is an essential factor that governs specific protein phosphorylation events that occur during prolonged HS in BY-2 cells. Notably, cAMP modulates HS-dependent phosphorylation of proteins that functions in mRNA processing, transcriptional control, vesicular trafficking, and cell cycle regulation and this is indicative for a systemic role of the messenger. In particular, changes of cAMP levels affect the phosphorylation status of highly conserved phosphosites in 19 RNA-binding proteins that are crucial during the reprogramming of the mRNA metabolism in response to HS. Furthermore, phosphorylation site motifs and molecular docking suggest that some proteins, including kinases and phosphatases, are conceivably able to directly interact with cAMP thus further supporting a regulatory role of cAMP in plant HS responses.

Determination of Differential Alternative Splicing Under Stress Conditions

Alternative splicing (AS) is an important mechanism contributing to stress-induced regulation of gene expression and proteome diversity. Massive sequencing technologies allow the identification of transcripts generated via stress-responsive AS, potentially important for adaptation to stress conditions. Several bioinformatics tools have been developed to identify differentially expressed alternative splicing events/transcripts from RNA-sequencing results. This chapter describes a detailed protocol for differential alternative splicing analysis using the rMATS tool. In addition, we provide guidelines for validation of the detected splice variants by qRT-PCR based on the obtained output files.

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.

RNA out of the mist

RNA has long been regarded primarily as the intermediate between genes and proteins. It was a surprise then to discover that eukaryotic genes are mosaics of mRNA sequences interrupted by large tracts of transcribed but untranslated sequences, and that multicellular organisms also express many long 'intergenic' and antisense noncoding RNAs (lncRNAs). The identification of small RNAs that regulate mRNA translation and half-life did not disturb the prevailing view that animals and plant genomes are full of evolutionary debris and that their development is mainly supervised by transcription factors. Gathering evidence to the contrary involved addressing the low conservation, expression, and genetic visibility of lncRNAs, demonstrating their cell-specific roles in cell and developmental biology, and their association with chromatin-modifying complexes and phase-separated domains. The emerging picture is that most lncRNAs are the products of genetic loci termed 'enhancers', which marshal generic effector proteins to their sites of action to control cell fate decisions during development.

From common gardens to candidate genes: exploring local adaptation to climate in red spruce

Local adaptation to climate is common in plant species and has been studied in a range of contexts, from improving crop yields to predicting population maladaptation to future conditions. The genomic era has brought new tools to study this process, which was historically explored through common garden experiments. In this study, we combine genomic methods and common gardens to investigate local adaptation in red spruce and identify environmental gradients and loci involved in climate adaptation. We first use climate transfer functions to estimate the impact of climate change on seedling performance in three common gardens. We then explore the use of multivariate gene-environment association methods to identify genes underlying climate adaptation, with particular attention to the implications of conducting genome scans with and without correction for neutral population structure. This integrative approach uncovered phenotypic evidence of local adaptation to climate and identified a set of putatively adaptive genes, some of which are involved in three main adaptive pathways found in other temperate and boreal coniferous species: drought tolerance, cold hardiness, and phenology. These putatively adaptive genes segregated into two 'modules' associated with different environmental gradients. This study nicely exemplifies the multivariate dimension of adaptation to climate in trees.

Alternative Polyadenylation Is a Novel Strategy for the Regulation of Gene Expression in Response to Stresses in Plants

Precursor message RNA requires processing to generate mature RNA. Cleavage and polyadenylation at the 3'-end in the maturation of mRNA is one of key processing steps in eukaryotes. The polyadenylation (poly(A)) tail of mRNA is an essential feature that is required to mediate its nuclear export, stability, translation efficiency, and subcellular localization. Most genes have at least two mRNA isoforms via alternative splicing (AS) or alternative polyadenylation (APA), which increases the diversity of transcriptome and proteome. However, most previous studies have focused on the role of alternative splicing on the regulation of gene expression. In this review, we summarize the recent advances concerning APA in the regulation of gene expression and in response to stresses in plants. We also discuss the mechanisms for the regulation of APA for plants in the adaptation to stress responses, and suggest that APA is a novel strategy for the adaptation to environmental changes and response to stresses in plants.

Genome-wide identification of MBD gene family members in Eleutherococcus senticosus with their expression motifs under drought stress and DNA demethylation

BACKGROUND

Methyl-binding domain (MBD) is a class of methyl-CpG-binding domain proteins that affects the regulation of gene expression through epigenetic modifications. MBD genes are not only inseparable from DNA methylation but have also been identified and validated in various plants. Although MBD is involved in a group of physiological processes and stress regulation in these plants, MBD genes in Eleutherococcus senticosus remain largely unknown.

RESULTS

Twenty EsMBD genes were identified in E. senticosus. Among the 24 chromosomes of E. senticosus, EsMBD genes were unevenly distributed on 12 chromosomes, and only one tandem repeat gene existed. Collinearity analysis showed that the fragment duplication was the main motif for EsMBD gene expansion. As the species of Araliaceae evolved, MBD genes also evolved and gradually exhibited different functional differentiation. Furthermore, cis-acting element analysis showed that there were numerous cis-acting elements in the EsMBD promoter region, among which light response elements and anaerobic induction elements were dominant. The expression motif analysis revealed that 60% of the EsMBDs were up-regulated in the 30% water content group.

CONCLUSIONS

By comparing the transcriptome data of different saponin contents of E. senticosus and integrating them with the outcomes of molecular docking analysis, we hypothesized that EsMBD2 and EsMBD5 jointly affect the secondary metabolic processes of E. senticosus saponins by binding to methylated CpG under conditions of drought stress. The results of this study laid the foundation for subsequent research on the E. senticosus and MBD genes.

Investigating nicotine pathway-related long non-coding RNAs in tobacco

Long non-coding RNAs (lncRNAs) are transcripts longer than 200 bp with low or no protein-coding ability, which play essential roles in various biological processes in plants. Tobacco is an ideal model plant for studying nicotine biosynthesis and metabolism, and there is little research on lncRNAs in this field. Therefore, how to take advantage of the mature tobacco system to profoundly investigate the lncRNAs involved in the nicotine pathway is intriguing. By exploiting 549 public RNA-Seq datasets of tobacco, 30,212 lncRNA candidates were identified, including 24,084 large intervening non-coding RNAs (lincRNAs), 5,778 natural antisense transcripts (NATs) and 350 intronic non-coding RNAs (incRNAs). Compared with protein-coding genes, lncRNAs have distinct properties in terms of exon number, sequence length, A/U content, and tissue-specific expression pattern. lincRNAs showed an asymmetric evolutionary pattern, with a higher proportion (68.71%) expressed from the Nicotiana sylvestris (S) subgenome. We predicted the potential cis/trans-regulatory effects on protein-coding genes. One hundred four lncRNAs were detected as precursors of 30 known microRNA (miRNA) family members, and 110 lncRNAs were expected to be the potential endogenous target mimics for 39 miRNAs. By combining the results of weighted gene co-expression network analysis with the differentially expressed gene analysis of topping RNA-seq data, we constructed a sub-network containing eight lncRNAs and 25 nicotine-related coding genes. We confirmed that the expression of seven lncRNAs could be affected by MeJA treatment and may be controlled by the transcription factor NtMYC2 using a quantitative PCR assay and gene editing. The results suggested that lncRNAs are involved in the nicotine pathway. Our findings further deepened the understanding of the features and functions of lncRNAs and provided new candidates for regulating nicotine biosynthesis in tobacco.

Identification of a Saltol-Independent Salinity Tolerance Polymorphism in Rice Mekong Delta Landraces and Characterization of a Promising Line, Doc Phung

The Mekong Delta River in Vietnam is facing salinity intrusion caused by climate change and sea-level rise that is severely affecting rice cultivation. Here, we evaluated salinity responses of 97 rice accessions (79 landraces and 18 improved accessions) from the Mekong Delta population by adding 100 mM NaCl to the nutrient solution for up to 20 days. We observed a wide distribution in salinity tolerance/sensitivity, with two major peaks across the 97 accessions when using the standard evaluation system (SES) developed by the International Rice Research Institute. SES scores revealed strong negative correlations (ranging from - 0.68 to - 0.83) with other phenotypic indices, such as shoot elongation length, root elongation length, shoot dry weight, and root dry weight. Mineral concentrations of Na⁺ in roots, stems, and leaves and Ca²⁺ in roots and stems were positively correlated with SES scores, suggesting that tolerant accessions lower their cation exchange capacity in the root cell wall. The salinity tolerance of Mekong Delta accessions was independent from the previously described salinity tolerance-related locus Saltol, which encodes an HKT1-type transporter in the salinity-tolerant cultivars Nona Bokra and Pokkali. Indeed, genome-wide association studies using SES scores and shoot dry weight ratios of the 79 accessions as traits identified a single common peak located on chromosome 1. This SNP did not form a linkage group with other nearby SNPs and mapped to the 3' untranslated region of gene LOC_Os01g32830, over 6.5 Mb away from the Saltol locus. LOC_Os01g32830 encodes chloroplast glycolate/glycerate translocator 1 (OsPLGG1), which is responsible for photorespiration and growth. SES and shoot dry weight ratios differed significantly between the two possible haplotypes at the causal SNP. Through these analyses, we characterize Doc Phung, one of the most salinity-tolerant varieties in the Mekong Delta population and a promising new genetic resource.

Chicken GLUT4 undergoes complex alternative splicing events and its expression in striated muscle changes dramatically during development

Glucose transporter protein 4 (GLUT4) plays an important role in regulating insulin-mediated glucose homeostasis in mammals. Until now, studies on GLUT4 have focused on mammals mostly, while chicken GLUT4 has been rarely investigated. In this study, chicken GLUT4 mRNA sequences were obtained by combining conventional amplification, 5'- and 3'- rapid amplification of cDNA ends technique (RACE), then bioinformatics analysis on its genomic structure, splicing pattern, subcellular localization prediction and homologous comparisons were carried out. In addition, the distribution of GLUT4 was detected by RT-qPCR in bird's liver and striated muscles (cardiac muscle, pectoralis and leg muscle) at different ages, including embryonic day 14 (E14), E19, 7-day-old (D7), D21 and D49 (n = 3-4). Results showed that chicken GLUT4 gene produced at least 14 transcripts (GenBank accession No: OP491293-OP491306) through alternative splicing and polyadenylation, which predicted encoding 12 types of amino acid (AA) sequences (with length ranged from 65 AA to 519 AA). These proteins contain typical major facilitator superfamily domain of glucose transporters with length variations, sharing a common sequence of 59 AA, and were predicted to have distinct subcellular localization. The dominant transcript (named as T1) consists of 11 exons with an open reading frame being predicted encoding 519 AA. In addition, analyzing on the spatio-temporal expression of chicken GLUT4 showed it dominantly expressed in pectoralis, leg muscles and cardiac muscle, and the mRNA level of chicken GLUT4 dramatically fluctuated with birds' development in cardiac muscle, pectoralis and leg muscles, with the level at D21 significantly higher than that at E14, E19, and D49 (P < 0.05). These data indicated that chicken GLUT4 undergoes complex alternative splicing events, and GLUT4 expression level in striated muscle was subjected to dynamic regulation with birds' development. Results indicate these isoforms may play overlapping and distinct roles in chicken.

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.

Analysis of mRNA-derived siRNAs in mutants of mRNA maturation and surveillance pathways in Arabidopsis thaliana

Defects in RNA maturation and RNA decay factors may generate substrates for the RNA interference machinery. This phenomenon was observed in plants where mutations in some RNA-related factors lead to the production of RNA-quality control small interfering RNAs and several mutants show enhanced silencing of reporter transgenes. To assess the potential of RNAi activation on endogenous transcripts, we sequenced small RNAs from a set of Arabidopsis thaliana mutants with defects in various RNA metabolism pathways. We observed a global production of siRNAs caused by inefficient pre-mRNA cleavage and polyadenylation leading to read-through transcription into downstream antisense genes. In addition, in the lsm1a lsm1b double mutant, we identified NIA1, SMXL5, and several miRNA-targeted mRNAs as producing siRNAs, a group of transcripts suggested being especially sensitive to deficiencies in RNA metabolism. However, in most cases, RNA metabolism perturbations do not lead to the widespread production of siRNA derived from mRNA molecules. This observation is contrary to multiple studies based on reporter transgenes and suggests that only a very high accumulation of defective mRNA species caused by specific mutations or substantial RNA processing defects trigger RNAi pathways.

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.

Immediate transcriptional responses of Arabidopsis leaves to heat shock

Plants have evolved efficient mechanisms for adapting to temperature fluctuations, known as heat stress response and heat stress memory. Although the transcriptional regulatory network of plant heat stress response has been established, little is known about the genome-wide transcriptional changes occurring within the first several minutes after heat shock. Here, we investigated the nascent RNA and mature messenger RNA (mRNA) from plant leaf tissues exposed to 5 min of heat shock treatment using global run-on sequencing and RNA sequencing methods. Only a small group of genes were up- or downregulated at both the nascent RNA and mRNA levels. Primed plants that were already exposed to mild heat stress exhibited a more drastic alteration at multiple transcriptional steps than naïve plants that had not experienced heat stress. Upon heat shock, we also observed the following: (i) engaged RNA polymerase II accumulated downstream of transcription start sites; (ii) 5' pausing release was a rate-limiting step for the induction of some heat shock protein genes; (iii) numerous genes switched transcription modes; (iv) pervasive read-through was induced at terminators; and (v) heat stress memory occurs at multiple steps of the transcription cycle, such as at Pol II recruitment, 5' pausing, elongation, and termination.

PLncDB V2.0: a comprehensive encyclopedia of plant long noncoding RNAs

Long noncoding RNAs (lncRNAs) are transcripts longer than 200 nucleotides with little or no protein coding potential. The expanding list of lncRNAs and accumulating evidence of their functions in plants have necessitated the creation of a comprehensive database for lncRNA research. However, currently available plant lncRNA databases have some deficiencies, including the lack of lncRNA data from some model plants, uneven annotation standards, a lack of visualization for expression patterns, and the absence of epigenetic information. To overcome these problems, we upgraded our Plant Long noncoding RNA Database (PLncDB, http://plncdb.tobaccodb.org/), which was based on a uniform annotation pipeline. PLncDB V2.0 currently contains 1 246 372 lncRNAs for 80 plant species based on 13 834 RNA-Seq datasets, integrating lncRNA information from four other resources including EVLncRNAs, RNAcentral and etc. Expression patterns and epigenetic signals can be visualized using multiple tools (JBrowse, eFP Browser and EPexplorer). Targets and regulatory networks for lncRNAs are also provided for function exploration. In addition, PLncDB V2.0 is hierarchical and user-friendly and has five built-in search engines. We believe PLncDB V2.0 is useful for the plant lncRNA community and data mining studies and provides a comprehensive resource for data-driven lncRNA research 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.

Arabidopsis SMN2/HEN2, Encoding DEAD-Box RNA Helicase, Governs Proper Expression of the Resistance Gene SMN1/RPS6 and Is Involved in Dwarf, Autoimmune Phenotypes of mekk1 and mpk4 Mutants

In Arabidopsis thaliana, a mitogen-activated protein kinase pathway, MEKK1-MKK1/MKK2-MPK4, is important for basal resistance and disruption of this pathway results in dwarf, autoimmune phenotypes. To elucidate the complex mechanisms activated by the disruption of this pathway, we have previously developed a mutant screening system based on a dwarf autoimmune line that overexpressed the N-terminal regulatory domain of MEKK1. Here, we report that the second group of mutants, smn2, had defects in the SMN2 gene, encoding a DEAD-box RNA helicase. SMN2 is identical to HEN2, whose function is vital for the nuclear RNA exosome because it provides non-ribosomal RNA specificity for RNA turnover, RNA quality control and RNA processing. Aberrant SMN1/RPS6 transcripts were detected in smn2 and hen2 mutants. Disease resistance against Pseudomonas syringae pv. tomato DC3000 (hopA1), which is conferred by SMN1/RPS6, was decreased in smn2 mutants, suggesting a functional connection between SMN1/RPS6 and SMN2/HEN2. We produced double mutants mekk1smn2 and mpk4smn2 to determine whether the smn2 mutations suppress the dwarf, autoimmune phenotypes of the mekk1 and mpk4 mutants, as the smn1 mutations do. As expected, the mekk1 and mpk4 phenotypes were suppressed by the smn2 mutations. These results suggested that SMN2 is involved in the proper function of SMN1/RPS6. The Gene Ontology enrichment analysis using RNA-seq data showed that defense genes were downregulated in smn2, suggesting a positive contribution of SMN2 to the genome-wide expression of defense genes. In conclusion, this study provides novel insight into plant immunity via SMN2/HEN2, an essential component of the nuclear RNA exosome.

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.

The transcription and export complex THO/TREX contributes to transcription termination in plants

Transcription termination has important regulatory functions, impacting mRNA stability, localization and translation potential. Failure to appropriately terminate transcription can also lead to read-through transcription and the synthesis of antisense RNAs which can have profound impact on gene expression. The Transcription-Export (THO/TREX) protein complex plays an important role in coupling transcription with splicing and export of mRNA. However, little is known about the role of the THO/TREX complex in the control of transcription termination. In this work, we show that two proteins of the THO/TREX complex, namely TREX COMPONENT 1 (TEX1 or THO3) and HYPER RECOMBINATION1 (HPR1 or THO1) contribute to the correct transcription termination at several loci in Arabidopsis thaliana. We first demonstrate this by showing defective termination in tex1 and hpr1 mutants at the nopaline synthase (NOS) terminator present in a T-DNA inserted between exon 1 and 3 of the PHO1 locus in the pho1-7 mutant. Read-through transcription beyond the NOS terminator and splicing-out of the T-DNA resulted in the generation of a near full-length PHO1 mRNA (minus exon 2) in the tex1 pho1-7 and hpr1 pho1-7 double mutants, with enhanced production of a truncated PHO1 protein that retained phosphate export activity. Consequently, the strong reduction of shoot growth associated with the severe phosphate deficiency of the pho1-7 mutant was alleviated in the tex1 pho1-7 and hpr1 pho1-7 double mutants. Additionally, we show that RNA termination defects in tex1 and hpr1 mutants leads to 3’UTR extensions in several endogenous genes. These results demonstrate that THO/TREX complex contributes to the regulation of transcription termination.

Production of messenger RNAs (mRNAs) involves numerous steps including initiation of transcription, elongation, splicing, termination, as well as export out of the nucleus. All these steps are highly coordinated and failure in any steps has a profound impact on the level and identity of mRNAs produced. The THO/TREX protein complex is associated with nascent RNAs and contributes to several mRNA biogenesis steps, including splicing and export. However, the contribution of the THO/TREX complex to mRNA termination was poorly defined. We have identified a role for two components of the THO/TREX complex, namely the proteins TEX1 and HPR1, in the control of transcription termination in the plant Arabidopsis thaliana. We show that the tex1 and hpr1 mutants have defects in terminating mRNA at the nopaline synthase (NOS) terminator found in a T-DNA insertion mutant leading to the transcriptional read-through pass the NOS terminator. We also show that tex1 and hpr1 mutants have defects in mRNA termination at several endogenous genes, leading to the production of 3’UTR extensions. Together, these results highlight a role for the THO/TREX complex in mRNA termination.

Metabolic Labeling of RNAs Uncovers Hidden Features and Dynamics of the Arabidopsis Transcriptome

Transcriptome analysis by RNA sequencing (RNA-seq) has become an indispensable research tool in modern plant biology. Virtually all RNA-seq studies provide a snapshot of the steady state transcriptome, which contains valuable information about RNA populations at a given time but lacks information about the dynamics of RNA synthesis and degradation. Only a few specialized sequencing techniques, such as global run-on sequencing, have been used to provide information about RNA synthesis rates in plants. Here, we demonstrate that RNA labeling with the modified, nontoxic uridine analog 5-ethynyl uridine (5-EU) in Arabidopsis (Arabidopsis thaliana) seedlings provides insight into plant transcriptome dynamics. Pulse labeling with 5-EU revealed nascent and unstable RNAs, RNA processing intermediates generated by splicing, and chloroplast RNAs. Pulse-chase experiments with 5-EU allowed us to determine RNA stabilities without the need for chemical transcription inhibitors such as actinomycin and cordycepin. Inhibitor-free, genome-wide analysis of polyadenylated RNA stability via 5-EU pulse-chase experiments revealed RNAs with shorter half-lives than those reported after chemical inhibition of transcription. In summary, our results indicate that the Arabidopsis nascent transcriptome contains unstable RNAs and RNA processing intermediates and suggest that polyadenylated RNAs have low stability in plants. Our technique lays the foundation for easy, affordable, nascent transcriptome analysis and inhibitor-free analysis of RNA stability in plants.

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.

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

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

UTR-Dependent Control of Gene Expression in Plants

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

Defective XRN3-mediated transcription termination in Arabidopsis affects the expression of protein-coding genes

Arabidopsis thaliana contains two nuclear XRN2/3 5'-3' exonucleases that are homologs of yeast and human Rat1/Xrn2 proteins involved in the processing and degradation of several classes of nuclear RNAs and in transcription termination of RNA polymerase II. Using strand-specific short read sequencing we show that knockdown of XRN3 leads to an altered expression of hundreds of genes and the accumulation of uncapped and polyadenylated read-through transcripts generated by inefficiently terminated Pol II. Our data support the notion that XRN3-mediated changes in the expression of a subset of genes are caused by upstream read-through transcription and these effects are enhanced by RNA-mRNA chimeras generated in xrn3 plants. In turn, read-through transcripts that are antisense to downstream genes may trigger production of siRNA. Our results highlight the importance of XRN3 exoribonuclease in Pol II transcription termination in plants and show that disturbance in this process may significantly alter gene expression.