- Research article
- Open Access
Landraces of snake melon, an ancient Middle Eastern crop, reveal extensive morphological and DNA diversity for potential genetic improvement
© The Author(s). 2018
- Received: 2 December 2017
- Accepted: 30 April 2018
- Published: 23 May 2018
The Correction to this article has been published in BMC Genetics 2018 19:44
Snake melon (Cucumis melo var. flexuosus, “Faqqous”) is a traditional and ancient vegetable in the Mediterranean area. A collection of landraces from 42 grower fields in Israel and Palestinian territories was grown and characterized in a “Common Garden” rain-fed experiment, at the morphological-horticultural and molecular level using seq-DArT markers.
The different landraces (“populations”) showed extensive variation in morphology and quantitative traits such as yield and femaleness, and clustered into four horticultural varieties. Yield was assessed by five harvests along the season, with middle harvests producing the highest yields. Yield correlated with early vigor, and with femaleness, but not with late vigor. At the molecular level, 2784 SNP were produced and > 90% were mapped to the melon genome. Populations were very polymorphic (46–72% of the markers biallelic in a 4 individuals sample), and observed heterozygosity was higher than the expected, suggesting gene flow among populations and extensive cross pollination among individuals in the field. Genetic distances between landraces were significantly correlated with the geographical distance between collecting sites, and with long term March precipitation average; variation in yield correlated with April temperature maxima.
The extensive variation suggests that selection of local snake melon could result in yield improvement. Correlations between traits and climatic variables could suggest local adaptation of landraces to the diverse environment in which they evolved. This study stresses the importance of preserving this germplasm, and its potential for breeding better snake melons as an heirloom crop in our region.
- Cucumis melo
- Snake melon
- Population structure
- Genetic diversity
Local landraces of crop plants represent an invaluable genetic resource for breeding in a changing environment. Landraces exhibit fine adaptation to the specific environment in which they have evolved under domestication. They often harbor rich genetic diversity, important for traditional agriculture sustainability, that has been eroded in modern cultivars [1, 2]. It is therefore important to generate germplasm collections, describe landraces thoroughly and evaluate their diversity.
Israel is located on the border of three different phyto-geographical zones: the Mediterranean, Saharo-Arabian, and Irano-Turanian zones, respectively . Over very short distances there are notable differences in rainfall, temperature regimes and soil types. In addition, yearly fluctuations in climatic parameters are large, and such instability has been proposed to favor genetic diversity, compared to very stable environments [4, 5]. In addition to natural factors, domestic landraces are also affected by active human agricultural practices, that directly shape plant phenology and anatomy, and exert selective forces that in the long run will determine heritable differences between populations [6–9]. These can be uncovered by growing populations from diverse environments in a “common garden” experiment (e.g., ). Being part of the Fertile Crescent, our region is the cradle where many crops such as durum wheat, barley, pea, lentil have been domesticated . They continue to evolve under domestication to present days, and many wild relatives being part of a crop’s gene pool still grow in nature.
Cucumis melo L. is perhaps the most genetically diverse crop among the Cucurbitaceae. In the Western world, sweet dessert melons are best known. They include Cucumis melo var. cantalupensis, having climacteric, aromatic fruits, and C. melo var. inodorus, with non climacteric, long keeping fruits, but the species also includes wild accessions (C. melo var. agrestis) and a few non-sweet cultivated varieties ( and refs. therein). Among these, Cucumis melo var. flexuosus, the snake melon, probably represents an early stage in melon domestication, where fruits have become non-bitter, but not yet sweet, offering a refreshing vegetable that can be eaten fresh, pickled or cooked at the immature fruit stage, like squash and cucumber ( and refs therein; [14, 15]). Later towards maturation they become inedible. Plants of snake melon are monoecious, ovaries are pubescent, and the fruits, as well as the stems, petioles and flower pedicel are elongated. The Arabic name for the local snake melon, Faqqous, originated from the Hebrew term “Piquous” (of Greek origin), indicating removal of fruit soft hair before consumption. Faqqous was known to ancient Hebrews as Qishut, mentioned in the bible as a vegetable eaten in ancient Egypt, and in the Talmud as a common cucurbit vegetable, whereas the cucumber (Cucumis sativus) arrived to the region much later, and Today’s squash is of course an acquisition from the New World ( and refs therein). This implies that Faqqous had a very long, probably uninterrupted history of cultivation in our region, rendering it an interesting germplasm resource, whose genetic composition could reflect its ancient and local origins.
At present, snake melon represents a traditional vegetable crop, mainly grown by Israeli Arab farmers in the Lower Galilee and in the Palestinian villages of Judea, Samaria and Jordan Valley under rain-fed conditions. In the framework of an Israeli-Palestinian collaboration, a collection of Faqqous accession was assembled . In the present study, samples of the collection were sown in a common plot in Israel, to attain a full description of their phonological and horticultural diversity. Such analysis could suggest factors and traits that could limit economic productivity and could be improved by future breeding. We also probed the genotypic diversity of the collection and assessed population structure, using the seq-DArT method [17, 18]. The data collected were correlated with climatic and geographical parameters, providing insight on the possible forces that have shaped snake melon landraces in this small but climatically diverse region.
Set up of “common garden” experiment
To set up the collection, a single open-pollinated fruit from approximately 10 individual plants of a given agricultural field was sampled and its seeds represented an accession. Accessions from the same field were evenly mixed and represented a population.
Morphological, phenological and horticultural traits were scored individually for each plant.
Early Vigor (coverage) of plants at the vegetative phase was evaluated 48 days after sowing. Each plant was scored twice by independent persons. We performed a qualitative visual score, where 1 represented the smallest plants and 5 – vigorous, largest plant (Fig. 1b). 2). We also counted all branches that exceeded 20 cm length and recorded the length of longest branch, as additional measures of early vegetative vigor. Since these measures were positively correlated we further analyzed only the visual scores.
Female flowering. We scored the earliest presence of open female flowers at three dates, starting 54 d post sowing, to attain good estimation of female-earliness. We also counted the female flowers aged from 2 days before anthesis to two days post anthesis in a continuous 11 day window in the beginning of female flowering, as a quantitative estimate of the total number of female flowers that develop on each plant.
Leaf morphology. Leaves were sorted in 5 different shapes (Fig. 1c), distinguishable by the lobes of varying depth.
Yield measurement and fruit traits. Immature fruit of snake melon are harvested continuously, every 2–3 days. We performed five “sample harvests” to provide yield estimate. Between these “snapshot harvests” the farmer made additional commercial harvests that were not recorded. Fruit longer than 10 cm were picked, and their number and total weight were recorded for each plant.
Fruit morphology and quality. In the two of the experimental harvests, the fruit from each plant were described. Two fruit from each plant (the largest and the smallest) were analyzed. Morphological traits were recorded as follows:
Color. 1 - light green, 2 - medium green, 3 - dark green, 4 - marbled (finely dotted over a gray background). Dots: present or absent. Straightness: 1 - straight, 2 - bent, 3 - winding. Fruit shape was determined as the length to diameter ratio that remains rather constant along fruit development. Diameter was measured in the middle of the fruit. Ribs/furrowing of the rind: 1 - deep, Sahouri like (few well demarcated, parallel furrows); 2 - like 1, but shallow; 3 - Baladi like (dense, less ordered wrinkles), deep furrows; 4 - like 3, but shallow (Figs. 1a, and 3).
For taste and texture evaluation, at each harvest two fruit each plant were assayed by a team of three participants. Bitterness: present/absent. Taste: 1 - tasty, 3 - not tasty. Texture: 1 - hard and crispy, 3 - soft. Hollow fruit (space between pericarp and placenta): present/absent (Fig. 1d).
Late vigor. Later in the season (73 d - 94 d after sowing) we scored six times the general vigor of the plant, until the field finally collapsed in July. We used a visual scale from 0 – dead plant to 5 – green and vigorous, the area under the plot providing an integrated, inverse measure of the senescence process.
Genotyping by diversity Array technology (seqDArT)
Young leaves or side shoot tips were collected for DNA extraction from individual plants of two of the replicates, refrigerated in the field and stored frozen in the deep freezer. A total of 94 DNA samples from 23 populations × 4 individuals and two outgroup samples were prepared by the CTAB method , and sent for seqDArT genotyping at DArT, Canberra http://www.diversityarrays.com/dart). Genotyping protocols were as detailed by Al-Beyroutiova et al. , except for replacing HpaII with MseI digestion. Amplicons containing SNP were located on the melon genome . Genotyping data was analyzed using the GenAlEx 6.502 software .
Determination of soil type, soil properties and climatic data
Soil samples were collected from each collection site and visually inspected. Soil type was determined according to the detailed soils map of Israel based soil surveys by the Agricultural Research Organization and Soil Conservation and Drainage Department (http://www.govmap.gov.il; . Calcium carbonate content of the soil was determined according to Rowell  using a devise modified by M. Sharabani, Geography and Environment Dept., Bar Ilan University. The % hygroscopic water content was determined by weighing the air-dried soil (Sw), then oven-drying for 24 h at 105 °C (Sd), calculated as 100 × (Sw-Sd)/Sd).
Climatic data based on yearly averages between 1960 and 1990 were taken from http://www.bioclime.org. Annual evaporation data was received from Dr. Noam Halfon, Israel Service of Meteorology.
The following traits were statistically analyzed as specified below: taste, texture and proportion of hollow fruits; early, late and total yield; early and late vigor; number of female flowers and days to first female flower; fruit taste, texture and hollow fruit. In order to select the proper statistical tests for each trait, we checked the Normality of distribution of the traits using the Kolmogorov-Smirnov test on the ~ 564 individual plants in the experimental field. Traits exhibiting significant deviation from normality were transformed using logarithmic and square-root transformation, respectively, and tested for normality. In case transformation resulted in normal distribution, we used the transformed data for further testing.
For traits with normal distribution, parametric tests were used: Pearson correlation, Tuckey test, 1-way ANOVA followed by Tuckey post hoc test for differences among groups). For traits with non-Gaussian distribution, non-parametric tests were used: Spearman correlation, Kruskal-Wallis rank sum test for differences among groups, followed by Mann Whitney pairwise comparisons, corrected for multiple comparisons. The tests were run using R (The R core development team; https://www.r-project.org/) and SPSS (IBM SPSS Statistics - Essentials for Python 23.0, SAS Institute, Cary NC).
All the above ten traits had a significant deviation from normality (not shown), and logarithmic transformation did not improve the distribution. Upon square-root transformation (and block-correction, see below), Early Yield, Late Yield, Total Yield and Female-Flower Number attained normal distribution (p = 0.83, 0.12, 0.43, 0.34, respectively) in the Kolmogorov-Smirnov test, while the distribution of the remaining five traits remained non-normal (p < 0.05).
Correction of block effect
We checked whether the three replicates (blocks), located along a moderate slope, could differ in plant performance. We used 1-way ANOVA for early and total yield (after sqrt transformation), and Kruskal-Wallis test for the other five traits. Six traits showed a significant block-effect (p < 0.05), and the seventh trait, late vigor (“senescence”), was near-significant (p = 0.08). To compensate for such effect we performed a correction, adding/subtracting the difference of the block average from the whole field average; the remainder of the analysis was carried out using block-corrected values.
Dendrograms of accessions
A dendrogram based on morphological traits was plotted using the Hierarchical cluster analysis, “between group linkage” method. Distance values were computed by dissimilarity squared Euclidian distances, rescaled between 0 and 1 (SPSS, IBM SPSS Statistics - Essentials for Python 23.0, SAS Institute, Cary NC). A Neighbor Joining tree was inferred from the populations SNP allele frequencies using the DIST.GENPOP function in ADEGENET R package, to calculate Nei genetic distances and NJ function in APE R package.
Correlation between climatic data and plant data
The Random forest algorithm implemented in the Boruta R package was used to explore the most important variables of climatic data that are correlated with quantitative morphological traits, coefficients of variation (CV) of the traits, and molecular diversity. All climatic variables were used in one regression model. Spearman test for non-parametric correlation was used to estimate the goodness of fit (r) and the direction of the correlation.
Evaluation of a snake melon collection in a common garden set up
Twelve plants of each population, (3 replicates × 4 plants, ca. 560 plants in total), were sown in a common field in Sandala, Yezre’el Valley, following the traditional, rain-fed agricultural practice used for snake melon, as described in the Methods. Phenotypic analysis, conducted by scoring each individual plant, was divided into morphological traits, mostly scored as discrete states, fruit quality parameters of interest to consumers (measured semi-quantitatively), and horticultural traits measured quantitatively.
Variation in morphological traits
The morphological-qualitative traits were as follows: fruit shape (L/W ratio), fruit color, presence of dots, rind furrows, fruit straightness, leaf shape. Scoring was done as detailed in the Methods section.
Additional file 1: Table S1 provides the scoring data regarding these traits. For each trait, it provides the segregation of the trait-states among individual plants.
Fruit dimensions constantly change due to rapid growth, but the length:width (l:w) ratios are rather constant from anthesis to ripening [25, 26], and provide a solid morphological character. When we plotted the l:w ratios from smaller and larger fruits in the same individual plant, we obtained a very high correlation (R2 = 0.97). The ratio (in addition to the furrowing pattern) readily distinguishes Sahouri (4.1 average ratio) from Baladi (5.3 ratio), while some exotic snake melons such as FLI have very elongated (15 ratio), snake-like fruits that are not found in our region. When we classified leaves by their shape (Fig. 1c), we noticed that the prominent leaf types were 3 and 4, with all Sahouri populations displaying state 4 and Baladi- state 3. FLI had Type 5 leaves, whereas leaf shapes 1 and 2 that are common in sweet melons, were not found consistently in our collection. Dotted fruits were quite rare, and bending was more common in Baladi than in Sahouri.
Fruit quality evaluation
Four fruit quality traits, namely the presence of bitterness, general taste (of non-bitter fruits), presence of a hollow fruit and fruit texture were evaluated and scored as described in the Methods section. Snake melon is eaten immature about one week post anthesis, and mature fruits become inedible. We tasted, in replicated trials, smaller and larger fruits of each plant (at 6–7 dpa and 8–9 dpa, respectively), the younger ones are preferred by consumers while the larger stage is still acceptable.
Bitterness was uncommon and few bitter fruits were found in 11 of the populations (mostly a single plant/fruit out of ~ 12), representing all four horticultural varieties. Populations AB60 from Bardala (N. Jordan Valley), QG23 (Qalqilia area) BB56 (Bet Lehem) and HS43 (Hebron area) stood out, with ~ 25% of the plants bearing bitter fruits.
Separation of the mesocarp from the placenta (Fig. 1d) is frequent in mature fruit, but is less desirable in immature fruit at the commercial stage. The proportion of hollow fruit among the smaller fruit of the four varieties ranged between 9 and 17%, and quickly increased to 44–64% in the larger fruit with higher ratios in SG landraces (Fig. 5c). Among the individual populations, BB55 (SG), HH49 (BW), JZ03 (BW) had the smallest proportions of hollow fruit.
Quantitative phenological/horticultural traits
The most important trait for future improvement of Faqqous is yield, a complex trait that is probably affected by numerous environmental and plant developmental factors. We compared the average values and extent of variation within the 47 “populations” and four horticultural varieties, for seven selected traits related to plant performance in general and fruit yield in particular. These included early plant vigor, the number of days until the appearance of the first female flower, number of female flowers in a determined time window, early, late and total fruit yield, and late vigor, measured towards the end of the season (an inverse measure of senescence).
Male flowers reached anthesis six weeks after germination, and one week later the female flowers, that are fewer and borne singly on side branches, first appeared. Plants that exhibited fast vegetative growth have the potential to attain female flowering earlier, and those having a stronger or longer female phase could have a larger yield. The number of days from sowing to first female flower ranged between 54 and 61; among the four varieties BG was earliest and SW the latest (Fig. 7b). However, variation within populations was large and individual plants could differ by 4 days. When we compare the number of female flowers that opened in an 11 day window, large variation was apparent, ranging between two flowers in SW population HS43 from Suba, Hebron and 16 flowers in BG population TB16 from Tell (Fig. 7c). Among varieties, BG displayed more profuse flowering and SW had the fewest female flowers. Because the time window selected was at the onset of flowering, an inverse correlation between the two femaleness parameters was expected; measuring an additional time window later in season would have been desirable but would interfere with fruit picking and yield measurement.
Picking of marketable fruit (10–30 cm long) was done every 3 days, starting after 55 days, and lasted for ~ 40 days. Because of the indeterminate growth habit and continuous requirement to harvest immature fruit (non-harvested fruit would inhibit further fruit set), it was not technically possible to measure the actual yield of each plant. Instead, we performed five “snapshots” or sample-harvests, in which fruits of individual plants were picked, counted and weighed, and, in-between, the farmer performed additional harvests and “cleaned” the plants. Total Yield is the sum of all five experimental harvests, and the actual yield along the season can be estimated by multiplying the total experimental yield by 2. The average experimental yield in our experiment was 1.4 kg/plant (representing a full yield estimate of 2.8 kg) and ranged between 0.8 and 2.3 kg/plants. Figure 7d presents the total yield (summed over the five experimental harvests) of 47 populations and four varieties. The highest yielding populations were SD30 from Dear Baloot and MNZ from Kafr Manda, and the lowest yielding were FLT from Italy and JM05 from Meslyeh. Among the varieties, BW ranked significantly higher and BG – lower.
Later in the season, we scored plant vigor using a qualitative scale (see Methods section), and integrated observations taken on 6 different dates (6 June - 5 July) into an “area under survival curve” value, where higher values represent higher late-vigor. Fig. 7e plots such values for each population and for the four Faqqous varieties. Late vigor scores varied from 27 to 54 units, Green Sahouri senesced earlier than the other varieties, while White Balladi exhibited, on average, prolonged vegetation.
Correlations among traits
Correlation (or lack of correlation) between quantitative traits could provide useful insight on factors that determine such traits, and generate ideas regarding future genetic improvement, based on the large phenotypic variation apparent in our collection. Spearman correlations were computed between all pairs of seven traits, using the individual scores from ca. 500 plants.
Early vigor was strongly correlated (0.65, p < 10-4) with the number of female flowers that was measured, as explained above, in an early time window, and with the days to first female flower (− 0.67, p < 10-4). Female flower number was correlated with early yield (0.37, p < 10-4; Fig. 8d). The correlations of early vigor with early and total yield had smaller, but statistically significant, correlation coefficients (0.28, p = 0.0004 and 0.27, p = 0.001, respectively). This suggests that while early vigor and establishment contribute to yield, plants can later catch up and compensate for this parameter. On the other hand, we saw no correlation between late vigor, i.e., whether a plant remained green longer, and its total yield, or even the late yield (Fig. 8e). It could be that the late-senescing plant produced a low yield in the earlier (or even later) phases, such that late senescence did not correlate with higher productivity. The lack of negative correlations could suggest that variation in both traits (early vigor and late vigor, and possibly prolonged productivity) could be combined by breeding.
Genetic diversity revealed by seqDArT
Seq-DArT (Diversity Array Technology , is a high-throughput genotyping method that generates sequence-based markers. These include co-dominant single nucleotide polymorphism (SNP), as well as dominant markers representing presence/absence of amplicons (“in silico DArTs”). In order to obtain a broad picture of the genetic structure in our snake melon collection, we selected a subset of 23 populations that represent the four horticultural varieties and the different geographical regions, and sampled four individual plants/population, in addition to two individuals of the out-group varieties, FLT (a Faqqous inbred genotype from Italy) and Vedrantais (a sweet melon inbred cultivar).
Distribution of DArT markers on the melon genome. Chromosome length in cM is according to Perin et al. 
Map Length (cM)
Based on DArT genotyping, we asked how much molecular variation is present in the collection, and how is it partitioned within individuals (in the form of heterozygosity), between individuals in a population, and among populations. Is heterozygosity higher or lower than expected, and what can be learned from this on population structure and reproductive strategy/ history of Fakus in our region? Does the molecular data divide the collection according to population, horticultural variety, geographic region?
Statistics of genetic variation 23 melon populations as measured by 2748 high quality SNP loci, GenAlex software
The Fixation index (“F”; Table 3) is calculated as 1-Ho/He. Most populations had a negative index, meaning that heterozygosity is higher than expected by Hardy-Weinberg equilibrium. Possible reasons for such situation are discussed below.
Statistically significant correlations between population means, coefficient of variation (CV) of quantitative traits, or genetic variation, and long-term climatic variables, determined by random forest feature selection and Spearman correlation
Early TotalYield CV
Early TotalYield Mean
Female Fl. Date Mean
Late Vigor Mean
The coefficients of variation of both the total yield and the early yield (two highly correlated traits) were negatively and strongly correlated (rho = − 0.5) with the maximum temperature in April. This parameter ranged between 30 °C in the Jordan valley and 22–23 °C in most regions. This could imply that extreme environmental conditions tend to reduce genetic variation, as shown by Peleg et al. .
More correlations were found between plant traits and climatic variables, but not soil or geographic parameters (Table 4). Early vigor was negatively correlated with the mean temperature in March, suggesting that cool springs could have selected quicker establishment of the plants. Early yield, on the other hand, positively correlated with evaporation in March. Early-season evaporation could reflect windy and/ or warm weather and this could favor early and higher yielding and perhaps a shorter life cycle. Late vigor, i.e., maintaining green foliage late in season, was very strongly correlated (Rho = 0.72) with minimum temperature in June. Such minima range between 14 °C in the Hebron area and 22 °C in the Jordan and Northern valleys. In the former areas, characterized by continental climate, plants could have evolved a shorter life cycle and senesced faster than in areas with warmer nights.
Finally, molecular variation, expressed as unbiased expected heterozygosity, was positively correlated with March precipitation, which range between 23 mm in the dry Jordan valley and 100 mm in Ramalla. In our Mediterranean, rain-limited area, less extreme conditions apparently favored genetic variation.
Snake melon (Cucumis melo var. flexuosus) is an ancient crop in Israel and Palestinian territories as well as other Mediterranean countries. In a previous study, a collection of landraces of Faqqous, the Palestinian snake melon, was assembled and a preliminary description provided based on observations in situ, at collection sites in Israel and the West Bank . Faqqous is an important summer vegetable grown without irrigation in Arab villages, and has not been subjected to modern breeding practices [14, 28]. During the collection trips we learned that each farmer selects a few open-pollinated fruits to save seeds for the following season. We hypothesized that such germplasm harbors significant genetic variation generated by long-term evolution under domestication, and by exposure to the large variation in temperatures and precipitations found in this small but very diverse region. We also hypothesized that such variation could be further selected by future breeding programs, to increase Faqqous yields and tolerance to biotic and a-biotic stress.
Here we report the characterization of 42 landraces from Judea, Samaria, the Lower Galilee and Jordan Valley in a common garden experiment conducted in Sandala, Yezre’el Valley. Our purpose was to describe the collection in a comparative manner, assess its phenotypic variation and agronomic potential and derive genetic and agronomic insights for future conservation and breeding. Such an experiment complements in situ observations, where each variety is assessed in the site to which it has adapted; the common garden provides a uniform environment (sowing time, soil-water content, simultaneous assessment of parameters), that is essential to characterize the collection and select items from it. The same collection will need to be analyzed in additional environments to provide a full picture of genetic x environment effects on important traits.
Cucumis melo is a highly diverse species, traditionally divided in two subspecies: ssp. agrestis, with wild accessions having inedible, often bitter fruits, and ssp. melo, further divided to six varieties: cantalupensis, inodorus, conomon, chito, dudaim and flexuosus [12, 29]. The first two varieties consist of sweet dessert melons (climacteric and long-keeping, respectively), popular in the western world since the Roman period. The others are mostly non-sweet varieties, but in the Far East, parallel domestication gave rise to sweet cultivars of the Makuwa group. Stepansky et al.  used molecular markers on a 60-accessions core-collection that broadly supported such intraspecific classification, and placed snake melons closer to the sweet melon, in concordance with their Middle-Eastern origin.
It has been thought that Cucumis melo speciated in Africa, but this view has been recalled by Sebastian et al. , based on analysis of a diversified array of Asian and Australian species, from which both melon (2n = 24) and cucumber (2n = 14) must have originated. The first domestication site of melon is unclear, and multiple domestication events are likely. It has been suggested that melons were domesticated in two stages, first as a non-sweet vegetable whose fruit is eaten shortly after anthesis, and only later as a sweet fruit eaten mature. Faqqous and a few other melon varieties (e.g., conomon) represent the first stage of domestication that required removal/mutation of bitterness genes. Cucurbit fruits produce the triterpenoid cucurbitacin, and interestingly such trait persists sporadically in the domesticated germplasm and needs to be selected against in crops such as cucumber and watermelon. Bitterness is affected by the environment, and bitter progeny may arise from crossing two non-bitter parents. In cucumber, it is controlled by several interacting genes [31, 32] that determine whether the foliage and/or the fruits are bitter. In our collection, 23% of the populations included bitter-fruited individuals, despite the obvious selection that had to be exerted against such deterrent trait. This is an obvious target for future improvement. The reasons for persistence of bitterness are presently unknown, perhaps poor heritability, or lack of awareness by growers who select the seeds. Growers’ folklore maintains that bitterness comes from pollen of nearby Echballium elaterium, a distantly related cucurbit weed.
A more realistic concern about Faqqous inter-crossing relates to the novel trend of consuming, as “Faqqous surrogate”, fruits of a sweet melon variety at the immature stage. As a result, Faqqous landraces are grown side by side with sweet melons and can cross via honeybees, resulting in genetic contamination of the traditional landraces. This further highlights the need for collecting and preserving the Faqqous germplasm.
Studies on melon diversity are numerous, but those focusing on snake melons are very few. Indian snake melons were described by Pandey et al. ; a Jordanian collection of 8 accessions was described by Abdel-Ghani & Mahadeen , who analyzed morphological traits and random amplified polymorphic DNA (RAPD) markers. Our study, with 46 populations × 12 plants/population provides a more in-depth evaluation of snake melons from a given area, with good genome sampling by a modern genotyping method.
When examining morphological traits in a qualitative manner, landraces from Israel/Palestinian territory fall into four main horticultural varieties: slender, wrinkled Baladi versus stout, straight-furrowed Sahouri, each of these including dark and light-colored accessions. Some populations appeared more variable that others and contained atypical fruits, indicating gene-mixing between populations. When examining fruit-quality traits, i.e., occurrence of bitterness, crispy texture, good flavor and hollow fruit, we observed significant variation within populations that could result from genetic variation, or from low heritability and difficulty in the objective assessment of such traits. Nevertheless, a few populations stood out as being better or worse than the majority (Fig. 5), and these should be rechecked in pre-breeding screens for stability of the trait. The correlation between good taste and crispiness or hollowness is not simple (Fig. 6), but hollow and softer fruits were often graded less tasty.
Yield improvement is an important breeding goal for 'Faqqous’ and progress towards such goal starts by evaluating diversity in yield, and in traits that could relate to yield. The randomized block experiment in a single environment allowed careful scoring of vegetative and reproductive growth of each individual plant. We detected significant variation in total yield (composed of early and late yields), in the early and late vigor representing, respectively, the establishment phase in young plants and the senescing phase at the end of the season. We also measured the degree of femaleness using two different parameters. Due to the high intra-population variability (later corroborated by our molecular data), differences among populations were statistically significant only between the populations from the two extremes of the spectrum (Fig. 7). Repeating such experiments with larger samples and fewer populations should confirm the heritability of such differences. In the present experiment, however, the large variation and statistical significance between extremes support the likely existence of useful genetic variation in the collection, both within and between populations.
We hypothesized that populations that can quickly establish a vegetative biomass and/or maintain a productive, green biomass later in the season (excelling in early and/or late vigor) would have higher yields. We also expected to find positive correlations between the femaleness parameters and yield, since in traditional varieties of cucurbit crops that have not been subject to modern breeding, female flowers are limiting, and represent an important yield component .
We plotted the Spearman correlations between quantitative traits, based on > 500 individual plant measurements (Fig. 8). Such plots indicated statistically significant correlations between early vigor, femaleness and yield. There was also large variation regarding performance late in the season, when some accessions stayed green for an extended period, when others had already collapsed. However, this potentially important trait was not statistically correlated with yield. This could suggest that different accessions employ different pathways/ adaptation strategies that could involve earlier or delayed fruiting, or even different ways to utilize soil water and resist stress. The large variation and lack of correlation between such traits actually indicate that exceptional genotypes are present in population that combine, for example, early vigor and late vigor, or high femaleness and high yield, and there is room for genetic amelioration by selecting superior trait-combinations.
Genotyping germplasm collections with markers can eliminate duplications and misnaming. More importantly, the distribution of diversity among accessions and collection sites can help understand population structure, and direct further collection efforts and germplasm maintenance (e.g., [35–37]. Genotyping core collections using molecular markers has been facilitated by novel SNP discovery techniques, based on next generation sequencing methods (NGS; [38, 39]. These provide thorough genotyping with thousands of SNP per individual genotype, allowing association studies that locate important genes; as a result, genotyped core collections become a permanent, invaluable resource not only for conservation but also for genetic research. In this study, we used the seqDArT genotyping method that has been developed from a microarray-based platform into an NGS-based one . We generated, in a cost-effective manner, 2460 best-quality SNPs, without missing points in the entire 94 sample-set. Over 90% of the SNP could be aligned, and were evenly distributed, on the melon genome; a larger number of dominant markers was derived as well.
Of the 2460 polymorphic markers, a very high proportion, 46–72%, were polymorphic within each population, considering the small sampling size (n = 4). AMOVA analysis and unbiased heterozygosity calculations show that heterozygosity is extremely high, 87% of the molecular variance being attributed to individual’s heterozygosity. This leaves small proportions of variation within and between populations (Fig. 10), and translates in higher-than-theoretically expected heterozygosity, and negative Fixation Indices, in 18 of 23 populations (Table 3). This is rather surprising, in view of the repetitive selection by the farmers who collect a very small number of fruits per field (~ 10, each with ~ 300 seeds, from a field of 1000 plants, to produce the 3000 seeds he may need for the next season). This could be expected to progressively reduce heterozygosity to low levels. Thus, high heterozygosity and negative F index could result from positive selection for heterozygous loci, or from very low levels of self-pollination of individual plants. The latter could be the likely result of the monoecious sex habit of Faqqous (separate male and female flowers on same plant), and the very small number of female flowers found at a given time point on a single plant, favoring cross pollination between plants by honey bees and other pollinators. These pollinators may transfer pollen from adjacent fields, that will increase the diversity within individuals and decrease the diversity between populations, as was observed in the AMOVA. A Neighbor-Joining tree based on the same data also reflected the high within-population variation leading to less demarcated separation between populations (Fig. 10b). Still, the best tree obtained faithfully separated the four horticultural varieties, even better than the morphological tree (Fig. 4), that was based on few traits, compared to the extensive genomic sampling of the molecular tree. This confirms that the four main horticultural varieties have been selected over a rather long period and reflect true genetic divergence; AWG, a commercial Green Baladi Faqous recently bred in Natufa valley, clusters with White Baladi from the same area, from which it could have been derived. Also in the Sahouri group, the green accessions seem to stem out of the larger White Sahouri cluster. In the morphological tree, fruit color greatly influences tree structure, being a clear-cut trait, and the two green varieties cluster together against the two white ones (with few exceptions).
In a study of Coffea genetic diversity using DArT seq, Observed Heterozygosity values were much lower (0.124; Garavito et al.  compared to our average of 0.364, and F values were positive. In Cicer, He and uHe values (based on DArT markers) were low in the cultivated accessions (~ 0.15), and reached ~ 0.3 in wild accessions .
Populations are seldom reproductively isolated, and gene flow either by pollen or seed exchange can provide mixing. On the other hand, the positive correlation between geographic distance and genetic variation (Fig. 11) supports local evolution of landraces, and their adaptation over a long period of local selection, to the diverse environments found in Israel over short distances, e.g., 170–600 mm annual precipitation. PCoA analysis divided the samples in a way that largely reflects the geographic map, including the separation between Israel and the West Bank territories that have been administratively separated for the last decades. Moreover, the fact that a few climatic parameters were correlated with genetic variation, with specific quantitative traits, or with their coefficients of variation (Table 4), also supports a long-term local evolution of landraces in response to climatic factors. The high correlation of genetic distance with geographic distance suggests that gene flow is mediated by local agents such as pollinators and not by long distance agents such as human trade.
Landraces in the surveyed area fall into four morphologically distinct horticultural varieties that are also genetically distinct. The extensive phenological-horticultural variation found in this collection, both within and between populations, suggests that selection for yield and related traits, such as femaleness and vigor, could result in yield improvement. Populations taken from individual grower field are highly heterozygous and polymorphic, in agreement with high rates of out-crossing between monoecious plants and significant gene flow between sites. Observed correlations between quantitative traits and climatic variables could suggest local adaptation of landraces to the diverse environment in which they evolved, and support our view of snake melon as an ancient link in melon domestication in the Near East. Our data stress the importance of preserving such germplasm, and demonstrate its potential for breeding better snake melons as an heirloom crop in our region.
This research was funded by USAID-MERC (US Agency for International Development - Middle East Regional Cooperation) project TA-MOU-12-M32–016. The funding served for performing the collection, generating and analyzing experimental data.
Availability of data and materials
Full phenotypic and genotypic data files that served to generate this study are deposited in the Israeli Plant Gene Bank database http://igb.agri.gov.il/web/ by clicking “view information” and “attached files” for items number 300031–300089 in the Bank catalogue. Seed samples for research purpose will be available from the Israeli Gene Bank and from the corresponding author.
SO performed most of the field and laboratory work and statistical analyses in this study. YK and JB supervised and performed statistical analyses. SP directed the soil analyses. HS planned and supervised population genetics and climate × variation analyses. RPT was the principal investigator in the project and led the writing of this manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Not applicable. Study has not involved humans, animals, wild plants, endangered plants, or engineered plants. Work with cultivated plants involved standard agronomic practices requiring no special permissions. Seeds for this study were donated by local farmers for research purpose, requiring no special permissions. Seed samples were deposited in the Israeli Plant Gene Bank as detailed below.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Clawson DL. Harvest security and intraspecific diversity in traditional tropical agriculture. Econ Bot. 1984;39:56–67.View ArticleGoogle Scholar
- Pasam RK, Sharma R, Walther A, Ouml Zkan H, Graner A, Kilian B. Genetic diversity and population sructure in a legacy collection of spring barley landraces adapted to a wide range of climates. PLoS One. 2014;9:e116164.View ArticlePubMedPubMed CentralGoogle Scholar
- Danin A, Plitmann U. Revision of the plant geographical territories of Israel and Sinai. Plant Syst Evol. 1987;156:43–53.View ArticleGoogle Scholar
- Westberg E, Ohali S, Shevelevich A, Fine P, Barazani O. Environmental effects on molecular and phenotypic variation in populations of Eruca sativa across a steep climatic gradient. Ecol Evol. 2013;3:2471–84.View ArticlePubMedPubMed CentralGoogle Scholar
- Peleg Z, Saranga Y, Krugman T, Abbo S, Nevo E, Fahima T. Allelic diversity associated with aridity gradient in wild emmer wheat populations. Plant Cell Environ. 2008;31:39–49.PubMedGoogle Scholar
- Martinez-Ainsworth NE, Tenaillon MI. Superheroes and masterminds of plant domestication. Comptes Rendus Biologies. 2016;339:268–73.View ArticlePubMedGoogle Scholar
- Abbo S, Pinhasi van-Oss R, Gopher A, Saranga Y, Ofner I, Peleg Z: Plant domestication versus crop evolution: a conceptual framework for cereals and grain legumes. Trends Plant Sci 2014, 19:351–360.Google Scholar
- Abbo S, Gopher A. Near eastern plant domestication: a history of thought. Trends Plant Sci. 2017;22:491–511.View ArticlePubMedGoogle Scholar
- Vergeer P, Kunin WE. Life history variation in Arabidopsis lyrata across its range: effects of climate, population size and herbivory. Oikos. 2011;120:979–90.View ArticleGoogle Scholar
- Aronson J, Kigel J, Shmida A, Klein J. Adaptive phenology of desert and Mediterranean populations of annual plants grown with and without water stress. Oecologia. 1992;89:17–26.View ArticlePubMedGoogle Scholar
- Zohary D. Wild genetic-resources of crops in Israel. Isr J Bot. 1983;32:97–127.Google Scholar
- Stepansky A, Kovalski I, Perl-Treves R. Intraspecific classification of melons (Cucumis melo L.) in view of their phenotypic and molecular variation. Plant Syst Evol. 1999;217:313–32.View ArticleGoogle Scholar
- Paris HS. Semitic-language records of snake melons (Cucumis melo, Cucurbitaceae) in the medieval period and the "piqqus" of the "faqqous". Genet Resour Crop Evol. 2012;59:31–8.View ArticleGoogle Scholar
- Mendİ YY, Eldoğan S, Gutakev R. Regeneration and histological analysis of snake melon (Cucumis melo var. flexuosus ( L .) Naudin ) by direct organogenesis. Turk J Agric For. 2010;34:309–17.Google Scholar
- Hammer K, Hanelt P, Perrino P. Carosello and the taxonomy of Cucumis melo L. especially of its vegetable races. Kulturpflanze. 1986;34:249–59.View ArticleGoogle Scholar
- Ali-Shtayeh MS, Jamous RM, Shtaya MJ, Mallah OB, Eid IS, SYA Z. Morphological characterization of snake melon (Cucumis melo var. flexuosus) populations from Palestine. Genetic Resources and Crop Evolution. 2017;64:7–22.View ArticleGoogle Scholar
- Steane DA, Potts BM, McLean E, Collins L, Prober SM, Stock WD, Vaillancourt RE, Byrne M. Genome-wide scans reveal cryptic population structure in a dry-adapted eucalypt. Tree Genet Genomes. 2015;11Google Scholar
- Sanchez-Sevilla JF, Horvath A, Botella MA, Gaston A, Folta K, Kilian A, Denoyes B, Amaya I. Diversity arrays technology (DArT) marker platforms for diversity analysis and linkage mapping in a complex crop, the octoploid cultivated strawberry (Fragaria x ananassa). PLoS One. 2015;10Google Scholar
- Rogers SO, Bendich AJ. Extraction of DNA from milligram amounts of fresh, herbarium and mummified plant-tissues. Plant Mol Biol. 1985;5:69–76.View ArticlePubMedGoogle Scholar
- Al-Beyroutiova M, Sabo M, Sleziak P, Dusinsky R, Erik Bircak E, Hauptvogel P, Kilian A, Svec M. Evolutionary relationships in the genus Secale revealed by DArTseq DNA polymorphism. Plant Syst Evol. 2016;8:1083–91.View ArticleGoogle Scholar
- Garcia-Mas J, Benjak A, Sanseverino W, Bourgeois M, Mir G, Gonzalez VM, Henaff E, Camara F, Cozzuto L, Lowy E, et al. The genome of melon (Cucumis melo L.). Proc Natl Acad Sci U S A. 2012;109:11872–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Peakall R, Smouse PE. GenAlEx 6.5: genetic analysis in excel. Population genetic software for teaching and research-an update. Bioinformatics. 2012;28:2537–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Koyumdjisky H, Dan J, Soriano S, Nissim S. Selected soil profiles from Israel (in Hebrew). Bet Dagan: Volcani Center; 1988.Google Scholar
- Rowell DL. Soil science methods and applications. USA: Longman Scientific and Technical; 1994.Google Scholar
- Harel-Beja R, Tzuri G, Portnoy V, Lotan-Pompan M, Lev S, Cohen S, Dai N, Yeselson L, Meir A, Libhaber SE, et al. A genetic map of melon highly enriched with fruit quality QTLs and EST markers, including sugar and carotenoid metabolism genes. Theoret Appl Genetics. 2010;121:511–33.View ArticleGoogle Scholar
- Monforte AJ, Eduardo I, Abad S, Arus P. Inheritance mode of fruit traits in melon: Heterosis for fruit shape and its correlation with genetic distance. Euphytica. 2005;144:31–8.View ArticleGoogle Scholar
- Perin C, Hagen S, De Conto V, Katzir N, Danin-Poleg Y, Portnoy V, Baudracco-Arnas S, Chadoeuf J, Dogimont C, Pitrat M. A reference map of Cucumis melo based on two recombinant inbred line populations. Theor Appl Genet. 2002;104:1017–34.View ArticlePubMedGoogle Scholar
- Pandey S, Dhillon NPS, Sureja AK, Singh D, Malik AA. Hybridization for increased yield and nutritional content of snake melon (Cucumis melo L. var. flexuosus). Plant Genetic Res-Characterization and Utilization. 2010;8:127–31.View ArticleGoogle Scholar
- Munger HM, Robinson RW. Nomenclature of Cucumis melo L. Cucurbit Genet Coop Reports. 1991;14:43–4.Google Scholar
- Sebastian P, Schaefer H, Telford IRH, Renner SS. 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 Nat Acad Sci USA. 2010;107:14269–73.View ArticlePubMedGoogle Scholar
- Zhang SP, Miao H, Sun RF, Wang XW, Huang SW, Wehner TC, Gu XF. Localization of a new gene for bitterness in cucumber. J Heredity. 2013;104:134–9.View ArticleGoogle Scholar
- Shang Y, Ma YS, Zhou Y, Zhang HM, Duan LX, Chen HM, Zeng JG, Zhou Q, Wang SH, Gu WJ, et al. Biosynthesis, regulation, and domestication of bitterness in cucumber. Science. 2014;346:1084–8.View ArticlePubMedGoogle Scholar
- Abdel-ghani AH, Mahadeen A. Genetic variation in snake melon (Cucumis melo var. flexuosus ) populations from Jordan using morphological traits and RAPDs. Jordan Journal of Agricultural Sciences. 2014;10:96–119.View ArticleGoogle Scholar
- Rajagopalan PA, Perl-Treves R. Improved cucumber transformation by a modified explant dissection and selection protocol. Hortscience. 2005;40:431–5.Google Scholar
- Roorkiwal M, von Wettberg EJ, Upadhyaya HD, Warschefsky E, Rathore A, Varshney RK. Exploring germplasm diversity to understand the domestication process in Cicer spp. using SNP and DArT Markers. PLoS One. 2014:9.Google Scholar
- Grzebelus D, Iorizzo M, Senalik D, Ellison S, Cavagnaro P, Macko-Podgorni A, Heller-Uszynska K, Kilian A, Nothnagel T, Allender C, et al. Diversity, genetic mapping, and signatures of domestication in the carrot (Daucus carota L.) genome, as revealed by diversity arrays technology (DArT) markers. Mol Breeding. 2014;33:625–37.View ArticleGoogle Scholar
- Baloch FS, Alsaleh A, Shahid MQ, Ciftci V, ESdM L, Aasim M, Nadeem MA, Aktas H, Ozkan H, Hatipoglu R. A Whole genome dartseq and SNP analysis for genetic diversity assessment in durum wheat from Central Fertile Crescent. PLoS One. 2017;12:e0167821.View ArticlePubMedPubMed CentralGoogle Scholar
- Millet EJ, Welcker C, Kruijer W, Negro S, Coupel-Ledru A, Nicolas SD, Laborde J, Bauland C, Praud S, Ranc N, et al. Genome-wide analysis of yield in Europe: allelic effects vary with drought and heat scenarios. Plant Physiol Plant Mol Biol. 2016;172:749–64.Google Scholar
- Gerard GS, Kobiljski B, Lohwasser U, Börner A, Simón MR. Genetic architecture of adult plant resistance to leaf rust in wheat association mapping panel. Plant Pathol. 2018;67:584–594.Google Scholar
- Garavito A, Montagnon C, Guyot R, Bertrand B. Identification by the DArTseq method of the genetic origin of the Coffea canephora cultivated in Vietnam and Mexico. BMC Plant Biol. 2016;16Google Scholar