Repeats of Unusual Size in Plant Mitochondrial Genomes: Identification, Incidence and Evolution

Assembly and Comparative Analysis of the Complete Mitochondrial Genome of Hippophae salicifolia

This study reports the first complete mitochondrial genome assembly of Hippophae salicifolia, an ecologically and economically important plant endemic to the Himalayas. The 475,105 bp genome has a 44.80% GC content and an overall AT bias, comprising 74 genes (37 protein-coding, 31 tRNA, three rRNA, and three pseudogenes). We identified extensive repetitive elements, including 188 SSRs, 20 tandem repeats, and 455 dispersed repeats, and explored their potential roles in genome evolution. Codon usage analysis showed a bias for codons ending in A or U, while RNA editing analysis revealed 415 sites that mostly convert hydrophilic to hydrophobic amino acids. Phylogenetic and collinearity analyses clarified evolutionary relationships within Hippophae and uncovered genome rearrangements. In addition, extensive gene transfer was detected between the mitochondrial and chloroplast genomes. Ka/Ks and nucleotide diversity analyses indicate that most genes are under purifying selection, with some possibly undergoing positive selection. Overall, these findings enhance our understanding of the structural and evolutionary features of the H. salicifolia mitochondrial genome and provide valuable insights for the genetic improvement and conservation of Hippophae species.

Assembly and analysis of the complete mitochondrial genome of Leonurus japonicus (Lamiaceae)

Leonurus japonicus Houtt. (L. japonicus), as an important plant resource with both ornamental and medicinal value, has now spread worldwide and is widely studied. Currently, its chromosomal genome and chloroplast genome (cpDNA) have been reported, but the mitochondrial genome (mtDNA) has not yet been explored. In this study, we extracted DNA from fresh leaves of L. japonicus and performed sequencing and assembly of its mtDNA using both second-generation and third-generation sequencing technologies. The complete mtDNA of L. japonicus is 382,905 bp in length, with a GC content of 45.13%. This genome includes 15 tRNA genes, 32 protein-coding genes (PCGs), and 4 rRNA genes. In this mtDNA genome, we predicted a total of 480 RNA editing sites among the 32 PCGs. Subsequently, we conducted analyses on repetitive sequences, organelle genome sequence migration, and Relative Synonymous Codon Usage (RSCU). There are 28 homologous sequence fragments between the mtDNA and cpDNA of L. japonicus, which are related to the migration of 10 mtDNA genes. The RSCU analysis predicted 28 high-frequency codons, most of which prefer to end with A/U. Selection pressure analysis indicated that the Ka/Ks ratio for the majority of PCGs is less than 1, suggesting they are highly conserved during evolutionary processes. Phylogenetic results from 24 species indicate that the genera Leonurus and Scutellaria within the Lamiaceae family have the closest relationships. In summary, we have successfully assembled the complete mtDNA of L. japonicus by integrating second-generation and third-generation sequencing data for the first time. Subsequent multi-faceted analyses have allowed us to gain deeper insights into the numerous features of this genome, providing important reference data for the molecular genetics, dynamic evolution, and species identification of this plant. This work promotes the conservation and development of this important resource of medicinal and edible plants.

Unveiling the mitochondrial genome of Salvia splendens insights into the evolutionary traits within the genus Salvia

Previously, we resolved the complete sequences of the mitochondrial genomes (mitogenome) of two Salvia species (S. miltiorrhiza and S. officinalis). The major configurations of these two species were two circular chromosomes. In this study, we further studied the mitogenome of a floral species of Salvia (Salvia splendens) to understand the diversity and evolution of the Salvia mitogenomes. We sequenced the total DNAs of S. splendens using the Nanopore and Illumina platforms and assembled the mitogenome using a hybrid assembly strategy. The major configurations of the S. splendens were two circular chromosomes with lengths of 182,239 and 165,055 bp. There were 32 protein-coding genes (PCGs), three rRNA genes, and 18 tRNA genes annotated in the S. splendens mitogenome. We found 56 pairs of repetitive sequences in the S. splendens mitogenome. Three of them (R01, 04, and 07) could mediate recombination, whose products could be identified by the mapping of Nanopore reads, PCR amplifications, and Sanger sequencing of the PCR products. 457 RNA editing sites were identified in the S. splendens mitochondrial RNAs when comparing the RNA-seq data with their corresponding DNA templates. We showed that S. splendens was a sister taxon to S. miltiorrhiza based on the mitogenomes, consistent with the phylogeny determined with the plastome sequences. Crucially, we developed 12 mitochondrial markers sourced from mitochondrial intron regions to facilitate the identification of three Salvia species. Our study offers a comprehensive view of the structure of the Salvia mitogenomes and provides robust mitochondrial markers for Salvia species identification.

Decoding the mitogenome of rosemary (Salvia rosmarinus): insights into genome evolution, structural dynamics and prospects for mitochondrial engineering

BACKGROUND

Rosemary (Salvia rosmarinus), an aromatic evergreen shrub of the Salvia (Lamiaceae), is native to the Mediterranean region, thriving in rocky or arid soils. Widely used in food, pharmaceuticals, and cosmetics, its clonal reproduction poses significant challenges for breeding and germplasm innovation. While mitogenome engineering holds promise for introducing heritable mutations, incomplete mitogenome information for rosemary has hindered such efforts. This study addresses this gap by assembling and analyzing the complete mitogenome of S. rosmarinus, focusing on its structure, repetitive sequences, RNA editing events, intracellular gene transfer (IGT), and phylogenetic relationships.

RESULTS

The S. rosmarinus mitogenome spans 384,113 bp with a GC content of 44.8%, containing 34 unique protein-coding genes and 114 simple sequence repeats. Comparative analysis revealed 28 homologous segments shared between the mitogenome and plastome, totaling 18,675 bp in length. Furthermore, homologous fragments between nuclear and organellar genomes were identified, including 1,069,255 bp of organelle-derived sequences in the nuclear genome, with 194,689 bp from nuclear plastid DNA transfers (NUPTs) and 15,192 bp from nuclear mitochondrial DNA transfers (NUMTs). NUPTs were more abundant and contributed more significantly to the total length. Synteny analysis of eight Lamiales species revealed extensive mitogenomic recombination and structural rearrangements. These findings highlight the dynamic nature of mitogenomes, offering insights into genome evolution and supporting future breeding programs to enhance the genetic diversity and adaptability of S. rosmarinus.

CONCLUSIONS

This study provides the first complete mitogenome of S. rosmarinus, revealing dispersed repeats, RNA editing, and horizontal gene transfer between the nuclear and organelle genomes. The mitogenome exhibits a typical circular structure with evidence of frequent recombination, providing valuable insights into Salvia mitochondrial genetics, genome evolution, and molecular biology.

Disruption of recombination machinery alters the mutational landscape in plant organellar genomes

Land plant organellar genomes have extremely low rates of point mutation yet also experience high rates of recombination and genome instability. Characterizing the molecular machinery responsible for these patterns is critical for understanding the evolution of these genomes. While much progress has been made toward understanding recombination activity in land plant organellar genomes, the relationship between recombination pathways and point mutation rates remains uncertain. The organellar-targeted mutS homolog MSH1 has previously been shown to suppress point mutations as well as non-allelic recombination between short repeats in Arabidopsis thaliana. We therefore implemented high-fidelity Duplex Sequencing to test if other genes that function in recombination and maintenance of genome stability also affect point mutation rates. We found small to moderate increases in the frequency of single nucleotide variants (SNVs) and indels in mitochondrial and/or plastid genomes of A. thaliana mutant lines lacking radA, recA1, or recA3. In contrast, osb2 and why2 mutants did not exhibit an increase in point mutations compared to wild-type (WT) controls. In addition, we analyzed the distribution of SNVs in previously generated Duplex Sequencing data from A. thaliana organellar genomes and found unexpected strand asymmetries and large effects of flanking nucleotides on mutation rates in WT plants and msh1 mutants. Finally, using long-read Oxford Nanopore sequencing, we characterized structural variants in organellar genomes of the mutant lines and show that different short repeat sequences become recombinationally active in different mutant backgrounds. Together, these complementary sequencing approaches shed light on how recombination may impact the extraordinarily low point mutation rates in plant organellar genomes.

Complete analysis and phylogenetic analysis of Polygonatum sibiricum mitochondria

In this project, we studied the complete mitogenome of the liliaceae medicinal plant Polygonatum sibiricum. The genome is represented by a circular ring molecule with a length of 691,910 bp and a GC content of 46.33%. Mitochondrial genome composition is slightly biased towards A+T, with AT accounting for 53.67%, and AT skewness slightly positive (0.092%). The complete mitogenome has a total of sixty-three unique genes, including thirty-nine protein-coding genes, twenty-one transfer RNAs (tRNAs) and three ribosomal RNAs (rRNAs). We examined codon use, repeat sequence, RNA editing in the mitogenome of P. sibiricum, and elucidated species classification based on phylogenetic trees of mitogenome of twenty-three species. Our results provide comprehensive information on the mitogenome of P. sibiricum and show for the first time the evolutionary relationship between the mitogenome of P. sibiricum and Chlorophytum comosum in the Asparagales family.

Comparative analysis of the complete mitochondrial genomes of Firmiana danxiaensis and F. kwangsiensis (Malvaceae), two endangered Firmiana species in China

MAIN CONCLUSION

We reported the mitogenomes of F. danxiaensis and F. kwangsiensis for the first time. Mitogenome size and structure differ significantly between them. Firmiana danxiaensis and F. kwangsiensis belong to the Firmiana genus and are distributed in the Danxia and Karst regions of southern China. Both species have been designated as endangered. Currently, the chloroplast genomes of F. danxiaensis and F. kwangsiensis have been sequenced, but the mitochondrial genome (mitogenome) of these two species has not been reported. To further understand the mitogenome characteristics, evolution, and phylogeny of F. danxiaensis and F. kwangsiensis, we assembled the mitogenomes of these two species based on a combination of Illumina and Nanopore sequencing methods. The mitogenome of F. danxiaensis exhibits a branching structure consisting of nine circular molecules with a total length of 938,890 bp, while the F. kwangsiensis has a circular structure with a length of 736,334 bp. Compared to F. kwangsiensis, F. danxiaensis has more tRNA genes, SSRs, tandem repeats, and dispersed repeats, while the codon use patterns are similar in these two species. There were 24 and 23 homologous sequences between mitogenome and chloroplast genome of F. danxiaensis and F. kwangsiensis, accounting for 0.37% and 0.49% of the mitogenome, respectively. In addition, the Ka/Ks ratio and the nucleic acid diversity analysis revealed that most of the mitochondria protein-coding genes in F. danxiaensis and F. kwangsiensis are highly conserved and may have undergone purifying selection. Furthermore, the collinear and comparative analysis showed that extensive genomic rearrangement events existed among the Malvaceae species. Lastly, a phylogenetic tree based on shared mitochondrial PCGs of 29 species revealed that F. danxiaensis and F. kwangsiensis form a sister group with high support values. Overall, the current study reports two mitogenomes (F. danxiaensis and F. kwangsiensis) in the Firmiana genus for the first time, which will help enhance comprehension of the mitogenome evolutionary patterns within Firmiana and promote the evolutionary and comparative genomic analyses within Malvaceae species.

Complete mitochondrial genomes of the hemiparasitic genus Cymbaria (Orobanchaceae): insights into repeat-mediated recombination, phylogenetic relationships, and horizontal gene transfer

BACKGROUND

The Orobanchaceae family is widely recognized as an exemplary model system for examining the evolutionary dynamics of parasitic plants. However, reports on the mitochondrial genome (mitogenome) of the hemiparasitic tribe Cymbarieae are currently lacking. Here, we sequenced, assembled and characterized the complete mitogenome of the genus Cymbaria L. sensu stricto (C. mongolica and C. daurica).

RESULTS

A total of 51 unique mitochondrial genes, including 33 protein-coding genes, three rRNA genes, and 15 tRNA genes, are shared by the mitogenomes of the two hemiparasitic plants, exhibiting the gene content characteristic of autotrophic plants. The mitogenomes of C. mongolica and C. daurica are characterized by a pentacyclic chromosome structure (their major conformation), with lengths of 1,576,465 bp and 1,539,836 bp, respectively. Moreover, we identified and validated the presence of four minor conformations mediated by four pairs of large repeats (> 1000 bp in size) in C. mongolica and eight minor conformations mediated by six large repeats in C. daurica. We further explored codon usage, RNA editing sites, selective pressure, and nucleotide diversity in two Cymbaria mitogenomes. Phylogenetic analyses of 26 species of Lamiales revealed that the two Cymbaria species form a sister clade to the other lineages of Orobanchaceae. Extensive mitogenomic rearrangements are also observed between Cymbaria and five closely related species. Although we identified mitochondrial plastid sequences in the Cymbaria mitogenomes, The mitochondrial plastid sequences (MTPTs) in their mitogenomes represent only 2.37% and 1.74%, respectively. Additionally, there is minimal evidence of intracellular and horizontal gene transfer, with only a few genes (rpl22, rps3, and ycf2) showing low bootstrap support (BS ≤ 70%) for the relationships with the potential host plants Allium mongolicum, Leymus chinensis, and Saposhnikovia divaricata, respectively.

CONCLUSIONS

We reported the mitochondrial genome in hemiparasitic Cymbaria species for the first time, which are characterized by multiple repeat-mediated recombination and little to no intracellular and horizontal gene transfer. Our findings provide valuable genetic insights for further studies on the mitogenome evolution of hemiparasitic plants.

Comparative genomics and phylogenetic analysis of mitochondrial genomes of Neocinnamomum

BACKGROUND

Neocinnamomum plants are considered a promising feedstock for biodiesel in China, due to the richness in long-chain fatty acids (LCFAs) found in their seeds. However, the mitochondrial genome (mitogenome) of this genus has not yet been systematically described, and the exploration of species relationships within this genus using mitogenome sequences is also an uncharted territory. This has hindered our understanding of mitogenome diversity and the evolutionary relationships within Neocinnamomum.

RESULTS

In this study, a total of 24 individuals representing seven distinct taxa from the genus Neocinnamomum were subjected to Illumina sequencing, and the species N. delavayi was sequenced using Oxford Nanopore sequencing technology. We successfully assembled the mitogenome of N. delavayi, which is 778,066 bp in size and exhibits a single circular structure. The analysis identified 659 dispersed repeats, 211 simple sequence repeats (SSRs), and 30 tandem repeats within the mitogenome. Additionally, 37 homologous fragments, totaling 9929 bp, were found between the mitogenome and the plastid genome (plastome). The codons of 41 protein-coding genes (PCGs) had a preference for ending in A/T, and the codon usage bias of the majority of these genes was influenced by natural selection pressures. Comparative genomic analysis revealed low collinearity and significant gene rearrangements between species. Phylogenetic analysis resulted in the classification of Neocinnamomum into six distinct clades, contradicting previous findings which based on complete plastomes and nuclear ribosomal cistron (nrDNA). In the PCGs of 24 individuals, 86 mutation events were identified, which included three indels and 83 SNPs. Notably, the ccmC gene underwent positive selection in pairwise comparisons of three species pairs. Furthermore, 748 RNA editing sites were predicted within the PCGs of the N. delavayi mitogenome.

CONCLUSIONS

This study enriches our knowledge of the mitogenomes in the family Lauraceae, and provides valuable data and a foundation for genomic evolution research, genetic resource conservation, and molecular breeding in Neocinnamomum.

Genomic insights and phylogenetics of the mitochondrial genome of Cryptocarya

The tropical genus Cryptocarya is known for its valuable timber and the constituents in these trees show potential for medicinal properties. However, the phylogenetic relationships among species in Asia remain unclear. Here, we report the first mitochondrial genome for Cryptocarya kwangtungensis, consisting of 758,020 bp, including 43 protein-coding genes, 23 tRNA genes and three rRNA genes, with 234 simple sequence repeats, and 1275 dispersed repeats, 35 homologous DNA fragments between the mitogenome and the plastome. Comparative genomic analysis indicated frequent recombination events among the sequences of five magnoliids mitogenomes and only five conserved clustered genes. Further phylogenetic analyses based on 91 mitochondrial regions and nuclear ribosomal cistron sequences of 21 species compound three well-resolved congruent groups for the Cryptocarya species in Asia, both of which support the genus divide into three clades, suggesting that the mitogenome sequence can provide strongly supported relationships within the genus in the family Lauraceae.

Assembly and comparative analysis of the complete mitochondrial genome of Daedaleopsissinensis (Polyporaceae, Basidiomycota), contributing to understanding fungal evolution and ecological functions

Daedaleopsissinensis is a crucial wood-decaying fungus with significant lignocellulose-degrading ability, which plays a vital role in the material cycle and energy flow of forest ecosystems. However, the mitochondrial genome of D.sinensis has not yet been revealed. In the present study, the complete mitochondrial genome of D.sinensis was assembled and compared with related species. The mitochondrial genome spans 69,155 bp and has a GC content of 25.0%. It comprises 15 protein-coding genes (PCGs), 26 transfer RNA genes, two ribosomal RNA genes and one DNA polymerase gene (dpo). Herein, we characterised and analysed the codon preferences, variation and evolution of PCGs, repeats, intron dynamics, as well as RNA editing events in the D.sinensis mitochondrial genome. Further, a phylogenetic analysis of D.sinensis and the other 86 Basidiomycota species was performed using mitochondrial genome data. The results revealed that four species, D.confragosa, D.sinensis, D.nitida and Fomesfomentarius, were grouped in a closely-related cluster with high support values, indicating that a close phylogenetic relationship existed between Daedaleopsis and Fomes. This study reported on the initial assembly and annotation of the mitochondrial genome of D.sinensis, which greatly improved the knowledge of the fungus. These results contribute to the limited understanding of the mitochondrial repository of wood-decaying fungi, thereby laying the foundation for subsequent research on fungal evolution and ecological functions.

Recombination and retroprocessing in broomrapes reveal a universal roadmap for mitochondrial evolution in heterotrophic plants

The altered life history strategies of heterotrophic organisms often leave a profound genetic footprint on energy metabolism related functions. In parasitic plants, the reliance on host-derived nutrients and loss of photosynthesis in holoparasites have led to highly degraded to absent plastid genomes, but its impact on mitochondrial genome (mitogenome) evolution has remained controversial. By examining mitogenomes from 45 Orobanchaceae species including three independent transitions to holoparasitism and key evolutionary intermediates, we identified measurable and predictable genetic alterations in genomic shuffling, RNA editing, and intracellular (IGT) and horizontal gene transfer (HGT) en route to a nonphotosynthetic lifestyle. In-depth comparative analyses revealed DNA recombination and repair processes, especially RNA-mediated retroprocessing, as significant drivers for genome structure evolution. In particular, we identified a novel RNA-mediated IGT and HGT mechanism, which has not been demonstrated in cross-species and inter-organelle transfers. Based on this, we propose a generalized dosage effect mechanism to explain the biased transferability of plastid DNA to mitochondria across green plants, especially in heterotrophic lineages like parasites and mycoheterotrophs. Evolutionary rates scaled with these genomic changes, but the direction and strength of selection varied substantially among genes and clades, resulting in high contingency in mitochondrial genome evolution. Finally, we describe a universal roadmap for mitochondrial evolution in heterotrophic plants where increased recombination and repair activities, rather than relaxed selection alone, lead to differentiated genome structure compared to free-living species.

Insights into the nuclear-organelle DNA integration in Cicuta virosa (Apiaceae) provided by complete plastid and mitochondrial genomes

BACKGROUND

Gene transfer between the organelles and the nucleus plays a central role in shaping plant genome evolution. The identification and analysis of nuclear DNA of plastid (NUPTs) and mitochondrial (NUMTs) origins are important for exploring the extent of intracellular DNA transfer in genomes.

RESULTS

We report the complete plastid and mitochondrial genomes (plastome and mitogenome) of Cicuta virosa (Apiaceae) as well as a draft nuclear genome using high-fidelity (HiFi) PacBio sequencing technologies. The C. virosa plastome (154,449 bp) is highly conserved, with a quadripartite structure, whereas the mitogenome (406,112 bp) exhibits two chromosomes (352,718 bp and 53,394 bp). The mitochondrial-encoded genes (rpl2, rps14, rps19, and sdh3) were successfully transferred to the nuclear genome. Our findings revealed extensive DNA transfer from organelles to the nucleus, with 6,686 NUPTs and 6,237 NUMTs detected, covering nearly the entire plastome (99.93%) and a substantial portion of the mitogenome (77.04%). These transfers exhibit a range of sequence identities (80-100%), suggesting multiple transfer events over evolutionary timescales. Recent DNA transfer between organelles and the nucleus is more frequent in mitochondria than that in plastids.

CONCLUSIONS

This study contributes to the understanding of ongoing genome evolution in C. virosa and underscores the significance of the organelle-nuclear genome interplay in plant species. Our findings provide valuable insights into the evolutionary processes that shape organelle genomes in Apiaceae, with implications for broader plant genome evolution.

Assembly and characterization of the complete mitogenome of Bauhinia purpurea (Leguminosae)

OBJECTIVE

Mitochondrial genome sequences are very useful for understanding the mitogenome evolution itself and reconstructing phylogeny of different plant lineages. Bauhinia purpurea, a species from the legume family Leguminosae, is widely distributed in South China and has high ornamental value. Here, we sequenced and assembled the mitogenome of B. purpurea to provide a useful genetic resource for further evolutionary studies.

DATA DESCRIPTION

We assembled and characterized the complete mitogenome of B. purpurea based on Illumina sequence data. The mitogenome size was 525,727 bp, and its GC content was 45.38%. A total of 35 protein-coding genes, 16 tRNA genes, and 3 rRNA genes were identified in the mitogenome. We also identified 124 pairs of repeats and 6 mitogenome sequences of plastid origin (MTPTs). These MTPTs range from 108 bp to 751 bp, covering 0.65% of the mitogenome.

Comprehensive analysis of 111 Pleuronectiformes mitochondrial genomes: insights into structure, conservation, variation and evolution

BACKGROUND

Pleuronectiformes, also known as flatfish, are important model and economic animals. However, a comprehensive genome survey of their important organelles, mitochondria, has been limited. Therefore, we aim to analyze the genomic structure, codon preference, nucleotide diversity, selective pressure and repeat sequences, as well as reconstruct the phylogenetic relationship using the mitochondrial genomes of 111 flatfish species.

RESULTS

Our analysis revealed a conserved gene content of protein-coding genes and rRNA genes, but varying numbers of tRNA genes and control regions across species. Various gene rearrangements were found in flatfish species, especially for the rearrangement of nad5-nad6-cytb block in Samaridae family, the swapping rearrangement of nad6 and cytb gene in Bothidae family, as well as the control region translocation and tRNA-Gln gene inversion in the subfamily Cynoglossinae, suggesting their unique evolutionary history and/or functional benefit. Codon usage showed obvious biases, with adenine being the most frequent nucleotide at the third codon position. Nucleotide diversity and selective pressure analysis suggested that different protein-coding genes underwent varying degrees of evolutionary pressure, with cytb and cox genes being the most conserved ones. Phylogenetic analysis using both whole mitogenome information and concatenated independently aligned protein-coding genes largely mirrored the taxonomic classification of the species, but showed different phylogeny. The identification of simple sequence repeats and various long repetitive sequences provided additional complexity of genome organization and offered markers for evolutionary studies and breeding practices.

CONCLUSIONS

This study represents a significant step forward in our comprehension of the flatfish mitochondrial genomes, providing valuable insights into the structure, conservation and variation within flatfish mitogenomes, with implications for understanding their evolutionary history, functional genomics and fisheries management. Future research can delve deeper into conservation biology, evolutionary biology and functional usages of variations.

Mitochondrial genome structural variants and candidate cytoplasmic male sterility-related gene in sugarcane

BACKGROUND

Sugarcane is a crucial crop for both sugar and bioethanol production. The nobilization breeding and utilization of wild germplasm have significantly enhanced its productivity. However, the pollen sterility in Saccharum officinarum restricts its role to being a female parent in crosses with Saccharum spontaneum during nobilization breeding, resulting in a narrow genetic basis for modern sugarcane cultivars. Mitochondria, often referred to as the intracellular "energy factories", provide energy for plant life activities, and are also implicated in cytoplasmic male sterility (CMS).

RESULTS

We performed mitochondrial genome assembly and structural analysis of two Saccharum founding species. We discovered that the proportions of repeat sequences are the primary factor contributing to the variations in mitochondrial genome structure and size between the two Saccharum species. Heterologous expression of the mitochondrial chimeric gene ORF113, which is highly expressed in male-sterile S. officinarum flowers, significantly inhibits growth and ATP synthesis in yeast cells, making it a key candidate CMS-related gene in sugarcane. Furthermore, we developed two co-dominant simple sequence repeat (SSR) markers based on the mitochondrial genome, which can effectively distinguish the cytoplasmic types of the two Saccharum species.

CONCLUSION

In this study, we identified structural variants and developed SSR molecular markers in the mitochondrial genomes of both S. officinarum and S. spontaneum. We also identified a novel mitochondrial chimeric ORF as a key candidate CMS-related gene. These findings offer valuable insights into variety identification, genetic resource development, and cross-breeding strategies in sugarcane.

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.

The evolutionarily diverged single-stranded DNA-binding proteins SSB1/SSB2 differentially affect the replication, recombination and mutation of organellar genomes in Arabidopsis thaliana

Single-stranded DNA-binding proteins (SSBs) play essential roles in the replication, recombination and repair processes of organellar DNA molecules. In Arabidopsis thaliana, SSBs are encoded by a small family of two genes (SSB1 and SSB2). However, the functional divergence of these two SSB copies in plants remains largely unknown, and detailed studies regarding their roles in the replication and recombination of organellar genomes are still incomplete. In this study, phylogenetic, gene structure and protein motif analyses all suggested that SSB1 and SSB2 probably diverged during the early evolution of seed plants. Based on accurate long-read sequencing results, ssb1 and ssb2 mutants had decreased copy numbers for both mitochondrial DNA (mtDNA) and plastid DNA (ptDNA), accompanied by a slight increase in structural rearrangements mediated by intermediate-sized repeats in mt genome and small-scale variants in both genomes. Our findings provide an important foundation for further investigating the effects of DNA dosage in the regulation of mutation frequencies in plant organellar genomes.

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.

Organelle genome assembly, annotation, and comparative analyses of Typha latifolia and T. domingensis: two keystone species for wetlands worldwide

Typha is a cosmopolitan aquatic plant genus that includes species with widespread global distributions. In previous studies, a revised molecular phylogeny was inferred using seven plastid loci from nine Typha species across different geographic regions. By utilizing complete organellar genomes, we aim to provide a more comprehensive dataset that offers a robust phylogenetic signal for resolving Typha species evolutionary relationships. Here, we assembled T. latifolia and T. domingensis mitochondrial genomes (mitogenomes) using a combination of short-read and long-read data (PacBio, ONT). The mitogenomes of both species are assembled into single circular molecules of 395,136 bp and 395,140 bp in length, respectively, with a similar GC content of 46.7%. A total of 39 protein-coding genes, 17 tRNA genes, and 3 rRNA genes were annotated in both mitogenomes. The plastid genomes (plastomes) of both species possess typical quadripartite structures observed across most plants, with sizes of 161,545 bp and 161,230 bp. The overall average GC content of the plastomes of both species was 36.6%. The comparative analysis of the plastome and mitogenome revealed that 12 mitogenome DNA fragments share similar sequences with in the repeat regions of the corresponding plastomes, suggesting a past transfer of repeat regions into the mitogenome. Additionally, the mitogenomes of the two Typha species exhibited high sequence conservation with several syntenic blocks. Phylogenetic analysis of the organellar genomes of the two Typha species and 10 related species produced congruent phylogenetic trees. The availability of these organellar genomes from two Typha species provide valuable genetic resources for studying the evolution of Typhaceae and will improve taxonomic classifications within the family.

The first complete mitochondrial genome of sesame (Sesamum indicum L.)

Sesame (Sesamum indicum L.), an important oilseed crop, has garnered considerable interest. The nuclear and chloroplast genomes of sesame have been extensively applied to sesame genetics and genomics research. The mitochondrial (mt) genome of sesame, however, has not been sequenced and annotated. In order to solve this issue, we reconstructed the first mt genome of sesame using third-generation sequencing data. The sesame mt genome was 724,998 bp in size and had 22 circular chromosomes. A total of 66 genes were annotated, including 37 protein-coding genes, 26 transfer RNAs, and three ribosomal RNAs. We investigated the codon usage patterns, simple sequence repeats, long tandem repeats, and dispersed repeats of the sesame mt genome. Furthermore, we investigated the DNA transfer from chloroplast to mitochondrion and compared the sesame mt genome to two other Lamiales mt genomes. Given the economic importance of this crop, our presented sesame mt genome is a valuable genomic resource and will allow for more comprehensive studies on sesame and related species.

Mitochondrial Genome Assembly and Comparative Analysis of Chionanthus Retusus (Oleaceae)

Background/Objectives: Chionanthus retusus Lindl. & Paxton is an ornamental tree species native to North China. Research on the mitochondrial genome can elucidate the evolution and biological characteristics of C. retusus and better protect this important species. Methods and Results: This work aimed to clarify the evolutionary and phylogenetic links by sequencing, assembling, annotating, and analyzing the entire mitochondrial genome of C. retusus. The single-loop structure that made up the mitochondrial genome had a total length of 657,640 bp and a GC content of 44.52%. In total, 37 unique protein-coding genes, 20 tRNA genes, and 3 rRNA genes were identified. Numerous repeat sequences and migrating fragments of chloroplast sequences were found. Using the mitochondrial protein-coding genes to construct evolutionary trees, it was found that the closest relative of C. retusus is C. rupicola (Lingelsh.) Kiew. Conclusions: This research represents the first comprehensive set of data on the mitochondrial genome of an ancient (>500 yr) C. retusus specimen. In addition to elucidating the biological characteristics of C. retusus. The findings contribute to the Oleaceae mitochondrial genome database and offer valuable insights for future studies in molecular breeding, evolutionary biology, and genetic diversity conservation.

Assembly and analysis of the complete mitochondrial genome of Carya illinoinensis to provide insights into the conserved sequences of tRNA genes

Carya illinoinensis is an economically important nut tree, and its chloroplast (cp.) genome has been reported; however, its mitochondrial (mt) genome remains unknown. In the present study, we assembled the first mt genome of C. illinoinensis. The circular mt genome of C. illinoinensis is 495,205 bp long, with 37 protein-coding genes(PCGs), 24 tRNA genes, and 3 rRNA genes. All the tRNAs could be folded into typical cloverleaf secondary structures, with lengths of 58-88 bp. A conserved U-U-C-x-A-x2 consensus nucleotide sequence was discovered in the Ψ-loops of tRNA sequences. In addition, 447 dispersed repeats were detected, as well as found 482 RNA editing sites and 9,960 codons in the mt genome. Furthermore, a total of 27 DNA sequences with a length of 43,277 bp were transferred from the cp. to the mt genome, and eight integrated cp-derived genes (trnL-CAA, trnV-GAC, trnD-GUC, trnW-CCA, trnN-GUU, trnH-GUG, trnM-CAU, and rps7) were identified. We also obtained 1,086 hits, including 364.023 kp of nuclear genome sequences, that were transferred to the mt genome. To determine the evolutionary position of C. illinoinensis, we conducted a phylogenetic analysis of the mitogenomes of C. illinoinensis and 14 other taxa. The results strongly suggested that C. illinoinensis and Fagus sylvatica formed a single clade with 100% bootstrap support. This study sequenced comprehensive data on the C. illinoinensis mitochondrial genome and provided insights into the conserved sequences of tRNA genes, which could facilitate evolutionary research in other Carya trees in the future.

The complete mitochondrial genome assembly of Capsicum pubescens reveals key evolutionary characteristics of mitochondrial genes of two Capsicum subspecies

BACKGROUND

Pepper (Capsicum pubescens), one of five domesticated pepper species, has unique characteristics, such as numerous hairs on the epidermis of its leaves and stems, black seeds, and vibrant purple flowers. To date, no studies have reported on the complete assembly of the mitochondrial genome (mitogenome) of C. pubescens. Understanding the mitogenome is crucial for further research on C. pubescens.

RESULTS

In our study, we successfully assembled the first mitogenome of C. pubescens, which was assigned the GenBank accession number OP957066. This mitogenome has a length of 454,165 bp and exhibits the typical circular structure observed in most mitogenomes. We annotated a total of 70 genes, including 35 protein-coding genes (PCGs), 30 tRNA genes, 3 rRNA genes, and 2 pseudogenes. Compared to the other three pepper mitogenomes (KJ865409, KJ865410, and MN196478), C. pubescens OP957066 exhibited four unique PCGs (atp4, atp8, mttB, and rps1), while two PCGs (rpl10 and rps3) were absent. Notably, each of the three pepper mitogenomes from C. annuum (KJ865409, KJ865410, and MN196478) experienced the loss of four PCGs (atp4, atp8, mttB, and rps1). To further explore the evolutionary relationships, we reconstructed a phylogenetic tree using the mitogenomes of C. pubescens and fourteen other species. Structural comparison and synteny analysis of the above four pepper mitogenomes revealed that C. pubescens shares high sequence similarity with KJ865409 and that C. pubescens has rearranged with the other three pepper mitogenomes. Interestingly, we observed 72 similar sequences between the mitochondrial and chloroplast genomes, which accounted for 12.60% of the mitogenome, with a total length of 57,207 bp. These sequences encompassed 12 tRNA genes and the rRNA gene (rrn18). Remarkably, selective pressure analysis suggested that the nad5 gene underwent obvious positive selection. Furthermore, a single-base mutation in three genes (nad1, nad2, and nad4) resulted in an amino acid change.

CONCLUSION

This study provides a high-quality mitogenome of pepper, providing valuable molecular data for future investigations into the exchange of genetic information between pepper organelle genomes.

Mitochondrial Genome Assembly and Structural Characteristics Analysis of Gentiana rigescens

Gentiana rigescens, an alpine plant with significant medicinal value, possesses a complex genetic background. However, comprehensive genomic research on G. rigescens is still lacking, particularly concerning its organelle genome. In this study, G. rigescens was studied to sequence the mitochondrial genome (mitogenome) and ascertain the assembly, informational content, and developmental expression of the mitogenome. The mitogenome of G. rigescens was 393,595 bp in length and comprised four circular chromosomes ranging in size from 6646 bp to 362,358 bp. The GC content was 43.73%. The mitogenome featured 30 distinct protein-coding genes, 26 tRNA genes, and 3 rRNA genes. The mitogenome of G. rigescens also revealed 70 SSRs, which were mostly tetra-nucleotides. In addition, 48 homologous fragments were found between the mitogenome and the chloroplast genome, with the longest measuring 23,330 bp. The documentation of the mitochondrial genome of G. rigescens is instrumental in advancing the understanding of its physiological development. Decoding the G. rigescens mitogenome will offer valuable genetic material for phylogenetic research on Gentianaceae and enhance the use of species germplasm resources.

Repeat-mediated recombination results in Complex DNA structure of the mitochondrial genome of Trachelospermum jasminoides

BACKGROUND

Trachelospermum jasminoides has medicinal and ornamental value and is widely distributed in China. Although the chloroplast genome has been documented, the mitochondrial genome has not yet been studied.

RESULTS

The mitochondrial genome of T. jasminoides was assembled and functionally annotated using Illumina and nanopore reads. The mitochondrial genome comprises a master circular molecular structure of 605,764 bp and encodes 65 genes: 39 protein-coding genes, 23 transfer RNA (tRNA) genes and 3 ribosomal RNA genes. In addition to the single circular conformation, we found many alternative conformations of the T. jasminoides mitochondrial genome mediated by 42 repetitive sequences. Six repetitive sequences (DRS01-DRS06) were supported by nanopore long reads, polymerase chain reaction (PCR) amplifications, and Sanger sequencing of the PCR products. Eleven homologous fragments were identified by comparing the mitochondrial and chloroplast genome sequences, including three complete tRNA genes. Moreover, 531 edited RNA sites were identified in the protein-coding sequences based on RNA sequencing data, with nad4 having the highest number of sites (54).

CONCLUSION

To our knowledge, this is the first description of the mitochondrial genome of T. jasminoides. Our results demonstrate the existence of multiple conformations. These findings lay a foundation for understanding the genetics and evolutionary dynamics of Apocynaceae.

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.

Assembly and comparative analysis of the first complete mitochondrial genome of Salix psammophila, a good windbreak and sand fixation shrub

Salix psammophila, commonly known as the sandlive willow, is a vital shrub species within the Salicaceae family, particularly significant for its ecological role in regions susceptible to desertification and sandy soils. In this study, we assembled the complete S. psammophila mitochondrial genome using Pacbio HiFi third-generation sequencing data. The genome was found to be a typical single circular structure, with a total length of 715,555 bp and a GC content of 44.89%. We annotated 33 unique protein-coding genes (PCGs), which included 24 core mitochondrial genes and 9 variable genes, as well as 18 tRNA genes (5 of which were multicopy genes) and 3 rRNA genes. Comparative analysis of the PCGs from the mitochondrial genomes of S. psammophila, Populus deltoides, Populus simonii, Salix wilsonii, and Salix suchowensis revealed that these genes are relatively conserved within the Salicaceae family, with variability primarily occurring in the ribosomal protein genes. The absence of the rps14, which encodes a ribosomal protein, may have played a role in the evolution of stress tolerance in Salicaceae plants. Additionally, we identified 232 SSRs, 19 tandem repeat sequences, and 236 dispersed repeat sequences in the S. psammophila mitochondrial genome, with palindromic and forward repeats being the most abundant. The longest palindromic repeat measured 260 bp, while the longest forward repeat was 86,068 bp. Furthermore, 324 potential RNA editing sites were discovered, all involving C-to-U edits, with the nad4 having the highest number of edits. These findings provide valuable insights into the phylogenetic and genetic research of Salicaceae plants.

Characterization of the complete mitogenome of Tiarella polyphylla, commonly known as Asian foamflower: insights into the multi-chromosomes structure and DNA transfers

BACKGROUND

Tiarella polyphylla D. Don has been traditionally used to cure asthma and skin eruptions. However, the sequence and the structure of the mitogenome of T. polyphylla remained elusive, limiting the genomic and evolution analysis based on the mitogenome.

RESULTS

Using a combination of Illumina and Nanopore sequencing reads, we de novo assembled the complete mitogenome of T. polyphylla. In addition to unveiling the major configuration of the T. polyphylla mitogenome was three circular chromosomes with lengths of 430,435 bp, 126,943 bp, and 55,269 bp, we revealed five (R01-R05) and one (R06) repetitive sequence could mediate the intra- and inter-chromosomal recombination, respectively. Furthermore, we identified 208 short and 25 long tandem segments, seven cp-derived mtDNAs, 106 segments of mtDNAs transferred to the nuclear genome, and 653 predicted RNA editing sites. Based on the sequence of the mitogenomes, we obtained the resolved phylogeny of the seven Saxifragales species.

CONCLUSIONS

These results presented the mitogenome features and expanded its potential applications in phylogenetics, species identification, and cytoplasmic male sterility (CMS) in the future.

Mitogenomes comparison of 3 species of Asparagus L shedding light on their functions due to domestication and adaptative evolution

BACKGROUND

Asparagus L., widely distributed in the old world is a genus under Asparagaceae, Asparagales. The species of the genus were mainly used as vegetables, traditional medicines as well as ornamental plants. However, the evolution and functions of mitochondrial (Mt) genomes (mitogenomes) remains largely unknown. In this study, the typical herbal medicine A. taliensis and ornamental plant A. setaceus were used to assemble and annotate the mitogenomes, and the resulting mitogenomes were further compared with published mitogenome of A. officinalis for the analysis of their functions in the context of domestication and adaptative evolution.

RESULTS

The mitochondrial genomes of both A. taliensis and A. setaceus were assembled as complete circular ones. The phylogenetic trees based on conserved protein-coding genes of Mt genomes and whole chloroplast (Cp) genomes showed that, the phylogenetic relationship of the sampled 13 species of Asparagus L. were not exactly consistent. The collinear analyses between the nuclear (Nu) and Mt genomes confirmed the existence of mutual horizontal genes transfers (HGTs) between Nu and Mt genomes within these species. Based on RNAseq data, the Mt RNA editing were predicted and atp1 and ccmB RNA editing of A. taliensis were further confirmed by DNA sequencing. Simultaneously homologous search found 5 Nu coding gene families including pentatricopeptide-repeats (PPRs) involved in Mt RNA editing. Finally, the Mt genome variations, gene expressions and mutual HGTs between Nu and Mt were detected with correlation to the growth and developmental phenotypes respectively. The results suggest that, both Mt and Nu genomes co-evolved and maintained the Mt organella replication and energy production through TCA and oxidative phosphorylation .

CONCLUSION

The assembled and annotated complete mitogenomes of both A. taliensis and A. setaceus provide valuable information for their phylogeny and concerted action of Nu and Mt genomes to maintain the energy production system of Asparagus L. in the context of domestication and adaptation to environmental niches.

Complete mitochondrial genome of Agropyron cristatum reveals gene transfer and RNA editing events

BACKGROUND

As an important forage in arid and semi-arid regions, Agropyron cristatum provides livestock with exceptionally high nutritional value. Additionally, A. cristatum exhibits outstanding genetic characteristics to endure drought and disease. Therefore, rich genetic diversity serves as a cornerstone for the improvement of major food crops. The purposes of this study were to systematically describe mitogenome of A.cristatum and preliminarily analyze its internal variations.

RESULT

The A. cristatum mitogenome was a single-ring molecular structure of 381,065 bp that comprised 52 genes, including 35 protein-coding, 3 rRNA and 14 tRNA genes. Among these, two pseudoprotein-coding genes and multiple copies of tRNA genes were observed. A total of 320 repetitive sequences was found to cover more than 10% of the mitogenome (105 simple sequences, 185 dispersed and 30 tandem repeats), which led to a large number of fragment rearrangements in the mitogenome of A. cristatum. Leucine was the most frequent amino acid (n = 1087,10.8%) in the protein-coding genes of A. cristatum mitogenome, and the highest usage codon was ATG (initiation codon). The number of A/T changes at the third base of the codon was much higher than that of G/C. Among 23 PCGs, the range of Pi values is from 0.0021 to 0.0539, with an average of 0.013. Additionally, 81 RNA editing sites were predicted, which were considerably fewer than those reported in other plant mitogenomes. Most of the RNA editing site base positions were concentrated at the first and second codon bases, which were C to T transitions. Moreover, we identified 95 sequence fragments (total length of 34, 343 bp) that were transferred from the chloroplast to mitochondria genes, introns, and intergenic regions. The stability of the tRNA genes was maintained during this process. Selection pressure analysis of 23 protein-coding genes shared by 15 Poaceae plants, showed that most genes were subjected to purifying selection during evolution, whereas rps4, cob, mttB, and ccmB underwent positive selection in different plants. Finally, a phylogenetic tree was constructed based on 22 plant mitogenomes, which showed that Agropyron plants have a high degree of independent heritability in Triticeae.

CONCLUSION

The findings of this study provide new data for a better understanding of A. cristatum genes, and demonstrate that mitogenomes are suitable for the study of plant classifications, such as those of Agropyron. Moreover, it provides a reference for further exploration of the phylogenetic relationships within Agropyron species, and establishes a theoretical basis for the subsequent development and utilization of A. cristatum plant germplasm resources.

Horizontally transferred mitochondrial DNA tracts become circular by microhomology-mediated repair pathways

The holoparasitic plant Lophophytum mirabile exhibits remarkable levels of mitochondrial horizontal gene transfer (HGT). Gathering comparative data from other individuals and host plants can provide insights into the HGT process. We sequenced the mitochondrial genome (mtDNA) from individuals of two species of Lophophytum and from mimosoid hosts. We applied a stringent phylogenomic approach to elucidate the origin of the whole mtDNAs, estimate the timing of the transfers, and understand the molecular mechanisms involved. Ancestral and recent HGT events replaced and enlarged the multichromosomal mtDNA of Lophophytum spp., with the foreign DNA ascending to 74%. A total of 14 foreign mitochondrial chromosomes originated from continuous regions in the host mtDNA flanked by short direct repeats. These foreign tracts are circularized by microhomology-mediated repair pathways and replicate independently until they are lost or they eventually recombine with other chromosomes. The foreign noncoding chromosomes are variably present in the population and likely evolve by genetic drift. We present the 'circle-mediated HGT' model in which foreign mitochondrial DNA tracts become circular and are maintained as plasmid-like molecules. This model challenges the conventional belief that foreign DNA must be integrated into the recipient genome for successful HGT.

Accumulation of Large Lineage-Specific Repeats Coincides with Sequence Acceleration and Structural Rearrangement in Plantago Plastomes

Repeats can mediate rearrangements and recombination in plant mitochondrial genomes and plastid genomes. While repeat accumulations are linked to heightened evolutionary rates and complex structures in specific lineages, debates persist regarding the extent of their influence on sequence and structural evolution. In this study, 75 Plantago plastomes were analyzed to investigate the relationships between repeats, nucleotide substitution rates, and structural variations. Extensive repeat accumulations were associated with significant rearrangements and inversions in the large inverted repeats (IRs), suggesting that repeats contribute to rearrangement hotspots. Repeats caused infrequent recombination that potentially led to substoichiometric shifting, supported by long-read sequencing. Repeats were implicated in elevating evolutionary rates by facilitating localized hypermutation, likely through DNA damage and repair processes. This study also observed a decrease in nucleotide substitution rates for loci translocating into IRs, supporting the role of biased gene conversion in maintaining lower substitution rates. Combined with known parallel changes in mitogenomes, it is proposed that potential dysfunction in nuclear-encoded genes associated with DNA replication, recombination, and repair may drive the evolution of Plantago organellar genomes. These findings contribute to understanding how repeats impact organellar evolution and stability, particularly in rapidly evolving plant lineages.

The Complete Mitogenome of Apostasia fujianica Y.Li & S.Lan and Comparative Analysis of Mitogenomes across Orchidaceae

Apostasia fujianica belongs to the genus Apostasia and is part of the basal lineage in the phylogenetic tree of the Orchidaceae. Currently, there are only ten reported complete mitochondrial genomes in orchids, which greatly hinders the understanding of mitochondrial evolution in Orchidaceae. Therefore, we assembled and annotated the mitochondrial genome of A. fujianica, which has a length of 573,612 bp and a GC content of 44.5%. We annotated a total of 44 genes, including 30 protein-coding genes, 12 tRNA genes, and two rRNA genes. We also performed relative synonymous codon usage (RSCU) analysis, repeat sequence analysis, intergenomic transfer (IGT) analysis, and Ka/Ks analysis for A. fujianica and conducted RNA editing site analysis on the mitochondrial genomes of eight orchid species. We found that most protein-coding genes are under purifying selection, but nad6 is under positive selection, with a Ka/Ks value of 1.35. During the IGT event in A. fujianica's mitogenome, the trnN-GUU, trnD-GUC, trnW-CCA, trnP-UGG, and psaJ genes were identified as having transferred from the plastid to the mitochondrion. Compared to other monocots, the family Orchidaceae appears to have lost the rpl10, rpl14, sdh3, and sdh4 genes. Additionally, to further elucidate the evolutionary relationships among monocots, we constructed a phylogenetic tree based on the complete mitogenomes of monocots. Our study results provide valuable data on the mitogenome of A. fujianica and lay the groundwork for future research on genetic variation, evolutionary relationships, and breeding of Orchidaceae.

De Novo Hybrid Assembly Unveils Multi-Chromosomal Mitochondrial Genomes in Ludwigia Species, Highlighting Genomic Recombination, Gene Transfer, and RNA Editing Events

Biological invasions have been identified as the fifth cause of biodiversity loss, and their subsequent dispersal represents a major ecological challenge. The aquatic invasive species Ludwigia grandiflora subsp. hexapetala (Lgh) and Ludwigia peploides subsp. montevidensis (Lpm) are largely distributed in aquatic environments in North America and in Europe. However, they also present worrying terrestrial forms that are able to colonize wet meadows. To comprehend the mechanisms of the terrestrial adaptation of Lgh and Lpm, it is necessary to develop their genomic resources, which are currently poorly documented. We performed de novo assembly of the mitogenomes of Lgh and Lpm through hybrid assemblies, combining short reads (SR) and/or long reads (LR) before annotating both mitogenomes. We successfully assembled the mitogenomes of Lgh and Lpm into two circular molecules each, resulting in a combined total length of 711,578 bp and 722,518 bp, respectively. Notably, both the Lgh and Lpm molecules contained plastome-origin sequences, comprising 7.8% of the mitochondrial genome length. Additionally, we identified recombinations that were mediated by large repeats, suggesting the presence of multiple alternative conformations. In conclusion, our study presents the first high-quality mitogenomes of Lpm and Lgh, which are the only ones in the Myrtales order found as two circular molecules.

Plant organellar genomes: much done, much more to do

Plastids and mitochondria are the only organelles that possess genomes of endosymbiotic origin. In recent decades, advances in sequencing technologies have contributed to a meteoric rise in the number of published organellar genomes, and have revealed greatly divergent evolutionary trajectories. In this review, we quantify the abundance and distribution of sequenced plant organellar genomes across the plant tree of life. We compare numerous genomic features between the two organellar genomes, with an emphasis on evolutionary trajectories, transfers, the current state of organellar genome editing by transcriptional activator-like effector nucleases (TALENs), transcription activator-like effector (TALE)-mediated deaminase, and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein (Cas), as well as genetic transformation. Finally, we propose future research to understand these different evolutionary trajectories, and genome-editing strategies to promote functional studies and eventually improve organellar genomes.

A High-Quality Assembly and Comparative Analysis of the Mitogenome of Actinidia macrosperma

The mitochondrial genome (mitogenome) of Actinidia macrosperma, a traditional medicinal plant within the Actinidia genus, remains relatively understudied. This study aimed to sequence the mitogenome of A. macrosperma, determining its assembly, informational content, and developmental expression. The results revealed that the mitogenome of A. macrosperma is circular, spanning 752,501 bp with a GC content of 46.16%. It comprises 63 unique genes, including 39 protein-coding genes (PCGs), 23 tRNA genes, and three rRNA genes. Moreover, the mitogenome was found to contain 63 SSRs, predominantly mono-nucleotides, as well as 25 tandem repeats and 650 pairs of dispersed repeats, each with lengths equal to or greater than 60, mainly comprising forward repeats and palindromic repeats. Moreover, 53 homologous fragments were identified between the mitogenome and chloroplast genome (cp-genome), with the longest segment measuring 4296 bp. This study represents the initial report on the mitogenome of the A. macrosperma, providing crucial genetic materials for phylogenetic research within the Actinidia genus and promoting the exploitation of species genetic resources.

Mitochondrial genome variation and intergenomic sequence transfers in Hevea species

Among the Hevea species, rubber tree (Hevea brasiliensis) is the most important source of natural rubber. In previous studies, we sequenced the complete nuclear and chloroplast genomes of Hevea species, providing an invaluable resource for studying their phylogeny, disease resistance, and breeding. However, given that plant mitochondrial genomes are more complex and more difficult to assemble than that of the other organelles, little is known about their mitochondrial genome, which limits the comprehensive understanding of Hevea genomic evolution. In this study, we sequenced and assembled the mitochondrial genomes of four Hevea species. The four mitochondrial genomes had consistent GC contents, codon usages and AT skews. However, there were significant differences in the genome lengths and sequence repeats. Specifically, the circular mitochondrial genomes of the four Hevea species ranged from 935,732 to 1,402,206 bp, with 34-35 unique protein-coding genes, 35-38 tRNA genes, and 6-13 rRNA genes. In addition, there were 17,294-46,552 bp intergenomic transfer fragments between the chloroplast and mitochondrial genomes, consisting of eight intact genes (psaA, rrn16S, tRNA-Val, rrn5S, rrn4.5S, tRNA-Arg, tRNA-Asp, and tRNA-Asn), intergenic spacer regions and partial gene sequences. The evolutionary position of Hevea species, crucial for understanding its adaptive strategies and relation to other species, was verified by phylogenetic analysis based on the protein-coding genes in the mitochondrial genomes of 21 Malpighiales species. The findings from this study not only provide valuable insights into the structure and evolution of the Hevea mitochondrial genome but also lay the foundation for further molecular, evolutionary studies, and genomic breeding studies on rubber tree and other Hevea species, thereby potentially informing conservation and utilization strategies.

Complete organelle genomes of Korean fir, Abies koreana and phylogenomics of the gymnosperm genus Abies using nuclear and cytoplasmic DNA sequence data

Abies koreana E.H.Wilson is an endangered evergreen coniferous tree that is native to high altitudes in South Korea and susceptible to the effects of climate change. Hybridization and reticulate evolution have been reported in the genus; therefore, multigene datasets from nuclear and cytoplasmic genomes are needed to better understand its evolutionary history. Using the Illumina NovaSeq 6000 and Oxford Nanopore Technologies (ONT) PromethION platforms, we generated complete mitochondrial (1,174,803 bp) and plastid (121,341 bp) genomes from A. koreana. The mitochondrial genome is highly dynamic, transitioning from cis- to trans-splicing and breaking conserved gene clusters. In the plastome, the ONT reads revealed two structural conformations of A. koreana. The short inverted repeats (1186 bp) of the A. koreana plastome are associated with different structural types. Transcriptomic sequencing revealed 1356 sites of C-to-U RNA editing in the 41 mitochondrial genes. Using A. koreana as a reference, we additionally produced nuclear and organelle genomic sequences from eight Abies species and generated multiple datasets for maximum likelihood and network analyses. Three sections (Balsamea, Momi, and Pseudopicea) were well grouped in the nuclear phylogeny, but the phylogenomic relationships showed conflicting signals in the mitochondrial and plastid genomes, indicating a complicated evolutionary history that may have included introgressive hybridization. The obtained data illustrate that phylogenomic analyses based on sequences from differently inherited organelle genomes have resulted in conflicting trees. Organelle capture, organelle genome recombination, and incomplete lineage sorting in an ancestral heteroplasmic individual can contribute to phylogenomic discordance. We provide strong support for the relationships within Abies and new insights into the phylogenomic complexity of this genus.

Plasticity of repetitive sequences demonstrated by the complete mitochondrial genome of Eucalyptus camaldulensis

The tree Eucalyptus camaldulensis is a ubiquitous member of the Eucalyptus genus, which includes several hundred species. Despite the extensive sequencing and assembly of nuclear genomes from various eucalypts, the genus has only one fully annotated and complete mitochondrial genome (mitogenome). Plant mitochondria are characterized by dynamic genomic rearrangements, facilitated by repeat content, a feature that has hindered the assembly of plant mitogenomes. This complexity is evident in the paucity of available mitogenomes. This study, to the best of our knowledge, presents the first E. camaldulensis mitogenome. Our findings suggest the presence of multiple isomeric forms of the E. camaldulensis mitogenome and provide novel insights into minor rearrangements triggered by nested repeat sequences. A comparative sequence analysis of the E. camaldulensis and E. grandis mitogenomes unveils evolutionary changes between the two genomes. A significant divergence is the evolution of a large repeat sequence, which may have contributed to the differences observed between the two genomes. The largest repeat sequences in the E. camaldulensis mitogenome align well with significant yet unexplained structural variations in the E. grandis mitogenome, highlighting the adaptability of repeat sequences in plant mitogenomes.

Mitochondrial genome features and systematic evolution of diospyros kaki thunb 'Taishuu'

BACKGROUND

'Taishuu' has a crisp texture, abundant juice, and sweet flavor with hints of cantaloupe. The availability of mitochondrial genome data of Diospyros species is far from the known number of species.

RESULTS

The sequencing data were assembled into a closed circular mitochondrial chromosome with a 421,308 bp length and a 45.79% GC content. The mitochondrial genome comprised 40 protein-coding, 24 tRNA, and three rRNA genes. The most common codons for arginine (Arg), proline (Pro), glycine (Gly), tryptophan (Trp), valine (Val), alanine (Ala), and leucine (Leu) were AGA, CCA, GGA, UGG, GUA, GCA, and CUA, respectively. The start codon for cox1 and nad4L protein-coding genes was ACG (ATG), whereas the remaining protein-coding genes started with ATG. There are four types of stop codons: CGA, TAA, TAG, and TGA, with TAA being the most frequently used stop codon (45.24%). In the D. kaki Thunb. 'Taishuu' mitochondrial genome, a total of 645 repeat sequences were identified, including 125 SSRs, 7 tandem repeats, and 513 dispersed repeats. Collinearity analysis revealed a close relationship between D. kaki Thunb. 'Taishuu' and Diospyros oleifera, with conserved homologous gene fragments shared among these species in large regions of the mitochondrial genome. The protein-coding genes ccmB and nad4L were observed to undergo positive selection. Analysis of homologous sequences between chloroplasts and mitochondria identified 28 homologous segments, with a total length of 24,075 bp, accounting for 5.71% of the mitochondrial genome. These homologous segments contain 8 annotated genes, including 6 tRNA genes and 2 protein-coding genes (rrn18 and ccmC). There are 23 homologous genes between chloroplasts and nuclei. Mitochondria, chloroplasts, and nuclei share two homologous genes, which are trnV-GAC and trnW-CCA.

CONCLUSION

In conclusion, a high-quality chromosome-level draft genome for D. kaki was generated in this study, which will contribute to further studies of major economic traits in the genus Diospyros.

Genomic insights into the clonal reproductive Opuntia cochenillifera: mitochondrial and chloroplast genomes of the cochineal cactus for enhanced understanding of structural dynamics and evolutionary implications

Background

The cochineal cactus (Opuntia cochenillifera), notable for its substantial agricultural and industrial applications, predominantly undergoes clonal reproduction, which presents significant challenges in breeding and germplasm innovation. Recent developments in mitochondrial genome engineering offer promising avenues for introducing heritable mutations, potentially facilitating selective sexual reproduction through the creation of cytoplasmic male sterile genotypes. However, the lack of comprehensive mitochondrial genome information for Opuntia species hinders these efforts. 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 the O. cochenillifera using DNBSEQ and Nanopore platforms. The mitochondrial genome was then assembled using a hybrid assembly strategy using Unicycler software. We found that the mitochondrial genome of O. cochenillifera has a length of 1,156,235 bp, a GC content of 43.06%, and contains 54 unique protein-coding genes and 346 simple repeats. Comparative genomic analysis revealed 48 homologous fragments shared between mitochondrial and chloroplast genomes, with a total length of 47,935 bp. Additionally, the comparison of mitochondrial genomes from four Cactaceae species highlighted their dynamic nature and frequent mitogenomic reorganizations.

Conclusion

Our study provides a new perspective on the evolution of the organelle genome and its potential application in genetic breeding. These findings offer valuable insights into the mitochondrial genetics of Cactaceae, potentially facilitating future research and breeding programs aimed at enhancing the genetic diversity and adaptability of O. cochenillifera by leveraging its unique mitochondrial genome characteristics.

Characterizing complete mitochondrial genome of Aquilegia amurensis and its evolutionary implications

BACKGROUND

Aquilegia is a model system for studying the evolution of adaptive radiation. However, very few studies have been conducted on the Aquilegia mitochondrial genome. Since mitochondria play a key role in plant adaptation to abiotic stress, analyzing the mitochondrial genome may provide a new perspective for understanding adaptive evolution.

RESULTS

The Aquilegia amurensis mitochondrial genome was characterized by a circular chromosome and two linear chromosomes, with a total length of 538,736 bp; the genes included 33 protein-coding genes, 24 transfer RNA (tRNA) genes and 3 ribosomal RNA (rRNA) genes. We subsequently conducted a phylogenetic analysis based on single nucleotide polymorphisms (SNPs) in the mitochondrial genomes of 18 Aquilegia species, which were roughly divided into two clades: the European-Asian clade and the North American clade. Moreover, the genes mttB and rpl5 were shown to be positively selected in European-Asian species, and they may help European and Asian species adapt to environmental changes.

CONCLUSIONS

In this study, we assembled and annotated the first mitochondrial genome of the adaptive evolution model plant Aquilegia. The subsequent analysis provided us with a basis for further molecular studies on Aquilegia mitochondrial genomes and valuable information on adaptive evolution in Aquilegia.

High-Quality Assembly and Analysis of the Complete Mitogenomes of German Chamomile (Matricaria recutita) and Roman Chamomile (Chamaemelum nobile)

German chamomile (Matricaria chamomilla L.) and Roman chamomile (Chamaemelum nobile) are the two well-known chamomile species from the Asteraceae family. Owing to their essential oils and higher medicinal value, these have been cultivated widely across Europe, Northwest Asia, North America, and Africa. Regarding medicinal applications, German chamomile is the most commonly utilized variety and is frequently recognized as the "star among medicinal species". The insufficient availability of genomic resources may negatively impact the progression of chamomile industrialization. Chamomile's mitochondrial genome is lacking in extensive empirical research. In this study, we achieved the successful sequencing and assembly of the complete mitochondrial genome of M. chamomilla and C. nobile for the first time. An analysis was conducted on codon usage, sequence repeats within the mitochondrial genome of M. chamomilla and C. nobile. The phylogenetic analysis revealed a consistent positioning of M. chamomilla and C. nobile branches within both mitochondrial and plastid-sequence-based phylogenetic trees. Furthermore, the phylogenetic analysis also showed a close relationship between M. chamomilla and C. nobile within the clade comprising species from the Asteraceae family. The results of our analyses provide valuable resources for evolutionary research and molecular barcoding in chamomile.

Multi-Chromosomal mitochondrial genome of medicinal plant Acorus tatarinowii (Acoraceae): Firstly reported from Acorales Order

Acorus tatarinowii Schott (A. tatarinowii), a well-known traditional Chinese medicinal plant renowned for its high medicinal value, but its mitochondrial genome (mitogenome) is still unexplored. In this study, we meticulously assembled the complete mitochondrial genome of A. tatarinowii using a combination of Illumina short reads and Oxford Nanopore long reads. Our findings revealed that A. tatarinowii possesses a complex chromosomal structural mitogenome, comprising two linear chromosomes and seven circular chromosomes. This mitogenome spans 1.81 Mb in length with a GC content of 38.29 %. Notably, it contained 24 unique mitochondrial core genes, seven unique variable genes, 17 tRNA genes, and three rRNA genes. Analyses of codon usage, most protein-coding genes (PCGs) exhibited a common codon usage preference, with RSCU values greater than 1, and the codon with the highest RSCU value was UAA(End, 1.90). We conducted a thorough analysis of repeat sequences, the distribution of repetitive sequences in nine mitochondrial chromosomes showed distinct patterns. Moreover, we identified 82 and 12 homologous fragments by comparing the sequences of chloroplast and nuclear genomes to the A. tatarinowii mitogenome, respectively. Lastly, We predicted a total of 234 potential RNA editing sites in 28 unique PCGs and discovered that the nad4 gene has been edited the most often, at 26 times. Our results contribute to the enrichment of mitochondrial genome resources for Acoraceae, and the mitogenome also can be used as a reference for other species.

Remarkable mitochondrial genome heterogeneity in Meniocus linifolius (Brassicaceae)

KEY MESSAGE

Detailed analyses of 16 genomes identified a remarkable acceleration of mutation rate, hence mitochondrial sequence and structural heterogeneity, in Meniocus linifolius (Brassicaceae). The powerhouse, mitochondria, in plants feature high levels of structural variation, while the encoded genes are normally conserved. However, the substitution rates and spectra of mitochondria DNA within the Brassicaceae, a family with substantial scientific and economic importance, have not been adequately deciphered. Here, by analyzing three newly assembled and 13 known mitochondrial genomes (mitogenomes), we report the highly variable genome structure and mutation rates in Brassicaceae. The genome sizes and GC contents are 196,604 bp and 46.83%, 288,122 bp and 44.79%, and 287,054 bp and 44.93%, for Meniocus linifolius (Mli), Crucihimalaya lasiocarpa (Cla), and Lepidium sativum (Lsa), respectively. In total, 29, 33, and 34 protein-coding genes (PCGs) and 14, 18, and 18 tRNAs are annotated for Mli, Cla, and Lsa, respectively, while all mitogenomes contain one complete circular molecule with three rRNAs and abundant RNA editing sites. The Mli mitogenome features four conformations likely mediated by the two pairs of long repeats, while at the same time seems to have an unusual evolutionary history due to higher GC content, loss of more genes and sequences, but having more repeats and plastid DNA insertions. Corroborating with these, an ambiguous phylogenetic position with long branch length and elevated synonymous substitution rate in nearly all PCGs are observed for Mli. Taken together, our results reveal a high level of mitogenome heterogeneity at the family level and provide valuable resources for further understanding the evolutionary pattern of organelle genomes in Brassicaceae.

Comparative mitochondrial genome analysis reveals a candidate ORF for cytoplasmic male sterility in tropical onion

Cytoplasmic male sterility (CMS) has been widely exploited for hybrid seed production in onions (Allium cepa L.). In contrast to long-day onion cultivars, short-day onion has not yet been investigated for mitochondrial genome structure and DNA rearrangements associated with CMS activity. Here, we report the 3,16,321 bp complete circular mitochondrial genome of tropical onion CMS line (97A). Due to the substantial number of repetitive regions, the assembled mitochondrial genome of maintainer line (97B) remained linear with 15 scaffolds. Additionally, 13 and 20 chloroplast-derived fragments with a size ranging from 143 to 13,984 bp and 153-17,725 bp were identified in the 97A and 97B genomes, respectively. Genome annotation revealed 24 core protein-coding genes along with 24 and 28 tRNA genes in the mitochondrial genomes of 97A and 97B, respectively. Furthermore, comparative genome analysis of the 97A and 97B mitochondrial genomes showed that gene content was almost similar except for the chimeric ORF725 gene which is the extended form of the COX1 gene.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03850-2.

Complete mitochondrial genome of the endangered Prunus pedunculata (Prunoideae, Rosaceae) in China: characterization and phylogenetic analysis

Introduction

Prunus pedunculata (Prunoideae: Rosaceae), a relic shrub with strong resistance and multiple application values, is endangered in China. Extensive research had been devoted to gene expression, molecular markers, plastid genome analysis, and genetic background investigations of P. pedunculata. However, the mitochondrial genome of this species has not been systematically described, owing to the complexity of the plant mitogenome.

Methods

In the present research, the complete mitochondrial genome of P. pedunculata was assembled, annotated, and characterized. The genomic features, gene content and repetitive sequences were analyzed. The genomic variation and phylogenetic analysis have been extensively enumerated.

Results and discussion

The P. pedunculata mitogenome is a circular molecule with a total length of 405,855 bp and a GC content of 45.63%, which are the smallest size and highest GC content among the known Prunus mitochondrial genomes. The mitogenome of P. pedunculata encodes 62 genes, including 34 unique protein-coding genes (PCGs, excluding three possible pseudogenes), three ribosomal RNA genes, and 19 transfer RNA genes. The mitogenome is rich in repetitive sequences, counting 112 simple sequence repeats, 15 tandem repeats, and 50 interspersed repetitive sequences, with a total repeat length of 11,793 bp, accounting for 2.91% of the complete genome. Leucine (Leu) was a predominant amino acid in PCGs, with a frequency of 10.67%, whereas cysteine (Cys) and tryptophan (Trp) were the least adopted. The most frequently used codon was UUU (Phe), with a relative synonymous codon usage (RSCU) value of 1.12. Selective pressure was calculated based on 20 shared PCGs in the mitogenomes of the 32 species, most of which were subjected to purifying selection (Ka/Ks < 1), whereas ccmC and ccmFn underwent positive selection. A total of 262 potential RNA editing sites in 26 PCGs were identified. Furthermore, 56 chloroplast-derived fragments were ascertained in the mitogenome, ranging from 30 to 858 bp, and were mainly located across IGS (intergenic spacer) regions or rRNA genes. These findings verify the occurrence of intracellular gene transfer events from the chloroplast to the mitochondria. Furthermore, the phylogenetic relationship of P. pedunculata was supported by the mitogenome data of 30 other taxa of the Rosaceae family. Understanding the mitochondrial genome characteristics of P. pedunculata is of great importance to promote comprehension of its genetic background and this study provides a basis for the genetic breeding of Prunus.

Intrageneric structural variation in organelle genomes from the genus Dystaenia (Apiaceae): genome rearrangement and mitochondrion-to-plastid DNA transfer

Introduction

During plant evolution, intracellular DNA transfer (IDT) occurs not only from organelles to the nucleus but also between organelles. To further comprehend these events, both organelle genomes and transcriptomes are needed.

Methods

In this study, we constructed organelle genomes and transcriptomes for two Dystaenia species and described their dynamic IDTs between their nuclear and mitochondrial genomes, or plastid and mitochondrial genomes (plastome and mitogenome).

Results and Discussion

We identified the putative functional transfers of the mitochondrial genes 5' rpl2, rps10, rps14, rps19, and sdh3 to the nucleus in both Dystaenia species and detected two transcripts for the rpl2 and sdh3 genes. Additional transcriptomes from the Apicaceae species also provided evidence for the transfers and duplications of these mitochondrial genes, showing lineage-specific patterns. Intrageneric variations of the IDT were found between the Dystaenia organelle genomes. Recurrent plastid-to-mitochondrion DNA transfer events were only identified in the D. takeshimana mitogenome, and a pair of mitochondrial DNAs of plastid origin (MIPTs) may generate minor alternative isoforms. We only found a mitochondrion-to-plastid DNA transfer event in the D. ibukiensis plastome. This event may be linked to inverted repeat boundary shifts in its plastome. We inferred that the insertion region involved an MIPT that had already acquired a plastid sequence in its mitogenome via IDT. We propose that the MIPT acts as a homologous region pairing between the donor and recipient sequences. Our results provide insight into the evolution of organelle genomes across the family Apiaceae.

Complete mitogenome assembly of Selenicereus monacanthus revealed its molecular features, genome evolution, and phylogenetic implications

BACKGROUND

Mitochondria are the powerhouse of the cell and are critical for plant growth and development. Pitaya (Selenicereus or Hylocereus) is the most important economic crop in the family Cactaceae and is grown worldwide, however its mitogenome is unreported.

RESULTS

This study assembled the complete mitogenome of the red skin and flesh of pitaya (Selenicereus monacanthus). It is a full-length, 2,290,019 bp circular molecule encoding 59 unique genes that only occupy 2.17% of the entire length. In addition, 4,459 pairs of dispersed repeats (≥ 50 bp) were identified, accounting for 84.78% of the total length, and three repeats (394,588, 124,827, and 13,437 bp) mediating genomic recombination were identified by long read mapping and Sanger sequencing. RNA editing events were identified in all 32 protein-coding genes (PCGs), among which four sites (nad1-2, nad4L-2, atp9-copy3-223, and ccmFC-1309) were associated with the initiation or termination of PCGs. Seventy-eight homologous fragments of the chloroplast genome were identified in the mitogenome, the longest having 4,523 bp. In addition, evolutionary analyses suggest that S. monacanthus may have undergone multiple genomic reorganization events during evolution, with the loss of at least nine PCGs (rpl2, rpl10, rps2, rps3, rps10, rps11, rps14, rps19, and sdh3).

CONCLUSIONS

This study revealed the genetic basis of the S. monacanthus mitogenome, and provided a scientific basis for further research on phenotypic traits and germplasm resource development.