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

Assembly and analysis of the complete mitochondrial genome of an endemic Camellia species of China, Camellia tachangensis

BACKGROUND

Camellia tachangensis F. C. Zhang is an endemic Camellia species of the junction of Yunnan, Guizhou and Guangxi Provinces in China. It is characterized by a primitive five-chambered ovary morphology and serves as the botanical source of the renowned "Pu'an Red Tea". Unfortunately, the populations of the species have declined due to the destruction of their habitats by human activities. The lack of mitochondrial genomic resources has hindered research into molecular breeding and phylogenetic evolution of C. tachangensis.

RESULT

In this study, we had sequenced, assembled, and annotated the mitochondrial genome of C. tachangensis to reveal its genetic characteristics and phylogenetic relation with other Camellia species. The assembly result indicated that the mitochondrial genome sequence of C. tachangensis was 746,931 bp (GC content = 45.86%). It consisted of one multibranched sequence (Chr1) and one circular sequence (Chr2), with Chr1 capable of producing 7 substructures. The comparative analysis of the mitochondrial and chloroplast DNA of C. tachangensis revealed 23 pairs of chloroplast homologous fragments, with 10 fully preserved tRNA genes within them. Interspecies comparison of Ka/Ks ratios revealed that mutations in mitochondrial protein-coding genes (PCGs) of C. tachangensis were predominantly shaped by purifying selection throughout its evolution (Ka/Ks < 1). The mitochondrial CDS-based phylogenetic tree indicated that within the Camellia lineage, C. tachangensis was phylogenetically independent of the species of sections Oleifera, Camellia, Heterogenea, and Chrysantha. However, it also did not support the clustering of C. tachangensis with certain variants of C. sinensis, due to the extremely low support (BS = 22, PP = 0.41). Meanwhile, the chloroplast PCG-based phylogenetic analysis revealed that C. tachangensis formed a strongly supported basal clade (BS = 100, PP = 1.00), alongside C. makuanica (NC_087766), C. taliensis (NC_022264), and C. gymnogyna (NC_039626).

CONCLUSIONS

Our study deciphered the mitochondrial genome and its multibranched structure of C. tachangensis. These findings not only enhanced our comprehension of the complexity and diversity of mitochondrial genome structures in Camellia species, but also established a foundational genetic data framework for future research on molecular breeding programs and phylogenetic relationship involving C. tachangensis and its related species.

De novo RNA base editing in plant organelles with engineered synthetic P-type PPR editing factors

In plant mitochondria and chloroplasts, cytidine-to-uridine RNA editing is necessary for the production of functional proteins. While natural PLS-type PPR proteins are specialized in this process, synthetic PPR proteins offer significant potential for targeted RNA editing. In this study, we engineered chimeric editing factors by fusing synthetic P-type PPR guides with the DYW cytidine deaminase domain of a moss mitochondrial editing factor, PPR56. These designer PPR editors (dPPRe) elicited efficient and precise de novo RNA editing in Escherichia coli as well as in the chloroplasts and mitochondria of Nicotiana benthamiana. Chloroplast transcriptome-wide analysis of the most efficient dPPRe revealed minimal off-target effects, with only three nontarget C sites edited due to sequence similarity with the intended target. This study introduces a novel and precise method for RNA base editing in plant organelles, paving the way for new approaches in gene regulation applicable to plants and potentially other organisms.

Mitochondrial genome assembly of the Chinese endemic species of Camellia luteoflora and revealing its repetitive sequence mediated recombination, codon preferences and MTPTs

Camellia luteoflora Y.K. Li ex Hung T. Chang & F.A. Zeng belongs to the Camellia L. genus (Theaceae Mirb.). As an endemic, rare, and critically endangered species in China, it holds significant ornamental and economic value, garnering global attention due to its ecological rarity. Despite its conservation importance, genomic investigations on this species remain limited, particularly in organelle genomics, hindering progress in phylogenetic classification and population identification. In this study, we employed high-throughput sequencing to assemble the first complete mitochondrial genome of C. luteoflora and reannotated its chloroplast genome. Through integrated bioinformatics analyses, we systematically characterized the mitochondrial genome's structural organization, gene content, interorganellar DNA transfer, sequence variation, and evolutionary relationships.Key findings revealed a circular mitochondrial genome spanning 587,847 bp with a GC content of 44.63%. The genome harbors70 unique functional genes, including 40 protein-coding genes (PCGs), 27 tRNA genes, and 3 rRNA genes. Notably, 9 PCGs contained 22 intronic regions. Codon usage analysis demonstrated a pronounced A/U bias in synonymous codon selection. Structural features included 506 dispersed repeats and 240 simple sequence repeats. Comparative genomics identified 19 chloroplast-derived transfer events, contributing 29,534 bp (3.77% of total mitochondrial DNA). RNA editing prediction revealed 539 C-to-T conversion events across PCGs. Phylogenetic reconstruction using mitochondrial PCGs positioned C. luteoflora in closest evolutionary proximity to Camellia sinensis var. sinensis. Selection pressure analysis (Ka/Ks ratios < 1 for 11 PCGs) and nucleotide diversity assessment (Pi values: 0-0.00711) indicated strong purifying selection and low sequence divergence.This study provides the first comprehensive mitochondrial genomic resource for C. luteoflora, offering critical insights for germplasm conservation, comparative organelle genomics, phylogenetic resolution, and evolutionary adaptation studies in Camellia species.

The first complete mitochondrial genome assembly and comparative analysis of the fern Blechnaceae family: Blechnopsis orientalis

Introduction

Blechnopsis orientalis (L.) C. Presl is a medicinal and edible fern species belonging to the Blechnaceae family. Currently, the complete mitochondrial genome of B. orientalis, as well as those of other Blechnaceae species, remains unreported, and studies on fern mitochondrial genome are limited.

Methods

In this study, the B. orientalis mitochondrial genome was sequenced using both Nanopore PromethION and Illumina NovaSeq 6000 platforms. Genome annotation was performed using MITOFY and MFANNOT, with structural visualization via OGDRAW. In-depth analyses were conducted, including assessments of non-synonymous/synonymous mutation ratios (Ka/Ks), codon usage bias, repeat sequence identification, RNA editing site prediction, collinearity, and the identification of homologous fragments between chloroplast and mitochondrial genomes. Finally, we employed both the maximum likelihood (ML) and Bayesian (BI) methods to analyze the phylogenetic relationships among B. orientalis and nine other fern and lycophyte species.

Results

The mitochondrial genome of B. orientalis has a complex structure comprising 80 contigs, with a total length of 501,663 bp and a GC content of 48.53%. A total of 179 genes were identified, including 40 protein-coding genes (PCGs), 98 tRNA genes, 40 rRNA genes, and one pseudogene (rps11). Phylogenetic analysis based on PCGs from both chloroplast genome and mitochondrial genome aligned with the relationships described in the Pteridophyte Phylogeny Group I (PPG I) system. Further comparison with mitochondrial genome of ten other reported fern and lycophyte species revealed that the mitochondrial genome PCGs in these plants are highly conserved, despite significant genome rearrangements among mitochondrial genome.

Discussion

The findings of this study provide valuable insights into the evolutionary analysis of B. orientalis and contribute to understanding the characteristics and evolutionary relationships of mitochondrial genome in ferns and lycophytes.

Mitogenome of Uncaria rhynchophylla: genome structure, characterization, and phylogenetic relationships

BACKGROUND

Uncaria rhynchophylla is listed in the Chinese pharmacopoeia as one of the five botanical sources of the traditional medicine Gou-Teng, which has been utilized for the treatment of mental and cardiovascular disorders. This particular species is well-known in China for its hook-like structures originating from the leaf axils. Despite available reports on its chloroplast genome, there persists a notable lack of understanding concerning the structural variations and evolution of its mitochondrial genome. This knowledge gap hinders our ability to fully comprehend its genetic attributes.

RESULTS

We successfully assembled the mitochondrial genome of U. rhynchophylla by seamlessly integrating Illumina short reads with Nanopore long reads, resulting in a non-circular genome comprising 1 circular contig and 2 linear contigs. The total length of this genome is 421,660 bp, encompassing 36 PCGs. The identification of 4 distinct pairs of repeats has unveiled their pivotal role in repeat-mediated recombination. Of the 28 homologous fragments derived from chloroplasts, the majority were observed to have been transferred from the inverted repeat (IR) regions of the chloroplast genome to the mitochondrial genome. The mitochondrial DNA provides a distinctive resolution for the positioning of several species within the Gentianales phylogenetic framework, which remains unresolved by chloroplast DNA.

CONCLUSION

By utilizing a newly assembled, high-quality mitochondrial genome of U. rhynchophylla, we have elucidated its intricate genomic structure, distinctive sequence characteristics, and potential for phylogenetic analysis. These findings mark significant strides in advancing our comprehension of the genetics of Uncaria.

De novo assembly and characterization of the complete mitochondrial genome of Phellodendron amurense reveals three repeat-mediated recombination

Phellodendron amurense Rupr., a rare herb renowned for its medicinal and ecological significance, has remained genetically unexplored at the mitochondrial level until now. This study presents the first-ever systematic assembly and annotation of the complete mitochondrial genome of P. amurense, achieved through a hybrid strategy combining Illumina and Nanopore sequencing data. The mitochondrial genome spans 566,285 bp with a GC content of 45.51 %, structured into two circular molecules. Our comprehensive analysis identified 32 protein-coding genes (PCGs), 33 tRNA genes, and 3 rRNA genes, alongside 181 simple sequence repeats, 19 tandem repeats, and 310 dispersed repeats. Notably, multiple genome conformations were predicted due to repeat-mediated homologous recombination. Additionally, we assembled the chloroplast genome, identifying 21 mitochondrial plastid sequences that provide insights into organelle genome interactions. A total of 380 RNA-editing sites within the mitochondrial PCGs were predicted, enhancing our understanding of gene regulation and function. Phylogenetic analysis using mitochondrial PCGs from 30 species revealed evolutionary relationships, confirming the homology between P. amurense and Citrus species. This foundational study offers a valuable genetic resource for the Rutaceae family, facilitating further research into genetic evolution and molecular diversity in plant mitochondrial genomes.

Plastome data provides new insights into population differentiation and evolution of Ginkgo in the Sichuan Basin of China

BACKGROUND

Ginkgo biloba L., an iconic living fossil, challenges traditional views of evolutionary stasis. While nuclear genomic studies have revealed population structure across China, the evolutionary patterns reflected in maternally inherited plastomes remain unclear, particularly in the Sichuan Basin - a potential glacial refugium that may have played a crucial role in Ginkgo's persistence.

RESULTS

Analysis of 227 complete plastomes, including 81 newly sampled individuals from the Sichuan Basin, revealed three distinct maternal lineages differing from known nuclear genome patterns. We identified 170 sequence variants and extensive RNA editing (235 sites) with a bias toward hydrophobic amino acid conversions, suggesting active molecular evolution. A previously undocumented haplotype (IIA2), predominant in western Sichuan Basin populations, showed close genetic affinity with rare refugial haplotypes. Western populations exhibited higher haplotypic diversity and distinctive genetic structure, supporting the basin's role as both glacial refugium and corridor for population expansion. Ancient trees (314-784 years) provided evidence for interaction between natural processes and historical human dispersal in shaping current genetic patterns.

CONCLUSIONS

Our findings demonstrate substantial genetic diversity within Sichuan Basin Ginkgo populations and reveal dynamic molecular evolution through plastome variation and RNA editing patterns, challenging the notion of evolutionary stasis in this living fossil. This study provides crucial genomic resources for understanding Ginkgo's evolution and informs conservation strategies for this endangered species.

Comparative analysis of mitochondrial genomes of Stemona tuberosa lour. reveals heterogeneity in structure, synteny, intercellular gene transfer, and RNA editing

BACKGROUND

Stemona tuberosa, a vital species in traditional Chinese medicine, has been extensively cultivated and utilized within its natural distribution over the past decades. While the chloroplast genome of S. tuberosa has been characterized, its mitochondrial genome (mitogenome) remains unexplored.

RESULTS

This paper details the assembly of the complete S. tuberosa mitogenome, achieved through the integration of Illumina and Nanopore sequencing technologies. The assembled mitogenome is 605,873 bp in size with a GC content of 45.63%. It comprises 66 genes, including 38 protein-coding genes, 25 tRNA genes, and 3 rRNA genes. Our analysis delved into codon usage, sequence repeats, and RNA editing within the mitogenome. Additionally, we conducted a phylogenetic analysis involving S. tuberosa and 17 other taxa to clarify its evolutionary and taxonomic status. This study provides a crucial genetic resource for evolutionary research within the genus Stemona and other related genera in the Stemonaceae family.

CONCLUSION

Our study provides the inaugural comprehensive analysis of the mitochondrial genome of S. tuberosa, revealing its unique multi-branched structure. Through our investigation of codon usage, sequence repeats, and RNA editing within the mitogenome, coupled with a phylogenetic analysis involving S. tuberosa and 17 other taxa, we have elucidated its evolutionary and taxonomic status. These investigations provide a crucial genetic resource for evolutionary research within the genus Stemona and other related genera in the Stemonaceae family.

De novo assembly of the mitochondrial genome of Glycyrrhiza glabra and identification of two types of homologous recombination configurations caused by repeat sequences

BACKGROUND

Glycyrrhiza glabra, which is widely used in medicine and therapy, is known as the 'king of traditional Chinese medicine'. In this study, we successfully assembled and annotated the mitochondrial and chloroplast genomes of G. glabra via high-throughput sequencing technology, combining the advantages of short-read (Illumina) and long-read (Oxford Nanopore) sequencing.

RESULTS

We revealed the ring structure of the mitochondrial genome, which spans 421,293 bp with 45.1% GC content and 56 annotated genes. Notably, we identified 514 repetitive sequences, including 123 Simple sequence repeats (SRs), 3 Tndem sequence repeats (TSRs), and 388 Dispersed sequence repeats (DSRs). We identified 79 out of the 388 DSRs as potentially involved in homologous recombination. We identified five forward repeats and four palindromic repeats that facilitate homologous recombination and induce alterations in the mitochondrial genome structure. We corroborated this finding via polymerase chain reaction (PCR). Furthermore, we identified chloroplast-derived sequence fragments within the mitochondrial genome, offering novel insights into the evolutionary history of plant mitochondrial genomes. We predicted 460 potential RNA editing sites, primarily involving cytosine-to-uracil transitions. This study reveals the complexity of repetitive sequence-mediated homologous recombination in the mitochondrial genome of G. glabra and provides new insights into its structure, function, and evolution.

CONCLUSIONS

These findings have important implications for conservation biology, population genetics, and evolutionary studies, underscoring the role of repetitive sequences in genome dynamics and highlighting the need for further research on mitochondrial genome evolution and function in plants.

ncRNA Editing: Functional Characterization and Computational Resources

Non-coding RNAs (ncRNAs) play crucial roles in gene expression regulation, translation, and disease development, including cancer. They are classified by size in short and long non-coding RNAs. This chapter focuses on the functional implications of adenosine-to-inosine (A-to-I) RNA editing in both short (e.g., miRNAs) and long ncRNAs. RNA editing dynamically alters the sequence and structure of primary transcripts, impacting ncRNA biogenesis and function. Notable findings include the role of miRNA editing in promoting glioblastoma invasiveness, characterizing RNA editing hotspots across cancers, and its implications in thyroid cancer and ischemia. This chapter also highlights bioinformatics resources and next-generation sequencing (NGS) technologies that enable comprehensive ncRNAome studies and genome-wide RNA editing detection. Dysregulation of RNA editing machinery has been linked to various human diseases, emphasizing the potential of RNA editing as a biomarker and therapeutic target. This overview integrates current knowledge and computational tools for studying ncRNA editing, providing insights into its biological significance and clinical applications.

Unraveling the organellar genomic landscape of the therapeutic and entheogenic plant Mimosa tenuiflora: insights into genetic, structural, and evolutionary dynamics

Mimosa tenuiflora, popularly known as "Jurema-Preta", is a perennial tree or shrub native to the tropical regions of the Americas, particularly among Afro-Brazilian and Indigenous Brazilian communities. Known for producing N,N-Dimethyltryptamine, a psychedelic compound with profound psychological effects, Jurema-Preta has been studied for its therapeutic potential in mental health. This study offers a comprehensive analysis of the plastid (ptDNA) and mitochondrion (mtDNA) genomes of M. tenuiflora. The 165,639 bp ptDNA sequence features the classical quadripartite structure with 130 protein-coding genes. Comparative genomics among Mimosa species shows high sequence identity in protein-coding genes, with variation in the rpoC1, clpP, ndhA, and ycf1 genes. The ptDNA junctions display distinct features, such as the deletion of the rpl22 gene, and specific simple sequence repeats highlight genetic variation and unique motifs as valuable genetic markers for population studies. Phylogenetic analysis places M. tenuiflora in the Caesalpinioideae, closely related to M. pigra and M. pudica. The 617,839 bp mtDNA sequence exhibits a complex structure with multiple genomic arrangements due to large repeats, encoding 107 protein-coding genes, including the ptDNA petG and psaA genes, and non-retroviral RNA mitoviruses sequences. Comparative analysis across Fabaceae species reveals limited conservation, emphasizing the dynamic nature of plant mitochondrial genomes. The genomic characterization of M. tenuiflora enhances understanding of its evolutionary dynamics, providing insights for population studies and potential applications in ethnopharmacology and conservation.

Deciphering the multi- partite mitochondrial genome of Crataegus pinnatifida: insights into the evolution and genetics of cultivated Hawthorn

Flowering plant (angiosperm) mitochondrial genomes are remarkably dynamic in their structures. We present the complete mitochondrial genome of hawthorn (Crataegus pinnatifida Bunge), a shrub that bears fruit and is celebrated for its extensive medicinal history. We successfully assembled the hawthorn mitogenome utilizing the PacBio long-read sequencing technique, which yielded 799,862 reads, and the Illumina novaseq6000 sequencing platform, which producing 6.6 million raw paired reads. The C. pinnatifida mitochondria sequences encompassed a total length of 440,295 bp with a GC content of 45.42%. The genome annotates 54 genes, including 34 that encode proteins, 17 that encode tRNA, and three genes for rRNA. A fascinating interplay was observed between the chloroplast and mitochondrial genomes, which share 17 homologous sequences sequences that rotal 1,933 bp. A total of 134 SSRs, 22 tandem repeats and 42 dispersed repeats were identified in the mitogenome. Four conformations of C. pinnatifida mitochondria sequences recombination were verified through PCR experiments and Sanger sequencing, and C. pinnatifida mitogenome is more likely to be assembled into three circular-mapping chromosomes. All the RNA editing sites that were identified C-U edits, which predominantly occurred at the first and second positions of the codons. Phylogenetic and collinearity analyses identified the evolutionary trajectory of C. pinnatifida, which reinforced the genetic identity of the hawthorn section. This unveiling of the unique multi-partite structure of the hawthorn mitogenome offers a foundational reference for future study into the evolution and genetics of C. pinnatifida.

Comparative analysis of the complete mitogenomes of Camellia sinensis var. sinensis and C. sinensis var. assamica provide insights into evolution and phylogeny relationship

Introduction

Among cultivated tea plants (Camellia sinensis), only four mitogenomes for C. sinensis var. assamica (CSA) have been reported so far but none for C. sinensis var. sinensis (CSS). Here, two mitogenomes of CSS (CSSDHP and CSSRG) have been sequenced and assembled.

Methods

Using a combination of Illumina and Nanopore data for the first time. Comparison between CSS and CSA mitogenomes revealed a huge heterogeneity.

Results

The number of the repetitive sequences was proportional to the mitogenome size and the repetitive sequences dominated the intracellular gene transfer segments (accounting for 88.7%- 92.8% of the total length). Predictive RNA editing analysis revealed that there might be significant editing in NADH dehydrogenase subunit transcripts. Codon preference analysis showed a tendency to favor A/T bases and T was used more frequently at the third base of the codon. ENc plots analysis showed that the natural selection play an important role in shaping the codon usage bias, and Ka/Ks ratios analysis indicated Nad1 and Sdh3 genes may have undergone positive selection. Further, phylogenetic analysis shows that six C. sinensis clustered together, with the CSA and CSS forming two distinct branches, suggesting two different evolutionary pathway.

Discussion

Altogether, this investigation provided an insight into evolution and phylogeny relationship of C. sinensis mitogenome, thereby enhancing comprehension of the evolutionary patterns within C. sinensis species.

Comparative analysis of the mitochondrial genomes of four Dendrobium species (Orchidaceae) reveals heterogeneity in structure, synteny, intercellular gene transfer, and RNA editing

The genus Dendrobium, part of the Orchidaceae family, encompasses species of significant medicinal, nutritional, and economic value. However, many Dendrobium species are threatened by environmental stresses, low seed germination rates, and overharvesting. Mitochondria generate the energy necessary for various plant life activities. Despite their importance, research on the mitochondrial genomes of Dendrobium species is currently limited. To address this gap, we performed a comprehensive genetic analysis of four Dendrobium species-D. flexicaule, D. nobile, D. officinale, and D. huoshanense-focusing on their mitochondrial and chloroplast genomes to elucidate their genetic architecture and support conservation efforts. We utilized advanced sequencing technologies, including Illumina for high-throughput sequencing and Nanopore for long-read sequencing capabilities. Our findings revealed the multichromosomal mitochondrial genome structures, with total lengths ranging from 596,506 bp to 772,523 bp. The mitochondrial genomes contained 265 functional genes, including 64-69 protein-coding genes, 23-28 tRNA genes, and 3 rRNA genes. We identified 647 simple sequence repeats (SSRs) and 352 tandem repeats, along with 440 instances of plastid-to-mitochondrial gene transfer. Additionally, we predicted 2,023 RNA editing sites within the mitochondrial protein-coding genes, predominantly characterized by cytosine-to-thymine transitions. Comparative analysis of mitochondrial DNA across the species highlighted 25 conserved genes, with evidence of positive selection in five genes: ccmFC, matR, mttB, rps2, and rps10. Phylogenetic assessments suggested a close sister relationship between D. nobile and D. huoshanense, and a similar proximity between D. officinale and D. flexicaule. This comprehensive genomic study provides a critical foundation for further exploration into the genetic mechanisms and biodiversity of Dendrobium species, contributing valuable insights for their conservation and sustainable utilization.

Conquering new grounds: plant organellar C-to-U RNA editing factors can be functional in the plant cytosol

Plant mitochondrial and chloroplast transcripts are subject to numerous events of specific cytidine-to-uridine (C-to-U) RNA editing to correct genetic information. Key protein factors for this process are specific RNA-binding pentatricopeptide repeat (PPR) proteins, which are encoded in the nucleus and post-translationally imported into the two endosymbiotic organelles. Despite hundreds of C-to-U editing sites in the plant organelles, no comparable editing has been found for nucleo-cytosolic mRNAs raising the question why plant RNA editing is restricted to chloroplasts and mitochondria. Here, we addressed this issue in the model moss Physcomitrium patens, where all PPR-type RNA editing factors comprise specific RNA-binding and cytidine deamination functionalities in single proteins. To explore whether organelle-type RNA editing can principally also take place in the plant cytosol, we expressed PPR56, PPR65 and PPR78, three editing factors recently shown to also function in a bacterial setup, together with cytosolic co-transcribed native targets in Physcomitrium. While we obtained unsatisfying results upon their constitutive expression, we found strong cytosolic RNA editing under hormone-inducible expression. Moreover, RNA-Seq analyses revealed varying numbers of up to more than 900 off-targets in other cytosolic transcripts. We conclude that PPR-mediated C-to-U RNA editing is not per se incompatible with the plant cytosol but that its limited target specificity has restricted its occurrence to the much less complex transcriptomes of mitochondria and chloroplast in the course of evolution.

Integration of large and diverse angiosperm DNA fragments into Asian Gnetum mitogenomes

BACKGROUND

Horizontal gene transfer (HGT) events have rarely been reported in gymnosperms. Gnetum is a gymnosperm genus comprising 25‒35 species sympatric with angiosperms in West African, South American, and Southeast Asian rainforests. Only a single acquisition of an angiosperm mitochondrial intron has been documented to date in Asian Gnetum mitogenomes. We wanted to develop a more comprehensive understanding of frequency and fragment length distribution of such events as well as their evolutionary history in this genus.

RESULTS

We sequenced and assembled mitogenomes from five Asian Gnetum species. These genomes vary remarkably in size and foreign DNA content. We identified 15 mitochondrion-derived and five plastid-derived (MTPT) foreign genes. Our phylogenetic analyses strongly indicate that these foreign genes were transferred from diverse eudicots-mostly from the Rubiaceae genus Coptosapelta and ten genera of Malpighiales. This indicates that Asian Gnetum has experienced multiple independent HGT events. Patterns of sequence evolution strongly suggest DNA-mediated transfer between mitochondria as the primary mechanism giving rise to these HGT events. Most Asian Gnetum species are lianas and often entwined with sympatric angiosperms. We therefore propose that close apposition of Gnetum and angiosperm stems presents opportunities for interspecific cell-to-cell contact through friction and wounding, leading to HGT.

CONCLUSIONS

Our study reveals that multiple HGT events have resulted in massive amounts of angiosperm mitochondrial DNA integrated into Asian Gnetum mitogenomes. Gnetum and its neighboring angiosperms are often entwined with each other, possibly accounting for frequent HGT between these two phylogenetically remote lineages.

Comparative chloroplast genomes study of five officinal Ardisia Species: Unraveling interspecific diversity and evolutionary insights in Ardisia

Ardisia S.W. (Primulaceae), naturally distributed in tropical and subtropical regions, has edible and medicinal values and is prevalent in clinical and daily use in China. More genetic information for distinct species delineation is needed to support the development and utilization of the genus Ardisia. We sequenced, annotated, and compared the chloroplast genomes of five Ardisia species: A. brunnescens, A. pusilla, A. squamulosa, A. crenata, and A. brevicaulis in this study. We found a typical quadripartite structure in all five chloroplast genomes, with lengths ranging from 155,045 to 156,943 bp. Except for A. pusilla, which lacked the ycf15 gene, the other four Ardisia species contained 114 unique genes, including 79 protein-coding genes, 30 tRNAs, and four rRNAs. In addition, the rps19 pseudogene gene was present only in A. brunnescens. Five highly variable DNA barcodes were identified for five Ardisia species, including trnT-GGU-psbD, trnT-UGU-trnL-UAA, rps4-trnT-UGU, rpl32-trnL-UAG, and rpoB-trnC-GAA. The RNA editiing sites of protein-coding genes in the five Ardisia plastome were characterized and compared, and 274 (A. crenata)-288 (A. brevicaulis) were found. The results of the phylogenetic analysis were consistent with the morphological classification. Sequence alignment and phylogenetic analysis showed that ycf15 genes were highly divergent in Primulaceae. Reconstructions of ancestral character states indicated that leaf margin morphology is critical for classifying the genus Ardisia, with a rodent-like character being the most primitive. These results provide valuable information on the taxonomy and evolution of Ardisia plants.

Assembly and comparative analysis of the first complete mitochondrial genome of Setaria italica

MAIN CONCLUSION

In this study, we assembled the first complete mitochondrial genome of Setaria italica and confirmed the multi-branched architecture. The foxtail millet (Setaria italica) holds significant agricultural importance, particularly in arid and semi-arid regions. It plays a pivotal role in diversifying dietary patterns and shaping planting strategies. Although the chloroplast genome of S. italica has been elucidated in recent studies, the complete mitochondrial genome remains largely unexplored. In this study, we employed PacBio HiFi sequencing platforms to sequence and assemble the complete mitochondrial genome. The mitochondrial genome spans a total length of 446,614 base pairs and harbors a comprehensive set of genetic elements, including 33 unique protein-coding genes (PCGs), encompassing 24 unique mitochondrial core genes and 9 variable genes, along with 20 transfer RNA (tRNA) genes and 3 ribosomal RNA (rRNA) genes. Our analysis of mitochondrial PCGs revealed a pronounced codon usage preference. For instance, the termination codon exhibits a marked preference for UAA, while alanine (Ala) exhibits a preference for GCU, and glutamine (Gln) favors CAA. Notably, the maximum Relative Synonymous Codon Usage (RSCU) values for cysteine (Cys) and phenylalanine (Phe) are both below 1.2, indicating a lack of strong codon usage preference for these amino acids. Phylogenetic analyses consistently place S. italica in close evolutionary proximity to Chrysopogon zizanioides, relative to other Panicoideae plants. Collinearity analysis showed that a total of 39 fragments were identified to display homology with both the mitochondrial and chloroplast genomes. A total of 417 potential RNA-editing sites were discovered across the 33 mitochondrial PCGs. Notably, all these editing events involved the conversion of cytosine (C) to uracil (U). Through the employment of PCR validation coupled with Sanger sequencing for the anticipated editing sites of these codons, RNA-editing events were conclusively identified at two specific loci: nad4L-2 and atp6-1030. The results of this study provide a pivotal foundation for advanced genomic breeding research in foxtail millet. Furthermore, they impart essential insights that will be instrumental for forthcoming investigations into the evolutionary and molecular dynamics of Panicoideae species.

Deep proteomics reveals incorporation of unedited proteins into mitochondrial protein complexes in Arabidopsis

The mitochondrial proteome consists of numerous types of proteins which either are encoded and synthesized in the mitochondria, or encoded in the cell nucleus, synthesized in the cytoplasm and imported into the mitochondria. Their synthesis in the mitochondria, but not in the nucleus, relies on the editing of the primary transcripts of their genes at defined sites. Here, we present an in-depth investigation of the mitochondrial proteome of Arabidopsis (Arabidopsis thaliana) and a public online platform for the exploration of the data. For the analysis of our shotgun proteomic data, an Arabidopsis sequence database was created comprising all available protein sequences from the TAIR10 and Araport11 databases, supplemented with sequences of proteins translated from edited and nonedited transcripts of mitochondria. Amino acid sequences derived from partially edited transcripts were also added to analyze proteins encoded by the mitochondrial genome. Proteins were digested in parallel with six different endoproteases to obtain maximum proteome coverage. The resulting peptide fractions were finally analyzed using liquid chromatography coupled to ion mobility spectrometry and tandem mass spectrometry. We generated a "deep mitochondrial proteome" of 4,692 proteins. 1,339 proteins assigned to mitochondria by the SUBA5 database (https://suba.live) accounted for >80% of the total protein mass of our fractions. The coverage of proteins by identified peptides was particularly high compared to single-protease digests, allowing the exploration of differential splicing and RNA editing events at the protein level. We show that proteins translated from nonedited transcripts can be incorporated into native mitoribosomes and the ATP synthase complex. We present a portal for the use of our data, based on "proteomaps" with directly linked protein data. The portal is available at www.proteomeexplorer.de.

Multiple factors interact in editing of PPR-E+-targeted sites in maize mitochondria and plastids

RNA cytidine-to-uridine editing is essential for plant organellar gene expression. Pentatricopeptide repeat (PPR)-E+ proteins have been proposed to bind to target sites and recruit the cytidine deaminase AtDYW2, facilitated by AtNUWA. Here we analyze the function of ZmNUWA, ZmDYW2A, and ZmDYW2B and their relationships with other editing factors in maize. The zmdyw2a and zmdyw2b single mutants are normal, but the zmdyw2a::zmdyw2b and zmnuwa mutants are severely arrested in seed development. ZmNUWA, ZmDYW2A, and ZmDYW2B are dual localized in mitochondria and plastids. Loss of ZmNUWA decreases the editing at 99 mitochondrial sites and 8 plastid sites. Surprisingly, loss of ZmDYW2A:ZmDYW2B affects almost the same set of sites targeted by PPR-E+ proteins. ZmNUWA interacts with ZmDYW2A and ZmDYW2B, suggesting that ZmNUWA recruits ZmDYW2A/2B in the editing of PPR-E+-targeted sites in maize. Further protein interaction analyses show that ZmNUWA and ZmDYW2A/2B interact with ZmMORF1, ZmMORF8, ZmMORF2, and ZmMORF9 and that ZmOZ1 interacts with ZmORRM1, ZmDYW2A, ZmDYW2B, ZmMORF8, and ZmMORF9. These results suggest that the maize mitochondrial PPR-E+ editosome contains PPR-E+, ZmDYW2A/2B, ZmNUWA, and ZmMORF1/8, whereas the plastid PPR-E+ editosome is composed of PPR-E+, ZmDYW2A/2B, ZmNUWA, ZmMORF2/8/9, ZmORRM1, and ZmOZ1.

PlantC2U: deep learning of cross-species sequence landscapes predicts plastid C-to-U RNA editing in plants

In plants, C-to-U RNA editing mainly occurs in plastid and mitochondrial transcripts, which contributes to a complex transcriptional regulatory network. More evidence reveals that RNA editing plays critical roles in plant growth and development. However, accurate detection of RNA editing sites using transcriptome sequencing data alone is still challenging. In the present study, we develop PlantC2U, which is a convolutional neural network, to predict plastid C-to-U RNA editing based on the genomic sequence. PlantC2U achieves >95% sensitivity and 99% specificity, which outperforms the PREPACT tool, random forests, and support vector machines. PlantC2U not only further checks RNA editing sites from transcriptome data to reduce possible false positives, but also assesses the effect of different mutations on C-to-U RNA editing based on the flanking sequences. Moreover, we found the patterns of tissue-specific RNA editing in the mangrove plant Kandelia obovata, and observed reduced C-to-U RNA editing rates in the cold stress response of K. obovata, suggesting their potential regulatory roles in plant stress adaptation. In addition, we present RNAeditDB, available online at https://jasonxu.shinyapps.io/RNAeditDB/. Together, PlantC2U and RNAeditDB will help researchers explore the RNA editing events in plants and thus will be of broad utility for the plant research community.

Detection and editing of the updated Arabidopsis plastid- and mitochondrial-encoded proteomes through PeptideAtlas

Arabidopsis (Arabidopsis thaliana) ecotype Col-0 has plastid and mitochondrial genomes encoding over 100 proteins. Public databases (e.g. Araport11) have redundancy and discrepancies in gene identifiers for these organelle-encoded proteins. RNA editing results in changes to specific amino acid residues or creation of start and stop codons for many of these proteins, but the impact of RNA editing at the protein level is largely unexplored due to the complexities of detection. Here, we assembled the nonredundant set of identifiers, their correct protein sequences, and 452 predicted nonsynonymous editing sites of which 56 are edited at lower frequency. We then determined accumulation of edited and/or unedited proteoforms by searching ∼259 million raw tandem MS spectra from ProteomeXchange, which is part of PeptideAtlas (www.peptideatlas.org/builds/arabidopsis/). We identified all mitochondrial proteins and all except 3 plastid-encoded proteins (NdhG/Ndh6, PsbM, and Rps16), but no proteins predicted from the 4 ORFs were identified. We suggest that Rps16 and 3 of the ORFs are pseudogenes. Detection frequencies for each edit site and type of edit (e.g. S to L/F) were determined at the protein level, cross-referenced against the metadata (e.g. tissue), and evaluated for technical detection challenges. We detected 167 predicted edit sites at the proteome level. Minor frequency sites were edited at low frequency at the protein level except for cytochrome C biogenesis 382 at residue 124 (Ccb382-124). Major frequency sites (>50% editing of RNA) only accumulated in edited form (>98% to 100% edited) at the protein level, with the exception of Rpl5-22. We conclude that RNA editing for major editing sites is required for stable protein accumulation.

Conservation of the moss RNA editing factor PPR78 despite the loss of its known cytidine-to-uridine editing sites is explained by a hidden extra target

Cytidine (C)-to-uridine (U) RNA editing in plant organelles relies on specific RNA-binding pentatricopeptide repeat (PPR) proteins. In the moss Physcomitrium patens, all such RNA editing factors feature a C-terminal DYW domain that acts as the cytidine deaminase for C-to-U conversion. PPR78 of Physcomitrium targets 2 mitochondrial editing sites, cox1eU755SL and rps14eU137SL. Remarkably, the latter is edited to highly variable degrees in different mosses. Here, we aimed to unravel the coevolution of PPR78 and its 2 target sites in mosses. Heterologous complementation in a Physcomitrium knockout line revealed that the variable editing of rps14eU137SL depends on the PPR arrays of different PPR78 orthologues but not their C-terminal domains. Intriguingly, PPR78 has remained conserved despite the simultaneous loss of editing at both known targets among Hypnales (feather mosses), suggesting it serves an additional function. Using a recently established RNA editing assay in Escherichia coli, we confirmed site-specific RNA editing by PPR78 in the bacterium and identified 4 additional off-targets in the bacterial transcriptome. Based on conservation profiles, we predicted ccmFNeU1465RC as a candidate editing target of PPR78 in moss mitochondrial transcriptomes. We confirmed editing at this site in several mosses and verified that PPR78 targets ccmFNeU1465RC in the bacterial editing system, explaining the conservation and functional adaptation of PPR78 during moss evolution.

RNA editing events and expression profiles of mitochondrial protein-coding genes in the endemic and endangered medicinal plant, Corydalis saxicola

Corydalis saxicola, an endangered medicinal plant endemic to karst habitats, is widely used in Traditional Chinese Medicine to treat hepatitis, abdominal pain, bleeding hemorrhoids and other conditions. However, to date, the mitochondrial (mt) genome of C. saxicola has not been reported, which limits our understanding of the genetic and biological mechanisms of C. saxicola. Here, the mt genome of C. saxicola was assembled by combining the Nanopore and Illumina reads. The mt genome of C. saxicola is represented by a circular chromosome which is 587,939 bp in length, with an overall GC content of 46.50%. 40 unique protein-coding genes (PCGs), 22 tRNA genes and three rRNA genes were identified. Codon usage of the PCGs was investigated and 167 simple sequence repeats were identified. Twelve homologous fragments were identified between the mt and ct genomes of C. saxicola, accounting for 1.04% of the entire mt genome. Phylogenetic examination of the mt genomes of C. saxicola and 30 other taxa provided an understanding of their evolutionary relationships. We also predicted 779 RNA editing sites in 40 C. saxicola mt PCGs and successfully validated 506 (65%) of these using PCR amplification and Sanger sequencing. In addition, we transcriptionally profiled 24 core mt PCGs in C. saxicola roots treated with different concentrations of CaCl2, as well as in other organs. These investigations will be useful for effective utilization and molecular breeding, and will also provide a reference for further studies of the genus Corydalis.

Assembly and comparative analysis of the first complete mitochondrial genome of a traditional Chinese medicine Angelica biserrata (Shan et Yuan) Yuan et Shan

Duhuo, a member of the Angelica family, is widely used to treat ailments such as rheumatic pain. It possesses a diverse array of bioactivities, including anti-tumor, anti-inflammatory, and analgesic properties, as recent pharmacological research has revealed. Nevertheless, the mtDNA of Angelica species remains relatively unexplored. To address this gap, we sequenced and assembled the mtDNA of A. biserrata to shed light on its genetic mechanisms and evolutionary pathways. Our investigation indicated a distinctive multi-branched conformation in the A. biserrata mtDNA. A comprehensive analysis of protein-coding sequences (PCGs) across six closely related species revealed the presence of 11 shared genes in their mitochondrial genomes. Intriguingly, positive selection emerged as a significant factor in the evolution of the atp4, matR, nad3, and nad7 genes. In addition, our data highlighted a recurring trend of homologous fragment migration between chloroplast and mitochondrial organelles. We identified 13 homologous fragments spanning both chloroplast and mitochondrial genomes. The phylogenetic tree established a close relationship between A. biserrata and Saposhnikovia divaricata. To sum up, our research would contribute to the application of population genetics and evolutionary studies in the genus Acanthopanax and other genera in the Araliaceae family.

The complete chloroplast genome sequence of Nepeta bracteata and comparison with congeneric species

Nepeta bracteata (N. bracteata) is an important medicinal plant used by Chinese ethnic minorities. However, the lack of knowledge regarding the chloroplast genome of N. bracteata has imposed current limitations on our study. Here, we used Next-generation sequencing to obtain the chloroplast genome of N. bracteata. The findings suggested that the 151,588 bp cp genome of N. bracteata comprises 130 genes, including 35 tRNA genes and 87 protein-coding genes. And its chloroplast genome exhibits a typical quadripartite structure, the largest single copy (LSC; 82,819 bp) and the smallest single copy (SSC; 17,557 bp) separate a pair of inverted repeats IR regions (IRa and IRb; 25,606 bp) from one another. Interestingly, palindromic repeats are more common, as shown by the examination of repetition. In the interim, 18 SSRs were discovered in the interim, the bulk of which were Adenine-Thymine (A-T) mononucleotides. Meanwhile, we compared it with five other species from the Nepeta genus. Five hypervariable areas were found by the study, including ndhH-rps15, accD-psal, ndhG-ndhl, trnH-GUG-psbA, and rpoC1-rpoB. Furthermore, the phylogenetic study revealed that N. bracteata and Nepeta stewartiana (N. stewartiana) were linked to each other most closely. In summary, our findings enrich the resources available for chloroplast genomes in the Nepeta genus. Moreover, these hypervariable regions have the potential to be developed into molecular markers, enabling the rapid identification of species within the Nepeta genus. Comparative analysis of species within the Nepeta genus can help enhance our study of their phylogenetic relationships, potential medicinal properties and bioprospecting.

Complete chloroplast genomes of Cerastium alpinum, C. arcticum and C. nigrescens: genome structures, comparative and phylogenetic analysis

The genus Cerastium includes about 200 species that are mostly found in the temperate climates of the Northern Hemisphere. Here we report the complete chloroplast genomes of Cerastium alpinum, C. arcticum and C. nigrescens. The length of cp genomes ranged from 147,940 to 148,722 bp. Their quadripartite circular structure had the same gene organization and content, containing 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. Repeat sequences varied from 16 to 23 per species, with palindromic repeats being the most frequent. The number of identified SSRs ranged from 20 to 23 per species and they were mainly composed of mononucleotide repeats containing A/T units. Based on Ka/Ks ratio values, most genes were subjected to purifying selection. The newly sequenced chloroplast genomes were characterized by a high frequency of RNA editing, including both C to U and U to C conversion. The phylogenetic relationships within the genus Cerastium and family Caryophyllaceae were reconstructed based on the sequences of 71 protein-coding genes. The topology of the phylogenetic tree was consistent with the systematic position of the studied species. All representatives of the genus Cerastium were gathered in a single clade with C. glomeratum sharing the least similarity with the others.

Chloroplast genome sequence of Nothapodytes nimmoniana (Icacinaceae), a vulnerable reservoir of camptothecin from Western Ghats

Nothapodytes nimmoniana belongs to family Icacinaceae and is a major source of compound Camptothecin. The global demand for Camptothecin has caused large-scale exploitation of N. nimmoniana from its wild habitat in Western Ghats of India, thereby making it vulnerable. The species is known to exhibit genetic diversity among the populations in Western Ghats. In this study, we report plastome sequence of N. nimmoniana, first for the genus. For the study, the species was collected from Western Ghats of Maharashtra. The plastome of N. nimmoniana was 150,726 bp in length and exhibited typical quadripartite structure with 83,771 bp LSC, 18,513 bp SSC and 24,221 IR region. The plastome was characterized by presence of 124 unique genes, 87 protein coding genes, 29 tRNA genes and four rRNA genes. Further, the plastome was compared with the available basal lamiid plastomes for gene order and composition. N. nimmoniana plastome exhibited SSC region in an inverted configuration. Phylogenomic study placed N. nimmoniana sister to Mappia mexicana. The SSR markers identified in this study, might help to distinguish genetically diverse populations, prioritizing the populations which need immediate conservation effects as well as for checking adulteration.

Genetic diversity of Coffea arabica L. mitochondrial genomes caused by repeat- mediated recombination and RNA editing

Background

Coffea arabica L. is one of the most important crops widely cultivated in 70 countries across Asia, Africa, and Latin America. Mitochondria are essential organelles that play critical roles in cellular respiration, metabolism, and differentiation. C. arabica's nuclear and chloroplast genomes have been reported. However, its mitochondrial genome remained unreported. Here, we intended to sequence and characterize its mitochondrial genome to maximize the potential of its genomes for evolutionary studies, molecular breeding, and molecular marker developments.

Results

We sequenced the total DNA of C. arabica using Illumina and Nanopore platforms. We then assembled the mitochondrial genome with a hybrid strategy using Unicycler software. We found that the mitochondrial genome comprised two circular chromosomes with lengths of 867,678 bp and 153,529 bp, encoding 40 protein-coding genes, 26 tRNA genes, and three rRNA genes. We also detected 270 Simple Sequence Repeats and 34 tandem repeats in the mitochondrial genome. We found 515 high-scoring sequence pairs (HSPs) for a self-to-self similarity comparison using BLASTn. Three HSPs were found to mediate recombination by the mapping of long reads. Furthermore, we predicted 472 using deep-mt with the convolutional neural network model. Then we randomly validated 90 RNA editing events by PCR amplification and Sanger sequencing, with the majority being non-synonymous substitutions and only three being synonymous substitutions. These findings provide valuable insights into the genetic characteristics of the C. arabica mitochondrial genome, which can be helpful for future study on coffee breeding and mitochondrial genome evolution.

Conclusion

Our study sheds new light on the evolution of C. arabica organelle genomes and their potential use in genetic breeding, providing valuable data for developing molecular markers that can improve crop productivity and quality. Furthermore, the discovery of RNA editing events in the mitochondrial genome of C. arabica offers insights into the regulation of gene expression in this species, contributing to a better understanding of coffee genetics and evolution.

Assembly and analysis of the mitochondrial genome of Prunella vulgaris

Prunella vulgaris (Lamiaceae) is widely distributed in Eurasia. Former studies have demonstrated that P. vulgaris has a wide range of pharmacological effects. Nevertheless, no complete P. vulgaris mitochondrial genome has been reported, which limits further understanding of the biology of P. vulgaris. Here, we assembled the first complete mitochondrial genome of P. vulgaris using a hybrid assembly strategy based on sequencing data from both Nanopore and Illumina platforms. Then, the mitochondrial genome of P. vulgaris was analyzed comprehensively in terms of gene content, codon preference, intercellular gene transfer, phylogeny, and RNA editing. The mitochondrial genome of P. vulgaris has two circular structures. It has a total length of 297, 777 bp, a GC content of 43.92%, and 29 unique protein-coding genes (PCGs). There are 76 simple sequence repeats (SSRs) in the mitochondrial genome, of which tetrameric accounts for a large percentage (43.4%). A comparative analysis between the mitochondrial and chloroplast genomes revealed that 36 homologous fragments exist in them, with a total length of 28, 895 bp. The phylogenetic analysis showed that P. vulgaris belongs to the Lamiales family Lamiaceae and P. vulgaris is closely related to Salvia miltiorrhiza. In addition, the mitochondrial genome sequences of seven species of Lamiaceae are unconservative in their alignments and undergo frequent genome reorganization. This work reports for the first time the complete mitochondrial genome of P. vulgaris, which provides useful genetic information for further Prunella studies.

Beyond a PPR-RNA recognition code: Many aspects matter for the multi-targeting properties of RNA editing factor PPR56

The mitochondrial C-to-U RNA editing factor PPR56 of the moss Physcomitrium patens is an RNA-binding pentatricopeptide repeat protein equipped with a terminal DYW-type cytidine deaminase domain. Transferred into Escherichia coli, PPR56 works faithfully on its two native RNA editing targets, nad3eU230SL and nad4eU272SL, and also converts cytidines into uridines at over 100 off-targets in the bacterial transcriptome. Accordingly, PPR56 is attractive for detailed mechanistic studies in the heterologous bacterial setup, allowing for scoring differential RNA editing activities of many target and protein variants in reasonable time. Here, we report (i) on the effects of numerous individual and combined PPR56 protein and target modifications, (ii) on the spectrum of off-target C-to-U editing in the bacterial background transcriptome for PPR56 and two variants engineered for target re-direction and (iii) on combinations of targets in tandem or separately at the 5'- and 3'-ends of large mRNAs. The latter experimentation finds enhancement of RNA editing at weak targets in many cases, including cox3eU290SF as a new candidate mitogenome target. We conclude that C-to-U RNA editing can be much enhanced by transcript features also outside the region ultimately targeted by PPRs of a plant editing factor, possibly facilitated by its enrichment or scanning along transcripts.

Comparative analysis of mitochondrial genomes of Schisandra repanda and Kadsura japonica

The family Schisandraceae is a basal angiosperm plant group distributed in East and Southeast Asia and includes many medicinal plant species such as Schisandra chinensis. In this study, mitochondrial genomes (mitogenomes) of two species, Schisandra repanda and Kadsura japonica, in the family were characterized through de novo assembly using sequencing data obtained with Oxford Nanopore and Illumina sequencing technologies. The mitogenomes of S. repanda were assembled into one circular contig (571,107 bp) and four linear contigs (10,898-607,430 bp), with a total of 60 genes: 38 protein-coding genes (PCGs), 19 tRNA genes, and 3 rRNA genes. The mitogenomes of K. japonica were assembled into five circular contigs (211,474-973,503 bp) and three linear contigs (8,010-72,712 bp), with a total of 66 genes: 44 PCGs, 19 tRNA genes, and 3 rRNA genes. The mitogenomes of the two species had complex structural features with high repeat numbers and chloroplast-derived sequences, as observed in other plant mitogenomes. Phylogenetic analysis based on PCGs revealed the taxonomical relationships of S. repanda and K. japonica with other species from Schisandraceae. Finally, molecular markers were developed to distinguish between S. repanda, K. japonica, and S. chinensis on the basis of InDel polymorphisms present in the mitogenomes. The mitogenomes of S. repanda and K. japonica will be valuable resources for molecular and taxonomic studies of plant species that belong to the family Schisandraceae.

Assembly and analysis of the first complete mitochondrial genome of Punica granatum and the gene transfer from chloroplast genome

Pomegranate (Punica granatum L.) is one of the oldest fruits with edible, medicinal and ornamental values. However, there is no report on the mitochondrial genome of pomegranate. In this study, the mitochondrial genome of P. granatum was sequenced, assembled and analyzed in detail, while the chloroplast genome was assembled using the same set of data. The results showed that the P. granatum mitogenome had a multi branched structure, using BGI + Nanopore mixed assembly strategy. The total genome length was 404,807 bp, with the GC content of 46.09%, and there were 37 protein coding genes, 20 tRNA genes and three rRNA genes. In the whole genome, 146 SSRs were identified. Besides, 400 pairs of dispersed repeats were detected, including 179 palindromic, 220 forward and one reverse. In the P. granatum mitochondrial genome, 14 homologous fragments of chloroplast genome were found, accounting for 0.54% of the total length. Phylogenetic analysis showed that among the published mitochondrial genomes of related genera, P. granatum had the closest genetic relationship with Lagerstroemia indica of Lythraceae. The 580 and 432 RNA editing sites were predicted on 37 protein coding genes of mitochondrial genome using BEDTools software and online website PREPACT respectively, but all were from C to U, of which ccmB and nad4 gene were most frequently edited, with 47 sites. This study provides a theoretical basis for understanding the evolution of higher plants, species classification and identification, and will also be useful for further utilization of pomegranate germplasm resources.

Mitogenomic Research of Silverleaf Sunflower (Helianthus argophyllus) and Its Interspecific Hybrids

Interspecific hybridization is widespread for sunflowers, both in wild populations and commercial breeding. One of the most common species that can efficiently cross with Helianthus annuus is the Silverleaf sunflower-Helianthus argophyllus. The current study carried out structural and functional organization analyses of mitochondrial DNA in H. argophyllus and the interspecific hybrid, H. annuus (VIR114A line) × H. argophyllus. The complete mitogenome of H. argophyllus counts 300,843 bp, has a similar organization to the mitogenome of cultivated sunflower, and holds SNPs typical for wild sunflowers. RNA editing analysis predicted 484 sites in H. argophyllus mitochondrial CDS. The mitochondrial genome of the H. annuus × H. argophyllus hybrid is identical to the maternal line (VIR114A). We expected that significant rearrangements in the mitochondrial DNA of the hybrid would occur, due to the frequent recombination. However, the hybrid mitogenome lacks rearrangements, presumably due to the preservation of nuclear-cytoplasmic interaction paths.

Assembly and Comparative Analysis of the Complete Mitochondrial Genome of Two Species of Calla Lilies (Zantedeschia, Araceae)

In this study, the mitochondrial genomes of two calla species, Zantedeschia aethiopica Spreng. and Zantedeschia odorata Perry., were assembled and compared for the first time. The Z. aethiopica mt genome was assembled into a single circular chromosome, measuring 675,575 bp in length with a 45.85% GC content. In contrast, the Z. odorata mt genome consisted of bicyclic chromosomes (chromosomes 1 and 2), measuring 719,764 bp and exhibiting a 45.79% GC content. Both mitogenomes harbored similar gene compositions, with 56 and 58 genes identified in Z. aethiopica and Z. odorata, respectively. Analyses of codon usage, sequence repeats, gene migration from chloroplast to mitochondrial, and RNA editing were conducted for both Z. aethiopica and Z. odorata mt genomes. Phylogenetic examination based on the mt genomes of these two species and 30 other taxa provided insights into their evolutionary relationships. Additionally, the core genes in the gynoecium, stamens, and mature pollen grains of the Z. aethiopica mt genome were investigated, which revealed maternal mitochondrial inheritance in this species. In summary, this study offers valuable genomic resources for future research on mitogenome evolution and the molecular breeding of calla lily.

Progress, challenge and prospect of plant plastome annotation

The plastome (plastid genome) represents an indispensable molecular data source for studying phylogeny and evolution in plants. Although the plastome size is much smaller than that of nuclear genome, and multiple plastome annotation tools have been specifically developed, accurate annotation of plastomes is still a challenging task. Different plastome annotation tools apply different principles and workflows, and annotation errors frequently occur in published plastomes and those issued in GenBank. It is therefore timely to compare available annotation tools and establish standards for plastome annotation. In this review, we review the basic characteristics of plastomes, trends in the publication of new plastomes, the annotation principles and application of major plastome annotation tools, and common errors in plastome annotation. We propose possible methods to judge pseudogenes and RNA-editing genes, jointly consider sequence similarity, customed algorithms, conserved domain or protein structure. We also propose the necessity of establishing a database of reference plastomes with standardized annotations, and put forward a set of quantitative standards for evaluating plastome annotation quality for the scientific community. In addition, we discuss how to generate standardized GenBank annotation flatfiles for submission and downstream analysis. Finally, we prospect future technologies for plastome annotation integrating plastome annotation approaches with diverse evidences and algorithms of nuclear genome annotation tools. This review will help researchers more efficiently use available tools to achieve high-quality plastome annotation, and promote the process of standardized annotation of the plastome.

Analysis of the Complete Mitochondrial Genome of the Bitter Gourd (Momordica charantia)

Bitter gourd (Momordica charantia L.) is a significant vegetable. Although it has a special bitter taste, it is still popular with the public. The industrialization of bitter gourd could be hampered by a lack of genetic resources. The bitter gourd's mitochondrial and chloroplast genomes have not been extensively studied. In the present study, the mitochondrial genome of bitter gourd was sequenced and assembled, and its substructure was investigated. The mitochondrial genome of bitter gourd is 331,440 bp with 24 unique core genes, 16 variable genes, 3 rRNAs, and 23 tRNAs. We identified 134 SSRs and 15 tandem repeats in the entire mitochondrial genome of bitter gourd. Moreover, 402 pairs of repeats with a length greater than or equal to 30 were observed in total. The longest palindromic repeat was 523 bp, and the longest forward repeat was 342 bp. We found 20 homologous DNA fragments in bitter gourd, and the summary insert length was 19,427 bp, accounting for 5.86% of the mitochondrial genome. We predicted a total of 447 potential RNA editing sites in 39 unique PCGs and also discovered that the ccmFN gene has been edited the most often, at 38 times. This study provides a basis for a better understanding and analysis of differences in the evolution and inheritance patterns of cucurbit mitochondrial genomes.

Complete chloroplast genome data of Shorea macrophylla (Engkabang): Structural features, comparative and phylogenetic analysis

Shorea macrophylla belongs to the Shorea genus under the Dipterocarpaceae family. It is a woody tree that grows in the rainforest in Southeast Asia. The complete chloroplast (cp) genome sequence of S. macrophylla is reported here. The genomic size of S. macrophylla is 150,778 bp and it possesses a circular structure with conserved constitute regions of large single copy (LSC, 83,681 bp) and small single copy (SSC, 19,813 bp) regions, as well as a pair of inverted repeats with a length of 23,642 bp. It has 112 unique genes, including 78 protein-coding genes, 30 tRNA genes, and four rRNA genes. The genome exhibits a similar GC content, gene order, structure, and codon usage when compared to previously reported chloroplast genomes from other plant species. The chloroplast genome of S. macrophylla contained 262 SSRs, the most prevalent of which was A/T, followed by AAT/ATT. Furthermore, the sequences contain 43 long repeat sequences, practically most of them are forward or palindrome type long repeats. The genome structure of S. macrophylla was compared to the genomic structures of closely related species from the same family, and eight mutational hotspots were discovered. The phylogenetic analysis demonstrated a close relationship between Shorea and Parashorea species, indicating that Shorea is not monophyletic. The complete chloroplast genome sequence analysis of S. macrophylla reported in this paper will contribute to further studies in molecular identification, genetic diversity, and phylogenetic research.

Rickettsial DNA and a trans-splicing rRNA group I intron in the unorthodox mitogenome of the fern Haplopteris ensiformis

Plant mitochondrial genomes can be complex owing to highly recombinant structures, lack of gene syntenies, heavy RNA editing and invasion of chloroplast, nuclear or even foreign DNA by horizontal gene transfer (HGT). Leptosporangiate ferns remained the last major plant clade without an assembled mitogenome, likely owing to a demanding combination of the above. We here present both organelle genomes now for Haplopteris ensiformis. More than 1,400 events of C-to-U RNA editing and over 500 events of reverse U-to-C edits affect its organelle transcriptomes. The Haplopteris mtDNA is gene-rich, lacking only the ccm gene suite present in ancestral land plant mitogenomes, but is highly unorthodox, indicating extraordinary recombinogenic activity. Although eleven group II introns known in disrupted trans-splicing states in seed plants exist in conventional cis-arrangements, a particularly complex structure is found for the mitochondrial rrnL gene, which is split into two parts needing reassembly on RNA level by a trans-splicing group I intron. Aside from ca. 80 chloroplast DNA inserts that complicated the mitogenome assembly, the Haplopteris mtDNA features as an idiosyncrasy 30 variably degenerated protein coding regions from Rickettiales bacteria indicative of heavy bacterial HGT on top of tRNA genes of chlamydial origin.

De Novo Assembly and Comparative Analysis of the Complete Mitochondrial Genome of Chaenomeles speciosa (Sweet) Nakai Revealed the Existence of Two Structural Isomers

As a valuable Chinese traditional medicinal species, Chaenomeles speciosa (Sweet) Nakai (C. speciosa) is a natural resource with significant economic and ornamental value. However, its genetic information is not well understood. In this study, the complete mitochondrial genome of C. speciosa was assembled and characterized to explore the repeat sequences, recombination events, rearrangements, and IGT, to predict RNA editing sites, and to clarify the phylogenetic and evolutionary relationship. The C. speciosa mitochondrial genome was found to have two circular chromosomes as its major conformation, with a total length of 436,464 bp and 45.2% GC content. The mitochondrial genome contained 54 genes, including 33 unique protein-coding genes, 18 tRNAs, and 3 rRNA genes. Seven pairs of repeat sequences involving recombination events were analyzed. Both the repeat pairs, R1 and R2, played significant roles in mediating the major and minor conformations. In total, 18 MTPTs were identified, 6 of which were complete tRNA genes. There were 454 RNA editing sites in the 33 protein-coding sequences predicted by the PREPACT3 program. A phylogenetic analysis based on 22 species of mitochondrial genomes was constructed and indicated highly conserved PCG sequences. Synteny analyses showed extensive genomic rearrangements in the mitochondrial genome of C. speciosa and closely related species. This work is the first to report the C. speciosa mitochondrial genome, which is of great significance for conducting additional genetic studies on this organism.

Multichromosomal Mitochondrial Genome of Paphiopedilum micranthum: Compact and Fragmented Genome, and Rampant Intracellular Gene Transfer

Orchidaceae is one of the largest families of angiosperms. Considering the large number of species in this family and its symbiotic relationship with fungi, Orchidaceae provide an ideal model to study the evolution of plant mitogenomes. However, to date, there is only one draft mitochondrial genome of this family available. Here, we present a fully assembled and annotated sequence of the mitochondrial genome (mitogenome) of Paphiopedilum micranthum, a species with high economic and ornamental value. The mitogenome of P. micranthum was 447,368 bp in length and comprised 26 circular subgenomes ranging in size from 5973 bp to 32,281 bp. The genome encoded for 39 mitochondrial-origin, protein-coding genes; 16 tRNAs (three of plastome origin); three rRNAs; and 16 ORFs, while rpl10 and sdh3 were lost from the mitogenome. Moreover, interorganellar DNA transfer was identified in 14 of the 26 chromosomes. These plastid-derived DNA fragments represented 28.32% (46,273 bp) of the P. micranthum plastome, including 12 intact plastome origin genes. Remarkably, the mitogenome of P. micranthum and Gastrodia elata shared 18% (about 81 kb) of their mitochondrial DNA sequences. Additionally, we found a positive correlation between repeat length and recombination frequency. The mitogenome of P. micranthum had more compact and fragmented chromosomes compared to other species with multichromosomal structures. We suggest that repeat-mediated homologous recombination enables the dynamic structure of mitochondrial genomes in Orchidaceae.

SlRIP1b is a global organellar RNA editing factor, required for normal fruit development in tomato plants

RNA editing in plant organelles involves numerous C-U conversions, which often restore evolutionarily conserved codons and may generate new translation initiation and termination codons. These RNA maturation events rely on a subset of nuclear-encoded protein cofactors. Here, we provide evidence of the role of SlRIP1b on RNA editing of mitochondrial transcripts in tomato (Solanum lycopersicum) plants. SlRIP1b is a RIP/MORF protein that was originally identified as an interacting partner of the organellar editing factor SlORRM4. Mutants of SlRIP1b, obtained by CRISPR/Cas9 strategy, exhibited abnormal carpel development and grew into fruit with more locules. RNA-sequencing revealed that SlRIP1b affects the C-U editing of numerous mitochondrial pre-RNA transcripts and in particular altered RNA editing of various cytochrome c maturation (CCM)-related genes. The slrip1b mutants display increased H₂ O₂ and aberrant mitochondrial morphologies, which are associated with defects in cytochrome c biosynthesis and assembly of respiratory complex III. Taken together, our results indicate that SlRIP1b is a global editing factor that plays a key role in CCM and oxidative phosphorylation system biogenesis during fruit development in tomato plants. These data provide important insights into the molecular roles of organellar RNA editing factors during fruit development.

Revisiting chloroplast genomic landscape and annotation towards comparative chloroplast genomes of Rhamnaceae

BACKGROUND

Massive parallel sequencing technologies have enabled the elucidation of plant phylogenetic relationships from chloroplast genomes at a high pace. These include members of the family Rhamnaceae. The current Rhamnaceae phylogenetic tree is from 13 out of 24 Rhamnaceae chloroplast genomes, and only one chloroplast genome of the genus Ventilago is available. Hence, the phylogenetic relationships in Rhamnaceae remain incomplete, and more representative species are needed.

RESULTS

The complete chloroplast genome of Ventilago harmandiana Pierre was outlined using a hybrid assembly of long- and short-read technologies. The accuracy and validity of the final genome were confirmed with PCR amplifications and investigation of coverage depth. Sanger sequencing was used to correct for differences in lengths and nucleotide bases between inverted repeats because of the homopolymers. The phylogenetic trees reconstructed using prevalent methods for phylogenetic inference were topologically similar. The clustering based on codon usage was congruent with the molecular phylogenetic tree. The groups of genera in each tribe were in accordance with tribal classification based on molecular markers. We resolved the phylogenetic relationships among six Hovenia species, three Rhamnus species, and two Ventilago species. Our reconstructed tree provides the most complete and reliable low-level taxonomy to date for the family Rhamnaceae. Similar to other higher plants, the RNA editing mostly resulted in converting serine to leucine. Besides, most genes were subjected to purifying selection. Annotation anomalies, including indel calling errors, unaligned open reading frames of the same gene, inconsistent prediction of intergenic regions, and misannotated genes, were identified in the published chloroplast genomes used in this study. These could be a result of the usual imperfections in computational tools, and/or existing errors in reference genomes. Importantly, these are points of concern with regards to utilizing published chloroplast genomes for comparative genomic analysis.

CONCLUSIONS

In summary, we successfully demonstrated the use of comprehensive genomic data, including DNA and amino acid sequences, to build a reliable and high-resolution phylogenetic tree for the family Rhamnaceae. Additionally, our study indicates that the revision of genome annotation before comparative genomic analyses is necessary to prevent the propagation of errors and complications in downstream analysis and interpretation.

RNA Editing in Chloroplast: Advancements and Opportunities

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

Plant mitochondrial RNA editing factors can perform targeted C-to-U editing of nuclear transcripts in human cells

RNA editing processes are strikingly different in animals and plants. Up to thousands of specific cytidines are converted into uridines in plant chloroplasts and mitochondria whereas up to millions of adenosines are converted into inosines in animal nucleo-cytosolic RNAs. It is unknown whether these two different RNA editing machineries are mutually incompatible. RNA-binding pentatricopeptide repeat (PPR) proteins are the key factors of plant organelle cytidine-to-uridine RNA editing. The complete absence of PPR mediated editing of cytosolic RNAs might be due to a yet unknown barrier that prevents its activity in the cytosol. Here, we transferred two plant mitochondrial PPR-type editing factors into human cell lines to explore whether they could operate in the nucleo-cytosolic environment. PPR56 and PPR65 not only faithfully edited their native, co-transcribed targets but also different sets of off-targets in the human background transcriptome. More than 900 of such off-targets with editing efficiencies up to 91%, largely explained by known PPR-RNA binding properties, were identified for PPR56. Engineering two crucial amino acid positions in its PPR array led to predictable shifts in target recognition. We conclude that plant PPR editing factors can operate in the entirely different genetic environment of the human nucleo-cytosol and can be intentionally re-engineered towards new targets.

DYW deaminase domain has a distinct preference for neighboring nucleotides of the target RNA editing sites

C-to-U RNA editing sites in plant organelles show a strong bias for neighboring nucleotides. The nucleotide upstream of the target cytidine is typically C or U, whereas A and G are less common and rare, respectively. In pentatricopeptide repeat (PPR)-type RNA editing factors, the PPR domain specifically binds to the 5' sequence of target cytidines, whereas the DYW domain catalyzes the C-to-U deamination. We comprehensively analyzed the effects of neighboring nucleotides of the target cytidines using an Escherichia coli orthogonal system. Physcomitrium PPR56 efficiently edited target cytidines when the nucleotide upstream was U or C, whereas it barely edited when the position was G or the nucleotide downstream was C. This preference pattern, which corresponds well with the observed nucleotide bias for neighboring nucleotides in plant organelles, was altered when the DYW domain of OTP86 or DYW1 was adopted. The PPR56 chimeric proteins edited the target sites even when the -1 position was G. Our results suggest that the DYW domain possesses a distinct preference for the neighboring nucleotides of the target sites, thus contributing to target selection in addition to the existing selection determined by the PPR domain.

The evolution of extremely diverged plastomes in Selaginellaceae (lycophyte) is driven by repeat patterns and the underlying DNA maintenance machinery

Two factors are proposed to account for the unusual features of organellar genomes: the disruptions of organelle-targeted DNA replication, repair, and recombination (DNA-RRR) systems in the nuclear genome and repetitive elements in organellar genomes. Little is known about how these factors affect organellar genome evolution. The deep-branching vascular plant family Selaginellaceae is known to have a deficient DNA-RRR system and convergently evolved organellar genomes. However, we found that the plastid genome (plastome) of Selaginella sinensis has extremely accelerated substitution rates, a low GC content, pervasive repeat elements, a dynamic network structure, and it lacks direct or inverted repeats. Unexpectedly, its organelle DNA-RRR system is short of a plastid-targeted Recombinase A1 (RecA1) and a mitochondrion-targeted RecA3, in line with other explored Selaginella species. The plastome contains a large collection of short- and medium-sized repeats. Given the absence of RecA1 surveillance, we propose that these repeats trigger illegitimate recombination, accelerated mutation rates, and structural instability. The correlations between repeat quantity and architectural complexity in the Selaginella plastomes support these conclusions. We, therefore, hypothesize that the interplay of the deficient DNA-RRR system and the high repeat content has led to the extraordinary divergence of the S. sinensis plastome. Our study not only sheds new light on the mechanism of plastome divergence by emphasizing the power of cytonuclear integration, but it also reconciles the longstanding contradiction on the effects of DNA-RRR system disruption on genome structure evolution.

Mitochondrial phylogenomics of Cuscuta (Convolvulaceae) reveals a potentially functional horizontal gene transfer from the host

Horizontal gene transfers (HGTs) from host or other organisms have been reported in mitochondrial genomes of parasitic plants. Genes transferred in this fashion have usually been found non-functional. Several examples of HGT from the mitochondrial genome of parasitic Cuscuta (Convolvulaceae) to its hosts have been reported, but not vice versa. Here we used 31 protein-coding mitochondrial genes to infer the phylogeny of Cuscuta, and compared it with previous nuclear and plastid phylogenetic estimates. We also investigated the presence of HGTs within these lineages. Unlike in plastid genomes, we did not find extensive gene loss in their mitochondrial counterparts. Our results reveal the first example of organellar HGT from host to Cuscuta. Mitochondrial atp1 genes of South African subgenus Pachystigma were inferred to be transferred from Lamiales, with high support. Moreover, the horizontally transferred atp1 gene has functionally replaced the native, vertically transmitted copy, has an intact open reading frame, and is under strong purifying selection, all of which suggests that this xenolog remains functional. The mitochondrial phylogeny of Cuscuta is generally consistent with previous plastid and nuclear phylogenies, except for the misplacement of Pachystigma when atp1 is included. This incongruence may be caused by the HGT mentioned above. No example of HGT was found within mitochondrial genes of three other, independently evolved parasitic lineages we sampled: Cassytha/Laurales, Krameria/Zygophyllales, and Lennooideae/Boraginales.

Plastid phylogenomic analyses of the Selaginella sanguinolenta group (Selaginellaceae) reveal conflict signatures resulting from sequence types, outlier genes, and pervasive RNA editing

Different from the generally conserved plastomes (plastid genomes) of most land plants, the Selaginellaceae plastomes exhibit dynamic structure, high GC content and high substitution rates. Previous plastome analyses identified strong conflict on several clades in Selaginella, however the factors causing the conflictions and the impact on the phylogenetic inference have not been sufficiently investigated. Here, we dissect the distribution of phylogenetic signals and conflicts in Selaginella sanguinolenta group, the plastome of which is DR (direct repeats) structure and with genome-wide RNA editing. We analyzed the data sets including 22 plastomes representing all species of the S. sanguinolenta group, covering the entire geographical distribution from the Himalayas to Siberia and the Russian Far East regions. We recovered four different topologies by applying multispecies coalescent (ASTRAL) and concatenation methods (IQ-TREE and RAxML) on four data sets of PC (protein-coding genes), NC (non-coding sequences), PCN (the concatenated PC and NC), and RC (predicted RNA editing sites "C" were corrected by "T"), respectively. Six monophyletic clades, S. nummularifolia clade, S. rossii clade, S. sajanensis clade, S. sanguinolenta I clade, S. sanguinolenta II clade, and S. sanguinolenta III clade, were consistently resolved and supported by the characteristics of GC content, RNA editing frequency, and gene content. However, the relationships among these clades varied across the four topologies. To explore the underlying causes of the uncertainty, we compared the phylogenetic signals of the four topologies. We identified that the sequence types (coding versus non-coding), outlier genes (genes with extremely high |ΔGLS| values), and C-to-U RNA editing frequency in the protein-coding genes were responsible for the unstable phylogenomic relationship. We further revealed a significant positive correlation between the |ΔGLS| values and the variation coefficient of the RNA editing number. Our results demonstrated that the coalescent method performed better than the concatenation method in overcoming the problems caused by outlier genes and extreme RNA editing events. Our study particularly focused on the importance of exploring the plastid phylogenomic conflicts and suggested conducting concatenated analyses cautiously when adopting organelle genome data.

Breaking the limits - multichromosomal structure of an early eudicot Pulsatilla patens mitogenome reveals extensive RNA-editing, longest repeats and chloroplast derived regions among sequenced land plant mitogenomes

BACKGROUND

The mitogenomes of vascular plants are one of the most structurally diverse molecules. In the present study we characterize mitogenomes of a rare and endangered species Pulsatilla patens. We investigated the gene content and its RNA editing potential, repeats distribution and plastid derived sequences.

RESULTS

The mitogenome structure of early divergent eudicot, endangered Pulsatilla patens does not support the master chromosome hypothesis, revealing the presence of three linear chromosomes of total length 986 613 bp. The molecules are shaped by the presence of extremely long, exceeding 87 kbp repeats and multiple chloroplast-derived regions including nearly complete inverted repeat. Since the plastid IR content of Ranunculales is very characteristic, the incorporation into mitogenome could be explained rather by intracellular transfer than mitochondrial HGT. The mitogenome contains almost a complete set of genes known from other vascular plants with exception of rps10 and sdh3, the latter being present but pseudogenized. Analysis of long ORFs enabled the identification of genes which are rarely present in plant mitogenomes, including RNA and DNA polymerases, albeit their presence even at species level is variable. Mitochondrial transcripts of P. patens were edited with a high frequency, which exceeded the level known in other analyzed angiosperms, despite the strict qualification criteria of counting the editing events and taking into analysis generally less frequently edited leaf transcriptome. The total number of edited sites was 902 and nad4 was identified as the most heavily edited gene with 65 C to U changes. Non-canonical, reverse U to C editing was not detected. Comparative analysis of mitochondrial genes of three Pulsatilla species revealed a level of variation comparable to chloroplast CDS dataset and much higher infrageneric differentiation than in other known angiosperm genera. The variation found in CDS of mitochondrial genes is comparable to values found among Pulsatilla plastomes. Despite the complicated mitogenome structure, 14 single copy regions of 329 kbp, not splitted by repeats or plastid-derived sequences (MTPT), revealed the potential for phylogenetic, phylogeographic and population genetics studies by revealing intra- and interspecific collinearity.

CONCLUSIONS

This study provides valuable new information about mitochondrial genome of early divergent eudicots, Pulsatilla patens, revealed multi-chromosomal structure and shed new light on mitogenomics of early eudicots.