Complete Chloroplast Genome Sequences of Kaempferia Galanga and Kaempferia Elegans: Molecular Structures and Comparative Analysis

Plastome structure, phylogeny and evolution of plastid genes in Reevesia (Helicteroideae, Malvaceae)

Reevesia is an eastern Asian-eastern North American disjunction genus in the family Malvaceae s.l. and comprises approximately 25 species. The relationships within the genus are not well understood. Here, 15 plastomes representing 12 Reevesia species were compared, with the aim of better understanding the species circumscription and phylogenetic relationships within the genus and among genera in the family Malvaceae s.l. The 11 newly sequenced plastomes range between 161,532 and 161, 945 bp in length. The genomes contain 114 unique genes, 18 of which are duplicated in the inverted repeats (IRs). Gene content of these plastomes is nearly identical. All the protein-coding genes are under purifying selection in the Reevesia plastomes compared. The top ten hypervariable regions, SSRs, and the long repeats identified are potential molecular markers for future population genetic and phylogenetic studies. Phylogenetic analysis based on the whole plastomes confirmed the monophyly of Reevesia and a close relationship with Durio (traditional Bombacaceae) in subfamily Helicteroideae, but not with the morphologically similar genera Pterospermum and Sterculia (both of traditional Sterculiaceae). Phylogenetic relationships within Reevesia suggested that two species, R. pubescens and R. thyrsoidea, as newly defined, are not monophyletic. Six taxa, R. membranacea, R. xuefengensis, R. botingensis, R. lofouensis, R. longipetiolata and R. pycnantha, are suggested to be recognized.

Characterization and Comparative Analysis of Complete Chloroplast Genomes of Four Bromus (Poaceae, Bromeae) Species

Bromus (Poaceae Bromeae) is a forage grass with high adaptability and ecological and economic value. Here, we sequenced Bromus ciliatus, Bromus benekenii, Bromus riparius, and Bromus rubens chloroplast genomes and compared them with four previously described species. The genome sizes of Bromus species ranged from 136,934 bp (Bromus vulgaris) to 137,189 bp (Bromus ciliates, Bromus biebersteinii), with a typical quadripartite structure. The studied species had 129 genes, consisting of 83 protein-coding, 38 tRNA-coding, and 8 rRNA-coding genes. The highest GC content was found in the inverted repeat (IR) region (43.85-44.15%), followed by the large single-copy (LSC) region (36.25-36.65%) and the small single-copy (SSC) region (32.21-32.46%). There were 33 high-frequency codons, with those ending in A/U accounting for 90.91%. A total of 350 simple sequence repeats (SSRs) were identified, with single-nucleotide repeats being the most common (61.43%). A total of 228 forward and 141 palindromic repeats were identified. No reverse or complementary repeats were detected. The sequence identities of all sequences were very similar, especially with respect to the protein-coding and inverted repeat regions. Seven highly variable regions were detected, which could be used for molecular marker development. The constructed phylogenetic tree indicates that Bromus is a monophyletic taxon closely related to Triticum. This comparative analysis of the chloroplast genome of Bromus provides a scientific basis for species identification and phylogenetic studies.

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.

The complete chloroplast genome of Mussaenda pubescens and phylogenetic analysis

The chloroplast (cp) genome sequence of Mussaenda pubescens, a promising resource that is used as a traditional medicine and drink, is important for understanding the phylogenetic relationships among the Mussaenda family and genetic improvement and reservation. This research represented the first comprehensive description of the morphological characteristics of M. pubescens, as well as an analysis of the complete cp genome and phylogenetic relationship. The results indicated a close relationship between M. pubescens and M. hirsutula based on the morphological characteristics of the flower and leaves. The cp was sequenced using the Illumina NovaSeq 6000 platform. The results indicated the cp genome of M. pubescens spanned a total length of 155,122 bp, including a pair of inverted repeats (IRA and IRB) with a length of 25,871 bp for each region, as well as a large single-copy (LSC) region and a small single-copy (SSC) region with lengths of 85,370 bp and 18,010 bp, respectively. The results of phylogenetic analyses demonstrated that species within the same genus displayed a tendency to group closely together. It was suggested that Antirhea, Cinchona, Mitragyna, Neolamarckia, and Uncaria might have experienced an early divergence. Furthermore, M. hirsutula showed a close genetic connection to M. pubescens, with the two species having partially overlapping distributions in China. This study presents crucial findings regarding the identification, evolution, and phylogenetic research on Mussaenda plants, specifically targeting M. pubescens.

The newly assembled chloroplast genome of Aeluropus littoralis: molecular feature characterization and phylogenetic analysis with related species

Aeluropus littoralis, a halophyte grass, is widely distributed from the Mediterranean to the Indian subcontinent through the Mongolian Gobi. This model halophyte has garnered increasing attention owing to its use as forage and its high tolerance to environmental stressors. The chloroplast genomes of many plants have been extensively examined for molecular, phylogenetic and transplastomic applications. However, no published research on the A. littoralis chloroplast (cp) genome was discovered. Here, the entire chloroplast genome of A. littoralis was assembled implementing accurate long-read sequences. The entire chloroplast genome, with an estimated length of 135,532 bp (GC content: 38.2%), has a quadripartite architecture and includes a pair of inverted repeat (IR) regions, IRa and IRb (21,012 bp each), separated by a large and a small single-copy regions (80,823 and 12,685 bp, respectively). The features of A. littoralis consist of 133 genes that synthesize 87 peptides, 38 transfer RNAs, and 8 ribosomal RNAs. Of these genes, 86 were unique, whereas 19 were duplicated in IR regions. Additionally, a total of forty-six simple sequence repeats, categorized into 32-mono, four-di, two-tri, and eight-tetranucleotides, were discovered. Furthermore, ten sets of repeats greater than 20 bp were located primarily in the LSC region. Evolutionary analysis based on chloroplast sequence data revealed that A. littoralis with A. lagopoides and A. sinensis belong to the Aeluropodinae subtribe, which is a sister to the Eleusininae in the tribe Cynodonteae and the subfamily Chloridoideae. This subfamily belongs to the PACMAD clade, which contains the majority of the C4 photosynthetic plants in the Poaceae. The newly constructed A. littoralis cp genome offers valuable knowledge for DNA barcoding, phylogenetic, transplastomic research, and other biological studies.

The complete chloroplast genome and phylogenetic analysis of Cyperus malaccensis Lam (Cyperaceae)

Cyperus malaccensis Lam is a perennial herbaceous plant that is distributed over a large area along the southern coast of China. Some plants of the Cyperaceae family are highly similar morphologically, which makes them difficult to classify and identify. In this study, the complete chloroplast genome of C. malaccensis was sequenced and assembled. The chloroplast genome is 186,098 bp long with a 33.18% content of GC. The structure of chloroplast genome includes a quadripartite structure that is composed of a pair of inverted repeats (IRs) of 37,434 bp, a small single copy (SSC) region of 10,296 bp, and a large single copy (LSC) region of 100,934 bp. The genome contains 141 genes, including 94 protein-coding genes, 39 tRNA genes, and 8 rRNA genes. A phylogenetic analysis showed that C. malaccensis is the most closely related to the congeneric species C. rotundus. These results enrich the genetic resources of the Cyperaceae and provide a molecular basis for further study on the phylogeny of this family.

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

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

Genomic exploration of Sesuvium sesuvioides: comparative study and phylogenetic analysis within the order Caryophyllales from Cholistan desert, Pakistan

BACKGROUND

The Aizoaceae family's Sesuvium sesuvioides (Fenzl) Verdc is a medicinal species of the Cholistan desert, Pakistan. The purpose of this study was to determine the genomic features and phylogenetic position of the Sesuvium genus in the Aizoaceae family. We used the Illumina HiSeq2500 and paired-end sequencing to publish the complete chloroplast sequence of S. sesuvioides.

RESULTS

The 155,849 bp length cp genome sequence of S. sesuvioides has a 36.8% GC content. The Leucine codon has the greatest codon use (10.6%), 81 simple sequence repetitions of 19 kinds, and 79 oligonucleotide repeats. We investigated the phylogeny of the order Caryophyllales' 27 species from 23 families and 25 distinct genera. The maximum likelihood tree indicated Sesuvium as a monophyletic genus, and sister to Tetragonia. A comparison of S. sesuvioides, with Sesuvium portulacastrum, Mesembryanthemum crystallinum, Mesembryanthemum cordifolium, and Tetragonia tetragonoides was performed using the NCBI platform. In the comparative investigation of genomes, all five genera revealed comparable cp genome structure, gene number and composition. All five species lacked the rps15 gene and the rpl2 intron. In most comparisons with S. sesuvioides, transition substitutions (Ts) were more frequent than transversion substitutions (Tv), producing Ts/Tv ratios larger than one, and the Ka/Ks ratio was lower than one. We determined ten highly polymorphic regions, comprising rpl22, rpl32-trnL-UAG, trnD-GUC-trnY-GUA, trnE-UUC-trnT-GGU, trnK-UUU-rps16, trnM-CAU-atpE, trnH-GUG-psbA, psaJ-rpl33, rps4-trnT-UGU, and trnF-GAA-ndhJ.

CONCLUSION

The whole S. sesuvioides chloroplast will be examined as a resource for in-depth taxonomic research of the genus when more Sesuvium and Aizoaceae species are sequenced in the future. The chloroplast genomes of the Aizoaceae family are well preserved, with little alterations, indicating the family's monophyletic origin. This study's highly polymorphic regions could be utilized to build realistic and low-cost molecular markers for resolving taxonomic discrepancies, new species identification, and finding evolutionary links among Aizoaceae species. To properly comprehend the evolution of the Aizoaceae family, further species need to be sequenced.

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 genomes provided insights into the evolution and species identification on the Datureae plants

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

Complete chloroplast genome sequence of Camellia sinensis: genome structure, adaptive evolution, and phylogenetic relationships

The chloroplast (cp) genome holds immense potential for a variety of applications including species identification, phylogenetic analysis, and evolutionary studies. In this study, we utilized Illumina NovaSeq 6000 to sequence the DNA of Camellia sinensis L. cultivar 'Zhuyeqi', followed by the assembly of its chloroplast genome using SPAdes v3.10.1, with subsequent analysis of its features and phylogenetic placement. The results showed that the cp genome of 'Zhuyeqi' was 157,072 bp, with a large single-copy region (LSC, 86,628 bp), a small single-copy region (SSC,18,282 bp), and two inverted repeat regions (IR, 26,081 bp). The total AT and GC contents of the cp genome of 'Zhuyeqi' were observed to be 62.21% and 37.29%, respectively. The cp genome encoded 135 unique genes, including 90 protein-coding genes (CDS), 37 tRNA genes, and 8 rRNA genes. Moreover, 31 codons and 247 simple sequence repeats (SSRs) were identified. The cp genomes of 'Zhuyeqi' were found to be relatively conserved, with conservation observed in the IR region, which showed no evidence of inversions or rearrangements. The five regions with the largest variations were identified, with four regions (rps12, rps19, rps16, and rpl33) located in the LSC region and one divergent region (trnI-GAU) in the IR region. Phylogenetic analysis revealed that Camellia sinensis (KJ996106.1) was closely related to 'Zhuyeqi', indicating a close phylogenetic relationship between these two species. These findings could provide important genetic information for further research into breeding of tea tree, phylogeny, and evolution of Camellia sinensis.

Comparative Phylogenetic Analysis for Aerides (Aeridinae, Orchidaceae) Based on Six Complete Plastid Genomes

Aerides Lour. (Orchidaceae, Aeridinae) is a group of epiphytic orchids with high ornamental value, mainly distributed in tropical and subtropical forests, that comprises approximately 20 species. The species are of great value in floriculture and garden designing because of their beautiful flower shapes and colors. Although the morphological boundaries of Aerides are clearly defined, the relationship between Aerides and other closely related genera is still ambiguous in terms of phylogeny. To better understand their phylogenetic relationships, this study used next-generation sequencing technology to investigate the phylogeny and DNA barcoding of this taxonomic unit using genetic information from six Aerides plastid genomes. The quadripartite-structure plastomes ranged from 147,244 bp to 148,391 bp and included 120 genes. Among them, 74 were protein coding genes, 38 were tRNA genes and 8 were rRNA genes, while the ndh genes were pseudogenized or lost. Four non-coding mutational hotspots (rpl20-rpl33, psbM, petB, rpoB-trnCGCA, Pi > 0.06) were identified. A total of 71-77 SSRs and 19-46 long repeats (>30 bp) were recognized in Aerides plastomes, which were mostly located in the large single-copy region. Phylogenetic analysis indicated that Aerides was monophylic and sister to Renanthera. Moreover, our results confirmed that six Aerides species can be divided into three major clades. These findings provide assistance for species identification and DNA barcoding investigation in Aerides, as well as contributes to future research on the phylogenomics of Orchidaceae.

The complete chloroplast genome sequence of Zingiber teres S. Q. Tong & Y. M. Xia (Zingiberaceae)

Here, the complete chloroplast genome sequence of Zingiber teres is described using MGI paired-end sequencing. The genome is 163,428 bp in length and contains a small single-copy region (SSC) of 15,782 bp, a large single-copy region (LSC) of 88,142 bp, and two inverted repeat (IR) regions of 29,752 bp. The overall GC content is 36.1%, and the GC content of the IR regions is 41.1%, which is higher than that of both the LSC region (33.8%) and SSC region (29.5%). The genome of Z. teres contains 133 complete genes, including 88 protein-coding genes (79 protein-coding gene species), 38 tRNA genes (28 tRNA species), and 8 rRNA genes (four rRNA species). Maximum likelihood phylogenetic analysis yielded a well-resolved tree of the genus Zingiber, and Z. teres and Zingiber mioga were sister species in this tree. The development of DNA barcodes could aid the identification of Zingiber species.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

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.

Plastomes of Garcinia mangostana L. and Comparative Analysis with Other Garcinia Species

The two varieties of mangosteen (Garcinia mangostana L.) cultivated in Malaysia are known as Manggis and Mesta. The latter is preferred for its flavor, texture, and seedlessness. Here, we report a complete plastome (156,580 bp) of the Mesta variety that was obtained through a hybrid assembly approach using PacBio and Illumina sequencing reads. It encompasses a large single-copy (LSC) region (85,383 bp) and a small single-copy (SSC) region (17,137 bp) that are separated by 27,230 bp of inverted repeat (IR) regions at both ends. The plastome comprises 128 genes, namely, 83 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The plastome of the Manggis variety (156,582 bp) obtained from reference-guided assembly of Illumina reads was found to be nearly identical to Mesta except for two indels and the presence of a single-nucleotide polymorphism (SNP). Comparative analyses with other publicly available Garcinia plastomes, including G. anomala, G. gummi-gutta, G. mangostana var. Thailand, G. oblongifolia, G. paucinervis, and G. pedunculata, found that the gene content, gene order, and gene orientation were highly conserved among the Garcinia species. Phylogenomic analysis divided the six Garcinia plastomes into three groups, with the Mesta and Manggis varieties clustered closer to G. anomala, G. gummi-gutta, and G. oblongifolia, while the Thailand variety clustered with G. pedunculata in another group. These findings serve as future references for the identification of species or varieties and facilitate phylogenomic analysis of lineages from the Garcinia genus to better understand their evolutionary history.

Characterization of the Plastid Genome of the Vulnerable Endemic Indosasa lipoensis and Phylogenetic Analysis

Indosasa lipoensis, an ornamental garden plant, belongs to the Indosasa genus of the subfamily Bambooaceae within Poaceae. Indosasa lipoensis is endangered and requires protection owing to its relatively narrow distribution area. Chloroplast (cp) genome offers a novel awareness of the evolutionary and genetic variation of higher plants. Herein, we assembled and elucidated the complete cp genome of I. lipoensis, and compared it with four previously published cp genomes from this genus. The I. lipoensis cp genome was 139,655 bp in size, with a typical quadripartite structure, encompassing a large single-copy region (LSC, 83,256 bp), a small single-copy region (SSC, 12,809 bp), and a pair of inverted repeat regions (IR, 21,795 bp). The cp genome consisted of 130 genes with 84 protein-coding genes (CDS), 38 tRNA genes, and 8 rRNA genes. The plastomes were highly conservative, compared to other bamboo species, and exhibited similar patterns of codon usage, number of repeat sequences, and expansion and contraction of the IR boundary. Five hypervariable hotspots were identified as potential DNA barcodes, namely rbcL, petA, petB, trnL-UAG, and ndhE-ndhI, respectively. Phylogenetic analysis based on the complete cp genomes revealed, with high resolution, that I. lipoensis and I. gigantea were most closely related. Overall, these results provided valuable characterization for the future conservation, genetic evaluation, and the breeding of I. lipoensis.

Revisiting chloroplast genomic landscape and annotation towards comparative chloroplast genomes of Rhamnaceae

BACKGROUND

Massive parallel sequencing technologies have enabled the elucidation of plant phylogenetic relationships from chloroplast genomes at a high pace. These include members of the family Rhamnaceae. The current Rhamnaceae phylogenetic tree is from 13 out of 24 Rhamnaceae chloroplast genomes, and only one chloroplast genome of the genus Ventilago is available. Hence, the phylogenetic relationships in Rhamnaceae remain incomplete, and more representative species are needed.

RESULTS

The complete chloroplast genome of Ventilago harmandiana Pierre was outlined using a hybrid assembly of long- and short-read technologies. The accuracy and validity of the final genome were confirmed with PCR amplifications and investigation of coverage depth. Sanger sequencing was used to correct for differences in lengths and nucleotide bases between inverted repeats because of the homopolymers. The phylogenetic trees reconstructed using prevalent methods for phylogenetic inference were topologically similar. The clustering based on codon usage was congruent with the molecular phylogenetic tree. The groups of genera in each tribe were in accordance with tribal classification based on molecular markers. We resolved the phylogenetic relationships among six Hovenia species, three Rhamnus species, and two Ventilago species. Our reconstructed tree provides the most complete and reliable low-level taxonomy to date for the family Rhamnaceae. Similar to other higher plants, the RNA editing mostly resulted in converting serine to leucine. Besides, most genes were subjected to purifying selection. Annotation anomalies, including indel calling errors, unaligned open reading frames of the same gene, inconsistent prediction of intergenic regions, and misannotated genes, were identified in the published chloroplast genomes used in this study. These could be a result of the usual imperfections in computational tools, and/or existing errors in reference genomes. Importantly, these are points of concern with regards to utilizing published chloroplast genomes for comparative genomic analysis.

CONCLUSIONS

In summary, we successfully demonstrated the use of comprehensive genomic data, including DNA and amino acid sequences, to build a reliable and high-resolution phylogenetic tree for the family Rhamnaceae. Additionally, our study indicates that the revision of genome annotation before comparative genomic analyses is necessary to prevent the propagation of errors and complications in downstream analysis and interpretation.

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.

The industrially important genus Kaempferia: An ethnopharmacological review

Kaempferia, a genus of the family Zingiberaceae, is widely distributed with more than 50 species which are mostly found throughout Southeast Asia. These plants have important ethnobotanical significance as many species are used in Ayurvedic and other traditional medicine preparations. This genus has received a lot of scholarly attention recently as a result of the numerous health advantages it possesses. In this review, we have compiled the scientific information regarding the relevance, distribution, industrial applications, phytochemistry, ethnopharmacology, tissue culture and conservation initiative of the Kaempferia genus along with the commercial realities and limitations of the research as well as missing industrial linkages followed by an exploration of some of the likely future promising clinical potential. The current review provides a richer and deeper understanding of Kaempferia, which can be applied in areas like phytopharmacology, molecular research, and industrial biology. The knowledge from this study can be further implemented for the establishment of new conservation strategies.

The antihyperglycemic effect of pulsed electric field-extracted polysaccharide of Kaempferia elegans officinale on streptozotocin induced diabetic mice

Kaempferia elegans polysaccharide (KEP) was extracted using a high-voltage pulsed electric field-assisted hot water method. Its physicochemical properties, in vitro activity and hypoglycemic effect was investigated. Experiments were undertaken with diabetic mice models and the potential mechanism of KEP to improve blood glucose levels was unveiled through measurements of relevant indicators in the serum and liver of the mice. Results showed that KEP is mainly composed of glucose, rhamnose, arabinose, and galactose. It has certain DPPH and ABTS free radical scavenging ability and good α-glucosidase inhibitory ability, indicating that KEP has the potential to improve blood glucose levels in diabetes patients. The experimental results of KEP treatment on mice showed that KEP could control the continuous increase of fasting blood glucose levels. The potential mechanisms behind this blood glucose level control composes of (1) increasing the glucokinase and C peptide levels and decreasing Glucose-6-phosphatase content for improving key enzyme activity in the glucose metabolism pathway. This promotes the consumption of blood glucose during glycolysis, thereby inhibiting the production of endogenous glucose in gluconeogenesis pathway; (2) reducing triglyceride, total cholesterol, low density lipoprotein cholesterol, and increasing high density lipoprotein cholesterol content, for regulating blood lipid indicators to normal levels; and (3) by improving the activities of catalase, glutathione peroxidase, and antioxidant enzymes superoxide dismutase for further improving the antioxidant defense system in the body to reduce blood glucose.

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.

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 of Boesenbergia rotunda and a comparative analysis with members of the family Zingiberaceae

Boesenbergia rotunda (L.) Mansf. is a medically important ginger species of the family Zingiberaceae but its genomic information on molecular phylogeny and identification is scarce. In this work, the chloroplast genome of B. rotunda was sequenced, characterized and compared to the other Zingiberaceae species to provide chloroplast genetic resources and to determine its phylogenetic position in the family. The chloroplast genome of B. rotunda was 163,817 bp in length and consisted of a large single-copy (LSC) region of 88,302 bp, a small single-copy (SSC) region of 16,023 bp and a pair of inverted repeats (IRA and IRB) of 29,746 bp each. The chloroplast genome contained 113 unique genes, including 79 protein-coding genes, 30 transfer RNA (tRNA) genes and four ribosomal RNA (rRNA) genes. Several genes had atypical start codons, while most amino acids exhibited biased usage of synonymous codons. Comparative analyses with various chloroplast genomes of Zingiberaceae taxa revealed several highly variable regions (psbK-psbI, trnT-GGU-psbD, rbcL-accD, ndhF-rpl32, and ycf1) in the LSC and SSC regions in the chloroplast genome of B. rotunda that could be utilized as molecular markers for DNA barcoding and species delimitation. Phylogenetic analyses based on shared protein-coding genes revealed that B. rotunda formed a distinct lineage with B. kingii Mood & L.M.Prince, in a subclade that also contained the genera Kaempferia and Zingiber. These findings constitute the first chloroplast genome information of B. rotunda that could be a reference for phylogenetic analysis and identification of genus Boesenbergia within the Zingiberaceae family.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s40415-022-00845-w.

The complete chloroplast genomes of three Hamamelidaceae species: Comparative and phylogenetic analyses

Hamamelidaceae is an important group that represents the origin and early evolution of angiosperms. Its plants have many uses, such as timber, medical, spice, and ornamental uses. In this study, the complete chloroplast genomes of Loropetalum chinense (R. Br.) Oliver, Corylopsis glandulifera Hemsl., and Corylopsis velutina Hand.-Mazz. were sequenced using the Illumina NovaSeq 6000 platform. The sizes of the three chloroplast genomes were 159,402 bp (C. glandulifera), 159,414 bp (C. velutina), and 159,444 bp (L. chinense), respectively. These chloroplast genomes contained typical quadripartite structures with a pair of inverted repeat (IR) regions (26,283, 26,283, and 26,257 bp), a large single-copy (LSC) region (88,134, 88,146, and 88,160 bp), and a small single-copy (SSC) region (18,702, 18,702, and 18,770 bp). The chloroplast genomes encoded 132-133 genes, including 85-87 protein-coding genes, 37-38 tRNA genes, and 8 rRNA genes. The coding regions were composed of 26,797, 26,574, and 26,415 codons, respectively, most of which ended in A/U. A total of 37-43 long repeats and 175-178 simple sequence repeats (SSRs) were identified, and the SSRs contained a higher number of A + T than G + C bases. The genome comparison showed that the IR regions were more conserved than the LSC or SSC regions, while the noncoding regions contained higher variability than the gene coding regions. Phylogenetic analyses revealed that species in the same genus tended to cluster together. Chunia Hung T. Chang, Mytilaria Lecomte, and Disanthus Maxim. may have diverged early and Corylopsis Siebold & Zucc. was closely related to Loropetalum R. Br. This study provides valuable information for further species identification, evolution, and phylogenetic studies of Hamamelidaceae plants.

Analysis of Chloroplast Genomes Provides Insights Into the Evolution of Agropyron

Plants of the Agropyron genus are important pasture resources, and they also play important roles in the ecological restoration. Chloroplast genomes are inherited from maternal parents, and they are important for studying species taxonomy and evolution. In this study, we sequenced the complete chloroplast genomes of five typical species of the Agropyron genus (eg., A. cristatum × A. desertorum Fisch. Schult, A. desertorum, A. desertorum Fisch. Schult. cv. Nordan, A. michnoi Roshev, and A. mongolicum Keng) using the Illumina NovaSeq platform. We found that these five chloroplast genomes exhibit a typical quadripartite structure with a conserved genome arrangement and structure. Their chloroplast genomes contain the large single-copy regions (LSC, 79,613 bp-79,634 bp), the small single-copy regions (SSC, 12,760 bp-12,768 bp), and the inverted repeat regions (IR, 43,060 bp-43,090 bp). Each of the five chloroplast genomes contains 129 genes, including 38 tRNA genes, eight rRNA genes, and 83 protein-coding genes. Among them, the genes trnG-GCC, matK, petL, ccsA, and rpl32 showed significant nucleotide diversity in these five species, and they may be used as molecular markers in taxonomic studies. Phylogenetic analysis showed that A. mongolicum is closely related to A. michnoi, while others have a closer genetic relationship with the Triticum genus.

Complete Chloroplast Genome Sequence of Fagus longipetiolata Seemen (Fagaceae): Genome Structure, Adaptive Evolution, and Phylogenetic Relationships

Fagus longipetiolata Seemen is a deciduous tree of the Fagus genus in Fagaceae, which is endemic to China. In this study, we successfully sequenced the cp genome of F. longipetiolata, compared the cp genomes of the Fagus genus, and reconstructed the phylogeny of Fagaceae. The results showed that the cp genome of F. longipetiolata was 158,350 bp, including a pair of inverted repeat (IRA and IRB) regions with a length of 25,894 bp each, a large single-copy (LSC) region of 87,671 bp, and a small single-copy (SSC) region of 18,891 bp. The genome encoded 131 unique genes, including 81 protein-coding genes, 37 transfer RNA genes (tRNAs), 8 ribosomal RNA genes (rRNAs), and 5 pseudogenes. In addition, 33 codons and 258 simple sequence repeats (SSRs) were identified. The cp genomes of Fagus were relatively conserved, especially the IR regions, which showed the best conservation, and no inversions or rearrangements were found. The five regions with the largest variations were the rps12, rpl32, ccsA, trnW-CCA, and rps3 genes, which spread over in LSC and SSC. The comparison of gene selection pressure indicated that purifying selection was the main selective pattern maintaining important biological functions in Fagus cp genomes. However, the ndhD, rpoA, and ndhF genes of F. longipetiolata were affected by positive selection. Phylogenetic analysis revealed that F. longipetiolata and F. engleriana formed a close relationship, which partially overlapped in their distribution in China. Our analysis of the cp genome of F. longipetiolata would provide important genetic information for further research into the classification, phylogeny and evolution of Fagus.

Comparative genomics and phylogenetic relationships of two endemic and endangered species (Handeliodendron bodinieri and Eurycorymbus cavaleriei) of two monotypic genera within Sapindales

BACKGROUND

Handeliodendron Rehder and Eurycorymbus Hand.-Mazz. are the monotypic genera in the Sapindaceae family. The phylogenetic relationship of these endangered species Handeliodendron bodinieri (Lévl.) Rehd. and Eurycorymbus cavaleriei (Lévl.) Rehd. et Hand.-Mazz. with other members of Sapindaceae s.l. is not well resolved. A previous study concluded that the genus Aesculus might be paraphyletic because Handeliodendron was nested within it based on small DNA fragments. Thus, their chloroplast genomic information and comparative genomic analysis with other Sapindaceae species are necessary and crucial to understand the circumscription and plastome evolution of this family.

RESULTS

The chloroplast genome sizes of Handeliodendron bodinieri and Eurycorymbus cavaleriei are 151,271 and 158,690 bp, respectively. Results showed that a total of 114 unique genes were annotated in H. bodinieri and E. cavaleriei, and the ycf1 gene contained abundant SSRs in both genomes. Comparative analysis revealed that gene content, PCGs, and total GC content were remarkably similar or identical within 13 genera from Sapindaceae, and the chloroplast genome size of four genera was generally smaller within the family, including Acer, Dipteronia, Aesculus, and Handeliodendron. IR boundaries of the H. bodinieri showed a significant contraction, whereas it presented a notable expansion in E. cavaleriei cp genome. Ycf1, ndhC-trnV-UAC, and rpl32-trnL-UAG-ccsA were remarkably divergent regions in the Sapindaceae species. Analysis of selection pressure showed that there are a few positively selected genes. Phylogenetic analysis based on different datasets, including whole chloroplast genome sequences, coding sequences, large single-copy, small single-copy, and inverted repeat regions, consistently demonstrated that H. bodinieri was sister to the clade consisting of Aesculus chinensis and A. wangii and strongly support Eurycorymbus cavaleriei as sister to Dodonaea viscosa.

CONCLUSION

This study revealed that the cp genome size of the Hippocastanoideae was generally smaller compared to the other subfamilies within Sapindaceae, and three highly divergent regions could be used as the specific DNA barcodes within Sapindaceae. Phylogenetic results strongly support that the subdivision of four subfamilies within Sapindaceae, and Handeliodendron is not nested within the genus Aesculus.

Comparative and Phylogenetic Analysis Based on the Chloroplast Genome of Coleanthus subtilis (Tratt.) Seidel, a Protected Rare Species of Monotypic Genus

Coleanthus subtilis (Tratt.) Seidel (Poaceae) is an ephemeral grass from the monotypic genus Coleanthus Seidl, which grows on wet muddy areas such as fishponds or reservoirs. As a rare species with strict habitat requirements, it is protected at international and national levels. In this study, we sequenced its whole chloroplast genome for the first time using the next-generation sequencing (NGS) technology on the Illumina platform, and performed a comparative and phylogenetic analysis with the related species in Poaceae. The complete chloroplast genome of C. subtilis is 135,915 bp in length, with a quadripartite structure having two 21,529 bp inverted repeat regions (IRs) dividing the entire circular genome into a large single copy region (LSC) of 80,100 bp and a small single copy region (SSC) of 12,757 bp. The overall GC content is 38.3%, while the GC contents in LSC, SSC, and IR regions are 36.3%, 32.4%, and 43.9%, respectively. A total of 129 genes were annotated in the chloroplast genome, including 83 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. The accD gene and the introns of both clpP and rpoC1 genes were missing. In addition, the ycf1, ycf2, ycf15, and ycf68 were pseudogenes. Although the chloroplast genome structure of C. subtilis was found to be conserved and stable in general, 26 SSRs and 13 highly variable loci were detected, these regions have the potential to be developed as important molecular markers for the subfamily Pooideae. Phylogenetic analysis with species in Poaceae indicated that Coleanthus and Phippsia were sister groups, and provided new insights into the relationship between Coleanthus, Zingeria, and Colpodium. This study presents the initial chloroplast genome report of C. subtilis, which provides an essential data reference for further research on its origin.

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.

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.

New Insight into Taxonomy of European Mountain Pines, Pinus mugo Complex, Based on Complete Chloroplast Genomes Sequencing

The Pinus mugo complex is a large group of closely related mountain pines, which are an important component of the ecosystems of the most important mountain ranges, such as the Alps, Carpathians and Pyrenees. The phylogenetic relationships between taxa in this complex have been under discussion for many years. Despite the use of many different approaches, they still need to be clarified and supplemented with new data, especially those obtained with high-throughput methods. Therefore, in this study, the complete sequences of the chloroplast genomes of the three most recognized members of the Pinus mugo complex, i.e., Pinus mugo, Pinus rotundata and Pinus uncinata, were sequenced and analyzed to gain new insight into their phylogenetic relationships. Comparative analysis of their complete chloroplast genome sequences revealed several mutational hotspots potentially useful for the genetic identification of taxa from the Pinus mugo complex. Phylogenetic inference based on sixteen complete chloroplast genomes of different coniferous representatives showed that pines from the Pinus mugo complex form one distinct monophyletic group. The results obtained in this study provide new and valuable omics data for further research within the European mountain pine complex. They also indicate which regions may be useful in the search for diagnostic DNA markers for the members of Pinus mugo complex and set the baseline in the conservation of genetic resources of its endangered taxa.

Complete Chloroplast Genome of Clethra fargesii Franch., an Original Sympetalous Plant from Central China: Comparative Analysis, Adaptive Evolution, and Phylogenetic Relationships

Clethra fargesii, an essential ecological and endemic woody plant of the genus Clethra in Clethraceae, is widely distributed in Central China. So far, there have been a paucity of studies on its chloroplast genome. In the present study, we sequenced and assembled the complete chloroplast genome of C. fargesii. We also analyzed the chloroplast genome features and compared them to Clethra delavayi and other closely related species in Ericales. The complete chloroplast genome is 157,486 bp in length, including a large single-copy (LSC) region of 87,034 bp and a small single-copy (SSC) region of 18,492 bp, separated by a pair of inverted repeat (IR) regions of 25,980 bp. The GC content of the whole genome is 37.3%, while those in LSC, SSC, and IR regions are 35.4%, 30.7%, and 43.0%, respectively. The chloroplast genome of C. fargesii encodes 132 genes in total, including 87 protein-coding genes (PCGs), 37 tRNA genes, and eight rRNA genes. A total of 26,407 codons and 73 SSRs were identified in C. fargesii chloroplast genome. Additionally, we postulated and demonstrated that the structure of the chloroplast genome in Clethra species may present evolutionary conservation based on the comparative analysis of genome features and genome alignment among eight Ericales species. The low Pi values revealed evolutionary conservation based on the nucleotide diversity analysis of chloroplast genome in two Clethra species. The low selection pressure was shown by a few positively selected genes by adaptive evolution analysis using 80 coding sequences (CDSs) of the chloroplast genomes of two Clethra species. The phylogenetic tree showed that Clethraceae and Ericaceae are sister clades, which reconfirm the previous hypothesis that Clethra is highly conserved in the chloroplast genome using 75 CDSs of chloroplast genome among 40 species. The genome information and analysis results presented in this study are valuable for further study on the intraspecies identification, biogeographic analysis, and phylogenetic relationship in Clethraceae.

Complete chloroplast genomes of Zingiber montanum and Zingiber zerumbet: Genome structure, comparative and phylogenetic analyses

Zingiber montanum (Z. montanum) and Zingiber zerumbet (Z. zerumbet) are important medicinal and ornamental herbs in the genus Zingiber and family Zingiberaceae. Chloroplast-derived markers are useful for species identification and phylogenetic studies, but further development is warranted for these two Zingiber species. In this study, we report the complete chloroplast genomes of Z. montanum and Z. zerumbet, which had lengths of 164,464 bp and 163,589 bp, respectively. These genomes had typical quadripartite structures with a large single copy (LSC, 87,856-89,161 bp), a small single copy (SSC, 15,803-15,642 bp), and a pair of inverted repeats (IRa and IRb, 29,393-30,449 bp). We identified 111 unique genes in each chloroplast genome, including 79 protein-coding genes, 28 tRNAs and 4 rRNA genes. We analyzed the molecular structures, gene information, amino acid frequencies, codon usage patterns, RNA editing sites, simple sequence repeats (SSRs) and long repeats from the two chloroplast genomes. A comparison of the Z. montanum and Z. zerumbet chloroplast genomes detected 489 single-nucleotide polymorphisms (SNPs) and 172 insertions/deletions (indels). Thirteen highly divergent regions, including ycf1, rps19, rps18-rpl20, accD-psaI, psaC-ndhE, psbA-trnK-UUU, trnfM-CAU-rps14, trnE-UUC-trnT-UGU, ccsA-ndhD, psbC-trnS-UGA, start-psbA, petA-psbJ, and rbcL-accD, were identified and might be useful for future species identification and phylogeny in the genus Zingiber. Positive selection was observed for ATP synthase (atpA and atpB), RNA polymerase (rpoA), small subunit ribosomal protein (rps3) and other protein-coding genes (accD, clpP, ycf1, and ycf2) based on the Ka/Ks ratios. Additionally, chloroplast SNP-based phylogeny analyses found that Zingiber was a monophyletic sister branch to Kaempferia and that chloroplast SNPs could be used to identify Zingiber species. The genome resources in our study provide valuable information for the identification and phylogenetic analysis of the genus Zingiber and family Zingiberaceae.

Complete Chloroplast Genome Sequence and Phylogenetic Inference of the Canary Islands Dragon Tree (Dracaena draco L.)

Dracaena draco, which belongs to the genus Dracaena, is an endemic succulent of the Canary Islands. Although it is one of the most popular and widely grown ornamental plants in the world, little is known about its genomic variability. Next generation sequencing, especially in combination with advanced bioinformatics analysis, is a new standard in taxonomic and phylogenetic research. Therefore, in this study, the complete D. draco chloroplast genome (cp) was sequenced and analyzed in order to provide new genomic information and to elucidate phylogenetic relationships, particularly within the genus Dracaena. The D. draco chloroplast genome is 155,422 bp, total guanine-cytosine (GC) content is 37.6%, and it has a typical quadripartite plastid genome structure with four separate regions, including one large single copy region of 83,942 bp length and one small single copy region of 18,472 bp length, separated by two inverted repeat regions, each 26,504 bp in length. One hundred and thirty-two genes were identified, 86 of which are protein-coding genes, 38 are transfer RNAs, and eight are ribosomal RNAs. Seventy-seven simple sequence repeats were also detected. Comparative analysis of the sequence data of various members of Asparagales revealed mutational hotspots potentially useful for their genetic identification. Phylogenetic inference based on 16 complete chloroplast genomes of Asparagales strongly suggested that Dracaena species form one monophyletic group, and that close relationships exist between D. draco, D. cochinchinensis and D. cambodiana. This study provides new and valuable data for further taxonomic, evolutionary and phylogenetic studies within the Dracaena genus.

Complete Chloroplast Genomes of Three Medicinal Alpinia Species: Genome Organization, Comparative Analyses and Phylogenetic Relationships in Family Zingiberaceae

Alpinia katsumadai (A. katsumadai), Alpinia oxyphylla (A. oxyphylla) and Alpinia pumila (A. pumila), which belong to the family Zingiberaceae, exhibit multiple medicinal properties. The chloroplast genome of a non-model plant provides valuable information for species identification and phylogenetic analysis. Here, we sequenced three complete chloroplast genomes of A. katsumadai, A. oxyphylla sampled from Guangdong and A. pumila, and analyzed the published chloroplast genomes of Alpinia zerumbet (A. zerumbet) and A. oxyphylla sampled from Hainan to retrieve useful chloroplast molecular resources for Alpinia. The five Alpinia chloroplast genomes possessed typical quadripartite structures comprising of a large single copy (LSC, 87,248-87,667 bp), a small single copy (SSC, 15,306-18,295 bp) and a pair of inverted repeats (IR, 26,917-29,707 bp). They had similar gene contents, gene orders and GC contents, but were slightly different in the numbers of small sequence repeats (SSRs) and long repeats. Interestingly, fifteen highly divergent regions (rpl36, ycf1, rps15, rpl22, infA, psbT-psbN, accD-psaI, petD-rpoA, psaC-ndhE, ccsA-ndhD, ndhF-rpl32, rps11-rpl36, infA-rps8, psbC-psbZ, and rpl32-ccsA), which could be suitable for species identification and phylogenetic studies, were detected in the Alpinia chloroplast genomes. Comparative analyses among the five chloroplast genomes indicated that 1891 mutational events, including 304 single nucleotide polymorphisms (SNPs) and 118 insertion/deletions (indels) between A. pumila and A. katsumadai, 367 SNPs and 122 indels between A. pumila and A. oxyphylla sampled from Guangdong, 331 SNPs and 115 indels between A. pumila and A. zerumbet, 371 SNPs and 120 indels between A. pumila and A. oxyphylla sampled from Hainan, and 20 SNPs and 23 indels between the two accessions of A. oxyphylla, were accurately located. Additionally, phylogenetic relationships based on SNP matrix among 28 whole chloroplast genomes showed that Alpinia was a sister branch to Amomum in the family Zingiberaceae, and that the five Alpinia accessions were divided into three groups, one including A. pumila, another including A. zerumbet and A. katsumadai, and the other including two accessions of A. oxyphylla. In conclusion, the complete chloroplast genomes of the three medicinal Alpinia species in this study provided valuable genomic resources for further phylogeny and species identification in the family Zingiberaceae.

The Chloroplast Genome of Carya illinoinensis: Genome Structure, Adaptive Evolution, and Phylogenetic Analysis

Research Highlights: For the first time, the complete chloroplast (cp) genome of Carya illinoinensis cv. ‘Pawnee’ was de novo assembled. Comprehensive analysis the cp genome of C. illinoinensis revealed potential cpDNA markers for intraspecies identification, genes involved in adaptation, and its phylogenetic position. Background and Objectives: C. illinoinensis is an economically important nut tree in the family Juglandaceae. Cp-derived markers are helpful for genetic research, but they still need to be developed in C. illinoinensis. Additionally, the adaptation and phylogenetic relationships of C. illinoinensis have not been revealed based on the complete cp genome. Materials and Methods: Chloroplast genomic DNA of C. illinoinensis cv. ‘Pawnee’ was extracted and subjected to Illumina sequencing. Results: The cp genome is 160,819 bp in size, exhibiting a typical quadripartite structure with a large single copy (LSC) of 90,022 bp, a small single copy (SSC) of 18,791 bp, and a pair of inverted repeats (IRA and IRB) regions of 26,003 bp each. The genome was predicted to encode 112 unique genes, including 79 protein-coding genes, 29 tRNAs, and four rRNAs, with 19 duplicates in the IR regions. In total, 213 SSRs and 44 long repeats were identified in the cp genome. A comparison of two different C. illinoinensis genotypes, ‘Pawnee’ and 87MX3-2.11, obtained 143 SNPs and 74 indels. The highly variable regions such as atpF, clpP, and ndhA genes, and matK-rps16, trnS-trnG, and trnT-psbD intergenic spacers might be helpful for future intraspecific identification. Positive selection was acting on the ccsA and rps12 cp genes based on the Ka/Ks ratios. Phylogenetic analysis indicated that C. illinoinensis forms a sister clade to Asian Carya species, represented by C. kweichowensis and Annamocarya sinensis. Conclusions: The genome information in our study will have significance for further research on the intraspecies identification and genetic improvement of C. illinoinensis.

Comparative analysis of chloroplast genomes for five Dicliptera species (Acanthaceae): molecular structure, phylogenetic relationships, and adaptive evolution

The genus Dicliptera (Justicieae, Acanthaceae) consists of approximately 150 species distributed throughout the tropical and subtropical regions of the world. Newly obtained chloroplast genomes (cp genomes) are reported for five species of Dilciptera (D. acuminata, D. peruviana, D. montana, D. ruiziana and D. mucronata) in this study. These cp genomes have circular structures of 150,689-150,811 bp and exhibit quadripartite organizations made up of a large single copy region (LSC, 82,796-82,919 bp), a small single copy region (SSC, 17,084-17,092 bp), and a pair of inverted repeat regions (IRs, 25,401-25,408 bp). Guanine-Cytosine (GC) content makes up 37.9%-38.0% of the total content. The complete cp genomes contain 114 unique genes, including 80 protein-coding genes, 30 transfer RNA (tRNA) genes, and four ribosomal RNA (rRNA) genes. Comparative analyses of nucleotide variability (Pi) reveal the five most variable regions (trnY-GUA-trnE-UUC, trnG-GCC, psbZ-trnG-GCC, petN-psbM, and rps4-trnL-UUA), which may be used as molecular markers in future taxonomic identification and phylogenetic analyses of Dicliptera. A total of 55-58 simple sequence repeats (SSRs) and 229 long repeats were identified in the cp genomes of the five Dicliptera species. Phylogenetic analysis identified a close relationship between D. ruiziana and D. montana, followed by D. acuminata, D. peruviana, and D. mucronata. Evolutionary analysis of orthologous protein-coding genes within the family Acanthaceae revealed only one gene, ycf15, to be under positive selection, which may contribute to future studies of its adaptive evolution. The completed genomes are useful for future research on species identification, phylogenetic relationships, and the adaptive evolution of the Dicliptera species.

Chloroplast Genome Sequence of Artemisia scoparia: Comparative Analyses and Screening of Mutational Hotspots

Artemisia L. is among the most diverse and medicinally important genera of the plant family Asteraceae. Discrepancies arise in the taxonomic classification of Artemisia due to the occurrence of multiple polyploidy events in separate lineages and its complex morphology. The discrepancies could be resolved by increasing the genomic resources. A. scoparia is one of the most medicinally important species in Artemisia. In this paper, we report the complete chloroplast genome sequence of Artemisia scoparia. The genome was 151,060 bp (base pairs), comprising a large single copy (82,834 bp) and small single copy (18,282 bp), separated by a pair of long inverted repeats (IRa and IRb: 24,972 bp each). We identified 114 unique genes, including four ribosomal RNAs, 30 transfer RNAs, and 80 protein-coding genes. We analysed the chloroplast genome features, including oligonucleotide repeats, microsatellites, amino acid frequencies, RNA editing sites, and codon usage. Transversion substitutions were twice as frequent as transition substitutions. Mutational hotspot loci included ccsA-ndhD, trnH-psbA, ndhG-ndhI, rps18-rpl20, and rps15-ycf1. These loci can be used to develop cost-effective and robust molecular markers for resolving the taxonomic discrepancies. The reconstructed phylogenetic tree supported previous findings of Artemisia as a monophyletic genus, sister to the genus Chrysanthemum, whereby A. scoparia appeared as sister to A. capillaris.

Plant microRNAs: biogenesis, gene silencing, web-based analysis tools and their use as molecular markers

MicroRNAs (miRNAs) are tiny (20-24 nt bp) regulatory non-protein-coding RNA molecules that have been extensively characterized and found important for many physiological and developmental processes. The miss-expression of miRNAs leads to various defects in plants. MicroRNAs repress gene expression by directing mRNA degradation or translational arrest. Several proteins such as PP43A, HYL1, DCL, HST are indispensable role players in promoting miRNA biogenesis in plants. During miRNA biogenesis, lariat RNAs are produced as by-products of pre-mRNA splicing which have a negative role in regulation of miRNA homeostasis. By acting as a decoy and by sequestering to the dicing complex, lariat RNA can prevent the processing of miRNAs. A number of bioinformatic tools with different methodologies are available to identify and validate miRNAs and their targets. Many miRNAs have been reported in different crops for different traits; however, no reports are available on their use in plant breeding. Recently, researchers have developed trait specific miRNA-based molecular markers (miRNA-SSRs/SNP) for many quantitative traits in different plant species. In the future, these molecular markers can be used for plant breeding programs. In this review, a comprehensive up-to-date information is provided on the bioinformatic tools used for analysis of plant miRNAs and their targets, the number of miRNAs, their biogenesis, gene silencing mechanism and miRNA-based molecular markers.

The Complete Chloroplast Genome of Trichopus zeylanicus, And Phylogenetic Analysis with Dioscoreales

We presents the first chloroplast genome from the genus Trichopus. Comparative analysis revealed that the IR regions are more conserved than the SC regions. Highly divergent sequence hot spots were identified, which could be used as molecular markers. Phylogenetic analysis gave insight into the evolutionary history of Trichopus zeylanicus. In this study, we determined the complete sequence of the chloroplast genome of an important, rare, and endangered medicinal plant, Trichopus zeylanicus. The analysis of the genome showed that the complete chloroplast genome of Trichopus zeylanicus is 153,497 bp in size, and has a quadripartite structure with a large single copy of 81,091 bp and a small single copy of 17,512 bp separated by inverted repeats of 27,447 bp. Sequence analysis revealed that the chloroplast genome encodes 112 unique genes, including 78 protein-coding genes, 30 tRNA genes, and four rRNA genes. We also identified 95 simple sequence repeats and 54 long repeats including 34 forward repeats, seven inverted repeats, nine palindromes, three reverse repeats, and one complementary repeat within the chloroplast genome of Trichopus zeylanicus. Whole chloroplast genome comparison with those of other Dioscoreales indicated that the inverted regions are more conserved than large single copy and small single copy regions. In the phylogenetic trees based on complete chloroplast genome and 78 shared chloroplast protein-coding genes in 15 monocot species, including 14 Dioscoreales, Trichopus zeylanicus formed a distinct clade. In summary, the first chloroplast genome from the genus Trichopus reported in this study gave a better insight into the phylogenetic relationships of different genera within the order Dioscoreales. Moreover, the present data will be a valuable chloroplast genomic resource for population genetics.

Complete chloroplast genome sequence and phylogenetic analysis of Spathiphyllum 'Parrish'

Spathiphyllum is a very important tropical plant used as a small, potted, ornamental plant in South China, with an annual output value of hundreds of millions of yuan. In this study, we sequenced and analyzed the complete nucleotide sequence of the Spathiphyllum 'Parrish' chloroplast genome. The whole chloroplast genome is 168,493 bp in length, and includes a pair of inverted repeat (IR) regions (IRa and IRb, each 31,600 bp), separated by a small single-copy (SSC, 15,799 bp) region and a large single-copy (LSC, 89,494 bp) region. Our annotation revealed that the S. 'Parrish' chloroplast genome contained 132 genes, including 87 protein coding genes, 37 transfer RNA genes, and 8 ribosomal RNA genes. In the repeat structure analysis, we detected 281 simple sequence repeats (SSRs) which included mononucleotides (223), dinucleotides (28), trinucleotides (12), tetranucleotides (11), pentanucleotides (6), and hexanucleotides (1), in the S. 'Parrish' chloroplast genome. In addition, we identified 50 long repeats, comprising 18 forward repeats, 13 reverse repeats, 17 palindromic repeats, and 2 complementary repeats. Single nucleotide polymorphism (SNP) and insertion/deletion (indel) analyses of the chloroplast genome of the S. 'Parrish' relative S. cannifolium revealed 962 SNPs in S. 'Parrish'. There were 158 indels (90 insertions and 68 deletions) in the S. 'Parrish' chloroplast genome relative to the S. cannifolium chloroplast genome. Phylogenetic analysis of five species found S. 'Parrish' to be more closely related to S. kochii than to S. cannifolium. This study identified the characteristics of the S. 'Parrish' chloroplast genome, which will facilitate species identification and phylogenetic analysis within the genus Spathiphyllum.

Characterization and phylogenetic analysis of the complete chloroplast genome of Curcuma longa (Zingiberaceae)

Curcuma longa, a well-known traditional medicinal plant in China, belongs to the genus Curcuma family Zingiberaceae. In this study, we firstly assembled the complete chloroplast genome of C. longa based on sequences from Illumina and PacBio sequencing platforms. We obtained the complete chloroplast genome with the total length of 162,176 bp. It consisted of a large single-copy region (LSC, 86,984 bp), a small single-copy region (SSC, 15,694 bp), and a pair of inverted repeats (IRs, 29,749 bp each). Sequence analyses indicated that the chloroplast genome contained 111 distinct genes including 79 protein-coding genes, 28 tRNA genes, and four rRNA genes. The nucleotide composition was asymmetric (31.62% A, 18.42% C, 17.79% G, 32.18% T) with an overall AT content of 63.80%. The AT contents of the LSC, SSC and IR regions were 66.00%, 70.35% and 58.85%, respectively. Sixteen genes owned a single intron, while another two genes had two introns. The phylogenetic analysis indicated that C. longa was closely related to species Curcuma roscoeana within the genus Curcuma in family Zingiberaceae.

Complete chloroplast genome sequence of Hedychium coronarium

The first complete chloroplast genome of Hedychium coronarium (Zingiberaceae) was reported in this study. The H. coronarium chloroplast genome was 163,949 bp in length and comprised a pair of inverted repeat (IR) regions of 29,780 bp each, a large single-copy (LSC) region of 88,581 bp and a small single-copy (SSC) region of 15,808 bp. It encoded 141 genes, including 87 protein-coding genes (79 PCG species), 46 tRNA genes (28 tRNA species), and eight rRNA genes (four rRNA species). The nucleotide composition was asymmetric (31.68% A, 18.35% C, 17.74% G, 32.23% T) with an overall AT content of 63.92%. Phylogenetic analysis showed that H. coronarium was classified into a monophyletic group within the genus Hedychium in family Zingiberaceae.

Complete chloroplast genome of the plant Stahlianthus Involucratus (Zingiberaceae)

The first complete chloroplast genome of Stahlianthus involucratus (Zingiberaceae) was reported in this study. The S. involucratus chloroplast genome was 163,300 bp in length and consisted of one large single copy (LSC) region of 87,498 bp, one small single copy (SSC) region of 15,568 bp, and a pair of inverted repeat (IR) regions 30,117 bp. It encoded 141 genes, including 87 protein-coding genes (79 PCG species), 46 tRNA genes (28 tRNA species) and 8 rRNA genes (4 rRNA species). The phylogenetic analysis based on single nucleotide polymorphisms strongly supported that S. involucratus, Curcuma roscoeana and Curcuma longa formed a cluster in group CurcumaII within family Zingiberaceae.