Complete sequencing of five araliaceae chloroplast genomes and the phylogenetic implications

Characterization and phylogenetic analysis of the complete chloroplast genome of Carpesium lipskyi (Asteraceae, Inuleae)

The species of Carpesium lipskyi C.Winkl. 1998 is an important traditional Chinese medicine in China. In this study, the complete chloroplast (cp) genome of C. lipskyi was determined and analyzed. The result showed that the complete cp genome of C. lipskyi was 151,244 bp in length, consisting of a large single-copy (LSC) region of 82,908 bp, a small single-copy (SSC) region of 18,430 bp, and a pair of inverted repeats (IRs) of 24,953 bp. The overall GC content of the C. lipskyi is 37.68%. The species of C. lipskyi possessed 127 genes, including 83 protein-coding genes, 36 transfer RNA genes, and eight ribosomal RNA genes. The present study found that Inula is sister groups with the closest genetic relationship. The obtained knowledge could provide useful information for future phylogenetic, taxonomic, and evolutionary studies on Inuleae.

Characterization and phylogenetic analysis of the chloroplast genome of Duhaldea cappa (Buch.-Ham. ex D.Don) Pruski & Anderb. (Asteraceae)

Duhaldea cappa, a valuable medicinal plant of genus Duhaldea in the tribe Inuleae, is predominantly found in China, Bhutan, India, Malaysia, Nepal, Pakistan, Thailand, and Vietnam. However, the genomic studies of Duhaldea cappa are limited. In this study, we successfully sequenced and assembled the complete chloroplast genome of Duhaldea cappa. The chloroplast genome is 150,819 bp in length with a 37.73% GC content. The chloroplast genome has a quadripartite structure, consisting of a large single-copy region of 82,731 bp, a small single-copy region of 18,168 bp, and a pair of inverted repeat sequences of 24,960 bp. The genome contains 133 genes. Among these genes, there are 88 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The phylogeny reconstructed from data of the complete chloroplast genome indicated that Duhaldea cappa is closely related to Pluchea indica in the tribe Inuleae. Analyzing and reporting the chloroplast genome of Duhaldea cappa will establish a solid theoretical and data foundation for the efficient development, conservation, and utilization of this plant species.

Evolution of the Araliaceae family involved rapid diversification of the Asian Palmate group and Hydrocotyle specific mutational pressure

The Araliaceae contain many valuable species in medicinal and industrial aspects. We performed intensive phylogenomics using the plastid genome (plastome) and 45S nuclear ribosomal DNA sequences. A total of 66 plastome sequences were used, 13 of which were newly assembled in this study, 12 from new sequences, and one from existing data. While Araliaceae plastomes showed conserved genome structure, phylogenetic reconstructions based on four different plastome datasets revealed phylogenetic discordance within the Asian Palmate group. The divergence time estimation revealed that splits in two Araliaceae subfamilies and the clades exhibiting phylogenetic discordances in the Asian Palmate group occurred at two climatic optima, suggesting that global warming events triggered species divergence, particularly the rapid diversification of the Asian Palmate group during the Middle Miocene. Nucleotide substitution analyses indicated that the Hydrocotyloideae plastomes have undergone accelerated AT-biased mutations (C-to-T transitions) compared with the Aralioideae plastomes, and the acceleration may occur in their mitochondrial and nuclear genomes as well. This implies that members of the genus Hydrocotyle, the only aquatic plants in the Araliaceae, have experienced a distinct evolutionary history from the other species. We also discussed the intercontinental disjunction in the genus Panax and proposed a hypothesis to complement the previously proposed hypothesis. Our results provide the evolutionary trajectory of Araliaceae and advance our current understanding of the evolution of Araliaceae species.

New Insights into Phylogenetic Relationship of Hydrocotyle (Araliaceae) Based on Plastid Genomes

Hydrocotyle, belonging to the Hydrocotyloideae of Araliaceae, consists of 95 perennial and 35 annual species. Due to the lack of stable diagnostic morphological characteristics and high-resolution molecular markers, the phylogenetic relationships of Hydrocotyle need to be further investigated. In this study, we newly sequenced and assembled 13 whole plastid genomes of Hydrocotyle and performed comparative plastid genomic analyses with four previously published Hydrocotyle plastomes and phylogenomic analyses within Araliaceae. The plastid genomes of Hydrocotyle exhibited typical quadripartite structures with lengths from 152,659 bp to 153,669 bp, comprising a large single-copy (LSC) region (83,958-84,792 bp), a small single-copy (SSC) region (18,585-18,768 bp), and a pair of inverted repeats (IRs) (25,058-25,145 bp). Each plastome encoded 113 unique genes, containing 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. Comparative analyses showed that the IR boundaries of Hydrocotyle plastomes were highly similar, and the coding and IR regions exhibited more conserved than non-coding and single-copy (SC) regions. A total of 2932 simple sequence repeats and 520 long sequence repeats were identified, with specificity in the number and distribution of repeat sequences. Six hypervariable regions were screened from the SC region, including four intergenic spacers (IGS) (ycf3-trnS, trnS-rps4, petA-psbJ, and ndhF-rpl32) and two coding genes (rpl16 and ycf1). Three protein-coding genes (atpE, rpl16, and ycf2) were subjected to positive selection only in a few species, implying that most protein-coding genes were relatively conserved during the plastid evolutionary process. Plastid phylogenomic analyses supported the treatment of Hydrocotyle from Apiaceae to Araliaceae, and topologies with a high resolution indicated that plastome data can be further used in the comprehensive phylogenetic research of Hydrocotyle. The diagnostic characteristics currently used in Hydrocotyle may not accurately reflect the phylogenetic relationships of this genus, and new taxonomic characteristics may need to be evaluated and selected in combination with more comprehensive molecular phylogenetic results.

Characterization and phylogenetic analysis of the complete chloroplast genome of Saussurea sagittifolia (Asteraceae, Cardueae)

The species of Saussurea sagittifolia Y. S. Chen & S. R. Yi belongs to the family Asteraceae (Cardueae). The complete chloroplast genome of S. sagittifolia was assembled and annotated for the first time in this study. The complete chloroplast genome of S. sagittifolia was 152,535 bp, including a large single-copy (LSC) region of 83,511 bp, a small single-copy (SSC) region of 18,632 bp, and a pair of inverted repeats (IRs) of 25,196 bp. The overall GC content of the chloroplast genome was 37.7%. The chloroplast genome encoded 131 genes, including 87 protein-coding genes, 36 tRNA genes, and eight rRNA genes. Phylogenetic analysis based on complete chloroplast sequences revealed that it related closely to Saussurea medusa.

Comparative and phylogenetic analysis of Asparagus meioclados Levl. and Asparagus munitus Wang et S. C. Chen plastomes and utility of plastomes mutational hotspots

Tiandong is a vital traditional Chinese herbal medicine. It is derived from the tuber root of the Asparagus cochinchinensis according to the Pharmacopoeia of the people's republic of China (2020 Edition). On account of the similar morphology, Asparagus meioclados and Asparagus munitus were used as Tian-Dong in southwest China. Chloroplast (cp) genomes are highly active genetic components of plants and play an extremely important role in improving the efficiency of the identification of plant species. To differentiate the medicinal plants belonging to the genus Asparagus, we sequenced and analyzed the complete plastomes (plastid genomes) of A. meioclados and A. munitus and obtained two plastomes whose length changed to 156,515 bp and 156,381 bp, respectively. A total of 111 unique genes have been detected in plastome, which included 78 protein-coding genes, 29 tRNA genes and 4 rRNA genes. In plastomes of A. meioclados and A. munitus, 14,685 and 14,987 codons were detected, among which 9942 and 10,207 had the relative synonymous codon usage (RSCU) values higher than 1, respectively. A. meioclados and A. munitus have 26 SSRs patterns, among which A. meioclados was 25 and A. munitus 21. The average Ka/Ks value was 0.36, and positive selection was detected in genes of the photosynthetic system (ndhF and rbcL) in Asparagus species. To perform the comparative analysis of plastomes, the two newly sequenced plastomes of the A. meioclados and A. munitus species were compared with that of A. cochinchinensis, and 12 hotspots, including 5 coding regions and 7 inter-genomic regions, were identified. Based on the whole plastome of Asparagus, 2 divergent hotspots (accD and rpl32-trnL-UAG) and 1 international barcode fragment (rbcL) were screened, which may be used as particular molecular markers for the identification of Asparagus species. In addition, we determined the phylogenetic relationship between A. meioclados and A. munitus in the genus Asparagus. This study enriches our knowledge of the molecular evolutionary relationships of the Asparagus genus and provides treasured data records for species identification, molecular breeding, and evolutionary analysis of this genus.

Comparative analyses and phylogenetic relationships of thirteen Pholidota species (Orchidaceae) inferred from complete chloroplast genomes

BACKGROUND

The orchid genus Pholidota Lindl. ex Hook. is economically important as some species has long been used in traditional medicine. However, the systematic status of the genus and intergeneric relationships inferred from previous molecular studies are unclear due to insufficient sampling and lack of informative sites. So far, only limited genomic information has been available. The taxonomy of Pholidota remains unresolved and somewhat controversial. In this study, the complete chloroplast (cp.) genomes of thirteen Pholidota species were sequenced and analyzed to gain insight into the phylogeny of Pholidota and mutation patterns in their cp. genomes.

RESULTS

All examined thirteen Pholidota cp. genomes exhibited typical quadripartite circular structures, with the size ranging from 158,786 to 159,781 bp. The annotation contained a total of 135 genes in each cp. genome, i.e., 89 protein-coding genes, 38 tRNA genes, and eight rRNA genes. The codon usage analysis indicated the preference of A/U-ending codons. Repeat sequence analysis identified 444 tandem repeats, 322 palindromic repeats and 189 dispersed repeats. A total of 525 SSRs, 13,834 SNPs and 8,630 InDels were detected. Six mutational hotspots were identified as potential molecular markers. These molecular markers and highly variable regions are expected to facilitate future genetic and genomic studies. Our phylogenetic analyses confirmed the polyphyletic status of the genus Pholidota, with species grouped into four main clades: Pholidota s.s. was resolved as the sister to a clade containing species of Coelogyne; the other two clades clustered together with species of Bulleyia and Panisea, respectively; species P. ventricosa was placed at the basal position, deviated from all other species.

CONCLUSION

This is the first study to comprehensively examine the genetic variations and systematically analyze the phylogeny and evolution of Pholidota based on plastid genomic data. These findings contribute to a better understanding of plastid genome evolution of Pholidota and provide new insights into the phylogeny of Pholidota and its closely related genera within the subtribe Coelogyninae. Our research has laid the foundation for future studies on the evolutionary mechanisms and classification of this economically and medicinally important genus.

Complete Chloroplast Genome Sequence of the Endemic and Endangered Plant Dendropanax oligodontus: Genome Structure, Comparative and Phylogenetic Analysis

Dendropanax oligodontus, which belongs to the family Araliaceae, is an endemic and endangered species of Hainan Island, China. It has potential economic and medicinal value owing to the presence of phenylpropanoids, flavonoids, triterpenoids, etc. The analysis of the structure and characteristics of the D. oligodontus chloroplast genome (cpDNA) is crucial for understanding the genetic and phylogenetic evolution of this species. In this study, the cpDNA of D. oligodontus was sequenced for the first time using next-generation sequencing methods, assembled, and annotated. We observed a circular quadripartite structure comprising a large single-copy region (86,440 bp), a small single-copy region (18,075 bp), and a pair of inverted repeat regions (25,944 bp). The total length of the cpDNA was 156,403 bp, and the GC% was 37.99%. We found that the D. oligodontus chloroplast genome comprised 131 genes, with 86 protein-coding genes, 8 rRNA genes, and 37 tRNAs. Furthermore, we identified 26,514 codons, 13 repetitive sequences, and 43 simple sequence repeat sites in the D. oligodontus cpDNA. The most common amino acid encoded was leucine, with a strong A/T preference at the third position of the codon. The prediction of RNA editing sites in the protein-coding genes indicated that RNA editing was observed in 19 genes with a total of 54 editing sites, all of which involved C-to-T transitions. Finally, the cpDNA of 11 species of the family Araliaceae were selected for comparative analysis. The sequences of the untranslated regions and coding regions among 11 species were highly conserved, and minor differences were observed in the length of the inverted repeat regions; therefore, the cpDNAs were relatively stable and consistent among these 11 species. The variable hotspots in the genome included clpP, ycf1, rnK-rps16, rps16-trnQ, atpH-atpI, trnE-trnT, psbM-trnD, ycf3-trnS, and rpl32-trnL, providing valuable molecular markers for species authentication and regions for inferring phylogenetic relationships among them, as well as for evolutionary studies. Evolutionary selection pressure analysis indicated that the atpF gene was strongly subjected to positive environmental selection. Phylogenetic analysis indicated that D. oligodontus and Dendropanax dentiger were the most closely related species within the genus, and D. oligodontus was closely related to the genera Kalopanax and Metapanax in the Araliaceae family. Overall, the cp genomes reported in this study will provide resources for studying the genetic diversity and conservation of the endangered plant D. oligodontus, as well as resolving phylogenetic relationships within the family.

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.

The entire chloroplast genome sequence of Asparagus cochinchinensis and genetic comparison to Asparagus species

Asparagus cochinchinensis is a traditional Chinese medicinal plant. The chloroplast (cp) genome study on A. cochinchinensis is poorly understood. In this research, we collected the data from the cp genome assembly and gene annotation of A. cochinchinensis, followed by further comparative analysis with six species in the genus Asparagus. The cp genome of A. cochinchinensis showed a circular quadripartite structure in the size of 157,095 bp, comprising a large single-copy (LSC), a small single-copy (SSC), and two inverted repeat (IR) regions. A total of 137 genes were annotated, consisting of 86 protein-coding genes, 8 ribosomal RNAs, 38 transfer RNAs, and 5 pseudo-genes. Forty scattered repetitive sequences and 247 simple sequence repeats loci were marked out. In addition, A/T-ending codons were shown to have a basis in the codon analysis. A cp genome comparative analysis revealed that a similar gene composition was detected in the IR and LSC/SSC regions with Asparagus species. Based on the complete cp genome sequence in Asparagaceae, the result showed that A. cochinchinensis was closely related to A. racemosus by phylogenetic analysis. Therefore, our study providing A. cochinchinensis genomic resources could effectively contribute to the phylogenetic analysis and molecular identification of the genus Asparagus.

The complete plastid genome of the endangered shrub Brassaiopsis angustifolia (Araliaceae): Comparative genetic and phylogenetic analysis

Brassaiopsis angustifolia K.M. Feng belongs to the family Araliaceae, and is an endangered shrub species in southwest China. Despite the importance of this species, the plastid genome has not been sequenced and analyzed. In this study, the complete plastid genome of B. angustifolia was sequenced, analyzed, and compared to the eight species in the Araliaceae family. Our study reveals that the complete plastid genome of B. angustifolia is 156,534 bp long, with an overall GC content of 37.9%. The chloroplast genome (cp) encodes 133 genes, including 88 protein-coding genes, 37 transfer RNA (tRNA) genes, and eight ribosomal RNA (rRNA) genes. All protein-coding genes consisted of 21,582 codons. Among the nine species of Araliaceae, simple sequence repeats (SSRs) and five large repeat sequences were identified with total numbers ranging from 37 to 46 and 66 to 78, respectively. Five highly divergent regions were successfully identified that could be used as potential genetic markers of Brassaiopsis and Asian Palmate group. Phylogenetic analysis of 47 plastomes, representing 19 genera of Araliaceae and two related families, was performed to reconstruct highly supported relationships for the Araliaceae, which highlight four well-supported clades of the Hydrocotyle group, Greater Raukaua group, Aralia-Panax group, and Asian Palmate group. The genus Brassaiopsis can be divided into four groups using internal transcribed spacer (ITS) data. The results indicate that plastome and ITS data can contribute to investigations of the taxonomy, and phylogeny of B. angustifolia. This study provides a theoretical basis for species identification and future biological research on resources of the genus Brassaiopsis.

Comparative chloroplast genome analyses of cultivated spinach and two wild progenitors shed light on the phylogenetic relationships and variation

Spinacia is a genus of important leafy vegetable crops worldwide and includes cultivated Spinacia oleracea and two wild progenitors, Spinacia turkestanica and Spinacia tetrandra. However, the chloroplast genomes of the two wild progenitors remain unpublished, limiting our knowledge of chloroplast genome evolution among these three Spinacia species. Here, we reported the complete chloroplast genomes of S. oleracea, S. turkestanica, and S. tetrandra obtained via Illumina sequencing. The three chloroplast genomes exhibited a typical quadripartite structure and were 150,739, 150,747, and 150,680 bp in size, respectively. Only three variants were identified between S. oleracea and S. turkestanica, whereas 690 variants were obtained between S. oleracea and S. tetrandra, strongly demonstrating the close relationship between S. turkestanica and S. oleracea. This was further supported by phylogenetic analysis. We reported a comprehensive variant dataset including 503 SNPs and 83 Indels using 85 Spinacia accessions containing 61 S. oleracea, 16 S. turkestanica, and eight S. tetrandra accessions. Thirteen S. oleracea accessions were derived through introgression from S. turkestanica that acts as the maternal parent. Together, these results provide a valuable resource for spinach breeding programs and improve our understanding of the phylogenetic relationships within Amaranthaceae.

Comparative Chloroplast Genome Analyses of the Winter-Blooming Eastern Asian Endemic Genus Chimonanthus (Calycanthaceae) With Implications For Its Phylogeny and Diversification

Chimonanthus of Calycanthaceae is a small endemic genus in China, with unusual winter-blooming sweet flowers widely cultivated for ornamentals and medicinal uses. The evolution of Chimonanthus plastomes and its phylogenetic relationships remain unresolved due to limited availability of genetic resources. Here, we report fully assembled and annotated chloroplast genomes of five Chimonanthus species. The chloroplast genomes of the genus (size range 153,010 – 153,299 bp) reveal high similarities in gene content, gene order, GC content, codon usage, amino acid frequency, simple sequence repeats, oligonucleotide repeats, synonymous and non-synonymous substitutions, and transition and transversion substitutions. Signatures of positive selection are detected in atpF and rpoB genes in C. campanulatus. The correlations among substitutions, InDels, and oligonucleotide repeats reveal weak to strong correlations in distantly related species at the intergeneric levels, and very weak to weak correlations among closely related Chimonanthus species. Chloroplast genomes are used to reconstruct a well-resolved phylogenetic tree, which supports the monophyly of Chimonanthus. Within Chimonanthus, C. praecox and C. campanulatus form one clade, while C. grammatus, C. salicifolius, C. zhejiangensis, and C. nitens constitute another clade. Chimonanthus nitens appears paraphyletic and is closely related to C. salicifolius and C. zhejiangensis, suggesting the need to reevaluate the species delimitation of C. nitens. Chimonanthus and Calycanthus diverged in mid-Oligocene; the radiation of extant Chimonanthus species was dated to the mid-Miocene, while C. grammatus diverged from other Chimonanthus species in the late Miocene. C. salicifolius, C. nitens(a), and C. zhejiangensis are inferred to have diverged in the Pleistocene of the Quaternary period, suggesting recent speciation of a relict lineage in the subtropical forest regions in eastern China. This study provides important insights into the chloroplast genome features and evolutionary history of Chimonanthus and family Calycanthaceae.

Phylogenomic and evolutionary dynamics of inverted repeats across Angelica plastomes

BACKGROUND

Angelica L. (family Apiaceae) is an economically important genus comprising ca. One hundred ten species. Angelica species are found on all continents of the Northern Hemisphere, and East Asia hosts the highest number of species. Morphological characters such as fruit anatomy, leaf morphology and subterranean structures of Angelica species show extreme diversity. Consequently, the taxonomic classification of Angelica species is complex and remains controversial, as the classifications proposed by previous studies based on morphological data and molecular data are highly discordant. In addition, the phylogenetic relationships of major clades in the Angelica group, particularly in the Angelica s. s. clade, remain unclear. Chloroplast (cp) genome sequences have been widely used in phylogenetic studies and for evaluating genetic diversity.

RESULTS

In this study, we sequenced and assembled 28 complete cp genomes from 22 species, two varieties and two cultivars of Angelica. Combined with 36 available cp genomes in GenBank from representative clades of the subfamily Apioideae, the characteristics and evolutionary patterns of Angelica cp genomes were studied, and the phylogenetic relationships of Angelica species were resolved. The Angelica cp genomes had the typical quadripartite structure including a pair of inverted repeats (IRs: 5836-34,706 bp) separated by a large single-copy region (LSC: 76,657-103,161 bp) and a small single-copy region (SSC: 17,433-21,794 bp). Extensive expansion and contraction of the IR region were observed among cp genomes of Angelica species, and the pattern of the diversification of cp genomes showed high consistency with the phylogenetic placement of Angelica species. Species of Angelica were grouped into two major clades, with most species grouped in the Angelica group and A. omeiensis and A. sinensis grouped in the Sinodielsia with Ligusticum tenuissimum.

CONCLUSIONS

Our results further demonstrate the power of plastid phylogenomics in enhancing the phylogenetic reconstructions of complex genera and provide new insights into plastome evolution across Angelica L.

Comparative and Phylogenetic Analysis of Complete Chloroplast Genomes in Eragrostideae (Chloridoideae, Poaceae)

Eragrostideae Stapf, the second-largest tribe in Chloridoideae (Poaceae), is a taxonomically complex tribe. In this study, chloroplast genomes of 13 Eragrostideae species were newly sequenced and used to resolve the phylogenetic relationships within Eragrostideae. Including seven reported chloroplast genomes from Eragrostideae, the genome structure, number and type of genes, codon usage, and repeat sequences of 20 Eragrostideae species were analyzed. The length of these chloroplast genomes varied from 130,773 bp to 135,322 bp. These chloroplast genomes showed a typical quadripartite structure, including a large single-copy region (77,993–80,643 bp), a small single-copy region (12,410–12,668 bp), and a pair of inverted repeats region (19,394–21,074 bp). There were, in total, 129–133 genes annotated in the genome, including 83–87 protein-coding genes, eight rRNA genes, and 38 tRNA genes. Forward and palindromic repeats were the most common repeat types. In total, 10 hypervariable regions (rpl22, rpoA, ndhF, matK, trnG–UCC-trnT–GGU, ndhF–rpl32, ycf4–cemA, rpl32–trnL–UAG, trnG–GCC–trnfM–CAU, and ccsA–ndhD) were found, which can be used as candidate molecular markers for Eragrostideae. Phylogenomic studies concluded that Enneapogon diverged first, and Eragrostis including Harpachne is the sister to Uniola. Furthermore, Harpachne harpachnoides is considered as a species of Eragrostis based on morphological and molecular evidence. In addition, the interspecies relationships within Eragrostis are resolved based on complete chloroplast genomes. This study provides useful chloroplast genomic information for further phylogenetic analysis of Eragrostideae.

Adaptation Evolution and Phylogenetic Analyses of Species in Chinese Allium Section Pallasia and Related Species Based on Complete Chloroplast Genome Sequences

The section Pallasia is one of the components of the genus Allium subgenus Allium (Amaryllidaceae), and species relationship in this section is still not resolved very well, which hinders further evolutionary and adaptive studies. Here, the complete chloroplast genomes of five sect. Pallasia species were reported, and a comparative analysis was performed with other three related Allium species. The genome size of the eight species ranged from 151,672 bp to 153,339 bp in length, GC content changed from 36.7% to 36.8%, and 130 genes (except Allium pallasii), 37 tRNA, and 8 rRNA were identified in each genome. By analyzing the IR/LSC and IR/SSC boundary, A. pallasii exhibited differences compared with other seven species. Phylogenetic analysis achieved high supports in each branch, seven of the eight Allium species cluster into a group, and A. pallasii exhibit a close relationship with A. obliquum. Higher pairwise Ka/Ks ratios were found in A. schoenoprasoides compared to A. caeruleum and A. macrostemon while a lower value of Ka/Ks ratios was detected between A. caeruleum and A. macrostemon. This study will be a great contribution to the future phylogenetic and adaptive research in Allium.

Reconstructing the Complex Evolutionary History of the Papuasian Schefflera Radiation Through Herbariomics

With its large proportion of endemic taxa, complex geological past, and location at the confluence of the highly diverse Malesian and Australian floristic regions, Papuasia - the floristic region comprising the Bismarck Archipelago, New Guinea, and the Solomon Islands - represents an ideal natural experiment in plant biogeography. However, scattered knowledge of its flora and limited representation in herbaria have hindered our understanding of the drivers of its diversity. Focusing on the woody angiosperm genus Schefflera (Araliaceae), we ask whether its morphologically defined infrageneric groupings are monophyletic, when these lineages diverged, and where (within Papuasia or elsewhere) they diversified. To address these questions, we use a high-throughput sequencing approach (Hyb-Seq) which combines target capture (with an angiosperm-wide bait kit targeting 353 single-copy nuclear loci) and genome shotgun sequencing (which allows retrieval of regions in high-copy number, e.g., organellar DNA) of historical herbarium collections. To reconstruct the evolutionary history of the genus we reconstruct molecular phylogenies with Bayesian inference, maximum likelihood, and pseudo-coalescent approaches, and co-estimate divergence times and ancestral areas in a Bayesian framework. We find strong support for most infrageneric morphological groupings, as currently circumscribed, and we show the efficacy of the Angiosperms-353 probe kit in resolving both deep and shallow phylogenetic relationships. We infer a sequence of colonization to explain the present-day distribution of Schefflera in Papuasia: from the Sunda Shelf, Schefflera arrived to the Woodlark plate (present-day eastern New Guinea) in the late Oligocene (when most of New Guinea was submerged) and, subsequently (throughout the Miocene), it migrated westwards (to the Maoke and Bird's Head Plates and thereon) and further diversified, in agreement with previous reconstructions.

New Insights Into the Plastome Evolution of the Millettioid/Phaseoloid Clade (Papilionoideae, Leguminosae)

The Millettioid/Phaseoloid (MP) clade from the subfamily Papilionoideae (Leguminosae) consists of six tribes and ca. 3,000 species. Previous studies have revealed some plastome structural variations (PSVs) within this clade. However, many deep evolutionary relationships within the clade remain unresolved. Due to limited taxon sampling and few genetic markers in previous studies, our understanding of the evolutionary history of this clade is limited. To address this issue, we sampled 43 plastomes (35 newly sequenced) representing all the six tribes of the MP clade to examine genomic structural variations and phylogenetic relationships. Plastomes of the species from the MP clade were typically quadripartite (size ranged from 140,029 to 160,040 bp) and contained 109-111 unique genes. We revealed four independent gene losses (ndhF, psbI, rps16, and trnS-GCU), multiple IR-SC boundary shifts, and six inversions in the tribes Desmodieae, Millettieae, and Phaseoleae. Plastomes of the species from the MP clade have experienced significant variations which provide valuable information on the evolution of the clade. Plastid phylogenomic analyses using Maximum Likelihood and Bayesian methods yielded a well-resolved phylogeny at the tribal and generic levels within the MP clade. This result indicates that plastome data is useful and reliable data for resolving the evolutionary relationships of the MP clade. This study provides new insights into the phylogenetic relationships and PSVs within this clade.

The chloroplasts genomic analyses of Rosa laevigata, R. rugosa and R. canina

BACKGROUND

Many species of the genus Rosa have been used as ornamental plants and traditional medicines. However, industrial development of roses is hampered due to highly divergent characteristics.

METHODS

We analyzed the chloroplast (cp) genomes of Rosa laevigata, R. rugosa and R. canina, including the repeat sequences, inverted-repeat (IR) contractions and expansions, and mutation sites.

RESULTS

The size of the cp genome of R. laevigata, R. rugosa and R. canina was between 156 333 bp and 156 533 bp, and contained 113 genes (30 tRNA genes, 4 rRNA genes and 79 protein-coding genes). The regions with a higher degree of variation were screened out (trnH-GUU, trnS-GCU, trnG-GCC, psbA-trnH, trnC-GCA,petN, trnT-GGU, psbD, petA, psbJ, ndhF, rpl32,psaC and ndhE). Such higher-resolution loci lay the foundation of barcode-based identification of cp genomes in Rosa genus. A phylogenetic tree of the genus Rosa was reconstructed using the full sequences of the cp genome. These results were largely in accordance with the current taxonomic status of Rosa.

CONCLUSIONS

Our data: (i) reveal that cp genomes can be used for the identification and classification of Rosa species; (ii) can aid studies on molecular identification, genetic transformation, expression of secondary metabolic pathways and resistant proteins; (iii) can lay a theoretical foundation for the discovery of disease-resistance genes and cultivation of Rosa species.

Comprehensive comparative analysis of chloroplast genomes from seven Panax species and development of an authentication system based on species-unique single nucleotide polymorphism markers

BACKGROUND

Panax species are important herbal medicinal plants in the Araliaceae family. Recently, we reported the complete chloroplast genomes and 45S nuclear ribosomal DNA sequences from seven Panax species, two (P . quinqu e folius and P . trifolius) from North America and five (P . ginseng, P . notoginseng, P . japonicus, P . vietnamensis, and P . stipuleanatus) from Asia.

METHODS

We conducted phylogenetic analysis of these chloroplast sequences with 12 other Araliaceae species and comprehensive comparative analysis among the seven Panax whole chloroplast genomes.

RESULTS

We identified 1,128 single nucleotide polymorphisms (SNP) in coding gene sequences, distributed among 72 of the 79 protein-coding genes in the chloroplast genomes of the seven Panax species. The other seven genes (including psaJ, psbN, rpl23, psbF, psbL, rps18, and rps7) were identical among the Panax species. We also discovered that 12 large chloroplast genome fragments were transferred into the mitochondrial genome based on sharing of more than 90% sequence similarity. The total size of transferred fragments was 60,331 bp, corresponding to approximately 38.6% of chloroplast genome. We developed 18 SNP markers from the chloroplast genic coding sequence regions that were not similar to regions in the mitochondrial genome. These markers included two or three species-specific markers for each species and can be used to authenticate all the seven Panax species from the others.

CONCLUSION

The comparative analysis of chloroplast genomes from seven Panax species elucidated their genetic diversity and evolutionary relationships, and 18 species-specific markers were able to discriminate among these species, thereby furthering efforts to protect the ginseng industry from economically motivated adulteration.

Development of molecular markers for invasive alien plants in Korea: a case study of a toxic weed, Cenchrus longispinus L., based on next generation sequencing data

BACKGROUND

Genomic data play an important role in plant research because of its implications in studying genomic evolution, phylogeny, and developing molecular markers. Although the information of invasive alien plants was collected, the genomic data of those species have not been intensively studied.

METHODS

We employ the next generation sequencing and PCR methods to explore the genomic data as well as to develop and test the molecular markers.

RESULTS

In this study, we characterize the chloroplast genomes (cpDNA) of Cenchrus longispinus and C. echinatus, of which the lengths are 137,144 and 137,131 bp, respectively. These two newly sequenced genomes include 78 protein-coding genes, 30 tRNA, and four rRNA. There are 56 simple single repeats and 17 forward repeats in the chloroplast genome of C. longispinus. Most of the repeats locate in non-coding regions. However, repeats can be found in infA, ndhD, ndhH, ndhK, psbC, rpl22, rpoC2, rps14, trnA-UGC, trnC-GCA, trnF-GAA, trnQ-UUG, trnS-UGA, trnS-GCU, and ycf15. The phylogenomic analysis revealed the monophyly of Cenchrus but not Panicum species in tribe Paniceae. The single nucleotide polymorphism sites in atpB, matK, and ndhD were successfully used for developing molecular markers to distinguish C. longispinus and related taxa. The simple PCR protocol for using the newly developed molecular markers was also provided.

Characterization of 20 complete plastomes from the tribe Laureae (Lauraceae) and distribution of small inversions

Lindera Thunb. (Lauraceae) consists of approximately 100 species, mainly distributed in the temperate and tropical regions of East Asia. In this study, we report 20 new, complete plastome sequences including 17 Lindera species and three related species, Actinodaphne lancifolia, Litsea japonica and Sassafras tzumu. The complete plastomes of Lindera range from 152,502 bp (L. neesiana) to 154,314 bp (L. erythrocarpa) in length. Eleven small inversion (SI) sites are documented among the plastomes. Six of the 11 SI sites are newly reported and they locate in rpoB-trnC, psbC-trnS, petA-psbJ, rpoA and ycf2 regions. The distribution patterns of SIs are useful for species identification. An average of 83 simple sequence repeats (SSRs) were detected in each plastome. The mono-SSRs accounted for 72.7% of total SSRs, followed by di- (12.4%), tetra- (9.4%), tri- (4.2%), and penta-SSRs (1.3%). Of these SSRs, 64.6% were distributed in an intergenic spacer (IGS) region. In addition, 79.8% of the SSRs are located in a large single copy (LSC) region. In contrast, almost no SSRs are distributed in inverted repeat (IR) regions. The SSR loci are useful to identifying species but the phylogenetic value is low because the majority of them show autapomorphic status or highly homoplastic characteristics. The nucleotide diversity (Pi) values also indicated the conserved nature of the IR region compared to LSC and small single copy (SSC) regions. Five spacer regions with high Pi values, trnH-psbA, petA-psbJ and ndhF-rpl32, rpl32-trnL and Ψycf1-ndhF, have a potential use for the molecular identification study of Lindera and related species. Lindera species form a paraphyletic group in the plastome tree because of the inclusion of related genera such as Actinodaphne, Laurus, Litsea and Neolitsea. A former member of tribe Laureae, Sassafras, forms a clade with the tribe Cinnamomeae. The SIs do not affect the phylogenetic relationship of Laureae. This result indicated that ancient plastome captures may have contribute to the mixed intergeneric relationship of Laureae. Alternatively, the result may indicate that the morphological characters defined the genera of Lauraceae originated for several times.

The complete chloroplast genome of Dendrobium nobile, an endangered medicinal orchid from north-east India and its comparison with related Dendrobium species

The medicinal orchid genus Dendrobium belonging to the Orchidaceae family is a huge genus comprising about 800-1,500 species. To better illustrate the species status in the genus Dendrobium, a comparative analysis of 33 available chloroplast genomes retrieved from NCBI RefSeq database was compared with that of the first complete chloroplast genome of D. nobile from north-east India based on next-generation sequencing methods (Illumina HiSeq 2500-PE150). Our results provide comparative chloroplast genomic information for taxonomical identification, alignment-free phylogenomic inference and other statistical features of Dendrobium plastomes, which can also provide valuable information on their mutational events and sequence divergence.

Comparative analysis of the complete chloroplast genome among Prunus mume, P. armeniaca, and P. salicina

Prunus mume Sieb. et Zucc., P. armeniaca L., and P. salicina L. are economically important fruit trees in temperate regions. These species are taxonomically perplexing because of shared interspecific morphological traits and variation, which are mainly attributed to hybridization. The chloroplast is cytoplasmically inherited and often used for evolutionary studies. We sequenced the complete chloroplast genomes of P. mume, P. armeniaca, and P. salicina using Illumina sequencing followed by de novo assembly. The three chloroplast genomes exhibit a typical quadripartite structure with conserved genome arrangement, structure, and moderate divergence. The lengths of the genomes are 157,815, 157,797, and 157,916 bp, respectively. The length of the large single-copy region (LSC) region is 86,113, 86,283, and 86,122 bp, and the length of the SSC region is 18,916, 18,734, and 19,028 bp; the IR region is 26,393, 26,390, and 26,383 bp, respectively. Each of the three chloroplast genomes encodes 133 genes, including 94 protein-coding, 31 tRNA, and eight rRNA genes. Differential gene analysis for the three species revealed that trnY-ATA is a unique gene in P. armeniaca; in contrast, the gene trnI-TAT is only present in P. mume and P. salicina, though the position of the gene in these chloroplast genomes differs. Further comparative analysis of the complete chloroplast genome sequences revealed that the ORF genes and the sequences of linked regions rps16 and atpA, atpH and atpI, trnc-GCA and psbD, ycf3 and atpB, and rpL32 and ndhD are significantly different and may be used as molecular markers in taxonomic studies. Phylogenetic evolution analysis of the three species suggests that P. mume has a closer genetic relationship to P. armeniaca than to P. salicina.

Analyzing and Characterizing the Chloroplast Genome of Salix wilsonii

Salix wilsonii is an important ornamental willow tree widely distributed in China. In this study, an integrated circular chloroplast genome was reconstructed for S. wilsonii based on the chloroplast reads screened from the whole-genome sequencing data generated with the PacBio RSII platform. The obtained pseudomolecule was 155,750 bp long and had a typical quadripartite structure, comprising a large single copy region (LSC, 84,638 bp) and a small single copy region (SSC, 16,282 bp) separated by two inverted repeat regions (IR, 27,415 bp). The S. wilsonii chloroplast genome encoded 115 unique genes, including four rRNA genes, 30 tRNA genes, 78 protein-coding genes, and three pseudogenes. Repetitive sequence analysis identified 32 tandem repeats, 22 forward repeats, two reverse repeats, and five palindromic repeats. Additionally, a total of 118 perfect microsatellites were detected, with mononucleotide repeats being the most common (89.83%). By comparing the S. wilsonii chloroplast genome with those of other rosid plant species, significant contractions or expansions were identified at the IR-LSC/SSC borders. Phylogenetic analysis of 17 willow species confirmed that S. wilsonii was most closely related to S. chaenomeloides and revealed the monophyly of the genus Salix. The complete S. wilsonii chloroplast genome provides an additional sequence-based resource for studying the evolution of organelle genomes in woody plants.

Sequencing and Structural Analysis of the Complete Chloroplast Genome of the Medicinal Plant Lycium chinense Mill

Lycium chinense Mill, an important Chinese herbal medicine, is widely used as a dietary supplement and food. Here the chloroplast (CP) genome of L. chinense was sequenced and analyzed, revealing a size of 155,756 bp and with a 37.8% GC content. The L. chinense CP genome comprises a large single copy region (LSC) of 86,595 bp and a small single copy region (SSC) of 18,209 bp, and two inverted repeat regions (IRa and IRb) of 25,476 bp separated by the single copy regions. The genome encodes 114 genes, 16 of which are duplicated. Most of the 85 protein-coding genes (CDS) had standard ATG start codons, while 3 genes including rps12, psbL and ndhD had abnormal start codons (ACT and ACG). In addition, a strong A/T bias was found in the majority of simple sequence repeats (SSRs) detected in the CP genome. Analysis of the phylogenetic relationships among 16 species revealed that L. chinense is a sister taxon to Lycium barbarum. Overall, the complete sequence and annotation of the L. chinense CP genome provides valuable genetic information to facilitate precise understanding of the taxonomy, species and phylogenetic evolution of the Solanaceae family.

The complete plastome of Panax stipuleanatus: Comparative and phylogenetic analyses of the genus Panax (Araliaceae)

Panax stipuleanatus (Araliaceae) is an endangered and medicinally important plant endemic to China. However, phylogenetic relationships within the genus Panax have remained unclear. In this study, we sequenced the complete plastome of P. stipuleanatus and included previously reported Panax plastomes to better understand the relationships between species and plastome evolution within the genus Panax. The plastome of P. stipuleanatus is 156,069 base pairs (bp) in length, consisting of a pair of inverted repeats (IRs, each 25,887 bp) that divide the plastome into a large single copy region (LSC, 86,126 bp) and a small single copy region (SSC, 8169 bp). The plastome contains 114 unigenes (80 protein-coding genes, 30 tRNA genes, and 4 rRNA genes). Comparative analyses indicated that the plastome gene content and order, as well as the expansion/contraction of the IR regions, are all highly conserved within Panax. No significant positive selection in the plastid protein-coding genes was observed across the eight Panax species, suggesting the Panax plastomes may have undergone a strong purifying selection. Our phylogenomic analyses resulted in a phylogeny with high resolution and supports for Panax. Nine protein-coding genes and 10 non-coding regions presented high sequence divergence, which could be useful for identifying different Panax species.

Complete Chloroplast Genome Sequence and Phylogenetic Analysis of Aster tataricus

We sequenced and analyzed the complete chloroplast genome of Aster tataricus (family Asteraceae), a Chinese herb used medicinally to relieve coughs and reduce sputum. The A. tataricus chloroplast genome was 152,992 bp in size, and harbored a pair of inverted repeat regions (IRa and IRb, each 24,850 bp) divided into a large single-copy (LSC, 84,698 bp) and a small single-copy (SSC, 18,250 bp) region. Our annotation revealed that the A. tataricus chloroplast genome contained 115 genes, including 81 protein-coding genes, 4 ribosomal RNA genes, and 30 transfer RNA genes. In addition, 70 simple sequence repeats (SSRs) were detected in the A. tataricus chloroplast genome, including mononucleotides (36), dinucleotides (1), trinucleotides (23), tetranucleotides (1), pentanucleotides (8), and hexanucleotides (1). Comparative chloroplast genome analysis of three Aster species indicated that a higher similarity was preserved in the IR regions than in the LSC and SSC regions, and that the differences in the degree of preservation were slighter between A. tataricus and A. altaicus than between A. tataricus and A. spathulifolius. Phylogenetic analysis revealed that A. tataricus was more closely related to A. altaicus than to A. spathulifolius. Our findings offer valuable information for future research on Aster species identification and selective breeding.

Positive Selection Driving Cytoplasmic Genome Evolution of the Medicinally Important Ginseng Plant Genus Panax

Panax L. (the ginseng genus) is a shade-demanding group within the family Araliaceae and all of its species are of crucial significance in traditional Chinese medicine. Phylogenetic and biogeographic analyses demonstrated that two rounds of whole genome duplications accompanying with geographic and ecological isolations promoted the diversification of Panax species. However, contributions of the cytoplasmic genomes to the adaptive evolution of Panax species remained largely uninvestigated. In this study, we sequenced the chloroplast and mitochondrial genomes of 11 accessions belonging to seven Panax species. Our results show that heterogeneity in nucleotide substitution rate is abundant in both of the two cytoplasmic genomes, with the mitochondrial genome possessing more variants at the total level but the chloroplast showing higher sequence polymorphisms at the genic regions. Genome-wide scanning of positive selection identified five and 12 genes from the chloroplast and mitochondrial genomes, respectively. Functional analyses further revealed that these selected genes play important roles in plant development, cellular metabolism and adaptation. We therefore conclude that positive selection might be one of the potential evolutionary forces that shaped nucleotide variation pattern of these Panax species. In particular, the mitochondrial genes evolved under stronger selective pressure compared to the chloroplast genes.

Phylogenomics and barcoding of Panax: toward the identification of ginseng species

BACKGROUND

The economic value of ginseng in the global medicinal plant trade is estimated to be in excess of US$2.1 billion. At the same time, the evolutionary placement of ginseng (Panax ginseng) and the complex evolutionary history of the genus is poorly understood despite several molecular phylogenetic studies. In this study, we use a full plastome phylogenomic framework to resolve relationships in Panax and to identify molecular markers for species discrimination.

RESULTS

We used high-throughput sequencing of MBD2-Fc fractionated Panax DNA to supplement publicly available plastid genomes to create a phylogeny based on fully assembled and annotated plastid genomes from 60 accessions of 8 species. The plastome phylogeny based on a 163 kbp matrix resolves the sister relationship of Panax ginseng with P. quinquefolius. The closely related species P. vietnamensis is supported as sister of P. japonicus. The plastome matrix also shows that the markers trnC-rps16, trnS-trnG, and trnE-trnM could be used for unambiguous molecular identification of all the represented species in the genus.

CONCLUSIONS

MBD2 depletion reduces the cost of plastome sequencing, which makes it a cost-effective alternative to Sanger sequencing based DNA barcoding for molecular identification. The plastome phylogeny provides a robust framework that can be used to study the evolution of morphological characters and biosynthesis pathways of ginsengosides for phylogenetic bioprospecting. Molecular identification of ginseng species is essential for authenticating ginseng in international trade and it provides an incentive for manufacturers to create authentic products with verified ingredients.

Integral Phylogenomic Approach over Ilex L. Species from Southern South America

The use of molecular markers with inadequate variation levels has resulted in poorly resolved phylogenetic relationships within Ilex. Focusing on southern South American and Asian species, we aimed at contributing informative plastid markers. Also, we intended to gain insights into the nature of morphological and physiological characters used to identify species. We obtained the chloroplast genomes of I.paraguariensis and I. dumosa, and combined these with all the congeneric plastomes currently available to accomplish interspecific comparisons and multilocus analyses. We selected seven introns and nine IGSs as variable non-coding markers that were used in phylogenomic analyses. Eight extra IGSs were proposed as candidate markers. Southern South American species formed one lineage, except for I. paraguariensis, I. dumosa and I. argentina, which occupied intermediate positions among sampled taxa; Euroasiatic species formed two lineages. Some concordant relationships were retrieved from nuclear sequence data. We also conducted integral analyses, involving a supernetwork of molecular data, and a simultaneous analysis of quantitative and qualitative morphological and phytochemical characters, together with molecular data. The total evidence tree was used to study the evolution of non-molecular data, evidencing fifteen non-ambiguous synapomorphic character states and consolidating the relationships among southern South American species. More South American representatives should be incorporated to elucidate their origin.

The complete chloroplast genome sequence of strawberry (Fragaria × ananassa Duch.) and comparison with related species of Rosaceae

Compared with other members of the family Rosaceae, the chloroplast genomes of Fragaria species exhibit low variation, and this situation has limited phylogenetic analyses; thus, complete chloroplast genome sequencing of Fragaria species is needed. In this study, we sequenced the complete chloroplast genome of F. × ananassa 'Benihoppe' using the Illumina HiSeq 2500-PE150 platform and then performed a combination of de novo assembly and reference-guided mapping of contigs to generate complete chloroplast genome sequences. The chloroplast genome exhibits a typical quadripartite structure with a pair of inverted repeats (IRs, 25,936 bp) separated by large (LSC, 85,531 bp) and small (SSC, 18,146 bp) single-copy (SC) regions. The length of the F. × ananassa 'Benihoppe' chloroplast genome is 155,549 bp, representing the smallest Fragaria chloroplast genome observed to date. The genome encodes 112 unique genes, comprising 78 protein-coding genes, 30 tRNA genes and four rRNA genes. Comparative analysis of the overall nucleotide sequence identity among ten complete chloroplast genomes confirmed that for both coding and non-coding regions in Rosaceae, SC regions exhibit higher sequence variation than IRs. The Ka/Ks ratio of most genes was less than 1, suggesting that most genes are under purifying selection. Moreover, the mVISTA results also showed a high degree of conservation in genome structure, gene order and gene content in Fragaria, particularly among three octoploid strawberries which were F. × ananassa 'Benihoppe', F. chiloensis (GP33) and F. virginiana (O477). However, when the sequences of the coding and non-coding regions of F. × ananassa 'Benihoppe' were compared in detail with those of F. chiloensis (GP33) and F. virginiana (O477), a number of SNPs and InDels were revealed by MEGA 7. Six non-coding regions (trnK-matK, trnS-trnG, atpF-atpH, trnC-petN, trnT-psbD and trnP-psaJ) with a percentage of variable sites greater than 1% and no less than five parsimony-informative sites were identified and may be useful for phylogenetic analysis of the genus Fragaria.

The complete chloroplast genome sequence of Aconitum coreanum and Aconitum carmichaelii and comparative analysis with other Aconitum species

Aconitum species (belonging to the Ranunculaceae) are well known herbaceous medicinal ingredients and have great economic value in Asian countries. However, there are still limited genomic resources available for Aconitum species. In this study, we sequenced the chloroplast (cp) genomes of two Aconitum species, A. coreanum and A. carmichaelii, using the MiSeq platform. The two Aconitum chloroplast genomes were 155,880 and 157,040 bp in length, respectively, and exhibited LSC and SSC regions separated by a pair of inverted repeat regions. Both cp genomes had 38% GC content and contained 131 unique functional genes including 86 protein-coding genes, eight ribosomal RNA genes, and 37 transfer RNA genes. The gene order, content, and orientation of the two Aconitum cp genomes exhibited the general structure of angiosperms, and were similar to those of other Aconitum species. Comparison of the cp genome structure and gene order with that of other Aconitum species revealed general contraction and expansion of the inverted repeat regions and single copy boundary regions. Divergent regions were also identified. In phylogenetic analysis, Aconitum species positon among the Ranunculaceae was determined with other family cp genomes in the Ranunculales. We obtained a barcoding target sequence in a divergent region, ndhC–trnV, and successfully developed a SCAR (sequence characterized amplified region) marker for discrimination of A. coreanum. Our results provide useful genetic information and a specific barcode for discrimination of Aconitum species.

Authentication Markers for Five Major Panax Species Developed via Comparative Analysis of Complete Chloroplast Genome Sequences

Ginseng represents a set of high-value medicinal plants of different species: Panax ginseng (Asian ginseng), Panax quinquefolius (American ginseng), Panax notoginseng (Chinese ginseng), Panax japonicus (Bamboo ginseng), and Panax vietnamensis (Vietnamese ginseng). Each species is pharmacologically and economically important, with differences in efficacy and price. Accordingly, an authentication system is needed to combat economically motivated adulteration of Panax products. We conducted comparative analysis of the chloroplast genome sequences of these five species, identifying 34-124 InDels and 141-560 SNPs. Fourteen InDel markers were developed to authenticate the Panax species. Among these, eight were species-unique markers that successfully differentiated one species from the others. We generated at least one species-unique marker for each of the five species, and any of the species can be authenticated by selection among these markers. The markers are reliable, easily detectable, and valuable for applications in the ginseng industry as well as in related research.

Evolution of the Araliaceae family inferred from complete chloroplast genomes and 45S nrDNAs of 10 Panax-related species

We produced complete sequences and conducted comparative analysis of the maternally inherited chloroplast (cp) genomes and bi-parentally inherited 45S nuclear ribosomal RNA genes (nrDNA) from ten Araliaceae species to elucidate the genetic diversity and evolution in that family. The cp genomes ranged from 155,993 bp to 156,730 bp with 97.1-99.6% similarity. Complete 45S nrDNA units were about 11 kb including a 5.8-kb 45S cistron. Among 79 cp protein-coding genes, 74 showed nucleotide variations among ten species, of which infA, rpl22, rps19 and ndhE genes showed the highest Ks values and atpF, atpE, ycf2 and rps15 genes showed the highest Ka/Ks values. Four genes, petN, psaJ, psbF, and psbN, related to photosynthesis and one gene, rpl23, related to the ribosomal large subunit remain conserved in all 10 Araliaceae species. Phylogenetic analysis revealed that the ten species could be resolved into two monophyletic lineages, the Panax-Aralia and the Eleutherococcus-Dendropanax groups, which diverged approximately 8.81-10.59 million years ago (MYA). The Panax genus divided into two groups, with diploid species including P. notoginseng, P. vietnamensis, and P. japonicus surviving in Southern Asia and a tetraploid group including P. ginseng and P. quinquefolius Northern Asia and North America 2.89-3.20 MYA.

Phylogenetic and paleobotanical evidence for late Miocene diversification of the Tertiary subtropical lineage of ivies (Hedera L., Araliaceae)

BACKGROUND

Hedera (ivies) is one of the few temperate genera of the primarily tropical Asian Palmate group of the Araliaceae, which extends its range out of Asia to Europe and the Mediterranean basin. Phylogenetic and phylogeographic results suggested Asia as the center of origin and the western Mediterranean region as one of the secondary centers of diversification. The bird-dispersed fleshy fruits of ivies suggest frequent dispersal over long distances (e.g. Macaronesian archipelagos), although reducing the impact of geographic barriers to gene flow in mainland species. Genetic isolation associated with geographic barriers and independent polyploidization events have been postulated as the main driving forces of diversification. In this study we aim to evaluate past and present diversification patterns in Hedera within a geographic and temporal framework to clarify the biogeographic history of the genus.

RESULTS

Phylogenetic (biogeographic, time divergence and diversification) and phylogeographic (coalescence) analyses using four DNA regions (nrITS, trnH-psbA, trnT-trnL, rpl32) revealed a complex spatial pattern of lineage divergence. Scarce geographic limitation to gene flow and limited diversification are observed during the early-mid Miocene, followed by a diversification rate increase related to geographic divergence from the Tortonian/Messinian. Genetic and palaeobotanical evidence points the origin of the Hedera clade in Asia, followed by a gradual E-W Asian extinction and the progressive E-W Mediterranean colonization. The temporal framework for the E Asia - W Mediterranean westward colonization herein reported is congruent with the fossil record. Subsequent range expansion in Europe and back colonization to Asia is also inferred. Uneven diversification among geographic areas occurred from the Tortonian/Messinian onwards with limited diversification in the newly colonized European and Asian regions. Eastern and western Mediterranean regions acted as refugia for Miocene and post-Miocene lineages, with a similar role as consecutive centers of centrifugal dispersal (including islands) and speciation.

CONCLUSIONS

The Miocene Asian extinction and European survival of Hedera question the general pattern of Tertiary regional extinction of temperate angiosperms in Europe while they survived in Asia. The Tortonian/Messinian diversification increase of ivies in the Mediterranean challenges the idea that this aridity period was responsible for the extinction of the Mediterranean subtropical Tertiary flora. Differential responses of Hedera to geographic barriers throughout its evolutionary history, linked to spatial isolation related to historical geologic and climatic constraints may have shaped diversification of ivies in concert with recurrent polyploidy.

A Dynamic Tandem Repeat in Monocotyledons Inferred from a Comparative Analysis of Chloroplast Genomes in Melanthiaceae

Chloroplast genomes (cpDNA) are highly valuable resources for evolutionary studies of angiosperms, since they are highly conserved, are small in size, and play critical roles in plants. Slipped-strand mispairing (SSM) was assumed to be a mechanism for generating repeat units in cpDNA. However, research on the employment of different small repeated sequences through SSM events, which may induce the accumulation of distinct types of repeats within the same region in cpDNA, has not been documented. Here, we sequenced two chloroplast genomes from the endemic species Heloniopsis tubiflora (Korea) and Xerophyllum tenax (USA) to cover the gap between molecular data and explore "hot spots" for genomic events in Melanthiaceae. Comparative analysis of 23 complete cpDNA sequences revealed that there were different stages of deletion in the rps16 region across the Melanthiaceae. Based on the partial or complete loss of rps16 gene in cpDNA, we have firstly reported potential molecular markers for recognizing two sections (Veratrum and Fuscoveratrum) of Veratrum. Melathiaceae exhibits a significant change in the junction between large single copy and inverted repeat regions, ranging from trnH_GUG to a part of rps3. Our results show an accumulation of tandem repeats in the rpl23-ycf2 regions of cpDNAs. Small conserved sequences exist and flank tandem repeats in further observation of this region across most of the examined taxa of Liliales. Therefore, we propose three scenarios in which different small repeated sequences were used during SSM events to generate newly distinct types of repeats. Occasionally, prior to the SSM process, point mutation event and double strand break repair occurred and induced the formation of initial repeat units which are indispensable in the SSM process. SSM may have likely occurred more frequently for short repeats than for long repeat sequences in tribe Parideae (Melanthiaceae, Liliales). Collectively, these findings add new evidence of dynamic results from SSM in chloroplast genomes which can be useful for further evolutionary studies in angiosperms. Additionally, genomics events in cpDNA are potential resources for mining molecular markers in Liliales.

Comparative Genomics and Phylogenomics of East Asian Tulips (Amana, Liliaceae)

The genus Amana Honda (Liliaceae), when it is treated as separate from Tulipa, comprises six perennial herbaceous species that are restricted to China, Japan and the Korean Peninsula. Although all six Amana species have important medicinal and horticultural uses, studies focused on species identification and molecular phylogenetics are few. Here we report the nucleotide sequences of six complete Amana chloroplast (cp) genomes. The cp genomes of Amana range from 150,613 bp to 151,136 bp in length, all including a pair of inverted repeats (25,629-25,859 bp) separated by the large single-copy (81,482-82,218 bp) and small single-copy (17,366-17,465 bp) regions. Each cp genome equivalently contains 112 unique genes consisting of 30 transfer RNA genes, four ribosomal RNA genes, and 78 protein coding genes. Gene content, gene order, AT content, and IR/SC boundary structure are nearly identical among all Amana cp genomes. However, the relative contraction and expansion of the IR/SC borders among the six Amana cp genomes results in length variation among them. Simple sequence repeat (SSR) analyses of these Amana cp genomes indicate that the richest SSRs are A/T mononucleotides. The number of repeats among the six Amana species varies from 54 (A. anhuiensis) to 69 (Amana kuocangshanica) with palindromic (28-35) and forward repeats (23-30) as the most common types. Phylogenomic analyses based on these complete cp genomes and 74 common protein-coding genes strongly support the monophyly of the genus, and a sister relationship between Amana and Erythronium, rather than a shared common ancestor with Tulipa. Nine DNA markers (rps15-ycf1, accD-psaI, petA-psbJ, rpl32-trnL, atpH-atpI, petD-rpoA, trnS-trnG, psbM-trnD, and ycf4-cemA) with number of variable sites greater than 0.9% were identified, and these may be useful for future population genetic and phylogeographic studies of Amana species.

Plastid phylogenomics and adaptive evolution of Gaultheria series Trichophyllae (Ericaceae), a clade from sky islands of the Himalaya-Hengduan Mountains

Gaultheria series Trichophyllae Airy Shaw is an angiosperm clade of high-alpine shrublets endemic to the Himalaya-Hengduan Mountains and characterized by recent species divergence and convergent character evolution that has until recently caused much confusion in species circumscription. Although multiple DNA sequence regions have been employed previously, phylogenetic relationships among species in the group have remained largely unresolved. Here we examined the effectiveness of the plastid genome for improving phylogenetic resolution within the G. series Trichophyllae clade. Plastid genomes of 31 samples representing all 19 recognized species of the series and three outgroup species were sequenced with Illumina Sequencing technology. Maximum likelihood (ML), maximum parsimony (MP) and Bayesian inference (BI) phylogenetic analyses were performed with various datasets, i.e., that from the whole plastid genome, coding regions, noncoding regions, large single-copy region (LSC) and inverted-repeat region a (IRa). The partitioned whole plastid genome with inverted-repeat region b (IRb) excluded was also analyzed with ML and BI. Tree topologies based on the whole plastid genome, noncoding regions, and LSC region datasets across all analyses, and that based on the partitioned dataset with ML and BI analyses, are identical and generally strongly supported. Gaultheria series Trichophyllae form a clade with three species and one variety that is sister to a clade of the remaining 16 species; the latter comprises seven main subclades. Interspecific relationships within the series are strongly supported except for those based on the coding-region and IRa-region datasets. Eight divergence hotspot regions, each possessing >5% percent variable sites, were screened across the whole plastid genome of the 28 individuals sampled in the series. Results of morphological character evolution reconstruction diagnose several clades, and a hypothesis of adaptive evolution for plant habit is postulated.

The complete chloroplast genome sequence of Hydrocotyle sibthorpioides (Apiales: araliaceae)

Hydrocotyle sibthorpioides Lam. is a widespread and important Chinese medicinal plant. In this study, we have sequenced the complete chloroplast genome of H. sibthorpioides, which is 152,880 bp in length with large (LSC 84,064 bp) and small (SSC 18,690 bp) single-copy regions separated by a pair of inverted repeats (IRs 25,063 bp) and contains 113 unique genes with 17 genes duplicated in the IR making a total of 130 genes. The phylogenetic analysis strongly supports the position of H. sibthorpioides and indicates it belongs to the Araliaceae family, potential to facilitate a better understanding of population and phylogenetic studies in Apiales.

The Chloroplast Genome of Utricularia reniformis Sheds Light on the Evolution of the ndh Gene Complex of Terrestrial Carnivorous Plants from the Lentibulariaceae Family

Lentibulariaceae is the richest family of carnivorous plants spanning three genera including Pinguicula, Genlisea, and Utricularia. Utricularia is globally distributed, and, unlike Pinguicula and Genlisea, has both aquatic and terrestrial forms. In this study we present the analysis of the chloroplast (cp) genome of the terrestrial Utricularia reniformis. U. reniformis has a standard cp genome of 139,725bp, encoding a gene repertoire similar to essentially all photosynthetic organisms. However, an exclusive combination of losses and pseudogenization of the plastid NAD(P)H-dehydrogenase (ndh) gene complex were observed. Comparisons among aquatic and terrestrial forms of Pinguicula, Genlisea, and Utricularia indicate that, whereas the aquatic forms retained functional copies of the eleven ndh genes, these have been lost or truncated in terrestrial forms, suggesting that the ndh function may be dispensable in terrestrial Lentibulariaceae. Phylogenetic scenarios of the ndh gene loss and recovery among Pinguicula, Genlisea, and Utricularia to the ancestral Lentibulariaceae cladeare proposed. Interestingly, RNAseq analysis evidenced that U. reniformis cp genes are transcribed, including the truncated ndh genes, suggesting that these are not completely inactivated. In addition, potential novel RNA-editing sites were identified in at least six U. reniformis cp genes, while none were identified in the truncated ndh genes. Moreover, phylogenomic analyses support that Lentibulariaceae is monophyletic, belonging to the higher core Lamiales clade, corroborating the hypothesis that the first Utricularia lineage emerged in terrestrial habitats and then evolved to epiphytic and aquatic forms. Furthermore, several truncated cp genes were found interspersed with U. reniformis mitochondrial and nuclear genome scaffolds, indicating that as observed in other smaller plant genomes, such as Arabidopsis thaliana, and the related and carnivorous Genlisea nigrocaulis and G. hispidula, the endosymbiotic gene transfer may also shape the U. reniformis genome in a similar fashion. Overall the comparative analysis of the U. reniformis cp genome provides new insight into the ndh genes and cp genome evolution of carnivorous plants from Lentibulariaceae family.

The Complete Chloroplast Genome of the Hare's Ear Root, Bupleurum falcatum: Its Molecular Features

Bupleurum falcatum, which belongs to the family Apiaceae, has long been applied for curative treatments, especially as a liver tonic, in herbal medicine. The chloroplast (cp) genome has been an ideal model to perform the evolutionary and comparative studies because of its highly conserved features and simple structure. The Apiaceae family is taxonomically close to the Araliaceae family and there have been numerous complete chloroplast genome sequences reported in the Araliaceae family, while little is known about the Apiaceae family. In this study, the complete sequence of the B. falcatum chloroplast genome was obtained. The full-length of the cp genome is 155,989 nucleotides with a 37.66% overall guanine-cytosine (GC) content and shows a quadripartite structure composed of three nomenclatural regions: a large single-copy (LSC) region, a small single-copy (SSC) region, and a pair of inverted repeat (IR) regions. The genome occupancy is 85,912-bp, 17,517-bp, and 26,280-bp for LSC, SSC, and IR, respectively. B. falcatum was shown to contain 111 unique genes (78 for protein-coding, 29 for tRNAs, and four for rRNAs, respectively) on its chloroplast genome. Genic comparison found that B. falcatum has no pseudogenes and has two gene losses, accD in the LSC and ycf15 in the IRs. A total of 55 unique tandem repeat sequences were detected in the B. falcatum cp genome. This report is the first to describe the complete chloroplast genome sequence in B. falcatum and will open up further avenues of research to understand the evolutionary panorama and the chloroplast genome conformation in related plant species.

Insight into infrageneric circumscription through complete chloroplast genome sequences of two Trillium species

Genomic events including gene loss, duplication, pseudogenization and rearrangement in plant genomes are valuable sources for exploring and understanding the process of evolution in angiosperms. The family Melanthiaceae is distributed in temperate regions of the Northern Hemisphere and divided into five tribes (Heloniadeae, Chionographideae, Xerophylleae, Melanthieae and Parideae) based on the molecular phylogenetic analyses. At present, complete chloroplast genomes of the Melanthiaceae have been reported from three species. In the previous genomic study of Liliales, atrnI-CAU gene duplication event was reported fromParis verticillata, a member of Parideae. To clarify the significant genomic events of the tribe Parideae, we analysed the complete chloroplast genome sequences of twoTrilliumspecies representing two subgenera:TrilliumandPhyllantherum InTrillium tschonoskii(subgenusTrillium), the circular double-stranded cpDNA sequence of 156 852 bp consists of two inverted repeat (IR) regions of 26 501 bp each, a large single-copy (LSC) region of 83 981 bp and a small single-copy (SSC) region of 19 869 bp. The chloroplast genome sequence ofT. maculatum(subgenusPhyllantherum) is 157 359 bp in length, consisting of two IRs (25 535 bp), one SSC (19 949 bp) and one LSC (86 340 bp), and is longer than that ofT. tschonoskii The results showed that the cpDNAs of Parideae are highly conserved across genome structure, gene order and contents. However, the chloroplast genome ofT. maculatumcontained a 3.4-kb inverted sequence betweenndhCandrbcLin the LSC region, and it was a unique feature for subgeneraPhyllantherum In addition, we found three different types of gene duplication in the intergenic spacer betweenrpl23andycf2containingtrnI-CAU, which were in agreement with the circumscription of subgenera and sections in Parideae excludingT. govanianum These genomic features provide informative molecular markers for identifying the infrageneric taxa ofTrilliumand improve our understanding of the evolution patterns of Parideae in Melanthiaceae.

High-throughput multiplex cpDNA resequencing clarifies the genetic diversity and genetic relationships among Brassica napus, Brassica rapa and Brassica oleracea

Brassica napus (rapeseed) is a recent allotetraploid plant and the second most important oilseed crop worldwide. The origin of B. napus and the genetic relationships with its diploid ancestor species remain largely unresolved. Here, chloroplast DNA (cpDNA) from 488 B. napus accessions of global origin, 139 B. rapa accessions and 49 B. oleracea accessions were populationally resequenced using Illumina Solexa sequencing technologies. The intraspecific cpDNA variants and their allelic frequencies were called genomewide and further validated via EcoTILLING analyses of the rpo region. The cpDNA of the current global B. napus population comprises more than 400 variants (SNPs and short InDels) and maintains one predominant haplotype (Bncp1). Whole-genome resequencing of the cpDNA of Bncp1 haplotype eliminated its direct inheritance from any accession of the B. rapa or B. oleracea species. The distribution of the polymorphism information content (PIC) values for each variant demonstrated that B. napus has much lower cpDNA diversity than B. rapa; however, a vast majority of the wild and cultivated B. oleracea specimens appeared to share one same distinct cpDNA haplotype, in contrast to its wild C-genome relatives. This finding suggests that the cpDNA of the three Brassica species is well differentiated. The predominant B. napus cpDNA haplotype may have originated from uninvestigated relatives or from interactions between cpDNA mutations and natural/artificial selection during speciation and evolution. These exhaustive data on variation in cpDNA would provide fundamental data for research on cpDNA and chloroplasts.

The complete plastid genome sequence of Panax notoginseng, a famous traditional Chinese medicinal plant of the family Araliaceae

We report complete nucleotide sequence of the Panax notoginseng chloroplast genome using next-generation sequencing technology. The genome consists of 156 324 bp containing a pair of inverted repeats (IRs) of 26 105 bp, which was separated by a large single-copy region and a small single-copy region of 86 082 bp and 18 032 bp, respectively. The P. notoginseng cp genome encodes 114 unigenes (80 protein-coding genes, 4 rRNA genes, and 30 tRNA genes), in which 18 are duplicated in the IR regions. The genic regions account for 51.1% of whole cp genome, and the GC content of the plastome was 38.1%. A phylogenomic analysis of the 10 complete chloroplast genomes from Araliaceae using Hydrocotyle verticillata outgroup showed that P. notoginseng is closely related to P. ginseng that belongs to the genus Panax.

The complete chloroplast genome of Eleutherococcus gracilistylus (W.W.Sm.) S.Y.Hu (Araliaceae)

Eleutherococcus gracilistylus is a plant species that is close to E. senticosus, a famous medicinal plant called Siberian ginseng. The complete chloroplast genome sequence of the E. gracilistylus was determined by de novo assembly using whole genome next generation sequences. The chloroplast genome of E. gracilistylus was 156 770 bp long and showed distinct four partite structures such as a large single copy region of 86 729 bp, a small single copy region of 18 175 bp, and a pair of inverted repeat regions of 25 933 bp. The overall GC contents of the genome sequence were 36.8%. The chloroplast genome of E. gracilistylus contains 79 protein-coding sequences, 30 tRNA genes, and four rRNA genes. The phylogenetic analysis with the reported chloroplast genomes confirmed close taxonomical relationship of E. gracilistylus with E. senticosus.

Development of chloroplast genomic resources for Cynara

In this study, new chloroplast (cp) resources were developed for the genus Cynara, using whole cp genomes from 20 genotypes, by means of high-throughput sequencing technologies. Our target species included seven globe artichokes, two cultivated cardoons, eight wild artichokes, and three other wild Cynara species (C. baetica, C. cornigera and C. syriaca). One complete cp genome was isolated using short reads from a whole-genome sequencing project, while the others were obtained by means of long-range PCR, for which primer pairs are provided here. A de novo assembly strategy combined with a reference-based assembly allowed us to reconstruct each cp genome. Comparative analyses among the newly sequenced genotypes and two additional Cynara cp genomes ('Brindisino' artichoke and C. humilis) retrieved from public databases revealed 126 parsimony informative characters and 258 singletons in Cynara, for a total of 384 variable characters. Thirty-nine SSR loci and 34 other INDEL events were detected. After data analysis, 37 primer pairs for SSR amplification were designed, and these molecular markers were subsequently validated in our Cynara genotypes. Phylogenetic analysis based on all cp variable characters provided the best resolution when compared to what was observed using only parsimony informative characters, or only short 'variable' cp regions. The evaluation of the molecular resources obtained from this study led us to support the 'super-barcode' theory and consider the total cp sequence of Cynara as a reliable and valuable molecular marker for exploring species diversity and examining variation below the species level.