Next Article in Journal
The Intestinal Microbiome and the Metabolic Syndrome—How Its Manipulation May Affect Metabolic-Associated Fatty Liver Disease (MAFLD)
Previous Article in Journal
Pathological Effects of SARS-CoV-2 Associated with Hematological Abnormalities
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
CIMB
Volume 45
Issue 9
10.3390/cimb45090454
cimb-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Six Express Sequence Tag–Simple Sequence Repeat Primers Reveal Genetic Diversity in the Cultivars of Three Zanthoxylum Species

by
Yangchuan Deng
1,
Zhoujian He
1,
Yanlin Li
1,
Meng Ye
1,*,† and
Li Xiang
2,*,†
1
College of Forestry, Sichuan Agricultural University, Huimin Road 211, Wenjiang District, Chengdu 611130, China
2
Institute of Chinese Materia Medica, China Academy of Chinese Medical Sciences, No. 16 Nanxiaojie, Dongzhimennei Ave., Beijing 100700, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Curr. Issues Mol. Biol. 2023, 45(9), 7183-7196; https://doi.org/10.3390/cimb45090454
Submission received: 28 July 2023 / Revised: 21 August 2023 / Accepted: 24 August 2023 / Published: 30 August 2023
(This article belongs to the Section Molecular Plant Sciences)
Browse Figures
Versions Notes

Abstract

:
Zanthoxylum (Sichuan pepper), with its rich cultivars, has long been widely cultivated in China for its unique seasoning and medicinal uses, but most of its cultivars have similar morphological characteristics. Therefore, we hypothesized that the genetic diversity of Zanthoxylum cultivars is low because of their apomixis and long cultivation history. In this study, we aimed to investigate the genetic diversity of three Zanthoxylum species on the cultivar level based on express sequence tag–simple sequence repeat (EST-SSR) primers. In total, 121 samples of three Zanthoxylum species (Z. bungeanum, Z. armatum and Z. piperitum) were collected from different areas in China for genetic diversity analysis. A total of six specificity and polymorphism EST-SSR primers, which we selected from among 120 primers based on two transcriptomes (Z. bungeanum, Z. armatum) in our earlier study, were used to evaluate genetic diversity based on capillary electrophoresis technology. The results of our analysis using the unweighted pair group method with arithmetic mean (UPGMA) indicated that most of the samples are clustered in one clade in the UPGMA dendrogram, and the average genetic distance was 0.6409. Principal component analysis (PCA) showed that Z. piperitum may have a closer genetic relationship with Z. bungeanum than with Z. armatum. An analysis of molecular variation (AMOVA) showed that the genetic variation mainly stemmed from individuals within populations; the genetic differentiation coefficient (PhiPT) was 0.429, the gene flow (Nm) between populations was 0.333, and the differences among populations were not significant (p > 0.001). For the intraspecific populations of ZB, the percentage of genetic variation was 53% among populations and 47% within populations, with non-significant differences between populations (p > 0.001). The genetic differentiation coefficient (PhiT) was 0.529, and the gene flow (Nm) was 0.223. For the intraspecific populations of ZA, the results indicated that the percentage of genetic variation was 29% among populations and 71% within populations, with non-significant differences between populations (p > 0.001); the genetic differentiation coefficient (PhiPT) was 0.293, and the gene flow (Nm) was 0.223. Through genetic structure analysis (GSA), we predicted that these 121 samples belonged to two optimal subgroups, which means that all the samples probably originated from two gene pools. Above all, this indicated that the genetic diversity of the 121 Zanthoxylum samples was relatively low at both the species and cultivar levels, a finding which was consistent with our initial assumptions. This study provides a reference, with molecular-level data, for the further identification of Zanthoxylum species.

1. Introduction

Sichuan pepper belongs to the genus Zanthoxylum (family Rutaceae), which includes a wide variety of shrubs and trees. To date, more than 250 species of Zanthoxylum have been identified, most of which are distributed in tropical and subtropical Oceania, Africa, and Asia [1,2,3]. In China, 41 species belonging to the genus Zanthoxylum have been identified, 25 of which are endemic [4]. In addition, several major species have been cultivated for their aromatic, edible, and medicinal properties [5,6,7,8]. Cultivars of Zanthoxylum species are currently being cultivated in China due to their economic and scientific value. Among these species, Zanthoxylum bungeanum Maxim. and Zanthoxylum armatum DC. are the main species of Sichuan pepper in mainland China [9]. However, an increasing number of cultivars of Z. bungeanum and Z. armatum are being domesticated in China. Some species and cultivars can be easily confused because of their similar morphological characteristics [10]. For instance, “Qinghuajiao” in trade is often regarded as Z. schinifolium, but according to the research, the majority of “Qinghuajiao” products are mostly Z. armatum [9,11]. Likewise, “Tengjiao”, “Jinyang Qinghuajiao” and “Jiuye Qinghuajiao” are all cultivars of Z. armatum, although they are generally treated as different species [12]. For most of the plants, their isolated geographical environments and different climates usually lead to consolidated differences between taxa. However, most Zanthoxylum species exhibit apomixis, in which differences between cultivars may be negligible. Therefore, we hypothesized that the genetic diversity of Zanthoxylum species is low on the cultivar level because of their apomixis and stable genotypes, although they have many different names and are distributed in several different areas [13,14].
Unfortunately, the published literature on the genetic diversity of Zanthoxylum species and the relationship between their morphological traits and genetic diversity is lacking. The existing studies have mainly focused on traditional uses such as medicine [15,16,17], phytochemical constituents [18,19,20] and antioxidant activities [21,22,23]. However, there are few studies on genetic diversity. Medhi et al. (2014) utilized random amplified polymorphic DNA (RAPD) and inter-simple sequence repeat (ISSR) analysis to explore the genetic diversity of three Zanthoxylum species, and the results revealed the presence of a significant variability within and between the different Zanthoxylum species [24]. Sun et al. (2010) investigated Z. schinifolium ecotypes using the common primers of the ITS region [25]. Feng et al. (2020) investigated the genetic diversity, population structure and evolutionary history of Zanthoxylum, and the results indicated that the Z. bungeanum found in China most likely originated from southeastern Gansu Province [26]. Kim et al. (2017) developed microsatellite markers as tools for the investigation of Z. schinifolium. They found 15 polymorphic and 3 monomorphic microsatellite markers for the genetic study of Z. schinifolium. Furthermore, 3 of the 11 cross-amplified primers could distinguish between two Zanthoxylum species [27]. Su et al. (2022) investigated the developmental mechanism of Z. bungeanum prickles via morphological and multiomic analyses. The results showed that nine differentially expressed genes related to prickle development were screened and validated [28]. Hence, it is highly important to develop a reliable method for detecting the genetic diversity of Zanthoxylum species and cultivars.
The EST-SSR marker with conserved sequences has been utilized to design primers with a combination of SSR and transcriptome techniques. This method can be used to conveniently and rapidly screen samples by relying on more accurate and specific primers for genetic diversity, classification and fingerprinting [29]. Moreover, this approach has previously been used for various species, such as wheat [30], grape [31], rice [32], soybean [33] and orange [34]. Such studies on Zanthoxyum have also been reported, in which similar techniques based on the transcriptome were applied, but with different aims [35,36,37]. Therefore, we sequenced transcriptomes of both Z. bungeanum and Z. armatum, hoping to screen EST-SSR primers for a genetic diversity analysis of the Zanthoxylum cultivars examined in our previously published article [38]. A total of 36.76 G high-quality clean data were screened for the subsequent analysis, while 64,944 and 75,669 unigenes were obtained from Z. bungeanum and Z. armatum, respectively. In total, 12,746 SSR loci were identified in 10,595 unigenes of Z. bungeanum, and 15,096 SSR loci were identified in 12,612 unigenes of Z. armatum. A total of 60 pairs of EST-SSR primers were randomly selected from both Z. armatum and Z. bungeanum, respectively. Six pairs showed significant specificity and polymorphism in 12 samples [38]. In this study, these six selected EST-SSR primers were used to investigate the genetic diversity of 121 Zanthoxylum cultivars based on the UPGMA, PCA, AMOVA and GSA. We aimed to reveal the genetic diversity among these Zanthoxylum cultivars through the EST-SSR method.

2. Materials and Methods

2.1. Materials

The 121 Zanthoxylum samples (dried pericarp or fresh leaves), which were collected from different provinces in China from 2015 to 2017, were tested for genetic diversity (Figure 1). These 121 samples included 67 Z. bungeanum (ZB) samples, 47 Z. armatum (ZA) samples and 7 Zanthoxylum piperitum De Candolle (ZP) samples [39] (Supplementary Table S1). All samples were authenticated by Professor Meng Ye, Professor of College of Forestry, Sichuan Agricultural University. In particular, samples ZB-HCDHP (Irr)-1, ZB-HCDHP (Irr)-2, and ZB-HCDHP (Irr)-3 (Accession: 9–11) were irradiated; sample ZB-X (MF) (Accession: 67) only had male flowers. All samples were selected at random with three replicates, frozen in liquid nitrogen and stored at −80 °C after being returned to the laboratory.

2.2. DNA Isolation, PCR Amplification and Capillary Electrophoresis Fluorescence

Tissues (100 mg) of the 121 Zanthoxylum samples were used to extract DNA in accordance with the protocol of the DNA Secure Plant Kit (Tiangen Biotech (Beijing) Co., Ltd., Beijing, China, DP320-02). The DNA was tested via 1% agarose gel electrophoresis under a voltage of 120 V for 15 min.
The 20 µL reaction mixtures used for PCR amplification (Bio-Rad C1000 Touch, Hercules, CA, USA) comprised the following components: 14.8 µL of ddH2O, 0.4 µL of dNTPs, 2 µL of buffer, 0.3 µL of forward primer (20 µM), 0.3 µL of reverse primer (20 µM), 0.2 µL of DNA and 0.2 µL of Taq polymerase. The PCR program used was as follows: predenaturation at 94 °C for 5 min; 35 cycles of denaturation at 94 °C for 30 s, renaturation at 55 °C for 35 s and extension at 72 °C for 40 s; and a final extension at 72 °C for 3 min, followed by storage at 4 °C.
After mixing the formamide with the molecular weight internal standard at a volume ratio of 100:1, 9 μL was taken and added to the upper sample plate, and then 1 μL of 10-fold-diluted PCR product was added. Then, capillary electrophoresis fluorescence was carried out using an ABI PRISM 3730XL sequencer.

2.3. Polymorphism Assessment of Primers

The polymorphism of the six pairs of primers, including Nei’s gene diversity (h) [40], Shannon’s information index (I) [41] and the percentage of polymorphic loci, was determined using POPGENE32 VERSION1.31 software. Nei’s standard genetic distance was determined using the following formula [40]:
J x = ∑ i = 1 r ∑ j = 1 k X i j 2 / n ,   J y = ∑ i = 1 r ∑ j = 1 k Y i j 2 / n ,   J y x = ∑ i = 1 r ∑ j = 1 k X i j Y i j / n
X i j , Y i j : gene frequency of the k-th allele on locus r in populations X and Y.

2.4. Genetic Distance Calculation, UPGMA and PCA

The raw data obtained using the ABI PRISM 3730XL sequencer were analyzed using Fragment (Plant) analysis VERSION 1.81 software in Gene Marker 2.2, and the fragment sizes obtained were analyzed by comparing and analyzing the position of the molecular weight internal standard in each lane with the position of the peak of each sample.
After the fragment sizes were confirmed, all the fragments were converted into a matrix of “0” and “1”, representing no production and substantial production, respectively. Based on the matrix of “0” and “1”, the genetic distances were calculated using NTSYSpc 2.11 software. The unweighted pair group method with arithmetic mean analysis (UPGMA) was performed using MEGA7.0 software. Principal component analysis (PCA) was performed using Canoco 4.5 software.

2.5. Analysis of Molecular Variance (AMOVA) and Genetic Structure Analysis (GSA)

Based on the fragment sizes of the PCR products, GenAlex 6.5 software [42] was used to perform an analysis of molecular variance (AMOVA) in order to calculate the genetic variation among the species and the coefficient of genetic differentiation (PhiPT), thus revealing the genetic component of the population. For the entire individual genetic structure analysis (GSA), a genetic structure calculation was performed using STRUCTURE VERSION 2.3.4. based on the Bayesian algorithm method [43,44]. At the beginning of the procedure, the number of Markov chain Monte Carlo (MCMC) iterations was set to 100,000 based on the admixture–mixture model, while it was set to 100,000 iterations during the burnin period and then 500,000 after the burnin period. Furthermore, the K-value was set at 1–3, and each K-value iteration was run 10 times. The STRUCTURE results were further processed using the Structure Harvester online tool and CLUMMP VERSION 1.1.2 software [45], and the final results were calculated and graphically presented using EXCEL 2019.

3. Results

3.1. Primers Showed Excellent Specificity and Polymorphism

In our previous study, six primer pairs showed significant specificity and polymorphism in 12 samples. Therefore, the same six pairs of EST-SSR primers were used in this trial to analyze 121 samples of Zanthoxylum species based on capillary electrophoresis technology. Two pairs of primers originated from Z. armatum, whereas the other four were derived from Z. bungeanum (Supplementary Table S2). We succeeded in amplifying 51 allelic loci, with a mean number of 8.5 polymorphic loci, using these primer pairs. The average number of amplified allelic genes (Na) was two, the average number of effective alleles (Ne) was 1.3323, the average Nei’s genetic diversity (h) was 0.2069, and the Shannon information index (I) was 0.3273. These values clearly indicated that the EST-SSR primers derived from the Zanthoxylum species themselves showed excellent specificity and polymorphism for the 121 samples used in this study. Among these primers, primer ZB26 (SSR sequence: 5′-(GGA)5(GGT)5)-3′) showed better characteristics. Specifically, its number of amplified allelic genes (Na) was two, its number of effective alleles (Ne) was 1.2738, the Nei’s genetic diversity (h) was 0.1608, the Shannon information index (I) was 0.2543, and the number of polymorphic loci was 13. Furthermore, polymorphisms of six primers for every single species were investigated. Specifically, for Z. bungeanum, the mean number of polymorphic loci was 6.8, the average number of amplified allelic genes (Na) was 1.8338, the average number of effective alleles (Ne) was 1.3280, the average Nei’s genetic diversity (h) was 0.1980, and the Shannon information index (I) was 0.3062. For Z. armatum, the mean number of polymorphic loci was 4.8, the average number of amplified allelic genes (Na) was 1.5733, the average number of effective alleles (Ne) was 1.2147, the average Nei’s genetic diversity (h) was 0.1305, and the Shannon information index (I) was 0.2059. For Z. piperitum, the average number of polymorphic loci was 3.7, the average number of amplified allelic genes (Na) was 1.3210, the average number of effective alleles (Ne) was 1.1630, the average Nei’s genetic diversity (h) was 0.0954, and the Shannon information index (I) was 0.1466.

3.2. Most Samples Are Clustered in Correct Clade in the UPGMA Dendrogram

As shown in Figure 2, on the basis of genetic distances and the UPGMA method, six pairs of primers were successfully used to cluster three species of the Zanthoxylum genus in this study, which also indicated that the EST-SSR primers originating from the Zanthoxylum genus themselves could distinguish between the Zanthoxylum species more clearly. In terms of clustering branches, the three species ZB, ZA and ZP were generally distinguished. Specifically, the whole samples were divided into three branches: one branch contained mostly ZB samples, one branch included all the ZA samples, and the final branch included all the ZP samples and coupled ZB samples [ZB-JLHJ, ZBYXGJ, ZBHYDHP(Sti), ZB-LWHJ(Sti) and ZB-HCDHP(LN)-1]. The genetic distances of the 121 samples used in this study ranged from 0 to 2.3811, with a mean of 0.6409. From an interspecies perspective, the average genetic distance between ZB and ZA was 0.7856, that between ZB and ZP was 0.8147, and that between ZA and ZP was 1.8888. Within a species, the average genetic distance was 0.5104 for all the ZB samples, 0.2333 for all the ZA samples, and 0.6481 for all the ZP samples. The genetic distances of some samples were even 0, especially for the ZA-TJ series. These data indicated that the genetic diversity of the Zanthoxylum samples in this study was not high, and the genetic differentiation of the intraspecific cultivars was relatively low, indicating that EST-SSR did not significantly discriminate between cultivars. All the ZP samples and several ZB samples (ZB-YXGJ, ZB-HCDHP(LN), ZB-LWHJ(STI)-1, ZB-HYDHP(STi) and ZB-(JLHJ)) were clustered close together, even within one branch, and the average genetic distance between the ZB and ZP samples was 0.8147, whereas that between the ZA and ZP samples was 1.1888. This indicated that ZP and ZB may have closer kinship with each other than with ZA, respectively.

3.3. The Distribution Areas of ZA Were Separated while ZP and ZB Were Well-Overlapped Based on PCA

As shown in Figure 3, the results revealed that the PCA method could clearly be used to investigate the relationships of ZA with ZB and ZP species. The distribution areas of ZP and ZB were relatively well-overlapped, suggesting that ZP and ZB may have a closer genetic relationship similar to that indicated by the results of the dendrogram. On the one hand, the ZB samples were more dispersed, and the ZA samples were more compactly distributed. On the other hand, the samples of both ZB and ZA were concentrated, and some samples overlapped, indicating that the intraspecific genetic diversity of both ZB and ZA was not high.

3.4. The Interspecific and Intraspecific Gene Flows Were Low

As shown in Table 1, the results of the AMOVA showed that the genetic variation among and within all the populations of Zanthoxylum in this study was 43% and 57%, respectively, and the differences among populations were not significant (p > 0.001). Furthermore, the genetic differentiation coefficient (PhiPT) was 0.429, and the gene flow between populations (Nm) was 0.333. For the intraspecific populations of ZB, the AMOVA showed that the percentage of genetic variation was 53% among populations and 47% within populations, with non-significant differences between populations (p > 0.001). The genetic differentiation coefficient (PhiT) was 0.529, and the gene flow (Nm) was 0.223. For the intraspecific populations of ZA, the results indicated that the percentage of genetic variation was 29% among populations and 71% within populations, with non-significant differences between populations (p > 0.001), that the genetic differentiation coefficient (PhiPT) was 0.293, and the gene flow (Nm) was 0.223. The percentage of genetic variation within populations was 71%, the difference between populations was highly insignificant (p > 0.001), the genetic differentiation coefficient (PhiPT) was 0.293, and the gene flow (Nm) was 0.605. Overall, the genetic variation of the Zanthoxylum populations in this study originated mainly within populations rather than among populations, and the interspecific gene flows and intraspecific gene flows were not high.

3.5. Three Zanthoxylum Species Divided into Two Genetic Structures

As shown in Figure 4, the GSA indicated that the optimal population of these 121 samples could be divided into two subgroups (the K-value was 2) (Figure 4a) and could more clearly distinguish between ZB, ZA and ZP. There was a greater similarity between ZP and ZB, which aligned with the results of the clustering analysis and PCA. In addition, some samples of ZB exhibited a degree of difference based on the GSA results, showing that the genetic diversity of some samples of ZB was higher than that of ZA. This may have been due to a higher number of samples and the distribution of sampling sites for ZB. In comparison, the similarity of the ZA samples was higher. There was less genetic differentiation within species, and the genetic diversity was lower. Among the seven ZP samples, the overall degree of difference was not significant, with only the ZP-PTSJ sample being slightly different from the other six samples. The other six samples consisted of three cultivars, with one sample of ZP-PTSJ, three samples of ZP-CCSJ and three samples of ZP-JLSJ. It is possible that the one sample of ZP-PTSJ differs from the three samples of the other two cultivars, but the presence of only one sample indicated that the findings could simply be due to an experimental error.

4. Discussion

4.1. EST-SSR Validation

Different kinds of molecular markers are now universally used to analyze genetic diversity for many species rapidly and efficiently. Examples include SSR in grape [46,47], chickpea [48], and Prunus armeniaca [49]; RAPD in grape [46], Apiaceae [50] and Piperaceae [51]; And CAPS in Citrus [52], cotton [53] and Camellia sinensis [54]. Although the species of Zanthoxylum are economically important trees in China [55], relatively few studies have focused on their genetic diversity [56]. In our published study, we sequenced the transcriptomes of Z. bungeanum and Z. armatum using the Illumina-Hiseq sequencing platform, which provided a good basis for the screening of EST-SSR markers and the detection of genetic diversity [38]. The benefit of obtaining EST-SSR primers using RNA-seq is that these markers have better specificity and polymorphism compared to conventional primers, because EST-SSR is derived from cDNA and is conserved [57]. These properties are considered ideal properties for molecular markers; thus, it is best to use EST-SSR for classification and identification. However, to date, only a few studies have detected the transcriptome of Zanthoxylum. Zhao et al. [58] reported the results of a germination analysis of Z. piperitum based on its transcriptome. Meanwhile, Feng et al. [37] employed a de novo transcriptome-developed SSR of Z. bungeanum. In addition, Fei et al. [14] investigated the mechanism of the apomixis of Z. bungeanum using an miRNA-based technique. However, the majority of previous studies only focused on one species, with no comparisons being performed. On the interspecies level, numerous published studies have proved that EST-SSRs have shown accurate and clear results with respect to the identity and genetic diversity of crop plants, such as Triticum [59], Hordeum vulgare [60] and Citrus [61]. On the cultivar level, this technique is more likely to fail to yield distinct results, as in the case of Zanthoxylum [62]. In this study, the six EST-SSR markers derived from two Zanthoxylum species showed significant availability and polymorphism for the three Zanthoxylum species samples, though it seemed that they were not able to definitively identify cultivars of Zanthoxylum clearly, suggesting that this technique may be more appropriate for species-level identification than cultivar-level identification. Additionally, the polymorphism for Z. piperitum was not as good as that for Z. bungeanum and Z. armatum. One possible reason is that there were no EST-SSR primers from Z. piperitum; another is that the number of samples was probably too low.

4.2. Genetic Diversity Assessments

In the present study, 121 samples were shown to belong to three species of the Zanthoxylum genus. Among the different analytical methods used in this study, the K-value analysis of the genetic structure, in particular, showed that the most suitable subgroup number for these samples was two. This indicated that the genetic diversity of the Zanthoxylum genus is not high. This finding may be attributed to several reasons. First, although the total number of samples in this study was large, many of the samples were cultivars. Second, only six primer pairs were used in this study, and while it was possible to clearly identify three Zanthoxylum species, the identification of intraspecific cultivars was relatively difficult, which indicated that the ability of the EST-SSR primers to identify intraspecific cultivars was somewhat limited. Third, the genus Zanthoxylum has a mainly reproductive mode of apomixis and thus has a stable phenotype [13,14]. Above all, the genetic differentiations were minimal, and the phenotypes were stable.
Additionally, we compared the morphological characteristics of Z. bungeanum, Z. armatum and Z. piperitum based on previously reported research (Supplementary Table S3) [9,63,64]. Z. piperitum is a special and stingless pepper that comes from Japan, which can be clearly differentiated from other Sichuan peppers [63,64]. There are few articles about Z. piperitum species regarding morphology and molecular biology, and there are no available EST databases. In this study, the ESR-SST primers for Z. bungeanum and Z. armatum were used to analyze the genetic diversity of Z. piperitum. The results indicated that Z. piperitum may be closer to Z. bungeanum than to Z. armatum. However, the lack of appropriate primers and the sample size may have had a certain impact on the results for Z. piperitum, and we will make improvements in future research.
Sample ZB-X (MF), a sample with only male flowers, did not show any significant differences from the rest of the samples in terms of clustering, PCA or genetic structure, probably because only one sample was analyzed. Three samples of ZB-HCDHP (Irr) were irradiated, and the clustering results revealed that these three samples did not show any significant differences when compared with the other samples. In addition, there was no genetic difference between the three samples. Furthermore, the PCA and GSA results revealed that these three samples showed no distinctive results. Irradiation mutation breeding is an approach used to induce gene mutations, promote gene recombination and improve the recombination rate [65,66,67]. Neither extraterrestrial radiation nor male flower production were found to play an important role in the genetic structure of these species, for which we could not detect differences based on the EST-SSR analysis. Since there were too few of the two specimen types, no statistically significant differences were identified between the individual samples. None of the ZB-DHP series samples, which were more widely distributed in the ZB branch, were particularly concentrated, suggesting the existence of some differences in ZB-DHP on the molecular level in different regions. In Figure 2, one of the clades included all the samples of ZP and several samples from ZB. Among them, ZB-JLHJ and ZB-YXGJY are both varieties with a long history of cultivation in China. However, because they have been cultivated in mountainous areas for a long time, there was less traffic, less foreign introduction and less exchange of genetic resources, so this may be the reason why they are not classified in the same branch as other ZB samples, or it may be the result of the chance of the experiments that and the results are not informative.
The ZA-Ws were wild Z. armatum samples, whereas the ZA-TJs were a Z. armatum cultivar. The genetic distance results and dendrogram demonstrated that these two types of samples were relatively genetically remote, indicating that the wild Z. armatum type and the cultivars used in this study differed to some extent. Meanwhile, the clustering, PCA and GSA results revealed that all the ZA-TJ samples showed an extremely close genetic relationship; the genetic distances of some samples were almost 0, indicating the very limited genetic differentiation of the cultivar. This may be because the ZA-TJ samples collected in this study were derived from Sichuan Province, where the sampling was performed within a small area. Thus, based on the genetic distance results, Zanthoxylum’s apomixis and the small sampling areas, it is likely that the samples may have originated from a single gene pool.
AMOVA is an important tool for studying genetic diversity among the populations of species. In this study, for the entire Zanthoxylum genus, interspecies variation accounted for 43% of the total variation, and intraspecies variation accounted for 57%, indicating that the population variation in the three Zanthoxylum samples in this study mainly originated from within the populations. Specifically, 53% of the inter-population variation was accounted for in ZB, suggesting that this variation mainly originated from among populations, which may be due to the fact that a larger number of ZB samples were collected in the present study and that the collection was performed using a wide range of regional sources. The genetic differentiation coefficient (PhiPT) and gene flow (Nm) are important indicators for analyzing the genetic structure of populations. According to Wright et al. (1990), if the genetic differentiation coefficient (PhiPT) ranges from 0 to 0.05, this indicates that no differentiation exists between populations; if it is 0.05–0.15, this indicates that there is moderate differentiation between populations; and if it is 0.15–0.25, this indicates high differentiation. If the PhiPT is greater than 0.25, this indicates extremely differentiated populations [68]. According to the results, it can be seen that the PhiPT was greater than 0.25, whether between the three species, within the ZB populations or within the ZA populations, which indicated that there was a great deal of genetic differentiation between the different Zanthoxylum species and also between the intraspecific populations in the current study. Ellstrand et al. suggested that the proportion of genetic variation among populations of asexual plants increases greatly when there are large barriers to gene flow among the populations [69]. Most of the Zanthoxylum species reproduce via apomixis, similar to asexual lineages; thus, it is reasonable to infer that the variation among populations is relatively large. According to Buso et al.’s research, a gene flow (Nm) > 1 results in a high degree of gene exchange and that of 0.250 < Nm < 0.999 results in a moderate degree of gene exchange [70]. None of the gene flow (Nm) values in this experiment were beyond 1; therefore, the degree of gene exchange was not high, probably because of the apomixis reproduction.
The Zanthoxylum market in China mostly consists of ZB and ZA species. Due to the different climatic conditions in China, Z. bungeanum is cultivated in the north and Z. armatum is cultivated in the south [9]. Therefore, as Z. bungeanum, Z. armatum and their cultivars are the main types of Sichuan pepper commercially grown in China, it is vital to make considerable efforts to maintain the diversity of Zanthoxylum with respect to characteristics such as its aroma, stinglessness and pungency.

5. Conclusions

In this research, we used six EST-SSR primers to detect the genetic diversity of three Zanthoxylum species, with 121 samples, using capillary electrophoresis technology. According to the UPGMA and PCA, the 121 samples were divided into three clades, and Z. piperitum and Z. bungeanum were more closely related to each other than either was to Z. armatum. On the basis of STRUCTURE, all the samples probably originated from two gene pools. Additionally, the AMOVA indicated that the differences among interspecies and intraspecies populations were not significant. Above all, the Zanthoxylum samples examined in this study showed relatively low genetic diversity at the cultivar level, which was consistent with our initial assumptions. However, we found that the EST-SSR primers can distinguish between the different species of Zanthoxylum but cannot clearly identify large numbers of cultivars; this topic requires further research. This study provides a reference, with molecular-level data, for further examination of the identity of Zanthoxylum.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cimb45090454/s1, Table S1: Table of 121 Zanthoxylum Samples information; Table S2: Table information of six pairs of EST-SSR primers and polymorphsim results; Table S3: Table of Morphology characteristics of Z. bungeanum, Z. armatum and Z. piperitum.

Author Contributions

M.Y. and L.X. conceived the experiments and helped in receiving funding. Y.D. performed the research and drafted the manuscript. Z.H. and Y.L. participated in the experiments. All authors have read and agreed to the published version of the manuscript.

Funding

Meng Ye: the Cultivation of Cultivars of Traditional Chinese Medicinal Herbs—Integrated Research and Industrialization Demonstration of Whole Industry Chain Technology of Z. armatum (No. 2018SZ0064); Li Xiang: the National Natural Science Foundation of China (U1812403-1); and the CACMS Innovation Fund (CI2021A05103).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are openly available with the article and Supplementary Materials.

Acknowledgments

We want to express our gratitude to Kunbin Cao (Sichuan Academy of Forestry), whose significant contributions greatly enriched this research. Throughout the course of our collaboration, we enjoyed a successful working relationship. Importantly, his generous support encompassed various aspects of this study, including experimental performance, data analysis, and refinement of the manuscript.

Conflicts of Interest

The authors declare that they have no competing financial interests or personal relationships that could have influenced the work reported in this manuscript.

Abbreviations

EST-SSRExpress sequence tag–simple sequence repeat
ZBZanthoxylum bungeanum
ZAZanthoxylum armatum
ZPZanthoxylum piperitum
RAPDRandom Amplified Polymorphic DNA
ISSRInter-Simple Sequence Repeat
MCMCMarkov Chain Monte Carlo
PCAPrincipal Component Analysis
UPGMAUnweighted Pair Group Method with Arithmetic Mean Analysis
GSAGenetic Structure Analysis
NRNon-Redundant Protein Sequence Database
COGClusters of Orthologous Group
KEGGKyoto Encyclopedia of Genes and Genomes
GOGene Ontology

References

  1. Asase, A.; Oppong-Mensah, G. Traditional antimalarial phytotherapy remedies in herbal markets in southern Ghana. J. Ethnopharmacol. 2009, 126, 492–499. [Google Scholar] [CrossRef] [PubMed]
  2. Massoma, L.; Gasco, M.; Rubio, J.; Yucra, S.; Sock, E.N.; Gonzales, G.F. Effect of the ethanolic extract from Fagara tessmannii on testicular function, sex reproductive organs and hormonal level in adult male rats. J. Androl. 2011, 43, 139–144. [Google Scholar] [CrossRef] [PubMed]
  3. Phuyal, N.; Jha, P.K.; Raturi, P.P.; Rajbhandary, S. Zanthoxylum armatum DC.: Current knowledge, gaps and opportunities in Nepal. J. Ethnopharmacol. 2019, 229, 326–341. [Google Scholar] [CrossRef] [PubMed]
  4. Zhang, D.X.; Hartley, T.G. Flora of China; Science Press: Beijing, China, 2008; Volume 11. [Google Scholar]
  5. Nurain, I.O.; Bewaji, C.O.; Johnson, J.S.; Davenport, R.D.; Zhang, Y. Potential of Three Ethnomedicinal Plants as Antisickling Agents. Mol. Pharm. 2017, 14, 172–182. [Google Scholar] [CrossRef] [PubMed]
  6. Yu, B.; Zhang, G.; Jin, L.; Zhang, B.; Yan, D.; Yang, H.; Ye, Z.; Ma, T. Inhibition of PAI-1 activity by toddalolactone as a mechanism for promoting blood circulation and removing stasis by Chinese herb Zanthoxylum nitidum var. tomentosum. Front. Pharm. 2017, 8, 489. [Google Scholar] [CrossRef] [PubMed]
  7. Akakpo-Akue, J.; Kplé, T.K.M.; Coulibaly, K.; Ahon, G.M.; Fofié, Y.; Yapo-Crezoit, A.; Zirihi, G.N.; Kra, A.K.M. Ethnobotanical study of medicinal plants used against sickle cell anaemi-a in the eastern part of the Côte d’Ivoire. J. Anim. Plant Sci. 2020, 45, 7839–7852. [Google Scholar] [CrossRef]
  8. Tatiana, K.M.K.; Akakpo-Akue, J.; Yvette, F.; Marcel, A.G.; Julien, G.K.; Mathieu, K.A.K.; Ibrahime, S.; Chiayé, C.A.Y. Sickling cells inhibition and radical scavenging activities of Zanthoxylum leprieurii’s (GUILL) bark extracts: Comparative study. J. Pharm. Res. Int. 2020, 32, 42–50. [Google Scholar] [CrossRef]
  9. Zhang, H.; Ye, M. Research Status on the Taxonomic and Component of Green Zanthoxylum bungeanum Maxim. North Hortic. 2010, 14, 918–924. [Google Scholar]
  10. Peterson, S.W. Identification of common Aspergillus species. Mycologist 2002, 17, 128. [Google Scholar] [CrossRef]
  11. Feng, S.T.; Yang, Z.; Liu, L.; Chen, N.; Hou, Y.; Wei, A.Z. Genetic diversity and relationships of wild and cultivated Zanthoxylum germplasms based on sequence-related amplified polymorphism (SRAP) markers. Genet. Resour. Crop Evol. 2015, 62, 1193–1204. [Google Scholar] [CrossRef]
  12. Zhu, L.; Chen, X.; Li, J.H.; Jiang, W.D.; Zhong, C.P.; Zhou, J.; Wu, C.J. Textual research on origin of Zanthoxylum armatum and consideration of its definition. Chin. Tradit. Herb. Drugs 2018, 4, 987–992. [Google Scholar] [CrossRef]
  13. Fei, X.; Shi, Q.; Qi, Y.; Wang, S.; Lei, Y.; Hu, H.; Liu, Y.; Yang, T.; Wei, A.Z. ZbAGL11, a class D MADS-box transcription factor of Zanthoxylum bungeanum, is involved in sporophytic apomixis. Hortic. Res. 2021, 8, 23. [Google Scholar] [CrossRef] [PubMed]
  14. Fei, X.; Lei, Y.; Qi, Y.; Wang, S.; Hu, H.; Wei, A. Small RNA sequencing provides candidate miRNA-target pairs for revealing the mechanism of apomixis in Zanthoxylum bungeanum. BMC Plant Biol. 2021, 21, 178. [Google Scholar] [CrossRef] [PubMed]
  15. Kigen, G.; Kipkore, W.; Wanjohi, B.; Haruki, B.; Kemboi, J. Medicinal plants used by traditional healers in Sangurur, Elgeyo Marakwet County, Kenya. Pharmacogn. Res. 2017, 9, 333–347. [Google Scholar] [CrossRef] [PubMed]
  16. Ameh, S.J.; Tarfa, F.D.; Ebeshi, B.U. Traditional herbal management of sickle cell anemia: Lessons from Nigeria. Anemia 2012, 2012, 607436. [Google Scholar] [CrossRef] [PubMed]
  17. Andima, M.; Coghi, P.; Yang, L.J.; Wong, V.K.W.; Ngule, C.M.; Heydenreich, M.; Ndakala, A.J.; Yenesew, A.; Derese, S. Antiproliferative activity of secondary metabolites from Zanthoxylum zanthoxyloides Lam: In vitro and in silico studies. Pharmacogn. Commun. 2019, 10, 44–51. [Google Scholar] [CrossRef]
  18. Tine, Y.; Yang, Y.; Renucci, F.; Costa, J.; Wélé, A.; Paolini, J. LC-MS/MS analysis of flavonoid compounds from Zanthoxylum zanthoxyloides extracts and their antioxidant activities. Nat. Prod. Commun. 2017, 12, 1865–1868. [Google Scholar] [CrossRef]
  19. Zeng, Q.; Wang, Z.J.; Chen, S.; Wang, H.; Xie, T.Z.; Xu, X.J.; Xiang, M.L.; Chen, Y.C.; Luo, X.D. Phytochemical and anti-MRSA constituents of Zanthoxylum nitidum. Biomed. Pharmacother. 2022, 148, 112758. [Google Scholar] [CrossRef]
  20. Ogunbolude, Y.; Ibrahim, M.; Elekofehinti, O.O.; Adeniran, A.; Abolaji, A.O.; Rocha, J.B.T.; Kamdem, J.P. Effects of Tapinanthus globiferus and Zanthoxylum zanthoxyloides extracts on human leukocytes in vitro. J. Intercult. Ethnopharmacol. 2014, 3, 167–172. [Google Scholar] [CrossRef]
  21. Djordjevic, V.B. Free Radicals in Cell Biology. Int. Rev. Cytol. 2004, 237, 57–89. [Google Scholar]
  22. Tan, B.L.; Norhaizan, M.E.; Liew, W.P.P.; Sulaiman, R.H. Antioxidant and oxidative stress: A mutual interplay in age-related diseases. Front. Pharmacol. 2018, 9, 1162. [Google Scholar] [CrossRef] [PubMed]
  23. Willcox, J.K.; Ash, S.L.; Catignani, G.L. Antioxidants and prevention of chronic disease. Crit. Rev. Food Sci. Nutr. 2004, 44, 275–295. [Google Scholar] [CrossRef] [PubMed]
  24. Medhi, K.; Sarmah, D.K.; Deka, M.; Bhau, B.S. High gene flow and genetic diversity in three economically important Zanthoxylum Spp. of Upper Brahmaputra Valley Zone of NE India using molecular markers. Meta Gene 2014, 2, 706–721. [Google Scholar] [CrossRef] [PubMed]
  25. Sun, Y.L.; Park, W.G.; Kwon, O.W.; Hong, S.W. Ribosomal DNA internal transcribed spacer 1, and internal transcribed spacer 2 regions as targets for molecular identification of medically important Zanthoxylum schinifolium. Afr. J. Biotechnol. 2010, 9, 4661–4673. [Google Scholar]
  26. Feng, S.; Liu, Z.; Hu, Y.; Tian, J.; Yang, T.; Wei, A. Genomic analysis reveals the genetic diversity, population structure, evolutionary history and relationships of Chinese pepper. Hortic. Res. 2020, 7, 158. [Google Scholar] [CrossRef] [PubMed]
  27. Kim, Y.M.; Jo, A.; Jeong, J.H.; Kwon, Y.R.; Kim, H.B. Development and characterization of microsatellite primers for Zanthoxylum schinifolium (Rutaceae). Appl. Plant Sci. 2017, 5, 1600145. [Google Scholar] [CrossRef] [PubMed]
  28. Su, K.; Sun, J.; Han, J.; Zheng, T.; Sun, B.; Liu, S. Combined morphological and multi-omics analyses to reveal the developmental mechanism of Zanthoxylum bungeanum prickles. Front. Plant Sci. 2022, 13, 950084. [Google Scholar] [CrossRef]
  29. Tangphatsornruang, S.; Sraphet, S.; Singh, R.; Okogbenin, E.; Triwitayakorn, K. Development of polymorphic markers from expressed sequence tags of Manihot esculenta Crantz. Mol. Ecol. Resour. 2008, 8, 682–685. [Google Scholar] [CrossRef]
  30. Holton, T.A.; Christopher, J.T.; McClure, L.; Harker, N.; Henry, R.J. Identification and mapping of polymorphic SSR markers from expressed gene sequences of barley and wheat. Mol. Breed. 2002, 9, 63–71. [Google Scholar] [CrossRef]
  31. Scott, K.D.; Eggler, P.; Seaton, G.; Rossetto, M.; Ablett, E.M.; Lee, L.S.; Henry, R.J. Analysis of SSRs derived from grape ESTs. Theor. Appl. Genet. 2000, 100, 723–726. [Google Scholar] [CrossRef]
  32. Kantety, R.V.; La Rota, M.; Matthews, D.E.; Sorrells, M.E. Data mining for simple sequence repeats in expressed sequence tags from barley, maize, rice, sorghum and wheat. Plant Mol. Biol. 2002, 48, 501–510. [Google Scholar] [CrossRef] [PubMed]
  33. Song, Q.J.; Marek, L.F.; Shoemaker, R.C.; Lark, K.G.; Concibido, V.C.; Delannay, X.; Specht, J.E.; Cregan, P.B. A new integrated genetic linkage map of the soybean. Theor. Appl. Genet. 2004, 109, 122–128. [Google Scholar] [CrossRef] [PubMed]
  34. Chen, X.C.; Ping, Z.; Young, A.C.; Shu, H.; Fred, G. Mining and characterizing microsatellites from citrus ESTs. Theor. Appl. Genet. 2006, 112, 1248–1257. [Google Scholar] [CrossRef] [PubMed]
  35. Hou, L.X.; Wei, A.Z.; Wang, L.H.; Liu, Y.L. Analysis of SSR Loci and Development of Molecular Markers in Zanthoxylum bungeanum Transcriptome. Chin. J. Agric. Biotech. 2018, 26, 1226–1236. [Google Scholar]
  36. Kim, J.A.; Roy, N.S.; Lee, I.H.; Choi, A.Y.; Choi, B.S.; Yu, Y.S.; Park, N.I.; Park, K.C.; Kim, S.; Yang, H.S.; et al. Genome-wide transcriptome profiling of the medicinal plant, Zanthoxylum planispinum, using a single-molecule direct RNA sequencing approach. Genomics 2019, 111, 973–979. [Google Scholar] [CrossRef] [PubMed]
  37. Feng, S.; Zhao, L.; Liu, Z.; Liu, Y.; Yang, T.; Wei, A. De novo transcriptome assembly of Zanthoxylum bungeanum using Illumina sequencing for evolutionary analysis and simple sequence repeat marker development. Sci. Rep. 2017, 7, 16754. [Google Scholar] [CrossRef] [PubMed]
  38. Deng, Y.C.; Xiang, L.; Su, Y.Y.; Yang, L.; Li, H.Y.; Luo, X.M.; Sun, W.; Ye, M. Development of EST-SSR markers of Zanthoxylum species based on RNA-sequencing. J. Northwest A F Univ. 2019, 4, 16–24. [Google Scholar] [CrossRef]
  39. Guo, W.Z.; Zhao, J.X.; Cao, J.H. Cultivation and management techniques of Zanthoxylum piperitum. North. Hortic 2012, 20, 36–37. [Google Scholar]
  40. Nei, M. Genetic Distance between Populations. Am. Nat. 1972, 106, 283–292. [Google Scholar] [CrossRef]
  41. Lewontin, R.C. Testing the Theory of Natural Selection. Nature 1972, 236, 181–182. [Google Scholar] [CrossRef]
  42. Peakall, R.; Peter, E.S. GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research—An update. Bioinformatics 2012, 28, 2537–2539. [Google Scholar] [CrossRef] [PubMed]
  43. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of Population Structure Using Multilocus Genotype Data. Genetics 2000, 155, 9197–9201. [Google Scholar] [CrossRef] [PubMed]
  44. Falush, D.; Stephens, M.; Pritchard, J.K. Inference of population structure using multilocus genotype data: Linked loci and correlated allele frequencies. Genetics 2003, 164, 1567–1587. [Google Scholar] [CrossRef] [PubMed]
  45. Evanno, G.S.; Regnaut, S.J.; Goudet, J. Detecting the number of clusters of individuals using the software Structure: A simulation study. Mol. Ecol. 2005, 14, 2611–2620. [Google Scholar] [CrossRef] [PubMed]
  46. Anuradha, U.; Kadam, U.S.; Priya, C.M.; Karibasappa, G.S. Microsatellite and RAPD analysis of grape (Vitis spp.) accessions and identification of duplicates/misnomers in germplasm collection. Indian J. Hortic. 2010, 67, 8–15. [Google Scholar] [CrossRef]
  47. Upadhyay, A.; Kadam, U.S.; Chacko, P.; Aher, L.; Karibasappa, G.M. Microsatellite analysis to differentiate clones of Thompson Seedless grapevine. Indian J. Hortic. 2010, 67, 260–263. [Google Scholar] [CrossRef]
  48. Mir, A.H.; Bhat, M.A.; Fayaz, H.; Wani, A.A.; Dar, S.A.; Maqbool, S.; Yasin, M.; Mir, J.I.; Khan, M.A.; Sofi, P.A.; et al. SSR markers in revealing extent of genetic diversity and phylogenetic relationships among chickpea core collection accessions for Western Himalayas. Mol. Biol. Rep. 2022, 49, 11469–11479. [Google Scholar] [CrossRef]
  49. Chen, J.; Liu, Q.; Lu, C.; Liu, Q.; Pan, J.; Zhang, J.; Dong, S. Genetic diversity of Prunus armeniaca L. var. ansu Maxim. germplasm revealed by simple sequence repeat (SSR) markers. PLoS ONE 2022, 17, e0269424. [Google Scholar] [CrossRef]
  50. Kazemeini, F.; Asri, Y.; Mostafavi, G.; Kalvandi, R.; Mehregan, I. RAPD-based evaluation of genetic diversity among populations of the Iranian endemic species Rhabdosciadium aucheri Boiss. (Apiaceae). Mol. Biol. Rep. 2020, 47, 9345–9352. [Google Scholar] [CrossRef]
  51. Sen, S.; Skaria, R.; Muneer, P.M.A. Genetic diversity analysis in Piper species (Piperaceae) using RAPD markers. Mol. Biotechnol. 2010, 46, 72–79. [Google Scholar] [CrossRef]
  52. Shahnazari, N.; Noormohammadi, Z.; Sheidai, M.; Koohdar, F. A new insight on genetic diversity of sweet oranges: CAPs-SSR and SSR markers. J. Genet. Eng. Biotechnol. 2022, 20, 105. [Google Scholar] [CrossRef] [PubMed]
  53. Chavhan, R.L.; Sable, S.; Narwade, A.V.; Hinge, V.R.; Kalbande, B.B.; Mukherjee, A.K.; Chakrabarty, P.K.; Kadam, U.S. Multiplex molecular marker-assisted analysis of significant pathogens of cotton (Gossypium sp.). Biocatal. Agric. Biotechnol. 2023, 47, 102557. [Google Scholar] [CrossRef]
  54. Hu, C.Y.; Tsai, Y.Z.; Lin, S.F. Development of STS and CAPS markers for variety identification and genetic diversity analysis of tea germplasm in Taiwan. Bot. Stud. 2014, 55, 12. [Google Scholar] [CrossRef] [PubMed]
  55. Choi, H.J. Evaluation of Antiviral Activity of Zanthoxylum Species against Picornaviruses. Osong Public Health Res. Perspect. 2016, 7, 400–403. [Google Scholar] [CrossRef] [PubMed]
  56. Alam, F.; Us Saqib, Q.N. Evaluation of Zanthoxylum armatum Roxb for in vitro biological activities. J. Tradit. Chin. Med. 2017, 7, 515–518. [Google Scholar] [CrossRef] [PubMed]
  57. Varshney, R.K.; Thiel, T.; Stein, N.; Langridge, P.; Graner, A. In silico analysis on frequency and distribution of microsatellites in ESTs of some cereal species. Cell. Mol. Biol. Lett. 2002, 7, 537–546. [Google Scholar] [CrossRef] [PubMed]
  58. Zhao, Y.C.; Liu, X.W.; Li, Y.; Zhao, D.G. Transcriptome Data Assembly and Gene Function Annotation of Zanthoxylum piperitum DC var. inerme Makino. Genom. Appl. Biol. 2016, 35, 1805–1819. [Google Scholar] [CrossRef]
  59. Eujayl, I.; Sorrells, M.; Baum, M.; Wolters, P.; Powell, W. Assessment of genotypic variation among cultivated durum wheat based on EST-SSRS and genomic SSRs. Euphytica 2001, 119, 39–43. [Google Scholar] [CrossRef]
  60. Thiel, T.; Michalek, W.; Varshney, R.; Graner, A. Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.). Theor. Appl. Genet. 2003, 106, 411–422. [Google Scholar] [CrossRef]
  61. Dong, J.; Guang, Y.Z. Analysis of Microsatellites in Citrus Unigenes. Act. Genet. Sin. 2006, 33, 345–353. [Google Scholar] [CrossRef]
  62. Zhao, L.L.; Feng, S.J.; Tian, J.Y.; Wei, A.Z.; Yang, T.X. Internal transcribed spacer 2 (ITS2) barcodes: A useful tool for identifying Chinese Zanthoxylum. Appl. Plant. Sci. 2018, 6, e01157. [Google Scholar] [CrossRef] [PubMed]
  63. Sun, Y.L.; Park, W.G.; Kwon, O.W.; Soon, K.H. Molecular identification of Zanthoxylum piperitum by the internal transcribed spacers as targets using newly designed species-specific primers. J. Biotechnol. 2010, 150, 505. [Google Scholar] [CrossRef]
  64. Wang, C.; Zhao, J.X.; Bi, J.; Wang, C.R.; Guo, W.Z.; Lin, Y. The Techniques of Fruitsetting Improvement of Zanthoxylu piperitum f. inerme Makino. North. Hortic. 2012, 26, 96–99. [Google Scholar] [CrossRef]
  65. Suprasanna, P.; Jain, S.M.; Ochatt, S.J.; Kulkarni, V.M.; Predieri, S. Applications of In Vitro Techniques in Mutation Breeding of Vegetatively Propagated Crops; CABI International: Wallingford, UK, 2012. [Google Scholar] [CrossRef]
  66. Beyaz, R.; Yildiz, M. The Use of Gamma Irradiation in Plant Mutation Breeding. In Plant Engineering; InTech: Tokyo, Japan, 2017. [Google Scholar] [CrossRef]
  67. Ahloowalia, B. Induced mutations-A new paradigm in plant breeding. Euphytica 2001, 118, 167–173. [Google Scholar] [CrossRef]
  68. Wright, S. Evolution in Mendelian populations. Bull. Math. Biol. 1990, 52, 241–295. [Google Scholar] [CrossRef] [PubMed]
  69. Ellstrand, N.C.; Roose, M.L. Patterns of genotypic diversity in clonal plant species. Am. J. Bot. 1987, 74, 123–131. [Google Scholar] [CrossRef]
  70. Buso, G.S.C.; Rangel, P.H.; Ferreira, M.E. Analysis of genetic variability of south American wild rice populations (Oryza glumaepatula) with isozymes and RAPD markers. Mol. Ecol. 1998, 7, 107–118. [Google Scholar] [CrossRef]
Cimb 45 00454 g001
Figure 1. Map showing the sites from which 121 Zanthoxylum samples were collected for genetic diversity analysis. (A): Distribution of sampling sites throughout China. Scale bar = 0–2160 km. (B): Distribution of sampling points in locally zoomed maps. Scale bar = 0–400 km. ZB (67), ZA (47), ZP (7). ZB: Zanthoxylum bungeanum. ZA: Zanthoxylum armatum. ZP: Zanthoxylum piperitum. Green pot: sample collection location. a: Three samples, Z. bungeanum (cultivars, accession: 9–11). b: Sixty-three samples, Z. bungeanum: (56 cultivars, accession: 12–67); Z. piperitum (7 cultivars, accession: 115–121). c: One sample, Z. bungeanum (cultivar, accession: 8). d: Fourteen samples, Z. armatum (14 cultivars and 3 wild-type, accession: 68–81). e: Two samples, Z. bungeanum (cultivars, accession: 4–5). f: Three samples, Z. armatum (wild-type, accession: 82–84). g: Ten samples, Z. armatum (cultivars, accession: 86–88, 98–104). h: Nine samples, Z. armatum (cultivars, accession: 89–97). i: Three samples, Z. bungeanum (2 cultivars, accession: 6–7); Z. armatum (1 cultivar, accession: 85). j: One sample, Z. bungeanum (cultivar, accession: 1). k: 1 sample, Z. bungeanum (cultivar, accession: 2). l: Two samples, Z. bungeanum (1 cultivar, accession: 3); Z. armatum (1 cultivar, accession: 113). m: One sample, Z. armatum (cultivar, accession: 112). n: Seven samples, Z. armatum (cultivars, accession: 105–111). o: One sample, Z. armatum (cultivar, accession: 114). All sample information is shown in Supplementary Table S1.
Figure 1. Map showing the sites from which 121 Zanthoxylum samples were collected for genetic diversity analysis. (A): Distribution of sampling sites throughout China. Scale bar = 0–2160 km. (B): Distribution of sampling points in locally zoomed maps. Scale bar = 0–400 km. ZB (67), ZA (47), ZP (7). ZB: Zanthoxylum bungeanum. ZA: Zanthoxylum armatum. ZP: Zanthoxylum piperitum. Green pot: sample collection location. a: Three samples, Z. bungeanum (cultivars, accession: 9–11). b: Sixty-three samples, Z. bungeanum: (56 cultivars, accession: 12–67); Z. piperitum (7 cultivars, accession: 115–121). c: One sample, Z. bungeanum (cultivar, accession: 8). d: Fourteen samples, Z. armatum (14 cultivars and 3 wild-type, accession: 68–81). e: Two samples, Z. bungeanum (cultivars, accession: 4–5). f: Three samples, Z. armatum (wild-type, accession: 82–84). g: Ten samples, Z. armatum (cultivars, accession: 86–88, 98–104). h: Nine samples, Z. armatum (cultivars, accession: 89–97). i: Three samples, Z. bungeanum (2 cultivars, accession: 6–7); Z. armatum (1 cultivar, accession: 85). j: One sample, Z. bungeanum (cultivar, accession: 1). k: 1 sample, Z. bungeanum (cultivar, accession: 2). l: Two samples, Z. bungeanum (1 cultivar, accession: 3); Z. armatum (1 cultivar, accession: 113). m: One sample, Z. armatum (cultivar, accession: 112). n: Seven samples, Z. armatum (cultivars, accession: 105–111). o: One sample, Z. armatum (cultivar, accession: 114). All sample information is shown in Supplementary Table S1.
Cimb 45 00454 g001
Cimb 45 00454 g002
Figure 2. Zanthoxylum clustering dendrogram of 121 samples based on the genetic distance and the UPGMA method using NTSYSpc 2.11 software and MEGA7.0 (tree scale: 0.1). ZB: Zanthoxylum bungeanum. ZA: Zanthoxylum armatum. ZP: Zanthoxylum piperitum. Red fonts: ZB samples in the clade. Yellow fonts: ZP samples in the clade.
Figure 2. Zanthoxylum clustering dendrogram of 121 samples based on the genetic distance and the UPGMA method using NTSYSpc 2.11 software and MEGA7.0 (tree scale: 0.1). ZB: Zanthoxylum bungeanum. ZA: Zanthoxylum armatum. ZP: Zanthoxylum piperitum. Red fonts: ZB samples in the clade. Yellow fonts: ZP samples in the clade.
Cimb 45 00454 g002
Cimb 45 00454 g003
Figure 3. Principal component analysis of 121 Zanthoxylum samples. ZB: Zanthoxylum bungeanum. ZA: Zanthoxylum armatum. ZP: Zanthoxylum piperitum. PC1: Principal component 1. PC2: Principal component 2.
Figure 3. Principal component analysis of 121 Zanthoxylum samples. ZB: Zanthoxylum bungeanum. ZA: Zanthoxylum armatum. ZP: Zanthoxylum piperitum. PC1: Principal component 1. PC2: Principal component 2.
Cimb 45 00454 g003
Cimb 45 00454 g004
Figure 4. Genetic structure and K-value analysis of all samples. (a) Reasonable number of subgroups for all samples according to K-value analysis. (b) Genetic structure analysis of all Zanthoxylum samples in this study. ZB: Zanthoxylum bungeanum. ZA: Zanthoxylum armatum. ZP: Zanthoxylum piperitum.
Figure 4. Genetic structure and K-value analysis of all samples. (a) Reasonable number of subgroups for all samples according to K-value analysis. (b) Genetic structure analysis of all Zanthoxylum samples in this study. ZB: Zanthoxylum bungeanum. ZA: Zanthoxylum armatum. ZP: Zanthoxylum piperitum.
Cimb 45 00454 g004
Table 1. The molecular variance (AMOVA) for both interspecies and intraspecies variations.
Table 1. The molecular variance (AMOVA) for both interspecies and intraspecies variations.
Source of VariationDegrees of Freedom (df)Statistics (SS)Mean Square (MS)Est. Var.Percentage of Variance Componentp-Value
Among Populations
(AP)		Interspecific populations of
ZB, ZA and ZP		2149.54374.7712.22443%p > 0.001
Intraspecific populations of ZB21130.7006.2241.72153%p > 0.001
Intraspecific populations of ZA929.3203.2580.52329%p > 0.001
Within Populations
(WP)		Interspecific populations of
ZB, ZA and ZP		118349.8542.9652.96557%\
Intraspecific populations of ZB3858.3331.5351.53547%\
Intraspecific populations of ZA2936.7061.2661.26671%\
TotalInterspecific populations of
ZB, ZA and ZP		120499.397\5.189100%\
Intraspecific populations of ZB59189.033\3.256100%\
Intraspecific populations of ZA3866.026\1.789100%\
Genetic differentiation coefficient
(PhiPT)		Interspecific populations of ZB, ZA and ZP0.429
Intraspecific populations of ZB0.529
Intraspecific populations of ZA0.293
Gene flow
(Nm)		Interspecific populations of ZB, ZA and ZP0.333
Intraspecific populations of ZB0.223
Intraspecific populations of ZA0.605
PhiPT = AP/(WP + AP) = AP/Total; Nm = [(1/PhiPT) − 1]/4; AP = Est. Var. Among Pops, WP = Est. Var. Within Pops.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Deng, Y.; He, Z.; Li, Y.; Ye, M.; Xiang, L. Six Express Sequence Tag–Simple Sequence Repeat Primers Reveal Genetic Diversity in the Cultivars of Three Zanthoxylum Species. Curr. Issues Mol. Biol. 2023, 45, 7183-7196. https://doi.org/10.3390/cimb45090454

AMA Style

Deng Y, He Z, Li Y, Ye M, Xiang L. Six Express Sequence Tag–Simple Sequence Repeat Primers Reveal Genetic Diversity in the Cultivars of Three Zanthoxylum Species. Current Issues in Molecular Biology. 2023; 45(9):7183-7196. https://doi.org/10.3390/cimb45090454

Chicago/Turabian Style

Deng, Yangchuan, Zhoujian He, Yanlin Li, Meng Ye, and Li Xiang. 2023. "Six Express Sequence Tag–Simple Sequence Repeat Primers Reveal Genetic Diversity in the Cultivars of Three Zanthoxylum Species" Current Issues in Molecular Biology 45, no. 9: 7183-7196. https://doi.org/10.3390/cimb45090454

APA Style

Deng, Y., He, Z., Li, Y., Ye, M., & Xiang, L. (2023). Six Express Sequence Tag–Simple Sequence Repeat Primers Reveal Genetic Diversity in the Cultivars of Three Zanthoxylum Species. Current Issues in Molecular Biology, 45(9), 7183-7196. https://doi.org/10.3390/cimb45090454

Article Metrics

Back to TopTop