Complete chloroplast genomes of all six Hosta species occurring in Korea: molecular structures, comparative, and phylogenetic analyses

Comparative chloroplast genome analyses of cultivated and wild Capsicum species shed light on evolution and phylogeny

The chloroplast (cp.) genome, also known as plastome, plays crucial roles in plant survival, adaptation, and evolution. The stable genetic structure of cp. genomes provides an ideal system for investigating species evolution. We sequenced three complete cp. genome sequences of Capsicum species and analyzed them using sequences of various Capsicum species retrieved from the NCBI database. The cp. genome of Capsicum species maintains a well-preserved quadripartite structure consisting of two inverted repeats (IRs) flanked by a large single copy (LSC) region and a small single copy (SSC) region. The sizes of cp. genome sequences ranged from 156,583 bp (C. lycianthoides) to 157,390 bp (C.pubescens). A total of 127-132 unique genes, including 83-87 protein-coding, 36-37 tRNA, and eight rRNA genes, were predicted. Comparison of cp. genomes of 10 Capsicum species revealed high sequence similarity in genome-wide organization and gene arrangements. Fragments of trnT-UGU/trnL-UAA, ccsA, ndhD, rps12, and ycf1 were identified as variable regions, and nucleotide variability of LSC and SSC was higher than that of IR. Phylogenetic speciation analysis showed that the major domesticated C. annuum species were the most extensively divergent species and closely related to C. tovarii and C. frutescens. Analysis of divergent times suggested that a substantial range of speciation events started occurring ~ 25.79 million years ago (Mya). Overall, comparative analysis of cp. genomes of Capsicum species not only offers new insights into their genetic variation and phylogenetic relationships, but also lays a foundation for evolutionary history, genetic diversity, conservation, and biological breeding of Capsicum species.

Plastome Evolution, Phylogenomics, and DNA Barcoding Investigation of Gastrochilus (Aeridinae, Orchidaceae), with a Focus on the Systematic Position of Haraella retrocalla

Gastrochilus is an orchid genus containing about 70 species in tropical and subtropical Asia with high morphological diversity. The phylogenetic relationships among this genus have not been fully resolved, and the plastome evolution has not been investigated either. In this study, five plastomes of Gastrochilus were newly reported, and sixteen plastomes of Gastrochilus were used to conduct comparative and phylogenetic analyses. Our results showed that the Gastrochilus plastomes ranged from 146,183 to 148,666 bp, with a GC content of 36.7-36.9%. There were 120 genes annotated, consisting of 74 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. No contraction and expansion of IR borders, gene rearrangements, or inversions were detected. Additionally, the repeat sequences and codon usage bias of Gastrochilus plastomes were highly conserved. Twenty hypervariable regions were selected as potential DNA barcodes. The phylogenetic relationships within Gastrochilus were well resolved based on the whole plastome, especially among main clades. Furthermore, both molecular and morphological data strongly supported Haraella retrocalla as a member of Gastrochilus (G. retrocallus).

Characterization, comparison, and phylogenetic analyses of chloroplast genomes of Euphorbia species

The genus Euphorbia (Euphorbiaceae) has near-cosmopolitan distribution and serves as a significant resource for both ornamental and medicinal purposes. Despite its economic importance, Euphorbia's taxonomy has long been challenged by the intricate nature of morphological traits exhibiting high levels of convergence. While molecular markers are essential for phylogenetic studies, their availability for Euphorbia has been limited. To address this gap, we conducted comparative analyses focusing on the chloroplast (CP) genomes of nine Euphorbia species, incorporating three newly sequenced and annotated accessions. In addition, phylogenetic informativeness and nucleotide diversity were computed to identify candidate markers for phylogenetic analyses among closely related taxa in the genus. Our investigation revealed relatively conserved sizes and structures of CP genomes across the studied species, with notable interspecific variations observed primarily in non-coding regions and IR/SC borders. By leveraging phylogenetic informativeness and nucleotide diversity, we identified rpoB gene as the optimal candidate for species delimitation and shallow-level phylogenetic inference within the genus. Through this comprehensive analysis of CP genomes across multiple taxa, our study sheds light on the evolutionary dynamics and taxonomic intricacies of Euphorbia, offering valuable insights into its CP genome evolution and taxonomy.

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.

Molecular phylogeny and taxonomy of Hosta (Asparagaceae) on Shikoku Island, Japan, including five new species, one new subspecies, and two new status assignments

Japan has 16 native species of the genus Hosta Tratt. (Asparagaceae). A recent study on Hosta based on field surveys and molecular phylogenetic analyses resulted in the discovery of six unknown taxa in Kochi Prefecture, Shikoku Island, southwestern Japan. We aimed to identify these unknown taxa. Therefore, we constructed a finely resolved phylogeny for 320 Hosta samples collected from the Honshu, Shikoku, and Kyushu Islands using multiplex inter-simple sequence repeat genotyping by sequencing (MIG-seq). Based on this phylogenetic analysis and related morphological observations, we describe five new species, H.longipedicellatasp. nov., H.minazukiflorasp. nov., H.polyneuronoidessp. nov., H.samukazemontanasp. nov., and H.takiminazukiflorasp. nov. and one new subspecies, H.takiminazukiflorasubsp.grandissubsp. nov. In addition, we propose two new status assignments, H.tardivasubsp.densinerviacomb. and stat. nov. and H.scabrinerviastat. nov. We also propose classifying H.kikutiivar.tosana as a species, H.tosana. Further studies that combine MIG-seq with careful morphological observations are needed for Hosta plants on all Japanese islands, which may result in the discovery of even more undescribed species.

Plastome evolution and phylogenomics of Trichosporeae (Gesneriaceae) with its morphological characters appraisal

Trichosporeae is the largest and most taxonomically difficult tribe of Gesneriaceae due to its diverse morphology. Previous studies have not clarified the phylogenetic relationships within this tribe on several DNA markers, including the generic relationships within subtribes. Recently, plastid phylogenomics have been successfully employed to resolve the phylogenetic relationships at different taxonomic levels. In this study, plastid phylogenomics were used to explore the relationships within Trichosporeae. Eleven plastomes of Hemiboea were newly reported. Comparative analyses, phylogeny and morphological character evolution within Trichosporeae were conducted on 79 species representing seven subtribes. The Hemiboea plastomes range from 152,742 bp to 153,695 bp in length. Within Trichosporeae, the sampled plastomes range from 152,196 bp to 156,614 bp and GC content from 37.2% to 37.8%. A total of 121-133 genes were annotated in each species, including 80-91 protein-coding genes, 34-37 tRNA genes, and 8 rRNA genes. The contraction and expansion of IR borders were not detected, and gene rearrangements and inversions did not occur. The 13 hypervariable regions were proposed as the potential molecular markers for species identification. A total of 24,299 SNPs and 3,378 indels were inferred, and most of the SNPs were functionally missense and silent variations. There were 1968 SSRs, 2055 tandem repeats and 2802 dispersed repeats. The RSCU and ENC values indicated that the codon usage pattern was conserved in Trichosporeae. Both the phylogenetic frameworks based on the whole plastome and 80 CDSs were basically concordant. The sister relationships between Loxocarpinae and Didymocarpinae were confirmed, and Oreocharis was a sister group of Hemiboea with high support. The morphological characters showed a complex evolutionary pattern of Trichosporeae. Our findings may contribute to future research on genetic diversity, morphological evolutionary patterns, and conservation of the tribe Trichosporeae.

Comparative and phylogenetic analysis of the complete chloroplast genomes of six Polygonatum species (Asparagaceae)

Polygonatum Miller belongs to the tribe Polygonateae of Asparagaceae. The horizontal creeping fleshy roots of several species in this genus serve as traditional Chinese medicine. Previous studies have mainly reported the size and gene contents of the plastomes, with little information on the comparative analysis of the plastid genomes of this genus. Additionally, there are still some species whose chloroplast genome information has not been reported. In this study, the complete plastomes of six Polygonatum were sequenced and assembled, among them, the chloroplast genome of P. campanulatum was reported for the first time. Comparative and phylogenetic analyses were then conducted with the published plastomes of three related species. Results indicated that the whole plastome length of the Polygonatum species ranged from 154,564 bp (P. multiflorum) to 156,028 bp (P. stenophyllum) having a quadripartite structure of LSC and SSC separated by two IR regions. A total of 113 unique genes were detected in each of the species. Comparative analysis revealed that gene content and total GC content in these species were highly identical. No significant contraction or expansion was observed in the IR boundaries among all the species except P. sibiricum1, in which the rps19 gene was pseudogenized owing to incomplete duplication. Abundant long dispersed repeats and SSRs were detected in each genome. There were five remarkably variable regions and 14 positively selected genes were identified among Polygonatum and Heteropolygonatum. Phylogenetic results based on chloroplast genome strongly supported the placement of P. campanulatum with alternate leaves in sect. Verticillata, a group characterized by whorled leaves. Moreover, P. verticillatum and P. cyrtonema were displayed as paraphyletic. This study revealed that the characters of plastomes in Polygonatum and Heteropolygonatum maintained a high degree of similarity. Five highly variable regions were found to be potential specific DNA barcodes in Polygonatum. Phylogenetic results suggested that leaf arrangement was not suitable as a basis for delimitation of subgeneric groups in Polygonatum and the definitions of P. cyrtonema and P. verticillatum require further study.

Decoding the complete chloroplast genome of Cissus quadrangularis: insights into molecular structure, comparative genome analysis and mining of mutational hotspot regions

Cissus quadrangularis L., a member of the Vitaceae family, is an important medicinal plant with widespread application in Indian traditional medicines. C. quadrangularis L. whole chloroplast genome of 160,404 bp was assembled using a genome skimming approach from the whole genome library. The assembled chloroplast genome contained a large single-copy region (88,987 bp), a small single-copy region (18,621 bp), and pairs of inverted repeat regions (26,398 bp). It also comprised 133 genes, including 37 tRNAs, eight rRNAs, and 88 protein-coding genes. Aside from that, we annotated three genes atpH, petB, and psbL, as well as one duplicated copy of the ycf1 gene in C. quadrangularis L. that had previously been missing from the annotation of compared Cissus chloroplast genomes. Five divergent hotspot regions such as petA_psbJ (0.1237), rps16_trnQ-UUG (0.0913), psbC_trnS-UGA (0.0847), rps15_ycf1 (0.0788), and rps2_rpoC2 (0.0788) were identified in the investigation that could aid in future species discrimination. Surprisingly, we found the overlapping genes ycf1 and ndhF on the IRb/SSC junction, rarely seen in angiosperms. The results of the phylogenetic study showed that the genomes of the Cissus species under study formed a single distinct clade. The detailed annotations given in this study could be useful in the future for genome annotations of Cissus species. The current findings of the study have the potential to serve as a useful resource for future research in the field of population genetics and the evolutionary relationships in the Cissus genus.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01312-w.

Plastome evolution in the East Asian lobelias (Lobelioideae) using phylogenomic and comparative analyses

Lobelia species, as rich source of the alkaloid lobeline which has been shown to have important biological activity, have been used in folk medicine throughout East Asia to treat various diseases. However, Lobelia is a complex and varied genus in East Asia and is thus difficult to identify. Genomic resources would aid identification, however the availability of such information is poor, preventing a clear understanding of their evolutionary history from being established. To close this gap in the available genomic data, in this study, 17 plastomes of East Asian lobelias were newly sequenced and assembled. Although the plastomes of Lobelia sect. Hypsela, L. sect. Speirema, and L. sect. Rhynchopetalum shared the gene structure, the inverted repeat (IR)/large single copy (LSC) boundaries, genome size, and the number of repeats were variable, indicating the non-conservative nature of plastome evolution within these sections. However, the genomes of the Lobelia sect. Delostemon and L. sect. Stenotium showed rearrangements, revealing that these two sections might have undergone different evolutionary histories. We assessed nine hotspot genes and 27-51 simple sequence repeat motifs, which will also serve as valuable DNA barcode regions in future population genetics studies and for the delineation of plant species. Our phylogenetic analysis resolved the evolutionary positions of the five sections in agreement with previous evolutionary trees based on morphological features. Although phylogenetic reconstruction of Lobelioideae based on the rpoC2 gene has rarely been performed, our results indicated that it contains a considerable amount of phylogenetic information and offers great promise for further phylogenetic analysis of Lobelioideae. Our site-specific model identified 173 sites under highly positive selections. The branch-site model exhibited 11 positive selection sites involving four genes in the East Asian branches. These four genes may play critical roles in the adaptation of East Asian taxa to diverse environments. Our study is the first to detect plastome organization, phylogenetic utility, and signatures of positive selection in the plastomes of East Asian lobelias, which will help to further advance taxonomic and evolutionary studies and the utilization of medicinal plant resources.

Complete Chloroplast Genome Sequences of Four Species in the Caladium Genus: Comparative and Phylogenetic Analyses

Caladiums are promising colorful foliage plants due to their dazzling colors of the leaves, veins, stripes, and patches, which are often cultivated in pots or gardens as decorations. Four wild species, including C. bicolor, C. humboldtii, C. praetermissum, and C. lindenii, were employed in this study, where their chloroplast (cp) genomes were sequenced, assembled, and annotated via high-throughput sequencing. The whole cp genome size ranged from 162,776 bp to 168,888 bp, and the GC contents ranged from 35.09% to 35.91%. Compared with the single large copy (LSC) and single small copy (SSC) regions, more conserved sequences were identified in the inverted repeat regions (IR). We further analyzed the different region borders of nine species of Araceae and found the expansion or contraction of IR/SSC regions might account for the cp genome size variation. Totally, 131 genes were annotated in the cp genomes, including 86 protein-coding genes (PCGs), 37 tRNAs, and eight rRNAs. The effective number of codons (ENC) values and neutrality plot analyses provided the foundation that the natural selection pressure could greatly affect the codon preference. The GC₃ content was significantly lower than that of GC₁ and GC₂, and codons ending with A/U had higher usage preferences. Finally, we conducted phylogenetic relationship analysis based on the chloroplast genomes of twelve species of Araceae, in which C. bicolor and C. humboldtii were grouped together, and C. lindenii was furthest from the other three Caladium species occupying a separate branch. These results will provide a basis for the identification, development, and utilization of Caladium germplasm.

Complete chloroplast genomes of Sorbus sensu stricto (Rosaceae): comparative analyses and phylogenetic relationships

BACKGROUND

Sorbus sensu stricto (Sorbus s.s.) is a genus with important economical values because of its beautiful leaves, and flowers and especially the colorful fruits. It belongs to the tribe Maleae of the family Rosaceae, and comprises about 90 species mainly distributed in China. There is on-going dispute about its infrageneric classification and species delimitation as the species are morphologically similar. With the aim of shedding light on the circumscription of taxa within the genus, phylogenetic analyses were performed using 29 Sorbus s.s. chloroplast (cp) genomes (16 newly sequenced) representing two subgenera and eight sections.

RESULTS

The 16 cp genomes newly sequenced range between 159,646 bp and 160,178 bp in length. All the samples examined and 22 taxa re-annotated in Sorbus sensu lato (Sorbus s.l.) contain 113 unique genes with 19 of these duplicated in the inverted repeat (IR). Six hypervariable regions including trnR-atpA, petN-psbM, rpl32-trnL, trnH-psbA, trnT-trnL and ndhC-trnV were screened and 44-53 SSRs and 14-31 dispersed repeats were identified as potential molecular markers. Phylogenetic analyses under ML/BI indicated that Sorbus s.l. is polyphyletic, but Sorbus s.s. and the other five segregate genera, Aria, Chamaemespilus, Cormus, Micromeles and Torminalis are monophyletic. Two major clades and four sub-clades resolved with full-support within Sorbus s.s. are not consistent with the existing infrageneric classification. Two subgenera, subg. Sorbus and subg. Albocarmesinae are supported as monophyletic when S. tianschanica is transferred to subg. Albocarmesinae from subg. Sorbus and S. hupehensis var. paucijuga transferred to subg. Sorbus from subg. Albocarmesinae, respectively. The current classification at sectional level is not supported by analysis of cp genome phylogeny.

CONCLUSION

Phylogenomic analyses of the cp genomes are useful for inferring phylogenetic relationships in Sorbus s.s. Though genome structure is highly conserved in the genus, hypervariable regions and repeat sequences used are the most promising molecule makers for population genetics, species delimitation and phylogenetic studies.

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.

Complete chloroplast genomes and comparative analyses of Hippeastrum ‘milady’, Hippeastrum albertii and Hippeastrum reticulatum (Amaryllidaceae)

Hippeastrum is a genus of ornamental plants with large, brightly colored flowers. Due to the very high seed-setting rate of the hybridization of Hippeastrum, the large population of hybrid progeny and the existence of superparent inheritance, it is difficult to trace the origin of the varieties collected from the market during breeding. In this study, we analyzed the chloroplast genomes of Hippeastrum ‘Milady’, H. alberti, and H. reticulatum using the Illumina NovaSeq sequencing platform and generated full-length sequences of 158,067, 158,067, and 158,522 bp, respectively. All three genomes had the typical tetrad structure. The large single copy, small single copy, and inverted repeat regions of H. reticulatum were observed to be respectively 277, 138, and 20 bp longer than the corresponding regions of H. ‘Milady’ and H. alberti. The results of comparative analysis of simple sequence repeats (SSRs), Ka/Ks ratios, codon preferences, and complete sequences of chloroplasts of these three taxa and 14 other plant species were as follows. First, the chloroplast genomes of H. ‘Milady’, H. alberti, and H. reticulatum contain 209, 209, and 211 SSR sites, respectively, most of which (123, 123, and 122, respectively) are single nucleotide repeats. Second, leucine, arginine, and serine are the most frequently used amino acids in the three chloroplast genomes. Third, H. ‘Milady’, H. alberti, and H. reticulatum are more closely related to Lycoris and Narcissus than to Allium and Agapanthus. Our results will provide information on the study of origins or relatedness of native species, and the identification of cultivars.

Comparative analysis of the complete plastid genomes in Prunus subgenus Cerasus (Rosaceae): Molecular structures and phylogenetic relationships

Prunus subgenus Cerasus (cherry) is an economically important group that distributed in temperate regions of the northern hemisphere. However, shared interspecific morphological traits and variability across taxa of Cerasus are among the impediments to taxonomic efforts to correctly delimit taxa. This is further complicated by a lack of genetic information on these taxa, with no focused genomic or phylogenetic studies being done on Cerasus. In this study, we conducted comparative analysis on the complete plastid genomes (plastomes) of 20 Cerasus species to gain a greater understanding of the attributes of the plastome of these taxa while helping resolve their phylogenetic placement in Prunus sensu lato and interspecific relationships within the subgenus. Our results displayed that (1) the plastomes of the 20 Cerasus species studied exhibited a typical quadripartite structure with conversed genome arrangement, structure, and moderate divergence. (2) The average size of complete plastomes for the Cerasus taxa studied was 157,861 bp, ranging from 157,458 to 158,024 bp. A total of 134 genes were annotated, including 86 protein-coding genes, 40 tRNAs, and 8 rRNAs across all species. In simple sequence repeat analysis, we found Cerasus had a comparable number of dispersed and tandem repeats to those identified in other angiosperm taxa, with only P. pseudocerasus found to contain trinucleotide repeats. Nucleotide diversity analysis revealed that the trnG-GCC gene and rpl32-trnL region had the highest Pi value showing potential as phylogenetic markers. (3) Two phylogenetic trees of the plastomes verified the monophyletic relationship of Cerasus and provided a more resolved species-level phylogeny. Our study provides detailed plastome information for exploring the phylogeny of subg. Cerasus taxa. We identified various types of repeats and nucleotide diversity hotspots, which can be a reference for species identification and reconstruction of phylogenetic relationships.

Comparative Chloroplast Genomics and Phylogenetic Analysis of Thuniopsis and Closely Related Genera within Coelogyninae (Orchidaceae)

The genus Thuniopsis was recently proposed for a rare orchid species T. cleistogama formerly classified in the genus Thunia. The relationships between Thuniopsis and its related genera have not yet been conclusively resolved. Recognition of the genus provides a new perspective to illustrate the morphological diversity and plastome evolution within Coelogyninae. In this study, we sequenced and assembled complete chloroplast (cp) genomes for three accessions of Thuniopsis cleistogama and two accessions of Thunia alba. A total of 135 genes were annotated for each cp genome, including 89 protein-coding genes, 38 tRNA genes, and eight rRNA genes. The ENC-plot and neutrality plot analyses revealed that natural selection dominated over mutation pressure in their evolutionary process. Specially, we found that selection played a vital role in shaping the codon usage in Thunia alba cp genome. General characteristics of the cp genomes were further analyzed and compared with those published plastomes of four other related species. Despite the conserved organization and structure, the whole individual cp genome size ranged from 158,394 bp to 159,950 bp. In all the examined plastomes, sequences in the inverted repeat (IR) regions were more conserved than those in the small single copy (SSC) and large single copy (LSC) regions. However, close examination identified contraction and expansion of the IR/SSC boundary regions, which might be the main reason for the cp genome size variation. Our comparative analysis of the cp genomes revealed that single nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) provided valuable information for identifying genetic variations within and among genera. Furthermore, sequence variations in the protein-coding regions were more conserved than those in the non-coding regions. We selected eight divergence hotspots with nucleotide sequence diversities (Pi values) higher than 0.08. Most of these polymorphisms were located in the intergenic regions. Phylogenomic analyses recovered largely congruent relationships among major clades and strongly supported the monophyly of Thuniopsis. The results obtained in this study can improve our understanding of the classification of this enigmatic genus. The chloroplast genomic data presented here provide valuable insights into the phylogeny and evolutionary patterns of the Coelogyninae as well as the orchids as a whole.

Disjunction and Vicariance Between East and West Asia: A Case Study on Euonymus sect. Uniloculares Based on Plastid Genome Analysis

Scientists have long been captivated by biogeographic disjunctions, and disjunctions between East Asia and North America have been particularly well-studied at the genus and family levels. By contrast, disjunctions between eastern and western Asia have received less attention. Euonymus L. is taxonomically divided into two sections based on the number of cells in anthers as follows: E. sect. Uniloculares has one-celled anthers and occurs mainly in Asia, whereas E. sect. Biloculares has two-celled anthers and is distributed globally. We used Illumina sequencing to investigate the genomes of four species in sect. Uniloculares. The chloroplast (cp) genomes are highly conserved (157,290-158,094 bp). Pseudogenisation of ndhF and intron loss in rps16 was detected. Based on the cp genomes of the four species of E. sect. Uniloculares, we propose a novel hypothesis of disjunction between eastern and western Asia. Biogeographic reconstruction and molecular dating revealed that sect. Uniloculares separated from its sect. Biloculares forebears 4.0 Mya during the Pliocene era. The radial diversification of sect. Uniloculares from East Asia and the establishment of the western Asian clade during the Pleistocene era (1.9 Mya) were the results of both dispersal and vicariance, making the section the youngest diverged clade conforming to age estimation. The centre of origin of sect. Uniloculares was determined to be in East Asia. Disjunctions and diversification between eastern and western Asia in sect. Uniloculares are thought to have been caused by changes in monsoon patterns, temperature variations, and the emergence of the Gobi Desert.

Complete chloroplast genome of Stephania tetrandra (Menispermaceae) from Zhejiang Province: insights into molecular structures, comparative genome analysis, mutational hotspots and phylogenetic relationships

BACKGROUND

The Stephania tetrandra S. Moore (S. tetrandra) is a medicinal plant belonging to the family Menispermaceae that has high medicinal value and is well worth doing further exploration. The wild resources of S. tetrandra were widely distributed in tropical and subtropical regions of China, generating potential genetic diversity and unique population structures. The geographical origin of S. tetrandra is an important factor influencing its quality and price in the market. In addition, the species relationship within Stephania genus still remains uncertain due to high morphological similarity and low support values of molecular analysis approach. The complete chloroplast (cp) genome data has become a promising strategy to determine geographical origin and understand species evolution for closely related plant species. Herein, we sequenced the complete cp genome of S. tetrandra from Zhejiang Province and conducted a comparative analysis within Stephania plants to reveal the structural variations, informative markers and phylogenetic relationship of Stephania species.

RESULTS

The cp genome of S. tetrandra voucher ZJ was 157,725 bp, consisting of a large single copy region (89,468 bp), a small single copy region (19,685 bp) and a pair of inverted repeat regions (24,286 bp each). A total of 134 genes were identified in the cp genome of S. tetrandra, including 87 protein-coding genes, 8 rRNA genes, 37 tRNA genes and 2 pseudogene copies (ycf1 and rps19). The gene order and GC content were highly consistent in the Stephania species according to the comparative analysis results, with the highest RSCU value in arginine (1.79) and lowest RSCU value in serine of S. tetrandra, respectively. A total of 90 SSRs have been identified in the cp genome of S. tetrandra, where repeats that consisting of A or T bases were much higher than that of G or C bases. In addition, 92 potential RNA editing sites were identified in 25 protein-coding genes, with the most predicted RNA editing sites in ndhB gene. The variations on length and expansion extent to the junction of ycf1 gene were observed between S. tetrandra vouchers from different regions, indicating potential markers for further geographical origin discrimination. Moreover, the values of transition to transversion ratio (Ts/Tv) in the Stephania species were significantly higher than 1 using Pericampylus glaucus as reference. Comparative analysis of the Stephania cp genomes revealed 5 highly variable regions, including 3 intergenic regions (trnH-psbA, trnD-trnY, trnP) and two protein coding genes (rps16 and ndhA). The identified mutational hotspots of Stephania plants exhibited multiple SNP sites and Gaps, as well as different Ka/Ks ratio values. In addition, five pairs of specific primers targeting the divergence regions were accordingly designed, which could be utilized as potential molecular markers for species identification, population genetic and phylogenetic analysis in Stephania species. Phylogenetic tree analysis based on the conserved chloroplast protein coding genes indicated a sister relationship between S. tetrandra and the monophyletic group of S. japonica and S. kwangsiensis with high support values, suggesting a close genetic relationship within Stephania plants. However, two S. tetrandra vouches from different regions failed to cluster into one clade, confirming the occurrences of genetic diversities and requiring further investigation for geographical tracing strategy.

CONCLUSIONS

Overall, we provided comprehensive and detailed information on the complete chloroplast genome and identified nucleotide diversity hotspots of Stephania species. The obtained genetic resource of S. tetrandra from Zhejiang Province would facilitate future studies in DNA barcode, species discrimination, the intraspecific and interspecific variability and the phylogenetic relationships of Stephania plants.

The complete chloroplast genome of Agave angustifolia

Agave angustifolia is commonly used for the production of bacanora, a kind of fermented and distilled beverage in Mexico. In the present study, we have successfully assembled its chloroplast genome. The full length of the genome is 157,274 bp with the GC content of 37.84%. There is a large single copy region (LSC) of 85,895 bp, a pair of inverted repeat regions (IR) of 26,575 bp and a small single copy region (SSC) of 18,229 bp in the genome. A total of 132 genes are annotated in the cp genome. Among these, there are 86 protein-coding genes, 38 tRNAs and 8 rRNAs. Phylogenetic analysis reveals that A. angustifolia is closely related with A. H11648.

Plastome Characterization and Phylogenomic Analysis Yield New Insights into the Evolutionary Relationships among the Species of the Subgenus Bryocles (Hosta; Asparagaceae) in East Asia

The genus Hosta, which has a native distribution in temperate East Asia and a number of species ranging from 23 to 40, represents a taxonomically important and ornamentally popular plant. Despite its taxonomic and horticultural importance, the genus Hosta has remained taxonomically challenging owing to insufficient diagnostic features, continuous morphological variation, and the process of hybridization and introgression, making species circumscription and phylogenetic inference difficult. In this study, we sequenced 11 accessions of Hosta plastomes, including members of three geographically defined subgenera, Hosta, Bryocles, and Giboshi, determined the characteristics of plastomes, and inferred their phylogenetic relationships. We found highly conserved plastomes among the three subgenera, identified several mutation hotspots that can be used as barcodes, and revealed the patterns of codon usage bias and RNA editing sites. Five positively selected plastome genes (rbcL, rpoB, rpoC2, rpl16, and rpl20) were identified. Phylogenetic analysis suggested (1) the earliest divergence of subg. Hosta, (2) non-monophyly of subg. Bryocles and its two sections (Lamellatae and Stoloniferae), (3) a sister relationship between H. sieboldiana (subg. Giboshi) and H. ventricosa (subg. Bryocles), and (4) reciprocally monophyletic and divergent lineages of H. capitata in Korea and Japan, requiring further studies of their taxonomic distinction.

Validation of Hostaalata (Asparagaceae) as a new species and its phylogenetic affinity

Molecular phylogenetic studies of Hostapulchella (Asparagaceae) and its relatives, which are native to Japan, have been conducted and resulted in a highly resolved phylogeny. Specifically, the relationship of H.pulchella to H.alata Hatusima, nom. nud. is investigated. These data include genome-wide SNPs obtained through conducting multiplexed ISSR genotyping by sequencing (MIG-seq). Based on these phylogenetic results, morphological observations, distribution, and differences in flowering periods of H.alata collections sympatric with H.pulchella, we find the two species closely related, but distinct. As such, we formally describe Hostaalata sp. nov. from the Oita Prefecture of Kyushu island, southwestern Japan.

Comparative chloroplast genome analysis of Impatiens species (Balsaminaceae) in the karst area of China: insights into genome evolution and phylogenomic implications

BACKGROUND

Impatiens L. is a genus of complex taxonomy that belongs to the family Balsaminaceae (Ericales) and contains approximately 1000 species. The genus is well known for its economic, medicinal, ornamental, and horticultural value. However, knowledge about its germplasm identification, molecular phylogeny, and chloroplast genomics is limited, and taxonomic uncertainties still exist due to overlapping morphological features and insufficient genomic resources.

RESULTS

We sequenced the chloroplast genomes of six different species (Impatiens chlorosepala, Impatiens fanjingshanica, Impatiens guizhouensis, Impatiens linearisepala, Impatiens loulanensis, and Impatiens stenosepala) in the karst area of China and compared them with those of six previously published Balsaminaceae species. We contrasted genomic features and repeat sequences, assessed sequence divergence and constructed phylogenetic relationships. Except for those of I. alpicola, I. pritzelii and I. glandulifera, the complete chloroplast genomes ranging in size from 151,366 bp (I. alpicola) to 154,189 bp (Hydrocera triflora) encoded 115 distinct genes [81 protein-coding, 30 transfer RNA (tRNA), and 4 ribosomal RNA (rRNA) genes]. Moreover, the characteristics of the long repeat sequences and simple sequence repeats (SSRs) were determined. psbK-psbI, trnT-GGU-psbD, rpl36-rps8, rpoB-trnC-GCA, trnK-UUU-rps16, trnQ-UUG, trnP-UGG-psaJ, trnT-UGU-trnL-UAA, and ycf4-cemA were identified as divergence hotspot regions and thus might be suitable for species identification and phylogenetic studies. Additionally, the phylogenetic relationships based on Maximum likelihood (ML) and Bayesian inference (BI) of the whole chloroplast genomes showed that the chloroplast genome structure of I. guizhouensis represents the ancestral state of the Balsaminaceae family.

CONCLUSION

Our study provided detailed information about nucleotide diversity hotspots and the types of repeats, which can be used to develop molecular markers applicable to Balsaminaceae species. We also reconstructed and analyzed the relationships of some Impatiens species and assessed their taxonomic statuses based on the complete chloroplast genomes. Together, the findings of the current study might provide valuable genomic resources for systematic evolution of the Balsaminaceae species.

The complete chloroplast genome of Agave fourcroydes

Agave fourcroydes (henequen) is the only cultivated Agave species in the Yucatan Peninsula, which is mainly used for fiber production. In the present study, we have successfully assembled the chloroplast (cp) genome of A. fourcroydes. The full length of the cp genome is 157,291 bp with a GC content at 37.8%. The cp genome is constructed with an inverted repeat region a (IRa) of 26,573 bp, a small single copy region (SSC) of 18,230 bp, an inverted repeat region b (IRb) of 26,573 bp and a large single copy region (LSC) of 85,915 bp. The annotation result reveals 132 genes on the cp genome, including 86 protein-coding genes, 38 tRNAs and 8 rRNAs. Phylogenetic tree reveals that A. fourcroydes is closely related with A. sisalana.

The complete chloroplast genome of Agave sisalana

Agave sisalana is one of the Agave cultivars for sisal fiber production around the tropical areas of the world. In the present study, we successfully sequenced and assembled its chloroplast genome. The full size of the genome is 157,268 bp with a GC content at 37.85%. The genome is constructed with a large single copy region (LSC, 85,894 bp), a pair of inverted repeat regions (IR, 26,573 bp) and a small single copy region (SSC, 18,228 bp). Besides, 86 protein-coding genes, 38 tRNAs and 8 rRNAs are annotated on the chloroplast genome. Phylogenetic result reveals that A. sisalana is closely related with A. americana and A. H11648.

Insights into phylogenetic relationships and genome evolution of subfamily Commelinoideae (Commelinaceae Mirb.) inferred from complete chloroplast genomes

Background

Commelinaceae (Commelinales) comprise 41 genera and are widely distributed in both the Old and New Worlds, except in Europe. The relationships among genera in this family have been suggested in several morphological and molecular studies. However, it is difficult to explain their relationships due to high morphological variations and low support values. Currently, many researchers have been using complete chloroplast genome data for inferring the evolution of land plants. In this study, we completed 15 new plastid genome sequences of subfamily Commelinoideae using the Mi-seq platform. We utilized genome data to reveal the structural variations and reconstruct the problematic positions of genera for the first time.

Results

All examined species of Commelinoideae have three pseudogenes (accD, rpoA, and ycf15), and the former two might be a synapomorphy within Commelinales. Only four species in tribe Commelineae presented IR expansion, which affected duplication of the rpl22 gene. We identified inversions that range from approximately 3 to 15 kb in four taxa (Amischotolype, Belosynapsis, Murdannia, and Streptolirion). The phylogenetic analysis using 77 chloroplast protein-coding genes with maximum parsimony, maximum likelihood, and Bayesian inference suggests that Palisota is most closely related to tribe Commelineae, supported by high support values. This result differs significantly from the current classification of Commelinaceae. Also, we resolved the unclear position of Streptoliriinae and the monophyly of Dichorisandrinae. Among the ten CDS (ndhH, rpoC2, ndhA, rps3, ndhG, ndhD, ccsA, ndhF, matK, and ycf1), which have high nucleotide diversity values (Pi > 0.045) and over 500 bp length, four CDS (ndhH, rpoC2, matK, and ycf1) show that they are congruent with the topology derived from 77 chloroplast protein-coding genes.

Conclusions

In this study, we provide detailed information on the 15 complete plastid genomes of Commelinoideae taxa. We identified characteristic pseudogenes and nucleotide diversity, which can be used to infer the family evolutionary history. Also, further research is needed to revise the position of Palisota in the current classification of Commelinaceae.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07541-1.

Two Korean Endemic Clematis Chloroplast Genomes: Inversion, Reposition, Expansion of the Inverted Repeat Region, Phylogenetic Analysis, and Nucleotide Substitution Rates

Previous studies on the chloroplast genome in Clematis focused on the chloroplast structure within Anemoneae. The chloroplast genomes of Cleamtis were sequenced to provide information for studies on phylogeny and evolution. Two Korean endemic Clematis chloroplast genomes (Clematis brachyura and C. trichotoma) range from 159,170 to 159,532 bp, containing 134 identical genes. Comparing the coding and non-coding regions among 12 Clematis species revealed divergent sites, with carination occurring in the petD-rpoA region. Comparing other Clematis chloroplast genomes suggested that Clematis has two inversions (trnH-rps16 and rps4), reposition (trnL-ndhC), and inverted repeat (IR) region expansion. For phylogenetic analysis, 71 protein-coding genes were aligned from 36 Ranunculaceae chloroplast genomes. Anemoneae (Anemoclema, Pulsatilla, Anemone, and Clematis) clades were monophyletic and well-supported by the bootstrap value (100%). Based on 70 chloroplast protein-coding genes, we compared nonsynonymous (dN) and synonymous (dS) substitution rates among Clematis, Anemoneae (excluding Clematis), and other Ranunculaceae species. The average synonymoussubstitution rates (dS)of large single copy (LSC), small single copy (SSC), and IR genes in Anemoneae and Clematis were significantly higher than those of other Ranunculaceae species, but not the nonsynonymous substitution rates (dN). This study provides fundamental information on plastid genome evolution in the Ranunculaceae.

Effect of chlorophyll biosynthesis-related genes on the leaf color in Hosta (Hosta plantaginea Aschers) and tobacco (Nicotiana tabacum L.)

BACKGROUND

'Regal Splendour' (Hosta variety) is famous for its multi-color leaves, which are useful resources for exploring chloroplast development and color changes. The expressions of chlorophyll biosynthesis-related genes (HrHEMA, HrPOR and HrCAO) in Hosta have been demonstrated to be associated with leaf color. Herein, we isolated, sequenced, and analyzed HrHEMA, HrPOR and HrCAO genes. Subcellular localization was also performed to determine the location of the corresponding enzymes. After plasmid construction, virus-induced gene silencing (VIGS) was carried out to reduce the expressions of those genes. In addition, HrHEMA-, HrPOR- and HrCAO-overexpressing tobacco plants were made to verify the genes function. Changes of transgenic tobacco were recorded under 2000 lx, 6000 lx and 10,000 lx light intensity. Additionally, the contents of enzyme 5-aminolevulinic acid (5-ALA), porphobilinogen (PBG), chlorophyll a and b (Chla and Chlb), carotenoid (Cxc), superoxide dismutase (SOD), peroxidase (POD), malondialdehyde (MDA), proline (Pro) and catalase (CAT) under different light intensities were evaluated.

RESULTS

The silencing of HrHEMA, HrPOR and HrCAO genes can induce leaf yellowing and chloroplast structure changes in Hosta. Specifically, leaves of Hosta with HrCAO silencing were the most affected, while those with HrPOR silencing were the least affected. Moreover, all three genes in tobacco were highly expressed, whereas no expression was detected in wild-type (WT). However, the sensitivities of the three genes to different light intensities were different. The highest expression level of HrHEMA and HrPOR was detected under 10,000 lx of illumination, while HrCAO showed the highest expression level under 6000 lx. Lastly, the 5-ALA, Chla, Cxc, SOD, POD, MDA, Pro and CAT contents in different transgenic tobaccos changed significantly under different light intensities.

CONCLUSION

The overexpression of these three genes in tobacco enhanced photosynthesis by accumulating chlorophyll content, but the influential level varied under different light intensities. Furthermore, HrHEMA-, HrPOR- and HrCAO- overexpressing in tobacco can enhance the antioxidant capacity of plants to cope with stress under higher light intensity. However, under lower light intensity, the antioxidant capacity was declined in HrHEMA-, HrPOR- and HrCAO- overexpressing tobaccos.

Phylogenomics With Hyb-Seq Unravels Korean Hosta Evolution

The genus Hosta (Agavoideae and Asparagaceae) is one of the most popular landscaping and ornamental plants native to temperate East Asia. Their popularity has led to extensive hybridization to develop various cultivars. However, their long history of hybridization, cultivation, and selection has brought about taxonomic confusion in the Hosta species delimitation along with their indistinguishable morphology. Here, we conducted the first broad phylogenetic analyses of Hosta species based on the most comprehensive genomic data set to date. To do so, we captured 246 nuclear gene sequences and plastomes from 55 accessions of Korean Hosta species using the Hyb-Seq method. As a result, this study provides the following novel and significant findings: (1) phylogenetic analyses of the captured sequences retrieved six species of Hosta in South Korea compared to five to eleven species based on the previous studies, (2) their phylogenetic relationships suggested that the large genome size was ancestral and the diversification of Korean Hosta species was accompanied by decreases in genome sizes, (3) comparison between nuclear genes and plastome revealed several introgressive hybridization events between Hosta species, and (4) divergence times estimated here showed that Hosta diverged 35.59 million years ago, while Korean Hosta species rapidly diversified during the late Miocene. Last, we explored whether these genomic data could be used to infer the origin of cultivars. In summary, this study provides the most comprehensive genomic resources to be used in phylogenetic, population, and conservation studies of Hosta, as well as for unraveling the origin of many cultivars.

Comparative Analysis of Complete Chloroplast Genome Sequences of Wild and Cultivated Bougainvillea (Nyctaginaceae)

Bougainvillea (Nyctaginaceae) is a popular ornamental plant group primarily grown for its striking colorful bracts. However, despite its established horticultural value, limited genomic resources and molecular studies have been reported for this genus. Thus, to address this existing gap, complete chloroplast genomes of four species (Bougainvillea glabra, Bougainvillea peruviana, Bougainvillea pachyphylla, Bougainvillea praecox) and one Bougainvillea cultivar were sequenced and characterized. The Bougainvillea cp genomes range from 153,966 bp to 154,541 bp in length, comprising a large single-copy region (85,159 bp-85,708 bp) and a small single-copy region (18,014 bp-18,078 bp) separated by a pair of inverted repeats (25,377-25,427 bp). All sequenced plastomes have 131 annotated genes, including 86 protein-coding, eight rRNA, and 37 tRNA genes. These five newly sequenced Bougainvillea cp genomes were compared to the Bougainvillea spectabilis cp genome deposited in GeBank. The results showed that all cp genomes have highly similar structures, contents, and organization. They all exhibit quadripartite structures and all have the same numbers of genes and introns. Codon usage, RNA editing sites, and repeat analyses also revealed highly similar results for the six cp genomes. The amino acid leucine has the highest proportion and almost all favored synonymous codons have either an A or U ending. Likewise, out of the 42 predicted RNA sites, most conversions were from serine (S) to leucine (L). The majority of the simple sequence repeats detected were A/T mononucleotides, making the cp genomes A/T-rich. The contractions and expansions of the IR boundaries were very minimal as well, hence contributing very little to the differences in genome size. In addition, sequence variation analyses showed that Bougainvillea cp genomes share nearly identical genomic profiles though several potential barcodes, such as ycf1, ndhF, and rpoA were identified. Higher variation was observed in both B. peruviana and B. pachyphylla cp sequences based on SNPs and indels analysis. Phylogenetic reconstructions further showed that these two species appear to be the basal taxa of Bougainvillea. The rarely cultivated and wild species of Bougainvillea (B. pachyphylla, B. peruviana, B. praecox) diverged earlier than the commonly cultivated species and cultivar (B. spectabilis, B. glabra, B. cv.). Overall, the results of this study provide additional genetic resources that can aid in further phylogenetic and evolutionary studies in Bougainvillea. Moreover, genetic information from this study is potentially useful in identifying Bougainvillea species and cultivars, which is essential for both taxonomic and plant breeding studies.

The complete chloroplast genome of agave hybrid 11648

Agave hybrid 11648 is the most widely cultivated agave variety for sisal fiber production around the world. In the present study, we have successfully sequenced the chloroplast genome of A. H11648. The complete chloroplast genome size is 157,274 bp in length with a GC content of 37.8%. The genome contains a large single copy region (LSC) of 85,896 bp, a small single copy region (SSC) of 18,230 bp, and a pair of inverted repeat regions (IR) of 26,574 bp. 121 genes are annotated in the chloroplast genome. The numbers of protein-coding, tRNA and rRNA genes are 99, 40 and 8, respectively. Phylogenetic tree reveals that A. H11648 is closely related to A. americana.

Comparative analysis of plastid genomes within the Campanulaceae and phylogenetic implications

The conflicts exist between the phylogeny of Campanulaceae based on nuclear ITS sequence and plastid markers, particularly in the subdivision of Cyanantheae (Campanulaceae). Besides, various and complicated plastid genome structures can be found in species of the Campanulaceae. However, limited availability of genomic information largely hinders the studies of molecular evolution and phylogeny of Campanulaceae. We reported the complete plastid genomes of three Cyanantheae species, compared them to eight published Campanulaceae plastomes, and shed light on a deeper understanding of the applicability of plastomes. We found that there were obvious differences among gene order, GC content, gene compositions and IR junctions of LSC/IRa. Almost all protein-coding genes and amino acid sequences showed obvious codon preferences. We identified 14 genes with highly positively selected sites and branch-site model displayed 96 sites under potentially positive selection on the three lineages of phylogenetic tree. Phylogenetic analyses showed that Cyananthus was more closely related to Codonopsis compared with Cyclocodon and also clearly illustrated the relationship among the Cyanantheae species. We also found six coding regions having high nucleotide divergence value. Hotpot regions were considered to be useful molecular markers for resolving phylogenetic relationships and species authentication of Campanulaceae.

The First Glimpse of Streptocarpus ionanthus (Gesneriaceae) Phylogenomics: Analysis of Five Subspecies' Chloroplast Genomes

Streptocarpus ionanthus (Gesneriaceae) comprise nine herbaceous subspecies, endemic to Kenya and Tanzania. The evolution of Str. ionanthus is perceived as complex due to morphological heterogeneity and unresolved phylogenetic relationships. Our study seeks to understand the molecular variation within Str. ionanthus using a phylogenomic approach. We sequence the chloroplast genomes of five subspecies of Str. ionanthus, compare their structural features and identify divergent regions. The five genomes are identical, with a conserved structure, a narrow size range (170 base pairs (bp)) and 115 unique genes (80 protein-coding, 31 tRNAs and 4 rRNAs). Genome alignment exhibits high synteny while the number of Simple Sequence Repeats (SSRs) are observed to be low (varying from 37 to 41), indicating high similarity. We identify ten divergent regions, including five variable regions (psbM, rps3, atpF-atpH, psbC-psbZ and psaA-ycf3) and five genes with a high number of polymorphic sites (rps16, rpoC2, rpoB, ycf1 and ndhA) which could be investigated further for phylogenetic utility in Str. ionanthus. Phylogenomic analyses here exhibit low polymorphism within Str. ionanthus and poor phylogenetic separation, which might be attributed to recent divergence. The complete chloroplast genome sequence data concerning the five subspecies provides genomic resources which can be expanded for future elucidation of Str. ionanthus phylogenetic relationships.