Deciphering the Multi-Chromosomal Mitochondrial Genome of Populus simonii

Characterize the Complete Mitogenome of Semiaquilegia guangxiensis and Assess the Efficiency of the Mitochondrial Genes in Ranunculales Phylogeny

Plant mitochondrial genome (mitogenome) has crucial functions underpinning survive, development, and reproduction of organisms. However, the complete mitogenomes have been far less assembled and annotated than plastomes and even nuclear genomes in plants, due to their highly frequent and long repeat sequences and genomic structural variations. These further hinder the understanding of the mitogenome evolution and restrict potential applications in phylogenetic analyses. In this study, we sequenced, assembled, and annotated the complete mitogenome of Semiaquilegia guangxiensis and explored its evolution and usefulness in phylogenetics. The results showed that the mitogenome was composed of two independent molecules, which had a total length of 522,981 bp, a GC content of 45.69%, and 58 genes including 34 protein-coding genes (PCGs), 21 tRNA genes, and three rRNA genes. A generalized codon usage preference was observed among the mitochondrial PCGs, and a total of 665 potential RNA editing sites were identified across the 34 mitochondrial PCGs, all of which were base C-to-U edits. Moreover, a large number of repetitive mitogenome sequences and chloroplast-sourced sequences transferred to the mitogenome were detected. The largest collinear block identified between S. guangxiensis and Paropyrum anemonoides was 4282 bp in length. The phylogenetic analyses based on the mitochondrial gene sequences resolved the phylogenetic relationships within Ranunculales, in which Semiaquilegia was close to Paropyrum. The 19 PCGs were ranked according to their efficiencies on phylogenetic resolution based on several metrics, and the combined metric suggested matR, rps3, and nad5 were the top three loci contributing most to phylogeny. As the first reported mitogenome in Semiaquilegia, our findings enrich the limited mitogenome library of plants, reecho the complex evolutionary dynamics of the mitogenome, and highlight the usefulness of mitochondrial gene sequences in phylogenetics.

Assembly and comparative analysis of the complete mitochondrial of Spodiopogon sagittifolius, an endemic and protective species from Yunnan, China

BACKGROUND

Spodiopogon sagittifolius, a C4 plant closely related to cultivated crops, is an edible resource and a Class II nationally protected species in China. Endemic to Yunnan, its populations are declining due to habitat destruction, highlighting its resource and conservation importance. Despite its significance, the molecular phylogenetic relationships and genetic mechanisms of adaptive evolution in the genus Spodiopogon remain poorly understood.

RESULTS

We successfully assembled and annotated the first mitochondrial genome of S. sagittifolius using HiFi sequencing and the PMAT tool. The genome is 500,699 bp in length with a GC content of 43.15%. Synteny and dN/dS analyses revealed structural and functional conservation of mitochondrial genomes in closely related species, with most protein-coding genes under purifying selection (dN/dS < 1). Notably, nad2 exhibited signs of positive selection (dN/dS = 1.49), indicating potential adaptive evolution. Extensive RNA editing events were detected across 27 protein-coding genes, predominantly involving C-to-U conversions, with synonymous mutations accounting for 49.65% of the edits. Strong codon usage bias favoring A/U-ending codons and the identification of repeat sequences suggest enhanced mitochondrial efficiency and stress adaptation. Phylogenetic analyses confirmed the taxonomic placement of S. sagittifolius within the Andropogoneae tribe.

CONCLUSIONS

This study provides the first insights into the mitochondrial genome evolution of S. sagittifolius, highlighting key features linked to stress tolerance and adaptive evolution. These findings establish a foundation for its conservation and potential domestication, with implications for crop improvement and ecological resilience.

De novo assembly and comparative analysis of cherry (Prunus subgenus Cerasus) mitogenomes

Prunus subgenus Cerasus (Mill) A. Gray, commonly known as cherries and cherry blossoms, possesses significant edible and ornamental value. However, the mitochondrial genomes (mitogenomes) of cherry species remain largely unexplored. Here, we successfully assembled the mitogenomes of five cherry species (P. campanulata, P. fruticosa, P. mahaleb, P. pseudocerasus, and P. speciosa), revealing common circular structures. The assembled mitogenomes exhibited sizes ranging from 383,398 bp to 447,498 bp, with GC content varying between 45.54% and 45.76%. A total of 62 to 69 genes were annotated, revealing variability in the copy number of protein-coding genes (PCGs) and tRNA genes. Mitogenome collinearity analysis indicated genomic rearrangements across Prunus species, driven by repetitive sequences, particularly dispersed repeats. Additionally, the five cherry species displayed highly conserved codon usage and RNA editing patterns, highlighting the evolutionary conservation of the mitochondrial PCGs. Phylogenetic analyses confirmed the monophyly of subg. Cerasus, although notable phylogenetic incongruences were observed between the mitochondrial and plastid datasets. These results provide significant genomic resources for forthcoming studies on the evolution and molecular breeding of cherry mitogenomes, enhancing the overall comprehension of mitogenome structure and evolution within Prunus.

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.

Comprehensive analysis of the mitochondrial genome of Iris domestica emphasizing multichromosomal organization and repeat-mediated homologous recombination

Background

Iris domestica is a perennial herb valued for both its ornamental and medicinal properties. Despite its significance, no comprehensive analysis of its mitochondrial genome has been previously reported. Plant mitochondrial genomes are known for their large size, structural complexity, and frequent recombination events. This study aims to provide the first complete assembly and characterization of the mitochondrial genome of I. domestica, with a focus on its structure, gene content, repeat elements, and RNA editing sites.

Results

We used GetOrganelle and Unicycler software to hybrid assemble Nanopore and Illumina data to obtain the mitochondrial genome of I. domestica. The mitochondrial genome of I. domestica consists of four contigs: contig1 (222,498 bp), contig2 (90,780 bp), contig3 (42,563 bp), and contig4 (39,247 bp). Two repeat sequences, R1 (7,784 bp) and R2 (3,519 bp), facilitate the conformation of three circular chromosomes, suggesting a complex multi-chromosomal structure. A total of 34 protein-coding genes, including 24 core genes and 10 non-core genes were identified. Analysis of tandem repeat elements revealed significant variability, with Chromosome 1 showing the highest diversity of SSRs and scattered repeats. Additionally, 20 homologous fragments were identified between the mitochondrial and chloroplast genomes, accounting for 1.10% of the mitochondrial genome. Phylogenetic analysis based on 24 conserved mitochondrial genes placed I. domestica in close relation to Iris domestica and Crocus sativus. Furthermore, 545 RNA editing sites were identified, with notable variations across genes, suggesting that RNA editing plays a significant role in regulating mitochondrial gene expression.

Conclusion

The complete assembly of the I. domestica mitochondrial genome reveals a complex multichromosomal structure characterized by recombination events. The high number of RNA editing sites and the presence of transferred plastid DNA highlight the dynamic nature of the genome, contributing to its adaptability and evolution. These findings provide a genetic foundation into the plant's medicinal properties, adaptive mechanisms, and potential for environmental resilience.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Comparative organelle genomics in Daphniphyllaceae reveal phylogenetic position and organelle structure evolution

The family Daphniphyllaceae has a single genus, and no relevant comparative phylogenetic study has been reported on it. To explore the phylogenetic relationships and organelle evolution mechanisms of Daphniphyllaceae species, we sequenced and assembled the chloroplast and mitochondrial genomes of Daphniphyllum macropodum. We also conducted comparative analyses of organelles in Daphniphyllaceae species in terms of genome structure, phylogenetic relationships, divergence times, RNA editing events, and evolutionary rates, etc. Results indicated differences in the evolutionary patterns of the plastome and mitogenome in D. macropodum. The plastome had a more conserved structure but a faster nucleotide substitution rate, and the mitogenome showed a more complex structure while the mitotic genome shows a more complex structure but a slower nucleotide substitution rate. We identified several unidirectional protein-coding gene transfer events from the plastome to the mitogenome based on homology analysis, but no transfer events occurred from the mitogenome to the plastome. Multiple TE fragments existed in organelle genomes, and two organelles showed different preferences for nuclear TE insertion types. The estimation of divergence time indicated that the differentiation of Daphniphyllaceae and Altingiaceae at around 29.86 Mya might be due to the dramatic uplift of Tibetan Plateau during the Oligocene. About 75% of codon changes in organelles were found to be hydrophilic to hydrophobic amino acids. The RNA editing in protein-coding transcripts is the result of amino acid changes to increase their hydrophobicity and conservation in alleles, which may contribute to the formation of functional 3D structures in proteins. This study would enrich genomic resources and provide valuable insights into the structural dynamics and molecular biology of Daphniphyllaceae species.

Comprehensive analysis of the multi-rings mitochondrial genome of Populus tomentosa

BACKGROUND

Populus tomentosa, known as Chinese white poplar, is indigenous and distributed across large areas of China, where it plays multiple important roles in forestry, agriculture, conservation, and urban horticulture. However, limited accessibility to the mitochondrial (mt) genome of P. tomentosa impedes phylogenetic and population genetic analyses and restricts functional gene research in Salicaceae family.

RESULTS

Single-molecule real-time (SMRT) sequencing technology was used to sequence, assemble, and annotate the mt genome of P. tomentosa. This genome has a complex structure composed of four circular molecules ranging from 153,004 to 330,873 base pairs (bp). Each of these four circular molecules contains unique gene sequences that constitute the mt genome of P. tomentosa. The mt genome comprises 69 functional genes, including 38 protein-coding genes (PCGs), 26 tRNA genes, and 5 rRNA genes. After removing duplications, 19 different tRNA coding genes remain, though only 10 amino acids can be recognized. The noncoding region constitutes 93.38% of the mt genome, comprising a large number of repetitive sequences, gene spacer regions, and insertion from chloroplast sequences. Specifically, 40 chloroplast-derived sequences, with a total length of 24,381 bp, were identified in P. tomentosa.

CONCLUSIONS

In the current study, the results provide mitochondrial genomic evidence for the maternal origin of P. tomentosa and enhance understanding of the gene dialog between organelle genomes, contributing to the conservation and utilization of the genetic resources of P. tomentosa.

Analysis of the complete mitochondrial genome of Panax quinquefolius reveals shifts from cis-splicing to trans-splicing of intron cox2i373

Panax quinquefolius is a perennial plant with medicinal values. In this study, we assembled the complete mitochondrial genome (mitogenome) of P. quinquefolius using PMAT assembler. The total length of P. quinquefolius mitogenome is 573,154 bp. We annotated a total of 34 protein-coding genes (PCGs), 35 tRNA genes, and 6 rRNA genes in this mitogenome. The analysis of repetitive elements shows that there are 153 SSRs, 24 tandem repeats and 242 pairs of dispersed repeats this mitogenome. Also, we found 24 homologous sequences with a total length of 64,070 bp among its mitogenome and plastome, accounting for 41.05 % of the plastome, and 11.18 % of the mitogenome, showing a remarkable frequent sequence dialogue between plastome and mitogenomes. Besides, a total of 583 C to U RNA editing sites on 34 PCGs of high confidence were predicted by using Deepred-mt. We also inferred the phylogenetic relationships of P. quinquefolius and other angiosperms based on mitochondrial PCGs. Finally, we observed a shift from cis- to trans-splicing in P. quinquefolius for two mitochondrial introns, namely cox2i373 and nad1i728, and a pair of 48 bp short repetitive sequences may be associated with the breaking and rearrangement of the cox2i373 intron. The fragmentation of the cox2i373 intron was further confirmed by our PCR amplification experiments. In summary, our report on the P. quinquefolius mitogenome provides a new perspective on the intron evolution of the mitogenome.

Comprehensive analysis of the first complete mitogenome and plastome of a traditional Chinese medicine Viola diffusa

BACKGROUND

Viola diffusa is used in the formulation of various Traditional Chinese Medicines (TCMs), including antiviral, antimicrobial, antitussive, and anti-inflammatory drugs, due to its richness in flavonoids and triterpenoids. The biosynthesis of these compounds is largely mediated by cytochrome P450 enzymes, which are primarily located in the membranes of mitochondria and the endoplasmic reticulum.

RESULTS

This study presents the complete assembly of the mitogenome and plastome of Viola diffusa. The circular mitogenome spans 474,721 bp with a GC content of 44.17% and encodes 36 unique protein-coding genes, 21 tRNA, and 3 rRNA. Except for the RSCU values of 1 observed for the start codon (AUG) and tryptophan (UGG), the mitochondrial protein-coding genes exhibited a codon usage bias, with most estimates deviating from 1, similar to patterns observed in closely related species. Analysis of repetitive sequences in the mitogenome demonstrated potential homologous recombination mediated by these repeats. Sequence transfer analysis revealed 24 homologous sequences shared between the mitogenome and plastome, including nine full-length genes. Collinearity was observed among Viola diffusa species within the other members of Malpighiales order, indicated by the presence of homologous fragments. The length and arrangement of collinear blocks varied, and the mitogenome exhibited a high frequency of gene rearrangement.

CONCLUSIONS

We present the first complete assembly of the mitogenome and plastome of Viola diffusa, highlighting its implications for pharmacological, evolutionary, and taxonomic studies. Our research underscores the multifaceted importance of comprehensive mitogenome analysis.

The mitochondrial genome of Lavandula angustifolia Mill. (Lamiaceae) sheds light on its genome structure and gene transfer between organelles

BACKGROUND

Lavandula angustifolia holds importance as an aromatic plant with extensive applications spanning the fragrance, perfume, cosmetics, aromatherapy, and spa sectors. Beyond its aesthetic and sensory applications, this plant offers medicinal benefits as a natural herbal remedy and finds use in household cleaning products. While extensive genomic data, inclusive of plastid and nuclear genomes, are available for this species, researchers have yet to characterize its mitochondrial genome. This gap in knowledge hampers deeper understanding of the genome organization and its evolutionary significance.

RESULTS

Through the course of this study, we successfully assembled and annotated the mitochondrial genome of L. angustifolia, marking a first in this domain. This assembled genome encompasses 61 genes, which comprise 34 protein-coding genes, 24 transfer RNA genes, and three ribosomal RNA genes. We identified a chloroplast sequence insertion into the mitogenome, which spans a length of 10,645 bp, accounting for 2.94% of the mitogenome size. Within these inserted sequences, there are seven intact tRNA genes (trnH-GUG, trnW-CCA, trnD-GUC, trnS-GGA, trnN-GUU, trnT-GGU, trnP-UGG) and four complete protein-coding genes (psbA, rps15, petL, petG) of chloroplast derivation. Additional discoveries include 88 microsatellites, 15 tandem repeats, 74 palindromic repeats, and 87 forward long repeats. An RNA editing analysis highlighted an elevated count of editing sites in the cytochrome c oxidase genes, notably ccmB with 34 editing sites, ccmFN with 32, and ccmC with 29. All protein-coding genes showed evidence of cytidine-to-uracil conversion. A phylogenetic analysis, utilizing common protein-coding genes from 23 Lamiales species, yielded a tree with consistent topology, supported by high confidence values.

CONCLUSIONS

Analysis of the current mitogenome resource revealed its typical circular genome structure. Notably, sequences originally from the chloroplast genome were found within the mitogenome, pointing to the occurrence of horizontal gene transfer between organelles. This assembled mitogenome stands as a valuable resource for subsequent studies on mitogenome structures, their evolution, and molecular biology.

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.

Complete mitochondrial genome assembly of Juglans regia unveiled its molecular characteristics, genome evolution, and phylogenetic implications

BACKGROUND

The Persian walnut (Juglans regia), an economically vital species within the Juglandaceae family, has seen its mitochondrial genome sequenced and assembled in the current study using advanced Illumina and Nanopore sequencing technology.

RESULTS

The 1,007,576 bp mitogenome of J. regia consisted of three circular chromosomes with a 44.52% GC content encoding 39 PCGs, 47 tRNA, and five rRNA genes. Extensive repetitive sequences, including 320 SSRs, 512 interspersed, and 83 tandem repeats, were identified, contributing to genomic complexity. The protein-coding sequences (PCGs) favored A/T-ending codons, and the codon usage bias was primarily shaped by selective pressure. Intracellular gene transfer occurred among the mitogenome, chloroplast, and nuclear genomes. Comparative genomic analysis unveiled abundant structure and sequence variation among J. regia and related species. The results of selective pressure analysis indicated that most PCGs underwent purifying selection, whereas the atp4 and ccmB genes had experienced positive selection between many species pairs. In addition, the phylogenetic examination, grounded in mitochondrial genome data, precisely delineated the evolutionary and taxonomic relationships of J. regia and its relatives. We identified a total of 539 RNA editing sites, among which 288 were corroborated by transcriptome sequencing data. Furthermore, expression profiling under temperature stress highlighted the complex regulation pattern of 28 differently expressed PCGs, wherein NADH dehydrogenase and ATP synthase genes might be critical in the mitochondria response to cold stress.

CONCLUSIONS

Our results provided valuable molecular resources for understanding the genetic characteristics of J. regia and offered novel perspectives for population genetics and evolutionary studies in Juglans and related woody species.

Cytoplasmic genomes of Jasminum sambac reveal divergent sub-mitogenomic conformations and a large nuclear chloroplast-derived insertion

BACKGROUND

Jasminum sambac, a widely recognized ornamental plant prized for its aromatic blossoms, exhibits three flora phenotypes: single-petal ("SP"), double-petal ("DP"), and multi-petal ("MP"). The lack of detailed characterization and comparison of J. sambac mitochondrial genomes (mitogenomes) hinders the exploration of the genetic and structural diversity underlying the varying floral phenotypes in jasmine accessions.

RESULTS

Here, we de novo assembled three mitogenomes of typical phenotypes of J. sambac, "SP", "DP", and "MP-hutou" ("HT"), with PacBio reads and the "HT" chloroplast (cp) genome with Illumina reads, and verified them with read mapping and fluorescence in situ hybridization (FISH). The three mitogenomes present divergent sub-genomic conformations, with two, two, and four autonomous circular chromosomes ranging in size from 35.7 kb to 405.3 kb. Each mitogenome contained 58 unique genes. Ribosome binding sites with conserved AAGAAx/AxAAAG motifs were detected upstream of uncanonical start codons TTG, CTG and GTG. The three mitogenomes were similar in genomic content but divergent in structure. The structural variations were mainly attributed to recombination mediated by a large (~ 5 kb) forward repeat pair and several short repeats. The three jasmine cp. genomes showed a well-conserved structure, apart from a 19.9 kb inversion in "HT". We identified a 14.3 kb "HT"-specific insertion on Chr7 of the "HT" nuclear genome, consisting of two 7 kb chloroplast-derived fragments with two intact ndhH and rps15 genes, further validated by polymerase chain reaction (PCR). The well-resolved phylogeny suggests faster mitogenome evolution in J. sambac compared to other Oleaceae species and outlines the mitogenome evolutionary trajectories within Lamiales. All evidence supports that "DP" and "HT" evolved from "SP", with "HT" being the most recent derivative of "DP".

CONCLUSION

The comprehensive characterization of jasmine organelle genomes has added to our knowledge of the structural diversity and evolutionary trajectories behind varying jasmine traits, paving the way for in-depth exploration of mechanisms and targeted genetic research.

De novo assembly and comprehensive analysis of the mitochondrial genome of Taxus wallichiana reveals different repeats mediate recombination to generate multiple conformations

Taxus plants are the exclusive source of paclitaxel, an anticancer drug with significant medicinal and economic value. Interspecies hybridization and gene introgression during evolution have obscured distinctions among Taxus species, complicating their phylogenetic classification. While the chloroplast genome of Taxus wallichiana, a widely distributed species in China, has been sequenced, its mitochondrial genome (mitogenome) remains uncharacterized.We sequenced and assembled the T. wallichiana mitogenome using BGI short reads and Nanopore long reads, facilitating comparisons with other gymnosperm mitogenomes. The T. wallichiana mitogenome spanning 469,949 bp, predominantly forms a circular configuration with a GC content of 50.51%, supplemented by 3 minor configurations mediated by one pair of LRs and two pairs of IntRs. It includes 32 protein-coding genes, 7 tRNA genes, and 3 rRNA genes, several of which exist in multiple copies.We detailed the mitogenome's structure, codon usage, RNA editing, and sequence migration between organelles, constructing a phylogenetic tree to elucidate evolutionary relationships. Unlike typical gymnosperm mitochondria, T. wallichiana shows no evidence of mitochondrial-plastid DNA transfer (MTPT), highlighting its unique genomic architecture. Synteny analysis indicated extensive genomic rearrangements in T. wallichiana, likely driven by recombination among abundant repetitive sequences. This study offers a high-quality T. wallichiana mitogenome, enhancing our understanding of gymnosperm mitochondrial evolution and supporting further cultivation and utilization of Taxus species.

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

Introduction

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

Methods

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

Results

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

Discussion

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

Comprehensive analysis of the Lycopodium japonicum mitogenome reveals abundant tRNA genes and cis-spliced introns in Lycopodiaceae species

Lycophytes and ferns represent one of the earliest-diverging lineages of vascular plants, with the Lycopodiaceae family constituting the basal clade among lycophytes. In this research, we successfully assembled and annotated the complete Lycopodium japonicum Thunb. (L. japonicum) mitochondrial genome (mitogenome) utilizing PacBio HiFi sequencing data, resulting in a single circular molecule with a size of 454,458 bp. 64 unique genes were annotated altogether, including 34 protein-coding genes, 27 tRNAs and 3 rRNAs. It also contains 32 group II introns, all of which undergo cis-splicing. We identified 195 simple sequence repeats, 1,948 dispersed repeats, and 92 tandem repeats in the L. japonicum mitogenome. Collinear analysis indicated that the mitogenomes of Lycopodiaceae are remarkably conserved compared to those of other vascular plants. We totally identified 326 RNA editing sites in 31 unique protein-coding genes with 299 sites converting cytosine to uracil and 27 sites the reverse. Notably, the L. japonicum mitogenome has small amounts foreign DNA from plastid or nuclear origin, accounting for only 2.81% of the mitogenome. The maximum likelihood phylogenetic analysis based on 23 diverse land plant mitogenomes and plastid genomes supports the basal position of lycophytes within vascular plants and they form a sister clade to all other vascular lineages, which is consistent with the PPG I classification system. As the first reported mitogenome of Lycopodioideae subfamily, this study enriches our understanding of Lycopodium mitogenomes, and sets the stage for future research on mitochondrial diversity and evolution within the lycophytes and ferns.

Assembly and comparative analysis of the complete mitochondrial genome of Fritillaria ussuriensis Maxim. (Liliales: Liliaceae), an endangered medicinal plant

BACKGROUND

Fritillaria ussuriensis is an endangered medicinal plant known for its notable therapeutic properties. Unfortunately, its population has drastically declined due to the destruction of forest habitats. Thus, effectively protecting F. ussuriensis from extinction poses a significant challenge. A profound understanding of its genetic foundation is crucial. To date, research on the complete mitochondrial genome of F. ussuriensis has not yet been reported.

RESULTS

The complete mitochondrial genome of F. ussuriensis was sequenced and assembled by integrating PacBio and Illumina sequencing technologies, revealing 13 circular chromosomes totaling 737,569 bp with an average GC content of 45.41%. A total of 55 genes were annotated in this mitogenome, including 2 rRNA genes, 12 tRNA genes, and 41 PCGs. The mitochondrial genome of F. ussuriensis contained 192 SSRs and 4,027 dispersed repeats. In the PCGs of F. ussuriensis mitogenome, 90.00% of the RSCU values exceeding 1 exhibited a preference for A-ended or U-ended codons. In addition, 505 RNA editing sites were predicted across these PCGs. Selective pressure analysis suggested negative selection on most PCGs to preserve mitochondrial functionality, as the notable exception of the gene nad3 showed positive selection. Comparison between the mitochondrial and chloroplast genomes of F. ussuriensis revealed 20 homologous fragments totaling 8,954 bp. Nucleotide diversity analysis revealed the variation among genes, and gene atp9 was the most notable. Despite the conservation of GC content, mitogenome sizes varied significantly among six closely related species, and colinear analysis confirmed the lack of conservation in their genomic structures. Phylogenetic analysis indicated a close relationship between F. ussuriensis and Lilium tsingtauense.

CONCLUSIONS

In this study, we sequenced and annotated the mitogenome of F. ussuriensis and compared it with the mitogenomes of other closely related species. In addition to genomic features and evolutionary position, this study also provides valuable genomic resources to further understand and utilize this medicinal plant.

Complete mitochondrial genome of Hippophae tibetana: insights into adaptation to high-altitude environments

Hippophae tibetana, belonging to the Elaeagnaceae family, is an endemic plant species of the Qinghai-Tibet Plateau, valued for its remarkable ecological restoration capabilities, as well as medicinal and edible properties. Despite being acknowledged as a useful species, its mitochondrial genome data and those of other species of the Elaeagnaceae family are lacking to date. In this study, we, for the first time, successfully assembled the mitochondrial genome of H. tibetana, which is 464,208 bp long and comprises 31 tRNA genes, 3 rRNA genes, 37 protein-coding genes, and 3 pseudogenes. Analysis of the genome revealed a high copy number of the trnM-CAT gene and a high prevalence of repetitive sequences, both of which likely contribute to genome rearrangement and adaptive evolution. Through nucleotide diversity and codon usage bias analyses, we identified specific genes that are crucial for adaptation to high-altitude conditions. Notably, genes such as atp6, ccmB, nad4L, and nad7 exhibited signs of positive selection, indicating the presence of unique adaptive traits for survival in extreme environments. Phylogenetic analysis confirmed the close relationship between the Elaeagnaceae family and other related families, whereas intergenomic sequence transfer analysis revealed a substantial presence of homologous fragments among the mitochondrial, chloroplast, and whole genomes, which may be linked to the high-altitude adaptation mechanisms of H. tibetana. The findings of this study not only enrich our knowledge of H. tibetana molecular biology but also advance our understanding of the adaptive evolution of plants on the Qinghai-Tibet Plateau. This study provides a solid scientific foundation for the molecular breeding, conservation, and utilization of H. tibetana genetic resources.

A Systematic Review and Developmental Perspective on Origin of CMS Genes in Crops

Cytoplasmic male sterility (CMS) arises from the incompatibility between the nucleus and cytoplasm as typical representatives of the chimeric structures in the mitochondrial genome (mitogenome), which has been extensively applied for hybrid seed production in various crops. The frequent occurrence of chimeric mitochondrial genes leading to CMS is consistent with the mitochondrial DNA (mtDNA) evolution. The sequence conservation resulting from faithfully maternal inheritance and the chimeric structure caused by frequent sequence recombination have been defined as two major features of the mitogenome. However, when and how these chimeric mitochondrial genes appear in the context of the highly conserved reproduction of mitochondria is an enigma. This review, therefore, presents the critical view of the research on CMS in plants to elucidate the mechanisms of this phenomenon. Generally, distant hybridization is the main mechanism to generate an original CMS source in natural populations and in breeding. Mitochondria and mitogenomes show pleomorphic and dynamic changes at key stages of the life cycle. The promitochondria in dry seeds develop into fully functioning mitochondria during seed imbibition, followed by massive mitochondria or mitogenome fusion and fission in the germination stage along with changes in the mtDNA structure and quantity. The mitogenome stability is controlled by nuclear loci, such as the nuclear gene Msh1. Its suppression leads to the rearrangement of mtDNA and the production of heritable CMS genes. An abundant recombination of mtDNA is also often found in distant hybrids and somatic/cybrid hybrids. Since mtDNA recombination is ubiquitous in distant hybridization, we put forward a hypothesis that the original CMS genes originated from mtDNA recombination during the germination of the hybrid seeds produced from distant hybridizations to solve the nucleo-cytoplasmic incompatibility resulting from the allogenic nuclear genome during seed germination.

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

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

Assembly and evolutionary analysis of the complete mitochondrial genome of Trichosanthes kirilowii, a traditional Chinese medicinal plant

Trichosanthes kirilowii (T. kirilowii) is a valuable plant used for both medicinal and edible purposes. It belongs to the Cucurbitaceae family. However, its phylogenetic position and relatives have been difficult to accurately determine due to the lack of mitochondrial genomic information. This limitation has been an obstacle to the potential applications of T. kirilowii in various fields. To address this issue, Illumina and Nanopore HiFi sequencing were used to assemble the mitogenome of T. kirilowii into two circular molecules with sizes of 245,700 bp and 107,049 bp, forming a unique multi-branched structure. The mitogenome contains 61 genes, including 38 protein-coding genes (PCGs), 20 tRNAs, and three rRNAs. Within the 38 PCGs of the T. kirilowii mitochondrial genome, 518 potential RNA editing sites were identified. The study also revealed the presence of 15 homologous fragments that span both the chloroplast and mitochondrial genomes. The phylogenetic analysis strongly supports that T. kirilowii belongs to the Cucurbitaceae family and is closely related to Luffa. Collinearity analysis of five Cucurbitaceae mitogenomes shows a high degree of structural variability. Interestingly, four genes, namely atp1, ccmFC, ccmFN, and matR, played significant roles in the evolution of T. kirilowii through selection pressure analysis. The comparative analysis of the T. kirilowii mitogenome not only sheds light on its functional and structural features but also provides essential information for genetic studies of the genus of Cucurbitaceae.

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.

Unraveling the complex evolutionary features of the Cinnamomum camphora mitochondrial genome

KEY MESSAGE

We reported the mitochondrial genome of Cinnamomum camphora for the first time, revealing frequent rearrangement events in the non-coding regions of Magnoliids mitochondrial genomes. As one of the representative species in the Lauraceae family of Magnoliids, Cinnamomum camphora holds significant economic and ecological value. In this study, the mitochondrial genome (mitogenome) of C. camphora was complete assembled and annotated using PacBio HiFi sequencing. The C. camphora mitogenome is characterized by a branch structure, spans 900,894 bp, and contains 43 protein-coding genes (PCGs), 24 tRNAs, and 3 rRNAs. Most of these PCGs are under purifying selection, with only two (ccmFc and rps7) exhibiting signs of positive selection. The C. camphora mitogenome contains numerous repetitive sequences and intracellular gene transfers, with a total of 36 mitochondrial plastid DNAs, amounting to a combined length of 23,816 bp. Comparative analysis revealed that the non-coding regions of Magnoliids mitogenomes have undergone frequent rearrangements during evolution, but the coding sequences remain highly conserved (more than 98% similarity for protein-coding sequences). Furthermore, a maximum-likelihood phylogenetic tree was reconstructed based on 25 PCGs from 23 plant mitogenomes. The analysis supports the closest relationship between C. camphora and C. chekiangense, consistent with the APG IV classification system. This study elucidates the unique evolutionary features of the C. camphora mitogenome, which will provide valuable insights into the study of genetics and evolution of the family Lauraceae.

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

MAIN CONCLUSION

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

Comprehensive analysis of the complete mitochondrial genome of Lilium tsingtauense reveals a novel multichromosome structure

KEY MESSAGE

Lilium tsingtauense mitogenome comprises 27 independent chromosome molecules, it undergoes frequent genomic recombination, and the rate of recombination and mutation between different repetitive sequences affects the formation of multichromosomal structures. Given the extremely large genome of Lily, which likely harbors additional genetic resources, it serves as an ideal material for studying the phylogenetic evolution of organisms. Although the Lilium chloroplast genome has been documented, the sequence of its mitochondrial genome (mitogenome) remains uncharted. Using BGI short reads and Nanopore long reads, we sequenced, assembled, and annotated the mitogenome of Lilium tsingtauense. This effort culminated in the characterization of Lilium's first complete mitogenome. Comparative analysis with other angiosperms revealed the unique multichromosomal structure of the L. tsingtauense mitogenome, spanning 1,125,108 bp and comprising 27 independent circular chromosomes. It contains 36 protein-coding genes, 12 tRNA genes, and 3 rRNA genes, with a GC content of 44.90%. Notably, three chromosomes in the L. tsingtauense mitogenome lack identifiable genes, hinting at the potential existence of novel genes and noncoding elements. The high degree of observed genome fragmentation implies frequent reorganization, with recombination and mutation rates among diverse repetitive sequences likely driving the formation of multichromosomal structures. Our comprehensive analysis, covering genome size, coding genes, structure, RNA editing, repetitive sequences, and sequence migration, sheds light on the evolutionary and molecular biology of multichromosomal mitochondria in Lilium. This high-quality mitogenome of L. tsingtauense not only enriches our understanding of multichromosomal mitogenomes but also establishes a solid foundation for future genome breeding and germplasm innovation in Lilium.

Assembly and comparative analysis of the complete multichromosomal mitochondrial genome of Cymbidium ensifolium, an orchid of high economic and ornamental value

BACKGROUND

Orchidaceae is one of the largest groups of angiosperms, and most species have high economic value and scientific research value due to their ornamental and medicinal properties. In China, Chinese Cymbidium is a popular ornamental orchid with high economic value and a long history. However, to date, no detailed information on the mitochondrial genome of any species of Chinese Cymbidium has been published.

RESULTS

Here, we present the complete assembly and annotation of the mitochondrial genome of Cymbidium ensifolium (L.) Sw. The mitogenome of C. ensifolium was 560,647 bp in length and consisted of 19 circular subgenomes ranging in size from 21,995 bp to 48,212 bp. The genome encoded 35 protein-coding genes, 36 tRNAs, 3 rRNAs, and 3405 ORFs. Repeat sequence analysis and prediction of RNA editing sites revealed a total of 915 dispersed repeats, 162 simple repeats, 45 tandem repeats, and 530 RNA editing sites. Analysis of codon usage showed a preference for codons ending in A/T. Interorganellar DNA transfer was identified in 13 of the 19 chromosomes, with plastid-derived DNA fragments representing 6.81% of the C. ensifolium mitochondrial genome. The homologous fragments of the mitochondrial genome and nuclear genome were also analysed. Comparative analysis showed that the GC content was conserved, but the size, structure, and gene content of the mitogenomes varied greatly among plants with multichromosomal mitogenome structure. Phylogenetic analysis based on the mitogenomes reflected the evolutionary and taxonomic statuses of C. ensifolium. Interestingly, compared with the mitogenomes of Cymbidium lancifolium Hook. and Cymbidium macrorhizon Lindl., the mitogenome of C. ensifolium lost 8 ribosomal protein-coding genes.

CONCLUSION

In this study, we assembled and annotated the mitogenome of C. ensifolium and compared it with the mitogenomes of other Liliidae and plants with multichromosomal mitogenome structures. Our findings enrich the mitochondrial genome database of orchid plants and reveal the rapid structural evolution of Cymbidium mitochondrial genomes, highlighting the potential for mitochondrial genes to help decipher plant evolutionary history.

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.

PMAT: an efficient plant mitogenome assembly toolkit using low-coverage HiFi sequencing data

Complete mitochondrial genomes (mitogenomes) of plants are valuable resources for nucleocytoplasmic interactions, plant evolution, and plant cytoplasmic male sterile line breeding. However, the complete assembly of plant mitogenomes is challenging due to frequent recombination events and horizontal gene transfers. Previous studies have adopted Illumina, PacBio, and Nanopore sequencing data to assemble plant mitogenomes, but the poor assembly completeness, low sequencing accuracy, and high cost limit the sampling capacity. Here, we present an efficient assembly toolkit (PMAT) for de novo assembly of plant mitogenomes using low-coverage HiFi sequencing data. PMAT has been applied to the de novo assembly of 13 broadly representative plant mitogenomes, outperforming existing organelle genome assemblers in terms of assembly accuracy and completeness. By evaluating the assembly of plant mitogenomes from different sequencing data, it was confirmed that PMAT only requires 1× HiFi sequencing data to obtain a complete plant mitogenome. The source code for PMAT is available at https://github.com/bichangwei/PMAT. The developed PMAT toolkit will indeed accelerate the understanding of evolutionary variation and breeding application of plant mitogenomes.

The first mitochondrial genome of Calophyllum soulattri Burm.f

Calophyllum soulattri Burm.f. is traditionally used to treat skin infections and reduce rheumatic pain, yet genetic and genomic studies are still limited. Here, we present the first complete mitochondrial genome of C. soulattri. It is 378,262 bp long with 43.97% GC content, containing 55 genes (30 protein-coding, 5 rRNA, and 20 tRNA). Repeat analysis of the mitochondrial genome revealed 194 SSRs, mostly mononucleotides, and 266 pairs of dispersed repeats ( ≥ 30 bp) that were predominantly palindromic. There were 23 homologous fragments found between the mitochondrial and plastome genomes. We also predicted 345 C-to-U RNA editing sites from 30 protein-coding genes (PCGs) of the C. soulatrii mitochondrial genome. These RNA editing events created the start codon of nad1 and the stop codon of ccmFc. Most PCGs of the C. soulattri mitochondrial genome underwent negative selection, but atp4 and ccmB experienced positive selection. Phylogenetic analyses showed C. soulattri is a sister taxon of Garcinia mangostana. This study has shed light on C. soulattri's evolution and Malpighiales' phylogeny. As the first complete mitochondrial genome in Calophyllaceae, it can be used as a reference genome for other medicinal plant species within the family for future genetic studies.

Assembly and characterization of the complete mitochondrial genome of Ventilago leiocarpa

KEY MESSAGE

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

The Complete Mitochondrial Genome of Paeonia lactiflora Pall. (Saxifragales: Paeoniaceae): Evidence of Gene Transfer from Chloroplast to Mitochondrial Genome

Paeonia lactiflora (P. lactiflora), a perennial plant renowned for its medicinal roots, provides a unique case for studying the phylogenetic relationships of species based on organelle genomes, as well as the transference of DNA across organelle genomes. In order to investigate this matter, we sequenced and characterized the mitochondrial genome (mitogenome) of P. lactiflora. Similar to the chloroplast genome (cpgenome), the mitogenome of P. lactiflora extends across 181,688 base pairs (bp). Its unique quadripartite structure results from a pair of extensive inverted repeats, each measuring 25,680 bp in length. The annotated mitogenome includes 27 protein-coding genes, 37 tRNAs, 8 rRNAs, and two pseudogenes (rpl5, rpl16). Phylogenetic analysis was performed to identify phylogenetic trees consistent with Paeonia species phylogeny in the APG Ⅳ system. Moreover, a total of 12 MTPT events were identified and 32 RNA editing sites were detected during mitogenome analysis of P. lactiflora. Our research successfully compiled and annotated the mitogenome of P. lactiflora. The study provides valuable insights regarding the taxonomic classification and molecular evolution within the Paeoniaceae family.

Insights into the multi-chromosomal mitochondrial genome structure of the xero-halophytic plant Haloxylon Ammodendron (C.A.Mey.) Bunge ex Fenzl

BACKGROUND

Haloxylon ammodendron holds significance as an ecological plant, showcasing remarkable adaptability to desert conditions, halophytic environments, and sand fixation. With its potential for carbon sequestration, it emerges as a promising candidate for environmental sustainability. Furthermore, it serves as a valuable C4 plant model, offering insights into the genetic foundations of extreme drought tolerance. Despite the availability of plastid and nuclear genomes, the absence of a mitochondrial genome (mitogenome or mtDNA) hinders a comprehensive understanding of its its mtDNA structure, organization, and phylogenetic implications.

RESULTS

In the present study, the mitochondrial genome of H. ammodendron was assembled and annotated, resulting in a multi-chromosomal configuration with two circular chromosomes. The mtDNA measured 210,149 bp in length and contained 31 protein-coding genes, 18 tRNA and three rRNA. Our analysis identified a total of 66 simple sequence repeats along with 27 tandem repeats, 312 forward repeats, and 303 palindromic repeats were found. Notably, 17 sequence fragments displayed homology between the mtDNA and chloroplast genome (cpDNA), spanning 5233 bp, accounting for 2.49% of the total mitogenome size. Additionally, we predicted 337 RNA editing sites, all of the C-to-U conversion type. Phylogenetic inference confidently placed H. ammodendron in the Amaranthacea family and its close relative, Suaeda glacum.

CONCLUSIONS

H. ammodendron mtDNA showed a multi-chromosomal structure with two fully circularized molecules. This newly characterized mtDNA represents a valuable resource for gaining insights into the basis of mtDNA structure variation within Caryophyllales and the evolution of land plants, contributing to their identification, and classification.

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.

Complete mitochondrial genome of Syzygium samarangense reveals genomic recombination, gene transfer, and RNA editing events

Wax apple (Syzygium samarangense) is a commercial fruit that belongs to one of the most species-rich tree genera in the world. We report here the first complete S. samarangense mitogenome obtained using a hybrid assembly strategy. The mitogenome was a 530,242 bp circular molecule encoding 61 unique genes accounting for 7.99% of the full-length genome. Additionally, 167 simple sequence repeats, 19 tandem repeats, and 529 pairs of interspersed repeats were identified. Long read mapping and Sanger sequencing revealed the involvement of two forward repeats (35,843 bp and 22,925 bp) in mediating recombination. Thirteen homologous fragments in the chloroplast genome were identified, accounting for 1.53% of the mitogenome, and the longest fragment was 2,432 bp. An evolutionary analysis showed that S. samarangense underwent multiple genomic reorganization events and lost at least four protein-coding genes (PCGs) (rps2, rps7, rps11, and rps19). A total of 591 RNA editing sites were predicted in 37 PCGs, of which nad1-2, nad4L-2, and rps10-2 led to the gain of new start codons, while atp6-1156, ccmFC-1315 and rps10-331 created new stop codons. This study reveals the genetic features of the S. samarangense mitogenome and provides a scientific basis for further studies of traits with an epistatic basis and for germplasm identification.

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.

Multiple configurations of the plastid and mitochondrial genomes of Caragana spinosa

MAIN CONCLUSION

In this study, we assembled the complete plastome and mitogenome of Caragana spinosa and explored the multiple configurations of the organelle genomes. Caragana spinosa belongs to the Papilionoidea subfamily and has significant pharmaceutical value. To explore the possible interaction between the organelle genomes, we assembled and analyzed the plastome and mitogenome of C. spinosa using the Illumina and Nanopore DNA sequencing data. The plastome of C. spinosa was 129,995 bp belonging to the inverted repeat lacking clade (IRLC), which contained 77 protein-coding genes, 29 tRNA genes, and four rRNA genes. The mitogenome was 378,373 bp long and encoded 54 unique genes, including 33 protein-coding, three ribosomal RNA (rRNA), and 18 transfer RNA (tRNA) genes. In addition to the single circular conformation, alternative conformations mediated by one and four repetitive sequences in the plastome and mitogenome were identified and validated, respectively. The inverted repeat (PDR12, the 12th dispersed repeat sequence in C. spinosa plastome) of plastome mediating recombinant was conserved in the genus Caragana. Furthermore, we identified 14 homologous fragments by comparing the sequences of mitogenome and plastome, including eight complete tRNA genes. A phylogenetic analysis of protein-coding genes extracted from the plastid and mitochondrial genomes revealed congruent topologies. Analyses of sequence divergence found one intergenic region, trnN-GUU-ycf1, exhibiting a high degree of variation, which can be used to develop novel molecular markers to distinguish the nine Caragana species accurately. This plastome and mitogenome of C. spinosa could provide critical information for the molecular breeding of C. spinosa and be used as a reference genome for other species of Caragana. In this study, we assembled the complete plastome and mitogenome of Caragana spinosa and explored the multiple configurations of the organelle genomes.

Comparative Analysis of the Complete Mitochondrial Genomes of Apium graveolens and Apium leptophyllum Provide Insights into Evolution and Phylogeny Relationships

The genus Apium, belonging to the family Apiaceae, comprises roughly 20 species. Only two species, Apium graveolens and Apium leptophyllum, are available in China and are both rich in nutrients and have favorable medicinal properties. However, the lack of genomic data has severely constrained the study of genetics and evolution in Apium plants. In this study, Illumina NovaSeq 6000 and Nanopore sequencing platforms were employed to identify the mitochondrial genomes of A. graveolens and A. leptophyllum. The complete lengths of the mitochondrial genomes of A. graveolens and A. leptophyllum were 263,017 bp and 260,164 bp, respectively, and contained 39 and 36 protein-coding genes, five and six rRNA genes, and 19 and 20 tRNA genes. Consistent with most angiosperms, both A. graveolens and A. leptophyllum showed a preference for codons encoding leucine (Leu). In the mitochondrial genome of A. graveolens, 335 SSRs were detected, which is higher than the 196 SSRs found in the mitochondrial genome of A. leptophyllum. Studies have shown that the most common RNA editing type is C-to-U, but, in our study, both A. graveolens and A. leptophyllum exhibited the U-C editing type. Furthermore, the transfer of the mitochondrial genomes of A. graveolens and A. leptophyllum into the chloroplast genomes revealed homologous sequences, accounting for 8.14% and 4.89% of the mitochondrial genome, respectively. Lastly, in comparing the mitochondrial genomes of 29 species, it was found that A. graveolens, A. leptophyllum, and Daucus carota form a sister group with a support rate of 100%. Overall, this investigation furnishes extensive insights into the mitochondrial genomes of A. graveolens and A. leptophyllum, thereby enhancing comprehension of the traits and evolutionary patterns within the Apium genus. Additionally, it offers supplementary data for evolutionary and comparative genomic analyses of other species within the Apiaceae family.

Factors contributing to mitogenome size variation and a recurrent intracellular DNA transfer in Melastoma

BACKGROUND

Mitogenome sizes of seed plants vary substantially even among closely related species, which are often related to horizontal or intracellular DNA transfer (HDT or IDT) events. However, the mechanisms of this size variation have not been well characterized.

RESULTS

Here we assembled and characterized the mitogenomes of three species of Melastoma, a tropical shrub genus experiencing rapid speciation. The mitogenomes of M. candidum (Mc), M. sanguineum (Ms) and M. dodecandrum (Md) were assembled to a circular mapping chromosome of 391,595 bp, 395,542 bp and 412,026 bp, respectively. While the mitogenomes of Mc and Ms showed good collinearity except for a large inversion of ~ 150 kb, there were many rearrangements in the mitogenomes between Md and either Mc or Ms. Most non-alignable sequences (> 80%) between Mc and Ms are from gain or loss of mitochondrial sequences. Whereas, between Md and either Mc or Ms, non-alignable sequences in Md are mainly chloroplast derived sequences (> 30%) and from putative horizontal DNA transfers (> 30%), and those in both Mc and Ms are from gain or loss of mitochondrial sequences (> 80%). We also identified a recurrent IDT event in another congeneric species, M. penicillatum, which has not been fixed as it is only found in one of the three examined populations.

CONCLUSIONS

By characterizing mitochondrial genome sequences of Melastoma, our study not only helps understand mitogenome size evolution in closely related species, but also cautions different evolutionary histories of mitochondrial regions due to potential recurrent IDT events in some populations or species.

Comprehensive analysis of the complete mitochondrial genomes of three Coptis species (C. chinensis, C. deltoidea and C. omeiensis): the important medicinal plants in China

Coptis plants (Ranunculaceae) contain high levels of isoquinoline alkaloids and have a long history of medicinal use. Coptis species are of great value in pharmaceutical industries and scientific research. Mitochondria are considered as one of the central units for receiving stress signals and arranging immediate responses. Comprehensive characterizations of plant mitogenomes are imperative for revealing the relationship between mitochondria, elucidating biological functions of mitochondria and understanding the environmental adaptation mechanisms of plants. Here, the mitochondrial genomes of C. chinensis, C. deltoidea and C. omeiensis were assembled through the Nanopore and Illumina sequencing platform for the first time. The genome organization, gene number, RNA editing sites, repeat sequences, gene migration from chloroplast to mitochondria were compared. The mitogenomes of C. chinensis, C. deltoidea and C. omeiensis have six, two, two circular-mapping molecules with the total length of 1,425,403 bp, 1,520,338 bp and 1,152,812 bp, respectively. The complete mitogenomes harbors 68-86 predicted functional genes including 39-51 PCGs, 26-35 tRNAs and 2-5 rRNAs. C. deltoidea mitogenome host the most abundant repeat sequences, while C. chinensis mitogenome has the largest number of transferred fragments from its chloroplasts. The large repeat sequences and foreign sequences in the mitochondrial genomes of Coptis species were related to substantial rearrangements, changes in relative position of genes and multiple copy genes. Further comparative analysis illustrated that the PCGs under selected pressure in mitochondrial genomes of the three Coptis species mainly belong to the mitochondrial complex I (NADH dehydrogenase). Heat stress adversely affected the mitochondrial complex I and V, antioxidant enzyme system, ROS accumulation and ATP production of the three Coptis species. The activation of antioxidant enzymes, increase of T-AOC and maintenance of low ROS accumulation in C. chinensis under heat stress were suggested as the factors for its thermal acclimation and normal growth at lower altitudes. This study provides comprehensive information on the Coptis mitogenomes and is of great importance to elucidate the mitochondrial functions, understand the different thermal acclimation mechanisms of Coptis plants, and breed heat-tolerant varieties.

When and why are mitochondria paternally inherited?

In contrast with nuclear genes that are passed on through both parents, mitochondrial genes are maternally inherited in most species, most of the time. The genetic conflict stemming from this transmission asymmetry is well-documented, and there is an abundance of population-genetic theory associated with it. While occasional or aberrant paternal inheritance occurs, there are only a few cases where exclusive paternal inheritance of mitochondrial genomes is the evolved state. Why this is remains poorly understood. By examining commonalities between species with exclusive paternal inheritance, we discuss what they may tell us about the evolutionary forces influencing mitochondrial inheritance patterns. We end by discussing recent technological advances that make exploring the causes and consequences of paternal inheritance feasible.

Comparative analysis of the organelle genomes of three Rhodiola species provide insights into their structural dynamics and sequence divergences

BACKGROUND

Plant organelle genomes are a valuable resource for evolutionary biology research, yet their genome architectures, evolutionary patterns and environmental adaptations are poorly understood in many lineages. Rhodiola species is a type of flora mainly distributed in highland habitats, with high medicinal value. Here, we assembled the organelle genomes of three Rhodiola species (R. wallichiana, R. crenulata and R. sacra) collected from the Qinghai-Tibet plateau (QTP), and compared their genome structure, gene content, structural rearrangements, sequence transfer and sequence evolution rates.

RESULTS

The results demonstrated the contrasting evolutionary pattern between plastomes and mitogenomes in three Rhodiola species, with the former possessing more conserved genome structure but faster evolutionary rates of sequence, while the latter exhibiting structural diversity but slower rates of sequence evolution. Some lineage-specific features were observed in Rhodiola mitogenomes, including chromosome fission, gene loss and structural rearrangement. Repeat element analysis shows that the repeats occurring between the two chromosomes may mediate the formation of multichromosomal structure in the mitogenomes of Rhodiola, and this multichromosomal structure may have recently formed. The identification of homologous sequences between plastomes and mitogenomes reveals several unidirectional protein-coding gene transfer events from chloroplasts to mitochondria. Moreover, we found that their organelle genomes contained multiple fragments of nuclear transposable elements (TEs) and exhibited different preferences for TEs insertion type. Genome-wide scans of positive selection identified one gene matR from the mitogenome. Since the matR is crucial for plant growth and development, as well as for respiration and stress responses, our findings suggest that matR may participate in the adaptive response of Rhodiola species to environmental stress of QTP.

CONCLUSION

The study analyzed the organelle genomes of three Rhodiola species and demonstrated the contrasting evolutionary pattern between plastomes and mitogenomes. Signals of positive selection were detected in the matR gene of Rhodiola mitogenomes, suggesting the potential role of this gene in Rhodiola adaptation to QTP. Together, the study is expected to enrich the genomic resources and provide valuable insights into the structural dynamics and sequence divergences of Rhodiola species.

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

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

Assembly and comparative analysis of the complete mitochondrial genome of Salix wilsonii using PacBio HiFi sequencing

Salix L. (willows) is one of the most taxonomically complex genera of flowering plants, including shrubs, tall trees, bushes, and prostrate plants. Despite the high species diversity, only five mitochondrial genomes (mitogenomes) have been released in this genus. Salix wilsonii is an important ornamental and economic willow tree in section Wilsonia of the genus Salix. In this study, the S. wilsonii mitogenome was assembled into a typical circular structure with a size of 711,456 bp using PacBio HiFi sequencing. A total of 58 genes were annotated in the S. wilsonii mitogenome, including 33 protein-coding genes (PCGs), 22 tRNAs, and 3 rRNAs. In the S. wilsonii mitogenome, four genes (mttB, nad3, nad4, and sdh4) were found to play important roles in its evolution through selection pressure analysis. Collinearity analysis of six Salix mitogenomes revealed high structural variability. To determine the evolutionary position of S. wilsonii, we conducted a phylogenetic analysis of the mitogenomes of S. wilsonii and 12 other species in the order Malpighiales. Results strongly supported the segregation of S. wilsonii and other five Salix species with 100% bootstrap support. The comparative analysis of the S. wilsonii mitogenome not only sheds light on the functional and structural features of S. wilsonii but also provides essential information for genetic studies of the genus Salix.

The chloroplast genome of Spiraea thunbergii (Rosaceae)

Spiraea thunbergii (S. thunbergia) is a very common ornamental shrub in the garden, with important horticultural and economic value. In this study, we assembled the complete chloroplast (cp) genome of S. thunbergia into a typical quadripartite structure. The genome size and GC content of the S. thunbergia cp genome are 155,922 bp and 36.76%, respectively. It contains 84 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The phylogenetic tree supported that S. thunbergii is closely related to Spiraea mongolica in the Rosaceae family. The study will provide significant genomic resource for elucidating the phylogenetic relationship of Spiraea.