1. Background
Prunus salicina L., belonging to the family of Rosaceae, are one of the most important economical fruit trees and are widely cultivated all over the world. They are mainly distributed in the northern hemisphere, especially in the temperate zone [Reference Chen, Chen and Chen1, Reference Herrera, Lora, Hormaza, Herrero and Rodrigo2]. China is one of the origin and distribution centers of Prunus L. species. Prunus L. species contain more than 430 species and are first segregated into six genera according to the morphology of fruit: Amygdalus L., Armeniaca Scop., Cerasus Mill., Laurocerasus, Padus Mill., Prunus species, and Tourn. ex Duh. However, phylogenetic analysis showed that Cerasus, Laurocerasus, and Padus were not monophyletic [Reference Wen, Berggren and Lee3, Reference Williams, Kubelik, Livak, Rafalski and Tingey4]. Besides, an increasing number of new cultivars from different countries result an important renewal of plant material worldwide [Reference Herrera, Lora, Hormaza, Herrero and Rodrigo2]. It is thus necessary to characterize genetic information of Prunus L. species to cultivate new breed with improved quality characteristics.
DNA polymorphism assay based on the amplification of random DNA segments with single primers of arbitrary nucleotide sequence has been widely used for genetic diversity analysis of species [Reference Pairon, Petitpierre and Campbell5]. Several studies have been devoted to the genetic diversity in Prunus L. species [Reference Fu, Li and Xu6–Reference Carrasco, González, Gebauer and García-González8]. Recently, a number of molecular marker techniques including random amplified polymorphic DNA (RAPD), simple sequence repeat (SSR), intersimple sequence repeat (ISSR), and amplified fragment length polymorphism (AFLP) have been developed and widely used in the identification of various organisms [Reference Williams, Kubelik, Livak, Rafalski and Tingey4, Reference Fu, Li and Xu6, Reference Godwin, Aitken and Smith9–Reference Liscum, Innis, Gelfand and Sninsky11]. Among these techniques, RAPD and ISSR methods are two PCR-based methods that require only small amounts of DNA sample without involving radioactive labels and therefore have been widely used for genetic characterization [Reference Verma, Ul Haq, Kachhwaha and Kothari12]. RAPD is a technique based on the amplification of the genomic DNA with either a single or multiple short oligonucleotide primers of an arbitrary or random sequence [Reference Verma, Ul Haq, Kachhwaha and Kothari12]. RAPD is simple, cost-efficient, and does not require DNA sequences before application [Reference Tripathi, Chouhan, Saini and Tiwari13]. ISSR is derived from SSR, which is more abundant, informative, highly polymorphic, and efficient [Reference Chaudhary, Kumar and Sharma14]. RAPD and ISSR methods have separately been used for genetic characterization in many species, such as Lonicera japonica Thunb. [Reference Fu, Yang, Khan and Mei15], synthetic hexaploid wheats [Reference Shakeel, Ilyas and Kazi16], Atractylodes lancea [Reference Kohjyouma, Nakajima, Namera, Shimizu, Mizukami and Kohda17], and Ocimum basilicum L. [Reference Giachino, Sönmez and Tonk18]. However, because of their advantages and disadvantages, more studies applied RAPD combined with ISSR to characterize the genetic variation of species, such as Litchi chinensis Sonn. [Reference Long, Cheng and Mei19], Allium species [Reference Mukherjee20], date palm [Reference El Sharabasy and Soliman21], and Cymbopogon [Reference Bishoyi, Sharma, Kavane and Geetha22]. However, only limited studies have been conducted to characterize the genetic relationships among different genus or cultivars of Prunus L. species [Reference Erturk, Ercisli, Maghradze, Orhan and Agar23–Reference Baránek, Raddová and Pidra25].
In this study, we applied the RAPD and ISSR technique for the genetic characterization of 25 P. salicina from China and other countries. This study may provide valuable insight into the genetic diversity of P. salicina L. and provide information to cultivate new breed with improved traits.
2. Methods
2.1. Plant Material Collection and DNA Extraction
This study included 25 P. salicina L. which were collected from 14 different regions from China (13 samples), Japan (4 samples), USA (7 samples), and New Zealand (1 sample) (Figure 1 and Table 1). Among them, P1, P2, and P3 are the three lines with different maturity of one cultivar. The flowers of the 25 P. salicina L. are shown in Figure 2.
The genomic DNA of 25 P. salicina L. was extracted from fresh leaves using a modified cetyl trimethylammonium bromide (CTAB) method as described previously [Reference Fu, Yang, Khan and Mei15, Reference Yang, Fu, Khan, Zeng and Fu26]. DNA integrity was checked by 0.8% agarose gel electrophoresis, and DNA purity was determined by the absorbance ratio at 260 nm : 280 nm on spectrophotometry. The final concentration of DNA samples was adjusted to 10 ng/µl for PCR and stored at −20°C until use.
2.2. Amplification of DNA by RAPD-PCR
The random RAPD primers were selected randomly for PCR amplification (Table 2). The PCR system in 10 μL volume contains 1 μL of 2.5 μmol/L primers, 1 μL of DNA template, 5 μL of 2 × PCR Taq Mastermix (TianGen Biotech Co. Ltd., Beijing), and 3 μL of deionized water. The PCR was executed on Applied Biosystems Veriti 96-Well Thermal Cycler (Thermo Fisher, USA) in the following procedure: initial denaturation at 95°C for 90 s, followed by 40 cycles of 40 s at 94°C, 60 s at 36°C, 90 s at 72°C, and final extension of 5 min at 72°C.
Note. R = (A/G), Y = (C/T), and D = (A/G/T); aaverage of the column.
2.3. ISSR Amplification
Fifteen ISSR primers were synthesized by Thermo Fisher (USA) (Table 2). ISSR amplification was performed in 10 μL reactions including 1 μL of 2.5 umol/L primers, 1 μL of DNA template, 5 μL of 29 PCR Taq Mastermix (TianGen Biotech Co. Ltd., Beijing), and 3 μL of deionized water. PCR was executed on Applied Biosystems Veriti 96-Well Thermal Cycler using the following procedure: initial denaturation at 95°C for 90 s, followed by 35 cycles of 30 s at 94°C, 30 s at 50°C, 90 s at 72°C, and final extension of 5 min at 72°C [Reference Fu, Yang, Khan and Mei15].
2.4. Agarose Gel Electrophoresis
The amplified PCR products were separated by electrophoresis on 1.8% agarose gel in 1 × TAE buffer. Gels were visualized by 0.5 g/ml ethidium bromide staining, and the images were documented using the ChemiDoc XR (Bio-Rad, USA). Bands that were unambiguous and reproducible in successive amplifications were selected for scoring.
2.5. Data Analysis
All PCRs were repeated five times for each of five samples. Bands in the gel profiles were scored as 1 for present and 0 for absent. The similarity matrix (SM) and the similarity index (SI) were calculated using SM coefficient in Numerical Taxonomy Multivariate Analysis System (NTSYS pc 2.1) software. The dendrogram based on the unweighted pair group method with arithmetic mean algorithm (UPGMA) was generated using the SAHN module in the NTSYS pc 2.1 software.
3. Results
3.1. Amplification of DNA by RAPD and ISSR
A total of nineteen RAPD primers and fifteen ISSR primers were used in this study for the evaluation of DNA polymorphism (Table 2). All RAPD primers and ISSR primers generated evaluable bands. Figure 3 shows the representative reproducible polymorphic amplification bands in these 25 samples generated from ISSR primer UBC807 and RAPD primer S201. For the RAPD primers, a total of 315 bands with an average of 16.58 bands per primer were obtained. Among these bands, 257 (81.59%) bands were polymorphic, and the approximate band size ranged from 200 bp to 2300 bp. The minimum number of bands was 10, which was generated by primer OPA-4 and the maximum was 21, which was produced by primer S43. The total number of polymorphic fragments ranged from 7 (primer OPA-4) to 18 (primer OPA-10). The average polymorphic fragments ratio (PFR) (in %) was 81.60% (min: 65%; max: 94.74%). The other information of the bands generated by RAPD primers, including polymorphism information content (PIC), resolving power (RP), effective multiplex ratio (EMR), and marker index (MI), are presented in Table 3.
Note. aAverage of the column. TF, total number of fragments; PF, number of polymorphic fragments; PFR, polymorphic fragments ratios (%); PIC, polymorphism information content; RP, resolving power; EMR, effective multiplex ratio; MI, marker index.
For the ISSR primers, a total of 204 bands with an average of 13.60 bands per primer were produced; of them, 179 (87.74%) were polymorphic. The approximate range of band size was 300 bp to 2500 bp (Table 4). The minimum number of bands was 8, which was yielded by primer UBC829, and the maximum was 19, which was produced by primer UBC807. The total number of PFs ranged from 8 (primer UBC829) to 15 (primers UBC807, UBC810, UBC846, and UBC881). The average PFR% was 87.80% (min: 69.23%; max: 100%). The other information of the bands generated by ISSR primers, including PIC, RP, EMR, and MI, are presented in Table 4.
Note. aAverage of the column. TF, total number of fragments; PF, number of polymorphic fragments; PFR, polymorphic fragments ratios (%); PIC, polymorphism information content; RP, resolving power; EMR, effective multiplex ratio; MI, marker index.
3.2. Genetic Distance and Cluster Analysis of RAPD and ISSR Markers
Based on the RAPD amplification profiles, cluster dendrogram was obtained using UPGMA (Figure 4). Since P1, P2, and P3 belong to one cultivar, we ignored their coefficients in the following analysis. The dendrogram showed that the similarity coefficients ranged from 0.584 to 0.860. In the RAPD-based dendrogram, the 25 P. salicina samples formed four clusters at a cutoff of 0.692. The similarity coefficient between sample P4 (Qiepili in Ningbo, Zhejiang, China) and P13 (Longyuanqiuli in Heilongjiang, China) was lowest (0.584), while that between sample P19 (Laroda in USA) and P20 (Red heart in USA) was highest (0.860) (Figure 4).
The ISSR analysis showed similar results to the RAPD analysis. The dendrogram showed that the similarity coefficients ranged from 0.558 to 0.892. In the ISSR-based dendrogram, the 25 P. salicina samples were divided into five clusters at a cutoff of 0.692. The similarity coefficient between sample P10 (Maihuangli in Anhui, China) and P13 (Longyuanqiuli in Heilongjiang, China) was lowest (0.558), while that between sample P10 (Maihuangli in Anhui, Chin) and P15 (Baili in Japan) was highest (0.892) (Figure 5).
3.3. Integrating Analysis of RAPD and ISSR Data
The dendrogram results of RAPD combined with ISSR showed that the similarity coefficients ranged from 0.597 to 0.865. Total 519 DNA fragments were yielded, of which 435 (84.7%) were polymorphic. The average number of PF per primer was 12.7. The mean PIC, RP, EMR, and MI values observed for all primers were 0.42, 17.77, 10.80, and 4.54, respectively (Table 5). The similarity coefficients between sample P10 (Maihuangli in Anhui, China) and P13 (Longyuanqiuli in Heilongjiang, China) was lowest (0.597), while that between sample P10 (Maihuangli in Anhui, China) and P15 (Baili in Japan) was highest (0.865) (Table 6).
Note. The bold values indicate the maximum and minimum genetic similarity values among the landraces.
3.4. Typical Band Patterns Amplified by ISSR and RAPD Markers
Sixteen primers, including 11 ISSR primers and 5 RAPD primers, could be used as the markers of molecular identification for 25 Prunus L. samples (Table 7). As shown in Table 7, UBC810, UBC834, and UBC836 could be considered as the markers of P1 (Zuili1 in Jiaxing, Zhejiang, China), P2 (Zuili2 in Jiaxing, Zhejiang, China), and P3 (Zuili3 in Jiaxing, Zhejiang, China). S17 could be considered as a marker of P4 (Qiepili in Ningbo, Zhejiang, China). UBC881 might be a marker of P5 (Jintangli in Zhoushan, Zhejiang, China). UBC847 was a marker of P6 (Furongli in Fujian, China). UBC847 and UBC855 could be used to distinguish P7 (Yuhuangli in Hubei, China). UBC848 could be considered as a marker of P8 (Jiuqianli in Guizhou, China). UBC857 might be a potential marker of P9 (Huangguli in Tongxiang, Zhejiang, China). RAPD-1 could be used as a marker of P12 (Niuxinli in Shandong, China). UBC889 could be considered as a marker of P16 (Akihime in Japan). S43 and S1403 might be the markers of P17 (Zhenzhuli in Japan). UBC829 might be a potential marker of P23 (Queen rose in USA). RAPD-5 also might be used a marker of P25 (Misili in New Zealand). The representative banding profiles obtained by ISSR primers UBC834, UBC847, UBC857, and RAPD primer S1403 are shown in Figure 6.
4. Discussion
Illustration of the genetic relationships or characterization of genetic diversity is important to provide genetic guidance for hybrid breeding. In this study, the genetic diversity and relationship among 25 P. salicina L. varieties were evaluated by RAPD and ISSR, respectively, and integrated. The Dice’s similarity coefficient of RAPD ranged from 0.584 to 0.860, and that of ISSR ranged from 0.558 to 0.892. Integrating analysis of RAPD and ISSR indicated the similarity coefficient varied from 0.597 to 0.865. The results indicated high diversity among the 25 varieties.
ISSR and RAPD were widely used for genetic diversity evaluations of Prunus L. species. Tian et al. used ISSR and RAPD for genetic diversity evaluations of 48 Prunus mira L. samples, the high levels of polymorphism, and the results imply that Tibet samples preserved higher genetic diversity and most genetic variations occurred [Reference Tian, Xing and Cao27]. However, the efficiency of RAPD markers and ISSR markers in detecting polymorphism is controversial. Tian et al. demonstrated that ISSR found 77.80% polymorphism, which is higher than that found by RAPD (72.73%). In the study of Kumar et al. the phylogenetic relationships of 36 locally grown P. armeniaca genotypes were analyzed using 20 RAPDs and 11 ISSRs markers. RAPD markers were found more efficient for polymorphism detection, as they detected 97.84% as compared to 96.5% for ISSR markers, and the pattern of clustering of the genotypes remained more or less the same in RAPD and combined data of RAPD + ISSR [Reference Kumar, Mishra, Singh, Kumar, Naik and Singh28]. In our study, the PFR% of RAPD primers was 81.60%, which is lower than that of ISSR primers (87.80%). Our results support the view that ISSR markers are more efficient than RAPD with regards to detecting polymorphism.
The RAPD results showed that the index of similarity coefficient between sample P4 (Qiepili in Ningbo, Zhejiang, China) and P13 (Longyuanqiuli in Heilongjiang, China) was lowest (0.584), while that between sample P19 (Laroda in USA) and P20 (Red heart in USA) was highest (0.860). However, the ISSR results showed that the index of similarity coefficient between sample P10 (Maihuangli in Anhui, Chin) and P13 (Longyuanqiuli in Heilongjiang, China) was lowest (0.558), while that between sample P10 (Maihuangli in Anhui, China) and P15 (Baili in Japan) was highest (0.892). In addition, the analysis of RAPD combined with ISSR showed that the similarity coefficient between sample P10 (Maihuangli in Anhui, China) and P13 (Longyuanqiuli in Heilongjiang, China) was lowest (0.597), while that between sample P10 (Maihuangli in Anhui, Chin) and P15 (Baili in Japan) was highest (0.865), which was consistent with the RAPD analysis. These findings demonstrated that the RAPD technique not only increased the resolution and yield but also was a reliable molecular tool for the genetic characterization of various organisms, which was reported in previous studies [Reference Fu, Li and Xu6, Reference Fu, Yang, Khan and Mei15]. Our RAPD and ISSR analysis showed potentiality to distinguish P. salicina L. from related genus or species.
5. Conclusion
In summary, our study indicates that the RAPD combined with ISSR techniques would be used for the genetic diversity, molecular-assisted breeding, and genetic characterization of P. salicina L. Our results might assist in parental gametophytes selection for hybrid breeding of P. salicina L.
Abbreviations
RAPD: Random amplified polymorphic DNA
SSR: Simple sequence repeat
ISSR: Intersimple sequence repeat
AFLP: Amplified fragment length polymorphism
CTAB: Cetyl trimethylammonium bromide
SM: Similarity matrix
SI: Similarity index
NTSYS: Numerical taxonomy multivariate analysis system
UPGMA: Unweighted pair group method with arithmetic mean algorithm
PIC: Polymorphism information content
RP: Resolving power
EMR: Effective multiplex ratio
MI: Marker index
TF: Total number of fragments
PF: Polymorphic fragments
PFR (%): Polymorphic fragments ratios (%).
Data Availability
The data that support the findings of this study are available on request to the corresponding author.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Authors’ Contributions
JL and GCG designed experiments. Bin acquired data. Bai and QHL analyzed and interpreted data. GCG obtained the funding. JL is a major contributor in drafting the manuscript. All authors read and approved the final version of the manuscript.
Acknowledgments
This work was supported in part by the National Natural Science Foundation of China (81703667), Basic Public Welfare Research Project of Zhejiang Province (LGN21C150004 and LGF21H280007), Jiaxing Science and Technology Project (2020AY10023), and Special Project of Jiaxing Science and Technology Commissioner (2021K117).