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.

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

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.

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.

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.

Plastid Phylogenomics Provide Evidence to Accept Two New Members of Ligusticopsis (Apiaceae, Angiosperms)

Peucedanum nanum and P. violaceum are recognized as members of the genus Peucedanum because of their dorsally compressed mericarps with slightly prominent dorsal ribs and narrowly winged lateral ribs. However, these species are not similar to other Peucedanum taxa but resemble Ligusticopsis in overall morphology. To check the taxonomic positions of P. nanum and P. violaceum, we sequenced their complete plastid genome (plastome) sequences and, together with eleven previously published Ligusticopsis plastomes, performed comprehensively comparative analyses. The thirteen plastomes were highly conserved and similar in structure, size, GC content, gene content and order, IR borders, and the patterns of codon bias, RNA editing, and simple sequence repeats (SSRs). Nevertheless, twelve mutation hotspots (matK, ndhC, rps15, rps8, ycf2, ccsA-ndhD, petN-psbM, psbA-trnK, rps2-rpoC2, rps4-trnT, trnH-psbA, and ycf2-trnL) were selected. Moreover, both the phylogenetic analyses based on plastomes and on nuclear ribosomal DNA internal transcribed spacer (ITS) sequences robustly supported that P. nanum and P. violaceum nested in Ligusticopsis, and this was further confirmed by the morphological evidence. Hence, transferring P. nanum and P. violaceum into Ligusticopsis genus is reasonable and convincing, and two new combinations are presented.

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.

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.

Polyphenolics, antioxidant characterization and DNA barcoding of Kala zeera [Bunium persicum (Boiss.) Fedtsch] through multiple barcode analysis to unravel best barcode combination

BACKGROUND

Kala zeera [Bunium persicum (Boiss.) Fedtsch] is one of the important spice crops of North Western Himalayas with lot of medicinal and culinary values. In spite of having great importance, this crop is under the threat of extinction due to loss of habitat and lack of awareness. The limited availability of the seeds has ultimately increased the economic value of this spice. The upmarket of Kala zeera leads to its adulteration with other black seeds and cumin seeds. The present investigation was undertaken to evaluate polyphenolics and antioxidant properties of Kala zeera genotypes collected from North Western Himalayas and to develop DNA barcodes that can ensure their purity and can also guide in conservation of selected Kala zeera germplasm lines.

METHODS AND RESULTS

Various locations of North Western Himalayas were explored for collecting 31 diverse germplasm lines of Kala zeera. The collected germplasm was maintained at our experimental stations during 2019-2020 and 2020-2021. These genotypes were evaluated for different seed traits and the methanolic extract from Kala zeera seeds was examined for total phenolic content, total flavonoid content, antioxidant activities by DPPH and FRAP. The results revealed significant variation in seed traits, polyphenolic content and antioxidant properties. 100 seed weight ranged from 0.05 to 0.35 g, TPC ranged from 7.5 to 22.56 mg/g, TFC ranged from 0.58 to 4.15 mg/g, antioxidant properties DPPH ranged from 168 to 624.4 μg/ml and FRAP ranged from 0.72 to 6.91 mg/g. Further, three different barcodes (ITS, rbcL and psbA-trnH) were used to reveal the authenticity of selected Kala zeera. MEGA 5 software was used for clustering and the barcodes did clustering based on geographical distribution of Kala zeera germplasm.

CONCLUSION

Based on molecular barcoding, best barcode combination was identified that may discriminate the Kala zeera germplasm vis-a-vis can authenticate their purity. Moreover, the identified DNA barcodes will have significant role in studying the evolutionary biology of Bunium species and will be important for designing a strategy to conserve the selected Kala zeera germplasm lines. The identified genotypes with high phenolic content and antioxidant activity can further be utilized in Kala zeera breeding programmes.

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

BACKGROUND

The Peucedanum genus is the backbone member of Apiaceae, with many economically and medically important plants. Although the previous studies on Peucedanum provide us with a good research basis, there are still unclear phylogenetic relationships and many taxonomic problems in Peucedanum, and a robust phylogenetic framework of this genus still has not been obtained, which severely hampers the improvement and revision of taxonomic system for this genus. The plastid genomes possessing more variable characters have potential for reconstructing a robust phylogeny in plants.

RESULTS

In the current study, we newly sequenced and assembled seven Peucedanum plastid genomes. Together with five previously published plastid genomes of Peucedanum, we performed a comprehensively comparative analyses for this genus. Twelve Peucedanum plastomes were similar in terms of genome structure, codon bias, RNA editing sites, and SSRs, but varied in genome size, gene content and arrangement, and border of SC/IR. Fifteen mutation hotspot regions were identified among plastid genomes that can serve as candidate DNA barcodes for species identification in Peucedanum. Our phylogenetic analyses based on plastid genomes generated a phylogeny with high supports and resolutions for Peucedanum that robustly supported the non-monophyly of genus Peucedanum.

CONCLUSION

The plastid genomes of Peucedanum showed both conservation and diversity. The plastid genome data were efficient and powerful for improving the supports and resolutions of phylogeny for the complex Peucedanum genus. In summary, our study provides new sights into the plastid genome evolution, taxonomy, and phylogeny for Peucedanum species.

Impact of Genomic and Transcriptomic Resources on Apiaceae Crop Breeding Strategies

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

Derived woodiness and annual habit evolved in African umbellifers as alternative solutions for coping with drought

BACKGROUND

One of the major trends in angiosperm evolution was the shift from woody to herbaceous habit. However, reversals known as derived woodiness have also been reported in numerous, distantly related clades. Among theories evoked to explain the factors promoting the evolution of derived woodiness are moderate climate theory and cavitation theory. The first assumes that woody habit evolves in response to mild climate allowing for prolonged life span, which in turn leads to bigger and woodier bodies. The second sees woodiness as a result of natural selection for higher cavitation resistance in seasonally dry environments. Here, we compare climatic niches of woody and herbaceous, mostly southern African, umbellifers from the Lefebvrea clade to assess whether woody taxa in fact occur in markedly drier habitats. We also calibrate their phylogeny to estimate when derived woodiness evolved. Finally, we describe the wood anatomy of selected woody and herbaceous taxa to see if life forms are linked to any particular wood traits.

RESULTS

The evolution of derived woodiness in chamaephytes and phanerophytes as well as the shifts to short-lived annual therophytes in the Lefebvrea clade took place at roughly the same time: in the Late Miocene during a trend of global climate aridification. Climatic niches of woody and herbaceous genera from the Cape Floristic Region overlap. There are only two genera with distinctly different climatic preferences: they are herbaceous and occur outside of the Cape Floristic Region. Therefore, studied herbs have an overall climatic niche wider than their woody cousins. Woody and herbaceous species do not differ in qualitative wood anatomy, which is more affected by stem architecture and, probably, reproductive strategy than by habit.

CONCLUSIONS

Palaeodrought was likely a stimulus for the evolution of derived woodiness in the Lefebvrea clade, supporting the cavitation theory. The concurrent evolution of short-lived annuals withering before summer exemplifies an alternative solution to the same problem of drought-induced cavitation. Changes of the life form were most likely neither spurred nor precluded by any qualitative wood traits, which in turn are more affected by internode length and probably also reproductive strategy.

A higher-level nuclear phylogenomic study of the carrot family (Apiaceae)

PREMISE

The carrot family (Apiaceae) comprises 466 genera, which include many well-known crops (e.g., aniseed, caraway, carrots, celery, coriander, cumin, dill, fennel, parsley, and parsnips). Higher-level phylogenetic relationships among subfamilies, tribes, and other major clades of Apiaceae are not fully resolved. This study aims to address this important knowledge gap.

METHODS

Target sequence capture with the universal Angiosperms353 probe set was used to examine phylogenetic relationships in 234 genera of Apiaceae, representing all four currently recognized subfamilies (Apioideae, Azorelloideae, Mackinlayoideae, and Saniculoideae). Recovered nuclear genes were analyzed using both multispecies coalescent and concatenation approaches.

RESULTS

We recovered hundreds of nuclear genes even from old and poor-quality herbarium specimens. Of particular note, we placed with strong support three incertae sedis genera (Platysace, Klotzchia, and Hermas); all three occupy isolated positions, with Platysace resolved as sister to all remaining Apiaceae. We placed nine genera (Apodicarpum, Bonannia, Grafia, Haplosciadium, Microsciadium, Physotrichia, Ptychotis, Tricholaser, Xatardia) that have never previously been included in any molecular phylogenetic study.

CONCLUSIONS

We provide support for the maintenance of the four existing subfamilies of Apiaceae, while recognizing that Hermas, Klotzschia, and the Platysace clade may each need to be accommodated in additional subfamilies (pending improved sampling). The placement of the currently apioid genus Phlyctidocarpa can be accommodated by the expansion of subfamily Saniculoideae, although adequate morphological synapomorphies for this grouping are yet to be defined. This is the first phylogenetic study of the Apiaceae using high-throughput sequencing methods and represents an unprecedented evolutionary framework for the group.

Backbone phylogeny and evolution of Apioideae (Apiaceae): New insights from phylogenomic analyses of plastome data

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

Functional characterization of cinnamate 4-hydroxylase from Helianthus annuus Linn using a fusion protein method

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

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

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.

Effects of contemporary shifts of range margins on patterns of genetic structure and mating system in two coastal plant species

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

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.

The Molecular and Structural Basis of O-methylation Reaction in Coumarin Biosynthesis in Peucedanum praeruptorum Dunn

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

Species delimitation and interspecific relationships of the endangered herb genus Notopterygium inferred from multilocus variations

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

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.

Isolation and cross-amplification of the first set of polymorphic microsatellite markers of two high-Andean cushion plants

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

Phylogeography of Libanotis buchtormensis (Umbelliferae) in Disjunct Populations along the Deserts in Northwest China

In Northwest China, aridification and desert expansion play significant roles in promoting desert plant diversification and speciation. However, to date, little is known about the effects of the desert barrier on the population structure of montane, non-desert species in the area. In this study, we sequenced chloroplast DNA regions (trnL–trnF and trnS–trnG) and a nuclear gene (rpb2) to investigate the population differentiation and phylogeographical history of Libanotis buchtormensis, a perennial montane species possessing a disjunct distribution at the periphery of the central desert. In total, 23 chloroplast haplotypes and 24 nuclear haplotypes were recovered from the 21 natural populations and six hebarium specimens. Phylogenetic analysis based on the combined plastid and nuclear dataset revealed two distinct lineages of L. buchtormensis, which inhabit the disjunct areas on both sides of the desert zone. The molecular dating analysis indicated that the divergence between the southeastern and the northwestern populations occurred in the middle Pleistocene, concomitantly with the desert expansion. The geographical vicariance likely contributed to the present disjunct distribution of L. buchtormensis across the deserts in Northwest China. Populations in the southeastern region may have migrated from the northwestern region, and seem to be a peripheral distribution of L. buchtormensis.

Selection of Reference Genes for Gene Expression Normalization in Peucedanum praeruptorum Dunn under Abiotic Stresses, Hormone Treatments and Different Tissues

Peucedanum praeruptorum Dunn is one of the main traditional Chinese medicines producing coumarins and plenty of literatures are focused on the biosynthesis of coumarins. Quantitative real-time reverse transcription PCR (qRT-PCR) is a widely used method in studying the biosynthesis pathway and the selection of reference genes plays a crucial role in accurate normalization. To facilitate biosynthesis study of coumarins, twelve candidate reference genes were selected from the transcriptome database of P. praeruptorum according to previous studies. Then, BestKeeper, geNoFrm and NormFinder were used for selecting stably expressed reference genes in different tissues and under various stress treatments. The results indicated that, among the twelve candidate reference genes, the SAND family protein (SAND), actin 2 (ACT2), ubiquitin-conjugating enzyme 9 (UBC9), protein phosphatase 2A gene (PP2A) and polypyrimidine tract-binding protein (PTBP1) were the most stable reference genes under different experimental treatments, while glyceraldehyde 3-phosphate dehydrogenase (GAPDH) and tubulin beta-6 (TUB6) were the least stable genes. In addition, the suitability of SAND, TIP41-like protein (TIP41), UBC9, ACT2, TUB6 and their combination as reference genes were confirmed by normalizing the expression of 1-aminocyclopropane-1-carboxylate oxidase (ACO) in different treatments. This work is the first survey of the stability of reference genes in P. praeruptorum and provides guidelines to obtain more accurate qRT-PCR results in P. praeruptorum and other plant species.

[ITS sequences variation and phylogenetic analysis on 31 geographical populations of Notopterygium incisum]

In this study, 31 Notopterygium incisum populations were analyzed using ITS sequences to investigate the genetic structure. The results showed that: the ITS region ranged in size from 634 to 635 bp and base composition was with high G + C content of 57.8%. Thirty-one polymorphic sites were detected from 402 sequences of 31 populations of N. incisum, and the proportion of polymorphic sites was 4.88%, in which parsimony informative sites were up to 12. And 31 haplotypes were identified based on these polymorphic sites. Molecular variance analysis (AMOVA) indicated that high genetic differentiation (57%) existed among population, and gene flow was low (N(m) = 0.38) among populations. Phylogenetic relationships of 31 haplotypes were analyzed using NJ method with N. forbesiias an out-group. Phylogenetic analysis showed that 31 haplotypes from different populations mixed together and did not form distinct geographically separated clades.

Evolution of substrate recognition sites (SRSs) in cytochromes P450 from Apiaceae exemplified by the CYP71AJ subfamily

BACKGROUND

Large proliferations of cytochrome P450 encoding genes resulting from gene duplications can be termed as 'blooms', providing genetic material for the genesis and evolution of biosynthetic pathways. Furanocoumarins are allelochemicals produced by many of the species in Apiaceaous plants belonging to the Apioideae subfamily of Apiaceae and have been described as being involved in the defence reaction against phytophageous insects.

RESULTS

A bloom in the cytochromes P450 CYP71AJ subfamily has been identified, showing at least 2 clades and 6 subclades within the CYP71AJ subfamily. Two of the subclades were functionally assigned to the biosynthesis of furanocoumarins. Six substrate recognition sites (SRS1-6) important for the enzymatic conversion were investigated in the described cytochromes P450 and display significant variability within the CYP71AJ subfamily. Homology models underline a significant modification of the accession to the iron atom, which might explain the difference of the substrate specificity between the cytochromes P450 restricted to furanocoumarins as substrates and the orphan CYP71AJ.

CONCLUSION

Two subclades functionally assigned to the biosynthesis of furanocoumarins and four other subclades were identified and shown to be part of two distinct clades within the CYP71AJ subfamily. The subclades show significant variability within their substrate recognition sites between the clades, suggesting different biochemical functions and providing insights into the evolution of cytochrome P450 'blooms' in response to environmental pressures.

Morphometrics of Daucus (Apiaceae): a counterpart to a phylogenomic study

PREMISE OF STUDY

Molecular phylogenetics of genome-scale data sets (phylogenomics) often produces phylogenetic trees with unprecedented resolution. A companion phylogenomics analysis of Daucus using 94 conserved nuclear orthologs supported many of the traditional species but showed unexpected results that require morphological analyses to help interpret them in a practical taxonomic context.

METHODS

We evaluated character state distributions, stepwise discriminant analyses, canonical variate analyses, and hierarchical cluster analyses from 40 morphological characters from 81 accessions of 14 taxa of Daucus and eight species in related genera in an experimental plot.

KEY RESULTS

Most characters showed tremendous variation with character state overlap across many taxa. Multivariate analyses separated the outgroup taxa easily from the Daucus ingroup. Concordant with molecular analyses, most species form phenetic groups, except the same taxa that are problematical in the molecular results: (1) the subspecies of D. carota, (2) D. sahariensis and D. syrticus, and (3) D. broteri and D. guttatus.

CONCLUSIONS

Phenetic analyses, in combination with molecular data, support many Daucus species, but mostly by overlapping ranges of size and meristic variation. The subspecies of D. carota are poorly separated morphologically, are paraphyletic, and all could be recognized at the subspecies rank under D. carota. Daucus sahariensis and D. syrticus are so similar morphologically that they could be placed in synonymy. Combined molecular and morphological data support three species in accessions previously identified as D. broteri and D. guttatus. Molecular and morphological results support the new combination Daucus carota subsp. capillifolius.

[Identification of peucedani radix, peucedani decursivi radix and its adulterants using ITS2 sequence]

In order to identify Peucedani Radix, Peucedani Decursivi Radix and their adulterants, the internal transcribed spacer 2 (ITS2) regions of Peucedani Radix, Peucedani Decursivi Radix and their adulterants were amplified and bidirectionally sequenced based on the Principles for Molecular Identification of Traditional Chinese Materia Medica Using DNA Barcoding, which has been promulgated by Chinese Pharmacopoeia Commission. Sequences were analyzed and assembled by Codon Code Aligner V3. 7.1. The relevant data were analyzed by MEGA 5. 0. Species identification analyses were performed by using the nearest distance methods and neighbor-joining (NJ) methods. The result showed that the ITS2 sequence lengths of Peucedani Radix were 229-230 bp and the average intra-specific genetic distances were 0.005. The ITS2 sequence lengths of Peucedani Decursivi Radix were 227 bp and the sequences contained no variation site. The average inter-specific K2P genetic distance of Peucedani Radix, Peucedani Decursivi Radix and their adulterants species were 0.044 and 0.065 respectively. The minimum inter-specific divergence is larger than the maximum intra-specific divergence of Peucedani Decursivi Radix. The nearest distance methods and NJ trees results indicated that Peucedani Radix, Peucedani Decursivi Radix and their adulterants species could be identification clearly. The ITS2 regions can stably and accurately distinguish Peucedani Radix, Peucedani Decursivi Radix and their adulterants.

[Identification of Peucedani Radix and its adulterants by DNA barcoding technique]

OBJECTIVE

To identify Peucedani Radix and its adulterants using DNA barcoding technique.

METHODS

Total genomic DNA was isolated from Peucedani Radix and its adulterants. Nuclear DNA ITS2 sequences were amplified and purified PCR products were sequenced. Sequence assembly and consensus sequence generation were performed using the CodonCode Aligner V3.0. The Kimura 2-Parameter(K2P) distances were calculated using software MEGA 4. 0. Identification analyses were performed using BLAST1, Nearest Distance and Neighbor-Joining (NJ) methods, and the secondary structure of the ITS2 sequence differences between species were analyzed.

RESULTS

Different samples of Peucedani Radix were gathered together and distinguished from its adulterants by NJ tree. The ITS2 secondary structure showed that Peucedani Radix could be differentiated obviously from its adulterants.

CONCLUSION

ITS2 sequence is able to identify Peucedani Radix and its adulterants correctly, which provides a scientific basis for fast and accurate identification of the herb.

[Stability and accuracy of the identification of Notopterygii Rhizoma et Radix using the ITS/ITS2 barcodes]

In this study, Notopterygii Rhizoma et Radix was used to verify the stability and accuracy of DNA barcodes in identification of Chinese materia medica for the first time. All genomic DNAs from thirty one samples were extracted. The ITS (internal transcribed spacer) regions were amplified and sequenced bi-directionally. Obtained sequences were assembled using the CodonCode Aligner. And the sequences of the ITS regions were aligned through Clustal-W and the genetic distances were computed using MEGA 5.0 in accordance with the kimura 2-parameter (K2P) model. The neighbor-joining (NJ) phylogenetic trees were constructed. The ITS2 regions were obtained by using the hidden Markov model (HMM)-based annotation methods from the ITS sequences. Results indicated that the lengths of ITS regions of Notopterygii Rhizoma et Radix were 603-604 bp, while the lengths of ITS2 regions were 228 bp. The haplotypes of ITS/ITS2 regions of Notopterygii Rhizoma et Radix were the same as those of the original plant leaves. The intra-specific genetic distances were smaller than inter-specific ones in ITS/ITS2 regions of Notopterygium incisum and N. franchetii. The NJ trees showed that N. incisum, N. franchetii and its adulterants can be easily differentiated according to their monophyly. Therefore, ITS/ITS2 regions as DNA barcodes can stably and accurately distinguish Notopterygii Rhizoma et Radix from its adulterants and could provide a new technique to ensure clinical safety in utilization of traditional Chinese medicines.

Fingerprinting the Asterid species using subtracted diversity array reveals novel species-specific sequences

Background

Asterids is one of the major plant clades comprising of many commercially important medicinal species. One of the major concerns in medicinal plant industry is adulteration/contamination resulting from misidentification of herbal plants. This study reports the construction and validation of a microarray capable of fingerprinting medicinally important species from the Asterids clade.

Methodology/Principal Findings

Pooled genomic DNA of 104 non-asterid angiosperm and non-angiosperm species was subtracted from pooled genomic DNA of 67 asterid species. Subsequently, 283 subtracted DNA fragments were used to construct an Asterid-specific array. The validation of Asterid-specific array revealed a high (99.5%) subtraction efficiency. Twenty-five Asterid species (mostly medicinal) representing 20 families and 9 orders within the clade were hybridized onto the array to reveal its level of species discrimination. All these species could be successfully differentiated using their hybridization patterns. A number of species-specific probes were identified for commercially important species like tea, coffee, dandelion, yarrow, motherwort, Japanese honeysuckle, valerian, wild celery, and yerba mate. Thirty-seven polymorphic probes were characterized by sequencing. A large number of probes were novel species-specific probes whilst some of them were from chloroplast region including genes like atpB, rpoB, and ndh that have extensively been used for fingerprinting and phylogenetic analysis of plants.

Conclusions/Significance

Subtracted Diversity Array technique is highly efficient in fingerprinting species with little or no genomic information. The Asterid-specific array could fingerprint all 25 species assessed including three species that were not used in constructing the array. This study validates the use of chloroplast genes for bar-coding (fingerprinting) plant species. In addition, this method allowed detection of several new loci that can be explored to solve existing discrepancies in phylogenetics and fingerprinting of plants.

Microsatellite development for an endangered riparian inhabitant, Lilaeopsis schaffneriana subsp. recurva (Apiaceae)

PREMISE OF THE STUDY

Microsatellite markers were developed and characterized to evaluate genetic diversity and population structure in Lilaeopsis schaffneriana subsp. recurva, an endangered species endemic to wetlands dispersed throughout southeastern Arizona, USA, and northern Sonora, Mexico.

METHODS AND RESULTS

Eight loci (one of which was monomorphic) were developed and characterized in 48 individuals from two populations. The total number of alleles was 35, ranging from one to 10 per locus. Many of the primers amplified in L. carolinensis, L. chinensis, L. masonii, L. occidentalis, L. schaffneriana subsp. schaffneriana, Oxypolis fendleri, and Eryngium lemmonii.

CONCLUSIONS

Development of these novel microsatellite loci will facilitate a deeper understanding of genetic diversity, mode of reproduction, and population structure not only in L. schaffneriana subsp. recurva, but also in apiaceous relatives.

Progenitor-derivative speciation in Pozoa (Apiaceae, Azorelloideae) of the southern Andes

BACKGROUND AND AIMS

Studies examining patterns and processes of speciation in South America are fewer than in North America and Europe. One of the least well documented processes has been progenitor-derivative speciation. A particularly instructive example occurs in the southern Andes in the genus Pozoa (Apiaceae, Azorelloideae), which consists of only two diploid outcrossing species, the widespread P. coriacea and the geographically and ecologically restricted P. volcanica. This paper tests the hypothesis that the latter species originated from the former through local geographical and ecological isolation by progenitor-derivative speciation.

METHODS

DNA sequences were analysed from Pozoa and the related South American genera Asteriscium, Eremocharis and Gymnophyton from non-coding regions of the plastid genome, ndhF-rpl32 and rpl32-trnL, plus incorporation of previously reported rpl16 intron and trnD-trnT intergenic spacer sequences. Amplified fragment length polymorphism (AFLP) data from 105 individuals in 21 populations throughout the entire range of distribution of the genus were used for estimation of genetic diversity, divergence and SplitsTree network analysis. Ecological factors, including habitat and associated species, were also examined.

KEY RESULTS

Pozoa coriacea is more similar genetically to the outgroup genera, Asteriscium and Eremocharis, than is P. volcanica. At the population level, only P. volcanica is monophyletic, whereas P. coriacea is paraphyletic. Analyses of genetic differentiation among populations and genetic divergence and diversity of the species show highest values in P. coriacea and clear reductions in P. volcanica. Pozoa coriacea occurs in several types of high elevation habitats, whereas P. volcanica is found only in newly formed open volcanic ash zones.

CONCLUSIONS

All facts support that Pozoa represents a good example of progenitor-derivative speciation in the Andes of southern South America.

Cryptic spatial aggregation of the cushion plant Azorella selago (Apiaceae) revealed by a multilocus molecular approach suggests frequent intraspecific facilitation under sub-Antarctic conditions

PREMISE OF THE STUDY

In abiotically severe habitats, intraspecific aggregations can increase species' fitness by ameliorating stressful environmental factors. However, the difficulty of identifying individual plants in some growth forms makes the measurements of intraspecific aggregation, and therefore the assessment of intraspecific facilitation, problematic. In this study, we examined the genotype composition within cushions of Azorella selago, a sub-Antarctic cushion plant, to investigate the potential extent of intraspecific facilitation.

METHODS

The study was performed on Marion Island, South Africa. Two to eight samples were collected from 42 A. selago cushions, comprising eight different growth forms. Samples were genotyped using seven microsatellite markers.

KEY RESULTS

We showed that all cushion shapes, with the exception of small cushions, may be comprised of more than one genetically distinct individual.

CONCLUSIONS

Under harsh sub-Antarctic conditions, intraspecific aggregation between A. selago individuals appears common and may be driven by the positive impacts of environmental amelioration.

RoBuST: an integrated genomics resource for the root and bulb crop families Apiaceae and Alliaceae

BACKGROUND

Root and bulb vegetables (RBV) include carrots, celeriac (root celery), parsnips (Apiaceae), onions, garlic, and leek (Alliaceae)--food crops grown globally and consumed worldwide. Few data analysis platforms are currently available where data collection, annotation and integration initiatives are focused on RBV plant groups. Scientists working on RBV include breeders, geneticists, taxonomists, plant pathologists, and plant physiologists who use genomic data for a wide range of activities including the development of molecular genetic maps, delineation of taxonomic relationships, and investigation of molecular aspects of gene expression in biochemical pathways and disease responses. With genomic data coming from such diverse areas of plant science, availability of a community resource focused on these RBV data types would be of great interest to this scientific community.

DESCRIPTION

The RoBuST database has been developed to initiate a platform for collecting and organizing genomic information useful for RBV researchers. The current release of RoBuST contains genomics data for 294 Alliaceae and 816 Apiaceae plant species and has the following features: (1) comprehensive sequence annotations of 3663 genes 5959 RNAs, 22,723 ESTs and 11,438 regulatory sequence elements from Apiaceae and Alliaceae plant families; (2) graphical tools for visualization and analysis of sequence data; (3) access to traits, biosynthetic pathways, genetic linkage maps and molecular taxonomy data associated with Alliaceae and Apiaceae plants; and (4) comprehensive plant splice signal repository of 659,369 splice signals collected from 6015 plant species for comparative analysis of plant splicing patterns.

CONCLUSIONS

RoBuST, available at http://robust.genome.com, provides an integrated platform for researchers to effortlessly explore and analyze genomic data associated with root and bulb vegetables.

A comparison of nrDNA ITS and ETS loci for phylogenetic inference in the Umbelliferae: an example from tribe Tordylieae

The Umbelliferae is a large and taxonomically complex family of flowering plants whose phylogenetic relationships, particularly at low taxonomic levels, are generally obscure based on current and widely used molecular markers. Thus, information on the phylogenetic utility of additional molecular markers at these levels is highly favorable. We investigate the utility of nuclear ribosomal DNA (nrDNA) external transcribed spacer (ETS) sequences for phylogenetic inference in Umbelliferae tribe Tordylieae, a group whose relationships have been previously difficult to resolve owing to low sequence variability, and compare the results to those obtained from the nrDNA internal transcribed spacer (ITS) region. We report that the ETS region evolves at a slightly faster rate and has a higher percentage of parsimony informative characters than that of ITS and all chloroplast DNA loci examined to date. The ETS region is a valuable phylogenetic marker in Umbelliferae for low level analysis, especially when used in combination with ITS.

[Genetic diversity of Changium smyrnioides based on SRAP]

OBJECTIVE

To study the genetic diversity of Changium smyrnioides and give a reference for utilization of the germplasm.

METHOD

Ten different populations of Ch. smyrnioides were analyzed by the approach of sequence-related amplified polymorphism (SRAP). Genetic similarity coefficient was calculated, and systematic relationships were constructed based on the UPGMA method.

RESULT

Seventeen primer pairs were selected from 160. A total of 363 bands were scored, 314 bands of them were polymorphic and the average was 18.47 polymorphic bands per primer pair, which were up to 86.50% polymorphic ratio. The results indicated that there was abundant genetic diversity among the tested materials. Genetic similarity coefficient was ranged from 0.4959 to 0.8182. Cluster analysis showed that ten different populations of Ch. smyrnioides could be distinguished into two groups.

CONCLUSION

High level genetic diversity was in different populations of Ch. smyrnioides, and genetic relationship was correlative to geographic position.

Cryopreservation of in vitro shoot apices of Glehnia littoralis - a medical plant

Cryopreservation was examined as a practical method for preserving the genetic resources of Glehnia littoralis Fr. Schm. a nearly exterminated medical plant. In vitro shoot apices were successfully cryopreserved using vitrification and encapsulation-dehydration. In vitrification, regrowth of apices loaded for 60 min with PVS2 (66.7 +/- 6.7 percent) was preferable to loading with PVS3 (50.0 +/- 5.8 percent). With encapsulation-dehydration, the best regrowth (43.3 +/- 3.3 percent) was achieved when the moisture content in the beads was reduced to 19 percent by drying with silica gel for 6 h. Increased regrowth of shoot tips cryopreserved by encapsulation-dehydration resulted from the addition of 0.1 mM acetylsalicylic acid to the loading solution (86.7 +/- 3.3 percent).

[Study of the Saposhnikovia divaricata with space mutagenesis by Fourier transform infrared spectroscopy]

The outer space Saposhnikovia divaricata of the fourth generation was studied and compared with the ground group and comparison group for the first time. The intensity of carbonyl absorption peak at 1640 cm(-1) is stronger than the other two, indicating that chromone content was enhanced obviously. The peaks at 2927 and 2856 cm(-1) were assigned to --CH2-- groups, the peak at 1050 cm(-1) was due to the C-O groups, and the intensities at 2927, 2856 and 1050 cm(-1) were stronger than other samples, indicating that polysaccharide and glycoside contents were enhanced obviously. The intensity of lactone absorption peak at 1743 cm(-1) was stronger than the ground group and weaker than the comparison, indicating that coumarin content was higher than ground group and lower than the comparison. The major components and the structures remained intact, and the effective component contents were enhanced in the outer space Saposhnikovia divaricata. Using FTIR can analyse and characterize the intrinsic quality of the outer space medicinal materials.

[Molecular characteristic marker method of 27 breeds of common Umbelliferae Chinese herb medicinal plants by sequencing rDNA]

OBJECTIVE

To probe a molecular marker method of accrediting fingerprinting of 27 kinds of common Umbelliferae Chinese herb medicinal plants by sequencing rDNA.

METHOD

The rDNA sequences of the 27 breeds of common Umbelliferae Chinese herb medicinal plants were amplified, and were digested by restriction endonuclease, and were seperated via polypropylene electrophoresis, at last 6 breeds of them were sequenced.

RESULTS

The rDNA sequence fragment we gained concluded ITS1, ITS2, 5.8S complete sequence and 18S, 26S part sequence. On the electrophoresis map of PCR products digested by restriction endonuclease MSP I, 27 breeds appeared 16 kinds of characteristic map, 11 of them differ from others; and PCR products digested by restriction endonuclease HaeIII, there appeared 5 kinds of characteristic map among 27 breeds, 3 of them differ from others. The sequenced result of 6 breeds showed genes whose length extented from 652bp to 656bp were acquired. These sequences of 3 breeds which showed the same electrophoresis map after digested by restriction endonuclease HaeIII exhibited great similarity according to similar phylogenetic tree constructed on rDNA sequence.

CONCLUSION

The rDNA sequence character is effective molecular marker for classifying the different Umbellerae Chinese herb medicinal plants. And the method of sequencing rDNA surpassed that of restriction fragment long polymorphism (RFLP).

[Study on molecular evidence of revision of taxonomic placement of Peucedanum decursivum]

OBJECTIVE

To explain the molecular evidence of revision of taxonomic placement of Peucedanum decursivum based on the nrDNA ITS sequence.

METHOD

PCR amplification, DNA sequencing and cladistic analysis.

RESULT

The ITS sequences and phylogenetic tree of 5 species of Angelica were and Peucedanum were acquired, in which 5 species were divided into 2 groups, Angelica group and Peucedanum group. P. decursivum was placed in the Angelica group.

CONCLUSION

P. decursivum belongs to genus Angelica. The scientific name of P. decursivum should be revised as A. decursivum. A. decursivum and P. praeruptorum should be used as crude drug respectively.

Single seed Raman measurements allow taxonomical discrimination of Apiaceae accessions collected in gene banks

NIR-FT-Raman spectroscopy was applied for a nondestructive analysis of single seeds (fruit mericarps) of 36 accessions belonging to various species of the Apiaceae family. Main seed components such as fatty acids, polysaccharides, proteins, and lignin were identified based on the obtained Raman spectra. Variation at the species and genus level was related to differences observed between spectra. The application of cluster analysis discriminated among most of the species evaluated and grouped them according to their taxonomical classification. The spectroscopically analyzed seeds germinated and developed into normal seedlings to demonstrate the additional advantage that Raman spectroscopy is nondestructive and can be applied to living seed without harm. These results indicate that Raman spectroscopy is a valuable tool for the rational evaluation and management of genetic resources in ex situ seed collections by providing useful information for taxonomical validation of the accessions.

Molecular evolution of flavonoid dioxygenases in the family Apiaceae

Plant species of the family Apiaceae are known to accumulate flavonoids mainly in the form of flavones and flavonols. Three 2-oxoglutarate-dependent dioxygenases, flavone synthase or flavanone 3 beta-hydroxylase and flavonol synthase are involved in the biosynthesis of these secondary metabolites. The corresponding genes were cloned recently from parsley (Petroselinum crispum) leaves. Flavone synthase I appears to be confined to the Apiaceae, and the unique occurrence as well as its high sequence similarity to flavanone 3beta-hydroxylase laid the basis for evolutionary studies. In order to examine the relationship of these two enzymes throughout the Apiaceae, RT-PCR based cloning and functional identification of flavone synthases I or flavanone 3beta-hydroxylases were accomplished from Ammi majus, Anethum graveolens, Apium graveolens, Pimpinella anisum, Conium maculatum and Daucus carota, yielding three additional synthase and three additional hydroxylase cDNAs. Molecular and phylogenetic analyses of these sequences were compatible with the phylogeny based on morphological characteristics and suggested that flavone synthase I most likely resulted from gene duplication of flavanone 3beta-hydroxylase, and functional diversification at some point during the development of the apiaceae subfamilies. Furthermore, the genomic sequences from Petroselinum crispum and Daucus carota revealed two introns in each of the synthases and a lack of introns in the hydroxylases. These results might be explained by intron losses from the hydroxylases occurring at a later stage of evolution.

[Analysis of population genetic structure and molecular identification of Changium smyrnioides and Chuanminshen violaceum with ISSR marker]

OBJECTIVE

To assess the population genetic diversity and genetic structure and screen species-specific bands for identification of Changium smyrnioides and Chuanminshen violaceum.

METHOD

Seven wild populations of Changium smyrnioides and one cultivated population of Chuanminshen violaceum were studied by ISSR analysis. The population genetic diversity and population genetic structure were assessed by using POPGENE software.

RESULT

A total of 152 ISSR markers were scored, among which 136 (90.8%) were polymorphic. The values of Gst tended to be high (mean Gst = 0.575). The level of genetic divesity of Changium smyrnioides (A = 1.272; P = 27.26%; I = 0.132; H = 0.087) was higher than that of Chuanminshen violaceum (A = 1.217; P = 21.7; I = 0.103; H = 0.067).

CONCLUSION

The genetic variation of Changium smyrnioides is high and the majority of genetic variation occur among populations. Substantial genetic divergence is shown by cluster analysis (UPGMA) to befound between Changium smyrnioides and Chuanminshen violaceum at DNA level. In addition, one species-specific marker has been obtained in Chuanminshen violaceum. The phylogenetic relationship of two species has also been discussed.