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

N-acetyltransferase 10 catalyzes RNA N4-acetylcytidine to regulate photosynthesis in rice

N4-acetylcytidine (ac4C) is a novel mRNA modification that enhances RNA stability and translation in mammals and plants. We previously identified ac4C in Arabidopsis, introduced by two homologs of human N-acetyltransferase 10 (NAT10). While ac4C influences leaf development in Arabidopsis, its role in rice is unclear. In this study, we identify OsNAT10 as the ac4C writer in rice. osnat10 mutants show developmental defects, including shorter roots, fewer tillers, and lower yield. Compared with wild type, ac4C-modified genes are less abundant in osnat10, particularly those related to photosynthesis. Additionally, osnat10 exhibits decreased photosynthetic capacity and reduced RNA stability and translation efficiency of ac4C target genes, like LIGHT-INDUCED RICE 1 (LIR1). Overexpressing OsLIR1 partially rescues osnat10 defects, underscoring OsNAT10's role in photosynthesis regulation. Our findings highlight ac4C's crucial function in photosynthesis and plant development, offering insights into epitranscriptomic modifications for crop improvement.

Quantitative RNA pseudouridine maps reveal multilayered translation control through plant rRNA, tRNA and mRNA pseudouridylation

Pseudouridine (Ψ) is the most abundant RNA modification, yet studies of Ψ have been hindered by a lack of robust methods to profile comprehensive Ψ maps. Here we utilize bisulfite-induced deletion sequencing to generate transcriptome-wide Ψ maps at single-base resolution across various plant species. Integrating ribosomal RNA, transfer RNA and messenger RNA Ψ stoichiometry with mRNA abundance and polysome profiling data, we uncover a multilayered regulation of translation efficiency through Ψ modifications. rRNA pseudouridylation could globally control translation, although the effects vary at different rRNA Ψ sites. Ψ in the tRNA T-arm loop shows strong positive correlations between Ψ stoichiometry and the translation efficiency of their respective codons. We observed a general inverse correlation between Ψ level and mRNA stability, but a positive correlation with translation efficiency in Arabidopsis seedlings. In conclusion, our study provides critical resources for Ψ research in plants and proposes prevalent translation regulation through rRNA, tRNA and mRNA pseudouridylation.

RNA Decay Assay: 5-Ethynyl-Uridine Labeling and Chasing

Eukaryotic RNA synthesis and degradation are intricately regulated, impacting on gene expression dynamics. RNA stability varies in individual transcripts and is modulated by trans-acting factors such as microRNAs, long noncoding RNAs, and RNA-binding proteins, which determine protein output and subsequent cellular processes. To measure RNA decay rate, accurate and reliable methodologies are essential in the field of RNA biology. Transcription inhibition and metabolic labeling enable comprehensive investigations on RNA decay, offering critical insights into dynamic regulation of RNA decay. Transcription shut-off has been employed widely by using various approaches, such as treatment with chemical inhibitors or generation of temperature-sensitive mutants of RNA polymerases. However, it has limitations, providing a static view and lacking real-time dynamics as well as precise measurement of decay rate. Metabolic labeling, a method of incorporating modified nucleotides into RNA transcripts, complements shut-off approaches, allowing selective monitoring of newly synthesized RNA and tracing decay intermediates. The purpose of the protocol described in this chapter is to assess the kinetics and statics of newly synthesized RNA and its decay by 5-ethynyl uridine labeling.

Dynamic changes in mRNA nucleocytoplasmic localization in the nitrate response of Arabidopsis roots

Nitrate is a nutrient and signal that regulates gene expression. The nitrate response has been extensively characterized at the organism, organ, and cell-type-specific levels, but intracellular mRNA dynamics remain unexplored. To characterize nuclear and cytoplasmic transcriptome dynamics in response to nitrate, we performed a time-course expression analysis after nitrate treatment in isolated nuclei, cytoplasm, and whole roots. We identified 402 differentially localized transcripts (DLTs) in response to nitrate treatment. Induced DLT genes showed rapid and transient recruitment of the RNA polymerase II, together with an increase in the mRNA turnover rates. DLTs code for genes involved in metabolic processes, localization, and response to stimulus indicating DLTs include genes with relevant functions for the nitrate response that have not been previously identified. Using single-molecule RNA FISH, we observed early nuclear accumulation of the NITRATE REDUCTASE 1 (NIA1) transcripts in their transcription sites. We found that transcription of NIA1, a gene showing delayed cytoplasmic accumulation, is rapidly and transiently activated; however, its transcripts become unstable when they reach the cytoplasm. Our study reveals the dynamic localization of mRNAs between the nucleus and cytoplasm as an emerging feature in the temporal control of gene expression in response to nitrate treatment in Arabidopsis roots.

Translational activation by a synthetic PPR protein elucidates control of psbA translation in Arabidopsis chloroplasts

Translation initiation on chloroplast psbA mRNA in plants scales with light intensity, providing its gene product, D1, as needed to replace photodamaged D1 in Photosystem II. The psbA translational activator HIGH CHLOROPHYLL FLUORESCENCE 173 (HCF173) has been hypothesized to mediate this regulation. HCF173 belongs to the short-chain dehydrogenase/reductase superfamily, associates with the psbA 5'-untranslated region (5'-UTR), and has been hypothesized to enhance translation by binding an RNA segment that would otherwise pair with and mask the ribosome binding region. To test these hypotheses, we examined whether a synthetic pentatricopeptide repeat (sPPR) protein can substitute for HCF173 when bound to the HCF173 binding site. We show that an sPPR designed to bind HCF173's footprint in the psbA 5'-UTR bound the intended site in vivo and partially substituted for HCF173 to activate psbA translation. However, sPPR-activated translation did not respond to light. These results imply that HCF173 activates translation, at least in part, by sequestering the RNA it binds to maintain an accessible ribosome binding region, and that HCF173 is also required to regulate psbA translation in response to light. Translational activation can be added to the functions that can be programmed with sPPR proteins for synthetic biology applications in chloroplasts.

Widespread position-dependent transcriptional regulatory sequences in plants

Much of what we know about eukaryotic transcription stems from animals and yeast; however, plants evolved separately for over a billion years, leaving ample time for divergence in transcriptional regulation. Here we set out to elucidate fundamental properties of cis-regulatory sequences in plants. Using massively parallel reporter assays across four plant species, we demonstrate the central role of sequences downstream of the transcription start site (TSS) in transcriptional regulation. Unlike animal enhancers that are position independent, plant regulatory elements depend on their position, as altering their location relative to the TSS significantly affects transcription. We highlight the importance of the region downstream of the TSS in regulating transcription by identifying a DNA motif that is conserved across vascular plants and is sufficient to enhance gene expression in a dose-dependent manner. The identification of a large number of position-dependent enhancers points to fundamental differences in gene regulation between plants and animals.

Study of the RNA splicing kinetics via in vivo 5-EU labeling

Splicing is an important step of gene expression in all eukaryotes. Splice sites might be used with different efficiency, giving rise to alternative splicing products. At the same time, splice sites might be used at a variable rate. We used 5-ethynyl uridine labeling to sequence a nascent transcriptome of HeLa cells and deduced the rate of splicing for each donor and acceptor splice site. The following correlation analysis showed a correspondence of primary transcript features with the rate of splicing. Some dependencies we revealed were anticipated, such as a splicing rate decrease with a decreased complementarity of the donor splice site to U1 and acceptor sites to U2 snRNAs. Other dependencies were more surprising, like a negative influence of a distance to the 5' end on the rate of the acceptor splicing site utilization, or the differences in splicing rate between long, short, and RBM17-dependent introns. We also observed a deceleration of last intron splicing with an increase of the distance to the poly(A) site, which might be explained by the cooperativity of the splicing and polyadenylation. Additional analysis of splicing kinetics of SF3B4 knockdown cells suggested the impairment of a U2 snRNA recognition step. As a result, we deconvoluted the effects of several examined features on the splicing rate into a single regression model. The data obtained here are useful for further studies in the field, as they provide general splicing rate dependencies as well as help to justify the existence of slowly removed splice sites.

TMK4-mediated FIP37 phosphorylation regulates auxin-triggered N6-methyladenosine modification of auxin biosynthetic genes in Arabidopsis

The dynamics of N6-methyladenosine (m6A) mRNA modification are tightly controlled by the m6A methyltransferase complex and demethylases. Here, we find that auxin treatment alters m6A modification on auxin-responsive genes. Mechanically, TRANSMEMBRANE KINASE 4 (TMK4), a component of the auxin signaling pathway, interacts with and phosphorylates FKBP12-INTERACTING PROTEIN 37 (FIP37), a core component of the m6A methyltransferase complex, in an auxin-dependent manner. Phosphorylation of FIP37 enhances its interaction with RNA, thereby increasing m6A modification on its target genes, such as NITRILASE 1 (NIT1), a gene involved in indole-3-acetic acid (IAA) biosynthesis. 1-Naphthalacetic acid (NAA) treatment accelerates the mRNA decay of NIT1, in a TMK4- and FIP37-dependent manner, which leads to inhibition of auxin biosynthesis. Our findings identify a regulatory mechanism by which auxin modulates m6A modification through the phosphorylation of FIP37, ultimately affecting mRNA stability and auxin biosynthesis in plants.

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.

Quantitative profiling of m6A at single base resolution across the life cycle of rice and Arabidopsis

N6-methyladenosine (m6A) plays critical roles in regulating mRNA metabolism. However, comprehensive m6A methylomes in different plant tissues with single-base precision have yet to be reported. Here, we present transcriptome-wide m6A maps at single-base resolution in different tissues of rice and Arabidopsis using m6A-SAC-seq. Our analysis uncovers a total of 205,691 m6A sites distributed across 22,574 genes in rice, and 188,282 m6A sites across 19,984 genes in Arabidopsis. The evolutionarily conserved m6A sites in rice and Arabidopsis ortholog gene pairs are involved in controlling tissue development, photosynthesis and stress response. We observe an overall mRNA stabilization effect by 3' UTR m6A sites in certain plant tissues. Like in mammals, a positive correlation between the m6A level and the length of internal exons is also observed in plant mRNA, except for the last exon. Our data suggest an active m6A deposition process occurring near the stop codon in plant mRNA. In addition, the MTA-installed plant mRNA m6A sites correlate with both translation promotion and translation suppression, depicting a more complicated regulatory picture. Our results therefore provide in-depth resources for relating single-base resolution m6A sites with functions in plants and uncover a suppression-activation model controlling m6A biogenesis across species.

Differences in transcription initiation directionality underlie distinctions between plants and animals in chromatin modification patterns at genes and cis-regulatory elements

Transcriptional initiation is among the first regulated steps controlling eukaryotic gene expression. High-throughput profiling of fungal and animal genomes has revealed that RNA Polymerase II often initiates transcription in both directions at the promoter transcription start site, but generally only elongates productively into the gene body. Additionally, Pol II can initiate transcription in both directions at cis-regulatory elements such as enhancers. These bidirectional RNA Polymerase II initiation events can be observed directly with methods that capture nascent transcripts, and they are also revealed indirectly by the presence of transcription-associated histone modifications on both sides of the transcription start site or cis-regulatory elements. Previous studies have shown that nascent RNAs and transcription-associated histone modifications in the model plant Arabidopsis thaliana accumulate mainly in the gene body, suggesting that transcription does not initiate widely in the upstream direction from genes in this plant. We compared transcription-associated histone modifications and nascent transcripts at both transcription start sites and cis-regulatory elements in A. thaliana, Drosophila melanogaster, and Homo sapiens. Our results provide evidence for mostly unidirectional RNA Polymerase II initiation at both promoters and gene-proximal cis-regulatory elements of A. thaliana, whereas bidirectional transcription initiation is observed widely at promoters in both D. melanogaster and H. sapiens, as well as cis-regulatory elements in Drosophila. Furthermore, the distribution of transcription-associated histone modifications around transcription start sites in the Oryza sativa (rice) and Glycine max (soybean) genomes suggests that unidirectional transcription initiation is the norm in these genomes as well. These results suggest that there are fundamental differences in transcriptional initiation directionality between flowering plant and metazoan genomes, which are manifested as distinct patterns of chromatin modifications around RNA polymerase initiation sites.

Dynamics of mRNA fate during light stress and recovery: from transcription to stability and translation

Transcript stability is an important determinant of its abundance and, consequently, translational output. Transcript destabilisation can be rapid and is well suited for modulating the cellular response. However, it is unclear the extent to which RNA stability is altered under changing environmental conditions in plants. We previously hypothesised that recovery-induced transcript destabilisation facilitated a phenomenon of rapid recovery gene downregulation (RRGD) in Arabidopsis thaliana (Arabidopsis) following light stress, based on mathematical calculations to account for ongoing transcription. Here, we test this hypothesis and investigate processes regulating transcript abundance and fate by quantifying changes in transcription, stability and translation before, during and after light stress. We adapt syringe infiltration to apply a transcriptional inhibitor to soil-grown plants in combination with stress treatments. Compared with measurements in juvenile plants and cell culture, we find reduced stability across a range of transcripts encoding proteins involved in RNA binding and processing. We also observe light-induced destabilisation of transcripts, followed by their stabilisation during recovery. We propose that this destabilisation facilitates RRGD, possibly in combination with transcriptional shut-off that was confirmed for HSP101, ROF1 and GOLS1. We also show that translation remains highly dynamic over the course of light stress and recovery, with a bias towards transcript-specific increases in ribosome association, independent of changes in total transcript abundance, after 30 min of light stress. Taken together, we provide evidence for the combinatorial regulation of transcription and stability that occurs to coordinate translation during light stress and recovery in Arabidopsis.

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.

Comparative analysis of nascent RNA sequencing methods and their applications in studies of cotranscriptional splicing dynamics

High-throughput detection of nascent RNA is critical for studies of transcription and much more challenging than that of mRNA. Recently, several massively parallel nascent RNA sequencing methods were established in eukaryotic cells. Here, we systematically compared 3 classes of methods on the same pure or crude nuclei preparations: GRO-seq for sequence nuclear run-on RNAs, pNET-seq for sequence RNA polymerase II-associated RNAs, and CB RNA-seq for sequence chromatin-bound (CB) RNAs in Arabidopsis (Arabidopsis thaliana). To improve the resolution of CB RNAs, 3'CB RNA-seq was established to sequence the 3' ends of CB RNAs. In addition, we modified pNET-seq to establish the Chromatin Native Elongation Transcript sequencing (ChrNET) method using chromatin as the starting material for RNA immunoprecipitation. Reproducibility, sensitivity and accuracy in detecting nascent transcripts, experimental procedures, and costs were analyzed, which revealed the strengths and weaknesses of each method. We found that pNET and GRO methods best detected active RNA polymerase II. CB RNA-seq is a simple and cost-effective alternative for nascent RNA studies, due to its high correlation with pNET-seq and GRO-seq. Compared with pNET, ChrNET has higher specificity for nascent RNA capture and lower sequencing cost. 3'CB is sensitive to transcription-coupled splicing. Using these methods, we identified 1,404 unknown transcripts, 4,482 unannotated splicing events, and 60 potential recursive splicing events. This comprehensive comparison of different nascent/chromatin RNA sequencing methods highlights the strengths of each method and serves as a guide for researchers aiming to select a method that best meets their study goals.

Differences in transcription initiation directionality underlie distinctions between plants and animals in chromatin modification patterns at genes and cis-regulatory elements

Transcriptional initiation is among the first regulated steps controlling eukaryotic gene expression. High-throughput profiling of fungal and animal genomes has revealed that RNA Polymerase II (Pol II) often initiates transcription in both directions at the promoter transcription start site (TSS), but generally only elongates productively into the gene body. Additionally, Pol II can initiate transcription in both directions at cis-regulatory elements (CREs) such as enhancers. These bidirectional Pol II initiation events can be observed directly with methods that capture nascent transcripts, and they are also revealed indirectly by the presence of transcription-associated histone modifications on both sides of the TSS or CRE. Previous studies have shown that nascent RNAs and transcription-associated histone modifications in the model plant Arabidopsis thaliana accumulate mainly in the gene body, suggesting that transcription does not initiate widely in the upstream direction from genes in this plant. We compared transcription-associated histone modifications and nascent transcripts at both TSSs and CREs in Arabidopsis thaliana, Drosophila melanogaster, and Homo sapiens. Our results provide evidence for mostly unidirectional Pol II initiation at both promoters and gene-proximal CREs of Arabidopsis thaliana, whereas bidirectional transcription initiation is observed widely at promoters in both Drosophila melanogaster and Homo sapiens, as well as CREs in Drosophila. Furthermore, the distribution of transcription-associated histone modifications around TSSs in the Oryza sativa (rice) and Glycine max (soybean) genomes suggests that unidirectional transcription initiation is the norm in these genomes as well. These results suggest that there are fundamental differences in transcriptional initiation directionality between flowering plant and metazoan genomes, which are manifested as distinct patterns of chromatin modifications around RNA polymerase initiation sites.

Plant organellar RNA maturation

Plant organellar RNA metabolism is run by a multitude of nucleus-encoded RNA-binding proteins (RBPs) that control RNA stability, processing, and degradation. In chloroplasts and mitochondria, these post-transcriptional processes are vital for the production of a small number of essential components of the photosynthetic and respiratory machinery-and consequently for organellar biogenesis and plant survival. Many organellar RBPs have been functionally assigned to individual steps in RNA maturation, often specific to selected transcripts. While the catalog of factors identified is ever-growing, our knowledge of how they achieve their functions mechanistically is far from complete. This review summarizes the current knowledge of plant organellar RNA metabolism taking an RBP-centric approach and focusing on mechanistic aspects of RBP functions and the kinetics of the processes they are involved in.

Application of nucleoside or nucleotide analogues in RNA dynamics and RNA-binding protein analysis

Cellular RNAs undergo dynamic changes during RNA biological processes, which are tightly orchestrated by RNA-binding proteins (RBPs). Yet, the investigation of RNA dynamics is hurdled by highly abundant steady-state RNAs, which make the signals of dynamic RNAs less detectable. Notably, the exert of nucleoside or nucleotide analogue-based RNA technologies has provided a remarkable platform for RNA dynamics research, revealing diverse unnoticed features in RNA metabolism. In this review, we focus on the application of two types of analogue-based RNA sequencing, antigen-/antibody- and click chemistry-based methodologies, and summarize the RNA dynamics features revealed. Moreover, we discuss emerging single-cell newly transcribed RNA sequencing methodologies based on nucleoside analogue labeling, which provides novel insights into RNA dynamics regulation at single-cell resolution. On the other hand, we also emphasize the identification of RBPs that interact with polyA, non-polyA RNAs, or newly transcribed RNAs and also their associated RNA-binding domains at genomewide level through ultraviolet crosslinking and mass spectrometry in different contexts. We anticipated that further modification and development of these analogue-based RNA and RBP capture technologies will aid in obtaining an unprecedented understanding of RNA biology. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry RNA Methods > RNA Analyses in Cells.

Monitoring mRNA Half-Life in Arabidopsis Using Droplet Digital PCR

mRNA decay is an important process in post-transcriptional regulation; in addition, it plays a crucial role in plant development and response to stress. The development of new tools to quantify mRNA decay intermediates is thus important to better characterize the dynamic of mRNA decay in various conditions. Here, we applied droplet digital PCR (ddPCR), a recent and precise PCR technology, to determine mRNA half-life in Arabidopsis seedlings. We demonstrated that ddPCR can correctly assess mRNA half-life from a wide variety of transcripts in a reproducible manner. We also demonstrated that thanks to multiplexing mRNA, the half-life of multiple transcripts can be followed in the same reaction. As ddPCR allows precise quantification, we proposed that this approach is highly suitable when a low amount of RNA is available; for the detection of many targets or for the analysis of lowly expressed transcripts.

Antisense Transcription in Plants: A Systematic Review and an Update on cis-NATs of Sugarcane

Initially, natural antisense transcripts (NATs, natRNAs, or asRNAs) were considered repressors; however, their functions in gene regulation are diverse. Positive, negative, or neutral correlations to the cognate gene expression have been noted. Although the first studies were published about 50 years ago, there is still much to be investigated regarding antisense transcripts in plants. A systematic review of scientific publications available in the Web of Science databases was conducted to contextualize how the studying of antisense transcripts has been addressed. Studies were classified considering three categories: “Natural antisense” (208), artificial antisense used in “Genetic Engineering” (797), or “Natural antisense and Genetic Engineering”-related publications (96). A similar string was used for a systematic search in the NCBI Gene database. Of the 1132 antisense sequences found for plants, only 0.8% were cited in PubMed and had antisense information confirmed. This value was the lowest when compared to fungi (2.9%), bacteria (2.3%), and mice (54.1%). Finally, we present an update for the cis-NATs identified in Saccharum spp. Of the 1413 antisense transcripts found in different experiments, 25 showed concordant expressions, 22 were discordant, 1264 did not correlate with the cognate genes, and 102 presented variable results depending on the experiment.

An atlas of plant full-length RNA reveals tissue-specific and monocots-dicots conserved regulation of poly(A) tail length

Poly(A) tail is a hallmark of eukaryotic messenger RNA and its length plays an essential role in regulating mRNA metabolism. However, a comprehensive resource for plant poly(A) tail length has yet to be established. Here, we applied a poly(A)-enrichment-free, nanopore-based method to profile full-length RNA with poly(A) tail information in plants. Our atlas contains over 120 million polyadenylated mRNA molecules from seven different tissues of Arabidopsis, as well as the shoot tissue of maize, soybean and rice. In most tissues, the size of plant poly(A) tails shows peaks at approximately 20 and 45 nucleotides, while the poly(A) tails in pollen exhibit a distinct pattern with strong peaks centred at 55 and 80 nucleotides. Moreover, poly(A) tail length is regulated in a gene-specific manner-mRNAs with short half-lives in general have long poly(A) tails, while mRNAs with long half-lives are featured with relatively short poly(A) tails that peak at ~45 nucleotides. Across species, poly(A) tails in the nucleus are almost twice as long as in the cytoplasm. Our comprehensive dataset lays the groundwork for future functional and evolutionary studies on poly(A) tail length regulation in plants.

Epigenetic regulation of thermomorphogenesis and heat stress tolerance

Many environmental conditions fluctuate and organisms need to respond effectively. This is especially true for temperature cues that can change in minutes to seasons and often follow a diurnal rhythm. Plants cannot migrate and most cannot regulate their temperature. Therefore, a broad array of responses have evolved to deal with temperature cues from freezing to heat stress. A particular response to mildly elevated temperatures is called thermomorphogenesis, a suite of morphological adaptations that includes thermonasty, formation of thin leaves and elongation growth of petioles and hypocotyl. Thermomorphogenesis allows for optimal performance in suboptimal temperature conditions by enhancing the cooling capacity. When temperatures rise further, heat stress tolerance mechanisms can be induced that enable the plant to survive the stressful temperature, which typically comprises cellular protection mechanisms and memory thereof. Induction of thermomorphogenesis, heat stress tolerance and stress memory depend on gene expression regulation, governed by diverse epigenetic processes. In this Tansley review we update on the current knowledge of epigenetic regulation of heat stress tolerance and elevated temperature signalling and response, with a focus on thermomorphogenesis regulation and heat stress memory. In particular we highlight the emerging role of H3K4 methylation marks in diverse temperature signalling pathways.

Synthesis of 4-thiouridines with prodrug functionalization for RNA metabolic labeling

Metabolic labeling has emerged as a powerful tool to endow RNA with reactive handles allowing for subsequent chemical derivatization and processing. Recently, thiolated nucleosides, such as 4-thiouridine (4sU), have attracted great interest in metabolic labeling-based RNA sequencing approaches (TUC-seq, SLAM-seq, TimeLapse-seq) to study cellular RNA expression and decay dynamics. For these and other applications (e.g. PAR-CLIP), thus far only the naked nucleoside 4sU has been applied. Here we examined the concept of derivatizing 4sU into a 5'-monophosphate prodrug that would allow for cell permeation and potentially improve labeling efficiency by bypassing the rate-limiting first step of 5' phosphorylation of the nucleoside into the ultimately bioactive 4sU triphosphate (4sUTP). To this end, we developed robust synthetic routes towards diverse 4sU monophosphate prodrugs. Using metabolic labeling assays, we found that most of the newly introduced 4sU prodrugs were well tolerated by the cells. One derivative, the bis(4-acetyloxybenzyl) 5'-monophosphate of 4sU, was also efficiently incorporated into nascent RNA.

The transcriptional response to oxidative stress is independent of stress-granule formation

Cells respond to stress with translational arrest, robust transcriptional changes, and transcription-independent formation of mRNP assemblies termed stress granules (SGs). Despite considerable interest in the role of SGs in oxidative, unfolded protein and viral stress responses, whether and how SGs contribute to stress-induced transcription have not been rigorously examined. To address this, we characterized transcriptional changes in Drosophila S2 cells induced by acute oxidative-stress and assessed how these were altered under conditions that disrupted SG assembly. Oxidative stress for 3 h predominantly resulted in induction or up-regulation of stress-responsive mRNAs whose levels peaked during recovery after stress cessation. The stress transcriptome is enriched in mRNAs coding for chaperones including HSP70s, small heat shock proteins, glutathione transferases, and several noncoding RNAs. Oxidative stress also induced cytoplasmic SGs that disassembled 3 h after stress cessation. As expected, RNAi-mediated knockdown of the conserved G3BP1/Rasputin protein inhibited SG assembly. However, this disruption had no significant effect on the stress-induced transcriptional response or stress-induced translational arrest. Thus SG assembly and stress-induced gene expression alterations appear to be driven by distinctive signaling processes. We suggest that while SG assembly represents a fast, transient mechanism, the transcriptional response enables a slower, longer-lasting mechanism for adaptation to and recovery from cell stress.

Conditional and tissue-specific approaches to dissect essential mechanisms in plant development

Reverse genetics approaches are routinely used to investigate gene function. However, mutations, especially in critical genes, can lead to pleiotropic effects as severe as lethality, thus limiting functional studies in specific contexts. Approaches that allow for modifications of genes or gene products in a specific spatial or temporal setting can overcome these limitations. The advent of CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats) technologies has not only revolutionized targeted genome modification in plants but also enabled new possibilities for inducible and tissue-specific manipulation of gene functions at the DNA and RNA levels. In addition, novel approaches for the direct manipulation of target proteins have been introduced in plant systems. Here, we review the current development in tissue-specific and conditional manipulation approaches at the DNA, RNA, and protein levels.

The life and death of RNA across temperatures

Temperature is an environmental condition that has a pervasive effect on cells along with all the molecules and reactions in them. The mechanisms by which prototypical RNA molecules sense and withstand heat have been identified mostly in bacteria and archaea. The relevance of these phenomena is, however, broader, and similar mechanisms have been recently found throughout the tree of life, from sex determination in reptiles to adaptation of viral RNA polymerases, to genetic disorders in humans. We illustrate the temperature dependence of RNA metabolism with examples from the synthesis to the degradation of mRNAs, and review recently emerged questions. Are cells exposed to greater temperature variations and gradients than previously surmised? How do cells reconcile the conflicting thermal stability requirements of primary and tertiary structures of RNAs? To what extent do enzymes contribute to the temperature compensation of the reaction rates in mRNA turnover by lowering the energy barrier of the catalyzed reactions? We conclude with the ecological, forensic applications of the temperature-dependence of RNA degradation and the biotechnological aspects of mRNA vaccine production.

Independent parental contributions initiate zygote polarization in Arabidopsis thaliana

Embryogenesis of flowering plants is initiated by polarization of the zygote, a prerequisite for correct axis formation in the embryo. The daughter cells of the asymmetric zygote division form the pro-embryo and the mostly extra-embryonic suspensor.¹ The suspensor plays a pivotal role in nutrient and hormone transport and rapid growth of the embryo.^2,3 Zygote polarization is controlled by a MITOGEN-ACTIVATING PROTEIN (MAP) kinase signaling pathway including the MAPKK kinase (MAP3K) YODA (YDA)⁴ and the upstream membrane-associated proteins BRASINOSTEROID SIGNALING KINASE 1 (BSK1) and BSK2.^5,6 Furthermore, suspensor development is controlled by cysteine-rich peptides of the EMBRYO SURROUNDING FACTOR 1 (ESF1) family.⁷ While they act genetically upstream of YDA, the corresponding receptor to perceive these potential ligands is unknown. In other developmental processes, such as stomata development, YDA activity is controlled by receptor kinases of the ERECTA family (ERf).^8-12 While the receptor kinases upstream of BSK1/2 in the embryo have so far not been identified,¹ YDA is in part activated by the sperm cell-derived BSK family member SHORT SUSPENSOR (SSP) that represents a naturally occurring, constitutively active variant of BSK1.^5,13 It has been speculated that SSP might be a paternal component of a parental tug-of-war controlling resource allocation toward the embryo.^2,13 Here, we show that in addition to SSP, the receptor kinase ERECTA plays a crucial role in zygote polarization as a maternally contributed part of the embryonic YDA pathway. We conclude that two independent parental contributions initiate zygote polarization and control embryo development.

Making sense of the natural antisense transcript puzzle

In plants, thousands of genes are associated with antisense transcription, which often produces noncoding RNAs. Although widespread, sense-antisense pairs have been implicated in a limited variety of functions in plants and are often thought to form extensive dsRNA stretches triggering gene silencing. In this opinion, we show that evidence does not support gene silencing as a major role for antisense transcription. In fact, it is more likely that antisense transcripts play diverse functions in gene regulation. We propose a general framework for the initial functional dissection of antisense transcripts, suggesting testable hypotheses relying on an experiment-based decision tree. By moving beyond the gene silencing paradigm, we argue that a broad and diverse role for natural antisense transcription will emerge.

Decoding co-/post-transcriptional complexities of plant transcriptomes and epitranscriptome using next-generation sequencing technologies

Next-generation sequencing (NGS) technologies - Illumina RNA-seq, Pacific Biosciences isoform sequencing (PacBio Iso-seq), and Oxford Nanopore direct RNA sequencing (DRS) - have revealed the complexity of plant transcriptomes and their regulation at the co-/post-transcriptional level. Global analysis of mature mRNAs, transcripts from nuclear run-on assays, and nascent chromatin-bound mRNAs using short as well as full-length and single-molecule DRS reads have uncovered potential roles of different forms of RNA polymerase II during the transcription process, and the extent of co-transcriptional pre-mRNA splicing and polyadenylation. These tools have also allowed mapping of transcriptome-wide start sites in cap-containing RNAs, poly(A) site choice, poly(A) tail length, and RNA base modifications. The emerging theme from recent studies is that reprogramming of gene expression in response to developmental cues and stresses at the co-/post-transcriptional level likely plays a crucial role in eliciting appropriate responses for optimal growth and plant survival under adverse conditions. Although the mechanisms by which developmental cues and different stresses regulate co-/post-transcriptional splicing are largely unknown, a few recent studies indicate that the external cues target spliceosomal and splicing regulatory proteins to modulate alternative splicing. In this review, we provide an overview of recent discoveries on the dynamics and complexities of plant transcriptomes, mechanistic insights into splicing regulation, and discuss critical gaps in co-/post-transcriptional research that need to be addressed using diverse genomic and biochemical approaches.

Epigenetic regulation of abiotic stress memory: maintaining the good things while they last

As sessile organisms, plants have evolved sophisticated ways to constantly gauge and adapt to changing environmental conditions including extremes that may be harmful to their growth and development and are thus perceived as stress. In nature, stressful events are often chronic or recurring and thus an initial stress may prime a plant to respond more efficiently to a subsequent stress event. An epigenetic basis of such stress memory was long postulated and in recent years it has been shown that this is indeed the case. High temperature stress has proven an excellent system to unpick the molecular basis of somatic stress memory, which includes histone modifications and nucleosome occupancy. This review discusses recent findings and pinpoints open questions in the field.

Acclimation in plants - the Green Hub consortium

Acclimation is the capacity to adapt to environmental changes within the lifetime of an individual. This ability allows plants to cope with the continuous variation in ambient conditions to which they are exposed as sessile organisms. Because environmental changes and extremes are becoming even more pronounced due to the current period of climate change, enhancing the efficacy of plant acclimation is a promising strategy for mitigating the consequences of global warming on crop yields. At the cellular level, the chloroplast plays a central role in many acclimation responses, acting both as a sensor of environmental change and as a target of cellular acclimation responses. In this Perspective article, we outline the activities of the Green Hub consortium funded by the German Science Foundation. The main aim of this research collaboration is to understand and strategically modify the cellular networks that mediate plant acclimation to adverse environments, employing Arabidopsis, tobacco (Nicotiana tabacum) and Chlamydomonas as model organisms. These efforts will contribute to 'smart breeding' methods designed to create crop plants with improved acclimation properties. To this end, the model oilseed crop Camelina sativa is being used to test modulators of acclimation for their potential to enhance crop yield under adverse environmental conditions. Here we highlight the current state of research on the role of gene expression, metabolism and signalling in acclimation, with a focus on chloroplast-related processes. In addition, further approaches to uncovering acclimation mechanisms derived from systems and computational biology, as well as adaptive laboratory evolution with photosynthetic microbes, are highlighted.

The TUTase URT1 connects decapping activators and prevents the accumulation of excessively deadenylated mRNAs to avoid siRNA biogenesis

Uridylation is a widespread modification destabilizing eukaryotic mRNAs. Yet, molecular mechanisms underlying TUTase-mediated mRNA degradation remain mostly unresolved. Here, we report that the Arabidopsis TUTase URT1 participates in a molecular network connecting several translational repressors/decapping activators. URT1 directly interacts with DECAPPING 5 (DCP5), the Arabidopsis ortholog of human LSM14 and yeast Scd6, and this interaction connects URT1 to additional decay factors like DDX6/Dhh1-like RNA helicases. Nanopore direct RNA sequencing reveals a global role of URT1 in shaping poly(A) tail length, notably by preventing the accumulation of excessively deadenylated mRNAs. Based on in vitro and in planta data, we propose a model that explains how URT1 could reduce the accumulation of oligo(A)-tailed mRNAs both by favoring their degradation and because 3' terminal uridines intrinsically hinder deadenylation. Importantly, preventing the accumulation of excessively deadenylated mRNAs avoids the biogenesis of illegitimate siRNAs that silence endogenous mRNAs and perturb Arabidopsis growth and development.

Barley Anther and Meiocyte Transcriptome Dynamics in Meiotic Prophase I

In flowering plants, successful germinal cell development and meiotic recombination depend upon a combination of environmental and genetic factors. To gain insights into this specialized reproductive development program we used short- and long-read RNA-sequencing (RNA-seq) to study the temporal dynamics of transcript abundance in immuno-cytologically staged barley (Hordeum vulgare) anthers and meiocytes. We show that the most significant transcriptional changes in anthers occur at the transition from pre-meiosis to leptotene-zygotene, which is followed by increasingly stable transcript abundance throughout prophase I into metaphase I-tetrad. Our analysis reveals that the pre-meiotic anthers are enriched in long non-coding RNAs (lncRNAs) and that entry to meiosis is characterized by their robust and significant down regulation. Intriguingly, only 24% of a collection of putative meiotic gene orthologs showed differential transcript abundance in at least one stage or tissue comparison. Argonautes, E3 ubiquitin ligases, and lys48 specific de-ubiquitinating enzymes were enriched in prophase I meiocyte samples. These developmental, time-resolved transcriptomes demonstrate remarkable stability in transcript abundance in meiocytes throughout prophase I after the initial and substantial reprogramming at meiosis entry and the complexity of the regulatory networks involved in early meiotic processes.

Comparative analysis of the transcriptomes of EDL, psoas, and soleus muscles from mice

BACKGROUND

Individual skeletal muscles have evolved to perform specific tasks based on their molecular composition. In general, muscle fibers are characterized as either fast-twitch or slow-twitch based on their myosin heavy chain isoform profiles. This approach made sense in the early days of muscle studies when SDS-PAGE was the primary tool for mapping fiber type. However, Next Generation Sequencing tools permit analysis of the entire muscle transcriptome in a single sample, which allows for more precise characterization of differences among fiber types, including distinguishing between different isoforms of specific proteins. We demonstrate the power of this approach by comparing the differential gene expression patterns of extensor digitorum longus (EDL), psoas, and soleus from mice using high throughput RNA sequencing.

RESULTS

EDL and psoas are typically classified as fast-twitch muscles based on their myosin expression pattern, while soleus is considered a slow-twitch muscle. The majority of the transcriptomic variability aligns with the fast-twitch and slow-twitch characterization. However, psoas and EDL exhibit unique expression patterns associated with the genes coding for extracellular matrix, myofibril, transcription, translation, striated muscle adaptation, mitochondrion distribution, and metabolism. Furthermore, significant expression differences between psoas and EDL were observed in genes coding for myosin light chain, troponin, tropomyosin isoforms, and several genes encoding the constituents of the Z-disk.

CONCLUSIONS

The observations highlight the intricate molecular nature of skeletal muscles and demonstrate the importance of utilizing transcriptomic information as a tool for skeletal muscle characterization.

Long non-coding RNAs in plants: emerging modulators of gene activity in development and stress responses

MAIN CONCLUSION

Long non-coding RNAs modulate gene activity in plant development and stress responses by various molecular mechanisms. Long non-coding RNAs (lncRNAs) are transcripts larger than 200 nucleotides without protein coding potential. Computational approaches have identified numerous lncRNAs in different plant species. Research in the past decade has unveiled that plant lncRNAs participate in a wide range of biological processes, including regulation of flowering time and morphogenesis of reproductive organs, as well as abiotic and biotic stress responses. LncRNAs execute their functions by interacting with DNA, RNA and protein molecules, and by modulating the expression level of their targets through epigenetic, transcriptional, post-transcriptional or translational regulation. In this review, we summarize characteristics of plant lncRNAs, discuss recent progress on understanding of lncRNA functions, and propose an experimental framework for functional characterization.