Assessment of the genetic diversity of the Tunisian citrus rootstock germplasm
© Snoussi et al; licensee BioMed Central Ltd. 2012
Received: 6 September 2011
Accepted: 19 March 2012
Published: 19 March 2012
Citrus represents a substantial income for farmers in the Mediterranean Basin. However, the Mediterranean citrus industry faces increasing biotic and abiotic constraints. Therefore the breeding and selection of new rootstocks are now of the utmost importance. In Tunisia, in addition to sour orange, the most widespread traditional rootstock of the Mediterranean area, other citrus rootstocks and well adapted to local environmental conditions, are traditionally used and should be important genetic resources for breeding. To characterize the diversity of Tunisian citrus rootstocks, two hundred and one local accessions belonging to four facultative apomictic species (C. aurantium, sour orange; C. sinensis, orange; C. limon, lemon; and C. aurantifolia, lime) were collected and genotyped using 20 nuclear SSR markers and four indel mitochondrial markers. Multi-locus genotypes (MLGs) were compared to references from French and Spanish collections.
The differentiation of the four varietal groups was well-marked. The groups displayed a relatively high allelic diversity, primarily due to very high heterozygosity. Sixteen distinct MLGs were identified. Ten of these were noted in sour oranges. However, the majority of the analysed sour orange accessions corresponded with only two MLGs, differentiated by a single allele, likely due to a mutation. The most frequent MLG is shared with the reference sour oranges. No polymorphism was found within the sweet orange group. Two MLGs, differentiated by a single locus, were noted in lemon. The predominant MLG was shared with the reference lemons. Limes were represented by three genotypes. Two corresponded to the 'Mexican lime' and 'limonette de Marrakech' references. The MLG of 'Chiiri' lime was unique.
The Tunisian citrus rootstock genetic diversity is predominantly due to high heterozygosity and differentiation between the four varietal groups. The phenotypic diversity within the varietal groups has resulted from multiple introductions, somatic mutations and rare sexual recombination events. Finally, this diversity study enabled the identification of a core sample of accessions for further physiological and agronomical evaluations. These core accessions will be integrated into citrus rootstock breeding programs for the Mediterranean Basin.
Worldwide production of citrus in 2009 reached greater than 120 million tons , making citrus the leading cultivated tree crop in the world. Oranges represent the majority of citrus production (54% in 2009) with over 67 million tons in 2009 . The other significant cultivated citrus are mandarins, lemons and grapefruits. The citrus production of the Mediterranean Basin is second only to Brazil. Cultivars are vegetatively propagated by bud-grafting onto seedling rootstocks. This ensures tree uniformity, early tree production, and tolerance to pathogens including Phytophthora sp., parasitic nematodes and viruses. Rootstocks are also significant in the adaptation of the tree to several abiotic constraints affecting the Mediterranean citriculture. These include water resource scarcity and soil salinity. Citrus rootstocks are generally apomictic and seed-propagated. Therefore, both scions and rootstocks are clonally propagated. Sour orange (C. aurantium L.), one of the most important citrus rootstocks in the world, is still the predominant rootstock on the southern edge of the Mediterranean Basin. Sour orange is tolerant to limestone, alkalinity and salinity and is resistant to Phytophthora sp. Furthermore, sour orange is widely compatible with scion varieties and confers good fruit quality. Unfortunately, combining sour orange with commercial citrus varieties results in trees that are sensitive to Citrus Tristeza Virus (CTV), causing the rapid decline and death of the trees grafted on sour orange rootstocks. As CTV spreads throughout the Mediterranean basin, citrus production on sour orange rootstock will soon be uneconomic. On the other hand, rootstocks selected for their resistance to CTV are not well adapted to other local constraints . Thus, there is an urgent need to diversify and select new citrus rootstocks exhibiting CTV resistance and adaptation to the regional abiotic stresses.
The most widely accepted taxonomic systems for Citrus are those of Swingle and Reece  and Tanaka , who recognized 16 and 162 species, respectively. Later, phylogenetic analysis by Scora  and Barrett and Rhodes  indicated only three true species within the cultivated citrus, i.e., C. medica L. (citron), C. reticulata Blanco (mandarin) and C. maxima (Burm.) Merr. (pummelo). Recent molecular studies have confirmed the central role of these three taxa [7–15] and concluded in favour of a fourth additional ancestral taxon, C. micrantha Wester [7, 8]. The secondary species, C. aurantium L. (sour orange), C. limon (L.) Burm.f. (lemon), C. aurantifolia (Christm.) Swing. (lime), C. sinensis (L.) Osb. (sweet orange) and C. paradisi Macfad. (grapefruit), were derived from hybridization amongst the true species [7, 9, 10, 16]. Despite significant phenotypical differentiation, all Citrus species and several related genera are sexually compatible [17, 18] and can be considered part of the same biological species. Most citrus species are characterized by facultative apomixis resulting from adventitious nucellar embryogenesis . This has deeply affected the intra- and inter-specific evolution of cultivated citrus. The total sexual compatibility within the genus permitted the occurrence of numerous inter-specific hybrids, which evolved through vegetative propagation due to their facultative apomixis. This led taxonomists to consider clonally propagated families of inter-specific origin as new species. As a result, Citrus taxonomy and systematics remain controversial . Citrus was domesticated in Southeast Asia, notably East India, North Burma and Southwest China, and then spread to other continents [5, 20]. Citrus introduction to Tunisia likely occurred during the 10th century, and the citrus industry was established in the beginning of the 20th century. Since 1934, export trade has undergone great expansion, and producers began producing the 'Maltaise demi-sanguine' orange. Today, citrus orchards occupy approximately 20,400 hectares, with citrus production oscillating between 210,000 and 300,000 tons during the last decade. Sour orange is still the most widely used traditional rootstock of Tunisian citriculture. Despite clonal propagation, the morphological variation exhibited by local accessions of sour orange is remarkable. Additional rootstock varieties and non-grafted citrus trees are still used, particularly in the oasis areas, and these trees are related to C. sinensis, C. limon and C. aurantifolia. The traditional Tunisian rootstocks are well adapted to the adverse regional soils and climatic conditions. Consequently, these rootstocks constitute a very important germplasm to be employed as an abiotic stress tolerance source for future rootstock breeding projects. Therefore, the activities of investigation, collection, preservation and characterization of this germplasm are a priority. The assessment of genetic diversity within the rootstock germplasm is a prerequisite to the optimization of its management at the national and regional level.
Simple sequence repeat (SSR) markers have gained considerable importance in plant genetics due to their many desirable genetic attributes, including high polymorphism, wide genomic distribution, co-dominant inheritance and reproducibility allowing networking activities, multi-allelic nature and chromosome-specific location. In citrus, nuclear SSR development from genomic libraries [21–24], ESTs [25, 26] and BACend sequences  has been important during the last decade. These markers have proved to be very useful for the evaluation of genetic diversity in citrus [13, 15] and for the analysis of the sexual or apomictic origin of plant seedlings [27–29].
The main objective of this study was to assess the genetic diversity of Tunisian citrus rootstock germplasm. Specifically, the following questions were addressed: (i) What are the extent and the structure of the Tunisian rootstocks genetic diversity? (ii) What are the origin and the extent of the genetic diversity within facultative apomictic species? (iii) How is the Tunisian germplasm diversity related to reference genotypes of the same varietal groups of certified Citrus collections? For this purpose, 201 individual from different regions of Tunisia were analysed, and a set of reference cultivars from Spanish and French germplasm banks were genotyped for 20 nuclear SSR markers and four Indels mitochondrial markers.
Genetic diversity of Tunisian germplasm
Allelic diversity and its organization
Genetic diversity of the Tunisian germplasm for 20 SSR markers
Genetic diversity of the four varietal groups prospected in Tunisia
Among the 201 analysed accessions, only 16 distinct MLGs were obtained for the 20 loci studied yielding a G/N ratio of 0.079. Considering each varietal group (Table 2), the number of MLGs identified from 169 sour orange, 8 lime, 12 lemon and 12 sweet orange accessions were 10, 3, 2 and 1, respectively. It is notable that the majority (160) of the sour orange accessions corresponded to two MLGs, 'sour orange A' and 'sour orange B', representing 128 and 32 accessions, respectively. These two genotypes differed only in a single allele at locus mCrCIR01C06. Sour orange C, representing two accessions, differed from sour orange A only by its homozygosity at one locus (mCrCIR02G12). Interestingly, five of the other sour orange MLGs, each corresponding to one accession (sour orange D, E, F, J, H), displayed approximately half of the heterozygosity of the primary sour orange MLG without additional alleles. These five accessions may have arisen from zygotic seedlings following self-pollination. The two last sour orange MLGs (G and I) displayed intermediate proportions of heterozygous loci (78.6% and 65%, respectively) and a high number of 'new' alleles (11 and 12 alleles, respectively) not shared with the sour orange groups A, B and C.
The four 'Chiiri' limes displayed the same genotype for the 20 analysed loci. The two 'Sweet' limes shared the second lime MLG, the third lime MLG being represented by 'Arbi' and 'Beldi lemon'. Eleven of the 12 lemon accessions were identical, and the last one differed only at a single locus.
Comparison of local germplasm with reference accessions
Different loci and alleles distinguishing local accessions from the references
Loci different from the corresponding reference MLG
Considered MLG reference
Other reference MLG sharing these alleles*
Sour Orange A (id ref)
Sour Orange B
Sour Orange C
Sour Orange D
Sour Orange E
Sour Orange F
Sour Orange J
Sour Orange H
Sour Orange G
mandarin (10), lime Rangpur (4), Orange (3), Rough lemon (4)...
Sour Orange I
Sweet orange (12), mandarin (5), pummelo (3)...
Lemon A (id ref)
orange(1), sweet lime(1), mandarin (1)
citron (8), lemon (6)
hybridization with citron
The clustering was generally consistent with varietal group classification. One cluster consisted of nine Tunisian sour orange MLGs and the three reference sour oranges. These three reference accessions and the most representative Tunisian MLG 'Sour orange A' were identical. The sour orange MLG G exhibited 14 loci which differed from the reference. Eleven alleles were not observed in the sour orange reference accessions. Ten of these alleles (91%) were found to be shared with 'Willow leaf' mandarin. Sour orange MLG I showed 12 alleles that were not observed in the sour orange references, although these 12 alleles are present in sweet oranges. These results suggest that G and I are hybrids of sour orange and mandarin, and of sour orange and sweet orange, respectively.
The sweet orange MLG from Tunisia and all sweet orange reference cultivars were identical. Furthermore, the MLG clustered with the 'Willow leaf' mandarin included as a C. deliciosa reference and with Pummelo 'Pink' included as the C. maxima reference.
All acid citrus (limes and lemons) clustered together. The 'lemon A' MLG, representing 11 of the 12 Tunisian lemon accessions, was identical to the two lemon references. The 'Lemon Beldi'/lime 'Arbi' MLG was identical to the 'Mexican' lime reference. Lime 'Chiiri' was closely related to these genotypes (bootstrap value of 50%) but differed from all the reference genotypes included in this study. Indeed, 'Chiiri' lime displayed 14 loci (50% of them were homozygous), which varied from the Mexican lime reference. This included nine alleles not observed in the Mexican lime, with eight (88.9%) of these being common to Corsican citron (Table 3). The 'sweet lime' MLG was identical to the 'limonette de Marrakech' genotype. When considering the results obtained from the Tunisian MLG and references (Table 1), the number of alleles for the 20 SSR locus increased from 120 for the Tunisian germplasm (six alleles per locus) to 162 (8.1 alleles per locus). All 120 alleles found in the Tunisian germplasm were also encountered in the reference set. However, 11 among the 16 local genotypes were not represented in the reference set.
Inter-specific differentiation is the main structuring factor of Tunisian germplasm diversity
Population genetic parameters and cluster analysis indicated a high genetic structuring of the citrus Tunisian germplasm. This was also characterised by a significant excess of heterozygosity, both in the whole sample and in each varietal group (sour orange, sweet orange, lime and lemon). These results evidenced the important differentiation between the four varietal group analysed (sweet orange, sour orange, lime and lemons) and the limitation of gene flow between these groups.
The set of 20 SSR markers used in this study was sufficient to obtain genetic aggregations, enabling the clear identification of varietal groups of citrus and producing results consistent with previous studies [7, 9, 15, 16]. The clustering of the Tunisian MLGs and reference cultivars confirmed the importance of inter-specific differentiation as a major element of the Tunisian germplasm variability. Indeed, the Tunisian accessions were split within three clusters. The first cluster included the majority of the Tunisian sour orange MLGs with the sour orange reference. All sweet orange accessions, including the Tunisian MLGs and reference, were identical and grouped in the second cluster, along with 'Willow Leaf' mandarin and 'Pink' pummelo. This grouping is in accordance with the origin of sweet orange, believed to have originated, as sour orange, from hybridizations between mandarin and pummelo gene pools . When analysing mitochondrial markers, sweet orange and sour orange were found to be identical to the pummelo mitotype (Figure 2). This supports the previous hypothesis that pummelo acted as the female parent of sweet and sour oranges [7, 8]. A third cluster grouped limes and lemons with citron and the other acid citrus species. All Tunisian lime and lemon were identical to the reference samples of the Spanish and French collection, with the exception of 'Chiiri' lime and the lemon MLG B. The currently accepted hypothesis is that lemon and lime are the result of secondary diversification by interspecific hybridization involving citron as one of the parents [7, 15]. Nicolosi et al.  and Gulsen and Roose  proposed that lemon arose via hybridization between sour orange and citron. The mitotype analysis presented here confirms the hypothesis that sour orange was the maternal parent [7, 8]. The numerous and complex origin of the lime cultivars will be discussed further with the source of intra-specific diversity in the Tunisian germplasm. All citrus species are sexually compatible [17, 18]. However, vegetative horticultural propagation methods and facultative apomixis of the considered genotypes have undoubtedly prevented genetic mixing in Tunisia. Therefore, the multiple introductions (from different highly heterozygous secondary species) and the limitation of sexual recombination underlie a genetic diversity organization mainly based on inter-specific differentiation and very high heterozygosity.
Intra-specific diversity and the source of variation in Tunisian germplasm
The observed heterozygosity within each varietal group is very high. This heterozygosity is associated with low inter-cultivar differentiation within species (low MLG/accession ratios). These results are consistent with the generally accepted inter-specific origin of these varietal groups [6, 7, 9, 10, 14, 15] and the predominant role of mutation or epigenetic variation in their inter-specific differentiation [10, 31].
Genotypic analysis of the 169 Tunisian sour orange accessions, at 20 microsatellite loci, revealed ten unique genotypes with important redundancies (G/N = 0.059). These redundancies were expected, as standard vegetatively-propagated rootstock varieties are distributed throughout the country. Thus, 162 accessions were found to correspond to three highly heterozygous, very closely clustered MLGs that differed only in a single allele. The most represented MLG (128 accessions) was identical to the three reference genotypes originating from other countries in the Mediterranean basin (Morocco, Spain, and Corsica). Within this group, different accessions were collected and named according to peculiar morphological traits. For example, the accessions 'sour orange Chiiri', 'sour orange with flattened fruits', 'sour orange with very small fruits and leaves', and sour oranges 'Arbi' and 'Souri' shared identical MLG A. It is suggested that the phenotypical variants observed between and within these three MLGs may be somatic mutants of the ancestral sour orange. The seven other sour orange MLGs corresponded to single accessions. All of them shared the same mitotype as the reference sour orange. Five of these MLGs exhibit reduced heterozygosity compared with the reference sour orange, as well as an absence of additional alleles when compared with the three predominant MLGs. It is highly probable that these MLGs are zygotic plants resulting from self-pollination within the predominant sour orange groups. The two remaining sour orange MLGs are likely sour orange × mandarin and sour orange × sweet orange hybrids. Six of the seven sour orange MLGs of zygotic origin were collected in Sbikha (Kairouan) from an orchard of mother trees used for seed propagation. All sour orange trees in this orchard were grown from seed. The other 22 plants analysed from this orchard genetically conformed to the predominant sour orange MLG. Thus, a rate of 21.4% zygotic plants can be estimated in this population. Facultative apomixis of citrus results from adventitious nucellar polyembryony in seeds also containing one zygotic embryo. The frequency of apomictic plants among seedlings is determined by the competition between nucellar embryos and the zygotic embryo during seed development and seed germination , as well as genotype and environment [19, 32]. Frost and Soost  reported zygotic rates of 15% for sour orange, while Moore and Castle  did not find any zygotics in sour orange seedlings. Observations from the Sbikha seed park affirm the need to control the genetic conformity of seedling-produced mother trees in a propagation scheme. It also underscores the biological potential of diversification through sexual recombination, even in areas where only facultative apomictic genotypes are present. However, among all the included accessions, these zygotic plants have been observed as unique accessions and in a very specific context. These results imply that nursery workers and growers exert highly efficient counter-selection against these zygotic off-types by visual observation.
The sweet orange accessions collected in Tunisia, although identified by phenotypic characters of the fruit (sweet, acid, red-coloured) or reported for their resistance to salinity (cv Meski), were identical for all analysed markers and highly heterozygous. Moreover, the sweet orange accessions also exhibited the same molecular profile as the five reference cultivars from the IVIA collection. This narrow genetic basis of the sweet orange accessions has been recognized in previous studies [13, 31, 35, 36]. As for sour orange, most authors believe that sweet orange had interspecific origin (between pummelo and mandarin gene pools), and the inter-varietal diversity within this species is attributed to somatic mutations [10, 15, 31, 37]. These mutations may alter horticultural characters, mostly fruit traits , and thus have been selected by man.
Lemons were highly heterozygous for the SSR loci analysed. The Tunisian accessions displayed two MLGs. The main one (A) was identical to 'Eureka' and 'Lisbon' lemon reference accessions, and the second displayed only one different SSR allele, along with a marked mammiform apex of a pyriform fruit. The latter type most probably resulted from mutations of the reference lemons.
Three different MLGs were identified for lime in Tunisia. Although easily identified by the farmers due to their phenotypic differences, the Tunisian 'Arbi' and 'Beldi' limes exhibited the same genotype as the 'Mexican lime' reference accession, of which they may be phenotypic variants. These two Tunisian limes are therefore part of C. aurantifolia Swing. (Tanaka classification). The sweet lime accessions did not differ from 'limonette de Marrakech', a lime from Morocco and are classified in C. limetta Risso (Tanaka classification). The very close phenotypic relationship of the Tunisian sweet lime and 'limonette de Marrakech' was mentioned by Hodgson . For this author, the only difference was the acidless flesh of the Tunisian lime while acidity was present in the 'limonette de Marrakech'. From our results, the Tunisian sweet lime is very likely an acidless mutant of an acid form such as 'limonette de Marrakech' as proposed by Hodgson. The 'Chiiri' lime MLG was not observed among the reference genotypes and clustered between the 'Mexican lime' genotype and the citron. It displayed a low heterozygosity level compared to the other citrus groups. A locus by locus comparison with 'Mexican lime' revealed a higher proportion of homozygous loci (50% versus 20%) and seven different alleles, six of which were shared with the reference citron. It is proposed that the Tunisian 'Chiiri' lime may have resulted from a cross between the 'Mexican' lime or one of its 'mutants' and citron. A recent study based on mitochondrial markers  suggested that limes have three different maternal origins: the C. maxima mitotype is shared with 'Marrakech limonette' and 'Palestine sweet lime', the C. micrantha mitotype is shared with 'Mexican' lime, and the acid mandarin (C. reshni) mitotype is shared with 'Rangpur' lime. In our analysis, the mitochondrial data confirmed the identities found at the nuclear level between limonette de Marrakech and Tunisian sweet limes as well as between Mexican lime and Beldi' lemon and 'Arbi' lime. 'Chiiri' lime also shared the same cytoplasm with 'Mexican' lime, which is consistent with Mexican lime being its maternal parent.
The analysis of Tunisian germplasm using microsatellites confirmed that the intra-specific diversification of the four apomictic species has resulted from three processes: (1) multiple introductions from diversified material, (2) mutation of local material, and (3) sexual recombination. The contributions of each of these processes varied depending on the species considered. For sweet oranges and lemons displaying very high heterozygosity and no or very little intra-group genetic diversity, sexual recombination can be discarded. For sour orange, sexual recombination was at the origin of seven accessions among 169 analysed (4%). For lime, the identities of two Tunisian MLGs with references from other countries points to a minimum of two introductions. The 'Chiiri' lime potentially arose from sexual events in Tunisia. However, it is difficult, without traceability of plant origin, to distinguish between multiple introductions of pre-differentiated materials or a de novo variation in Tunisia. These results show that, although rare, sexual reproduction occurs, and the progeny of self or inter-specific crossing may be selected by growers. However, it appears that mutation or epigenetic variation were the major factor of diversification within the different varietal groups, as previously underscored for other apomictic species [39, 40]. The contribution of somatic mutations to the evolution of other vegetatively propagated crops, such as grapes , olives , yams  and cassava , has also been demonstrated. Human selection of new phenotypes and further clonal propagation are also key factors underlying the relatively high inter-varietal morphological polymorphism [6, 10].
In Tunisia, sour orange is found in nearly all citrus growing areas and is well adapted to environmental conditions. Lemon and lime represent other species traditionally used as self-rooted trees or as rootstocks, mainly in southern Tunisia, where rainfall is only 120 to 200 mm per year and most irrigation water is brackish. In this region, irrigation water can contain up to 4-7 g/l of salt, especially during the summer. These rootstocks are well-adapted to local conditions and bear full flavoured fruits. Farmers consider the 'Chiiri' lime variety to be the best adapted citrus rootstock to local bioclimatic conditions (salinity and drought) and report that it has better water uptake ability in comparison to sour orange. This study revealed that Tunisian citrus rootstock diversity has resulted from multiple introductions, mutations and residual sexual recombination of apomictic genotypes, as well as the impact of human selection. Original germplasm has been identified, particularly mutants and zygotic plants (mostly from selfing) of sour oranges, phenotypical variants of Mexican limes and a probable hybrid between a Mexican lime type and citron ('Chiiri' lime). This variability may be exploited either directly by nursery men and growers or by breeders to combine resistance to biotic and abiotic stresses with improvements in horticultural behaviour through conventional breeding or biotechnological approaches, such as somatic hybridization [2, 45, 46]. This study indicates the value of an accurate physiological and agronomic evaluation of the different MLGs identified with particular emphasis for tolerance to biotic and abiotic stresses and agronomic behaviour under unfavourable environmental conditions.
Materials and methods
Primers of the 20 SSR loci used in the analysis
Primer sequence Forward
Primer sequence Reverse
Four mitochondrial indel markers (rrn5/rrn18-1, nad2/4-3, nad5/2-1, and nad7/1-2), developed by Froelicher et al. , were also used to analyse the maternal phylogeny of the Tunisian germplasm. These markers allowed a clear differentiation among the four basic taxa of cultivated citrus. Moreover they differentiated 3 mitotypes for limes and were therefore well adapted to the objectives of the study.
Molecular marker analysis
Initially, the 201 collected accessions were genotyped using 20 polymorphic SSR loci at CIRAD (Montpellier, France). Primers were radio-labelled with 33P prior to PCR amplification. The PCR reaction for SSR analysis was performed in a total volume of 20 μl containing 25 ng genomic DNA, 0.2 μM of each primer, 200 μM of each dNTP, 2 μl of 10X PCR buffer (100 mM Tris-HCl pH 8.3, 500 mM KCl, 15 mM MgCl2, 0.01% glycerol), and 0.75 units of Taq DNA polymerase. PCR reactions were performed using a MJ research thermal cycler model PTC-100™ and the following parameters: an initial denaturation at 94°C for 5 min, 35 amplification cycles of 30 s at 94°C, 1 min at 50-55°C (depending on the primer pair), 30 s at 72°C, and a final elongation at 72°C for 4 min. PCR products were visualized by running 6 μl of each sample on 5% polyacrylamide denaturing gels. Gels were dried and exposed to autoradiographic film. To minimize genotyping errors and for further comparison, all gels included a control individual (sour orange 'Granito') in four gel lanes and a ladder as a size standard. Unclear genotypes were systematically reanalysed.
Using the same 20 nuclear SSR markers and four Indels mitochondrial markers, the different multi-locus genotypes (MLGs) observed in the initial analysis were reanalysed along with the reference set of 23 accessions. These analyses were conducted using capillary electrophoresis (CEQ™ 8000 Genetic Analysis System) at IVIA (Valence, Spain) and repeated two times. The PCR reaction was performed using 2 ng/μl of template DNA, 0.2 mM of dNTPs, 0.2 μM of WellRED dye-labeled forward primer, 0.2 μM of reverse primer, 10X PCR buffer (Fermentas), 1.5 mM of MgCl2 and 0.5 U/μl of Taq DNA polymerase (Fermentas). The PCR conditions were as follows: 5 min at 94°C; 40 cycles of 30 s at 94°C, 30 s at 55°C (or 50°C depending on the primer) and 1 min at 72°C; and 4 min at 72°C.
Analysis of genetic data
For each primer, bands (PAGE) or peaks (capillary electrophoresis) were scored as allelic data. The allelic pattern of each accession was cross-checked using two readers.
where dij represents the dissimilarity between units i and j, L represents the number of loci, and ml represents the number of matching alleles between i and j for locus l. From the dissimilarity matrix, a Neighbor-Joining tree  was computed using the DARwin software version 5.0.158 (Dissimilarity Analysis and Representation for Windows, http://darwin.cirad.fr/darwin . Branch robustness was tested using 1000 bootstraps.
GENETIX software 4.05.2 http://www.genetix.univ-montp2.fr/genetix/genetix.htm  was used to estimate several parameters describing the genetic diversity and its organization within and between the four varietal groups encountered in Tunisia (sour orange, sweet orange, lemon and lime). The parameters estimated for each marker were the allele number (n), the mean number of alleles, allele frequencies, observed and expected mean heterozygosities (Ho and He, respectively). The number of distinct MLGs in each varietal group was identified from the whole set of markers.
This work was funded by the European INCO CIBEWU Project entitled Citrus rootstock breeding for efficient water and nutrient use, Contract no.: 015453, and benefited from the support of CIRAD, the Ministry of Agriculture and Agricultural Development Services in Tunisia. The authors thank all farmers and nursery gardeners for their collaboration.
- FAOSTAT. 2009, [http://faostat.fao.org/site/567/default.aspx]
- Dambier D, Benyahia H, Pensabene-Bellavia G, Aka Kacar Y, Froelicher Y, Belfalah Z, Beniken L, Handaji N, Printz B, Morillon R, Yesiloglu T, Navarro L, Ollitrault P: Somatic hybridization for Citrus rootstock breeding: an effective tool to solve some important issues of the Mediterranean citrus industry. Plant Cells report. 2011, 30: 883-900. 10.1007/s00299-010-1000-z.View ArticleGoogle Scholar
- Swingle WT, Reece PC: The botany of Citru and its wild relatives. The Citrus industry. volume 1. Edited by: Reuther W, Webber HJ, Batchelor DL. 1967, Berkeley, CA, USA: University of California Press, 190-430.Google Scholar
- Tanaka T: Fundamental discussion of Citru classification. Stud Citrol. 1977, 14: 1-6.Google Scholar
- Scora RW: On the history and origin of Citrus. Bull Torr Bot Club. 1975, 102: 369-375. 10.2307/2484763.View ArticleGoogle Scholar
- Barrett HC, Rhodes AM: A numerical taxonomic study of affinity relationships in cultivated Citru and its close relatives. Syst Bot. 1976, 1: 105-136. 10.2307/2418763.View ArticleGoogle Scholar
- Nicolosi E, Deng ZN, Gentile A, La Malfa S, Continella G, Tribulato E: Citrus phylogeny and genetic origin of important species as investigated by molecular markers. Theor Appl Genet. 2000, 100: 1155-1166. 10.1007/s001220051419.View ArticleGoogle Scholar
- Froelicher Y, Mouhaya W, Bassene JB, Costantino G, Kamiri M, Luro F, Morillon R, Ollitrault P: New universal mitochondrial PCR markers reveal new information on maternal citrus phylogeny. Tree Genet Genomes. 2011, 7 (1): 49-61. 10.1007/s11295-010-0314-x.View ArticleGoogle Scholar
- Federici CT, Fang DQ, Scora RW, Roose ML: Phylogenetic relationships within the genus Citru (Rutacea) and related genera as revealed by RFLP and RAPD analysis. Theor Appl Genet. 1998, 96: 812-822. 10.1007/s001220050807.View ArticleGoogle Scholar
- Ollitrault P, Jacquemond C, Dubois C, Luro F: Citrus. Genetic diversity of cultivated tropical plants. Edited by: Hamon P, Seguin M, Perrier X, Glaszmann J. 2003, Paris: Science Publishers and Cirad, 193-217.Google Scholar
- Pang XM, Hu CG, Deng XX: Phylogenetic relationships within Citru and its related genera as inferred from AFLP markers. Genet Resour Crop Evol. 2007, 54: 429-436. 10.1007/s10722-006-0005-5.View ArticleGoogle Scholar
- Yamamoto M, Kobayashi S, Nakamura Y, Yamada Y: Phylogenetic relationships of Citrus revealed by diversity of cytoplasmic genomes. Techniques on Gene Diagnosis and Breeding in fruit Trees. Edited by: Hayashi T, Omura M, Scott NS. 1993, Okitsu, Japan: Fruit Trees Research Station, 39-46.Google Scholar
- Luro F, Rist D, Ollitrault P: Evaluation of genetic relationships in Citrus Genus by means of Sequence Tagged Microsatellites. Proceedings of the International Symposium on Molecular Markers for Characterizing Genotypes and Identifying Cultivars in Horticulture February 2001; Montpellier, France. Acta Horticulturae (ISHS). 2001, 546: 237-242.Google Scholar
- Ollitrault F, Terol J, Pina JA, Navarro L, Talon M, Ollitrault P: Development of SSR markers from Citrus clementin (Rutaceae) BAC end sequences and interspecific transferability in Citrus. Am J Bot. 2010, 97 (11): 124-129. 10.3732/ajb.1000280.View ArticleGoogle Scholar
- Barkley N, Roose M, Krueger R, Frederici C: Assessing genetic diversity and population structure in a Citrus germplasm collection utilizing simple sequence repeat markers (SSRs). Theor Appl Genet. 2006, 112: 1519-1531. 10.1007/s00122-006-0255-9.View ArticlePubMedGoogle Scholar
- Herrero R, Asíns MJ, Pina JA, Carbonell EA, Navarro L: Genetic diversity in the orange subfamily Aurantioideae. II. Genetic relationships among genera and species. Theor Appl Genet. 1996, 93: 1327-1334. 10.1007/BF00223466.View ArticlePubMedGoogle Scholar
- Barrett HC: Intergeneric hybridization of citrus and other genera in citrus cultivar improvement. Proceedings of the second International Citrus Congress: 1977; Orlando, Florida. International society of Citriculture (ISC). 1977, 2: 586-589.Google Scholar
- Iwamasa M, Nito N, Ling JT: Intra- and intergeneric hybridization in the orange subfamily, Aurantioideae. Proceedings of the 6th International Citrus Congress: 1988; Tel Aviv, Israel. International society of Citriculture (ISC). 1988, Margraf Scientific Books, 1: 123-130.Google Scholar
- Kepiro JL, Roose ML: Nucellar embryony. Citrus Genetics, Breeding and Biotechnology. Edited by: Khan IA. 2007, UK, USA: CAB International publishers, 141-149.View ArticleGoogle Scholar
- Tanaka T: Species problem in Citrus (revisio aurantiacearum, IX). 1954, Ueno, Tokyo: Japanese Society for Promotion of ScienceGoogle Scholar
- Kijas JMH, Fowler JCS, Thomas MR: An evaluation of sequence-tagged microsatellite-site markers for genetic analysis within Citru and related species. Genome. 1995, 38: 349-355. 10.1139/g95-045.View ArticlePubMedGoogle Scholar
- Ahmad R, Struss D, Southwick MS: Development and characterization of microsatellite markers in Citru. J Am Soc Hortic Sci. 2003, 128 (4): 584-590.Google Scholar
- Novelli VM, Cristofani M, Souza AA, Machado MA: Development and characterization of polymorphic microsatellite markers for the sweet orange (Citrus sinensi L. Osbeck). Genet Mol Biol. 2006, 29 (1): 90-96. 10.1590/S1415-47572006000100018.View ArticleGoogle Scholar
- Froelicher Y, Dambier D, Bassene JB, Costantino G, Lotfy S, Didout C, Beaumont V, Brottier P, Risterucci AM, Luro F, Ollitrault P: Characterization of microsatellite markers in mandarin orange (Citrus reticulat Blanco). Mol Ecol Resour. 2008, 8: 119-122. 10.1111/j.1471-8286.2007.01893.x.View ArticlePubMedGoogle Scholar
- Chen C, Bowman KD, Choi YA, Dang PM, Rao MN, Huang S, Soneji JR, McCollum TG, Gmitter FG: EST-SSR genetic maps for Citrus sinensi and Poncirus trifoliat. Tree Genet Genomes. 2008, 4: 1-10.View ArticleGoogle Scholar
- Luro F, Costantino G, Terol J, Argout X, Allario T, Wincker P, Talon M, Ollitrault P, Morillon R: Transferability of the EST-SSRs developed on Nules clementine (Citrus clementin Hort ex Tan) to other Citru species and their effectiveness for genetic mapping. BMC Genomics. 2008, 9 (287): 1-13.Google Scholar
- Rao MN, Soneji JR, Chen C, Huang S, Gmitter FG: Characterization of zygotic and nucellar seedlings from sour orange-like citrus rootstock candidates using RAPD and EST-SSR markers. Tree Genet Genomes. 2008, 4 (1): 113-124.Google Scholar
- Ruiz C, Paz Breto M, Asíns MJ: A quick methodology to identify sexual seedlings in citrus breeding programs using SSR markers. Euphytica. 2000, 112 (1): 89-94. 10.1023/A:1003992719598.View ArticleGoogle Scholar
- Aleza P, Juarez J, Ollitrault P, Navarro L: Polyembryony in non-apomictic citrus genotypes. Ann Bot. 2010, 106 (4): 533-545. 10.1093/aob/mcq148.PubMed CentralView ArticlePubMedGoogle Scholar
- Gulsen O, Roose ML: Lemons: diversity and relationships with selected citrus genotypes as measured with nuclear genome markers. J Am Soc Hortic Sci. 2001, 126: 309-317.Google Scholar
- Herrero R, Asíns MJ, Carbonell EA, Navarro L: Genetic diversity in the orange subfamily Aurantioideae. I. Intraspecies and intragenus genetic variability. Theor Appl Genet. 1996, 92: 599-609. 10.1007/BF00224564.View ArticlePubMedGoogle Scholar
- Khan IA, Roose ML: Frequency and characteristics of nucellar and zygotic seedlings in three cultivars of trifoliate orange. J Am Soc Hortic Sci. 1988, 113: 105-110.Google Scholar
- Frost HB, Soost RK: Seed reproduction: development of gametes and embryos. The Citrus Industry. Volume 1. Edited by: Reuther W, Batchelor LD, Webber HJ. 1968, Berkeley, California: University of California Press, 290-324.Google Scholar
- Moore GA, Castle WS: Morphological and isozymic analysis of open-pollinated Citru rootstock populations. J Hered. 1988, 79: 59-63.Google Scholar
- Fang DQ, Roose ML: Identification of closely related citrus cultivars with inter-simple sequence repeat markers. Theor Appl Genet. 1997, 95: 408-417. 10.1007/s001220050577.View ArticleGoogle Scholar
- Targon MLPN, Machado MA, Coletta Filho HD, Cristofani M: Genetic polymorphism of sweet orange (Citrus sinensi L. Osbeck) varieties evaluated by random amplified polymorphic DNA. Acta Hortic. 2000, 535: 51-53.View ArticleGoogle Scholar
- Nicolosi E: Origin and taxonomy. Citrus Genetics, Breeding and Biotechnology. Edited by: Khan IA. 2007, UK, USA: CAB International publishers, 19-43.View ArticleGoogle Scholar
- Hodgson RW: Horticultural varieties of Citrus. The Citrus Industry. Volume 1. Edited by: Reuther W, Batchelor LD, Webber HJ. 1967, Berkeley, California: University of California Press, 431-591.Google Scholar
- HÖrandl E, Paun O: Patterns and sources of genetic diversity in apomictic plants: implications for evolutionary potentials. Apomixis: evolution, mechanisms and perspectives. Edited by: HÖrandl E, Grossniklaus U, Van Dijk PJ, Sharbel TF. 2007, Koenigstein. Germany, 169-194.Google Scholar
- Nybom H: Unique reproduction in dogroses (Ros sect. Canina) maintains successful and highly heterozygous genotypes. Apomixis: evolution, mechanisms and perspectives. Edited by: HÖrandl E, Grossniklaus U, Van Dijk PJ, Sharbel TF. 2007, Koenigstein, Germany: Regnum Vegetabile, 281-298.Google Scholar
- Crespan M: Evidence of the evolution of polymorphism of microsatellite markers in varieties of Vitis vinifer L. Theor Appl Genet. 2004, 108: 231-237. 10.1007/s00122-003-1419-5.View ArticlePubMedGoogle Scholar
- Cipriani G, Marrazzo MT, Marconi R, Cimato A: Microsatellie markers isolated in olive (Olea europae L.) are suitable for individual fingerprinting and reveal polymorphism within ancient cultivars. Theor Appl Genet. 2002, 104: 223-228. 10.1007/s001220100685.View ArticlePubMedGoogle Scholar
- Scarcelli N: Structure et dynamique de la diversité d'une plante cultivée à multiplication végétative: le cas des ignames au Bénin (Dioscorea sp.). PhD thesis. 2005, Montpellier II University, FranceGoogle Scholar
- Sardos J, McKey D, Duval MF, Malapa R, Noyer JL, Lebot V: Evolution of cassava (Manihot esculent Crantz) after recent introduction into a South Pacific Island system: the contribution of sex to the diversification of a clonally propagated crop. Genome. 2008, 51: 912-921. 10.1139/G08-080.View ArticlePubMedGoogle Scholar
- Ollitrault P, Dambier D, Froelicher Y, Carreel F, d'Hont A, Luro F, Bruyère S, Cabasson C, Lotfy S, Joumaa A, Vanel F, Maddi F, Treanton K, Grisoni M: Apport de l'hybridation somatique pour l'exploitation des ressources génétiques des agrumes. Cahiers Agricultures, Ressources génétiques. 2000, 9 (3): 223-236.Google Scholar
- Grosser JW, Gmitter FG: Protoplast fusion for production of tetraploïds and triploids: applications for scion and rootstock breeding in citrus. Plant Cell Tissue Organ Cult. 2010, 104 (3): 343-357.View ArticleGoogle Scholar
- Perrier X, Jacquemoud-Collet JP: DARwin software. 2006, [http://darwin.cirad.fr/darwin]Google Scholar
- Saitou N, Nei M: The Neighbor-Joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4 (4): 406-425.PubMedGoogle Scholar
- Perrier X, Flori A, Bonnot F: Data analysis methods. Genetic diversity of cultivated tropical plants. Edited by: Hamon P, Seguin M, Perrier X, Glaszmann JC. 1999, Montpellier: CIRAD, 43-76.Google Scholar
- Belkhir K, Borsa P, Chikhi L, Raufaste N, Bonhomme F: GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations. Laboratoire Génome, Populations, Interactions, CNRS UMR 5171, Université de Montpellier II. 1996, Montpellier (France), 2004Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.