Master graph: an essential integrated assembly model for the plant mitogenome based on a graph-based framework

Complete mitochondrial genome of Houttuynia cordata sheds light on the evolution of complex structural features in early angiosperms

BACKGROUND

Houttuynia cordata is a well-known medicinal and edible plant with important economic value. It belongs to the magnoliid clade, a basal lineage of angiosperms. Currently, the limited number of reported mitochondrial genomes from magnoliids severely restricts our understanding of early mitochondrial genome evolution. To address this gap, we present the complete mitochondrial genome of H. cordata.

RESULTS

The H. cordata mitochondrial genome spans 534,194 bp, and predominantly adopts a bicyclic structure among four possible configurations. It encodes 45 protein-coding genes (PCGs) (40 unique), 29 transfer RNAs (21 unique), and 3 ribosomal RNAs. Repeat analysis revealed significant differences in tandem repeats, dispersed repeats, and simple sequence repeat (SSR) types and densities between the two chromosomes, with Chromosome 1 exhibiting higher repeat density and diversity. Additionally, the H. cordata mitochondrial genome harbors 18 potential mitochondrial plastid DNAs (MTPTs), most of which originate from the inverted repeats (IR) regions of the plastid genome. Variant detection at the individual-level revealed significant heteroplasmy across most regions, while large repeat regions and MTPTs regions exhibited relative conservation. Phylogenetic and collinearity analyses further indicated substantial differences between the mitochondrial genome structure of H. cordata and 15 other magnoliid species, while conserved regions were maintained.

CONCLUSION

This study elucidates the unique evolutionary characteristics of the H. cordata mitochondrial genome and provides valuable insights into its genetics and evolution. Comparative analysis with other magnoliid species reveals both dynamic evolutionary changes and the conservation of mitochondrial genome structure and gene content in early angiosperms.

The characterisation of the complete mitochondrial genome of Polygonatum Kingianum reveals recombination mediated by repeats associated with DNA replication

Polygonatum kingianum Coll et Hemsl (P. kingianum) is recognized as one of the most valuable medicinal species, demonstrating a diverse range of bioactivities, including antidiabetic effects, reduction of blood lipids, anti-tumor properties, and hyperglycemic activity. However, the mitochondrial genome of P. kingianum has not yet been reported. Therefore, we sequenced and assembled the mitochondrial DNA (mtDNA) of P. kingianum. The results indicate that the mitochondrial genome possesses a multi-branched conformational structure. Further annotation of the mitochondrial genome reveals that it is 647,110 bp in length, containing a total of 35 distinct protein-coding genes, 22 tRNA genes, and 4 rRNA genes. The analysis of gene loss revealed that a total of seven genes were absent in the mitochondrial genome of P. kingianum. This includes six ribosomal protein genes (rpl2, rpl6, rpl10, rps7, rps8, rps11) and two succinate dehydrogenase genes (sdh3 and sdh4), as well as one pseudogene. Furthermore, Analysis of mitochondrial genome coding sequences of ten closely related species revealed a total of 24 shared genes. The Ka/Ks values of all shared coding protein genes are less than 1, indicating that these genes have undergone purification selection during evolution and their protein functions are relatively stable. In addition, codon studies on P. kingianum mitochondria showed 29 high-frequency codons (RSCU > 1). The preference for codon usage is influenced by mutations and natural selection, but mainly determined by natural selection. Furthermore, we identified 31 homologous fragments spanning both chloroplast and mitochondrial genomes and 585 C-to-U RNA editing sites were predicted in mitochondrial PCGs. The phylogenetic tree established a close relationship between P. kingianum and the 27 closely related species. To sum up, this study will contribute to the application of population genetics and evolutionary research in the genus Polygonatum and other genera in the Asparagales family.

The complete mitochondrial genome of Cenchrus fungigraminus indicates structural dynamics and sequence divergences in Poaceae family

Background

Cenchrus fungigraminus is a perennial grass species, which is used in the cultivation of edible fungi, helping to resolve land use conflicts between edible fungus cultivation and forest production. It can be also utilized for bioethanol production and fiberboard manufacturing. Currently, no report was released on the mitochondrial genome. Thus, the C. fungigraminus mitochondrial genome was assembled and annotated utilizing a combination of Illumina and PacBio sequencing data.

Results

The assembled genome measures 404,398 bp and exhibits a linear structure. The GC content is 43. 98%, and 64 characteristic genes were identified, containing 34 protein-coding, six rRNA, 22 tRNA and two pseudo-genes, with 17 genes containing 30 introns. The mitochondrial genome codon usage displayed a significant A/T preference. A total of 235 dispersed, 93 simple sequence, and 38 tandem repeats were detected. Furthermore, 34 sequence fragments (17,634 bp) were transferred from the chloroplast to the mitochondria, accounting for 4.36% of the complete mitochondrial genome. The gene rps2 was identified as having the highest polymorphism. Additionally, these genes (ccmB, nad2, and rps7) were influenced by positive selection. Among the reported Poaceae mitochondrial genomes, the genome sizes range from 381,065 to 704,100 bp, with GC content ranging from 42.27% to 44.43%. Three genes (sdh3, rps10, and rpl6) are missing in the Poaceae mitochondrial genome, indicating significant rearrangements in the overall genome structure and low homology. It was shown that C. fungigraminus is most closely related to Saccharum officinarum within the Poaceae family by Phylogenetic analysis. The interspecific Ka/Ks comparison indicated that 28 protein-coding genes have a Ka/Ks ratio < 1.

Conclusion

Therefore, the mitochondrial genome data of C. fungigraminus was enriched and a theoretical basis was provided for germplasm evaluation and new variety utilization.

De-novo assembly and comparative analysis of the complete mitogenome of traditional Chinese medicine Strobilanthes sarcorrhiza

BACKGROUND

Strobilanthes sarcorrhiza is a traditional medicinal plant known for its heat-clearing and kidney-nourishing properties. While its plastid genome has been reported, there is a scarcity of genetic information regarding its mitogenome, leading to unclear phylogenetic relationships. We sequenced and assembled the complete its mitogenome and conducted a series of genetic analyses in conjunction with the plastid genome to gain a better understanding of the species' genetic background.

RESULTS

The mitogenome comprised a linear structure spanning 617,134 bp. It included 35 protein-coding genes (PCGs), 19 transfer RNAs (tRNAs), and 3 ribosomal RNAs (rRNAs) that have been annotated. Additionally, 122 simple sequence repeats (SSRs) and 25 tandem repeats were identified. A total of 1,482 pairs of dispersed repeats were detected, which account for 17.58% of the entire mitogenome. Furthermore, 37 migration fragments between the mitochondrial and plastid genomes were discovered, consisting of 5 complete PCGs, 7 tRNAs, and 1 rRNA. Based on the analysis of 38 mitogenomes and 46 plastid genomes, the evolutionary relationship and phylogenetic position of S. sarcorrhiza were elucidated.

CONCLUSIONS

This study has, for the first time, provided insights into the mitochondrial genomic characteristics of S. sarcorrhiza and clarified its phylogenetic position. These findings offered significant insights for the future identification and classification of this genus, as well as for the genetic breeding of medicinal plants.

Organelle genomes of two Scaevola species, S. taccada and S. hainanensis, provide new insights into evolutionary divergence between Scaevola and its related species

Chloroplast and mitochondrial genomes harbor crucial information that can be utilized for elucidating plant evolution and environmental adaptation. The organellar genomic characteristics of Goodeniaceae, a sister family to Asteraceae, remain unexplored. Here, using a combination of short-read and long-read sequencing technologies, we successfully assembled the complete organellar genomes of two Goodeniaceae species native to China, Scaevola taccada and S. hainanensis. Chloroplast genome collinearity analysis revealed that Scaevola expanded its genome length through inverted repeat expansion and large single copy fragment duplication, resulting in 181,022 bp (S. taccada) and 182,726 bp (S. hainanensis), ~30 kb increase compared to its related species. Mitochondrial genomes of two Scaevola species exhibit multi-ring topology, forming dual mitochondrial chromosomes of 314,251 bp (S. taccada) and 276,175 bp (S. hainanensis). Sequence variation analysis demonstrated substantial chloroplast sequence divergence (Pi = 0.45) and an increase in gene copy number within the genus. Relative synonymous codon usage (RSCU) analysis revealed that Scaevola chloroplast has a higher bias for A/U-ending codons than mitochondria, with chloroplasts RSCU values ranging from 0.32 to 1.94, whereas mitochondrial RSCU values ranging from 0.38 to 1.62. Phylogenetic analyses support the monophyly of the Asteraceae-Goodeniaceae sister group, whereas the extended evolutionary branches of Scaevola, coupled with mitochondrial collinearity analysis, indicate rapid organellar genome evolution of Scaevola. Organellar-nuclear horizontal gene transfer analysis identified specific increased in the copy numbers of photosynthesis-related genes and chloroplast-nuclear transfer events in S. taccada. Our study not only provides insights for understanding environmental adaptation mechanisms of coastal plants, but also contributes to elucidating organellar genome evolution in Scaevola and Goodeniaceae.

Unequally Abundant Chromosomes and Unusual Collections of Transferred Sequences Characterize Mitochondrial Genomes of Gastrodia (Orchidaceae), One of the Largest Mycoheterotrophic Plant Genera

The mystery of genomic alternations in heterotrophic plants is among the most intriguing in evolutionary biology. Compared to plastid genomes (plastomes) with parallel size reduction and gene loss, mitochondrial genome (mitogenome) variation in heterotrophic plants remains underexplored in many aspects. To further unravel the evolutionary outcomes of heterotrophy, we present a comparative mitogenomic study with 13 de novo assemblies of Gastrodia (Orchidaceae), one of the largest fully mycoheterotrophic plant genera, and its relatives. Analyzed Gastrodia mitogenomes range from 0.56 to 2.1 Mb, each consisting of numerous, unequally abundant chromosomes or contigs. Size variation might have evolved through chromosome rearrangements followed by stochastic loss of "dispensable" chromosomes, with deletion-biased mutations. The discovery of a hyper-abundant (∼15 times intragenomic average) chromosome in two assemblies represents the hitherto most extreme copy number variation in any mitogenomes, with similar architectures discovered in two metazoan lineages. Transferred sequence contents highlight asymmetric evolutionary consequences of heterotrophy: despite drastically reduced intracellular plastome transfers convergent across heterotrophic plants, their rarity of horizontally acquired sequences sharply contrasts parasitic plants, where massive transfers from their hosts prevail. Rates of sequence evolution are markedly elevated but not explained by copy number variation, extending prior findings of accelerated molecular evolution from parasitic to heterotrophic plants. Putative evolutionary scenarios for these mitogenomic convergence and divergence fit well with the common (e.g. plastome contraction) and specific (e.g. host identity) aspects of the two heterotrophic types. These idiosyncratic mycoheterotrophs expand known architectural variability of plant mitogenomes and provide mechanistic insights into their content and size variation.

Genomic variation and evolutionary patterns in organelle genomes between annual and perennial Glycine species

BACKGROUND

The complexity of structural variations and long stretches of repetitive DNA make the analysis of plant mitochondrial genomes (mitogenomes) exceptionally challenging. A thorough investigation of plant mitogenomes is essential for uncovering the evolutionary processes of plant organelles and optimizing traits related to plant cellular metabolism. The genus Glycine includes groups with both perennial and annual life strategies, making it an ideal subject for studying the complexity and variations of plant mitogenomes during evolution across different life strategies.

RESULTS

Here, we assembled 20 complete mitochondrial and plastid genomes of Glycine accessions, including both annual and perennial species using the latest organelle genome assembly tool. Significant structural variations and differences in tRNA content were observed in the mitogenomes between the two life-history strategy subgenera, while protein-coding genes and rRNAs content were highly conserved. Distinct patterns of nuclear plastid DNAs and nuclear mitochondrial DNAs (NUPTs/NUMTs) were uncovered among annual and perennial species. Genes residing in NUMTs (NUMGs) showed a substantial presence in Glycine accessions, with annual soybeans exhibiting a higher proportion of protein-coding genes fully integrated into the nuclear genome. Phylogenetic analysis indicated a closely related evolutionary trajectory between mitochondrial and nuclear genomes in Glycine, providing supplementary evidence relevant to the evolutionary history of Glycine.

CONCLUSIONS

This study showed the structural variations and evolutionary patterns of mitochondrial genomes between annual and perennial Glycine species. These findings contribute to our understanding of plant organelle complexity, variation and history of intracellular genomic integration.

The Taihangia mitogenome provides new insights into its adaptation and organelle genome evolution in Rosaceae

MAIN CONCLUSION

We present the first Taihangia mitogenome, uncovering frequent rearrangements and significant length variation in Rosaceae, likely driven by hybridization and repeat content, alongside widespread mito-chloroplast phylogenetic conflicts. Taihangia, an ancient and endangered monotypic genus within the subfamily Rosoideae of the family Rosaceae, is endemic to cliffs and serves as an ideal material for studying the adaptations of cliff-dwelling plants and the evolutionary processes of the Rosaceae family. In this study, the mitogenome and plastome of T. rupestris var. ciliata were assembled, with lengths of 265,633 bp and 155,467 bp, both exhibiting typical circular structures. Positive selection was detected in the nad4L and sdh4 genes, likely playing a role in adaptation to harsh environments. Comparative genomic analysis indicated that repetitive sequences are likely the main contributors to genome size variation in Rosaceae and also influence horizontal gene transfer between organelle genomes. In T. rupestris var. ciliata, 20 mitochondrial plastid DNA sequences were identified, including 16 complete plastid genes. Moreover, frequent rearrangements were observed in the non-coding regions of mitogenome within the subfamily Rosoideae, potentially linked to the complex evolutionary history and the presence of repetitive sequences. In contrast, coding regions remained highly conserved (over 83% similarity) to maintain essential mitochondrial functions. Phylogenomic analysis of the two organelle genomes revealed conflicts in the phylogenetic relationships within Rosaceae, potentially due to the inconsistent mutation rates and frequent hybridization events in the evolutionary history of the family. In conclusion, the organelle genome analysis of Taihangia provides crucial genomic resources for understanding the evolution and adaptation of Rosaceae species.

Evolutionary dynamics of mitochondrial genomes and intracellular transfers among diploid and allopolyploid cotton species

BACKGROUND

Plant mitochondrial genomes (mitogenomes) exhibit extensive structural variation yet extremely low nucleotide mutation rates, phenomena that remain only partially understood. The genus Gossypium, a globally important source of cotton, offers a wealth of long-read sequencing resources to explore mitogenome and plastome variation and dynamics accompanying the evolutionary divergence of its approximately 50 diploid and allopolyploid species.

RESULTS

Here, we assembled 19 mitogenomes from Gossypium species, representing all genome groups (diploids A through G, K, and the allopolyploids AD) based on a uniformly applied strategy. A graph-based mitogenome assembly method revealed more alternative structural conformations than previously recognized, some of which confirmed the mitogenome structure reported in earlier studies on cotton. Using long-read data, we quantified alternative conformations mediated by recombination events between repeats, and phylogenetically informative structural variants were noted. Nucleotide substitution rate comparisons between coding and non-coding regions revealed low mutation rates across the entire mitogenome. Genome-wide mapping of nuclear organellar DNA transfers (NUOTs) in Gossypium revealed a nonrandom distribution of transfers in the nuclear genome. In cotton, the fate of NUOT events varied, with mitochondrion-to-nucleus transfer (NUMT) predominantly retained as short fragments in the nuclear genome, with more plastid sequences integrated into the nucleus. Phylogenetic relationships inferred using different data sets highlighted distinct evolutionary histories among these cellular compartments, providing ancillary evidence relevant to the evolutionary history of Gossypium.

CONCLUSIONS

A comprehensive analysis of organellar genome variation demonstrates complex structural variation and low mutation rates across the entire mitogenome and reveals the history of organellar genome transfer among the three genomes throughout the cotton genus. The findings enhance our general understanding of mitogenome evolution, comparative organellar and nuclear evolutionary rates, and the history of inter-compartment genomic integration.

TIPPo: A User-Friendly Tool for De Novo Assembly of Organellar Genomes with High-Fidelity Data

Plant cells have two major organelles with their own genomes: chloroplasts and mitochondria. While chloroplast genomes tend to be structurally conserved, the mitochondrial genomes of plants, which are much larger than those of animals, are characterized by complex structural variation. We introduce TIPPo, a user-friendly, reference-free assembly tool that uses PacBio high-fidelity long-read data and that does not rely on genomes from related species or nuclear genome information for the assembly of organellar genomes. TIPPo employs a deep learning model for initial read classification and leverages k-mer counting for further refinement, significantly reducing the impact of nuclear insertions of organellar DNA on the assembly process. We used TIPPo to completely assemble a set of 54 complete chloroplast genomes. No other tool was able to completely assemble this set. TIPPo is comparable with PMAT in assembling mitochondrial genomes from most species but does achieve even higher completeness for several species. We also used the assembled organelle genomes to identify instances of nuclear plastid DNA (NUPTs) and nuclear mitochondrial DNA (NUMTs) insertions. The cumulative length of NUPTs/NUMTs positively correlates with the size of the nuclear genome, suggesting that insertions occur stochastically. NUPTs/NUMTs show predominantly C:G to T:A changes, with the mutated cytosines typically found in CG and CHG contexts, suggesting that degradation of NUPT and NUMT sequences is driven by the known elevated mutation rate of methylated cytosines. Small interfering RNA loci are enriched in NUPTs and NUMTs, consistent with the RdDM pathway mediating DNA methylation in these sequences.

Unveiling the intricate structural variability induced by repeat-mediated recombination in the complete mitochondrial genome of Cuscuta gronovii Willd

Cuscuta gronovii Willd., a member of the Convolvulaceae family, is noted for its potential medicinal and nutritional benefits. In this study, we utilized a combination of Illumina and Oxford Nanopore sequencing technologies to successfully assemble the complete circular mitochondrial genome (mitogenome) of C. gronovii. The mitogenome, spanning 304,467 base pairs, includes 54 genes: 33 protein-coding genes, three ribosomal RNA (rRNA) genes, and 18 transfer RNA (tRNA) genes. Beyond its primary circular structure, we discovered and validated several alternative genomic conformations, driven by five specific repeat sequences. Three inverted repeats were found to initiate rearrangements, resulting in the creation of seven distinct chromosomal structures, while two direct repeats split a larger molecule into two subgenomic entities. We also mapped 421 RNA editing sites across the protein-coding sequences, influencing 33 protein-coding genes with varying distribution, particularly noting high frequencies in the nad4 and ccmB genes. Sixteen of these RNA editing sites were experimentally validated through PCR amplification and Sanger sequencing, confirming their presence with 100 % accuracy. This research not only introduces the first mitochondrial genome of C. gronovii but also highlights its complex conformational variability induced by repeat-mediated recombination, providing a valuable genomic resource for further molecular breeding efforts and phylogenetic evolution within the genus Cuscuta.

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.

Structural analysis of the mitochondrial genome of Santalum album reveals a complex branched configuration

BACKGROUND

Santalum album L. is an evergreen tree which is mainly distributes throughout tropical and temperate regions. And it has a great medicinal and economic value.

RESULTS

In this study, the complete mitochondrial genome of S. album were assembled and annotated, which could be descried by a complex branched structure consisting of three contigs. The lengths of these three contigs are 165,122 bp, 93,430 bp and 92,491 bp. We annotated 34 genes coding for proteins (PCGs), 26 tRNA genes, and 4 rRNA genes. The analysis of repeated elements shows that there are 89 SSRs and 242 pairs of dispersed repeats in S. album mitochondrial genome. Also we found 20 MTPTs among the chloroplast and mitochondria. The 20 MTPTs sequences span a combined length of 22,353 bp, making up 15.52 % of the plastome, 6.37 % of the mitochondrial genome. Additionally, by using the Deepred-mt tool, we found 628 RNA editing sites in 34 PCGs. Moreover, significant genomic rearrangement is observed between S. album and its associated mitochondrial genomes. Finally, based on mitochondrial genome PCGs, we deduced the phylogenetic ties between S. album and other angiosperms.

CONCLUSIONS

We reported the mitochondrial genome from Santalales for the first time, which provides a crucial genetic resource for our study of the evolution of mitochondrial genome.

Morphological Structure Identification, Comparative Mitochondrial Genomics and Population Genetic Analysis toward Exploring Interspecific Variations and Phylogenetic Implications of Malus baccata 'ZA' and Other Species

Malus baccata, a valuable germplasm resource in the genus Malus, is indigenous to China and widely distributed. However, little is known about the lineage composition and genetic basis of 'ZA', a mutant type of M. baccata. In this study, we compared the differences between 'ZA' and wild type from the perspective of morphology and ultrastructure and analyzed their chloroplast pigment content based on biochemical methods. Further, the complete mitogenome of M. baccata 'ZA' was assembled and obtained by next-generation sequencing. Subsequently, its molecular characteristics were analyzed using Geneious, MISA-web, and CodonW toolkits. Furthermore, by examining 106 Malus germplasms and 42 Rosaceae species, we deduced and elucidated the evolutionary position of M. baccata 'ZA', as well as interspecific variations among different individuals. In comparison, the total length of the 'ZA' mitogenome (GC content: 45.4%) is 374,023 bp, which is approximately 2.33 times larger than the size (160,202 bp) of the plastome (GC: 36.5%). The collinear analysis results revealed abundant repeats and genome rearrangements occurring between different Malus species. Additionally, we identified 14 plastid-driven fragment transfer events. A total of 54 genes have been annotated in the 'ZA' mitogenome, including 35 protein-coding genes, 16 tRNAs, and three rRNAs. By calculating nucleotide polymorphisms and selection pressure for 24 shared core mitochondrial CDSs from 42 Rosaceae species (including 'ZA'), we observed that the nad3 gene exhibited minimal variation, while nad4L appeared to be evolving rapidly. Population genetics analysis detected a total of 1578 high-quality variants (1424 SNPs, 60 insertions, and 94 deletions; variation rate: 1/237) among samples from 106 Malus individuals. Furthermore, by constructing phylogenetic trees based on both Malus and Rosaceae taxa datasets, it was preliminarily demonstrated that 'ZA' is closely related to M. baccata, M. sieversii, and other proximate species in terms of evolution. The sequencing data obtained in this study, along with our findings, contribute to expanding the mitogenomic resources available for Rosaceae research. They also hold reference significance for molecular identification studies as well as conservation and breeding efforts focused on excellent germplasms.

Graph-based mitochondrial genomes of three foundation species in the Saccharum genus

KEY MESSAGE

We reported the graph-based mitochondrial genomes of three foundation species (Saccharum spontaneum, S. robustum and S. officinarum) for the first time. The results revealed pan-structural variation and evolutionary processes in the mitochondrial genomes within Saccharum. Saccharum belongs to the Andropogoneae, and cultivars species in Saccharum contribute nearly 80% of sugar production in the world. To explore the genomic studies in Saccharum, we assembled 15 complete mitochondrial genomes (mitogenome) of three foundation species (Saccharum spontaneum, S. robustum and S. officinarum) using Illumina and Oxford Nanopore Technologies sequencing data. The mitogenomes of the three species were divided into a total of eight types based on contig numbers and linkages. All mitogenomes in the three species encoded 51 unique genes, including 32 protein-coding, 3 ribosomal RNA (rRNA) and 16 transfer RNA (tRNA) genes. The existence of long and short-repeat-mediated recombinations in the mitogenome of S. officinarum and S. robustum was revealed and confirmed through PCR validation. Furthermore, employing comparative genomics and phylogenetic analyses of the organelle genomes, we unveiled the evolutionary relationships and history of the major interspecific lineages in Saccharum genus. Phylogenetic analyses of homologous fragments between S. officinarum and S. robustum showed that S. officinarum and S. robustum are phylogenetically distinct and that they were likely parallel rather than domesticated. The variations between ancient (S. sinense and S. barberi) and modern cultivated species (S. hybrid) possibly resulted from hybridization involving different S. officinarum accessions. Lastly, this project reported the first graph-based mitogenomes of three Saccharum species, and a systematic comparison of the structural organization, evolutionary processes, and pan-structural variation of the Saccharum mitogenomes revealed the differential features of the Saccharum mitogenomes.

The complete mitochondrial genome of Pontederia crassipes: using HiFi reads to investigate genome recombination and gene transfer from chloroplast genome

Water hyacinth (Pontederia crassipes Mart.) is a monocotyledonous aquatic plant renowned for its rapid growth, extensive proliferation, biological invasiveness, and ecological resilience to variations in pH, nutrients, and temperature. The International Union for Conservation of Nature (IUCN) has listed P. crassipes among the top 100 invasive species. However, comprehensive genomic information, particularly concerning its mitochondrial genome (mitogenome), remains surprisingly limited. In this study, the complete mitogenome of P. crassipes was analyzed using bioinformatics approaches. The mitogenome is 399,263 bp long and contains 38 protein-coding genes (PCGs), 24 tRNA genes, and 3 rRNA genes. Sequence analysis revealed that the complete mitogenome of the species contains 3,289 dispersed repeats, and 765 RNA editing sites in protein-coding genes. The P. crassipes mitogenome possessed un-conserved structures, including extensive sequence transfer between its chloroplasts and mitochondria. Our study on the mitogenome of P. crassipes offers critical insights into its evolutionary patterns and phylogenetic relationships with related taxa. This research enhances our understanding of this invasive species, known for its significant biomass and rapid overgrowth in aquatic environments.

Insights into structure, codon usage, repeats, and RNA editing of the complete mitochondrial genome of Perilla frutescens (Lamiaceae)

Perilla frutescens (L.) Britton, a member of the Lamiaceae family, stands out as a versatile plant highly valued for its unique aroma and medicinal properties. Additionally, P. frutescens seeds are rich in Îś-linolenic acid, holding substantial economic importance. While the nuclear and chloroplast genomes of P. frutescens have already been documented, the complete mitochondrial genome sequence remains unreported. To this end, the sequencing, annotation, and assembly of the entire Mitochondrial genome of P. frutescens were hereby conducted using a combination of Illumina and PacBio data. The assembled P. frutescens mitochondrial genome spanned 299,551 bp and exhibited a typical circular structure, involving a GC content of 45.23%. Within the genome, a total of 59 unique genes were identified, encompassing 37 protein-coding genes, 20 tRNA genes, and 2 rRNA genes. Additionally, 18 introns were observed in 8 protein-coding genes. Notably, the codons of the P. frutescens mitochondrial genome displayed a notable A/T bias. The analysis also revealed 293 dispersed repeat sequences, 77 simple sequence repeats (SSRs), and 6 tandem repeat sequences. Moreover, RNA editing sites preferentially produced leucine at amino acid editing sites. Furthermore, 70 sequence fragments (12,680 bp) having been transferred from the chloroplast to the mitochondrial genome were identified, accounting for 4.23% of the entire mitochondrial genome. Phylogenetic analysis indicated that among Lamiaceae plants, P. frutescens is most closely related to Salvia miltiorrhiza and Platostoma chinense. Meanwhile, inter-species Ka/Ks results suggested that Ka/Ks < 1 for 28 PCGs, indicating that these genes were evolving under purifying selection. Overall, this study enriches the mitochondrial genome data for P. frutescens and forges a theoretical foundation for future molecular breeding research.

Analysis of the Rhodomyrtus tomentosa mitochondrial genome: Insights into repeat-mediated recombination and intra-cellular DNA transfer

Plant mitochondrial genomes participate in encoding proteins crucial to the major producers of ATP in the cell and replication and heredity of their own DNA. The sequences and structure of the plant mitochondrial genomes profoundly impact these fundamental processes, and studies of plant mitochondrial genomes are needed. We reported the complete sequences of the Rhodomyrtus tomentosa mitochondrial genome here, totaling 400,482 bp. Nanopore ONT reads and PCR amplification provided evidence for recombination mediated by the eight repeat pairs for the R. tomentosa mitochondrial genome. Thirty-eight genes were identified in the R. tomentosa mitochondrial genome. Comparative analyses of the mitochondrial genome and plastome and PCR amplification suggest that five fragments of mitochondrial plastid DNA were unfunctional sequences resulting from intracellular gene transfer. Phylogenetic analysis based on each and all of the 27 mitochondrial protein-coding genes of nine Myrtales species revealed that R. tomentosa always clustered with other species of Myrtaceae. This study uncovered the enormous complexity of the R. tomentosa mitochondrial genome, the active repeat-mediated recombinations, the presence of mitochondrial plastid DNAs, and the topological incongruence of Myrtales among the single-gene trees.

Comprehensive analysis of the mitochondrial genome of Rehmannia glutinosa: insights into repeat-mediated recombinations and RNA editing-induced stop codon acquisition

Rehmannia glutinosa is an economically significant medicinal plant. Yet, the structure and sequence of its mitochondrial genome has not been published, which plays a crucial role in evolutionary analysis and regulating respiratory-related macromolecule synthesis. In this study, the R. glutinosa mitogenome was sequenced employing a combination of Illumina short reads and Nanopore long reads, with subsequent assembly using a hybrid strategy. We found that the predominant configuration of the R. glutinosa mitogenome comprises two circular chromosomes. The primary structure of the mitogenome encompasses two mitochondrial chromosomes corresponding to the two major configurations, Mac1-1 and Mac1-2. The R. glutinosa mitogenome encoded an angiosperm-typical set of 24 core genes, nine variable genes, three rRNA genes, and 15 tRNA genes. A phylogenetic analysis using the 16 shared protein-coding genes (PCG) yielded a tree consistent with the phylogeny of Lamiales species and two outgroup taxa. Mapping RNA-seq data to the coding sequences (CDS) of the PCGs revealed 507 C-to-U RNA editing sites across 31 PCGs of the R. glutinosa mitogenome. Furthermore, one start codon (nad4L) and two stop codons (rpl10 and atp6) were identified as products of RNA editing events in the R. glutinosa mitogenome.

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.

The complete mitochondrial genome of Aglaia odorata, insights into its genomic structure and RNA editing sites

Aglaia odorata, native to Guangdong, Guangxi, and Hainan provinces in China, has long been utilized as an herbal remedy in ancient China. In this study, we assembled and annotated the complete mitochondrial genome (mitogenome) of A. odorata, which spans a total length of 537,321 bp. Conformation of the A. odorata recombination was verified through PCR experiments and Sanger sequencing. We identified and annotated 35 protein-coding genes (PCGs), 22 tRNA genes, and 3 rRNA genes within the mitogenome. Analysis of repeated elements revealed the presence of 192 SSRs, 29 pairs of tandem repeats, and 333 pairs of dispersed repeats in the A. odorata mitogenome. Additionally, we analyzed codon usage and mitochondrial plastid DNAs (MTPTs). Twelve MTPTs between the plastome and mitogenome of A. odorata were identified, with a combined length of 2,501 bp, accounting for 0.47% of the mitogenome. Furthermore, 359 high-confidence C to U RNA editing sites were predicted on PCGs, and four selected RNA editing sites were specially examined to verify the creation of start and/or stop codons. Extensive genomic rearrangement was observed between A. odorata and related mitogenomes. Phylogenetic analysis based on mitochondrial PCGs were conducted to elucidate the evolutionary relationships between A. odorata and other angiosperms.

Comprehensive Analysis of the Complete Mitochondrial Genome of Rehmannia chingii: An Autotrophic Species in the Orobanchaceae Family

Rehmannia chingii is an important medicinal plant with immense value in scientific research. However, its mitochondrial genome (mitogenome) has not yet been characterized. Herein, based on whole-genome Illumina short reads and PacBio HiFi reads, we obtained the complete mitogenome of R. chingii through a de novo assembly strategy. We carried out comparative genomic analyses and found that, in comparison with the plastid genome (plastome) showing a high degree of structural conservation, the R. chingii mitogenome structure is relatively complex, showing an intricate ring structure with 16 connections, owing to five repetitive sequences. The R. chingii mitogenome was 783,161 bp with a GC content of 44.8% and contained 77 genes, comprising 47 protein-coding genes (CDS), 27 tRNA genes, and 3 rRNA genes. We counted 579 RNA editing events in 47 CDS and 12,828 codons in all CDSs of the R. chingii mitogenome. Furthermore, 24 unique sequence transfer fragments were found between the mitogenome and plastome, comprising 8 mitogenome CDS genes and 16 plastome CDS genes, corresponding to 2.39% of the R. chingii mitogenome. Mitogenomes had shorter but more collinear regions, evidenced by a comparison of the organelles of non-parasitic R. chingii, hemiparasitic Pedicularis chinensis, and holoparasitic Aeginetia indica in the Orobanchaceae family. Moreover, from non-parasitic to holoparasitic species, the genome size in the mitogenomes of Orobanchaceae species did not decrease gradually. Instead, the smallest mitogenome was found in the hemiparasitic species P. chinensis, with a size of 225,612 bp. The findings fill the gap in the mitogenome research of the medicinal plant R. chingii, promote the progress of the organelle genome research of the Orobanchaceae family, and provide clues for molecular breeding.

The telomere-to-telomere (T2T) genome provides insights into the evolution of specialized centromere sequences in sandalwood

BACKGROUND

Sandalwood, a prized hemiparasitic plant, is highly sought in the commercial market because of its aromatic core materia. The structure and stability of the genome are instrumental in the rapid adaptation of parasitic plants to their surroundings. However, there is a conspicuous lack of research on the genomic-level adaptive evolution of sandalwood.

RESULTS

In this study, we assembled a gap-free telomere-to-telomere (T2T) reference genome for Santalum album using PacBio HiFi, Hi-C, and ultra-long ONT data. The T2T reference genome (Sal_t2t) encompassed annotations of 24,171 genes and 25.34% repetitive sequences, in addition to all 10 centromeres and 20 telomeres across the 10 chromosomes. The results revealed that the 3 distinct parasitic species of Santalales had diverse centromeric compositions. The Copia-type long terminal repeat transposon emerged as the most significant in the S. album genome, constituting the primary sequence of the centromere and influencing gene expression. Third, in sandalwood, the presence of Copia affected the size of the centromeres and, consequently, the genome size. Identification of the sandalwood T2T genome in this study also enabled the identification of more precise organelle transfer fragments.

CONCLUSIONS

Our research provides a sandalwood T2T genome, laying the groundwork for future investigations on the evolution of energy organs in parasitic plants. Moreover, it offers novel insights into the function and evolution of centromeres, as well as the mechanisms of adaptation and parasitism.

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.

Complete Mitochondrial Genome Assembly of an Upland Wild Rice Species, Oryza granulata and Comparative Mitochondrial Genomic Analyses of the Genus Oryza

Wild upland rice species, including Oryza granulata, possess unique characteristics that distinguish them from other Oryza species. For instance, O. granulata characteristically has a GG genome and is accordingly classified as a basal lineage of the genus Oryza. Here, we deployed a versatile hybrid approach by integrating Illumina and PacBio sequencing data to generate a high-quality mitochondrial genome (mitogenome) assembly for O. granulata. The mitogenome of O. granulata was 509,311 base pairs (bp) with sixty-seven genes comprising two circular chromosomes, five ribosomal RNA (rRNA) coding genes, twenty-five transfer RNA (tRNA) coding genes, and thirty-seven genes coding for proteins. We identified a total of 378 simple sequence repeats (SSRs). The genome also contained 643 pairs of dispersed repeats comprising 340 palindromic and 303 forward. In the O. granulata mitogenome, the length of 57 homologous fragments in the chloroplast genome occupied 5.96% of the mitogenome length. Collinearity analysis of three Oryza mitogenomes revealed high structural variability and frequent rearrangements. Phylogenetic analysis showed that, compared to other related genera, O. granulata had the closest genetic relationship with mitogenomes reported for all members of Oryza, and occupies a position at the base of the Oryza phylogeny. Comparative analysis of complete mitochondrial genome assemblies for Oryza species revealed high levels of mitogenomic diversity, providing a foundation for future conservation and utilization of wild rice biodiversity.

Assembly and comparative analysis of the complete mitochondrial genome of Viburnum chinshanense

BACKGROUND

Viburnum chinshanense is an endemic species found exclusively in the North-Central and South-Central regions of China. This species is a lush garden ornamental tree and is extensively utilized for vegetation restoration in rocky desertification areas.

RESULTS

In this study, we obtained 13.96 Gb of Oxford Nanopore data for the whole genome, and subsequently, by combining Illumina short-reads, we successfully assembled the complete mitochondrial genome (mitogenome) of the V. chinshanense using a hybrid assembly strategy. The assembled genome can be described as a circular genome. The total length of the V. chinshanense mitogenome measures 643,971 bp, with a GC content of 46.18%. Our annotation efforts have revealed a total of 39 protein-coding genes (PCGs), 28 tRNA genes, and 3 rRNA genes within the V. chinshanense mitogenome. The analysis of repeated elements has identified 212 SSRs, 19 long tandem repeat elements, and 325 pairs of dispersed repeats in the V. chinshanense mitogenome. Additionally, we have investigated mitochondrial plastid DNAs (MTPTs) and identified 21 MTPTs within the mitogenome and plastidial genome. These MTPTs collectively span a length of 9,902 bp, accounting for 1.54% of the mitogenome. Moreover, employing Deepred-mt, we have confidently predicted 623 C to U RNA editing sites across the 39 protein-coding genes. Furthermore, extensive genomic rearrangements have been observed between V. chinshanense and the mitogenomes of related species. Interestingly, we have also identified a bacterial-derived tRNA gene (trnC-GCA) in the V. chinshanense mitogenome. Lastly, we have inferred the phylogenetic relationships of V. chinshanense with other angiosperms based on mitochondrial PCGs.

CONCLUSIONS

This study marks the first report of a mitogenome from the Viburnum genus, offering a valuable genomic resource for exploring the evolution of mitogenomes within the Dipsacales order.

Comparative analysis of organellar genomes between diploid and tetraploid Chrysanthemum indicum with its relatives

Chrysanthemum indicum, a species native to Eastern Asia is well known as one of the progenitor species of the cultivated Chrysanthemum which is grown for its ornamental and medicinal value. Previous genomic studies on Chrysanthemum have largely ignored the dynamics of plastid genome (plastome) and mitochondria genome (mitogenome) evolution when analyzing this plant lineage. In this study, we sequenced and assembled the plastomes and mitogenomes of diploid and tetraploid C. indicum as well as the morphologically divergent variety C. indicum var. aromaticum. We used published data from 27 species with both plastome and mitogenome complete sequences to explore differences in sequence evolution between the organellar genomes. The size and structure of organellar genome between diploid and tetraploid C. indicum were generally similar but the tetraploid C. indicum and C. indicum var. aromaticum were found to contain unique sequences in the mitogenomes which also contained previously undescribed open reading frames (ORFs). Across Chrysanthemum mitogenome structure varied greatly but sequences transferred from plastomes in to the mitogenomes were conserved. Finally, differences observed between mitogenome and plastome gene trees may be the result of the difference in the rate of sequence evolution between genes in these two genomes. In total the findings presented here greatly expand the resources for studying Chrysanthemum organellar genome evolution with possible applications to conservation, breeding, and gene banking in the future.

The pan-plastome of tartary buckwheat (fagopyrum tataricum): key insights into genetic diversity and the history of lineage divergence

BACKGROUND

Tartary buckwheat (Fagopyrum tataricum) is an important food and medicine crop plant, which has been cultivated for 4000 years. A nuclear genome has been generated for this species, while an intraspecific pan-plastome has yet to be produced. As such a detailed understanding of the maternal genealogy of Tartary buckwheat has not been thoroughly investigated.

RESULTS

In this study, we de novo assembled 513 complete plastomes of Fagopyrum and compared with 8 complete plastomes of Fagopyrum downloaded from the NCBI database to construct a pan-plastome for F. tartaricum and resolve genomic variation. The complete plastomes of the 513 newly assembled Fagopyrum plastome sizes ranged from 159,253 bp to 159,576 bp with total GC contents ranged from 37.76 to 37.97%. These plastomes all maintained the typical quadripartite structure, consisting of a pair of inverted repeat regions (IRA and IRB) separated by a large single copy region (LSC) and a small single copy region (SSC). Although the structure and gene content of the Fagopyrum plastomes are conserved, numerous nucleotide variations were detected from which population structure could be resolved. The nucleotide variants were most abundant in the non-coding regions of the genome and of those the intergenic regions had the most. Mutational hotspots were primarily found in the LSC regions. The complete 521 Fagopyrum plastomes were divided into five genetic clusters, among which 509 Tartary buckwheat plastomes were divided into three genetic clusters (Ft-I/Ft-II/Ft-III). The genetic diversity in the Tartary buckwheat genetic clusters was the greatest in Ft-III, and the genetic distance between Ft-I and Ft-II was the largest. Based on the results of population structure and genetic diversity analysis, Ft-III was further subdivided into three subgroups Ft-IIIa, Ft-IIIb, and Ft-IIIc. Divergence time estimation indicated that the genera Fagopyrum and Rheum (rhubarb) shared a common ancestor about 48 million years ago (mya) and that intraspecies divergence in Tartary buckwheat began around 0.42 mya.

CONCLUSIONS

The resolution of pan-plastome diversity in Tartary buckwheat provides an important resource for future projects such as marker-assisted breeding and germplasm preservation.

Variation in mitogenome structural conformation in wild and cultivated lineages of sorghum corresponds with domestication history and plastome evolution

BACKGROUND

Mitochondria are organelles within eukaryotic cells that are central to the metabolic processes of cellular respiration and ATP production. However, the evolution of mitochondrial genomes (mitogenomes) in plants is virtually unknown compared to animal mitogenomes or plant plastids, due to complex structural variation and long stretches of repetitive DNA making accurate genome assembly more challenging. Comparing the structural and sequence differences of organellar genomes within and between sorghum species is an essential step in understanding evolutionary processes such as organellar sequence transfer to the nuclear genome as well as improving agronomic traits in sorghum related to cellular metabolism.

RESULTS

Here, we assembled seven sorghum mitochondrial and plastid genomes and resolved reticulated mitogenome structures with multilinked relationships that could be grouped into three structural conformations that differ in the content of repeats and genes by contig. The grouping of these mitogenome structural types reflects the two domestication events for sorghum in east and west Africa.

CONCLUSIONS

We report seven mitogenomes of sorghum from different cultivars and wild sources. The assembly method used here will be helpful in resolving complex genomic structures in other plant species. Our findings give new insights into the structure of sorghum mitogenomes that provides an important foundation for future research into the improvement of sorghum traits related to cellular respiration, cytonuclear incompatibly, and disease resistance.