The complete plastomes of thirteen Libanotis (Apiaceae, Apioideae) plants: comparative and phylogenetic analyses provide insights into the plastome evolution and taxonomy of Libanotis.

BACKGROUND: The genus Libanotis Haller ex Zinn, nom. cons., a contentious member of Apiaceae, encompasses numerous economically and medicinally significant plants, comprising approximately 30 species distributed across Eurasia. Despite many previous taxonomic insights into it, phylogenetic studies of the genus are still lacking. And the establishment of a robust phylogenetic framework remains elusive, impeding advancements and revisions in the taxonomic system for this genus. Plastomes with greater variability in their genetic characteristics hold promise for building a more robust Libanotis phylogeny. RESULTS: During our research, we sequenced, assembled, and annotated complete plastomes for twelve Libanotis species belong to three sections and two closely related taxa. We conducted a comprehensive comparative analysis through totally thirteen Libanotis plastomes for the genus, including an additional plastome that had been published. Our results suggested that Libanotis plastome was highly conserved between different subclades, while the coding regions were more conserved than the non-coding regions, and the IR regions were more conserved than the single copy regions. Nevertheless, eight mutation hotspot regions were identified among plastomes, which can be considered as candidate DNA barcodes for accurate species identification in Libanotis. The phylogenetic analyses generated a robustly framework for Libanotis and revealed that Libanotis was not a monophyletic group and their all three sections were polygenetic. Libanotis schrenkiana was sister to L. sibirica, type species of this genus, but the remainders scattered within Selineae. CONCLUSION: The plastomes of Libanotis exhibited a high degree of conservation and was effective in enhancing the support and resolution of phylogenetic analyses within this genus. Based on evidence from both phylogeny and morphology, we propose the recognition of "Libanotis sensu stricto" and provide taxonomic recommendations for other taxa that previously belonged to Libanotis. In conclusion, our study not only revealed the phylogenetic position and plastid evolution of Libanotis, but also provided new insights into the phylogeny of the family Apiaceae and phylogenetic relationships within the tribe Selineae.

Phylogeny and diversification of genus Sanicula L. (Apiaceae): novel insights from plastid phylogenomic analyses.

BACKGROUND: The genus Sanicula L. is a unique perennial herb that holds important medicinal values. Although the previous studies on Sanicula provided us with a good research basis, its taxonomic system and interspecific relationships have not been satisfactorily resolved, especially for those endemic to China. Moreover, the evolutionary history of this genus also remains inadequately understood. The plastid genomes possessing highly conserved structure and limited evolutionary rate have proved to be an effective tool for studying plant phylogeny and evolution. RESULTS: In the current study, we newly sequenced and assembled fifteen Sanicula complete plastomes. Combined with two previously reported plastomes, we performed comprehensively plastid phylogenomics analyses to gain novel insights into the evolutionary history of this genus. The comparative results indicated that the seventeen plastomes exhibited a high degree of conservation and similarity in terms of their structure, size, GC content, gene order, IR borders, codon bias patterns and SSRs profiles. Such as all of them displayed a typical quadripartite structure, including a large single copy region (LSC: 85,074-86,197 bp), a small single copy region (SSC: 17,047-17,132 bp) separated by a pair of inverted repeat regions (IRs: 26,176-26,334 bp). And the seventeen plastomes had similar IR boundaries and the adjacent genes were identical. The rps19 gene was located at the junction of the LSC/IRa, the IRa/SSC junction region was located between the trnN gene and ndhF gene, the ycf1 gene appeared in the SSC/IRb junction and the IRb/LSC boundary was located between rpl12 gene and trnH gene. Twelve specific mutation hotspots (atpF, cemA, accD, rpl22, rbcL, matK, ycf1, trnH-psbA, ycf4-cemA, rbcL-accD, trnE-trnT and trnG-trnR) were identified that can serve as potential DNA barcodes for species identification within the genus Sanicula. Furthermore, the plastomes data and Internal Transcribed Spacer (ITS) sequences were performed to reconstruct the phylogeny of Sanicula. Although the tree topologies of them were incongruent, both provided strong evidence supporting the monophyly of Saniculoideae and Apioideae. In addition, the sister groups between Saniculoideae and Apioideae were strongly suggested. The Sanicula species involved in this study were clustered into a clade, and the Eryngium species were also clustered together. However, it was clearly observed that the sections of Sanicula involved in the current study were not respectively recovered as monophyletic group. Molecular dating analysis explored that the origin of this genus was occurred during the late Eocene period, approximately 37.84 Ma (95% HPD: 20.33-52.21 Ma) years ago and the diversification of the genus was occurred in early Miocene 18.38 Ma (95% HPD: 10.68-25.28 Ma). CONCLUSION: The plastome-based tree and ITS-based tree generated incongruences, which may be attributed to the event of hybridization/introgression, incomplete lineage sorting (ILS) and chloroplast capture. Our study highlighted the power of plastome data to significantly improve the phylogenetic supports and resolutions, and to efficiently explore the evolutionary history of this genus. Molecular dating analysis explored that the diversification of the genus occurred in the early Miocene, which was largely influenced by the prevalence of the East Asian monsoon and the uplift of the Hengduan Mountains (HDM). In summary, our study provides novel insights into the plastome evolution, phylogenetic relationships, taxonomic framework and evolution of genus Sanicula.

Phylogeny and adaptive evolution of subgenus Rhizirideum (Amaryllidaceae, Allium) based on plastid genomes.

The subgenus Rhizirideum in the genus Allium consists of 38 species worldwide and forms five sections (A. sect. Rhizomatosa, A. sect. Tenuissima, A. sect. Rhizirideum, A. sect. Eduardia, and A. sect. Caespitosoprason), A. sect. Caespitosoprason being merged into A. sect. Rhizomatosa recently. Previous studies on this subgenus mainly focused on separate sections. To investigate the inter-section and inter-subgenera phylogenetic relationships and adaptive evolution of A. subg. Rhizirideum, we selected thirteen representative species, which cover five sections of this subgenus and can represent four typical phenotypes of it. We conducted the comparative plastome analysis with our thirteen plastomes. And phylogenetic inferences with CDSs and complete sequences of plastomes of our thirteen species and another fifty-four related species were also performed. As a result, the A. subg. Rhizirideum plastomes were relatively conservative in structure, IR/SC borders, codon usage, and repeat sequence. In phylogenetic results, the inter-subgenera relationships among A. subg. Rhizirideum and other genus Allium subgenera were generally similar to the previous reports. In contrast, the inter-section relationships within our subgenus A. subg. Rhizirideum were newly resolved in this study. A. sect. Rhizomatosa and A. sect. Tenuissima were sister branches, which were then clustered with A. sect. Rhizirideum and A. sect. Eduardia successively. However, Allium Polyrhizum Turcz. ex Regel, type species of A. sect. Caespitosoprason, was resolved as the basal taxon of A. subg. Rhizirideum. Allium siphonanthum J. M. Xu was also found in clade A. subg. Cyathophora instead of clade A. subg. Rhizirideum. The selective pressure analysis was also conducted, and most protein-coding genes were under purifying selection. At the same time, just one gene, ycf2, was found under positive selection, and another three genes (rbcL, ycf1a, ycf1b) presented relaxed selection, which were all involved in the photosynthesis. The low temperature, dry climate, and high altitude of the extreme habitats where A. subg. Rhizirideum species grow might impose intense natural selection forces on their plastome genes for photosynthesis. In summary, our research provides new insights into the phylogeny and adaptive evolution of A. subg. Rhizirideum. Moreover, we suggest that the positions of the A. subg. Rhizirideum species A. polyrhizum and A. siphonanthum should be reconsidered.

The plastid genome of twenty-two species from Ferula, Talassia, and Soranthus: comparative analysis, phylogenetic implications, and adaptive evolution.

BACKGROUND: The Ferula genus encompasses 180-185 species and is one of the largest genera in Apiaceae, with many of Ferula species possessing important medical value. The previous studies provided more information for Ferula, but its infrageneric relationships are still confusing. In addition, its genetic basis of its adaptive evolution remains poorly understood. Plastid genomes with more variable sites have the potential to reconstruct robust phylogeny in plants and investigate the adaptive evolution of plants. Although chloroplast genomes have been reported within the Ferula genus, few studies have been conducted using chloroplast genomes, especially for endemic species in China. RESULTS: Comprehensively comparative analyses of 22 newly sequenced and assembled plastomes indicated that these plastomes had highly conserved genome structure, gene number, codon usage, and repeats type and distribution, but varied in plastomes size, GC content, and the SC/IR boundaries. Thirteen mutation hotspot regions were detected and they would serve as the promising DNA barcodes candidates for species identification in Ferula and related genera. Phylogenomic analyses with high supports and resolutions showed that Talassia transiliensis and Soranthus meyeri were nested in the Ferula genus, and thus they should be transferred into the Ferula genus. Our phylogenies also indicated the monophyly of subgenera Sinoferula and subgenera Narthex in Ferula genus. Twelve genes with significant posterior probabilities for codon sites were identified in the positively selective analysis, and their function may relate to the photosystem II, ATP subunit, and NADH dehydrogenase. Most of them might play an important role to help Ferula species adapt to high-temperatures, strong-light, and drought habitats. CONCLUSION: Plastome data is powerful and efficient to improve the support and resolution of the complicated Ferula phylogeny. Twelve genes with significant posterior probabilities for codon sites were helpful for Ferula to adapt to the harsh environment. Overall, our study supplies a new perspective for comprehending the phylogeny and evolution of Ferula.

The complete plastomes of seven Peucedanum plants: comparative and phylogenetic analyses for the Peucedanum genus.

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.

The phylogeny of Seseli (Apiaceae, Apioideae): insights from molecular and morphological data.

BACKGROUND: The genus Seseli L., which consists of 125-140 species distributed in the Old World from western Europe and northwestern Africa to China and Japan, is one of the largest and most taxonomically difficult genera of Apiaceae Lindl. Although several previous studies have been conducted on Seseli based on limited morphological characteristics and molecular fragments, a robust and comprehensive phylogeny of Seseli remains elusive. Plastomes provide abundant genetic information and have been widely used in studying plant phylogeny and evolution. Consequently, we newly generated the complete plastomes of eleven Seseli taxa. We combined plastome data and morphological characteristics to investigate the phylogeny of Seseli. RESULTS: In our study, we observed that the genome length, gene numbers, IR/SC borders, and repeat composition of the eleven Seseli plastomes were variable. Several appropriate mutation hotspot regions may be developed as candidate DNA barcodes for evolution, phylogeny, and species identification of Seseli. The phylogenetic results identified that Seseli was not a monophyletic group. Moreover, the eleven newly sequenced Seseli taxa did not cluster with S. tortuosum (the type species of Seseli, belonging to the tribe Selineae), where S. delavayi clustered with Eriocycla belonging to the tribe Echinophoreae and the other ten belonged to Selineae. The comparative plastome and morphological characteristics analyses confirmed the reliability of the phylogenetic analyses and implied the complex evolution of Seseli. CONCLUSION: Combining molecular and morphological data is efficient and useful for studying the phylogeny of Seseli. We suggest that "a narrow sense" of Seseli will be meaningful for further study and the current taxonomic system of Seseli needs to be revised. In summary, our study can provide new insights into the phylogenetic relationships and taxonomic framework of Seseli.

Plastid phylogenomic insights into relationships, divergence, and evolution of Apiales.

MAIN CONCLUSION: Members of Apiales are monophyletic and radiated in the Late Cretaceous. Fruit morphologies are critical for Apiales evolution and negative selection and mutation pressure play important roles in environmental adaptation. Apiales include many foods, spices, medicinal, and ornamental plants, but the phylogenetic relationships, origin and divergence, and adaptive evolution remain poorly understood. Here, we reconstructed Apiales phylogeny based on 72 plastid genes from 280 species plastid genomes representing six of seven families of this order. Highly supported phylogenetic relationships were detected, which revealed that each family of Apiales is monophyletic and confirmed that Pennanticeae is a member of Apiales. Genera Centella and Dickinsia are members of Apiaceae, and the genus Hydrocotyle previously classified into Apiaceae is confirmed to belong to Araliaceae. Besides, coalescent phylogenetic analysis and gene trees cluster revealed ten genes that can be used for distinguishing species among families of Apiales. Molecular dating suggested that the Apiales originated during the mid-Cretaceous (109.51 Ma), with the families' radiation occurring in the Late Cretaceous. Apiaceae species exhibit higher differentiation compared to other families. Ancestral trait reconstruction suggested that fruit morphological evolution may be related to shifts in plant types (herbaceous or woody), which in turn is related to the distribution areas and species numbers. Codon bias and positive selection analyses suggest that negative selection and mutation pressure may play important roles in environmental adaptation of Apiales members. Our results improve the phylogenetic framework of Apiales and provide insights into the origin, divergence, and adaptive evolution of this order and its members.

Phylogeny and highland adaptation of Chinese species in Allium section Daghestanica (Amaryllidaceae) revealed by transcriptome sequencing.

Allium L. is one of the largest monocotyledonous genera with extensive distribution in the Northern Hemisphere. The fundamental phylogenies of Allium have been investigated using many morphological and molecular characters. However, the morphological characters may not agree with the molecular results in some Allium groups or sections (such as the Chinese Allium section Daghestanica), which may result in ambiguous species relationships and hinder further evolutionary and adaptive researches. Here, transcriptome sequences of the six Chinese endemics from Allium section Daghestanica were collected, with their single-copy genes (SCGs) were extracted. The interspecies relationships were analyzed using concatenation and coalescent methods. The branch-site model (BSM) was conducted to detect the positively selected genes (PSGs) in five highland species of this section. Based on 1644, 1281 and 1580 SCGs in flowers, leaves, and flowers-leaves combination respectively, a robust consistent and well-resolved phylogeny was generated from the concatenation method. Strong conflicts among individual gene trees were detected in the coalescent method, and morphological characters were incongruent with molecular relationships to some degree. Many PSGs were involved in responses of various stresses and stimuli (e.g. hypoxia, low temperature, aridity), DNA repair, metabolism, nutrient or energy intake, photosynthesis, and signal transduction. Our study revealed a clear interspecies relationship of Chinese endemics in Allium section Daghestanica and suggested that the discordance between morphological characters and molecular relationships might result from that the former are more susceptible to convergence compared with the latter. PSGs detected in our study may provide some insights into highland adaptation in Allium species.

Insights into phylogeny, age and evolution of Allium (Amaryllidaceae) based on the whole plastome sequences.

BACKGROUND AND AIMS: The genus Allium L., one of the largest monocotyledonous genera and one that includes many economically important crops with nutritional and medicinal value, has been the focus of classification or phylogeny studies for centuries. Recent studies suggested that the genus can be divided into 15 subgenera and 72 sections, which were further classified into three evolutionary lineages. However, the phylogenetic relationships reconstructed by one or two loci showed weaker support, especially for the third evolutionary lineage, which might not show the species relationships very clearly and could hinder further adaptive and evolutionary study. METHODS: In this study, a total of 39 complete chloroplast genomes of Allium (covering 12 Allium subgenera) were collected, and combining these with 125 species of plastomes from 19 other families of monocots, we reconstructed the phylogeny of the genus Allium, estimated the origin and divergence time of the three evolutionary lineages and investigated the adaptive evolution in this genus and related families. RESULTS: Our phylogenetic analysis confirmed the monophyly and three evolutionary lineages of Allium, while new species relationships were detected within the third evolutionary lineage. The divergence time of the three evolutionary lineages was estimated to be in the early Eocene to the middle Miocene, and numerous positive selected genes (PSGs) and PSGs with high average Ka/Ks values were found in Allium species. CONCLUSIONS: Our results detected a well-supported phylogenetic relationship of Allium. The PSGs and PSGs with high Ka/Ks values, as well as diversified morphologies, complicated chromosome characteristics and unique reproductive modes may play important roles in the adaptation and evolution of Allium species. This is the first study that conducted phylogenetic and evolutionary analyses on the genus Allium combined with the plastome and morphological and cytological data. We hope that this study can contribute to further analysis of Allium for other researchers.

A transcriptome-based study on the phylogeny and evolution of the taxonomically controversial subfamily Apioideae (Apiaceae).

BACKGROUND AND AIMS: A long-standing controversy in the subfamily Apioideae concerns relationships among the major lineages, which has prevented a comprehensive study of their fruits and evolutionary history. Here we use single copy genes (SCGs) generated from transcriptome datasets to generate a reliable species tree and explore the evolutionary history of Apioideae. METHODS: In total, 3351 SCGs were generated from 27 transcriptome datasets and one genome, and further used for phylogenetic analysis using coalescent-based methods. Fruit morphology and anatomy were studied in combination with the species tree. Eleven SCGs were screened out for dating analysis with two fossils selected for calibration. KEY RESULTS: A well-supported species tree was generated with a topology [Chamaesieae, (Bupleureae, (Pleurospermeae, (Physospermopsis Clade, (Group C, (Group A, Group B)))))] that differed from previous trees. Daucinae and Torilidinae were not in the tribe Scandiceae and existed as sister groups to the Acronema Clade. Five branches (I-V) of the species tree showed low quartet support but strong local posterior probabilities. Dating analysis suggested that Apioideae originated around 56.64 Mya (95 % highest posterior density interval, 45.18-73.53 Mya). CONCLUSIONS: This study resolves a controversial phylogenetic relationship in Apioideae based on 3351 SCGs and coalescent-based species tree estimation methods. Gene trees that contributed to the species tree may undergoing rapid evolutionary divergence and incomplete lineage sorting. Fruits of Apioideae might have evolved in two directions, anemochorous and hydrochorous, with epizoochorous as a derived mode. Molecular and morphological evidence suggests that Daucinae and Torilidinae should be restored to the tribe level. Our results provide new insights into the morphological evolution of this subfamily, which may contribute to a better understanding of species diversification in Apioideae. Molecular dating analysis suggests that uplift of the Qinghai-Tibetan Plateau (QTP) and climate changes probably drove rapid speciation and diversification of Apioideae in the QTP region.

The effect of Hengduan Mountains Region (HMR) uplift to environmental changes in the HMR and its eastern adjacent area: Tracing the evolutionary history of Allium section Sikkimensia (Amaryllidaceae).

Exploring the effects of orographic events and climatic shifts on geographic distribution of organism in the Hengduan Mountains Region (HMR) and its eastern adjacent area is crucial to the understanding of the environmental changes to organismal evolution. To gain further insight into these processes, we reconstruct evolutionary history of ten species in Allium section Sikkimensia, distributed across regions abovementioned. Using chloroplast and nuclear sequence variation of 79 populations of these ten Allium species with known morphological preferences, we elucidate the phylogenetic relationship, divergence time, ancestral area and genetic structures. Climatic variables analysis, Isolation by distance (IBD) and environment (IBE) and Species distribution modeling (SDM) were analyzed along different genetic clades. These analyses indicated that the initial split of Sikkimensia was triggered by climate changes following Qinghai-Tibet Plateau sensu lato (QTPsl) uplift during the late Miocene. Subsequently, divergences within lineage (lineage A)/among lineages (lineage C and D) in Sikkimensia may be induced by the intense uplift of the HMR around 3-4 Ma and abrupt intensifying of the Asian monsoon regimes. Furthermore, Sikkimensia populations exhibited lopsided demographic history in the Last Glacial Maximum (LGM), as was indicated by the expansion of their range in the QDM and contraction in the HMR. Our findings appear to suggest that the HMR uplift could have strengthened the orographic difference between the HMR and its eastern adjacent area and led to a colder climate in the HMR, while geological topography also played an important role for taxa to respond the climate change that had taken place in the HMR and its eastern adjacent area during the Pleistocene.

Molecular phylogenetics and historical biogeography of the tribe Lilieae (Liliaceae): bi-directional dispersal between biodiversity hotspots in Eurasia.

Background and Aims: The role played by the Qinghai-Tibet Plateau (QTP) in the organismal diversification and biogeography of plants in the Northern Hemisphere has attracted much attention from evolutionary biologists. Here we use tribe Lilieae (Liliaceae), including primarily temperate and alpine lineages with disjunct distributions in the North Temperate Zone, as a case study to shed light upon these processes. Methods: Using 191 taxa (five outgroup taxa) comprising more than 60 % of extant Lilieae species across the entire geographical range, we analyse phylogenetic relationships based on three plastid markers (matK, rbcL, rpl16) and nuclear ITS. Divergence time estimation and ancestral range reconstruction were further inferred. Key Results: The results support a monophyletic Lilieae divided into four clades. Lilium is nested within Fritillaria, which is paraphyletic and partitioned into two clades, New World and Old World, in the chloroplast DNA (cpDNA) analysis. Incongruences between the ITS and cpDNA trees may be explained by divergent ITS paralogues and hybridization. Lilieae originated around 40-49 (28-67) Mya and probably diversified in the QTP region with four major clades that were established during the Oligocene and the Early Miocene. Uplift of the QTP and climatic changes probably drove early diversification of Lilieae in the QTP region. A rapid radiation occurred during the Late Miocene and the Pleistocene, coinciding temporally with recent orogenic process in the QTP region and climatic oscillations. Several lineages dispersed out of the QTP. Conclusions: Lineage persistence and explosive radiation were important processes for establishing high species diversity of Lilieae in the QTP region. Both long-distance dispersal and migration across Beringia probably contributed to the modern distribution range of Lilieae. Our study shows that biotic interchanges between the QTP region and Irano-Turanian region and the Mediterranean Basin were bi-directional, suggesting the latter was a secondary centre of diversity.

Comparative Analysis of the Chloroplast Genomes of the Chinese Endemic Genus Urophysa and Their Contribution to Chloroplast Phylogeny and Adaptive Evolution.

Urophysa is a Chinese endemic genus comprising two species, Urophysa rockii and Urophysa henryi. In this study, we sequenced the complete chloroplast (cp) genomes of these two species and of their relative Semiquilegia adoxoides. Illumina sequencing technology was used to compare sequences, elucidate the intra- and interspecies variations, and infer the phylogeny relationship with other Ranunculaceae family species. A typical quadripartite structure was detected, with a genome size from 158,473 to 158,512 bp, consisting of a pair of inverted repeats separated by a small single-copy region and a large single-copy region. We analyzed the nucleotide diversity and repeated sequences components and conducted a positive selection analysis by the codon-based substitution on single-copy coding sequence (CDS). Seven regions were found to possess relatively high nucleotide diversity, and numerous variable repeats and simple sequence repeats (SSR) markers were detected. Six single-copy genes (atpA, rpl20, psaA, atpB, ndhI, and rbcL) resulted to have high posterior probabilities of codon sites in the positive selection analysis, which means that the six genes may be under a great selection pressure. The visualization results of the six genes showed that the amino acid properties across each column of all species are variable in different genera. All these regions with high nucleotide diversity, abundant repeats, and under positive selection will provide potential plastid markers for further taxonomic, phylogenetic, and population genetics studies in Urophysa and its relatives. Phylogenetic analyses based on the 79 single-copy genes, the whole complete genome sequences, and all CDS sequences showed same topologies with high support, and U. rockii was closely clustered with U. henryi within the Urophysa genus, with S. adoxoides as their closest relative. Therefore, the complete cp genomes in Urophysa species provide interesting insights and valuable information that can be used to identify related species and reconstruct their phylogeny.

The Chloroplast Genome of Lilium henrici: Genome Structure and Comparative Analysis.

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.

Molecular phylogeny, biogeography and ecological niche modelling of Cardiocrinum (Liliaceae): insights into the evolutionary history of endemic genera distributed across the Sino-Japanese floristic region.

BACKGROUND AND AIMS: The patterns of evolutionary assembly in the Sino-Japanese floristic region (SJFR) remain largely unknown due to a lack of integrative multidimensional studies throughout the region. To address this issue, we elucidated the evolutionary history of Cardiocrinum (Liliaceae), a genus containing four taxa distributed across the SJFR. METHODS: Fifty-four populations were sampled throughout the geographical range of Cardiocrinum to assess genetic structure, analyse phylogenetic relationships and reconstruct ancestral area based on six chloroplast DNA (cpDNA) fragments and three low copy nuclear genes (LCNG). Ecological niche modelling was used to examine the potential range shifts of Cardiocrinum in response to climatic change. KEY RESULTS: The molecular data showed high genetic similarity in the cpDNA (98·37 %) and LCNG (94·53 %) sequences. The biogeographical analyses revealed that the ancestor of Cardiocrinum diversified during the late Miocene (approx. 7·32 Mya) in Central China. The ancestor of the C. giganteum lineage dispersed westward to the Himalayas and south-west China with the split between C. giganteum and C. giganteum var. yunnanense occurring around 4·11 Mya consistent with the period of orogeny of the Hengduan Mountains. Some populations of the C. cathayanum lineage dispersed eastward to south Japan via the land bridge approx. 4·97 Mya, providing opportunities for allopatric speciation of C. cordatum The predicted suitable habitats of Cardiocrinum have become smaller and more fragmented since the Last Glacial Maximum. CONCLUSIONS: Our study provides evidence of a biogeographical pattern of dispersal from Central China to the Himalayas in the west and Japan in the east for genera distributed across the SJFR, and highlights that the orogeny of the Hengduan Mountains and fluctuations of the sea level of the East China Sea played important roles in promoting species divergence.

Phylogenetic reappraisal of Allium subgenus Cyathophora (Amaryllidaceae) and related taxa, with a proposal of two new sections.

The phylogeny of subgenus Cyathophora and representatives of its closely related taxa within Allium were reconstructed based on nrDNA ITS and two plastid fragments (trnL-F and rpl32-trnL). The constructed phylogenies indicated that subgenus Cyathophora was not monophyletic and to be split in three parts positioned in different clusters. Allium kingdonii was unequivocally placed within subgenus Amerallium and formed an immediate sister relationship with New World Amerallium clade, suggesting an unexpected intercontinental disjunct distribution. For another, Allium trifurcatum was firmly nested within subgenus Butomissa next to A. tuberosum and A. ramosum, but it is distinctly different morphologically from the latter by thinly leathery bulb tunics, uniovulate locule and obviously 3-cleft stigma. Based on the geographic features, morphological and molecular evidences, two new sections, Kingdonia X.J.He et D.Q.Huang for A. kingdonii and Trifurcatum X.J.He et D.Q.Huang for A. trifurcatum, were proposed. The remaining three species of subgenus Cyathophora formed a well-defined clade, and the phylogenetic relationships among them recovered were consistent with previous findings. In addition, A. weschniakowii and A. subtilissimum were proven to be a member of subgenera Rhizirideum sensu stricto (s. str.) and Cepa, respectively, rather than subgenera Cepa and Polyprason previously proposed. Section Rhizomatosa represented by A. caespitosum should be subsumed within section Caespitosoprason of subgenus Rhizirideum s. str.

Phylogeny and Taxonomic Revision of the Genus Melanosciadium (Apiaceae), Based on Plastid Genomes and Morphological Evidence.

Melanosciadium is considered a monotypic genus and is also endemic to the southwest of China. No detailed phylogenetic studies or plastid genomes have been identified in Melanosciadium. In this study, the plastid genome sequence and nrDNA sequence were used for the phylogenetic analysis of Melanosciadium and its related groups. Angelica tsinlingensis was previously considered a synonym of Hansenia forbesii. Similarly, Ligusticum angelicifolium was previously thought to be the genus Angelica or Ligusticopsis. Through field observations and morphological evidence, we believe that the two species are more similar to M. pimpinelloideum in leaves, umbel rays, and fruits. Meanwhile, we found a new species from Anhui Province (eastern China) that is similar to M. pimpinelloideum and have named it M. Jinzhaiensis. We sequenced and assembled the complete plastid genomes of these species and another three Angelica species. The genome comparison results show that M. pimpinelloideum, A. tsinlingensis, Ligusticum angelicifolium, and M. jinzhaiensis have similarities to each other in the plastid genome size, gene number, and length of the LSC and IR regions; the plastid genomes of these species are distinct from those of the Angelica species. In addition, we reconstruct the phylogenetic relationships using both plastid genome sequences and nrDNA sequences. The phylogenetic analysis revealed that A. tsinlingensis, M. pimpinelloideum, L. angelicifolium, and M. jinzhaiensis are closely related to each other and form a monophyletic group with strong support within the Selineae clade. Consequently, A. tsinlingensis and L. angelicifolium should be classified as members of the genus Melanosciadium, and suitable taxonomical treatments have been proposed. Meanwhile, a comprehensive description of the new species, M. jinzhaiensis, is presented, encompassing its habitat environment and detailed morphological traits.

Plastid phylogenomics contributes to the taxonomic revision of taxa within the genus Sanicula L. and acceptance of two new members of the genus.

Introduction: The genus Sanicula L. is a taxonomically complicated taxa within Apiaceae, as its high variability in morphology. Although taxonomists have performed several taxonomic revisions for this genus, the interspecific relationships and species boundaries have not been satisfactorily resolved, especially for those endemic to China. This study mainly focused on S. giraldii var. ovicalycina, S. tienmuensis var. pauciflora, and S. orthacantha var. stolonifera and also described two new members of the genus. Methods: We newly sequenced sixteen plastomes from nine Sanicula species. Combined with eleven plastomes previously reported by us and one plastome downloaded, we performed a comprehensively plastid phylogenomics analysis of 21 Sanicula taxa. Results and Discussion: The comparative results showed that 21 Sanicula plastomes in their structure and features were highly conserved and further justified that two new species were indeed members of Sanicula. Nevertheless, eleven mutation hotspot regions were still identified. Phylogenetic analyses based on plastome data and the ITS sequences strongly supported that these three varieties were clearly distant from three type varieties. The results implied that these three varieties should be considered as three independent species, which were further justified by their multiple morphological characters. Therefore, revising these three varieties into three independent species was reasonable and convincing. Moreover, we also identified and described two new Sanicula species (S. hanyuanensis and S. langaoensis) from Sichuan and Shanxi, China, respectively. Based on their distinct morphological characteristics and molecular phylogenetic analysis, two new species were included in Sanicula. In summary, our study impelled the revisions of Sanicula members and improved the taxonomic system of the genus.

Evolutionary events in Lilium (including Nomocharis, Liliaceae) are temporally correlated with orogenies of the Q-T plateau and the Hengduan Mountains.

The Hengduan Mountains (H-D Mountains) in China flank the eastern edge of the Qinghai-Tibet Plateau (Q-T Plateau) and are a center of great temperate plant diversity. The geological history and complex topography of these mountains may have prompted the in situ evolution of many diverse and narrowly endemic species. Despite the importance of the H-D Mountains to biodiversity, many uncertainties remain regarding the timing and tempo of their uplift. One hypothesis is that the Q-T Plateau underwent a final, rapid phase of uplift 8-7 million years ago (Mya) and that the H-D Mountains orogeny was a separate event occurring 4-3 Mya. To evaluate this hypothesis, we performed phylogenetic, biogeographic, divergence time dating, and diversification rate analyses of the horticulturally important genus Lilium, including Nomocharis. The Lilium-Nomocharis complex is distributed throughout the temperate Northern Hemisphere but is most diverse within the H-D Mountains and Q-T Plateau. Our matK and ITS phylogenies support previous studies showing that Nomocharis is nested within Lilium. However, we detected incongruence between the two gene trees which may result from hybridization. Dating analyses performed using the ITS dataset showed that the evolution of major lineages within Lilium-Nomocharis may be temporally coincident with Q-T Plateau uplift occurring 8-7 Mya and H-D Mountains uplift approximately 4-3 Mya. Our analyses of diversification times and rates among Lilium-Nomocharis clades are less conclusive. However, these do suggest high extinction rates among H-D Mountains lineages.

Backbone phylogeny and adaptive evolution of Pleurospermum s. l.: New insights from phylogenomic analyses of complete plastome data.

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.