Unveiling Chloroplast RNA Editing Events Using Next Generation Small RNA Sequencing Data

Deciphering the enigma of RNA editing in the ATP1_alpha subunit of ATP synthase in Triticum aestivum

There is evidence that RNA editing is related to plant cellular stress as well as electron transport organelles, such as mitochondria. The mitochondrial atp1 gene encodes the alpha-subunit of Atp synthase. Control as well as two periods of drought stress treatments were analyzed in the cDNAs generated from the mitochondrial atp1 gene of two cultivars of Triticum aestivum [Giza 168 (G168) and Gemmiza 10 (GM10)]. Following RNA-seq data assembly, atp1 cDNAs from the control (acc. no. OQ129415), 2-hour (acc. no. OQ129416), and 12-hour (acc. no. OQ129417) time points of the T. aestivum cultivar G168 were obtained. Control (acc. no. OQ129419), 2-hour (acc. no. OQ129420), and 12-hour (acc. no. OQ129421) samples all included reconstructed atp1 transcripts from Gemmiza 10. Atp1 transcripts were assembled using the wheat atp1 gene (acc. no. NC_036024). RNA-seq raw data was utilized to identify 11 RNA editing sites in atp1 in the tolerant cultivar Giza168 and 6 in the sensitive cultivar Gemmiza10. The significant difference in RNA editing observed between control and drought stress conditions in sites led to synonymous amino acids. This led to no change in tertiary structure between tolerant and sensitive cultivars. But the change was focused between produced protein and its correspondence sequence on DNA.

The whole chloroplast genome of Mesua ferrea: Insight into the dynamic pattern of evolution and its comparison with species from recently diverged families

Mesua ferrea is an important source of timber, oil and herbal medicines. In the present investigation, we assembled the whole chloroplast genome of M. ferrea of size 161.4 kb. The genome contained 86 protein-coding genes, 38 tRNAs, 8 rRNA genes and exhibited a characteristic quadripartite structural orientation, with two inverted repeats (27,614 bp) separated by an LSC (88,746 bp) region and an SSC (27,614 bp) (17,470 bp). Interestingly, no gene loss was identified in the M. ferrea genome, contrary to what has been observed in other Clusioid species. We compared the chloroplast genome of M. ferrea with the chloroplast genome of Bonnetia and Garcinia belonging to Bonnetiaceae and Clusiaceae families. Overall, the compared genomes possess a similar synteny of gene order except for a small inversion in Garcinia species. M. ferrea has the largest chloroplast genome size in Clusioid clade owing to the lengthening of the LSC, IR, and non-coding regions. Substantial differences were observed in population of simple sequence repeats (SSRs) and RNA editing sites among the studied genomes. A comparative assessment of chloroplast genomes revealed five highly divergence regions: rpl32, trnS-GCU_trnG-UCC, petN-psbM, psbZ_trnG-GCC and ccsA_ndhD among the analyzed sequences. Phylogenetic analyses and sequence homology search indicate that M. ferrea is closely related to the Garcinia species.

Comparative Chloroplast Genomes of Nicotiana Species (Solanaceae): Insights Into the Genetic Variation, Phylogenetic Relationship, and Polyploid Speciation

Nicotiana L. is a genus rich in polyploidy, which represents an ideal natural system for investigating speciation, biodiversity, and phytogeography. Despite a wealth of phylogenetic work on this genus, a robust evolutionary framework with a dated molecular phylogeny for the genus is still lacking. In this study, the 19 complete chloroplast genomes of Nicotiana species were assembled, and five published chloroplast genomes of Nicotiana were retrieved for comparative analyses. The results showed that the 24 chloroplast genomes of Nicotiana, ranging from 155,327 bp (N. paniculata) to 156,142 bp (N. heterantha) in size, exhibited typical quadripartite structure. The chloroplast genomes were rather conserved in genome structure, GC content, RNA editing sites, and gene content and order. The higher GC content observed in the IR regions could be a result of the presence of abundant rRNA and tRNA genes, which contained a relatively higher GC content. A total of seven hypervariable regions, as new molecular markers for phylogenetic analysis, were uncovered. Based on 78 protein-coding genes, we constructed a well-supported phylogenetic tree, which was largely in agreement with previous studies, except for a slight conflict in several sections. Chloroplast phylogenetic results indicated that the progenitors of diploid N. sylvestris, N. knightiana, and the common ancestor of N. sylvestris and N. glauca might have donated the maternal genomes of allopolyploid N. tabacum, N. rustica, and section Repandae, respectively. Meanwhile, the diploid section Noctiflorae lineages (N. glauca) acted as the most likely maternal progenitor of section Suaveolentes. Molecular dating results show that the polyploid events range considerably in ~0.12 million (section Nicotiana) to ~5.77 million (section Repandae) years ago. The younger polyploids (N. tabacum and N. rustica) were estimated to have arisen ~0.120 and ~0.186 Mya, respectively. The older polyploids (section Repandae and Suaveolentes) were considered to have originated from a single polyploid event at ~5.77 and ~4.49 Mya, respectively. In summary, the comparative analysis of chloroplast genomes of Nicotiana species has not only revealed a series of new insights into the genetic variation and phylogenetic relationships in Nicotiana but also provided rich genetic resources for speciation and biodiversity research in the future.

Chloroplast Genome of Lithocarpus dealbatus (Hook.f. & Thomson ex Miq.) Rehder Establishes Monophyletic Origin of the Species and Reveals Mutational Hotspots with Taxon Delimitation Potential

There is phylogenetic ambiguity in the genus Lithocarpus and subfamily Quercoideae (Family: Fagaceae). Lithocarpus dealbatus, an ecologically important tree, is the dominant species among the Quercoideae in India. Although several studies have been conducted on the species' regeneration and ecological and economic significance, limited information is available on its phylo-genomics. To resolve the phylogeny in Quercoideae, we sequenced and assembled the 161,476 bp chloroplast genome of L. dealbatus, which has a large single-copy section of 90,732 bp and a small single-copy region of 18,987 bp, separated by a pair of inverted repeat regions of 25,879 bp. The chloroplast genome contained 133 genes, of which 86 were protein-coding genes, 39 were transfer RNAs, and eight were ribosomal RNAs. Analysis of repeat elements and RNA editing sites revealed interspecific similarities within the Lithocarpus genus. DNA diversity analysis identified five highly diverged coding and noncoding hotspot regions in the four genera, which can be used as polymorphic markers for species/taxon delimitation across the four genera of Quercoideae viz., Lithocarpus, Quercus, Castanea, and Castanopsis. The chloroplast-based phylogenetic analysis among the Quercoideae established a monophyletic origin of Lithocarpus, and a closer evolutionary lineage with a few Quercus species. Besides providing insights into the chloroplast genome architecture of L. dealbatus, the study identified five mutational hotspots having high taxon-delimitation potential across four genera of Quercoideae.

The first report of RNA U to C or G editing in the mitochondrial NADH dehydrogenase subunit 5 (Nad5) transcript of wild barley

BACKGROUND

Nad dehydrogenase complex in mtDNA has a significant role in cellular respiration. One of the largest subunits in the complex is subunit 5 (Nad5).

METHODS AND RESULTS

Four cDNAs of the Hordeum vulgare subsp. spontaneum nad5 gene have been characterized and subjected to four phases of 0.5 M salinity, at 0 h (control, accession no. MT235236), after 2 h (acc. no. MT235237), after 12 h (acc. no. MT235238) and after 24 h (acc. no. MT235239). Utilizing raw data from RNA-seq, ten RNA editing sites were reported. Seven sites have common editing from C to U in positions (C1490, C1859, C1895, C1900, C1901, C1916, C1918). A rare editing event U to C was detected in two positions (U1650 and U1652) and a novel editing event U to G was for the first time in positions nad5-U231. The highest editing level was shown in 2 and 12 h after salinity exposure. After 24 h, these edits were disrupted, possibly due to the launch of the programed cell death mechanism. However, the RNA editing in positions U1650, U1652 and U231 was fixed at all exposure times.

CONCLUSIONS

Although study clarified the role of salinity stress in nad5 RNA editing sites, the main achievements are first report of U to G RNA editing in plants at position U231 and first report of U to C editing in the nad5 gene at U1650 and U1652.

The Algal Chloroplast as a Testbed for Synthetic Biology Designs Aimed at Radically Rewiring Plant Metabolism

Sustainable and economically viable support for an ever-increasing global population requires a paradigm shift in agricultural productivity, including the application of biotechnology to generate future crop plants. Current genetic engineering approaches aimed at enhancing the photosynthetic efficiency or composition of the harvested tissues involve relatively simple manipulations of endogenous metabolism. However, radical rewiring of central metabolism using new-to-nature pathways, so-called "synthetic metabolism", may be needed to really bring about significant step changes. In many cases, this will require re-programming the metabolism of the chloroplast, or other plastids in non-green tissues, through a combination of chloroplast and nuclear engineering. However, current technologies for sophisticated chloroplast engineering ("transplastomics") of plants are limited to just a handful of species. Moreover, the testing of metabolic rewiring in the chloroplast of plant models is often impractical given their obligate phototrophy, the extended time needed to create stable non-chimeric transplastomic lines, and the technical challenges associated with regeneration of whole plants. In contrast, the unicellular green alga, Chlamydomonas reinhardtii is a facultative heterotroph that allows for extensive modification of chloroplast function, including non-photosynthetic designs. Moreover, chloroplast engineering in C. reinhardtii is facile, with the ability to generate novel lines in a matter of weeks, and a well-defined molecular toolbox allows for rapid iterations of the "Design-Build-Test-Learn" (DBTL) cycle of modern synthetic biology approaches. The recent development of combinatorial DNA assembly pipelines for designing and building transgene clusters, simple methods for marker-free delivery of these clusters into the chloroplast genome, and the pre-existing wealth of knowledge regarding chloroplast gene expression and regulation in C. reinhardtii further adds to the versatility of transplastomics using this organism. Herein, we review the inherent advantages of the algal chloroplast as a simple and tractable testbed for metabolic engineering designs, which could then be implemented in higher plants.

Light/heat effects on RNA editing in chloroplast NADH-plastoquinone oxidoreductase subunit 2 (ndhB) gene of Calotropis (Calotropis procera)

Background

RNA editing is common in terrestrial plants, especially in mitochondria and chloroplast. In the photosynthesis process, NAD dehydrogenase plays a very important role. Subunit 2 of NADH-dehydrogenase is one of the major subunits in NAD dehydrogenase complex. Using desert plant Calotropis (Calotropis procera), this study focuses on the RNA editing activity of ndhB based on light time.

Results

NdhB (NADH-dehydrogenase subunit 2) gene accession no. MK144329 was isolated from Calotropis procera genomic data (PRJNA292713). Additionally, using RNA-seq data, the cDNA of the ndhB gene of C. procera was isolated at three daylight periods, i.e., dawn (accession no. MK165161), at midday (accession no. MK165160), and pre-dusk (accession no. MK165159). Seven RNA editing sites have been found in several different positions (nucleotide no. C467, C586, C611, C737, C746, C830, and C1481) within the ndhB coding region. The rate of these alterations was deferentially edited across the three daylight periods. RNA editing rate of ndhB gene was highest at dawn, (87.5, 79.6, 78.5, 76, 68.6, 39.3, and 96.9%, respectively), less in midday (74.8, 54.1, 62.6, 47.4, 45.5, 47.4, and 93.4%, respectively), and less at pre-dusk (67, 52.6, 56.9, 40.1, 40.7, 33.2, and 90%, respectively), also all these sites were validated by qRT-PCR.

Conclusion

The differential editing of chloroplast ndhB gene across light periods may be led to a somehow relations between the RNA editing and control of photosynthesis.

Plastid genomics of Nicotiana (Solanaceae): insights into molecular evolution, positive selection and the origin of the maternal genome of Aztec tobacco (Nicotiana rustica)

Species of the genus Nicotiana (Solanaceae), commonly referred to as tobacco plants, are often cultivated as non-food crops and garden ornamentals. In addition to the worldwide production of tobacco leaves, they are also used as evolutionary model systems due to their complex development history tangled by polyploidy and hybridization. Here, we assembled the plastid genomes of five tobacco species: N. knightiana, N. rustica, N. paniculata, N. obtusifolia and N. glauca. De novo assembled tobacco plastid genomes had the typical quadripartite structure, consisting of a pair of inverted repeat (IR) regions (25,323-25,369 bp each) separated by a large single-copy (LSC) region (86,510-86,716 bp) and a small single-copy (SSC) region (18,441-18,555 bp). Comparative analyses of Nicotiana plastid genomes with currently available Solanaceae genome sequences showed similar GC and gene content, codon usage, simple sequence and oligonucleotide repeats, RNA editing sites, and substitutions. We identified 20 highly polymorphic regions, mostly belonging to intergenic spacer regions (IGS), which could be suitable for the development of robust and cost-effective markers for inferring the phylogeny of the genus Nicotiana and family Solanaceae. Our comparative plastid genome analysis revealed that the maternal parent of the tetraploid N. rustica was the common ancestor of N. paniculata and N. knightiana, and the later species is more closely related to N. rustica. Relaxed molecular clock analyses estimated the speciation event between N. rustica and N. knightiana appeared 0.56 Ma (HPD 0.65-0.46). Biogeographical analysis supported a south-to-north range expansion and diversification for N. rustica and related species, where N. undulata and N. paniculata evolved in North/Central Peru, while N. rustica developed in Southern Peru and separated from N. knightiana, which adapted to the Southern coastal climatic regimes. We further inspected selective pressure on protein-coding genes among tobacco species to determine if this adaptation process affected the evolution of plastid genes. These analyses indicate that four genes involved in different plastid functions, including DNA replication (rpoA) and photosynthesis (atpB, ndhD and ndhF), came under positive selective pressure as a result of specific environmental conditions. Genetic mutations in these genes might have contributed to better survival and superior adaptations during the evolutionary history of tobacco species.

Identification of soybean trans-factors associated with plastid RNA editing sites

RNA editing is a posttranscriptional process that changes nucleotide sequences, among which cytosine-to-uracil by a deamination reaction can revert non-neutral codon mutations. Pentatricopeptide repeat (PPR) proteins comprise a family of RNA-binding proteins, with members acting as editing trans-factors that recognize specific RNA cis-elements and perform the deamination reaction. PPR proteins are classified into P and PLS subfamilies. In this work, we have designed RNA biotinylated probes based in soybean plastid RNA editing sites to perform trans-factor specific protein isolation. Soybean cis-elements from these three different RNA probes show differences in respect to other species. Pulldown samples were submitted to mass spectrometry for protein identification. Among detected proteins, five corresponded to PPR proteins. More than one PPR protein, with distinct functional domains, was pulled down with each one of the RNA probes. Comparison of the soybean PPR proteins to Arabidopsis allowed identification of the closest homologous. Differential gene expression analysis demonstrated that the PPR locus Glyma.02G174500 doubled its expression under salt stress, which correlates with the increase of its potential rps14 editing. The present study represents the first identification of RNA editing trans-factors in soybean. Data also indicated that potential multiple trans-factors should interact with RNA cis-elements to perform the RNA editing.

Salinity effects on nad3 gene RNA editing of wild barley mitochondria

Nad complex plays a very important role during cellular respiration. nad3 (nad dehydrogenase subunit 3) is one of the biggest subunits in this complex. Four cDNAs of nad3 gene were characterized in Hordeum vulgare subsp. spontaneum at exposed to four periods of 500 mM salinity, 0 h or control (accession no. MN066165), after 2 h (accession no. MN066166), after 12 h (accession no. MN066167) and after 24 h (accession no. MN066168) using RNA-seq raw data. Seventeen RNA editing sites were found in positions (or nucleotide nos. C5, C39, C44, C61, C62, C79, C80, C147, C185, C190, C191, C208, C209, C275, C317, C344, C349) within the nad3 coding region. These alterations represent differential editing at four exposure times. The maximum editing rate was revealed 2 and 12 h after salinity exposure. However, these edits were disrupted after 24 h probably due to the initiation of program cell death machinery. We found that RNA editing not only improved protein function but also may improve codon bias by altering the nucleotide without any change in amino acid. Characterization of pentatricopeptide repeat-containing protein At4g13650 (PPRSp1) in wild barley helped us to understand the behavior of editing sites C190 and C191 under salinity. Position - 6 in cis-element upstream editing sites of C155, C190 and C191 may be vital to the editing process in these sites by PPRSp1 protein. The differential editing of this gene under salinity led to a relationship between RNA editing and cellular respiration regulation.

Comparison of different annotation tools for characterization of the complete chloroplast genome of Corylus avellana cv Tombul

BACKGROUND

Several bioinformatics tools have been designed for assembly and annotation of chloroplast (cp) genomes, making it difficult to decide which is most useful and applicable to a specific case. The increasing number of plant genomes provide an opportunity to accurately obtain cp genomes from whole genome shotgun (WGS) sequences. Due to the limited genetic information available for European hazelnut (Corylus avellana L.) and as part of a genome sequencing project, we analyzed the complete chloroplast genome of the cultivar 'Tombul' with multiple annotation tools.

RESULTS

Three different annotation strategies were tested, and the complete cp genome of C. avellana cv Tombul was constructed, which was 161,667 bp in length, and had a typical quadripartite structure. A large single copy (LSC) region of 90,198 bp and a small single copy (SSC) region of 18,733 bp were separated by a pair of inverted repeat (IR) regions of 26,368 bp. In total, 125 predicted functional genes were annotated, including 76 protein-coding, 25 tRNA, and 4 rRNA unique genes. Comparative genomics indicated that the cp genome sequences were relatively highly conserved in species belonging to the same order. However, there were still some variations, especially in intergenic regions, that could be used as molecular markers for analyses of phylogeny and plant identification. Simple sequence repeat (SSR) analysis showed that there were 83 SSRs in the cp genome of cv Tombul. Phylogenetic analysis suggested that C. avellana cv Tombul had a close affinity to the sister group of C. fargesii and C. chinensis, and then a closer evolutionary relationship with Betulaceae family than other species of Fagales.

CONCLUSION

In this study, the complete cp genome of Corylus avellana cv Tombul, the most widely cultivated variety in Turkey, was obtained and annotated, and additionally phylogenetic relationships were predicted among Fagales species. Our results suggest a very accurate assembly of chloroplast genome from next generation whole genome shotgun (WGS) sequences. Enhancement of taxon sampling in Corylus species provide genomic insights into phylogenetic analyses. The nucleotide sequences of cv Tombul cp genomes can provide comprehensive genetic insight into the evolution of genus Corylus.

Potential of Transcript Editing Across Mitogenomes of Early Land Plants Shows Novel and Familiar Trends

RNA editing alters the identity of nucleotides in an RNA sequence so that the mature transcript differs from the template defined in the genome. This process has been observed in chloroplasts and mitochondria of both seed and early land plants. However, the frequency of RNA editing in plant mitochondria ranges from zero to thousands of editing sites. To date, analyses of RNA editing in mitochondria of early land plants have been conducted on a small number of genes or mitochondrial genomes of a single species. This study provides an overview of the mitogenomic RNA editing potential of the main lineages of these two groups of early land plants by predicting the RNA editing sites of 33 mitochondrial genes of 37 species of liverworts and mosses. For the purpose of the research, we newly assembled seven mitochondrial genomes of liverworts. The total number of liverwort genera with known complete mitogenome sequences has doubled and, as a result, the available complete mitogenome sequences now span almost all orders of liverworts. The RNA editing site predictions revealed that C-to-U RNA editing in liverworts and mosses is group-specific. This is especially evident in the case of liverwort lineages. The average level of C-to-U RNA editing appears to be over three times higher in liverworts than in mosses, while the C-to-U editing frequency of the majority of genes seems to be consistent for each gene across bryophytes.

Post-transcriptional regulation of the oxidative stress response in plants

Due to their sessile lifestyle, plants can be exposed to several kinds of stresses that will increase the production of reactive oxygen species (ROS), such as hydrogen peroxide, singlet oxygen, and hydroxyl radicals, in the plant cells and activate several signaling pathways that cause alterations in the cellular metabolism. Nevertheless, when ROS production outreaches a certain level, oxidative damage to nucleic acids, lipids, metabolites, and proteins will occur, finally leading to cell death. Until now, the most comprehensive and detailed readout of oxidative stress responses is undoubtedly obtained at the transcriptome level. However, transcript levels often do not correlate with the corresponding protein levels. Indeed, together with transcriptional regulations, post-transcriptional, translational, and/or post-translational regulations will shape the active proteome. Here, we review the current knowledge on the post-transcriptional gene regulation during the oxidative stress responses in planta.

Expected and unexpected evolution of plant RNA editing factors CLB19, CRR28 and RARE1: retention of CLB19 despite a phylogenetically deep loss of its two known editing targets in Poaceae

BACKGROUND

C-to-U RNA editing in mitochondria and chloroplasts and the nuclear-encoded, RNA-binding PPR proteins acting as editing factors present a wide field of co-evolution between the different genetic systems in a plant cell. Recent studies on chloroplast editing factors RARE1 and CRR28 addressing one or two chloroplast editing sites, respectively, found them strictly conserved among 65 flowering plants as long as one of their RNA editing targets remained present.

RESULTS

Extending the earlier sampling to 117 angiosperms with high-quality genome or transcriptome data, we find more evidence confirming previous conclusions but now also identify cases for expected evolutionary transition states such as retention of RARE1 despite loss of its editing target or the degeneration of CRR28 truncating its carboxyterminal DYW domain. The extended angiosperm set was now used to explore CLB19, an "E+"-type PPR editing factor targeting two chloroplast editing sites, rpoAeU200SF and clpPeU559HY, in Arabidopsis thaliana. We found CLB19 consistently conserved if one of the two targets was retained and three independent losses of CLB19 after elimination of both targets. The Ericales show independent regains of the ancestrally lost clpPeU559HY editing, further explaining why multiple-target editing factors are lost much more rarely than single target factors like RARE1. The retention of CLB19 despite loss of both editing targets in some Ericaceae, Apocynaceae and in Camptotheca (Nyssaceae) likely represents evolutionary transitions. However, the retention of CLB19 after a phylogenetic deep loss in the Poaceae rather suggests a yet unrecognized further editing target, for which we suggest editing event ndhAeU473SL.

CONCLUSION

Extending the scope of studies on plant organelle RNA editing to further taxa and additional nuclear cofactors reveals expected evolutionary transitions, strikingly different evolutionary dynamics for multiple-target editing factors like CLB19 and CRR28 and suggests additional functions for editing factor CLB19 among the Poaceae.