Search

Divergent patterns of phenotypic selection on floral traits can arise in response to interactions with functionally distinct pollinators. However, there are a limited number of studies that relate patterns of phenotypic selection on floral traits to variation in local pollinator assemblages in pollination-generalized plant species. We studied phenotypic selection on floral traits of Viscaria vulgaris, a plant that interacts with a broad range of diurnal and nocturnal pollinators, and related divergence in phenotypic selection on floral traits to the expected level of divergence in local pollinator assemblages. We detected phenotypic selection on floral traits involved in the attraction of pollinators and the mechanics of pollen removal and deposition, and demonstrated that floral traits are subject to spatiotemporal variation in the strength and direction of phenotypic selection. We revealed that diurnal and nocturnal pollinators, when considered in isolation, mediated divergent patterns of phenotypic selection on floral traits. Consistent with the Grant-Stebbins model, we observed that divergence in phenotypic selection on floral traits increased with the expected level of divergence in local pollinator assemblages. Thus, generalized plant-pollinator interactions can mediate phenotypic selection on floral traits, and distinct local pollinator assemblages can generate a geographic mosaic of divergent patterns of phenotypic selection. We underscore that these outcomes are not exclusive to specialized plant-pollinator interactions and can emerge at a local geographic scale., (© The Author(s) 2024. Published by Oxford University Press on behalf of the European Society of Evolutionary Biology.)

Background: Glehnia littoralis F. Schmidt ex Miq., an endangered plant species with significant medicinal, edible, and ecological value, is now a central concern for conservation and sustainable utilization. Investigating the physiological and ecological mechanisms leading to its endangerment and elucidating its genetic background constitutes the foundation for conducting in-depth research on G. littoralis., Results: Our observations have revealed a significant degree of floral sterility in wild populations of G. littoralis. The inflorescences of G. littoralis are classified into three types: completely fertile, completely sterile, and partially fertile compound umbels. Moreover, the flowers of G. littoralis can be categorized into fertile and sterile types. Sterile flowers exhibited abnormalities in the stigma, ovary, and ovules. This study is the first to discover that the presence or absence of a giant cell at the funiculus during the initiation of ovule primordium determines whether the flower can develop normally, providing cytological evidence for female sterility in G. littoralis. Conversely, both fertile and sterile flowers produced normally developed pollen. Field observations have suggested that robust plants bear more fertile umbels, while weaker ones have fewer or even no fertile umbels, indicating a close relationship between flower fertility and plant nutritional status. Our model correctly predicted that the eastern coastal regions of China, as well as prospective areas in Neimenggu and Sichuan, are suitable environments for its cultivation. Additionally, Using flow cytometry and genome survey, we estimated the genome size of G. littoralis to be 3.06 Gb and the heterozygosity to be 4.58%., Conclusion: The observations and findings presented in this study were expected to provide valuable insights for further conserving its genetic resources and sustainable utilization of G. littoralis., (© 2024. The Author(s).)

Background: The genus Hydrocotyle Tourn. ex L. is a key group for further study on the evolution of Apiales, comprising around 170 species globally. Previous studies mainly focused on separate sections and provided much information about this genus, but its infrageneric relationships are still confusing. In addition, the genetic basis of its adaptive evolution remains poorly understood. To investigate the phylogeny and evolution of the genus, we selected ten representative species covering two of three diversity distribution centers and exhibiting rich morphology diversity. Comparative plastome analysis was conducted to clarify the structural character of Hydrocotyle plastomes. Positive selection analyses were implemented to assess the evolution of the genus. Phylogenetic inferences with protein-coding sequences (CDS) of Hydrocotyle and 17 related species were also performed., Results: Plastomes within Hydrocotyle were generally conservative in structure, gene order, and size. A total of 14 regions (rps16-trnK, trnQ-rps16, atpI-atpH, trnC-petN-psbM, ycf3-trnS, accD-psaI-ycf4, petA-psbJ, rps12-rpl20, rpl16 intron, rps3-rpl16 intron, rps9-rpl22, ndhF-rpl32, ndhA intron, and ycf1a) were recognized as hotspot regions within the genus, which suggested to be promising DNA barcodes for global phylogenetic analysis of Hydrocotyle. The ycf15 gene was suggested to be a protein-coding gene for Hydrocotyle species, and it could be used as a DNA barcode to identify Hydrocotyle. In phylogenetic analysis, three monophyletic clades (Clade I, II, III) were identified with evidence of rapid radiation speciation within Clade I. The selective pressure analysis detected that six CDS genes (ycf1b, matK, atpF, accD, rps14, and psbB) of Hydrocotyle species were under positive selection. Within the genus, the last four genes were conservative, suggesting a relation to the unique evolution of the genus in Apiales. Seven genes (atpE, matK, psbH, ycf1a, ycf1b, rpoA, and ycf2) were detected to be under some degree of positive selection in different taxa within the genus Hydrocotyle, indicating their role in the adaptive evolution of species., Conclusions: Our study offers new insights into the phylogeny and adaptive evolution of Hydrocotyle. The plastome sequences could significantly enhance phylogenetic resolution and provide genomic resources and potential DNA markers useful for future studies of the genus., (© 2024. The Author(s).)

Complex coumarins (CCs) represent characteristic metabolites found in Apiaceae plants, possessing significant medical value. Their essential functional role is likely as protectants against pathogens and regulators responding to environmental stimuli. Utilizing genomes and transcriptomes from 34 Apiaceae plants, including our recently sequenced Peucedanum praeruptorum, we conduct comprehensive phylogenetic analyses to reconstruct the detailed evolutionary process of the CC biosynthetic pathway in Apiaceae. Our results show that three key enzymes - p-coumaroyl CoA 2'-hydroxylase (C2'H), C-prenyltransferase (C-PT), and cyclase - originated successively at different evolutionary nodes within Apiaceae through various means of gene duplications: ectopic and tandem duplications. Neofunctionalization endows these enzymes with novel functions necessary for CC biosynthesis, thus completing the pathway. Candidate genes are cloned for heterologous expression and subjected to in vitro enzymatic assays to test our hypothesis regarding the origins of the key enzymes, and the results precisely validate our evolutionary inferences. Among the three enzymes, C-PTs are likely the primary determinant of the structural diversity of CCs (linear/angular), due to divergent activities evolved to target different positions (C-6 or C-8) of umbelliferone. A key amino acid variation (Ala161/Thr161) is identified and proven to play a crucial role in the alteration of enzymatic activity, possibly resulting in distinct binding forms between enzymes and substrates, thereby leading to different products. In conclusion, this study provides a detailed trajectory for the establishment and evolution of the CC biosynthetic pathway in Apiaceae. It explains why only a portion, not all, of Apiaceae plants can produce CCs and reveals the mechanisms of CC structural diversity among different Apiaceae plants., (© 2024. The Author(s).)

Ten SSR markers based on transcriptome sequencing were employed to genotype 231 samples of G. littoralis subsp. littoralis (Apiaceae) from nine cultivated populations and seven wild populations, aiming to assess the genetic diversity and genetic structure, and elucidate the origin of the cultivated populations. Cultivated populations exhibited relatively high genetic diversity (h = 0.441, I = 0.877), slightly lower than that of their wild counterparts (h = 0.491, I = 0.930), likely due to recent domestication and ongoing gene flow between wild and cultivated germplasm. The primary cultivated population in Shandong have the crucial genetic status. A single origin of domestication was inferred through multiple analysis, and wild populations from Liaoning and Shandong are inferred to be potentially the ancestor source for the present cultivated populations. Phenotypic analysis revealed a relatively high heritability of root length across three growth periods (0.683, 0.284, 0.402), with significant correlations observed between root length and petiole length (Pearson correlation coefficient = 0.30, P<0.05), as well as between root diameter and leaf area (Pearson correlation coefficient = 0.36, P<0.01). These parameters can serve as valuable indicators for monitoring the developmental progress of medicinal plants during field management. In summary, this study can shed light on the intricate genetic landscape of G. littoralis subsp. littoralis, providing foundational insights crucial for conservation strategies, targeted breeding initiatives, and sustainable management practices in both agricultural and natural habitats., Competing Interests: The authors declare no competing interests., (Copyright: © 2024 Li et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Background: Tongoloa is a genus comprising approximately 20 species, primarily distributed in the mountainous regions of southwest China. The insufficiency of specimen materials and morphological similarities among species render it a taxonomically challenging genus within the Apiaceae family. To elucidate the phylogenetic relationships and taxonomy of Chinese Tongoloa, this study utilized a total of 115 nrITS sequences, including 47 recently obtained sequences, for phylogenetic reconstruction., Results: Phylogenetic relationships reconstructed from ITS sequences indicate that the East Asia Clade and the Komarovia Clade are sister groups, and Tongoloa belongs to the East Asia Clade. Species of Tongoloa are subdivided into 3 distinct groups, all sharing similar fruit morphologies and are clearly differentiated from related taxa. Several Tongoloa-like members classified under other genera are interpreted to be closely related to Tongoloa. Morphological and molecular data indicate that Tongoloa, Sinolimprichtia subclade and Chinese Trachydium subclade are separate yet genetically contiguous taxa. It is confirmed that Tongoloa zhongdianensis belongs to the Hymenidium Clade, while Sinocarum is classified within the Acronema Clade. Two new taxa are found in the Hengduan Mountains., Conclusion: Tongoloa is a genus within the East Asia Clade of Apiaceae, and the phylogeny reconstructed based on ITS sequences divides it into 3 main groups. By integrating fruit morphology and molecular phylogenetic analyses, we preliminary clarified the intricate taxonomic relationships among Tongoloa and related taxa. These results provide valuable opportunities for a deeper understanding of the phylogeny of Tongoloa., (© 2024. The Author(s).)

Background: Powdery mildew, caused by Eeysiphe heraclei, seriously threatens Heracleum moellendorffii Hance. Plant secondary metabolites are essential to many activities and are necessary for defense against biotic stress. In order to clarify the functions of these metabolites in response to the pathogen, our work concentrated on the variations in the accumulation of secondary metabolites in H. moellendorffii during E. heraclei infection., Results: Following E. heraclei infection, a significant upregulation of coumarin metabolites-particularly simple coumarins and associated genes was detected by RNA-seq and UPLC-MS/MS association analysis. Identifying HmF6'H1, a Feruloyl CoA 6'-hydroxylase pivotal in the biosynthesis of the coumarin basic skeleton through ortho-hydroxylation, was a significant outcome. The cytoplasmic HmF6'H1 protein was shown to be able to catalyze the ortho-hydroxylation of p-coumaroyl-CoA and caffeoyl-CoA, resulting in the formation of umbelliferone and esculetin, respectively. Over-expression of the HmF6'H1 gene resulted in increased levels of simple coumarins, inhibiting the biosynthesis of furanocoumarins and pyranocoumarins by suppressing PT gene expression, enhancing H. moellendorffii resistance to powdery mildew., Conclusions: These results established HmF6'H1 as a resistance gene aiding H. moellendorffii in combatting E. heraclei infection, offering additional evidence of feruloyl-CoA 6'-hydroxylase role in catalyzing various types of simple coumarins. Therefore, this work contributes to our understanding of the function of simple coumarins in plants' defense against powdery mildew infection., (© 2024. The Author(s).)

Background: Peucedani Radix, also known as "Qian-hu" is a traditional Chinese medicine derived from Peucedanum praeruptorum Dunn. It is widely utilized for treating wind-heat colds and coughs accompanied by excessive phlegm. However, due to morphological similarities, limited resources, and heightened market demand, numerous substitutes and adulterants of Peucedani Radix have emerged within the herbal medicine market. Moreover, Peucedani Radix is typically dried and sliced for sale, rendering traditional identification methods challenging., Materials and Methods: We initially examined and compared 104 commercial "Qian-hu" samples from various Chinese medicinal markets and 44 species representing genuine, adulterants or substitutes, utilizing the mini barcode ITS2 region to elucidate the botanical origins of the commercial "Qian-hu". The nucleotide signature specific to Peucedani Radix was subsequently developed by analyzing the polymorphic sites within the aligned ITS2 sequences., Results: The results demonstrated a success rate of 100% and 93.3% for DNA extraction and PCR amplification, respectively. Forty-five samples were authentic "Qian-hu", while the remaining samples were all adulterants, originating from nine distinct species. Peucedani Radix, its substitutes, and adulterants were successfully identified based on the neighbor-joining tree. The 24-bp nucleotide signature (5'-ATTGTCGTACGAATCCTCGTCGTC-3') revealed distinct differences between Peucedani Radix and its common substitutes and adulterants. The newly designed specific primers (PR-F/PR-R) can amplify the nucleotide signature region from commercial samples and processed materials with severe DNA degradation., Conclusions: We advocate for the utilization of ITS2 and nucleotide signature for the rapid and precise identification of herbal medicines and their adulterants to regulate the Chinese herbal medicine industry., (© 2024. The Author(s), under exclusive licence to Springer Nature B.V.)

Cryptotaenia japonica, a traditional medicinal and edible vegetable crops, is well-known for its attractive flavors and health care functions. As a member of the Apiaceae family, the evolutionary trajectory and biological properties of C. japonica are not clearly understood. Here, we first reported a high-quality genome of C. japonica with a total length of 427 Mb and N50 length 50.76 Mb, was anchored into 10 chromosomes, which confirmed by chromosome (cytogenetic) analysis. Comparative genomic analysis revealed C. japonica exhibited low genetic redundancy, contained a higher percentage of single-cope gene families. The homoeologous blocks, K s , and collinearity were analyzed among Apiaceae species contributed to the evidence that C. japonica lacked recent species-specific WGD. Through comparative genomic and transcriptomic analyses of Apiaceae species, we revealed the genetic basis of the production of anthocyanins. Several structural genes encoding enzymes and transcription factor genes of the anthocyanin biosynthesis pathway in different species were also identified. The CjANSa, CjDFRb, and CjF3H gene might be the target of Cjaponica&#95;2.2062 (bHLH) and Cjaponica&#95;1.3743 (MYB). Our findings provided a high-quality reference genome of C. japonica and offered new insights into Apiaceae evolution and biology., (© 2024 Society for Experimental Biology and John Wiley & Sons Ltd.)

Transcriptome sequencing was employed to mine the simple sequence repeat(SSR) locus information of Saposhnikovia divaricata and design specific primers, which aimed to provide a basis for the research on the genetic diversity of S. divaricata germplasm resources. The seed purity, 1 000-seed weight, germination rate, and seed vigor were determined. MISA was used to obtain the SSR locus information from 12 606 unigene longer than 1 kb in the transcriptome database. Forty-three pairs of SSR primers designed in Primer 3 were used to analyze the polymorphism of 28 S. divaricata samples of different sources. The results showed that there were differences in the seed purity, 1 000-seed weight, germination rate, vigor, and seed length and width among S. divaricata samples of different sources. Particularly, the germination rate and seed vigor had significant differences, and HB-ZJK1, NMG-CF4, NMG-BT, NMG-HLE1, and NMG-CF2 had significantly higher 1 000-seed weight, germination rate, and seed vigor than the samples of other sources. Among the 86 233 unigene, 12 606(14.62%) unigene contained 15 958 SSR loci, with one SSR locus every 5 009 bp on average. The SSR loci were mainly single nucleotide and dinucleotide repeats, which were dominated by G/C and TC/AG, respectively. All the primers were screened by using 28 S. divaricata sample from different habitats, and the primers corresponding to the amplification products with clear bands and stable polymorphism were obtained. The clustering results of the biological characteristics and genetic diversity of the 28 S. divaricata samples were basically consistent, and the samples of the same origin(HB-AG1, HB-AG2, HB-ZJK1, and HB-ZJK2) generally gathered together and had close genetic relationship. The SSRs in S. divaricata transcriptome has high frequency, rich types, and high polymorphism, which provides candidate molecular markers for the germplasm identification, genetic map construction, and molecular-assisted breeding.

Background: Traditional Chinese medicine has used Peucedanum praeruptorum Dunn (Apiaceae) for a long time. Various coumarins, including the significant constituents praeruptorin (A-E), are the active constituents in the dried roots of P. praeruptorum. Previous transcriptomic and metabolomic studies have attempted to elucidate the distribution and biosynthetic network of these medicinal-valuable compounds. However, the lack of a high-quality reference genome impedes an in-depth understanding of genetic traits and thus the development of better breeding strategies., Results: A telomere-to-telomere (T2T) genome was assembled for P. praeruptorum by combining PacBio HiFi, ONT ultra-long, and Hi-C data. The final genome assembly was approximately 1.798 Gb, assigned to 11 chromosomes with genome completeness >98%. Comparative genomic analysis suggested that P. praeruptorum experienced 2 whole-genome duplication events. By the transcriptomic and metabolomic analysis of the coumarin metabolic pathway, we presented coumarins' spatial and temporal distribution and the expression patterns of critical genes for its biosynthesis. Notably, the COSY and cytochrome P450 genes showed tandem duplications on several chromosomes, which may be responsible for the high accumulation of coumarins., Conclusions: A T2T genome for P. praeruptorum was obtained, providing molecular insights into the chromosomal distribution of the coumarin biosynthetic genes. This high-quality genome is an essential resource for designing engineering strategies for improving the production of these valuable compounds., (© The Author(s) 2024. Published by Oxford University Press GigaScience.)

Peucedanum praeruptorum Dunn, a traditional Chinese medicine rich in coumarin, belongs to the Apiaceae family. A high-quality assembled genome of P. praeruptorum is lacking, which has posed obstacles to functional identification and molecular evolution studies of genes associated with coumarin production. Here, a chromosome-scale reference genome of P. praeruptorum, an important medicinal and aromatic plant, was first sequenced and assembled using Oxford Nanopore Technologies and Hi-C sequencing. The final assembled genome size was 1.83 Gb, with a contig N50 of 11.12 Mb. The entire BUSCO evaluation and second-generation read comparability rates were 96.0 % and 99.31 %, respectively. Furthermore, 99.91 % of the genome was anchored to 11 pseudochromosomes. The comparative genomic study revealed the presence of 18,593 orthogroups, which included 476 species-specific orthogroups and 1211 expanded gene families. Two whole-genome duplication (WGD) events and one whole-genome triplication (WGT) event occurred in P. praeruptorum. In addition to the γ-WGT shared by core eudicots or most eudicots, the first WGD was shared by Apiales, while the most recent WGD was unique to Apiaceae. Our study demonstrated that WGD events that occurred in Apioideae highlighted the important role of tandem duplication in the biosynthesis of coumarins and terpenes in P. praeruptorum. Additionally, the expansion of the cytochrome P450 monooxygenase, O-methyltransferase, ATP-binding cassette (ABC) transporter, and terpene synthase families may be associated with the abundance of coumarins and terpenoids. Moreover, we identified >170 UDP-glucosyltransferase members that may be involved in the glycosylation post-modification of coumarins. Significant gene expansion was observed in the ABCG, ABCB, and ABCC subgroups of the ABC transporter family, potentially facilitating the transmembrane transport of coumarins after bolting. The P. praeruptorum genome provides valuable insights into the machinery of coumarin biosynthesis and enhances our understanding of Apiaceae evolution., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023. Published by Elsevier B.V.)

The genus Apium, belonging to the family Apiaceae, comprises roughly 20 species. Only two species, Apium graveolens and Apium leptophyllum, are available in China and are both rich in nutrients and have favorable medicinal properties. However, the lack of genomic data has severely constrained the study of genetics and evolution in Apium plants. In this study, Illumina NovaSeq 6000 and Nanopore sequencing platforms were employed to identify the mitochondrial genomes of A. graveolens and A. leptophyllum . The complete lengths of the mitochondrial genomes of A. graveolens and A. leptophyllum were 263,017 bp and 260,164 bp, respectively, and contained 39 and 36 protein-coding genes, five and six rRNA genes, and 19 and 20 tRNA genes. Consistent with most angiosperms, both A. graveolens and A. leptophyllum showed a preference for codons encoding leucine (Leu). In the mitochondrial genome of A. graveolens , 335 SSRs were detected, which is higher than the 196 SSRs found in the mitochondrial genome of A. leptophyllum . Studies have shown that the most common RNA editing type is C-to-U, but, in our study, both A. graveolens and A. leptophyllum exhibited the U-C editing type. Furthermore, the transfer of the mitochondrial genomes of A. graveolens and A. leptophyllum into the chloroplast genomes revealed homologous sequences, accounting for 8.14% and 4.89% of the mitochondrial genome, respectively. Lastly, in comparing the mitochondrial genomes of 29 species, it was found that A. graveolens , A. leptophyllum , and Daucus carota form a sister group with a support rate of 100%. Overall, this investigation furnishes extensive insights into the mitochondrial genomes of A. graveolens and A. leptophyllum , thereby enhancing comprehension of the traits and evolutionary patterns within the Apium genus. Additionally, it offers supplementary data for evolutionary and comparative genomic analyses of other species within the Apiaceae family.

Background: Peucedanum praeruptorum Dunn, a traditional Chinese herbal medicine, contains coumarin and volatile oil components that have clinical application value. However, early bolting often occurs in the medicinal materials of Apiaceae plants. The rhizomes of the medicinal parts are gradually lignified after bolting, resulting in a sharp decrease in the content of coumarins. At present, the link between coumarin biosynthesis and early bolting in P. praeruptorum has not been elucidated., Results: Combining the genome sequencing and the previous transcriptome sequencing results, we reanalyzed the differential transcripts of P. praeruptorum before and after bolting. A total of 62,088 new transcripts were identified, of which 31,500 were unknown transcripts. Functional classification and annotation showed that many genes were involved in the regulation of transcription, defense response, and carbohydrate metabolic processes. The main domains are the pentatricopeptide repeat, protein kinase, RNA recognition motif, leucine-rich repeat, and ankyrin repeat domains, indicating their pivotal roles in protein modification and signal transduction. Gene structure analysis showed that skipped exon (SE) was the most dominant alternative splicing, followed by the alternative 3' splice site (A3SS) and the alternative 5' splice site (A5SS). Functional enrichment of differentially expressed genes showed that these differentially expressed genes mainly include transmembrane transporters, channel proteins, DNA-binding proteins, polysaccharide-binding proteins, etc. In addition, genes involved in peroxisome, hexose phosphate pathway, phosphatidylinositol signaling system, and inositol phosphate metabolism pathway were greatly enriched. A protein-protein interaction network analysis discoverd 1,457 pairs of proteins that interact with each other. The expression levels of six UbiA genes, three UGT genes, and four OMT genes were higher during the bolting stage. This observation suggests their potential involvement in the catalytic processes of prenylation, glycosylation, and methylation of coumarins, respectively. A total of 100 peroxidase (PRX) genes were identified being involved in lignin polymerization, but only nine PRX genes were highly expressed at the bolting stage. It is worth noting that 73 autophagy-related genes (ATGs) were first identified from the KEGG pathway-enriched genes. Some ATGs, such as BHQH00009837, BHQH00013830, and novel8944, had higher expression levels after bolting., Conclusions: Comparative transcriptome analysis and large-scale genome screening provide guidance and new opinions for the identification of bolting-related genes in P. praeruptorum., (© 2023. BioMed Central Ltd., part of Springer Nature.)

Background: Sium L. (Apiaceae) is a small genus distributed primarily in Eurasia, with one species also occurring in North America. Recently, its circumscription has been revised to include 10 species, however, the phylogenetic relationships within its two inclusive clades were poorly supported or collapsed in previous studies based on nuclear ribosomal DNA ITS or cpDNA sequences. To identify molecular markers suitable for future intraspecific phylogeographic and population genetic studies, and to evaluate the efficacy of plastome in resolving the phylogenetic relationships of the genus, the complete chloroplast (cp) genomes of six Sium species were sequenced., Results: The Sium plastomes exhibited typical quadripartite structures of Apiaceae and most other higher plant plastid DNAs, and were relatively conserved in their size (153,029-155,006 bp), gene arrangement and content (with 114 unique genes). A total of 61-67 SSRs, along with 12 highly divergent regions (trnQ, trnG-atpA, trnE-trnT, rps4-trnT, accD-psbI, rpl16, ycf1-ndhF, ndhF-rpl32, rpl32-trnL, ndhE-ndhG, ycf1a and ycf1b) were discovered in the plastomes. No significant IR length variation was detected showing that plastome evolution was conserved within this genus. Phylogenomic analysis based on whole chloroplast genome sequences produced a highly resolved phylogenetic tree, in which the monophyly of Sium, as well as the sister relationship of its two inclusive clades were strongly supported., Conclusions: The plastome sequences could greatly improve phylogenetic resolution, and will provide genomic resources and potential markers useful for future studies of the genus., (© 2023. The Author(s).)

Background: The Sinodielsia clade of the subfamily Apioideae (Apiacieae) was established in 2008, and it is composed of 37 species from 17 genera. Its circumscription is still poorly delimited and unstable, and interspecific relationships in the clade lack comprehensive analysis. Chloroplast (cp.) genomes provide valuable and informative data sources for evolutionary biology and have been widely used in studies on plant phylogeny. To infer the phylogenetic history of the Sinodielsia clade, we assembled complete cp. genomes of 39 species and then performed phylogenetic analysis based on these cp. genome sequence data combined with 66 published cp. genomes from 16 genera relative to the Sinodielsia clade., Results: These 39 newly assembled genomes had a typical quadripartite structure with two inverted repeat regions (IRs: 17,599-31,486 bp) separated by a large single-copy region (LSC: 82,048-94,046 bp) and a small single-copy region (SSC: 16,343-17,917 bp). The phylogenetic analysis showed that 19 species were clustered into the Sinodielsia clade, and they were divided into two subclades. Six mutation hotspot regions were detected from the whole cp. genomes among the Sinodielsia clade, namely, rbcL-accD, ycf4-cemA, petA-psbJ, ycf1-ndhF, ndhF-rpl32 and ycf1, and it was found that ndhF-rpl32 and ycf1 were highly variable in the 105 sampled cp. genomes., Conclusion: The Sinodielsia clade was subdivided into two subclades relevant to geographical distributions, except for cultivated and introduced species. Six mutation hotspot regions, especially ndhF-rpl32 and ycf1, could be used as potential DNA markers in the identification and phylogenetic analyses of the Sinodielsia clade and Apioideae. Our study provided new insights into the phylogeny of the Sinodielsia clade and valuable information on cp. genome evolution in Apioideae., (© 2023. The Author(s).)

Hybridization is recognized as a major force in species evolution and biodiversity formation, generally leading to the origin and differentiation of new species. Multiple hybridization events cannot easily be reconstructed, yet they offer the potential to study a number of evolutionary processes. Here, we used nuclear expressed sequence tag-simple sequence repeat and large-scale single nucleotide polymorphism variation data, combined with niche analysis, to investigate the putative independent hybridization events in Notopterygium, a group of perennial herb plants endemic to China. Population genomic analysis indicated that the four studied species are genetically well-delimited and that N. forrestii and N. oviforme have originated by hybridization. According to Approximate Bayesian Computation, the best-fit model involved the formation of N. forrestii from the crossing of N. franchetii and N. incisum, with N. forrestii further backcrossing to N. franchetii to form N. oviforme. The niche analyses indicated that niche divergence [likely triggered by the regional climate changes, particularly the intensification of East Asian winter monsoon, and tectonic movements (affecting both Qinghai-Tibetan Plateau and Qinling Mountains)] may have promoted and maintained the reproductive isolation among hybrid species. N. forrestii shows ecological specialization with respect to their parental species, whereas N. oviforme has completely shifted its niche. These results suggested that the climate and environmental factors together triggered the two-step hybridization of the East Asia herb plants. Our study also emphasizes the power of genome-wide SNPs for investigating suspected cases of hybridization, particularly unravelling old hybridization events., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Inc. All rights reserved.)

The current study examined the phylogenetic pattern of medicinal species of the family Apiaceae based on flavonoid groups production, as well as the overall mechanism of the key genes involved in flavonol and flavone production. Thirteen species of the family Apiaceae were used, including Eryngium campestre from the subfamily Saniculoideae, as well as Cuminum cyminum, Carum carvi, Coriandrum sativum, Apium graveolens, Petroselinum crispum, Pimpinella anisum, Anethum graveolens, Foeniculum vulgare, Daucus carota, Ammi majus, Torilis arvensis, and Deverra tortuosa from the subfamily Apioideae. The seeds were cultivated, and the leaves were collected to estimate flavonoids and their groups, physiological factors, transcription levels of flavonol and flavone production-related genes. The phylogenetic relationship between the studied species was established using the L-ribosomal 16 (rpl16) chloroplast gene. The results revealed that the studied species were divided into two patterns: six plant species, E. campestre, C. carvi, C. sativum, P. anisum, An. graveolens, and D. carota, contained low content of flavonoids, while the other seven species had high content. This pattern of flavonoids production coincided with the phylogenetic relationships between the studied species. In contrast, the phylogeny of the flavonol and flavone synthase genes was incompatible with the quantitative production of their products. The study concluded that the increment in the production of flavonol depends on the high expression of chalcone synthase, chalcone isomerase, flavanone 3 hydroxylase, flavonol synthase, the increase of Abscisic acid, sucrose, and phenyl ammonia lyase, while flavone mainly depends on evolution and on the high expression of the flavone synthase gene., (© 2023. The Author(s).)

Ligusticopsis litangensis is identified and described as a cryptic species from Sichuan Province, China. Although the distribution of this cryptic species overlaps with that of Ligusticopsis capillacea and Ligusticopsis dielsiana , the morphological boundaries between them are explicit and have obviously distinguishable characters. The main distinguishing features of the cryptic species are as follows: long conical multi-branched roots, very short pedicels in compound umbels, unequal rays, oblong-globose fruits, 1-2 vittae per furrow and 3-4 vittae on the commissure. The above-mentioned features differ somewhat from other species within the genus Ligusticopsis , but generally coincide with the morphological boundaries defined for the genus Ligusticopsis . To determine the taxonomic position of L. litangensis , we sequenced and assembled the plastomes of L. litangensis and compared them with the plastomes of 11 other species of the genus Ligusticopsis . Notably, both phylogenetic analyses based on ITS sequences and the complete chloroplast genome robustly supported that three accessions of L. litangensis are monophyletic clade and then nested in Ligusticopsis genus. Moreover, the plastid genomes of 12 Ligusticopsis species, including the new species, were highly conserved in terms of gene order, gene content, codon bias, IR boundaries and SSR content. Overall, the integration of morphological, comparative genomic and phylogenetic evidence indicates that Ligusticopsis litangensis actually represents a new species.

Background: Drought is one of the crucial constraints limiting horticultural plant's production and development around the world. Pugionium cornutum is an annual or biennial xerophyte with strong environmental adaptability and drought resistance; however, the mechanisms with respect to response to drought stress remain largely unclear., Methods: After seedling emergence, the gravimetric method was used to control soil relative water content (SRWC). Drought stress was applied to the six-leaf stage P. cornutum seedlings. The soil water content of different drought stress levels (L) was controlled by gravimetric method as follows: control (L1): 70-75% SRWC; moderate drought level (L2): 40-45% SRWC; severe drought level (L3): 30-35% SRWC, and the water was added to different drought stress levels at about 18:00 p.m. every day. The experiment ended when the leaves of P. cornutum showed severe wilting (10-leaf stage). Samples were harvested and stored at -80 °C for physiological determination, and transcriptomic and proteomic sequencing., Results: Compared with L1, the leaves of P. cornutum seedlings were increasingly wilted after drought treatment; the SRWC of the drought-stress leaves decreased notably while the leaf water potential was rose; the proline, malondialdehyde (MDA) content increased with the continuous drought treatment but peroxidase (POD) activity decreased. Besides, 3,027 differential genes (DGs) and 196 differential proteins (DPs), along with 1,943 DGs and 489 DPs were identified in L2-L1 and L3-L1, respectively. The transcriptome and proteome integrated analysis manifested that only 30 and 70 were commonly regulated both in L2-L1 and L3-L1, respectively. Of which, 24 and 61 DGs or DPs showed the same trend including sHSPs, APX2, GSTU4, CML42, and POD, etc . However, most of DGs or DPs were regulated only at the transcriptome or proteome level mainly including genes encoding signal pathway (PYR1, PYLs, SnRK2J, PLC2, CDPK9/16/29, CML9, MAPKs), transcription factors (WRKYs, DREB2A, NAC055, NAC072, MYB and, HB7) and ion channel transporters (ALMT4, NHX1, NHX2 and TPK2). These genes or proteins were involved in multiple signaling pathways and some important metabolism processes, which offers valuable information on drought-responsive genes and proteins for further study in P. cornutum ., Competing Interests: The authors declare that they have no competing interests., (© 2023 Wang et al.)

Background: Water scarcity has become one of the most prevalent environmental factors adversely affecting plant growth and development. Different species have developed multiple ways of drought resistance. Saposhnikovia divaricata is a commonly used traditional herb in East Asia. However, limited information is available on the drought response of this herb and further clarification of underlying molecular mechanism remains a challenge., Methods and Results: In this study, a comparative transcriptome analysis was firstly conducted to identify the major pathways and candidate genes involved in the drought adaptive response of S. divaricata. The seedlings of S. divaricata were subjected to progressive drought by withholding water for 16 days followed by 8 days of rehydration. Transcriptome analysis identified a total of 89,784 annotated unigenes. The number of differentially expressed genes (DEGs) gradually increased with the deepening of drought and decreased after rehydration. Gene Ontology enrichment analysis and Kyoto Encyclopedia of Genes and Genomes pathway analysis suggested genes related to oxidoreductase activity, carbohydrate metabolism, plant hormone signaling pathway and secondary metabolism were important in drought response of S. divaricata. Specific genes involved in reactive oxygen species scavenging system (POD, Cu/Zn-SOD, APX), abscisic acid and jasmonic acid signaling pathway (PYL4, PP2Cs, JAR1, JAZ) and phenylpropanoid biosynthesis (4CL, CCR, CAD) underwent dynamic alterations under drought and rehydration. Finally, the expression pattern of 12 selected DEGs from the transcriptomic profiling was validated by real-time quantitative PCR., Conclusion: Our study laid a foundation for understanding the stress response of S. divaricata and other Apiaceae family plant at molecular level., (© 2023. The Author(s), under exclusive licence to Springer Nature B.V.)

Early bolting of Peucedanum praeruptorum Dunn severely affects its quality. In this study, we compared with the root structure of P. praeruptorum and its four coumarins content between early bolting (CT) and unbolting (WT) at different growth stages. We found that the proportion of area outside the root cambium (Rs) was higher in the WT plants than in the CT plants and correlated positively with the proximity to the root tip. Furthermore, the content of all four coumarins was also higher in the WT plants relative to the CT plants. In addition, we identified 15,524 differentially expressed genes (DEGs) between the two plant varieties. 11 DEGs are involved in the photoperiod and gibberellin pathways that regulate early bolting and 24 genes involved in coumarins biosynthesis were also identified. Nevertheless, early bolting of P. praeruptorum does affect its quality formation, and further studies are needed to confirm its mechanism., Competing Interests: Declaration of Competing Interest The authors declare no conflict of interests., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Apiaceae plants are used as medicinal herbs, pesticides, spices, and vegetables; thus, accurately identifying Apiaceae species is important. The grassland ecosystem of Heilongjiang Province in northern China has huge reserves of wild Apiaceae plants, but few reports have systematically documented their diversity. In this study, 275 Apiaceae plants of 23 species in 18 genera were collected from this area. We identified Apiaceae species by using nuclear internal transcribed spacer (ITS/ITS2) and psbA-trnH (chloroplast non-coding region) sequences based on experimental data. The identification efficiency of ITS, ITS2 and psbA-trnH sequences was determined and evaluated by sequence alignment and analysis, intraspecific and interspecific genetic distance analyses, and phylogenetic tree construction. ITS, ITS2 could distinguish 21 species from 17 genera of Apiaceae with good identification effect. When identifying species in the Apiaceae family, ITS2 can be used as the core barcode and psbA-trnH can be used as the supplementary barcode. These results can enrich the reference Apiaceae DNA barcode database., (© 2022. The Author(s), under exclusive licence to Springer Nature B.V.)

Peucedanum nanum and P. violaceum are recognized as members of the genus Peucedanum because of their dorsally compressed mericarps with slightly prominent dorsal ribs and narrowly winged lateral ribs. However, these species are not similar to other Peucedanum taxa but resemble Ligusticopsis in overall morphology. To check the taxonomic positions of P. nanum and P. violaceum , we sequenced their complete plastid genome (plastome) sequences and, together with eleven previously published Ligusticopsis plastomes, performed comprehensively comparative analyses. The thirteen plastomes were highly conserved and similar in structure, size, GC content, gene content and order, IR borders, and the patterns of codon bias, RNA editing, and simple sequence repeats (SSRs). Nevertheless, twelve mutation hotspots ( mat K, ndh C, rps 15, rps 8, ycf 2, ccs A- ndh D, pet N- psb M, psb A- trn K, rps 2- rpo C2, rps 4- trn T, trn H- psb A, and ycf 2- trn L) were selected. Moreover, both the phylogenetic analyses based on plastomes and on nuclear ribosomal DNA internal transcribed spacer (ITS) sequences robustly supported that P. nanum and P. violaceum nested in Ligusticopsis , and this was further confirmed by the morphological evidence. Hence, transferring P. nanum and P. violaceum into Ligusticopsis genus is reasonable and convincing, and two new combinations are presented.

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., (© 2022. The Author(s).)

Saposhnikovia divaricata is a commonly used bulk medicinal plant. To explore the key enzyme genes and their expression in the biosynthesis of chromone and coumarin, the key active components, we carried out transcriptome sequencing(Illumina HiSeq) and bioinformatics analysis for the 1-year-old(S1) and 2-year-old(S2) plants of S. divaricata. A total of 40.8 Gb data was obtained. After the sequence assembly via Trinity, 110 732 transcripts and 86 233 unigenes were obtained, which were aligned and annotated with NR, Swiss-Prot, GO, KEGG, and PFAM. Daucus carota and S. divaricata had the highest sequence homology. KEGG pathway enrichment showed that the differentially expressed genes were mainly enriched in plant hormone signal transduction, phenylpropanoid biosynthesis, and flavonoid biosynthesis pathways. A total of 27 differentially expressed unigenes, including 13 enzyme genes, were identified in the pathways related to the synthesis of active ingredients in S. divaricata. Compared with S1 plant, S2 plant showed up-regulated expression of PAL, BGL, C4H, 4CL, CYP98A, CSE, REF, and CCoAOMT and down-regulated expression of CHS, CAD, and COMT. HCT and POD had both up-regulated and down-regulated unigenes. Among them, PAL, C4H, 4CL, BGL, and CHS can be used as candidate genes for the synthesis of the active ingredients in S. divaricata. The four key enzyme genes were verified by RT-qPCR, which showed the results consistent with transcriptome sequencing. This study enriches the genetic information of S. divaricata and provides support for the identification of candidate genes in the biosynthesis of secondary metabolites.

Background: Saposhnikovia divaricata (Turcz.) Schischk. is a perennial herb whose dried roots are commonly used as a source of traditional medicines. To elucidate the organelle-genome-based phylogeny of Saposhnikovia species and the transfer of DNA between organelle genomes, we sequenced and characterised the mitochondrial genome (mitogenome) of S. divaricata., Results: The mitogenome of S. divaricata is a circular molecule of 293,897 bp. The nucleotide composition of the mitogenome is as follows: A, 27.73%; T, 27.03%; C, 22.39%; and G, 22.85. The entire gene content is 45.24%. A total of 31 protein-coding genes, 20 tRNAs and 4 rRNAs, including one pseudogene (rpl16), were annotated in the mitogenome. Phylogenetic analysis of the organelle genomes from S. divaricata and 10 related species produced congruent phylogenetic trees. Selection pressure analysis revealed that most of the mitochondrial genes of related species are highly conserved. Moreover, 2 and 46 RNA-editing sites were found in the chloroplast genome (cpgenome) and mitogenome protein-coding regions, respectively. Finally, a comparison of the cpgenome and the mitogenome assembled from the same dataset revealed 10 mitochondrial DNA fragments with sequences similar to those in the repeat regions of the cpgenome, suggesting that the repeat regions might be transferred into the mitogenome., Conclusions: In this study, we assembled and annotated the mitogenome of S. divaricata. This study provides valuable information on the taxonomic classification and molecular evolution of members of the family Apiaceae., (© 2022. The Author(s).)

Background: The genus Daucus (Apiaceae) comprises about 40 wild species and the cultivated carrot, a crop of great economic and nutritional importance. The rich genetic diversity of wild Daucus species makes them a valuable gene pool for carrot improvement breeding programs. Therefore, it is essential to have good knowledge of the genome structure and relationships among wild Daucus species. To broaden such knowledge, in this research, the nuclear DNA content for 14 Daucus accessions and four closely related species was estimated by flow cytometry and their pollen morphology was analyzed by light and scanning electron microscopy (SEM)., Results: The flow cytometric analysis showed a 3.2-fold variation in the mean 2C values among Daucus taxa, ranging from 0.999 (D. carota subsp. sativus) to 3.228 pg (D. littoralis). Among the outgroup species, the mean 2C values were 1.775-2.882 pg. The pollen grains of Daucus were tricolporate, mainly prolate or perprolate (rarely) in shape, and mainly medium or small (rarely) in size (21.19-40.38 µm), whereas the outgroup species had tricolporate, perprolate-shaped, and medium-sized (26.01-49.86 µm) pollen grains. In the studied taxa, SEM analysis revealed that exine ornamentation was striate, rugulate, perforate, or the ornamentation pattern was mixed. At the time of shedding, all pollen grains were three-celled, as evidenced by DAPI staining. We also found high positive correlations between the length of the polar axis (P) and the length of the equatorial diameter (E) of pollen grains, as well as between P and P/E. However, when comparing cytogenetic information with palynological data, no significant correlations were observed., Conclusions: This study complements the information on the nuclear DNA content in Daucus and provides comprehensive knowledge of the pollen morphology of its taxa. These findings may be important in elucidating the taxonomic relationships among Daucus species and can help in the correct identification of gene bank accessions. In a broader view, they could also be meaningful for the interpretation of evolutionary trends in the genus., (© 2022. The Author(s).)

Background: Kala zeera [Bunium persicum (Boiss.) Fedtsch] is one of the important spice crops of North Western Himalayas with lot of medicinal and culinary values. In spite of having great importance, this crop is under the threat of extinction due to loss of habitat and lack of awareness. The limited availability of the seeds has ultimately increased the economic value of this spice. The upmarket of Kala zeera leads to its adulteration with other black seeds and cumin seeds. The present investigation was undertaken to evaluate polyphenolics and antioxidant properties of Kala zeera genotypes collected from North Western Himalayas and to develop DNA barcodes that can ensure their purity and can also guide in conservation of selected Kala zeera germplasm lines., Methods and Results: Various locations of North Western Himalayas were explored for collecting 31 diverse germplasm lines of Kala zeera. The collected germplasm was maintained at our experimental stations during 2019-2020 and 2020-2021. These genotypes were evaluated for different seed traits and the methanolic extract from Kala zeera seeds was examined for total phenolic content, total flavonoid content, antioxidant activities by DPPH and FRAP. The results revealed significant variation in seed traits, polyphenolic content and antioxidant properties. 100 seed weight ranged from 0.05 to 0.35 g, TPC ranged from 7.5 to 22.56 mg/g, TFC ranged from 0.58 to 4.15 mg/g, antioxidant properties DPPH ranged from 168 to 624.4 μg/ml and FRAP ranged from 0.72 to 6.91 mg/g. Further, three different barcodes (ITS, rbcL and psbA-trnH) were used to reveal the authenticity of selected Kala zeera. MEGA 5 software was used for clustering and the barcodes did clustering based on geographical distribution of Kala zeera germplasm., Conclusion: Based on molecular barcoding, best barcode combination was identified that may discriminate the Kala zeera germplasm vis-a-vis can authenticate their purity. Moreover, the identified DNA barcodes will have significant role in studying the evolutionary biology of Bunium species and will be important for designing a strategy to conserve the selected Kala zeera germplasm lines. The identified genotypes with high phenolic content and antioxidant activity can further be utilized in Kala zeera breeding programmes., (© 2022. The Author(s), under exclusive licence to Springer Nature B.V.)

Saposhnikovia divaricata , a well-known Chinese medicinal herb, is the sole species under the genus Saposhnikovia of the Apiaceae subfamily Apioideae Drude. However, information regarding its genetic diversity and evolution is still limited. In this study, the first complete chloroplast genome (cpDNA) of wild S. divaricata was generated using de novo sequencing technology. Similar to the characteristics of Ledebouriella seseloides , the 147,834 bp-long S. divaricata cpDNA contained a large single copy, a small single copy, and two inverted repeat regions. A total of 85 protein-coding, 8 ribosomal RNA, and 36 transfer RNA genes were identified. Compared with five other species, the non-coding regions in the S. divaricata cpDNA exhibited greater variation than the coding regions. Several repeat sequences were also discovered, namely, 33 forward, 14 reverse, 3 complement, and 49 microsatellite repeats. Furthermore, phylogenetic analysis using 47 cpDNA sequences of Apioideae members revealed that L. seseloides and S. divaricata clustered together with a 100% bootstrap value, thereby supporting the validity of renaming L. seseloides to S. divaricata at the genomic level. Notably, S. divaricata was most closely related to Libanotis buchtormensis , which contradicts previous reports. Therefore, these findings provide a valuable foundation for future studies on the genetic diversity and evolution of S. divaricata .

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., (© 2022. The Author(s).)

Objective: Angelica heterocarpa (Apiaceae) is a wild endemic French species with special conservation interest in the European Union. It belongs to Angelica complex genus which is widespread throughout the north temperate zone, and is sympatric with other congeneric species. The objective of this work is to develop and characterize microsatellite markers as a new tool for understanding the ecology and evolution of Angelica species complex., Results: We identified simple sequence repeat (SSR) regions in a microsatellite-enriched library from A. heterocarpa and A. sylvestris. All 16 selected SSR regions were found to amplify in these species and were highly polymorphic. Marker transferability was validated in A. razulii and A. archangelica. These markers will help us to better understand the evolutionary dynamic between rare endemics and widespread sister species, and be useful for conservation of the endemic species. Moreover, they can provide new tools for studying the numerous traditional medicinal herbs of the Angelica genus., (© 2022. The Author(s).)

The Apiaceae taxon is one of the most important families of flowering plants and includes thousands of species used for food, flavoring, fragrance, medical and industrial purposes. This study had the specific intent of reviewing the main genomics and transcriptomic data available for this family and their use for the constitution of new varieties. This was achieved starting from the description of the main reproductive systems and barriers, with particular reference to cytoplasmic (CMS) and nuclear (NMS) male sterility. We found that CMS and NMS systems have been discovered and successfully exploited for the development of varieties only in Foeniculum vulgare , Daucus carota , Apium graveolens and Pastinaca sativa; whereas, strategies to limit self-pollination have been poorly considered. Since the constitution of new varieties benefits from the synergistic use of marker-assisted breeding in combination with conventional breeding schemes, we also analyzed and discussed the available SNP and SSR marker datasets (20 species) and genomes (8 species). Furthermore, the RNA-seq studies aimed at elucidating key pathways in stress tolerance or biosynthesis of the metabolites of interest were limited and proportional to the economic weight of each species. Finally, by aligning 53 plastid genomes from as many species as possible, we demonstrated the precision offered by the super barcoding approach to reconstruct the phylogenetic relationships of Apiaceae species. Overall, despite the impressive size of this family, we documented an evident lack of molecular data, especially because genomic and transcriptomic resources are circumscribed to a small number of species. We believe that our contribution can help future studies aimed at developing molecular tools for boosting breeding programs in crop plants of the Apiaceae family.

Background: One of the major trends in angiosperm evolution was the shift from woody to herbaceous habit. However, reversals known as derived woodiness have also been reported in numerous, distantly related clades. Among theories evoked to explain the factors promoting the evolution of derived woodiness are moderate climate theory and cavitation theory. The first assumes that woody habit evolves in response to mild climate allowing for prolonged life span, which in turn leads to bigger and woodier bodies. The second sees woodiness as a result of natural selection for higher cavitation resistance in seasonally dry environments. Here, we compare climatic niches of woody and herbaceous, mostly southern African, umbellifers from the Lefebvrea clade to assess whether woody taxa in fact occur in markedly drier habitats. We also calibrate their phylogeny to estimate when derived woodiness evolved. Finally, we describe the wood anatomy of selected woody and herbaceous taxa to see if life forms are linked to any particular wood traits., Results: The evolution of derived woodiness in chamaephytes and phanerophytes as well as the shifts to short-lived annual therophytes in the Lefebvrea clade took place at roughly the same time: in the Late Miocene during a trend of global climate aridification. Climatic niches of woody and herbaceous genera from the Cape Floristic Region overlap. There are only two genera with distinctly different climatic preferences: they are herbaceous and occur outside of the Cape Floristic Region. Therefore, studied herbs have an overall climatic niche wider than their woody cousins. Woody and herbaceous species do not differ in qualitative wood anatomy, which is more affected by stem architecture and, probably, reproductive strategy than by habit., Conclusions: Palaeodrought was likely a stimulus for the evolution of derived woodiness in the Lefebvrea clade, supporting the cavitation theory. The concurrent evolution of short-lived annuals withering before summer exemplifies an alternative solution to the same problem of drought-induced cavitation. Changes of the life form were most likely neither spurred nor precluded by any qualitative wood traits, which in turn are more affected by internode length and probably also reproductive strategy., (© 2021. The Author(s).)

Traditional phylogenies inferred from chloroplast DNA fragments have not obtained a well-resolved evolutionary history for the backbone of Apioideae, the largest subfamily of Apiaceae. In this study, we applied the genome skimming approach of next-generation sequencing to address whether the lack of resolution at the tip of the Apioideae phylogenetic tree is due to limited information loci or the footprint of ancient radiation. A total of 90 complete chloroplast genomes (including 23 newly sequenced genomes and covering 20 major clades of Apioideae) were analyzed (RAxML and MrBayes) to provide a phylogenomic reconstruction of Apioideae. Dating analysis was also implemented using BEAST to estimate the origin and divergence time of the major clades. As a result, the early divergences of Apioideae have been clarified but the relationship among its distally branching clades (Group A) was only partially resolved, with short internal branches pointing to an ancient radiation scenario. Four major clades, Tordyliinae I, Pimpinelleae I, Apieae and Coriandreae, were hypothesized to have originated from chloroplast capture events induced by early hybridization according to the incongruence between chloroplast-based and nrDNA-based phylogenetic trees. Furthermore, the variable and nested distribution of junction positions of LSC (Large single copy region) and IR B (inverted repeat region B) in Group A may reflect incomplete lineage sorting within this group, which possibly contributed to the unclear phylogenetic relationships among these clades inferred from plastome data. Molecular clock analysis revealed the chloroplast capture events mainly occurred during the middle to late Miocene, providing a geological and climate context for the evolution of Apioideae., (Published by Elsevier Inc.)

Premise: The carrot family (Apiaceae) comprises 466 genera, which include many well-known crops (e.g., aniseed, caraway, carrots, celery, coriander, cumin, dill, fennel, parsley, and parsnips). Higher-level phylogenetic relationships among subfamilies, tribes, and other major clades of Apiaceae are not fully resolved. This study aims to address this important knowledge gap., Methods: Target sequence capture with the universal Angiosperms353 probe set was used to examine phylogenetic relationships in 234 genera of Apiaceae, representing all four currently recognized subfamilies (Apioideae, Azorelloideae, Mackinlayoideae, and Saniculoideae). Recovered nuclear genes were analyzed using both multispecies coalescent and concatenation approaches., Results: We recovered hundreds of nuclear genes even from old and poor-quality herbarium specimens. Of particular note, we placed with strong support three incertae sedis genera (Platysace, Klotzchia, and Hermas); all three occupy isolated positions, with Platysace resolved as sister to all remaining Apiaceae. We placed nine genera (Apodicarpum, Bonannia, Grafia, Haplosciadium, Microsciadium, Physotrichia, Ptychotis, Tricholaser, Xatardia) that have never previously been included in any molecular phylogenetic study., Conclusions: We provide support for the maintenance of the four existing subfamilies of Apiaceae, while recognizing that Hermas, Klotzschia, and the Platysace clade may each need to be accommodated in additional subfamilies (pending improved sampling). The placement of the currently apioid genus Phlyctidocarpa can be accommodated by the expansion of subfamily Saniculoideae, although adequate morphological synapomorphies for this grouping are yet to be defined. This is the first phylogenetic study of the Apiaceae using high-throughput sequencing methods and represents an unprecedented evolutionary framework for the group., (© 2021 The Authors. American Journal of Botany published by Wiley Periodicals LLC on behalf of Botanical Society of America.)

Rhabosciadium aucheri is an Iranian endemic herbaceous species that grows in the west, center, and south regions of Iran. In the present study, genetic variation of 70 individuals belonging to seven natural populations of four provinces was investigated using Random Amplified Polymorphic DNA (RAPD) markers. Ten out of twenty-two RAPD primers employed in this study, generated 110 highly amplified and reproducible loci and a mean of 11.1 bands per primer and 48.13% of polymorphism was obtained. According to our results, the primer OPA10 presented the highest effective number of alleles, Shannon's index, and genetic diversity. The highest value of genetic identity (0.916) was determined between Hamadan, Nahavand and Hamadan, Alvand Mts. populations and the highest genetic distance (0.277) was observed between Hamadan, Asadabad and Kurdistan, Qorveh populations. Therefore, there is an obvious correlation between genetic diversity and geographical distribution. PCA was obtained based on RAPD molecular data and Neighbor Joining (NJ) dendrogram was provided successively. Similar results were attained employing UPGMA and Neighbor Joining dendrograms, supported by PCA ordination plot. Overall, almost moderate level of polymorphism was obtained by RAPD molecular markers at the population level that shows a reasonable amount of intraspecific variability. The reason might be due to the low level of gene flow between populations that could give rise to high genetic differentiation.

Sunflower (Helianthus annuus L.) is an important oil crop, the secondary metabolites of it include many compounds such as flavonoids and lignin. However, the research on the biosynthesis of phenolic compounds in sunflowers is still scarce. Cinnamate 4-hydroxylase (C4H) belongs to the cytochrome P450-dependent monooxygenase family and is involved in the synthesis of many phenolic compounds, but C4H in sunflowers has not yet been cloned and functionally characterized. In this study, we screened three C4H genes from the sunflower transcriptome and genomic databases, named HaC4H1, HaC4H2, and, HaC4H3, respectively. In heterologous expression experiments, we had improved a method from previous studies by the addition of restriction sites to make it easier to express multiple C4H functions and suitable for in vitro activity verification. HaC4Hs without the N-terminal membrane anchor region was fused with a redox partner of Arabidopsis thaliana cytochrome P450 enzyme (CYP450) by the method and functionally expressed in E. coli and the results showed that these three enzymes catalyzed the formation of p-coumaric acid. To further investigate whether our fusion protein approach is applicable to other C4Hs, we used this method to explore the functions of C4H from Peucedanum praeruptorum and Angelica decursiva, and they can also convert trans-cinnamic acid to p-coumaric acid. The gene expression profile showed that all three HaC4H genes showed the highest transcription levels in the roots and might be up-regulated by MeJA. In summary, these results reveal the function of HaC4Hs in sunflower and provide a simpler way to explore C4H and even other cytochrome P450 enzymes in prokaryotic expression systems., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2020 Elsevier B.V. All rights reserved.)

Key Message: Psoralen synthase and angelicin synthase responsible for the formation of psoralen and angelicin in Peucedanum praeruptorum Dunn were identified and functionally characterized, respectively. Furanocoumarins were reported to possess several activities such as anticancer, anti-inflammatory and neuroprotective, and function as phytotoxin and allelochemical in plants. Furanocoumarins are the main bioactive ingredient in P. praeruptorum which is a commonly used traditional Chinese medicine. Phenylalanine ammonia lyase (PAL), 4-coumarate: CoA ligase (4CL), p-coumaroyl CoA 2'-hyfroxylase (C2'H) were cloned previously to elucidate the biosynthetic mechanism of coumarin lactone ring. However, the genes involved in complex coumarins in P. praeruptorum have not been explored. Herein, putative psoralen synthase CYP71AJ49 and angelicin synthase CYP71AJ51 were cloned from P. praeruptorum. In vivo and in vitro yeast assays were conducted to confirm their activities. Furthermore, the results of High Performance Liquid Chromatography-Electrospray Ionization Mass Spectrometry (HPLC-ESI-MS) verified that CYP71AJ49 catalyzed the conversion of marmesin to psoralen, and CYP71AJ51 catalyzed columbianetin to angelicin. Subsequently, the expression profile showed that CYP71AJ49 and CYP71AJ51 were easily affected by environmental conditions, especially UV and temperature. The genes tissue-specific expression and compounds tissue-specific distribution pattern indicated the existence of substance transport in P. praeruptorum. Phylogenetic analysis was conducted with 27 CYP71AJs, CYP71AJ49 and CYP71AJ51 were classified in I-4 and I-2, respectively. These results provide further insight to understand the biosynthetic mechanism of complex coumarins.

Glehnia littoralis is an important medicinal halophyte-the dried root of which is used as Chinese herbal medicine. However, the use, selection and stability of reference genes are rarely verified in studies of G. littoralis, which hampers investigation of its salt tolerance and metabolism. In this study, we selected 13 candidate reference genes from the transcriptome data of G. littoralis-serine/threonine-protein phosphatase PP2A (PP2A), polyubiquitin 10 (UBQ10), actin (ACT), elongation factor 1-α (EF1-α), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), α-tubulin (α-TUB), β-tubulin (β-TUB), polypyrimidine tract-binding protein 1 (PTBP1), expressed protein 1 (EXP1), expressed protein 2 (EXP2), TIP41-like (TIP41), SAND family (SAND), and cyclophilin 2 (CYP2), and used qRT-PCR to analyse their expression levels in roots of G. littoralis treated with NaCl, polyethylene glycol (PEG), abscisic acid (ABA), and methyl jasmonate (MeJA), as well as in various organs of G. littoralis. The ΔCt, geNorm, NormFinder, and BestKeeper algorithms were used to assess the expression stability of the candidate reference genes and the results were then used to generate a comprehensive rank list with the RankAggreg R package. The most stable reference genes for normalisation were EXP1 and PP2A in response to NaCl, EXP2 and PP2A in response to ABA, CYP2 and α-TUB in response to MeJA, and ACT and EXP1 in the PEG and the organ subsets. GAPDH, β-TUB, and UBQ10 exhibited low stability and so were unsuitable for normalisation. This study is the first systematic analysis of candidate reference genes in G. littoralis and will facilitate further investigation of normalisation of gene expression in G. littoralis.

Species' geographical ranges are often restricted due to niche limitation resulting in geographical isolation and reduced population size at range margins. Under the "abundant center" paradigm, static marginal populations are thus expected to show higher genetic differentiation and lower genetic diversity than core populations. Low mate availability may also drive shifts toward higher propensity for selfing in geographically marginal populations. However, these predictions remain to be validated for contemporary range shifts occurring under current environmental change. This study is devoted to bridging this gap and assesses the spatial patterns of genetic structure and mating system across the geographical range of two coastal plant species characterized by contrasting contemporary range dynamics: the receding myrmecochorous Dune pansy (Viola tricolor subsp. curtisii) and the widespread expanding hydrochorous Rock samphire (Crithmum maritimum). Both species exhibited high propensity for selfing, with indications of inbreeding depression acting at early life stages. In Dune pansy, a biogeographical break was observed between core and marginal populations, with trailing-edge populations showing higher levels of genetic differentiation, reduced genetic diversity, and higher levels of selfing estimated through progeny arrays. In contrast, genetic structuring was weak in Rock samphire and no clear spatial trends were observed in genetic diversity nor in mating system, likely the result of efficient long-distance seed dispersal by sea-surface currents. Our study highlights that key species differences in life-history traits related to dispersal and/or mate limitation modify the expectations of genetic diversity loss and mating system shift in contemporary range-expanding populations, as compared with historical core populations.

The endangered medicinal plant Glehnia littoralis is one of the important natural source of furanocoumarin, which has been used as mucolytic, antitussive, antitumour and antibacterial. However, the genetic information of furanocoumarin biosynthesis in G. littoralis is scarce at present. The objective of this study was to mine the putative candidate genes involved in the biosynthesis pathwayof furanocoumarin and provide references for gene identification, and functional genomics of G. littoralis . We carried out the transcriptome analysis of leaves and roots in G. littoralis , which provided a dataset for gene mining. Psoralen, imperatorin and isoimperatorin were detected in G. littoralis by high performance liquid chromatography analysis. Candidate key genes were mined based on the annotations and local BLAST with homologous sequences using BioEdit software. The relative expression of genes was analysed using quantitative real-time polymerase chain reaction. Further, the CYP450 genes were mined using phylogenetic analyses using MEGA 6.0 software. Atotal of 156,949 unigenes were generated, of which 9021 were differentially-expressed between leaves and roots. A total of 82 unigenes encoding eight enzymes in furanocoumarin biosynthetic pathway were first obtained. Seven genes that encoded key enzymes in the downstream furanocoumarin biosynthetic pathway and expressed more in roots than leaves were screened. Twenty-six candidate CYP450 unigenes expressed abundantly in roots and were chiefly concentrated in CYP71, CYP85 and CYP72 clans. Finally, we filtered 102 differentially expressed transcription factors (TFs) unigenes. The transcriptome of G. littoralis was characterized which would help to elucidate the furanocoumarin biosynthetic pathway in G. littoralis and provide an invaluable resource for further study of furanocoumarin.

The distribution of genetic variation in species is governed by factors that act differently across spatial scales. To tease apart the contribution of different processes, especially at intermediate spatial scales, it is useful to study simple ecosystems such as those on sub-Antarctic oceanic islands. In this study, we characterize spatial genetic patterns of two keystone plant species, Azorella selago on sub-Antarctic Marion Island and Azorella macquariensis on sub-Antarctic Macquarie Island. Although both islands experience a similar climate and have a similar vegetation structure, they differ significantly in topography and geological history. We genotyped six microsatellites for 1,149 individuals from 123 sites across Marion Island and 372 individuals from 42 sites across Macquarie Island. We tested for spatial patterns in genetic diversity, including correlation with elevation and vegetation type, and clines in different directional bearings. We also examined genetic differentiation within islands, isolation-by-distance with and without accounting for direction, and signals of demographic change. Marion Island was found to have a distinct northwest-southeast divide, with lower genetic diversity and more sites with a signal of population expansion in the northwest. We attribute this to asymmetric seed dispersal by the dominant northwesterly winds, and to population persistence in a southwestern refugium during the Last Glacial Maximum. No apparent spatial pattern, but greater genetic diversity and differentiation between sites, was found on Macquarie Island, which may be due to the narrow length of the island in the direction of the dominant winds and longer population persistence permitted by the lack of extensive glaciation on the island. Together, our results clearly illustrate the implications of island shape and geography, and the importance of direction-dependent drivers, in shaping spatial genetic structure., (© 2019 John Wiley & Sons Ltd.)

Three Apiaceae species Ledebouriella seseloides , Peucedanum japonicum , and Glehnia littoralis are used as Asian herbal medicines, with the confusingly similar common name "Bang-poong". We characterized the complete chloroplast (cp) genomes and 45S nuclear ribosomal DNA (45S nrDNA) sequences of two accessions for each species. The complete cp genomes of G. littoralis , L. seseloides , and P. japonicum were 147,467, 147,830, and 164,633 bp, respectively. Compared to the other species, the P. japonicum cp genome had a huge inverted repeat expansion and a segmental inversion. The 45S nrDNA cistron sequences of the three species were almost identical in size and structure. Despite the structural variation in the P. japonicum cp genome, phylogenetic analysis revealed that G. littoralis diverged 5-6 million years ago (Mya), while P. japonicum diverged from L. seseloides only 2-3 Mya. Abundant copy number variations including tandem repeats, insertion/deletions, and single nucleotide polymorphisms, were found at the interspecies level. Intraspecies-level polymorphism was also found for L. seseloides and G. littoralis . We developed nine PCR barcode markers to authenticate all three species. This study characterizes the genomic differences between L. seseloides , P. japonicum , and G. littoralis ; provides a method of species identification; and sheds light on the evolutionary history of these three species.

Species identification and discrimination is the basis of biodiversity research. In general, it is considered that numerous nucleotide variations (e.g., whole chloroplast genomes) can identify species with higher resolution than a few loci, e.g., partial chloroplast or nuclear gene fragments. In this study, we tested this hypothesis by sampling population genetics samples of the endangered herb genus Notopterygium. We sequenced the complete plastomes, five nuclear gene regions, three chloroplast DNA fragments, and a nuclear internal transcribed spacer (nrITS) region for 18 populations sampled throughout most of the geographic ranges of all six Notopterygium species. Species identification analysis showed that four DNA barcodes (matK, rbcL, trnS-trnG, and nrITS) and/or combinations of these markers achieved Notopterygium species discrimination at higher resolution than the general plastomes and nuclear gene sequences. In particular, nrITS had the highest discriminatory power among all of the individual markers. Molecular data sets and morphological evidence indicated that all six Notopterygium species could be reclassified unambiguously to four putative species clades. N. oviforme and N. franchetii had the closest relationship. Molecular dating showed that the origin and divergence of Notopterygium species was significantly associated with geological and climatic fluctuations during the middle of the Pliocene. In conclusion, our results suggest that a few nucleotide variations can achieve species discrimination with higher resolution than numerous plastomes and general nuclear gene fragments when discerning related Notopterygium species., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Methoxylated coumarins represent a large proportion of officinal value coumarins while only one enzyme specific to bergaptol O -methylation (BMT) has been identified to date. The multiple types of methoxylated coumarins indicate that at least one unknown enzyme participates in the O -methylation of other hydroxylated coumarins and remains to be identified. Combined transcriptome and metabonomics analysis revealed that an enzyme similar to caffeic acid O -methyltransferase (COMT-S, S is short for similar) was involved in catalyzing all the hydroxylated coumarins in Peucedanum praeruptorum . However, the precise molecular mechanism of its substrate heterozygosis remains unsolved. Pursuing this question, we determined the crystal structure of COMT-S to clarify its substrate preference. The result revealed that Asn132, Asp271, and Asn325 govern the substrate heterozygosis of COMT-S. A single mutation, such as N132A, determines the catalytic selectivity of hydroxyl groups in esculetin and also causes production differences in bergapten. Evolution-based analysis indicated that BMT was only recently derived as a paralogue of caffeic acid O -methyltransferase (COMT) via gene duplication, occurring before the Apiaceae family divergence between 37 and 100 mya. The present study identified the previously unknown O -methylation steps in coumarin biosynthesis. The crystallographic and mutational studies provided a deeper understanding of the substrate preference, which can be used for producing specific O -methylation coumarins. Moreover, the evolutionary relationship between BMT and COMT-S was clarified to facilitate understanding of evolutionary events in the Apiaceae family.

The species of the aromatic plant family Apiaceae are mainly used as spices and foods, but the family also includes medicinal and some poisonous plant species. Due to the similar chemical compounds or aroma and morphology, the poisonous species are often mistaken for the edible aromatic species. It is therefore imperative to correctly identify the species present at the initial raw stage samples to ensure product safety and efficacy. At the molecular level, plant species can be identified using DNA loci either from nuclear or plastid genome with easily available universal oligonucleotides, a technique called DNA barcoding. However, this is possible when single-species plant material is present but may not work on a mixture of plants species. Another disadvantage is that using universal oligonucleotides is of limited help, especially if the adulterating material is present in low quantities. On the other hand, if using the species-specific oligonucleotides, only single specific adulterating plant material could be detected and, consequently, the unexpected adulterants may go undetected. Therefore, in the current work, four degenerated oligonucleotides from ITS1 and ITS2 regions of the nuclear genome were designed that can bind to a variety of Apiaceae genera only and not to other genera belonging to different plant families. These family-specific oligonucleotides were able to amplify a diagnostic PCR product from 16 Apiaceae species that, upon sequencing, revealed the identity of the plant it was derived from. The size of these products is around 140 bp for ITS1 and approximately 80 bp for the ITS2 region. The size range of the amplified products falls in the category of a desired mini-barcode size to be used for damaged/fragmented DNA and next generation sequencing., Competing Interests: The authors declare no conflict of interest., (Georg Thieme Verlag KG Stuttgart · New York.)

Background: Bunium persicum commonly called as Kala zeera, a very high value herbaceous spice used for medicinal purposes is often adulterated with Cuminum cyminum or Safed zeera, a closely related species. Lack of distinctive morphological features makes the identification of genuine kala zeera from its adulterant difficult, the problem is even exaggerated in case of powdered material., Methodology: Genomic DNA was extracted from all the plant materials by using CTAB-SDS method (Möller et al., 1992) with slight modifications. On the basis of reproducibility and high amplification ability, four universal barcoding loci viz. ITS2, rbcL-a, mat K and psbA-trnH and a specific locus Cum were used in the present study. The amplified PCR products were sequenced bidirectionally and assembled to obtain contigs. The sequences thus obtained were aligned using MUSCLE algorithm (Edgar, 2004) and information pertaining to conserved/ variable/ parsimony informative sites, number of transitions, transversions and Indels was obtained after analyzing the sequences., Results and Conclusion: Among the tested barcoding loci, psbA-trnH has proven to be best barcode in authentication of kala zeera as its amplification and sequencing success was high and it showed the presence of polymorphic sites to detect interspecific variation. This barcode could differentiate between safed zeera and kala zeera in a single reaction, simultaneously., (Copyright © 2018 Elsevier GmbH. All rights reserved.)

In the southern Andes mountains (27-39◦S) Azorella madreporica and Laretia acaulis , two Apiaceae cushion plant species commonly known as yaretas , conform a well-established altitudinal vegetation belt along the lower Andean zone. These species have been considered as fundamental components of several ecological dynamics within their communities; however, high mountain ecosystems are increasingly threatened worldwide by natural and anthropogenic pressures and the southern Andes are not the exception. Recognizing that genetic information is crucial for the success of any conservation or restoration initiative inwild populations, we developed and cross-amplified 28 specifically designed microsatellite markers (14 in A. madreporica and 14 in L. acaulis ), and also tested the cross amplification of 25 markers from the related species Azorella selago . In a region which is particularly vulnerable to global change trends, this new polymorphic microsatellite loci will be useful in the study of the genetic diversity of these high-mountain cushion plants, which are pivotal in the structuring of their native ecosystems.