Next Article in Journal
Heterozygosity of ALG9 in Association with Autosomal Dominant Polycystic Liver Disease
Previous Article in Journal
Novel HSPG2 Gene Mutation Causing Schwartz–Jampel Syndrome in a Moroccan Family: A Literature Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Population Genetic Analysis of Paris polyphylla var. yunnanensis Based on cpDNA Fragments

1
The Chongqing Engineering Laboratory for Green Cultivation and Deep Processing of Three Gorges Reservoir Area’s Medicinal Herbs, College of Food and Biology Engineering, Chongqing Three Gorges University, Chongqing 404120, China
2
College of Pharmacy, Dali University, Dali 671000, China
3
Research Center of Biomass Clean Utilization, Beijing Key Laboratory of Lignocellulosic Chemistry, College of Materials Science and Technology, Beijing Forestry University, Beijing 100083, China
*
Authors to whom correspondence should be addressed.
Genes 2023, 14(9), 1754; https://doi.org/10.3390/genes14091754
Submission received: 21 July 2023 / Revised: 23 August 2023 / Accepted: 25 August 2023 / Published: 2 September 2023
(This article belongs to the Special Issue Molecular Genetics and Resources of Medicinal Plants)

Abstract

:
Paris polyphylla var. yunnanensis is a well-known medicinal plant that is mainly distributed in Southwest China; however, its genetic diversity and biodiversity processes are poorly understood. In this study, the sequences of cpDNA trnL-trnF fragments of 15 wild populations and 17 cultivated populations of P. polyphylla var. yunnanensis were amplified, sequenced, and aligned to study the population genetics of this species. Genetic diversity was analyzed based on nucleotide diversity, haplotype diversity, Watterson diversity, population-level diversity, and species-level genetic diversity. Genetic structure and genetic differentiation were explored using haplotype distribution maps and genetic distance matrices. A total of 15 haplotypes were identified in the 32 populations of P. polyphylla var. yunnanensis. Five unique haplotypes were identified from the fourteen haplotypes of the cultivated populations, while only one unique haplotype was identified from the ten haplotypes of the wild populations. The haplotype richness and genetic diversity of the cultivated populations were higher than those of the wild populations (HT = 0.900 vs. 0.861). In addition, there were no statistically significant correlations between geographic distance and genetic distance in the cultivated populations (r = 0.16, p > 0.05), whereas there was a significant correlation between geographical distance and genetic structure in the wild populations (r = 0.32, p > 0.05), indicating that there was a geographical and genetic connection between the wild populations. There was only 2.5% genetic variation between the wild populations and cultivated populations, indicating no obvious genetic differentiation between the wild and cultivated populations. Overall, the genetic background of the cultivated populations was complex, and it was hypothesized that the unique haplotypes and higher diversity of the cultivated populations were caused by the mixed provenance of the cultivated populations.

1. Introduction

P. polyphylla var. yunnanensis is a perennial herb whose stem is generally 20–30 cm in length and purplish red in color, and the plant has a membranous leaf sheath at the base that embraces the stem. The pedicel protrudes from the top of the stem and bears a terminal flower. The flowering period is from June to July, and the fruiting period is from September to October. This plant is mainly found in Southwest China (provinces of Sichuan, Guizhou, Yunnan, and Guangxi); prefers acidic or slightly acidic soil with rich humus layers; and is suitable for growing in evergreen broad-leaved forests, bamboo forests, or shrublands with temperatures in the range of 16–28 °C and elevations ranging from 1400 m to 3100 m [1]. Pharmacological studies of P. polyphylla var. yunnanensis have shown that the main functional components of this plant are Paris saponin and diosgenin, which have demonstrated antitumor, antibacterial, antiviral, and hemostatic functions [2,3]. P. polyphylla var. yunnanensis has high medicinal and economic value and is the natural raw material that comprises 81 different Chinese patent medicines, such as “Yunnan Baiyao” and “Gongxuening Capsule” [4], corroborating its role in the development of traditional Chinese medicines (TCMs) for tumor treatment, immunomodulation, and others [5]. P. polyphylla var. yunnanensis is listed in the Pharmacopoeia of the People’s Republic of China as one of the source plants of Paris polyphylla. However, wild P. polyphylla var. yunnanensis suffers from slow growth, and it has been overexploited, making it difficult to obtain sufficient numbers of seedlings. Many farms collect wild P. polyphylla var. yunnanensis and cultivate its samples to meet the growing market demand. Despite its popularity, few studies have analyzed the genetic diversity of wild and cultivated P. polyphylla var. yunnanensis, which is important to its future survival. As a consequence, it is urgent to understand the genetic diversity and genetic structure of P. polyphylla var. yunnanensis.
The harsh and unstable climate of the Quaternary Ice Age caused large-scale migration of organisms that led to disastrous survival rates. After the Ice Age, the glacial refugium provided a refuge for organisms to survive and ultimately rediffuse to other geographical regions because of its unique geographical conditions. P. polyphylla var. yunnanensis also underwent this migration, ultimately surviving and differentiating well in Southwest China. As a result, the genetic diversity and structure of P. polyphylla var. yunnanensis was deeply affected by the unique geographical conditions of Southwest China. Some studies have sought to identify the molecular markers that elucidate the genetic diversity of P. polyphylla var. yunnanensis. A total of 153 individuals from six populations (three wild populations and three cultivated populations) in the Sichuan and Yunnan provinces (China) were investigated based on ISSR molecular markers. The genetic diversity of the cultivated populations was higher than that of the wild populations (0.153 vs. 0.151), and there was no obvious genetic differentiation between the cultivated populations and the wild populations [6]. The genetic diversity of 62 individuals from the Yunnan province was explored based on SSR molecular markers. Ten pairs of primers were used as single-sequence-repeat (SSR) markers, and a heterozygosity range of 0.790–0.976 was observed following the analysis [7]. SSR molecular markers were also utilized to analyze the genetic diversity of 115 samples from five populations in Yunnan province, and the as-expected heterozygosity was 0.7744 at the species level and 0.6548 at the population level, indicating a high genetic diversity level of P. polyphylla var. Yunnanensis [8]. Using amplified fragment length polymorphism (AFLP) molecular markers, the genetic diversity at the species level and population level was 0.2768 and 0.1821, respectively, indicating that the 15 populations of P. polyphylla var. Yunnanensis in the Yunnan province had a low genetic diversity [9]. AFLPs were used in the genetic analysis of 15 wild and 17 cultivated populations of P. polyphylla var. yunnanensis [10], which revealed that the cultivated populations had higher genetic diversity than the wild populations at the species level (HE = 0.2636 vs. 0.2616, respectively).
Chloroplast genome (cpDNA) is maternally inherited in most angiosperms, while it is paternally inherited in most gymnosperms, meaning it undergoes almost no recombination. Compared to nuclear genomic markers inherited from both parents, cpDNA markers have the ability to reveal genetic differences in small population sizes. Therefore, cpDNA markers serve as an effective tool for studying phylogeography and elucidating the migration and diffusion routes of species in shelters and post-glacial periods [11,12]. Although some population genetic studies of P. polyphylla var. yunnanensis had been reported, earlier studies that focused on the Yunnan province had a narrow research area and only studied a small number of populations, which did not effectively reflect the genetic diversity and genetic structure of P. polyphylla var. yunnanensis at the entire species level. In addition, information regarding the phylogenetic structure of P. polyphylla var. yunnanensis based on cpDNA and nuclear genes has not yet been reported.
The purpose of the present study was to analyze the genetic diversity and structure of 15 wild populations and 17 cultivated populations of P. polyphylla var. yunnanensis using cpDNA fragments (trnL–trnF).

2. Materials and Methods

2.1. Materials and Reagents

In the trnL–trnF experiment using cpDNA fragments, the as-used materials were the same as AFLP markers [10]. There were 364 effective individuals from 32 populations, including 15 wild and 17 cultivated populations.

2.2. PCR Amplification and Sequencing of cpDNA trnL–trnF Fragments

2.2.1. Total DNA Extraction and Purification

The DNA from approximately 0.5 g of dried P. polyphylla var. yunnanensis leaves from each sample was extracted using cetyltrimethylammonium bromide (CTAB) according to a previously published (and improved) method [13].

2.2.2. PCR Amplification and Sequencing of cpDNA trnL-trnF Fragments

PCR amplification and sequencing were performed on 364 individuals using screened primers of cpDNA trnL-trnF (base sequence of C: 5′-CGAAATCGGTAGACGCTACG-3′, F: 5′-ATTTGAACTGGTGACACGAG-3′). PCR amplification of the trnL-trnF fragments was performed in a 20 µL reaction system using a K960 PCR apparatus with the following thermocycling program: 95 °C for 5 min, 95 °C for 1 min, 58 °C for 1 min, 72 °C for 1 min, and 33 cycles at 72 °C for 10 min; afterward, the PCR products were stored at 4 °C. The PCR products were then electrophoresed over a 1% agarose gel, and the clear bands with no impurities were selected and sent to Beijing Boyoushun Biotechnology Co., Ltd. (Beijing, China) for sequencing.

2.3. DNA Sequence Splicing and Comparison

SeqMan 17.1(Lasergene Genomics software package, DNASTAR) was used to check the peak map of all individual sequenced sequences. The positive and negative primers were spliced together and saved in the “.seq” format. Then, the chloroplast haplotypes were determined based on the variations in the sequences; the haplotypes in the population were named “H1”, “H2”, “H3”.... Then, all haplotype sequences saved in the “.seq” format were imported into SeqMan again and saved as a FASTA file for comparison and sequencing (the FASTA files of the sequences of all individuals in the wild populations and cultivated populations were formatted the same way). Finally, the FASTA files were imported into the software for comparison, and the matched sequences were saved as FASTA and NEXUS files for subsequent sequence alignment analysis using MEGA-X [14].

2.4. Genetic Diversity Analysis

Genetic diversity parameters based on the DNA sequences included nucleotide diversity (π), haplotype diversity (Hd), Watterson diversity (θ), population-level genetic diversity (HS), and species-level genetic diversity (HT) [15,16]. To determine π, two sequences were randomly selected from the samples, and the average nucleotide variation was calculated based on the average unit site. Hd referred to the frequency of two different haplotype sequences randomly selected from the samples. HS and HT referred to the haplotype diversity at the population level and species level, respectively. The parameters π and Hd were statistically calculated using analysis of molecular variance (AMOVA) tests, while HS and HT were calculated using Permut v.1.2.1 with 1000 permutations [17,18]. The genetic diversities of the wild populations, cultivated populations, and all populations together were calculated.

2.5. Genetic Structure and Genetic Differentiation Analysis

The haplotype distribution maps of the wild and cultivated populations were edited using the ArcGIS software (v. 10.5) to visualize the distribution of the haplotypes in the wild and cultivated populations. The genetic distances between the wild and cultivated populations were calculated using MEGA-X based on Kimura’s 2-parameter model, which sorted out the geographic distance and genetic distance matrix, and the data were analyzed using GenAlEx (v. 6.503) to perform Mantel tests [19,20]. The AMOVA test was conducted using Arlequin (v. 3.5) to analyze the genetic differentiation of the wild and cultivated populations [18], with 1000 substitutions being used for the significance test. The spatial genetic structure of all populations, including the wild and cultivated populations, was analyzed using the SAMOVA 1.0 (Spatial Analysis of Molecular Variance) software. The software clustered the populations, which were similar in geography and genetics, into groups according to the haplotype composition and geographical distribution of the populations. The K values were simulated from 2 to 10 and annealed 100 times. Each K value had a corresponding intergroup genetic differentiation value (FCT); when the FCT value reached a maximum or stabilized at a specific value without a single population, the corresponding K value was considered the best group.
The GST and NST of the wild and cultivated populations were calculated using Permut v. 1.2.1 (1000 permutations). GST is the genetic differentiation coefficient of the frequency of haplotypes (range of 0–1). The smaller the coefficient, the lower the degree of genetic differentiation between populations, indicating that genetic variation mainly exists within each population. NST not only considers the frequency of the haplotypes but also the genetic differentiation coefficient related to the genetic distance between haplotypes. In many cases, these two values are not equal. One of the reasons for the difference may be the genetic and geographical structure of the populations. Therefore, the magnitude of the difference could be used to judge the genetic and geographical relationship between populations. When the value of NST is significantly greater than the value of GST, there is a clear differentiation in the geographical genetic structure between populations. However, when the value of NST is equal to or less than the value of GST, there is no geographical relationship and no differentiation in genetic structure between populations.

3. Results

3.1. Genetic Diversity Analysis of P. polyphylla var. yunnanensis

In this study, the cpDNA trnL–trnF intergenic spacer was used to classify 15 wild populations and 17 cultivated populations of P. polyphylla var. yunnanensis. Based on the sequence alignment, 15 haplotypes were detected in these 32 populations of P. polyphylla var. yunnanensis. Five unique haplotypes were identified from the fourteen haplotypes of the cultivated populations, while only one unique haplotype was identified from the ten haplotypes of the wild populations (Table S1). The wild populations of P. polyphylla var. yunnanensis displayed high genetic diversity at the species level (HT = 0.861) (Table 1), while the genetic diversity at the population level was low (HS = 0.135), resulting in a high genetic differentiation among populations (NST = 0.919 and GST = 0.843). The cultivated populations showed the same trend as the wild populations of P. polyphylla var. yunnanensis (H T = 0.900, HS = 0.222, NST = 0.815, and GST = 0.754). All populations of P. polyphylla var. yunnanensis displayed high genetic diversity at the species level (HT = 0.866), while the genetic diversity at the population level was low (HS = 0.181), resulting in a high genetic differentiation among populations (GST = 0.791 and NST = 0.860). Therefore, the total genetic diversity, haplotype diversity, and nucleotide diversity results indicated that the genetic diversity of the wild populations was slightly lower than the diversity of the cultivated populations, suggesting that artificial domestication in different planting regions increased the diversity of P. polyphylla var. yunnanensis.

3.2. Genetic Structure and Genetic Differentiation Analysis of P. polyphylla var. yunnanensis

In a previous study, the haplotypes determined based on the genetic sequences of the trnL–trnF fragments were used to generate geographical distribution maps of the haplotypes of wild populations [21]. The Guizhou, Central Yunnan, and Western Yunnan provinces had their own haplotypes without crossing. There was a geographical genetic structure among the wild populations. Additionally, the results of the Mantel test analysis showed that there was a significant correlation between the genetic geographical structure of the wild populations (r = 0.3178, p < 0.05) (Figure 1). Similarly, the genetic differentiation coefficients, GST and NST, also indicated the existence of a genetic geographical structure. In the wild populations, the GST was 0.843. Out of 1000 replacements, the NST value for 950 of the replacements was 0.902. The higher value of NST compared to the GST value (p < 0.05) indicated that haplotypes with a similar genetic distance appeared in the same population or populations with a similar geographic distance. Furthermore, out of 1000 replacements, the NST value for 990 of the replacements was 0.913. In this case, the value of NST was also higher than that of GST (p < 0.01), indicating a significant geographical genetic structure among the wild populations (Table 2).
According to the haplotype geographical distribution map of the cultivated populations [21], there was a crossing haplotype (H3) between Songming of Kunming from Yunnan province and Huidong of Sichuan province. According to the geographical structure of the wild populations, there was a geographical isolation between the population in Central Yunnan province and that in Huidong of Sichuan province, and the haplotype of the cultivated populations appeared as a cross haplotype in these two areas. Therefore, there was no obvious geographical genetic structure. This inference was verified by the Mantel test of the cultivated populations, the results of which showed that there was no geographical genetic structure among the cultivated populations (r = 0.1629, p > 0.05) (Figure 1). At the same time, the GST value of the cultivation populations was 0.754. In the 1000-replacement test, the value of NST was 0.821 after 950 replacements (Table 2), which was higher than that of GST (p > 0.05), indicating no obvious geographical genetic structure between the populations of this group. Furthermore, after 990 replacements, the value of NST was 0.832, which was also higher than that of GST (p > 0.01) (Table 2), indicating no obvious geographical genetic structure between the cultivated populations.
Following the SAMOVA analysis of all wild populations and cultivated populations (Figure 2), when the K value increased from 2 to 10, the corresponding FCT value did not stabilize. A maximum FCT value indicated that there was a single population group, which meant that there was no optimal group. Following the AMOVA analysis, genetic differentiation in the wild populations and cultivated populations (76.66% and 78.86%) was observed. However, there were no genetic differences between the wild populations and cultivated populations (2.5%, p > 0.05) (Table 3). In terms of the relative proportions of different chloroplast haplotypes found in the wild and cultivated P. polyphylla var. yunnanensis populations, there were nine haplotypes (H1, H2, H3, H4, H5, H6, H8, H9, and H11) shared between the wild populations and the cultivated populations, one unique haplotype (H13) in the wild populations, and five unique haplotypes (H7, H10, H12, H14, and H15) in the cultivated populations (Figure 3 and Table 4).

4. Discussion

4.1. Genetic Diversity of P. polyphylla var. yunnanensis Based on cpDNA

A higher genetic diversity of plants means they can more strongly adapt to environmental changes, which is important for their survival. Genetic diversity is affected by a number of natural and human factors, such as geographical distribution, habitat destruction, and overexploitation [22]. Because chloroplast genome (cpDNA) is inherited maternally without undergoing genetic recombination, chloroplast haplotype remains almost unchanged between generations. Therefore, cpDNA is an excellent tool to study the genetic diversity of plants. In this study, there was a moderate level of cpDNA diversity (HT = 0.861), which was higher than that of Cedrella odorata (HT = 0.700) [23], Juniperus przewalskii (HT = 0.700) [24], Alnus glutinosa (HT = 0.773), and Cyclobalanopsis glauca (HT = 0.681) [25]. Petit et al. [26] determined the average cpDNA diversity of 170 plants (HT = 0.67) as measured using different molecular markers. However, Cycas taitungensis (HT = 0.998) [27], Cunninghamia lanceolate (HT = 0.952) [28], and Scutellaria baicalensis (HT = 0.888) all had high cpDNA diversity [29]. In this study, the moderate level of genetic diversity of wild P. polyphylla var. yunnanensis as determined by sequencing the cpDNA trnL–trnF intergenic spacing was consistent with the diversity measured using AFLP markers [18]. However, the genetic diversity of P. polyphylla var. yunnanensis was still higher than that of many other endangered plants, indicating that being endangered was not due to its own biological characteristics but rather by overexploitation or other factors.

4.2. Comparison of the Genetic Diversity between the Cultivated and Wild Populations

In this study, we identified 15 haplotypes in P. polyphylla var. yunnanensis, encompassing 10 haplotypes in the wild populations and 14 haplotypes in the cultivated populations (Table S1). Among them, there were nine haplotypes shared between the wild group and the cultivated group, one unique haplotype in the wild group, and five unique haplotypes in the cultivated group (Figure 3). Moreover, the genetic diversity of the cultivated populations was higher than that of the wild populations (HT = 0.900 vs. 0.861). According to the UPGMA tree of all individuals in Yunnan constructed based on AFLP markers [10], the genetic background of the cultivated P. polyphylla var. yunnanensis populations was relatively complex. Therefore, both the haplotype richness and genetic diversity of the cultivated populations were higher than those of the wild populations. However, the correlation between geographical distance and genetic distance was not strong in the cultivation populations, negating the existence of a geographical genetic structure and suggesting that the unique haplotype of the cultivation populations might have been obtained due to the mixed provenance of the cultivation populations.

4.3. Selection of a Conservation Area for P. polyphylla var. yunnanensis

Wild P. polyphylla var. yunnanensis was considered an endangered species in the industrial age. Currently, P. polyphylla var. yunnanensis is classified as a national second-class endangered medicinal plant, and the development and sustainable applications of its wild resources are of far-reaching significance in order to formulate effective protection strategies. We found that P. polyphylla var. yunnanensis was widely distributed throughout Southwest China, but there were few wild individuals collected in some areas. For example, only six and four individuals were collected in Longli county (LL) and Qishe town (XY), respectively, in Xingyi city. Seven and nine individuals were collected in Yimen County and Changning County (CN), respectively, in Yunnan Province. The genetic diversity of 15 wild populations from three provinces was analyzed and compared by sequencing cpDNA fragments, the results of which showed that the wild populations in the Guizhou and Sichuan provinces had high genetic diversity; however, more studies need to be carried out in various populations.
Previously published evolution network diagrams of wild P. polyphylla var. yunnanensis showed that the H1 and H6 haplotypes are ancient haplotypes [10,21]. The genealogical analysis of wild P. polyphylla var. yunnanensis in this study indicated that Guizhou and Western Yunnan provinces might have been two independent glacial refugia, suggesting that priority should be given to these areas to protect wild P. polyphylla var. yunnanensis in these regions and strictly prohibit the excavation of land harboring wild P. polyphylla var. yunnanensis to ensure the protection of its habitat.
P. polyphylla var. yunnanensis is an important naturally growing TCM. In recent years, it has been rapidly cultivated in large quantities to meet the world’s demand for traditional Chinese medicine [30]. It takes as long as ten years for wild individuals to grow from seed germination to harvest. To protect wild resources and ensure sustainable utilization, artificial cultivation has been an effective alternative method [31]. Traditional cultivation methods entail the collection of seeds from wild populations and then planting them in similar habitats, which is an effective way to maintain the gene pool of medicinal plants [6,32]. From our survey, this traditional approach is used in only some small remote mountainous areas, while rare sites may be preserved in local distribution, such as Huidong and Sichuan (HD). Therefore, it is equally important to establish a wild P. polyphylla var. yunnanensis germplasm resource nursery for systematic relocation protection while taking local protection measures. Although we observed that these populations had high genetic diversity, hybrid provenances might produce an uneven quality of P. polyphylla var. yunnanensis. Therefore, to obtain high-quality raw materials for the development of the heavy building industry, the screening and purification of high-quality seed sources should be carried out in conjunction with relocation protection.

5. Conclusions

In summary, the genetic diversity and structure of 15 wild populations and 17 cultivated populations of P. polyphylla were analyzed by sequencing the cpDNA trnL-trnF fragments and comparing the results with those of AFLP markers. Based on the analysis of the trnL-trnF fragments, the wild populations of P. polyphylla had a moderate level of genetic diversity (HT = 0.861), which was consistent with the results of similar analyses using AFLP markers. A total of 15 haplotypes were identified in the 32 populations of Paris polyphylla; five unique haplotypes were identified in the cultivated populations, while only one unique haplotype was identified in the wild populations. Moreover, the genetic diversity of the cultivated populations was higher than that of the wild populations (HT = 0.900 vs. 0.861), which was also consistent with the results of the AFLP marker analysis. Combining the results shown in the UPGMA tree of the cultivated population grown in Yunnan province established based on AFLP markers with the results of the Mantel test of the cultivation group based on trnL-trnF fragment sequences, it was speculated that the unique haplotype of the cultivation group might have been obtained due to the mixed provenance of the cultivation populations.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes14091754/s1, Table S1: Collecting information of 15 wild populations and 17 cultivated populations of P. polyphylla var. yunnanensis.

Author Contributions

Conceptualization, D.W., Y.H. and L.R., methodology, Y.H. and L.R.; data curation and statistical analysis, H.D., M.M. and N.Z.; investigation, D.W. and Y.H.; writing—original draft preparation, D.W., Y.H., L.R. and H.D.; writing—review and editing, H.D., J.Q., M.M. and N.Z.; supervision, N.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Chongqing Science and Technology Commission (cstc2018jcyjAX0267), and the Scientific and Technological Research Program of Chongqing Municipal Education Commission (KJQN202201205).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data supporting the findings of this study are included in this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. IUCN. IUCN Red List Categories and Criteria v. 3.1; IUCN Species Survival Commission: Gland, Switzerland; Cambridge, UK, 2001; p. 35. [Google Scholar]
  2. Jin, L.Y.; Lu, T.X.; Qin, X.J.; Ni, W.; Yan, H.; Chen, Y.; Liu, H.Y. Two new highly oxygenated spirostanol saponins from Paris polyphylla var. stenophylla. Nat. Prod. Bioprospecting 2016, 6, 205–210. [Google Scholar] [CrossRef] [PubMed]
  3. Yu, M.J.; Xiang, H.Y.; Zhang, S.S.; Huang, Z.J.; Liu, Y. Contrastive Study on content and anti-oxidation activity of total saponins from cultivated and wild paris plants in Southwest China. Med. Plant 2017, 8, 7–11. [Google Scholar]
  4. Qi, J.J.; Zheng, N.; Zhang, B.; Sun, P.; Hu, S.N.; Xu, W.J.; Li, X.E. Mining genes involved in the stratification of Paris polyphylla seeds using high-throughput embryo transcriptome sequencing. BMC Genom. 2013, 14, 358. [Google Scholar] [CrossRef] [PubMed]
  5. He, J.; Zhang, S.; Wang, H.; Chen, C.X.; Chen, S.F. Advances in studies on and uses of Paris polyphylla var. yunnanensis (Trilliaceae). Acta Bot. Yunnanica 2006, 28, 271–276. [Google Scholar]
  6. He, J.; Wang, H.; Li, D.Z.; Chen, S.F. Genetic diversity of Paris polyphylla var. yunnanensis, a traditional Chinese medicinal Herb, Detected by ISSR Markers. Planta Medica 2007, 73, 1316–1321. [Google Scholar] [CrossRef]
  7. Song, Y.; Li, M.F.; Xu, J.; Zhao, Z.; Chen, N.Z. Polymorphic microsatellite markers in the traditional Chinese medicinal plant Paris polyphylla var. yunnanensis. Genet. Mol. Res. GMR 2015, 14, 9939–9942. [Google Scholar] [CrossRef]
  8. Chen, Z.; Tian, B.; Cai, C.T. Genetic diversity of Paris polyphylla var. yunnanensis by SSR marker. Chin. Tradit. Herb. Drugs 2017, 48, 1834–1838. [Google Scholar]
  9. Li, Y.B.; Wang, L.; He, Z.J.; Xie, S.Q.; Liao, Q.H.; Yang, S.C.; Liu, T. AFLP Analysis on Genetic Diversity of P. polyphylla var. yunnanensis. Bot. Res. 2018, 7, 45–53. [Google Scholar]
  10. Huang, Y.; Zhou, N.; Yang, M.; Shen, Y.X.; Zhang, D.Q. A comparative study of the population genetics of wild and cultivated populations of Paris polyphylla var. yunnanensis based on amplified fragment length polymorphism markers. Ecol. Evol. 2019, 9, 10707–10722. [Google Scholar] [CrossRef]
  11. Palme, A.E. Chloroplast DNA variation, postglacial recolonization and hybridization in hazel, Corylus avellana. Mol. Ecol. 2002, 11, 1769–1779. [Google Scholar] [CrossRef]
  12. Petit, R.J.; Hampe, A.; Cheddadi, R. Climate changes and tree phylogeography in the Mediterranean. Taxon 2005, 54, 877–885. [Google Scholar] [CrossRef]
  13. Williams, J.G.; Kubelik, A.R.; Livak, K.J.; Rafalski, J.A.; Tingey, S.V. DNA polymorphisms amplified by arbitrary primers are useful as genetic markers. Nucleic Acids Res. 1990, 18, 6531–6535. [Google Scholar] [CrossRef]
  14. Kumar, S.; Stecher, G.; Li, M.; Stecher, G.; Li, M.; Knyaz, G.; Tamura, K. MEGA X: Molecular evoluti-onary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  15. Watterson, G.A. On the number of segregating sites in genetical models without recombination. Theor. Popul. Biol. 1975, 7, 256–276. [Google Scholar] [CrossRef]
  16. Nei, M.; Li, W.H. Mathematical model for studying genetic variation in terms of restriction endonucleases. Proc. Natl. Acad. Sci. USA 1979, 76, 5269–5273. [Google Scholar] [CrossRef]
  17. Pons, O.; Petit, R.J. Measuring and testing genetic differentiation with ordered versus unordered alleles. Genetics 1996, 144, 1237–1245. [Google Scholar] [CrossRef] [PubMed]
  18. Excoffier, L.; Lischer, H.E. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 2010, 10, 564–567. [Google Scholar] [CrossRef]
  19. Mantel, N. The detection of disease clustering and a generalized regression approach. Cancer Res. 1967, 27, 209–220. [Google Scholar] [PubMed]
  20. Peakall, R.; Smouse, P.E. GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics 2012, 28, 2537–2539. [Google Scholar] [CrossRef]
  21. Zhao, J.J.; Huang, Y.; Zhang, D.Q.; Zhou, N. Phylogeography of Paris poliphylla var. yunnanensis based on chloroplast gene trnL-trnF sequences. China J. Chin. Mater. Medica 2021, 46, 1094–1101. [Google Scholar]
  22. Hamrick, J.L.; Godt, M.J.W.; Sherman-Broyles, S.L. Factors Influencing Levels of Genetic Diversity in Woody Plant Species. New For. 1992, 6, 95–124. [Google Scholar] [CrossRef]
  23. Cavers, S.; Navarro, C.; Lowe, A.J. A combination of molecular markers identifies evolutionarily significant units in Cedrela odorata L. (Meliaceae) in Costa Rica. Conserv. Genet. 2003, 4, 571–580. [Google Scholar] [CrossRef]
  24. Zhang, Q.; Chiang, T.; George, M.; Liu, J.; Abbott, R. Phylogeography of the Qinghai-Tibetan Plateau endemic Juniperus przewalskii (Cupressaceae) inferred from chloroplast DNA sequence variation. Mol. Ecol. 2010, 14, 3513–3524. [Google Scholar] [CrossRef] [PubMed]
  25. Huang, H.; Han, X.; Kang, L.; Raven, P.; Jackson, P.W.; Chen, Y. Conserving native plants in China. Science 2002, 297, 935–936. [Google Scholar] [CrossRef] [PubMed]
  26. Petit, R.J.; Duminil, J.; Fineschi, S.; Hampe, A.; Salvini, D.; Vendramin, G.G. Comparative organization of chloroplast, mitochondrial and nuclear diversity in plant populations. Mol. Ecol. 2005, 14, 689–701. [Google Scholar] [CrossRef] [PubMed]
  27. Huang, S.; Chiang, Y.C.; Schaal, B.A.; Chou, C.H.; Chiang, T.Y. Organelle DNA phylogeography of Cycas taitungensis, a relict species in Taiwan. Mol. Ecol. 2001, 10, 2669–2681. [Google Scholar] [CrossRef] [PubMed]
  28. Lu, S.Y.; Peng, C.I.; Cheng, Y.P.; Hong, K.H.; Chiang, T.Y. Chloroplast DNA phylogeography of Cunninghamia konishii (Cupressaceae), an endemic conifer of Taiwan. Genome 2001, 44, 797–807. [Google Scholar] [CrossRef] [PubMed]
  29. Yuan, Q.J.; Zhang, Z.Y.; Hu, J.; Guo, L.P.; Shao, A.J.; Huang, L.Q. Impacts of recent cultivation on genetic diversity pattern of a medicinal plant, Scutellaria baicalensis (Lamiaceae). BMC Genet. 2010, 11, 29. [Google Scholar] [CrossRef] [PubMed]
  30. Hamilton, A. Medicinal plants, conservation and livelihoods. Biodivers. Conserv. 2004, 13, 1477–1517. [Google Scholar] [CrossRef]
  31. Mullen, S.P. Conservation and the Genetics of Populations; Blackwell Publishing: Hoboken, NJ, USA, 2007. [Google Scholar]
  32. He, J.S.; Chen, L.; Si, Y.; Huang, B.; Ban, X.Q.; Wang, Y.W. Population structure and genetic diversity distribution in wild and cultivated populations of the traditional Chinese medicinal plant Magnolia officinalis subsp. biloba (Magnoliaceae). Genetica 2009, 135, 233. [Google Scholar] [CrossRef]
Figure 1. Mantel test analysis: (a) Mantel test of wild populations and (b) Mantel test of cultivated populations.
Figure 1. Mantel test analysis: (a) Mantel test of wild populations and (b) Mantel test of cultivated populations.
Genes 14 01754 g001
Figure 2. Cluster analysis diagram of the wild (a) and cultivated (b) populations based on SAMOVA analysis. Note: FCT means the intergroup genetic differentiation value.
Figure 2. Cluster analysis diagram of the wild (a) and cultivated (b) populations based on SAMOVA analysis. Note: FCT means the intergroup genetic differentiation value.
Genes 14 01754 g002
Figure 3. The relative proportions of different chloroplast haplotypes found in the wild and cultivated P. polyphylla var. yunnanensis populations. H1–H15: haplotype 1–haplotype 15; H in bold refers to the unique haplotypes of the wild and cultivated P. polyphylla var. yunnanensis populations. The data in parentheses indicate the number of plants.
Figure 3. The relative proportions of different chloroplast haplotypes found in the wild and cultivated P. polyphylla var. yunnanensis populations. H1–H15: haplotype 1–haplotype 15; H in bold refers to the unique haplotypes of the wild and cultivated P. polyphylla var. yunnanensis populations. The data in parentheses indicate the number of plants.
Genes 14 01754 g003
Table 1. Genetic diversity and genetic differentiation of P. polyphylla var. yunnanensis based on cpDNA (trnL-trnF).
Table 1. Genetic diversity and genetic differentiation of P. polyphylla var. yunnanensis based on cpDNA (trnL-trnF).
ParametersWild PopulationsCultivated PopulationsTotal Populations
Haplotype101415
HS0.135 (0.0489)0.222 (0.0557)0.181 (0.0376)
HT0.861 (0.0479)0.900 (0.0298)0.866 (0.0269)
GST0.843 (0.0558)0.754 (0.0618)0.791 (0.0422)
NST0.919 (0.0450)0.815 (0.0627)0.860 (0.0399)
Notes: HS, haplotype diversity at the population level; HT, haplotype diversity at the species level; GST and NST, genetic differentiation obtained using different calculation methods. The data in parentheses indicate nucleotide diversity π.
Table 2. Genetic differentiation coefficients (NST and GST) for the cpDNA of P. polyphylla var. yunnanensis determined by conducting 1000 substitution U-tests.
Table 2. Genetic differentiation coefficients (NST and GST) for the cpDNA of P. polyphylla var. yunnanensis determined by conducting 1000 substitution U-tests.
PermutationsWild Populationsp-ValueCultivated Populationsp-Value
NSTGSTNSTGST
9500.9020.843<0.050.8210.754>0.05
9900.9130.843<0.010.8320.754>0.01
Notes: GST and NST: genetic differentiation coefficients obtained using different calculation methods.
Table 3. Analysis of molecular variance (AMOVA) of the 32 populations of P. polyphylla var. yunnanensis based on the cpDNA sequences.
Table 3. Analysis of molecular variance (AMOVA) of the 32 populations of P. polyphylla var. yunnanensis based on the cpDNA sequences.
Source of Varianced.f.SSDVCPV (%)Fixation Indicesp-Value
All populations without hierarchy
Among populations31122.2420.33976.660FST = 0.767<0.001
Within populations33234.3040.10323.340
Total336156.5470.443
Wild populations vs. cultivated populations
Among groups12.3830.0112.500FCT = −0.025>0.05
Among populations within groups30119.8590.34478.86FSC = 0.769<0.001
Within populations33234.3040.10320.64FST = 0.764<0.001
Total363156.5470.459
Notes: FST represents the inter-population genetic differentiation; FCT represents the intergroup genetic differentiation.
Table 4. Sample and geographical information of the 15 wild populations and 17 cultivated populations of P. polyphylla var. yunnanensis.
Table 4. Sample and geographical information of the 15 wild populations and 17 cultivated populations of P. polyphylla var. yunnanensis.
Population IDNHaplotypesLocalityLatitude (N)/Longitude (E)Altitude (m)Sample No.
Wild populations
W-LL6H1 (6)Longli County, Qiannan Prefecture, Guizhou Province26°27′47″, 106°59′33″10802017CL-003
W-XY4H1 (4)Qishe Town, Xingyi City, Guizhou Province25°00′40″, 104°49′10″17532017CL-041
W-HD13H3 (2), H4 (11)Laojuntan Township, Huidong County, Sichuan Province26°23′32″, 102°57′55″2205HY16080401
W-YM7H5 (7)Yimen County, Yuxi City, Yunnan Province24°58′38″, 102°12′52″1893HY16080901
W-CX10H5 (10)Lvhe Town, Chuxiong City, Yunnan Province25°07′55″, 101°22′30″1871HY16081401
W-LY11H6 (11)Puyao Town, Longyang District, Baoshan City, Yunnan Province25°02′09″, 099°04′108″2298HY16081801
W-CN9H9 (9)Changning County, Baoshan City, Yunnan Province24°94′20″, 099°56′52″2074HY16081901
W-YP10H6 (6), H8 (4)Yongping County, Dali Bai Autonomous Prefecture, Yunnan Province25°21′27″, 099°23′14″1925HY16082001
W-YL11H8 (11)Yunlong County, Dali Bai Autonomous Prefecture, Yunnan Province25°34′57″, 099°07′29″2268HY16082101
W-LJ7H6 (6), H11 (1)Yulong Naxi Autonomous County, Lijiang City, Yunnan Province27°01′94″, 100°22′01″3200HY16082201
W-LG13H1 (9), H13 (4)Longgong Town, Xixiu District, Anshun City, Guizhou Province26°05′42″, 105°52′43″1178HY16071301
W-XX12H1 (11), H13 (1)Xixiu District, Anshun City, Guizhou Province26°15′50″, 106°00′35″1410HY16071501
W-QZ10H1 (10)Weicheng Town, Qingzhen City, Guizhou Province26°44′43″, 106°22′57″1363HY16071601
W-XR12H1 (2), H2 (10)Xingren County, Xingyi City, Guizhou Province25°32′46″, 105°27′35″1515HY16072001
W-GD6H5 (6)Guandu District, Kunming City, Yunnan Province24°59′19″, 102°58′43″2308HY16080701
Cultivated populations
C-LL20H1 (20)Longli County, Qiannan Prefecture, Guizhou Province26°27′47″, 106°59′33″1080HY16071703
C-XY5H1 (4), H2 (1)Qishe Town, Xingyi City, Guizhou Province25°00′40″, 104°49′10″1753HY16072103
C-HD14H3 (4), H4 (10)Laojuntan Township, Huidong County, Sichuan Province26°23′43″, 102°58′06″2073HY16080301
C-YM12H5 (12)Yimen County, Yuxi City, Yunnan Province24°58′38″, 102°12′52″1893HY16080903
C-CX8H5 (8)Lvhe Town, Chuxiong City, Yunnan Province25°07′55″, 101°22′30″1871HY16081403
C-LY20H7 (3), H8 (17)Puyao Town, Longyang District, Baoshan City, Yunnan Province25°02′09″, 099°04′08″2298HY16081601
C-CN20H6 (5), H9 (15)Changning County, Baoshan City, Yunnan Province24°94′20″, 099°56′52″2074HY16081903
C-YP15H6 (13), H8(3)Yongping County, Dali Bai Autonomous Prefecture, Yunnan Province25°21′08″, 099°23′06″1829HY16082003
C-YL12H8 (11), H10 (1)Yunlong County, Dali Bai Autonomous Prefecture, Yunnan Province25°34′57″, 099°07′29″2172HY16082005
C-LJ15H6 (8), H10(1), H12 (2), H13 (4)Yulong Naxi Autonomous County, Lijiang City, Yunnan Province27°01′94″, 100°22′01″3200HY16082203
C-ZJ5H1 (5)Zhijin County, Anshun City, Guizhou Province26°48′33″, 105°38′57″1369HY16071201
C-ZY5H1 (5)Ziyun County, Anshun City, Guizhou Province25°59′01″, 106°04′53″1148HY16071401
C-SM14H3 (1), H14(13)Yanglin Town, Songming County, Kunming City, Yunnan Province25°09′25″, 103°02′31″1996HY16080801
C-BC15H6 (10), H8(1), H9(1), H15(3)Paicai Village, Longyang District, Baoshan City, Yunnan Province25°11′50″, 099°19′15″2321HY16081603
C-LX20H6 (20)Longxin Township, Longling County, Baoshan City, Yunnan Province24°32′12″, 098°46′54″1842HY16081701
C-MS9H6 (3), H9 (6)Mang City, Dehong Prefecture, Yunnan Province24°29′24″, 098°20′11″1936HY16081703
C-JC14H10 (14)Yangcen Township, Jianchuan County, Dali Bai Autonomous Prefecture, Yunnan Province26°48′42″, 099°80′75″2944HY16082301
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, D.; Huang, Y.; Rui, L.; Du, H.; Qi, J.; Ma, M.; Zhou, N. Population Genetic Analysis of Paris polyphylla var. yunnanensis Based on cpDNA Fragments. Genes 2023, 14, 1754. https://doi.org/10.3390/genes14091754

AMA Style

Wang D, Huang Y, Rui L, Du H, Qi J, Ma M, Zhou N. Population Genetic Analysis of Paris polyphylla var. yunnanensis Based on cpDNA Fragments. Genes. 2023; 14(9):1754. https://doi.org/10.3390/genes14091754

Chicago/Turabian Style

Wang, Dan, Yu Huang, Lu Rui, Huihui Du, Junsheng Qi, Mingguo Ma, and Nong Zhou. 2023. "Population Genetic Analysis of Paris polyphylla var. yunnanensis Based on cpDNA Fragments" Genes 14, no. 9: 1754. https://doi.org/10.3390/genes14091754

APA Style

Wang, D., Huang, Y., Rui, L., Du, H., Qi, J., Ma, M., & Zhou, N. (2023). Population Genetic Analysis of Paris polyphylla var. yunnanensis Based on cpDNA Fragments. Genes, 14(9), 1754. https://doi.org/10.3390/genes14091754

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop