Assembly and comparative analysis of complete mitochondrial genome sequence of an economic plant Salix suchowensis

Assemble and comparative analysis of the mitochondrial genome of Rhododendron delavayi: Insights into phylogenetic relationships and genomic variations

Rhododendron delavayi, a notable ornamental plant primarily found in regions of China like Yunnan and Guizhou provinces, holds substantial horticultural value. To elucidate the systematic phylogenetic relationships and organelle genomic differences within R. delavayi and related Rhododendron species, we conducted sequencing and assembly of the complete mitochondrial genome of R. delavayi. The full-length mitochondrial genome of it was a singular circular molecule spanning 1,009,263 bp, comprising 53 protein-coding genes, including 18 transfer RNA (tRNA) genes, 3 ribosomal RNA (rRNA) genes, and 32 protein-coding genes. A total of 1,182 simple sequence repeats (SSRs) loci were identified in the R. delavayi mitochondrial genome, primarily consisting of single nucleotide, dinucleotide, and trinucleotide repeats. Nucleotide diversity analysis highlighted five genes (atp6, atp9, cox2, nad1, and rpl10) with the highest diversity within the mitochondrial genomes of Rhododendron genus. Comparative analysis of the mitochondrial genome of R. delavayi with those of four other Rhododendron species indicated complex rearrangements in 21 genes, including rps4, nad6, rps3, atp6, cob, atp9, nad7, among others. The mitochondrial phylogenetic tree revealed a close relationship between R. delavayi and R. decorum, forming a sister clade to R. × pulchrum and R. simsii. Furthermore, 126 plastid-to-mitochondrial gene transfers in R. delavayi were identified, ranging from 30 bp to 19,385 bp. These fragments collectively constituted 47.54 % and 9.52 % of the chloroplast and mitochondrial genomes (202,169 bp), respectively. Complex mitochondrial-to-mitochondrial transfers were also observed, with 843 identified fragments totaling 312,036 bp (30.92 % of the mitochondrial genome). Segments exceeding 10 kb may mediate homologous recombination within the mitochondrial molecules. Remarkably, our study underscores that the mitochondrial genome of R. delavayi was the largest reported within the Rhododendron genus to date. The intricate rearrangements observed in the mitochondrial genomes of Rhododendron species, alone with the identification of five potential molecular marker sites, provided valuable insights for species classification and parentage identification within the Rhododendron genus.

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.

De novo assembly of the complete mitochondrial genome of pepino (Solanum muricatum) using PacBio HiFi sequencing: insights into structure, phylogenetic implications, and RNA editing

BACKGROUND

Solanum muricatum is an emerging horticultural fruit crop with rich nutritional and antioxidant properties. Although the chromosome-scale genome of this species has been sequenced, its mitochondrial genome sequence has not been reported to date.

RESULTS

PacBio HiFi sequencing was used to assemble the circular mitogenome of S. muricatum, which was 433,466 bp in length. In total, 38 protein-coding, 19 tRNA, and 3 rRNA genes were annotated. The reticulate mitochondrial conformations with multiple junctions were verified by polymerase chain reaction, and codon usage, sequence repeats, and gene migration from chloroplast to mitochondrial genome were determined. A collinearity analysis of eight Solanum mitogenomes revealed high structural variability. Overall, 585 RNA editing sites in protein coding genes were identified based on RNA-seq data. Among them, mttB was the most frequently edited (52 times), followed by ccmB (46 times). A phylogenetic analysis based on the S. muricatum mitogenome and those of 39 other taxa (including 25 Solanaceae species) revealed the evolutionary and taxonomic status of S. muricatum.

CONCLUSIONS

We provide the first report of the assembled and annotated S. muricatum mitogenome. This information will help to lay the groundwork for future research on the evolutionary biology of Solanaceae species. Furthermore, the results will assist the development of molecular breeding strategies for S. muricatum based on the most beneficial agronomic traits of this species.

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.

Assembly and analysis of the complete mitochondrial genome of the Chinese wild dwarf almond (Prunus tenella)

Background: The wild dwarf almond (Prunus tenella) is one of the national key grade II-protected wild plants in China. It is a relic deciduous forest species from the middle Eocene of the ancient Mediterranean Sea and is also known as a "living fossil of plants." It is distributed in Southeast Europe, West Asia, Central Asia, Siberia, and Xinjiang (Tacheng) and other areas of China. The plant grows on arid slopes, steppes, depressions, and valleys at an altitude of 1,200 m. The seeds of wild dwarf almonds are frost resistant and contain oil and bitter lentil glycosides, which possess medicinal value. Additionally, the seeds of wild dwarf almonds can be used as the original material for breeding new varieties of almonds and obtain ornamental flowers and trees. Results: The complete mitochondrial genome of P. tenella was sequenced and assembled using two sequencing platforms, namely, Illumina Novaseq6000 and Oxford Nanopore PromethION. The assembled genome was 452,158-bp long with a typical loop structure. The total number of A, T, C, and G bases in the genome was 122,066 (26.99%), 124,114 (27.45%), 103,285 (22.84%), and 102,693 (22.71%), respectively, with a GC content of 45.55%. A total of 63 unique genes, including 36 protein-coding genes, 24 tRNA genes, and 3 rRNA genes, were identified in the genome. Furthermore, codon usage, sequence duplication, RNA editing, and mitochondrial and chloroplast DNA fragment transfer events in the genome were analyzed. A phylogenetic tree was also constructed using 30 protein-coding genes that are common to the mitochondrial genomes of 24 species, which indicated that the genome of wild lentils is highly conserved with those of apples and pears belonging to Rosaceae. Conclusion: Assembly and annotation of the P. tenella mitochondrial genome provided comprehensive information about the mitochondrial genome of wild dwarf almonds, This study provides information on the mitochondrial genome of Prunus species and serves as a reference for further evolutionary studies on wild dwarf almonds.

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Assembly and comparative analysis of the complete mitochondrial genome of Isopyrum anemonoides (Ranunculaceae)

Ranunculaceae is a large family of angiosperms comprising 2500 known species-a few with medicinal and ornamental values. Despite this, only two mitochondrial genomes (mitogenomes) of the family have been released in GenBank. Isopyrum anemonoides is a medicinal plant belonging to the family Ranunculaceae, and its chloroplast genome has recently been reported; however, its mitogenome remains unexplored. In this study, we assembled and analyzed the complete mitochondrial genome of I. anemonoides and performed a comparative analysis against different Ranunculaceae species, reconstructing the phylogenetic framework of Isopyrum. The circular mitogenome of I. anemonoides has a length of 206,722 bp, with a nucleotide composition of A (26.4%), T (26.4%), C (23.6%), and G (23.6%), and contains 62 genes, comprising 37 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, and three ribosomal RNA (rRNA) genes. Abundantly interspersed repetitive and simple sequence repeat (SSR) loci were detected in the I. anemonoides mitogenome, with tetranucleotide repeats accounting for the highest proportion of SSRs. By detecting gene migration, we observed gene exchange between the chloroplast and mitogenome in I. anemonoides, including six intact tRNA genes, six PCG fragments, and fragments from two rRNA genes. Comparative mitogenome analysis of three Ranunculaceae species indicated that the PCG contents were conserved and the GC contents were similar. Selective pressure analysis revealed that only two genes (nad1 and rpl5) were under positive selection during their evolution in Ranunculales, and two specific RNA editing sites (atp6 and mttB) were detected in the I. anemonoides mitogenome. Moreover, a phylogenetic analysis based on the mitogenomes of I. anemonoides and the other 15 taxa accurately reflected the evolutionary and taxonomic status of I. anemonoides. Overall, this study provides new insights into the genetics, systematics, and evolution of mitochondrial evolution in Ranunculaceae, particularly I. anemonoides.

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.

The complete mitochondrial genome of Amorphophallus albus and development of molecular markers for five Amorphophallus species based on mitochondrial DNA

Introduction

Amorphophallus albus is an herbaceous, cormous, perennial plant used as a food source and traditional medicine in Asia.

Methods

In this study, we assembled and annotated the complete mitochondrial genome (mitogenome) of A. albus. Then we analyzed the repeated elements and mitochondrial plastid sequences (MTPTs), predicted RNA editing sites in mitochondrial protein-coding genes (PCGs). Lastly, we inferred the phylogenetic relationships of A. albus and other angiosperms based on mitochondrial PCGs, and designed two molecular markers based on mitochondrial DNA.

Results and discussion

The complete mitogenome of A. albus consists of 19 circular chromosomes. And the total length of A. albus mitogenome is 537,044 bp, with the longest chromosome measuring 56,458 bp and the shortest measuring 12,040 bp. We identified and annotated a total of 36 protein-coding genes (PCGs), 21 tRNA genes, and 3 rRNA genes in the mitogenome. Additionally, we analyzed mitochondrial plastid DNAs (MTPTs) and identified 20 MTPTs between the two organelle genomes, with a combined length of 22,421 bp, accounting for 12.76% of the plastome. Besides, we predicted a total of 676 C to U RNA editing sites on 36 protein-coding genes of high confidence using Deepred-mt. Furthermore, extensive genomic rearrangement was observed between A. albus and the related mitogenomes. We conducted phylogenetic analyses based on mitochondrial PCGs to determine the evolutionary relationships between A. albus and other angiosperms. Finally, we developed and validated two molecular markers, Ai156 and Ai976, based on two intron regions (nad2i156 and nad4i976) respectively. The discrimination success rate was 100 % in validation experiments for five widely grown konjac species. Our results reveal the multi-chromosome mitogenome of A. albus, and the developed markers will facilitate molecular identification of this genus.

The mitochondrial genome of the diploid oat Avena longiglumis

BACKGROUND

Avena longiglumis Durieu (2n = 2x = 14) is a wild relative of cultivated oat (Avena sativa, 2n = 6x = 42) with good agronomic and nutritional traits. The plant mitochondrial genome has a complex organization and carries genetic traits of value in exploiting genetic resources, not least male sterility alleles used to generate F₁ hybrid seeds. Therefore, we aim to complement the chromosomal-level nuclear and chloroplast genome assemblies of A. longiglumis with the complete assembly of the mitochondrial genome (mitogenome) based on Illumina and ONT long reads, comparing its structure with Poaceae species.

RESULTS

The complete mitochondrial genome of A. longiglumis can be represented by one master circular genome being 548,445 bp long with a GC content of 44.05%. It can be represented by linear or circular DNA molecules (isoforms or contigs), with multiple alternative configurations mediated by long (4,100-31,235 bp) and medium (144-792 bp) size repeats. Thirty-five unique protein-coding genes, three unique rRNA genes, and 11 unique tRNA genes are identified. The mitogenome is rich in duplications (up to 233 kb long) and multiple tandem or simple sequence repeats, together accounting for more than 42.5% of the total length. We identify homologous sequences between the mitochondrial, plastid and nuclear genomes, including the exchange of eight plastid-derived tRNA genes, and nuclear-derived retroelement fragments. At least 85% of the mitogenome is duplicated in the A. longiglumis nuclear genome. We identify 269 RNA editing sites in mitochondrial protein-coding genes including stop codons truncating ccmFC transcripts.

CONCLUSIONS

Comparative analysis with Poaceae species reveals the dynamic and ongoing evolutionary changes in mitochondrial genome structure and gene content. The complete mitochondrial genome of A. longiglumis completes the last link of the oat reference genome and lays the foundation for oat breeding and exploiting the biodiversity in the genus.

Apostasia Mitochondrial Genome Analysis and Monocot Mitochondria Phylogenomics

Apostasia shenzhenica belongs to the subfamily Apostasioideae and is a primitive group located at the base of the Orchidaceae phylogenetic tree. However, the A. shenzhenica mitochondrial genome (mitogenome) is still unexplored, and the phylogenetic relationships between monocots mitogenomes remain unexplored. In this study, we discussed the genetic diversity of A. shenzhenica and the phylogenetic relationships within its monocotyledon mitogenome. We sequenced and assembled the complete mitogenome of A. shenzhenica, resulting in a circular mitochondrial draft of 672,872 bp, with an average read coverage of 122× and a GC content of 44.4%. A. shenzhenica mitogenome contained 36 protein-coding genes, 16 tRNAs, two rRNAs, and two copies of nad4L. Repeat sequence analysis revealed a large number of medium and small repeats, accounting for 1.28% of the mitogenome sequence. Selection pressure analysis indicated high mitogenome conservation in related species. RNA editing identified 416 sites in the protein-coding region. Furthermore, we found 44 chloroplast genomic DNA fragments that were transferred from the chloroplast to the mitogenome of A. shenzhenica, with five plastid-derived genes remaining intact in the mitogenome. Finally, the phylogenetic analysis of the mitogenomes from A. shenzhenica and 28 other monocots showed that the evolution and classification of most monocots were well determined. These findings enrich the genetic resources of orchids and provide valuable information on the taxonomic classification and molecular evolution of monocots.

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

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

The garden asparagus (Asparagus officinalis L.) mitochondrial genome revealed rich sequence variation throughout whole sequencing data

Garden asparagus (Asparagus officinalis L.) is a horticultural crop with high nutritional and medical value, considered an ideal plant for sex determination research among many dioecious plants, whose genomic information can support genetic analysis and breeding programs. In this research, the entire mitochondrial genome of A. officinalis was sequenced, annotated and assembled using a mixed Illumina and PacBio data. The garden asparagus circular mitochondrial genome measures 492,062 bp with a GC value of 45.9%. Thirty-six protein-coding genes, 17 tRNA and 6 rRNA genes were annotated, among which 8 protein-coding genes contained 16 introns. In addition, 254 SSRs with 10 complete tandem repeats and 293 non-tandem repeats were identified. It was found that the codons of edited sites located in the amino acids showed a leucine-formation trend, and RNA editing sites mainly caused the mutual transformation of amino acids with the same properties. Furthermore, 72 sequence fragments accounting for 20,240 bp, presentating 4.11% of the whole mitochondrial genome, were observed to migrate from chloroplast to mitochondrial genome of A. officinalis. The phylogenetic analysis showed that the closest genetic relationship between A. officinalis with onion (Allium cepa) inside the Liliaceae family. Our results demonstrated that high percentage of protein-coding genes had evolutionary conservative properties, with Ka/Ks values less than 1. Therefore, this study provides a high-quality garden asparagus mitochondrial genome, useful to promote better understanding of gene exchange between organelle genomes.

Comparative analysis of mitochondrial genomes of two alpine medicinal plants of Gentiana (Gentianaceae)

Gentiana crassicaulis and G. straminea are alpine plants of Gentiana with important medicinal value and complex genetic backgrounds. In this study, the mitochondrial genomes (mtDNAs) of these two species were sequenced. The mtDNAs of G. crassicaulis and G. straminea are 368,808 and 410,086 bp long, respectively, 52 and 49 unique genes are annotated in the two species, and the gene arrangement varies widely. Compared to G. crassicaulis, G. straminea loses three effective genes, namely atp6, trnG-GCC and trnV-GAC. As a pseudogene, the atp6 gene of G. straminea is incomplete, which is rare in higher plants. We detected 1696 and 1858 pairs of long repeats and 213 SSRs and 250 SSs in the mtDNAs of G. crassicaulis and G. straminea, respectively. There are 392 SNPs and 18 InDels between the two genomes, and syntenic sequence and structural variation analysis show low collinearity between the two genomes. Chloroplast DNA transferring to mtDNA is observed in both species, and 46,511 and 55,043 bp transferred segments containing three tRNA genes are identified, respectively. Comparative analysis of mtDNAs of G. crassicaulis, G. straminea and four species of Gentianales determined 18 core genes, and there is no specific gene in G. crassicaulis and G. straminea. The phylogenetic tree based on mtDNAs places Gentianaceae in a branch of Gentianales. This study is the first to analyze the mtDNAs of Gentianaceae, which could provide information for analysis of the structure of mtDNAs of higher plants and phylogenetic research of Gentianaceae and Gentianales.

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.

A Comparative Genomics Approach for Analysis of Complete Mitogenomes of Five Actinidiaceae Plants

Actinidiaceae, an economically important plant family, includes the Actinidia, Clematoclethra and Saurauia genus. Kiwifruit, with remarkably high vitamin C content, is an endemic species widely distributed in China with high economic value. Although many Actinidiaceae chloroplast genomes have been reported, few complete mitogenomes of Actinidiaceae have been studied. Here, complete circular mitogenomes of the four kiwifruit species and Saurauia tristyla were assembled. Codon usage, sequence repeats, RNA editing, gene transfers, selective pressure, and phylogenetic relationships in the four kiwifruit species and S. tristyla were comparatively analyzed. This research will contribute to the study of phylogenetic relationships within Actiniaceae and molecular barcoding in kiwifruit.

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

BACKGROUND

Bupleurum chinense(B. chinense) is a plant that is widely distributed globally and has strong pharmacological effects. Though the chloroplast(cp) genome of B. chinense has been studied, no reports regarding the mitochondrial(mt) genome of B. chinense have been published yet.

RESULTS

The mt genome of B.chinense was assembled and functionally annotated. The circular mt genome of B. chinense was 435,023 bp in length, and 78 genes, including 39 protein-coding genes, 35 tRNA genes, and 4 rRNA genes, were annotated. Repeat sequences were analyzed and sites at which RNA editing would occur were predicted. Gene migration was observed to occur between the mt and cp genomes of B. chinense via the detection of homologous gene fragments. In addition, the sizes of plant mt genomes and their GC content were analyzed and compared. The sizes of mt genomes of plants varied greatly, but their GC content was conserved to a greater extent during evolution. Ka/Ks analysis was based on code substitutions, and the results showed that most of the coding genes were negatively selected. This indicates that mt genes were conserved during evolution.

CONCLUSION

In this study, we assembled and annotated the mt genome of the medicinal plant B. chinense. Our findings provide extensive information regarding the mt genome of B. chinense, and help lay the foundation for future studies on the genetic variations, phylogeny, and breeding of B. chinense via an analysis of the mt genome.

Comparative analysis of the chloroplast and mitochondrial genomes of Saposhnikovia divaricata revealed the possible transfer of plastome repeat regions into the mitogenome

Background

Saposhnikovia divaricata (Turcz.) Schischk. is a perennial herb whose dried roots are commonly used as a source of traditional medicines. To elucidate the organelle-genome-based phylogeny of Saposhnikovia species and the transfer of DNA between organelle genomes, we sequenced and characterised the mitochondrial genome (mitogenome) of S. divaricata.

Results

The mitogenome of S. divaricata is a circular molecule of 293,897 bp. The nucleotide composition of the mitogenome is as follows: A, 27.73%; T, 27.03%; C, 22.39%; and G, 22.85. The entire gene content is 45.24%. A total of 31 protein-coding genes, 20 tRNAs and 4 rRNAs, including one pseudogene (rpl16), were annotated in the mitogenome. Phylogenetic analysis of the organelle genomes from S. divaricata and 10 related species produced congruent phylogenetic trees. Selection pressure analysis revealed that most of the mitochondrial genes of related species are highly conserved. Moreover, 2 and 46 RNA-editing sites were found in the chloroplast genome (cpgenome) and mitogenome protein-coding regions, respectively. Finally, a comparison of the cpgenome and the mitogenome assembled from the same dataset revealed 10 mitochondrial DNA fragments with sequences similar to those in the repeat regions of the cpgenome, suggesting that the repeat regions might be transferred into the mitogenome.

Conclusions

In this study, we assembled and annotated the mitogenome of S. divaricata. This study provides valuable information on the taxonomic classification and molecular evolution of members of the family Apiaceae.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08821-0.

Complex Physical Structure of Complete Mitochondrial Genome of Quercus acutissima (Fagaceae): A Significant Energy Plant

Quercus acutissima Carruth. is a Chinese important energy plant with high ecological and economic values. While the species chloroplast genome has been reported, its mitochondrial genome (mitogenome) is still unexplored. Here, we assembled and annotated the Q. acutissima mitogenome, and we compared its characteristic differences with several closely related species. The Q. acutissima mitogenome’s main structure is branched with three distinguished contigs (linear molecule 1, circular molecule 2, and circular molecule 3) with 448,982 bp total length and 45.72% GC content. The mitogenome contained 51 genes, including 32 protein-coding, 16 tRNA and 3 rRNA genes. We examined codon usage, repeated sequences, genome recombination, chloroplast to mitochondrion DNA transformation, RNA editing, and synteny in the Q. acutissima mitogenome. Phylogenetic trees based on 29 species mitogenomes clarified the species classification. Our results provided comprehensive information of Q. acutissima mitogenome, and they are expected to provide valuable information for Fagaceae evolutionary biology and to promote the species germplasm utilization.

Deciphering the Multi-Chromosomal Mitochondrial Genome of Populus simonii

Mitochondria, inherited maternally, are energy metabolism organelles that generate most of the chemical energy needed to power cellular various biochemical reactions. Deciphering mitochondrial genome (mitogenome) is important for elucidating vital activities of species. The complete chloroplast (cp) and nuclear genome sequences of Populus simonii (P. simonii) have been reported, but there has been little progress in its mitogenome. Here, we assemble the complete P. simonii mitogenome into three circular-mapping molecules (lengths 312.5, 283, and 186 kb) with the total length of 781.5 kb. All three molecules of the P. simonii mitogenome had protein-coding capability. Whole-genome alignment analyses of four Populus species revealed the fission of poplar mitogenome in P. simonii. Comparative repeat analyses of four Populus mitogenomes showed that there were no repeats longer than 350 bp in Populus mitogenomes, contributing to the stability of genome sizes and gene contents in the genus Populus. As the first reported multi-circular mitogenome in Populus, this study of P. simonii mitogenome are imperative for better elucidating their biological functions, replication and recombination mechanisms, and their unique evolutionary trajectories in Populus.

De novo assembly of the complete mitochondrial genome of sweet potato (Ipomoea batatas [L.] Lam) revealed the existence of homologous conformations generated by the repeat-mediated recombination

Sweet potato (Ipomoea batatas [L.] Lam) is an important food crop, an excellent fodder crop, and a new type of industrial raw material crop. The lack of genomic resources could affect the process of industrialization of sweet potato. Few detailed reports have been completed on the mitochondrial genome of sweet potato. In this research, we sequenced and assembled the mitochondrial genome of sweet potato and investigated its substructure. The mitochondrial genome of sweet potato is 270,304 bp with 23 unique core genes and 12 variable genes. We detected 279 pairs of repeat sequences and found that three pairs of direct repeats could mediate the homologous recombination into four independent circular molecules. We identified 70 SSRs in the whole mitochondrial genome of sweet potato. The longest dispersed repeat in mitochondrial genome was a palindromic repeat with a length of 915 bp. The homologous fragments between the chloroplast and mitochondrial genome account for 7.35% of the mitochondrial genome. We also predicted 597 RNA editing sites and found that the rps3 gene was edited 54 times, which occurred most frequently. This study further demonstrates the existence of multiple conformations in sweet potato mitochondrial genomes and provides a theoretical basis for the evolution of higher plants and cytoplasmic male sterility breeding.

Assembly and comparative analysis of the first complete mitochondrial genome of Acer truncatum Bunge: a woody oil-tree species producing nervonic acid

BACKGROUND

Acer truncatum (purpleblow maple) is a woody tree species that produces seeds with high levels of valuable fatty acids (especially nervonic acid). The species is admired as a landscape plant with high developmental prospects and scientific research value. The A. truncatum chloroplast genome has recently been reported; however, the mitochondrial genome (mitogenome) is still unexplored.

RESULTS

We characterized the A. truncatum mitogenome, which was assembled using reads from PacBio and Illumina sequencing platforms, performed a comparative analysis against different species of Acer. The circular mitogenome of A. truncatum has a length of 791,052 bp, with a base composition of 27.11% A, 27.21% T, 22.79% G, and 22.89% C. The A. truncatum mitogenome contains 62 genes, including 35 protein-coding genes, 23 tRNA genes and 4 rRNA genes. We also examined codon usage, sequence repeats, RNA editing and selective pressure in the A. truncatum mitogenome. To determine the evolutionary and taxonomic status of A. truncatum, we conducted a phylogenetic analysis based on the mitogenomes of A. truncatum and 25 other taxa. In addition, the gene migration from chloroplast and nuclear genomes to the mitogenome were analyzed. Finally, we developed a novel NAD1 intron indel marker for distinguishing several Acer species.

CONCLUSIONS

In this study, we assembled and annotated the mitogenome of A. truncatum, a woody oil-tree species producing nervonic acid. The results of our analyses provide comprehensive information on the A. truncatum mitogenome, which would facilitate evolutionary research and molecular barcoding in Acer.

The complete chloroplast genome of Camellia fluviatilis (Theaceae), a wild oil-Camellia species

Camellia fluviatilis is an important shrub producing edible seed oil, which is widely cultivated in South China. In this study, the complete chloroplast genome was sequenced and analyzed based on the Illumina HiSeq platform. The results showed that the complete chloroplast genome is 157,041 bp with 37.29% GC content, including a large single copy (LSC) region of 86,718 bp, a small single copy (SSC) region of 18,293 bp, and a pair of inverted repeats (IRs) regions of 26,015 bp. There are 128 genes in the chloroplast genome of C. fluviatilis, including 83 protein-coding genes, 8 ribosomal RNAs, and 37 transfer RNAs. Phylogenetic analysis revealed that C. fluviatilis is closely related to C. lanceoleosa, indicating that both belong to the Sect. Paracamellia Sealy.

The complete chloroplast genome of Camellia grijsii 'zhenzhucha', a variant cultivar with floral aroma

Camellia grijsii 'zhenzhucha' is a variant cultivar from Camellia grijsii, which is also called Camellia grijsii 'juhuacha'.C. grijsii 'zhenzhucha' is an ornamental shrub with a floral aroma, and is oftenused in landscape. To provide genetic information for genetic research, we have sequenced and assembled the complete chloroplast (cp) genome of C. grijsii 'zhenzhucha' based on the Illumina Hiseq platform. The assembled complete cp genome of C. grijsii 'zhenzhucha' was 161,478 bp in length with 37.24% GC, including a large single-copy (LSC) region of 59,942 bp, a small single-copy (SSC) region of 17,294 bp, and a pair of inverted repeats (IRs) of 20,293 bp. The cp genome was annotated with 130 functional genes, consisting of 81 protein-coding genes, 45 transporter RNAs, and 4 ribosomal RNAs. To obtain the phylogeny relationship, the cp genome of C. grijsii 'zhenzhucha'has been compared with other Camellia species, and the results indicate that C. grijsii 'zhenzhucha' is closely related to C. grijsii. This study provides fundamental information of C. grijsii 'zhenzhucha' cp genome, and has an important reference value for the evolutionary analysis.

Assembly of the complete mitochondrial genome of an endemic plant, Scutellaria tsinyunensis, revealed the existence of two conformations generated by a repeat-mediated recombination

We assembled the complete mitochondrial genome of Scutellaria tsinyunensis in this study. Repeat-mediated recombination resulted in the formation of two conformations of the mitochondrial genome in S. tsinyunensis. Scutellaria tsinyunensis belongs to the family Lamiaceae, distributed only in the Jinyun Mountain, Chongqing, China. As a valuable endemic and small population species, it is regarded as a natural resource potentially with significant economic and ecological importance. In this study, we assembled a complete and gap-free mitochondrial genome of S. tsinyunensis. This genome had a length of 354,073 bp and the base composition of the genome was A (27.44%), T (27.30%), C (22.58%), and G (22.68%). This genome encodes 59 genes, including 32 protein-coding genes, 24 tRNA genes, and 3 rRNA genes. The Sanger sequencing and Oxford Nanopore sequencing confirmed a pair of direct repeats had mediated genome recombination, resulting in the formation of two conformations. The gene conversation between plastome and mitochondrial genome was also observed in S. tsinyunensis by detecting gene migration, including six tRNA genes (namely, trnW-CCA, trnI-CAU, trnH-UUU, trnD-GUC, trnN-GUU, and trnM-CAU), five protein-coding gene fragments, and the fragments from 2 rRNA genes. Moreover, the dN/dS analysis revealed the atp9 gene had undergone strong negative selection, and four genes (atp4, mttB, ccmFc, and ccmB) probably had undergone positive selection during evolution in Lamiales. This work reported the first mitochondrial genome of S. tsinyunensis, which could be used as a reference genome for the important medicinal plants of the genus Scutellaria, and also provide much-desired information for molecular breeding.

The complete chloroplast genome of Camellia grijsii, an ornamental shrub with floral aroma

Camellia grijsii is an ornamental shrub with a floral aroma, which is widely cultivated and used for landscaping in China. To obtain the genetic information of C. grijsii, we have sequenced and assembled the complete chloroplast (cp) genome based on the Illumina Hiseq platform. The total genome size is 161,078 bp in length with 37.18% GC, which contains a large single copy (LSC, 84,645 bp) region, a small single copy (SSC, 15,772 bp) region, and a pair of inverted repeat (IRs, 30,330 bp) regions. It is composed of 81 protein-coding genes, 4 ribosomal RNAs, and 43 transfer RNAs. The cp genome of C. grijsii has also been compared with other species of Camellia, and the results showed that the C. grijsii and the C. grandbibracteata are closely related. This study provides the complete cp genome of C. grijsii and has an important reference value for the evolutionary analysis.

Mitochondrial genome of Salix cardiophylla and its implications for infrageneric division of the genus of Salix

Salix cardiophylla was a member of the genus of Salix in family Salicaceae with unique morphological traits, and once recognized as a separate genus, Toisusu Kimura. Here, we sequenced and assembled the complete mitochondrial genome of S. cardiophylla, which was 735,173 bp in length, including 56 genes, 28 protein-coding genes, 3 rRNA genes, 25 tRNA genes, and one large inverted repeat regions with length of 13,603 bp. Phylogenetic analysis based on 26 mitochondrial CDS confirmed that S. cardiophylla is a member of Salix, and support its merge into Salix in aspect of our new insights on mitogenome phylogenomics.

Multicentric origin and diversification of atp6-orf79-like structures reveal mitochondrial gene flows in Oryza rufipogon and Oryza sativa

Cytoplasmic male sterility (CMS) is a widely used genetic tool in modern hybrid rice breeding. Most genes conferring rice gametophytic CMS are homologous to orf79 and co-transcribe with atp6. However, the origin, differentiation and flow of these mitochondrial genes in wild and cultivated rice species remain unclear. In this study, we performed de novo assembly of the mitochondrial genomes of 221 common wild rice (Oryza rufipogon Griff.) and 369 Asian cultivated rice (Oryza sativa L.) accessions, and identified 16 haplotypes of atp6-orf79-like structures and 11 orf79 alleles. These homologous structures were classified into 4 distinct groups (AO-I, AO-II, AO-III and AO-IV), all of which were observed in O. rufipogon but only AO-I was detected in O. sativa, causing a decrease in the frequency of atp6-orf79-like structures from 19.9% to 8.1%. Phylogenetic and biogeographic analyses revealed that the different groups of these gametophytic CMS-related genes in O. rufipogon evolved in a multicentric pattern. The geographical origin of the atp6-orf79-like structures was further traced back, and a candidate region in north-east of Gangetic Plain on the Indian Peninsula (South Asia) was identified as the origin centre of AO-I. The orf79 alleles were detected in all three cytoplasmic types (Or-CT0, Or-CT1 and Or-CT2) of O. rufipogon, but only two alleles (orf79a and orf79b) were observed in Or-CT0 type of O. sativa, while no orf79 allele was found in other types of O. sativa. Our results also revealed that the orf79 alleles in cultivated rice originated from the wild rice population in South and South-East Asia. In addition, strong positive selection pressure was detected on the sequence variations of orf79 alleles, and a special evolutionary strategy was noted in these gametophytic CMS-related genes, suggesting that their divergence could be beneficial to their survival in evolution.

The complete mitochondrial genome of Salix polaris, a specie in harsh arctic environment

Salix polaris is a dwarf species of the genus Salix which distribute in harsh arctic environment. In this study, we sequenced and assembled whole mitochondrial genome of S. polaris for the first time. The complete mitochondrial genome sequence of this species is a circular molecule of 562920 bp in size, encoding 23 CDS, 21 tRNA genes, and 3 rRNA genes. Mitogenome based phylogenetic tree showed that Salix polaris clustered in the robust clade consist of all Salix species sampled in this study. The complete mitochondrial genome of Salix polaris adds new data into the whole mitochondrial genome of Salix, which is still scarce so far, and we anticipate the accumulation of mitogenome of Salix will facilitate phylogenetic study for this systematic difficult genus.

Characterization and Analysis of the Mitochondrial Genome of Common Bean (Phaseolus vulgaris) by Comparative Genomic Approaches

The common bean (Phaseolus vulgaris) is a major source of protein and essential nutrients for humans. To explore the genetic diversity and phylogenetic relationships of P. vulgaris, its complete mitochondrial genome (mitogenome) was sequenced and assembled. The mitogenome is 395,516 bp in length, including 31 unique protein-coding genes (PCGs), 15 transfer RNA (tRNA) genes, and 3 ribosomal RNA (rRNA) genes. Among the 31 PCGs, four genes (mttB, nad1, nad4L, and rps10) use ACG as initiation codons, which are altered to standard initiation codons by RNA editing. In addition, the termination codon CGA in the ccmFC gene is converted to UGA. Selective pressure analysis indicates that the ccmB, ccmFC, rps1, rps10, and rps14 genes were under evolutionary positive selection. The proportions of five amino acids (Phe, Leu, Pro, Arg, and Ser) in the whole amino acid profile of the proteins in each mitogenome can be used to distinguish angiosperms from gymnosperms. Phylogenetic analyses show that P. vulgaris is evolutionarily closer to the Glycininae than other leguminous plants. The results of the present study not only provide an important opportunity to conduct further genomic breeding studies in the common bean, they also provide valuable information for future evolutionary and molecular studies of leguminous plants.

The complete mitochondrial genome of Taxus cuspidata (Taxaceae): eight protein-coding genes have transferred to the nuclear genome

BACKGROUND

Gymnosperms represent five of the six lineages of seed plants. However, most sequenced plant mitochondrial genomes (mitogenomes) have been generated for angiosperms, whereas mitogenomic sequences have been generated for only six gymnosperms. In particular, complete mitogenomes are available for all major seed plant lineages except Conifer II (non-Pinaceae conifers or Cupressophyta), an important lineage including six families, which impedes a comprehensive understanding of the mitogenomic diversity and evolution in gymnosperms.

RESULTS

Here, we report the complete mitogenome of Taxus cuspidata in Conifer II. In comparison with previously released gymnosperm mitogenomes, we found that the mitogenomes of Taxus and Welwitschia have lost many genes individually, whereas all genes were identified in the mitogenomes of Cycas, Ginkgo and Pinaceae. Multiple tRNA genes and introns also have been lost in some lineages of gymnosperms, similar to the pattern observed in angiosperms. In general, gene clusters could be less conserved in gymnosperms than in angiosperms. Moreover, fewer RNA editing sites were identified in the Taxus and Welwitschia mitogenomes than in other mitogenomes, which could be correlated with fewer introns and frequent gene losses in these two species.

CONCLUSIONS

We have sequenced the Taxus cuspidata mitogenome, and compared it with mitogenomes from the other four gymnosperm lineages. The results revealed the diversity in size, structure, gene and intron contents, foreign sequences, and mutation rates of gymnosperm mitogenomes, which are different from angiosperm mitogenomes.

Characterization and phylogenetic significance of the complete chloroplast genome of Camellia Kissii, an economic crop for producing oil

Camellia kissii is cultivated for a long time as an oil crop for edible and industrial oils, and has the functions of high oil production rate and unique health care. The complete chloroplast (cp) genome sequence of C. kissii is 156,961 bp in length with GC content of 39.29%. It presents a quadrate structure, including a large single-copy (LSC) region (86,640 bp), a small single-copy (SSC) region (18,399 bp), and a pair of inverted repeats (IRs) (25,961 bp). Meanwhile, 15 complete chloroplast genome of Camellia was aligned to explore the phylogenetic significance of Camellia. And the genetic relationship between Camellia kissii and Camellia huana was found to be closest.

Potential of Transcript Editing Across Mitogenomes of Early Land Plants Shows Novel and Familiar Trends

RNA editing alters the identity of nucleotides in an RNA sequence so that the mature transcript differs from the template defined in the genome. This process has been observed in chloroplasts and mitochondria of both seed and early land plants. However, the frequency of RNA editing in plant mitochondria ranges from zero to thousands of editing sites. To date, analyses of RNA editing in mitochondria of early land plants have been conducted on a small number of genes or mitochondrial genomes of a single species. This study provides an overview of the mitogenomic RNA editing potential of the main lineages of these two groups of early land plants by predicting the RNA editing sites of 33 mitochondrial genes of 37 species of liverworts and mosses. For the purpose of the research, we newly assembled seven mitochondrial genomes of liverworts. The total number of liverwort genera with known complete mitogenome sequences has doubled and, as a result, the available complete mitogenome sequences now span almost all orders of liverworts. The RNA editing site predictions revealed that C-to-U RNA editing in liverworts and mosses is group-specific. This is especially evident in the case of liverwort lineages. The average level of C-to-U RNA editing appears to be over three times higher in liverworts than in mosses, while the C-to-U editing frequency of the majority of genes seems to be consistent for each gene across bryophytes.

Decoding and analysis of organelle genomes of Indian tea (Camellia assamica) for phylogenetic confirmation

The NCBI database has >15 chloroplast (cp) genome sequences available for different Camellia species but none for C. assamica. There is no report of any mitochondrial (mt) genome in the Camellia genus or Theaceae family. With the strong believes that these organelle genomes can play a great tool for taxonomic and phylogenetic analysis, we successfully assembled and analyzed cp and mt genome of C. assamica. We assembled the complete mt genome of C. assamica in a single circular contig of 707,441 bp length comprising of a total of 66 annotated genes, including 35 protein-coding genes, 29 tRNAs and two rRNAs. The first ever cp genome of C. assamica resulted in a circular contig of 157,353 bp length with a typical quadripartite structure. Phylogenetic analysis based on these organelle genomes showed that C. assamica was closely related to C. sinensis and C. leptophylla. It also supports Caryophyllales as Superasterids.

Pacbio Sequencing Reveals Identical Organelle Genomes between American Cranberry (Vaccinium macrocarpon Ait.) and a Wild Relative

Breeding efforts in the American cranberry (Vaccinium macrocarpon Ait.), a North American perennial fruit crop of great importance, have been hampered by the limited genetic and phenotypic variability observed among cultivars and experimental materials. Most of the cultivars commercially used by cranberry growers today were derived from a few wild accessions bred in the 1950s. In different crops, wild germplasm has been used as an important genetic resource to incorporate novel traits and increase the phenotypic diversity of breeding materials. Vaccinium microcarpum (Turcz. ex Rupr.) Schmalh. and V. oxycoccos L., two closely related species, may be cross-compatible with the American cranberry, and could be useful to improve fruit quality such as phytochemical content. Furthermore, given their northern distribution, they could also help develop cold hardy cultivars. Although these species have previously been analyzed in diversity studies, genomic characterization and comparative studies are still lacking. In this study, we sequenced and assembled the organelle genomes of the cultivated American cranberry and its wild relative, V. microcarpum. PacBio sequencing technology allowed us to assemble both mitochondrial and plastid genomes at very high coverage and in a single circular scaffold. A comparative analysis revealed that the mitochondrial genome sequences were identical between both species and that the plastids presented only two synonymous single nucleotide polymorphisms (SNPs). Moreover, the Illumina resequencing of additional accessions of V. microcarpum and V. oxycoccos revealed high genetic variation in both species. Based on these results, we provided a hypothesis involving the extension and dynamics of the last glaciation period in North America, and how this could have shaped the distribution and dispersal of V. microcarpum. Finally, we provided important data regarding the polyploid origin of V. oxycoccos.

The plastid and mitochondrial genomes of Eucalyptus grandis

BACKGROUND

Land plant organellar genomes have significant impact on metabolism and adaptation, and as such, accurate assembly and annotation of plant organellar genomes is an important tool in understanding the evolutionary history and interactions between these genomes. Intracellular DNA transfer is ongoing between the nuclear and organellar genomes, and can lead to significant genomic variation between, and within, species that impacts downstream analysis of genomes and transcriptomes.

RESULTS

In order to facilitate further studies of cytonuclear interactions in Eucalyptus, we report an updated annotation of the E. grandis plastid genome, and the second sequenced and annotated mitochondrial genome of the Myrtales, that of E. grandis. The 478,813 bp mitochondrial genome shows the conserved protein coding regions and gene order rearrangements typical of land plants. There have been widespread insertions of organellar DNA into the E. grandis nuclear genome, which span 141 annotated nuclear genes. Further, we identify predicted editing sites to allow for the discrimination of RNA-sequencing reads between nuclear and organellar gene copies, finding that nuclear copies of organellar genes are not expressed in E. grandis.

CONCLUSIONS

The implications of organellar DNA transfer to the nucleus are often ignored, despite the insight they can give into the ongoing evolution of plant genomes, and the problems they can cause in many applications of genomics. Future comparisons of the transcription and regulation of organellar genes between Eucalyptus genotypes may provide insight to the cytonuclear interactions that impact economically important traits in this widely grown lignocellulosic crop species.

Assembly of the Mitochondrial Genome in the Campanulaceae Family Using Illumina Low-Coverage Sequencing

Platycodongrandiflorus (balloon flower) and Codonopsislanceolata (bonnet bellflower) are important herbs used in Asian traditional medicine, and both belong to the botanical family Campanulaceae. In this study, we designed and implemented a de novo DNA sequencing and assembly strategy to map the complete mitochondrial genomes of the first two members of the Campanulaceae using low-coverage Illumina DNA sequencing data. We produced a total of 28.9 Gb of paired-end sequencing data from the genomic DNA of P.grandiflorus (20.9 Gb) and C.lanceolata (8.0 Gb). The assembled mitochondrial genome of P.grandiflorus was found to consist of two circular chromosomes; the master circle contains 56 genes, and the minor circle contains 42 genes. The C.lanceolata mitochondrial genome consists of a single circle harboring 54 genes. Using a comparative genome structure and a pattern of repeated sequences, we show that the P.grandiflorus minor circle resulted from a recombination event involving the direct repeats of the master circle. Our dataset will be useful for comparative genomics and for evolutionary studies, and will facilitate further biological and phylogenetic characterization of species in the Campanulaceae.

The complete mitochondrial genome sequence of an alpine plant Arabis alpina

Arabis alpina (A. alpina) is an arctic-alpine flowering plant in the family Brassicaceae, naturally growing in the tundra of arctic regions and in mountains at southern latitudes. In this study, we first report the assembly of the complete A. alpina mitochondrial (mt) genome using the next-generation sequencing technologies. The A. alpina mt circular genome is 323,159 bp in length and contains 33 protein-coding genes, 18 tRNA genes and 3 rRNA genes. To analyze the phylogenic and evolutional relationship of A. alpina, a neighbour-joining phylogenetic tree was reconstructed based on the mt genome of A. alpina and other 27 plants. The complete A. alpina mt genome will be helpful in population studies or investigations of mt functions of these alpine plant species.

Towards a comprehensive picture of C-to-U RNA editing sites in angiosperm mitochondria

Our understanding of the dynamic and evolution of RNA editing in angiosperms is in part limited by the few editing sites identified to date. This study identified 10,217 editing sites from 17 diverse angiosperms. Our analyses confirmed the universality of certain features of RNA editing, and offer new evidence behind the loss of editing sites in angiosperms. RNA editing is a post-transcriptional process that substitutes cytidines (C) for uridines (U) in organellar transcripts of angiosperms. These substitutions mostly take place in mitochondrial messenger RNAs at specific positions called editing sites. By means of publicly available RNA-seq data, this study identified 10,217 editing sites in mitochondrial protein-coding genes of 17 diverse angiosperms. Even though other types of mismatches were also identified, we did not find evidence of non-canonical editing processes. The results showed an uneven distribution of editing sites among species, genes, and codon positions. The analyses revealed that editing sites were conserved across angiosperms but there were some species-specific sites. Non-synonymous editing sites were particularly highly conserved (~ 80%) across the plant species and were efficiently edited (80% editing extent). In contrast, editing sites at third codon positions were poorly conserved (~ 30%) and only partially edited (~ 40% editing extent). We found that the loss of editing sites along angiosperm evolution is mainly occurring by replacing editing sites with thymidines, instead of a degradation of the editing recognition motif around editing sites. Consecutive and highly conserved editing sites had been replaced by thymidines as result of retroprocessing, by which edited transcripts are reverse transcribed to cDNA and then integrated into the genome by homologous recombination. This phenomenon was more pronounced in eudicots, and in the gene cox1. These results suggest that retroprocessing is a widespread driving force underlying the loss of editing sites in angiosperm mitochondria.

Organellar genome assembly methods and comparative analysis of horticultural plants

Although organellar genomes (including chloroplast and mitochondrial genomes) are smaller than nuclear genomes in size and gene number, organellar genomes are very important for the investigation of plant evolution and molecular ecology mechanisms. Few studies have focused on the organellar genomes of horticultural plants. Approximately 1193 chloroplast genomes and 199 mitochondrial genomes of land plants are available in the National Center for Biotechnology Information (NCBI), of which only 39 are from horticultural plants. In this paper, we report an innovative and efficient method for high-quality horticultural organellar genome assembly from next-generation sequencing (NGS) data. Sequencing reads were first assembled by Newbler, Amos, and Minimus software with default parameters. The remaining gaps were then filled through BLASTN search and PCR. The complete DNA sequence was corrected based on Illumina sequencing data using BWA (Burrows-Wheeler Alignment tool) software. The advantage of this approach is that there is no need to isolate organellar DNA from total DNA during sample preparation. Using this procedure, the complete mitochondrial and chloroplast genomes of an ornamental plant, Salix suchowensis, and a fruit tree, Ziziphus jujuba, were identified. This study shows that horticultural plants have similar mitochondrial and chloroplast sequence organization to other seed plants. Most horticultural plants demonstrate a slight bias toward A+T rich features in the mitochondrial genome. In addition, a phylogenetic analysis of 39 horticultural plants based on 15 protein-coding genes showed that some mitochondrial genes are horizontally transferred from chloroplast DNA. Our study will provide an important reference for organellar genome assembly in other horticultural plants. Furthermore, phylogenetic analysis of the organellar genomes of horticultural plants could accurately clarify the unanticipated relationships among these plants.