Background: The Peucedanum genus is the backbone member of Apiaceae, with many economically and medically important plants. Although the previous studies on Peucedanum provide us with a good research basis, there are still unclear phylogenetic relationships and many taxonomic problems in Peucedanum, and a robust phylogenetic framework of this genus still has not been obtained, which severely hampers the improvement and revision of taxonomic system for this genus. The plastid genomes possessing more variable characters have potential for reconstructing a robust phylogeny in plants. Results: In the current study, we newly sequenced and assembled seven Peucedanum plastid genomes. Together with five previously published plastid genomes of Peucedanum, we performed a comprehensively comparative analyses for this genus. Twelve Peucedanum plastomes were similar in terms of genome structure, codon bias, RNA editing sites, and SSRs, but varied in genome size, gene content and arrangement, and border of SC/IR. Fifteen mutation hotspot regions were identified among plastid genomes that can serve as candidate DNA barcodes for species identification in Peucedanum. Our phylogenetic analyses based on plastid genomes generated a phylogeny with high supports and resolutions for Peucedanum that robustly supported the non-monophyly of genus Peucedanum. Conclusion: The plastid genomes of Peucedanum showed both conservation and diversity. The plastid genome data were efficient and powerful for improving the supports and resolutions of phylogeny for the complex Peucedanum genus. In summary, our study provides new sights into the plastid genome evolution, taxonomy, and phylogeny for Peucedanum species.

Lilium henrici Franchet, which belongs to the family Liliaceae, is an endangered plant native to China. The wild populations of L. henrici have been largely reduced by habitat degradation or loss. In our study, we determined the whole chloroplast genome sequence for L. henrici and compared its structure with other Lilium (including Nomocharis) species. The chloroplast genome of L. henrici is a circular structure and 152,784 bp in length. The large single copy and small single copy is 82,429 bp and 17,533 bp in size, respectively, and the inverted repeats are 26,411 bp in size. The L. henrici chloroplast genome contains 116 different genes, including 78 protein coding genes, 30 tRNA genes, 4 rRNA genes, and 4 pseudogenes. There were 51 SSRs detected in the L. henrici chloroplast genome sequence. Genic comparison among L. henrici with other Lilium (including Nomocharis) chloroplast genomes shows that the sequence lengths and gene contents show little variation, the only differences being in three pseudogenes. Phylogenetic analysis revealed that N. pardanthina was a sister species to L. henrici. Overall, this study, providing L. henrici genomic resources and the comparative analysis of Lilium chloroplast genomes, will be beneficial for the evolutionary study and phylogenetic reconstruction of the genus Lilium, molecular barcoding in population genetics.

Exploring the effects of orographic events and climatic shifts on the geographic distribution of organisms in the Himalayas-Hengduan Mountains (HHM) region and Qinghai-Tibetan Plateau (QTP) is crucial to understand the impact of environmental changes on organism evolution. To gain further insight into these processes, we reconstructed the evolutionary history of nine Chamaesium species distributed across the HHM and QTP regions. In total, 525 individuals from 56 populations of the nine species were analyzed based on three maternally inherited chloroplast fragments (rpl16, trnT-trnL, and trnQ-rps16) and one nuclear DNA region (internal transcribed spacer, ITS). Fifty-two chloroplast DNA (cpDNA) and 47 ITS haplotypes were identified in nine species. All of the cpDNA and ITS haplotypes were species-specific. Phylogenetic analysis suggested that all nine species form a monophyletic clade with high support. Dating analysis and ancestral area reconstruction revealed that the ancestral group of Chamaesium originated in the southern Himalayan region at the beginning of the Paleogene (60.85 Ma). The nine species of Chamaesium then separated well during the last 25 million years started in Miocene. Our maxent modeling indicated the broad-scale distributions of all nine species remained fairly stable from LIG to the present and predicted that it will remain stable into the future. The initial split of Chamaesium was triggered by climate changes following the collision of the Indian plate with the Eurasia plate during the Eocene. Subsequently, divergences within Chamaesium may have been induced by the intense uplift of the QTP, the onset of the monsoon system, and Central Asian aridification. Long evolutionary history, sexual reproduction, and habitat fragmentation could contribute to the high level of genetic diversity of Chamaesium. The higher genetic differentiation among Chamaesium populations may be related to the drastic changes of the external environment in this region and limited seed/pollen dispersal capacity.

Fritillaria is a genus that has important medicinal and horticultural values. The study involved the most comprehensive chloroplast genome samples referring to Old and New World clades of Fritillaria for marker selection and phylogenetic studies. We reported and compared eleven newly sequenced whole-plastome sequences of Fritillaria which proved highly similar in overall size (151,652–152,434 bp), genome structure, gene content, and order. Comparing them with other species of Liliales (6 out of 10 families) indicated the same similarity but showed some structural variations due to the contraction or expansion of the inverted repeat (IR) regions. A/T mononucleotides, palindromic, and forward repeats were the most common types. Six hypervariable regions (rps16-trnQ, rbcL-accD, accD-psaI, psaJ-rpl33, petD-rpoA, and rpl32-trnL) were discovered based on 26 Fritillaria whole-plastomes to be potential molecular markers. Based on the plastome data that were collected from 26 Fritillaria and 21 Lilium species, a phylogenomic study was carried out with three Cardiocrinum species as outgroups. Fritillaria was sister to Lilium with a high support value, and the interspecies relationships within subgenus Fritillaria were resolved very well. The six hypervariable regions can be used as candidate DNA barcodes of Fritillaria and the phylogenomic framework can guide extensive genomic sampling for further phylogenetic analyses.

Pleurospermum is a taxonomically challenging taxon of Apiaceae, as its circumscription and composition remain controversial for morphological similarities with several related genera, leading to a dispute between Pleurospermum in the broad sense and strict sense. While evidence from previous molecular studies recognized plural branching lineages within the Pleurospermum s. l., it did not support the latest delimitation of Pleurospermum s. str. by only two closely related northern species. So far, no proper delimitation for Pleurospermum has come up, and many of the plural taxa in Pleurospermum s. l. remain unresolved, which may be due to poor phylogenetic resolution yielded barely from ITS sequences. Herein, we newly assembled 40 complete plastomes from 36 species of Pleurospermum s. l. and related genera, 34 of which were first reported and generated a well-resolved backbone phylogeny in a framework of the subfamily Apioideae. From the phylogeny with greatly improved resolution, a total of six well-supported monophyletic lineages within Pleurospermum s. l. were recognized falling in different major clades of Apioideae. Combining morphological characteristics with phylogenetic inference, we suggested to re-delimit the Pleurospermum s. str. by introducing nine species mainly from the Himalayan regions and proposed its boundary features; the remaining species were suggested to be excluded from Pleurospermum to incorporate into their more related taxa being revealed. On this basis, the plastome comparison revealed not only the high conservatism but also the mild differences among lineages in plastome structure and gene evolution. Overall, our study provided a backbone phylogeny essential for further studies of the taxonomically difficult taxa within Pleurospermum s. l.

Hansenia Turcz., Haplosphaera Hand.-Mazz. and Sinodielsia H.Wolff are three Apiaceae genera endemic to the Hengduan Mountains and the Himalayas, which usually inhabit elevations greater than 2000 m. The phylogenetic relationships between and within the genera were uncertain, especially the placement of Hap. himalayensis and S. microloba. Therefore, we aimed to conduct comparative (simple sequence repeat (SSR) structure, codon usage bias, nucleotide diversity (Pi) and inverted repeat (IR) boundaries) and phylogenetic analyses of Hansenia, Haplosphaera and Sinodielsia (also compared with Chamaesium and Bupleurum) to reduce uncertainties in intergeneric and interspecific relationships. We newly assembled eight plastid genomes from Hansenia, Haplosphaera and Sinodielsia species, and analyzed them with two plastid genomes from GenBank of Hap. phaea,S. yunnanensis. Phylogenetic analyses used these ten genomes and another 22 plastid genome sequences of Apiaceae. We found that the newly assembled eight genomes ranged from 155,435 bp to 157,797 bp in length and all had a typical quadripartite structure. Fifty-five to 75 SSRs were found in Hansenia, Haplosphaera and Sinodielsia species, and the most abundant SSR was mononucleotide, which accounted for 58.47% of Hansenia, 60.21% of Haplosphaera and 48.01% of Sinodielsia. There was no evident divergence of codon usage frequency between the three genera, where codons ranged from 21,134 to 21,254. The Pi analysis showed that trnE(UUC)-trnT(GGU), trnH(GUG)-psbA and trnE(UUC)-trnT(GGU) spacer regions had the highest Pi values in the plastid genomes of Hansenia (0.01889), Haplosphaera (0.04333) and Sinodielsia (0.01222), respectively. The ndhG-ndhI spacer regions were found in all three genera to have higher diversity values (Pi values: 0.01028–0.2), and thus may provide potential DNA barcodes in phylogenetic analysis. IR boundary analysis showed that the length of rps19 and ycf1 genes entering IRs were usually stable in the same genus. Our phylogenetic tree demonstrated that Hap. himalayensis is sister to Han. weberbaueriana; meanwhile, Haplosphaera and Hansenia are nested together in the East Asia clade, and S. microloba is nested within individuals of S. yunnanensis in the Acronema clade. This study will enrich the complete plastid genome dataset of the Apiaceae genera and has provided a new insight into phylogeny reconstruction using complete plastid genomes of Hansenia, Haplosphaera and Sinodielsia.