Development and characterization of 16 novel microsatellite markers by Transcriptome sequencing for Angelica dahurica and test for cross-species amplification

Molecular evidence provides new insights into the evolutionary origin of an ancient traditional Chinese medicine, the domesticated "Baizhi"

Introduction

"Baizhi" is a famous herbal medicine in China, and it includes four landraces named as 'Hangbaizhi', 'Chuanbaizhi', 'Qibaizhi', and 'Yubaizhi'. Long-term artificial selection had caused serious degradation of these germplasms. Determining the wild progenitor of the landraces would be benefit for their breed improvements. Previous studies have suggested Angelica dahurica var. dahurica, A. dahurica var. formosana, or A. porphyrocaulis as potential candidates, but the conclusion remains uncertain, and their phylogenetic relationships are still in controversy.

Methods

In this study, the genetic variation and phylogenetic analyses of these species and four landraces were conducted on the basis of both the nrITS and plastome datasets.

Results

Genetic variation analysis showed that all 8 population of four landraces shared only one ITS haplotype, meanwhile extremely low variation occurred within 6 population at plastid genome level. Both datasets supported the four landraces might be originated from a single wild germplasm. Phylogenetic analyses with both datasets revealed largely consistent topology using Bayesian inference and Maximum likelihood methods. Samples of the four landraces and all wild A. dahurica var. dahurica formed a highly supported monophyletic clade, and then sister to the monophyly clade comprised by samples of A. porphyrocaulis, while four landraces were clustered into one clade, which further clustered with a mixed branches of A. porphyrocaulis and A. dahurica var. dahurica to form sister branches for plastid genomes. Furthermore, the monophyletic A. dahurica var. formosana was far distant from the A. dahurica var. dahurica-"Baizhi" clade in Angelica phylogeny. Such inferences was also supported by the evolutionary patterns of nrITS haplotype network and K2P genetic distances. The outcomes indicated A. dahurica var. dahurica is most likely the original plant of "Baizhi".

Discussion

Considering of phylogenetic inference and evolutionary history, the species-level status of A. dahurica var. formosana should be accepted, and the taxonomic level and phylgenetic position of A. porphyrocaulis should be further confirmed. This study preliminarily determined the wild progenitor of "Baizhi" and clarified the phylogenetic relationships among A. dahurica var. dahurica, A. dahurica var. formosana and A. porphyrocaulis, which will provide scientific guidance for wild resources protections and improvement of "Baizhi".

Genetic variation and structure of endemic and endangered wild celery (Kelussia odoratissima Mozaff.) quantified using novel microsatellite markers developed by next-generation sequencing

Kelussia odoratissima Mozaff. (Apiaceae) is a native plant that has been traditionally consumed in Iran's food and pharmaceutical industries. Overharvesting of the taxon, especially at the beginning of the growing season, due to its considerable medicinal and economic value, is believed to be the main reason for the extirpating of this plant. The consequences of the severe anthropogenic impacts on the genetic diversity of populations are poorly known. In order to investigate the level of genetic variation and patterns of the genetic structure of K. odoratissima, we developed novel microsatellite markers using the 454 Roche next-generation sequencing (NGS) platform for the first time. Out of 1,165 microsatellite markers bioinformatically confirmed, twenty-five were tested, of which 23 were used to screen genetic variation across 12 natural populations. Our results showed that the average number of alleles per locus and the polymorphic information content (PIC) were 10.87 (range 7 to 27), and 0.81 (range 0.67 to 0.94), respectively. The mean observed and expected heterozygosities (± SD) across all populations were 0.80 ± 0.31 and 0.72 ± 0.14, respectively. The average pairwise FST among the populations was 0.37 (range 0.04 to 0.81). Bayesian and distance-based clustering, and principal coordinate analyses revealed at least four major genetic clusters. Although high level of structure can be explained by landscape topography and geographic distance, presence of admixed populations can be associated to seed or pollen dispersal. Contrary to expectations, the high level of genetic variation and lack of inbreeding suggest that overexploitation has not yet significantly purged the allelic variability within the natural populations in protected areas.

Transcriptome analysis provides insights into coumarin biosynthesis in the medicinal plant Angelica dahurica cv. Yubaizhi

Angelica dahurica roots have a long history of use in traditional Chinese medicine due to their high coumarin content. To address the increasing demand for these roots, a synthetic biology approach has been proposed. Nevertheless, our comprehension of coumarin biosynthesis and its regulation remains limited. In this study, we utilized Hiseq2500 sequencing to analyze the transcriptomes of A. dahurica at different growth stages while concurrently quantifying coumarin content. Differentially expressed gene (DEG) analysis was employed to identify key genes involved in coumarin and terpenoid backbone biosynthesis. Weighted gene co-expression network analysis (WGCNA) was applied to identify gene modules strongly associated with coumarin content, elucidating the regulatory relationships between transcription factors (TFs) and pathway genes. Furthermore, KEGG enrichment analysis was used to explore essential pathways governing coumarin biosynthesis, with the identification of hub genes. Our results indicated that total coumarin content was highest in the roots, followed by leaves and stems, across all three developmental stages. Transcriptome analysis identified a total of 92,478 genes, among which 215 and 30 genes were implicated in coumarin and terpenoid backbone biosynthesis, respectively. Within the 73 identified gene modules by WGCNA, three modules-namely aquamarine1 (comprising two OMTs, one CSE, one AACT, one HDS, two PSs, one 2OGO, four UGTs, and seven CYP450s), darkmagenta (containing one UGT and one HDR), and navajowhite2 (consisting of one HCT, three UGTs, one CYP71A25, one OMT, one CSE, one HDS, and one PT)-were strongly associated with imperatorin, oxypeucedanin, and isoimperatorin content, respectively. KEGG enrichment analysis highlighted significant enrichment of cytochrome P450, transporter, and ubiquitin system pathways. Moreover, TF-gene regulatory analysis unveiled the complexity of coumarin biosynthesis, with 17 TF families regulating 17 genes in the aquamarine1 module, 8 TF families regulating 2 genes in the darkmagenta module, and 8 TF families regulating 7 genes in the navajowhite2 module. These comprehensive findings provide valuable insights into coumarin biosynthesis in A. dahurica, facilitating future research and potential applications in traditional Chinese medicine and synthetic biology strategies.

The novel developed microsatellite markers revealed potential hybridization among Cymbidium species and the interspecies sub-division of C. goeringii and C. ensifolium

BACKGROUND

Orchids (Cymbidium spp.) exhibit significant variations in floral morphology, pollinator relations, and ecological habitats. Due to their exceptional economic and ornamental value, Cymbidium spp. have been commercially cultivated for centuries. SSR markers are extensively used genetic tools for biology identification and population genetics analysis.

RESULT

In this study, nine polymorphic EST-SSR loci were isolated from Cymbidium goeringii using RNA-Seq technology. All nine SSR loci showed transferability in seven other congeneric species, including 51 cultivars. The novel SSR markers detected inter-species gene flow among the Cymbidium species and intra-species sub-division of C. goeringii and C. ensifolium, as revealed by neighborhood-joining and Structure clustering analyses.

CONCLUSION

In this study, we developed nine microsatellites using RNA-Seq technology. These SSR markers aided in detecting potential gene flow among Cymbidium species and identified the intra-species sub-division of C. goeringii and C. ensifolium.

High genetic diversity and low population differentiation of a medical plant Ficus hirta Vahl., uncovered by microsatellite loci: implications for conservation and breeding

BACKGROUND

Wuzhimaotao (Radix Fici Hirtae) originates from the dry root of Ficus hirta (Moraceae), which is widely known as a medical and edible plant distributed in South China. As the increasing demand for Wuzhimaotao, the wild F. hirta has been extremely reduced during the past years. It is urgent to protect and rationally develop the wild resources of F. hirta for its sustainable utilization. However, a lack of genetic background of F. hirta makes it difficult to plan conservation and breeding strategies for this medical plant. In the present study, a total of 414 accessions of F. hirta from 7 provinces in southern China were evaluated for the population genetics using 9 polymorphic SSR markers.

RESULTS

A mean of 17.1 alleles per locus was observed. The expected heterozygosity (He) varied from 0.142 to 0.861 (mean = 0.706) in nine SSR loci. High genetic diversity (H_e = 0.706, ranged from 0.613 to 0.755) and low genetic differentiation among populations (G'_ST = 0.147) were revealed at population level. In addition, analysis of molecular variance (AMOVA) indicated that the principal molecular variance existed within populations (96.2%) was significantly higher than that among populations (3.8%). Meanwhile, the three kinds of clustering methods analysis (STRUCTURE, PCoA and UPGMA) suggested that the sampled populations were clustered into two main genetic groups (K = 2). Mantel test showed a significant correlation between geographic and genetic distance among populations (R² = 0.281, P < 0.001). Pollen flow, seed flow and/or geographical barriers might be the main factors that formed the current genetic patterns of F. hirta populations.

CONCLUSIONS

This is a comprehensive study of genetic diversity and population structure of F. hirta in southern China. We revealed the high genetic diversity and low population differentiation in this medicinal plant and clarified the causes of its current genetic patterns. Our study will provide novel insights into the exploitation and conservation strategies for F. hirta.

Limited genetic diversity and high differentiation in Angelica dahurica resulted from domestication: insights to breeding and conservation

BACKGROUND

Angelica dahurica belongs to the Apiaceae family, whose dry root is a famous traditional Chinese medicine named as "Bai zhi". There are two cultivars (A. dahurica cv. 'Hangbaizhi' and A. dahurica cv. 'Qibaizhi'), which have been domesticated for thousands of years. Long term artificial selection has led to great changes in root phenotypes of the two cultivars, and also decreased their adaptability to environment. We proposed hypothesis that the cultivars may have lost some of the genetic diversity found in the wild species and may be highly differentiated from the latter during the domestication process. However, few studies have been carried out on how domestication affected the genetic variation of this species. Here, we accessed the levels of genetic variation and differentiation within and between wild A. dahurica populations and two cultivars using 12 microsatellite markers.

RESULTS

The results revealed that the genetic diversity of the cultivars was much lower than that of wild A. dahurica, and A. dahurica cv. 'Qibaizhi' had lower genetic diversity compared to A. dahurica cv. 'Hangbaizhi'. AMOVA analysis showed significant genetic differentiation between the wild and cultivated A. dahurica populations, and between A. dahurica cv. 'Hangbaizhi' and A. dahurica cv. 'Qibaizhi'. Results from Bayesian, UPGMA, NJ and PcoA clustering analysis indicated that all 15 populations were assigned to two genetic clusters corresponding to the wild and cultivated populations. Bayesian clustering analysis further divided the cultivated populations into two sub-clusters corresponding to the two cultivars.

CONCLUSIONS

Our study suggests that the domestication process is likely the major factor resulting in the loss of genetic diversity in cultivated A. dahurica populations and in significant genetic differentiation from the wild populations due to founder effect and/or artificially directional selections. This large-scale analysis of population genetics could provide valuable information for genetic resources conservation and breeding programs of Angelica dahurica.