Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Genetic Diversity and Population Structure of Cucumber (Cucumis sativus L.)

  • Jing Lv ,

    Contributed equally to this work with: Jing Lv, Jianjian Qi, Qiuxiang Shi, Di Shen, Shengping Zhang

    Current address: Tobacco Research Institute of Chinese Academy of Agricultural Sciences/Key Laboratory of Tobacco Biology and Processing, Ministry of Agriculture, Qingdao, China

    Affiliation Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Sino-Dutch Joint Lab of Horticultural Genomics, Opening Lab of Genetic Improvement of Agricultural Crops of Ministry of Agriculture, Beijing, China

  • Jianjian Qi ,

    Contributed equally to this work with: Jing Lv, Jianjian Qi, Qiuxiang Shi, Di Shen, Shengping Zhang

    Affiliation Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Sino-Dutch Joint Lab of Horticultural Genomics, Opening Lab of Genetic Improvement of Agricultural Crops of Ministry of Agriculture, Beijing, China

  • Qiuxiang Shi ,

    Contributed equally to this work with: Jing Lv, Jianjian Qi, Qiuxiang Shi, Di Shen, Shengping Zhang

    Affiliation Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Sino-Dutch Joint Lab of Horticultural Genomics, Opening Lab of Genetic Improvement of Agricultural Crops of Ministry of Agriculture, Beijing, China

  • Di Shen ,

    Contributed equally to this work with: Jing Lv, Jianjian Qi, Qiuxiang Shi, Di Shen, Shengping Zhang

    Affiliation Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Sino-Dutch Joint Lab of Horticultural Genomics, Opening Lab of Genetic Improvement of Agricultural Crops of Ministry of Agriculture, Beijing, China

  • Shengping Zhang ,

    Contributed equally to this work with: Jing Lv, Jianjian Qi, Qiuxiang Shi, Di Shen, Shengping Zhang

    Affiliation Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Sino-Dutch Joint Lab of Horticultural Genomics, Opening Lab of Genetic Improvement of Agricultural Crops of Ministry of Agriculture, Beijing, China

  • Guangjin Shao,

    Affiliation Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Sino-Dutch Joint Lab of Horticultural Genomics, Opening Lab of Genetic Improvement of Agricultural Crops of Ministry of Agriculture, Beijing, China

  • Hang Li,

    Affiliations Middlebury College, Middlebury, Vermont, United States of America, High School Affiliated to Renmin University of China, Beijing, China

  • Zhanyong Sun,

    Affiliation East-West Seed International Ltd., Nanning, Guangxi, China

  • Yiqun Weng,

    Affiliation United States Department of Agriculture (USDA), ARS, Vegetable Crops Research Unit, Department of Horticulture, University of Wisconsin, Madison, Wisconsin, United States of America

  • Yi Shang,

    Affiliation Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Sino-Dutch Joint Lab of Horticultural Genomics, Opening Lab of Genetic Improvement of Agricultural Crops of Ministry of Agriculture, Beijing, China

  • Xingfang Gu,

    Current address: State Key Lab Crop Stress Biol Arid Areas, Coll Plant Protect, Northwest A&F University, Yangling, Shaanxi, China

    Affiliation Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Sino-Dutch Joint Lab of Horticultural Genomics, Opening Lab of Genetic Improvement of Agricultural Crops of Ministry of Agriculture, Beijing, China

  • Xixiang Li,

    Affiliation Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Sino-Dutch Joint Lab of Horticultural Genomics, Opening Lab of Genetic Improvement of Agricultural Crops of Ministry of Agriculture, Beijing, China

  • Xiaoguo Zhu,

    Affiliation Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Sino-Dutch Joint Lab of Horticultural Genomics, Opening Lab of Genetic Improvement of Agricultural Crops of Ministry of Agriculture, Beijing, China

  • Jinzhe Zhang,

    Current address: National Laboratory for Protein Engineering and Plant Genetic Engineering, Peking-Yale Joint Research Center for Plant Molecular Genetics and AgroBiotechnology, College of Life Sciences, Peking University, Beijing, China

    Affiliation Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Sino-Dutch Joint Lab of Horticultural Genomics, Opening Lab of Genetic Improvement of Agricultural Crops of Ministry of Agriculture, Beijing, China

  • Robbert van Treuren,

    Affiliation Centre for Genetic Resources, The Netherlands, Wageningen University and Research Centre, Wageningen, The Netherlands

  • Willem van Dooijeweert,

    Affiliation Centre for Genetic Resources, The Netherlands, Wageningen University and Research Centre, Wageningen, The Netherlands

  • Zhonghua Zhang,

    Affiliation Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Sino-Dutch Joint Lab of Horticultural Genomics, Opening Lab of Genetic Improvement of Agricultural Crops of Ministry of Agriculture, Beijing, China

  •  [ ... ],
  • Sanwen Huang

    huangsanwen@caas.net.cn

    Affiliation Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Sino-Dutch Joint Lab of Horticultural Genomics, Opening Lab of Genetic Improvement of Agricultural Crops of Ministry of Agriculture, Beijing, China

  • [ view all ]
  • [ view less ]

Abstract

Knowing the extent and structure of genetic variation in germplasm collections is essential for the conservation and utilization of biodiversity in cultivated plants. Cucumber is the fourth most important vegetable crop worldwide and is a model system for other Cucurbitaceae, a family that also includes melon, watermelon, pumpkin and squash. Previous isozyme studies revealed a low genetic diversity in cucumber, but detailed insights into the crop's genetic structure and diversity are largely missing. We have fingerprinted 3,342 accessions from the Chinese, Dutch and U.S. cucumber collections with 23 highly polymorphic Simple Sequence Repeat (SSR) markers evenly distributed in the genome. The data reveal three distinct populations, largely corresponding to three geographic regions. Population 1 corresponds to germplasm from China, except for the unique semi-wild landraces found in Xishuangbanna in Southwest China and East Asia; population 2 to Europe, America, and Central and West Asia; and population 3 to India and Xishuangbanna. Admixtures were also detected, reflecting hybridization and migration events between the populations. The genetic background of the Indian germplasm is heterogeneous, indicating that the Indian cucumbers maintain a large proportion of the genetic diversity and that only a small fraction was introduced to other parts of the world. Subsequently, we defined a core collection consisting of 115 accessions and capturing over 77% of the SSR alleles. Insight into the genetic structure of cucumber will help developing appropriate conservation strategies and provides a basis for population-level genome sequencing in cucumber.

Introduction

Commonly known as cucurbits or gourds, the botanical family Cucurbitaceae includes a number of cultivated species of global or local economical importance [1]. Cucumber (Cucumis sativus L.) is the fourth most important vegetable worldwide [2]. As the first in cucurbits, the cucumber genome sequence elucidated chromosomal evolution in the genus Cucumis and afforded novel insights into several important biological processes such as biosynthesis of cucurbitacin and “fresh green” odor. Cucumber is being developed as a new model species in plant biology due to its small number of genes, rich diversity of sex expression, suitability for vascular biology studies, short life cycle (three months from seed to seed), and accumulating resources in genetics [3], [4] and genomics [5], [6].

The rapid advance of Next Generation Sequencing (NGS) technologies makes it affordable to re-sequence multiple genotypes of a given species to generate a haplotype map that displays the genome-wide patterns of genetic variation at a single base resolution [7], [8], [9]. However, such population sequencing requires a well-designed sampling of the species' genetic diversity. To date, little is known about the genetic diversity and structure in cucumber.

Cucumber is indigenous to India [10] and likely originated from the foothills of the Himalayan Mountain, where its only two botanical varieties both were discovered, namely the domesticated cucumber C.s. var. sativus and the wild cucumber C.s. var. hardwickii (Royle) Alef. The cultivation of cucumber seems to have spread rapidly from India to Western Asia, and then to Southern Europe. It has been reported that cucumber was introduced into China through the Silk Route by the diplomat Zhang Qian (164 B.C.–114 B.C.) during the Han Dynasty and subsequently spread to East Asia [11]. A unique semi-wild landrace, named Xishuangbanna Gourd, was discovered in the Prefecture Xishuangbanna of the Province Yunnan in Southwest China. To date, its relationship to other phylogeographic populations is unknown.

Genetic diversity in the U.S. cucumber collection (∼1,000 Plant Introductions, or PIs) maintained by the National Plant Germplasm System (NPGS) was investigated in a series of studies by analysis of variation at 21 isozyme loci [12], [13], [14], [15]. Overall, these studies revealed low diversity (∼2.2 alleles per polymorphic locus) in domesticated cucumber. Possibly due to the lack of polymorphism, the isozyme data did not allow consistent clustering of the PIs. Consequently, the genetic structure of domesticated cucumber remained largely unknown.

In this study we assembled a mega-collection consisting of all available cucumber accessions from the national germplasm centres of China, the Netherlands and U.S. We fingerprinted a total of 3,342 accessions with 23 highly polymorphic Simple Sequence Repeat (SSR) markers. The SSR data were used to assess overall genetic diversity and to reveal a clear population structure in cucumber. In addition, the data allowed the construction of a core set of 115 accessions capturing over 77% of the total allelic diversity observed.

Materials and Methods

Materials

A mega-collection of 3,342 cucumber accessions was used in this study. A total of 1,692 accessions were collected from the China National Vegetable Germplasm Bank located at the Institute of Vegetables and Flowers at the Chinese Academy of Agricultural Sciences (IVF-CAAS), 883 accessions from the U.S. National Plant Germplasm System (NPGS, http://www.ars-grin.gov/npgs/), 759 accessions from the Centre for Genetic Resources, the Netherlands (CGN, http://www.cgn.wur.nl), 2 accessions from the IPK (http://www.ipk-gatersleben.de/), and 6 accessions from the East-West Seed International Ltd. The mega-collection included 235 accessions originating from India, 1,361 from China (excluding Xishuangbanna), 208 from Xishuangbanna, 209 from Southeast Asia, 141 from Western Asia, 802 from Europe, 176 from the Americas, 50 from Africa and 8 from Oceania. Additionally, 152 accessions of unknown origin were included in the study. Passport information about the mega-collection is listed in Table S1.

SSR analyses

DNA extraction.

For most of the study accessions multiple self-crossing was performed during gene bank curation, and previous studies on cucumber germplasm have revealed limited intra-accession variation. To increase the efficiency in analysing the large number of accessions, leaves from five randomly chosen plants per accession were bulked into a single sample. Genomic DNA was isolated from about one gram of leaf tissue from each bulk using a CTAB-based method [16].

SSR markers.

In previous studies [3], [4], we determined the polymorphism level of 995 SSR markers using a selected panel of 11 cucumber genotypes that represent six market types worldwide: ‘Chinese Long’, Southern China type, Xishuangbanna type, European greenhouse type, American slicing type, and Japanese type. For the present study, we selected 23 highly polymorphic SSR markers that are evenly distributed across the cucumber genetic map (Table 1). Each selected SSR marker is able to produce a clear, stable, and single polymerase chain reaction (PCR) product. On each chromosome arm, there is at least one SSR (Figure 1). The genetic distance between two neighbouring markers ranges from7.1 cM to 53.7 cM with an average of 29.4 cM.

thumbnail
Figure 1. Position of the 23 SSR loci on the cucumber genetic map [18].

https://doi.org/10.1371/journal.pone.0046919.g001

thumbnail
Table 1. Characteristics of the 23 investigated SSR markers and the diversity detected in 3,342 cucumber accessions.

https://doi.org/10.1371/journal.pone.0046919.t001

PCR amplification.

The 23 SSR markers were analyzed using four sets of multiplex PCR reactions. Each multiplex was carefully assembled according to the compatibility of the SSRs during PCR and the molecular size of their amplicons. The forward primer of the SSR markers was labelled with one of the four fluorescent dyes, carboxy fluorescein (FAM), carboxytetramethylrhodamine (TAMRA), hexachloro-6-carboxyfluorescein (HEX) and ROX (carboxy-X-rhodamine). Multiplex PCRs were performed in a 11 µl volume containing approximately 20–50 ng template DNA, 0.4 pmol of each primer, and 5 µl Multiplex PCR Master Mix (QIAGEN Multiplex PCR Kit, Qiagen). Reactions were performed in an ABI 9700 thermocycler with an initial denaturation step of 15 min at 95°C, followed by 35 cycles of 95°C for 30 s, 57°C for 90 s and 72°C for 70 s, and a final extension at 60°C for 30 min. For PCR fragment size determinations, 0.25 µl of an internal size standard (Liz-500, LIZ) was mixed with 1 µl of diluted PCR product (1/100) and 9 µl formamide. The mixture was heated at 94°C for 3 minutes and then cooled within icy water. Electrophoresis was carried out on an ABI 3130×l genetic analyzer (Applied Biosystems). Analyses were performed using the Genemapper 4.0 software.

Data analysis

Genetic structure.

The model-based program STRUCTURE [17] was used to infer population structure and to identify hybrid forms. In order to identify the number of populations (K) capturing the major structure in the data, we used a burn-in period of 10,000 Markov Chain Monte Carlo iterations and 100,000 run length, an admixture model following Hardy-Weinberg equilibrium and correlated allele frequencies as well as independent loci for each run. Three independent runs were performed for each simulated value of K, ranging from 1 to 20. Subsequently, an Pr(X|K) index with respect to each K was used to calculate ΔK using the formula described by Evanno et al. [18]. The optimal K depends on the first peak of ÄΔK = |L″(K)|/s[Pr(x|k)], where (|L″(K)| denotes the absolute value of the second order rate of change of Pr(X|K), and s[Pr(x|k)] the standard deviation of the Pr(X|K)).

Genetic diversity assessment.

For the entire collection as well as the subpopulations, allele richness, gene diversity and the number of unique alleles were calculated with PowerMarker version 3.25 [19], using genetic distances calculated by the CS Chord 1967 method [20]. Tree topologies were constructed based on the neighbour-joining method with MEGA [21]. For each subpopulation inferred from STRUCTURE, the observed and effective number of alleles, the observed and expected heterozygosity, as well as Shannon's information index were calculated by POPGENE. As a measure of genetic differentiation between pairs of sub-populations, we computed FST statistics with the AMOVA means in Arlequin 3.11 [22]. Principal Component Analysis (PCA) in NTSYSpc2.11 was used to analyse genetic relationships among accessions and to determine the optimal number of clusters in the study.

Core collection sampling.

For sampling core collections, we used the Maximization (M) algorithm implemented in MSTRAT software version 4.1 [23]. For evaluation of core collection's minimal size and for individual sampling of the collections, 20 replicates of 20 iterations for each replicate were performed.

Results

Molecular fingerprinting of the mega-collection

We fingerprinted 3,342 cucumber accessions using 23 SSR markers. In total, 94.9% of PCR reactions yielded clear peaks in electrophoresis and generated 72,960 data points. The remaining 5.1% resulted in multiple bands with low quality and consequently were treated as missing values in subsequent analyses. Out of the 72,960 data points, 61,976 (84.9%) corresponded to single peaks that were regarded as homozygous for the given SSR locus, while 8,929 (12.2%) corresponded to double peaks that were regarded as heterozygous. Only 2,055 (2.8%) of the data points corresponded to more than two clear peaks, of which only the two highest peaks were recorded that were regarded as heterozygous. This simplified scoring may cause under-estimation of rare alleles.

The 23 SSR markers generated a total of 316 alleles in the mega-collection (Table 1). Out of the 316 SSR alleles, 64 displayed a frequency of more than 5% in the total sample and hence were classified as ‘common’ alleles, while another 64 displayed frequencies between 1% and 5% that were denoted as ‘less common’ alleles. Out of the 188 remaining alleles, 121 were denoted as ‘rare alleles’ with frequencies between 0.1% and 1% and 67 as ‘very rare’ alleles with frequencies smaller than 0.1%. That 59% of the observed alleles showed a frequency of less than 1% in the mega-collection suggested that this collection represents the bio-diversity of cucumber broadly.

Genetic structure in cucumber

The population structure of the cucumber mega-collection was firstly inferred using STRUCTURE 2.3.1 [18]. The Evanno et al. (2005) correction showing only peak of ΔK, for K = 3, suggested the presence of three main populations in cucumbers (Figure 2A). The classification of accessions into populations by the model-based method was shown in Figure 2B and Table S1. In total 2,932 accessions (87.7%) were assigned to one of the three populations, where more than 70% of their inferred ancestry was derived from one of the model-based populations. Population 1, 2, and 3 (P1, P2, and P3) consisted of 1399, 1129, and 404 accessions, respectively. The remaining 410 accessions (12.3%) were categorized as having admixed ancestry, including 152 admixtures between P1 and P2 (P1P2), 163 between P1 and P3 (P1P3), and 95 between P2 and P3 (P2P3).

thumbnail
Figure 2. Model-based populations in the cucumber mega-population.

(A) Delta K values for different numbers of populations assumed (K) in the STRUCURE analysis. (B) Classification of 3,342 accessions into three populations using STRUCTURE 2.3.1. The distribution of the accessions to different populations is indicated by the color code (P1: red, P2: yellow, P3: blue). Numbers on the y-axis show the subgroup membership, and the x-axis shows the different accession. (C) Unrooted Neighbor-Joining (NJ) tree of the 3,342 accessions. Bootstrap values are indicated at the major branches. (D) Principal Component Analysis (PCA) of the 1,026 accessions without missing SSR fingerprinting data. Color codes for each population and admixture are consistent in Figure 2C, 2D, Figure 3 and Figure 4C.

https://doi.org/10.1371/journal.pone.0046919.g002

The unrooted neighbour-joining tree (Figure 2C) is consistent with the aforementioned model-based population structure. Accessions of the three populations are largely, but not completely, separated. Admixtures are mostly located in between two populations.

We also performed Principal Component Analysis (PCA) on 1,026 accessions without any missing SSR data (Figure 2D). This analysis largely supported the separation of the accessions into three populations. P3 can be divided into two subgroups, one consisting mainly of accessions from India and the other with accessions from Xishuangbanna. Most admixtures appear in between the two populations.

In summary, the model-based ancestry analysis, the phylogenetic tree and the PCA strongly supported that cucumber has three well-differentiated genetic populations and admixtures.

Geographical distribution of cucumber populations

The classification of populations appeared highly correlated with the geographical distribution of the cucumber accessions (Table 2, Figure 3). P1 includes cucumber germplasm from North and East China, Japan and South Korea, P2 includes accessions collected from Central and West Asia, Europe, America and Africa, while P3 comprises 404 accessions mainly from India and Xishuangbanna. Xishuangbanna is geographically close to the Indo-Gangetic plain where cucumber originated. Also the distribution of the admixtures appeared related to geography. About a quarter of the P1P2 admixtures (36 of 152) originated from East Europe including the Russian Federation and the Former Soviet Union. The P1P3 admixtures are mainly from Southwest China and India, and the P2P3 admixtures mainly from India. The admixtures are most likely due to natural or artificial hybridization.

thumbnail
Figure 3. Geographical distribution of the investigated cucumber accessions.

Color codes denote the populations shown in Figure 2.

https://doi.org/10.1371/journal.pone.0046919.g003

thumbnail
Table 2. Geographic origin of the accessions assigned by the software STRUCURE to populations (P1–P3) and admixtures (P1P2–P2P3).

https://doi.org/10.1371/journal.pone.0046919.t002

In PCA analysis (Figure 2D), P3 can be divided into two subgroups, one with accessions from India mainly and the other with accessions from Xishuangbanna. With more markers, the subgroups might be differentiated into two separate populations.

Genetic diversity assessment

AMOVA analysis of the populations P1, P2 and P3 indicated that 68.1% of the variation was due to differences within populations, while 31.9% was due to differences among populations (Table 3). Pairwise estimates of FST using AMOVA indicated a high degree of differentiation between the three model-based populations with values ranging from 0.30 to 0.33 (between P1 and P2: 0.33, P1 and P3: 0.30; P2 and P3: 0.31).

thumbnail
Table 3. Analysis of molecular variance (AMOVA) among populations and within populations.

https://doi.org/10.1371/journal.pone.0046919.t003

Genetic parameters for all accessions and the three model based populations are given in Table 4. For all accessions genotyped, the average number of alleles per locus was 13.70, ranging from 6 (SSR05012) to 25 (SSR13787) (Table 1). The Expected number of alleles was much smaller than that of observed, indicating the existence of few high frequency alleles in cucumber genome. The overall genetic diversity of cucumber was similar to tomato (0.2–0.58) [24] and rice (0.16–0.79) [25], but lower than maize (0.54–0.82) [26].

thumbnail
Table 4. Summary of SSR diversity parameters of the model-based populations and admixtures.

https://doi.org/10.1371/journal.pone.0046919.t004

Among the three model based populations, P3 possessed the highest diversity (He = 0.52), though there was no significant differences in Na and Ne. Even with the least number of accessions involved, P3 displayed the largest number of population specific alleles (n = 40). This supports the notion that cucumber was originated from India [10].

For all populations, the observed heterozygosity was lower than expected heterozygosity, which may be caused by inbreeding during seed propagation by germplasm banks.

Core Collection

The Maximization or M strategy algorithm was used to select the core set of the mega-population. The M curve reached the plateau when approximately 115 accessions were selected (Figure 4A). Taking geographic distribution, existing phenotype data, and population structure into account, we selected 115 accessions to define the core collection (Table S2). The core collection captures 77.2% of the alleles discovered in the mega-collection, including all the 64 common alleles, 61 of the 64 less common alleles, 94 of the 121 rare alleles and 25 of the 67 very rare alleles (Figure 4B). The core collection has higher effective number of alleles, higher expected heterozygosity, and higher gene diversity than the entire collection (Table 5). In general, the topology of the unrooted neighbour-joining tree of the core collection (Figure 4C) resembles the tree of the mega-population, while better separating the populations. The P2P3 admixtures are clustered with Indian P3 accessions, indicating their evolutionary relationship.

thumbnail
Figure 4. Development and evaluation of a core collection.

(A) The minimized number of accessions for a core collection capturing the majority of alleles by random sampling and the Mstrat strategy. (B) Number of alleles observed in the total sample and those captured in the core collection, also presented separately for alleles occurring in different frequency. (C) NJ tree of the 115 core accessions. Values at the branches denote bootstrap values.

https://doi.org/10.1371/journal.pone.0046919.g004

thumbnail
Table 5. Genetic diversity parameters of the entire sample and those of the core collection based on the 23 investigated SSR loci.

https://doi.org/10.1371/journal.pone.0046919.t005

Discussion

In this study, we carefully selected 23 highly polymorphic SSR markers and used them to fingerprint a large cucumber collection. The dataset allowed us to visualize the population structure of cucumber and to assemble a core collection of 115 accessions that cover over 77% of the observed SSR alleles.

Out of the 316 SSR alleles scored, 59% and 80% are present in less than 1% and 5% of the mega-collection, respectively. As the level of gene flow is negatively related to the number of rare alleles [27], the abundance of rare alleles in the mega-collection indicates that there has been limited genetic exchange among sub-populations in cucumber.

The results also supports that independent migrations of cucumber have occurred from India to other parts of the world. The migration route to Europe and then to the Americas and Africa are clearly through Central and West Asian countries. Cucumber accessions from China and other East Asian countries are genetically differentiated from those in Central and West Asia. It is well documented in Chinese history that cucumber was introduced from there through the Silk Route by the diplomat Zhang Qian (164 B.C.–114 B.C.) [11]. In addition, cucumber cultivation in heated houses during the winter season was recorded in poems from the Tang Dynasty (618–907). These records support that the crop has been intensively cultivated and therefore selected by Chinese farmers for more than one thousand years, which likely caused the genetic differentiation of Chinese (East Asian) cucumbers related to other populations.

Cucumber plants in India and Xishuangbanna bear short, ovary to round fruits with sparse spines (Figure S1). The occurrence of black spines, an unfavourable market trait, within Indian resources, is more frequent than the other two populations. Cucumbers should have undergone differential human selection when introduced to other parts of the world. For instance, in Europe and America, short-fruited gherkins and long fruited Dutch cucumbers were both developed to satisfy different market need (Figure S1). To enhance total yield, gynoecious varieties were widely adopted by western farmers as they could bear more fruits [28]. But in East Asia where cucumber is often served for direct consumption, characteristics such as proper appearance, crispness and “fresh green” flavour outweigh other traits (Figure S1). As a result, the ‘Chinese Long’ type had been developed to meet with that specific demand.

The Xishuangbanna cucumbers are affiliated to the Indian population. As semi-wild landraces, they apparently lost the unfavourable fruit bitterness but still possess the dependence of short daylight for flowering and fruit setting. The latter trait was discarded in most cucumber cultivars. Accessions from Xishuangbanna are genetically close to each other, suggesting that a single dispersal occurred from India. It is worth mentioning that among the African accessions, six lines are present from Zambia and Zimbabwe classified as P3 and admixture P2P3 (CG7743–CG7745, CG7747–CG7749, Table S1). The accessions were initially misclassified as melon and likely represent another dispersal from India.

The study also has potential implications for cucumber breeding. Currently, most crosses have been made within populations and less frequently among populations. The large differentiation among populations indicates that each population may possess its own alleles and haplotypes. Therefore, crosses between populations will broaden the genetic diversity within current breeding programs and may result in new types of heterosis. This study provides a sound basis for further characterization of the biodiversity of this important vegetable. For example, we are currently re-sequencing the core collection, enabling in-depth discovery of the genetic variation in cucumber.

Supporting Information

Table S1.

Detail information of cucumber resources.

https://doi.org/10.1371/journal.pone.0046919.s001

(XLSX)

Figure S1.

Pictures of typical cucumber accessions from the three model based populations. CG1149, “Chinese Long”; CG4357, Southern China type; CG5786, Dutch greenhouse; CG6600, Pickling; CG9191, Xishuangbanna; CG8039, Indian

https://doi.org/10.1371/journal.pone.0046919.s003

(PDF)

Acknowledgments

We are grateful to Dr. Thomas Städler at ETH Zurich and to Prof. Susanne Renner at University of Munich for their critical reading. The authors wish to acknowledge NPGS (http://www.ars-grin.gov/npgs/), CGN (http://www.cgn.wur.nl), National Mid-term Genebank for Vegetables (Beijing, China), and IPK (http://www.ipk-gatersleben.de/) for kindly providing cucumber accessions for this study.

Author Contributions

Conceived and designed the experiments: SH. Performed the experiments: QS JQ DS SZ GS HL JL YS XZ JZ. Analyzed the data: JL JQ QS ZZ RT WD ZS YW XG XL. Contributed reagents/materials/analysis tools: DS SZ ZS YW. Wrote the paper: JL JQ ZZ SH.

References

  1. 1. Robinson RW, Decker-Walters D (1997) Cucurbits: Cab international.
  2. 2. Tatlioglu T (1993) Cucumber Cucumis sativus L. Genetic improvement of vegetable crops. Tarrytown: Pergamon Press Ltd 197–234.
  3. 3. Ren Y, Zhang Z, Liu J, Staub JE, Han Y, et al. (2009) An integrated genetic and cytogenetic map of the cucumber genome. PLoS One 4: e5795.
  4. 4. Han Y, Zhang Z, Liu C, Liu J, Huang S, et al. (2009) Centromere repositioning in cucurbit species: Implication of the genomic impact from centromere activation and inactivation. Proc Natl Acad Sci USA 106: 14937.
  5. 5. Huang S, Li R, Zhang Z, Li L, Gu X, et al. (2009) The genome of the cucumber, Cucumis sativus L. Nat Genet 41: 1275–1281.
  6. 6. Guo S, Zheng Y, Joung JG, Liu S, Zhang Z, et al. (2010) Transcriptome sequencing and comparative analysis of cucumber flowers with different sex types. BMC Genomics 11: 384.
  7. 7. Durbin RM, Altshuler DL, Abecasis GR, Bentley DR, Chakravarti A, et al. (2010) A map of human genome variation from population-scale sequencing. Nature 467: 1061–1073.
  8. 8. Weigel D, Mott R (2009) The 1001 genomes project for Arabidopsis thaliana. Genome Biol 10: 107.
  9. 9. Huang X, Wei X, Sang T, Zhao Q, Feng Q, et al. (2010) Genome-wide association studies of 14 agronomic traits in rice landraces. Nat Genet 42: 961–967.
  10. 10. Sebastian P, Schaefer H, Telford IRH, Renner SS (2010) Cucumber (Cucumis sativus) and melon (C. melo) have numerous wild relatives in Asia and Australia, and the sister species of melon is from Australia. Proc Natl Acad Sci USA 107: 14269–14273.
  11. 11. Li Shizhen Compendium of Materia Medica. Volume 28: Chapter vegetables, part III: cucurbits 14–15.
  12. 12. Meglic V, Serquen F, Staub JE (1996) Genetic diversity in cucumber (Cucumis sativus L.): I. A reevaluation of the US germplasm collection. Genet Resour Crop Evol 43: 533–546.
  13. 13. Meglic V, Staub JE (1996) Genetic diversity in cucumber (Cucumis sativus L.): II. An evaluation of selected cultivars released between 1846 and 1978. Genet Resour Crop Evol 43: 547–558.
  14. 14. Staub JE, Serquen FC, Mccreight JD (1997) Genetic diversity in cucumber (Cucumis sativus L.): III. An evaluation of Indian germplasm. Genet Resour Crop Evol 44: 315–326.
  15. 15. Staub JE, Serquen FC, Horejsi T, Chen J (1999) Genetic diversity in cucumber (Cucumis sativus L.): IV. An evaluation of Chinese germplasm1. Genet Resour Crop Evol 46: 297–310.
  16. 16. Stewart CN Jr, Via LE (1993) A rapid CTAB DNA isolation technique useful for RAPD fingerprinting and other PCR applications. BioTechniques 14: 748–749.
  17. 17. Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155: 945–959.
  18. 18. Evanno G, Regnaut S, Goudet J (2005) Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol 14: 2611–2620.
  19. 19. Liu K, Muse SV (2005) PowerMarker: an integrated analysis environment for genetic marker analysis. Bioinformatics 21: 2128–2129.
  20. 20. Cavalli-Sforza LL, Edwards AWF (1967) Phylogenetic analysis models and estimation procedures. Am J Hum Genet 19: 233–257.
  21. 21. Tamura K, Dudley J, Nei M, Kumar S (2007) MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Bio Evol 24: 1596–1599.
  22. 22. Excoffier L, Laval G, Schneider S (2005) Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evol Bioinform online 1: 47–50.
  23. 23. Gouesnard B, Bataillon T, Decoux G, Rozale C, Schoen D, et al. (2001) MSTRAT: An algorithm for building germplasm core collections by maximizing allelic or phenotypic richness. J Hered 92: 93–94.
  24. 24. Ranc N, Muños S, Santoni S, Causse M (2008) A clarified position for solanum lycopersicum var. cerasiforme in the evolutionary history of tomatoes (solanaceae). BMC plant biology 8: 130.
  25. 25. Zhou H, Xie Z, Ge S (2003) Microsatellite analysis of genetic diversity and population genetic structure of a wild rice (Oryza rufipogon Griff.) in China. Theor Appl Genet 107: 332–339.
  26. 26. Liu K, Goodman M, Muse S, Smith JS, Buckler E, et al. (2003) Genetic structure and diversity among maize inbred lines as inferred from DNA microsatellites. Genetics 165: 2117–2128.
  27. 27. Slatkin M (1985) Rare alleles as indicators of gene flow. Evolution 39: 53–65.
  28. 28. Dabholkar AR (2006) General Plant Breeding. New Delhi: Concept Publishing Company. 318 p.