Next Article in Journal
mRNA Expression and DNA Methylation Analysis of Serotonin Receptor 2A (HTR2A) in the Human Schizophrenic Brain
Next Article in Special Issue
Applying Human ADAR1p110 and ADAR1p150 for Site-Directed RNA Editing—G/C Substitution Stabilizes GuideRNAs against Editing
Previous Article in Journal
Genetic Regulation of Virulence and Antibiotic Resistance in Acinetobacter baumannii
Previous Article in Special Issue
RNA Editing and Its Molecular Mechanism in Plant Organelles
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 8
Issue 1
10.3390/genes8010013
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification and Analysis of RNA Editing Sites in the Chloroplast Transcripts of Aegilops tauschii L.

by
Mengxing Wang
,
Hui Liu
,
Lingqiao Ge
,
Guangwei Xing
,
Meng Wang
,
Song Weining
* and
Xiaojun Nie
*
State Key Laboratory of Crop Stress Biology in Arid Areas, College of Agronomy, Northwest Agriculture and Forestry University, Yangling 712100, Xianyang, China
*
Authors to whom correspondence should be addressed.
Genes 2017, 8(1), 13; https://0-doi-org.brum.beds.ac.uk/10.3390/genes8010013
Submission received: 23 September 2016 / Revised: 9 December 2016 / Accepted: 20 December 2016 / Published: 30 December 2016
(This article belongs to the Special Issue RNA Editing)
Browse Figures
Versions Notes

Abstract

:
RNA editing is an important way to convert cytidine (C) to uridine (U) at specific sites within RNA molecules at a post-transcriptional level in the chloroplasts of higher plants. Although it has been systematically studied in many plants, little is known about RNA editing in the wheat D genome donor Aegilops tauschii L. Here, we investigated the chloroplast RNA editing of Ae. tauschii and compared it with other wheat relatives to trace the evolution of wheat. Through bioinformatics prediction, a total of 34 C-to-U editing sites were identified, 17 of which were validated using RT-PCR product sequencing. Furthermore, 60 sites were found by the RNA-Seq read mapping approach, 24 of which agreed with the prediction and six were validated experimentally. The editing sites were biased toward tCn or nCa trinucleotides and 5′-pyrimidines, which were consistent with the flanking bases of editing sites of other seed plants. Furthermore, the editing events could result in the alteration of the secondary structures and topologies of the corresponding proteins, suggesting that RNA editing might impact the function of target genes. Finally, comparative analysis found some evolutionarily conserved editing sites in wheat and two species-specific sites were also obtained. This study is the first to report on RNA editing in Aegilops tauschii L, which not only sheds light on the evolution of wheat from the point of view of RNA editing, but also lays a foundation for further studies to identify the mechanisms of C-to-U alterations.

1. Introduction

RNA editing is a post-transcriptional process for converting cytidine (C) to uridine (U), or U to C, at specific sites within RNA molecules to alter the identity of nucleotides between RNA and genomic DNA, serving as a mechanism to correct missense mutations of genes at the RNA level, and also to enrich genetic information [1,2]. It can act upon transcripts in the cell nucleus, cytoplasm, mitochondria, and chloroplast (cp), which has been widely documented in various organisms, including eukaryotes, prokaryotes, archaebacteria, viruses, etc. [3,4,5]. Since the first editing event of the land plant plastome was found in the mRNA transcript of the rpl2 gene in the maize chloroplast genome in 1991 [6], extensive RNA editing sites have been detected in a variety of higher plants, such as tobacco [7,8], Arabidopsis [9], pea [10], and tomato [11], and some of them have been validated by experiments [12]. Hornwort and fern exhibit hundreds of C-to-U as well as U-to-C editing sites [13,14]; however, in angiosperms, approximately 30 to 40 C-U editing sites are found in the chloroplasts [15]. In other land plant clades, editing frequencies vary widely, with no editing in marchantiid liverworts [16]. Although a few cases of editing within the noncoding regions or introns have been reported, most of the known editing events occur in exons, restoring corresponding amino acids to a more conserved direction [17,18]. Additionally, editing usually occurs in the second or first base of codons, resulting in the conversion of hydrophilic amino acid to hydrophobic, in most cases [19,20]. As an essential step of RNA maturation in plant organelles, the molecular mechanisms of RNA editing have been extensively studied [21]. It is widely accepted that PPR (pentatricopeptide repeat) proteins are the key factor to regulate RNA editing through binding to the target transcripts and then altering the RNA sequence turnover, processing, and translation [22,23,24].
Aegilops tauschii L., also known as Tausch’s goatgrass, is an annual grass belonging to the Triticeae tribe of the Poaceae family, which is found all over the world, including in China, Central and Western Asia, the Caucasus, Eastern Europe, and America. Ae. tauschii donated the D genome of hexaploid wheat (Triticum aestivum L.) by a spontaneous hybridization of the wild diploid grass Ae. tauschii with the tetraploid wheat Triticum turgidum about 8000–10,000 years ago in the Fertile Crescent [25,26], which holds promise for the underlying processes of wheat formation and evolution. Extensive studies have been conducted to reveal the roles of Ae. tauschii in the evolutionary origin of polyploid wheat from different viewpoints, such as complete genome sequencing [27], physical mapping [28], and plastome sequences [29], as well as population genetic analyses [30]. However, the RNA editing sites in Ae. tauschii have not been well understood up to now. To trace the evolutionary process of polyploid wheat from the perspective of chloroplast RNA editing, we investigated and compared the chloroplast RNA editing patterns of Ae. tauschii together with other wheat relative species in this study. Firstly, the RNA editing sites in 76 chloroplast transcripts of Ae. tauschii were identified using bioinformatics prediction, combined with the RNA-Seq read mapping approach. Then, these predicted editing sites were further validated through PCR and RT-PCR. Furthermore, the alterations of the secondary structures of these proteins were investigated to understand the potential biological functions of these editing events. Finally, the chloroplast RNA editing sites in Ae. tauschii were compared with those of common wheat, wheat A and B genome progenitors, as well as five other Poaceae species to provide information regarding the origin and evolution of chloroplast RNA editing in wheat and beyond.

2. Materials and Methods

2.1. Prediction of RNA Editing Sites

All of the 76 protein-coding genes in the plastome of Ae. tauschii (AL8/78) (GenBank Accession: KJ614412.1) were downloaded from the NCBI (National Center for Biotechnology Information) database (https://www.ncbi.nlm.nih.gov) according to their annotation information [29]. Potential RNA editing sites in these genes were predicted by the Predictive RNA Editor for Plants (PREP)-Cp web server (http://prep.unl.edu/cgi-bin/cp-input.pl) [31] with a cutoff value of 0.8. Bioedit [32] and Microsoft Excel tools were used to format the required input file.

2.2. Editing Site Detection by Read Mapping and SNP (Single Nucleotide Polymorphism) Calling

RNA-Seq data from eight tissues of Ae. tauschii (AL8/78) were downloaded from the Sequence Read Archive (SRA) database with the accession numbers of SRX209402~SRX209425 (Table S1) [27]. The RNA-Seq reads were mapped to the cp genome of Ae. tauschii (AL8/78) using the Tophat tool [33]. SNPs were called by SAMtools to detect potential changes of C to U in the sense-strand and G (guanine) to A (adenine) in the antisense-strand of the cp genome. The mapped reads were subsequently visualized in the Integrative Genomics Viewer [34]. Only the DNA:RNA mismatches that had at least 5% editing efficiency, and more than five individual read supports were kept. RNA editing efficiency was counted by edited reads divided by total mapped reads [2].

2.3. Protein Structure Analysis before and after Editing

CFSSP (Chou and Fasman Secondary Structure Prediction Server) (http://www.biogem.org/tool/chou-fasman/) [35,36] was used to analyze the protein secondary structures of amino acid sequences translated from edited and unedited RNA. Prediction of transmembrane segments before and after editing, was performed using the Phobius tool (http://phobius.sbc.su.se/) [37].

2.4. Comparative Analysis of RNA Editing Sites in Cp Protein-Coding Genes among Poaceae Species

Previous study has been conducted to reveal the origin and evolution of RNA editing in cereal mitochondria [38]. Similarly, we compared and analyzed RNA editing in chloroplasts of Poaceae species. Firstly, the RNA editing sites in cp protein-coding genes of Ae. tauschii were compared to those of Hordeum vulgare [39], Lolium perenne [40], Oryza sativa [41], Saccharum officinarum [42], and Zea mays [17] to identify the patterns of RNA editing in Poaceae plastomes. Furthermore, all of the chloroplast protein-coding genes of Triticum urartu, Aegilops speltoides, and T. turgidum, as well as T. aestivum (Table 1) [29], were downloaded and potential RNA editing sites in these genes were predicted using the methods described above, and then the differences among these species were investigated.

2.5. Validation of Partial RNA Editing Sites

Seeds of Ae. tauschii (AL8/78) were put on moistened filter paper in a dish and then cultured in a growth chamber under the condition of 20 °C and a natural light/dark cycle. One-week-old seedlings were harvested for DNA and RNA isolation. Genomic DNA was extracted, as described previously [43]. Total RNA was isolated from the seedlings using the E.Z.N.A.® Plant RNA Kit (Omega, Norcross, GA, USA) and then treated with HiScriptTM Q RT SuperMix (+gDNA wiper) (Vazyme, Nanjing, China) to remove genomic DNA and to synthesize cDNAs, according to the manufacturer’s instructions.
The extracted DNA and obtained cDNAs were used as PCR templates for the validation of specific RNA editing sites, namely 18 sites of six genes, selected randomly from the prediction, as well as six sites in six genes identified by RNA-seq read mapping. In total, 11 specific primers (Table S2) were designed according to the candidate regions (Table S3). The PCR reaction (20 μL) contained 10 μM of each primer, 1 × EsTaq Mix (Cwbiotech, Beijing, China) and up to 1000 ng of template. The reaction proceeded through 35 cycles after an initial denaturation for 2 min at 94 °C, and was followed by a final extension for 5 min at 72 °C. One thermo cycle was 94 °C for 30 s, annealing for 30 s, and 72 °C for extension. The annealing temperature and extension time were adjusted to the Tm value of the primer pair and the length of product, respectively. The resulting products were separated by 1% agarose gel electrophoresis and were purified with a Universal DNA Purification Kit (Tiangen, Beijing, China), following the manufacturer’s instructions. The homogeneous products of petL, ndhK, rpoA, rps3, ndhD, ndhG, and rpoB were sequenced directly from both directions by Sangon Biotech Corporation (Shanghai, China). The products of ndhB, atpA, matK, ndhA, and rpoB were subsequently ligated into pMD18-T Vector (Takara, Shiga, Japan) and transformed into Escherichia coli DH5α. Two positive clones were selected for each fragment to sequence (both directions for ndhB and only one direction for others) by Sangon Biotech Corporation. Both DNA and cDNA sequences were spliced using the ContigExpress tool implemented in Vector NTI suite 6.0 (Invitrogen, Waltham, MA, USA), and the RNA editing sites were validated by aligning cDNAs with corresponding genomic DNA sequences using the ClustalW program (http://www.ebi.ac.uk/Tools/msa/clustalw2/) [44].

3. Results and Discussion

3.1. Prediction of RNA Editing Sites

In total, 34 editing sites present in 15 chloroplast protein-coding genes (Table 2) were predicted by the PREP-Cp program in Ae. tauschii, all of which were C-to-U conversions. Among them, ndh genes showed the most abundant editing sites with a total of 15, followed by rpo genes with a value of 10, and ycf with two editing events. Furthermore, 30 out of 34 (88.24%) editing events were observed in the second position of the codon and the rest were in the first position.

3.2. Editing Site Detection Using Read Mapping and SNP Calling

A total of 60 sites were identified as RNA editing sites by read mapping and SNP calling, 24 of which agreed with the prediction (Table 3). Among them, only 11 sites were found to be completely edited, and the remaining sites were partially edited, accounting for 81% (49 of 60), which was lower than the 100% value of A. thaliana [45]. The abundant, partially edited transcripts, in which not all of the C targets had been converted to Us, might be an intermediate phase in the editing process [15] and might result in more than one protein product encoded by the same gene, suggesting that chloroplast RNA editing was a post-transcriptional process [4]. Furthermore, 20 sites were found to be localized in the intergenic region, presumably affecting 5′ or 3′ UTRs, and two were located within the intron of rps16.
It has been demonstrated that the editing sites in the chloroplasts of seed plants preferred 5′-T–A-3′ flanking bases and/or a 5′-pyrimidine due to the genomic context with lower substitution rates [46,47,48]. Thus, we further analyzed the flanking bases of all 70 editing sites identified via the bioinformatics prediction and RNA-Seq data. Results revealed that 31.43% (22 of 70) of sites showed a tCa bias and 81.43% (57 of 70) possessed tCn or nCa trinucleotides (Figure 1). Additionally, a 5′-pyrimidine was observed in 61 sites, accounting for 87.14%, of which thymine accounted for 57.14% and cytosine accounted for 30.00%. Furthermore, all the editing events in the protein-coding genes could result in eight types of amino acid transitions (Figure 2). Among them, 27 out of 40 editing events (67.5%) led to a conversion from hydrophilic amino acid to hydrophobic amino acid (Figure 2, red bar), with only one changing from hydrophobic amino acid to hydrophilic amino acid (P→S) (Figure 2, blue bar), and the remaining 12 without hydrophobicity exchanges (Figure 2, green bar). Moreover, most of the amino acid conversions affected by RNA editing were serine to leucine (17), followed by proline to leucine (11). This was most likely due to serine being encoded by tCn and proline by the cCn codon. The low re-mutation biased genomic context could result in the conservative (functionally similar) amino acid exchange, probably without negative effects on protein function [49].
Additionally, six silent editing sites with no impact on the encoded amino acid were also detected, including pistil-specific atpA-234 and rpl16-36, spike-specific psbB-414, and rps3-30 edited in all tissues. Editing events at sites rps18-448 and ndhF-1891 generated stop codons, resulting in the premature termination of transcription, which only occurred in pistils with relatively low efficiency (9%). Furthermore, editing events at ndhB-586, ndhB-611, ndhB-704, and ndhB-1481 on the sense-strand, as well as ndhB-836, ndhB-149, ndhA-1070, and ndhA-563 on the antisense-strand, were found to only occur in leaves, while rps4-305 was specifically edited in seeds.

3.3. Validation of the RNA Editing Sites by RT-PCR Analysis

To validate the predicted RNA editing sites, 24 sites of 11 genes were randomly selected for PCR and RT-PCR analyses. By comparing cDNAs with genomic sequences, 23 sites in 11 genes were proven to be true (approximately 95.83%). Only rpoB-617 was not validated, as nucleotide C was present in both the mRNA and genomic DNA (Figure S1). In addition to false positives resulting from the bioinformatics prediction, one possible explanation was that rpoB-617 was edited in a tissue- or stage-specific manner, or modulated under different environmental conditions, which were not included in current study.
Generally, there are three methods to identify RNA editing sites, namely bioinformatics prediction and RNA-Seq, as well as molecular cloning. Molecular cloning is a common, widely used and accurate method to identify RNA editing in different species. Kugita et al. [50] verified a total of 509 C-to-U and 433 U-to-C conversions in hornwort chloroplasts and showed that RNA editing in hornwort chloroplasts made more than half the genes function; Jiang et al. [51] identified 54 editing sites in the cotton chloroplast genome, which was the highest editing frequency reported in angiosperms. However, this method requires a series of time-consuming and tedious experimental procedures, making it not suitable for genome-wide identification. In light of the evolutionary conservation, bioinformatics prediction tools have been developed to identify RNA editing, such as PREP-Cp and CURE (Cytidine to Uridine Recognizing Editor) [52]. In the cp genome of Cycas taitungensis, 85 editing sites were found in 25 transcripts using a combination of prediction by CURE and experimental determination [53]. In coconut cpDNA, RT-PCR analysis confirmed editing at 64 of the 83 PREP-Cp predicted RNA editing sites out of 27 genes [54]. With the advent of high-throughput sequencing, RNA-Seq technology provides an opportunity to investigate RNA editing sites at the transcriptome level. By read mapping and SNP calling, SNPs that emerged repeatedly and existed stable could be interpreted as edited sites. It has been gradually used to identify RNA editing on a large scale [2,55]. However, edited targets can be easily neglected, or SNPs can be falsely interpreted as true edited targets using this method [55]. Thus, to efficiently and accurately identify RNA editing sites in Aegilops tauschii L. chloroplasts, we performed identification through a combination of these three methods, which provide a systematic and comprehensive RNA editing profile of the transcriptome of Aegilops tauschii L. chloroplasts.

3.4. Impact of RNA Editing on Protein Structure

To preliminarily understand the biological functions of RNA editing, the secondary structures of these proteins, encoded by 15 genes with editing events, were analyzed using CFSSP and Phobius tools. Results showed that petB, rpoA, and rps8 possessed the same secondary structures, before and after editing, while editing events of the remaining 12 genes could alter the composition of the secondary structures around the editing sites (Table S4). More than half of the edited genes preferred to form new a α-helix and β-sheet, and also more than half of them underwent a decrease of turn structures at up- and down-stream regions around the editing codon (Table S5).
In the cp genome of Chlamydomonas reinhardtii, no editing was observed at position 204 of the petB gene encoding a cytochrome b6 subunit after introducing a proline residue in place of a leucine residue [56]. This gene modification led to a defective assembly of cytochrome b6f complexes, consistent with a block in the photosynthetic electron transfer so that the transformants displayed non-phototrophic phenotypes. The 204 position of petB in C. reinhardtii corresponded to the petB-662 in Ae. tauschii (AL8/78), which was in the proline codon CCA at the DNA level and was subsequently edited to the leucine codon CTA at the RNA level. Although the editing event at petB-662 had no impact on the hydrophobicity of the amino acids, the composition of the secondary structures or the number of transmembrane structures, the transmembrane region at codons 201 to 220 expanded to the C terminal of the protein (Figure 3A). We thus speculated that the editing event at petB-662 contributed to maintaining the functional activity of cytochrome b6f in Ae. tauschii.
There was no change in the number of transmembrane structures of ndhA, before and after editing, which encoded NADH dehydrogenase subunit A. Nevertheless, a contraction in the region of codons 296 to 323 was observed and the transmembrane region at codons 335 to 353 shifted to the C terminal (Figure 3B). The edited codons also restored universally conserved amino acids in the ndhA-encoded peptides of other chloroplast species, as was found in maize [57]. In addition, ndhB, with all sites edited, revealed a new transmembrane region at codons 449 to 510 (Figure 3C). NdhB editing site III was reported to remain unedited in the non-photosynthetic mutants of tobacco, as well as during the etiolated stage of the seedling development of maize, whereas it was fully edited in plants with an intact photosynthetic apparatus, indicating that the editing events in ndhB of Ae. tauschii might occur with a similar developmental stage-dependent manner, suggesting an unexpected link between ndhB editing and photosynthesis and light-induced chloroplast development [58].

3.5. Comparative Analysis of Chloroplast RNA Editing Sites among Poaceae Species

A total of 31 conserved RNA editing sites in 15 chloroplast protein-coding genes were systematically compared in Ae. tauschii, H. vulgare, L. perenne, O. sativa, S. officinarum, and Z. mays, 14 sites of which were found to be fully or partially edited (Table 4), establishing an evolutionarily conserved feature among these six species. No editing events occurred at ndhA-1, ndhK-1, rpl2-1 and rps14-1, and the absence of edited rpl2-1 was specific to Ae. tauschii. NdhK-1 was not edited in Ae. tauschii, and nucleotide C was present at the genomic DNA level. For the other three sites, nucleotide T already existed at the genomic DNA level. These sites had probably re-mutated to T at the DNA level to remedy detrimental genomic mutations that occurred after the water-to-land transition of land plants [59]. In terms of ndhB-1, ndhB-5 and ndhB-8, Ae. tauschii, H. vulgare, L. perenne and O. sativa shared edited ndhB-5 and ndhB-8, whereas ndhB-1 was only edited in Triticeae crops, suggesting that patterns of RNA editing were in accordance with the phylogenetic relationship (Figure S2) and there was a specific evolutionary force to drive the RNA editing events in these species.
Furthermore, to trace the evolution information of polyploid wheat from the chloroplast RNA editing point of view, we investigated and compared the chloroplast RNA editing patterns between Ae. tauschii and other wheat relatives including T. urartu, Ae. speltoides and T. turgidum. In total, 35 editing sites were detected in Ae. speltoides, T. turgidum, and T. aestivum, while only 33 sites were found in T. urartu (Tables S6–S9). Moreover, common wheat and its genome donors shared the 33 sites in 15 protein-coding genes of T. urartu (Figure 4), indicating that a highly conserved characteristic of RNA editing was present during the evolution of common wheat. It is worth noting that atpB-35 and ndhB-737 were edited in Ae. speltoides and T. turgidum compared to T. urartu (Figure 4A). As atpB-35 and ndhB-737 in T. urartu, as well as atpB-35 in Ae. tauschii, were not edited, the encoded T in the DNA level might be a result of re-mutation to repair the deficient transcripts that arose from mutations as described previously. If we posited that editing sites remained stable when the genome progenitors evolved to common wheat, undergoing the natural cross and chromosome doubling, atpB-35 and ndhB-737 of T. turgidum might come from Ae. speltoides. As predicted, the editing event at atpB-35 was not observed in Ae. tauschii but occurred in T. aestivum (Figure 4B,C), suggesting that atpB-35 of T. aestivum came from T. turgidum and originated from Ae. speltoides. However, the editing event at ndhB-737 also occurred in Ae. tauschii. Thus, further study is needed to reveal the origin of ndhB-737 in common wheat.

4. Conclusions

Here, we identified the RNA editing sites in the Ae. tauschii chloroplast genome through bioinformatics prediction, combined with RNA-Seq mapping and molecular cloning. A total of 34 editing sites in 15 protein-coding genes were predicted by PREP-Cp, of which 23 sites in 12 genes were validated by RNA-Seq mapping and 17 sites in six genes by molecular cloning. Additionally, RNA-Seq mapping detected 60 sites, seven of which were confirmed experimentally. Ae. tauschii shared similar patterns and features of RNA editing with other land plants and editing events could result in the alteration of the secondary structures of the encoded protein, which might have an impact on biological functions. Finally, comparative analysis was performed to identify the evolutionarily conserved or species-specific sites. Two sites, including atpB-35 and ndhB-737, were found to underlie some origin and evolutionary information of common wheat. This study, for the first time, reported on RNA editing in Aegilops tauschii L. and compared the RNA editing pattern among wheat relative species, which not only shed light on the evolution of wheat from the point of view of RNA editing, but also lay the foundation for further studies to identify the mechanisms of RNA editing in wheat and beyond.

Supplementary Materials

The following are available online at www.mdpi.com/2073-4425/8/1/13/s1, Figure S1: RNA editing sites confirmed experimentally, except rpoB-617 which was not edited. Arrow points to the editing site; Figure S2: The phylogeny of the compared six Poaceae species, Ae. tauschii, H. vulgare, L. perenne, O. sativa, S. officinarum and Z. mays. The phylogenetic tree of the compared six Poaceae species was reconstructed using maximum likelihood (ML) based on the RNA editing sites’ composition; Table S1: Information of RNA-Seq data of Ae. tauschii (AL8/78); Table S2: Primer information used in this study; Table S3: RNA editing sites validated by each primer pair; Table S4: Quantity of protein secondary structures before and after editing; Table S5: Conversion of protein secondary structure before and after editing; Table S6: RNA editing sites in chloroplast genes of T. aestivum cv CS (TA3008) predicted by PREP-Cp; Table S7: RNA editing sites in chloroplast genes of T. urartu (PI428335) predicted by PREP-Cp; Table S8: RNA editing sites in chloroplast genes of Ae. speltoides (AE918, TA1796 and PI487232) predicted by PREP-Cp; Table S9: RNA editing sites in chloroplast genes of T. turgidum ssp. dicoccoides (TA0073, TA0060 and TA1133) predicted by PREP-Cp.

Acknowledgments

We thank Le Wang for good advice and technical support. This work was supported by the National Natural Science Foundation of China (Grant No. 31401373 and 31561143005), and partially supported by the Open Project of the State Key Laboratory of Crop Stress Biology in Arid Areas, Northwest Agriculture and Forestry University (Grant No. CSBAA2014002).

Author Contributions

Xiaojun Nie conceived and designed the experiments. Hui Liu and Lingqiao Ge performed the experiments. Mengxing Wang, Guangwei Xing and Meng Wang analyzed the data. Mengxing Wang and Xiaojun Nie wrote the paper. Song Weining contributed reagents/materials/analysis tools. Xiaojun Nie and Song Weining supervised the entire study.

Conflicts of Interest

The authors declare no conflicts of interest. The funding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. Wang, C.; Xu, J.-R.; Liu, H. A-to-I RNA editing independent of adars in filamentous fungi. RNA Biol. 2016, 13, 940–945. [Google Scholar] [CrossRef] [PubMed]
  2. Wang, W.; Zhang, W.; Wu, Y.; Maliga, P.; Messing, J. RNA editing in chloroplasts of Spirodela polyrhiza, an aquatic monocotelydonous species. PLoS ONE 2015, 10, e0140285. [Google Scholar] [CrossRef] [PubMed]
  3. Brennicke, A.; Marchfelder, A.; Binder, S. RNA editing. FEMS Microbiol. Rev. 1999, 23, 297–316. [Google Scholar] [CrossRef] [PubMed]
  4. Hanson, M.R.; Sutton, C.; Luis, B. Plant organelle gene expression: Altered by RNA editing. Trends Plant Sci. 1996, 1, 57–64. [Google Scholar] [CrossRef]
  5. Liu, H.; Wang, Q.; He, Y.; Chen, L.; Hao, C.; Jiang, C.; Li, Y.; Dai, Y.; Kang, Z.; Xu, J.-R. Genome-wide A-to-I RNA editing in fungi independent of ADAR enzymes. Genome Res. 2016, 26, 499–509. [Google Scholar] [CrossRef] [PubMed]
  6. Hoch, B.; Maier, R.M.; Appel, K.; Igloi, G.L.; Kossel, H. Editing of a chloroplast mRNA by creation of an initiation codon. Nature 1991, 353, 178–180. [Google Scholar] [CrossRef] [PubMed]
  7. Hirose, T.; Kusumegi, T.; Tsudzuki, T.; Sugiura, M. RNA editing sites in tobacco chloroplast transcripts: Editing as a possible regulator of chloroplast RNA polymerase activity. Mol. Gen. Genet. MGG 1999, 262, 462–467. [Google Scholar] [CrossRef] [PubMed]
  8. Sasaki, T.; Yukawa, Y.; Miyamoto, T.; Obokata, J.; Sugiura, M. Identification of RNA editing sites in chloroplast transcripts from the maternal and paternal progenitors of tobacco (Nicotiana tabacum): Comparative analysis shows the involvement of distinct trans-factors for ndhB editing. Mol. Biol. Evol. 2003, 20, 1028–1035. [Google Scholar] [CrossRef] [PubMed]
  9. Lutz, K.A.; Maliga, P. Lack of conservation of editing sites in mRNAs that encode subunits of the NAD(P)H dehydrogenase complex in plastids and mitochondria of Arabidopsis thaliana. Curr. Genet. 2001, 40, 214–219. [Google Scholar] [CrossRef] [PubMed]
  10. Inada, M.; Sasaki, T.; Yukawa, M.; Tsudzuki, T.; Sugiura, M. A systematic search for RNA editing sites in pea chloroplasts: An editing event causes diversification from the evolutionarily conserved amino acid sequence. Plant Cell Physiol. 2004, 45, 1615–1622. [Google Scholar] [CrossRef] [PubMed]
  11. Kahlau, S.; Aspinall, S.; Gray, J.C.; Bock, R. Sequence of the tomato chloroplast DNA and evolutionary comparison of solanaceous plastid genomes. J. Mol. Evol. 2006, 63, 194–207. [Google Scholar] [CrossRef] [PubMed]
  12. Sasaki, T.; Yukawa, Y.; Wakasugi, T.; Yamada, K.; Sugiura, M. A simple in vitro RNA editing assay for chloroplast transcripts using fluorescent dideoxynucleotides: Distinct types of sequence elements required for editing of NDH transcripts. Plant J. 2006, 47, 802–810. [Google Scholar] [CrossRef] [PubMed]
  13. Duff, R.J.; Moore, F.B.-G. Pervasive RNA editing among hornwort rbcL transcripts except leiosporoceros. J. Mol. Evol. 2005, 61, 571–578. [Google Scholar] [CrossRef] [PubMed]
  14. Wolf, P.G.; Rowe, C.A.; Hasebe, M. High levels of RNA editing in a vascular plant chloroplast genome: Analysis of transcripts from the fern Adiantum capillus-veneris. Gene 2004, 339, 89–97. [Google Scholar] [CrossRef] [PubMed]
  15. Stern, D.B.; Goldschmidt-Clermont, M.; Hanson, M.R. Chloroplast RNA metabolism. Annu. Rev. Plant Biol. 2010, 61, 125–155. [Google Scholar] [CrossRef] [PubMed]
  16. Freyer, R.; Kiefer-Meyer, M.-C.; Kössel, H. Occurrence of plastid RNA editing in all major lineages of land plants. Proc. Natl. Acad. Sci. USA 1997, 94, 6285–6290. [Google Scholar] [CrossRef] [PubMed]
  17. Maier, R.M.; Neckermann, K.; Igloi, G.L.; Kössel, H. Complete sequence of the maize chloroplast genome: Gene content, hotspots of divergence and fine tuning of genetic information by transcript editing. J. Mol. Biol. 1995, 251, 614–628. [Google Scholar] [CrossRef] [PubMed]
  18. Tsudzuki, T.; Wakasugi, T.; Sugiura, M. Comparative analysis of RNA editing sites in higher plant chloroplasts. J. Mol. Evol. 2001, 53, 327–332. [Google Scholar] [CrossRef] [PubMed]
  19. Shikanai, T. RNA editing in plant organelles: Machinery, physiological function and evolution. Cell. Mol. Life Sci. CMLS 2006, 63, 698–708. [Google Scholar] [CrossRef] [PubMed]
  20. Smith, H.C.; Gott, J.M.; Hanson, M.R. A guide to RNA editing. RNA 1997, 3, 1105–1123. [Google Scholar] [PubMed]
  21. Barkan, A.; Small, I. Pentatricopeptide repeat proteins in plants. Annu. Rev. Plant Biol. 2014, 65, 415–442. [Google Scholar] [CrossRef] [PubMed]
  22. Kotera, E.; Tasaka, M.; Shikanai, T. A pentatricopeptide repeat protein is essential for RNA editing in chloroplasts. Nature 2005, 433, 326–330. [Google Scholar] [CrossRef] [PubMed]
  23. Takenaka, M.; Zehrmann, A.; Verbitskiy, D.; Härtel, B.; Brennicke, A. RNA editing in plants and its evolution. Annu. Rev. Genet. 2013, 47, 335–352. [Google Scholar] [CrossRef] [PubMed]
  24. Zehrmann, A.; Verbitskiy, D.; Härtel, B.; Brennicke, A.; Takenaka, M. RNA editing competence of trans-factor MEF1 is modulated by ecotype-specific differences but requires the DYW domain. FEBS Lett. 2010, 584, 4181–4186. [Google Scholar] [CrossRef] [PubMed]
  25. Dubcovsky, J.; Dvorak, J. Genome plasticity a key factor in the success of polyploid wheat under domestication. Science 2007, 316, 1862–1866. [Google Scholar] [CrossRef] [PubMed]
  26. Salamini, F.; Ozkan, H.; Brandolini, A.; Schafer-Pregl, R.; Martin, W. Genetics and geography of wild cereal domestication in the near east. Nat. Rev. Genet. 2002, 3, 429–441. [Google Scholar] [PubMed]
  27. Jia, J.; Zhao, S.; Kong, X.; Li, Y.; Zhao, G.; He, W.; Appels, R.; Pfeifer, M.; Tao, Y.; Zhang, X.; et al. Aegilops tauschii draft genome sequence reveals a gene repertoire for wheat adaptation. Nature 2013, 496, 91–95. [Google Scholar] [CrossRef] [Green Version]
  28. Luo, M.-C.; Gu, Y.Q.; You, F.M.; Deal, K.R.; Ma, Y.; Hu, Y.; Huo, N.; Wang, Y.; Wang, J.; Chen, S.; et al. A 4-gigabase physical map unlocks the structure and evolution of the complex genome of Aegilops tauschii, the wheat D-genome progenitor. Proc. Natl. Acad. Sci. USA 2013, 110, 7940–7945. [Google Scholar]
  29. Gornicki, P.; Zhu, H.; Wang, J.; Challa, G.S.; Zhang, Z.; Gill, B.S.; Li, W. The chloroplast view of the evolution of polyploid wheat. New Phytol. 2014, 204, 704–714. [Google Scholar]
  30. Wang, J.; Luo, M.-C.; Chen, Z.; You, F.M.; Wei, Y.; Zheng, Y.; Dvorak, J. Aegilops tauschii single nucleotide polymorphisms shed light on the origins of wheat D-genome genetic diversity and pinpoint the geographic origin of hexaploid wheat. New Phytol. 2013, 198, 925–937. [Google Scholar] [CrossRef] [PubMed]
  31. Mower, J.P. The PREP suite: Predictive rna editors for plant mitochondrial genes, chloroplast genes and user-defined alignments. Nucleic Acids Res. 2009, 37, W253–W259. [Google Scholar] [CrossRef] [PubMed]
  32. Hall, T.A. Bioedit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 1999, 41, 95–98. [Google Scholar]
  33. Trapnell, C.; Pachter, L.; Salzberg, S.L. Tophat: Discovering splice junctions with RNA-seq. Bioinformatics 2009, 25, 1105–1111. [Google Scholar] [CrossRef] [PubMed]
  34. Thorvaldsdóttir, H.; Robinson, J.T.; Mesirov, J.P. Integrative Genomics Viewer (IGV): High-performance genomics data visualization and exploration. Brief. Bioinform. 2013, 14, 178–192. [Google Scholar] [CrossRef] [PubMed]
  35. Chou, P.Y.; Fasman, G.D. Conformational parameters for amino acids in helical, β-sheet, and random coil regions calculated from proteins. Biochemistry 1974, 13, 211–222. [Google Scholar] [CrossRef] [PubMed]
  36. Chou, P.Y.; Fasman, G.D. Prediction of protein conformation. Biochemistry 1974, 13, 222–245. [Google Scholar] [CrossRef] [PubMed]
  37. Käll, L.; Krogh, A.; Sonnhammer, E.L.L. A combined transmembrane topology and signal peptide prediction method. J. Mol. Biol. 2004, 338, 1027–1036. [Google Scholar] [CrossRef] [PubMed]
  38. Tsunewaki, K.; Matsuoka, Y.; Yamazaki, Y.; Ogihara, Y. Evolutionary dynamics of wheat mitochondrial gene structure with special remarks on the origin and effects of rna editing in cereals. Genes Genet. Syst. 2008, 83, 301–320. [Google Scholar] [CrossRef] [PubMed]
  39. Middleton, C.P.; Senerchia, N.; Stein, N.; Akhunov, E.D.; Keller, B.; Wicker, T.; Kilian, B. Sequencing of chloroplast genomes from wheat, barley, rye and their relatives provides a detailed insight into the evolution of the Triticeae tribe. PLoS ONE 2014, 9, e85761. [Google Scholar] [CrossRef] [PubMed]
  40. Diekmann, K.; Hodkinson, T.R.; Wolfe, K.H.; van den Bekerom, R.; Dix, P.J.; Barth, S. Complete chloroplast genome sequence of a major allogamous forage species, perennial ryegrass (Lolium perenne L.). DNA Res. 2009, 16, 165–176. [Google Scholar] [CrossRef] [PubMed]
  41. Corneille, S.; Lutz, K.; Maliga, P. Conservation of rna editing between rice and maize plastids: Are most editing events dispensable? Mol. Gen. Genet. MGG 2000, 264, 419–424. [Google Scholar] [CrossRef] [PubMed]
  42. Calsa Júnior, T.; Carraro, D.M.; Benatti, M.R.; Barbosa, A.C.; Kitajima, J.P.; Carrer, H. Structural features and transcript-editing analysis of sugarcane (Saccharum officinarum L.) chloroplast genome. Curr. Genet. 2004, 46, 366–373. [Google Scholar] [CrossRef] [PubMed]
  43. Song, W.; Henry, R.J. Molecular analysis of the DNA polymorphism of wild barley (Hordeum spontaneum) germplasm using the polymerase chain reaction. Genet. Resour. Crop Evol. 1995, 42, 273–280. [Google Scholar] [CrossRef]
  44. Larkin, M.A.; Blackshields, G.; Brown, N.P.; Chenna, R.; McGettigan, P.A.; McWilliam, H.; Valentin, F.; Wallace, I.M.; Wilm, A.; Lopez, R.; et al. Clustal W and clustal X version 2.0. Bioinformatics 2007, 23, 2947–2948. [Google Scholar] [CrossRef] [PubMed]
  45. Ruwe, H.; Castandet, B.; Schmitz-Linneweber, C.; Stern, D.B. Arabidopsis chloroplast quantitative editotype. FEBS Lett. 2013, 587, 1429–1433. [Google Scholar] [CrossRef] [PubMed]
  46. Morton, R.B. The role of context-dependent mutations in generating compositional and codon usage bias in grass chloroplast DNA. J. Mol. Evol. 2003, 56, 616–629. [Google Scholar] [CrossRef] [PubMed]
  47. Morton, R.B.; Oberholzer, M.V.; Clegg, T.M. The influence of specific neighboring bases on substitution bias in noncoding regions of the plant chloroplast genome. J. Mol. Evol. 1997, 45, 227–231. [Google Scholar] [CrossRef] [PubMed]
  48. Tillich, M.; Lehwark, P.; Morton, B.R.; Maier, U.G. The evolution of chloroplast RNA editing. Mol. Biol. Evol. 2006, 23, 1912–1921. [Google Scholar] [CrossRef] [PubMed]
  49. Grosche, C.; Funk, H.T.; Maier, U.G.; Zauner, S. The chloroplast genome of Pellia endiviifolia: Gene content, RNA-editing pattern, and the origin of chloroplast editing. Genome Biol. Evol. 2012, 4, 1349–1357. [Google Scholar] [CrossRef] [PubMed]
  50. Kugita, M.; Yamamoto, Y.; Fujikawa, T.; Matsumoto, T.; Yoshinaga, K. RNA editing in hornwort chloroplasts makes more than half the genes functional. Nucleic Acids Res. 2003, 31, 2417–2423. [Google Scholar] [CrossRef] [PubMed]
  51. Jiang, Y.; Fan, S.L.; Song, M.Z.; Yu, J.N.; Yu, S.X. Identification of rna editing sites in cotton (Gossypium hirsutum) chloroplasts and editing events that affect secondary and three-dimensional protein structures. Genet. Mol. Res. GMR 2012, 11, 987–1001. [Google Scholar] [CrossRef] [PubMed]
  52. Du, P.; Li, Y. Prediction of C-to-U RNA editing sites in plant mitochondria using both biochemical and evolutionary information. J. Theor. Biol. 2008, 253, 579–586. [Google Scholar] [CrossRef] [PubMed]
  53. Chen, H.; Deng, L.; Jiang, Y.; Lu, P.; Yu, J. RNA editing sites exist in protein-coding genes in the chloroplast genome of Cycas taitungensis. J. Integr. Plant Biol. 2011, 53, 961–970. [Google Scholar] [CrossRef] [PubMed]
  54. Huang, Y.-Y.; Matzke, A.J.M.; Matzke, M. Complete sequence and comparative analysis of the chloroplast genome of coconut palm (Cocos nucifera). PLoS ONE 2013, 8, e74736. [Google Scholar] [CrossRef] [PubMed]
  55. Timme, R.E.; Kuehl, J.V.; Boore, J.L.; Jansen, R.K. A comparative analysis of the Lactuca and Helianthus (Asteraceae) plastid genomes: Identification of divergent regions and categorization of shared repeats. Am. J. Bot. 2007, 94, 302–312. [Google Scholar] [CrossRef] [PubMed]
  56. Zito, F.; Kuras, R.; Choquet, Y.; Kössel, H.; Wollman, F.-A. Mutations of cytochrome b6 in Chlamydomonas reinhardtii disclose the functional significance for a proline to leucine conversion by petB editing in maize and tobacco. Plant Mol. Biol. 1997, 33, 79–86. [Google Scholar] [CrossRef] [PubMed]
  57. Maier, R.M.; Hoch, B.; Zeltz, P.; Kössel, H. Internal editing of the maize chloroplast ndha transcript restores codons for conserved amino acids. Plant Cell 1992, 4, 609–616. [Google Scholar] [CrossRef] [PubMed]
  58. Karcher, D.; Bock, R. The amino acid sequence of a plastid protein is developmentally regulated by RNA editing. J. Biol. Chem. 2002, 277, 5570–5574. [Google Scholar] [CrossRef] [PubMed]
  59. Maier, U.G.; Bozarth, A.; Funk, H.T.; Zauner, S.; Rensing, S.A.; Schmitz-Linneweber, C.; Börner, T.; Tillich, M. Complex chloroplast RNA metabolism: Just debugging the genetic programme? BMC Biol. 2008, 6, 1–9. [Google Scholar] [CrossRef] [PubMed]
Genes 08 00013 g001 550
Figure 1. Nucleotide context of editing sites in Ae. tauschii (AL8/78). pos −1 and pos +1 indicate the direct upstream and downstream nucleotide of the edited sites (Ed).
Figure 1. Nucleotide context of editing sites in Ae. tauschii (AL8/78). pos −1 and pos +1 indicate the direct upstream and downstream nucleotide of the edited sites (Ed).
Genes 08 00013 g001
Genes 08 00013 g002 550
Figure 2. Types and quantities of amino acid exchanges arising from chloroplast RNA editing in Ae. tauschii (AL8/78). The red bars represent a conversion of amino acid from hydrophilic to hydrophobic, while blue is for the opposite direction and green is for no hydrophobicity exchanges.
Figure 2. Types and quantities of amino acid exchanges arising from chloroplast RNA editing in Ae. tauschii (AL8/78). The red bars represent a conversion of amino acid from hydrophilic to hydrophobic, while blue is for the opposite direction and green is for no hydrophobicity exchanges.
Genes 08 00013 g002
Genes 08 00013 g003 550
Figure 3. Transmembrane topology of petB, ndhA and ndhB before and after editing. (A) The transmembrane region of petB expanded to the C terminal of the protein after editing; (B) The transmembrane region of ndhA displayed a contraction and shift after editing; (C) NdhB showed a larger transmembrane region after editing.
Figure 3. Transmembrane topology of petB, ndhA and ndhB before and after editing. (A) The transmembrane region of petB expanded to the C terminal of the protein after editing; (B) The transmembrane region of ndhA displayed a contraction and shift after editing; (C) NdhB showed a larger transmembrane region after editing.
Genes 08 00013 g003
Genes 08 00013 g004 550
Figure 4. Comparison of RNA editing sites in chloroplast genes of T. aestivum and its genome donors. (A) Comparison among T. urartu, Ae. speltoides and T. turgidum; (B) Comparison among T. turgidum, Ae. tauschii and T. aestivum; (C) Comparison among T. urartu, Ae. speltoides, Ae. tauschii and T. aestivum. Tu indicates T. urartu (PI428335); As indicates Ae. speltoides (AE918, TA1796 and PI487232); At indicates Ae. tauschii (AL8/78); Tt indicates T. turgidum ssp. dicoccoides (TA0073, TA0060 and TA1133); and Ta indicates T. aestivum cv CS (TA3008). T. urartu, Triticum urartu L. Ae. speltoides, Aegilops speltoides L. T. turgidum, Triticum turgidum L. Ae. tauschii, Aegilops tauschii L. T. aestivum, Triticum aestivum L.
Figure 4. Comparison of RNA editing sites in chloroplast genes of T. aestivum and its genome donors. (A) Comparison among T. urartu, Ae. speltoides and T. turgidum; (B) Comparison among T. turgidum, Ae. tauschii and T. aestivum; (C) Comparison among T. urartu, Ae. speltoides, Ae. tauschii and T. aestivum. Tu indicates T. urartu (PI428335); As indicates Ae. speltoides (AE918, TA1796 and PI487232); At indicates Ae. tauschii (AL8/78); Tt indicates T. turgidum ssp. dicoccoides (TA0073, TA0060 and TA1133); and Ta indicates T. aestivum cv CS (TA3008). T. urartu, Triticum urartu L. Ae. speltoides, Aegilops speltoides L. T. turgidum, Triticum turgidum L. Ae. tauschii, Aegilops tauschii L. T. aestivum, Triticum aestivum L.
Genes 08 00013 g004
Table
Table 1. Information of complete chloroplast genome for comparative analysis of RNA editing sites in protein-coding genes among Poaceae species.
Table 1. Information of complete chloroplast genome for comparative analysis of RNA editing sites in protein-coding genes among Poaceae species.
Material NameGenBank Accession
Ae. tauschii (AL8/78)KJ614412.1
T. urartu (PI428335)KJ614411.1
Ae. speltoides ssp. ligustica (AE918)KJ614404.1
Ae. speltoides ssp. ligustica (TA1796)KJ614405.1
Ae. speltoides ssp. speltoides (PI487232)KJ614406.1
T. turgidum ssp. dicoccoides (TA0073)KJ614400.1
T. turgidum ssp. dicoccoides (TA0060)KJ614401.1
T. turgidum ssp. dicoccoides (TA1133)KJ614402.1
T. aestivum ssp. aestivum cv Chinese Spring (CS) (TA3008)KJ614396.1
Table
Table 2. RNA editing sites in chloroplast protein-coding genes of Ae. tauschii (AL8/78) predicted by PREP-Cp.
Table 2. RNA editing sites in chloroplast protein-coding genes of Ae. tauschii (AL8/78) predicted by PREP-Cp.
GeneNucleotide PositionAmino Acid PositionCodon ConversionAmino Acid Conversion
ndhA473 a158tCa>tTaS>L
563 a188tCa>tTaS>L
1070357tCt>tTtS>F
atpA1148 a383tCa>tTaS>L
atpB1487496tCg>tTgS>L
matK1261 a421Cat>TatH>Y
ndhB149a50tCa>tTaS>L
467 a156cCa>cTaP>L
586 a196Cat>TatH>Y
611 a204tCa>tTaS>L
704 a235tCc>tTcS>F
737 a246cCa>cTaP>L
830 a277tCa>tTaS>L
836 a279tCa>tTaS>L
1481 a494cCa>cTaP>L
ndhD878 a293tCa>tTaS>L
ndhF6221tCa>tTaS>L
1487496aCg>aTgT>M
petB662221cCa>cTaP>L
rpoC22009670cCa>cTaP>L
2030677cCa>cTaP>L
2158720Ccc>TccP>S
30021001cCg>cTgP>L
40311344tCg>tTgS>L
rpl26221aCt>aTtT>I
rpl20308103tCa>tTaS>L
rpoA1009337Ctc>TtcL>F
rpoB467 a156tCg>tTgS>L
545 a182tCa>tTaS>L
560 a187tCg>tTgS>L
617 b206cCg>cTgP>L
rps818261tCa>tTaS>L
ycf34415tCc>tTcS>F
19164aCg>aTgT>M
Capitals in column Codon Conversion indicate target nucleotides. a Editing sites with experimental verification. b Unedited by experimental verification.
Table
Table 3. Cytidine (C) to uridine (U) editing sites in chloroplast genome of Ae. tauschii (AL8/78) validated by RNA-Seq data.
Table 3. Cytidine (C) to uridine (U) editing sites in chloroplast genome of Ae. tauschii (AL8/78) validated by RNA-Seq data.
Genome PositionGene PositionCodonAmino AcidLeafPistilRootSeedSheathSpikeStamenStem
1 b23451123412312341231123
Sense-strand
3090~psbK 22
3551~psbK 7
8259psbI~psbD 35 c36373637 4744 33293233 2522
8481psbI~psbD 9
11,311psbC~psbZ 21 1920 18 10 10 9
20,093rpoB-467 atCgS>L 3828 66757890 72 32
20,171rpoB-545 atCaS>L 48 70868878 71
20,186rpoB-560 atCgS>L 57 73868085 69
29,181rpoC2-4031 atCgS>L94++9198+ +827160739799939792 +9299 98
32,404atpH~atpF 1818 15
34,440atpA-234ggCG>G 11
35,354atpA-1148 atCaS>L9998999999 839498919491879690949695 9471939594
42,415atpA~trnS-GGA 14
57,960ycf4~cemA 525067485512 37422731 23 26
62,758petL-56 dcCaP>L8375627760 58 3630 83
65,352rps18-448CgaR>Stop 9
68,447psbB-414atCL>L 17
71,835petB-662 acCaP>L9698979890 43 27 34 665826 13
127,814ndhB-149 atCaS>L94 96
128,251ndhB-586 aCatH>Y9588919393
128,276ndhB-611 atCaS>L95 939594
Genome PositionGene PositionCodonAmino AcidLeafPistilRootSeedSheathSpikeStamenStem
123451123412312341231123
Sense-strand
128,369ndhB-704 atCcS>F 92
129,858ndhB-1481 acCaP>L969999 97
Antisense-strand
1956matK-1261 aCatH>Y664541695213 36 6955386955 503661 71
5100rps16-intron 191722232289545150 65 3633323783 8360667060
5179rps16-intron 151523
36,673rps14~psaB 4545454445 20 2014 1923 8
42,982ycf3-191 aaCgT>M79838281806344554539 497069676364 746472 76
43,880ycf3-44 atCcS>F9292969792407168702646 7070698242 5043697774
44,727ycf3~rps4 786871756538 41 35293519 2021
44,868ycf3~rps4 11
45,244rps4-305tCaS>L 20
49,393ndhK-125 dcCaP>L9998999999+ 70+71 96979494+ 938995 89
49,858ndhC-13CacH>Y9898989999758683636386 86928591 568083 63
61,027psbL-111ttCF>F778383814976 47453329 56
65,383trnP-GGG~rpl20 24
65,632rpl20-308 atCaS>L 11 24 60386375 56
67,968clpP~psbN 31
73,402psbN~rpoA 4550515145 35 1912132120 1514 22
73,980rpoA-527 dtCcS>F3034302828 231831357964728184819084758541878286
74,553rpoA~rps11 414444454121 1911291512222120221818 1817232020
76,074rps8-182 atCaS>L899089909096838683849696959797979793959396959595
77,285rpl16-36ccCP>P 18
78,313rpl16~rps3 92
Genome PositionGene PositionCodonAmino AcidLeafPistilRootSeedSheathSpikeStamenStem
123451123412312341231123
Antisense-strand
79,174rps3-30 dttCF>F625762665977586671697882636668697168636362565748
81,816rpl2~rpl23 82 878888
86,095ndhB-836 atCaS>L 99
87,176ndhB-467 acCaP>L9995969798 96 89
87,494ndhB-149 atCaS>L 95
96,705rps7~trnN-GUU 6 12
101,658ndhF-1891CagQ>Stop 9
103,487ndhF-62 atCaS>L4543444951 28 39394318 211744 37
106,742ndhD-1398 datCL>L3628373437 23
107,262ndhD-878 a,dtCaS>L9998999998 + 8995 8981897784 779489+83
109,393ndhE~ndhG 32
109,684ndhG-347 dcCaP>L9897979896 65596472 4078 65
110,034ndhG~ndhI 3827222819
110,946ndhA-1070 atCtS>F75 6492
111,453ndhA-563 atCaS>L13
112,572ndhA-473 atCaS>L + 4533
No. of site3529333335279101611167828282829133192520921
Capitals in column Codon indicate target nucleotides. a Editing sites consistent with the prediction. b Biological replicates of each tissue. c This value indicates editing efficiency (%). d Editing sites validated by direct sequencing of PCR products. +, Complete editing (editing efficiency, 100%). Blank, no editing/not yet determined.
Table
Table 4. Comparative analysis of RNA editing sites in chloroplast genes among six Poaceae species.
Table 4. Comparative analysis of RNA editing sites in chloroplast genes among six Poaceae species.
GeneSiteCodon PositionEdited CodonAmino Acid ConversionAe. tauschiiH. vulgareL. perenneO. sativaS. officinarumZ. mays
matK1420 (421)CatH→Y+++ a+
rpoB1156tCaS→L+++ a+ a++
2182tCaS→L+++ a+ a++
3187tCgS→L+++ a+ a uCa++
4206cCgP→L+−−−−+
rps14127tCaS→L(−) (−)++ a+
atpA1383tCaS→L+ ++++
ycf3115tCcS→F+ +(−)−+
262 (64)aCgT→M+ ++ a++ a
ndhG1116cCaP→L+ a ++ (−)
rpl201103tCaS→L+ + a(−)+a+
psbL137ttCF→F+ a + a(−)
rps8161tCaS→L+ ++++
rpl211aCgT→M(−)++ a+ a+a+
ndhB150tCaS→L+++(−)(−)(−)
2156cCaP→L++++++
3196CatH→Y++++++
4204tCaS→L++++ a++
5235tCcS→F++++(−)(−)
6246cCaP→L++++++
7277tCaS→L++++++
8279tCaS→L++++(−)(−)
9494cCaP→L++++ a++
ndhF121tCaS→L+ +++ a+
ndhD1295 (293)tCaS→L+++ a+++
ndhA117tCaS→L(−)++(−)++
2158tCaS→L++++++
3188tCaS→L++++++
4357tCcS→F++(−)+++
ndhK12gtCV→V− +
243(42)cCaP→L+ +
Capitals in column Edited Codon indicate target nucleotides. “+”, complete editing; “−”, no editing and C encoded in DNA; “(−)”, no editing and T encoded in DNA; “+ a”, partial editing; Blank, no information available/not yet determined. Ae. tauschii, Aegilops tauschii L. H. vulgare, Hordeum vulgare L. L. perenne, Lolium perenne L. O. sativa, Oryza sativa L. S. officinarum, Saccharum officinarum L. Z. mays, Zea mays L.

Share and Cite

MDPI and ACS Style

Wang, M.; Liu, H.; Ge, L.; Xing, G.; Wang, M.; Weining, S.; Nie, X. Identification and Analysis of RNA Editing Sites in the Chloroplast Transcripts of Aegilops tauschii L. Genes 2017, 8, 13. https://0-doi-org.brum.beds.ac.uk/10.3390/genes8010013

AMA Style

Wang M, Liu H, Ge L, Xing G, Wang M, Weining S, Nie X. Identification and Analysis of RNA Editing Sites in the Chloroplast Transcripts of Aegilops tauschii L. Genes. 2017; 8(1):13. https://0-doi-org.brum.beds.ac.uk/10.3390/genes8010013

Chicago/Turabian Style

Wang, Mengxing, Hui Liu, Lingqiao Ge, Guangwei Xing, Meng Wang, Song Weining, and Xiaojun Nie. 2017. "Identification and Analysis of RNA Editing Sites in the Chloroplast Transcripts of Aegilops tauschii L." Genes 8, no. 1: 13. https://0-doi-org.brum.beds.ac.uk/10.3390/genes8010013

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop