A single origin and high genetic diversity of cultivated medicinal herb Glehnia littoralis subsp. littoralis (Apiaceae) deciphered by SSR marker and phenotypic analysis.

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.

The gradual establishment of complex coumarin biosynthetic pathway in Apiaceae.

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.

Plastome structure and phylogenetic relationships of genus Hydrocotyle (apiales): provide insights into the plastome evolution of Hydrocotyle.

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.

Morphology and ITS sequences provide insights into the phylogeny of Tongoloa (Apiaceae) from China.

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.

The telomere-to-telomere (T2T) genome of Peucedanum praeruptorum Dunn provides insights into the genome evolution and coumarin biosynthesis.

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.

Combined metabolome and transcriptome reveal HmF6'H1 regulating simple coumarin accumulation against powdery mildew infection in Heracleum moellendorffii Hance.

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.

Molecular authentication of commercial "Qian-hu" through the integration of nrDNA internal transcribed spacer 2 and nucleotide signature.

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.

[Development and application of SSR markers of Saposhnikovia divaricata based on transcriptome].

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.

Identification of Angelica acutiloba, A. sinensis, and other Chinese medicinal Apiaceae plants by DNA barcoding.

Crude drug Angelicae acutilobae radix is one of the most important crude drugs in Japanese traditional medicine and is used mainly for the treatment of gynecological disorders. In the listing in the Japanese Pharmacopoeia XVIII, Angelicae acutilobae radix is defined as the root of Angelica acutiloba (Apiaceae), which has long been produced on an industrial scale in Japan. With the aging of farmers and depopulation of production areas, the domestic supply has recently declined and the majority of the supply is now imported from China. Due to having only slightly different morphological and chemical characteristics for the Apiaceae roots used to produce dried roots for Chinese medicines, the plant species originating the crude drug Apiaceae roots may be incorrectly identified. In particular, Angelicae sinensis radix, which is widely used in China, and Angelicae acutilobae radix are difficult to accurately identify by morphology and chemical profiles. Thus, in order to differentiate among Angelicae acutilobae radix and other radixes originated from Chinese medicinal Apiaceae plants, we established DNA markers. Using DNA sequences for the chloroplast psbA-trnH intergenic spacer and nuclear internal transcribed spacer regions, Angelicae acutilobae radix and other Chinese Apiaceae roots, including Angelicae sinensis radix, can be definitively identified.

The high-quality genome of Cryptotaenia japonica and comparative genomics analysis reveals anthocyanin biosynthesis in Apiaceae.

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, Ks, 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_2.2062 (bHLH) and Cjaponica_1.3743 (MYB). Our findings provided a high-quality reference genome of C. japonica and offered new insights into Apiaceae evolution and biology.

A chromosome-scale genome of Peucedanum praeruptorum provide insights into Apioideae evolution and medicinal ingredient biosynthesis.

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.

Comparative Analysis of the Complete Mitochondrial Genomes of Apium graveolens and Apium leptophyllum Provide Insights into Evolution and Phylogeny Relationships.

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.

Genome sequencing-based transcriptomic analysis reveals novel genes in Peucedanum praeruptorum.

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.

Plastome evolution in the genus Sium (Apiaceae, Oenantheae) inferred from phylogenomic and comparative analyses.

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.

Chloroplast genome characteristics and phylogeny of the sinodielsia clade (apiaceae: apioideae).

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.

Plastid Phylogenomic Analyses Reveal a Cryptic Species of Ligusticopsis (Apiaceae, Angiosperms).

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.

Understanding the response in Pugionium cornutum (L.) Gaertn. seedling leaves under drought stress using transcriptome and proteome integrated analysis.

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.

Molecular phylogenetic study of flavonoids in medicinal plants: a case study family Apiaceae.

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.

Comparative transcriptome analysis of Saposhnikovia divaricata to reveal drought and rehydration adaption strategies.

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.

The hybridization origin of the Chinese endemic herb genus Notopterygium (Apiaceae): Evidence from population genomics and ecological niche analysis.

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.