The complete chloroplast genome sequence of the medicinal plant Salvia miltiorrhiza

Comparative and Phylogenetic Analysis of Complete Chloroplast Genomes of Five Mangifera Species

Background/Objectives: Mango, which is known as the "King of Tropical Fruits", is an evergreen plant belonging to the Anacardiaceae family. It belongs to the genus Mangifera, which comprises 69 species of plants found in tropical and subtropical regions, including India, Indonesia, the Malay Peninsula, Thailand, and South China. However, research on the structural information of complete chloroplast genomes of Mangifera is limited. Methods: The rapid advancement of high-throughput sequencing technology enables the acquisition of the entire chloroplast (cp) genome sequence, providing a molecular foundation for phylogenetic research. This work sequenced the chloroplast genomes of six Mangifera samples, performed a comparative analysis of the cp genomes, and investigated the evolutionary relationships within the Mangifera genus. Results: All six Mangifera samples showed a single circular molecule with a quadripartite structure, ranging from 157,604 bp to 158,889 bp in length. The number of RNA editing sites ranged from 60 to 61, with ndhB exhibiting the highest number of RNA editing sites across all species. Seven genes-namely, atpB, cemA, clpP, ndhD, petB, petD, and ycf15-exhibited a Ka/Ks value > 1, suggesting they may be under positive selection. Phylogenetic analysis revealed that Mangifera siamensis showed a close relationship between Mangifera indica and Mangifera sylvatica. Conclusions: Our comprehensive analysis of the whole cp genomes of the five Mangifera species offers significant insights regarding their phylogenetic reconstruction. Moreover, it elucidates the evolutionary processes of the cp genome within the Mangifera genus.

Decoding the chloroplast genomes of five Iranian Salvia species: insights into genomic structure, phylogenetic relationships, and molecular marker development

BACKGROUND

The genus Salvia, a prominent member of the Lamiaceae family, is renowned for its ecological, medicinal, and economic significance. Despite its importance, molecular data, particularly chloroplast (cp.) genome information, remain scarce for many native Iranian Salvia species. In this study, we sequenced and analyzed the complete cp. genomes of five Iranian Salvia species (S. aethiopis, S. sclarea, S. glutinosa, S. verticillata, and S. officinalis) to elucidate their genomic structure, evolutionary relationships, and potential for biotechnological applications.

RESULTS

The cp. genomes of the five Salvia species exhibited a conserved quadripartite structure, with sizes ranging from 151,163 to 151,662 bp, and a GC content of 38%. Each genome contained 132 or 131 genes, comprising 86 or 87 protein-coding, 8 rRNA, and 37 tRNA genes, with duplications in rpl2, rpl23, and rps12. Minor variations in gene content were observed, such as the absence of trnS-CGA in S. glutinosa. Comparative analysis of IR boundaries showed subtle expansions in S. officinalis and S. sclarea, while S. glutinosa remained stable. Trans-splicing of the rps12 gene was observed in all species, with complex structures in S. glutinosa and S. sclarea. Codon usage analysis revealed a preference for A/U-ending codons, with S. verticillata displaying unique patterns. Nucleotide diversity (Pi) identified highly variable regions, such as rpl14-rpl16 and psbK-psbI, as potential molecular markers. Phylogenetic analysis resolved distinct clades, with S. aethiopis and S. sclarea forming a close group, S. glutinosa clustering with S. chanryoenica, and S. officinalis showing genetic homogeneity with Mediterranean species. S. verticillata exhibited an earlier divergence, highlighting the genus's evolutionary complexity.

CONCLUSIONS

This study provides critical genomic resources for species identification, phylogenetic studies, and the development of molecular markers, facilitating the conservation of native Salvia species and their utilization in breeding programs for medicinal and aromatic traits.

Comparative Chloroplast Genomics of Actinidia deliciosa Cultivars: Insights into Positive Selection and Population Evolution

The chloroplast genome, as an important evolutionary marker, can provide a new breakthrough direction for the population evolution of plant species. Actinidia deliciosa represents one of the most economically significant and widely cultivated fruit species in the genus Actinidia. In this study, we sequenced and analyzed the complete chloroplast genomes of seven cultivars of Actinidia. deliciosa to detect the structural variation and population evolutionary characteristics. The total genome size ranged from 156,404 bp (A. deliciosa cv. Hayward) to 156,495 bp (A. deliciosa cv. Yate). A total of 321 simple sequence repeats (SSRs) and 1335 repetitive sequences were identified. Large-scale repeat sequences may facilitate indels and substitutions, molecular variations in A. deliciosa varieties' chloroplast genomes. Additionally, four polymorphic chloroplast DNA loci (atpF-atpH, atpH-atpI, atpB, and accD) were detected, which could potentially provide useful molecular genetic markers for further population genetics studies within A. deliciosa varieties. Site-specific selection analysis revealed that six genes (atpA, rps3, rps7, rpl22, rbcL, and ycf2) underwent protein sequence evolution. These genes may have played key roles in the adaptation of A. deliciosa to various environments. The population evolutionary analysis suggested that A. deliciosa cultivars were clustered into an individual evolutionary branch with moderate-to-high support values. These results provided a foundational genomic resource that will be a major contribution to future studies of population genetics, adaptive evolution, and genetic improvement in Actinidia.

New insights of advanced biotechnological engineering strategies for tanshinone biosynthesis in Salvia miltiorrhiza

Salvia miltiorrhiza Bunge, a well-known traditional Chinese herbal medicine, has been served as not only medicine for human ailments, but also health care products. As one of major bioactive ingredients, tanshinones are widely used to treat cardiovascular and cerebrovascular diseases, and also possess different pharmacological activities including anti-tumor, anti-inflammatory, anti-fibrotic and others. However, the content of tanshinones is relatively low in S. miltiorrhiza plants. Recently, multiple biotechnological strategies have been applied to improve tanshinone production. In this review, advances in bioactivities, biosynthesis pathway and regulation, transcriptional regulatory network, epigenetic modification and synthetic biology are summarized, and future perspectives are discussed, which will help develop high-quality S. miltiorrhiza resources.

Comprehensive Chloroplast Genomic Insights into Amaranthus: Resolving the Phylogenetic and Taxonomic Status of A. powellii and A. bouchonii

Amaranthus, a genus in Amaranthaceae, is divided into three subgenera-Amaranthus, Acnida, and Albersia-and contains approximately 70 to 80 species. Understanding its phylogenetic relationships is essential for species classification, genetic diversity assessment, and evolutionary studies. This knowledge is vital for improving Amaranthus utilization in crop improvement and managing the ecological impacts of invasive weeds. In this study, we analyzed the chloroplast genomes of 27 Amaranthus species across all three subgenera to characterize their genomic features and construct a comprehensive phylogenetic tree. Our aim was to elucidate the phylogenetic relationships within the genus and evaluate interspecific affinities among the subgenera. We also addressed the taxonomic ambiguity surrounding A. bouchonii and A. powellii to determine their distinct species within the genus. Chloroplast genome sizes ranged from 149,949 to 150,818 bp, with GC content varying between 36.52% and 36.63%. Comparative structural analyses confirmed highly conserved quadripartite structures, gene content, and organization, comprising 87 protein-coding genes, 37 tRNAs, and 8 rRNAs. Repeat and codon usage analyses revealed conserved repeat patterns and a preference for codons ending in A or U. Selection pressure analysis indicated a predominantly purifying selection, with matK showing signs of positive selection, particularly in A. spinosus. Phylogenetic analysis of 80 protein-coding genes confirmed the monophyly of subgenus Amaranthus but found Alberisa and Acnida to be paraphyletic. Despite their morphological similarity, A. bouchonii and A. powellii were placed in separate clades within subgenus Amaranthus, with A. bouchonii clustering with A. retroflexus, and A. powellii aligning with the A. hybridus complex. Additionally, we identified 16 variable regions as potential molecular markers for species identification. Our study provides the most comprehensive Amaranthus chloroplast genome dataset to date, offering new insights into its evolutionary relationships and valuable genomic resources for taxonomy, germplasm management, and invasive risk assessment.

Identification and management of a novel Danshen leaf anthracnose caused by Colletotrichum karstii in Salvia miltiorrhiza Bunge in China

Danshen (Salvia miltiorrhiza Bunge), a member of the genus Salvia within the Lamiaceae family, holds significant economic and medicinal value. Regrettably, the emergence of a novel leaf anthracnose in 2020 has significantly impacted its cultivation, leading to decreased yield and compromised quality. This newly identified pathogen was meticulously isolated from affected leaves, employing meticulous single conidia isolation techniques. Subsequent confirmation of pathogenicity was achieved through strict adherence of Koch's postulates. To ensure precise identification, morphological characteristics were supplemented with tandem sequence analysis targeting the rDNA internal transcribed spacer (ITS), β-tubulin (TUB), and histone (His3) regions. Combining molecular biology techniques with morphological observation and Koch's postulates, the pathogen was conclusively identified as Colletotrichum karstii. Further investigations focused on understanding the environmental factors influencing the mycelial growth and sporulation of the pathogen. The optimum temperature for the growth of C.karstii is 25°C, the suitable light conditions are 12h light/12h dark or 24h dark, and the suitable pH is 5 to 9. Utilizing BIOLOG phenotypic analysis technique, the metabolic utilization of carbon and nitrogen sources by the pathogen was assessed across different temperatures (20°C, 25°C, and 30°C). Results indicated the highest utilization rates at 25°C, particularly for arbutin and L-tryptophan. Lastly, the efficacy of 15 chemical fungicides and six botanical fungiticide against C. karstii was evaluated in vitro, revealing fluazinam as the most potent inhibitor against mycelial growth with EC50 of 0.0725 mg/mL for mycelium and 0.0378 mg/mL for spore germination, respectively. The 1 % osthole emulsion in water was found to have the strongest inhibitory effect on the growth of mycelium, with an EC50 value of 4.8984 µg/mL. Spore germination was most strongly inhibited by the 80 % ethylicin EC, which had an EC50 value of 0.5541 µg/mL. This study represents the first documentation of C. karstii as a causative agent of anthrax in Danshen, underscoring the significance of these findings for agricultural management and disease control strategies.

Analysis of microbial community composition and diversity in the rhizosphere of Salvia miltiorrhiza at different growth stages

Salvia miltiorrhiza is a kind of medicinal plant with various pharmacological activities. Few studies on the composition and diversity of rhizosphere microbial communities at different growth stages have been conducted on Salvia miltiorrhiz; in particular, salviorrhiza grows in soil that has been continuously planted for 3 years. The purpose of this study was to understand the changes of soil physicochemical properties of Salvia miltiorrhiza at different growth stages, and to study the composition and diversity of rhizosphere microbial community at different growth stages. Illumina NovaSeq sequencing technology was used to analyze the bacterial 16S rRNA gene and the fungal ITS region in the rhizosphere soil of Salvia miltiorrhiza at different growth stages. The results showed that the dominant bacterial phyla in the Salvia miltiorrhiza rhizosphere were Proteobacteria, Bacteroidetes, Acidobacteria, Firmicutes, Actinobacteria, and Chloroflexi. The dominant fungal phyla were Ascomycota, Mortierellomycota, Basidiomycota, and Rozellomycota. During the growth of Salvia miltiorrhiza, the physical and chemical properties of soil changed. As the Salvia miltiorrhiza grew, the content of available phosphorus, available potassium, pH, nitrate nitrogen, and ammonium nitrogen significantly decreased. Ammonium nitrogen and nitrate nitrogen had a greater impact on the bacterial community structure in the rhizosphere than on the fungal community structure. The work was to reveal differences in the rhizosphere bacterial and fungal community structure during different growth stages of Salvia miltiorrhiza, further understand the changes of rhizosphere microbial ecological characteristics and soil physicochemical properties during the cultivation of Salvia miltiorrhiza.

Deciphering the complete chloroplast genome sequence of Meconopsis torquata Prain: Insights into genome structure, comparative analysis and phylogenetic relationship

In the present study, we have characterized the complete chloroplast (Cp) genome of Meconopsis torquata Prain (family Papaveraceae), revealing the plastome size of 153,290 bp, and a GC content of 38.72 %. The cp genome features the typical circular quadripartite structure found in flowering plants, including a pair of inverted repeat regions (25,816 bp), isolated by a small single-copy region (17,740 bp) and a large single-copy (83,918 bp). Genome annotation revealed 132 genes: 87 protein-coding genes, 37 tRNAs and eight rRNAs. This comparative study demonstrated that the genome structure, gene number and GC ratio are consistent with several other cp genomes of Meconopsis and Papaver genera. A total of 120 SSRs were detected in the plastome, the majority (111) of which were mononucleotide repeats. Among the longer repeats, palindromic sequences were most common, followed by forward, reverse, and complement repeats. The whole genome alignment revealed the conserved nature of the inverted repeat region over single-copy zones. Nucleotide diversity unveiled hypervariable sites (ycf1, rps16, accD, atpB and psbD) in both the small and large single-copy regions, which could be useful for designing molecular markers for taxonomic identification. Phylogenetic analysis revealed a close alliance of M. torquata with other Meconopsis species, such as M. pinnatifolia and M. paniculata, with strong bootstrap support. Molecular dating suggests that M. torquata originated during the Tortonian age of the Miocene epoch of the Cenozoic era. These findings provide valuable insights for biological research, especially in understanding the genetic and evolutionary divergence within the Papaveraceae family.

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.

The complete chloroplast genome sequences of nine melon varieties (Cucumis melo L.): lights into comparative analysis and phylogenetic relationships

Melon (Cucumis melo L.) is one of the most extensively grown horticulture crops of the world. Based on the morphological characters, melon was formerly divided into two subspecies, Cucumis melo ssp. melo and C. melo ssp. agrestis. However, the present methods are still inadequate to distinguish between them. The phylogenetic analysis based on chloroplast genome sequences could provide essential evidence for the classification of melon varieties. We sequenced the chloroplast genomes of nine different melon varieties by the Illumina Hiseq and performed bioinformatic analyses including repeat element analysis, genome comparison and phylogenetic analysis. The results showed that the melon chloroplast genome has a typical quadripartite structure that was conserved across the analyzed sequences. Its length ranges between 155, 558 and 156, 569 bp, with a total GC content varying from 36.7% to 37%. We found 127-132 genes in melon chloroplast genomes, including 85-87 protein-coding regions, 34-37 tRNA and 6-8 rRNA genes. The molecular structure, gene order, content, codon usage, long repeats, and simple sequence repeats (SSRs) were mostly conserved among the nine sequenced genomes. Phylogenetic analysis showed that the chloroplast genome could clearly distinguish between C. melo ssp. melo and C. melo ssp. agrestis. This study not only provides valuable knowledge on melon chloroplasts, but also offers a theoretical basis and technical support for the genetic breeding of melons.

Genetic characterization and phylogenetic analysis of the Nigella sativa (black seed) plastome

In this study, the complete plastome sequence of Nigella sativa (black seed), was analyzed for the first time. The plastome spans approximately 154,120 bp, comprising four sections: the Large Single-Copy (LSC) (85,538 bp), the Small Single-Copy (SSC) (17,984 bp), and two Inverted Repeat (IR) regions (25,299 bp). A comparative study of N. sativa's plastome with ten other species from various genera in the Ranunculaceae family reveals substantial structural variations. The contraction of the inverted repeat region in N. sativa influences the boundaries of single-copy regions, resulting in a shorter plastome size than other species. When comparing the plastome of N. sativa with those of its related species, significant divergence is observed, particularly except for N. damascena. Among these, the plastome of A. glaucifolium displays the highest average pairwise sequence divergence (0.2851) with N. sativa, followed by A. raddeana (0.2290) and A. coerulea (0.1222). Furthermore, the study identified 12 distinct hotspot regions characterized by elevated Pi values (> 0.1). These regions include trnH-GUG-psbA, matK-trnQ-UUG, psbK-trnR-UCU, atpF-atpI, rpoB-psbD, ycf3-ndhJ, ndhC-cemA, petA-psaJ, trnN-GUU-ndhF, trnV-GAC-rps12, ycf2-trnI-CAU, and ndhA-ycf1. Approximately, 24 tandem and 48 palindromic and forward repeats were detected in N. sativa plastome. The analysis revealed 32 microsatellites with the majority being mononucleotide repeats. In the N. sativa plastome, phenylalanine had the highest number of codons (1982 codons), while alanine was the least common amino acid with 260 codons. A phylogenetic tree, constructed using protein-coding genes, revealed a distinct monophyletic clade comprising N. sativa and N. damascene, closely aligned with the Cimicifugeae tribe and exhibiting robust support. This plastome provides valuable genetic information for precise species identification, phylogenetic resolution, and evolutionary studies of N. sativa.

The chloroplast genomes of two medicinal species (Veronica anagallis-aquatica L. and Veronica undulata Wall.) and its comparative analysis with related Veronica species

Veronica anagallis-aquatica L. and Veronica undulata Wall. are widely used ethnomedicinal plants in China. The two species have different clinical efficacies, while their extremely similar morphology and unclear interspecific relationship make it difficult to accurately identify them, leading to increased instances of mixed usage. This article reports on the complete chloroplast genomes sequence of these two species and their related Veronica species to conduct a comparative genomics analysis and phylogenetic construction. The results showed that the chloroplast (cp) genomes of Veronica exhibited typical circular quadripartite structures, with total lengths of 149,386 to 152,319 base pairs (bp), and GC content of 37.9 to 38.1%, and the number of genes was between 129-134. The total number of simple sequence repeats (SSRs) in V. anagallis-aquatica and V. undulata is 37 and 36, while V. arvensis had the highest total number of 56, predominantly characterized by A/T single bases. The vast majority of long repeat sequence types are forward repeats and palindromic repeats. Selective Ka/Ks values showed that three genes were under positive selection. Sequence differences often occur in the non-coding regions of the large single-copy region (LSC) and small single-copy region (SSC), with the lowest sequence variation in the inverted repeat regions (IR). Seven highly variable regions (trnT-GGU-psbD, rps8-rpl16, trnQ-UUG, trnN-GUU-ndhF, petL, ycf3, and ycf1) were detected, which may be potential molecular markers for identifying V. anagallis-aquatica and V. undulata. The phylogenetic tree indicates that there is a close genetic relationship between the genera Veronica and Neopicrorhiza, and V. anagallis-aquatica and V. undulata are sister groups. The molecular clock analysis results indicate that the divergence time of Veronica may occur at ∼ 9.09 Ma, and the divergence time of these two species occurs at ∼ 0.48 Ma. It is speculated that climate change may be the cause of Veronica species diversity and promote the radiation of the genus. The chloroplast genome data of nine Veronica specie provides important insights into the characteristics and evolution of the chloroplast genome of this genus, as well as the phylogenetic relationships of the genus Veronica.

Plastid genomes provide insights into the phylogeny and chloroplast evolution of the paper daisy tribe Gnaphalieae (Asteraceae)

Chloroplast genomes, as an essential source of phylogenetic information, are increasingly utilized in the evolutionary study of angiosperms. Gnaphalieae is a medium-sized tribe of the sunflower family of Asteraceae, with about 2,100 species in 178 genera distributed in temperate habitats worldwide. There has been considerable progress in our understanding of their phylogenetic evolution using both nuclear and chloroplast sequences, but no focus on chloroplast genomic data. In this study, we performed sequencing, assembly, and annotation of 16 representative chloroplast genomes from all the major lineages of Gnaphalieae. Our results showed that the plastomes exhibited a typical circular tetrad structure with similar genomic structure gene content. But there were differences in genome size, SSRs, and codon usage within the tribe. Phylogenetic analysis revealed Relhania clade is the earliest diverged lineages with the Lasiopogon clade and the Gnaphalium s.s. clade diverged subsequently. The core group includes FLAG clade sister to the HAP and Australasian group. Compared with the outgroup species, chloroplast genome size of the FLAG clade is much reduced whereas those of Australasian, HAP, Gnaphalium s.s., Lasiopogon and Relhania clades are relatively expanded. Insertions and deletions in the intergenic regions associated with repetitive sequence variations are supposed to be the main factor leading to length variations in the chloroplast genomes of Gnaphalieae. The comparative analyses of chloroplast genomes would provide useful implications into understanding the taxonomic and evolutionary history of Gnaphalieae.

The complete chloroplast genome of Lagochilus ilicifolius Bunge ex Bentham, Labiat. Gen. 1834 (Lamiaceae) and its phylogenetic analysis

Lagochilus ilicifolius Bunge ex Bentham, Labiat. Gen is a perennial herb with much-branched stems native to Nei Mongol, Ningxia, Gansu, N Shaanxi. It can be used clinically as a hemostatic agent. The chloroplast genome length is 151,466 bp. It contained two inverted repeat regions of 25,660 bp each, a large single-copy region of length 82,504 bp, and a small single-copy region of length 17,642 bp. Also, the GC content is 38.6%. There were 133 genes annotated, including 88 known protein-coding genes, 37 tRNAs, and eight rRNAs. The phylogenetic tree was constructed using Bayesian method for plastome data of 29 species. The entire chloroplast genome of L. ilicifolius within the Lamiaceae is the first to reveal genetic taxonomy at the molecular level, and the new phylogenetic tree data can be used for future evolutionary studies.

Comparative genomics and phylogenomics of the genus Glycyrrhiza (Fabaceae) based on chloroplast genomes

Glycyrrhiza (Fabaceae) species are rich in metabolites and widely used in medicine. Research on the chloroplast genome of Glycyrrhiza is important for understanding its phylogenetics, biogeography, genetic diversity, species identification, and medicinal properties. In this study, comparative genomics and phylogenomics of Glycyrrhiza were analyzed based on the chloroplast genome. The chloroplast genomes of six Glycyrrhiza species were obtained using various assembly and annotation tools. The final assembled chloroplast genome sizes for the six Glycyrrhiza species ranged from 126,380 bp to 129,115 bp, with a total of 109-110 genes annotated. Comparative genomics results showed that the chloroplast genomes of Glycyrrhiza showed typically lacking inverted repeat regions, and the genome length, structure, GC content, codon usage, and gene distribution were highly similar. Bioinformatics analysis revealed the presence of 69-96 simple sequence repeats and 61-138 long repeats in the chloroplast genomes. Combining the results of mVISTA and nucleotide diversity, four highly variable regions were screened for species identification and relationship studies. Selection pressure analysis indicated overall purifying selection in the chloroplast genomes of Glycyrrhiza, with a few positively selected genes potentially linked to environmental adaptation. Phylogenetic analyses involving all tribes of Fabaceae with published chloroplast genomes elucidated the evolutionary relationships, and divergence time estimation estimated the chronological order of species differentiations within the Fabaceae family. The results of phylogenetic analysis indicated that species from the six subfamilies formed distinct clusters, consistent with the classification scheme of the six subfamilies. In addition, the inverted repeat-lacking clade in the subfamily Papilionoideae clustered together, and it was the last to differentiate. Co-linear analysis confirmed the conserved nature of Glycyrrhiza chloroplast genomes, and instances of gene rearrangements and inversions were observed in the subfamily Papilionoideae.

The complete plastome sequences of invasive weed Parthenium hysterophorus: genome organization, evolutionary significance, structural features, and comparative analysis

Parthenium hysterophorus, a globally widespread weed, poses a significant threat to agricultural ecosystems due to its invasive nature. We investigated the chloroplast genome of P. hysterophorus in this study. Our analysis revealed that the chloroplast genome of P. hysterophorus spans a length of 151,881 base pairs (bp). It exhibits typical quadripartite structure commonly found in chloroplast genomes, including inverted repeat regions (IR) of 25,085 bp, a small single copy (SSC) region of 18,052 bp, and a large single copy (LSC) region of 83,588 bp. A total of 129 unique genes were identified in P. hysterophorus chloroplast genomes, including 85 protein-coding genes, 36 tRNAs, and eight rRNAs genes. Comparative analysis of the P. hysterophorus plastome with those of related species from the tribe Heliantheae revealed both conserved structures and intriguing variations. While many structural elements were shared among the species, we identified a rearrangement in the large single-copy region of P. hysterophorus. Moreover, our study highlighted notable gene divergence in several specific genes, namely matK, ndhF, clpP, rps16, ndhA, rps3, and ndhD. Phylogenetic analysis based on the 72 shared genes placed P. hysterophorus in a distinct clade alongside another species, P. argentatum. Additionally, the estimated divergence time between the Parthenium genus and Helianthus (sunflowers) was approximately 15.1 million years ago (Mya). These findings provide valuable insights into the evolutionary history and genetic relationships of P. hysterophorus, shedding light on its divergence and adaptation over time.

Chloroplast genome of Justicia procumbens: genomic features, comparative analysis, and phylogenetic relationships among Justicieae species

Justicia procumbens L. is a traditional medicinal plant that is widely distributed in China. However, little is known about the genetic diversity and evolution of this genus, and no genomic studies have been carried out on J. procumbens previously. In this study, we aimed to assemble and annotate the first complete chloroplast genome (cpDNA) of J. procumbens and compare it with all previously published cpDNAs within the tribe Justicieae. Genome structure and comparative and phylogenetic analyses were performed. The 150,454 bp-long J. procumbens cpDNA has a circular and quadripartite structure consisting of a large single copy, a small single copy, and two inverted repeat regions. It contains 133 genes, of which 88 are protein-coding genes, 37 are tRNA genes, and eight are rRNA genes. Twenty-four simple sequence repeats (SSRs) and 81 repeat sequences were identified. Comparative analyses with other Justicieae species revealed that the non-coding regions of J. procumbens cpDNA showed greater variation than did the coding regions. Moreover, phylogenetic analysis based on 14 cpDNA sequences from Justicieae species showed that J. procumbens and J. flava were most closely related. This study provides valuable genetic information to support further research on the genetic diversity and evolutionary development of the tribe Justicieae.

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

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

Comparative chloroplast genomics of 24 species shed light on the genome evolution and phylogeny of subtribe Coelogyninae (Orchidaceae)

BACKGROUND

The orchids of the subtribe Coelogyninae are among the most morphologically diverse and economically important groups within the subfamily Epidendroideae. Previous molecular studies have revealed that Coelogyninae is an unambiguously monophyletic group. However, intergeneric and infrageneric relationships within Coelogyninae are largely unresolved. There has been long controversy over the classification among the genera within the subtribe.

RESULTS

The complete chloroplast (cp.) genomes of 15 species in the subtribe Coelogyninae were newly sequenced and assembled. Together with nine available cp. genomes in GenBank from representative clades of the subtribe, we compared and elucidated the characteristics of 24 Coelogyninae cp. genomes. The results showed that all cp. genomes shared highly conserved structure and contained 135 genes arranged in the same order, including 89 protein-coding genes, 38 tRNAs, and eight rRNAs. Nevertheless, structural variations in relation to particular genes at the IR/SC boundary regions were identified. The diversification pattern of the cp. genomes showed high consistency with the phylogenetic placement of Coelogyninae. The number of different types of SSRs and long repeats exhibited significant differences in the 24 Coelogyninae cp. genomes, wherein mononucleotide repeats (A/T), and palindromic repeats were the most abundant. Four mutation hotspot regions (ycf1a, ndhF-rp132, psaC-ndhE, and rp132-trnL) were determined, which could serve as effective molecular markers. Selection pressure analysis revealed that three genes (ycf1a, rpoC2 and ycf2 genes) might have experienced apparent positive selection during the evolution. Using the alignments of whole cp. genomes and protein-coding sequences, this study presents a well-resolved phylogenetic framework of Coelogyninae.

CONCLUSION

The inclusion of 55 plastid genome data from a nearly complete generic-level sampling provide a comprehensive view of the phylogenetic relationships among genera and species in subtribe Coelogyninae and illustrate the diverse genetic variation patterns of plastid genomes in this species-rich plant group. The inferred relationships and informally recognized major clades within the subtribe are presented. The genetic markers identified here will facilitate future studies on the genetics and phylogeny of subtribe Coelogyninae.

An insight on the complete chloroplast genome of Gomphocarpus siniacus and Duvalia velutina, Asclepiadoideae (Apocynaceae)

Abstract We studied the complete chloroplast genome of Gomphocarpus siniacus and Duvalia velutina from Asclepiadoideae subfamily; due to their medicinal importance and distribution worldwide their interest became high. In this study we analyzed the complete chloroplast genomes of G. siniacus and D. velutina using Illumina sequencing technology. The sequences were compared with the other species from Apocynaceae family. The complete genome of G. siniacus is 162,570 bp while D. velutina has154, 478 bp in length. Both genomes consist of 119 genes; encode 31 tRNA genes, and eight rRNA genes. Comparative studies of the two genomes showed variations in SSR markers in which G. siniacus possesses 223 while D. velutina has 186. This could be used for barcoding in order to aid in easy identification of the species. Phylogenetic analysis on the other hand reaffirms the tribal position of G. siniacus in Asclepiadeae and D. velutina in Ceropegieae. These findings could be used in subsequent research studies of angiosperms identification, genetic engineering, herb genomics and phylogenomic studies of Apocynaceae family.

Phylogeny and Evolution of Cocconeiopsis (Cocconeidaceae) as Revealed by Complete Chloroplast and Mitochondrial Genomes

The genus Cocconeiopsis was separated from Navicula, but its systematic position is in debate. We sequenced the complete chloroplast and mitochondrial genome of Cocconeidaceae for the first time with Cocconeiopsis kantsiensis and investigated its phylogeny and evolutionary history. Results showed that the plastid genome was 140,415 bp long with 167 genes. The mitochondrial genome was 43,732 bp long with 66 genes. Comparative analysis showed that the plastid genome structure of C. kantsiensis was most similar to those of three Navicula species and Halamphora americana, and its size was significantly smaller than that of a monoraphid species. Its mitochondrial genome was similar to that of related species except for Phaeodactylum tricornutum. The multigene phylogeny reconstruction showed that Cocconeiopsis was sister to Didymosphenia but distant from Naviculaceae. The two-gene phylogenetic analysis containing 255 species showed Cocconeiopsis was sister to Cocconeis, and distant from Naviculaceae as well. Divergence time estimation indicates the common ancestor of cocconeid species occurred about 62.8 Ma and Cocconeiopsis diverged with monoraphid Cocconeis about 58.9 Ma. Our results support the assignment of Cocconeiopsis to Cocconeidaceae and that monoraphid cocconeids were likely evolved from the lineage of Cocconeiopsis.

Complete Chloroplast Genomes of Pterodon emarginatus Vogel and Pterodon pubescens Benth: Comparative and Phylogenetic Analyses

Background

The species Pterodon emarginatus and P. pubescens, popularly known as white sucupira or faveira, are native to the Cerrado biome and have the potential for medicinal use and reforestation. They are sister species with evolutionary proximity.

Objective

Considering that the chloroplast genome exhibits a conserved structure and genes, the analysis of its sequences can contribute to the understanding of evolutionary, phylogenetic, and diversity issues.

Methods

The chloroplast genomes of P. emarginatus and P. pubescens were sequenced on the Illumina MiSeq platform. The genomes were assembled based on the de novo strategy. We performed the annotation of the genes and the repetitive regions of the genomes. The nucleotide diversity and phylogenetic relationships were analyzed using the gene sequences of these species and others of the Leguminosae family, whose genomes are available in databases.

Results

The complete chloroplast genome of P. emarginatus is 159,877 bp, and that of P. pubescens is 159,873 bp. The genomes of both species have circular and quadripartite structures. A total of 127 genes were predicted in both species, including 110 single-copy genes and 17 duplicated genes in the inverted regions. 141 microsatellite regions were identified in P. emarginatus and 140 in P. pubescens. The nucleotide diversity estimates of the gene regions in twenty-one species of the Leguminosae family were 0.062 in LSC, 0.086 in SSC, and 0.036 in IR. The phylogenetic analysis demonstrated the proximity between the genera Pterodon and Dipteryx, both from the clade Dipterygeae. Ten pairs of primers with potential for the development of molecular markers were designed.

Conclusion

The genetic information obtained on the chloroplast genomes of P. emarginatus and P. pubescens presented here reinforces the similarity and evolutionary proximity between these species, with a similarity percentage of 99.8%.

Comparative analysis of chloroplast genome and new insights into phylogenetic relationships of Ajuga and common adulterants

Introduction

The potential contamination of herbal medicinal products poses a significant concern for consumer health. Given the limited availability of genetic information concerning Ajuga species, it becomes imperative to incorporate supplementary molecular markers to enhance and ensure accurate species identification.

Methods

In this study, the chloroplast (cp) genomes of seven species of the genus Ajuag were sequenced, de novo assembled and characterized.

Results

exhibiting lengths ranging from 150,342 bp to 150,472 bp, encompassing 86 - 88 protein-coding genes (PCGs), 35 - 37 transfer RNA, and eight ribosomal RNA. The repetitive sequences, codon uses, and cp genomes of seven species were highly conserved, and PCGs were the reliable molecular markers for investigating the phylogenetic relationship within the Ajuga genus. Moreover, four mutation hotspot regions (accD-psaI, atpH-atpI, ndhC-trnV(UAC), and ndhF-rpl23) were identified within cp genomes of Ajuga, which could help distinguish A. bracteosa and its contaminants. Based on cp genomes and PCGs, the phylogenetic tree preliminary confirmed the position of Ajuga within the Lamiaceae family. It strongly supported a sister relationship between Subsect. Genevense and Subsect. Biflorae, suggesting the merger of Subsect. Biflorae and Subsect. Genevenses into one group rather than maintaining separate categorizations. Additionally, molecular clock analysis estimated the divergence time of Ajuga to be around 7.78 million years ago.

Discussion

The species authentication, phylogeny, and evolution analyses of the Ajuga species may benefit from the above findings.

Comparative chloroplast genomes provided insights into the evolution and species identification on the Datureae plants

Generally, chloroplast genomes of angiosperms are always highly conserved but carry a certain number of variation among species. In this study, chloroplast genomes of 13 species from Datureae tribe that are of importance both in ornamental gardening and medicinal usage were studied. In addition, seven chloroplast genomes from Datureae together with two from Solanaceae species retrieved from the National Center for Biotechnology Information (NCBI) were integrated into this study. The chloroplast genomes ranged in size from 154,686 to 155,979 and from 155,497 to 155,919 bp for species of Datura and Brugmansia, respectively. As to Datura and Brugmansia, a total of 128 and 132 genes were identified, in which 83 and 87 protein coding genes were identified, respectively; Furthermore, 37 tRNA genes and 8 rRNA genes were both identified in Datura and Brugmansia. Repeats analysis indicated that the number and type varied among species for Simple sequence repeat (SSR), long repeats, and tandem repeats ranged in number from 53 to 59, 98 to 99, and 22 to 30, respectively. Phylogenetic analysis based on the plastid genomes supported the monophyletic relationship among Datura and Brugmansia and Trompettia, and a refined phylogenic relationships among each individual was resolved. In addition, a species-specific marker was designed based on variation spot that resulted from a comparative analysis of chloroplast genomes and verified as effective maker for identification of D. stramonium and D. stramonium var. inermis. Interestingly, we found that 31 genes were likely to be under positive selection, including genes encoding ATP protein subunits, photosystem protein subunit, ribosome protein subunits, NAD(P)H dehydrogenase complex subunits, and clpP, petB, rbcL, rpoCl, ycf4, and cemA genes. These genes may function as key roles in the adaption to diverse environment during evolution. The diversification of Datureae members was dated back to the late Oligocene periods. These chloroplast genomes are useful genetic resources for taxonomy, phylogeny, and evolution for Datureae.

Comparative Analysis of Plastomes in Elsholtzieae: Phylogenetic Relationships and Potential Molecular Markers

The Elsholtzieae, comprising ca. 7 genera and 70 species, is a small tribe of Lamiaceae (mint family). Members of Elsholtzieae are of high medicinal, aromatic, culinary, and ornamentals value. Despite the rich diversity and value of Elsholtzieae, few molecular markers or plastomes are available for phylogenetics. In the present study, we employed high-throughput sequencing to assemble two Mosla plastomes, M. dianthera and M. scabra, for the first time, and compared with other plastomes of Elsholtzieae. The plastomes of Elsholtzieae exhibited a quadripartite structure, ranging in size from 148,288 bp to 152,602 bp. Excepting the absence of the pseudogene rps19 in Elsholtzia densa, the exhaustive tally revealed the presence of 132 genes (113 unique genes). Among these, 85 protein-coding genes (CDS), 37 tRNA genes, 8 rRNA genes, and 2 pseudogenes (rps19 and ycf1) were annotated. Comparative analyses showed that the plastomes of these species have minor variations at the gene level. Notably, the E. eriostchya plastid genome exhibited increased GC content regions in the LSC and SSC, resulting in an increased overall GC content of the entire plastid genome. The E. densa plastid genome displayed modified boundaries due to inverted repeat (IR) contraction. The sequences of CDS and intergenic regions (IGS) with elevated variability were identified as potential molecular markers for taxonomic inquiries within Elsholtzieae. Phylogenetic analysis indicated that four genera formed monophyletic entities, with Mosla and Perilla forming a sister clade. This clade was, in turn, sister to Collinsonia, collectively forming a sister group to Elsholtzia. Both CDS, and CDS + IGS could construct a phylogenetic tree with stronger support. These findings facilitate species identification and DNA barcoding investigations in Elsholtzieae and provide a foundation for further exploration and resource utilization within this tribe.

Comparative and phylogenetic analysis of complete chloroplast genome sequences of Salvia regarding its worldwide distribution

Salvia is widely used as medicine, food, and ornamental plants all over the world, with three main distribution centers, the Central and western Asia/Mediterranean (CAM), the East Aisa (EA), and the Central and South America (CASA). Along with its large number of species and world-wide distribution, Salvia is paraphyletic with multiple diversity. Chloroplast genomes (CPs) are useful tools for analyzing the phylogeny of plants at lower taxonomic levels. In this study, we reported chloroplast genomes of five species of Salvia and performed phylogenetic analysis with current available CPs of Salvia. Repeated sequence analysis and comparative analysis of Salvia CPs were also performed with representative species from different distribution centers. The results showed that the genetic characters of the CPs are related to the geographic distribution of plants. Species from CAM diverged first to form a separate group, followed by species from EA, and finally species from CASA. Larger variations of CPs were observed in species from CAM, whereas more deficient sequences and less repeated sequences in the CPs were observed in species from CASA. These results provide valuable information on the development and utilization of the worldwide genetic resources of Salvia.

Comparative chloroplast genomics and insights into the molecular evolution of Tanaecium (Bignonieae, Bignoniaceae)

Species of Tanaecium (Bignonieae, Bignoniaceae) are lianas distributed in the Neotropics and centered in the Amazon. Members of the genus exhibit exceptionally diverse flower morphology and pollination systems. Here, we sequenced, assembled, and annotated 12 complete and four partial chloroplast genomes representing 15 Tanaecium species and more than 70% of the known diversity in the genus. Gene content and order were similar in all species of Tanaecium studied, with genome sizes ranging between 158,470 and 160,935 bp. Tanaecium chloroplast genomes have 137 genes, including 80-81 protein-coding genes, 37 tRNA genes, and four rRNA genes. No rearrangements were found in Tanaecium plastomes, but two different patterns of boundaries between regions were recovered. Tanaecium plastomes show nucleotide variability, although only rpoA was hypervariable. Multiple SSRs and repeat regions were detected, and eight genes were found to have signatures of positive selection. Phylogeny reconstruction using 15 Tanaecium plastomes resulted in a strongly supported topology, elucidating several relationships not recovered previously and bringing new insights into the evolution of the genus.

Comparative chloroplast genome analysis of Sambucus L. (Viburnaceae): inference for phylogenetic relationships among the closely related Sambucus adnata Wall. ex DC Sambucus javanica Blume

Sambucus L. is found in the family Viburnaceae (syn. Adoxaceae) and encompasses approximately 29 accepted species. The complex morphology of these species has caused continued confusion concerning their nomenclature, classification, and identification. Despite previous attempts to resolve taxonomic complexities in the Sambucus genus, there are still unclear phylogenetic relationships among several species. In this study, the newly obtained plastome of Sambucus williamsii Hance. as well as the populations of Sambucus canadensis L., Sambucus javanica Blume, and Sambucus adnata Wall. ex DC were sequenced, and their sizes, structural similarity, gene order, gene number, and guanine-cytosine (GC) contents were analyzed. The phylogenetic analyses were conducted using the whole chloroplast genomes and protein-coding genes (PCGs). The findings revealed that the chloroplast genomes of Sambucus species exhibited typical quadripartite double-stranded DNA molecules. Their lengths ranged from 158,012 base pairs (bp) (S. javanica) to 158,716 bp (S. canadensis L). Each genome comprised a pair of inverted repeats (IRs), which separated the large single-copy (LSC) and small single-copy (SSC) regions. In addition, the plastomes contained 132 genes, encompassing 87 protein-coding, 37 tRNA, and four rRNA genes. In the simple sequence repeat (SSR) analysis, A/T mononucleotides had the highest proportion, with the most repetitive sequences observed in S. williamsii. The comparative genome analyses showed high similarities in structure, order, and gene contents. The hypervariable regions in the studied chloroplast genomes were trnT-GGU, trnF-GAA, psaJ, trnL-UAG, ndhF, and ndhE, which may be used as candidate barcodes for species discrimination in Sambucus genus. Phylogenetic analyses supported the monophyly of Sambucus and revealed the separation of S. javanica and S. adnata populations. Sambucus chinensis Lindl. was nested within S. javanica in the same clade, collaborating their conspecific treatment. These outcomes indicate that the chloroplast genome of Sambucus plants is a valuable genetic resource for resolving taxonomic discrepancies at the lower taxonomic levels and can be applied in molecular evolutionary studies.

Comparative chloroplast genome analysis of Sambucus L. (Viburnaceae): inference for phylogenetic relationships among the closely related Sambucus adnata Wall. ex DC Sambucus javanica Blume

Sambucus L. is found in the family Viburnaceae (syn. Adoxaceae) and encompasses approximately 29 accepted species. The complex morphology of these species has caused continued confusion concerning their nomenclature, classification, and identification. Despite previous attempts to resolve taxonomic complexities in the Sambucus genus, there are still unclear phylogenetic relationships among several species. In this study, the newly obtained plastome of Sambucus williamsii Hance. as well as the populations of Sambucus canadensis L., Sambucus javanica Blume, and Sambucus adnata Wall. ex DC were sequenced, and their sizes, structural similarity, gene order, gene number, and guanine–cytosine (GC) contents were analyzed. The phylogenetic analyses were conducted using the whole chloroplast genomes and protein-coding genes (PCGs). The findings revealed that the chloroplast genomes of Sambucus species exhibited typical quadripartite double-stranded DNA molecules. Their lengths ranged from 158,012 base pairs (bp) (S. javanica) to 158,716 bp (S. canadensis L). Each genome comprised a pair of inverted repeats (IRs), which separated the large single-copy (LSC) and small single-copy (SSC) regions. In addition, the plastomes contained 132 genes, encompassing 87 protein-coding, 37 tRNA, and four rRNA genes. In the simple sequence repeat (SSR) analysis, A/T mononucleotides had the highest proportion, with the most repetitive sequences observed in S. williamsii. The comparative genome analyses showed high similarities in structure, order, and gene contents. The hypervariable regions in the studied chloroplast genomes were trnT-GGU, trnF-GAA, psaJ, trnL-UAG, ndhF, and ndhE, which may be used as candidate barcodes for species discrimination in Sambucus genus. Phylogenetic analyses supported the monophyly of Sambucus and revealed the separation of S. javanica and S. adnata populations. Sambucus chinensis Lindl. was nested within S. javanica in the same clade, collaborating their conspecific treatment. These outcomes indicate that the chloroplast genome of Sambucus plants is a valuable genetic resource for resolving taxonomic discrepancies at the lower taxonomic levels and can be applied in molecular evolutionary studies.

Complete chloroplast genomes of five Cuscuta species and their evolutionary significance in the Cuscuta genus

BACKGROUND

Cuscuta, a parasitic plant species in the Convolvulaceae family, grows in many countries and regions. However, the relationship between some species is still unclear. Therefore, more studies are needed to assess the variation of the chloroplast (cp) genome in Cuscuta species and their relationship with subgenera or sections, thus, providing important information on the evolution of Cuscuta species.

RESULTS

In the present study, we identified the whole cp genomes of C. epithymum, C. europaea, C. gronovii, C. chinensis and C. japonica, and then constructed a phylogenetic tree of 23 Cuscuta species based on the complete genome sequences and protein-coding genes. The complete cp genome sequences of C. epithymum and C. europaea were 96,292 and 97,661 bp long, respectively, and lacked an inverted repeat region. Most cp genomes of Cuscuta spp. have tetragonal and circular structures except for C. epithymum, C. europaea, C. pedicellata and C. approximata. Based on the number of genes and the structure of cp genome and the patterns of gene reduction, we found that C. epithymum and C. europaea belonged to subgenus Cuscuta. Most of the cp genomes of the 23 Cuscuta species had single nucleotide repeats of A and T. The inverted repeat region boundaries among species were similar in the same subgenera. Several cp genes were lost. In addition, the numbers and types of the lost genes in the same subgenus were similar. Most of the lost genes were related to photosynthesis (ndh, rpo, psa, psb, pet, and rbcL), which could have gradually caused the plants to lose the ability to photosynthesize.

CONCLUSION

Our results enrich the data on cp. genomes of genus Cuscuta. This study provides new insights into understanding the phylogenetic relationships and variations in the cp genome of Cuscuta species.

Molecular structure, comparative and phylogenetic analysis of the complete chloroplast genome sequences of weedy rye Secale cereale ssp. segetale

The complete chloroplast genome of Secale cereale ssp. segetale (Zhuk.) Roshev. (Poaceae: Triticeae) was sequenced and analyzed to better use its genetic resources to enrich rye and wheat breeding. The study was carried out using the following methods: DNA extraction, sequencing, assembly and annotation, comparison with other complete chloroplast genomes of the five Secale species, and multigene phylogeny. As a result of the study, it was determined that the chloroplast genome is 137,042 base pair (bp) long and contains 137 genes, including 113 unique genes and 24 genes which are duplicated in the IRs. Moreover, a total of 29 SSRs were detected in the Secale cereale ssp. segetale chloroplast genome. The phylogenetic analysis showed that Secale cereale ssp. segetale appeared to share the highest degree of similarity with S. cereale and S. strictum. Intraspecific diversity has been observed between the published chloroplast genome sequences of S. cereale ssp. segetale. The genome can be accessed on GenBank with the accession number (OL688773).

Genetic relationship of Pleione based on the chloroplast genome

Pleione (Orchidaceae) is a very famous horticultural plant with a high international reputation for its unique flower shape and abundant color. The small difference in morphological characteristics among Pleione species caused by weak reproductive isolation and easy hybridization makes the taxonomic status of individual species very confusing. Chloroplast (cp) genome information is of great significance for the study of plant phylogeny and taxonomy. In this study, the cp genomes of Pleione were sequenced and the complete cp structure and sequence characteristics of 19 species were compared and analyzed. The cp genome of Pleione species exhibited a conserved tetrad structure and each species encoded 135 protein-coding genes, 38 tRNA and 8 RNA genes. The cp genome sizes of 19 Pleione were 157964-159269 bp and the length of LSC, SSC and IR were 85953-87172 bp, 18499-18712 bp, 26459-26756 bp, respectively. Palindromic and forward repeats accounted for a high proportion and the SSRs were mainly mononucleotide repeats in Pleione. Analysis of chloroplast sequence differences indicated that the differences in coding regions were smaller than those in non-coding regions, and the variation in LSC and SSC regions was greater than that in IR regions. Phylogenetic analysis showed that all Pleione species inferred from the cp genome were clustered together and received high support. However, the genetic relationship of Pleione plants is different from the current update system of this genus. Therefore, the demarcation of Pleione interspecific relationships still needs further investigation due to the lack of sufficient evidence. The cp genome serves as valuable information for the identification of Pleione species and the study of phylogenetic relationships.

A near complete genome assembly of chia assists in identification of key fatty acid desaturases in developing seeds

Chia is an annual crop whose seeds have the highest content of α-linolenic acid (ALA) of any plant known to date. We generated a high-quality assembly of the chia genome using circular consensus sequencing (CCS) of PacBio. The assembled six chromosomes are composed of 21 contigs and have a total length of 361.7 Mb. Genome annotation revealed a 53.5% repeat content and 35,850 protein-coding genes. Chia shared a common ancestor with Salvia splendens ~6.1 million years ago. Utilizing the reference genome and two transcriptome datasets, we identified candidate fatty acid desaturases responsible for ALA biosynthesis during chia seed development. Because the seed of S. splendens contains significantly lower proportion of ALA but similar total contents of unsaturated fatty acids, we suggest that strong expression of two ShFAD3 genes are critical for the high ALA content of chia seeds. This genome assembly will serve as a valuable resource for breeding, comparative genomics, and functional genomics studies of chia.

Application of third-generation sequencing to herbal genomics

There is a long history of traditional medicine use. However, little genetic information is available for the plants used in traditional medicine, which limits the exploitation of these natural resources. Third-generation sequencing (TGS) techniques have made it possible to gather invaluable genetic information and develop herbal genomics. In this review, we introduce two main TGS techniques, PacBio SMRT technology and Oxford Nanopore technology, and compare the two techniques against Illumina, the predominant next-generation sequencing technique. In addition, we summarize the nuclear and organelle genome assemblies of commonly used medicinal plants, choose several examples from genomics, transcriptomics, and molecular identification studies to dissect the specific processes and summarize the advantages and disadvantages of the two TGS techniques when applied to medicinal organisms. Finally, we describe how we expect that TGS techniques will be widely utilized to assemble telomere-to-telomere (T2T) genomes and in epigenomics research involving medicinal plants.

The complete chloroplast genome sequence of the medicinal plant, Dracocephalum rupestre (Lamiaceae)

Dracocephalum rupestre (tribe Mentheae; Lamiaceae) is a perennial herb from Korea and China with high ornamental and medicinal value. Here, we report its complete chloroplast genome to provide insight into the phylogenetic relationships of Dracocephalum. The genome is 151,230 bp long, with two inverted repeat regions (25,643 bp each) that separate a large single-copy region (82,536 bp) and a small single-copy region (17,408 bp). It contains 133 genes that encode 88 proteins, eight rRNAs, and 37 tRNAs. The maximum-likelihood phylogenetic analysis strongly supported Dracocephalum monophyly, showing that the genus forms a sister group with the subtribe Menthinae in the tribe Mentheae.

The complete chloroplast genome of Crataegus scabrifolia (Franch.) Rehd (Rosaceae), a medicinal and edible plant in Southwest China

Crataegus scabrifolia (Franch.) Rehd is a medicinal and edible plant in Southwest China. The chloroplast genome of C. scabrifolia was analyzed by high-throughput sequencing technology, and its genetic relationship to related species was discussed. The chloroplast genome is 159,637 bp long, with two inverted repeat (IR) regions (26,384 bp each) that separate a large single-copy (LSC) region (87,730 bp) and a small single-copy (SSC) region (19,139 bp). A total of 127 genes were annotated, including 83 protein-coding genes, 8 rRNA genes, and 36 tRNA genes. The phylogenetic tree shows that C. hupehensis is closely related to C. scabrifolia with strong bootstrap support.

A reappraisal of the phylogeny and historical biogeography of Sparganium (Typhaceae) using complete chloroplast genomes

BACKGROUND

Sparganium (Typhaceae) is a widespread temperate genus of ecologically important aquatic plants. Previous reconstructions of the phylogenetic relationships among Sparganium species are incompletely resolved partly because they were based on molecular markers comprising < 7,000 bp. Here, we sequenced and assembled the complete chloroplast genomes from 19 Sparganium samples representing 15 putative species and three putative subspecies in order to explore chloroplast genome evolution in this genus, clarify taxonomic lineages, estimate the divergence times of Sparganium species, and reconstruct aspects of the biogeographic history of the genus.

RESULTS

The 19 chloroplast genomes shared a conserved genome structure, gene content, and gene order. Our phylogenomic analysis presented a well-resolved phylogeny with robust support for most clades. Non-monophyly was revealed in three species: S. erectum, S. eurycarpum, and S. stoloniferum. Divergence time estimates suggest that the two subgenera of Sparganium split from each other ca. 30.67 Ma in the middle Oligocene. The subgenus Xanthosparganium diversified in the late Oligocene and Miocene, while the subgenus Sparganium diversified in the late Pliocene and Pleistocene. Ancestral area reconstruction suggested that the two subgenera may have originated in East Eurasia and North America.

CONCLUSION

The non-monophyletic nature of three putative species underscores the necessity of taxonomic revision for Sparganium: S. stoloniferum subsp. choui may be more appropriately identified as S. choui, and subspecies of S. erectum may be in fact distinct species. The estimated diversification times of the two subgenera correspond to their species and nucleotide diversities. The likely ancestral area for most of subgenus Xanthosparganium was East Eurasia and North America from where it dispersed into West Eurasia and Australia. Most of subgenus Sparganium likely originated in North America and then dispersed into Eurasia. Our study demonstrates some of the ways in which complete chloroplast genome sequences can provide new insights into the evolution, phylogeny, and biogeography of the genus Sparganium.

Evidence from Phylogenomics and Morphology Provide Insights into the Phylogeny, Plastome Evolution, and Taxonomy of Kitagawia

Kitagawia Pimenov is one of the segregate genera of Peucedanum sensu lato within the Apiaceae. The phylogenetic position and morphological delimitation of Kitagawia have been controversial. In this study, we used plastid genome (plastome) and nuclear ribosomal DNA (nrDNA) sequences to reconstruct the phylogeny of Kitagawia, along with comparative plastome and morphological analyses between Kitagawia and related taxa. The phylogenetic results identified that all examined Kitagawia species were divided into Subclade I and Subclade II within the tribe Selineae, and they were all distant from the representative members of Peucedanum sensu stricto. The plastomes of Kitagawia and related taxa showed visible differences in the LSC/IRa junction (JLA) and several hypervariable regions, which separated Subclade I and Subclade II from other taxa. Fruit anatomical and micromorphological characteristics, as well as general morphological characteristics, distinguished the four Kitagawia species within Subclade I from Subclade II and other related genera. This study supported the separation of Kitagawia from Peucedanum sensu lato, confirmed that Kitagawia belongs to Selineae, and two species (K. praeruptora and K. formosana) within Subclade II should be placed in a new genus. We believe that the "core" Kitagawia should be limited to Subclade I, and this genus can be distinguished by the association of a series of morphological characteristics. Overall, our study provides new insights into the phylogeny, plastome evolution, and taxonomy of Kitagawia.

The entire chloroplast genome sequence of Asparagus setaceus (Kunth) Jessop: Genome structure, gene composition, and phylogenetic analysis in Asparagaceae

Asparagus setaceus (Kunth) Jessop is a horticultural plant of the genus Asparagus. Herein, the whole chloroplast (cp) genome of A. setaceus was sequenced with PacBio and Illumina sequencing systems. The cp genome shows a characteristic quadripartite structure with 158,076 bp. In total, 135 genes were annotated, containing 89 protein-coding, 38 tRNA, and 8 rRNA genes. Contrast with the previous cp genome of A. setaceus registered in NCBI, we identified 7 single-nucleotide polymorphisms and 15 indels, mostly situated in noncoding areas. Meanwhile, 36 repeat structures and 260 simple sequence repeats were marked out. A bias for A/T-ending codons was shown in this cp genome. Furthermore, we predicted 78 RNA-editing sites in 29 genes, which were all for C-to-U transitions. And it was also proven that positive selection was exerted on the rpoC1 gene of A. setaceus with the K a/K s data. Meanwhile, a conservative gene order and highly similar sequences of protein-coding genes were revealed within Asparagus species. Phylogenetic tree analysis indicated that A. setaceus was a sister to Asparagus cochinchinensis. Taken together, our released genome provided valuable information for the gene composition, genetics comparison, and the phylogeny studies of A. setaceus.

Comparative analysis of medicinal plant Isodon rubescens and its common adulterants based on chloroplast genome sequencing

Isodon rubescens (Hemsley) H. Hara is the source of Donglingcao under the monograph Rabdosiae Rubescentis Herba in Chinese Pharmacopoeia. In the local marketplace, this medicine can be accidentally contaminated, deliberately substituted, or mixed with other related species. The contaminants of herbal products are a threat to consumer safety. Due to the scarcity of genetic information on Isodon plants, more molecular markers are needed to avoid misidentification. In the present study, the complete chloroplast (cp) genome of seven species of Isodon was sequenced, de novo assembled and characterized. The cp genomes of these species universally exhibited a conserved quadripartite structure, i.e., two inverted repeats (IRs) containing most of the ribosomal RNA genes and two unique regions (large single copy and small single copy). Moreover, the genome structure, codon usage, and repeat sequences were highly conserved and showed similarities among the seven species. Five highly variable regions (trnS-GCU-trnT-CGU, atpH-atpI, trnE-UUC-trnT-GGU, ndhC-trnM-CAU, and rps15-ycf1) might be potential molecular markers for identifying I. rubescens and its contaminants. These findings provide valuable information for further species identification, evolution, and phylogenetic research of Isodon.

De Novo Hybrid Assembly of the Salvia miltiorrhiza Mitochondrial Genome Provides the First Evidence of the Multi-Chromosomal Mitochondrial DNA Structure of Salvia Species

Salvia miltiorrhiza has been an economically important medicinal plant. Previously, an S. miltiorrhiza mitochondrial genome (mitogenome) assembled from Illumina short reads, appearing to be a single circular molecule, has been published. Based on the recent reports on the plant mitogenome structure, we suspected that this conformation does not accurately represent the complexity of the S. miltiorrhiza mitogenome. In the current study, we assembled the mitogenome of S. miltiorrhiza using the PacBio and Illumina sequencing technologies. The primary structure of the mitogenome contained two mitochondrial chromosomes (MC1 and MC2), which corresponded to two major conformations, namely, Mac1 and Mac2, respectively. Using two approaches, including (1) long reads mapping and (2) polymerase chain reaction amplification followed by Sanger sequencing, we observed nine repeats that can mediate recombination. We predicted 55 genes, including 33 mitochondrial protein-coding genes (PCGs), 3 rRNA genes, and 19 tRNA genes. Repeat analysis identified 112 microsatellite repeats and 3 long-tandem repeats. Phylogenetic analysis using the 26 shared PCGs resulted in a tree that was congruent with the phylogeny of Lamiales species in the APG IV system. The analysis of mitochondrial plastid DNA (MTPT) identified 16 MTPTs in the mitogenome. Moreover, the analysis of nucleotide substitution rates in Lamiales showed that the genes atp4, ccmB, ccmFc, and mttB might have been positively selected. The results lay the foundation for future studies on the evolution of the Salvia mitogenome and the molecular breeding of S. miltiorrhiza.

Comparative Genomics and Phylogenetic Analysis of the Chloroplast Genomes in Three Medicinal Salvia Species for Bioexploration

To systematically determine their phylogenetic relationships and develop molecular markers for species discrimination of Salvia bowleyana, S. splendens, and S. officinalis, we sequenced their chloroplast genomes using the Illumina Hiseq 2500 platform. The chloroplast genomes length of S. bowleyana, S. splendens, and S. officinalis were 151,387 bp, 150,604 bp, and 151,163 bp, respectively. The six genes ndhB, rpl2, rpl23, rps7, rps12, and ycf2 were present in the IR regions. The chloroplast genomes of S. bowleyana, S. splendens, and S. officinalis contain 29 tandem repeats; 35, 29, 24 simple-sequence repeats, and 47, 49, 40 interspersed repeats, respectively. The three specific intergenic sequences (IGS) of rps16-trnQ-UUG, trnL-UAA-trnF-GAA, and trnM-CAU-atpE were found to discriminate the 23 Salvia species. A total of 91 intergenic spacer sequences were identified through genetic distance analysis. The two specific IGS regions (trnG-GCC-trnM-CAU and ycf3-trnS-GGA) have the highest K2p value identified in the three studied Salvia species. Furthermore, the phylogenetic tree showed that the 23 Salvia species formed a monophyletic group. Two pairs of genus-specific DNA barcode primers were found. The results will provide a solid foundation to understand the phylogenetic classification of the three Salvia species. Moreover, the specific intergenic regions can provide the probability to discriminate the Salvia species between the phenotype and the distinction of gene fragments.

Complete chloroplast genome sequencing and comparative analysis of threatened dragon trees Dracaena serrulata and Dracaena cinnabari

Dracaena (Asparagaceae family) tree is famous for producing "dragon blood"—a bioactive red-colored resin. Despite its long history of use in traditional medicine, little knowledge exists on the genomic architecture, phylogenetic position, or evolution. Hence, in this study, we sequenced the whole chloroplast (cp) genomes of D. serrulata and D. cinnabari and performed comparative genomics of nine genomes of the genus Dracaena. The results showed that the genome sizes range from 155,055 (D. elliptica) to 155,449 (D. cochinchinensis). The cp genomes of D. serrulata and D. cinnabari encode 131 genes, each including 85 and 84 protein-coding genes, respectively. However, the D. hokouensis had the highest number of genes (133), with 85 protein coding genes. Similarly, about 80 and 82 repeats were identified in the cp genomes of D. serrulata and D. cinnabari, respectively, while the highest repeats (103) were detected in the cp genome of D. terniflora. The number of simple sequence repeats (SSRs) was 176 and 159 in D. serrulata and D. cinnabari cp genomes, respectively. Furthermore, the comparative analysis of complete cp genomes revealed high sequence similarity. However, some sequence divergences were observed in accD, matK, rpl16, rpoC2, and ycf1 genes and some intergenic spacers. The phylogenomic analysis revealed that D. serrulata and D. cinnabari form a monophyletic clade, sister to the remaining Dracaena species sampled in this study, with high bootstrap values. In conclusion, this study provides valuable genetic information for studying the evolutionary relationships and population genetics of Dracaena, which is threatened in its conservation status.

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.

Complete Chloroplast Genomes of 14 Subspecies of D. glomerata: Phylogenetic and Comparative Genomic Analyses

Orchardgrass (Dactylis glomerata L.) is a species in the Gramineae family that is highly important economically and valued for its role in ecology. However, the phylogeny and taxonomy of D. glomerata are still controversial based on current morphological and molecular evidence. The study of chloroplast (cp) genomes has developed into a powerful tool to develop molecular markers for related species and reveal the relationships between plant evolution and phylogenetics. In this study, we conducted comparative genomic analyses and phylogenetic inferences on 14 cp genomes of D. glomerata originating from the Mediterranean and Eurasia. The genome size ranged from 134,375 bp to 134,993 bp and exhibited synteny of gene organization and order. A total of 129–131 genes were identified, including 85–87 protein coding genes, 38 tRNA genes and 8 rRNA genes. The cp sequences were highly conserved, and key sequence variations were detected at the junctions of inverted repeats (IRs)/small single–copy (SSC) regions. Moreover, nine highly variable regions were identified among the subspecies based on a sequence divergence analysis. A total of 285 RNA editing sites were detected that were relevant to 52 genes, where rpoB exhibited the most abundant RNA editing sites. The phylogenetic analysis revealed that all Dactylis subspecies clustered into a monophyletic group and most branches provided a high support bootstrap. The main divergence time of D. glomerata was dated to the Miocene era, and this could have been due to changes in the climate. These findings will provide useful insights for further studies on phylogeny, the identification of subspecies and the development of hypotheses for the evolutionary history of the genus Dactylis and of the Gramineae family.

Newly reported chloroplast genome of Sinosenecio albonervius Y. Liu & Q. E. Yang and comparative analyses with other Sinosenecio species

Background

Sinosenecio B. Nordenstam (Asteraceae) currently comprises 44 species. To investigate the interspecific relationship, several chloroplast markers, including ndhC-trnV, rpl32-trnL, matK, and rbcL, are used to analyze the phylogeny of Sinosenecio. However, the chloroplast genomes of this genus have not been thoroughly investigated. We sequenced and assembled the Sinosenecio albonervius chloroplast genome for the first time. A detailed comparative analysis was performed in this study using the previously reported chloroplast genomes of three Sinosenecio species.

Results

The results showed that the chloroplast genomes of four Sinosenecio species exhibit a typical quadripartite structure. There are equal numbers of total genes, protein-coding genes and RNA genes among the annotated genomes. Per genome, 49–56 simple sequence repeats and 99 repeat sequences were identified. Thirty codons were identified as RSCU values greater than 1 in the chloroplast genome of S. albonervius based on 54 protein-coding genes, indicating that they showed biased usage. Among 18 protein-coding genes, 46 potential RNA editing sites were discovered. By comparing these chloroplast genomes' structures, inverted repeat regions and coding regions were more conserved than single-copy and non-coding regions. The junctions among inverted repeat and single-copy regions showed slight difference. Several hot spots of genomic divergence were detected, which can be used as new DNA barcodes for species identification. Phylogenetic analysis of the whole chloroplast genome showed that the four Sinosenecio species have close interspecific relationships.

Conclusions

The complete chloroplast genome of Sinosenecio albonervius was revealed in this study, which included a comparison of Sinosenecio chloroplast genome structure, variation, and phylogenetic analysis for related species. These will help future research on Sinosenecio taxonomy, identification, origin, and evolution to some extent.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08872-3.

Revealing the Complete Chloroplast Genome of an Andean Horticultural Crop, Sweet Cucumber (Solanum muricatum), and Its Comparison with Other Solanaceae Species

Sweet cucumber (Solanum muricatum) sect. Basarthrum is a neglected horticultural crop native to the Andean region. It is naturally distributed very close to other two Solanum crops of high importance, potatoes, and tomatoes. To date, molecular tools for this crop remain undetermined. In this study, the complete sweet cucumber chloroplast (cp) genome was obtained and compared with seven Solanaceae species. The cp genome of S. muricatum was 155,681 bp in length and included a large single copy (LSC) region of 86,182 bp and a small single-copy (SSC) region of 18,360 bp, separated by a pair of inverted repeats (IR) regions of 25,568 bp. The cp genome possessed 87 protein-coding genes (CDS), 37 transfer RNA (tRNA) genes, eight ribosomal RNA (rRNA) genes, and one pseudogene. Furthermore, 48 perfect microsatellites were identified. These repeats were mainly located in the noncoding regions. Whole cp genome comparative analysis revealed that the SSC and LSC regions showed more divergence than IR regions. Similar to previous studies, our phylogenetic analysis showed that S. muricatum is a sister species to members of sections Petota + Lycopersicum + Etuberosum. We expect that this first sweet cucumber chloroplast genome will provide potential molecular markers and genomic resources to shed light on the genetic diversity and population studies of S. muricatum, which will allow us to identify varieties and ecotypes. Finally, the features and the structural differentiation will provide us with information about the genes of interest, generating tools for the most precise selection of the best individuals of sweet cucumber, in less time and with fewer resources.

The complete chloroplast genomes of Tetrastigma hemsleyanum (Vitaceae) from different regions of China: molecular structure, comparative analysis and development of DNA barcodes for its geographical origin discrimination

BACKGROUND

Tetrastigma hemsleyanum is a valuable traditional Chinese medicinal plant widely distributed in the subtropical areas of China. It belongs to the Cayratieae tribe, family Vitaceae, and exhibited significant anti-tumor and anti-inflammatory activities. However, obvious differences were observed on the quality of T. hemsleyanum root from different regions, requiring the discrimination strategy for the geographical origins.

RESULT

This study characterized five complete chloroplast (cp) genomes of T. hemsleynum samples from different regions, and conducted a comparative analysis with other representing species from family Vitaceae to reveal the structural variations, informative markers and phylogenetic relationships. The sequenced cp genomes of T. hemsleyanum exhibited a conserved quadripartite structure with full length ranging from 160,124 bp of Jiangxi Province to 160,618 bp of Zhejiang Province. We identified 112 unique genes (80 protein-coding, 28 tRNA and 4 rRNA genes) in the cp genomes of T. hemsleyanum with highly similar gene order, content and structure. The IR contraction/expansion events occurred on the junctions of ycf1, rps19 and rpl2 genes with different degrees, causing the differences of genome sizes in T. hemsleyanum and Vitaceae plants. The number of SSR markers discovered in T. hemsleyanum was 56-57, exhibiting multiple differences among the five geographic groups. Phylogenetic analysis based on conserved cp genome proteins strongly grouped the five T. hemsleyanum species into one clade, showing a sister relationship with T. planicaule. Comparative analysis of the cp genomes from T. hemsleyanum and Vitaceae revealed five highly variable spacers, including 4 intergenic regions and one protein-coding gene (ycf1). Furthermore, five mutational hotspots were observed among T. hemsleyanum cp genomes from different regions, providing data for designing DNA barcodes trnL and trnN. The combination of molecular markers of trnL and trnN clustered the T. hemsleyanum samples from different regions into four groups, thus successfully separating specimens of Sichuan and Zhejiang from other areas.

CONCLUSION

Our study obtained the chloroplast genomes of T. hemsleyanum from different regions, and provided a potential molecular tracing tool for determining the geographical origins of T. hemsleyanum, as well as important insights into the molecular identification approach and and phylogeny in Tetrastigma genus and Vitaceae family.