The primary transcriptome of barley chloroplasts: numerous noncoding RNAs and the dominating role of the plastid-encoded RNA polymerase

Catalyzing plant science research with RNA-seq

Next generation DNA sequencing technologies are driving increasingly rapid, affordable and high resolution analyses of plant transcriptomes through sequencing of their associated cDNA (complementary DNA) populations; an analytical platform commonly referred to as RNA-sequencing (RNA-seq). Since entering the arena of whole genome profiling technologies only a few years ago, RNA-seq has proven itself to be a powerful tool with a remarkably diverse range of applications, from detailed studies of biological processes at the cell type-specific level, to providing insights into fundamental questions in plant biology on an evolutionary time scale. Applications include generating genomic data for heretofore unsequenced species, thus expanding the boundaries of what had been considered "model organisms," elucidating structural and regulatory gene networks, revealing how plants respond to developmental cues and their environment, allowing a better understanding of the relationships between genes and their products, and uniting the "omics" fields of transcriptomics, proteomics, and metabolomics into a now common systems biology paradigm. We provide an overview of the breadth of such studies and summarize the range of RNA-seq protocols that have been developed to address questions spanning cell type-specific-based transcriptomics, transcript secondary structure and gene mapping.

RNA processing and decay in plastids

Plastids were derived through endosymbiosis from a cyanobacterial ancestor, whose uptake was followed by massive gene transfer to the nucleus, resulting in the compact size and modest coding capacity of the extant plastid genome. Plastid gene expression is essential for plant development, but depends on nucleus-encoded proteins recruited from cyanobacterial or host-cell origins. The plastid genome is heavily transcribed from numerous promoters, giving posttranscriptional events a critical role in determining the quantity and sizes of accumulating RNA species. The major events reviewed here are RNA editing, which restores protein conservation or creates correct open reading frames by converting C residues to U, RNA splicing, which occurs both in cis and trans, and RNA cleavage, which relies on a variety of exoribonucleases and endoribonucleases. Because the RNases have little sequence specificity, they are collectively able to remove extraneous RNAs whose ends are not protected by RNA secondary structures or sequence-specific RNA-binding proteins (RBPs). Other plastid RBPs, largely members of the helical-repeat superfamily, confer specificity to editing and splicing reactions. The enzymes that catalyze RNA processing are also the main actors in RNA decay, implying that these antagonistic roles are optimally balanced. We place the actions of RBPs and RNases in the context of a recent proteomic analysis that identifies components of the plastid nucleoid, a protein-DNA complex with multiple roles in gene expression. These results suggest that sublocalization and/or concentration gradients of plastid proteins could underpin the regulation of RNA maturation and degradation.

Contradiction between plastid gene transcription and function due to complex posttranscriptional splicing: an exemplary study of ycf15 function and evolution in angiosperms

Plant chloroplast genes are usually co-transcribed while its posttranscriptional splicing is fairly complex and remains largely unsolved. On basis of sequencing the three complete Camellia (Theaceae) chloroplast genomes for the first time, we comprehensively analyzed the evolutionary patterns of ycf15, a plastid gene quite paradoxical in terms of its function and evolution, along the inferred angiosperm phylogeny. Although many species in separate lineages including the three species reported here contained an intact ycf15 gene in their chloroplast genomes, the phylogenetic mixture of both intact and obviously disabled ycf15 genes imply that they are all non-functional. Both intracellular gene transfer (IGT) and horizontal gene transfer (HGT) failed to explain such distributional anomalies. While, transcriptome analyses revealed that ycf15 was transcribed as precursor polycistronic transcript which contained ycf2, ycf15 and antisense trnL-CAA. The transcriptome assembly was surprisingly found to cover near the complete Camellia chloroplast genome. Many non-coding regions including pseudogenes were mapped by multiple transcripts, indicating the generality of pseudogene transcriptions. Our results suggest that plastid DNA posttranscriptional splicing may involve complex cleavage of non-functional genes.

Plastid non-coding RNAs: emerging candidates for gene regulation

Recent advances in transcriptomics and bioinformatics, specifically strand-specific RNA sequencing, have allowed high-throughput, comprehensive detection of low-abundance transcripts typical of the non-coding RNAs studied in bacteria and eukaryotes. Before this, few plastid non-coding RNAs (pncRNAs) had been identified, and even fewer had been investigated for any functional role in gene regulation. Relaxed plastid transcription initiation and termination result in full transcription of both chloroplast DNA strands. Following this, post-transcriptional processing produces a pool of metastable RNA species, including distinct pncRNAs. Here we review pncRNA biogenesis and possible functionality, and speculate that this RNA class may have an underappreciated role in plastid gene regulation.

Effects of reduced chloroplast gene copy number on chloroplast gene expression in maize

Chloroplasts and other members of the plastid organelle family contain a small genome of bacterial ancestry. Young chloroplasts contain hundreds of genome copies, but the functional significance of this high genome copy number has been unclear. We describe molecular phenotypes associated with mutations in a nuclear gene in maize (Zea mays), white2 (w2), encoding a predicted organellar DNA polymerase. Weak and strong mutant alleles cause a moderate (approximately 5-fold) and severe (approximately 100-fold) decrease in plastid DNA copy number, respectively, as assayed by quantitative PCR and Southern-blot hybridization of leaf DNA. Both alleles condition a decrease in most chloroplast RNAs, with the magnitude of the RNA deficiencies roughly paralleling that of the DNA deficiency. However, some RNAs are more sensitive to a decrease in genome copy number than others. The rpoB messenger RNA (mRNA) exhibited a unique response, accumulating to dramatically elevated levels in response to a moderate reduction in plastid DNA. Subunits of photosynthetic enzyme complexes were reduced more severely than were plastid mRNAs, possibly because of impaired translation resulting from limiting ribosomal RNA, transfer RNA, and ribosomal protein mRNA. These results indicate that chloroplast genome copy number is a limiting factor for the expression of a subset of chloroplast genes in maize. Whereas in Arabidopsis (Arabidopsis thaliana) a pair of orthologous genes function redundantly to catalyze DNA replication in both mitochondria and chloroplasts, the w2 gene is responsible for virtually all chloroplast DNA replication in maize. Mitochondrial DNA copy number was reduced approximately 2-fold in mutants harboring strong w2 alleles, suggesting that w2 also contributes to mitochondrial DNA replication.

Rogowsky P. PPR8522 encodes a chloroplast-targeted pentatricopeptide repeat protein necessary for maize embryogenesis and vegetative development

The pentatricopeptide repeat (PPR) domain is an RNA binding domain allowing members of the PPR superfamily to participate in post-transcriptional processing of organellar RNA. Loss of PPR8522 from maize (Zea mays) confers an embryo-specific (emb) phenotype. The emb8522 mutation was isolated in an active Mutator (Mu) population and co-segregation analysis revealed that it was tightly linked to a MuDR insertion in the first exon of PPR8522. Independent evidence that disruption of PPR8522 caused the emb phenotype was provided by fine mapping to a region of 116kb containing no other gene than PPR8522 and complementation of the emb8522 mutant by a PPR8522 cDNA. The deduced PPR8522 amino acid sequence of 832 amino acids contains 10 PPR repeats and a chloroplast target peptide, the function of which was experimentally demonstrated by transient expression in Nicotiana benthamiana. Whereas mutant endosperm is apparently normal, mutant embryos deviate from normal development as early as 3 days after pollination, are reduced in size, exhibit more or less severe morphological aberrations depending on the genetic background, and generally do not germinate. The emb8522 mutation is the first to associate the loss of a PPR gene with an embryo-lethal phenotype in maize. Analyses of mutant plantlets generated by embryo-rescue experiments indicate that emb8522 also affects vegetative plant growth and chloroplast development. The loss of chloroplast transcription dependent on plastid-encoded RNA polymerase is the likely cause for the lack of an organized thylakoid network and an albino, seedling-lethal phenotype.

Positive regulation of psbA gene expression by cis-encoded antisense RNAs in Synechocystis sp. PCC 6803

The D1 protein of photosystem II in the thylakoid membrane of photosynthetic organisms is encoded by psbA genes, which in cyanobacteria occur in the form of a small gene family. Light-dependent up-regulation of psbA gene expression is crucial to ensure the proper replacement of the D1 protein. To gain a high level of gene expression, psbA transcription can be enhanced by several orders of magnitude. Recent transcriptome analyses demonstrated a high number of cis-encoded antisense RNAs (asRNAs) in bacteria, but very little is known about their possible functions. Here, we show the presence of two cis-encoded asRNAs (PsbA2R and PsbA3R) of psbA2 and psbA3 from Synechocystis sp. PCC 6803. These asRNAs are located in the 5' untranslated region of psbA2 and psbA3 genes. Their expression becomes up-regulated by light and down-regulated by darkness, similar to their target mRNAs. In the PsbA2R-suppressing strain [PsbA2R(-)], the amount of psbA2 mRNA was only about 50% compared with the control strain. Likewise, we identified a 15% lowered activity of photosystem II and a reduced amount of the D1 protein in PsbA2R(-) compared with the control strain. The function of PsbA2R in the stabilization of psbA2 mRNA was shown from in vitro RNase E assay when the AU box and the ribosome-binding site in the 5' untranslated region of psbA2 mRNA were both covered by PsbA2R. These results add another layer of complexity to the mechanisms that contribute to psbA gene expression and show PsbA2R as a positively acting factor to achieve a maximum level of D1 synthesis.

Motif analysis unveils the possible co-regulation of chloroplast genes and nuclear genes encoding chloroplast proteins

Chloroplasts play critical roles in land plant cells. Despite their importance and the availability of at least 200 sequenced chloroplast genomes, the number of known DNA regulatory sequences in chloroplast genomes are limited. In this paper, we designed computational methods to systematically study putative DNA regulatory sequences in intergenic regions near chloroplast genes in seven plant species and in promoter sequences of nuclear genes in Arabidopsis and rice. We found that -35/-10 elements alone cannot explain the transcriptional regulation of chloroplast genes. We also concluded that there are unlikely motifs shared by intergenic sequences of most of chloroplast genes, indicating that these genes are regulated differently. Finally and surprisingly, we found five conserved motifs, each of which occurs in no more than six chloroplast intergenic sequences, are significantly shared by promoters of nuclear-genes encoding chloroplast proteins. By integrating information from gene function annotation, protein subcellular localization analyses, protein-protein interaction data, and gene expression data, we further showed support of the functionality of these conserved motifs. Our study implies the existence of unknown nuclear-encoded transcription factors that regulate both chloroplast genes and nuclear genes encoding chloroplast protein, which sheds light on the understanding of the transcriptional regulation of chloroplast genes.

The antisense RNA As1_flv4 in the Cyanobacterium Synechocystis sp. PCC 6803 prevents premature expression of the flv4-2 operon upon shift in inorganic carbon supply

The functional relevance of natural cis-antisense transcripts is mostly unknown. Here we have characterized the association of three antisense RNAs and one intergenically encoded noncoding RNA with an operon that plays a crucial role in photoprotection of photosystem II under low carbon conditions in the cyanobacterium Synechocystis sp. PCC 6803. Cyanobacteria show strong gene expression dynamics in response to a shift of cells from high carbon to low levels of inorganic carbon (C_i), but the regulatory mechanisms are poorly understood. Among the most up-regulated genes in Synechocystis are flv4, sll0218, and flv2, which are organized in the flv4-2 operon. The flavodiiron proteins encoded by this operon open up an alternative electron transfer route, likely starting from the Q_B site in photosystem II, under photooxidative stress conditions. Our expression analysis of cells shifted from high carbon to low carbon demonstrated an inversely correlated transcript accumulation of the flv4-2 operon mRNA and one antisense RNA to flv4, designated as As1_flv4. Overexpression of As1_flv4 led to a decrease in flv4-2 mRNA. The promoter activity of as1_flv4 was transiently stimulated by C_i limitation and negatively regulated by the AbrB-like transcription regulator Sll0822, whereas the flv4-2 operon was positively regulated by the transcription factor NdhR. The results indicate that the tightly regulated antisense RNA As1_flv4 establishes a transient threshold for flv4-2 expression in the early phase after a change in C_i conditions. Thus, it prevents unfavorable synthesis of the proteins from the flv4-2 operon.

Transcriptional organization of the large and the small ATP synthase operons, atpI/H/F/A and atpB/E, in Arabidopsis thaliana chloroplasts

The ATP synthase is a ubiquitous enzyme which is found in bacteria and eukaryotic organelles. It is essential in the photosynthetic and respiratory processes, by transforming the electrochemical proton gradient into ATP energy via proton transport across the membranes. In Escherichia coli, the atp genes coding for the subunits of the ATP synthase enzyme are grouped in the same transcriptional unit, while in higher plants the plastid atp genes are organized into a large (atpI/H/F/A) and a small (atpB/E) atp operon. By using the model plant Arabidopsis thaliana, we have investigated the strategy evolved in chloroplasts to overcome the physical separation of the atp gene clusters and to coordinate their transcription. We show that all the identified promoters in the two atp operons are PEP dependent and require sigma factors for specific recognition. Our results indicate that transcription of the two atp operons is initiated by at least one common factor, the essential SIG2 factor. Our data show that SIG3 and SIG6 also participate in transcription initiation of the large and the small atp operon, respectively. We propose that SIG2 might be the factor responsible for coordinating the basal transcription of the plastid atp genes and that SIG3 and SIG6 might serve to modulate plastid atp expression with respect to physiological and environmental conditions. However, we observe that in the sigma mutants (sig2, sig3 and sig6) the deficiency in the recognition of specific atp promoters is largely balanced by mRNA stabilization and/or by activation of otherwise silent promoters, indicating that the rate-limiting step for expression of the atp operons is mostly post-transcriptional.

RIP1, a member of an Arabidopsis protein family, interacts with the protein RARE1 and broadly affects RNA editing

Transcripts of plant organelle genes are modified by cytidine-to-uridine (C-to-U) RNA editing, often changing the encoded amino acid predicted from the DNA sequence. Members of the PLS subclass of the pentatricopeptide repeat (PPR) motif-containing family are site-specific recognition factors for either chloroplast or mitochondrial C targets of editing. However, other than PPR proteins and the cis-elements on the organelle transcripts, no other components of the editing machinery in either organelle have previously been identified. The Arabidopsis chloroplast PPR protein Required for AccD RNA Editing 1 (RARE1) specifies editing of a C in the accD transcript. RARE1 was detected in a complex of >200 kDa. We immunoprecipitated epitope-tagged RARE1, and tandem MS/MS analysis identified a protein of unknown function lacking PPR motifs; we named it RNA-editing factor interacting protein 1 (RIP1). Yeast two-hybrid analysis confirmed RIP1 interaction with RARE1, and RIP1-GFP fusions were found in both chloroplasts and mitochondria. Editing assays for all 34 known Arabidopsis chloroplast targets in a rip1 mutant revealed altered efficiency of 14 editing events. In mitochondria, 266 editing events were found to have reduced efficiency, with major loss of editing at 108 C targets. Virus-induced gene silencing of RIP1 confirmed the altered editing efficiency. Transient introduction of a WT RIP1 allele into rip1 improved the defective RNA editing. The presence of RIP1 in a protein complex along with chloroplast editing factor RARE1 indicates that RIP1 is an important component of the RNA editing apparatus that acts on many chloroplast and mitochondrial C targets.

Reverse protection assay: a tool to analyze transcriptional rates from individual promoters

Transcriptional activity of entire genes in chloroplasts is usually assayed by run-on analyses. To determine not only the overall intensity of transcription of a gene, but also the rate of transcription from a particular promoter, we created the Reverse RNase Protection Assay (RePro): in-organello run-on transcription coupled to RNase protection to define distinct transcript ends during transcription. We demonstrate successful application of RePro in plastid promoter analysis and transcript 3' end processing.