Skip to main content
Download PDF

Genome-wide development of simple sequence repeats markers and genetic diversity analysis of chayote

BMC Plant Biology volume 24, Article number: 603 (2024) Cite this article

Abstract

Background

Chayote is a high economic crop in the Cucurbitaceae family, playing an important role in food production, disease treatment and the production of degradable materials in industries. Due to the harsh environment, such as high temperature, drought and frost, some chayote resources are gradually disappearing. It is crucial to collect, characterize, and conserve chayote resources. However, the genetic diversity of chayote resources in China has not been studied so far.

Results

In this study, we collected 35 individuals of chayote from 14 provinces in China. Subsequently, we found 363,156 SSR motifs from the chayote genome and designed 57 pairs of SSR primers for validation. Out of these, 48 primer pairs successfully amplified bands, with 42 of them showing polymorphism. These 42 primer pairs detected a total of 153 alleles, averaging 3.64 alleles per locus. The polymorphic information content ranged from 0.03 to 0.78, with an average value of 0.41, indicating a high level of polymorphism. Based on the analysis using STRUCTURE, PCoA, and UPGMA methods, the 35 chayote individuals were divided into two major clusters. Through further association analysis, 7 significantly associated SSR markers were identified, including four related to peel color and three related to spine.

Conclusions

These molecular markers will contribute to the analysis of genetic diversity and genetic breeding improvement of chayote in the future.

Peer Review reports

Introduction

Chayote (Sechium edule) is a perennial climbing plant belonging to the Cucurbitaceae family. It originated in Mexico and is now also cultivated in southwest China and several coastal provinces [1]. Chayote was introduced to China during the late Ming Dynasty (1368–1644) or early Qing Dynasty (1644–1912). The exact timeline is debated, but it is generally believed that Portuguese or Spanish traders brought the plant to Asia. It quickly adapted to the subtropical and tropical climates of southern China. In China, chayote is known as ‘fó shǒu guā’. It has become a staple in southern Chinese cuisine, appreciated for its mild flavor and versatility in dishes [2]. The fruits are primarily white or green and may have a variable number of spines on their surface [3]. Chayote is a dual-use plant. Its fruit is rich in vitamins, protein, cucurbitacin, phenolic compounds, and other nutrients beneficial to humans, including trace elements like zinc and selenium [4]. Additionally, chayote has been noted for its potential benefits in treating cardiovascular diseases, lowering blood pressure, and its possible anticancer properties [5]. In recent years, chayote extracts have been explored as a source for new biodegradable plastic products in various industries [6]. Additionally, its root tubers and tendrils are edible [7]. The root tubers are rich in starch and can serve as a substitute for potatoes and sweet potatoes as a primary source of starch [8].

The precocious germination characteristics of chayote make traditional conservation methods infeasible, thereby limiting efforts to study and protect these resources [9]. Therefore, implementing field protection measures is essential for the conservation of chayote. Additionally, chayote, an underutilized crop in the Cucurbitaceae family, typically grows in poor soil conditions [10]. However, the intensification of frost, high temperatures, drought and root diseases in recent years has led to a gradual decline in wild chayote resources. In some areas, such as Mexico and Costa Rica, these resources have decreased by half [11]. The Germplasm National Bank of Sechium edule (BANGESe) has been established in Mexico, with most accessions originating from backyard areas. Reproductive isolation of most wild chayote materials has been influenced by consumer preferences [12]. To date, China has not carried out a systematic germplasm collection of chayote. Although chayote germplasm resources are gradually decreasing, collecting them helps maintain genetic diversity. Protecting these genetic resources is crucial for future crop improvement and ensuring food security. Thus, it is imperative to accelerate the collection of chayote germplasm resources in China [13].

Studying the genetic diversity of species aids in exploring plant evolution history, adaptation potential, and promoting genetic resource conservation. Research has shifted from enzyme-based to various DNA-based molecular markers. SSR molecular markers are considered ideal due to their co-dominance, wide distribution, high mutation rate, and high polymorphism [14]. In addition, SSR molecular markers have been widely used in Cucurbitaceae plants, such as melon [15], bitter gourd [16], cucumber [17] and bottle gourd [18]. The genetic diversity of 42 chayote from Costa Rica was analyzed using 8 isozymes, revealing 35 different polymorphic genotypes, indicating high genetic diversity [19]. With the wide application of molecular markers, 36 chayote samples from India were analyzed using 12 directed amplification of minisatellite-region DNA (DAMD) molecular markers, and 97 polymorphic bands were identified, which were finally divided into three groups due to different living environment conditions [11]. Similarly, 28 pairs of RAPD primers and 5 pairs of ISSR primers were used to analyze 74 chayote in the North Eastern Hill (NEH) region of India, and it was found that the germplasms were not grouped according to their geographical distribution. Due to the varying availability periods of chayote fruits across different regions and the fact that chayote is propagated through its fruits, traders are able to transfer popular local varieties from one planting area to another with ease [10]. Furthermore, 223 pairs of AFLP molecular markers were utilized to analyze 133 accessions of chayote resources from Mexico, revealing a high degree of genetic polymorphism between cultivated germplasm and their wild ancestors [12]. Moreover, the study showed that eight pairs of polymorphic SSR primers were screened from 20 chayote in Mexico by using 11 SSR markers developed, which will be helpful to characterize the population genetic structure of chayote varieties in Mexico and Latin America and decipher the dynamics of their genetic diversity [20]. Subsequently, the researchers used 6 pairs of primers to analyze the genetic diversity of 21 chayote samples from Japan [21].

In this study, chayote resources from the main provinces in China were collected and preserved. Additionally, we identified a rich set of Simple Sequence Repeats (SSR) molecular markers from the chayote genome data and developed reproducible polymorphic SSR markers to assess the genetic diversity and genetic structure of the collected wild germplasm resources. These findings will serve as valuable references for the conservation of chayote resources, genetic improvement, and molecular marker-assisted breeding.

Results

SSR analysis in chayote genome

A total of 363,156 SSRs, with a relative abundance of 597.12/Mb, were identified in the genome sequence (606.42 Mb) of S. edule. Additionally, SSRs were found to have a total length of 5,518,289 bp and an average length of 15.2 bp. They exhibited a relative density of 9,073.52 bp/Mb. The distribution of SSR motifs with higher frequency on 14 chromosomes is summarized in Fig. 1. The information of the SSR loci’ motifs and their positions on chromosomes was also counted (Supplementary Table 2). We have designed primers for all SSR loci and successfully designed a total of 185,137 pairs of primers. We have conducted statistical analysis on the sequence and expected PCR fragment size of these primers (Supplementary Table 3).

Fig. 1
figure 1

Overall view of SSR motif distribution in the chayote reference genome. A chayote chromosomes; B gene density; C-H Mono-, Di-, Tri-, Tetra-, Penta-, and Hexa- repeats, respectively

Full size image

Among all identified SSRs, the number decreased as repeat units increased, with the largest number of repeats at 10 and 11 (Fig. 2A; Supplementary Table 4). The mononucleotide repeats were the most abundant in all identified nucleotides (247,753, 68.22%), followed by dinucleotide repeats (84,342, 23.22%) and trinucleotide repeats (27,074, 7.46%). The number and proportion of SSRs with long motifs were significantly reduced. The repeats of tetranucleotide, pentanucleotide and hexanucleotide were 2,349 (0.65%), 1,020 (0.28%) and 618 (0.17%), respectively (Fig. 2A). The number of SSRs was positively correlated with chromosome length (r = 0.885, p < 0.001), with the most SSRs on chromosome 1 and the fewest on chromosome 14 (Fig. 2B). The most abundant repetitive motif in mononucleotides was the A motif (98.8%). Among dinucleotides, AT and AG repeat motifs were the most abundant, accounting for 49.2% and 25.48%. The most abundant motifs were AAT in trinucleotides (54.16%), AAAT in tetranucleotides (53.23%), AAAAT in pentanucleotides (41.67%), and AAATTT in hexanucleotides (40%) (Fig. 2C).

Fig. 2
figure 2

Distribution characteristics of SSR in the chayote genome. A Number of different SSR types; B Correlation between SSR number and chromosome length; C Frequency distribution of different SSR motif types

Full size image

Comparative genomics analysis between chayote and other cucurbits species

We compared and analyzed the collinearity of universal SSR markers between chayote and different cucurbit species (Fig. 3). The results showed that chayote and melon had the most common SSR markers, with 1,673 (Fig. 3B), followed by cucumber with 1,264 (Fig. 3A), and pumpkin with 1,061(Fig. 3C). Watermelon had the fewest common markers, with 988 (Fig. 3D). In the collinearity between chayote and cucumber, the most common SSR markers were found on chromosome 11 of chayote, with a total of 166 interspecific SSR markers, followed by 149 SSR markers on chromosome 5. The minimum number of SSR markers found on chromosome 9 and chromosome 2 was 17 and 18, respectively (Fig. 3A; Supplementary Table 5). Among them, chromosome 6 of chayote has the least corresponding relationship with chromosomes 3, 5, and 7 of cucumber, while chromosomes 1, 3–5, 7–8, and 10–14 are collinear with all 7 chromosomes of cucumber (Supplementary Table 5). Among the collinearity between chayote and melon, the most collinearity occurred on chromosome 8 of chayote, with 163 shared SSR markers, while only 39 shared SSR markers were found on chromosome 2. The chromosomes 4, 7, and 14 of the chayote and 12 pairs of chromosomes in the melon exhibit collinearity (Supplementary Table 6). While chayote chromosomes 7 and 8 share 193 and 161 SSR markers with 20 chromosomes of pumpkin, respectively, and chayote chromosome 10 only has a minimum of 10 SSR markers with pumpkin chromosomes 5 and 9 (Supplementary Table 7). There were 159 and 20 common SSR markers on chromosomes 12 and 9 of chayote, which were the most and least common with watermelon (Supplementary Table 8).

Fig. 3
figure 3

Collinearity analysis of common SSR molecular markers in chayote and cucurbits species crops. A Sechium edule with Cucumis sativus. B Sechium edule with Cucumis melo. C Sechium edule with Cucurbita moschata. D Sechium edule with Citrullus lanatus. Black represents the chromosome of chayote

Full size image

Analysis of genetic diversity

We randomly selected 57 pairs of SSR primers distributed on 14 chromosomes of chayote for PCR amplification. The PCR amplification results showed that 48 out of 57 primer pairs amplified bands, of which 42 primer pairs had polymorphic bands. These 42 primer pairs detected a total of 153 alleles in 35 chayote resources, with an average amplification of 3.64 loci per primer pair (Supplementary Table 9). The observed variation in the number of alleles (Na) ranges from 2 to 9, with SSR61 having the highest Na value of 9. The variation range of effective allele number (Ne) ranged between 1.03 to 5.23, with an average value of 2.289. The variation ranges of observed heterozygosity (Ho) and expected heterozygosity (He) were 0.17 to 0.97 and 0.03 to 0.83, respectively. In addition, the Shannon Information Index (I) ranged from 0.08 to 1.80, with a mean of 0.85. The Nei's gene diversity index (H) ranged from 0.03 to 0.81, with a mean of 0.48. The variation range of polymorphism information content (PIC) was 0.03 to 0.78, with a mean of 0.41. Among them, 14 SSR primers were highly polymorphic (PIC ≥ 0.5), 20 SSR primers were moderately polymorphic (0.25 < PIC < 0.5), and 8 primers were low polymorphic (PIC ≤ 0.25) (Table 1).

Table 1 Genetic diversity parameters of 42 SSR markers in 35 chayote materials
Full size table

In order to determine whether the difference in the bands of polyacrylamide gel electrophoresis results is due to the difference in the number of SSR repeats. We selected SSR267366 molecular marker and randomly recovered the target bands of PCR products from four materials (No. 1, 12, 27 and 28) for sequencing. The results showed that No. 1, No. 12 and No. 27 had the same number of AG repeat motifs, while No. 28 had fewer repeat motifs than materials 1, 12 and 27 (Fig. 4B). We observed the same results on polyacrylamide gel (Fig. 4A). Multiple sequence alignment further confirmed that all four materials contained conserved AG motifs (35–60 bp), and the differences in bands were due to differences in the number of SSR motif repeats among different chayote resources (Fig. 4C).

Fig. 4
figure 4

Representative SSR polyacrylamide gel electrophoresis and sanger sequencing. A Results of PCR amplification profile of 35 chayote germplasm resources using SSR151039 and SSR267366. M: DNA Ladder; 1 to 35: materials of the 35 chayote individuals listed in Table 3. B Sequencing peak map of SSR267366 marker in chayote with material numbers 1, 12, 27 and 28. C: Multiple sequence alignment results of SSR267366 marker in chayote with material numbers 1, 12, 27 and 28. SSR repeat motifs were marked with black lines

Full size image

Cluster analysis based on SSR markers

The genetic similarity coefficient of 35 germplasms was calculated based on SSR data, and UPGMA cluster analysis was performed using genetic distance (Fig. 5). The genetic similarity coefficients ranged from 0.45 to 0.97, with a variance of 0.52. Materials numbered 14 and 15 had the highest similarity coefficient of 0.97, indicating the closest genetic relationship. At a genetic coefficient of 0.476, the 35 materials were divided into two major clusters. The first cluster was divided into two sub-clusters. The first sub-cluster included No.12 chayote from Hubei Province, and the remaining 9 materials from Sichuan Province. Additionally, materials numbered 22, 25, and 27, all with white fruit skin, formed a separate small cluster. The second sub-cluster included 6 chayote materials (numbered 29–34) from Yunnan Province. The remaining 19 materials formed the second major cluster.

Fig. 5
figure 5

Dendrogram of 35 chayote individuals based on UPGMA

Full size image

Principal coordinate analysis

Principal coordinate analysis (PCoA) based on genetic distance was performed to observe the genetic relationship among 35 chayote genotypes (Fig. 6). The first and second axes explained 16.82% and 13.82% of the variation, respectively. The 35 chayote samples were divided into two clusters. The 35 chayote samples were divided into two clusters. The first cluster included 16 samples (1, 12, and 21–34), while the second cluster included 19 samples (2–11, 13–20, and 35).

Fig. 6
figure 6

Principal coordinates analysis of 35 chayote based on SSR data

Full size image

Population structure analysis

Structure software was used to infer population structure of the 35 chayote samples (Fig. 7). The results show that when K = 2, Δ K has a clear highest peak. Therefore, the 35 materials were divided into two main clusters. The population structure classification results were consistent with the PCoA classification results. Among the 16 materials in the first cluster, 9 were sourced from Sichuan Province (numbered 1 and 21–28), 6 were sourced from Yunnan Province (numbered 29–34), and 1 was sourced from Hubei Province (numbered 12). The 19 materials for the second cluster were sourced from 11 provinces, respectively.

Fig. 7
figure 7

The population structure of 35 chayote germplasms. A Line chart of K value and ΔK value based on SSR markers; B Population structure analysis of 35 chayote germplasm materials (K = 2)

Full size image

Trait (peel color and spine)—SSR association analysis

We conducted an association analysis between fruit peel color and spine with SSR molecular markers by taking Q (population structure) and K (relative kinship) as covariables using MLM and GLM models (Table 2). GLM model showed that three SSR markers were significantly associated with peel color (P < 0.001) and explained 44.71–46.05% of phenotypic variation. SSR108420, SSR185330 and SSR247694 were significantly associated with spine (P < 0.05), which explained the phenotypic variation of 11.94, 16.08 and 11.75% respectively. In MLM model, SSR332709 was significantly associated with peel color (P < 0.05), which explained 23.33% phenotypic variation. SSR108420 and SSR185330 also detected significant associated with spine in MLM model (P < 0.05).

Table 2 Results of association analysis between SSR markers and fruit traits
Full size table

Discussion

Chayote is an edible and medicinal plant distributed in southwest China. Currently, the planting area of chayote is gradually decreasing, leading to the loss of some germplasm resources. It is also an important but underutilized Cucurbitaceae crop [2]. Therefore, collecting and analyzing chayote germplasm resources is crucial for their protection and utilization. This study collected germplasm resources from the main chayote growing areas in China, which are the most reported current germplasm resources of chayote. These germplasm resources will provide considerable opportunities for genetic improvement and breeding in the future.

SSRs are reliable genetic markers for genetic diversity research and breeding. Few studies have developed and applied SSR markers to chayote. Initially, only 11 SSR markers were developed from AC and GA-enriched genomic libraries [20]. Subsequently, six SSR markers were used to analyze the genetic diversity of chayote [21]. With the release of chayote genome, it is possible to develop SSRs at the whole genome level [22]. In this study, a total of 363,156 SSRs were identified in chayote, more than 296,360 SSRs in watermelon (97,103), 254,243 SSRs in melon (DHL92), and 188,091 SSRs in bitter gourd (Dali-11). [16]. At the same time, the SSR marker density of chayote is 597.12 loci/Mb, which is lower than that of bitter gourd (639.73 loci/Mb), cucumber (658.72 loci/Mb), melon (624.78 loci/Mb) and watermelon (834.24 loci/Mb) [16]. This situation may be due to the difference of species, as well as the different SSR identification software and parameters used by developers. The parameters set by SSR of these Cucurbitaceae species are that trinucleotide (Tri-), tetranucleotide (Tetra-), pentanucleotide (Penta-) and hexanucleotide (Hexa-) are repeated four times, while the standard set in this study is at least 5 times. Previous studies have shown that the more repetitions, the higher the polymorphism [23]. We identified 124,177, 221,986, 260,094, and 163,710 SSR motifs in the genomes of cucumber, melon, pumpkin, and watermelon, respectively, using the same parameters. All of these numbers are lower than the number of SSR motifs found in the chayote genome. The SSR marker density in chayote is higher than that in cucumber (549.03 loci/Mb) but lower than that in melon (620.53 loci/Mb), pumpkin (722.88 loci/Mb), and watermelon (622.53 loci/Mb). This situation may be attributed to differences in genome size directly affecting marker density, uneven distribution of SSR markers within the genome, variations in the number and types of repeat sequences among different plants, evolutionary and breeding histories influencing SSR marker distribution, and varying quality of genome sequencing and annotation across crops.

Previous studies have shown that SSR will be enriched at the end of chromosome [24], and SSR motifs of A and AT account for a higher proportion in mononucleotide and dinucleotide [25]. These characteristics can be explained by the fact that there are fewer genes in the centromere region and the downstream genes are rich in poly(A) structure [26]. In this study, the proportion of mononucleotide and dinucleotide in chayote was up to 91.44%, which was consistent with the reported Cucurbitaceae crops [16]. Notably, the number of SSR motifs of A and AT is much higher than that of C and GC. This is consistent with many reported plants, such as Arabidopsis thaliana [27], Akebia trifoliata [28], Monochasma savatieri [29]. This may be attributed to the long evolutionary history of dicotyledons and the higher frequency of AT motifs [30]. Moreover, we conducted comparative genomic analysis of chayote and cucurbit species to explore the collinearity of shared SSR markers among different species. The collinearity results revealed that the genomes of chayote, cucumber, melon, pumpkin, and watermelon were relatively conserved. Comparative genomics has important value in studying the evolutionary process between different cucurbit species.

Studying the genetic diversity of species is not only helpful to explore its adaptability and evolution history, but also provides an important theoretical basis for the protection and utilization of species [31]. In this study, 48 pairs of polymorphic primers were screened from 57 pairs of SSR primers designed in the genome of chayote. The average values of observed heterozygosity and expected heterozygosity among 35 chayote germplasm were 0.52 and 0.48, respectively, and the average value of Shannon information index was 0.85. Previous studies showed that the average values of Ho and He in 21 samples of chayote were 0.29 and 0.41 respectively, and the average value of Shannon information index was 0.61 [21]. The results showed that the genetic differentiation and polymorphism of chayote germplasm resources collected in this study were high. Moreover, the observed heterozygosity is slightly higher than the expected heterozygosity, indicating that there is a surplus of homozygotes among chayote resources. The level of genetic diversity of species is closely related to their mating system, and the observed heterozygosity of mainly outcrossing species is obviously less than the expected heterozygosity [32]. Although chayote is an insect-borne plant, natural geographical barriers such as mountains and lakes lead to inbreeding among individuals of chayote, resulting in too many homozygotes. The high expected heterozygosity among chayote germplasm resources indicates that these resources are rich in genetic diversity, which can be properly protected.

Current studies indicate that factors influencing species genetic diversity include reproductive mode, natural selection, gene flow, genetic drift, and gene mutation [33]. Therefore, genetic diversity analysis is crucial for accurately understanding the genetic relationships between germplasm [34]. In this study, 35 chayote germplasm resources were divided into two clusters based on the three kinds of analysis of Structure, PCoA, and UPGMA. UPGMA clustering analysis divided 35 materials into two clusters, with the first cluster mainly including materials from Yunnan and Sichuan provinces. Notably, three portions of chayote with white peel color were grouped into a small group. Materials from similar sources were also grouped separately into a small group. The structure analysis showed that the materials in Sichuan and Yunnan were divided into one cluster, and the remaining cluster mainly included coastal provinces and some provinces in the central region. The chayote materials showed a certain regional distribution pattern. PCoA analysis once again obtained the same results. In general, geographical isolation may be the main factor affecting the genetic differentiation of chayote population, leading to a high level of genetic differentiation and variation of the species.

Peel color and spine are important agronomic and commercial characteristics that affect the appearance, quality, and value of fruits. They also play a crucial role in plant defense against biotic and abiotic stresses [35]. Chayote fruit is mainly white or green, and the presence or absence of spines is also an important characteristic [4]. Association analysis can combine polymorphic allelic variation sites with target phenotypic data, and can directly identify gene sites closely related to phenotypic variation. However, the accuracy of correlation analysis results is influenced by many factors [36]. The population structure is an important factor affecting the results of association analysis. The mixing of subgroups makes the LD intensity of the whole population increase, which may lead to the association between unlinked gene polymorphism sites and traits, thus obtaining false positive results. Therefore, the simpler the population structure, the more reliable the association analysis results [37]. In this study, 35 chayote materials can be divided into two clusters, and the population structure is relatively simple, which ensures the accuracy of the correlation analysis results to some extent. Moreover, GLM and MLM models are selected for correlation analysis to ensure the reliability of the correlation results. A total of 4 SSR loci were detected which were significantly associated with peel color and 3 with spine, among which SSR108420 and SSR185330 were detected in both GLM and MLM models. At present, the research on correlation analysis of chayote has not been reported, so the SSR markers significantly correlated in this study need further experiments to verify the accuracy of the results.

Conclusion

In this study, 363,156 SSR motifs were found at the genome level of chayote for the first time. Among them, 48 highly polymorphic SSR markers divided 35 chayote resources into two clusters. Geographic isolation plays an important role in genetic differentiation between populations. Moreover, 4 and 3 SSR markers were significantly associated with peel color and spine traits. These SSR markers and population clustering results will play an important role in the future genetic improvement and breeding plans of chayote.

Materials and methods

Plant materials and DNA extraction

A total of 35 chayote genotypes collected from 14 provinces in China were used in this study (Table 3, Fig. 8). There were 9 accessions from different regions in Sichuan Province, 6 from Yunnan Province, 3 from Guizhou, Guangxi, and Jiangsu provinces, 2 from Anhui and Shandong provinces, and 1 from Fujian, Guangdong, Hunan, Hubei, Shaanxi, Shanghai, and Zhejiang provinces. All materials were planted in the modern agricultural research and development base of Sichuan Agricultural University (Chongzhou, Sichuan). The leaves of three samples from each material were mixed in equal quantities and stored at -80℃ for DNA extraction. Genomic DNA was extracted by CTAB method [38]. The quality and concentration of DNA were determined by Nanodrop One ultra-micro ultraviolet spectrophotometer, and the final concentration of DNA was adjusted to 100 ng/μL for subsequent experiments.

Table 3 Details and sources of the 35 chayote samples analyzed in this study
Full size table
Fig. 8
figure 8

Phenotypic photos of 35 chayote fruits germplasm resources

Full size image

Genome-wide SSR mining and primer development

Chayote genome data was sourced from the CuGenDBv2 genome website (http://cucurbitgenomics.org/v2/). Whole-genome SSRs of chayote were detected using Krait software [39], with minimum repeat numbers of ten for mononucleotides (Mono-) and six for dinucleotides (Di-). Five repeats were set for trinucleotides (Tri-), tetranucleotides (Tetra-), pentanucleotides (Penta-), and hexanucleotides (Hexa-). 57 SSR markers distributed across 14 chromosomes of the chayote genome were randomly selected, and primers were designed using Primer 5.0 software (Supplementary Table 1). The specific parameters were: PCR product length of 100–250 bp, primer length of 18–24 bp, annealing temperature of 50–60℃, and GC content of 40–60% [40].

Comparative genome between chayote and other cucurbits species

To analyze the collinearity of the developed SSR markers in the genomes of Cucurbitaceae species, we conducted a BLAST search on the genomes of cucumbers, cantaloupes, pumpkins, and watermelons in the Cucurbitaceae family using all identified repetitive motifs from the chayote genome, with an expected value (E-value) of 10. We allowed primers with a maximum of five nucleotide mismatches at the 5' end, but no mismatches at the 3' end, and a total matching homology of at least 90%. We removed the duplicated parts and the comparison on the scaffolds from the comparison results. Finally, we analyzed the collinearity of SSR markers between species using Circos v0.55 software.

PCR amplification and polyacrylamide gel electrophoresis

We used a 20 μL reaction mixture for PCR amplification, specifically including 10 μL 2 × SanTaq PCR Mix (Shanghai Biotech. Shanghai, China), 1 μL primer and genomic DNA, 7 μL ddH2O. The PCR amplification reaction procedure is as follows: 94℃ pre-denaturation for 5 min; The following 35 cycles: 94℃ denaturation for 30 s, 55℃ annealing for 30 s, and 72 ℃ extension for 30 s; Finally, extension at 72℃ for 10 min. PCR products were detected by using 10% polyacrylamide gel containing 1 × TBE buffer, and electrophoresis was carried out at 180 V for 2 h. After electrophoresis, silver staining was performed and photos were recorded [41].

Cloning and sequencing of PCR products

In order to confirm whether the difference in the banding pattern is caused by the difference in the number of SSR repeats, we randomly selected an SSR marker to randomly cut the amplified PCR products of four chayote resources from the 10% polyacrylamide gel, and recycled the target DNA through the EZ-10 Spin Column DNA PAGE Gel Extraction Kit (Sangon Biotech, China), then connected it to the clone T vector to transform E. coli, and finally selected a single clone for Sanger sequencing. Multi-sequence alignment was carried out by using MEGA software to confirm the difference of SSR motifs numbers.

Data analysis

The basic parameters of each SSR marker were counted by Power Marker (v3.25) software, including observed allele number (Na), effective allele number (Ne), observed heterozygosity (Ho), expected heterozygosity (He), Shannon's information index (I), Nei's gene diversity index (H) and polymorphism information index (PIC) [42]. Principal component analysis was conducted using GenAIEx (v6.51b) software. Structure (v2.3.4) software was used to analyze the population genetic structure. Set the ancestral value from 1 to 10, and assume that all loci are independent. Set the initial uncounted iteration of MCMC to 10,000 times, and then set the uncounted iteration to 100,000 times and 10 repetitions. Determine the best k value by using the online evaluation software Structure Harvester (Web v0.6.94) [43]. The software NTsys-pc (v2.10 e) was used for cluster analysis among varieties based on unweighted pair group with arithmetic average (UPGMA). The relationship between SSR markers and the phenotypic traits of chayote peel color and spine was analyzed by using the mixed linear model (MLM) and general linear model (GLM) of Tassel (v2.1) software [44].

Availability of data and materials

All data generated or analyzed during this study are included in the manuscript and its Supplemental Table 1 to Supplemental Table 8. The materials used during the current study are available from the corresponding authors. The SSR makers of chayote were deposited and are available at https://github.com/shaobocheng/SSR-maker.

References

  1. Vieira EF, Pinho O, Ferreira IMPLVO, Delerue-Matos C. Chayote (Sechium edule): A review of nutritional composition, bioactivities and potential applications. Food Chem. 2019;275:557–68.

    Article  CAS  PubMed  Google Scholar 

  2. Pu Y, Luo Q, Wen L, Li Y, Meng P, Wang X, Tan G. Origin, evolution, breeding, and omics of chayote, an important Cucurbitaceae vegetable crop. Front Plant Sci. 2021;12:739091.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Cadena-Iñiguez J, Avendaño-Arrazate CH, Soto-Hernández M, Ruiz-Posadas LM, Aguirre-Medina JF, Arévalo-Galarza L. Infraspecific variation of Sechium edule (Jacq.) Sw. in the state of Veracruz, Mexico. Genetic Resources and Crop Evolution. 2008;55(6):835–47.

    Article  Google Scholar 

  4. Díaz-de-Cerio E, Verardo V, Fernández-Gutiérrez A, Gómez-Caravaca AM. New insight into phenolic composition of chayote (Sechium edule (Jacq.) Sw.). Food Chemistry. 2019;295:514–9.

    Article  PubMed  Google Scholar 

  5. Trejo-Moreno C, Castro-Martínez G, Méndez-Martínez M, Jiménez-Ferrer JE, Pedraza-Chaverri J, Arrellín G, Zamilpa A, Medina-Campos ON, Lombardo-Earl G, Barrita-Cruz GJ, Hernández B, Ramírez CC, Santana MA, Fragoso G, Rosas G. Acetone fraction from Sechium edule (Jacq.) S.w. edible roots exhibits anti-endothelial dysfunction activity. J J Ethnopharmacol. 2018;220:75–86.

    Article  CAS  PubMed  Google Scholar 

  6. Wu H, Wang J, Li T, Lei Y, Peng L, Chang J, Li S, Yuan X, Zhou M, Zhang Z. Effects of cinnamon essential oil-loaded Pickering emulsion on the structure, properties and application of chayote tuber starch-based composite films. Int J Biol Macromol. 2023;240:124444.

    Article  CAS  PubMed  Google Scholar 

  7. Cheng S, Jia M, Su L, Liu X, Chu Q, He Z, Zhou X, Lu W, Jiang C. Genome-wide identification of the MADS-box gene family during male and female flower development in Chayote (Sechium edule). Int J Mol Sci. 2023;24(7):6114.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Su L, Cheng S, Liu Y, Xie Y, He Z, Jia M, Zhou X, Zhang R, Li C. Transcriptome and metabolome analysis provide new insights into the process of Tuberization of Sechium edule roots. Int J Mol Sci. 2022;23(12):6390.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Ramírez-Rodas YC, Arévalo-Galarza MD, Cadena-Iñiguez J, Soto-Hernández RM, Peña-Valdivia CB, Guerrero-Analco JA. Chayote Fruit (Sechium edule var. virens levis) development and the effect of growth regulators on seed germination. Plants. 2023;12(1):108.

    Article  Google Scholar 

  10. Verma VK, Pandey A, Jha AK, Ngachan SV. Genetic characterization of chayote [Sechium edule (Jacq.) Swartz.] landraces of North Eastern Hills of India and conservation measure. Physiology and Molecular Biology of Plants. 2017;23(4):911–24.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Jain JR, Timsina B, Satyan KB, Manohar SH. A comparative assessment of morphological and molecular diversity among Sechium edule (Jacq.) Sw. accessions in India. 3 Biotech. 2017;7(2):106.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Cadena-Iñiguez J, Avendaño-Arrazate CH, Arévalo-Galarza MD, Cisneros-Solano VM, Ruiz-Posadas LD, Aguirre-Medina JF, Watanabe K, Machida-Hirano R, Barrera-Guzmán LA. Varietal descriptors for the distinction of underutilized varieties of Sechium edule (Jacq) Swartz. Plants. 2022;11(23):3309.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Cheng S, Liu Y, Su L, Liu X, Chu Q, He Z, Zhou X, Lu W, Jiang C, Zheng W. Physiological, anatomical and quality indexes of root tuber formation and development in chayote (Sechium edule). BMC Plant Biol. 2023;23(1):413.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Wang R, Zhong Y, Hong W, Luo H, Li D, Zhao L, Zhang H, Wang J. Genetic diversity evaluation and core collection construction of pomegranate (Punica granatum L.) using genomic SSR markers. Scientia Horticulturae. 2023;319:112192.

    Article  CAS  Google Scholar 

  15. Zhang J, Yang J, Lv Y, Zhang X, Xia C, Zhao H, Wen C. Genetic diversity analysis and variety identification using SSR and SNP markers in melon. BMC Plant Biol. 2023;23(1):39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Cui J, Cheng J, Nong D, Peng J, Hu Y, He W, Zhou Q, Dhillon NPS, Hu K. Genome-wide analysis of simple sequence repeats in bitter gourd (Momordica charantia). Front Plant Sci. 2017;8:1103.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Fukino N, Yoshioka Y, Kubo N, Hirai M, Sugiyama M, Sakata Y, Matsumoto S. Development of 101 novel SSR markers and construction of an SSR-based genetic linkage map in cucumber (Cucumis sativus L.). Breed Sci. 2008;58(4):475–83.

    Article  CAS  Google Scholar 

  18. Bonthala B, Abdin MZ, Arya L, Pandey CD, Sharma V, Yadav P, Verma M. Genome-wide SSR markers in bottle gourd: development, characterization, utilization in assessment of genetic diversity of National Genebank of India and synteny with other related cucurbits. J Appl Genet. 2022;63(2):237–63.

    Article  CAS  PubMed  Google Scholar 

  19. Abdelnour A, Rocha OJ. Genetic characterization of a collection of chayote, Sechium edule (Jacq.) Swartz, in Costa Rica by using isozyme markers. Genetic Resources and Crop Evolution. 2008;55(1):163–70.

    Article  Google Scholar 

  20. Machida-Hirano R, Cortés-Cruz M, Amaro González B, Íñiguez J, Shirata K, Watanabe K. Isolation and characterization of novel microsatellite markers in chayote (Sechium edule (Jacq.) Sw.). Am J Plant Sci. 2015;6:2033–41.

    Article  CAS  Google Scholar 

  21. Shi M, Wang Y, Olvera-Vazquez SG, Iñiguez JC, Thein MS, Watanabe KN. Comparison of chayote (Sechium edule (Jacq.) Sw.) accessions from Mexico, Japan, and Myanmar using reproductive characters and microsatellite markers. Plants. 2023;12(3):476.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Fu A, Wang Q, Mu J, Ma L, Wen C, Zhao X, Gao L, Li J, Shi K, Wang Y, Zhang X, Zhang X, Wang F, Grierson D, Zuo J. Combined genomic, transcriptomic, and metabolomic analyses provide insights into chayote (Sechium edule) evolution and fruit development. Horticulture Res. 2021;8(1):35.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Bansal S, Kumar A, Lone AA, Khan MH, Malhotra EV, Singh R. Development of novel genome-wide simple sequence repeats (SSR) markers in Bunium persicum. Ind Crops Prod. 2022;178:114625.

    Article  CAS  Google Scholar 

  24. Liu L, Wang H, Li Y, Chen S, Wu M, Dou M, Qi Y, Fang J, Zhang J. Genome-wide development of interspecific microsatellite markers for Saccharum officinarum and Saccharum spontaneum. J Integr Agric. 2022;21(11):3230–44.

    Article  CAS  Google Scholar 

  25. Shi Z, Zhao W, Li Z, Kang D, Ai P, Ding H, Wang Z. Development and validation of SSR markers related to flower color based on full-length transcriptome sequencing in Chrysanthemum. Sci Rep. 2022;12(1):22310.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Oko Y, Ito N, Sakamoto T. The mechanisms and significance of the positional control of centromeres and telomeres in plants. J Plant Res. 2020;133(4):471–8.

    Article  CAS  PubMed  Google Scholar 

  27. Lawson MJ, Zhang L. Distinct patterns of SSR distribution in the Arabidopsis thalianaand rice genomes. Genome Biol. 2006;7(2):R14.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Zhong S, Chen W, Yang H, Shen J, Ren T, Li Z, Tan F, Luo P. Characterization of microsatellites in the Akebia trifoliata genome and their transferability and development of a whole set of effective, polymorphic, and physically mapped simple sequence repeat markers. Front Plant Sci. 2022;13:860101.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Yang W, Bai Z, Wang F, Zou M, Wang X, Xie J, Zhang F. Analysis of the genetic diversity and population structure of Monochasma savatieri Franch. ex Maxim using novel EST-SSR markers. BMC Genomics. 2022;23(1):597.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Cavagnaro PF, Senalik DA, Yang L, Simon PW, Harkins TT, Kodira CD, Huang S, Weng Y. Genome-wide characterization of simple sequence repeats in cucumber (Cucumis sativus L.). BMC Genomics. 2010;11(1):569.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Sood T, Sood S, Sood VK, Badiyal A, Anuradha, Kapoor S, Sood V, Kumar N. Characterisation of bell pepper (Capsicum annuum L. var. grossum Sendt.) accessions for genetic diversity and population structure based on agro-morphological and microsatellite markers. Sci Hortic. 2023;321:112308.

    Article  CAS  Google Scholar 

  32. Mutegi E, Snow AA, Rajkumar M, Pasquet R, Ponniah H, Daunay M-C, Davidar P. Genetic diversity and population structure of wild/weedy eggplant (Solanum insanum, Solanaceae) in southern India: Implications for conservation. Am J Bot. 2015;102(1):140–8.

    Article  PubMed  Google Scholar 

  33. Ivandic V, Hackett CA, Nevo E, Keith R, Thomas WTB, Forster BP. Analysis of simple sequence repeats (SSRs) in wild barley from the Fertile Crescent: associations with ecology, geography and flowering time. Plant Mol Biol. 2002;48(5):511–27.

    Article  CAS  PubMed  Google Scholar 

  34. Nordlander J, Skytte af Sätra J, Mattisson H, Udovychenko K, Lushpigan O, Lose L, Naess H, Leino M, Hjalmarsson I, Garkava-Gustavsson L. Genetic diversity in gooseberry (Ribes uva-crispa), as estimated with SSR markers. Scientia Horticulturae. 2022;306:111438.

    Article  CAS  Google Scholar 

  35. Che G, Zhang X. Molecular basis of cucumber fruit domestication. Curr Opin Plant Biol. 2019;47:38–46.

    Article  CAS  PubMed  Google Scholar 

  36. Luan M, Liu C, Wang X, Xu Y, Sun Z, Chen J. SSR markers associated with fiber yield traits in ramie (Boehmeria nivea L. Gaudich). Industrial Crops and Products. 2017;107:439–45.

    Article  CAS  Google Scholar 

  37. Liao G, Li Z, Huang C, Zhong M, Tao J, Qu X, Chen L, Xu X. Genetic diversity of inner quality and SSR association analysis of wild kiwifruit (Actinidia eriantha). Sci Hortic. 2019;248:241–7.

    Article  CAS  Google Scholar 

  38. Edwards K, Johnstone C, Thompson C. A simple and rapid method for the preparation of plant genomic DNA for PCR analysis. Nucleic Acids Res. 1991;19(6):1349–1349.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Du L, Zhang C, Liu Q, Zhang X, Yue B. Krait: an ultrafast tool for genome-wide survey of microsatellites and primer design. Bioinformatics. 2018;34(4):681–3.

    Article  CAS  PubMed  Google Scholar 

  40. Guo Y, Liu L, Yue Y, Fan X, Teng W, Zhang H, Gao K, Guan J, Chang Z, Teng K. Development of SSR markers based on transcriptome sequencing and verification of their conservation across species of ornamental pennisetum rich. (Poaceae). Agronomy. 2022;12(7):1683.

    Article  CAS  Google Scholar 

  41. Lu X, Adedze YMN, Chofong GN, Gandeka M, Deng Z, Teng L, Zhang X, Sun G, Si L, Li W. Identification of high-efficiency SSR markers for assessing watermelon genetic purity. J Genet. 2018;97(5):1295–306.

    Article  CAS  PubMed  Google Scholar 

  42. Zhu H, Guo L, Song P, Luan F, Hu J, Sun X, Yang L. Development of genome-wide SSR markers in melon with their cross-species transferability analysis and utilization in genetic diversity study. Mol Breeding. 2016;36(11):153.

    Article  Google Scholar 

  43. Yin J, Zhao H, Wu X, Ma Y, Zhang J, Li Y, Shao G, Chen H, Han R, Xu Z. SSR marker based analysis for identification and of genetic diversity of non-heading Chinese cabbage varieties. Front Plant Sci. 2023;14:1112748.

    Article  PubMed  PubMed Central  Google Scholar 

  44. An X, Luo X, Liu T, Li W, Zou L. Development and application of fruit color-related expressed sequence tag-simple sequence repeat markers in Abelmoschus esculentus on the basis of Transcriptome sequencing. Front Plant Sci. 2022;13:907895.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Statement

We have permission to collect all plant materials used in the study.

Funding

This study was financially supported by the Second Tibetan Plateau Scientific Expedition and Research (2019QZKK0303) and Ya 'an Science and Technology Bureau.

Author information

Author notes

  1. Shaobo Cheng, Lihong Su and Xin Guo contributed equally to this work and share the first authorship.

Authors and Affiliations

  1. College of Horticulture, Sichuan Agricultural University, Chengdu, 611130, PR China

    Shaobo Cheng, Lihong Su, Xin Guo, Dalong Shao, Yanmei Qin, Xuanxuan Liu, Qianwen Chu, Xiaoting Zhou & Zhongqun He

Authors
  1. Shaobo Cheng
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Lihong Su
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Xin Guo
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Dalong Shao
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Yanmei Qin
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Xuanxuan Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Qianwen Chu
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Xiaoting Zhou
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Zhongqun He
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

ZH conceived the study and designed the experiments. SC, LS, XG, DS, YQ, XL, QC and XZ carried out the experiments and analyzed the data. SC wrote the manuscript. All authors read and approved the final manuscript

Corresponding author

Correspondence to Zhongqun He.

Ethics declarations

Ethics approval and consent to participate

We declare that all methods in this study were carried out in accordance with relevant guidelines and regulations.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cheng, S., Su, L., Guo, X. et al. Genome-wide development of simple sequence repeats markers and genetic diversity analysis of chayote. BMC Plant Biol 24, 603 (2024). https://doi.org/10.1186/s12870-024-05317-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-024-05317-9

Keywords

BMC Plant Biology

ISSN: 1471-2229

Contact us