Quantification of co-transcriptional splicing from RNA-Seq data

Truncating the spliceosomal 'rope protein' Prp45 results in Htz1 dependent phenotypes

Spliceosome assembly contributes an important but incompletely understood aspect of splicing regulation. Prp45 is a yeast splicing factor which runs as an extended fold through the spliceosome, and which may be important for bringing its components together. We performed a whole genome analysis of the genetic interaction network of the truncated allele of PRP45 (prp45(1-169)) using synthetic genetic array technology and found chromatin remodellers and modifiers as an enriched category. In agreement with related studies, H2A.Z-encoding HTZ1, and the components of SWR1, INO80, and SAGA complexes represented prominent interactors, with htz1 conferring the strongest growth defect. Because the truncation of Prp45 disproportionately affected low copy number transcripts of intron-containing genes, we prepared strains carrying intronless versions of SRB2, VPS75, or HRB1, the most affected cases with transcription-related function. Intron removal from SRB2, but not from the other genes, partly repaired some but not all the growth phenotypes identified in the genetic screen. The interaction of prp45(1-169) and htz1Δ was detectable even in cells with SRB2 intron deleted (srb2Δi). The less truncated variant, prp45(1-330), had a synthetic growth defect with htz1Δ at 16°C, which also persisted in the srb2Δi background. Moreover, htz1Δ enhanced prp45(1-330) dependent pre-mRNA hyper-accumulation of both high and low efficiency splicers, genes ECM33 and COF1, respectively. We conclude that while the expression defects of low expression intron-containing genes contribute to the genetic interactome of prp45(1-169), the genetic interactions between prp45 and htz1 alleles demonstrate the sensitivity of spliceosome assembly, delayed in prp45(1-169), to the chromatin environment.

Insplico: effective computational tool for studying splicing order of adjacent introns genome-wide with short and long RNA-seq reads

Although splicing occurs largely co-transcriptionally, the order by which introns are removed does not necessarily follow the order in which they are transcribed. Whereas several genomic features are known to influence whether or not an intron is spliced before its downstream neighbor, multiple questions related to adjacent introns' splicing order (AISO) remain unanswered. Here, we present Insplico, the first standalone software for quantifying AISO that works with both short and long read sequencing technologies. We first demonstrate its applicability and effectiveness using simulated reads and by recapitulating previously reported AISO patterns, which unveiled overlooked biases associated with long read sequencing. We next show that AISO around individual exons is remarkably constant across cell and tissue types and even upon major spliceosomal disruption, and it is evolutionarily conserved between human and mouse brains. We also establish a set of universal features associated with AISO patterns across various animal and plant species. Finally, we used Insplico to investigate AISO in the context of tissue-specific exons, particularly focusing on SRRM4-dependent microexons. We found that the majority of such microexons have non-canonical AISO, in which the downstream intron is spliced first, and we suggest two potential modes of SRRM4 regulation of microexons related to their AISO and various splicing-related features. Insplico is available on gitlab.com/aghr/insplico.

Profiling lariat intermediates reveals genetic determinants of early and late co-transcriptional splicing

Long introns with short exons in vertebrate genes are thought to require spliceosome assembly across exons (exon definition), rather than introns, thereby requiring transcription of an exon to splice an upstream intron. Here, we developed CoLa-seq (co-transcriptional lariat sequencing) to investigate the timing and determinants of co-transcriptional splicing genome wide. Unexpectedly, 90% of all introns, including long introns, can splice before transcription of a downstream exon, indicating that exon definition is not obligatory for most human introns. Still, splicing timing varies dramatically across introns, and various genetic elements determine this variation. Strong U2AF2 binding to the polypyrimidine tract predicts early splicing, explaining exon definition-independent splicing. Together, our findings question the essentiality of exon definition and reveal features beyond intron and exon length that are determinative for splicing timing.

A simple and robust method for isolating and analyzing chromatin-bound RNAs in Arabidopsis

BACKGROUND

Chromatin-bound RNAs are the primary product of transcription that undergo on-chromatin processing such as capping, splicing, and polyadenylation. These processing steps then determine the fate of RNAs. Albeit its vital importance, a simple and robust method for isolating different fractions of chromatin-bound RNAs is missing in plants.

RESULT

Here, we describe our updated method and the associated step-by-step protocol for chromatin-bound RNAs isolation in A. thaliana. The chromatin-bound RNAs isolation is based on the 1 M UREA wash that removes the majority of non-chromatin-associated proteins from the nucleus, as previously developed in mammalian cells. On-demand, the isolated chromatin-bound RNAs can be either used directly for gene-specific analysis or subject to further rRNA removal and also the optional polyadenylated RNA removal, followed by high-throughput sequencing. Detailed protocols for these procedures are also provided. Comparison of sequencing results of chromatin-bound RNAs with and without polyadenylated RNA removal revealed that a small fraction of CB-RNAs is polyadenylated but not yet fully spliced, representing RNA-processing intermediate on-chromatin.

CONCLUSION

This optimized chromatin-bound RNAs purification method is simple and robust and can be used to study transcription and its-coupled RNA processing in plants.

Tracking pre-mRNA maturation across subcellular compartments identifies developmental gene regulation through intron retention and nuclear anchoring

Steps of mRNA maturation are important gene regulatory events that occur in distinct cellular locations. However, transcriptomic analyses often lose information on the subcellular distribution of processed and unprocessed transcripts. We generated extensive RNA-seq data sets to track mRNA maturation across subcellular locations in mouse embryonic stem cells, neuronal progenitor cells, and postmitotic neurons. We find disparate patterns of RNA enrichment between the cytoplasmic, nucleoplasmic, and chromatin fractions, with some genes maintaining more polyadenylated RNA in chromatin than in the cytoplasm. We bioinformatically defined four regulatory groups for intron retention, including complete cotranscriptional splicing, complete intron retention in the cytoplasmic RNA, and two intron groups present in nuclear and chromatin transcripts but fully excised in cytoplasm. We found that introns switch their regulatory group between cell types, including neuronally excised introns repressed by polypyrimidine track binding protein 1 (PTBP1). Transcripts for the neuronal gamma-aminobutyric acid (GABA) B receptor, 1 (Gabbr1) are highly expressed in mESCs but are absent from the cytoplasm. Instead, incompletely spliced Gabbr1 RNA remains sequestered on chromatin, where it is bound by PTBP1, similar to certain long noncoding RNAs. Upon neuronal differentiation, Gabbr1 RNA becomes fully processed and exported for translation. Thus, splicing repression and chromatin anchoring of RNA combine to allow posttranscriptional regulation of Gabbr1 over development. For this and other genes, polyadenylated RNA abundance does not indicate functional gene expression. Our data sets provide a rich resource for analyzing many other aspects of mRNA maturation in subcellular locations and across development.

Long-Read cDNA Sequencing Enables a "Gene-Like" Transcript Annotation of Transposable Elements

Transcript-based annotations of genes facilitate both genome-wide analyses and detailed single-locus research. In contrast, transposable element (TE) annotations are rudimentary, consisting of information only on TE location and type. The repetitiveness and limited annotation of TEs prevent the ability to distinguish between potentially functional expressed elements and degraded copies. To improve genome-wide TE bioinformatics, we performed long-read sequencing of cDNAs from Arabidopsis (Arabidopsis thaliana) lines deficient in multiple layers of TE repression. These uniquely mapping transcripts were used to identify the set of TEs able to generate polyadenylated RNAs and create a new transcript-based annotation of TEs that we have layered upon the existing high-quality community standard annotation. We used this annotation to reduce the bioinformatic complexity associated with multimapping reads from short-read RNA sequencing experiments, and we show that this improvement is expanded in a TE-rich genome such as maize (Zea mays). Our TE annotation also enables the testing of specific standing hypotheses in the TE field. We demonstrate that inaccurate TE splicing does not trigger small RNA production, and the cell more strongly targets DNA methylation to TEs that have the potential to make mRNAs. This work provides a transcript-based TE annotation for Arabidopsis and maize, which serves as a blueprint to reduce the bioinformatic complexity associated with repetitive TEs in any organism.

The changing paradigm of intron retention: regulation, ramifications and recipes

Intron retention (IR) is a form of alternative splicing that has long been neglected in mammalian systems although it has been studied for decades in non-mammalian species such as plants, fungi, insects and viruses. It was generally assumed that mis-splicing, leading to the retention of introns, would have no physiological consequence other than reducing gene expression by nonsense-mediated decay. Relatively recent landmark discoveries have highlighted the pivotal role that IR serves in normal and disease-related human biology. Significant technical hurdles have been overcome, thereby enabling the robust detection and quantification of IR. Still, relatively little is known about the cis- and trans-acting modulators controlling this phenomenon. The fate of an intron to be, or not to be, retained in the mature transcript is the direct result of the influence exerted by numerous intrinsic and extrinsic factors at multiple levels of regulation. These factors have altered current biological paradigms and provided unexpected insights into the transcriptional landscape. In this review, we discuss the regulators of IR and methods to identify them. Our focus is primarily on mammals, however, we broaden the scope to non-mammalian organisms in which IR has been shown to be biologically relevant.

Gallop Racing Shifts Mature mRNA towards Introns: Does Exercise-Induced Stress Enhance Genome Plasticity?

Physical exercise is universally recognized as stressful. Among the "sport species", the horse is probably the most appropriate model for investigating the genomic response to stress due to the homogeneity of its genetic background. The aim of this work is to dissect the whole transcription modulation in Peripheral Blood Mononuclear Cells (PBMCs) after exercise with a time course framework focusing on unexplored regions related to introns and intergenic portions. PBMCs NGS from five 3 year old Sardinian Anglo-Arab racehorses collected at rest and after a 2000 m race was performed. Apart from differential gene expression ascertainment between the two time points the complexity of transcription for alternative transcripts was identified. Interestingly, we noted a transcription shift from the coding to the non-coding regions. We further investigated the possible causes of this phenomenon focusing on genomic repeats, using a differential expression approach and finding a strong general up-regulation of repetitive elements such as LINE. Since their modulation is also associated with the "exonization", the recruitment of repeats that act with regulatory functions, suggesting that there might be an active regulation of this transcriptional shift. Thanks to an innovative bioinformatic approach, our study could represent a model for the transcriptomic investigation of stress.

The Features and Regulation of Co-transcriptional Splicing in Arabidopsis

Precursor mRNA (pre-mRNA) splicing is essential for gene expression in most eukaryotic organisms. Previous studies from mammals, Drosophila, and yeast show that the majority of splicing events occurs co-transcriptionally. In plants, however, the features of co-transcriptional splicing (CTS) and its regulation still remain largely unknown. Here, we used chromatin-bound RNA sequencing to study CTS in Arabidopsis thaliana. We found that CTS is widespread in Arabidopsis seedlings, with a large proportion of alternative splicing events determined co-transcriptionally. CTS efficiency correlated with gene expression level, the chromatin landscape and, most surprisingly, the number of introns and exons of individual genes, but is independent of gene length. In combination with enhanced crosslinking and immunoprecipitation sequencing analysis, we further showed that the hnRNP-like proteins RZ-1B and RZ-1C promote efficient CTS globally through direct binding, frequently to exonic sequences. Notably, this general effect of RZ-1B/1C on splicing promotion is mainly observed at the chromatin level, not at the mRNA level. RZ-1C promotes CTS of multiple-exon genes in association with its binding to regions both proximal and distal to the regulated introns. We propose that RZ-1C promotes efficient CTS of genes with multiple exons through cooperative interactions with many exons, introns, and splicing factors. Our work thus reveals important features of CTS in plants and provides methodologies for the investigation of CTS and RNA-binding proteins in plants.

MeCP2 Represses the Rate of Transcriptional Initiation of Highly Methylated Long Genes

Mutations in the methyl-DNA-binding repressor protein MeCP2 cause the devastating neurodevelopmental disorder Rett syndrome. It has been challenging to understand how MeCP2 regulates transcription because MeCP2 binds broadly across the genome and MeCP2 mutations are associated with widespread small-magnitude changes in neuronal gene expression. We demonstrate here that MeCP2 represses nascent RNA transcription of highly methylated long genes in the brain through its interaction with the NCoR co-repressor complex. By measuring the rates of transcriptional initiation and elongation directly in the brain, we find that MeCP2 has no measurable effect on transcriptional elongation, but instead represses the rate at which Pol II initiates transcription of highly methylated long genes. These findings suggest a new model of MeCP2 function in which MeCP2 binds broadly across highly methylated regions of DNA, but acts at transcription start sites to attenuate transcriptional initiation.

Long-read sequencing of nascent RNA reveals coupling among RNA processing events

Pre-mRNA splicing is accomplished by the spliceosome, a megadalton complex that assembles de novo on each intron. Because spliceosome assembly and catalysis occur cotranscriptionally, we hypothesized that introns are removed in the order of their transcription in genomes dominated by constitutive splicing. Remarkably little is known about splicing order and the regulatory potential of nascent transcript remodeling by splicing, due to the limitations of existing methods that focus on analysis of mature splicing products (mRNAs) rather than substrates and intermediates. Here, we overcome this obstacle through long-read RNA sequencing of nascent, multi-intron transcripts in the fission yeast Schizosaccharomyces pombe. Most multi-intron transcripts were fully spliced, consistent with rapid cotranscriptional splicing. However, an unexpectedly high proportion of transcripts were either fully spliced or fully unspliced, suggesting that splicing of any given intron is dependent on the splicing status of other introns in the transcript. Supporting this, mild inhibition of splicing by a temperature-sensitive mutation in prp2, the homolog of vertebrate U2AF65, increased the frequency of fully unspliced transcripts. Importantly, fully unspliced transcripts displayed transcriptional read-through at the polyA site and were degraded cotranscriptionally by the nuclear exosome. Finally, we show that cellular mRNA levels were reduced in genes with a high number of unspliced nascent transcripts during caffeine treatment, showing regulatory significance of cotranscriptional splicing. Therefore, overall splicing of individual nascent transcripts, 3′ end formation, and mRNA half-life depend on the splicing status of neighboring introns, suggesting crosstalk among spliceosomes and the polyA cleavage machinery during transcription elongation.

Evaluation of two main RNA-seq approaches for gene quantification in clinical RNA sequencing: polyA+ selection versus rRNA depletion

To allow efficient transcript/gene detection, highly abundant ribosomal RNAs (rRNA) are generally removed from total RNA either by positive polyA+ selection or by rRNA depletion (negative selection) before sequencing. Comparisons between the two methods have been carried out by various groups, but the assessments have relied largely on non-clinical samples. In this study, we evaluated these two RNA sequencing approaches using human blood and colon tissue samples. Our analyses showed that rRNA depletion captured more unique transcriptome features, whereas polyA+ selection outperformed rRNA depletion with higher exonic coverage and better accuracy of gene quantification. For blood- and colon-derived RNAs, we found that 220% and 50% more reads, respectively, would have to be sequenced to achieve the same level of exonic coverage in the rRNA depletion method compared with the polyA+ selection method. Therefore, in most cases we strongly recommend polyA+ selection over rRNA depletion for gene quantification in clinical RNA sequencing. Our evaluation revealed that a small number of lncRNAs and small RNAs made up a large fraction of the reads in the rRNA depletion RNA sequencing data. Thus, we recommend that these RNAs are specifically depleted to improve the sequencing depth of the remaining RNAs.

The circadian dynamics of small nucleolar RNA in the mouse liver

The circadian regulation of gene expression allows plants and animals to anticipate predictable environmental changes. While the influence of the circadian clock has recently been shown to extend to ribosome biogenesis, the dynamics and regulation of the many small nucleolar RNA that are required in pre-ribosomal RNA folding and modification are unknown. Using a novel computational method, we show that 18S and 28S pre-rRNA are subject to circadian regulation in a nuclear RNA sequencing time course. A population of snoRNA with circadian expression is identified that is functionally associated with rRNA modification. More generally, we find the abundance of snoRNA known to modify 18S and 28S to be inversely correlated with the abundance of their target. Cyclic patterns in the expression of a number of snoRNA indicate a coordination with rRNA maturation, potentially through an upregulation in their biogenesis, or their release from mature rRNA at the end of the previous cycle of rRNA maturation, in antiphase with the diurnal peak in pre-rRNA. Few cyclic snoRNA have cyclic host genes, indicating the action of regulatory mechanisms in addition to transcriptional activation of the host gene. For highly expressed independently transcribed snoRNA, we find a characteristic RNA polymerase II and H3K4me3 signature that correlates with mean snoRNA expression over the day.

Dynamic RNA-protein interactions underlie the zebrafish maternal-to-zygotic transition

During the maternal-to-zygotic transition (MZT), transcriptionally silent embryos rely on post-transcriptional regulation of maternal mRNAs until zygotic genome activation (ZGA). RNA-binding proteins (RBPs) are important regulators of post-transcriptional RNA processing events, yet their identities and functions during developmental transitions in vertebrates remain largely unexplored. Using mRNA interactome capture, we identified 227 RBPs in zebrafish embryos before and during ZGA, hereby named the zebrafish MZT mRNA-bound proteome. This protein constellation consists of many conserved RBPs, some of which are potential stage-specific mRNA interactors that likely reflect the dynamics of RNA-protein interactions during MZT. The enrichment of numerous splicing factors like hnRNP proteins before ZGA was surprising, because maternal mRNAs were found to be fully spliced. To address potentially unique roles of these RBPs in embryogenesis, we focused on Hnrnpa1. iCLIP and subsequent mRNA reporter assays revealed a function for Hnrnpa1 in the regulation of poly(A) tail length and translation of maternal mRNAs through sequence-specific association with 3' UTRs before ZGA. Comparison of iCLIP data from two developmental stages revealed that Hnrnpa1 dissociates from maternal mRNAs at ZGA and instead regulates the nuclear processing of pri-mir-430 transcripts, which we validated experimentally. The shift from cytoplasmic to nuclear RNA targets was accompanied by a dramatic translocation of Hnrnpa1 and other pre-mRNA splicing factors to the nucleus in a transcription-dependent manner. Thus, our study identifies global changes in RNA-protein interactions during vertebrate MZT and shows that Hnrnpa1 RNA-binding activities are spatially and temporally coordinated to regulate RNA metabolism during early development.

Workflow for Genome-Wide Determination of Pre-mRNA Splicing Efficiency from Yeast RNA-seq Data

Pre-mRNA splicing represents an important regulatory layer of eukaryotic gene expression. In the simple budding yeast Saccharomyces cerevisiae, about one-third of all mRNA molecules undergo splicing, and splicing efficiency is tightly regulated, for example, during meiotic differentiation. S. cerevisiae features a streamlined, evolutionarily highly conserved splicing machinery and serves as a favourite model for studies of various aspects of splicing. RNA-seq represents a robust, versatile, and affordable technique for transcriptome interrogation, which can also be used to study splicing efficiency. However, convenient bioinformatics tools for the analysis of splicing efficiency from yeast RNA-seq data are lacking. We present a complete workflow for the calculation of genome-wide splicing efficiency in S. cerevisiae using strand-specific RNA-seq data. Our pipeline takes sequencing reads in the FASTQ format and provides splicing efficiency values for the 5′ and 3′ splice junctions of each intron. The pipeline is based on up-to-date open-source software tools and requires very limited input from the user. We provide all relevant scripts in a ready-to-use form. We demonstrate the functionality of the workflow using RNA-seq datasets from three spliceosome mutants. The workflow should prove useful for studies of yeast splicing mutants or of regulated splicing, for example, under specific growth conditions.

Perfect timing: splicing and transcription rates in living cells

An important step toward understanding gene regulation is the elucidation of the time necessary for the completion of individual steps. Measurement of reaction rates can reveal potential nodes for regulation. For example, measurements of in vivo transcription elongation rates reveal regulation by DNA sequence, gene architecture, and chromatin. Pre-mRNA splicing is regulated by transcription elongation rates and vice versa, yet the rates of RNA processing reactions remain largely elusive. Since the 1980s, numerous model systems and approaches have been used to determine the precise timing of splicing in vivo. Because splicing can be co-transcriptional, the position of Pol II when splicing is detected has been used as a proxy for time by some investigators. In addition to these 'distance-based' measurements, 'time-based' measurements have been possible through live cell imaging, metabolic labeling of RNA, and gene induction. Yet splicing rates can be convolved by the time it takes for transcription, spliceosome assembly and spliceosome disassembly. The variety of assays and systems used has, perhaps not surprisingly, led to reports of widely differing splicing rates in vivo. Recently, single molecule RNA-seq has indicated that splicing occurs more quickly than previously deduced. Here we comprehensively review these findings and discuss evidence that splicing and transcription rates are closely coordinated, facilitating the efficiency of gene expression. On the other hand, introduction of splicing delays through as yet unknown mechanisms provide opportunity for regulation. More work is needed to understand how cells optimize the rates of gene expression for a range of biological conditions. WIREs RNA 2017, 8:e1401. doi: 10.1002/wrna.1401 For further resources related to this article, please visit the WIREs website.

Identification of Coilin Mutants in a Screen for Enhanced Expression of an Alternatively Spliced GFP Reporter Gene in Arabidopsis thaliana

Coilin is a marker protein for subnuclear organelles known as Cajal bodies, which are sites of various RNA metabolic processes including the biogenesis of spliceosomal small nuclear ribonucleoprotein particles. Through self-associations and interactions with other proteins and RNA, coilin provides a structural scaffold for Cajal body formation. However, despite a conspicuous presence in Cajal bodies, most coilin is dispersed in the nucleoplasm and expressed in cell types that lack these organelles. The molecular function of coilin, particularly of the substantial nucleoplasmic fraction, remains uncertain. We identified coilin loss-of-function mutations in a genetic screen for mutants showing either reduced or enhanced expression of an alternatively spliced GFP reporter gene in Arabidopsis thaliana The coilin mutants feature enhanced GFP fluorescence and diminished Cajal bodies compared with wild-type plants. The amount of GFP protein is several-fold higher in the coilin mutants owing to elevated GFP transcript levels and more efficient splicing to produce a translatable GFP mRNA. Genome-wide RNA-sequencing data from two distinct coilin mutants revealed a small, shared subset of differentially expressed genes, many encoding stress-related proteins, and, unexpectedly, a trend toward increased splicing efficiency. These results suggest that coilin attenuates splicing and modulates transcription of a select group of genes. The transcriptional and splicing changes observed in coilin mutants are not accompanied by gross phenotypic abnormalities or dramatically altered stress responses, supporting a role for coilin in fine tuning gene expression. Our GFP reporter gene provides a sensitive monitor of coilin activity that will facilitate further investigations into the functions of this enigmatic protein.

Lessons from non-canonical splicing

Recent improvements in experimental and computational techniques that are used to study the transcriptome have enabled an unprecedented view of RNA processing, revealing many previously unknown non-canonical splicing events. This includes cryptic events located far from the currently annotated exons and unconventional splicing mechanisms that have important roles in regulating gene expression. These non-canonical splicing events are a major source of newly emerging transcripts during evolution, especially when they involve sequences derived from transposable elements. They are therefore under precise regulation and quality control, which minimizes their potential to disrupt gene expression. We explain how non-canonical splicing can lead to aberrant transcripts that cause many diseases, and also how it can be exploited for new therapeutic strategies.

Splicing of Nascent RNA Coincides with Intron Exit from RNA Polymerase II

Protein-coding genes in eukaryotes are transcribed by RNA polymerase II (Pol II) and introns are removed from pre-mRNA by the spliceosome. Understanding the time lag between Pol II progression and splicing could provide mechanistic insights into the regulation of gene expression. Here, we present two single-molecule nascent RNA sequencing methods that directly determine the progress of splicing catalysis as a function of Pol II position. Endogenous genes were analyzed on a global scale in budding yeast. We show that splicing is 50% complete when Pol II is only 45 nt downstream of introns, with the first spliced products observed as introns emerge from Pol II. Perturbations that slow the rate of spliceosome assembly or speed up the rate of transcription caused splicing delays, showing that regulation of both processes determines in vivo splicing profiles. We propose that matched rates streamline the gene expression pathway, while allowing regulation through kinetic competition.

De novo Sequencing of the Leaf Transcriptome Reveals Complex Light-Responsive Regulatory Networks in Camellia sinensis cv. Baijiguan

Tea plants (Camellia sinensis L.) possess high genetic diversity that is important for breeding. One cultivar, Baijiguan, exhibits a yellow leaf phenotype, reduced chlorophyll (Chl) content, and aberrant chloroplast structures under high light intensity. In contrast, under low light intensity, the flush shoot from Baijiguan becomes green, the Chl content increases significantly, and the chloroplasts exhibit normal structures. To understand the underlying molecular mechanisms for these observations, we performed de novo transcriptome sequencing and digital gene expression (DGE) profiling using Illumina sequencing technology. De novo transcriptome assembly identified 88,788 unigenes, including 1652 transcription factors from 25 families. In total, 1993 and 2576 differentially expressed genes (DEGs) were identified in Baijiguan plants treated with 3 and 6 days of shade, respectively. Gene Ontology (GO) and pathway enrichment analyses indicated that the DEGs are predominantly involved in the ROS scavenging system, chloroplast development, photosynthetic pigment synthesis, secondary metabolism, and circadian systems. The light-responsive gene POR (protochlorophyllide oxidoreductase) and transcription factor HY5 were identified. Quantitative real-time PCR (qRT-PCR) analysis of 20 selected DEGs confirmed the RNA-sequencing (RNA-Seq) results. Overall, these findings suggest that high light intensity inhibits the expression of photosystem II 10-kDa protein (PsbR) in Baijiguan, thus affecting PSII stability, chloroplast development and chlorophyll biosynthesis.

Transcriptome-wide RNA processing kinetics revealed using extremely short 4tU labeling

BACKGROUND

RNA levels detected at steady state are the consequence of multiple dynamic processes within the cell. In addition to synthesis and decay, transcripts undergo processing. Metabolic tagging with a nucleotide analog is one way of determining the relative contributions of synthesis, decay and conversion processes globally.

RESULTS

By improving 4-thiouracil labeling of RNA in Saccharomyces cerevisiae we were able to isolate RNA produced during as little as 1 minute, allowing the detection of nascent pervasive transcription. Nascent RNA labeled for 1.5, 2.5 or 5 minutes was isolated and analyzed by reverse transcriptase-quantitative polymerase chain reaction and RNA sequencing. High kinetic resolution enabled detection and analysis of short-lived non-coding RNAs as well as intron-containing pre-mRNAs in wild-type yeast. From these data we measured the relative stability of pre-mRNA species with different high turnover rates and investigated potential correlations with sequence features.

CONCLUSIONS

Our analysis of non-coding RNAs reveals a highly significant association between non-coding RNA stability, transcript length and predicted secondary structure. Our quantitative analysis of the kinetics of pre-mRNA splicing in yeast reveals that ribosomal protein transcripts are more efficiently spliced if they contain intron secondary structures that are predicted to be less stable. These data, in combination with previous results, indicate that there is an optimal range of stability of intron secondary structures that allows for rapid splicing.