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.

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.

AbstractEukaryotic genomes give rise to thousands of long non-coding RNAs (lncRNAs), yet the purpose of lncRNAs remains largely enigmatic. Functional characterization of lncRNAs is challenging due to multiple orthogonal hypothesis for molecular activities of lncRNA loci. Here, we identified afloweringassociatedintergeniclncRNA (FLAIL) that represses flowering inArabidopsis. An allelic series offlailloss-of-function mutants generated by CRISPR/Cas9 and T-DNA mutagenesis showed an early flowering phenotype. Gene expression analyses inflailmutants revealed differentially expressed genes linked to the regulation of flowering. A genomic rescue fragment ofFLAILintroduced inflailmutants complemented gene expression defects and early flowering, consistent withtrans-acting effects of theFLAILRNA. Knock-down ofFLAILRNA levels using the artificial microRNA approach revealed an early flowering phenotype shared with genomic mutations, indicating atrans-acting role ofFLAILRNA in the repression of flowering time. Genome-wide detection ofFLAIL-DNA interactions by ChIRP-seq suggested thatFLAILmay directly bind genomic regions.FLAILbound to genes involved in regulation of flowering that were differentially expressed inflail, consistent with the interpretation ofFLAILas atrans-acting lncRNA directly shaping gene expression. Our findings highlightFLAILas atrans-acting lncRNA that affects flowering inArabidopsis, likely through mediating transcriptional regulation of genes directly bound byFLAIL.

Nucleosome-depleted regions (NDRs) at gene promoters support initiation of RNA polymerase II transcription. Interestingly, transcription often initiates in both directions, resulting in an mRNA and a divergent non-coding (DNC) transcript of unclear purpose. Here, we characterized the genetic architecture and molecular mechanism of DNC transcription in budding yeast. Using high-throughput reverse genetic screens based on quantitative single-cell fluorescence measurements, we identified the Hda1 histone deacetylase complex (Hda1C) as a repressor of DNC transcription. Nascent transcription profiling showed a genome-wide role of Hda1C in repression of DNC transcription. Live-cell imaging of transcription revealed that mutations in the Hda3 subunit increased the frequency of DNC transcription. Hda1C contributed to decreased acetylation of histone H3 in DNC transcription regions, supporting DNC transcription repression by histone deacetylation. Our data support the interpretation that DNC transcription results as a consequence of the NDR-based architecture of eukaryotic promoters, but that it is governed by locus-specific repression to maintain genome fidelity.

Nucleosome-depleted regions (NDRs) at gene promoters support initiation of RNA Polymerase II transcription. Interestingly, transcription often initiates in both directions, resulting in an mRNA, and a divergent non-coding (DNC) transcript with an unclear purpose. Here, we characterized the genetic architecture and molecular mechanism of DNC transcription in budding yeast. We identified the Hda1 histone deacetylase complex (Hda1C) as a repressor of DNC in high-throughput reverse genetic screens based on quantitative single-cell fluorescence measurements. Nascent transcription profiling showed a genome-wide role of Hda1C in DNC repression. Live-cell imaging of transcription revealed that Hda1C reduced the frequency of DNC transcription. Hda1C contributed to decreased acetylation of histone H3 in DNC regions, supporting DNC repression by histone deacetylation. Our data support the interpretation that DNC results as a consequence of the NDR-based architecture of eukaryotic promoters, but that it is governed by locus-specific repression to maintain genome fidelity.

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 as a proof-of-principle. A comparison to the existing transcriptome annotations revealed that our gene model is more accurate and comprehensive than the two most commonly used community gene models, TAIR10 and Araport11. In particular, we identify thousands of transient transcripts missing from the existing annotations. Our new annotation promises to improve the quality of A.thaliana 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.

Eukaryotic genomes give rise to thousands of long non-coding RNAs (lncRNAs), yet the purpose of lncRNAs remains largely enigmatic. Functional characterization of lncRNAs is challenging due to multiple orthogonal hypothesis for molecular activities of lncRNA loci. Here, we identified a fl owering a ssociated i ntergenic l ncRNA ( FLAIL ) that represses flowering in Arabidopsis . An allelic series of flail loss-of-function mutants generated by CRISPR/Cas9 and T-DNA mutagenesis showed an early flowering phenotype. Gene expression analyses in flail mutants revealed differentially expressed genes linked to the regulation of flowering. A genomic rescue fragment of FLAIL introduced in flail mutants complemented gene expression defects and early flowering, consistent with trans- acting effects of the FLAIL RNA. Knock-down of FLAIL RNA levels using the artificial microRNA approach revealed an early flowering phenotype shared with genomic mutations, indicating a trans- acting role of FLAIL RNA in the repression of flowering time. Genome-wide detection of FLAIL -DNA interactions by ChIRP-seq suggested that FLAIL may directly bind genomic regions. FLAIL bound to genes involved in regulation of flowering that were differentially expressed in flail , consistent with the interpretation of FLAIL as a trans- acting lncRNA directly shaping gene expression. Our findings highlight FLAIL as a trans- acting lncRNA that affects flowering in Arabidopsis , likely through mediating transcriptional regulation of genes directly bound by FLAIL .### Competing Interest StatementThe authors have declared no competing interest.* amiRNAs : artificial microRNAs ChIRP : chromatin isolation by RNA purification CIR1 : CIRCADIAN 1 CNIT : Coding-NonCoding Identifying Tool CO : CONSTANS CPC2 : Coding Potential Calculator DEGs : differentially expressed genes FLAIL : flowering associated intergenic lncRNA FLC : FLOWERING LOCUS C LAC8 : LACCASE 8 lncRNAs : long non-coding RNAs LUC : Luciferase NMT3 : PHOSPHOETHANOLAMINE N-METHYLTRANSFERASE 3 sORFs : small open reading frames PCP : PORCUPINE PIF4 : PHYTOCHROME INTERACTING FACTOR 4 plaNET-seq : plant Native Elongating Transcripts sequencing PLL25 : PECTIN LYASE-LIKE 25 RNAPII : RNA polymerase II SAMs : shoot apical meristems SME : Sm protein E1 TIF-seq : Transcript Isoform sequencing TSS-seq : Transcription Start Site sequencing UBQ : UBIQUITIN WMD : Web MicroRNA Designer

ABSTRACT In animals, transcription by RNA polymerase II initiates bidirectionally from gene promoters to produce pre-mRNAs on the forward strand and promoter upstream transcripts (PROMPTs) on the reverse strand. PROMPTs are rapidly degraded by the nuclear exosome. Similarly, active enhancer regions in animals initiate transcription of exosome-sensitive enhancer RNAs (eRNAs). Previous studies based on nascent RNA approaches concluded that Arabidopsis thaliana does not produce PROMPTs. Here, we used steady-state RNA sequencing methods in mutants defective in nuclear RNA decay, including by the exosome, to reassess the existence of PROMPTs and eRNAs in A. thaliana . While PROMPTs are overall rare in A. thaliana , about 100 clear cases of exosome-sensitive PROMPTs and 113 loci producing eRNA-like transcripts were identified. In addition, we found ∼200 transcription start sites within 3’-UTR-encoding regions that produce unspliced exosome-sensitive antisense RNAs covering much of the cognate pre-mRNA. A typical representative of this class of RNAs is the previously characterized non-coding RNA controlling the expression of the key seed dormancy regulator, DELAY OF GERMINATION1 . Exosome-sensitive antisense RNAs are overrepresented in transcription factor genes, suggesting a potential for widespread control of gene expression. Lastly, we assess the use of alternative promoters in A. thaliana and compare the accuracy of existing TSS annotations.

RNA polymerase II (RNAPII) transcription is crucial for gene expression. RNAPII density peaks at gene boundaries, associating these key regions for gene expression control with limited RNAPII movement. The connections between RNAPII transcription speed and gene regulation in multicellular organisms are poorly understood. Here, we directly modulate RNAPII transcription speed by point mutations in the second largest subunit of RNAPII in Arabidopsis thaliana. A RNAPII mutation predicted to decelerate transcription is inviable, while accelerating RNAPII transcription confers phenotypes resembling auto-immunity. Nascent transcription profiling revealed that RNAPII complexes with accelerated transcription clear stalling sites at both gene ends, resulting in read-through transcription. The accelerated transcription mutant NRPB2-Y732F exhibits increased association with 5' splice site (5'SS) intermediates and enhanced splicing efficiency. Our findings highlight potential advantages of RNAPII stalling through local reduction in transcription speed to optimize gene expression for the development of multicellular organisms.

Progression of RNA polymerase II (RNAPII) transcription relies on the appropriately positioned activities of elongation factors. The resulting profile of factors and chromatin signatures along transcription units provides a "positional information system" for transcribing RNAPII. Here, we investigate a chromatin-based mechanism that suppresses intragenic initiation of RNAPII transcription. We demonstrate that RNAPII transcription across gene promoters represses their function in plants. This repression is characterized by reduced promoter-specific molecular signatures and increased molecular signatures associated with RNAPII elongation. The conserved FACT histone chaperone complex is required for this repression mechanism. Genome-wide Transcription Start Site (TSS) mapping reveals thousands of discrete intragenic TSS positions in fact mutants, including downstream promoters that initiate alternative transcript isoforms. We find that histone H3 lysine 4 mono-methylation (H3K4me1), an Arabidopsis RNAPII elongation signature, is enriched at FACT-repressed intragenic TSSs. Our analyses suggest that FACT is required to repress intragenic TSSs at positions that are in part characterized by elevated H3K4me1 levels. In sum, conserved and plant-specific chromatin features correlate with the co-transcriptional repression of intragenic TSSs. Our insights into TSS repression by RNAPII transcription promise to inform the regulation of alternative transcript isoforms and the characterization of gene regulation through the act of pervasive transcription across eukaryotic genomes.

Temperature profoundly affects the kinetics of biochemical reactions, yet how large molecular complexes such as the transcription machinery accommodate changing temperatures to maintain cellular function is poorly understood. Here, we developed plant native elongating transcripts sequencing (plaNET-seq) to profile genome-wide nascent RNA polymerase II (RNAPII) transcription during the cold-response of Arabidopsis thaliana with single-nucleotide resolution. Combined with temporal resolution, these data revealed transient genome-wide reprogramming of nascent RNAPII transcription during cold, including characteristics of RNAPII elongation and thousands of non-coding transcripts connected to gene expression. Our results suggest a role for promoter-proximal RNAPII stalling in predisposing genes for transcriptional activation during plant-environment interactions. At gene 3’-ends, cold initially facilitated transcriptional termination by limiting the distance of read-through transcription. Within gene bodies, cold reduced the kinetics of co-transcriptional splicing leading to increased intragenic stalling. Our data resolved multiple distinct mechanisms by which temperature transiently altered the dynamics of nascent RNAPII transcription and associated RNA processing, illustrating potential biotechnological solutions and future focus areas to promote food security in the context of a changing climate.

Most DNA in the genomes of higher organisms does not code for proteins. RNA Polymerase II (Pol II) transcribes non-coding DNA into long non-coding RNAs (lncRNAs), but biological roles of lncRNA are unclear. We find that mutations in the yeast lncRNA CUT60 result in poor growth. Defective termination of CUT60 transcription causes read-through transcription across the ATP16 gene promoter. Read-through transcription localizes chromatin signatures associated with Pol II elongation to the ATP16 promoter. The act of Pol II elongation across this promoter represses functional ATP16 expression by a Transcriptional Interference (TI) mechanism. Atp16p function in the mitochondrial ATP-synthase complex promotes mitochondrial DNA stability. ATP16 repression by TI through inefficient termination of CUT60 therefore triggers mitochondrial genome loss. Our results expand the functional and mechanistic implications of non-coding DNA in eukaryotes by highlighting termination of nuclear lncRNA transcription as mechanism to stabilize an organellar genome.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Many studies assess epigenetic signatures in bodily fluids and make conclusions about epigenetic patterns in inaccessible tissues of interest such as brain, thus neglecting the tissue specificity of epigenetic marks.