Simplified map showing the geographic distribution of the 28 sampling sites of Cymbidium tortisepalum in Yunnan Province, China. Blue marks represent populations in the northern region of three parallel rivers in the western Yunnan Province (NW); green marks markers represent those in the middle region of three parallel rivers in the western Yunnan Province (MD), and red marks represent those in the southern region of three parallel rivers in the western Yunnan Province. The abbreviations of the sites refer to Table 1.
Fig. 2.
Geographic distribution of the chloroplast DNA haplotypes of wild Cymbidium tortisepalum in Yunnan Province, China. Different colors in the label box represent different haplotypes, and the pie charts represent the distribution proportion of each haplotype in the populations.
Fig. 3.
Neighbor-joining trees obtained by analysis of 11 chloroplast DNA haplotypes of wild Cymbidium tortisepalum in Yunnan Province, China. Hap1–Hap11 represent 11 haplotypes, respectively. Bootstrap values (>50%) are based on 1000 replicates and are indicated above branches.
Fig. 4.
The haplotype network of 11 chloroplast DNA haplotypes detected in 28 populations of wild Cymbidium tortisepalum in Yunnan Province, China. The size of circles corresponds to the frequency of each haplotype, H1–H11, represent 11 haplotypes, respectively, with each solid line representing one step mutation.
Fig. 5.
Mismatch distribution analysis based on chloroplast DNA of wild Cymbidium tortisepalum in Yunnan Province (A), in northern region of three parallel rivers in the western Yunnan Province (B), in middle region of three parallel rivers in the western Yunnan Province (C), and in southern region of three parallel rivers in the western Yunnan Province (D). Solid lines indicate expected value (Exp), dotted lines reflect observed value (Obs).
Genetic Diversity and Population Structure Analysis of Wild Cymbidium tortisepalum Based on Chloroplast DNA in Yunnan Province of China
Authors:
,
,
, and
Click on author name to view affiliation information
Cymbidium tortisepalum is a primary orchid species in Yunnan Province, China, and has an extremely high ornamental and economic value. To reveal the levels and distribution of genetic variation and structure of wild C. tortisepalum resources, sequence variations of six chloroplast DNA intergenic spacers (psbM-trnD, trnV-trnA, accD-psal, rrn23, trnk-rps16, and ycf1) were analyzed in 404 wild individuals from 28 populations in the three river area in Yunnan Province, China. The results showed that the six chloroplast DNA sequences were aligned with 61 polymorphic sites, including 50 indels and 11 haplotypes in 404 individuals, which revealed a low level of genetic diversity (total genetic diversity = 0.240, and the average value of nucleotide diversity = 0.00024). In addition, a fairly low genetic differentiation [coefficients for genetic differentiation among populations (GST) = 0.099, number of substitution (NST) = 0.081] was found among the studied populations, and NST value was less than GST, which indicated that no significant phylogeographic structure existed in those populations. Furthermore, analysis of molecular variance revealed that great genetic variance (91%) came from individuals within the populations, which indicated that there was no clear genetic differentiation among populations. On the basis of these findings, a conservation plan was proposed to sample or preserve fewer populations but with more individuals from each population.
Cymbidium tortisepalum is a perennial herb, which belongs to Orchidaceae family (Du Puy and Cribb, 1988). Wild C. tortisepalum mainly grows in the river valleys or under secondary forest of Pinus yunnanensis, Abies fabri, and Quercus aquifolioides in the middle of the mountains, with a range of 1000 to 2700 m, 70% to 80% shading intensities, and 60% to 80% air humidity, and where the soil is moist all year with a high nutrient content and pH between 5.5 and 6.5 (Xue et al., 2007). Because this plant has an extremely narrow geographic spread, high dependence on ecological environment, and poor natural reproduction ability, its wild resources are rare, being scattered only within Taiwan Province, the southwest of Sichuan Province, and the northwest of Yunnan Province, China (Lin et al., 2016; Peng, 2003). It is primarily distributed in the three parallel river region (Jinsha, Lancang, and Nujiang rivers) of the northwestern Yunnan Province, at latitude 24°N to 28.4°N and longitude 98°E to 101°E. C. tortisepalum is famous for its graceful shape, varied flowers, and beautiful color, which gives it a high ornamental and economic value, and it has become one of the most sought-after orchid species in Yunnan Province (Du Puy and Cribb, 1988). The C. tortisepalum industry has reached a value of more than 100 million yuan in the Yunnan (Qi, 2004). With increased market, the price of one wild C. tortisepalum can reach up to tens of thousands and even millions of yuan. Consequently, the wild resources of C. tortisepalum in Yunnan Province are being diminished by overexploitation and habitat loss (Rao et al., 2008; Xue et al., 2007). Therefore, the research on germplasm resources, and the evaluation of the genetic diversity and structure of wild C. tortisepalum in Yunnan Province, could provide valuable knowledge to the protection of wild resources. However, current studies of C. tortisepalum, both within China and abroad, are mainly focused on cultivation management (He et al., 2019), tissue culture, and rapid propagation (Huan et and Tanake, 2004; Teixeira da Silva et al., 2006), and cross-breeding (Li et al., 2009, 2010; Zhao et al., 2019). Reports of genetic diversity and genetic structure of wild C. tortisepalum are limited to morphological diversity (Suo et al., 2016; Ye et al., 2015), nuclear simple sequence repeats [nSSR (Zhao et al., 2017)] and direct amplification of length polymorphism [DALP (Jia et al., 2012)]. Moreover, population samples in other studies represented less than 10% of the wild C. tortisepalum found in Yunnan Province and were lower than the sample size in this study. Furthermore, there is no report about the genetic diversity and structure of wild C. tortisepalum base on chloroplast DNA (cpDNA). Therefore, it is highly important to further explore the genetic diversity of populations and the genetic structure of wild C. tortisepalum based on cpDNA in Yunnan Province to adopt appropriate references for the formation of effective conservation strategies.
In our study, 404 wild individuals from 28 populations in Yunnan Province were used as experimental materials, and their genetic structure and their lineage differentiation were investigated by cpDNA sequence analysis. The purpose of our study is to provide a theoretical basis for the formation of a protection strategy and to provide technical support for the industrial development of wild C. tortisepalum in Yunnan Province of China.
Methods
Experimental material.
On the basis of two investigations of the wild C. tortisepalum resources in Yunnan Province before sampling, 28 sampling sites located in the western area of the Yunnan Province were determined. Approximately 15 to 24 individuals were selected from each population with a distance of at least 20 m. Finally, 404 individuals from 28 populations of wild C. tortisepalum were collected between 2010 and 2011. Afterward, 1 g of leaves were sampled from each plant, sealed in a bag, and placed in a sealing box with desiccant. Prepared materials were stored at –80 °C in the Southwest China Germplasm Bank of Wild Species (Kunming, China). The collection area name, latitude, longitude, altitude, and habitat were recorded for each sample. Detailed locations and geographic distribution of the populations are shown in Table 1 and Fig. 1.
Fig. 1.Simplified map showing the geographic distribution of the 28 sampling sites of Cymbidium tortisepalum in Yunnan Province, China. Blue marks represent populations in the northern region of three parallel rivers in the western Yunnan Province (NW); green marks markers represent those in the middle region of three parallel rivers in the western Yunnan Province (MD), and red marks represent those in the southern region of three parallel rivers in the western Yunnan Province. The abbreviations of the sites refer to Table 1.
Citation: Journal of the American Society for Horticultural Science 146, 6; 10.21273/JASHS05046-21
Table 1.Details of sampling sites of 28 populations of wild Cymbidium tortisepalum in Yunnan Province, China. From 2010 to 2011, we collected 404 individuals from 28 naturally distributed populations of C. tortisepalum in northwestern Yunnan. Approximately 15 to 24 individuals were selected from each population with a distance at least 20 m.
Total DNA extraction.
Genomic DNA was extracted from a single individual using the cetyltrimethylammonium bromide (CTAB) method (Doyle, 1986) with some modifications (Tang, 2014). The quality and quantity of the DNA were checked on 1% agarose gels. The extracted DNA was diluted to 50 ng·μL−1 and stored at –20 °C for future use.
Polymerase chain reaction amplification and sequencing.
On the basis of complete chloroplast genome sequence of the genus C. tortisepalum by our research group (Yang et al., 2013), we selected 43 regions for designing polymerase chain reaction (PCR) primers. Of these, six chloroplast fragments had variable sites (psbM-trnD, trnV-trnA, accD-psal, rrn23, trnk-rps16, and ycf1) among 10 tested samples in the pilot study, which were finally used to sequence all samples. Sequences were deposited at GenBank with accession numbers KY621375-KY621467.
PCR amplification was performed using a thermal cycler (GeneAmp PCR System 9700; PerkinElmer, Waltham, MA). Primer information is shown in Table 2. PCR reactions were 25 μL in volume and contained 2.5 μL 10 × PCR buffer, 2.5 μL MgCl2 (25 mm), 2.0 μL dNTP mixture (2.5 mm), 1.0 μL primer F (5 μM), 1.0 μL primer R (5 μM) DNA, 0.125 μL Taq polymerase (5 U.μL−1), 2 μL Template DNA (50 ng·μL−1), and then adding double-distilled water (ddH2O) to a total volume of 25 μL. The cycling parameters were as follows: 3 min of denaturation at 94 °C; followed by 35 cycles of 1 min denaturation at 94 °C, 1.5 min annealing at 56 °C, and 2 min extension at 72 °C; with a final extension of 72 °C for 10 min. PCR products were visualized using gel electrophoresis on an agarose gel (0.8% to 1%). Sequencing of PCR products was accomplished using a DNA analyzer (ABI 3730; Thermo Fisher Scientific, Waltham, MA).
Table 2.Primer sequences for the six genetic loci of wild Cymbidium tortisepalum in Yunnan Province, China. In our previous study (Yang et al., 2013), 43 primer sequences were designed based on complete chloroplast genome sequence and six chloroplast fragments had variable sites among 10 tested samples in the pilot study, which were finally used to sequence all samples. The sequence has been deposited at GenBank with accession numbers KY621375-KY621467.
Data analysis.
Sequencher version 4.1.4 software (Larkin et al., 2014) was used for sequence splicing. Geneious version 4.8 (Drummond et al., 2011) was used for sequence alignment and calibration. DnaSP version 5.0 software (Librado and Rozas, 2009) was used to perform haplotype statistics and sites with gaps were not considered in our analysis. ArcMap software (version 9.3; ESRI, Redlands, CA) was adopted to construct the geographic information map of all haplotypes within each population (Gmbh, 2004). Network version 4.6 software (Forster et al., 2007) was used to construct a haplotype network diagram. Arlequin version 3.5 (Excoffier and Lischer, 2010) was used to detect the genetic structure and perform analysis of molecular variance (AMOVA). GenALEx version 6.5 (Peakall and Smouse, 2006) was used to calculate the geographic distance, and the Mantel test (Mantel, 1967) was used to test the correlation between the two matrices. Genetic diversity (HS), the total genetic diversity (HT), and the genetic differentiation coefficients (GST and NST) were calculated using Permut version 2.0 software (Pons and Petit, 1996), GST and NST were compared via the U statistical method (1000 repeated permutation test) (Nei, 1987; Tajima, 1993). We used MEGA5 (Tamura et al., 2011) to construct a phylogenetic tree with neighbor-joining (NJ) method. Arlequin version 3.5 (Excoffier and Lischer, 2010) was used to perform mismatch analyze and calculate the Tajima’s D (Tajima, 1993, 1996) and Fu’s parameters (Fu, 1997) based on the infinite alleles model. DnaSP version 5.0 (Librado et al., 2009) was employed to construct mismatch distribution diagram. Our results suggested that although these methods and indices were commonly to infer population expansion, they vary in sensitivity and in their adaptability to data.
Results
Sequence characteristics and haplotype distribution.
A sequence with a total length of 2702 bp was obtained from the six chloroplast fragments of the 404 individuals of wild C. tortisepalum from the Yunnan Province. Sixty-one mutation sites were detected, which included 50 indels and 11 haplotypes based on 11 base substitution sites. The distribution of each haplotype with the population is shown in Table 3 and Fig. 2.
Fig. 2.Geographic distribution of the chloroplast DNA haplotypes of wild Cymbidium tortisepalum in Yunnan Province, China. Different colors in the label box represent different haplotypes, and the pie charts represent the distribution proportion of each haplotype in the populations.
Citation: Journal of the American Society for Horticultural Science 146, 6; 10.21273/JASHS05046-21
Table 3.Distribution frequencies of the 11 haplotypes in the 28 populations of wild Cymbidium tortisepalum in Yunnan Province, China.
Among the six sequences, the sequence length of the chloroplast trnk-rps16 fragment was the longest, which was 786 bp, and two base substitutions were detected. The sequence length of the chloroplast trnV-trnA fragment was the second longest at 704 bp, with two base substitutions. The sequence length of the chloroplast fragments psbM-trnD, accD-psal, rrn23, and ycfl, were 171, 283, 486, and 272 bp, respectively. There were only three, one, two, and one polymorphic loci detected in these fragments, respectively.
Nine of the 28 wild populations of C. tortisepalum sampled contained only one haplotype (32% of the population is within a single haplotype group). The remaining 19 populations showed polymorphisms, which included two to four haplotypes. Haplotypes H6, H7, H10, and H11 existed within only one individual, which were private haplotypes belonging to some of the populations. H6 and H7, two kinds of private haplotypes, both were distributed in the Luzhang (LZ) population. The haplotypes H10 and H11 were only found in the Pupiao (PP) population. Among the 11 haplotypes, H1 was the most widely distributed, and was present in all 28 populations. The H3 haplotype was the second most widely distributed.
Genetic diversity and genetic structure.
Results revealed that the genetic diversity of wild C. tortisepalum at the species level (HT = 0.204, π = 0.00024) and population level were low. In regard to genetic differentiation, there was no significant differentiation among populations (GST = 0.099, NST = 0.081). The higher GST value than NST value, indicates that the probability of the presence of haplotypes with close relation in the same or similar populations is low and there is no significant phylogeographic structure within the wild C. tortisepalum populations (Tables 4 and 5).
Table 4.Total genetic diversity of wild Cymbidium tortisepalum based on chloroplast DNA in Yunnan Province, China. Average genetic diversity within populations (HS), the total genetic diversity (HT), the genetic differentiation coefficients (GST) and number of substitution (NST) were calculated using Permut version 2.0 software (Pons and Petit, 1996), GST and NST were compared via the U statistical method [1000 repeated permutation test (Nei, 1987; Tajima, 1993)].
Table 5.Genetic diversity and neutral test of wild Cymbidium tortisepalum based on chloroplast DNA in Yunnan Province, China.
AMOVA revealed that the greatest genetic variance (≈91%) came from individuals within the population of wild C. tortisepalum, and the remaining 9.001% of the genetic variation existed among the populations, which indicated there is no clear genetic differentiation among populations (Table 6).
Table 6.Analyses of molecular variance of wild Cymbidium tortisepalum based on chloroplast DNA in Yunnan Province, China.
The results of the Mantel test analysis showed that the genetic distance of the wild C. tortisepalum had no significant correlation with geographical distance (r = 0.279, P = 0.020). The NJ adjacency tree constructed by the nonforeign group is shown in Fig. 3. Results revealed that 11 haplotypes of the wild C. tortisepalum formed six major branches with a 75% bootstrap value and six haplotypes (H2, H3, H6, H9, H10, and H11) constituted an independent branch. In addition, the remaining five haplotypes formed the collateral independent branch. The chloroplast haplotype network diagram constructed using the network of wild C. tortisepalum is shown in Fig. 4. All the haplotypes constituted a relatively simple linear relationship.
Fig. 3.Neighbor-joining trees obtained by analysis of 11 chloroplast DNA haplotypes of wild Cymbidium tortisepalum in Yunnan Province, China. Hap1–Hap11 represent 11 haplotypes, respectively. Bootstrap values (>50%) are based on 1000 replicates and are indicated above branches.
Citation: Journal of the American Society for Horticultural Science 146, 6; 10.21273/JASHS05046-21
Fig. 4.The haplotype network of 11 chloroplast DNA haplotypes detected in 28 populations of wild Cymbidium tortisepalum in Yunnan Province, China. The size of circles corresponds to the frequency of each haplotype, H1–H11, represent 11 haplotypes, respectively, with each solid line representing one step mutation.
Citation: Journal of the American Society for Horticultural Science 146, 6; 10.21273/JASHS05046-21
We examined the expansion history of wild C. tortisepalum using the tools of mismatch analysis and neutral test (Table 5, Fig. 5). According to the results of mismatch analysis, the Harpending’s raggedness (HRag) index of C. tortisepalum changed as a curve with a single peak both in the overall level or the three division areas and the observed values were in good agreement with the expected values, which indicated that the populations might have experienced rapid expansion, or it may have been caused by high mobility among neighboring groups accompanied with rapid population expansion. However, this inference was not supported by the neutral test because no significant differences were found among Tajima’s D value and Fu and Li’s value.
Fig. 5.Mismatch distribution analysis based on chloroplast DNA of wild Cymbidium tortisepalum in Yunnan Province (A), in northern region of three parallel rivers in the western Yunnan Province (B), in middle region of three parallel rivers in the western Yunnan Province (C), and in southern region of three parallel rivers in the western Yunnan Province (D). Solid lines indicate expected value (Exp), dotted lines reflect observed value (Obs).
Citation: Journal of the American Society for Horticultural Science 146, 6; 10.21273/JASHS05046-21
Population genetic diversity and genetic differentiation.
Chloroplast DNA is characterized by the small genome, maternal inheritance of angiosperm, and small selection pressure and rapid variation of its nonencoding area (Wolfe et al., 1989). Therefore, it is highly suitable for studies on natural population genetic structure, the system evolution of plants. This study reports the largest survey of molecular diversity of native C. tortisepalumto date, with 404 individiuals from 28 populations sequenced at six cpDNA fragments. The total genetic diversity of wild C. tortisepalum was low. The reason for this result may be firstly related to the properties of the chloroplast genome. The chloroplast genome is a haploidisation of the cytoplasmic genome, and its evolutionary rate is low. Compared with the nuclear genome, the effective population size reflected by the chloroplast genome is much smaller, which is about half of the nuclear genes (Birky, 1988). In addition, long-term human activities that exert pressure on the population cause a sharp reduction in the number of groups in the distribution of the population. According to our investigations in recent years, the distribution of wild C. tortisepalum in Dali, Nujiang, Baoshan, and other places has been in declining state due to overutilization of the species. The effective population has been severely reduced and the random genetic drift intensified, which has caused the loss of genetic diversity (Tang, 2014). Furthermore, as the effective size of the populations decrease, the effects of random genetic drift and the bottleneck effect on chloroplast genome are more severe than those on the nuclear genome. The loss of genetic diversity, reflected by the chloroplast genome, is also amplified (Jordano and Godoy, 2000; Xun et al., 2011).
Alternatively, the results of this study may be related to the narrow geographic distribution and breeding system of this species. Wild resources of C. tortisepalum are mainly concentrated in the Dali Bai Autonomous Region, the Nujiang Lisu Nationality Autonomous Prefecture, Baoshan City, Lijiang City, and the Diqing Tibetan Autonomous Prefecture in the Yunnan Province. In general, species with a small distribution tend to have lower genetic diversity than those with a wide range distribution.
AMOVA of the cpDNA sequence results showed that the genetic differentiation of wild C. tortisepalum primarily existed within the populations (91%), and genetic differentiation among populations (FSC = 8.4%) and divergence among all populations (FST = 0.09001) were relatively low, which was consistent with the theory that 90% of variation exists within populations of cross-pollinated plants (Hamrickn and Godt,1990). A similar result of 89.25% molecular variation within populations was obtained in our previous investigation involving AMOVA of SSR data of these plants (Tang, 2014; Zhao et al., 2017), which, combined with the present results, indicate that the genetic background of wild C. tortisepalum is relatively simple, and there are no complex lineage differentiations. Low genetic differentiation of wild C. tortisepalum may be due to its unique life history, its highly specialized floral traits, its long evolutionary history, and severe bottlenecks, which the species may have experienced during the last ice age, as well as the fragmentation of the population distribution region with effects on seed and pollen flow.
The nucleotide polymorphism (Hd= 0.241) and haplotype polymorphism (π = 0.24 × 102) of C. tortisepalum at the species level were low. The genetic diversity of this species was primarily concentrated within populations (HS = 0.216), which was consistent with the results of a study by Jia et al., (2012). At the species level, the NST value was less than the GST value (P > 0.05), which indicated that the probability of closely related haplotypes in the same or similar populations was low and that- there was no clear spectrum structure of the 28 populations in this study. This result was consistent with that of our previous study on SSR (Tang, 2014), which was in agreement with the theory that there is no significant correlation between geographic distance and genetic distance among outcrossing species, especially species with a small geographic distribution (Willyard et al., 2007).
The level of genetic diversity within a species depends on the breeding system, especially the number of effective populations and their variation among generations, which has significance to the future evolution and survival of the species. Endemic, rare, or endangered species are generally considered to have low levels of genetic diversity, due to low population numbers, the isolation of populations, and their adaptation to a single habitat (Barrett and Kohn, 1991). The genetic diversity of the endangered and their nonendangered species were compared by Frankham (1995) and Spielman et al. (2004), which demonstrated that the genetic diversity of most endangered species was lower than that of corresponding nonendangered species. However, results of another study revealed that the genetic diversity level was still high for some endangered species, such as Iris aphylla (Wroblewska et al., 2003), Suzukia luchuensis (Maki et al., 2003), and Sarracenia leucophylla (Wang et al., 2004).
Chloroplast haplotype analysis.
When interpreting the haplotype network diagram, the newly derived haplotype is often located at the tip of the clade, while the ancestral haplotypes are usually found in interior node (Crandall and Templeton, 1993; Schaal et al., 1998; Templeton, 1998). The most widely distributed haplotype, H1, located in the central position, was the ancestral haplotype, and H5, H8, H4, and the unique haplotypes H6 and H7, located at the end of the network, were the newly derived haplotypes.
Moreover, the haplotype diversity of the Bady (BD) population located in Badi Township in Diqing Weixi County was the highest. The number of haplotypes (four) was the highest in the Changxin (CX) population located in Dali and the LZ population located in the NuJiang region. The private haplotypes H10 and H11 found in PP and Taiping (TP) populations located in Pupiao and Taiping town of Baoshan Longyang District, respectively. Populations in a refuge should have a higher genetic diversity than those in a nonrefuge, and they should contain the ancestral (or ancient) haplotype (Hewitt, 1999, 2004; Lascoux et al., 2004; Petit et al., 2005). It can be inferred that BD, CX, Lushui (LS), and the towns of PP and TP in Baoshan may have been a refuge for wild C. tortisepalum during the last ice age, and these areas may have richer resources of wild C. tortisepalum.
Conservation and future studies on wild C. tortisepalum.
Orchidaceae is one of the largest families of flowering plants, although most are rare and endangered plants. All wild orchid species are included in Appendix II of the Convention on International Trade in Endangered Species of Wild Fauna and Flora. The main goal for protecting an endangered species is to conserve its genetic diversity because the species may only survive and adapt to the environmental changes by maintaining sufficient genetic variation within and among populations (Barrett and Kohn, 1991). The results of this study showed that the total genetic diversity and genetic deviation of wild C. tortisepalum, as analyzed by cpDNA, were low. Therefore, this species may have a relatively low tolerance to ecologic or random climatic events and natural phenomena. On the basis of our investigation, we considered that the main cause in impending danger to wild C. tortisepalum resources was not due to its genetic basis and biological characteristics but to habitat fragmentation and the rapid decrease in the number of populations and the number of individuals within populations, which is caused by habitat destruction and overutilization by humans. We found that the survey region for wild C. tortisepalum had experienced severe habitat destruction, which made it difficult to find wild C. tortisepalum resources in its reported distribution areas. Most of the populations were small, and the individuals in the population were scattered. With time, the loss of genetic diversity would further exacerbate the endangered status of this species. Therefore, to maintain the long-term survival of this species, the protection of the habitat of wild C. tortisepalum is critical. As possible refuges for resources of wild C. tortisepalum, BD, CX, LS, PP, and TP towns should be given priority for protection and future studies.
The results of this study also showed that the vast majority of genetic variation in wild C. tortisepalum was located within populations. The haplotype frequency had a large difference within populations, and low frequency haplotypes accounted for the majority of the populations. The low frequency of haplotypes was likely caused by decreases in the populations and random genetic drift and also the result of declining genetic diversity within populations. Therefore, we suggest that the number of populations could be reduced and that individuals from each population should be increased for conservation or sampling. In particular, this protection strategy should be preferred in populations with high genetic diversity and low frequency haplotype groups.
Future studies should be based on the investigation of large-scale resources, with the collection of as many representative population samples as possible from areas with less human activity and more simple ecological climates, such as the Diqing Tibetan Autonomous Region. The diversity and genetic differentiation of wild C. tortisepalum should be studied by the method of single-copy nuclear gene sequence and cpDNA sequences of a certain length, which would provide a more reliable basis for further elucidation of the genetic diversity and the population genetic differentiation patterns of wild C. tortisepalum.
Received: 20 Jan 2021
Accepted: 11 Aug 2021
Published online: 30 Sept 2021
Published print: 01 Nov 2021
Fig. 1.
Simplified map showing the geographic distribution of the 28 sampling sites of Cymbidium tortisepalum in Yunnan Province, China. Blue marks represent populations in the northern region of three parallel rivers in the western Yunnan Province (NW); green marks markers represent those in the middle region of three parallel rivers in the western Yunnan Province (MD), and red marks represent those in the southern region of three parallel rivers in the western Yunnan Province. The abbreviations of the sites refer to Table 1.
Fig. 2.
Geographic distribution of the chloroplast DNA haplotypes of wild Cymbidium tortisepalum in Yunnan Province, China. Different colors in the label box represent different haplotypes, and the pie charts represent the distribution proportion of each haplotype in the populations.
Fig. 3.
Neighbor-joining trees obtained by analysis of 11 chloroplast DNA haplotypes of wild Cymbidium tortisepalum in Yunnan Province, China. Hap1–Hap11 represent 11 haplotypes, respectively. Bootstrap values (>50%) are based on 1000 replicates and are indicated above branches.
Fig. 4.
The haplotype network of 11 chloroplast DNA haplotypes detected in 28 populations of wild Cymbidium tortisepalum in Yunnan Province, China. The size of circles corresponds to the frequency of each haplotype, H1–H11, represent 11 haplotypes, respectively, with each solid line representing one step mutation.
Fig. 5.
Mismatch distribution analysis based on chloroplast DNA of wild Cymbidium tortisepalum in Yunnan Province (A), in northern region of three parallel rivers in the western Yunnan Province (B), in middle region of three parallel rivers in the western Yunnan Province (C), and in southern region of three parallel rivers in the western Yunnan Province (D). Solid lines indicate expected value (Exp), dotted lines reflect observed value (Obs).
Barrett,S.C.H. & Kohn,J.R.1991Genetic and evolutionary consequences of small population size in plants: implications for conservation130Genetics and conservation of rare plants.Oxford University PressNew York, NY
Birky,C.W.1988Evolution and variation in plant chloroplast and mitochondrial genomes2354Gottlieb,L.D. & Jain,S.K.Plant evolutionary biology.Chapman and HallLondon, UK
Crandall,K.A. & Templeton,A.R.1993Empirical tests of some predictions from coalescent theory with application to intraspecific phylogeny reconstructionGenetics1343959969https://doi.org/10.1101/gad.7.7b.145910.1093/genetics/134.3.959
Excoffier,L. & Lischer,H.E.L.2010Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and WindowsMol. Ecol. Resour.103564567https://doi.org/10.1111/j.1755-0998.2010.02847.x10.1111/j.1755-0998.2010.02847.x
Fu,Y.X1997Statistical tests of neutrality of mutations against population growth, hitchhiking and background selectionGenetics1472915925https://doi.org/10.1017/S001667239700296610.1093/genetics/147.2.915
Hewitt,G.M1999Post-glacial re-colonization of European biotaBiol. J. Linn. Soc. Lond.681-287112https://doi.org/10.1006/bijl.1999.033210.1111/j.1095-8312.1999.tb01160.x
Hewitt,G.M2004Genetic consequence of climate oscillations in the QuaternaryPhilos. Trans. R. Soc. Lond. Ser. B Biol. Sci.3591442183195https://doi.org/10.1098/rstb.2003.138810.1098/rstb.2003.1388
Huan,L.V.T. & Tanake,M.2004Effects of red and blue light-emitting diodes on callus induction, callus proliferation, and protocorm-like body formation from callus in Cymbidium orchidEnviron. Control Biol.4215764https://doi.org/10.2525/ecb1963.42.5710.2525/ecb1963.42.57
Jordano,P. & Godoy,J.A.2000RAPD variation and population genetic structure in Prunus mahaleb (Rosaceae), an animal-dispersed treeMol. Ecol.9912931305https://doi.org/10.1046/j.1365-294x.2000.01009.x10.1046/j.1365-294x.2000.01009.x
Larkin,M.A., Blackshields,G., Brown,N.P., Chenna,R., Mcgettigan,P.A. & Mcwilliam,H.2014Clustal W and Clustal X version 2.0Bioinformatics232129472948https://doi.org/10.1093/bioinformatics/btm404
Lascoux,M., Palme,A.E., Cheddadi,R. & Latta,R.G.2004Impact of ice ages on the genetic structure of tree and shrubsPhilos. Trans. R. Soc. Lond. Ser. B-Biol. Sci.3591442197207https://doi.org/10.1098/rstb.2003.139010.1098/rstb.2003.1390
Li,H.F., Zhao,Z.L. & Liu,G.M.2009Study on effects of different transplantation medium on test-tube plantlet of Cymbidium torliseplumJ. Anhui Agr. Sci.3799509951https://doi.org/CNKI: SUN:AHNY.0.2009-21-055
Librado,P. & Rozas,J.2009DnaSP v5: A software for comprehensive analysis of DNA polymorphism dataBioinformatics251114511452https://doi.org/10.1093/bioinformatics/btp18710.1093/bioinformatics/btp187
Lin,T.P., Liu,H.Y., Hsieh,C.F. & Wang,K.H.2016Complete list of the native orchids of Taiwan and their type informationTaiwania6178126https://doi.org/10.6165/tai.2016.61.78
Liu,Y., Yang,S.X., Ji,P.Z. & Gao,L.Z.2012Phylogeography of Camellia taliensis (Theaceae) inferred from chloroplast and nuclear DNA: Insights into evolutionary history and conservationBMC Evol. Biol.1292https://doi.org/10.1186/1471-2148-12-9210.1186/1471-2148-12-92
Maki,M., Yamashiro,T. & Matsumura,S.2003High levels of genetic diversity in island populations of the island endemic Suzukia luchuensis (Labiatae)Heredity91300306https://doi.org/10.1038/sj.hdy.680031110.1038/sj.hdy.6800311
Ning,H.J., Shao,L. & Dai,S.L.2020The complete chloroplast genome of Cymbidium longibracteatum (Orchidaceae) and phylogenetic analysisMitochondrial DNA B Resour.5329482949https://doi.org/10.1080/23802359.2020.1780968
Peakall,R. & Smouse,P.E.2006Genelex6: Genetic analysis in Excel. Population genetic software for teaching and researchMol. Ecol. Notes6288295https://doi.org/10.1111/j.1471-8286.2005.01155.x10.1111/j.1471-8286.2005.01155.x
Pons,O. & Petit,R.J.1996Measuring and testing genetic differentiation with ordered versus unordered allelesGenetics144312371245https://doi.org/10.1093/genetics/144.3.123710.1093/genetics/144.3.1237
Qi,F.Y2004Conservation and development strategy for Cymbidium lianpan resource of WeixiForest Invent. Plann.292105108https://doi.org/10.3969/j.issn.1671-3168.2004.02.027
Shaw,J., Lickey,B.E., Schilling,E.E. & Small,R.L.2007Comparison of whole chloroplast genome sequences to choose non-coding regions for phylogenetic studies in angiosperms: The tortoise and the hare IIIAmer. J. Bot.943275288https://doi.org/10.3732/ajb.94.3.27510.3732/ajb.94.3.275
Spielman,D.B., Brook,W. & Frankham,R.2004Most species are not driven to extinction before genetic factors impact themProc. Natl. Acad. Sci. USA1011562115624https://doi.org/10.1073/pnas.0403809101
Sun,Y., Chen,Y., Lin,H. & Ma,X.2019The complete chloroplast genome of a Cymbidium tortisepalum (Orchidaceae) male mutantMitochondrial DNA B Resour.4240874088https://doi.org/0.1080/23802359.2019.169194810.1080/23802359.2019.1691948
Suo,M., Bi,Y.F., Wang,H.Z., Tang,M. & Ye,G.2016Relationship between morphology and elevation in native Cymbidium tortisepalum in Yunnan Province, ChinaActa Ecol. Sin.361961706177https://doi.org/10.5846/stxb201507261561
Tajima,F1993Simple methods for testing the molecular evolutionary clock hypothesisGenetics1352599607https://doi.org/10.1101/gad.7.10.204810.1093/genetics/135.2.599
Tajima,F1996The amount of DNA polymorphism maintained in a finite population when the neutral mutation rate varies among sitesGenetics14331475https://doi.org/10.1007/BF00120227
Tang,M2014Wild population genetic structure and cultivars molecular identification of Cymbidium tortisepalum in Yunnan ProvincePh.D. Diss.Yunnan Agr. Univ.Kumming, China
Templeton,A.R1998Nested clade analyses of phylogeographic data: Testing hypotheses about gene flow and population historyMol. Ecol.74381397https://doi.org/10.1046/j.1365-294x.1998.00308.x10.1046/j.1365-294x.1998.00308.x
Wroblewska,A., Brzosko,E., Czarnecka,B. & Nowosielski,J.2003High levels of genetic diversity in populations of Iris aphylla L. (Iridaceae), an endangered species in PolandBiol. J. Linn. Soc. Lond.1426572https://doi.org/10.1046/j.1095-8339.2003.00162.x
Xue,R.G., He,S.X., Li,Z.G., He,J.W., Yang,S.H. & Li,Y.2007A review of conservation status of Cymbidium tortisepalum and the utilizationForest Invent. Plann.3239497https://doi.org/10.3969/j.issn.1671-3168.2007.03.024
Xun,G., Luan,S.S., Hung,K.H., Hwang,C.C., Lin,C.J., Chiang,Y.C. & Chiang,T.Y.2011Population structure of Nouelia insignis (Asteraceae), an endangered species in southwestern China, based on chloroplast DNA sequences: Recent demographic shrinkingJ. Plant Res.1242221230https://doi.org/10.1007/s10265-010-0363-010.1007/s10265-010-0363-0
Yang,J.B., Tang,M., Li,H.T., Zhang,Z.R. & Li,D.Z.2013Complete chloroplast genome of the genus Cymbidium: Lights into the species identification, phylogenetic implications and population genetic analysesBMC Evol. Biol.1384https://doi.org/10.1186/1471-2148-13-8410.1186/1471-2148-13-84
Ye,G., Tang,M., Wang,H.Z., Suo,M., Che,W.G., Chen,X. & Bi,Y.F.2015Relationship between morphological differentiation and elevation of wild Cymbidium tortisepalum in Diqing of Yunnan ProvinceJ. Yunnan Agr. Univ.3018489https://doi.org/10.16211/j.issn.1004-390X(n).2015.01.014
Zhao,Y., Li,D. & Liu,T.2019Pollination-Induced transcriptome and phylogenetic analysis in Cymbidium tortisepalum (Orchidaceae)Russ. J. Plant Physiol.66618627https://doi.org/10.1134/S102144371904017410.1134/S1021443719040174
Zhao,Y., Tang,M. & Bi,Y.F.2017Nuclear genetic diversity and population structure of a vulnerable and endemic orchid (Cymbidium tortisepalum) in Northwestern Yunnan, ChinaScientia Hort.2192230https://doi.org/10.1016/j.scienta.2017.02.03310.1016/j.scienta.2017.02.033
This study was supported by the National Science Foundation of China: Genetic Response Mechanism of Main Ecological Factors for Populations of Wild Cymbidium Tortisepalum in Yunnan Province (no. 31260487) and the Germplasm Bank of Wild Species in Southwest China.
Simplified map showing the geographic distribution of the 28 sampling sites of Cymbidium tortisepalum in Yunnan Province, China. Blue marks represent populations in the northern region of three parallel rivers in the western Yunnan Province (NW); green marks markers represent those in the middle region of three parallel rivers in the western Yunnan Province (MD), and red marks represent those in the southern region of three parallel rivers in the western Yunnan Province. The abbreviations of the sites refer to Table 1.
Fig. 2.
Geographic distribution of the chloroplast DNA haplotypes of wild Cymbidium tortisepalum in Yunnan Province, China. Different colors in the label box represent different haplotypes, and the pie charts represent the distribution proportion of each haplotype in the populations.
Fig. 3.
Neighbor-joining trees obtained by analysis of 11 chloroplast DNA haplotypes of wild Cymbidium tortisepalum in Yunnan Province, China. Hap1–Hap11 represent 11 haplotypes, respectively. Bootstrap values (>50%) are based on 1000 replicates and are indicated above branches.
Fig. 4.
The haplotype network of 11 chloroplast DNA haplotypes detected in 28 populations of wild Cymbidium tortisepalum in Yunnan Province, China. The size of circles corresponds to the frequency of each haplotype, H1–H11, represent 11 haplotypes, respectively, with each solid line representing one step mutation.
Fig. 5.
Mismatch distribution analysis based on chloroplast DNA of wild Cymbidium tortisepalum in Yunnan Province (A), in northern region of three parallel rivers in the western Yunnan Province (B), in middle region of three parallel rivers in the western Yunnan Province (C), and in southern region of three parallel rivers in the western Yunnan Province (D). Solid lines indicate expected value (Exp), dotted lines reflect observed value (Obs).