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.

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.

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.

Comparative analysis of root anatomical structure, chemical components and differentially expressed genes between early bolting and unbolting in Peucedanum praeruptorum Dunn.

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.

Identification of Apiaceae using ITS, ITS2 and psba-trnH barcodes.

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.

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

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

[Transcriptome profiling of Saposhnikovia divaricata growing for different years and mining of key genes in active ingredient biosynthesis].

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.