Triploid and tetraploid citrus hybrids were obtained via 2x X 2x and 2x X 4x sexual hybridizations using diploid “Eureka Frost” and “Fino” lemon genotypes as female parents pollinated with diploid “Fortune” mandarin (C. clementina x C. tangerina) and C. ichangensis Swing and tetraploid C. macrophylla Wester. Flowers in pre-anthesis were emasculated, pollinated, and enclosed with a cloth bag. A total of 115 “Eureka Frost” lemon flowers were pollinated, including 55 with “Fortune” mandarin (named EuFor) and 60 with C. ichangensis (named EuIch), while 15 “Fino” lemon flowers were pollinated with tetraploid C. macrophylla (named FinMac). The detailed methods used for plant recovery via in vitro embryo rescue and ploidy level analysis via flow cytometry can be found in Aleza et al. (2010a; 2010b; 2012a; 2012b).

The male and female parents and 48 hybrids were genotyped using 34 molecular markers (16 Simple Sequence Repeats [SSRs] and 18 Single Nucleotide Polymorphisms [SNPs]) showing heterozygosity for the lemon genotypes and polymorphism with the male parents. These markers are distributed across all LGs of the clementine genetic map (Ollitrault et al., 2012a). Detailed information about the markers is provided in Table 1.

Genomic DNA was isolated using a Plant DNeasy kit from Qiagen Inc. (Valencia, CA, USA) following the manufacturer's protocol. PCR amplifications using 16 SSR markers were performed using a Thermocycler rep gradient S (Eppendorf®) in a 10 μL final volume containing 0.8 U of Taq DNA polymerase (Fermentas®), 2 ng/μL citrus DNA, 0.2 mM welled (Sigma®) dye-labeled forward primer, 0.2 mM non-dye-labeled reverse primer, 0.2 mM of each dNTP, 10 × PCR buffer, and 1.5 mM MgCl₂. The PCR protocol was as follows: denaturation at 94°C for 5 min followed by 40 cycles of 30 s at 94°C, 1 min at 50 or 55°C, and 45 s at 72°C; and a final elongation step of 4 min at 72°C. Capillary electrophoresis was carried out using a CEN™ 8,000 Genetic Analysis System (Beckman Coulter Inc.). The PCR products were initially denatured at 90°C for 2 min, injected at 2 kV for 30 s, and separated at 6 kV for 35 min. Alleles were sized based on a DNA size standard (400 bp). Genome Lab™ Gap v.10.0 genetic analysis software was used for data collection. Allele dosage was calculated using the MAC-PR (microsatellite DNA allele counting-peak ratio) method (Esselink et al., 2004), validated in citrus by Cuenca et al. (2011).

Triploid and tetraploid hybrids were also genotyped with 18 SNP markers using KASPar™ technology by LGC Genomics (www.lgcgenomics.com). The KASPar Genotyping System is a competitive, allele-specific dual Förster Resonance Energy Transfer (FRET)-based assay for SNP genotyping. Primers were directly designed by LGC Genomics Company based on the SNP locus-flanking sequence (~50 nt on each side of the SNP). SNP genotyping was performed using the KASPar technique. A detailed description of specific conditions and reagents can be found in Cuppen (2007). Identification of allele doses in heterozygous triploid and tetraploid hybrids was carried out based on the relative allele signals, as described by Cuenca et al. (2013a) and Aleza et al. (2015).

For triploid and tetraploid hybrids, the 2n gamete origin was determined by identifying the parent that passed double genetic information onto the hybrid. Markers with total differentiation between the parents (A₁A₁ x A₂A₂A₂A₂, A₁A₂ x A₃A₃A₃A₃, and A₁A₂ x A₃A₃A₄A₄ in 2x X 4x crosses) for tetraploids and (A₁A₁ x A₂A₂, A₁A₂ x A₃A₃, and A₁A₂ x A₃A₄ in 2x X 2x crosses) for triploids were the best allelic configurations, as described by Aleza et al. (2015) and Cuenca et al. (2015). Indeed, conclusive results can be obtained using only one marker, as was the case for FinMac hybridization using the JK-TAA41 SSR marker. However, for EuFor and EuIch hybridizations, more than one marker had to be analyzed to observe both alleles from the female parent at least once for each hybrid. The SSRs JK-TAA1, JK-TAA41, and MEST131 were used for EuFor hybridization, and JK-TAA1, JK-TAA41, MEST001, and Ci02B07 were used for EuIch.

Once the female origin of the diploid gamete was demonstrated, inference of the allelic configurations of the 2n gametes from hybrid genotyping was performed as described by Cuenca et al. (2011). In the case of FinMac tetraploid hybridization, for the A₁A₂ × A₃A₃A₃A₃ and A₁A₂ × A₃A₃A₄A₄ allelic configurations, the genotype of the unreduced gamete was deduced directly from observation of both A₁ and A₂ alleles in the tetraploid hybrids. However, when the male and female parents shared one allele (A₁A₂ × A₁A₁A₁A₁ and A₁A₂ × A₁A₁A₃A₃), for the tetraploid hybrids that inherited the common allele (A₁), inference of the unreduced female gamete structure was carried out based on the estimated allele dosage in the tetraploid hybrid.

In the case of triploid hybrids obtained from EuFor and EuIch hybridizations, for A₁A₂ x A₃A₃ and A₁A₂ x A₃A₄, the genotype of the 2n gamete was deducted directly from the triploid hybrid structure. When the male and female genitors shared one allele (A₁A₂ x A₂A₂ and A₁A₂ x A₂A₃), the 2n female gamete structure for the triploid hybrids with a common allele from the male genitor was inferred from the estimated allele dosage in the triploid hybrid.

For the EuFor and EuIch progenies, nine SSR and SNP molecular markers within 20 cM of the centromere (Aleza et al., 2015) located in all nine LGs of the clementine genetic map (Ollitrault et al., 2012a) were genotyped to determine the mechanism of 2n gamete formation for each population. The molecular markers used included MEST001, JK-CAC15, 3p11355960, mCrCIR07D06, 5p22687304, 6p7496245, mCrCIR03B07, 8p18684429, and Ci07C09 for EuFor hybridization and MEST001, JK-CAC15, 3p11355960, mCrCIR07D06, 5p22687304, CiC4356-06, mCrCIR03B07, 8p18684429, and 9p4699283 for EuIch hybridization. For FinMac hybridization, seven molecular markers distributed in seven LGs were used, including MEST001, JK-CAC15, CiC4240-04, LapXcF238, mCrCIR03B07, 8P16570424, and CIBE3966.

To distinguish between the SDR and FDR hypotheses, the maximum-likelihood method based on the LOD score test described by Cuenca et al. (2015) was employed. LODs >2 were considered to be significant for SDR, those < −2 were considered to be significant for FDR, and those between 2 and −2 were considered not to be significant. To compare the SDR hypothesis with the PRD hypothesis using LOD scores, we considered the minimum value of 66% of PHR as the theoretical value for the PRD hypothesis.

Additionally, a set of six SSR and SNP molecular markers distributed along LG1 were used to analyze PHR evolution, including SSR markers MEST001, mCrCIR06B05, CIBE6126, and MEST431 and SNP markers CiC5950-02 and CiC2110-02. Moreover, a complementary experiment was performed to differentiate between PMD and SDR mechanisms using 11 telomeric molecular markers in LG2 to LG9. These included SSR markers mCrCIR03C08, JK-TAA41, MEST488 and mCrCIR03G05 and SNP markers CiC4876-07, CiC3674-02, CiC5842-02, CiC1459-02, CiC3712-01, 3p165889, and 8p2427684. The marker positions and distances to the centromeres of each LG are shown in Table 1.

Taking into account the centromere position, three-point linkage mapping was performed to estimate chiasma interference for each chromosome arm of chromosome I. The centromere was used as the first point, and two markers were selected on each arm (MEST001 and MEST431 on one arm and mCrCIR06B05 and CIBE6126 on the other arm). The chromosome interference coefficient (IC) is defined as follows (as per Griffiths et al., 1996): (1)IC=1-[rdrCM1·rM1M2] Where r_CM1 indicates the observed recombination rate (heterozygous to homozygous and vice versa) between the centromere and locus 1; r_M1M2, the observed recombination between locus1 and 2; and rd, the observed rate of double recombination between the centromere and locus 2.

For sexual hybridizations between “Eureka Frost” lemon as the female parent and “Fortune” mandarin and C. ichangensis as the male parents, the average fruit set was 45.5 and 36.7%, respectively (Table 2), yielding 250 and 464 seeds, respectively, from both hybridizations. We classified the seeds by size, since, according to Aleza et al. (2010a), seed size is highly correlated to ploidy level. While small seeds are expected to contain triploid embryos, tetraploids are generally observed in normal size seeds. Thus, we selected 45 and 40 small seeds from the EuFor and EuIch hybridizations, respectively, for plant regeneration by embryo rescue.

From the 45 small seeds obtained in the EuFor hybridization, 54 embryos were cultured in vitro, with an average of 1.2 embryos per seed, indicating a low rate of polyembryony in “Eureka” lemon. Of the 53 plantlets recovered, 32 were diploid and 21 triploid. All 40 small seeds recovered from the EuIch hybridization contained only a single embryo. Of the 35 plants regenerated, 21 were diploid and 14 were triploid. For the FinMac 2x X 4x sexual hybridization, the average fruit set was 53.3%, and 36 normal seeds were obtained according to the size classification of Aleza et al. (2012b). Of the 36 plants recovered, 23 were triploid and 13 were tetraploid (Table 2).

To determine which parent passed double genetic information onto the hybrids, we genotyped triploid hybrids recovered from the 2x X 2x hybridizations using markers that displayed total allelic differentiation between “Eureka Frost” lemon and the male parents, “Fortune” mandarin and C. ichangensis (Figure 1): SSR markers JK-TAA1, JK-TAA41, and MEST131 for the EuFor hybridization and SSR markers JK-TAA1, JK-TAA41, MEST001, and Ci02B07 for the EuIch hybridization. Genetic analysis enabled us to unequivocally identify the hybrid origins of all triploid plants, except for one plant from the EuFor sexual hybridization and four from the EuIch sexual hybridization, which were rejected since they could have originated from autopollination of the female parents. Genetic analysis showed that “Eureka Frost” lemon produced the 2n gametes for all triploid hybrids, as shown in Figure 1.

For the tetraploid hybrids, the JK-TAA41 SSR marker displayed total allelic differentiation between “Fino” lemon and tetraploid C. macrophylla, allowing us to conclude that all plants were hybrids and that “Fino” lemon produced the 2n gametes (Figure 1). Analysis of the genetic origins of the 23 triploid plants recovered from this 2x X 4x hybridization showed that, as expected, they were obtained from the union of a normal reduced haploid female gamete and a normal reduced diploid pollen gamete, as previously observed in other citrus species (Aleza et al., 2012a).

Lemon hybrids were obtained from 2n gametes at a frequency of 4.9% for “Eureka Frost” and 8.3% for “Fino”. Geraci et al. (1975) reported frequencies of 1 and 5% for triploid hybrids assumed to be obtained through unreduced gametes of “Lisbon” and “Eureka” lemons, respectively. Pérez-Tornero et al. (2012) obtained triploid hybrids at a frequency of 5.8 to 8.6% in hybridizations between diploid plants of “Verna” as the female parent and “Fino” as the male parent. In mandarins, greater differences between genotypes have been observed, ranging from <1% for clementines to over 22% for “Sukega” and “Ortanique” tangor (Esen and Soost, 1971; Wakana et al., 1982; Ollitrault et al., 2008; Aleza et al., 2010a; Xie et al., 2014).

The frequency of 2n gametes was shown to be genotype-dependent in citrus and in other herbaceous and woody plants such as Brassica, potato, and peach (Dermen, 1938; Mok et al., 1975; Ollitrault et al., 2008; Aleza et al., 2010a; Mason et al., 2011; Younis et al., 2014). This hypothesis is supported by the genetic improvement of unreduced gamete rates for Trifolium (frequencies increased from 0.04 to 47%) and Medicago sativa (from 9 to 78%) in only three generations of recurrent selection (Gallais, 2003).

In the current study, we observed a rate of 4.9% 2n gametes in the 2x X 2x hybridizations (EuFor and EuIch), whereas, in the 2x X 4x hybridization (FinMac), the percentage was higher (8.3%). These differences might be due to a genotypic effect of the parents, but are more likely due to the modification of the embryo/endosperm ploidy level ratio in interploid hybridizations. Esen and Soost (1971) reported that, in diploid plants, when an unreduced gamete is pollinated with normal reduced pollen, the embryo/endosperm ploidy level ratio (3/5) is less favorable for embryo development than that for normal diploid embryos (2/3), whereas the pollination of a 2n female gamete with diploid pollen in 2x X 4x sexual hybridizations provides the correct embryo/endosperm ploidy level ratio (4/6 = 2/3), leading to normal seed development. Therefore, 2x X 4x hybridization appears to be a more favorable situation for revealing unreduced gametes via the development of tetraploid embryos in normal seeds.

To determine the mechanism leading to unreduced gamete formation, we used nine unlinked molecular markers localized in the nine LGs for EuFor and EuIch and seven markers in seven different LGs for FinMac to perform a LOD score test for SDR/FDR and SDR/PRD probability ratios for all genotypes analyzed (Tables 3, 4, 5). The analysis of six markers covering LG1 and additional telomeric loci allowed us to distinguish between SDR and PMD when the inferred gametes were totally homozygous for the centromeric loci.

For the EuFor hybridization, 20 triploid hybrids were genotyped using nine centromeric loci found in all LGs. Ten of the inferred 2n gametes were totally homozygous for these markers. However, all displayed at least one heterozygous marker when six markers covering LG1 were analyzed, allowing the PMD hypothesis to be rejected for all inferred 2n gametes. For the SDR/FDR hypothesis test at the individual level, 19 inferred 2n gametes displayed LOD values >2 (ranging from 12.05 to 15.22; Table 3). For the same 19 gametes, the LOD values for SDR/PRD were also >2. Therefore, these 19 plants were considered to have originated from SDR. One plant obtained negative LODs of −4.52 and −6.86 for the SDR/FDR and SDR/PRD hypotheses, respectively, suggesting that this plant is of FDR or PRD origin. At the population level, the LOD values were 267.82 and 57.03 for the SDR/FDR and SDR/PRD hypotheses, respectively, revealing a high rate of SDR.

For EuIch hybridization, 10 triploid hybrids were genotyped with nine centromeric markers located on all LGs. Two inferred 2n gametes were totally homozygous for these markers, but at least one heterozygous locus was observed for each 2n gamete in the complementary analysis of PHR along the LG1, thus discarding the PMD hypothesis. At the individual level, eight plants displayed LOD values >2 for SDR/FDR (from 8.69 to 14.53), rejecting the FDR hypothesis (Table 4). Among them, seven displayed a LOD >2 for SDR/PRD (ranging from 2.13 to 3.86) and were considered to have arisen from SDR. The LOD value for the remaining 2n gamete was 0.55, suggesting that this 2n gamete had arisen from SDR rather than PRD, but, since this value is below our threshold, this result is not conclusive. Two plants produced negative LOD values (< −2) in both the SDR/FDR and SDR/PRD tests, suggesting that they originated by FDR or PRD. The population LODs were 80.21 and 2.77 for SDR/FDR and SDR/PRD respectively, confirming the predominance of the SDR mechanism.

For FinMac, 13 tetraploid hybrids were genotyped with seven centromeric markers (LGs 1, 2, 4, 6, 7, 8, and 9). Six inferred 2n gametes were totally homozygous for these markers (Table 5). Among these, two unreduced gametes (from FinMac 12 and FinMac 13) remained totally homozygous after analyzing six markers covering LG1 and were subjected to additional analysis to distinguish between the SDR and PMD hypothesis. The 11 2n gametes with at least one heterozygous locus produced LOD values >2 for SDR/FDR, rejecting the FDR hypothesis. Among these, four displayed LOD values of 2.81 for the SDR/PRD test and were therefore considered to have arisen from SDR. The seven remaining 2n gametes displayed positive values ranging from 0.52 to 1.91. These gametes had a higher probability of arising from SDR than from PRD, but this result is not conclusive because the values are below our threshold. The population LOD values were 78.84 and 19.81 for SDR/FDR and SDR/PRD, respectively, again confirming the prevalence of SDR. The seven 2n gametes with inconclusive individual LODs display a population LOD of 43.12 and 8.56 for SDR/FDR and SDR/PRD, respectively. It is therefore highly probable that they also arose from SDR.

To validate, at the population level, the finding that 38 2n gametes were derived by SDR (as determined by individual LOD analysis) and to distinguish between SDR and PRD for the eight gametes with inconclusive individual LODs, we compared the PHR patterns of the two set of gametes in LG1. For this analysis, we used four SSR markers (CIBE6126, mCrCIR06B05, MEST001, and MEST431) and two SNP markers (CiC2110-02 and CiC5950-02) (Figure 2) mapped in LG1 (Figures 3, 4).

For the conclusive SDR 2n gametes, the PHR values in LG1 (Figure 3) decreased from 67% for the telomeric marker CIBE6126 to 3% for the centromeric marker mCrCIR06B05 and progressively increased to 77% when moving toward the other telomeric marker, MEST431. The average PHR value was 42%. For the eight inconclusive 2n gametes, the same PHR pattern was observed: the lowest value was obtained for the centromeric marker mCrCIR06B05 (0%) and the highest for the telomeric markers (63% for CiC2110-02 in one telomere and 75% for MEST431 in the other). The average PHR for these eight gametes was 46% (Figure 3). These PHR patterns totally fit the profile for SDR. The average PHR value over the two sets of 2n gametes was 43%. Various studies have indicated that the global restitution of heterozygosity is expected to be near 80% for FDR and 40% for SDR, assuming a random distribution of heterozygous loci along the chromosomes (Peloquin, 1983; Hutten et al., 1994; Carputo et al., 2003). Both the patterns along LG1 and the average PHR values comply with the SDR hypothesis. Therefore, we conclude that the eight 2n gametes of indeterminate origin identified from the individual LOD (SDR/PRD) analysis also originated from SDR. Under this conclusion, the PHR pattern in LG 1 is very similar for “Eureka Frost” and “Fino” lemon SDR 2n gamete populations (Figure 4).

We performed additional analyses of the two inferred 2n gametes (FinMac 12 and FinMac 13 tetraploid plants) fully homozygous for the seven centromeric markers and the six markers of LG1. Fully homozygous 2n female gametes for centromeric loci can originate through SDR or PMD, with different consequences for the genetic structures of 2n gametes. Bastiaanssen et al. (1998) defined two conditions that are necessary to conclude that PMD rather than SDR has occurred, i.e., 100% homozygosity for all genotyped loci and the occurrence of recombination between homozygous alleles in the same LG. Therefore, we genotyped FinMac 12 and FinMac 13 using 11 telomeric loci found in different LGs to provide genetic evidence for a particular PMD mechanism. The average distance from these markers to their corresponding centromere is 53.22 cM (ranging from 25.32 to 89.59 cM). Both plants were homozygous for all molecular markers analyzed. Furthermore, C. limon is a direct hybrid between two genetically distant genotypes, C. aurantium and C. medica (Nicolosi et al., 2000; Curk et al., 2016), and the specific origins of the homozygous alleles can easily be distinguished. We found that some homozygous markers of the same LG were inherited from the C. aurantium ancestor and the others from C. medica. For example, multilocus analyses of the homozygous alleles in LG1 (Figure 5) revealed interspecific recombination in the two plants with alternation of homozygosity originated from C. aurantium and C. medica. Consequently, according to Bastiaanssen et al. (1998), the observation of 100% homozygosity and recombination between C. aurantium and C. medica along the same LG provides evidence discarding the SDR mechanism and leads us to conclude that these two 2n gametes originated through PMD. To our knowledge, this is the first report of the identification of a new mechanism, Post-Meiotic genome Doubling, leading to 2n ovule gametes in citrus, and specifically in lemon.

On the whole, we conclude that 38 (88%) of the 2n gametes analyzed had arisen from SDR, three (7%) from FDR or PRD, and two (5%) from PMD. At the population level, SDR appears to be by far the most common mechanism for 2n ovule formation in both C. limon genotypes, “Eureka Frost” and “Fino”. This is the first report of the production of a large number of lemon progenies from 2n gametes produced by different mechanisms of unreduced ovule gametes. Luro et al. (2004), Aleza et al. (2015), and Cuenca et al. (2015) also found that SDR was the predominant mechanism leading to 2n megagametophyte production in mandarins. Among the 19 mandarins investigated, the authors concluded that only 1.1 and 2.9% of plants were recovered from FDR in the “Ellendale” and “Fortune” mandarins, respectively. The coexistence of SDR and FDR has been recently observed in unreduced pollen gametes by Rouiss et al. (2017). 53 plants were obtained from 2n pollen gametes produced by a diploid hybrid between clementine and sweet orange. FDR was the predominant mechanism (77%) and SDR was the mechanism for the remaining plants (23%). In addition, FDR was the main mechanism for 2n female gamete production in “Femminello” lemon (Ferrante et al., 2010). These results are questionable because the authors used only a few molecular markers and lacked previous information about centromere location and the relative distances between the markers and the centromeres. With the recent location of centromeres in the citrus genetic map (Ollitrault et al., 2012a; Aleza et al., 2015), the markers used by Ferrante et al. (2010), JK-TAA1, JK-TAA15, JK-TAA41, and NB-GT03, are located at 87.29, 59.07, 74.99, and 50.47 cM from the centromere of the LGs 6, 1, 2, and 8 respectively, being mostly telomeric, and therefore the high PHR values obtained in their study can fit both SDR or FDR mechanisms.

At the methodological level, we demonstrated the power of using two complementary approaches, namely, analysis of the PHR pattern in one LG with the maximum-likelihood method proposed by Cuenca et al. (2015). Considering only centromeric loci, different mechanisms can lead to the same homozygous patterns. Therefore, analyzing the heterozygosity restitution pattern along LGs at the individual level is a useful approach for distinguishing between SDR and PMD, since, under this mechanism, the heterozygosity restitution value is zero for all markers in all LGs. After LOD analysis at the individual level, this method is used to analyze PHR patterns at the population level to distinguish between SDR and PRD when individual LODs are under the threshold required to obtain conclusive results. When enough number of individuals is analyzed, this technique should also be utilized to distinguish between FDR and PRD. With FDR-2n gametes, heterozygosity restitution varies from 100% in centromeric loci to close to 66% in telomeric areas under the non-interference model (Cuenca et al., 2011), whereas, with PRD, heterozygosity restitution is expected to be very similar along the entire chromosome.

Crossover interference ensures the appropriate distribution of crossovers along the chromosome, since one crossover reduces the likelihood of other crossovers occurring nearby (Youds et al., 2010). The analysis of crossover rates (Table 6) for both arms of chromosome I revealed the presence of up to four crossovers on one arm and three on the other arm. In addition, three complementary crossovers (double crossing over involving four chromatids) were observed as a result of phase-changing between two homozygous markers. Similarly, Cuenca et al. (2011) and Aleza et al. (2015) detected up to four crossovers on one arm and complementary crossovers in “Fortune” mandarin and C. clementina. We estimated the IC for each chromosome arm, finding partial interference in both arms (IC = 0.27 and 0.44). Such variation in interference values between both arms has also been observed in other citrus species, ranging from 0.82 to 0.48 for “Fina” clementine on LG 1 (Aleza et al., 2015) and 0.73 to 0.53 for “Fortune” mandarin on LG 2 (Cuenca et al., 2011). Variation in the level of interference between different parts of the genome has also been observed in Arabidopsis (Drouaud et al., 2007), humans (Lian et al., 2008), and mice (Broman et al., 2002).

Sexual polyploidization via 2n gametes and interploid sexual hybridizations using tetraploid parents (doubled diploids) are the main strategies used to produce triploid citrus hybrids (Ollitrault et al., 2008; Aleza et al., 2010b, 2012a,b, 2016; Navarro et al., 2015). These different strategies and the different meiotic behaviors result in different genetic structures in the diploid gametes and, consequently, the resulting triploid progenies. The three hybrids obtained via FDR or PRD 2n gametes have a higher rate of heterozygosity than hybrids obtained via SDR. By contrast, the two plants obtained by PMD transmit 0% of PHR (Bastiaanssen et al., 1998). Therefore, such a mechanism generally promotes inbreeding in the hybrid progenies (Tai, 1986; Gallais, 2003). However, these lines constitute interesting parentals to be used as test lines in inheritance studies (Bastiaanssen et al., 1998).

In addition, the mechanism that generates the 2n gametes affects the breeding efficiency for a character in relation to the genetic distance to the centromeres of the major genes controlling this character. For instance, Cuenca et al. (2013b, 2016) found that resistance to Alternaria brown-spot fungal disease is a recessive trait controlled by a single locus located at 10.5 cM from the centromere of chromosome III. Therefore, in crosses between a heterozygous parent producing diploid gametes and a resistant genotype, PMD is the most favorable mechanism (50% of resistant hybrids), followed by SDR (40%). Under FDR, only 5% of the hybrids will be resistant. For diploid gametes produced by the doubled-diploid genotype or resulting from PRD, the rates of resistant hybrids should vary from 16% (tetrasomic segregation) to 0% (disomic segregation) according to the preferential pairing behavior.

The aim of some lemon-breeding programs is to produce new lemon-like types of fruit, which essentially involves 2x X 4x crosses using diploid lemons as female parents and more or less complex hybrids as tetraploid parents (Recupero et al., 2005; Viloria and Grosser, 2005). This approach is used in an attempt to solve some of the problems caused by the low genetic variation of C. limon, although relatively few tetraploids are available. This approach has allowed for the selection and protection of the triploid “Lemox”, a hybrid between a diploid female complex hybrid, and tetraploid lemon (Recupero et al., 2005). “Lemox” produces quality fruits resembling lemons with high tolerance to Mal secco. The 2n lemon gametes will be very useful for producing new lemon-like seedless citrus types via 2x X 2x hybridizations, thereby dramatically increasing the gene pool of genotypes that could be used as parents. Furthermore, the production of 2n gametes has been investigated in a small number of lemon genotypes. Evaluating the many existing lemon genotypes may result in the detection of specific genotypes that produce higher rates of 2n gametes and (eventually) genotypes with different ratios of FDR and SDR 2n gametes, which will increase the efficiency of breeding programs.

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.