Comparative Analyses of Chloroplast Genomes of Cucurbitaceae Species: Lights into Selective Pressures and Phylogenetic Relationships

The complete chloroplast genome sequence of Lemna turionifera (Araceae)

Lemna turionifera is native to North America and northern Asia, with significant potential for industrial wastewater remediation. The complete nucleotide sequence of the L. turionifera chloroplast genome (cpDNA) was determined. The cpDNA is a circular molecule of 166,606 bp and containing a pair of inverted repeats (IRs) measuting 31,663 bp each. These IRs are flanked by a small single-copy region of 13,542 bp and a large single-copy region of 89,738 bp. The chloroplast genome of L. turionifera consisted of 112 unique genes, including 78 protein-encoding genes, 30 tRNA genes, and four rRNA genes. The phylogenetic analysis utilizing cpDNA provided a well-supported resolution of the relationships among subfamilies within the Araceae family. Our findings indicated a close relationship between L. turionifera and a clade consisting of L. minor, L. japonica, and L. gibba. The availability of the complete chloroplast genome sequence of L. turionifera presents valuable data for future phylogenetic investigations within the Lemnaceae family.

The First Complete Chloroplast Genome of Cordia monoica: Structure and Comparative Analysis

Cordia monoica is a member of the Boraginaceae family. This plant is widely distributed in tropical regions and has a great deal of medical value as well as economic importance. In the current study, the complete chloroplast (cp) genome of C. monoica was sequenced, assembled, annotated, and reported. This circular chloroplast genome had a size of 148,711 bp, with a quadripartite structure alternating between a pair of repeated inverted regions (26,897-26,901 bp) and a single copy region (77,893 bp). Among the 134 genes encoded by the cp genome, there were 89 protein-coding genes, 37 transfer RNA (tRNA) genes, and 8 ribosomal RNA (rRNA) genes. A total of 1387 tandem repeats were detected, with the hexanucleotides class making up 28 percent of the repeats. Cordia monoica has 26,303 codons in its protein-coding regions, and leucine amino acid was the most frequently encoded amino acid in contrast to cysteine. In addition, 12 of the 89 protein-coding genes were found to be under positive selection. The phyloplastomic taxonomical clustering of the Boraginaceae species provides further evidence that chloroplast genome data are reliable not only at family level but also in deciphering the phylogeny at genus level (e.g., Cordia).

Comparative chloroplast genome analysis of seven extant Citrullus species insight into genetic variation, phylogenetic relationships, and selective pressure

Citrullus ecirrhosus, Citrullus rehmii, and Citrullus naudinianus are three important related wild species of watermelon in the genus Citrullus, and their morphological differences are clear, however, their chloroplast genome differences remain unknown. This study is the first to assemble, analyze, and publish the complete chloroplast genomes of C. ecirrhosus, C. rehmii, and C. naudinianus. A comparative analysis was then conducted among the complete chloroplast genomes of seven extant Citrullus species, and the results demonstrated that the average genome sizes of Citrullus is 157,005 bp, a total of 130-133 annotated genes were identified, including 8 rRNA, 37 tRNA and 85-88 protein-encoding genes. Their gene content, order, and genome structure were similar. However, noncoding regions were more divergent than coding regions, and rps16-trnQ was a hypervariable fragment. Thirty-four polymorphic SSRs, 1,271 SNPs and 234 INDELs were identified. Phylogenetic trees revealed a clear phylogenetic relationship of Citrullus species, and the developed molecular markers (SNPs and rps16-trnQ) could be used for taxonomy in Citrullus. Three genes (atpB, clpP1, and rpoC2) were identified to undergo selection and would promote the environmental adaptation of Citrullus.

The Complete Chloroplast Genomes of Gynostemma Reveal the Phylogenetic Relationships of Species within the Genus

Gynostemma is an important medicinal and food plant of the Cucurbitaceae family. The phylogenetic position of the genus Gynostemma in the Cucurbitaceae family has been determined by morphology and phylogenetics, but the evolutionary relationships within the genus Gynostemma remain to be explored. The chloroplast genomes of seven species of the genus Gynostemma were sequenced and annotated, of which the genomes of Gynostemma simplicifolium, Gynostemma guangxiense and Gynostemma laxum were sequenced and annotated for the first time. The chloroplast genomes ranged from 157,419 bp (Gynostemma compressum) to 157,840 bp (G. simplicifolium) in length, including 133 identical genes: 87 protein-coding genes, 37 tRNA genes, eight rRNA genes and one pseudogene. Phylogenetic analysis showed that the genus Gynostemma is divided into three primary taxonomic clusters, which differs from the traditional morphological classification of the genus Gynostemma into the subgenus Gynostemma and Trirostellum. The highly variable regions of atpH-atpL, rpl32-trnL, and ccsA-ndhD, the repeat unilts of AAG/CTT and ATC/ATG in simple sequence repeats (SSRs) and the length of overlapping regions between rps19 and inverted repeats(IRb) and between ycf1 and small single-copy (SSC) were found to be consistent with the phylogeny. Observations of fruit morphology of the genus Gynostemma revealed that transitional state species have independent morphological characteristics, such as oblate fruit and inferior ovaries. In conclusion, both molecular and morphological results showed consistency with those of phylogenetic analysis.

Comparative karyotype analysis of eight Cucurbitaceae crops using fluorochrome banding and 45S rDNA-FISH

To have an insight into the karyotype variation of eight Cucurbitaceae crops including Cucumissativus Linnaeus, 1753, Cucumismelo Linnaeus, 1753, Citrulluslanatus (Thunberg, 1794) Matsumura et Nakai, 1916, Benincasahispida (Thunberg, 1784) Cogniaux, 1881, Momordicacharantia Linnaeus, 1753, Luffacylindrica (Linnaeus, 1753) Roemer, 1846, Lagenariasicerariavar.hispida (Thunberg, 1783) Hara, 1948 and Cucurbitamoschata Duchesne ex Poiret, 1819, well morphologically differentiated mitotic metaphase chromosomes were prepared using the enzymatic maceration and flame-drying method, and the chromosomal distribution of heterochromatin and 18S-5.8S-26S rRNA genes (45S rDNA) was investigated using sequential combined PI and DAPI (CPD) staining and fluorescence in situ hybridization (FISH) with 45S rDNA probe. Detailed karyotypes were established using the dataset of chromosome measurements, fluorochrome bands and rDNA FISH signals. Four karyotype asymmetry indices, CVCI, CVCL, MCA and Stebbins' category, were measured to elucidate the karyological relationships among species. All the species studied had symmetrical karyotypes composed of metacentric and submetacentric or only metacentric chromosomes, but their karyotype structure can be discriminated by the scatter plot of MCA vs. CVCL. The karyological relationships among these species revealed by PCoA based on x, 2n, TCL, MCA, CVCL and CVCI was basically in agreement with the phylogenetic relationships revealed by DNA sequences. CPD staining revealed all 45S rDNA sites in all species, (peri)centromeric GC-rich heterochromatin in C.sativus, C.melo, C.lanatus, M.charantia and L.cylindrica, terminal GC-rich heterochromatin in C.sativus. DAPI counterstaining after FISH revealed pericentromeric DAPI+ heterochromatin in C.moschata. rDNA FISH detected two 45S loci in five species and five 45S loci in three species. Among these 45S loci, most were located at the terminals of chromosome arms, and a few in the proximal regions. In C.sativus, individual chromosomes can be precisely distinguished by the CPD band and 45S rDNA signal patterns, providing an easy method for chromosome identification of cucumber. The genome differentiation among these species was discussed in terms of genome size, heterochromatin, 45S rDNA site, and karyotype asymmetry based on the data of this study and previous reports.

The complete plastome sequence of Momordica cochinchinensis (Cucurbitaceae)

Momordica cochinchinensis (Lour.) Spreng. is an important medicinal plant that is used to treat various diseases in South and Southeast Asia. In this study, the complete plastome of M. cochinchinensis was sequenced and found to exhibit a total length of 158,955 bp, with a large single copy (LSC) region of 87,924 bp and a small single copy (SSC) region of 18,479 bp, as well as with two inverted repeats (IRs) that were both 26,726 bp in length. In total, 129 genes were detected, comprising 86 protein-encoding genes, 8 ribosomal RNA (rRNA) genes, and 35 transfer RNA (tRNA) genes. Furthermore, the inferred phylogenetic tree confirmed that M. cochinchinensis belongs to the genus Momordica in the Cucurbitaceae family. The research results will be used for authenticating M. cochinchinensis plant materials and for analyzing the genetic diversity and phylogenetic relationships in Momordica.

Diploid chromosome-level reference genome and population genomic analyses provide insights into Gypenoside biosynthesis and demographic evolution of Gynostemma pentaphyllum (Cucurbitaceae)

Gynostemma pentaphyllum (Thunb.) Makino is a perennial creeping herbaceous plant in the family Cucurbitaceae, which has great medicinal value and commercial potential, but urgent conservation efforts are needed due to the gradual decreases and fragmented distribution of its wild populations. Here, we report the high-quality diploid chromosome-level genome of G. pentaphyllum obtained using a combination of next-generation sequencing short reads, Nanopore long reads, and Hi-C sequencing technologies. The genome is anchored to 11 pseudo-chromosomes with a total size of 608.95 Mb and 26 588 predicted genes. Comparative genomic analyses indicate that G. pentaphyllum is estimated to have diverged from Momordica charantia 60.7 million years ago, with no recent whole-genome duplication event. Genomic population analyses based on genotyping-by-sequencing and ecological niche analyses indicated low genetic diversity but a strong population structure within the species, which could classify 32 G. pentaphyllum populations into three geographical groups shaped jointly by geographic and climate factors. Furthermore, comparative transcriptome analyses showed that the genes encoding enzyme involved in gypenoside biosynthesis had higher expression levels in the leaves and tendrils. Overall, the findings obtained in this study provide an effective molecular basis for further studies of demographic genetics, ecological adaption, and systematic evolution in Cucurbitaceae species, as well as contributing to molecular breeding, and the biosynthesis and biotransformation of gypenoside.

Evolutionary Comparison of the Complete Chloroplast Genomes in Convallaria Species and Phylogenetic Study of Asparagaceae

The genus Convallaria (Asparagaceae) comprises three herbaceous perennial species that are widely distributed in the understory of temperate deciduous forests in the Northern Hemisphere. Although Convallaria species have high medicinal and horticultural values, studies related to the phylogenetic analysis of this genus are few. In the present study, we assembled and reported five complete chloroplast (cp) sequences of three Convallaria species (two of C. keiskei Miq., two of C. majalis L., and one of C. montana Raf.) using Illumina paired-end sequencing data. The cp genomes were highly similar in overall size (161,365–162,972 bp), and all consisted of a pair of inverted repeats (IR) regions (29,140–29,486 bp) separated by a large single-copy (LSC) (85,183–85,521 bp) and a small single-copy (SSC) region (17,877–18,502 bp). Each cp genome contained the same 113 unique genes, including 78 protein-coding genes, 30 transfer RNA genes, and 4 ribosomal RNA genes. Gene content, gene order, AT content and IR/SC boundary structure were nearly identical among all of the Convallaria cp genomes. However, their lengths varied due to contraction/expansion at the IR/LSC borders. Simple sequence repeat (SSR) analyses indicated that the richest SSRs are A/T mononucleotides. Three highly variable regions (petA-psbJ, psbI-trnS and ccsA-ndhD) were identified as valuable molecular markers. Phylogenetic analysis of the family Asparagaceae using 48 cp genome sequences supported the monophyly of Convallaria, which formed a sister clade to the genus Rohdea. Our study provides a robust phylogeny of the Asparagaceae family. The complete cp genome sequences will contribute to further studies in the molecular identification, genetic diversity, and phylogeny of Convallaria.

Comparative Analysis of the Complete Chloroplast Genomes in Allium Section Bromatorrhiza Species (Amaryllidaceae): Phylogenetic Relationship and Adaptive Evolution

With the development of molecular sequencing approaches, many taxonomic and phylogenetic problems of the genus Allium L. have been solved; however, the phylogenetic relationships of some subgenera or sections, such as section Bromatorrhiza, remain unresolved, which has greatly impeded our full understanding of the species relationships among the major clades of Allium. In this study, the complete chloroplast (cp) genomes of nine species in the Allium sect. Bromatorrhiza were determined using the Illumina paired-end sequencing, the NOVOPlasty de novo assembly strategy, and the PGA annotation method. The results showed that the cp genome exhibited high conservation and revealed a typical circular tetrad structure. Among the sect. Bromatorrhiza species, the gene content, SSRs, codon usage, and RNA editing site were similar. The genome structure and IR regions’ fluctuation were investigated while genes, CDSs, and non-coding regions were extracted for phylogeny reconstruction. Evolutionary rates (Ka/Ks values) were calculated, and positive selection analysis was further performed using the branch-site model. Five hypervariable regions were identified as candidate molecular markers for species authentication. A clear relationship among the sect. Bromatorrhiza species were detected based on concatenated genes and CDSs, respectively, which suggested that sect. Bromatorrhiza is monophyly. In addition, there were three genes with higher Ka/Ks values (rps2, ycf1, and ycf2), and four genes (rpoC2, atpF, atpI, and rpl14) were further revealed to own positive selected sites. These results provide new insights into the plastome component, phylogeny, and evolution of Allium species.

The complete chloroplast genome of critically endangered Chimonobambusa hirtinoda (Poaceae: Chimonobambusa) and phylogenetic analysis

Chimonobambusa hirtinoda, a threatened species, is only naturally distributed in Doupeng Mountain, Duyun, Guizhou, China. Next-generation sequencing (NGS) is used to obtain the complete chloroplast (cp) genome sequence of C. hirtinoda. The sequence was assembled and analyzed for phylogenetic and evolutionary studies. Additionally, we compared the cp genome of C. hirtinoda with previously published Chimonobambusa species. The cp genome of C. hirtinoda has a total length of 139, 561 bp and 38.90% GC content. This genome included a large single -copy (LSC) region of 83, 166 bp, a small single-copy (SSC) region of 20, 811 bp and a pair of inverted repeats of 21,792 bp each. We discovered 130 genes in the cp genome, including 85 protein-coding genes, 37 tRNA, and 8 rRNA genes. A total of 48 simple sequence repeats (SSRs) were detected. The A/U preference of the third nucleotide in the cp genome of C. hirtinoda was obtained by measuring the codon usage frequency of amino acids. Furthermore, phylogenetic analysis using complete cp sequences and matK gene revealed a genetic relationship within the Chimonobambusa genus. This study reported the chloroplast genome of the C. hirtinoda.

Complete chloroplast genome sequences of three aroideae species (Araceae): lights into selective pressure, marker development and phylogenetic relationships

Background

Colocasia gigantea, Caladium bicolor and Xanthosoma sagittifolium are three worldwide famous ornamental and/or vegetable plants in the Araceae family, these species in the subfamily Aroideae are phylogenetically perplexing due to shared interspecific morphological traits and variation.

Result

This study, for the first time ever, assembled and analyzed complete chloroplast genomes of C. gigantea, C. bicolor and X. sagittifolium with genome sizes of 165,906 bp, 153,149 bp and 165,169 bp in length, respectively. The genomes were composed of conserved quadripartite circular structures with a total of 131 annotated genes, including 8 rRNA, 37 tRNA and 86 protein-coding genes. A comparison within Aroideae showed seven protein-coding genes (accD, ndhF, ndhK, rbcL, rpoC1, rpoC2 and matK) linked to environmental adaptation. Phylogenetic analysis confirmed a close relationship of C. gigantea with C. esculenta and S. colocasiifolia, and the C. bicolor with X. sagittifolium. Furthermore, three DNA barcodes (atpH-atpI + psaC-ndhE, atpH-atpI + trnS-trnG, atpH-atpI + psaC-ndhE + trnS-trnG) harbored highly variable regions to distinguish species in Aroideae subfamily.

Conclusion

These results would be beneficial for species identification, phylogenetic relationship, genetic diversity, and potential of germplasm resources in Aroideae.

Supplementary Information

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

Characterization of the complete chloroplast genome sequence of Gynostemma microspermum (Cucurbitaceae)

Gynostemma microspermum C. Y. Wu et S. K. Chen is an endemic creeping herbaceous species mainly distributed in dense forests on limestone in northwestern China. Here, the complete chloroplast genome sequence of G. microspermum was obtained by Illumina pair-end sequencing. The circular complete chloroplast genome of G. microspermum is 158,692 bp in length and contains a large single copy region (87,452 bp), a small single copy region (19,068 bp) and two short inverted repeat regions (26,086 bp). The genome sequence encodes 133 genes including 87 protein-coding genes, 37 transfer RNA genes, 8 ribosomal RNA genes and 1 pseudogene. The maximum likelihood (ML) phylogeny estimation shows that G. microspermum is sister to all other analyzed species of the genus Gynostemma with high bootstrap support.

Plastid Phylogenomic Data Offers Novel Insights Into the Taxonomic Status of the Trichosanthes kirilowii Complex (Cucurbitaceae) in South Korea

Trichosanthes is a genus in Cucurbitaceae comprising 90-100 species. Trichosanthes species are valuable as herbaceous medicinal ingredients. The fruits, seeds, and roots of species such as T. kirilowii and T. rosthornii are used in Korean traditional herbal medicines. T. rosthornii is only found in China, whereas in South Korea two varieties, T. kirilowii var. kirilowii and T. kirilowii var. japonica, are distributed. T. kirilowii var. kirilowii and T. kirilowii var. japonica have different fruit and leaf shapes but are recognized as belonging to the same species. Furthermore, although its members have herbal medicine applications, genomic information of the genus is still limited. The broad goals of this study were (i) to evaluate the taxonomy of Trichosanthes using plastid phylogenomic data and (ii) provide molecular markers specific for T. kirilowii var. kirilowii and T. kirilowii var. japonica, as these have differences in their pharmacological effectiveness and thus should not be confused and adulterated. Comparison of five Trichosanthes plastid genomes revealed locally divergent regions, mainly within intergenic spacer regions (trnT-UGU-trnL-UAA: marker name Tri, rrn4.5-rrn5: TRr, trnE-UUC-trnT-GGU: TRtt). Using these three markers as DNA-barcodes for important herbal medicine species in Trichosanthes, the identity of Trichosanthes material in commercial medicinal products in South Korea could be successfully determined. Phylogenetic analysis of the five Trichosanthes species revealed that the species are clustered within tribe Sicyoeae. T. kirilowii var. kirilowii and T. rosthornii formed a clade with T. kirilowii var. japonica as their sister group. As T. kirilowii in its current circumscription is paraphyletic and as the two varieties can be readily distinguished morphologically (e.g., in leaf shape), T. kirilowii var. japonica should be treated (again) as an independent species, T. japonica.

Uncovering dynamic evolution in the plastid genome of seven Ligusticum species provides insights into species discrimination and phylogenetic implications

Ligusticum L., one of the largest members in Apiaceae, encompasses medicinally important plants, the taxonomic statuses of which have been proved to be difficult to resolve. In the current study, the complete chloroplast genomes of seven crucial plants of the best-known herbs in Ligusticum were presented. The seven genomes ranged from 148,275 to 148,564 bp in length with a highly conserved gene content, gene order and genomic arrangement. A shared dramatic decrease in genome size resulted from a lineage-specific inverted repeat (IR) contraction, which could potentially be a promising diagnostic character for taxonomic investigation of Ligusticum, was discovered, without affecting the synonymous rate. Although a higher variability was uncovered in hotspot divergence regions that were unevenly distributed across the chloroplast genome, a concatenated strategy for rapid species identification was proposed because separate fragments inadequately provided variation for fine resolution. Phylogenetic inference using plastid genome-scale data produced a concordant topology receiving a robust support value, which revealed that L. chuanxiong had a closer relationship with L. jeholense than L. sinense, and L. sinense cv. Fuxiong had a closer relationship to L. sinense than L. chuanxiong, for the first time. Our results not only furnish concrete evidence for clarifying Ligusticum taxonomy but also provide a solid foundation for further pharmaphylogenetic investigation.

Plastomes of eight Ligusticum species: characterization, genome evolution, and phylogenetic relationships

BACKGROUND

The genus Ligusticum consists of approximately 60 species distributed in the Northern Hemisphere. It is one of the most taxonomically difficult taxa within Apiaceae, largely due to the varied morphological characteristics. To investigate the plastome evolution and phylogenetic relationships of Ligusticum, we determined the complete plastome sequences of eight Ligusticum species using a de novo assembly approach.

RESULTS

Through a comprehensive comparative analysis, we found that the eight plastomes were similar in terms of repeat sequence, SSR, codon usage, and RNA editing site. However, compared with the other seven species, L. delavayi exhibited striking differences in genome size, gene number, IR/SC borders, and sequence identity. Most of the genes remained under the purifying selection, whereas four genes showed relaxed selection, namely ccsA, rpoA, ycf1, and ycf2. Non-monophyly of Ligusticum species was inferred from the plastomes and internal transcribed spacer (ITS) sequences phylogenetic analyses.

CONCLUSION

The plastome tree and ITS tree produced incongruent tree topologies, which may be attributed to the hybridization and incomplete lineage sorting. Our study highlighted the advantage of plastome with mass informative sites in resolving phylogenetic relationships. Moreover, combined with the previous studies, we considered that the current taxonomy system of Ligusticum needs to be improved and revised. In summary, our study provides new insights into the plastome evolution, phylogeny, and taxonomy of Ligusticum species.

Chloroplast Genome Sequences and Comparative Analyses of Combretaceae Mangroves with Related Species

In the Combretaceae family, only two species of Lumnitzera and one species of Laguncularia belong to mangroves. Among them, Lumnitzera littorea (Jack) Voigt. is an endangered mangrove plant in China for the limited occurrence and seed abortion. In contrast, Lumnitzera racemosa Willd. is known as the most widespread mangrove plant in China. Laguncularia racemosa C. F. Gaertn., an exotic mangrove in China, has the fast growth and high adaptation ability. To better understand the phylogenetic positions of these mangroves in Combretaceae and in Myrtales and to provide information for studies on evolutionary adaptation for intertidal habitat, the complete chloroplast (cp) genomes of Lu. racemosa and La. racemosa were sequenced. Furthermore, we present here the results from the assembly and annotation of the two cp genomes, which were further subjected to the comparative analysis with Lu. littorea cp genomes we published before and other eleven closely related species within Myrtales. The chloroplast genomes of the three Combretaceae mangrove species: Lu. littorea, Lu. racemosa, and La. racemosa are 159,687 bp, 159,473 bp, and 158,311 bp in size. All three cp genomes host 130 genes including 85 protein-coding genes, 37 tRNAs, and 4 rRNAs. A comparative analysis of those three genomes revealed the high similarity of genes in coding-regions and conserved gene order in the IR and LSC/SSC regions. The differences between Lumnitzera and Laguncularia cp genomes are the locations of rps19 and rpl2 genes in the IR/SC boundary regions. Investigating the effects of selection events on shared protein-coding genes showed a relaxed selection had acted on the ycf2, ycf1, and matK genes of Combretaceae mangroves compared to the nonmangrove species Eucalyptus aromaphloia. The phylogenetic analysis based on the whole chloroplast genome sequence with one outgroup species strongly supported three Combretaceae mangroves together with other two Combretaceae species formed a cluster in Combretaceae. This study is the first report on the comparative analysis of three Combretaceae mangrove chloroplast genomes, which will provide the significant information for understanding photosynthesis and evolution in Combretaceae mangrove plants.

Complete chloroplast genomes of two Siraitia Merrill species: Comparative analysis, positive selection and novel molecular marker development

Siraitia grosvenorii fruit, known as Luo-Han-Guo, has been used as a traditional Chinese medicine for many years, and mogrosides are its primary active ingredients. Unfortunately, Siraitia siamensis, its wild relative, might be misused due to its indistinguishable appearance, not only threatening the reliability of the medication but also partly exacerbating wild resource scarcity. Therefore, high-resolution genetic markers must be developed to discriminate between these species. Here, the complete chloroplast genomes of S. grosvenorii and S. siamensis were assembled and analyzed for the first time; they were 158,757 and 159,190 bp in length, respectively, and possessed conserved quadripartite circular structures. Both contained 134 annotated genes, including 8 rRNA, 37 tRNA and 89 protein-coding genes. Twenty divergences (Pi > 0.03) were found in the intergenic regions. Nine protein-coding genes, accD, atpA, atpE, atpF, clpP, ndhF, psbH, rbcL, and rpoC2, underwent selection within Cucurbitaceae. Phylogenetic relationship analysis indicated that these two species originated from the same ancestor. Finally, four pairs of molecular markers were developed to distinguish the two species. The results of this study will be beneficial for taxonomic research, identification and conservation of Siraitia Merrill wild resources in the future.