Bioinformatic Analysis of Chloroplast Gene Expression and RNA Posttranscriptional Maturations Using RNA Sequencing

The Use of Nanopore Sequencing to Analyze the Chloroplast Transcriptome Part II: Bioinformatic Analyzes and Virtual RNA Blots

Nanopore sequencing of full-length cDNAs offers unprecedented details of the plastid RNA metabolism. After the generation of the nanopore reads, several bioinformatic steps are required to analyze the data. In this chapter, we describe in a few simple command lines the processing and mapping of the reads as well as the generation of virtual Northern blots as a simple and familiar way to visualize Nanopore sequencing data.

The Use of Nanopore Sequencing to Analyze the Chloroplast Transcriptome Part I: Library Preparation

Global understanding of plastid gene expression has always been impaired by its complexity. RNA splicing, editing, and intercistronic processing create multiple transcripts isoforms that can hardly be resolved using traditional molecular biology techniques. During the last decade, the wide adoption of RNA-seq-based techniques has, however, allowed an unprecedented understanding of all the different steps of chloroplast gene expression, from transcription to translation. Current strategies are nonetheless unable to identify and quantify full length transcripts isoforms, a limitation that can now be overcome using Nanopore Sequencing. We here provide a complete protocol to produce, from total leaf RNA, cDNA libraries ready for Nanopore sequencing. While most Nanopore protocols take advantage of the mRNA polyA tail we here first ligate an RNA adapter to the 3' ends of the RNAs and use it to initiate the template switching reverse transcription. The cDNA is then prepared and indexed for use with the regular Oxford Nanopore v14 chemistry. This protocol is of particular interest to researchers willing to simultaneously study the multiple post-transcriptional processes prevalent in the chloroplast.

Response of the organellar and nuclear (post)transcriptomes of Arabidopsis to drought

Plants have evolved sophisticated mechanisms to cope with drought, which involve massive changes in nuclear gene expression. However, little is known about the roles of post-transcriptional processing of nuclear or organellar transcripts and how meaningful these changes are. To address these issues, we used RNA-sequencing after ribosomal RNA depletion to monitor (post)transcriptional changes during different times of drought exposure in Arabidopsis Col-0. Concerning the changes detected in the organellar transcriptomes, chloroplast transcript levels were globally reduced, editing efficiency dropped, but splicing was not affected. Mitochondrial transcripts were slightly elevated, while editing and splicing were unchanged. Conversely, alternative splicing (AS) affected nearly 1,500 genes (9% of expressed nuclear genes). Of these, 42% were regulated solely at the level of AS, representing transcripts that would have gone unnoticed in a microarray-based approach. Moreover, we identified 927 isoform switching events. We provide a table of the most interesting candidates, and as proof of principle, increased drought tolerance of the carbonic anhydrase ca1 and ca2 mutants is shown. In addition, altering the relative contributions of the spliced isoforms could increase drought resistance. For example, our data suggest that the accumulation of a nonfunctional FLM (FLOWERING LOCUS M) isoform and not the ratio of FLM-ß and -δ isoforms may be responsible for the phenotype of early flowering under long-day drought conditions. In sum, our data show that AS enhances proteome diversity to counteract drought stress and represent a valuable resource that will facilitate the development of new strategies to improve plant performance under drought.

Full Length Transcriptome Highlights the Coordination of Plastid Transcript Processing

Plastid gene expression involves many post-transcriptional maturation steps resulting in a complex transcriptome composed of multiple isoforms. Although short-read RNA-Seq has considerably improved our understanding of the molecular mechanisms controlling these processes, it is unable to sequence full-length transcripts. This information is crucial, however, when it comes to understanding the interplay between the various steps of plastid gene expression. Here, we describe a protocol to study the plastid transcriptome using nanopore sequencing. In the leaf of Arabidopsis thaliana, with about 1.5 million strand-specific reads mapped to the chloroplast genome, we could recapitulate most of the complexity of the plastid transcriptome (polygenic transcripts, multiple isoforms associated with post-transcriptional processing) using virtual Northern blots. Even if the transcripts longer than about 2500 nucleotides were missing, the study of the co-occurrence of editing and splicing events identified 42 pairs of events that were not occurring independently. This study also highlighted a preferential chronology of maturation events with splicing happening after most sites were edited.

Systematic sequencing of chloroplast transcript termini from Arabidopsis thaliana reveals >200 transcription initiation sites and the extensive imprints of RNA-binding proteins and secondary structures

Chloroplast transcription requires numerous quality control steps to generate the complex but selective mixture of accumulating RNAs. To gain insight into how this RNA diversity is achieved and regulated, we systematically mapped transcript ends by developing a protocol called Terminome-seq. Using Arabidopsis thaliana as a model, we catalogued >215 primary 5′ ends corresponding to transcription start sites (TSS), as well as 1628 processed 5′ ends and 1299 3′ ends. While most termini were found in intergenic regions, numerous abundant termini were also found within coding regions and introns, including several major TSS at unexpected locations. A consistent feature was the clustering of both 5′ and 3′ ends, contrasting with the prevailing description of discrete 5′ termini, suggesting an imprecision of the transcription and/or RNA processing machinery. Numerous termini correlated with the extremities of small RNA footprints or predicted stem-loop structures, in agreement with the model of passive RNA protection. Terminome-seq was also implemented for pnp1–1, a mutant lacking the processing enzyme polynucleotide phosphorylase. Nearly 2000 termini were altered in pnp1–1, revealing a dominant role in shaping the transcriptome. In summary, Terminome-seq permits precise delineation of the roles and regulation of the many factors involved in organellar transcriptome quality control.

Plant Ribonuclease J: An Essential Player in Maintaining Chloroplast RNA Quality Control for Gene Expression

RNA quality control is an indispensable but poorly understood process that enables organisms to distinguish functional RNAs from nonfunctional or inhibitory ones. In chloroplasts, whose gene expression activities are required for photosynthesis, retrograde signaling, and plant development, RNA quality control is of paramount importance, as transcription is relatively unregulated. The functional RNA population is distilled from this initial transcriptome by a combination of RNA-binding proteins and ribonucleases. One of the key enzymes is RNase J, a 5′→3′ exoribonuclease and an endoribonuclease that has been shown to trim 5′ RNA termini and eliminate deleterious antisense RNA. In the absence of RNase J, embryo development cannot be completed. Land plant RNase J contains a highly conserved C-terminal domain that is found in GT-1 DNA-binding transcription factors and is not present in its bacterial, archaeal, and algal counterparts. The GT-1 domain may confer specificity through DNA and/or RNA binding and/or protein–protein interactions and thus be an element in the mechanisms that identify target transcripts among diverse RNA populations. Further understanding of chloroplast RNA quality control relies on discovering how RNase J is regulated and how its specificity is imparted.

The Analysis of the Editing Defects in the dyw2 Mutant Provides New Clues for the Prediction of RNA Targets of Arabidopsis E+-Class PPR Proteins

C to U editing is one of the post-transcriptional steps which are required for the proper expression of chloroplast and mitochondrial genes in plants. It depends on several proteins acting together which include the PLS-class pentatricopeptide repeat proteins (PPR). DYW2 was recently shown to be required for the editing of many sites in both organelles. In particular almost all the sites associated with the E+ subfamily of PPR proteins are depending on DYW2, suggesting that DYW2 is required for the function of E+-type PPR proteins. Here we strengthened this link by identifying 16 major editing sites controlled by 3 PPR proteins: OTP90, a DYW-type PPR and PGN and MEF37, 2 E+-type PPR proteins. A re-analysis of the DYW2 editotype showed that the 49 sites known to be associated with the 18 characterized E+-type PPR proteins all depend on DYW2. Considering only the 288 DYW2-dependent editing sites as potential E+-type PPR sites, instead of the 795 known editing sites, improves the performances of binding predictions systems based on the PPR code for E+-type PPR proteins. However, it does not compensate for poor binding predictions.