Transcript isoform sequencing reveals widespread promoter-proximal transcriptional termination in Arabidopsis

Cytoplasmic binding partners of the Integrator endonuclease INTS11 and its paralog CPSF73 are required for their nuclear function

INTS11 and CPSF73 are metal-dependent endonucleases for Integrator and pre-mRNA 3'-end processing, respectively. Here, we show that the INTS11 binding partner BRAT1/CG7044, a factor important for neuronal fitness, stabilizes INTS11 in the cytoplasm and is required for Integrator function in the nucleus. Loss of BRAT1 in neural organoids leads to transcriptomic disruption and precocious expression of neurogenesis-driving transcription factors. The structures of the human INTS9-INTS11-BRAT1 and Drosophila dIntS11-CG7044 complexes reveal that the conserved C terminus of BRAT1/CG7044 is captured in the active site of INTS11, with a cysteine residue directly coordinating the metal ions. Inspired by these observations, we find that UBE3D is a binding partner for CPSF73, and UBE3D likely also uses a conserved cysteine residue to directly coordinate the active site metal ions. Our studies have revealed binding partners for INTS11 and CPSF73 that behave like cytoplasmic chaperones with a conserved impact on the nuclear functions of these enzymes.

Enhancers associated with unstable RNAs are rare in plants

Unstable transcripts have emerged as markers of active enhancers in vertebrates and shown to be involved in many cellular processes and medical disorders. However, their prevalence and role in plants is largely unexplored. Here, we comprehensively captured all actively initiating (nascent) transcripts across diverse crops and other plants using capped small (cs)RNA sequencing. We discovered that unstable transcripts are rare in plants, unlike in vertebrates, and when present, often originate from promoters. In addition, many 'distal' elements in plants initiate tissue-specific stable transcripts and are likely bona fide promoters of as-yet-unannotated genes or non-coding RNAs, cautioning against using reference genome annotations to infer putative enhancer sites. To investigate enhancer function, we integrated data from self-transcribing active regulatory region (STARR) sequencing. We found that annotated promoters and other regions that initiate stable transcripts, but not those marked by unstable or bidirectional unstable transcripts, showed stronger enhancer activity in this assay. Our findings underscore the blurred line between promoters and enhancers and suggest that cis-regulatory elements can encompass diverse structures and mechanisms in eukaryotes, including humans.

A data-driven genome annotation approach for cassava

Genome annotation files play a critical role in dictating the quality of downstream analyses by providing essential predictions for gene positions and structures. These files are pivotal in decoding the complex information encoded within DNA sequences. Here, we generated experimental data resolving RNA 5'- and 3'-ends as well as full-length RNAs for cassava TME12 sticklings in ambient temperature and cold. We used these data to generate genome annotation files using the TranscriptomeReconstructoR (TR) tool. A careful comparison to high-quality genome annotations suggests that our new TR genome annotations identified additional genes, resolved the transcript boundaries more accurately and identified additional RNA isoforms. We enhanced existing cassava genome annotation files with the information from TR that maintained the different transcript models as RNA isoforms. The resultant merged annotation was subsequently utilized for comprehensive analysis. To examine the effects of genome annotation files on gene expression studies, we compared the detection of differentially expressed genes during cold using the same RNA-seq data but alternative genome annotation files. We found that our merged genome annotation that included cold-specific TR gene models identified about twice as many cold-induced genes. These data indicate that environmentally induced genes may be missing in off-the-shelf genome annotation files. In conclusion, TR offers the opportunity to enhance crop genome annotations with implications for the discovery of differentially expressed candidate genes during plant-environment interactions.

RNAirport: a deep neural network-based database characterizing representative gene models in plants

A 5'-leader, known initially as the 5'-untranslated region, contains multiple isoforms due to alternative splicing (aS) and alternative transcription start site (aTSS). Therefore, a representative 5'-leader is demanded to examine the embedded RNA regulatory elements in controlling translation efficiency. Here, we develop a ranking algorithm and a deep-learning model to annotate representative 5'-leaders for five plant species. We rank the intra-sample and inter-sample frequency of aS-mediated transcript isoforms using the Kruskal-Wallis test-based algorithm and identify the representative aS-5'-leader. To further assign a representative 5'-end, we train the deep-learning model 5'leaderP to learn aTSS-mediated 5'-end distribution patterns from cap-analysis gene expression data. The model accurately predicts the 5'-end, confirmed experimentally in Arabidopsis and rice. The representative 5'-leader-contained gene models and 5'leaderP can be accessed at RNAirport (http://www.rnairport.com/leader5P/). The Stage 1 annotation of 5'-leader records 5'-leader diversity and will pave the way to Ribo-Seq open-reading frame annotation, identical to the project recently initiated by human GENCODE.

Transcript elongation by RNA polymerase II in plants: factors, regulation and impact on gene expression

Transcriptional elongation by RNA polymerase II (RNAPII) through chromatin is a dynamic and highly regulated step of eukaryotic gene expression. A combination of transcript elongation factors (TEFs) including modulators of RNAPII activity and histone chaperones facilitate efficient transcription on nucleosomal templates. Biochemical and genetic analyses, primarily performed in Arabidopsis, provided insight into the contribution of TEFs to establish gene expression patterns during plant growth and development. In addition to summarising the role of TEFs in plant gene expression, we emphasise in our review recent advances in the field. Thus, mechanisms are presented how aberrant intragenic transcript initiation is suppressed by repressing transcriptional start sites within coding sequences. We also discuss how transcriptional interference of ongoing transcription with neighbouring genes is prevented. Moreover, it appears that plants make no use of promoter-proximal RNAPII pausing in the way mammals do, but there are nucleosome-defined mechanism(s) that determine the efficiency of mRNA synthesis by RNAPII. Accordingly, a still growing number of processes related to plant growth, development and responses to changing environmental conditions prove to be regulated at the level of transcriptional elongation.

Size fractionated NET-Seq reveals a conserved architecture of transcription units around yeast genes

Genomes from yeast to humans are subject to pervasive transcription. A single round of pervasive transcription is sufficient to alter local chromatin conformation, nucleosome dynamics and gene expression, but is hard to distinguish from background signals. Size fractionated native elongating transcript sequencing (sfNET-Seq) was developed to precisely map nascent transcripts independent of expression levels. RNAPII-associated nascent transcripts are fractionation into different size ranges before library construction. When anchored to the transcription start sites (TSS) of annotated genes, the combined pattern of the output metagenes gives the expected reference pattern. Bioinformatic pattern matching to the reference pattern identified 9542 transcription units in Saccharomyces cerevisiae, of which 47% are coding and 53% are noncoding. In total, 3113 (33%) are unannotated noncoding transcription units. Anchoring all transcription units to the TSS or polyadenylation site (PAS) of annotated genes reveals distinctive architectures of linked pairs of divergent transcripts approximately 200nt apart. The Reb1 transcription factor is enriched 30nt downstream of the PAS only when an upstream (TSS -60nt with respect to PAS) noncoding transcription unit co-occurs with a downstream (TSS +150nt) coding transcription unit and acts to limit levels of upstream antisense transcripts. The potential for extensive transcriptional interference is evident from low abundance unannotated transcription units with variable TSS (median -240nt) initiating within a 500nt window upstream of, and transcribing over, the promoters of protein-coding genes. This study confirms a highly interleaved yeast genome with different types of transcription units altering the chromatin landscape in distinctive ways, with the potential to exert extensive regulatory control.

Population-level annotation of lncRNAs in Arabidopsis reveals extensive expression variation associated with transposable element-like silencing

Long noncoding RNAs (lncRNAs) are understudied and underannotated in plants. In mammals, lncRNA loci are nearly as ubiquitous as protein-coding genes, and their expression is highly variable between individuals of the same species. Using Arabidopsis thaliana as a model, we aimed to elucidate the true scope of lncRNA transcription across plants from different regions and study its natural variation. We used transcriptome deep sequencing data sets spanning hundreds of natural accessions and several developmental stages to create a population-wide annotation of lncRNAs, revealing thousands of previously unannotated lncRNA loci. While lncRNA transcription is ubiquitous in the genome, most loci appear to be actively silenced and their expression is extremely variable between natural accessions. This high expression variability is largely caused by the high variability of repressive chromatin levels at lncRNA loci. High variability was particularly common for intergenic lncRNAs (lincRNAs), where pieces of transposable elements (TEs) present in 50% of these lincRNA loci are associated with increased silencing and variation, and such lncRNAs tend to be targeted by the TE silencing machinery. We created a population-wide lncRNA annotation in Arabidopsis and improve our understanding of plant lncRNA genome biology, raising fundamental questions about what causes transcription and silencing across the genome.

LncRNAs exert indispensable roles in orchestrating the interaction among diverse noncoding RNAs and enrich the regulatory network of plant growth and its adaptive environmental stress response

With the advent of advanced sequencing technologies, non-coding RNAs (ncRNAs) are increasingly pivotal and play highly regulated roles in the modulation of diverse aspects of plant growth and stress response. This includes a spectrum of ncRNA classes, ranging from small RNAs to long non-coding RNAs (lncRNAs). Notably, among these, lncRNAs emerge as significant and intricate components within the broader ncRNA regulatory networks. Here, we categorize ncRNAs based on their length and structure into small RNAs, medium-sized ncRNAs, lncRNAs, and circle RNAs. Furthermore, the review delves into the detailed biosynthesis and origin of these ncRNAs. Subsequently, we emphasize the diverse regulatory mechanisms employed by lncRNAs that are located at various gene regions of coding genes, embodying promoters, 5'UTRs, introns, exons, and 3'UTR regions. Furthermore, we elucidate these regulatory modes through one or two concrete examples. Besides, lncRNAs have emerged as novel central components that participate in phase separation processes. Moreover, we illustrate the coordinated regulatory mechanisms among lncRNAs, miRNAs, and siRNAs with a particular emphasis on the central role of lncRNAs in serving as sponges, precursors, spliceosome, stabilization, scaffolds, or interaction factors to bridge interactions with other ncRNAs. The review also sheds light on the intriguing possibility that some ncRNAs may encode functional micropeptides. Therefore, the review underscores the emergent roles of ncRNAs as potent regulatory factors that significantly enrich the regulatory network governing plant growth, development, and responses to environmental stimuli. There are yet-to-be-discovered roles of ncRNAs waiting for us to explore.

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.

Long-read direct RNA sequencing reveals epigenetic regulation of chimeric gene-transposon transcripts in Arabidopsis thaliana

Transposable elements (TEs) are accumulated in both intergenic and intragenic regions in plant genomes. Intragenic TEs often act as regulatory elements of associated genes and are also co-transcribed with genes, generating chimeric TE-gene transcripts. Despite the potential impact on mRNA regulation and gene function, the prevalence and transcriptional regulation of TE-gene transcripts are poorly understood. By long-read direct RNA sequencing and a dedicated bioinformatics pipeline, ParasiTE, we investigated the transcription and RNA processing of TE-gene transcripts in Arabidopsis thaliana. We identified a global production of TE-gene transcripts in thousands of A. thaliana gene loci, with TE sequences often being associated with alternative transcription start sites or transcription termination sites. The epigenetic state of intragenic TEs affects RNAPII elongation and usage of alternative poly(A) signals within TE sequences, regulating alternative TE-gene isoform production. Co-transcription and inclusion of TE-derived sequences into gene transcripts impact regulation of RNA stability and environmental responses of some loci. Our study provides insights into TE-gene interactions that contributes to mRNA regulation, transcriptome diversity, and environmental responses in plants.

RNA degradome analysis reveals DNE1 endoribonuclease is required for the turnover of diverse mRNA substrates in Arabidopsis

In plants, cytoplasmic mRNA decay is critical for posttranscriptionally controlling gene expression and for maintaining cellular RNA homeostasis. Arabidopsis DCP1-ASSOCIATED NYN ENDORIBONUCLEASE 1 (DNE1) is a cytoplasmic mRNA decay factor that interacts with proteins involved in mRNA decapping and nonsense-mediated mRNA decay (NMD). There is limited information on the functional role of DNE1 in RNA turnover, and the identities of its endogenous targets are unknown. In this study, we utilized RNA degradome approaches to globally investigate DNE1 substrates. Monophosphorylated 5' ends, produced by DNE1, should accumulate in mutants lacking the cytoplasmic exoribonuclease XRN4, but be absent from DNE1 and XRN4 double mutants. In seedlings, we identified over 200 such transcripts, most of which reflect cleavage within coding regions. While most DNE1 targets were NMD-insensitive, some were upstream ORF (uORF)-containing and NMD-sensitive transcripts, indicating that this endoribonuclease is required for turnover of a diverse set of mRNAs. Transgenic plants expressing DNE1 cDNA with an active-site mutation in the endoribonuclease domain abolished the in planta cleavage of transcripts, demonstrating that DNE1 endoribonuclease activity is required for cleavage. Our work provides key insights into the identity of DNE1 substrates and enhances our understanding of DNE1-mediated mRNA decay.

Cotranscriptional RNA processing and modification in plants

The activities of RNA polymerases shape the epigenetic landscape of genomes with profound consequences for genome integrity and gene expression. A fundamental event during the regulation of eukaryotic gene expression is the coordination between transcription and RNA processing. Most primary RNAs mature through various RNA processing and modification events to become fully functional. While pioneering results positioned RNA maturation steps after transcription ends, the coupling between the maturation of diverse RNA species and their transcription is becoming increasingly evident in plants. In this review, we discuss recent advances in our understanding of the crosstalk between RNA Polymerase II, IV, and V transcription and nascent RNA processing of both coding and noncoding RNAs.

LncRNA FLAIL affects alternative splicing and represses flowering in Arabidopsis

How the noncoding genome affects cellular functions is a key biological question. A particular challenge is to distinguish the effects of noncoding DNA elements from long noncoding RNAs (lncRNAs) that coincide at the same loci. Here, we identified the flowering-associated intergenic lncRNA (FLAIL) in Arabidopsis through early flowering flail mutants. Expression of FLAIL RNA from a different chromosomal location in combination with strand-specific RNA knockdown characterized FLAIL as a trans-acting RNA molecule. FLAIL directly binds to differentially expressed target genes that control flowering via RNA-DNA interactions through conserved sequence motifs. FLAIL interacts with protein and RNA components of the spliceosome to affect target mRNA expression through co-transcriptional alternative splicing (AS) and linked chromatin regulation. In the absence of FLAIL, splicing defects at the direct FLAIL target flowering gene LACCASE 8 (LAC8) correlated with reduced mRNA expression. Double mutant analyses support a model where FLAIL-mediated splicing of LAC8 promotes its mRNA expression and represses flowering. Our study suggests lncRNAs as accessory components of the spliceosome that regulate AS and gene expression to impact organismal development.

Plant terminators: the unsung heroes of gene expression

To be properly expressed, genes need to be accompanied by a terminator, a region downstream of the coding sequence that contains the information necessary for the maturation of the mRNA 3' end. The main event in this process is the addition of a poly(A) tail at the 3' end of the new transcript, a critical step in mRNA biology that has important consequences for the expression of genes. Here, we review the mechanism leading to cleavage and polyadenylation of newly transcribed mRNAs and how this process can affect the final levels of gene expression. We give special attention to an aspect often overlooked, the effect that different terminators can have on the expression of genes. We also discuss some exciting findings connecting the choice of terminator to the biogenesis of small RNAs, which are a central part of one of the most important mechanisms of regulation of gene expression in plants.

TFIIS Is Crucial During Early Transcript Elongation for Transcriptional Reprogramming in Response to Heat Stress

In addition to the stage of transcriptional initiation, the production of mRNAs is regulated during elongation. Accordingly, the synthesis of mRNAs by RNA polymerase II (RNAPII) in the chromatin context is modulated by various transcript elongation factors. TFIIS is an elongation factor that stimulates the transcript cleavage activity of RNAPII to reactivate stalled elongation complexes at barriers to transcription including nucleosomes. Since Arabidopsis tfIIs mutants grow normally under standard conditions, we have exposed them to heat stress (HS), revealing that tfIIs plants are highly sensitive to elevated temperatures. Transcriptomic analyses demonstrate that particularly HS-induced genes are expressed at lower levels in tfIIs than in wildtype. Mapping the distribution of elongating RNAPII uncovered that in tfIIs plants RNAPII accumulates at the +1 nucleosome of genes that are upregulated upon HS. The promoter-proximal RNAPII accumulation in tfIIs under HS conditions conforms to that observed upon inhibition of the RNAPII transcript cleavage activity. Further analysis of the RNAPII accumulation downstream of transcriptional start sites illustrated that RNAPII stalling occurs at +1 nucleosomes that are depleted in the histone variant H2A.Z upon HS. Therefore, assistance of early transcript elongation by TFIIS is required for reprogramming gene expression to establish plant thermotolerance.

Intergenic splicing-stimulated transcriptional readthrough is suppressed by nonsense-mediated mRNA decay in Arabidopsis

Recent emerging evidence has shown that readthrough transcripts (RTs), including polycistronic mRNAs, are also transcribed in eukaryotes. However, the post-transcriptional regulation for these remains to be elucidated. Here, we identify 271 polycistronic RT-producing loci in Arabidopsis. Increased accumulation of RTs is detected in the nonsense-mediated mRNA decay (NMD)-deficient mutants compared with wild type, and the second open reading frames (ORFs) of bicistronic mRNAs are rarely translated in contrast to the first ORFs. Intergenic splicing (IS) events which occur between first and second genes are seen in 158 RTs. Splicing inhibition assays suggest that IS eliminates the chance of transcription termination at the polyadenylation sites of the first gene and promotes accumulation of RTs. These results indicate that RTs arise from genes whose transcription termination is relatively weak or attenuated by IS, but NMD selectively degrades them. Ultimately, this report presents a eukaryotic strategy for RNA metabolism.

Comprehensive Transcriptome Analysis Reveals Genome-Wide Changes Associated with Endoplasmic Reticulum (ER) Stress in Potato (Solanum tuberosum L.)

We treated potato (Solanum tuberosum L.) plantlets with TM and performed gene expression studies to identify genome-wide changes associated with endoplasmic reticulum (ER) stress and the unfolded protein response (UPR). An extensive network of responses was identified, including chromatin remodeling, transcriptional reprogramming, as well as changes in the structural components of the endomembrane network system. Limited genome-wide changes in alternative RNA splicing patterns of protein-coding transcripts were also discovered. Significant changes in RNA metabolism, components of the translation machinery, as well as factors involved in protein folding and maturation occurred, which included a broader set of genes than expected based on Arabidopsis research. Antioxidant defenses and oxygen metabolic enzymes are differentially regulated, which is expected of cells that may be experiencing oxidative stress or adapting to protect proteins from oxidation. Surges in protein kinase expression indicated early signal transduction events. This study shows early genomic responses including an array of differentially expressed genes that have not been reported in Arabidopsis. These data describe novel ER stress responses in a solanaceous host.

The TRIPLE PHD FINGERS proteins are required for SWI/SNF complex-mediated +1 nucleosome positioning and transcription start site determination in Arabidopsis

Eukaryotes have evolved multiple ATP-dependent chromatin remodelers to shape the nucleosome landscape. We recently uncovered an evolutionarily conserved SWItch/Sucrose Non-Fermentable (SWI/SNF) chromatin remodeler complex in plants reminiscent of the mammalian BAF subclass, which specifically incorporates the MINUSCULE (MINU) catalytic subunits and the TRIPLE PHD FINGERS (TPF) signature subunits. Here we report experimental evidence that establishes the functional relevance of TPF proteins for the complex activity. Our results show that depletion of TPF triggers similar pleiotropic phenotypes and molecular defects to those found in minu mutants. Moreover, we report the genomic location of MINU2 and TPF proteins as representative members of this SWI/SNF complex and their impact on nucleosome positioning and transcription. These analyses unravel the binding of the complex to thousands of genes where it modulates the position of the +1 nucleosome. These targets tend to produce 5′-shifted transcripts in the tpf and minu mutants pointing to the participation of the complex in alternative transcription start site usage. Interestingly, there is a remarkable correlation between +1 nucleosome shift and 5′ transcript length change suggesting their functional connection. In summary, this study unravels the function of a plant SWI/SNF complex involved in +1 nucleosome positioning and transcription start site determination.

Inositol hexakisphosphate is required for Integrator function

Integrator is a multi-subunit protein complex associated with RNA polymerase II (Pol II), with critical roles in noncoding RNA 3'-end processing and transcription attenuation of a broad collection of mRNAs. IntS11 is the endonuclease for RNA cleavage, as a part of the IntS4-IntS9-IntS11 Integrator cleavage module (ICM). Here we report a cryo-EM structure of the Drosophila ICM, at 2.74 Å resolution, revealing stable association of an inositol hexakisphosphate (IP6) molecule. The IP6 binding site is located in a highly electropositive pocket at an interface among all three subunits of ICM, 55 Å away from the IntS11 active site and generally conserved in other ICMs. We also confirmed IP6 association with the same site in human ICM. IP6 binding is not detected in ICM samples harboring mutations in this binding site. Such mutations or disruption of IP6 biosynthesis significantly reduced Integrator function in snRNA 3'-end processing and mRNA transcription attenuation. Our structural and functional studies reveal that IP6 is required for Integrator function in Drosophila, humans, and likely other organisms.

Transcription-coupled and epigenome-encoded mechanisms direct H3K4 methylation

Mono-, di-, and trimethylation of histone H3 lysine 4 (H3K4me1/2/3) are associated with transcription, yet it remains controversial whether H3K4me1/2/3 promote or result from transcription. Our previous characterizations of Arabidopsis H3K4 demethylases suggest roles for H3K4me1 in transcription. However, the control of H3K4me1 remains unexplored in Arabidopsis, in which no methyltransferase for H3K4me1 has been identified. Here, we identify three Arabidopsis methyltransferases that direct H3K4me1. Analyses of their genome-wide localization using ChIP-seq and machine learning reveal that one of the enzymes cooperates with the transcription machinery, while the other two are associated with specific histone modifications and DNA sequences. Importantly, these two types of localization patterns are also found for the other H3K4 methyltransferases in Arabidopsis and mice. These results suggest that H3K4me1/2/3 are established and maintained via interplay with transcription as well as inputs from other chromatin features, presumably enabling elaborate gene control.

PAMP-triggered genetic reprogramming involves widespread alternative transcription initiation and an immediate transcription factor wave

Immune responses triggered by pathogen-associated molecular patterns (PAMPs) are key to pathogen defense, but drivers and stabilizers of the growth-to-defense genetic reprogramming remain incompletely understood in plants. Here, we report a time-course study of the establishment of PAMP-triggered immunity (PTI) using cap analysis of gene expression. We show that around 15% of all transcription start sites (TSSs) rapidly induced during PTI define alternative transcription initiation events. From these, we identify clear examples of regulatory TSS change via alternative inclusion of target peptides or domains in encoded proteins, or of upstream open reading frames in mRNA leader sequences. We also find that 60% of PAMP response genes respond earlier than previously thought. In particular, a cluster of rapidly and transiently PAMP-induced genes is enriched in transcription factors (TFs) whose functions, previously associated with biological processes as diverse as abiotic stress adaptation and stem cell activity, appear to converge on growth restriction. Furthermore, examples of known potentiators of PTI, in one case under direct mitogen-activated protein kinase control, support the notion that the rapidly induced TFs could constitute direct links to PTI signaling pathways and drive gene expression changes underlying establishment of the immune state.

Catalytic activities, molecular connections, and biological functions of plant RNA exosome complexes

RNA exosome complexes provide the main 3'-5'-exoribonuclease activities in eukaryotic cells and contribute to the maturation and degradation of virtually all types of RNA. RNA exosomes consist of a conserved core complex that associates with exoribonucleases and with multimeric cofactors that recruit the enzyme to its RNA targets. Despite an overall high level of structural and functional conservation, the enzymatic activities and compositions of exosome complexes and their cofactor modules differ among eukaryotes. This review highlights unique features of plant exosome complexes, such as the phosphorolytic activity of the core complex, and discusses the exosome cofactors that operate in plants and are dedicated to the maturation of ribosomal RNA, the elimination of spurious, misprocessed, and superfluous transcripts, or the removal of mRNAs cleaved by the RNA-induced silencing complex and other mRNAs prone to undergo silencing.

Transcriptomic Complexity in Strawberry Fruit Development and Maturation Revealed by Nanopore Sequencing

The use of alternative transcription start or termination sites (aTSS or aTTS) as well as alternative splicing (AS) produce diverse transcript isoforms, playing indispensable roles in the plant development and environmental adaptations. Despite the advances in the finding of the genome-wide alternatively spliced genes in strawberry, it remains unexplored how AS responds to the developmental cues and what relevance do these outcomes have to the gene function. In this study, we have systematically investigated the transcriptome complexity using long-read Oxford Nanopore Technologies along the four successive developmental stages. The full-length cDNA sequencing results unraveled thousands of previously unexplored transcript isoforms raised from aTSS, aTTS, and AS. The relative contributions of these three processes to the complexity of strawberry fruit transcripts were compared. The aTSS and aTTS were more abundant than the AS. Differentially expressed transcripts unraveled the key transitional role of the white fruit stage. Isoform switches of transcripts from 757 genes were observed. They were associated with protein-coding potential change and domain gain or loss as the main consequences. Those genes with switched isoforms take part in the key processes of maturation in the late stages. A case study using yeast two hybrid analysis supported the functional divergence of the two isoforms of the B-box protein 22. Our results provided a new comprehensive overview of the dynamic transcriptomic landscape during strawberry fruit development and maturation.

Landscape of transcription termination in Arabidopsis revealed by single-molecule nascent RNA sequencing

BACKGROUND

The dynamic process of transcription termination produces transient RNA intermediates that are difficult to distinguish from each other via short-read sequencing methods.

RESULTS

Here, we use single-molecule nascent RNA sequencing to characterize the various forms of transient RNAs during termination at genome-wide scale in wildtype Arabidopsis and in atxrn3, fpa, and met1 mutants. Our data reveal a wide range of termination windows among genes, ranging from ~ 50 nt to over 1000 nt. We also observe efficient termination before downstream tRNA genes, suggesting that chromatin structure around the promoter region of tRNA genes may block pol II elongation. 5' Cleaved readthrough transcription in atxrn3 with delayed termination can run into downstream genes to produce normally spliced and polyadenylated mRNAs in the absence of their own transcription initiation. Consistent with previous reports, we also observe long chimeric transcripts with cryptic splicing in fpa mutant; but loss of CG DNA methylation has no obvious impact on termination in the met1 mutant.

CONCLUSIONS

Our method is applicable to establish a comprehensive termination landscape in a broad range of species.

Cryptic promoter activation occurs by at least two different mechanisms in the Arabidopsis genome

In gene-trap screening of plant genomes, promoterless reporter constructs are often expressed without trapping of annotated gene promoters. The molecular basis of this phenomenon, which has been interpreted as the trapping of cryptic promoters, is poorly understood. Here, we found that cryptic promoter activation occurs by at least two different mechanisms using Arabidopsis gene-trap lines in which a firefly luciferase (LUC) open reading frame (ORF) without an apparent promoter sequence was expressed from intergenic regions: one mechanism is 'cryptic promoter capturing', in which the LUC ORF captured pre-existing promoter-like chromatin marked by H3K4me3 and H2A.Z, and the other is 'promoter de novo origination', in which the promoter chromatin was newly formed near the 5' end of the inserted LUC ORF. The latter finding raises a question as to how the inserted LUC ORF sequence is involved in this phenomenon. To examine this, we performed a model experiment with chimeric LUC genes in transgenic plants. Using Arabidopsis psaH1 promoter-LUC constructs, we found that the functional core promoter region, where transcription start sites (TSSs) occur, cannot simply be determined by the upstream nor core promoter sequences; rather, its positioning proximal to the inserted LUC ORF sequence was more critical. This result suggests that the insertion of the coding sequence alters the local distribution of TSSs in the plant genome. The possible impact of the two types of cryptic promoter activation mechanisms on plant genome evolution and endosymbiotic gene transfer is discussed.

G-Quadruplex in Gene Encoding Large Subunit of Plant RNA Polymerase II: A Billion-Year-Old Story

G-quadruplexes have long been perceived as rare and physiologically unimportant nucleic acid structures. However, several studies have revealed their importance in molecular processes, suggesting their possible role in replication and gene expression regulation. Pathways involving G-quadruplexes are intensively studied, especially in the context of human diseases, while their involvement in gene expression regulation in plants remains largely unexplored. Here, we conducted a bioinformatic study and performed a complex circular dichroism measurement to identify a stable G-quadruplex in the gene RPB1, coding for the RNA polymerase II large subunit. We found that this G-quadruplex-forming locus is highly evolutionarily conserved amongst plants sensu lato (Archaeplastida) that share a common ancestor more than one billion years old. Finally, we discussed a new hypothesis regarding G-quadruplexes interacting with UV light in plants to potentially form an additional layer of the regulatory network.

Long Noncoding RNAs in Plants

Plants have an extraordinary diversity of transcription machineries, including five nuclear DNA-dependent RNA polymerases. Four of these enzymes are dedicated to the production of long noncoding RNAs (lncRNAs), which are ribonucleic acids with functions independent of their protein-coding potential. lncRNAs display a broad range of lengths and structures, but they are distinct from the small RNA guides of RNA interference (RNAi) pathways. lncRNAs frequently serve as structural, catalytic, or regulatory molecules for gene expression. They can affect all elements of genes, including promoters, untranslated regions, exons, introns, and terminators, controlling gene expression at various levels, including modifying chromatin accessibility, transcription, splicing, and translation. Certain lncRNAs protect genome integrity, while others respond to environmental cues like temperature, drought, nutrients, and pathogens. In this review, we explain the challenge of defining lncRNAs, introduce the machineries responsible for their production, and organize this knowledge by viewing the functions of lncRNAs throughout the structure of a typical plant gene.

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.

Winter fields antisense RNAs to kick off flowering

This Outlook discusses the finding by Zhao et al. demonstrating that COOLAIR-mediated FLC silencing is induced by the first seasonal frost in the field, and this acts as a key molecular indicator during autumn for winter arrival.

FLOWERING LOCUS C (FLC), a MADS-box transcription factor, plays a major role in determining flowering time in Arabidopsis. In this issue of Genes & Development, Zhao and colleagues (pp. 888–898) elucidate the role of COOLAIR antisense noncoding RNAs in FLC regulation through field trials and laboratory experiments. COOLAIR-mediated FLC silencing is induced by the first seasonal frost in the field and thus acts as a key molecular indicator during autumn for winter arrival.

TrancriptomeReconstructoR: data-driven annotation of complex transcriptomes

BACKGROUND

The quality of gene annotation determines the interpretation of results obtained in transcriptomic studies. The growing number of genome sequence information calls for experimental and computational pipelines for de novo transcriptome annotation. Ideally, gene and transcript models should be called from a limited set of key experimental data.

RESULTS

We developed TranscriptomeReconstructoR, an R package which implements a pipeline for automated transcriptome annotation. It relies on integrating features from independent and complementary datasets: (i) full-length RNA-seq for detection of splicing patterns and (ii) high-throughput 5' and 3' tag sequencing data for accurate definition of gene borders. The pipeline can also take a nascent RNA-seq dataset to supplement the called gene model with transient transcripts. We reconstructed de novo the transcriptional landscape of wild type Arabidopsis thaliana seedlings and Saccharomyces cerevisiae cells as a proof-of-principle. A comparison to the existing transcriptome annotations revealed that our gene model is more accurate and comprehensive than the most commonly used community gene models, TAIR10 and Araport11 for A.thaliana and SacCer3 for S.cerevisiae. In particular, we identify multiple transient transcripts missing from the existing annotations. Our new annotations promise to improve the quality of A.thaliana and S.cerevisiae genome research.

CONCLUSIONS

Our proof-of-concept data suggest a cost-efficient strategy for rapid and accurate annotation of complex eukaryotic transcriptomes. We combine the choice of library preparation methods and sequencing platforms with the dedicated computational pipeline implemented in the TranscriptomeReconstructoR package. The pipeline only requires prior knowledge on the reference genomic DNA sequence, but not the transcriptome. The package seamlessly integrates with Bioconductor packages for downstream analysis.

Using TIF-Seq2 to investigate association between 5´ and 3´mRNA ends

The development of high-throughput technologies has revealed pervasive transcription in all genomes that have been investigated so far. This has uncovered a highly interleaved transcriptome organization involving thousands of overlapping coding and non-coding RNA isoforms that challenge our traditional definitions of genes and functional regions of the genome. In this chapter, we discuss the application of an improved Transcript Isoform Sequencing approach (TIF-Seq2) able to concurrently determine the start and end sites of individual RNA molecules. We exemplify its use for the investigation of the human transcriptome and show how it is especially well suited to discriminate between overlapping molecules and accurately define their boundaries.

Natural temperature fluctuations promote COOLAIR regulation of FLC

In this study, Zhao et al. set out to characterize how plants respond to cold through regulation of FLC expression. Using genetics and genomics approaches, the authors reveal how natural temperature fluctuations promote COOLAIR regulation of FLC, with the first autumn frost acting as a key indicator of autumn/winter arrival.

Plants monitor many aspects of their fluctuating environments to help align their development with seasons. Molecular understanding of how noisy temperature cues are registered has emerged from dissection of vernalization in Arabidopsis, which involves a multiphase cold-dependent silencing of the floral repressor locus FLOWERING LOCUS C (FLC). Cold-induced transcriptional silencing precedes a low probability PRC2 epigenetic switching mechanism. The epigenetic switch requires the absence of warm temperatures as well as long-term cold exposure. However, the natural temperature inputs into the earlier transcriptional silencing phase are less well understood. Here, through investigation of Arabidopsis accessions in natural and climatically distinct field sites, we show that the first seasonal frost strongly induces expression of COOLAIR, the antisense transcripts at FLC. Chamber experiments delivering a constant mean temperature with different fluctuations showed the freezing induction of COOLAIR correlates with stronger repression of FLC mRNA. Identification of a mutant that ectopically activates COOLAIR revealed how COOLAIR up-regulation can directly reduce FLC expression. Consistent with this, transgenes designed to knockout COOLAIR perturbed the early phase of FLC silencing. However, all transgenes designed to remove COOLAIR resulted in increased production of novel convergent FLC antisense transcripts. Our study reveals how natural temperature fluctuations promote COOLAIR regulation of FLC, with the first autumn frost acting as a key indicator of autumn/winter arrival.

Chromatin-associated SUMOylation controls the transcriptional switch between plant development and heat stress responses

The post-translational protein modification known as SUMOylation has conserved roles in the heat stress responses of various species. The functional connection between the global regulation of gene expression and chromatin-associated SUMOylation in plant cells is unknown. Here, we uncovered a genome-wide relationship between chromatin-associated SUMOylation and transcriptional switches in Arabidopsis thaliana grown at room temperature, exposed to heat stress, and exposed to heat stress followed by recovery. The small ubiquitin-like modifier (SUMO)-associated chromatin sites, characterized by whole-genome ChIP-seq, were generally associated with active chromatin markers. In response to heat stress, chromatin-associated SUMO signals increased at promoter-transcriptional start site regions and decreased in gene bodies. RNA-seq analysis supported the role of chromatin-associated SUMOylation in transcriptional activation during rapid responses to high temperature. Changes in SUMO signals on chromatin were associated with the upregulation of heat-responsive genes and the downregulation of growth-related genes. Disruption of the SUMO ligase gene SIZ1 abolished SUMO signals on chromatin and attenuated rapid transcriptional responses to heat stress. The SUMO signal peaks were enriched in DNA elements recognized by distinct groups of transcription factors under different temperature conditions. These observations provide evidence that chromatin-associated SUMOylation regulates the transcriptional switch between development and heat stress response in plant cells.

Terminus-Associated Non-coding RNAs: Trash or Treasure?

3′ untranslated regions (3′ UTRs) of protein-coding genes are well known for their important roles in determining the fate of mRNAs in diverse processes, including trafficking, stabilization, translation, and RNA–protein interactions. However, non-coding RNAs (ncRNAs) scattered around 3′ termini of the protein-coding genes, here referred to as terminus-associated non-coding RNAs (TANRs), have not attracted wide attention in RNA research. Indeed, whether TANRs are transcriptional noise, degraded mRNA products, alternative 3′ UTRs, or functional molecules has remained unclear for a long time. As a new category of ncRNAs, TANRs are widespread, abundant, and conserved in diverse eukaryotes. The biogenesis of TANRs mainly follows the same promoter model, the RNA-dependent RNA polymerase activity-dependent model, or the independent promoter model. Functional studies of TANRs suggested that they are significantly involved in the versatile regulation of gene expression. For instance, at the transcriptional level, they can lead to transcriptional interference, induce the formation of gene loops, and participate in transcriptional termination. Furthermore, at the posttranscriptional level, they can act as microRNA sponges, and guide cleavage or modification of target RNAs. Here, we review current knowledge of the potential role of TANRs in the modulation of gene expression. In this review, we comprehensively summarize the current state of knowledge about TANRs, and discuss TANR nomenclature, relation to ncRNAs, cross-talk biogenesis pathways and potential functions. We further outline directions of future studies of TANRs, to promote investigations of this emerging and enigmatic category of RNA.