The Complete Chloroplast Genome Sequences of 14 Curcuma Species: Insights Into Genome Evolution and Phylogenetic Relationships Within Zingiberales

The chloroplast genome of Chrozophora sabulosa Kar. & Kir. and its exploration in the evolutionary position uncertainty of genus Chrozophora

Chrozophora sabulosa Kar. & Kir. is a biennial herbaceous plant that belongs to the Euphorbiaceae family and has medicinal properties. This research aimed to identify the genetic characteristics and phylogenetic position of the Chrozophora genus within the Euphorbiaceae family. The evolutionary position of the Chrozophora genus was previously unknown due to insufficient research. Therefore, to determine the evolutionary link between C. sabulosa and other related species, we conducted a study using the NGS Illumina platform to sequence the C. sabulosa chloroplast (cp.) genome. The study results showed that the genome was 156,488 bp in length. It had a quadripartite structure consisting of two inverted repeats (IRb and IRa) of 24,649-bp, separated by an 87,696-bp LSC region and a 19,494-bp SSC region. The CP genome contained 113 unique genes, including four rRNA genes, 30 tRNA genes, and 79 CDS genes. In the second copy of the inverted repeat, there were 18 duplicated genes. The C. sabulosa lacks the petD, petB, rpl2, and rps16 intron. The analysis of simple sequence repeats (SSRs) revealed 93 SSR loci of 22 types and 78 oligonucleotide repeats of four kinds. The phylogenetic investigation showed that the Chrozophora genus evolved paraphyletically from other members of the Euphorbiaceae family. To support the phylogenetic findings, we selected species from the Euphorbiaceae and Phyllanthaceae families to compare with C. sabulosa for Ks and Ka substitution rates, InDels investigation, IR contraction and expansion, and SNPs analysis. The results of these comparative studies align with the phylogenetic findings. We identified six highly polymorphic regions shared by both families, which could be used as molecular identifiers for the Chrozophora genus (rpl33-rps18, rps18-rpl20, rps15-ycf1, ndhG-ndhI, psaI-ycf4, petA-psbJ). The cp. genome sequence of C. sabulosa reveals the evolution of plastid sequences in Chrozophora species. This is the first time the cp. genome of a Chrozophora genus has been sequenced, serving as a foundation for future sequencing of other species within the Chrozophoreae tribe and facilitating in-depth taxonomic research. The results of this research will also aid in identifying new Chrozophora species.

Characterization and phylogenetic analysis of the complete chloroplast genome of Curcuma comosa and C. latifolia

Members of the Curcuma genus, a crop in the Zingiberaceae, are widely utilized rhizomatous herbs globally. There are two distinct species, C. comosa Roxb. and C. latifolia Roscoe, referred to the same vernacular name "Wan Chak Motluk" in Thai. C. comosa holds economic importance and is extensively used as a Thai traditional medicine due to its phytoestrogenic properties. However, its morphology closely resembles that of C. latifolia, which contains zederone, a compound known for its hepatotoxic effects. They are often confused, which may affect the quality, efficacy and safety of the derived herbal materials. Thus, DNA markers were developed for discriminating C. comosa from C. latifolia. This study focused on analyzing core DNA barcode regions, including rbcL, matK, psbA-trnH spacer and ITS2, of the authentic C. comosa and C. latifolia species. As a result, no variable nucleotides in core DNA barcode regions were observed. The complete chloroplast (cp) genome was introduced to differentiate between the two species. The comparison revealed that the cp genomes of C. comosa and C. latifolia were 162,272 and 162,289 bp, respectively, with a total of 133 identified genes. The phylogenetic analysis revealed that C. comosa and C. latifolia exhibited a very close relationship with other Curcuma species. The cp genome of C. comosa and C. latifolia were identified for the first time, providing valuable insights for species identification and evolutionary research within the Zingiberaceae family.

Elucidating phylogenetic relationships within the genus Curcuma through the comprehensive analysis of the chloroplast genome of Curcuma viridiflora Roxb. 1810 (Zingiberaceae)

Curcuma viridiflora Roxb., a plant species of significant pharmaceutical interest, has been the subject of limited chloroplast genomic research. In this study, we present the sequencing and assembly of the C. viridiflora chloroplast genome, which is characterized by a circular chromosome spanning 162,212 base pairs and a GC content of 36.20%. The genome encodes 87 protein-coding genes (PCGs), 38 transfer RNA (tRNA) genes, and eight ribosomal RNA (rRNA) genes. A phylogenetic analysis was conducted, incorporating eight related species, and based on the complete chloroplast genome and protein-coding DNA sequences of six related taxa within the genus. Outgroup species Zingiber zerumbet and Zingiber officinale were also included in the analysis. The results indicate a close relationship between C. viridiflora and Curcuma phaeocaulis, Curcuma sichuanensis, and Curcuma yunnanensis. This study provides the first chloroplast genome of C. viridiflora, thereby contributing a valuable genomic resource for future research on medicinal plants within the Curcuma genus.

Comparative Analysis of Chloroplast Pan-Genomes and Transcriptomics Reveals Cold Adaptation in Medicago sativa

Alfalfa (Medicago sativa) is a perennial forage legume that is widely distributed all over the world; therefore, it has an extremely complex genetic background. Though population structure and phylogenetic studies have been conducted on a large group of alfalfa nuclear genomes, information about the chloroplast genomes is still lacking. Chloroplast genomes are generally considered to be conservative and play an important role in population diversity analysis and species adaptation in plants. Here, 231 complete alfalfa chloroplast genomes were successfully assembled from 359 alfalfa resequencing data, on the basis of which the alfalfa chloroplast pan-genome was constructed. We investigated the genetic variations of the alfalfa chloroplast genome through comparative genomic, genetic diversity, phylogenetic, population genetic structure, and haplotype analysis. Meanwhile, the expression of alfalfa chloroplast genes under cold stress was explored through transcriptome analysis. As a result, chloroplast genomes of 231 alfalfa lack an IR region, and the size of the chloroplast genome ranges from 125,192 bp to 126,105 bp. Using population structure, haplotypes, and construction of a phylogenetic tree, it was found that alfalfa populations could be divided into four groups, and multiple highly variable regions were found in the alfalfa chloroplast genome. Transcriptome analysis showed that tRNA genes were significantly up-regulated in the cold-sensitive varieties, while rps7, rpl32, and ndhB were down-regulated, and the editing efficiency of ycf1, ycf2, and ndhF was decreased in the cold-tolerant varieties, which may be due to the fact that chloroplasts store nutrients through photosynthesis to resist cold. The huge number of genetic variants in this study provide powerful resources for molecular markers.

Phylogenetics Study to Compare Chloroplast Genomes in Four Magnoliaceae Species

Magnoliaceae, a family of perennial woody plants, contains several endangered species whose taxonomic status remains ambiguous. The study of chloroplast genome information can help in the protection of Magnoliaceae plants and confirmation of their phylogenetic relationships. In this study, the chloroplast genomes were sequenced, assembled, and annotated in Woonyoungia septentrionalis and three Michelia species (Michelia champaca, Michelia figo, and Michelia macclurei). Comparative analyses of genomic characteristics, repetitive sequences, and sequence differences were performed among the four Magnoliaceae plants, and phylogenetic relationships were constructed with twenty different magnolia species. The length of the chloroplast genomes varied among the four studied species ranging from 159,838 bp (Woonyoungia septentrionalis) to 160,127 bp (Michelia macclurei). Four distinct hotspot regions were identified based on nucleotide polymorphism analysis. They were petA-psbJ, psbJ-psbE, ndhD-ndhE, and rps15-ycf1. These gene fragments may be developed and utilized as new molecular marker primers. By using Liriodendron tulipifera and Liriodendron chinense as outgroups reference, a phylogenetic tree of the four Magnoliaceae species and eighteen other Magnoliaceae species was constructed with the method of Shared Coding Sequences (CDS). Results showed that the endangered species, W. septentrionalis, is relatively genetically distinct from the other three species, indicating the different phylogenetic processes among Magnoliaceae plants. Therefore, further genetic information is required to determine the relationships within Magnoliaceae. Overall, complete chloroplast genome sequences for four Magnoliaceae species reported in this paper have shed more light on phylogenetic relationships within the botanical group.

The complete chloroplast genome of Hedychium flavum Roxb., an ornamental, edible and medicinal plant

Hedychium flavum Roxb. 1820 is a perennial herb mainly distributed in China, India, Myanmar and Thailand with ornamental, edible and medicinal value. It is extensively cultivated as a source of aromatic essential oils, ornamental plant, food flavorings and vegetables, and folk medicine. In this study, we sequence the complete chloroplast genome of H. flavum by de novo assembly. The assembled genome has a typical quadripartite circular structure with 163,909 bp in length, containing a large single-copy region (LSC, 88,589 bp), a small single-copy region (SSC, 15,762 bp), and two inverted repeat regions (IRs, 29,779 bp). The cp genome contains 133 genes, including 87 protein-coding genes, 38 tRNA genes and 8 rRNA genes. Phylogenetic analysis based on the complete cp genome shows a close affinity of H. flavum and H. neocarneum with 100% bootstrap support. This study will provide useful genetic resource for further phylogenetic analysis of the genus Hedychium and Zingiberaceae.

Comparative Chloroplast Genomics of 21 Species in Zingiberales with Implications for Their Phylogenetic Relationships and Molecular Dating

Zingiberales includes eight families and more than 2600 species, with many species having important economic and ecological value. However, the backbone phylogenetic relationships of Zingiberales still remain controversial, as demonstrated in previous studies, and molecular dating based on chloroplast genomes has not been comprehensively studied for the whole order. Herein, 22 complete chloroplast genomes from 21 species in Zingiberales were sequenced, assembled, and analyzed. These 22 genomes displayed typical quadripartite structures, which ranged from 161,303 bp to 163,979 bp in length and contained 111-112 different genes. The genome structures, gene contents, simple sequence repeats, long repeats, and codon usage were highly conserved, with slight differences among these genomes. Further comparative analysis of the 111 complete chloroplast genomes of Zingiberales, including 22 newly sequenced ones and the remaining ones from the national center for biotechnology information (NCBI) database, identified three highly divergent regions comprising ccsA, psaC, and psaC-ndhE. Maximum likelihood and Bayesian inference phylogenetic analyses based on chloroplast genome sequences found identical topological structures and identified a strongly supported backbone of phylogenetic relationships. Cannaceae was sister to Marantaceae, forming a clade that was collectively sister to the clade of (Costaceae, Zingiberaceae) with strong support (bootstrap (BS) = 100%, and posterior probability (PP) = 0.99-1.0); Heliconiaceae was sister to the clade of (Lowiaceae, Strelitziaceae), then collectively sister to Musaceae with strong support (BS = 94-100%, and PP = 0.93-1.0); the clade of ((Cannaceae, Marantaceae), (Costaceae, Zingiberaceae)) was sister to the clade of (Musaceae, (Heliconiaceae, (Lowiaceae, Strelitziaceae))) with robust support (BS = 100%, and PP = 1.0). The results of divergence time estimation of Zingiberales indicated that the crown node of Zingiberales occurred approximately 85.0 Mya (95% highest posterior density (HPD) = 81.6-89.3 million years ago (Mya)), with major family-level lineages becoming from 46.8 to 80.5 Mya. These findings proved that chloroplast genomes could contribute to the study of phylogenetic relationships and molecular dating in Zingiberales, as well as provide potential molecular markers for further taxonomic and phylogenetic studies of Zingiberales.

Comparative chloroplast genome analysis of four Polygonatum species insights into DNA barcoding, evolution, and phylogeny

The Polygonatum genus represents a perennial herb with the Liliaceae family, boasting substantial economic and medicinal significance. The majority of Polygonatum plants exhibit notable similarity while lacking distinctive identifying characteristics, thus resulting in the proliferation of adulterated medicinal materials within the market. Within this study, we conducted an in-depth analysis of the complete chloroplast (cp) genomes of four Polygonatum plants and compared them with four closely akin species. The primary objectives were to unveil structural variations, species divergence, and the phylogenetic interrelations among taxa. The cp genomes of the four Polygonatum species were typified by a conventional quadripartite structure, incorporating a large single copy region (LSC), a small single copy region (SSC), and a pair of inverted repeat regions. In total, we annotated a range of 131 to 133 genes, encompassing 84 to 86 protein-coding genes, 38 transfer RNA (tRNA) genes, 8 ribosomal RNA (rRNA) genes, and 0 to 2 pseudogenes (ycf1, infA). Our comparative analyses unequivocally revealed a remarkable consistency in gene order and GC content within the Polygonatum genus. Furthermore, we predicted a potential 59 to 64 RNA editing sites distributed across 22 protein-coding genes, with the ndhB gene exhibiting the most prominent propensity for RNA editing sites, boasting a tally of 15 sites. Notably, six regions of substantial potential variability were ascertained, characterized by elevated Pi values. Noteworthy, molecular markers for species identification, population genetic scrutiny, and phylogenetic investigations within the genus were identified in the form of the psaJ-rpl33 and trnS + trnT-psaD barcodes. The resultant phylogenetic tree unequivocally depicted the formation of a monophyletic clade comprising species within the evolutionary framework of Liliaceae, demonstrating closer evolutionary affinities with Maianthemum, Dracaeneae, and Asparageae. This comprehensive compendium of findings collectively contributes to the advancement of molecular species identification, elucidation of phylogenetic interrelationships, and the establishment of DNA barcodes tailored to the Polygonatum species.

Plastid phylogenomics and species discrimination in the "Chinese" clade of Roscoea (Zingiberaceae)

Roscoea is an alpine or subalpine genus from the pan-tropical family Zingiberaceae, which consists of two disjunct groups in geography, namely the "Chinese" clade and the "Himalayan" clade. Despite extensive research on the genus, Roscoea species remain poorly defined and relationships between these species are not well resolved. In this study, we used plastid genomes of nine species and one variety to resolve phylogenetic relationships within the "Chinese" clade of Roscoea and as DNA super barcodes for species discrimination. We found that Roscoea plastid genomes ranged in length from 163,063 to 163,796 bp, and encoded 113 genes, including 79 protein-coding genes, 30 tRNA genes, four rRNA genes. In addition, expansion and contraction of the IR regions showed obvious infraspecific conservatism and interspecific differentiation. Plastid phylogenomics revealed that species belonging to the "Chinese" clade of Roscoea can be divided into four distinct subclades. Furthermore, our analysis supported the independence of R. cautleoides var. pubescens, the recovery of Roscoea pubescens Z.Y. Zhu, and a close relationship between R. humeana and R. cautloides. When we used the plastid genome as a super barcode, we found that it possessed strong discriminatory power (90%) with high support values. Intergenic regions provided similar resolution, which was much better than that of protein-coding regions, hypervariable regions, and DNA universal barcodes. However, plastid genomes could not completely resolve Roscoea phylogeny or definitively discriminate species. These limitations are likely related to the complex history of Roscoea speciation, poorly defined species within the genus, and the maternal inheritance of plastid genomes.

Comparative Plastomes of Curcuma alismatifolia (Zingiberaceae) Reveal Diversified Patterns among 56 Different Cut-Flower Cultivars

Curcuma alismatifolia (Zingiberaceae) is an ornamental species with high economic value due to its recent rise in popularity among floriculturists. Cultivars within this species have mixed genetic backgrounds from multiple hybridization events and can be difficult to distinguish via morphological and histological methods alone. Given the need to improve identification resources, we carried out the first systematic study using plastomic data wherein genomic evolution and phylogenetic relationships from 56 accessions of C. alismatifolia were analyzed. The newly assembled plastomes were highly conserved and ranged from 162,139 bp to 164,111 bp, including 79 genes that code for proteins, 30 tRNA genes, and 4 rRNA genes. The A/T motif was the most common of SSRs in the assembled genomes. The Ka/Ks values of most genes were less than 1, and only two genes had Ka/Ks values above 1, which were rps15 (1.15), and ndhl (1.13) with petA equal to 1. The sequence divergence between different varieties of C. alismatifolia was large, and the percentage of variation in coding regions was lower than that in the non-coding regions. Such data will improve cultivar identification, marker assisted breeding, and preservation of germplasm resources.

Phylogenomics and plastome evolution of Indigofera (Fabaceae)

Introduction

Indigofera L. is the third largest genus in Fabaceae and includes economically important species that are used for indigo dye-producing, medicinal, ornamental, and soil and water conservation. The genus is taxonomically difficult due to the high level of overlap in morphological characters of interspecies, fewer reliability states for classification, and extensive adaptive evolution. Previous characteristic-based taxonomy and nuclear ITS-based phylogenies have contributed to our understanding of Indigofera taxonomy and evolution. However, the lack of chloroplast genomic resources limits our comprehensive understanding of the phylogenetic relationships and evolutionary processes of Indigofera.

Methods

Here, we newly assembled 18 chloroplast genomes of Indigofera. We performed a series of analyses of genome structure, nucleotide diversity, phylogenetic analysis, species pairwise Ka/Ks ratios, and positive selection analysis by combining with allied species in Papilionoideae.

Results and discussion

The chloroplast genomes of Indigofera exhibited highly conserved structures and ranged in size from 157,918 to 160,040 bp, containing 83 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Thirteen highly variable regions were identified, of which trnK-rbcL, ndhF-trnL, and ycf1 were considered as candidate DNA barcodes for species identification of Indigofera. Phylogenetic analysis using maximum likelihood (ML) and Bayesian inference (BI) methods based on complete chloroplast genome and protein-coding genes (PCGs) generated a well-resolved phylogeny of Indigofera and allied species. Indigofera monophyly was strongly supported, and four monophyletic lineages (i.e., the Pantropical, East Asian, Tethyan, and Palaeotropical clades) were resolved within the genus. The species pairwise Ka/Ks ratios showed values lower than 1, and 13 genes with significant posterior probabilities for codon sites were identified in the positive selection analysis using the branch-site model, eight of which were associated with photosynthesis. Positive selection of accD suggested that Indigofera species have experienced adaptive evolution to selection pressures imposed by their herbivores and pathogens. Our study provided insight into the structural variation of chloroplast genomes, phylogenetic relationships, and adaptive evolution in Indigofera. These results will facilitate future studies on species identification, interspecific and intraspecific delimitation, adaptive evolution, and the phylogenetic relationships of the genus Indigofera.

Widely Targeted Volatilomics and Metabolomics Analysis Reveal the Metabolic Composition and Diversity of Zingiberaceae Plants

Zingiberaceae plants are widely used in the food and pharmaceutical industries; however, research on the chemical composition and interspecific differences in the metabolome and volatilome of Zingiberaceae plants is still limited. In this study, seven species of Zingiberaceae plants were selected, including Curcuma longa L., Zingiber officinale Rosc., Alpinia officinarum Hance, Alpinia tonkinensis Gagnep, Amomum tsaoko Crevost et Lemarie, Alpinia hainanensis K. Schum. and Amomum villosum Lour. Myristica fragrans Houtt. was also selected due to its flavor being similar to that of the Zingiberaceae plant. The metabolome and volatilome of selected plants were profiled by widely targeted approaches; 542 volatiles and 738 non-volatile metabolites were detected, and β-myrcene, α-phellandrene and α-cadinene were detected in all the selected plants, while chamigren, thymol, perilla, acetocinnamone and cis-α-bisabolene were exclusively detected in certain Zingiberaceae plants. Differential analysis showed that some terpenoids, such as cadalene, cadalene-1,3,5-triene, cadalene-1,3,8-triene and (E)-β-farnesene, and some lipids, including palmitic acid, linoleic acid and oleic acid were amongst the most varied compounds in Zingiberaceae plants. In conclusion, this study provided comprehensive metabolome and volatilome profiles for Zingiberaceae plants and revealed the metabolic differences between these plants. The results of this study could be used as a guide for the nutrition and flavor improvement of Zingiberaceae plants.

Cross-Compatibility in Interspecific Hybridization of Different Curcuma Accessions

Curcuma is extensively cultivated as a medicinal and ornamental plant in tropical and subtropical regions. Due to the bright bract color, distinctive inflorescence and long blooming period, it has become a new favorite in terms of the urban landscape, potted flowers and cut flowers. However, little research on breeding new cultivars using traditional plant breeding methods is available on the genus Curcuma. In the present study, pollen viability and stigma receptivity evaluation were performed, and the genetic relationship of 38 Curcuma accessions was evaluated, then 5 C. alismatifolia Gagnep. (Ca), 2 C. hybrid (Ch), 2 C. sparganiifolia Gagnep. cultivars and 4 Curcuma native species were selected as parents for subsequent interspecific cross-breeding. A total of 132 reciprocal crosses were carried out for interspecific hybridization, including 70 obverse and 62 inverse crosses. Obvious discrepancies among fruit-setting rates were manifested in different combinations and in reciprocal crosses. Results showed that the highest fruit-setting rate (87.5%) was observed in the Ca combinations. There were 87 combinations with a fruit-setting rate of 0%, which meant nearly 65.9% was incompatible. We concluded that C. alismatifolia 'Siam Shadow' (Ch34) was suitable as a male parent and C. petiolata Roxb. (Cpet) was suitable as a female parent to improve the fruit-setting rates. The maximum number of seeds per fruit (45.4) was obtained when C. alismatifolia 'Chiang Mai Pink' (Ca01) was used as a female parent followed by C. attenuata Wall. ex Baker (Catt) (42.8) and C. alismatifolia 'Splash' (Ca63) (39.6) as male parents. The highest germination rate was observed for the Ca group followed by Catt and C. sparganiifolia 'Maetang Sunrise' (Csms). The germination rates of Ca accessions ranged from 58.2% (C. alismatifolia 'Siam Scarlet' (Ca06) as a male parent) to 89.3% (C. alismatifolia 'Sitone' (Ca10) as a male parent) with an average value of 74.0%. Based on the results of hybrid identification, all the individuals from the four combinations exhibited paternal-specific bands, indicating that the true hybrid rates of crossings were 100%. Our results would facilitate the interspecific hybridization and introduction of genetic variation from wild species into the cultivars in Curcuma in the future, which could be helpful in realizing the sustainable application in urban green areas.

The population genomic analyses of chloroplast genomes shed new insights on the complicated ploidy and evolutionary history in Fragaria

The genus Fragaria consists of a rich diversity of ploidy levels with diploid (2x), tetraploid (4x), pentaploid (5x), hexaploidy (6x), octoploid (8x) and decaploid (10x) species. Only a few studies have explored the origin of diploid and octoploid strawberry, and little is known about the roles of tetraploidy and hexaploidy during the evolution of octoploid strawberry. The chloroplast genome is usually a stable circular genome and is widely used in investigating the evolution and matrilineal identification. Here, we assembled the chloroplast genomes of F. x ananassa cv. 'Benihoppe' (8x) using Illumina and HiFi data seperately. The genome alignment results showed that more InDels were located in the chloroplast genomes based on the PacBio HiFi data than Illumina data. We obtain highly accurate chloroplast genomes assembled through GetOrganelle using Illumina reads. We assembled 200 chloroplast genomes including 198 Fragaria (21 species) and 2 Potentilla samples. Analyses of sequence variation, phylogenetic and PCA analyses showed that Fragaria was divided into five groups. F. iinumae, F. nilgerrensis and all octoploid accessions formed Group A, C and E separately. Species native to western China were clustered into Group B. Group D consisted of F. virdis, F. orientalis, F. moschata, and F. vesca. STRUCTURE and haplotype network confirmed that the diploid F. vesca subsp. bracteata was the last maternal donator of octoploid strawberry. The dN/dS ratio estimated for the protein-coding genes revealed that genes involved in ATP synthase and photosystem function were under positive selection. These findings demonstrate the phylogeny of totally 21 Fragaria species and the origin of octoploid species. F. vesca was the last female donator of octoploid, which confirms the hypothesis that the hexaploid species F. moschata may be an evolutionary intermediate between the diploids and wild octoploid species.

Chloroplast genomes of four Carex species: Long repetitive sequences trigger dramatic changes in chloroplast genome structure

The chloroplast genomes of angiosperms usually have a stable circular quadripartite structure that exhibits high consistency in genome size and gene order. As one of the most diverse genera of angiosperms, Carex is of great value for the study of evolutionary relationships and speciation within its genus, but the study of the structure of its chloroplast genome is limited due to its highly expanded and restructured genome with a large number of repeats. In this study, we provided a more detailed account of the chloroplast genomes of Carex using a hybrid assembly of second- and third-generation sequencing and examined structural variation within this genus. The study revealed that chloroplast genomes of four Carex species are significantly longer than that of most angiosperms and are characterized by high sequence rearrangement rates, low GC content and gene density, and increased repetitive sequences. The location of chloroplast genome structural variation in the species of Carex studied is closely related to the positions of long repeat sequences; this genus provides a typical example of chloroplast structural variation and expansion caused by long repeats. Phylogenetic relationships constructed based on the chloroplast protein-coding genes support the latest taxonomic system of Carex, while revealing that structural variation in the chloroplast genome of Carex may have some phylogenetic significance. Moreover, this study demonstrated a hybrid assembly approach based on long and short reads to analyze complex chloroplast genome assembly and also provided an important reference for the analysis of structural rearrangements of chloroplast genomes in other taxa.

Complete chloroplast genomes provide insights into evolution and phylogeny of Zingiber (Zingiberaceae)

BACKGROUND

The genus Zingiber of the Zingiberaceae is distributed in tropical, subtropical, and in Far East Asia. This genus contains about 100-150 species, with many species valued as important agricultural, medicinal and horticultural resources. However, genomic resources and suitable molecular markers for species identification are currently sparse.

RESULTS

We conducted comparative genomics and phylogenetic analyses on Zingiber species. The Zingiber chloroplast genome (size range 162,507-163,711 bp) possess typical quadripartite structures that consist of a large single copy (LSC, 86,986-88,200 bp), a small single copy (SSC, 15,498-15,891 bp) and a pair of inverted repeats (IRs, 29,765-29,934 bp). The genomes contain 113 unique genes, including 79 protein coding genes, 30 tRNA and 4 rRNA genes. The genome structures, gene contents, amino acid frequencies, codon usage patterns, RNA editing sites, simple sequence repeats and long repeats are conservative in the genomes of Zingiber. The analysis of sequence divergence indicates that the following genes undergo positive selection (ccsA, ndhA, ndhB, petD, psbA, psbB, psbC, rbcL, rpl12, rpl20, rpl23, rpl33, rpoC2, rps7, rps12 and ycf3). Eight highly variable regions are identified including seven intergenic regions (petA-pabJ, rbcL-accD, rpl32-trnL-UAG, rps16-trnQ-UUG, trnC-GCA-psbM, psbC-trnS-UGA and ndhF-rpl32) and one genic regions (ycf1). The phylogenetic analysis revealed that the sect. Zingiber was sister to sect. Cryptanthium rather than sect. Pleuranthesis.

CONCLUSIONS

This study reports 14 complete chloroplast genomes of Zingiber species. Overall, this study provided a solid backbone phylogeny of Zingiber. The polymorphisms we have uncovered in the sequencing of the genome offer a rare possibility (for Zingiber) of the generation of DNA markers. These results provide a foundation for future studies that seek to understand the molecular evolutionary dynamics or individual population variation in the genus Zingiber.

Characterization of the complete plastid genome of Korean endemic, Ajuga spectabilis Nakai (Lamiaceae)

Ajuga spectabilis Nakai is a Korean endemic species in Lamiaceae. In spite of its importance, genomic studies are not performed on this species. Here, we report the complete plastid genome sequences of A. spectabilis, which will provide valuable information for its natural conservation and future studies for the plastid genome evolution. The plastid genome is 150,417 bp in length, containing a large single-copy region (LSC) of 82,140 bp and a small single-copy (SSC) region of 17,165 bp which are separated by a pair of inverted repeats (IR) of 25,556 bp. It encodes 113 genes, including 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. The overall GC content is 38.3%, and those in the LSC, SSC, and IR regions are 36.4%, 32.2%, and 43.3%, respectively, which is consistent with other Ajuga species. Our phylogenetic analysis revealed that A. spectabilis formed a close relationship with A. ciliata and A. decumbens.

Complete chloroplast genome sequences of Myristicaceae species with the comparative chloroplast genomics and phylogenetic relationships among them

BACKGROUND

Myristicaceae was widly distributed from tropical Asia to Oceania, Africa, and tropical America. There are 3 genera and 10 species of Myristicaceae present in China, mainly distributed in the south of Yunnan Province. Most research on this family focuses on fatty acids, medicine, and morphology. Based on the morphology, fatty acid chemotaxonomy, and a few of molecular data, the phylogenetic position of Horsfieldia pandurifolia Hu was controversial.

RESULTS

In this study, the chloroplast genomes of two Knema species, Knema globularia (Lam.) Warb. and Knema cinerea (Poir.) Warb., were characterized. Comparing the genome structure of these two species with those of other eight published species, including three Horsfieldia species, four Knema species, and one Myristica species, it was found that the chloroplast genomes of these species were relatively conserved, retaining the same gene order. Through sequence divergence analysis, there were 11 genes and 18 intergenic spacers were subject to positive selection, which can be used to analyze the population genetic structure of this family. Phylogenetic analysis showed that all Knema species were clustered in the same group and formed a sister clade with Myristica species support by both high maximum likelihood bootstrap values and Bayesian posterior probabilities; among Horsfieldia species, Horsfieldia amygdalina (Wall.) Warb., Horsfieldia kingii (Hook.f.) Warb., Horsfieldia hainanensis Merr. and Horsfieldia tetratepala C.Y.Wu. were grouped together, but H. pandurifolia formed a single group and formed a sister clade with genus Myristica and Knema. Through the phylogenetic analysis, we support de Wilde' view that the H. pandurifolia should be separated from Horsfieldia and placed in the genus Endocomia, namely Endocomia macrocoma subsp. prainii (King) W.J.de Wilde.

CONCLUSION

The findings of this study provide a novel genetic resources for future research in Myristicaceae and provide a molecular evidence for the taxonomic classification of Myristicaceae.

Comparative chloroplast genome analyses of Amomum: insights into evolutionary history and species identification

BACKGROUND

Species in genus Amomum always have important medicinal and economic values. Classification of Amomum using morphological characters has long been a challenge because they exhibit high similarity. The main goals of this study were to mine genetic markers from cp genomes for Amomum species identification and discover their evolutionary history through comparative analysis.

RESULTS

Three species Amomum villosum, Amomum maximum and Amomum longipetiolatum were sequenced and annotated for the complete chloroplast (cp) genomes, and the cp genomes of A. longipetiolatum and A. maximum were the first reported. Three cp genomes exhibited typical quadripartite structures with 163,269-163,591 bp in length. Each genome encodes 130 functional genes including 79 protein-coding, 26 tRNAs and 3 rRNAs genes. 113-152 SSRs and 99 long repeats were identified in the three cp genomes. By designing specific primers, we amplified the highly variable loci and the mined genetic marker ccsA exhibited a relatively high species identification resolution in Amomum. The nonsynonymous and synonymous substitution ratios (Ka/Ks) in Amomum and Alpinia showed that most genes were subjected to a purifying selection. Phylogenetic analysis revealed the evolutionary relationships of Amomum and Alpinia species and proved that Amomum is paraphyletic. In addition, the sequenced sample of A. villosum was found to be a hybrid, becoming the first report of natural hybridization of this genus. Meanwhile, the high-throughput sequencing-based ITS2 analysis was proved to be an efficient tool for interspecific hybrid identification and with the help of the chloroplast genome, the hybrid parents can be also be determined.

CONCLUSION

The comparative analysis and mined genetic markers of cp genomes were conducive to species identification and evolutionary relationships of Amomum.

Complete chloroplast genome sequences of Phlomis fruticosa and Phlomoides strigosa and comparative analysis of the genus Phlomis sensu lato (Lamiaceae)

The taxonomic terms "Phlomis" and "Phlomoides" had been used to describe two sections within the genus Phlomis belonging to the Lamiaceae family. Recently, phylogenetic analyses using molecular markers showed that Phlomis and Phlomoides formed two monophyletic clades, and thus they are generally accepted as separate genera. In this study, we assembled the complete chloroplast genome of Phlomis fruticosa, which is the first reported chloroplast genome belonging to Phlomis genus, as well as the complete chloroplast genome of Phlomoides strigosa belonging to Phlomoides genus. The results showed that the length of chloroplast genome was 151,639 bp (Phlomis fruticosa) and 152,432 bp (Phlomoides strigosa), with conserved large single copy regions, small single copy regions, and inverted repeat regions. 121 genes in Phlomis fruticosa and 120 genes in Phlomoides strigosa were annotated. The chloroplast genomes of Phlomis fruticosa, Phlomoides strigosa, and three reported Phlomoides species, as well as those of 51 species from the Lamiaceae family, which covered 12 subfamilies, were subjected to phylogenetic analyses. The Phlomis and Phlomoides species were split into two groups, which were well supported by both maximum likelihood and Bayesian inference tree analyses. Our study provided further evidence to recognize Phlomis and Phlomoides as independent genera.

A chromosome-scale genome assembly of turmeric provides insights into curcumin biosynthesis and tuber formation mechanism

Curcuma longa, known as the 'golden spice' and 'life spice', is one of the most commonly utilized spices in the world and also has medicinal, cosmetic, dye and flavoring values. Herein, we present the chromosomal-level genome for turmeric to explore the differences between tubers and rhizomes in the regulation of curcumin biosynthesis and the mechanism of tuber formation. We assembled the turmeric genome into 21 pseudochromosomes using Pacbio long reads complemented with Hi-C technologies, which has a total length of 1.11 Gb with scaffold N50 of 50.12 Mb and contains 49,612 protein-coding genes. Genomic evolutionary analysis indicated that turmeric and ginger have shared a recent WGD event. Contraction analysis of gene families showed possible roles for transcription factors, phytohormone signaling, and plant-pathogen interactions associated genes in adaptation to harsh environments. Transcriptomic data from tubers at different developmental stages indicated that candidate genes related to phytohormone signaling and carbohydrate metabolic responses may be associated with the induction of tuber formation. The difference in curcumin content between rhizomes and tubers reflected the remodeling of secondary metabolites under environmental stress, which was associated with plant defense in response to abiotic stresses. Overall, the availability of the C. longa genome provides insight into tuber formation and curcumin biosynthesis in turmeric as well as facilitating the understanding of other Curcuma species.

The Complete Chloroplast Genome of Endangered Species Stemona parviflora: Insight into the Phylogenetic Relationship and Conservation Implications

Stemona parviflora is an endangered species, narrowly endemic to Hainan and Southwest Guangdong. The taxonomic classification of S. parviflora remains controversial. Moreover, studying endangered species is helpful for current management and conservation. In this study, the first complete chloroplast genome of S. parviflora was assembled and compared with other Stemona species. The chloroplast genome size of S. parviflora was 154,552 bp, consisting of 87 protein-coding genes, 38 tRNA genes, 8 rRNA genes, and one pseudogene. The ψycf1 gene was lost in the cp genome of S. sessilifolia, but it was detected in four other species of Stemona. The inverted repeats (IR) regions have a relatively lower length variation compared with the large single copy (LSC) and small single copy (SSC) regions. Long repeat sequences and simple sequence repeat (SSR) were detected, and most SSR were distributed in the LSC region. Codon usage bias analyses revealed that the RSCU value of the genus Stemona has almost no difference. As with most angiosperm chloroplast genomes, protein-coding regions were more conservative than the inter-gene spacer. Seven genes (atpI, ccsA, cemA, matK, ndhA, petA, and rpoC1) were detected under positive selection in different Stemona species, which may result from adaptive evolution to different habitats. Phylogenetic analyses show the Stemona cluster in two main groups; S. parviflora were closest to S. tuberosa. A highly suitable region of S. parviflora was simulated by Maxent in this study; it is worth noting that the whole territory of Taiwan has changed to a low fitness area and below in the 2050 s, which may not be suitable for the introduction and cultivation of S. parviflora. In addition, limited by the dispersal capacity of S. parviflora, it is necessary to carry out artificial grafts to expand the survival areas of S. parviflora. Our results provide valuable information on characteristics of the chloroplast genome, phylogenetic relationships, and potential distribution range of the endangered species S. parviflora.

Structural Characterization of the Acer ukurunduense Chloroplast Genome Relative to Related Species in the Acer Genus

Acer ukurunduense refers to a deciduous tree distributed in Northeast Asia and is a widely used landscaping tree species. Although several studies have been conducted on the species’ ecological and economic significance, limited information is available on its phylo-genomics. Our study newly constitutes the complete chloroplast genome of A. ukurunduense into a 156,645-bp circular DNA, which displayed a typical quadripartite structure. In addition, 133 genes were identified, containing 88 protein-coding genes, 37 tRNA genes, and eight rRNA genes. In total, 107 simple sequence repeats and 49 repetitive sequences were observed. Thirty-two codons indicated that biased usages were estimated across 20 protein-coding genes (CDS) in A. ukurunduense. Four hotspot regions (trnK-UUU/rps16, ndhF/rpl32, rpl32/trnL-UAG, and ycf1) were detected among the five analyzed Acer species. Those hotspot regions may be useful molecular markers and contribute to future population genetics studies. The phylogenetic analysis demonstrated that A. ukurunduense is most closely associated with the species of Sect. Palmata. A. ukurunduense and A. pubipetiolatum var. pingpienense diverged in 22.11 Mya. We selected one of the hypervariable regions (trnK-UUU/rps16) to develop a new molecular marker and designed primers and confirmed that the molecular markers could accurately discriminate five Acer species through Sanger sequencing. By sequencing the cp genome of A. ukurunduense and comparing it with the relative species of Acer, we can effectively address the phylogenetic problems of Acer at the species level and provide insights into future research on population genetics and genetic diversity.

Plastid genome evolution of a monophyletic group in the subtribe Lauriineae (Laureae, Lauraceae)

Litsea, a non-monophyletic group of the tribe Laureae (Lauraceae), plays important roles in the tropical and subtropical forests of Asia, Australia, Central and North America, and the islands of the Pacific. However, intergeneric relationships between Litsea and Laurus, Lindera, Parasassafras and Sinosassafras of the tribe Laureae remain unresolved. In this study, we present phylogenetic analyses of seven newly sequenced Litsea plastomes, together with 47 Laureae plastomes obtained from public databases, representing six genera of the Laureae. Our results highlight two highly supported monophyletic groups of Litsea taxa. One is composed of 16 Litsea taxa and two Lindera taxa. The 18 plastomes of these taxa were further compared for their gene structure, codon usage, contraction and expansion of inverted repeats, sequence repeats, divergence hotspots, and gene evolution. The complete plastome size of newly sequenced taxa varied between 152,377 bp (Litsea auriculata) and 154,117 bp (Litsea pierrei). Seven of the 16 Litsea plastomes have a pair of insertions in the IRa (trnL-trnH) and IRb (ycf2) regions. The 18 plastomes of Litsea and Lindera taxa exhibit similar gene features, codon usage, oligonucleotide repeats, and inverted repeat dynamics. The codons with the highest frequency among these taxa favored A/T endings and each of these plastomes had nine divergence hotspots, which are located in the same regions. We also identified six protein coding genes (accD, ndhJ, rbcL, rpoC2, ycf1 and ycf2) under positive selection in Litsea; these genes may play important roles in adaptation of Litsea species to various environments.

Molecular evolution and phylogenetic relationships of Ligusticum (Apiaceae) inferred from the whole plastome sequences

Background

The genus Ligusticum belongs to Apiaceae, and its taxonomy has long been a major difficulty. A robust phylogenetic tree is the basis of accurate taxonomic classification of Ligusticum. We herein used 26 (including 14 newly sequenced) plastome-scale data to generate reliable phylogenetic trees to explore the phylogenetic relationships of Chinese Ligusticum.

Results

We found that these plastid genomes exhibited diverse plastome characteristics across all four currently identified clades in China, while the plastid protein-coding genes were conserved. The phylogenetic analyses by the concatenation and coalescent methods obtained a more robust molecular phylogeny than prior studies and showed the non-monophyly of Chinese Ligusticum. In the concatenation-based phylogeny analyses, the two datasets yielded slightly different topologies that may be primarily due to the discrepancy in the number of variable sites.

Conclusions

Our plastid phylogenomics analyses emphasized that the current circumscription of the Chinese Ligusticum should be reduced, and the taxonomy of Ligusticum urgently needs revision. Wider taxon sampling including the related species of Ligusticum will be necessary to explore the phylogenetic relationships of this genus. Overall, our study provided new insights into the taxonomic classification of Ligusticum and would serve as a framework for future studies on taxonomy and delimitation of Ligusticum from the perspective of the plastid genome.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02010-z.

Comparative Analyses of Chloroplast Genomes Provide Comprehensive Insights into the Adaptive Evolution of Paphiopedilum (Orchidaceae)

An elucidation of how the selection pressures caused by habitat environments affect plant plastid genomes and lead to the adaptive evolution of plants, is a very intense area of research in evolutionary biology. The genus Paphiopedilum is a predominant group of orchids that includes over 66 species with high horticultural and ornamental value. However, owing to the destructive exploitation and habitat deterioration of wild germplasm resources of Paphiopedilum, it needs more molecular genetic resources and studies on this genus. The chloroplast is cytoplasmically inherited and often used in evolutionary studies. Thus, for this study, we newly sequenced, assembled and annotated five chloroplast genomes of the Paphiopedilum species. The size of these genomes ranged from 155,886 bp (P. henryanum) to 160,503 bp (P. ‘GZSLKY’ Youyou) and they contained 121–122 genes, which consisted of 76 protein coding genes, eight ribosomal RNAs, and 37–38 transfer RNAs. Combined with the other 14 Paphiopedilum species, the characteristics of the repeat sequences, divergent hotspot regions, and the condo usage bias were evaluated and identified, respectively. The gene transfer analysis showed that some fragments of the ndh and ycf gene families were shared by both the chloroplast and nucleus. Although the genomic structure and gene content was conserved, there was a significant boundary shift caused by the inverted repeat (IR) expansion and small single copy (SSC) contraction. The lower GC content and loss of ndh genes could be the result of adaptive evolutionary responses to its unique habitats. The genes under positive selection, including accD, matK, psbM, rpl20, rps12, ycf1, and ycf2 might be regarded as potential candidate genes for further study, which significantly contribute to the adaptive evolution of Paphiopedilum.

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.

Evaluation of curcuminoids, physiological adaptation, and growth of Curcuma longa under water deficit and controlled temperature

Turmeric (Curcuma longa L.; Zingiberaceae), an economically important crop and a major spice in Indian cuisine, produces natural yellow color (curcumin) as well as curcuminoids which are widely utilized in traditional and modern medicinal practices. During the turmeric culture, the fluctuations of precipitation and seasonal changes in the whole life cycle play a major role, especially water shortage and decreasing temperature (in winter season), leading to rhizome dormancy under extreme weather conditions. The objective of this investigation was to understand how the water deficit and reduced temperature affect turmeric growth, physiological adaptation, quantity, and quality of turmeric rhizomes. Four-month-old turmeric plants were subjected to four treatments, namely normal temperature and well-watered (RT-WW), or water-deficit (RT-WD) conditions in the greenhouse, 25 °C controlled temperature and well-watered (CT-WW), or water-deficit (CT-WD) conditions in glasshouse. Leaf osmotic potential considerably declined in 30 days CT-WD treatment, leading to chlorophyll degradation by 26.04%, diminution of maximum quantum yield of PSII (F_v/F_m) by 23.50%, photon yield of PSII (Φ_PSII) by 29.01%, and reduction of net photosynthetic rate (P_n) by 89.39% over CT-WW (control). After 30 days water withholding, fresh- and dry-weights of rhizomes of turmeric plants grown under CT-WD declined by 30-50% when compared with RT-WW conditions. Subsequently, curcuminoid content was reduced by 40% over RT-WW plants (control), whereas transcriptional expression levels of curcuminoids-related genes (CURS1, CURS2, CURS3, and DCS) were upregulated in CT-WD conditions. In summary, the water withholding and controlled temperature (constant at 25 °C day/night) negatively affected turmeric plants as abiotic stresses tend to limit overall plant growth performances and curcuminoid yield.

Leaf epidermal micromorphology of Zingiber (Zingiberaceae) from China and its systematic significance

Leaf epidermal characteristics are important for phylogenetic and taxonomic studies of many plants, but there is currently insufficient such data for this application in Zingiber species. Therefore, the leaf epidermal micromorphology of 22 species in three sections of Zingiber was investigated by light microscopy and scanning electron microscopy. Differences between various taxonomic groups of Zingiberaceae were also compared to assess their phylogenetic and taxonomic significance. As in other genera of Zingiberaceae, the epidermal cells in both the adaxial and abaxial epidermis of Zingiber species were found to be hexagonal or polygonal, with non-sinuous anticlinal walls that are arranged parallel to leaf veins. Tetracytic stomata are mostly randomly distributed in the intercostal regions of both surfaces and are more common on the abaxial surface. The stomatal density of the species in sect. Pleuranthesis is significantly lower than that in sects. Zingiber and Cryptanthium. There are two types of trichome in Zingiber: so-called "delicate" trichomes are present in most species, while "stout" trichomes with a swollen base are only found in Z.corallinum and Z.montanum. Oil cells occur in both epidermal layers of some species in sects. Zingiber and Cryptanthium, but only in the abaxial epidermis of Z.ellipticum in sect. Pleuranthesis. Crystals are found in the abaxial epidermis only in all species, but are present in both epidermal layers of Z.corallinum and Z.montanum. Although the epidermal morphology is similar in most Zingiber species, stomatal density, type of trichome and distribution of oil cells and crystals offer valuable information for the systematic and taxonomic studies in this genus.

Incompatibility Phylogenetic Signals between Double-Digest Restriction Site-Associated DNA Sequencing and Plastid Genomes in Chinese Curcuma (Zingiberaceae)—A Recent Qinghai–Tibetan Plateau Diversification Genera

Curcuma is of high economic value, credited to its medicinal, edible, and ornamental properties, which possess all signatures of adaptability, and rapid radiation, especially species of Curcuma (Chinese Curcuma, a recent Qinghai–Tibetan Plateau diversification genera) scattered in China. However, little is known about the incongruent phylogenetic signals within this genera from different inheritance patterns that will militate against the further development of this genera. In this research, we applied complete chloroplast genome data together with double-digest restriction site-associated DNA sequencing data (ddRAD-seq) strategy to investigate phylogenetic signals of Chinese Curcuma species, clustering using two RAD analysis pipelines (STACKS and pyRAD). Phylogenetic trees were obtained from each locus based on the maximum likelihood (ML) and multispecies coalescent (BEAST) methods. For visual comparison, multi-method and different datasets were used to infer the phylogeny. We discovered inconsistent relationships for the Chinese Curcuma with varying degrees of support using different methods and datasets.

Whence Came These Plants Most Foul? Phylogenomics and Biogeography of Lowiaceae (Zingiberales)

Lowiaceae (order Zingiberales) is a small family of forest herbs in Southeast Asia. All species belong to the genus Orchidantha. They are known for possessing orchid-like flowers that are smelly, apparently mimicking dead animals, feces, or mushrooms. Little is known of the biogeographic patterns or character evolution of the family. We sampled the family extensively, including many recently discovered species, and reconstructed the phylogeny of the family using HybSeq with Lowiaceae-specific RNA baits. Our phylogenetic reconstructions confirm that the family is most closely related to Strelitziaceae, and that species with dark, foul-smelling flowers form a grade in which a clade of species with paler flowers are embedded. The pale-flowered species produce a distinct odor, resembling edible mushrooms. Apart from a single species, the species from Borneo form a clade, and the same is true for Indochinese species. The remaining species form a more widespread clade. A biogeographic analysis shows that the distribution of Lowiaceae can explained by vicariance and gradual dispersal from a shared ancestral range of Borneo and Indochina. There is no evidence of long-distance dispersal, only a later extension in distribution to Peninsular Malaysia which coincides with the presence of a land bridge. Different directions of spread are possible, but none require long-distance dispersal. The results are consistent with the geological history of Southeast Asia. In particular, the relatively early isolation between Indochina and Borneo could be explained by the presence of a sea barrier that developed 10–15 MYA, and the continuous movement of plant species between Borneo and Peninsular Malaysia could be explained by a land bridge that existed until c. 5 MYA. The lack of an extensive land bridge with a suitable habitat may explain the absence of this genus from Sumatra and other Indonesian islands aside from Borneo. The strict reliance on a continuous habitat for the range expansion of Lowiaceae can be explained by their fruits and seeds, which lack obvious adaptations for long-distance dispersal. The inability to disperse to new areas may also explain why the extant species have very restricted distributions.

The Complete Chloroplast Genome Sequences of Eight Fagopyrum Species: Insights Into Genome Evolution and Phylogenetic Relationships

Buckwheat (Fagopyrum genus, Polygonaceae), is an annual or perennial, herbaceous or semi-shrub dicotyledonous plant. There are mainly three cultivated buckwheat species, common buckwheat (Fagopyrum esculentum) is widely cultivated in Asia, Europe, and America, while Tartary buckwheat (F. tataricum) and F. cymosum (also known as F. dibotrys) are mainly cultivated in China. The genus Fagopyrum is taxonomically confusing due to the complex phenotypes of different Fagopyrum species. In this study, the chloroplast (cp) genomes of three Fagopyrum species, F. longistylum, F. leptopodum, F. urophyllum, were sequenced, and five published cp genomes of Fagopyrum were retrieved for comparative analyses. We determined the sequence differentiation, repeated sequences of the cp genomes, and the phylogeny of Fagopyrum species. The eight cp genomes ranged, gene number, gene order, and GC content were presented. Most of variations of Fagopyrum species cp genomes existed in the LSC and SSC regions. Among eight Fagopyrum chloroplast genomes, six variable regions (ndhF-rpl32, trnS-trnG, trnC, trnE-trnT, psbD, and trnV) were detected as promising DNA barcodes. In addition, a total of 66 different SSR (simple sequence repeats) types were found in the eight Fagopyrum species, ranging from 8 to 16 bp. Interestingly, many SSRs showed significant differences especially in some photosystem genes, which provided valuable information for understanding the differences in light adaptation among different Fagopyrum species. Genus Fagopyrum has shown a typical branch that is distinguished from the Rumex, Rheum, and Reynoutria, which supports the unique taxonomic status in Fagopyrum among the Polygonaceae. In addition, phylogenetic analysis based on the cp genomes strongly supported the division of eight Fagopyrum species into two independent evolutionary directions, suggesting that the separation of cymosum group and urophyllum group may be earlier than the flower type differentiation in Fagopyrum plants. The results of the chloroplast-based phylogenetic tree were further supported by the matK and Internal Transcribed Spacer (ITS) sequences of 17 Fagopyrum species, which may help to further anchor the taxonomic status of other members in the urophyllum group in Fagopyrum. This study provides valuable information and high-quality cp genomes for identifying species and evolutionary analysis for future Fagopyrum research.

Molecular evolution of chloroplast genomes in subfamily Zingiberoideae (Zingiberaceae)

BACKGROUND

Zingiberoideae is a large and diverse subfamily of the family Zingiberaceae. Four genera in subfamily Zingiberoideae each possess 50 or more species, including Globba (100), Hedychium (> 80), Kaempferia (50) and Zingiber (150). Despite the agricultural, medicinal and horticultural importance of these species, genomic resources and suitable molecular markers for them are currently sparse.

RESULTS

Here, we have sequenced, assembled and analyzed ten complete chloroplast genomes from nine species of subfamily Zingiberoideae: Globba lancangensis, Globba marantina, Globba multiflora, Globba schomburgkii, Globba schomburgkii var. angustata, Hedychium coccineum, Hedychium neocarneum, Kaempferia rotunda 'Red Leaf', Kaempferia rotunda 'Silver Diamonds' and Zingiber recurvatum. These ten chloroplast genomes (size range 162,630-163,968 bp) possess typical quadripartite structures that consist of a large single copy (LSC, 87,172-88,632 bp), a small single copy (SSC, 15,393-15,917 bp) and a pair of inverted repeats (IRs, 29,673-29,833 bp). The genomes contain 111-113 different genes, including 79 protein coding genes, 28-30 tRNAs and 4 rRNA genes. The dynamics of the genome structures, gene contents, amino acid frequencies, codon usage patterns, RNA editing sites, simple sequence repeats and long repeats exhibit similarities, with slight differences observed among the ten genomes. Further comparative analysis of seventeen related Zingiberoideae species, 12 divergent hotspots are identified. Positive selection is observed in 14 protein coding genes, including accD, ccsA, ndhA, ndhB, psbJ, rbcL, rpl20, rpoC1, rpoC2, rps12, rps18, ycf1, ycf2 and ycf4. Phylogenetic analyses, based on the complete chloroplast-derived single-nucleotide polymorphism data, strongly support that Globba, Hedychium, and Curcuma I + "the Kaempferia clade" consisting of Curcuma II, Kaempferia and Zingiber, form a nested evolutionary relationship in subfamily Zingiberoideae.

CONCLUSIONS

Our study provides detailed information on ten complete Zingiberoideae chloroplast genomes, representing a valuable resource for future studies that seek to understand the molecular evolutionary dynamics in family Zingiberaceae. The identified divergent hotspots can be used for development of molecular markers for phylogenetic inference and species identification among closely related species within four genera of Globba, Hedychium, Kaempferia and Zingiber in subfamily Zingiberoideae.

Authentication of Zingiber Species Based on Analysis of Metabolite Profiles

Zingiber corallinum and Zingiber montanum, which belong to the Zingiberaceae family, are traditional Chinese folk medicinal herbs in Guizhou and Yunnan Province of China. They share great similarities in morphology, chemical constituent, and DNA barcoding sequence. The taxonomy of the two Zingiber species is controversial and discrimination of traditional Chinese medicines directly affects the pharmacological and clinical effects. In the present study, we performed a systemic analysis of "super-barcode" and untargeted metabolomics between Z. corallinum and Z. montanum using chloroplast (cp) genome sequencing and gas chromatography-mass spectrometry (GC-MS) analysis. Comparison and phylogenetic analysis of cp genomes of the two Zingiber species showed that the cp genome could not guarantee the accuracy of identification. An untargeted metabolomics strategy combining GC-MS with chemometric methods was proposed to distinguish the Zingiber samples of known variety. A total of 51 volatile compounds extracted from Z. corallinum and Z. montanum were identified, and nine compounds were selected as candidate metabolic markers to reveal the significant difference between Z. corallinum and Z. montanum. The performance of the untargeted metabolomic approach was verified with unknown Zingiber samples. Although the cp genomes could not be used to identify Zingiber species in this study, it will still provide a valuable genomics resource for population studies in the Zingiberaceae family, and the GC-MS based metabolic fingerprint is more promising for species identification and safe application of Z. corallinum and Z. montanum.

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 Comparative Analyses of Six Complete Chloroplast Genomes of Morphologically Diverse Chenopodium album L. (Amaranthaceae) Collected in Korea

Chenopodium album sensu stricto belonging to C. album aggregate is an annual cosmopolitan weed displaying the diversity of morphologies. We completed the six chloroplast genomes of C. album s. str. collected in Korea to understand the relationship between the diversity of chloroplast genomes and their morphological variations. All six C. album chloroplast genomes have a typical quadripartite structure with length ranging from 151,906 bp to 152,199 bp, similar to the previously sequenced C. album chloroplast genome (NC_034950). In total, 56 single nucleotide polymorphisms (SNPs) and 26 insertion and deletion (INDEL) regions (308 bp in total) were identified from the six chloroplast genomes, presenting a low level of intraspecific variations in comparison to the other angiosperm species. 376 normal simple sequence repeats were identified in all seven C. album chloroplast genomes. The phylogenetic analysis based on all available complete Amaranthaceae chloroplast genomes presents phylogenetic positions of six C. album samples as well as correlation with one of C. album morphological features. Our results provide the way to investigate intraspecific features of C. album chloroplast genomes and also the insights of understanding various intraspecific characteristics including morphological features.