ChloroSeq, an Optimized Chloroplast RNA-Seq Bioinformatic Pipeline, Reveals Remodeling of the Organellar Transcriptome Under Heat Stress

A prion-like domain is required for phase separation and chloroplast RNA processing during cold acclimation in Arabidopsis

Arabidopsis (Arabidopsis thaliana) plants can produce photosynthetic tissue with active chloroplasts at temperatures as low as 4°C, and this process depends on the presence of the nuclear-encoded, chloroplast-localized RNA-binding protein CP29A. In this study, we demonstrate that CP29A undergoes phase separation in vitro and in vivo in a temperature-dependent manner, which is mediated by a prion-like domain (PLD) located between the two RNA recognition motif domains of CP29A. The resulting droplets display liquid-like properties and are found near chloroplast nucleoids. The PLD is required to support chloroplast RNA splicing and translation in cold-treated tissue. Together, our findings suggest that plant chloroplast gene expression is compartmentalized by inducible condensation of CP29A at low temperatures, a mechanism that could play a crucial role in plant cold resistance.

Multi-scale analysis of heat stress acclimation in Arabidopsis seedlings highlights the primordial contribution of energy-transducing organelles

Much progress has been made in understanding the molecular mechanisms of plant adaptation to heat stress. However, the great diversity of models and stress conditions, and the fact that analyses are often limited to a small number of approaches, complicate the picture. We took advantage of a liquid culture system in which Arabidopsis seedlings are arrested in their development, thus avoiding interference with development and drought stress responses, to investigate through an integrative approach seedlings' global response to heat stress and acclimation. Seedlings perfectly tolerate a noxious heat shock (43°C) when subjected to a heat priming treatment at a lower temperature (38°C) the day before, displaying a thermotolerance comparable to that previously observed for Arabidopsis. A major effect of the pre-treatment was to partially protect energy metabolism under heat shock and favor its subsequent rapid recovery, which was correlated with the survival of seedlings. Rapid recovery of actin cytoskeleton and mitochondrial dynamics were another landmark of heat shock tolerance. The omics confirmed the role of the ubiquitous heat shock response actors but also revealed specific or overlapping responses to priming, heat shock, and their combination. Since only a few components or functions of chloroplast and mitochondria were highlighted in these analyses, the preservation and rapid recovery of their bioenergetic roles upon acute heat stress do not require extensive remodeling of the organelles. Protection of these organelles is rather integrated into the overall heat shock response, thus allowing them to provide the energy required to elaborate other cellular responses toward acclimation.

Long-read RNA-Seq for the discovery of long noncoding and antisense RNAs in plant organelles

Plant organelle transcription has been studied for decades. As techniques advanced, so did the fields of mitochondrial and plastid transcriptomics. The current view is that organelle genomes are pervasively transcribed, irrespective of their size, content, structure, and taxonomic origin. However, little is known about the nature of organelle noncoding transcriptomes, including pervasively transcribed noncoding RNAs (ncRNAs). Next-generation sequencing data have uncovered small ncRNAs in the organelles of plants and other organisms, but long ncRNAs remain poorly understood. Here, we argue that publicly available third-generation long-read RNA sequencing data from plants can provide a fine-tuned picture of long ncRNAs within organelles. Indeed, given their bloated architectures, plant mitochondrial genomes are well suited for studying pervasive transcription of ncRNAs. Ultimately, we hope to showcase this new avenue of plant research while also underlining the limitations of the proposed approach.

Long-read RNA sequencing can probe organelle genome pervasive transcription

40 years ago, organelle genomes were assumed to be streamlined and, perhaps, unexciting remnants of their prokaryotic past. However, the field of organelle genomics has exposed an unparallel diversity in genome architecture (i.e. genome size, structure, and content). The transcription of these eccentric genomes can be just as elaborate - organelle genomes are pervasively transcribed into a plethora of RNA types. However, while organelle protein-coding genes are known to produce polycistronic transcripts that undergo heavy posttranscriptional processing, the nature of organelle noncoding transcriptomes is still poorly resolved. Here, we review how wet-lab experiments and second-generation sequencing data (i.e. short reads) have been useful to determine certain types of organelle RNAs, particularly noncoding RNAs. We then explain how third-generation (long-read) RNA-Seq data represent the new frontier in organelle transcriptomics. We show that public repositories (e.g. NCBI SRA) already contain enough data for inter-phyla comparative studies and argue that organelle biologists can benefit from such data. We discuss the prospects of using publicly available sequencing data for organelle-focused studies and examine the challenges of such an approach. We highlight that the lack of a comprehensive database dedicated to organelle genomics/transcriptomics is a major impediment to the development of a field with implications in basic and applied science.

Uncovering the photosystem I assembly pathway in land plants

Photosystem I (PSI) is one of two large pigment-protein complexes responsible for converting solar energy into chemical energy in all oxygenic photosynthetic organisms. The PSI supercomplex consists of the PSI core complex and peripheral light-harvesting complex I (LHCI) in eukaryotic photosynthetic organisms. However, how the PSI complex assembles in land plants is unknown. Here we describe PHOTOSYSTEM I BIOGENESIS FACTOR 8 (PBF8), a thylakoid-anchored protein in Arabidopsis thaliana that is required for PSI assembly. PBF8 regulates two key consecutive steps in this process, the building of two assembly intermediates comprising eight or nine subunits, by interacting with PSI core subunits. We identified putative PBF8 orthologues in charophytic algae and land plants but not in Cyanobacteria or Chlorophyta. Our data reveal the major PSI assembly pathway in land plants. Our findings suggest that novel assembly mechanisms evolved during plant terrestrialization to regulate PSI assembly, perhaps as a means to cope with terrestrial environments.

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.

RNA Editing in Chloroplast: Advancements and Opportunities

Many eukaryotic and prokaryotic organisms employ RNA editing (insertion, deletion, or conversion) as a post-transcriptional modification mechanism. RNA editing events are common in these organelles of plants and have gained particular attention due to their role in the development and growth of plants, as well as their ability to cope with abiotic stress. Owing to rapid developments in sequencing technologies and data analysis methods, such editing sites are being accurately predicted, and many factors that influence RNA editing are being discovered. The mechanism and role of the pentatricopeptide repeat protein family of proteins in RNA editing are being uncovered with the growing realization of accessory proteins that might help these proteins. This review will discuss the role and type of RNA editing events in plants with an emphasis on chloroplast RNA editing, involved factors, gaps in knowledge, and future outlooks.

Does IR-loss promote plastome structural variation and sequence evolution?

Plastids are one of the main distinguishing characteristics of the plant cell. The plastid genome (plastome) of most autotrophic seed plants possesses a highly conserved quadripartite structure containing a large single-copy (LSC) and a small single-copy (SSC) region separated by two copies of the inverted repeat (termed as IR_A and IR_B). The IRs have been inferred to stabilize the plastid genome via homologous recombination-induced repair mechanisms. IR loss has been documented in seven autotrophic flowering plant lineages and two autotrophic gymnosperm lineages, and the plastomes of these species (with a few exceptions) are rearranged to a great extent. However, some plastomes containing normal IRs also show high structural variation. Therefore, the role of IRs in maintaining plastome stability is still controversial. In this study, we first integrated and compared genome structure and sequence evolution of representative plastomes of all nine reported IR-lacking lineages and those of their closest relative(s) with canonical inverted repeats (CRCIRs for short) to explore the role of the IR in maintaining plastome structural stability and sequence evolution. We found the plastomes of most IR-lacking lineages have experienced significant structural rearrangement, gene loss and duplication, accumulation of novel small repeats, and acceleration of synonymous substitution compared with those of their CRCIRs. However, the IR-lacking plastomes show similar structural variation and sequence evolution rate, and even less rearrangement distance, dispersed repeat number, tandem repeat number, indels frequency and GC3 content than those of IR-present plastomes with variation in Geraniaceae. We argue that IR loss is not a driver of these changes but is instead itself a consequence of other processes that more broadly shape both structural and sequence-level plastome evolution.

Organellar transcripts dominate the cellular mRNA pool across plants of varying ploidy levels

Mitochondrial and plastid functions depend on coordinated expression of proteins encoded by genomic compartments that have radical differences in copy number of organellar and nuclear genomes. In polyploids, doubling of the nuclear genome may add challenges to maintaining balanced expression of proteins involved in cytonuclear interactions. Here, we use ribo-depleted RNA sequencing (RNA-seq) to analyze transcript abundance for nuclear and organellar genomes in leaf tissue from four different polyploid angiosperms and their close diploid relatives. We find that even though plastid genomes contain <1% of the number of genes in the nuclear genome, they generate the majority (69.9 to 82.3%) of messenger RNA (mRNA) transcripts in the cell. Mitochondrial genes are responsible for a much smaller percentage (1.3 to 3.7%) of the leaf mRNA pool but still produce much higher transcript abundances per gene compared to nuclear genome. Nuclear genes encoding proteins that functionally interact with mitochondrial or plastid gene products exhibit mRNA expression levels that are consistently more than 10-fold lower than their organellar counterparts, indicating an extreme cytonuclear imbalance at the RNA level despite the predominance of equimolar interactions at the protein level. Nevertheless, interacting nuclear and organellar genes show strongly correlated transcript abundances across functional categories, suggesting that the observed mRNA stoichiometric imbalance does not preclude coordination of cytonuclear expression. Finally, we show that nuclear genome doubling does not alter the cytonuclear expression ratios observed in diploid relatives in consistent or systematic ways, indicating that successful polyploid plants are able to compensate for cytonuclear perturbations associated with nuclear genome doubling.

Context-Dependent Mutation Dynamics, Not Selection, Explains the Codon Usage Bias of Most Angiosperm Chloroplast Genes

Two competing proposals about the degree to which selection affects codon usage of angiosperm chloroplast genes are examined. The first, based on observations that codon usage does not match expectations under the naïve assumption that base composition will be identical at all neutral sites, is that selection plays a significant role. The second is that codon usage is determined almost solely by mutation bias and drift, with selection influencing only one or two highly expressed genes, in particular psbA. First it is shown that, as a result of an influence of neighboring base composition on mutation dynamics, compositional biases are expected to be widely divergent at different sites in the absence of selection. The observed mutation properties are then used to predict expected neutral codon usage biases and to show that observed deviations from the naïve expectations are in fact expected given the context-dependent mutational dynamics. It is also shown that there is a match between the observed and expected codon usage when context effects are taken into consideration, with psbA being a notable exception. Overall, the data support the model that selection is not a widespread factor affecting the codon usage of angiosperm chloroplast genes and highlight the need to have an accurate model of mutational dynamics.

CAF Proteins Help SOT1 Regulate the Stability of Chloroplast ndhA Transcripts

Protein-mediated RNA stabilization plays profound roles in chloroplast gene expression. Genetic studies have indicated that chloroplast ndhA transcripts, encoding a key subunit of the NADH dehydrogenase-like complex that mediates photosystem I cyclic electron transport and facilitates chlororespiration, are stabilized by PPR53 and its orthologs, but the underlying mechanisms are unclear. Here, we report that CHLOROPLAST RNA SPLICING 2 (CRS2)-ASSOCIATED FACTOR (CAF) proteins activate SUPPRESSOR OF THYLAKOID FORMATION 1 (SOT1), an ortholog of PPR53 in Arabidopsis thaliana, enhancing their affinity for the 5' ends of ndhA transcripts to stabilize these molecules while inhibiting the RNA endonuclease activity of the SOT1 C-terminal SMR domain. In addition, we established that SOT1 improves the splicing efficiency of ndhA by facilitating the association of CAF2 with the ndhA intron, which may be due to the SOT1-mediated stability of the ndhA transcripts. Our findings shed light on the importance of PPR protein interaction partners in moderating RNA metabolism.

Pre-mRNA alternative splicing as a modulator for heat stress response in plants

The molecular responses of plants to the important abiotic stress, heat stress (HS), have been extensively studied at the transcriptional level. Alternative splicing (AS) is a post-transcriptional regulatory process in which an intron-containing gene can generate more than one mRNA variant. The impact of HS on the pre-mRNA splicing process has been reported in various eukaryotes but seldom discussed in-depth, especially in plants. Here, we review AS regulation in response to HS in different plant species. We discuss potential molecular mechanisms controlling heat-inducible AS regulation in plants and hypothesize that AS regulation participates in heat-priming establishment and HS memory maintenance. We propose that the pre-mRNA splicing variation is an important regulator of plant HS responses (HSRs).

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.

OMICs, Epigenetics, and Genome Editing Techniques for Food and Nutritional Security

The incredible success of crop breeding and agricultural innovation in the last century greatly contributed to the Green Revolution, which significantly increased yields and ensures food security, despite the population explosion. However, new challenges such as rapid climate change, deteriorating soil, and the accumulation of pollutants require much faster responses and more effective solutions that cannot be achieved through traditional breeding. Further prospects for increasing the efficiency of agriculture are undoubtedly associated with the inclusion in the breeding strategy of new knowledge obtained using high-throughput technologies and new tools in the future to ensure the design of new plant genomes and predict the desired phenotype. This article provides an overview of the current state of research in these areas, as well as the study of soil and plant microbiomes, and the prospective use of their potential in a new field of microbiome engineering. In terms of genomic and phenomic predictions, we also propose an integrated approach that combines high-density genotyping and high-throughput phenotyping techniques, which can improve the prediction accuracy of quantitative traits in crop species.

Acclimation in plants - the Green Hub consortium

Acclimation is the capacity to adapt to environmental changes within the lifetime of an individual. This ability allows plants to cope with the continuous variation in ambient conditions to which they are exposed as sessile organisms. Because environmental changes and extremes are becoming even more pronounced due to the current period of climate change, enhancing the efficacy of plant acclimation is a promising strategy for mitigating the consequences of global warming on crop yields. At the cellular level, the chloroplast plays a central role in many acclimation responses, acting both as a sensor of environmental change and as a target of cellular acclimation responses. In this Perspective article, we outline the activities of the Green Hub consortium funded by the German Science Foundation. The main aim of this research collaboration is to understand and strategically modify the cellular networks that mediate plant acclimation to adverse environments, employing Arabidopsis, tobacco (Nicotiana tabacum) and Chlamydomonas as model organisms. These efforts will contribute to 'smart breeding' methods designed to create crop plants with improved acclimation properties. To this end, the model oilseed crop Camelina sativa is being used to test modulators of acclimation for their potential to enhance crop yield under adverse environmental conditions. Here we highlight the current state of research on the role of gene expression, metabolism and signalling in acclimation, with a focus on chloroplast-related processes. In addition, further approaches to uncovering acclimation mechanisms derived from systems and computational biology, as well as adaptive laboratory evolution with photosynthetic microbes, are highlighted.

Noncoding RNA: An Insight into Chloroplast and Mitochondrial Gene Expressions

Regulation of gene expression in any biological system is a complex process with many checkpoints at the transcriptional, post-transcriptional and translational levels. The control mechanism is mediated by various protein factors, secondary metabolites and a newly included regulatory member, i.e., noncoding RNAs (ncRNAs). It is known that ncRNAs modulate the mRNA or protein profiles of the cell depending on the degree of complementary and context of the microenvironment. In plants, ncRNAs are essential for growth and development in normal conditions by controlling various gene expressions and have emerged as a key player to guard plants during adverse conditions. In order to have smooth functioning of the plants under any environmental pressure, two very important DNA-harboring semi-autonomous organelles, namely, chloroplasts and mitochondria, are considered as main players. These organelles conduct the most crucial metabolic pathways that are required to maintain cell homeostasis. Thus, it is imperative to explore and envisage the molecular machineries responsible for gene regulation within the organelles and their coordination with nuclear transcripts. Therefore, the present review mainly focuses on ncRNAs origination and their gene regulation in chloroplasts and plant mitochondria.

Computational Detection of Plant RNA Editing Events

Computers are able to systematically exploit RNA-seq data allowing us to efficiently detect RNA editing sites in a genome-wide scale. This chapter introduces a very flexible computational framework for detecting RNA editing sites in plant organelles. This framework comprises three major steps: RNA-seq data processing, RNA read alignment, and RNA editing site detection. Each step is discussed in sufficient detail to be implemented by the reader. As a study case, the framework will be used with publicly available sequencing data to detect C-to-U RNA editing sites in the coding sequences of the mitochondrial genome of Nicotiana tabacum.

Pentatricopeptide repeat protein PHOTOSYSTEM I BIOGENESIS FACTOR2 is required for splicing of ycf3

To gain a better understanding of the molecular mechanisms of photosystem I (PSI) biogenesis, we characterized the Arabidopsis thaliana photosystem I biogenesis factor 2 (pbf2) mutant, which lacks PSI complex. PBF2 encodes a P-class pentatricopeptide repeat (PPR) protein. In the pbf2 mutants, we observed a striking decrease in the transcript level of only one gene, the chloroplast gene ycf3, which is essential for PSI assembly. Further analysis of ycf3 transcripts showed that PBF2 is specifically required for the splicing of ycf3 intron 1. Computational prediction of binding sequences and electrophoretic mobility shift assays reveal that PBF2 specifically binds to a sequence in ycf3 intron 1. Moreover, we found that PBF2 interacted with two general factors for group II intron splicing CHLOROPLAST RNA SPLICING2-ASSOCIATED FACTOR1 (CAF1) and CAF2, and facilitated the association of these two factors with ycf3 intron 1. Our results suggest that PBF2 is specifically required for the splicing of ycf3 intron 1 through cooperating with CAF1 and CAF2. Our results also suggest that additional proteins are required to contribute to the specificity of CAF-dependent group II intron splicing.

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.

Mutation of DELAYED GREENING1 impairs chloroplast RNA editing at elevated ambient temperature in Arabidopsis

Chloroplasts are important for plant growth and development. RNA editing in chloroplast converts cytidines (Cs) to uridines (Us) at specific transcript positions and provides a correction mechanism to restore conserved codons or creates start or stop codons. However, the underlined molecular mechanism is not yet fully understood. In the present study, we identified a thermo-sensitive mutant in leaf color 1 (tsl1) and found that TSL1 is allelic to DELAYED GREENING 1 (DG1). The missense mutation of DG1 in tsl1 mutant confers a high temperature sensitivity and impaired chloroplast development at an elevated ambient temperature in Arabidopsis. Subsequent analysis showed that chloroplast RNA editing at several sites including accD-1568, ndhD-2, and petL-5 is impaired in tsl1 mutant plants grown at an elevated temperature. DG1 interacts with MORF2 and other proteins such as DYW1 and DYW2 involved in chloroplast RNA editing. In vitro RNA electrophoretic mobility shift assay demonstrated that DG1 binds to RNA targets such as accD, ndhD, and petL. Thus, our results revealed that DG1 is important for maintaining chloroplast mRNA editing in Arabidopsis.

ECuADOR-Easy Curation of Angiosperm Duplicated Organellar Regions, a tool for cleaning and curating plastomes assembled from next generation sequencing pipelines

Background

With the rapid increase in availability of genomic resources offered by Next-Generation Sequencing (NGS) and the availability of free online genomic databases, efficient and standardized metadata curation approaches have become increasingly critical for the post-processing stages of biological data. Especially in organelle-based studies using circular chloroplast genome datasets, the assembly of the main structural regions in random order and orientation represents a major limitation in our ability to easily generate “ready-to-align” datasets for phylogenetic reconstruction, at both small and large taxonomic scales. In addition, current practices discard the most variable regions of the genomes to facilitate the alignment of the remaining coding regions. Nevertheless, no software is currently available to perform curation to such a degree, through simple detection, organization and positioning of the main plastome regions, making it a time-consuming and error-prone process. Here we introduce a fast and user friendly software ECuADOR, a Perl script specifically designed to automate the detection and reorganization of newly assembled plastomes obtained from any source available (NGS, sanger sequencing or assembler output).

Methods

ECuADOR uses a sliding-window approach to detect long repeated sequences in draft sequences, which then identifies the inverted repeat regions (IRs), even in case of artifactual breaks or sequencing errors and automates the rearrangement of the sequence to the widely used LSC–Irb–SSC–IRa order. This facilitates rapid post-editing steps such as creation of genome alignments, detection of variable regions, SNP detection and phylogenomic analyses.

Results

ECuADOR was successfully tested on plant families throughout the angiosperm phylogeny by curating 161 chloroplast datasets. ECuADOR first identified and reordered the central regions (LSC–Irb–SSC–IRa) for each dataset and then produced a new annotation for the chloroplast sequences. The process took less than 20 min with a maximum memory requirement of 150 MB and an accuracy of over 99%.

Conclusions

ECuADOR is the sole de novo one-step recognition and re-ordination tool that provides facilitation in the post-processing analysis of the extra nuclear genomes from NGS data. The program is available at https://github.com/BiodivGenomic/ECuADOR/.

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.

Methods and Tools for Plant Organelle Genome Sequencing, Assembly, and Downstream Analysis

Organelles play an important role in a eukaryotic cell. Among them, the two organelles, chloroplast and mitochondria, are responsible for the critical function of photosynthesis and aerobic respiration. Organellar genomes are also very important for plant systematic studies. Here we have described the methods for isolation of the mitochondrial and plastid DNA and its subsequent sequencing with the help of NGS technology. We have also discussed in detail the various tools available for assembly, annotation, and visualization of the organelle genome sequence.

Variation in plastid genomes in the gynodioecious species Silene vulgaris

BACKGROUND

Gynodioecious species exist in two sexes - male-sterile females and hermaphrodites. Male sterility in higher plants often results from mitonuclear interaction between the CMS (cytoplasmic male sterility) gene(s) encoded by mitochondrial genome and by nuclear-encoded restorer genes. Mitochondrial and nuclear-encoded transcriptomes in females and hermaphrodites are intensively studied, but little is known about sex-specific gene expression in plastids. We have compared plastid transcriptomes between females and hermaphrodites in two haplotypes of a gynodioecious species Silene vulgaris with known CMS candidate genes.

RESULTS

We generated complete plastid genome sequences from five haplotypes S. vulgaris including the haplotypes KRA and KOV, for which complete mitochondrial genome sequences were already published. We constructed a phylogenetic tree based on plastid sequences of S. vulgaris. Whereas lowland S. vulgaris haplotypes including KRA and KOV clustered together, the accessions from high European mountains diverged early in the phylogram. S. vulgaris belongs among Silene species with slowly evolving plastid genomes, but we still detected 212 substitutions and 112 indels between two accessions of this species. We estimated elevated Ka/Ks in the ndhF gene, which may reflect the adaptation of S. vulgaris to high altitudes, or relaxed selection. We compared depth of coverage and editing rates between female and hermaphrodite plastid transcriptomes and found no significant differences between the two sexes. We identified 51 unique C to U editing sites in the plastid genomes of S. vulgaris, 38 of them in protein coding regions, 2 in introns, and 11 in intergenic regions. The editing site in the psbZ gene was edited only in one of two plastid genomes under study.

CONCLUSIONS

We revealed no significant differences between the sexes in plastid transcriptomes of two haplotypes of S. vulgaris. It suggests that gene expression of plastid genes is not affected by CMS in flower buds of S. vulgaris, although both sexes may still differ in plastid gene expression in specific tissues. We revealed the difference between the plastid transcriptomes of two S. vulgaris haplotypes in editing rate and in the coverage of several antisense transcripts. Our results document the variation in plastid genomes and transcriptomes in S. vulgaris.

Complete mitochondrial genomes from transcriptomes: assessing pros and cons of data mining for assembling new mitogenomes

Thousands of eukaryotes transcriptomes have been generated, mainly to investigate nuclear genes expression, and the amount of available data is constantly increasing. A neglected but promising use of this large amount of data is to assemble organelle genomes. To assess the reliability of this approach, we attempted to reconstruct complete mitochondrial genomes from RNA-Seq experiments of Reticulitermes termite species, for which transcriptomes and conspecific mitogenomes are available. We successfully assembled complete molecules, although a few gaps corresponding to tRNAs had to be filled manually. We also reconstructed, for the first time, the mitogenome of Reticulitermes banyulensis. The accuracy and completeness of mitogenomes reconstruction appeared independent from transcriptome size, read length and sequencing design (single/paired end), and using reference genomes from congeneric or intra-familial taxa did not significantly affect the assembly. Transcriptome-derived mitogenomes were found highly similar to the conspecific ones obtained from genome sequencing (nucleotide divergence ranging from 0% to 3.5%) and yielded a congruent phylogenetic tree. Reads from contaminants and nuclear transcripts, although slowing down the process, did not result in chimeric sequence reconstruction. We suggest that the described approach has the potential to increase the number of available mitogenomes by exploiting the rapidly increasing number of transcriptomes.

MORF9 Functions in Plastid RNA Editing with Tissue Specificity

RNA editing in plant mitochondria and plastids converts specific nucleotides from cytidine (C) to uridine (U). These editing events differ among plant species and are relevant to developmental stages or are impacted by environmental conditions. Proteins of the MORF family are essential components of plant editosomes. One of the members, MORF9, is considered the core protein of the editing complex and is involved in the editing of most sites in chloroplasts. In this study, the phenotypes of a T-DNA insertion line with loss of MORF9 and of the genetic complementation line of Arabidopsis were analyzed, and the editing efficiencies of plastid RNAs in roots, rosette leaves, and flowers from the morf9 mutant and the wild-type (WT) control were compared by bulk-cDNA sequencing. The results showed that most of the known MORF9-associated plastid RNA editing events in rosette leaves and flowers were similarly reduced by morf9 mutation, with the exception that the editing rate of the sites ndhB-872 and psbF-65 declined in the leaves and that of ndhB-586 decreased only in the flowers. In the roots, however, the loss of MORF9 had a much lower effect on overall plastid RNA editing, with nine sites showing no significant editing efficiency change, including accD-794, ndhD-383, psbZ-50, ndhF-290, ndhD-878, matK-706, clpP1-559, rpoA-200, and ndhD-674, which were reduced in the other tissues. Furthermore, we found that during plant aging, MORF9 mRNA level, but not the protein level, was downregulated in senescent leaves. On the basis of these observations, we suggest that MORF9-mediated RNA editing is tissue-dependent and the resultant organelle proteomes are pertinent to the specific tissue functions.

Integrated Genomic and Transcriptomic Analysis of the Peridinin Dinoflagellate Amphidinium carterae Plastid

The plastid genomes of peridinin-containing dinoflagellates are highly unusual, possessing very few genes, which are located on small chromosomal elements termed "minicircles". These minicircles may contain genes, or no recognisable coding information. Transcripts produced from minicircles may undergo unusual processing events, such as the addition of a 3' poly(U) tail. To date, little is known about the genetic or transcriptional diversity of non-coding sequences in peridinin dinoflagellate plastids. These sequences include empty minicircles, and regions of non-coding DNA in coding minicircles. Here, we present an integrated plastid genome and transcriptome for the model peridinin dinoflagellate Amphidinium carterae, identifying a previously undescribed minicircle. We also profile transcripts covering non-coding regions of the psbA and petB/atpA minicircles. We present evidence that antisense transcripts are produced within the A. carterae plastid, but show that these transcripts undergo different end cleavage events from sense transcripts, and do not receive 3' poly(U) tails. The difference in processing events between sense and antisense transcripts may enable the removal of non-coding transcripts from peridinin dinoflagellate plastid transcript pools.

Transcription is a major driving force for plastid genome instability in Arabidopsis

Though it is an essential process, transcription can be a source of genomic instability. For instance, it may generate RNA:DNA hybrids as the nascent transcript hybridizes with the complementary DNA template. These hybrids, called R-loops, act as a major cause of replication fork stalling and DNA breaks. In this study, we show that lowering transcription and R-loop levels in plastids of Arabidopsis thaliana reduces DNA rearrangements and mitigates plastid genome instability phenotypes. This effect can be observed on a genome-wide scale, as the loss of the plastid sigma transcription factor SIG6 prevents DNA rearrangements by favoring conservative repair in the presence of ciprofloxacin-induced DNA damage or in the absence of plastid genome maintenance actors such as WHY1/WHY3, RECA1 and POLIB. Additionally, resolving R-loops by the expression of a plastid-targeted exogenous RNAse H1 produces similar results. We also show that highly-transcribed genes are more susceptible to DNA rearrangements, as increased transcription of the psbD operon by SIG5 correlates with more locus-specific rearrangements. The effect of transcription is not specific to Sigma factors, as decreased global transcription levels by mutation of heat-stress-induced factor HSP21, mutation of nuclear-encoded polymerase RPOTp, or treatment with transcription-inhibitor rifampicin all prevent the formation of plastid genome rearrangements, especially under induced DNA damage conditions.

RNA editing in plants: A comprehensive survey of bioinformatics tools and databases

RNA editing is a widespread epitranscriptomic mechanism by which primary RNAs are specifically modified through insertions/deletions or nucleotide substitutions. In plants, RNA editing occurs in organelles (plastids and mitochondria), involves the cytosine to uridine modification (rarely uridine to cytosine) within protein-coding and non-protein-coding regions of RNAs and affects organelle biogenesis, adaptation to environmental changes and signal transduction. High-throughput sequencing technologies have dramatically improved the detection of RNA editing sites at genomic scale. Consequently, different bioinformatics resources have been released to discovery and/or collect novel events. Here, we review and describe the state-of-the-art bioinformatics tools devoted to the characterization of RNA editing in plant organelles with the aim to improve our knowledge about this fascinating but yet under investigated process.

Multifarious Evolutionary Pathways of a Nuclear RNA Editing Factor: Disjunctions in Coevolution of DOT4 and Its Chloroplast Target rpoC1eU488SL

Nuclear-encoded pentatricopeptide repeat (PPR) proteins are site-specific factors for C-to-U RNA editing in plant organelles coevolving with their targets. Losing an editing target by C-to-T conversion allows for eventual loss of its editing factor, as recently confirmed for editing factors CLB19, CRR28, and RARE1 targeting ancient chloroplast editing sites in flowering plants. Here, we report on alternative evolutionary pathways for DOT4 addressing rpoC1eU488SL, a chloroplast editing site in the RNA polymerase β' subunit mRNA. Upon loss of rpoC1eU488SL by C-to-T conversion, DOT4 got lost multiple times independently in angiosperm evolution with intermediate states of DOT4 orthologs in various stages of degeneration. Surprisingly, we now also observe degeneration and loss of DOT4 despite retention of a C in the editing position (in Carica, Coffea, Vicia, and Spirodela). We find that the cytidine remains unedited, proving that DOT4 was not replaced by another editing factor. Yet another pathway of DOT4 evolution is observed among the Poaceae. Although the rpoC1eU488SL edit has been lost through C-to-T conversion, DOT4 orthologs not only remain conserved but also have their array of PPRs extended by six additional repeats. Here, the loss of the ancient target has likely allowed DOT4 to adapt for a new function. We suggest rps3 antisense transcripts as previously demonstrated in barley (Hordeum vulgare) arising from promotor sequences newly emerging in the rpl16 intron of Poaceae as a new candidate target for the extended PPR stretch of DOT4. Altogether, DOT4 and its target show more flexible pathways for evolution than the previously explored editing factors CLB19, CRR28, and RARE1. Certain plant clades (e.g., Amaranthus, Vaccinium, Carica, the Poaceae, Fabales, and Caryophyllales) show pronounced dynamics in the evolution of editing sites and corresponding factors.

Divergence of RNA editing among Arabidopsis species

RNA editing altered the RNA sequence by replacing the C nucleotide to U in the organellar genomes of plants. RNA editing status sometimes differed among distant species. The pattern of conservation and variation of RNA editing status made it possible to evaluate evolutionary mechanisms impacting functional aspects of RNA editing. In this study, divergence of RNA editing in the chloroplast genome among Arabidopsis species was analyzed to determine 9 losses and 1 gain in RNA editing. All changes in A. thaliana lineage resulted from changes to the chloroplast genome sequence, whereas changes in the A. lyrata / halleri lineage were possibly due to exclusive changes in the nuclear editing factors. One loss of RNA editing in A. lyrata was caused by a deficiency in the PPR gene OTP80. The changes in RNA editing occurred approximately every two million years and were not observed at functionally important sites. These results highlight the conserved nature of RNA editing status suggesting the importance of RNA editing during evolution.

Sequencing the Plastid Genome of Giant Ragweed (Ambrosia trifida, Asteraceae) From a Herbarium Specimen

We report the first plastome sequence of giant ragweed (Ambrosia trifida); with this new genome information, we assessed the phylogeny of Asteraceae and the transcriptional profiling against glyphosate resistance in giant ragweed. Assembly and genic features show a normal angiosperm quadripartite plastome structure with no signatures of deviation in gene directionality. Comparative analysis revealed large inversions across the plastome of giant ragweed and the previously sequenced members of the plant family. Asteraceae plastid genomes contain two inversions of 22.8 and 3.3 kb; the former is located between trnS-GCU and trnG-UCC genes, and the latter between trnE-UUC and trnT-GGU genes. The plastid genome sequences of A. trifida and the related species, Ambrosia artemisiifolia, are identical in gene content and arrangement, but they differ in length. The phylogeny is well-resolved and congruent with previous hypotheses about the phylogenetic relationship of Asteraceae. Transcriptomic analysis revealed divergence in the relative expressions at the exonic and intronic levels, providing hints toward the ecological adaptation of the genus. Giant ragweed shows various levels of glyphosate resistance, with introns displaying higher expression patterns at resistant time points after the assumed herbicide treatment.

Complete chloroplast genome sequence of Betula platyphylla: gene organization, RNA editing, and comparative and phylogenetic analyses

BACKGROUND

Betula platyphylla is a common tree species in northern China that has high economic and medicinal value. Our laboratory has been devoted to genome research on B. platyphylla for approximately 10 years. As primary organelle genomes, the complete genome sequences of chloroplasts are important to study the divergence of species, RNA editing and phylogeny. In this study, we sequenced and analyzed the complete chloroplast (cp) genome sequence of B. platyphylla.

RESULTS

The complete cp genome of B. platyphylla was 160,518 bp in length, which included a pair of inverted repeats (IRs) of 26,056 bp that separated a large single copy (LSC) region of 89,397 bp and a small single copy (SSC) region of 19,009 bp. The annotation contained a total of 129 genes, including 84 protein-coding genes, 37 tRNA genes and 8 rRNA genes. There were 3 genes using alternative initiation codons. Comparative genomics showed that the sequence of the Fagales species cp genome was relatively conserved, but there were still some high variation regions that could be used as molecular markers. The IR expansion event of B. platyphylla resulted in larger cp genomes and rps19 pseudogene formation. The simple sequence repeat (SSR) analysis showed that there were 105 SSRs in the cp genome of B. platyphylla. RNA editing sites recognition indicated that at least 80 RNA editing events occurred in the cp genome. Most of the substitutions were C to U, while a small proportion of them were not. In particular, three editing loci on the rRNA were converted to more than two other bases that had never been reported. For synonymous conversion, most of them increased the relative synonymous codon usage (RSCU) value of the codons. The phylogenetic analysis suggested that B. platyphylla had a closer evolutionary relationship with B. pendula than B. nana.

CONCLUSIONS

In this study, we not only obtained and annotated the complete cp genome sequence of B. platyphylla, but we also identified new RNA editing sites and predicted the phylogenetic relationships among Fagales species. These findings will facilitate genomic, genetic engineering and phylogenetic studies of this important species.

CaHSL1 Acts as a Positive Regulator of Pepper Thermotolerance Under High Humidity and Is Transcriptionally Modulated by CaWRKY40

Pepper (Capsicum annuum) is an economically important vegetable and heat stress can severely impair pepper growth, development, and productivity. The molecular mechanisms underlying pepper thermotolerance are therefore important to understand but remain elusive. In the present study, we characterized the function of CaHSL1, encoding a HAESA-LIKE (HSL) receptor-like protein kinase (RLK), during the response of pepper to high temperature and high humidity (HTHH). CaHSL1 exhibits the typical structural features of an arginine-aspartate RLK. Transient overexpression of CaHSL1 in the mesophyll cells of Nicotiana benthamiana showed that CaHSL1 localizes throughout the cell, including the plasma membrane, cytoplasm, and the nucleus. CaHSL1 was significantly upregulated by HTHH or the exogenous application of abscisic acid but not by R. solanacearum inoculation. However, CaHSL1 was downregulated by exogenously applied salicylic acid, methyl jasmonate, or ethephon. Silencing of CaHSL1 by virus-induced gene silencing significantly was reduced tolerance to HTHH and downregulated transcript levels of an associated gene CaHSP24. In contrast, transient overexpression of CaHSL1 enhanced the transcript abundance of CaHSP24 and increased tolerance to HTHH, as manifested by enhanced optimal/maximal photochemical efficiency of photosystem II in the dark (Fv/Fm) and actual photochemical efficiency of photosystem II in the light. In addition, CaWRKY40 targeted the promoter of CaHSL1 and induced transcription during HTHH but not in response to R. solanacearum. All of these results suggest that CaHSL1 is directly modulated at the transcriptional level by CaWRKY40 and functions as a positive regulator in the response of pepper to HTHH.

Unraveling chloroplast transcriptomes with ChloroSeq, an organelle RNA-Seq bioinformatics pipeline

Online sequence repositories are teeming with RNA sequencing (RNA-Seq) data from a wide range of eukaryotes. Although most of these data sets contain large numbers of organelle-derived reads, researchers tend to ignore these data, focusing instead on the nuclear-derived transcripts. Consequently, GenBank contains massive amounts of organelle RNA-Seq data that are just waiting to be downloaded and analyzed. Recently, a team of scientists designed an open-source bioinformatics program called ChloroSeq, which systemically analyzes an organelle transcriptome using RNA-Seq. The ChloroSeq pipeline uses RNA-Seq alignment data to deliver detailed analyses of organelle transcriptomes, which can be fed into statistical software for further analysis and for generating graphical representations of the data. In addition to providing data on expression levels via coverage statistics, ChloroSeq can examine splicing efficiency and RNA editing profiles. Ultimately, ChloroSeq provides a well-needed avenue for researchers of all stripes to start exploring organelle transcription and could be a key step toward a more thorough understanding of organelle gene expression.

Comparative assessment of chloroplast transcriptional responses highlights conserved and unique patterns across Triticeae members under salt stress

Chloroplast functional genomics, in particular understanding the chloroplast transcriptional response is of immense importance mainly due to its role in oxygenic photosynthesis. As a photosynthetic unit, its efficiency and transcriptional activity is directly regulated by reactive oxygen species during abiotic and biotic stress and subsequently affects carbon assimilation, and plant biomass. In crops, understanding photosynthesis is crucial for crop domestication by identifying the traits that could be exploited for crop improvement. Transcriptionally and translationally active chloroplast plays a key role by regulating the PSI and PSII photo-reaction centres, which ubiquitously affects the light harvesting. Using a comparative transcriptomics mapping approach, we identified differential regulation of key chloroplast genes during salt stress across Triticeae members with potential genes involved in photosynthesis and electron transport system such as CytB6f. Apart from differentially regulated genes involved in PSI and PSII, we found widespread evidence of intron splicing events, specifically uniquely spliced petB and petD in Triticum aestivum and high proportion of RNA editing in ndh genes across the Triticeae members during salt stress. We also highlight the role and differential regulation of ATP synthase as member of CF₀CF₁ and also revealed the effect of salt stress on the water-splitting complex under salt stress. It is worthwhile to mention that the observed conserved down-regulation of psbJ across the Triticeae is limiting the assembly of water-splitting complexes and thus making the BEP clade Triticeae members more vulnerable to high light during the salt stress. Comparative understanding of the chloroplast transcriptional dynamics and photosynthetic regulation will improve the approaches for improved crop domestication.

A novel chloroplast gene reported for flagellate plants

PREMISE OF THE STUDY

Gene space in plant plastid genomes is well characterized and annotated, yet we discovered an unrecognized open reading frame (ORF) in the fern lineage that is conserved across flagellate plants.

METHODS

We initially detected a putative uncharacterized ORF by the existence of a highly conserved region between rps16 and matK in a series of matK alignments of leptosporangiate ferns. We mined available plastid genomes for this ORF, which we now refer to as ycf94, to infer evolutionary selection pressures and assist in functional prediction. To further examine the transcription of ycf94, we assembled the plastid genome and sequenced the transcriptome of the leptosporangiate fern Adiantum shastense Huiet & A.R. Sm.

KEY RESULTS

The ycf94 predicted protein has a distinct transmembrane domain but with no sequence homology to other proteins with known function. The nonsynonymous/synonymous substitution rate ratio of ycf94 is on par with other fern plastid protein-encoding genes, and additional homologs can be found in a few lycophyte, moss, hornwort, and liverwort plastid genomes. Homologs of ycf94 were not found in seed plants. In addition, we report a high level of RNA editing for ycf94 transcripts-a hallmark of protein-coding genes in fern plastomes.

CONCLUSIONS

The degree of sequence conservation, together with the presence of a distinct transmembrane domain and RNA-editing sites, suggests that ycf94 is a protein-coding gene of functional significance in ferns and, potentially, bryophytes and lycophytes. However, the origin and exact function of this gene require further investigation.

Editing of Chloroplast rps14 by PPR Editing Factor EMB2261 Is Essential for Arabidopsis Development

RNA editing in plastids is known to be required for embryogenesis, but no single editing event had been shown to be essential. We show that the emb2261-2 mutation is lethal through a failure to express an editing factor that specifically recognizes the rps14-2 site. EMB2261 was predicted to bind the cis-element upstream of the rps14-2 site and genetic complementation with promoters of different strength followed by RNA-seq analysis was conducted to test the correlation between rps14-2 editing and EMB2261 expression. Rps14-2 is the only editing event in Arabidopsis chloroplasts that correlates with EMB2261 expression. Sequence divergence between the cis-element and the EMB2261 protein sequence in plants where rps14-2 editing is not required adds support to the association between them. We conclude that EMB2261 is the specificity factor for rps14-2 editing. This editing event converts P51 in Rps14 to L51, which is conserved among species lacking RNA editing, implying the importance of the editing event to Rps14 function. Rps14 is an essential ribosomal subunit for plastid translation, which, in turn, is essential for Arabidopsis embryogenesis.

A Guide to the Chloroplast Transcriptome Analysis Using RNA-Seq

Since its first use in plants in 2007, high-throughput RNA sequencing (RNA-Seq) has generated a vast amount of data for both model and nonmodel species. Organellar transcriptomes, however, are virtually always overlooked at the data analysis step. We therefore developed ChloroSeq, a bioinformatic pipeline aimed at facilitating the systematic analysis of chloroplast RNA metabolism, and we provide here a step-by-step user's manual. Following the alignment of quality-controlled data to the genome of interest, ChloroSeq measures genome expression level along with splicing and RNA editing efficiencies. When used in combination with the Tuxedo suite (TopHat and Cufflinks), ChloroSeq allows the simultaneous analysis of organellar and nuclear transcriptomes, opening the way to a better understanding of nucleus-organelle cross talk. We also describe the use of R commands to produce publication-quality figures based on ChloroSeq outputs. The effectiveness of the pipeline is illustrated through analysis of an RNA-Seq dataset covering the transition from growth to maturation to senescence of Arabidopsis thaliana leaves.

Pervasive, Genome-Wide Transcription in the Organelle Genomes of Diverse Plastid-Bearing Protists

Organelle genomes are among the most sequenced kinds of chromosome. This is largely because they are small and widely used in molecular studies, but also because next-generation sequencing technologies made sequencing easier, faster, and cheaper. However, studies of organelle RNA have not kept pace with those of DNA, despite huge amounts of freely available eukaryotic RNA-sequencing (RNA-seq) data. Little is known about organelle transcription in nonmodel species, and most of the available eukaryotic RNA-seq data have not been mined for organelle transcripts. Here, we use publicly available RNA-seq experiments to investigate organelle transcription in 30 diverse plastid-bearing protists with varying organelle genomic architectures. Mapping RNA-seq data to organelle genomes revealed pervasive, genome-wide transcription, regardless of the taxonomic grouping, gene organization, or noncoding content. For every species analyzed, transcripts covered ≥85% of the mitochondrial and/or plastid genomes (all of which were ≤105 kb), indicating that most of the organelle DNA—coding and noncoding—is transcriptionally active. These results follow earlier studies of model species showing that organellar transcription is coupled and ubiquitous across the genome, requiring significant downstream processing of polycistronic transcripts. Our findings suggest that noncoding organelle DNA can be transcriptionally active, raising questions about the underlying function of these transcripts and underscoring the utility of publicly available RNA-seq data for recovering complete genome sequences. If pervasive transcription is also found in bigger organelle genomes (>105 kb) and across a broader range of eukaryotes, this could indicate that noncoding organelle RNAs are regulating fundamental processes within eukaryotic cells.

Two interacting PPR proteins are major Arabidopsis editing factors in plastid and mitochondria

RNA editing is converting hundreds of cytosines into uridines during organelle gene expression of land plants. The pentatricopeptide repeat (PPR) proteins are at the core of this posttranscriptional RNA modification. Even if a PPR protein defines the editing site, a DYW domain of the same or another PPR protein is believed to catalyze the deamination. To give insight into the organelle RNA editosome, we performed tandem affinity purification of the plastidial CHLOROPLAST BIOGENESIS 19 (CLB19) PPR editing factor. Two PPR proteins, dually targeted to mitochondria and chloroplasts, were identified as potential partners of CLB19. These two proteins, a P-type PPR and a member of a small PPR-DYW subfamily, were shown to interact in yeast. Insertional mutations resulted in embryo lethality that could be rescued by embryo-specific complementation. A transcriptome analysis of these complemented plants showed major editing defects in both organelles with a very high PPR type specificity, indicating that the two proteins are core members of E+-type PPR editosomes.

ORRM5, an RNA recognition motif-containing protein, has a unique effect on mitochondrial RNA editing

Plants have an RNA editing mechanism that prevents deleterious organelle mutations from resulting in impaired proteins. A typical flowering plant modifies about 40 cytidines in chloroplast transcripts and many hundreds of cytidines in mitochondrial transcripts. The plant editosome, the molecular machinery responsible for this process, contains members of several protein families, including the organelle RNA recognition motif (ORRM)-containing family. ORRM1 and ORRM6 are chloroplast editing factors, while ORRM2, ORRM3, and ORRM4 are mitochondrial editing factors. Here we report the identification of organelle RRM protein 5 (ORRM5) as a mitochondrial editing factor with a unique mode of action. Unlike other ORRM editing factors, the absence of ORRM5 in orrm5 mutant plants results in an increase of the editing extent in 14% of the mitochondrial sites surveyed. The orrm5 mutant also exhibits a reduced splicing efficiency of the first nad5 intron and slower growth and delayed flowering time. ORRM5 contains an RNA recognition motif (RRM) and a glycine-rich domain at the C terminus. The RRM provides the editing activity of ORRM5 and is able to complement the splicing but not the morphological defects.

The Roles of Mutation, Selection, and Expression in Determining Relative Rates of Evolution in Mitochondrial versus Nuclear Genomes

Eukaryotes rely on proteins encoded by the nuclear and mitochondrial (mt) genomes, which interact within multisubunit complexes such as oxidative-phosphorylation enzymes. Although selection is thought to be less efficient on the asexual mt genome, in bilaterian animals the ratio of nonsynonymous to synonymous substitutions (ω) is lower in mt- compared with nuclear-encoded OXPHOS subunits, suggesting stronger effects of purifying selection in the mt genome. Because high levels of gene expression constrain protein sequence evolution, one proposed resolution to this paradox is that mt genes are expressed more highly than nuclear genes. To test this hypothesis, we investigated expression and sequence evolution of mt and nuclear genes from 84 diverse eukaryotes that vary in mt gene content and mutation rate. We found that the relationship between mt and nuclear ω values varied dramatically across eukaryotes. In contrast, transcript abundance is consistently higher for mt genes than nuclear genes, regardless of which genes happen to be in the mt genome. Consequently, expression levels cannot be responsible for the differences in ω Rather, 84% of the variance in the ratio of ω values between mt and nuclear genes could be explained by differences in mutation rate between the two genomes. We relate these findings to the hypothesis that high rates of mt mutation select for compensatory changes in the nuclear genome. We also propose an explanation for why mt transcripts consistently outnumber their nuclear counterparts, with implications for mitonuclear protein imbalance and aging.