[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.

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.

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.

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.

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.

[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.

Combining genome size and pollen morphology data to study species relationships in the genus Daucus (Apiaceae).

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.

The Chloroplast Genome of Wild Saposhnikovia divaricata: Genomic Features, Comparative Analysis, and Phylogenetic Relationships.

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.

Development and characterization of microsatellite markers for the French endemic Angelica heterocarpa (Apiaceae) and congeneric sympatric species.

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.

Two CYP71AJ enzymes function as psoralen synthase and angelicin synthase in the biosynthesis of furanocoumarins in Peucedanum praeruptorum Dunn.

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.

Selection and Validation of Reference Genes for Normalisation of Gene Expression in Glehnia littoralis.

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.

Mining genes associated with furanocoumarin biosynthesis in an endangered medicinal plant, Glehnia littoralis.

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.

Dynamic Chloroplast Genome Rearrangement and DNA Barcoding for Three Apiaceae Species Known as the Medicinal Herb "Bang-Poong".

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.

Identification of Species in the Aromatic Spice Family Apiaceae Using DNA Mini-barcodes.

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.

DNA barcoding for specific and sensitive detection of Cuminum cyminum adulteration in Bunium persicum.

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.

Comparative transcriptome analysis to identify putative genes involved in thymol biosynthesis pathway in medicinal plant Trachyspermum ammi L.

Thymol, as a dietary monoterpene, is a phenol derivative of cymene, which is the major component of the essential oil of Trachyspermum ammi (L.). It shows multiple biological activities: antifungal, antibacterial, antivirus and anti-inflammatory. T. ammi, commonly known as ajowan, belongs to Apiaceae and is an important medicinal seed spice. To identify the putative genes involved in thymol and other monoterpene biosynthesis, we provided transcriptomes of four inflorescence tissues of two ajowan ecotypes, containing different thymol yield. This study has detected the genes encoding enzymes for the go-between stages of the terpenoid biosynthesis pathways. A large number of unigenes, differentially expressed between four inflorescence tissues of two ajowan ecotypes, was revealed by a transcriptome analysis. Furthermore, differentially expressed unigenes encoding dehydrogenases, transcription factors, and cytochrome P450s, which might be associated with terpenoid diversity in T. ammi, were identified. The sequencing data obtained in this study formed a valuable repository of genetic information for an understanding of the formation of the main constituents of ajowan essential oil and functional analysis of thymol-specific genes. Comparative transcriptome analysis led to the development of new resources for a functional breeding of ajowan.

Transcriptomic analysis identifies differentially expressed genes (DEGs) associated with bolting and flowering in Saposhnikovia divaricata.

Saposhnikovia divaricata is a valuable Chinese medicinal herb; the transformation from vegetative growth to reproductive growth may lead to the decrease of its pharmacological activities. Therefore, the study of bolting and flowering for Saposhnikovia divaricata is warranted. The present study aimed to reveal differentially expressed genes (DEGs) and regularity of expression during the bolting and flowering process, and the results of this study might provide a theoretical foundation for the suppression of early bolting for future research and practical application. Three sample groups, early flowering, flower bud differentiation, and late flowering (groups A, B, and C, respectively) were selected. Transcriptomic analysis identified 67, 010 annotated unigenes, among which 50, 165 were differentially expressed including 16, 108 in A vs B, and 17, 459 in B vs C, respectively. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway functional classification analysis were performed on these differentially expressed genes, and five important pathways were significantly impacted (P ≤ 0.01): plant circadian rhythm, other glycan degradation, oxidative phosphorylation, plant hormone signal transduction, and starch and sucrose metabolism. Plant hormone signal transduction might play an important role in the bolting and flowering process. The differentially expressed indole-3-acetic acid (IAA) gene showed significant down-regulation during bolting and flowering, while the transport inhibitor response 1 (TIR1) gene showed no significant change during the bolting process. The expression of flowering related genes FLC, LYF, and AP1 also showed a greater difference at different development stages. In conclusion, we speculate that the decrease in auxin concentration is not caused by the degrading effect of TIR1 but by an alternative mechanism.

Comparative Transcriptome Analysis Reveals Adaptive Evolution of Notopterygium incisum and Notopterygium franchetii, Two High-Alpine Herbal Species Endemic to China.

The extreme conditions (e.g., cold, low oxygen, and strong ultraviolet radiation) of the high mountains provide an ideal natural laboratory for studies on speciation and the adaptive evolution of organisms. Up to now, few genome/transcriptome-based studies have been carried out on how plants adapt to conditions at extremely high altitudes. Notopterygium incisum and Notopterygium franchetii (Notopterygium, Apiaceae) are two endangered high-alpine herbal plants endemic to China. To explore the molecular genetic mechanisms of adaptation to high altitudes, we performed high-throughput RNA sequencing (RNA-seq) to characterize the transcriptomes of the two species. In total, more than 130 million sequence reads, 81,446 and 63,153 unigenes with total lengths of 86,924,837 and 62,615,693 bp, were generated for the two herbal species, respectively. OrthoMCL analysis identified 6375 single-copy orthologous genes between N. incisum and N. franchetii. In total, 381 positively-selected candidate genes were identified for both plants by using estimations of the non-synonymous to synonymous substitution rate. At least 18 of these genes potentially participate in RNA splicing, DNA repair, glutathione metabolism and the plant-pathogen interaction pathway, which were further enriched in various functional gene categories possibly responsible for environment adaptation in high mountains. Meanwhile, we detected various transcription factors that regulated the material and energy metabolism in N. incisum and N. franchetii, which probably play vital roles in the tolerance to stress in surroundings. In addition, 60 primer pairs based on orthologous microsatellite-containing sequences between the both Notopterygium species were determined. Finally, 17 polymorphic microsatellite markers (SSR) were successfully characterized for the two endangered species. Based on these candidate orthologous and SSR markers, we detected that the adaptive evolution and species divergence of N. incisum and N. franchetii were significantly associated with the extremely heterogeneous environments and climatic oscillations in high-altitude areas. This work provides important insights into the molecular mechanisms of adaptation to high-altitudes in alpine herbal plants.