Assembly and comparative analysis of the complete mitochondrial genome of Bupleurum chinense DC

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

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

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

MAIN CONCLUSION

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

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

Introduction

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

Methods

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

Results

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

Discussion

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

Assembly and comparative analysis of the complete mitochondrial genome of Lactuca sativa var. ramosa Hort

Lettuce (Lactuca sativa var. ramosa Hort) is an important leaf vegetable that widely cultivates due to its high-quality, short growth cycle, and less diseases. L. sativa var. ramosa Hort belongs to Asteraceae family and its evolutionary relationships with related species of Asteraceae are not completely assessed based on genome sequences. Here, we assembled the whole mitochondrial (mt) genome of L. sativa var. ramosa Hort, and performed a comparative with other related species. The L. sativa var. ramosa Hort mt genome has a typical circular structure with a length of 363,324 bp, within GC content accounted for 45.35%. In total of 71 genes, comprising 35 protein-coding genes (PCGs), 6 rRNAs, 28 tRNAs, and 2 pseudogenes were annotated. Codon preference, RNA-editing sites, repetitive sequences, and genes migrating from chloroplast (cp) to mt genomes were investigated in the L. sativa var. ramosa Hort mt genome. Nucleotide diversity (Pi) showed that the L. sativa var. ramosa Hort mt genome was relatively conserved. A Bayesian phylogenetic tree showed that L. sativa var. ramosa Hort was closely to L. sativa var. capitata and L. virosa, which belonged to the Lactuca genus in the Asteraceae family. Our findings will provide useful information to explore genetic variation, genetic diversity, and molecular breeding on the Lactuca genus.

Mitochondrial Genome Characteristics and Comparative Genomic Analysis of Spartina alterniflora

The mitochondrial genome of Spartina alterniflora, an invasive species with significant ecological and economic impacts, was analyzed to provide a theoretical basis for understanding its phylogenetic relationships and molecular biology. Mitochondrial genome sequences of S. alterniflora and 23 related species from NCBI were utilized for bioinformatics and comparative genomic analyses. A sliding window analysis identified three genes (rps2, atp9, and nad6) as potential DNA barcodes for species identification. Intracellular gene transfer (IGT) events between mitochondrial and chloroplast genome were detected, highlighting the dynamic nature of genomic evolution. A selective pressure analysis revealed that most protein-coding genes (PCGs) underwent purifying selection (Ka/Ks < 1), while the nad2 and ccmB genes showed signs of positive selection pressure (Ka/Ks > 1), indicating their role in adaptation. A phylogenetic analysis demonstrated a close relationship between S. alterniflora and Eleusine indica, supported by a collinearity analysis, which suggests environmental convergence. This study provides novel insights into the structural and evolutionary characteristics of the S. alterniflora mitochondrial genome, offering valuable genomic resources for future research on invasive species management and evolutionary biology.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Complete mitochondrial genome assembly and comparative analysis of Colocasia esculenta

Colocasia esculenta ranks as the fifth most important tuber crop and is known for its high nutritional and medicinal value. However, there is no research on its mitochondrial genome, hindering in-depth exploration of its genomic resources and genetic relationships. Using second- and third-generation sequencing technologies, we assembled and annotated the mitogenome of C. esculenta. Its mitogenome mainly consists of five circular DNA molecules, with a total length of 594,811 bp and a GC content of 46.25%. A total of 55 genes, 157 simple sequence repeats, 29 tandem repeat sequences, 202 dispersed repeat sequences, and 625 RNA editing sites were detected. Most protein-coding genes use ATG as the start codon, and the third position of the codon tends to be A or T (U). GAA, AUU, and UUU are the most common codons in C. esculenta mitochondria. Finally, based on 28 representative plant species, a phylogenetic tree was constructed, revealing a close relationship between C. esculenta and Araceae. This study provides comprehensive information on C. esculenta, laying a foundation for crop genetics and molecular breeding.

Mitochondrial Genome Insights into Evolution and Gene Regulation in Phragmites australis

As a globally distributed perennial Gramineae, Phragmites australis can adapt to harsh ecological environments and has significant economic and environmental values. Here, we performed a complete assembly and annotation of the mitogenome of P. australis using genomic data from the PacBio and BGI platforms. The P. australis mitogenome is a multibranched structure of 501,134 bp, divided into two circular chromosomes of 325,493 bp and 175,641 bp, respectively. A sequence-simplified succinate dehydrogenase 4 gene was identified in this mitogenome, which is often translocated to the nuclear genome in the mitogenomes of gramineous species. We also identified tissue-specific mitochondrial differentially expressed genes using RNAseq data, providing new insights into understanding energy allocation and gene regulatory strategies in the long-term adaptive evolution of P. australis mitochondria. In addition, we studied the mitogenome features of P. australis in more detail, including repetitive sequences, gene Ka/Ks analyses, codon preferences, intracellular gene transfer, RNA editing, and multispecies phylogenetic analyses. Our results provide an essential molecular resource for understanding the genetic characterisation of the mitogenome of P. australis and provide a research basis for population genetics and species evolution in Arundiaceae.

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.

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

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

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Complete Mitogenome Assembly and Comparative Analysis of Vaccinium bracteatum (Ericaceae), a Rich Source of Health-Promoting Molecules

Vaccinium bracteatum is a valuable plant used both as food and medicine in China, but low production limits the development of its industry. As such, it is important to develop genetic resources for the high-value species for preservation of wild populations and utilization. The complete chloroplast and nuclear genomes have already been available; however, its mitogenome has not yet been characterized. Here, the V. bracteatum mitogenome was assembled using HiFi reads, and a comparative analysis was conducted. The mitogenome was a circular sequence of 708,384 bp with a GC content of 45.28%, in which 67 genes were annotated, including 36 protein-coding genes, 26 tRNA genes, 3 rRNA genes, and 2 pseudogenes. Overall, 370 dispersed repeats, 161 simple repeats, and 42 tandem repeats were identified, and 360 RNA editing sites were predicted. There was extensive DNA migration among the three genomes. In addition, most of the protein-coding genes underwent purifying selection throughout evolution, and the nucleotide diversity was highly variable. In addition, comparative analysis indicated that the sizes, structures, and gene contents of the mitogenomes differed significantly, but the GC contents and functional genes were relatively conserved among the Ericales species. Mitogenome-based phylogenetic analysis indicated the precise. evolutionary and taxonomic status of V. bracteatum. The complete mitogenome represents the last link of the reference genome of V. bracteatum and lays the foundation for effective utilization and molecular breeding of this plant.

Assembly and comparative analysis of the multichromosomal mitochondrial genome of globally endangered seagrass Halophila beccarii

BACKGROUND

Halophila beccarii is one of the oldest two generations of seagrass plants and one of the 10 species of seagrass currently at risk of extinction worldwide. Therefore, how to effectively protect the H. beccarii resources from extinction is a huge challenge. Molecular biology research can provide a scientific basis for species conservation. So far, there has been no detailed analysis of the mitochondrial genome of the genus Halophila.

RESULTS

The mitochondrial genome of H. beccarii was assembled into 28 circular chromosomes, ranging in length from 41,738 bp to 104,744 bp, with a total length of 1,964,072 bp and a GC content of 46.71%. It contains 39 genes, including 26 protein coding genes, 10 tRNA genes, and 3 rRNA genes. Repeat sequence analysis and prediction of RNA editing sites revealed a total of 850 dispersed repeats, 1,205 simple repeats, 61 tandem repeats, and 120 RNA editing sites. Analysis of codon usage indicates that codons ending in A/U are preferred. Gene migration between the mitochondrial genome and the chloroplast genome was observed through homologous fragment detection. In addition, Ka/Ks analysis showed that most protein coding genes in the mitochondrial genome experienced negative selection, while only the nad3 gene experienced potential positive selection in most Alismatales. Nucleotide polymorphism analysis revealed variations in each gene, with rpl10 being the most significant. In addition, comparative analysis shows that the GC content is conserved, but there are significant differences in the size and structure of mitochondrial genomes among different species of Alismatales. The phylogenetic analysis based on the mitochondrial genome reflects the exact evolutionary and taxonomic status of H. beccarii.

CONCLUSION

In this study, we sequenced and annotated the mitochondrial genome of H. beccarii, and compared it with the mitochondrial genomes of other plants in Alismatales. Our findings enrich the mitogenome database of seagrass plants and highlight the potential for mitochondrial genes to help decipher plant evolutionary history.

Comparative analysis of mitochondrial genomes of invasive weed Mikania micrantha and its indigenous congener Mikania cordata

Mikania micrantha and Mikania cordata are two distinct species in China. The former is notorious as one of the top 100 worst invasive species, whereas the latter is an indigenous species harmless to native plants or the environment. They form an ideal congener pair for comparative studies aimed at deeply understanding the invasion mechanisms of the exotic weed. In this study, we have assembled and annotated the mitogenomes of both species using Illumina and PacBio sequencing data and compared their characteristic differences. The complete mitogenome of M. micrantha is a double-stranded DNA with a length of 336,564 bp, while the mitogenome of M. cordata exhibits a branching structure, consisting of two small circular molecules and six linear molecules, with a combined length totaling 335,444 bp. Compared to M. cordata, M. micrantha has less SSRs, tandem repeats, dispersed repeats, mitochondrial protein coding genes (PCGs). The two plants show similar codon usage patterns. This comparative study has revealed the structure and function of the mitogenomes of the two species and laid a solid foundation for investigating the effects of gene loss and duplication on the development of invasive traits in M. micrantha.

Assembly and comparative analysis of the complete mitochondrial genome of Vaccinium carlesii Dunn

Vaccinium L. is an important fruit tree with nutritional, medicinal, and ornamental values. However, the mitochondrial (mt) genome of Vaccinium L. remains largely unexplored. Vaccinium carlesii Dunn is an endemic wild resource in China, which is crucial for blueberry breeding. The V. carlesii mt genomes were sequenced using Illumina and Nanopore, which total length was 636,904 bp with 37 protein coding genes, 20 tRNA genes, and three rRNA genes. We found four pairs of long repeat fragments homologous recombination mediated the generation of substructures in the V. carlesii mt genome. We predicted 383 RNA editing sites, all converting cytosine (C) to uracil (U). According to the phylogenetic analysis, V. carlesii and V. macrocarpon of the Ericaceae exhibited the closest genetic relationship. This study provides a theoretical basis for understanding the evolution of higher plants, species classification and identification, and will also be useful for further utilization of Vaccinium germplasm resources.

Assembly and Comparative Analysis of the Complete Mitochondrial Genome of Saussurea inversa (Asteraceae)

Saussurea inversa is a perennial herb used in traditional Chinese medicine and is effective against rheumatoid arthritis. In this study, we sequenced the complete mitochondrial (mt) genome of S. inversa (GenBank accession number: ON584565.1). The circular mt genome of S. inversa was 335,372 bp in length, containing 62 genes, including 33 mRNAs, 22 tRNAs, 6 rRNAs, and 1 pseudogene, along with 1626 open reading frames. The GC content was 45.14%. Predictive analysis revealed substantial RNA editing, with ccmFn being the most abundantly edited gene, showing 36 sites. Gene migration between the mt and chloroplast (cp) genomes of S. inversa was observed through the detection of homologous gene fragments. Phylogenetic analysis revealed that S. inversa was clustered with Arctium tomentosum (Asteraceae). Our findings provide extensive information regarding the mt genome of S. inversa and help lay the foundation for future studies on its genetic variations, phylogeny, and breeding via the analysis of the mt genome.

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.

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.

The Whole Mitochondrial Genome Sequence of Dendrobium loddigesii Rolfe, an Endangered Orchid Species in China, Reveals a Complex Multi-Chromosome Structure

Dendrobium loddigesii is a precious traditional Chinese medicine with high medicinal and ornamental value. However, the characterization of its mitochondrial genome is still pending. Here, we assembled the complete mitochondrial genome of D. loddigesii and discovered that its genome possessed a complex multi-chromosome structure. The mitogenome of D. loddigesii consisted of 17 circular subgenomes, ranging in size from 16,323 bp to 56,781 bp. The total length of the mitogenome was 513,356 bp, with a GC content of 43.41%. The mitogenome contained 70 genes, comprising 36 protein-coding genes (PCGs), 31 tRNA genes, and 3 rRNA genes. Furthermore, we detected 403 repeat sequences as well as identified 482 RNA-editing sites and 8154 codons across all PCGs. Following the sequence similarity analysis, 27 fragments exhibiting homology to both the mitogenome and chloroplast genome were discovered, accounting for 9.86% mitogenome of D. loddigesii. Synteny analysis revealed numerous sequence rearrangements in D. loddigesii and the mitogenomes of related species. Phylogenetic analysis strongly supported that D. loddigesii and D. Amplum formed a single clade with 100% bootstrap support. The outcomes will significantly augment the orchid mitochondrial genome database, offering profound insights into Dendrobium's intricate mitochondrial genome architecture.

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.

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.

Assembly and comparative analysis of the complete mitochondrial genome of Brassica rapa var. Purpuraria

BACKGROUND

Purple flowering stalk (Brassica rapa var. purpuraria) is a widely cultivated plant with high nutritional and medicinal value and exhibiting strong adaptability during growing. Mitochondrial (mt) play important role in plant cells for energy production, developing with an independent genetic system. Therefore, it is meaningful to assemble and annotate the functions for the mt genome of plants independently. Though there have been several reports referring the mt genome of in Brassica species, the genome of mt in B. rapa var. purpuraria and its functional gene variations when compared to its closely related species has not yet been addressed.

RESULTS

The mt genome of B. rapa var. purpuraria was assembled through the Illumina and Nanopore sequencing platforms, which revealed a length of 219,775 bp with a typical circular structure. The base composition of the whole B. rapa var. purpuraria mt genome revealed A (27.45%), T (27.31%), C (22.91%), and G (22.32%). 59 functional genes, composing of 33 protein-coding genes (PCGs), 23 tRNA genes, and 3 rRNA genes, were annotated. The sequence repeats, codon usage, RNA editing, nucleotide diversity and gene transfer between the cp genome and mt genome were examined in the B. rapa var. purpuraria mt genome. Phylogenetic analysis show that B. rapa var. Purpuraria was closely related to B. rapa subsp. Oleifera and B. juncea. Ka/Ks analysis reflected that most of the PCGs in the B. rapa var. Purpuraria were negatively selected, illustrating that those mt genes were conserved during evolution.

CONCLUSIONS

The results of our findings provide valuable information on the B.rapa var. Purpuraria genome, which might facilitate molecular breeding, genetic variation and evolutionary researches for Brassica species in the future.

Assembly and comparative analysis of the complete mitochondrial genome of white towel gourd (Luffa cylindrica)

Mitochondria play an important role in the energy production of plant cells through independent genetic systems. This study has aimed to assemble and annotate the functions of the mitochondrial (mt) genome of Luffa cylindrica. The mt genome of L. cylindrica contained two chromosomes with lengths of 380,879 bp and 67,982 bp, respectively. Seventy-seven genes including 39 protein-coding genes, 34 tRNA genes, 3 rRNA genes, and 1 pseudogene, were identified. About 90.63% of the codons ended with A or U bases, and 98.63% of monomers contained A/T, which contributed to the high A/T content (55.91%) of the complete mt genome. Six genes (ATP8, CCMFC, NAD4, RPL10, RPL5 and RPS4) showed positive selection. Phylogenetic analysis indicates that L. cylindrica is closely related to L. acutangula. The present results provide the mt genome of L. cylindrica, which may facilitate possible genetic variation, evolutionary, and molecular breeding studies of L. cylindrica.

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Assembly and comparative analysis of the complete mitochondrial genome of Trigonella foenum-graecum L

BACKGROUND

Trigonella foenum-graecum L. is a Leguminosae plant, and the stems, leaves, and seeds of this plant are rich in chemical components that are of high research value. The chloroplast (cp) genome of T. foenum-graecum has been reported, but the mitochondrial (mt) genome remains unexplored.

RESULTS

In this study, we used second- and third-generation sequencing methods, which have the dual advantage of combining high accuracy and longer read length. The results showed that the mt genome of T. foenum-graecum was 345,604 bp in length and 45.28% in GC content. There were 59 genes, including: 33 protein-coding genes (PCGs), 21 tRNA genes, 4 rRNA genes and 1 pseudo gene. Among them, 11 genes contained introns. The mt genome codons of T. foenum-graecum had a significant A/T preference. A total of 202 dispersed repetitive sequences, 96 simple repetitive sequences (SSRs) and 19 tandem repetitive sequences were detected. Nucleotide diversity (Pi) analysis counted the variation in each gene, with atp6 being the most notable. Both synteny and phylogenetic analyses showed close genetic relationship among Trifolium pratense, Trifolium meduseum, Trifolium grandiflorum, Trifolium aureum, Medicago truncatula and T. foenum-graecum. Notably, in the phylogenetic tree, Medicago truncatula demonstrated the highest level of genetic relatedness to T. foenum-graecum, with a strong support value of 100%. The interspecies non-synonymous substitutions (Ka)/synonymous substitutions (Ks) results showed that 23 PCGs had Ka/Ks < 1, indicating that these genes would continue to evolve under purifying selection pressure. In addition, setting the similarity at 70%, 23 homologous sequences were found in the mt genome of T. foenum-graecum.

CONCLUSIONS

This study explores the mt genome sequence information of T. foenum-graecum and complements our knowledge of the phylogenetic diversity of Leguminosae plants.

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

Introduction

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

Methods

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

Results and Discussion

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

Assembly and analysis of the complete mitochondrial genome of Forsythia suspensa (Thunb.) Vahl

BACKGROUND

Forsythia suspensa (Thunb.) Vahl is a valuable ornamental and medicinal plant. Although the nuclear and chloroplast genomes of F. suspensa have been published, its complete mitochondrial genome sequence has yet to be reported. In this study, the genomic DNA of F. suspensa yellowish leaf material was extracted, sequenced by using a mixture of Illumina Novaseq6000 short reads and Oxford Nanopore PromethION long reads, and the sequencing data were assembled and annotated.

RESULT

The F. suspensa mitochondrial genome was obtained in the length of 535,692 bp with a circular structure, and the GC content was 44.90%. The genome contains 60 genes, including 36 protein-coding genes, 21 tRNA genes, and three rRNA genes. We further analyzed RNA editing of the protein-coding genes, relative synonymous codon usage, and sequence repeats based on the genomic data. There were 25 homologous sequences between F. suspensa mitochondria and chloroplast genome, which involved the transfer of 8 mitochondrial genes, and 9473 homologous sequences between mitochondrial and nuclear genomes. Analysis of the nucleic acid substitution rate, nucleic acid diversity, and collinearity of protein-coding genes of the F. suspensa mitochondrial genome revealed that the majority of genes may have undergone purifying selection, exhibiting a slower rate of evolution and a relatively conserved structure. Analysis of the phylogenetic relationships among different species revealed that F. suspensa was most closely related to Olea europaea subsp. Europaea.

CONCLUSION

In this study, we sequenced, assembled, and annotated a high-quality F. suspensa mitochondrial genome. The results of this study will enrich the mitochondrial genome data of Forsythia, lay a foundation for the phylogenetic development of Forsythia, and promote the evolutionary analysis of Oleaceae species.

Multichromosomal mitochondrial genome of Punica granatum: comparative evolutionary analysis and gene transformation from chloroplast genomes

BACKGROUND

Punica granatum is a fundamentally important fruit tree that has important economic, medicinal and ornamental properties. At present, there are few reports on the mitochondrial genome of pomegranate. Hence, in this study the P. granatum mitogenome was sequenced and assembled to further understanding of organization, variation, and evolution of mitogenomes of this tree species.

RESULTS

The genome structure was multi-chromosomes with seven circular contigs, measuring 382,774 bp in length with a 45.91% GC content. It contained 74 genes, including 46 protein-coding genes, 25 tRNA genes, and three rRNA genes. There were 188 pairs of dispersed repeats with lengths of 30 or greater, primarily consisting of reverse complementary repeats. The mitogenome analysis identified 114SSRs and 466 RNA editing sites. Analyses of codon usage, nucleotide diversity and gene migration from chloroplast to mitochondrial were also conducted. The collinear and comparative analysis of mitochondrial structures between P. granatum and its proximal species indicated that P. granatum 'Taishanhong' was closely related to P. granatum 'Qingpitian' and Lagerstroemia indica. Phylogenetic examination based on the mitogenome also confirmed the evolutionary relationship.

CONCLUSION

The results offered crucial information on the evolutionary biology of pomegranate and highlighted ways to promote the utilization of the species' germplasm.

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.

Assembly and phylogenetic analysis of the mitochondrial genome of endangered medicinal plant Huperzia crispata

Huperzia crispata is a traditional Chinese herb plant and has attracted special attention in recent years for its products Hup A can serve as an acetylcholinesterase inhibitor (AChEI). Although the chloroplast (cp) genome of H. crispata has been studied, there are no reports regarding the Huperzia mitochondrial (mt) genome since the previously reported H. squarrosa has been revised as Phlegmariurus squarrosus. The mt genome of H. crispata was sequenced using a combination of long-read nanopore and Illumina sequencing platforms. The entire H. crispata mt genome was assembled in a circular with a length of 412,594 bp and a total of 91 genes, including 45 tRNAs, 6 rRNAs, 37 protein-coding genes (PCGs), and 3 pseudogenes. Notably, the rps8 gene was present in P. squarrosus and a pseudogene rps8 was presented in H. crispata, which was lacking in most of Pteridophyta and Gymnospermae. Intron-encoded maturase (mat-atp9i85 and mat-cobi787) genes were present in H. crispata and P. squarrosus, but lost in other examined lycophytes, ferns, and Gymnospermae plants. Collinearity analysis showed that the mt genome of H. crispata and P. squarrossus is highly conservative compared to other ferns. Relative synonymous codon usage (RSCU) analysis showed that the amino acids most frequently found were phenylalanine (Phe) (4.77%), isoleucine (Ile) (4.71%), lysine (Lys) (4.26%), while arginine (Arg) (0.32%), and histidine (His) (0.42%) were rarely found. Simple sequence repeats (SSR) analysis revealed that a total of 114 SSRs were identified in the mt genome of H. crispata and account for 0.35% of the whole mt genome. Monomer repeats were the majority types of SSRs and represent 91.89% of the total SSRs. In addition, a total of 1948 interspersed repeats (158 forward, 147 palindromic, and 5 reverse repeats) with a length ranging from 30 bp to 14,945 bp were identified in the H. crispata mt genome and the 30-39-bp repeats were the most abundant type. Gene transfer analysis indicated that a total of 12 homologous fragments were discovered between the cp and mt genomes of H. crispata, accounting for 0.93% and 2.48% of the total cp and mt genomes, respectively. The phylogenetic trees revealed that H. crispata was the sister of P. squarrosus. The Ka/Ks analysis results suggested that most PCGs, except atp6 gene, were subject to purification selection during evolution. Our study provides extensive information on the features of the H. crispata mt genome and will help unravel evolutionary relationships, and molecular identification within lycophytes.

Assembly and comparative analysis of the complete mitochondrial genome of Ilex metabaptista (Aquifoliaceae), a Chinese endemic species with a narrow distribution

BACKGROUND

Ilex metabaptista is a woody tree species with strong waterlogging tolerance and is also admired as a landscape plant with high development prospects and scientific research value. Unfortunately, populations of this species have declined due to habitat loss. Thus, it is a great challenge for us to efficiently protect I. metabaptista resources from extinction. Molecular biology research can provide the scientific basis for the conservation of species. However, the study of I. metabaptista genetics is still in its infancy. To date, no mitochondrial genome (mitogenome) in the genus Ilex has been analysed in detail.

RESULTS

The mitogenome of I. metabaptista was assembled based on the reads from Illumina and Nanopore sequencing platforms; it was a typical circular DNA molecule of 529,560 bp with a GC content of 45.61% and contained 67 genes, including 42 protein-coding genes, 22 tRNA genes, and 3 rRNA genes. Repeat sequence analysis and prediction of RNA editing sites revealed a total of 286 dispersed repeats, 140 simple repeats, 18 tandem repeats, and 543 RNA editing sites. Analysis of codon usage showed that codons ending in A/T were preferred. Gene migration was observed to occur between the mitogenome and chloroplast genome via the detection of homologous fragments. In addition, Ka/Ks analysis revealed that most of the protein-coding genes in the mitogenome had undergone negative selection, and only the ccmB gene had undergone potential positive selection in most asterids. Nucleotide polymorphism analysis revealed the variation in each gene, with atp9 being the most notable. Furthermore, comparative analysis showed that the GC contents were conserved, but the sizes and structure of mitogenomes varied greatly among asterids. Phylogenetic analysis based on the mitogenomes reflected the exact evolutionary and taxonomic status of I. metabaptista.

CONCLUSION

In this study, we sequenced and annotated the mitogenome of I. metabaptista and compared it with the mitogenomes of other asterids, which provided essential background information for further understanding of the genetics of this plant and helped lay the foundation for future studies on molecular breeding of I. metabaptista.

Complete mitochondrial genome assembly and comparison of Camellia sinensis var. Assamica cv. Duntsa

Camellia sinensis var. Assamica cv. Duntsa (C.duntsa), a valuable Theaceae from Hunan Province, has been looked at as a precious tea resource by local farmers because of its economic and ecological value. Genomics study on C.duntsa is essential for the domestication and enhancement of tea tree varieties. In the present study, we used a hybrid approach based on Illumina and PacBio data to sequence and assemble the mitochondrial genome of C.duntsa. The mitochondrial genome of C.duntsa was estimated to be 1,081,996 base pairs (bp) and eighty-one genes consisting of one pseudogene, three ribosomal RNA (rRNA) genes, thirty transfer RNA (tRNA) genes, and forty-seven protein-coding genes (PCGs). Tetramer repetitions made up 43.90% of simple sequence repeats (SSRs). The codon usage bias of the Theaceae mitochondrial gene atp9 was altered by mutation, but the codon usage of other genes was shaped by natural selection. Besides, there are eighteen gene-containing homologous regions between the chloroplast and mitochondrial genomes of C. duntsa.Some genomes including atp8, cox1, cox3, nad7, nad9, rpl16, rpl2, rps19, rps4, and sdh4 are absent in the mitochondrial genome of several Theaceae plant. However, C. duntsa maintains these genes integrity and functionality. Another gene, rps16, is either lacking from the mitochondrial genome of C. duntsa or is present as a pseudogene. C. duntsa and C. sinensis (OM809792) are very similar, as shown by a collinear match across four species of Theaceae; the most conservative genes are nad5, atp9, cox2, rps3, trnA-TGC, trnI-GAT, rrn18, trnV-GAC, and ccmFN. Similarly, the genome's phylogenetic trees revealed that C. duntsa was the sister species to C. sinensis. The results confirmed that the C. duntsa and C. sinensis (OM809792) mitochondrial genome underwent gene rearrangement.In general, our results shows that genomic information from organelles can help us understand plant phylogeny and can also be used to make molecular markers and study how genetic traits change over time. Our research will contribute to the population genetics and evolution of tea plant.