Characterization of Lebanese Germplasm of Snake Melon ( Cucumis melo subsp. melo var. ﬂexuosus ) Using Morphological Traits and SSR Markers

: Snake melon ( Cucumis melo subsp. melo L. var. ﬂexuosus (L.) Naudin) is an ancient and traditional crop in the Mediterranean region. Nevertheless, there has been poor interest in assessing snake melon germplasm where its genetic resources have not been surveyed before despite their potential in adaptation to environmental changes. In this study, we assess the genetic diversity of snake melon landraces collected from di ﬀ erent Lebanese regions at both morphological and molecular levels. Morphological characterization using a set of 18 descriptors revealed an important phenotypic variability among the landraces studied. Principle component analysis indicated that fruit hair and its consistency, fruit size, and skin color pattern were good criteria for discriminating among landraces. Based on the scatter plot diagram, landraces of snake melon formed ﬁve di ﬀ erent groups with one being deﬁned as typical var. ﬂexuosus . Ten simple sequence repeat (SSR) markers were used for the molecular characterization. Fifty-six di ﬀ erent alleles were detected, with an average of 5.6 alleles per locus. Polymorphism information content of SSR markers ranged from 0.06 to 0.84 (average 0.38). Cluster analysis based on molecular markers showed high genetic diversity and divided the landraces into ﬁve distinct genetic groups, conﬁrming thereby the morphological variability. Findings of this study indicate a signiﬁcant diversity for the Lebanese snake melon germplasm that must be further conserved and considered in improvement programs of this ancient crop.


Introduction
Snake melon or Armenian cucumber, also known as snake cucumber (Cucumis melo subsp. melo L. var. flexuosus (L.) Naudin; 2n = 2x = 24), is a member of the cucurbits family. It is an ancient crop that has been enjoyed and esteemed in Mediterranean antiquity by several civilizations [1]. Snake melon is a rain-fed crop where it flourishes in sunny and warm climates and it is commonly grown in the Levant area, Asia Minor, and North Africa. Snake melon fruits are harvested immature and green and consumed when fresh or pickled [2,3]. It earned its common name "snake melon" from its very long, slender, curved, and twisted fruits [4]. However, it is recognized around the world by different

Morphological Characterization
For the morphological evaluation 10 seeds from each accession were sown in a compost media. At the 2-leaf stage, seedlings were transplanted to the field. Characterization of the snake melon accessions was done at the production stage based on a list of descriptors previously developed by IPGRI, Stepansky et al., and Youssif et al. [17,25,26]. For each accession, a set of 15 immature fresh fruits were randomly taken five days after fruit set. In total, 18 traits were studied including 12 qualitative and 6 quantitative ones ( Table 2): pedicel separation and length, predominant and secondary skin color along with the color pattern covering the skin, presence or absence of fruit hair and their consistency, skin glossiness, fruit shape, texture and weight, ribs appearance and number, flesh color, sweetness, fruit length from stem end to blossom, width at the broadest point, and seed weight per fruit. Also, flower sex expression was recorded for each accession.   Collected accessions were submitted to both morphological and molecular assessments that have been carried on at the experimental field and laboratories of Warsaw University of Life Science, thus within the in-kind collaboration and documented exchange of seed plant genetic resources between Lebanon and Poland.

Morphological Characterization
For the morphological evaluation 10 seeds from each accession were sown in a compost media. At the 2-leaf stage, seedlings were transplanted to the field. Characterization of the snake melon accessions was done at the production stage based on a list of descriptors previously developed by IPGRI, Stepansky et al., and Youssif et al. [17,25,26]. For each accession, a set of 15 immature fresh fruits were randomly taken five days after fruit set. In total, 18 traits were studied including 12 qualitative and 6 quantitative ones ( Table 2): pedicel separation and length, predominant and secondary skin color along with the color pattern covering the skin, presence or absence of fruit hair and their consistency, skin glossiness, fruit shape, texture and weight, ribs appearance and number, flesh color, sweetness, fruit length from stem end to blossom, width at the broadest point, and seed weight per fruit. Also, flower sex expression was recorded for each accession.

DNA Isolation
Seedlings of snake cucumber accessions were grown in greenhouse under controlled conditions. At the 2-3 leaf stage, young leaf samples were collected from four seedlings from each accession and directly frozen in liquid nitrogen. Genomic DNA was extracted from the leaf samples by using DNeasy Plant Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer's protocol. The quality and quantity of the isolated DNA was checked spectrophotometrically using Nanodrop 2000 (TermoFisher TM , Waltham, MA, USA). All DNA samples were diluted using distilled water to equal concentration and stored at −20 • C.

SSR Analysis
Ten simple sequence repeat (SSR) primers previously used in similar analyses and/or developed for melon were selected for this study [27][28][29]. (Thermo Fisher Scientific, USA). The amplification was performed in a Master Gradient thermal cycler (Eppendorf, Hamburg, Germany). The cycling conditions were as follows: an initial denaturation at 94 • C for 2 min; 35 cycles of 2 min at 94 • C of denaturation, 1 min of primer annealing at 55 • C, and 2 min of extension at 72 • C; the program ended with a 10-min final extension step at 72 • C. PCR products were stored at 4 • C prior to analysis. After amplification, 6 µL of loading buffer (95% formamide, 10 mM EDTA (pH 8.0), 0.025% xylene cyanol, and 0.025% bromophenol blue) was added to each reaction tube. The samples were heat denatured for 5 min at 95 • C and quickly cooled on ice. To detect polymorphism, 1 µL of the amplified product was loaded on a 6% denaturation polyacrylamide gel that had been preheated for 30 min and visualized by silver staining. DNA fragment size was estimated using a 50-350 DNA ladder [30,31].

Statistical Analysis
For morphological characterization and with the aim to compare together qualitative and quantitative data, a numerical transformation was applied to all traits. For quantitative traits, arithmetic mean, standard deviation (s.d.) and coefficient of variation were calculated. The descriptive statistics were later used to classify the accessions into three classes according to Zar [32]; small (<Mean − s.d.), medium (>Mean − s.d.) and large (>Mean + s.d.). The Shannon-Weaver diversity index (H ) for qualitative and quantitative traits was estimated as follows [33]: where, pi is the proportion of individuals in the i th class of the n class trait and n is the number of phenotypic classes for a given trait.
Principal component analysis (PCA) was performed to generate a cluster diagram. Eigen value and contribution percentage of the first three principal component axes were calculated using the correlation matrix among 15 descriptors for 66 genotypes. All computations were performed using R software version 3.6.1 [9].
Bands revealed by SSR markers were scored as absent (0) or present (1), and then data were used to generate a pair-wise similarity matrix using Jaccard's coefficient [34]. The unweighted pair-group method with arithmetic mean (UPGMA) was employed to create the clustering dendrogram [35]. Polymorphism information content (PIC) values were calculated according to Smith et al. [36], using the algorithm for all primer combinations as follows:

General Status of Collected Accessions
A total of 70 snake melon accessions were collected during this survey. Figure 1 and Table 1 show the distribution of the sites visited during the collection missions. Eight accessions were found in Akkar, 10 in North Lebanon, 13 in Mount Lebanon, 11 in Southern Lebanon, and 28 in the Bekaa region. Sources of plant material were obtained from farmers and commercial agricultural stores. The accessions were collected from rural areas, where farmers apply the seed saving system and plant small or vast surfaces of snake melon for their specific uses. When accessions were collected from local farmers' seed stores, efforts were made to obtain seeds from local conservation systems rather than improved materials from breeding programs. For the 70 collected accessions, only eight vernacular names were found referring mainly to the color of the skin or to the shape. About 38 of the 70 studied accessions were named 'Abyad' (meaning white), 24 were named 'Akhdar' (meaning green), two were named 'Azrak' (meaning blue), and two were named 'Mshakal' (meaning mixed colors). In association to the fruit shape, four accessions were found namely 'Mdalaa' (meaning ribbed), 'Hayye' (meaning snake), 'Tawile' (meaning long), and 'Chikawi' (meaning thorny), represented by one accession each. It is worthy to note that there was no apparent link between the given names for accessions and their geographic distribution. However, the most diversity of vernacular names was observed in the Bekaa area (five names), followed by Mount Lebanon (four names).

Morphological Characterization
Over the 18 observed traits (Table 2), three were monomorphic and 15 were polymorphic. Studied landraces were found to be monecious and monomorphic in fruit skin glossiness and texture. They had various types of fruits but mostly with ribbed and dull skin ( Figure 2).
in Akkar, 10 in North Lebanon, 13 in Mount Lebanon, 11 in Southern Lebanon, and 28 in the Bekaa region. Sources of plant material were obtained from farmers and commercial agricultural stores. The accessions were collected from rural areas, where farmers apply the seed saving system and plant small or vast surfaces of snake melon for their specific uses. When accessions were collected from local farmers' seed stores, efforts were made to obtain seeds from local conservation systems rather than improved materials from breeding programs.
For the 70 collected accessions, only eight vernacular names were found referring mainly to the color of the skin or to the shape. About 38 of the 70 studied accessions were named 'Abyad' (meaning white), 24 were named 'Akhdar' (meaning green), two were named 'Azrak' (meaning blue), and two were named 'Mshakal' (meaning mixed colors). In association to the fruit shape, four accessions were found namely 'Mdalaa' (meaning ribbed), 'Hayye' (meaning snake), 'Tawile' (meaning long), and 'Chikawi' (meaning thorny), represented by one accession each. It is worthy to note that there was no apparent link between the given names for accessions and their geographic distribution. However, the most diversity of vernacular names was observed in the Bekaa area (five names), followed by Mount Lebanon (four names).

Morphological Characterization
Over the 18 observed traits (Table 2), three were monomorphic and 15 were polymorphic. Studied landraces were found to be monecious and monomorphic in fruit skin glossiness and texture. They had various types of fruits but mostly with ribbed and dull skin ( Figure 2).  Descriptive statistics including mean, minimum, maximum, and coefficient of variation (CV) were recorded for each quantitative trait measured (Table 3). Among all quantitative traits recorded, pedicel length, fruit number of ribs, sweetness, size, length/width ratio, and seed weight, showed relatively low CV values (<45%) ranging from 21.19% to 32.31%. All quantitative traits showed medium to high H diversity index ranging from 0.51 for fruit number of ribs and 0.99 for pedicel length. The frequency distribution of polymorphic phenotypic traits for the nine qualitative traits is shown in Figure 3. Two phenotypic classes were observed according to pedicel separation from fruit, primary fruit skin color, secondary fruit skin color pattern, fruit hair, and fruit flesh color. Easy pedicel separation, white primary color, the absence of fruit skin color pattern and pubescence, and green flesh color were more predominant among all accessions.
shown in Figure 3. Two phenotypic classes were observed according to pedicel separation from fruit, primary fruit skin color, secondary fruit skin color pattern, fruit hair, and fruit flesh color. Easy pedicel separation, white primary color, the absence of fruit skin color pattern and pubescence, and green flesh color were more predominant among all accessions.
On the other hand, three phenotypic classes were generated by secondary skin fruit color, fruit shape, hair consistency, and fruit ribs appearance. The white secondary color, semi-curved fruit shape, absence in hair consistency, and strong appearance of fruit ribs were predominantly observed in the studied accessions.
The H' diversity index among classes of the studied traits was calculated ( Figure 3). Overall, the qualitative traits showed intermediate to high diversity indices. The highest H' was observed with fruit shape (0.96) followed by hair consistency (0.92), while the lowest diversity among trait classes was obtained with fruit flesh color (0.57).

Principal Component Analysis
Principal component analysis (PCA) involving 15 polymorphic morphological traits (nine qualitative and six quantitative descriptors) was performed to identify the most discriminating characteristics ( Table 4). The first three components present 51.51% of the total variation. The first component is characterized by a percentage of variation of 28.39% and is mainly represented by fruit hair and its consistency, skin primary and secondary colors and their respective pattern, and fruit ribs appearance. The second component explains 12.42% of the total variation and includes principally fruit size, flesh color and length to width ratio. The third component was characterized by a percentage of 10.69% and is dominated by fruit shape and sweetness.
The first and second components distribute the 70 Lebanese snake melon accessions into four main groups (Figure 4). All accessions of group I have large size fruits with hairy skin, while accessions of group II have small to medium size fruits with hairy fruits. However, group III have glabrous fruits with a striped skin color pattern, and accessions of group IV have glabrous fruits with the absence of fruit skin color pattern. Within group IV, a small subgroup of seven accessions (B5, B7, B13, B14, B18, B34, and B41), all coming from Bekaa province, could be distinguished by their long and slender fruits, strong fruit ribs, white flesh color, and slight fruit sweetness.  On the other hand, three phenotypic classes were generated by secondary skin fruit color, fruit shape, hair consistency, and fruit ribs appearance. The white secondary color, semi-curved fruit shape, absence in hair consistency, and strong appearance of fruit ribs were predominantly observed in the studied accessions.
The H diversity index among classes of the studied traits was calculated ( Figure 3). Overall, the qualitative traits showed intermediate to high diversity indices. The highest H was observed with fruit shape (0.96) followed by hair consistency (0.92), while the lowest diversity among trait classes was obtained with fruit flesh color (0.57).

Principal Component Analysis
Principal component analysis (PCA) involving 15 polymorphic morphological traits (nine qualitative and six quantitative descriptors) was performed to identify the most discriminating characteristics ( Table 4). The first three components present 51.51% of the total variation. The first component is characterized by a percentage of variation of 28.39% and is mainly represented by fruit hair and its consistency, skin primary and secondary colors and their respective pattern, and fruit ribs appearance. The second component explains 12.42% of the total variation and includes principally fruit size, flesh color and length to width ratio. The third component was characterized by a percentage of 10.69% and is dominated by fruit shape and sweetness. The first and second components distribute the 70 Lebanese snake melon accessions into four main groups ( Figure 4). All accessions of group I have large size fruits with hairy skin, while accessions of group II have small to medium size fruits with hairy fruits. However, group III have glabrous fruits with a striped skin color pattern, and accessions of group IV have glabrous fruits with the absence of fruit skin color pattern. Within group IV, a small subgroup of seven accessions (B5, B7, B13, B14, B18, B34, and B41), all coming from Bekaa province, could be distinguished by their long and slender fruits, strong fruit ribs, white flesh color, and slight fruit sweetness.
Also, it was noticed that the majority of the accessions collected from South Lebanon are present in group I, while accessions from North Lebanon and Mount Lebanon were mainly scattered in groups II and III. Moreover, a large number of the studied accessions coming from Bekaa and Akkar provinces are grouped in cluster IV.

Genetic Similarity
With the 10 SSRs markers used, a total of 56 alleles were detected among the 70 Lebanese snake melon accessions ( Table 5). The number of alleles identified by each primer pair ranged from 2 in GMGS2-3 to 11 in CMTA134a with an average alleles number of 5.71. Size of allele fragments ranged from 59 to 340 bp. Moreover, all SSR markers were found to be polymorphic (polymorphism rate: 100%) among snake melon accessions. PIC values for SSRs ranged from 0.06 with GMGS2-3 to 0.84 with CMTA134a and the average PIC value for all the studied primer pairs was 0.38.

Genetic Similarity
With the 10 SSRs markers used, a total of 56 alleles were detected among the 70 Lebanese snake melon accessions ( Table 5). The number of alleles identified by each primer pair ranged from 2 in GMGS2-3 to 11 in CMTA134a with an average alleles number of 5.71. Size of allele fragments ranged from 59 to 340 bp. Moreover, all SSR markers were found to be polymorphic (polymorphism rate: 100%) among snake melon accessions. PIC values for SSRs ranged from 0.06 with GMGS2-3 to 0.84 with CMTA134a and the average PIC value for all the studied primer pairs was 0.38. In all accessions tested more than one allele was observed in at least three SSR loci analyzed. Generally multiple alleles were observed in four or five SSR loci (in 36 and 18 accessions, respectively). The individual SSR alleles occurred in between one and 67 of the tested accessions. Overall, 11 alleles were rare (occurred in less than 5% of the accessions), 27 were common (in 5 to 25% of accessions) and 18 alleles were frequent (over 50% of accessions).

Discussion
This is the first report on the assessment of the genetic diversity of the Lebanese germplasm of snake melon by means of both morphological descriptors and SSR markers. The morphological characterization of the Lebanese snake melon accessions was done using 18 quantitative and qualitative morphological descriptors. All accessions were monoecious and had a dull and ribbed fruit skin. Overall, all traits revealed intermediate to high H' diversity index ranging from 0.57 to 0.96 for the qualitative traits and between 0.51 and 0.99 for quantitative traits. These relatively high diversity indices are mainly due to the existence of more than two phenotypic classes as the case in fruit hair consistency, or to the almost even distribution of individuals in the different classes of a studied trait as the case in fruit shape. The principle component analysis (PCA) indicated that fruit pubescence and its consistency, fruit size, along with skin color pattern were the most relevant traits to explain the variability within the studied accessions. The percent of the total variability found in our study was 51.51% as accounted by the first three components. This result was higher than that previously obtained in Jordanian snake melon accessions (39.9%) [13]. However, higher percentages were observed for Palestinian snake melon landraces (69.9%) [9] and in Tunisian snake melon accessions (89.8%) [37].
The PCA scatter diagram clustered the accessions into four main groups according to the presence (group I and II) or absence (group III and IV) of hair. However, the seven accessions collected from farmers in the Bekaa area present in the small subgroup of the major group IV produced elongated, ribbed, hairless, white skinned and flesh fruits, with a Brix value less than 4.2°, and length-to-width ratio higher than 8. These characters agree well with the description of snake melon varieties by Pitrat et al. [4]. Thus, these accessions are considered to be typical and original flexuous landraces. It has become very clear that the method of seed-saving practiced by local farmers The second main group was also portioned into two subgroups as the second node. Fifteen accessions coming from Bekaa, Mount Lebanon, South and North Lebanon, namely ML3, ML4, B6, B32, B18, B25, ML20, SL22, B35, ML12, B30, ML19, NL49, SL16, NL21, were found in the first subgroup G3, of which B35, ML12, and B30 accessions were clustered very closely together having genetic similarity of 0.94. The second subgroup G4 was subdivided into two main sub-groups G4A and G4B of the third node containing the remaining studied accessions. The first sub-group G4A constituted of SL15, B28, ML1, and ML2 accessions, and the second sub-group G4B constituted from 21 accessions, B27, NL46, ML53, B31, A55, SL24, NL47, NL48, NL52, B9, B10, B14, SL37, ML43, B62, SL23, A57, B67, ML59, A68, and A70.

Discussion
This is the first report on the assessment of the genetic diversity of the Lebanese germplasm of snake melon by means of both morphological descriptors and SSR markers. The morphological characterization of the Lebanese snake melon accessions was done using 18 quantitative and qualitative morphological descriptors. All accessions were monoecious and had a dull and ribbed fruit skin. Overall, all traits revealed intermediate to high H diversity index ranging from 0.57 to 0.96 for the qualitative traits and between 0.51 and 0.99 for quantitative traits. These relatively high diversity indices are mainly due to the existence of more than two phenotypic classes as the case in fruit hair consistency, or to the almost even distribution of individuals in the different classes of a studied trait as the case in fruit shape. The principle component analysis (PCA) indicated that fruit pubescence and its consistency, fruit size, along with skin color pattern were the most relevant traits to explain the variability within the studied accessions. The percent of the total variability found in our study was 51.51% as accounted by the first three components. This result was higher than that previously obtained in Jordanian snake melon accessions (39.9%) [13]. However, higher percentages were observed for Palestinian snake melon landraces (69.9%) [9] and in Tunisian snake melon accessions (89.8%) [37].
The PCA scatter diagram clustered the accessions into four main groups according to the presence (group I and II) or absence (group III and IV) of hair. However, the seven accessions collected from farmers in the Bekaa area present in the small subgroup of the major group IV produced elongated, ribbed, hairless, white skinned and flesh fruits, with a Brix value less than 4.2 • , and length-to-width ratio higher than 8. These characters agree well with the description of snake melon varieties by Pitrat et al. [4]. Thus, these accessions are considered to be typical and original flexuous landraces. It has become very clear that the method of seed-saving practiced by local farmers since antiquity has maintained genetically diverse crops where farmers continuously selected better snake melon features to overcome the change in the agroclimatic conditions and/or to meet the consumers' preferences. As to the remaining 63 studied accessions growing in the different provinces, the large variability revealed within this set is probably the result of outcrossing nature of the snake melon. Similarly, Jordanian snake melon germplasm exhibited a genetic variation among populations for only 18.52% of the total variation, while more than 80% occurred within populations due to the outcrossing nature of this crop [13].
In addition, it was noticed that accessions coming from the same province have similar morphological characteristics. Majority of accessions collected from Bekaa and Akkar were characterized by their white color and hairless skin, while a significant number of accessions growing in the North, South, and Mount Lebanon are glabrous. Although the consumer's preference of snake melon differs among Lebanese regions, accessions from different geographic sources were found to be clustered in the same group. This could be the result of seeds dissemination that is actively practiced for landraces through the informal farmers' network and the traditional seed stores [21].
The molecular analysis of the 70 snake melon accessions performed through the 10 SSR markers generated a total of 56 alleles with two to 11 alleles per primer pair and an average number of 5.6 alleles per primer pair. This result is in line with the one previously obtained for Turkish snake melon accessions by Solmaz et al. [15], while Henane et al. [18] had a much lower value of 2.54 alleles per primer when studying the diversity of the Tunisian snake melon genotypes. Also, our average number of alleles per primer was higher than the one obtained by Ning et al. [20] through 36 SSR markers (4.03 alleles/primer) used to characterize the Chinese snake melon and its relationship with melon germplasms of diverse origins. Conversely, higher diversity was detected in melon germplasms for which the primer pairs were initially developed. Monforte et al.
[38] identified 6.3 alleles on a collection of 27 accessions of wild and cultivated melons. Likewise, Kacar et al. [19] reported 6.55 alleles on 96 Turkish melon genotypes belonging to various subspecies. On the other hand, the polymorphism percentage obtained here with the Lebanese snake melon accessions was 100%, similarly to the findings of Solmaz et al. [15] for the Turkish snake melon germplasm and the ones of Tzitzikas et al. [16] for the Greek/Cypriot germplasm. Nevertheless, Tzitzikas et al. [16] reported a polymorphism information content (PIC) average value of 0.44, which is higher than the 0.38 value obtained here. This difference might be due to the different germplasms studied, or to the selected set of SSR markers used that might be more informative when compared to our study. Also, the multiple number of alleles per locus observed in this study is suggested to be linked to the heterogeneity and/or heterozygosity of the plants sampled to represent the respective accessions that may be caused by the outcrossing nature of this crop [13].
In addition, the polymorphism per locus in the used markers varied between our results and previous findings. For instance, in marker GMAGN68, 11 different alleles were detected by Solmaz et al. [15], while only seven alleles were observed in our study, of which four alleles (168, 180, 184, and 187 bp) were found to be the same in both studies. In marker GMGS2-3 there were six different alleles observed while assessing the Turkish germplasm, while only two alleles were detected in our study. However, these two alleles were found at higher molecular weights with 320 and 315 bp as compared to the Turkish alleles ranging between 190 and 221 bp. Marker CMTA134a identified six alleles while assessing the Greek/Cypriot germplasm versus a total of 11 different alleles detected in our study. These findings show that the Lebanese snake melon germplasm have some distinctive characters as compared to other germplasms in the region.
The SSR markers used here allowed us to differentiate the studied Lebanese snake melon accessions into five main clusters. The similarity among accessions according to Jaccard's similarity distance coefficient ranged from 0.25 to 0.96. This proved the existence of wide diversity among Lebanese snake melon accessions studied.
Although both SSR markers and morphological descriptors revealed considerable variation among the Lebanese snake melon accessions, the clustering generated by SSRs was poorly connected to the one made by the morphological traits from one side, and to the geographic origin of studied accessions from the other side. Such weak connection could be due to the absence of linkage between SSR markers chosen for this study and loci that control the traits [38]. Similar results were also observed in the reports of Kacar et al. [19] and Solmaz et al. [15] indicating no correlation between molecular clustering and major morphological traits from one side and geographical origin of melon and snake melon genotypes from the other side. It is worthy to note that Kacar et al. [19] analyzed 12 snake melon accessions, along with other melon accessions, that were fairly clustered together but did not have any relation to geographic origin of the studied accessions, and Solmaz et al. [15] did not identify the geographic region and source of the used Turkish snake melon accessions.
Nevertheless, the three cases of similar accessions commonly named 'Abyad' as revealed by SSR analysis at 0.96 level of similarity, had also similar fruit characteristics and were very closely related in the morphological scatter plot. Specifically, the similar ML12, B30, and B35 accessions had white, semi curved, hairless, medium length fruits. B63 and B66 had white, elongated, hairless, short fruits, and SL40, and A69 had white, curved, hairless, long fruits. It is worthy to note that these 'Abyad' accessions that sound similar are sold in farmers' seed stores and disseminated in different agroclimatic areas in terms of altitude and precipitations where they are cultivated under different agricultural practices.
The crop vegetable diversification in the Mediterranean area is a complex matter [39]. The morphological assessment showed that some classes among studied characteristics were dominant. These include the white primary (71%) and secondary (74%) color, absence of pubescence (57%), the elongated to semi curved fruit shape (88%) and the medium sized fruits with respect to weight (63%) and dimensions (71%). Interestingly, the majority of the traits observed in this study are of various interest to consumers, growers, and breeders. For instance, some desirable phenotypic criteria of the fruit including elongated shape, sweet taste, white skin, and crisp or soft texture are critical for fresh snake melon consumption, while growers are mainly interested in marketable characters related to fruit dimensions and weight [26]. Breeders tend to search in the existing germplasm for desirable traits to improve the fruit appearance and plant resistance to biotic and abiotic stresses [40].
The overall analysis proves the existence of a wide diversity in plant material that may have important implications for genetic resources management. This high diversity could be linked to the old history of snake melon cultivation in Lebanon and around the Mediterranean [1]. Moreover, cultural practices like seed conservation and the necessity to protect rain-fed crops may have increased diversification of snake melon

Conclusions
This study highlighted the diversity of the Lebanese germplasm of snake melon as revealed by both morphological descriptors and SSR markers. Snake melon accessions assessed here represent a