Chenxing Zhang1,2,3
Chunfeng Jia4Xinru Liu1,2,3Hanqing Zhao1,2,3
Lu Hou1,2,3
Meng Li1,2,3Binbin Cui4*
Yingyue Li1,2,3*- 1National Engineering Research Center of Tree Breeding and Ecological Restoration, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- 2Key Laboratory of Genetics and Breeding in Forest Trees and Ornamental Plants, Ministry of Education, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- 3The Tree and Ornamental Plant Breeding and Biotechnology Laboratory of National Forestry and Grassland Administration, Beijing Forestry University, Beijing, China
- 4College of Biochemistry and Environmental Engineering, Baoding University, Baoding, China
Elsholtzia stauntonii Benth. (Lamiaceae) is an economically important ornamental, medicinal and aromatic plant species. To meet the increasing market demand for E. stauntonii, it is necessary to assess genetic diversity within the species to accelerate the process of genetic improvement. Analysis of the transferability of simple sequence repeat (SSR) markers from related species or genera is a fast and economical method to evaluate diversity, and can ensure the availability of molecular markers in crops with limited genomic resources. In this study, the cross-genera transferability of 497 SSR markers selected from other members of the Lamiaceae (Salvia L., Perilla L., Mentha L., Hyptis Jacq., Leonurus L., Pogostemon Desf., Rosmarinus L., and Scutella L.) to E. stauntonii was 9.05% (45 primers). Among the 45 transferable markers, 10 markers revealed relatively high polymorphism in E. stauntonii. The genetic variation among 825 individuals from 18 natural populations of E. stauntonii in Hebei Province of China was analyzed using the 10 polymorphic SSR markers. On the basis of the SSR data, the average number of alleles (NA), expected heterozygosity (HE), and Shannon’s information index (I) of the 10 primers pairs were 7.000, 0.478, and 0.688, respectively. Lower gene flow (Nm = 1.252) and high genetic differentiation (Fst = 0.181) were detected in the populations. Analysis of molecular variance (AMOVA) revealed that most of the variation (81.47%) was within the populations. Integrating the results of STRUCTURE, UPGMA (Unweighted Pair Group Method with Arithmetic Mean) clustering, and principal coordinate analysis, the 825 samples were grouped into two clusters associated with geographical provenance (southwestern and northeastern regions), which was consistent with the results of a Mantel test (r = 0.56, p < 0.001). Overall, SSR markers developed in related genera were effective to study the genetic structure and genetic diversity in geographical populations of E. stauntonii. The results provide a theoretical basis for conservation of genetic resources, genetic improvement, and construction of a core collection for E. stauntonii.
Introduction
Plants provide the fundamental support system for life on Earth and plant diversity is a key property of ecosystems (Moreira et al., 2016). Lamiaceae, with more than 7,000 species around the world, is a large family of flowering plants and an important part of plant diversity (Harley et al., 2004). Some species of the family have been commercially exploited worldwide (e.g., Lavandula angustifolia) (Pistelli et al., 2017), while many others are known locally for their ornamental, biological, and medicinal properties (e.g., Hedeoma patens) (Leyva-Lopez et al., 2016). Elsholtzia stauntonii Benth., a perennial deciduous subshrub of the Lamiaceae family, is widely distributed in China, especially in Hebei Province (Editorial committee of the Flora of China, and Chinese Academy of Sciences, 1977). Elsholtzia stauntonii is aromatic and its essential oils show insecticidal activity against certain pests of stored products (Lv et al., 2012). Consequently, E. stauntonii is used in novel biological insecticides (Liang et al., 2021). With its dense flowers, late flowering and bright colors, E. stauntonii is a popular bedding plant cultivated in public gardens and a wild nectar plant (Tian et al., 2007). Its branches and leaves are also used as a traditional Chinese medicinal material for treating common colds, toothache, and digestive system diseases (Zheng et al., 1999; Pekhova et al., 2020). Furthermore, the species is of ecological importance owing to its strong tolerance of drought and barren environments (Si and Peng, 2017). Previous studies on E. stauntonii have mainly focused on its morphological traits (Zhou et al., 2016), cultivation and propagation (Zhou et al., 2014), and extraction and identification of volatile oils (Guo et al., 2012; Xing et al., 2019). Little information is available on the genome and genetic diversity of this species. Studies of genetic diversity have been performed on many Lamiaceae species to clarify the genetic relationships among species (Bahadirli and Ayanoglu, 2020) and to identify molecular markers associated with agronomic traits for marker-assisted breeding (Lim et al., 2021). Therefore, with recognition of the value of this important but neglected species, it is crucial to evaluate the genetic variation in E. stauntonii.
Understanding the genetic diversity and structure of germplasm resources is a prerequisite for formulating effective breeding and conservation strategies (Gepts, 2006; Engelhardt et al., 2014; Vanavermaete et al., 2021). In recent years, with advances in molecular biological technologies, molecular marker-based analysis of genetic variation has avoided the limitations of morphological and biochemical indicators that are readily affected by environmental factors, and are more accurate for evaluation of population genetic parameters (Agarwal et al., 2008; Nadeem et al., 2018). Among the various types of molecular markers, microsatellites (or simple sequence repeats; SSRs) have become the preferred choice because of their high degree of reproducibility, ability to identify high levels of genetic polymorphism, codominant inheritance, and abundance in plant genomes (Powell et al., 1996; Grover and Sharma, 2016; Jose et al., 2018). Currently, SSR markers have been widely used in assessments of population genetic diversity (Gadissa et al., 2018), cultivar identification (Pinto et al., 2018; Sulu et al., 2020), and association mapping (Yang et al., 2020) as effective tools to describe genetic variation in plants. Using available genomic or expressed sequence tag databases (Vieira et al., 2016), a large number of SSR markers have been developed in Lamiaceae species, such as Scutellaria baicalensis (Zhang et al., 2014), Perilla frutescens (Ha et al., 2021), and Salvia splendens (Jiao et al., 2020). However, owing to the lack of genetic information, to date no SSR markers are available for evaluation of the genetic diversity within E. stauntonii, which limits its breeding and commercial development.
An alternative approach to overcome the time-consuming and costly development of SSR primers for target species is to exploit the transferability of SSR primers between related species or genera (Peakall et al., 1998; Pan et al., 2018). The transfer of SSR markers has been successful in many cases, as documented among Lamiaceae species (Karaca et al., 2013), from Gossypium hirsutum and Corchorus olitorius to Malvaceae species (Satya et al., 2016), among Allium species and among members of the Alliaceae (Barboza et al., 2018), and from Ricinus communis to Euphorbiaceae species (Dharajiya et al., 2020).
In this study, cross-genera amplification of SSR markers in Lamiaceae species was used to enrich the genetic database of E. stauntonii. On this basis, the genetic diversity and population structure of 825 E. stauntonii individuals sampled from 18 regions in Hebei Province, China, which is an important center in the distribution range of the species, were successfully evaluated using 10 highly polymorphic cross-transferable SSR markers. The results provide insight into the genetic resources of E. stauntonii, and represent a foundation for the effective management and genetic improvement of the species.
Materials and Methods
Plant Materials
A field survey of extant E. stauntonii populations in Hebei Province, China was conducted from July to September 2020. Samples of 825 mature and disease-free individuals were randomly collected from 18 populations in the natural distribution area; the sample size per population ranged from 43 to 50 individuals (Table 1). To ensure sufficient geographical representation of the samples, the distance between sampled individuals was more than 50 m. The GPS coordinates were recorded for each population. The altitude of all sampled populations ranged from 540 to 1,200 m, with an average altitude of 937 m. Fresh leaves were desiccated in silica gel and stored at –20°C until use.
Table 1. Location and sampling site characteristics of 18 Elsholtzia stauntonii geographical populations (825 individuals).
DNA Extraction and Genotyping With Simple Sequence Repeats Markers
Total genomic DNA of all samples was extracted using a Plant Genomic DNA Kit (Tiangen Biotech, Beijing, China). The DNA quality was evaluated by 1.0% agarose gel electrophoresis. DNA quantification was conducted with a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, United States) and the DNA concentration was adjusted to 20 ng/μL.
In total, 497 SSR markers from eight genera of Lamiaceae were used to perform cross-amplification tests on 16 randomly selected E. stauntonii samples from different populations (Table 2). Primers were synthesized by RuiBiotech (Beijing, China). To allow fluorescent labeling, the universal M13 sequence was added to the 5′ end of all forward primers. The PCR mixture was prepared in accordance with the method of Zhang et al. (2015) with modifications. Briefly, the PCR mixture (20 μl) contained 10 μl of 2 × Taq PCR Mix, 2 μl (40 ng) genomic DNA, 4 μl (4 pmol) fluorescent-dye-labeled M13 primer, and 4 μl (4 pmol) mixed complementary forward and reverse primers. The temperature profile used was as follows: an initial step of 5 min at 94°C, followed by 35 cycles of 30 s at 94°C, 30 s at 54°C, and 1 min at 72°C, then eight cycles of 30 s at 94°C, 30 s at 53°C, and 30 s at 72°C, and a final step of 10 min at 72°C. Capillary electrophoresis of the PCR products was performed using an ABI 3730xl DNA Analyzer (Tiangen Biotech, Beijing, China). GeneMarker version 2.2.0 software (Holland and Parson, 2011) was used to analyze the fragment size of the products and determine whether the primers detected polymorphisms.
Table 2. SSR markers from eight genera in the Lamiaceae used to evaluate cross-genera transferability to Elsholtzia stauntonii.
Data Analysis
The SSR markers that successfully amplified single and specific bands within the expected product size range were considered to be transferable. The percentage transferability for each donor genus was calculated as (number of SSRs transferred/total number of SSRs screened) × 100. Similarly, the percentage transferability of SSR markers classified by repeat motifs (di-/tri-/tetra-/penta-/hexa-/complex-nucleotide) was estimated. Population diversity and population structure of E. stauntonii were assessed based on the transferred microsatellite loci.
For each SSR locus, important indices of genetic diversity were calculated using GenAlEx version 6.5 (Peakall and Smouse, 2012) and Microsatellite-Toolkit (Park, 2001) software. The indices comprised number of alleles (NA), number of effective alleles (NE), Shannon’s information index (I), observed heterozygosity (HO), expected heterozygosity (HE), inbreeding coefficient (Fis), gene flow (Nm), and genetic differentiation index (Fst). The polymorphic information content (PIC) value and Hardy–Weinberg equilibrium (HWE) for the loci were calculated using PowerMarker version 3.25 (Liu and Muse, 2005).
Analysis of molecular variance (AMOVA) was conducted using Arlequin version 3.5 (Excoffier and Lischer, 2010) to estimate the variance components of genetic variation among and within populations. Bayesian analysis was performed to evaluate the population genetic structure and detect the most likely number of population genetic clusters of E. stauntonii using STRUCTURE version 2.3.4 (Pritchard et al., 2000). For each simulated value of K (range from 1 to 18), 15 independent runs were performed with a burn-in period of 200,000 iterations followed by 1,000,000 Markov chain Monte Carlo repetitions. The ΔK method (Evanno et al., 2005) was implemented in the Structure Harvester (Earl and VonHoldt, 2012) program based on the STRUCTURE results to determine the optimum K value. The percentage membership of each individual in every cluster (Q value) was determined; an individual with a Q value higher than 0.80 was considered to have a single genetic component (a pure individual). Nei’s genetic distance (DA) (Nei et al., 1983) among the 825 individuals and 18 populations was calculated using PowerMarker version 3.25 (Liu and Muse, 2005). To explore the genetic relationships among all samples, the DA matrix was used to construct a dendrogram by hierarchical clustering with the unweighted pair group method with arithmetic mean (UPGMA) using PowerMarker version 3.25 (Liu and Muse, 2005). The dendrograms were visualized, manipulated, and annotated with the Interactive Tree of Life online tool (Letunic and Bork, 2007). Based on the standardized covariance of genetic distance, a principal coordinate analysis (PCoA) was conducted with GenAlEx version 6.5 (Peakall and Smouse, 2012). A clustering heatmap was prepared with Microsoft Excel based on the pairwise Nei’s genetic distance between populations and the Fst values.
The geographic distance between the sampling locations was calculated with Geographic Distance Matrix Generator version 1.2.3 software1. The Mantel test was performed to analyze the correlation between geographic and genetic distance of different populations using the “ggpubr2” and “ggplot2” (Ginestet, 2011) packages in R version 4.1.1 (Ihaka and Gentleman, 1996). The map was created using the “sf” (Pebesma, 2018), “ggplot2” (Ginestet, 2011), and “ggspatial3” packages in R version 4.1.1.
Results
Transferability of Simple Sequence Repeats Markers
Of the 497 SSRs mined, 45 primers (9.05%) successfully amplified genomic DNA of E. stauntonii and produced PCR products of the expected size (Supplementary Table 1). Transferability of the SSR markers varied among eight genera of Lamiaceae: 26.67% of the SSR markers of Hyptis Jacq., 25.00% of Leonurus L., 16.67% of Pogostemon Desf., 13.28% of Perilla L., 8.33% of Rosmarinus L., 6.73% of Salvia L., and 4.00% of Mentha L. were transferable to E. stauntonii, whereas no SSR marker was transferable from Scutellaria L. (Figure 1). The transferability of SSR markers based on simple di-/tri-nucleotide repeat motifs (7.89–9.65%) was higher than that based on tetra-/penta-/hexa-nucleotide repeat motifs (Figure 2). The SSR markers based on complex (simple imperfect, compound perfect/imperfect) nucleotide repeats also showed high transferability (24.00%).
Figure 1. Confamiliar transferability of SSR markers to Elsholtzia stauntonii from other genera of Lamiaceae.
Figure 2. Transferability of SSR markers to Elsholtzia stauntonii based on SSR repeat motif length. The repeat motif of eighteen SSRs was not available.
Polymorphism Analysis Based on Transferable Simple Sequence Repeats Markers
Among the 45 transferable SSR markers, 35 markers were monomorphic (33 primers) or showed very low levels of polymorphism (two primers) based on capillary electrophoresis results (Supplementary Table 1). Ten (2.01%) primer pairs showed higher polymorphism and were used for analysis of genetic diversity in E. stauntonii (Table 3).
Table 3. Characterization of 10 polymorphic SSR primer pairs.
In total, 70 alleles were amplified by the 10 polymorphic SSR markers among 825 E. stauntonii individuals, with a mean of 7.000 observed alleles per locus, and ranging in length from 132 to 350 bp (Table 3). The number of alleles per marker varied from 2 (ssps344, GBPF203, and C3787) to 24 (ssps251). Among the 10 markers, the observed heterozygosity (HO) of eight markers was lower than the expected heterozygosity (HE), and the average inbreeding coefficient (Fis) was positive. The PIC values ranged from 0.191 to 0.868, with a mean of 0.429. Eight SSR markers showed moderate or high polymorphism levels (PIC > 0.25) and two markers (ssps344 and KNUPF67) showed low polymorphism levels (PIC < 0.25) among the tested E. stauntonii populations. All SSR genotyping data for the 10 loci showed strongly significant deviation from the HWE (p < 0.01) (Table 3).
Population Genetic Diversity and Differentiation
In general, the difference in genetic diversity among the 18 geographical populations of E. stauntonii was observed (Table 4). The genetic diversity indices NA and NE of individual populations ranged from 2.400 (XW) to 3.900 (XM), and from 1.406 (ZQ) to 2.532 (CX), with means of 3.267 and 1.971, respectively. The average HO and HE were 0.359 and 0.388, ranging from 0.253 (XW) to 0.478 (CX), and from 0.260 (ZQ) to 0.506 (CX), respectively. The I value ranged from 0.450 (ZQ) to 0.887 (CX), with an average of 0.688. The highest level of genetic diversity was detected in CX, whereas XW and ZQ showed the lowest genetic diversity (Table 4).
Table 4. Genetic diversity for 10 polymorphic SSR markers in 18 natural populations of Elsholtzia stauntonii (total of 825 individuals).
The average Fst, a measure of the degree of genetic differentiation among populations, was 0.181 and the average Nm was 1.252 (Table 4). These results were consistent with those of AMOVA, which revealed that 18.53% of the genetic variation occurred among the populations and 81.47% was present within the populations (Table 5).
Table 5. Analysis of molecular variance (AMOVA) for 825 Elsholtzia stauntonii individuals clustered in 18 populations.
The pairwise Fst ranged from 0.020 (between populations CY and ZC) to 0.199 (between populations CX and ZQ) (Figure 3). The pairwise Nei’s genetic distance ranged from 0.034 to 0.359, with the highest value (0.359) between the CX and ZQ populations and the second highest (0.328) between the CX and XW populations (Figure 3). In general, populations that were geographically closer had lower genetic distances (SZ–XM, 52.43 km, 0.034), whereas higher genetic distances were observed for populations that were farther apart (CX–XW, 490.81 km, 0.328). Mantel test results revealed a strongly significant correlation between genetic and geographic distance among the 18 natural populations (r = 0.56, p < 0.001) (Figure 4).
Figure 3. Clustering heatmap of pairwise genetic differentiation index (Fst) (values below the diagonal), and pairwise Nei’s genetic distance (values above the diagonal) among 18 populations of Elsholtzia stauntonii.
Figure 4. Mantel test results of the correlation between Nei’s genetic distance and geographical distance among 18 populations of Elsholtzia stauntonii.
Genetic Relationship and Population Structure
To determine the genetic structure of E. stauntonii populations, the co-ancestral relationship among populations was analyzed based on a Bayesian assignment model. The Q value of 680 individuals, accounting for 82.42% of the total samples, was greater than 0.8, indicating that most individuals had a relatively simple genetic background (Supplementary Table 2). STRUCTURE analysis revealed a distinct peak for ΔK at K = 2, indicating that the E. stauntonii samples could be grouped into two major clusters; these two subpopulations roughly corresponded to the southwestern and northeastern regions of Hebei Province (Figures 5A,B, 6). The four populations (ZQ, ZF, BY, and BL) situated at the junction of the two regions contained abundant genetic variation and the lineages were obviously mixed (Figure 6). These populations contained genetic information of dual origins, with an average assignment of 49.7% to Cluster I (red) and 50.3% to Cluster II (green), whereas the other populations were relatively independent (Figure 6 and Supplementary Table 3).
Figure 5. Population structure of 825 Elsholtzia stauntonii individuals based on 10 SSR markers. (A) Delta K values from STRUCTURE analysis of the E. stauntonii individuals. (B) Histogram from STRUCTURE analysis for the model with K = 2.
Figure 6. Geographical distribution of 18 populations of Elsholtzia stauntonii and lineage contribution to gene pools. The pie diagram was drawn based on the results of STRUCTURE analysis with K = 2 and represents the average proportion of cluster membership of all individuals in each population.
The UPGMA dendrogram grouped the 18 E. stauntonii populations into three major clusters. The populations XW, HS, XM, and SZ formed one cluster, the populations BF, SP, BZ, BL, BY, ZF, and ZQ formed a separate cluster, and the remaining populations (ZC, CY, CC, CF, CL, QD, and CX) formed a third group (Figure 3). However, the UPGMA tree of 825 individuals contained clusters of individuals from different populations and the samples were grouped into two major clades (Figure 7A). Consistent with the population dendrogram, the individuals of seven populations (ZC, CY, CC, CF, CL, QD, and CX) in the northeastern region were merged into a cluster and showed a broadly mixed ancestry, whereas the majority of individuals from the 11 populations in the southwestern region were grouped in the same cluster. In general, genetic proximity was observed between neighboring populations, whereas few genotypes shared an ancestry with individuals in geographically distant populations (Figure 7A).
Figure 7. (A) UPGMA tree constructed based on Nei’s genetic distance of 825 Elsholtzia stauntonii individuals. (B) Principal coordinate analysis (PCoA) of 825 individuals from 18 populations.
Principal coordinate analysis was performed to visualize the multidimensional relationships among individuals based on genetic distance (Figure 7B). The closer they were on the PCoA graph, the less genetically differentiated they were, such as the HS, XM and SZ populations. Principal coordinates (PC) 1 and 2 explained 25.51% and 23.83% of the variance, respectively, and a combined 49.34% of the total variation. The PCoA analysis separated the populations into two clusters, which further confirmed the presence of genetic structure among E. stauntonii populations (Figure 7B).
Discussion
Cross-Genera Transferability of Microsatellite Markers
Cross-taxa transferability enables the development of microsatellite markers, based on markers from related species, at low cost (Kalia et al., 2011; Bernardes et al., 2021). In the present study, the transferability observed (9.05%) was similar to the average of approximately 10% reported in cross-genera transferability studies of eudicots between 1997 and mid-2006 (Barbará et al., 2007). However, the percentage transferability determined here was lower than the cross-genera amplification percentage observed in some families, such as Bignoniaceae (40.58%) (Kalia et al., 2020) and Cactaceae (35.16%) (Bombonato et al., 2019), and much lower than that between species within the same genus (Miranda et al., 2020; Pern et al., 2020; Li et al., 2021). The transferability of SSR markers between species or genera is determined by the conservation of DNA sequences and the stability of primer binding sites in flanking regions of SSRs during evolution (Ellegren, 2000; Saeed et al., 2016). A higher cross-taxa amplification percentage of SSR markers is expected between species with a close phylogenetic relationship, especially when both donor and target species are congeneric (Lesser et al., 2012; Karaca et al., 2013; Bharti et al., 2018). Therefore, it was speculated that E. stauntonii might have a relatively distant genetic relationship with its donor species (Jarne and Lagoda, 1996; Li et al., 2017). Nevertheless, although the primary contributing factor is phylogenetic proximity, transferability may also be influenced by other factors, such as the size and complexity of the relevant genome and whether the microsatellite is sited within a coding region (Oliveira et al., 2006).
The SSR transferability is also affected by the base repeat type of SSR core units (Gao et al., 2013). In general, SSR markers with a trinucleotide repeat motif show higher transferability than those with dinucleotide-type repeat motifs (Ellis and Burke, 2007; Harijan et al., 2017). A similar trend in heterologous SSR marker transferability has been reported for Euphorbiaceae, with the transfer rate of simple trinucleotide repeats being 11.44% higher than that of dinucleotide repeats (Dharajiya et al., 2020). One reason might be that the selection of frameshift mutations limits the fixation of microsatellites whose motif lengths are not divisible by three (e.g., di-/tetra-/pentanucleotide repeats), whereas fixation of mutant microsatellite whose motif lengths are divisible by three (tri-/hexanucleotide repeats) is not affected by differential selection pressures in coding and non-coding regions (Metzgar et al., 2000; Dutta et al., 2011; Shivakumar et al., 2017). Moreover, microsatellites based on a complex nucleotide sequence show a lower mutation rate and higher transferability than simple repeats (Jin et al., 1996; Kalia et al., 2020). Given that more than 63% of the SSR markers used in the present study were di-nucleotide motif repeats, this may be an additional reason for the low transferability. Regardless, the 10 polymorphic transferable SSR markers identified in this study could be used to accelerate analysis of the genetic diversity of E. stauntonii and related species.
Genetic Diversity and Population Structure Analysis
All 10 transferable SSR primers used in the present study deviated significantly from the HWE (p < 0.01), which might reflect that the samples were collected from natural populations affected by complex environmental factors and human activities (Wang and Shete, 2012; Peng et al., 2021). The investigation found that E. stauntonii mainly resided beside rivers and valleys at lower altitudes, and its living environment was easily disturbed by tourism, grazing, and farmland planting. The mean PIC was 0.429, indicating that the genetic diversity of E. stauntonii germplasm in Hebei Province was generally moderate (0.25 < PIC < 0.5) (Botstein et al., 1980), which was lower than that reported for Scutellaria baicalensis (0.72) (Qi et al., 2015) and Perilla frutescens (0.582) (Park et al., 2019). The genetic diversity possessed by a species is determined by its evolutionary and ecological history (Aboukhalid et al., 2017). For E. stauntonii, the wide original geographic distribution theoretically contributes to the high level of genetic diversity detected. The moderate genetic diversity observed in this study might be limited by the number of markers, the test method, or the regions of sample collection (Budak et al., 2004; Serrote et al., 2020). Furthermore, it is generally accepted that markers of high transferability are less polymorphic among closely related species (Ni et al., 2015).
Genetic variation of a population can be measured by the number of alleles and heterozygosity (Huang et al., 2020; Savić et al., 2021). The mean NA (7.000) and NE (1.971) observed in the present study were slightly higher than the values (6.45 and 1.67, respectively) for 12 Plectranthus edulis populations evaluated using 20 SSR markers (Gadissa et al., 2018), but lower than the values reported for Salvia splendens (NA = 11.829, NE = 3.261) (Jiao et al., 2020) and Juglans nigra (NA = 9.99, NE = 7.51) (Schneider et al., 2019). Compared with previous studies on Lamiaceae species, E. stauntonii (HO = 0.360) has higher genetic diversity than Salvia fruticosa (HO = 0.307) (Radosavljević et al., 2015) and Hedeoma piperita (HO = 0.3) (Herrera-Arroyo et al., 2020), but a higher level of genetic diversity has been observed in Rosmarinus officinalis (HO = 0.63) (Segarra-Moragues et al., 2015). These results also indicated that the tested E. stauntonii samples showed moderate genetic diversity. The observed heterozygosity in the present study was lower than the expected heterozygosity (HE = 0.478), indicating that the number of heterozygotes was lower than the predicted value (García-Arias et al., 2018). The average inbreeding coefficient (Fis = 0.192) in the population was positive and high, which similarly suggested there was an overabundance of homozygotes (Liu et al., 2019; Bairu et al., 2021). The reproductive characteristics of species might be an important factor in heterozygote deficiency in populations; for example, the combination of hermaphroditic flowers and profuse flowering of E. stauntonii may lead to improved levels of inbreeding which can reduce the genetic diversity of this widely distributed species (Zigene et al., 2021).
Genetic differentiation is an important index of the genetic structure of populations within a species (Wang, 2020). In the current study, gene flow of E. stauntonii (Nm = 1.252) was moderately low, lower than that reported for Eragrostis tef (Nm = 4.742) (Tadesse et al., 2021) and Mucuna pruriens (Nm = 7.48) (Lepcha and Sathyanarayana, 2021), but it was consistent with the high genetic differentiation (Fst = 0.181) (Balloux and Lugon-Moulin, 2002). The low gene flow, which might be caused by the small sampling range and habitat fragmentation, increased the within-population genetic variation and promoted population differentiation (Young et al., 1996; Li et al., 2014; Peng et al., 2021). Notably, among the 18 populations, the greatest genetic differentiation was observed between the CX and ZQ populations, which may be considered as potential parents for hybrid breeding. Furthermore, AMOVA revealed that the variation among E. stauntonii populations (18.53%) was higher than that of many other plant species, such as Origanum vulgare subsp. hirtum (8%) (Alekseeva et al., 2021) and Allium sativum L. (12.29%) (Poljuha et al., 2021), suggesting that disruption of gene flow between the populations played a role in determining the genetic differentiation of populations. Similarly, the results of the Mantel test (r = 0.56, p < 0.001) showed that increase in geographic distance reduced gene flow. Therefore, geographic distance may be the principal factor for genetic differentiation in E. stauntonii.
In the UPGMA clustering analysis, a certain deviation was observed in the clustering of individuals and of populations, which might be due to the high degree of gene exchange between populations. Overall, the clustering results observed with STRUCTURE, PCoA, and UPGMA were essentially similar, indicating that there might be two gene pools of E. stauntonii. Examination of the topography and geomorphology in Hebei Province indicated that the genetic structure of E. stauntonii was delimited according to geographical location and topographic barriers: one genealogical lineage each was located near the Taihang Mountains and the Yanshan Mountains. Previous studies have shown that geographical barriers play a crucial role in preventing plant dispersal (Fan et al., 2013; Zhao and Gong, 2015; Zeng et al., 2018), for example, the Tanaka Line in southwest China has a marked effect on the divergence and evolution of plant species (Ju et al., 2018). In this study, populations ZQ and CX on both sides of the mountain ranges presented the greatest Fst and genetic distance values, although their geographic distance was relatively close (194.14 km apart from each other), which also indicated the role of geographical barriers in genetic differentiation. The Taihang and Yanshan mountains may be responsible for maintaining the major southwest–northeast split of E. stauntonii populations by hindering animal-mediated gene exchange. Thus, natural adaptation and geographical barriers may explain the differentiation of the two regions. Moreover, in the present analyses, each group was more or less intermingled with lineages of the other group, indicating that there was a certain level of gene flow between the populations, which has also been reported in Tectona grandis (Balakrishnan et al., 2021). The levels of genetic admixture differed substantially from population to population and seemed to be associated with the geographic location of the population, which indicated that interbreeding was an ongoing process, especially for the ZF, ZQ, BY, and BL populations (Mangaravite et al., 2016). These four populations had relatively high gene introgression, indicating that these populations had relatively rich genetic backgrounds and might be potential centers of diversification in this species. From the perspective of cultivar development and conservation, these populations may be valuable because they may contain novel alleles.
Conclusion
In this study, SSR markers from seven genera of Lamiaceae were transferable to E. stauntonii. Analysis based on 10 transferable SSR markers showed that, although the genetic diversity in 18 geographical populations of E. stauntonii was slightly lower than that of certain other species in the Lamiaceae, the genetic diversity was generally rich and the populations showed a high degree of genetic differentiation. STRUCTURE, PCoA, and UPGMA clustering analyses of all individuals resolved two genetic clusters, which might reflect geographical barriers to gene flow. By screening the study materials with superior traits and rich genetic variation, some individuals in the present experiment are expected to be utilized in the future improvement programmes of E. stauntonii to meet the increasing demand from the floriculture, essential oil, pharmaceutical, and other industries. The information obtained from the analysis of genetic diversity and population structure provides a foundation for the conservation, management, and utilization of E. stauntonii.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author Contributions
YL and BC conceived and designed the experiments and obtained the funding. CZ, XL, and HZ conducted the experiments. CJ and BC contributed to the investigation. CZ, LH, and ML collected and analyzed the data. CZ wrote the manuscript. YL provided suggestions and comments for the manuscript. All authors read and approved the final manuscript.
Funding
This work was financially supported by the Key Research and Development Program of Hebei Province (20326332D).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.903674/full#supplementary-material
Footnotes
- ^ http://biodiversityinformatics.amnh.org/opensource/gdmg/index.php
- ^ https://cran.r-project.org/package=ggpubr
- ^ https://cran.r-project.org/package=ggspatial
References
Aboukhalid, K., Machon, N., Lambourdiere, J., Abdelkrim, J., Bakha, M., Douaik, A., et al. (2017). Analysis of genetic diversity and population structure of the endangered Origanum compactum from Morocco, using SSR markers: implication for conservation. Biol. Conserv. 212, 172–182. doi: 10.1016/j.biocon.2017.05.030
Agarwal, M., Shrivastava, N., and Padh, H. (2008). Advances in molecular marker techniques and their applications in plant sciences. Plant Cell Rep. 27, 617–631. doi: 10.1007/s00299-008-0507-z
Alekseeva, M., Zagorcheva, T., Rusanova, M., Rusanov, K., and Atanassov, I. (2021). Genetic and flower volatile diversity in natural populations of Origanum vulgare subsp. hirtum (Link) Ietsw. in Bulgaria: toward the development of a core collection. Front. Plant Sci. 12:679063. doi: 10.3389/fpls.2021.679063
Bahadirli, N. P., and Ayanoglu, F. (2020). Genetic diversity of Salvia species from Turkey assessed by microsatellite markers. J. Appl. Res. Med. Aromat. Plants 20:100281. doi: 10.1016/j.jarmap.2020.100281
Bairu, M. W., Amelework, A. B., and Coetzer, W. G. (2021). Genetic diversity and population structure of six South African Acacia mearnsii breeding populations based on SSR markers. J. Plant Res. 134, 1243–1252. doi: 10.1007/s10265-021-01331-2
Balakrishnan, S., Dev, S. A., Sakthi, A. R., Vikashini, B., Bhasker, T. R., Magesh, N. S., et al. (2021). Gene-ecological zonation and population genetic structure of Tectona grandis L.f. in India revealed by genome-wide SSR markers. Tree Genet. Genomes 17:33. doi: 10.1007/s11295-021-01514-x
Balloux, F., and Lugon-Moulin, N. (2002). The estimation of population differentiation with microsatellite markers. Mol. Ecol. 11, 155–165. doi: 10.1046/j.0962-1083.2001.01436.x
Barbará, T., Palma-Silva, C., Paggi, G. M., Bered, F., Fay, M. F., and Lexer, C. (2007). Cross-species transfer of nuclear microsatellite markers: potential and limitations. Mol. Ecol. 16, 3759–3767. doi: 10.1111/j.1365-294X.2007.03439.x
Barboza, K., Beretta, V., Kozub, P. C., Salinas, C., Morgenfeld, M. M., Galmarini, C. R., et al. (2018). Microsatellite analysis and marker development in garlic: distribution in EST sequence, genetic diversity analysis, and marker transferability across Alliaceae. Mol. Genet. Genomics 293, 1091–1106. doi: 10.1007/s00438-018-1442-5
Bernardes, V., Murakami, D. M., Bizão, N., Souza, T. N., Silva, M. J., and Telles, M. P. C. (2021). Transferability and characterization of microsatellite markers from Byrsonima cydoniifolia A. Juss. (MALPIGHIACEAE) in seven related taxa from Cerrado biome reveal genetic relationships. Mol. Biol. Rep. 48, 4039–4046. doi: 10.1007/s11033-021-06411-z
Bharti, R., Kumar, S., and Parekh, M. J. (2018). Development of genomic simple sequence repeat (gSSR) markers in cumin and their application in diversity analyses and cross-transferability. Ind. Crops Prod. 111, 158–164. doi: 10.1016/j.indcrop.2017.10.018
Blank, A. F., Jesus, A. S., Santos, C. P., Grando, C., Pinheiro, J. B., Zucchi, M. I., et al. (2014). Development and characterization of novel microsatellite markers in Hyptis pectinate (Lamiaceae). Genet. Mol. Res. 13, 10173–10176. doi: 10.4238/2014.december.4.11
Bombonato, J. R., Bonatelli, I. A. S., Silva, R. G. A., Moraes, E. M., Zappi, D. C., Taylor, N. P., et al. (2019). Cross-genera SSR transferability in cacti revealed by a case study using Cereus (Cereeae, Cactaceae). Genet. Mol. Biol. 42, 87–94. doi: 10.1590/1678-4685-GMB-2017-0293
Botstein, D., White, R. L., Skolnick, M., and Davis, R. W. (1980). Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am. J. Hum. Genet. 32, 314–331.
Budak, H., Shearman, R. C., Parmaksiz, I., and Dweikat, I. (2004). Comparative analysis of seeded and vegetative biotype buffalograsses based on phylogenetic relationship using ISSRs, SSRs, RAPDs, and SRAPs. Theor. Appl. Genet. 109, 280–288. doi: 10.1007/s00122-004-1630-z
Dharajiya, D. T., Shah, A., Galvadiya, B. P., Patel, M. P., Srivastava, R., Pagi, N. K., et al. (2020). Genome-wide microsatellite markers in castor (Ricinus communis L.): identification, development, characterization, and transferability in Euphorbiaceae. Ind. Crops Prod. 151:112461. doi: 10.1016/j.indcrop.2020.112461
Dutta, S., Kumawat, G., Singh, B. P., Gupta, D. K., Singh, S., Dogra, V., et al. (2011). Development of genic-SSR markers by deep transcriptome sequencing in pigeonpea [Cajanus cajan (L.) Millspaugh]. BMC Plant Biol. 11:17. doi: 10.1186/1471-2229-11-17
Earl, D. A., and VonHoldt, B. M. (2012). STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour. 4, 359–361. doi: 10.1007/s12686-011-9548-7
Editorial committee of the Flora of China, and Chinese Academy of Sciences (1977). “Elsholtzia(Lamiaceae),” in Flora of China, Vol. 66, eds Z. Y. Wu and X. W. Li (Beijing: Science Press). 304–348.
Ellegren, H. (2000). Microsatellite mutations in the germline: implications for evolutionary inference. Trends Genet. 16, 551–558. doi: 10.1016/S0168-9525(00)02139-9
Ellis, J. R., and Burke, J. M. (2007). EST-SSRs as a resource for population genetic analyses. Heredity 99, 125–132. doi: 10.1038/sj.hdy.6801001
Engelhardt, K. A. M., Lloyd, M. W., and Neel, M. C. (2014). Effects of genetic diversity on conservation and restoration potential at individual, population, and regional scales. Biol. Conserv. 179, 6–16. doi: 10.1016/j.biocon.2014.08.011
Evanno, G., Regnaut, S., and Goudet, J. (2005). Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol. Ecol. 14, 2611–2620. doi: 10.1111/j.1365-294X.2005.02553.x
Excoffier, L., and Lischer, H. E. L. (2010). Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 10, 564–567. doi: 10.1111/j.1755-0998.2010.02847.x
Fan, D. M., Yue, J. P., Nie, Z. L., Li, Z. M., Comes, H. P., and Sun, H. (2013). Phylogeography of Sophora davidii (Leguminosae) across the ‘Tanaka-Kaiyong Line’, an important phytogeographic boundary in Southwest China. Mol. Ecol. 22, 4270–4288. doi: 10.1111/mec.12388
Gadissa, F., Tesfaye, K., Dagne, K., and Geleta, M. (2018). Genetic diversity and population structure analyses of Plectranthus edulis (Vatke) Agnew collections from diverse agro-ecologies in Ethiopia using newly developed EST-SSRs marker system. BMC Genet. 19:92. doi: 10.1186/s12863-018-0682-z
Gao, C. H., Ren, X. D., Mason, A. S., Li, J., Wang, W., Xiao, M. L., et al. (2013). Revisiting an important component of plant genomes: microsatellites. Funct. Plant Biol. 40, 645–661. doi: 10.1071/FP12325
García-Arias, F. L., Osorio-Guarín, J. A., and Zarantes, V. M. N. (2018). Association study reveals novel genes related to yield and quality of fruit in cape gooseberry (Physalis peruviana L.). Front. Plant Sci. 9:362. doi: 10.3389/fpls.2018.00362
Gepts, P. (2006). Plant genetic resources conservation and utilization. Crop Sci. 46, 2278–2292. doi: 10.2135/cropsci2006.03.0169gas
Ginestet, C. (2011). Ggplot2: elegant graphics for data analysis. J. R. Stat. Soc. A Stat. 174, 245–245. doi: 10.1111/j.1467-985X.2010.00676_9.x
Grover, A., and Sharma, P. C. (2016). Development and use of molecular markers: past and present. Crit. Rev. Biotechnol. 36, 290–302. doi: 10.3109/07388551.2014.959891
Guo, Z. Q., Liu, Z. Z., Wang, X. H., Liu, W. R., Jiang, R., Cheng, R. Y., et al. (2012). Elsholtzia: phytochemistry and biological activities. Chem. Cent. J. 6:147. doi: 10.1186/1752-153X-6-147
Ha, Y. J., Sa, K. J., and Lee, J. K. (2021). Identifying SSR markers associated with seed characteristics in Perilla (Perilla frutescens L.). Physiol. Mol. Biol. Plants 27, 93–105. doi: 10.1007/s12298-021-00933-3
Harijan, Y., Nishanth, G. K., Katageri, I. S., and Khadi, B. M. (2017). Transferability of heterologous SSR markers to cotton genotypes. Electron. J. Plant Breed. 8, 379–384. doi: 10.5958/0975-928X.2017.00057.6
Harley, R. M., Atkins, S., Budantsev, A. L., Cantino, P. D., Conn, B. J., Grayer, R., et al. (2004). “Labiatae,” in Families and Genera of Vascular Plants, Vol. 7, eds K. Kubitzki and J. W. Kadereit (Berlin: Springer). 167–275.
Herrera-Arroyo, M. L., Rico, Y., and Bedolla-García, B. Y. (2020). Morphotype divergence and genetic diversity of Hedeoma piperita Benth. in western Mexico. Mol. Biol. Rep. 47, 8925–8934. doi: 10.1007/s11033-020-05946-x
Holland, M. M., and Parson, W. (2011). GeneMarker® HID: a reliable software tool for the analyses of forensic STR data. J. Forensic Sci. 56, 29–35. doi: 10.1111/j.1556-4029.2010.01565.x
Huang, C. J., Chu, F. H., Huang, Y. S., Hung, Y. M., Tseng, Y. H., Pu, C. E., et al. (2020). Development and technical application of SSR-based individual identification system for Chamaecyparis taiwanensis against illegal logging convictions. Sci. Rep. 10:22095. doi: 10.1038/s41598-020-79061-z
Ihaka, R., and Gentleman, R. (1996). R: a language for data analyses and graphics. J. Comput. Graph. Stat. 5, 299–314. doi: 10.1080/10618600.1996.10474713
Jarne, P., and Lagoda, P. J. (1996). Microsatellites, from molecules to populations and back. Trends Ecol. Evol. 11, 424–429. doi: 10.1016/0169-5347(96)10049-5
Jiao, S. Q., Dong, A. X., Shi, T. L., Liu, H., Porth, I., Xin, H. B., et al. (2020). Development of a large gene-associated SSR marker set and in-depth genetic characterization in scarlet sage. Front. Genet. 11:504. doi: 10.3389/fgene.2020.00504
Jin, L., Macaubas, C., Hallmayer, J., Kimura, A., and Mignot, E. (1996). Mutation rate varies among alleles at a microsatellite locus: phylogenetic evidence. Proc. Natl. Acad. Sci. U.S.A. 93, 15285–15288. doi: 10.1073/pnas.93.26.15285
Jose, A. G., Concepción, M., and Francisco, M. (2018). Trends in plant research using molecular markers. Planta 247, 543–557. doi: 10.1007/s00425-017-2829-y
Ju, M. M., Fu, Y., Zhao, G. F., He, C. Z., Li, Z. H., and Tian, B. (2018). Effects of the Tanaka Line on the genetic structure of Bombax ceiba (Malvaceae) in dry-hot valley areas of southwest China. Ecol. Evol. 8, 3599–3608. doi: 10.1002/ece3.3888
Kalia, R. K., Chhajer, S., and Pathak, R. (2020). Cross genera transferability of microsatellite markers from other members of family Bignoniaceae to Tecomella undulata (Sm.) Seem. Acta Physiol. Plant. 42:151. doi: 10.1007/s11738-020-03138-5
Kalia, R. K., Rai, M. K., Kalia, S., Singh, R., and Dhawan, A. K. (2011). Microsatellite markers: an overview of the recent progress in plants. Euphytica 177, 309–334. doi: 10.1007/s10681-010-0286-9
Karaca, M., Ince, A. G., Aydin, A., and Ay, S. T. (2013). Cross-genera transferable e-microsatellite markers for 12 genera of the Lamiaceae family. J. Sci. Food Agric. 93, 1869–1879. doi: 10.1002/jsfa.5982
Krak, K., Vit, P., Douda, J., and Mandak, B. (2020). Development of 18 microsatellite markers for Salvia pratensis. Appl. Plant Sci. 8, e11316. doi: 10.1002/aps3.11316
Kwon, S. J., Lee, J. K., Kim, N. S., Yu, J. W., Dixit, A., Cho, E. G., et al. (2005). Isolation and characterization of microsatellite markers in Perilla frutescens Brit. Mol. Ecol. Notes 5, 455–457. doi: 10.1111/j.1471-8286.2005.00901.x
Lepcha, P., and Sathyanarayana, N. (2021). Variability for seed-based economic traits and genetic diversity analysis in Mucuna pruriens population of northeast India. Agric. Res. 1–11. doi: 10.1007/s40003-021-00568-6
Lesser, M. R., Parchman, T. L., and Buerkle, C. A. (2012). Cross-species transferability of SSR loci developed from transciptome sequencing in lodgepole pine. Mol. Ecol. Resour. 12, 448–455. doi: 10.1111/j.1755-0998.2011.03102.x
Letunic, I., and Bork, P. (2007). Interactive Tree of Life (iTOL): an online tool for phylogenetic tree display and annotation. Bioinformatics 23, 127–128. doi: 10.1093/bioinformatics/btl529
Leyva-Lopez, N., Nair, V., Bang, W. Y., Cisneros-Zevallos, L., and Heredia, J. B. (2016). Protective role of terpenes and polyphenols from three species of Oregano (Lippia graveolens, Lippia palmeri and Hedeoma patens) on the suppression of lipopolysaccharide-induced inflammation in RAW 264.7 macrophage cells. J. Ethnopharmacol. 187, 302–312. doi: 10.1016/j.jep.2016.04.051
Li, C. H., Zheng, Y. Q., Liu, Y., Lin, F. R., and Huang, P. (2021). Development of genomic SSR for the subtropical hardwood tree Dalbergia hupeana and assessment of their transferability to other related species. Forests 12, 804–804. doi: 10.3390/f12060804
Li, K. Q., Chen, L., Feng, Y. H., Yao, J. X., Li, B., Xu, M., et al. (2014). High genetic diversity but limited gene flow among remnant and fragmented natural populations of Liriodendron chinense Sarg. Biochem. Syst. Ecol. 54, 230–236. doi: 10.1016/j.bse.2014.01.019
Li, P., Qi, Z. C., Liu, L. X., Ohi-Toma, T., Lee, J., Hsieh, T. H., et al. (2017). Molecular phylogenetics and biogeography of the mint tribe Elsholtzieae (Nepetoideae, Lamiaceae), with an emphasis on its diversification in East Asia. Sci. Rep. 7:2057. doi: 10.1038/s41598-017-02157-6
Liang, J. Y., Ning, A. Q., Lu, P. Y., An, Y., Wang, Z. L., Zhang, J., et al. (2021). Biological activities and synergistic effects of Elsholtziastauntonii essential oil from flowers and leaves and their major constituents against Tribolium castaneum. Eur. Food Res. Technol. 247, 2609–2619. doi: 10.1007/s00217-021-03829-4
Lim, S. E., Sa, K. J., and Lee, J. K. (2021). Bulk segregant analysis identifies SSR markers associated with leaf and seed related traits in Perilla crop (Perilla frutescens L.). Genes Genom. 43, 323–332. doi: 10.1007/s13258-021-01056-5
Liu, K., and Muse, S. V. (2005). PowerMarker: an integrated analyses environment for genetic marker analyses. Bioinformatics 21, 2128–2129. doi: 10.1093/bioinformatics/bti282
Liu, Y. L., Geng, Y. P., Song, M. L., Zhang, P. F., Hou, J. L., and Wang, W. Q. (2019). Genetic structure and diversity of glycyrrhiza populations based on transcriptome SSR markers. Plant Mole. Biol. Rep. 37, 401–412. doi: 10.1007/s11105-019-01165-2
Lv, J. H., Su, X. H., and Zhong, J. J. (2012). Fumigant activity of Elsholtzia stauntonii extract against Lasioderma serricorne. S. Afr. J. Sci. 108, 77–79. doi: 10.4102/sajs.v108i7/8.556
Ma, S. J., Sa, K. J., Hong, T. K., and Lee, J. K. (2017). Genetic diversity and population structure analysis in Perilla frutescens from Northern areas of China based on simple sequence repeats. Genet. Mol. Res. 16:16039746. doi: 10.4238/gmr16039746
Mangaravite, E., Vinson, C. C., Rody, H. V. S., Carniello, M. A., Silva, R. S., and Oliveira, L. O. (2016). Contemporary patterns of genetic diversity of Cedrela fissilis offer insight into the shaping of seasonal forests in eastern South America. Am. J. Bot. 103, 307–316. doi: 10.3732/ajb.1500370
Metzgar, D., Bytof, J., and Wills, C. (2000). Selection against frameshift mutations limits microsatellite expansion in coding DNA. Genome Res. 10, 72–80. doi: 10.1101/gr.10.1.72
Miranda, K. M. C., Guimarães, R. A., Silva, M. J., Oliveira, P. R. O., Ribeiro, T. G., Mendes, T. P., et al. (2020). Cross-amplification and characterization of microsatellite markers in species of Manihot Mill. (Euphorbiaceae) endemic to the Brazilian Cerrado. Acta Bot. Bras. 34, 772–777. doi: 10.1590/0102-33062019abb0374
Moreira, X., Abdala-Roberts, L., Rasmann, S., Castagneyrol, B., and Mooney, K. A. (2016). Plant diversity effects on insect herbivores and their natural enemies: current thinking, recent findings, and future directions. Curr. Opin. Insect Sci. 14, 1–7. doi: 10.1016/j.cois.2015.10.003
Nadeem, M. A., Nawaz, M. A., Shahid, M. Q., Dogan, Y., Comertpay, G., Yildiz, M., et al. (2018). DNA molecular markers in plant breeding: current status and recent advancements in genomic selection and genome editing. Biotechnol. Biotechnol. Equip. 32, 261–285. doi: 10.1080/13102818.2017.1400401
Nei, M., Tajima, F., and Tateno, Y. (1983). Accuracy of estimated phylogenetic trees from molecular data. J. Mol. Evol. 19, 153–170. doi: 10.1007/BF01840887
Ni, Z. X., Bai, T. D., Cai, H., Chen, S. F., and Xu, L. (2015). The Transferability of Pinus massoniana SSR in other Pinus species. Mol. Plant Breed. 13, 2811–2817. doi: 10.13271/j.mpb.013.002811
Oliveira, E. J., Pádua, J. G., Zucchi, M. I., Vencovsky, R., and Vieira, M. L. C. (2006). Origin, evolution and genome distribution of microsatellites. Genet. Mol. Biol. 29, 294–307. doi: 10.1590/S1415-47572006000200018
Pan, L., Huang, T., Yang, Z., Tang, L., Cheng, Y., Wang, J., et al. (2018). EST-SSR marker characterization based on RNA-sequencing of Lolium multiflorum and cross transferability to related species. Mol. Breed. 38:80. doi: 10.1007/s11032-018-0775-4
Park, D. H., Sa, K. J., Lim, S. E., Ma, S. J., and Lee, J. K. (2019). Genetic diversity and population structure of Perilla frutescens collected from Korea and China based on simple sequence repeats (SSRs). Genes Genom. 41, 1329–1340. doi: 10.1007/s13258-019-00860-4
Park, S. D. E. (2001). Trypanotolerance in West African Cattle and the Population Genetic Effects of Selection. Dissertation’s thesis. Dublin: University of Dublin.
Peakall, R., Gilmore, S., Keys, W., Morgante, M., and Rafalski, A. (1998). Cross-species amplification of soybean (Glycine max) simple sequence repeats (SSRs) within the genus and other legume genera: implications for the transferability of SSRs in plants. Mol. Biol. Evol. 15, 1275–1287. doi: 10.1093/oxfordjournals.molbev.a025856
Peakall, R., and Smouse, P. E. (2012). GenAlEx 6.5: genetic analysis in Excel. population genetic software for teaching and research—an update. Bioinformatics 28, 2537–2539. doi: 10.1093/bioinformatics/bts460
Pebesma, E. (2018). Simple features for R: standardized support for spatial vector data. R J. 10, 439–446. doi: 10.32614/RJ-2018-009
Pekhova, O. A., Timasheva, L. A., Danilova, I. L., and Belova, I. V. (2020). Accumulation of biologically active substances in plants of Elsholtziastauntonii Benth. grown in the foothill zone of the Crimea. Agrar. Bullet. Urals. 11, 76–84. doi: 10.32417/1997-4868-2020-202-11-76-84
Peng, J. Y., Shi, C., Wang, D. W., Li, S. Z., Zhao, X. L., Duan, A. N., et al. (2021). Genetic diversity and population structure of the medicinal plant Docynia delavayi (Franch.) Schneid revealed by transcriptome-based SSR markers. J. Appl. Res. Med. Aromat. Plants 21:100294. doi: 10.1016/j.jarmap.2021.100294
Pern, Y. C., Lee, S. Y., Ng, W. L., and Mohamed, R. (2020). Cross-amplification of microsatellite markers across agarwood-producing species of the Aquilarieae tribe (Thymelaeaceae). 3 Biotech 10:103. doi: 10.1007/s13205-020-2072-2
Pinto, M. V., Poornima, H. S., Sivaprasad, V., and Naik, V. G. (2018). A new set of mulberry-specific SSR markers for application in cultivar identification and DUS testing. J. Genet. 97, e31–e37. doi: 10.1007/s12041-018-0900-5
Pistelli, L., Najar, B., Giovanelli, S., Lorenzini, L., Tavarini, S., and Angelini, L. G. (2017). Agronomic and phytochemical evaluation of lavandin and lavender cultivars cultivated in the Tyrrhenian area of Tuscany (Italy). Ind. Crops Prod. 109, 37–44. doi: 10.1016/j.indcrop.2017.07.041
Poljuha, D., Franić, M., Kralj, I., Weber, T., Šatović, Z., Ban, D., et al. (2021). Genetic diversity and structure analysis of Croatian garlic collection assessed by SSR markers. Folia Hortic. 33, 157–171. doi: 10.2478/fhort-2021-0011
Powell, W., Morgante, M., Andre, C., Hanafey, M., Vogel, J., Tingey, S., et al. (1996). The comparison of RFLP, RAPD, AFLP and SSR (microsatellite) markers for germplasm analysis. Mol. Breed. 2, 225–238. doi: 10.1016/j.scienta.2011.06.015
Pritchard, J. K., Stephens, M., and Donnelly, P. (2000). Inference of population structure using multilocus genotype data. Genetics 155, 945–959. doi: 10.1093/genetics/155.2.945
Qi, L. J., Long, P., Jiang, C., Yuan, Y., and Huang, L. Q. (2015). Development of microsatellites and genetic diversity analysis of Scutellaria baicalensis Georgi using genomic-SSR markers. Acta Pharm. Sin. 50, 500–505. doi: 10.16438/j.0513-4870.2015.04.002
Radosavljević, I., Jakse, J., Javornik, B., Satovic, Z., and Liber, Z. (2011). New microsatellite markers for Salvia officinalis (Lamiaceae) and cross-amplification in closely related species. Am. J. Bot. 98, e316–e318. doi: 10.3732/ajb.1000462
Radosavljević, I., Satovic, Z., Jakse, J., Javornik, B., Greguraš, D., Jug-Dujaković, M., et al. (2012). Development of new microsatellite markers for Salvia officinalis L. and its potential use in conservation-genetic studies of narrow sndemic Salvia brachyodon Vandas. Int. J. Mol. Sci. 13, 12082–12093. doi: 10.3390/ijms130912082
Radosavljević, I., Satovic, Z., and Liber, Z. (2015). Causes and consequences of contrasting genetic structure in sympatrically growing and closely related species. AoB Plants 7:106. doi: 10.1093/aobpla/plv106
Ren, P. (2012). The research on genetic diversity of Leonurus japonicus germplasm resources in Henan province. [master’s thesis]. Zhengzhou (HN): Henan Agricultural University.
Sa, K. J., Choi, I. Y., Park, K. C., and Lee, J. K. (2018). Genetic diversity and population structure among accessions of Perilla frutescens (L.) Britton in East Asia using new developed microsatellite markers. Genes Genom. 40, 1319–1329. doi: 10.1007/s13258-018-0727-8
Sa, K. J., Choi, I. Y., Ueno, M., and Lee, J. K. (2015). Genetic diversity and population structure in cultivated and weedy types of Perilla in East Asia and other countries as revealed by SSR marker. Hortic. Environ. Biotechnol. 56, 524–534. doi: 10.1007/s13580-015-0039-8
Sa, K. J., Lim, S. E., Choi, I. Y., Park, K. C., and Lee, J. K. (2019). Development and characterization of new microsatellite markers for Perilla frutescens (L.) Britton. Am. J. Plant Sci. 10, 1623–1630. doi: 10.4236/ajps.2019.109115
Saeed, A. F., Wang, R. Z., and Wang, S. H. (2016). Microsatellites in pursuit of microbial genome evolution. Front. Microbiol. 6:1462. doi: 10.3389/fmicb.2015.01462
Sandes, S. S., Pinheiro, J. B., Zucchi, M. I., Monteiro, M., Arrigoni-Blank, M. F., and Blank, A. F. (2013). Development and characterization of microsatellite primers in Pogostemon cablin (Lamiaceae). Genet. Mol. Res. 12, 2837–2840. doi: 10.4238/2013.August.8.4
Satya, P., Paswan, P. K., Ghosh, S., Majumdar, S., and Ali, N. (2016). Confamiliar transferability of simple sequence repeat (SSR) markers from cotton (Gossypium hirsutum L.) and jute (Corchorus olitorius L.) to twenty two Malvaceous species. 3 Biotech 6:65. doi: 10.1007/s13205-016-0392-z
Savić, A., Pipan, B., Vasić, M., and Meglič, V. (2021). Genetic diversity of common bean (Phaseolus vulgaris L.) germplasm from Serbia, as revealed by single sequence repeats (SSR). Sci. Hortic. 288:110405. doi: 10.1016/j.scienta.2021.110405
Schneider, S. J., Hwang, A. Y., Land, S. D., Chen, L. L., Thomas, A. L., and Hwang, C. F. (2019). Genetic diversity of ten black walnut (Juglans nigra L.) cultivars and construction of a mapping population. Tree Genet. Genomes 15:62. doi: 10.1007/s11295-019-1369-y
Segarra-Moragues, J. G., and Gleiser, G. (2009). Isolation and characterisation of di and tri nucleotide microsatellite loci in Rosmarinus officinalis (Lamiaceae), using enriched genomic libraries. Conserv. Genet. 10, 571–575. doi: 10.1007/s10592-008-9572-7
Segarra-Moragues, J. G., Marco, Y. C., Castellanos, M. C., Molina, M. J., and García-Fayos, P. (2015). Ecological and historical determinants of population genetic structure and diversity in the Mediterranean shrub Rosmarinus officinalis (Lamiaceae). Bot. J. Linn. Soc. 180, 50–63. doi: 10.1111/boj.12353
Serrote, C. M. L., Reiniger, L. R. S., Silva, K. B., Rabaiolli, S. M. D., and Stefanel, C. M. (2020). Determining the polymorphism information content of a molecular marker. Gene 726:144175. doi: 10.1016/j.gene.2019.144175
Shivakumar, M. S., Ramesh, S., Mohan Rao, A., Udaykumar, H. R., and Keerthi, C. M. (2017). Cross legume species/genera transferability of SSR markers and their utility in assessing transferability of SSR markers and their utility in assessing polymorphism among advanced breeding lines in Dolichos bean (Lablab purpureus L.). Int. J. Curr. Microbiol. App. Sci. 6, 656–668. doi: 10.20546/ijcmas.2017.608.083
Si, S. X., and Peng, Z. F. (2017). Introduction and domestication and application of Elsholtzia stauntonii in the middle and lower reaches of the Yellow River Plain. J. Anhui Agri. Sci. 45:214. doi: 10.13989/j.cnki.0517-6611.2017.04.053
Sulu, G., Kacar, Y. A., Polat, I., Kitapci, A., Turgutoglu, E., Simsek, O., et al. (2020). Identification of genetic diversity among mutant lemon and mandarin varieties using different molecular markers. Turk. J. Agric. For. 44, 465–478. doi: 10.3906/tar-1909-67
Tadesse, M., Kebede, M., and Girma, D. (2021). Genetic diversity of Tef [Eragrostis tef (Zucc.) Trotter] as revealed by microsatellite markers. Int. J. Genomics 2021:6672397. doi: 10.1155/2021/6672397
Tian, Y. L., Li, J. Q., Wang, W. H., Shi, A. P., and Wang, H. L. (2007). The germplasm resource and ornamental character of the Lamiaceae in Beijing. J. Beijing Univ. Agric. 22, 41–43. doi: 10.13473/j.cnki.issn.1002-3186.2007.03.006
Vanavermaete, D., Fostier, J., Maenhout, S., and De Baets, B. (2021). Deep scoping: a breeding strategy to preserve, reintroduce and exploit genetic variation. Theor. Appl. Genet. 134, 3845–3861. doi: 10.1007/s00122-021-03932-w
Vieira, M. L. C., Santini, L., Diniz, A. L., and Munhoz, C. D. F. (2016). Microsatellite markers: what they mean and why they are so useful. Genet. Mol. Biol. 39, 312–328. doi: 10.1590/1678-4685-GMB-2016-0027
Vining, K. J., Pandelova, I., Hummer, K., Bassil, N., Contreras, R., Neill, K., et al. (2019). Genetic diversity survey of Mentha aquatica L. and Mentha suaveolens Ehrh., mint crop ancestors. Genet. Resour. Crop Evol. 66, 825–845. doi: 10.1007/s10722-019-00750-4
Wang, J., and Shete, S. (2012). Testing departure from Hardy-Weinberg proportions. Methods Mol. Biol. 850, 77–102. doi: 10.1007/978-1-61779-555-8_6
Wang, S. Q. (2020). Genetic diversity and population structure of the endangered species Paeonia decomposita endemic to China and implications for its conservation. BMC Plant Biol. 20:510. doi: 10.1186/s12870-020-02682-z
Wang, X. P., Wen, H., Shang, Z. W., Yang, S., Xu, J., and Shen, Q. (2017). Genetic relationship and genetic diversity of 13 Perilla frutescens varieties based on SSR markers. Guizhou Agri. Sci. 45, 103–106.
Xing, X. Q., Li, M. D., Gong, Y. M., Niu, B., Hao, N., Tian, Y. L., et al. (2019). Comparative analysis of volatile components in essential oil of Elsholtziastauntonii and different organs. South Chin. Agri. 13, 150–155. doi: 10.19415/j.cnki.1673-890x.2019.02.077
Yang, Y. Y., He, R. Q., Zheng, J., Hu, Z. H., Wu, J., and Leng, P. S. (2020). Development of EST-SSR markers and association mapping with floral traits in Syringa oblata. BMC Plant Biol. 20:436. doi: 10.1186/s12870-020-02652-5
Young, A., Boyle, T., and Brown, T. (1996). The population genetic consequences of habitat fragmentation for plants. Trends Ecol. Evol. 11, 413–418. doi: 10.1016/0169-5347(96)10045-8
Zeng, Y. F., Zhang, J. G., Abuduhamiti, B., Wang, W. T., and Jia, Z. Q. (2018). Phylogeographic patterns of the desert poplar in Northwest China shaped by both geology and climatic oscillations. BMC Ecol. Evol. 18:75. doi: 10.1186/s12862-018-1194-1
Zhang, W. J., Yuan, Q. J., Jiang, D., Zhang, Y. Q., and Huang, L. Q. (2014). Development and characterisation of microsatellite markers for the medicinal plant Scutellaria baicalensis (Lamiaceae). Biochem. Syst. Ecol. 54, 267–271. doi: 10.1016/j.bse.2014.01.005
Zhang, Z. Y., Cui, B. B., Mao, J. F., Pang, X. M., and Li, Y. Y. (2015). Novel polymorphic EST-derived microsatellite markers for the red-listed five needle pine, Pinus dabeshanensis. Conserv. Genet. Resour. 7, 191–192. doi: 10.1007/s12686-014-0329-y
Zhao, Y. J., and Gong, X. (2015). Genetic divergence and phylogeographic history of two closely related species (Leucomeris decora and Nouelia insignis) across the’ Tanaka Line’ in southwest China. BMC Evol. Biol. 15:134. doi: 10.1186/s12862-015-0374-5
Zheng, S., Kang, S., Shen, Y., and Sun, L. (1999). Three new c-methylated flavones from Elsholtziastauntonii. Planta Med. 65, 173–175. doi: 10.1055/s-2006-960459
Zhou, X. M., Li, B. Y., and Wang, Y. J. (2014). Cultivation and management of Elsholtziastauntonii. Chin. Flowers Hortic. 24, 40–41.
Zhou, X. M., Liu, H. C., and Li, B. Y. (2016). Type and distribution of trichomes on the leaf epidermis of Elsholtzia stautonii and the secreting process of the glandular trichomes. Acta Hortic. Sin. 43, 1555–1565. doi: 10.16420/j.issn.0513-353x.2016-0154
Keywords: cross-transferability, Elsholtzia stauntonii Benth., simple sequence repeat (SSR) markers, genetic diversity, population structure
Citation: Zhang C, Jia C, Liu X, Zhao H, Hou L, Li M, Cui B and Li Y (2022) Genetic Diversity Study on Geographical Populations of the Multipurpose Species Elsholtzia stauntonii Using Transferable Microsatellite Markers. Front. Plant Sci. 13:903674. doi: 10.3389/fpls.2022.903674
Received: 24 March 2022; Accepted: 21 April 2022;
Published: 12 May 2022.
Edited by:
Shunli Wang, Chinese Academy of Agricultural Sciences (CAAS), ChinaReviewed by:
Sezai Ercisli, Atatürk University, TurkeyZhang Qingyu, Northwest A&F University, China
Copyright © 2022 Zhang, Jia, Liu, Zhao, Hou, Li, Cui and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yingyue Li, yingyueli@bjfu.edu.cn; Binbin Cui, cbb0508@163.com