Research | Open | Published:
Analysis of genetic diversity of Tunisian caprifig (Ficus carica L.) accessions using simple sequence repeat (SSR) markers
Hereditasvolume 152, Article number: 1 (2015)
The common fig (Ficus carica L.) is a gynodioecious species with two sexual forms: male trees (caprifigs) with male and female flowers and female trees that produce only female flowers that will result in the edible fig syconium. In this study the genetic diversity of 20 Tunisian accessions of caprifig is analyzed using SSR markers previously developed for this crop.
The results revealed that the 13 pairs of primers used amplified a total of 37 alleles in the accessions studied. The number of alleles per locus ranged from two to six, with a mean value of 2.85 alleles per locus. Observed and expected heterozygosities showed mean values of 0.33 and 0.29 respectively. UPGMA cluster analysis and Principal Component Analysis grouped the caprifig accessions analyzed in three groups.
The results obtained show a low genetic diversity in the Tunisian accessions of caprifig studied and, in spite of analyzing samples from different geographic regions, no clear groupings based on geographical origin are observed suggesting widespread exchange of caprifig plant material through vegetative propagation among different areas in Tunisia.
The common fig (Ficus carica L.; 2n = 2x = 26)  belongs to the Moraceae, a family with over 1400 species distributed in about 40 genera. The genus Ficus L. contains about 750 species of woody trees, epiphytes and shrubs, mainly of tropical and subtropical distribution, divided into six subgenera [2, 3] that share a unique inflorescence, the syconium. Common fig seems to be originated from Southern Arabia and the eastern part of the Mediterranean regions and, together with the grapevine and the olive, is considered one of the three classical fruit trees associated with the beginning of horticulture in the Mediterranean Basin [4–6] domesticated at a very early stage contemporarily with cereal crops .
Ficus carica L. is a gynodioecious species with two sexual forms: male trees (caprifigs) that produce syconia with separate male and female flowers and female trees that produce syconia with only female flowers that will develop into edible seeded figs if pollinated. Since only male trees produce pollen the common fig is functionally a dioecious species. Three types of female figs are cultivated : the common-type that develops fruit parthenocarpically without pollination and can produce one (unifera varieties) or two (bifera varieties) crops, the Smyrna-type that requires pollination with pollen from caprifigs, and the San Pedro-type that produces a first crop parthenocarpically (breba) and a second crop (fig) only after pollination with pollen from caprifigs. Pollination in the genus Ficus is dependent on the coevolution of Ficus species with pollinating wasps of the family Agaonidae . In the case of the common fig, pollination (caprification) is performed by a specific pollinating insect, Blastophaga psenes L. The caprifigs produce fruits in three crop cycles during each growing season, each harboring the larvae, pupae and, temporarily, the adults of the pollinating Blastophaga : “profichi” that ripen in early summer, “mammoni” that ripen in autumn and “mamme” that overwinter on the tree and mature in spring. Only the "profichi” carry pollen and are used for pollination. Pollination occurs naturally when female and caprifig trees are present in the same orchard or when caprifig branches with flowers are placed close to female fig trees. The main problems of caprification in Tunisia, which are common to other fig producing areas, are the disruption of the cycle of Blastophaga psenes in cold mountainous areas and the unavailability of mature caprifig “profichi” when female figs are receptive.
In Tunisia, as in other Mediterranean countries , fig has been traditionally cultivated since ancient times in diverse edaphoclimatic conditions, in association with date palms in the south or olive trees in other regions of the country resulting in a high number of local varieties and frequent exchange of varieties among different regions [11, 12]. The denominations of the cultivars are usually based on the color, size and time of fruit ripening or geographical origin resulting in confusion in nomenclature. Hence, appropriate characterization and differentiation among cultivars is necessary to optimize fig germplasm management and conservation, jeopardized by intensive urbanization, cultivation of selected clonal varieties or biotic and abiotic stresses. This genetic erosion is even more important in caprifigs since they do not produce edible fruits and, consequently, the pressure to conserve genetic resources is lower. However, artificial caprification is still a common practice in Tunisian fig cultivation and, consequently, caprifigs are used by farmers in order to obtain edible fig production .
Various studies have reported the use of morphological traits [11–19] and isozyme markers [20, 21] for fig characterization. However, these parameters are influenced by environmental conditions and the phenological status of the trees. To overcome these difficulties, different molecular tools such as RAPDs [13, 20–27] ISSRs [6, 19, 24, 26–28], AFLPs [20, 29, 30], RFLPs  or SSRs [6, 10, 26, 27, 30–41] have been used for fig germplasm characterization and diversity analyses. However, most of those studies include mainly female fig cultivars and molecular characterization and diversity studies in caprifigs are very scarce.
The main objective of this study was to characterize and evaluate the genetic diversity of 20 Tunisian caprifig accessions using SSR markers in order to develop strategies to preserve the endangered genetic resources of this species.
This study was carried out on 20 Tunisian caprifig accessions originated in different geographic regions and with different phenotypic traits (Table 1 and Fig. 1). All the accessions were planted in the fig germplasm collection of the Arid Land Institute of Médenine established in El Gordhab, Tataouine in Southeastern Tunisia. The plant material was propagated by hardwood cuttings.
Genomic DNA extraction and PCR amplification
Genomic DNA was extracted from lyophilized young leaves following a CTAB-based method optimized previously for fig . After purification, DNA quantity and quality were checked using a Nanodrop ND-1000 UV-visible spectrophotometer, diluted to 10 ng μl−1 with modified TE buffer (1 M Tris–HCl pH 8.0; 0.5 M EDTA) and used for PCR amplification.
A set of 13 SSR loci (MFC1, MFC2, MFAC4, LMFC15, LMFC18, LMFC21, LMFC23, LMFC24, LMFC27, LMFC28, LMFC30, LMFC31 and LMFC32), previously developed in fig [31, 33] were used in 15 μl reactions containing 16 mM (NH4)2SO4, 67 mM Tris–HCl, pH 8.8, 0.01 % Tween 20, 2 mM MgCl2, 0.1 mM each dNTP, 0.4 μM each primer, 25 ng genomic DNA and 0.5 Units of BioTaq DNA polymerase (Bioline, London, UK). Amplifications were performed in a thermocycler (Bio-Rad Laboratories, Hercules, CA, USA) using the following temperature profile: an initial step of 1 min at 94 °C, 35 cycles of 30 s at 94 °C, 30 s at 55 °C and 1 min at 72 °C, and a final step of 5 min at 72 °C. The amplification products were resolved using a CEQTM 2000XL capillary DNA analysis system (Beckman Coulter, Fullerton, CA, USA).
Different parameters of genetic diversity were estimated: number of alleles per locus (A), allelic frequencies, observed heterozygosity (Ho), expected heterozygosity (He), Wright’s fixation index and probability of identity (PI). Those parameters were only computed using the SSRs that produced amplification of a single locus (e.g. one or two alleles amplified in each of the samples analyzed). The computations were performed with the program POPGENE 1.32  and IDENTITY 1.0 . Genetic relationships within the accessions studied were calculated using UPGMA cluster analysis of the similarity matrix obtained from the proportion of shared amplification fragments  using the program NTSYS 2.11 (Exeter Software, Stauket, NY). The cophenetic correlation coefficient was estimated by comparing with the Mantel test the cophenetic matrix obtained from the dendrogram with the original similarity matrix. The bootstrap values were obtained using 2000 replicates with TREECON 1.3b . Principal Component Analysis (PCA) was performed using NTSYS 2.11(Exeter Software, Stauket, NY).
Results and discussion
Genetic polymorphism and SSR patterns
The 13 SSR loci produced successful and repeatable amplification fragments in all the 20 caprifig accessions analyzed, resulting in a total of 37 alleles ranging from two (LMFC32, LMFC15, LMFC21, LMFC31, LMFC18, LMFC27 and LMFC23) to six (LMFC30) alleles per locus (Table 2), with an average of 2.85 alleles per locus and amplification fragment sizes between 120 and 278 bp. Some accessions showed the amplification of more than two alleles (three alleles in LMFC30, MFC1, MFC4 and four alleles in LMFC28) suggesting the probable amplification of more than one locus. Similar results for those 4 loci were obtained by Giraldo et al.  analyzing an ex situ germplasm collection of 209 fig accessions. The remaining 9 SSRs showed one or two bands per genotype, suggesting the amplification of a single locus. For these 9 loci, the accessions studied were considered homozygous or heterozygous when one or two fragments were present per locus respectively . Genetic diversity was studied with the 9 loci that produced one or two alleles per locus. For these loci, allelic frequencies varied from 0.03 to 0.98 with a mean of 0.45 (data not shown). Furthermore, 15 % of the alleles studied were considered rare (p < 0.05) and fixed (p > 0.9). Additionally, six alleles were found only in one accession (allele 205 of LMFC32 is present only in ‘Jrani’; allele 278 of LMFC24 is present only in ‘Dhokkar2’; allele 271 of LMFC21 is present only in ‘Dhokkar3’; allele 256 of LMFC30 is present only in ‘Dhokkar1’ and allele 193 of LMFC 28 is present in ‘Bouharrag2’).
Observed heterozygosity ranged from 0.05 (LMFC32 and LMFC21) to 0.80 (LMFC31) with a mean of 0.33 (Table 2). The expected heterozygosity ranged from 0.05 (LMFC32 and LMFC21) to 0.50 (LMFC31) with a mean of 0.29 (Table 2). A heterozygote excess (Hobs > Hexp) was observed for the LMFC15, LMFC31, LMFC18, LMFC27 loci whereas a deficit of heterozygosity (Hobs < Hexp) was observed in LMFC24 and MFC2 loci (Table 2). For the loci LMFC 32, LMFC21 and LMFC23 observed and expected heterozygosities were similar (Table 2). The maximum value of the probability of identity (0.91) was detected in LMFC32 and LMFC21 and the minimum (0.45) in MFC2 (Table 2).
The diversity parameters obtained in this work are mostly below the range of those reported for fig using microsatellites in previous works that mainly analyzed female cultivars [6, 10, 26, 27, 30, 41]. Although most of those works included the analysis of a higher number of female genotypes from different origins, the results obtained suggest a low genetic diversity in the Tunisian caprifig accessions studied (see below).
Cluster and principal component analyses
Among all the possible UPGMA dendrograms generated, that with the highest cophenetic correlation (r = 0.76) between the cophenetic coefficient and the similarity matrix was chosen. The UPGMA dendrogram obtained showed three main groups among the genotypes analyzed (Fig. 2) with three undistinguishable accessions that represent a case of synonymy (Assafri from Kerkennah, Beldi and Dhokkar6 from Zarzis). The first group (A) includes 7 accessions (Magouli1, Bouharrag1, Khadhouri, Magouli2, Dhokkar4, Bithri2 and Bouharrag2), the second (B) contains 10 accessions (Jrani, Dhokkar7, Dhokkar1, Assafri, Beldi, Dhokkar6, Bouharrag3, Bithri1, Dhokkar3 and Dhokkar2) and the third (C) includes 3 accessions (Limi, Sawoudi and Tebessi). Several homonymies were detected in the genotypes analyzed since genotypes having the same name in the same or different locations are genetically different. This includes all the different samples analyzed of Magouli (2 samples), Bouharrag (3 samples), Bithri (2 samples) and Dhokkar (6 samples). Among all the nodes obtained only three groups have bootstrap values higher than 50 %. These groups include the three undistinguishable accessions (Assafri, Beldi and Dhokkar6) (bootstrap of 95 %), those three accessions and Bouharrag3 (bootstrap of 71 %) and the accessions Limi and Sawoudi (bootstrap of 62 %).
Principal Component Analysis shows that the three first principal components explain 50.50 % of the total variability. The contributions of PCA1, PCA2 and PCA3 were 21.33, 19.26 and 9.92 % respectively. Fig. 3 shows the distribution of accessions according to the first two components (PCA1 and PCA 2) in which three groups that correspond to those found in the UPGMA analysis can be clearly identified.
Some of the main groups obtained in the UPGMA and PCA studies are correlated with the geographic origin of the genotypes. Thus the group C in Figs. 2 and 3 corresponds to three cultivars (Limi, Tebessi and Sawoudi) from the same location in Southwestern Tunisia (Kébéli) and the group of three undistinguishable cultivars (Assafri, Beldi and Dhokkar6), which are probable synonymies, were originated in Eastern Tunisia. Six of the accessions, even from different geographic origins, wear the name ‘Dhokkar’; although this could be a case of homonymy, it has to be considered that ‘dhokkar’ is the common word for caprifig in Arabic. Regarding the rest of the accessions analyzed, although some of the ones that cluster together (such as Dhokkar2 and Dhokkar3) have a common geographical origin, no clear geographical groupings can be found for most of them. Similar results have also been obtained with common type figs in Tunisia using different types of molecular markers [24, 28] These results could suggest a common genetic base for most Tunisian caprifigs that can be explained by the easy vegetative propagation of the crop that allows exchange of plant material between different regions. This contrasts with a higher genetic diversity found in Tunisian common figs by characterizing cultivated and wild figs .
The present study provides the first molecular database of Tunisian caprifig accessions using SSR markers. The set of SSR markers used indicates that the genetic diversity between accessions studied is overall narrow and that no clear relationship is found between geographical origin and genetic composition suggesting exchange of caprifig genetic material among different regions. Additional studies with caprifigs from other countries should be performed in order to have a clear picture on overall caprifig genetic diversity and optimize collaborative caprifig genetic resource management and conservation.
Weiblen GD. Phylogenetic relationships of functionally dioecious Ficus (Moraceae) based on ribosomal DNA sequences and morphology. Am J Bot. 2000;87:1342–57.
Berg CC. Flora Malesiana precursor for the treatment of Moraceae 1: The main subdivision of Ficus: the subgenera. Blumea. 2003;48:167–78.
Ronsted N, Weiblen GD, Clement WL, Zerega NJC, Savolainen V. Reconstructing the phylogeny of figs (Ficus, Moraceae) to reveal the history of the fig pollination mutualism. Symbiosis. 2008;45:45–55.
Zohary D, Spiegel-Roy P. Beginnings of fruit growing in the old world. Science. 1975;187:319–27.
Janick J. The origin of fruits, fruit growing, and fruit breeding. Plant Breed Rev. 2005;25:255–320.
Khadari B, Oukabli A, Ater M, Mamouni A, Roger JP, Kjellberg F. Molecular characterization of Moroccan fig germplasm using intersimple sequence repeat and simple sequence repeat markers to establish a reference collection. HortSci. 2005;40:29–32.
Kislev ME, Hartmann A, Bar-Yosef O. Early domesticated fig in the Jordan Valley. Science. 2006;312:1372–4.
Flaishman M, Rodov V, Stover E. The fig: botany, horticulture and breeding. Hortic Rev. 2008;34:113–96.
Herre EA, Jander KC, Machado CA. Evolutionary ecology of figs and their associates: Recent progress and outstanding puzzles. Ann Rev Ecol Evol Syst. 2008;39:439–58.
Achtak H, Oukabli A, Ater M, Santoni S, Kjellberg F, Khadari B. Microsatellite markers as reliable tools for fig cultivar identification. J Amer Soc Hort Sci. 2009;134:624–31.
Chatti K, Salhi-Hannachi A, Mars M, Marrakchi M, Trifi M. Analyse de la diversité de cultivars tunisiens de figuier (Ficus carica L.) par les caractères morphologiques. Fruits. 2004;59:49–61.
Mars M, Chatti K, Saddoud A, Salhi-Hannachi A, Trifi M, Marrakchi M. Fig cultivation and genetic resources in Tunisia, an overview. Acta Hortic. 2008;798:27–32.
Papadopoulou K, Ehaliotis C, Tourna M, Kastanis P, Karydis I, Zervakis G. Genetic relatedness among dioecious Ficus carica L. cultivars by random amplified polymorphic DNA analysis, and evaluation of agronomic and morphological characters. Genetics. 2002;114:183–94.
Caliskan O, Polat AA. Fruit characteristics of fig cultivars and genotypes grown in Turkey. Sci Hortic. 2008;115:360–7.
Aljane F, Ferchichi A. Assessment of genetic diversity among some southern Tunisian fig (Ficus carica L.) cultivars based on morphological and chemical characters. Acta Bot Gallica. 2010;157:171–82.
Podgornik M, Vul I, Vrhovnik I, Mavsar DB. A survey and morphological evaluation of fig (Ficus carica L.) genetic resources from Slovenia. Sci Hortic. 2010;125:380–9.
Giraldo E, Lόpez-Corrales M, Hormaza JI. Selection of the most discriminating morphological qualitative variables for characterization of fig germplasm. J Amer Soc Hort Sci. 2010;135:240–9.
Aljane F, Nahdi S, Essid A. Genetic diversity of some accessions of Tunisian fig tree (Ficus carica L.) based in morphological and chemical traits. J Nat Prod Plant Resour. 2012;2:350–9.
Almajali D, Abdel-Ghani AH, Migdadi H. Evaluation of genetic diversity among Jordanian fig germplasm accessions by morphological traits and ISSR markers. Sci Hortic. 2012;147:8–19.
Cabrita LF, Aksoy U, Hepaksoy S, Leitao JM. Suitability of isozyme, RAPD and AFLP markers to assess genetic differences and relatedness among fig (Ficus carica L.) clones. Sci Hortic. 2001;87:261–73.
Chessa I, Nieddu G. Analysis of diversity in the fruit tree genetic resources from a Mediterranean island. Genet Resour Crop Evol. 2005;52:267–76.
Khadari B, Lashermes P, Kjellberg F. RAPD fingerprints for identification and genetic characterization of fig (Ficus carica L.) genotypes. J Genet Breed. 1995;49:77–86.
Galderisi U, Cipollaro M, Di Bernardo G, De Masi L, Galano G, Cascino A. Identification of the edible fig “Bianco del Cilento” by random amplified polymorphic DNA analysis. HortSci. 1999;34:1263–5.
Salhi-Hannachi A, Chatti K, Mars M, Marrakchi M. Trifi M (2005) Comparative analysis of genetic diversity in two Tunisian collections of fig cultivars based on random amplified polymorphic DNA and inter simple sequence repeats fingerprints. Genet Resour Crop Evol. 2005;52:563–73.
Akbulut M, Ercisli S, Karlidag H. RAPD-based study of genetic variation and relationships among wild fig genotypes in Turkey. Genet Mol Res. 2009;8:1109–15.
Ikegami H, Nogata H, Hirashima K, Awamura M, Nakahara T. Analysis of genetic diversity among European and Asian fig varieties (Ficus carica L.) using ISSR, RAPD, and SSR markers. Genet Resour Crop Evol. 2009;56:201–9.
Chatti K, Baraket G, Ben Abdelkrim A, Saddoud O, Mars M, Trifi M, et al. Development of molecular tools for characterization and genetic diversity analysis in Tunisian fig (Ficus carica) cultivars. Biochem Genet. 2010;48:789–806.
Salhi-Hannachi A, Trifi M, Zehdi S, Hadfi J, Mars M, Rhouma A, et al. Inter simple sequence repeat fingerprints to assess genetic diversity in Tunisian fig (Ficus carica L.). Genet Resour Crop Evol. 2004;51:269–75.
Baraket G, Chatti K, Saddoud O, Mars M, Marrakchi M, Trifi M, et al. Genetic analysis of Tunisian fig (Ficus carica L.) cultivars using amplified fragment length polymorphism (AFLP) markers. Sci Hortic. 2009;120:487–92.
Baraket G, Chatti K, Saddoud O, Ben Abdelkarim A, Mars M, Trifi M, et al. Comparative assessment of SSR and AFLP markers for evaluation of genetic diversity and conservation of fig, Ficus carica L., genetic resources in Tunisia. Plant Molec Biol Rep. 2011;29:171–84.
Khadari B, Hochu I, Santoni S, Kjellberg F. Identification and characterisation of microsatellite loci in the common fig (Ficus carica L.) and representative species of genus Ficus. Mol Ecol Notes. 2001;1:191–3.
Saddoud O. Salhi- Hannachi A, Chatti K, Mars M Rhouma A, Marrakchi M, Trifi M. Tunisian fig (Ficus carica L.) genetic diversity and cultivar characterization using microsatellite markers. Fruits. 2005;60:143–53.
Giraldo E, Viruel MA, Lopez-Corrales M, Hormaza JI. Characterisation and cross-species transferability of microsatellites in the common fig (Ficus carica L.). J Hort Sci Biotech. 2005;80:217–24.
Bandelj D, Javornik B, Jakse J. Development of microsatellite markers in the common fig, Ficus carica L. Mol Ecol Notes. 2007;7:1311–4.
Giraldo E, Lopez-Corrales M, Hormaza JI. Optimization of the management of an ex-situ germplasm bank in common fig with SSRs. J Amer Soc Hort Sci. 2008;133:69–77.
Achtak H, Ater M, Oukabli A, Santoni S, Kjellberg F, Khadari B. Traditional agroecosystems as conservatories and incubators of cultivated plant varietal diversity: the case of fig (Ficus carica L.) in Morocco. BMC Plant Biol. 2010;10:28.
Aradhya MK, Stover E, Velasco D, Koehmstedt A. Genetic structure and differentiation in cultivated fig (Ficus carica L.). Genetica. 2010;138:681–94.
Do Val ADB, Souza CS, Ferreira EA, Salgado SML, Pasqual M, Cançado GMA. Evaluation of genetic diversity in fig accessions by using microsatellite markers. Genet Molec Res. 2013;12:1383–91.
Perez-Jiménez M, López B, Dorado GM, Pujadas-Salvá A, Guzmán G, Hernandez P. Analysis of genetic diversity of southern Spain fig tree (Ficus carica L.) and reference materials as a tool for breeding and conservation. Hereditas. 2012;149:108–13.
Balas FC, Osuna MD, Domínguez G, Pérez-Gragera F, López-Corrales M. Ex situ conservation of underutilised fruit tree species: establishment of a core collection for Ficus carica L. using microsatellite markers (SSRs). Tree Genet Genom. 2014;10:703–10.
Abdelkrim AB, Baraket G, Essalouh L, Achtak H, Khadari B, Salhi-Hannachi A. Use of morphological traits and microsatellite markers to characterize the Tunisian cultivated and wild figs (Ficus carica L.). Biochem Syst Ecol. 2015;59:209–19.
Yeh FC, Young RC, Timothy B, Boyle TBJ, Ye ZH, Mao JX. Popgene, the user-friendly shareware for population genetic analysis. Canada: Molecular Biology and Biotechnology Center, University of Alberta; 1997.
Wagner H, Sefc K (1999) Identity 1.0 Centre for Applied Genetics, Univ. Agric. Sci., Vienna (http://www.uni-graz.at/~sefck/manual.pdf); last accessed May 2015.
Nei M, Li WH. Mathematical model for studying genetic variation in terms of restriction endonucleases. Proc Natl Acad Sci U S A. 1979;76:5269–73.
Van de Peer Y, De Watchter R. TREECON for Windows: a software package for the construction and drawing of evolutionary trees for the Microsoft Windows environment. Comput Appl Biosci. 1994;10:569–70.
Callen DF, Thompson AD, Shen Y, Phillips HA, Richards RI, Mulley JC, et al. Incidence and origin of “null” alleles in the (AC)n microsatellite markers. Am J Human Genet. 1993;52:922–7.
Financial support for this work was provided by the Tunisian Education Ministry, the University of Monastir and the Spanish Ministry of Economy and Innovation (Project Grant AGL2013-43732-R). We acknowledge support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI). We gratefully acknowledge the farmers in ‘CFPA’ ‘El Gordhab’, Tatouine and the Institute of Arid Region, Médenine Tunisia for their cooperation and Y. Verdún with help with the molecular analyses.
The authors declare that they have no competing interests.
AE, FA, AF and JIH designed research; AE performed research; AE, FA, AF and JIH analyzed data and wrote the paper. All authors read and approved the final manuscript.