Editing site analysis in a gymnosperm mitochondrial genome reveals similarities with angiosperm mitochondrial genomes

Assembly and characterization of the complete mitochondrial genome of Ventilago leiocarpa

KEY MESSAGE

We reported the mitochondrial genome of Ventilago leiocarpa for the first time. Two and one sites lead to the generation of stop and stat codon through editing were verified. Ventilago leiocarpa, a member of the Rhamnaceae family, is frequently utilized in traditional medicine due to the medicinal properties of its roots. In this study, we successfully assembled the mitogenome of V. leiocarpa using both BGI short reads and Nanopore long reads. This mitogenome has a total length of 331,839 bp. The annotated results showed 36 unique protein-coding, 16 tRNA and 3 rRNA genes in this mitogenome. Furthermore, we confirmed the presence of a branched structure through the utilization of long reads mapping, PCR amplification, and Sanger sequencing. Specifically, the ctg1 can form a single circular molecule or combine with ctg4 to form a linear molecule. Likewise, ctg2 can form a single circular molecule or can be connected to ctg4 to form a linear molecule. Subsequently, through a comparative analysis of the mitogenome and cpgenome sequences, we identified ten mitochondrial plastid sequences (MTPTs), including two complete protein-coding genes and five complete tRNA genes. The existence of MTPTs was verified by long reads. Colinear analysis showed that the mitogenomes of Rosales were highly divergent in structure. Finally, we identified 545 RNA editing sites involving 36 protein-coding genes by Deepred-mt. To validate our findings, we conducted PCR amplification and Sanger sequencing, which confirmed the generation of stop codons in atp9-223 and rps10-391, as well as the generation of a start codon in nad4L-2. This project reported the complex structure and RNA editing event of the V. Leiocarpa mitogenome, which will provide valuable information for the study of mitochondrial gene expression.

Complete sequence and comparative analysis of the mitochondrial genome of the rare and endangered Clematis acerifolia, the first clematis mitogenome to provide new insights into the phylogenetic evolutionary status of the genus

Clematis is one of the large worldwide genera of the Ranunculaceae Juss. Family, with high ornamental and medicinal value. China is the modern distribution centre of Clematis with abundant natural populations. Due to the complexity and high morphological diversity of Clematis, the genus is difficult to classify systematically, and in particular, the phylogenetic position of the endangered Clematis acerifolia is highly controversial. The use of the mitochondrial complete genome is a powerful molecular method that is frequently used for inferring plants phylogenies. However, studies on Clematis mitogenome are rare, thus limiting our full understanding of its phylogeny and genome evolution. Here, we sequenced and annotated the C. acerifolia mt genome using Illumina short- and Nanopore long-reads, characterized the species first complete mitogenome, and performed a comparative phylogenetic analysis with its close relatives. The total length of the C. acerifolia mitogenome is 698,247 bp and the main structure is multi-branched (linear molecule 1 and circular molecule 2). We annotated 55 genes, including 35 protein-coding, 17 tRNA, and 3 rRNA genes. The C. acerifolia mitogenome has extremely unconserved structurally, with extensive sequence transfer between the chloroplast and mitochondrial organelles, sequence repeats, and RNA editing. The phylogenetic position of C. acerifolia was determined by constructing the species mitogenome with 24 angiosperms. Further, our C. acerifolia mitogenome characteristics investigation included GC contents, codon usage, repeats and synteny analysis. Overall, our results are expected to provide fundamental information for C. acerifolia mitogenome evolution and confirm the validity of mitochondrial analysis in determining the phylogenetic positioning of Clematis plants.

The complete mitochondrial genome of Taxus cuspidata (Taxaceae): eight protein-coding genes have transferred to the nuclear genome

BACKGROUND

Gymnosperms represent five of the six lineages of seed plants. However, most sequenced plant mitochondrial genomes (mitogenomes) have been generated for angiosperms, whereas mitogenomic sequences have been generated for only six gymnosperms. In particular, complete mitogenomes are available for all major seed plant lineages except Conifer II (non-Pinaceae conifers or Cupressophyta), an important lineage including six families, which impedes a comprehensive understanding of the mitogenomic diversity and evolution in gymnosperms.

RESULTS

Here, we report the complete mitogenome of Taxus cuspidata in Conifer II. In comparison with previously released gymnosperm mitogenomes, we found that the mitogenomes of Taxus and Welwitschia have lost many genes individually, whereas all genes were identified in the mitogenomes of Cycas, Ginkgo and Pinaceae. Multiple tRNA genes and introns also have been lost in some lineages of gymnosperms, similar to the pattern observed in angiosperms. In general, gene clusters could be less conserved in gymnosperms than in angiosperms. Moreover, fewer RNA editing sites were identified in the Taxus and Welwitschia mitogenomes than in other mitogenomes, which could be correlated with fewer introns and frequent gene losses in these two species.

CONCLUSIONS

We have sequenced the Taxus cuspidata mitogenome, and compared it with mitogenomes from the other four gymnosperm lineages. The results revealed the diversity in size, structure, gene and intron contents, foreign sequences, and mutation rates of gymnosperm mitogenomes, which are different from angiosperm mitogenomes.

Comparative RNA editing profile of mitochondrial transcripts in cytoplasmic male sterile and fertile pigeonpea reveal significant changes at the protein level

RNA editing is a process which leads to post-transcriptional alteration of the nucleotide sequence of the corresponding mRNA molecule which may or may not lead to changes at the protein level. Apart from its role in providing variability at the transcript and protein levels, sometimes, such changes may lead to abnormal expression of the mitochondrial gene leading to a cytoplasmic male sterile phenotype. Here we report the editing status of 20 major mitochondrial transcripts in both male sterile (AKCMS11) and male fertile (AKPR303) pigeonpea genotypes. The validation of the predicted editing sites was done by mapping RNA-seq reads onto the amplified mitochondrial genes, and 165 and 159 editing sites were observed in bud tissues of the male sterile and fertile plant respectively. Among the resulting amino acid alterations, the most frequent one was the conversion of hydrophilic amino acids to hydrophobic. The alterations thus detected in our study indicates differential editing, but no major change in terms of the abnormal protein structure was detected. However, the above investigation provides an insight into the behaviour of pigeonpea mitochondrial genome in native and alloplasmic state and could hold clues in identification of editing factors and their role in adaptive evolution in pigeonpea.

Plant organelle RNA editing and its specificity factors: enhancements of analyses and new database features in PREPACT 3.0

Background

Gene expression in plant chloroplasts and mitochondria is affected by RNA editing. Numerous C-to-U conversions, accompanied by reverse U-to-C exchanges in some plant clades, alter the genetic information encoded in the organelle genomes. Predicting and analyzing RNA editing, which ranges from only few sites in some species to thousands in other taxa, is bioinformatically demanding.

Results

Here, we present major enhancements and extensions of PREPACT, a WWW-based service for analysing, predicting and cataloguing plant-type RNA editing. New features in PREPACT’s core include direct GenBank accession query input and options to restrict searches to candidate U-to-C editing or to sites where editing has been documented previously in the references. The reference database has been extended by 20 new organelle editomes. PREPACT 3.0 features new modules “EdiFacts” and “TargetScan”. EdiFacts integrates information on pentatricopeptide repeat (PPR) proteins characterized as site-specific RNA editing factors. PREPACT’s editome references connect into EdiFacts, linking editing events to specific co-factors where known. TargetScan allows position-weighted querying for sequence motifs in the organelle references, optionally restricted to coding regions or sequences around editing sites, or in queries uploaded by the user. TargetScan is mainly intended to evaluate and further refine the proposed PPR-RNA recognition code but may be handy for other tasks as well. We present an analysis for the immediate sequence environment of more than 15,000 documented editing sites finding strong and different bias in the editome data sets.

Conclusions

We exemplarily present the novel features of PREPACT 3.0 aimed to enhance the analyses of plant-type RNA editing, including its new modules EdiFacts integrating information on characterized editing factors and TargetScan aimed to analyse RNA editing site recognition specificities.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2244-9) contains supplementary material, which is available to authorized users.

Towards a comprehensive picture of C-to-U RNA editing sites in angiosperm mitochondria

Our understanding of the dynamic and evolution of RNA editing in angiosperms is in part limited by the few editing sites identified to date. This study identified 10,217 editing sites from 17 diverse angiosperms. Our analyses confirmed the universality of certain features of RNA editing, and offer new evidence behind the loss of editing sites in angiosperms. RNA editing is a post-transcriptional process that substitutes cytidines (C) for uridines (U) in organellar transcripts of angiosperms. These substitutions mostly take place in mitochondrial messenger RNAs at specific positions called editing sites. By means of publicly available RNA-seq data, this study identified 10,217 editing sites in mitochondrial protein-coding genes of 17 diverse angiosperms. Even though other types of mismatches were also identified, we did not find evidence of non-canonical editing processes. The results showed an uneven distribution of editing sites among species, genes, and codon positions. The analyses revealed that editing sites were conserved across angiosperms but there were some species-specific sites. Non-synonymous editing sites were particularly highly conserved (~ 80%) across the plant species and were efficiently edited (80% editing extent). In contrast, editing sites at third codon positions were poorly conserved (~ 30%) and only partially edited (~ 40% editing extent). We found that the loss of editing sites along angiosperm evolution is mainly occurring by replacing editing sites with thymidines, instead of a degradation of the editing recognition motif around editing sites. Consecutive and highly conserved editing sites had been replaced by thymidines as result of retroprocessing, by which edited transcripts are reverse transcribed to cDNA and then integrated into the genome by homologous recombination. This phenomenon was more pronounced in eudicots, and in the gene cox1. These results suggest that retroprocessing is a widespread driving force underlying the loss of editing sites in angiosperm mitochondria.

RNA Editing and Its Molecular Mechanism in Plant Organelles

RNA editing by cytidine (C) to uridine (U) conversions is widespread in plant mitochondria and chloroplasts. In some plant taxa, “reverse” U-to-C editing also occurs. However, to date, no instance of RNA editing has yet been reported in green algae and the complex thalloid liverworts. RNA editing may have evolved in early land plants 450 million years ago. However, in some plant species, including the liverwort, Marchantia polymorpha, editing may have been lost during evolution. Most RNA editing events can restore the evolutionarily conserved amino acid residues in mRNAs or create translation start and stop codons. Therefore, RNA editing is an essential process to maintain genetic information at the RNA level. Individual RNA editing sites are recognized by plant-specific pentatricopeptide repeat (PPR) proteins that are encoded in the nuclear genome. These PPR proteins are characterized by repeat elements that bind specifically to RNA sequences upstream of target editing sites. In flowering plants, non-PPR proteins also participate in multiple RNA editing events as auxiliary factors. C-to-U editing can be explained by cytidine deamination. The proteins discovered to date are important factors for RNA editing but a bona fide RNA editing enzyme has yet to be identified.

Comparative analysis of 11 Brassicales mitochondrial genomes and the mitochondrial transcriptome of Brassica oleracea

To elucidate the evolution of mitochondrial genomic diversity within a single order of angiosperms, we sequenced seven Brassicales genomes and the transcriptome of Brassica oleracea. In the common ancestor of Brassicaceae, several genes of known function were lost and the ccmFN gene was split into two independent genes, which also coincides with a trend of genome reduction towards the smallest sequenced angiosperm genomes of Brassica. For most ORFs of unknown function, the lack of conservation throughout Brassicales and the generally low expression and absence of RNA editing in B. oleracea argue against functionality. However, two chimeric ORFs were expressed and edited in B. oleracea, suggesting a potential role in cytoplasmic male sterility in certain nuclear backgrounds. These results demonstrate how frequent shifts in size, structure, and content of plant mitochondrial genomes can occur over short evolutionary time scales.

RNA editing in plant mitochondria—connecting RNA target sequences and acting proteins

RNA editing changes several hundred cytidines to uridines in the mRNAs of mitochondria in flowering plants. The target cytidines are identified by a subtype of PPR proteins characterized by tandem modules which each binds with a specific upstream nucleotide. Recent progress in correlating repeat structures with nucleotide identities allows to predict and identify target sites in mitochondrial RNAs. Additional proteins have been found to play a role in RNA editing; their precise function still needs to be elucidated. The enzymatic activity performing the C to U reaction may reside in the C-terminal DYW extensions of the PPR proteins; however, this still needs to be proven. Here we update recent progress in understanding RNA editing in flowering plant mitochondria.

The "fossilized" mitochondrial genome of Liriodendron tulipifera: ancestral gene content and order, ancestral editing sites, and extraordinarily low mutation rate

BACKGROUND

The mitochondrial genomes of flowering plants vary greatly in size, gene content, gene order, mutation rate and level of RNA editing. However, the narrow phylogenetic breadth of available genomic data has limited our ability to reconstruct these traits in the ancestral flowering plant and, therefore, to infer subsequent patterns of evolution across angiosperms.

RESULTS

We sequenced the mitochondrial genome of Liriodendron tulipifera, the first from outside the monocots or eudicots. This 553,721 bp mitochondrial genome has evolved remarkably slowly in virtually all respects, with an extraordinarily low genome-wide silent substitution rate, retention of genes frequently lost in other angiosperm lineages, and conservation of ancestral gene clusters. The mitochondrial protein genes in Liriodendron are the most heavily edited of any angiosperm characterized to date. Most of these sites are also edited in various other lineages, which allowed us to polarize losses of editing sites in other parts of the angiosperm phylogeny. Finally, we added comprehensive gene sequence data for two other magnoliids, Magnolia stellata and the more distantly related Calycanthus floridus, to measure rates of sequence evolution in Liriodendron with greater accuracy. The Magnolia genome has evolved at an even lower rate, revealing a roughly 5,000-fold range of synonymous-site divergence among angiosperms whose mitochondrial gene space has been comprehensively sequenced.

CONCLUSIONS

Using Liriodendron as a guide, we estimate that the ancestral flowering plant mitochondrial genome contained 41 protein genes, 14 tRNA genes of mitochondrial origin, as many as 7 tRNA genes of chloroplast origin, >700 sites of RNA editing, and some 14 colinear gene clusters. Many of these gene clusters, genes and RNA editing sites have been variously lost in different lineages over the course of the ensuing ∽200 million years of angiosperm evolution.

The "fossilized" mitochondrial genome of Liriodendron tulipifera: ancestral gene content and order, ancestral editing sites, and extraordinarily low mutation rate

BACKGROUND

The mitochondrial genomes of flowering plants vary greatly in size, gene content, gene order, mutation rate and level of RNA editing. However, the narrow phylogenetic breadth of available genomic data has limited our ability to reconstruct these traits in the ancestral flowering plant and, therefore, to infer subsequent patterns of evolution across angiosperms.

RESULTS

We sequenced the mitochondrial genome of Liriodendron tulipifera, the first from outside the monocots or eudicots. This 553,721 bp mitochondrial genome has evolved remarkably slowly in virtually all respects, with an extraordinarily low genome-wide silent substitution rate, retention of genes frequently lost in other angiosperm lineages, and conservation of ancestral gene clusters. The mitochondrial protein genes in Liriodendron are the most heavily edited of any angiosperm characterized to date. Most of these sites are also edited in various other lineages, which allowed us to polarize losses of editing sites in other parts of the angiosperm phylogeny. Finally, we added comprehensive gene sequence data for two other magnoliids, Magnolia stellata and the more distantly related Calycanthus floridus, to measure rates of sequence evolution in Liriodendron with greater accuracy. The Magnolia genome has evolved at an even lower rate, revealing a roughly 5,000-fold range of synonymous-site divergence among angiosperms whose mitochondrial gene space has been comprehensively sequenced.

CONCLUSIONS

Using Liriodendron as a guide, we estimate that the ancestral flowering plant mitochondrial genome contained 41 protein genes, 14 tRNA genes of mitochondrial origin, as many as 7 tRNA genes of chloroplast origin, >700 sites of RNA editing, and some 14 colinear gene clusters. Many of these gene clusters, genes and RNA editing sites have been variously lost in different lineages over the course of the ensuing ∽200 million years of angiosperm evolution.

PREPACT 2.0: Predicting C-to-U and U-to-C RNA Editing in Organelle Genome Sequences with Multiple References and Curated RNA Editing Annotation

RNA editing is vast in some genetic systems, with up to thousands of targeted C-to-U and U-to-C substitutions in mitochondria and chloroplasts of certain plants. Efficient prognoses of RNA editing in organelle genomes will help to reveal overlooked cases of editing. We present PREPACT 2.0 (http://www.prepact.de) with numerous enhancements of our previously developed Plant RNA Editing Prediction & Analysis Computer Tool. Reference organelle transcriptomes for editing prediction have been extended and reorganized to include 19 curated mitochondrial and 13 chloroplast genomes, now allowing to distinguish RNA editing sites from “pre-edited” sites. Queries may be run against multiple references and a new “commons” function identifies and highlights orthologous candidate editing sites congruently predicted by multiple references. Enhancements to the BLASTX mode in PREPACT 2.0 allow querying of complete novel organelle genomes within a few minutes, identifying protein genes and candidate RNA editing sites simultaneously without prior user analyses.

PPR2263, a DYW-Subgroup Pentatricopeptide repeat protein, is required for mitochondrial nad5 and cob transcript editing, mitochondrion biogenesis, and maize growth

RNA editing plays an important role in organelle gene expression in various organisms, including flowering plants, changing the nucleotide information at precise sites. Here, we present evidence that the maize (Zea mays) nuclear gene Pentatricopeptide repeat 2263 (PPR2263) encoding a DYW domain-containing PPR protein is required for RNA editing in the mitochondrial NADH dehydrogenase5 (nad5) and cytochrome b (cob) transcripts at the nad5-1550 and cob-908 sites, respectively. Its putative ortholog, MITOCHONDRIAL EDITING FACTOR29, fulfills the same role in Arabidopsis thaliana. Both the maize and the Arabidopsis proteins show preferential localization to mitochondria but are also detected in chloroplasts. In maize, the corresponding ppr2263 mutation causes growth defects in kernels and seedlings. Embryo and endosperm growth are reduced, leading to the production of small but viable kernels. Mutant plants have narrower and shorter leaves, exhibit a strong delay in flowering time, and generally do not reach sexual maturity. Whereas mutant chloroplasts do not have major defects, mutant mitochondria lack complex III and are characterized by a compromised ultrastructure, increased transcript levels, and the induction of alternative oxidase. The results suggest that mitochondrial RNA editing at the cob-908 site is necessary for mitochondrion biogenesis, cell division, and plant growth in maize.

RNA editing sites exist in protein-coding genes in the chloroplast genome of Cycas taitungensis

RNA editing is a post-transcriptional process that results in modifications of ribonucleotides at specific locations. In land plants editing can occur in both mitochondria and chloroplasts and most commonly involves C-to-U changes, especially in seed plants. Using prediction and experimental determination, we investigated RNA editing in 40 protein-coding genes from the chloroplast genome of Cycas taitungensis. A total of 85 editing sites were identified in 25 transcripts. Comparison analysis of the published editotypes of these 25 transcripts in eight species showed that RNA editing events gradually disappear during plant evolution. The editing in the first and third codon position disappeared quicker than that in the second codon position. ndh genes have the highest editing frequency while serine and proline codons were more frequently edited than the codons of other amino acids. These results imply that retained RNA editing sites have imbalanced distribution in genes and most of them may function by changing protein structure or interaction. Mitochondrion protein-coding genes have three times the editing sites compared with chloroplast genes of Cycas, most likely due to slower evolution speed.

The RNA editing pattern of cox2 mRNA is affected by point mutations in plant mitochondria

The mitochondrial transcriptome from land plants undergoes hundreds of specific C-to-U changes by RNA editing. These events are important since most of them occur in the coding region of mRNAs. One challenging question is to understand the mechanism of recognition of a selected C residue (editing sites) on the transcript. It has been reported that a short region surrounding the target C forms the cis-recognition elements, but individual residues on it do not play similar roles for the different editing sites. Here, we studied the role of the −1 and +1 nucleotide in wheat cox2 editing site recognition using an in organello approach. We found that four different recognition patterns can be distinguished: (a) +1 dependency, (b) −1 dependency, (c) +1/−1 dependency, and (d) no dependency on nearest neighbor residues. A striking observation was that whereas a 23 nt cis region is necessary for editing, some mutants affect the editing efficiency of unmodified distant sites. As a rule, mutations or pre-edited variants of the transcript have an impact on the complete set of editing targets. When some Cs were changed into Us, the remaining editing sites presented a higher efficiency of C-to-U conversion than in wild type mRNA. Our data suggest that the complex response observed for cox2 mRNA may be a consequence of the fate of the transcript during mitochondrial gene expression.

A unique transcriptome: 1782 positions of RNA editing alter 1406 codon identities in mitochondrial mRNAs of the lycophyte Isoetes engelmannii

The analysis of the mitochondrial DNA of Isoetes engelmannii as a first representative of the lycophytes recently revealed very small introns and indications for extremely frequent RNA editing. To analyze functionality of intron splicing and the extent of RNA editing in I. engelmannii, we performed a comprehensive analysis of its mitochondrial transcriptome. All 30 groups I and II introns were found to be correctly removed, showing that intron size reduction does not impede splicing. We find that mRNA editing affects 1782 sites, which lead to a total of 1406 changes in codon meanings. This includes the removal of stop codons from 23 of the 25 mitochondrial protein encoding genes. Comprehensive sequence analysis of multiple cDNAs per locus allowed classification of partially edited sites as either inefficiently edited but relevant or as non-specifically edited at mostly low frequencies. Abundant RNA editing was also found to affect tRNAs in hitherto unseen frequency, taking place at 41 positions in tRNA-precursors, including the first identification of U-to-C exchanges in two tRNA species. We finally investigated the four group II introns of the nad7 gene and could identify 27 sites of editing, most of which improve base pairing for proper secondary structure formation.