Original Research ARTICLE
Analysis of Crossover Events and Allele Segregation Distortion in Interspecific Citrus Hybrids by Single Pollen Genotyping
- 1Centro de Citricultura y Producción Vegetal, Instituto Valenciano de Investigaciones Agrarias, Valencia, Spain
- 2Concordia Agricultural Experiment Station, National Agricultural Technology Institute, Entre Ríos, Argentina
- 3Department of Breeding Research, Leibniz Institute of Plant Genetics and Crop Plant Research (IPK), Seeland, Germany
- 4Institute of Agricultural and Nutritional Sciences, Faculty of Natural Sciences III, Martin Luther University Halle-Wittenberg, Halle (Saale), Germany
In citrus, a classical method of studying crossovers and segregation distortion (SD) is the genetic analysis of progenies. A new strategy combining fluorescence-activated cell sorting and whole genome amplification of haploid pollen nuclei with a large set of molecular markers, offers the opportunity to efficiently determine the frequency of crossovers and the identification of SD without the need to generate segregating populations. Here we have analyzed meiotic crossover events in a pollen nuclei population from “Eureka” lemon and the allelic SD was evaluated in a pollen nuclei population from a clementine × sweet orange hybrid (“CSO”). Data obtained from the “CSO” pollen nuclei population were compared to those obtained from genotyping of a segregating population (“RTSO”) arising from a hand-made sexual hybridization between diploid non apomictic selected tangor (mandarin × sweet orange; “RTO” tangor) as female parent pollinated with “CSO” tangor as male parent. The analysis of crossovers rates on chromosome 1 revealed the presence of up to five crossovers events on one arm and four on the corresponding other arm, with an average of 1.97 crossovers per chromosome while no crossover events were observed in five “Eureka” lemon pollen nuclei. The rate of SD observed in “CSO” pollen nuclei (13.8%) was slightly lower than that recovered in the “RTSO” population (20.7%). In the pollen nuclei population, SD was found on linkage group (LG) 2, while the “RTSO” population showed SD on LGs 2 and 7. Potential male gametic selection mechanisms were distinguished in pollen grains, while in the population, mechanisms of gametophytic selection and/or zygotic selection were observed. This methodology is a very useful tool to facilitate research focused on the reproductive biology of citrus and study the mechanisms that affect crossovers and SD.
Most plant species that reproduce sexually alternate the life cycle between a diploid sporophytic phase and a reduced gametophytic haploid phase (Xu et al., 2013). Meiosis is the process in which the number of chromosomes is halved, leading to haploid gametes (Zamariola et al., 2014; Bomblies et al., 2015). The formation of crossovers between homologous parental chromosomes occurs in the prophase of meiosis I (Hunter, 2015; Lambing et al., 2017; Li et al., 2017), and consists in the exchange of genetic information between two non-sister homologous chromatids (Zamariola et al., 2014; Lambing et al., 2017). Crossovers are fundamental for chromosome segregation in most eukaryotes, but play also an important role in evolution and speciation of plants (Li et al., 2017), as they generate novel genetic combinations (Morgan et al., 2017). Many studies have demonstrated that the distribution of meiotic crossover events along chromosomes is non-random in plants, varying locally and displaying several hotspots (Mézard, 2006; Choi et al., 2008; Lambing et al., 2017).
On the other hand, allelic segregation distortion (SD), that is, the deviation of the allelic segregation ratio from the expected Mendelian-ratio, is increasingly recognized as an additional powerful evolutionary force (Dai et al., 2017). SD may result mainly from non-random segregation of gametes during meiosis (due to several causes such as the presence of deleterious alleles or gamete competition), post-meiotic gamete dysfunction or differential gamete success (such as differential pollen tube growth, pollen–pistil interactions, or partial incompatibility) and differential zygotic fitness (Li et al., 2014; Bodénès et al., 2016; Fishman and McIntosh, 2019; Seymour et al., 2019).
Many studies have been conducted in crossover formation and allelic SD in several plant species, such as Arabidopsis spp., wheat, barley, maize, tomato, and rice (Zamariola et al., 2014; Wang and Copenhaver, 2018). In citrus species, belonging to the family Rutaceae, subfamily Aurantioideae and with nine (x = 9) as basic chromosome number (Krug, 1943), the analysis of the crossover rate was determined in mandarin and lemon by analyzing progenies with molecular markers. Results showed that up to four crossover events per chromosome arm can be found (Cuenca et al., 2011; Aleza et al., 2015; Rouiss et al., 2017b). In addition, skewed markers appear to be concentrated in specific areas, usually named SD regions, which are different for the different genotypes analyzed (Ollitrault et al., 2012a), and could contain deleterious genes or be associated with hybridization incompatibility (Zhou et al., 2015). In addition to parental genotype, differential aptitude of gamete genotypes, direction of the crosses, and regulatory gene interactions can also contribute to the high level of SD in citrus (Bernet et al., 2010; Curtolo et al., 2017). Such SD was proposed to result from gametic selection rather than zygotic selection for the analyzed crosses (Bernet et al., 2010; Ollitrault et al., 2012a).
A classical approach to locate crossovers is analyzing the co-inheritance of linked heterozygous markers through meiosis in segregating populations, which requires precise genetic maps and the availability of known genome structures (Choi and Henderson, 2015). Strategies for identifying whether gametic or zygotic selection (or both) is causing SD are based on genotyping segregating populations from reciprocal and mixed crosses (Ruiz and Asins, 2003; Bernet et al., 2010; Ollitrault et al., 2012a; Xu et al., 2013; Bodénès et al., 2016; Fishman and McIntosh, 2019; Seymour et al., 2019). Alternative methods involve cytogenetic analyses based on recombination nodule maps, chiasma counting or C-band maps (Mézard, 2006; Li et al., 2017).
One alternative strategy for investigating both meiotic recombination patterns and allelic SD is to perform direct genotyping of the gametes, rather than individuals from segregating populations. This methodology paves the way for new information into the genetic basis of non-Mendelian inheritance in contrast to population-based analyses (Choi et al., 2008; Dreissig et al., 2015, 2017; Sun et al., 2019).
The use of the fluorescence activated cell sorting (FACS) technique (Galbraith, 2010) coupled with whole genome amplification (WGA) allows the analysis of individual nuclei with a high number of molecular markers (Dreissig et al., 2015, 2017). In this context, genotyping of individual pollen grains nuclei can be useful for the determination of the male parent haplotypes, the evaluation of meiotic recombination and potential allelic SD (Mase et al., 2014; Dreissig et al., 2017). This strategy also allows performing studies on the genetic structures of pollen grain populations as compared with those originated at the plant level, without interferences due to a potential partial cross-incompatibility or gamete competition (Gu et al., 2013). In citrus, the strategy of ploidy level analysis, sorting, and genotyping of single pollen grains has been previously assessed for several diploid, triploid, and tetraploid genotypes (Garavello et al., 2019), with successful results when applied to haploid pollen grains.
In the current study, we used a strategy combining flow-sorting of pollen nuclei with WGA and genotyping of individualized nuclei by simple sequence repeats (SSR) and single nucleotide polymorphism (SNP) molecular markers. The meiotic crossover events were investigated in a pollen grain population derived from “Eureka” lemon [Citrus limon (L.) Burm. f.] and the allelic SD was evaluated in a pollen grain population derived from a clementine × sweet orange hybrid (C. clementina × C. sinensis; hereafter referred to as “CSO” tangor). Data obtained from pollen grain population were compared to those obtained from genotyping of a segregating population recovered from a hand-made hybridization between diploid non apomictic selected tangor (C. reticulata × C. sinensis) as female parent pollinated with “CSO” tangor as male parent.
Materials and Methods
Plant material was sampled from the parental collection of our triploid breeding program (Navarro et al., 2015) carried out in the Instituto Valenciano de Investigaciones Agrarias (IVIA) located at Moncada (Valencia, Spain).
For crossover analysis, we isolated and whole-genome amplified pollen nuclei from the diploid “Eureka” lemon. Lemon resulted from an ancient direct interspecific hybridization between sour orange (C. aurantium L.) as female parent and citron (C. medica L.) as male parent (Nicolosi et al., 2000; Froelicher et al., 2011; Garcia-Lor et al., 2013b; Curk et al., 2016; Wu et al., 2018). Since sour orange and citron are genetically distant species, the identification of alleles from each parent is greatly facilitated (Wu et al., 2018).
On the other hand, for the SD analysis, we used a progeny of 86 diploid hybrids arising from a hand-made cross between diploid non apomictic selected tangor (C. reticulata × C. sinensis; hereafter referred to as “RTO” tangor) as female parent pollinated with “CSO” tangor as male parent, hereafter referred as “RTSO” population. This hybridization is of great interest in citrus breeding programs to combine desired attributes from both parents such as late maturing from “RTO” and anthocyanin content from “CSO.” Isolated and whole-genome amplified pollen nuclei from “CSO” tangor were genotyped to be compared with the diploid hybrids belonging to the “RTSO” population. Further information about recovery and ploidy level analysis of regenerated diploid plants can be found in Aleza et al. (2009b).
DNA Extraction From Leaves
Genomic DNA was extracted from “RTO” and “CSO” tangors, “RTSO” diploid progeny, “Eureka” lemon, sour orange and citron, using a Plant DNeasy kit from Qiagen Inc. (Valencia, CA, United States) following the manufacturer’s protocol and measured using a spectrophotometer (NanoDrop 2000C, Thermo Fisher). The samples were diluted with sterile water (Sigma-Aldrich, Co., United Kingdom) at a concentration of 10 ng/μl and stored at −20°C until use.
Fluorescence Activated Cell Sorting-Based Isolation of Single Pollen Nuclei and Whole Genome Amplification
For pollen grain recovery, between 40 and 50 flowers of the “CSO” tangor and “Eureka” lemon were collected in pre-anthesis during the spring of 2018. The anther’s were removed from the flowers with forceps and were dried in Petri dishes over silica gel in a desiccator at room temperature until the anthers opened after 1 or 2 days. Then, dehiscent anthers were selected under magnifying glass, discarding those that were not fully open, and they were sealed with parafilm and stored at −20°C until use. To isolate pollen nuclei we used the methodology described by Kron and Husband (2012) and validated in citrus by Garavello et al. (2019). Briefly, 4–5 dehisced anthers were vortexed in a small tube with 300 μl of nuclei isolation buffer in order to recover all pollen grains. Afterward, the suspension was filtered using the pre-filter and the bursting filter (CellTrics® filters, Partec®). With the help of a plastic rod, collected pollen grains over the bursting filter were rubbed against the filter and washed with nuclei isolation buffer, twice. Subsequently, DAPI staining solution (1.5 μg/ml) was added to the suspension and incubated for 10 min. Stained suspension were run in a BD Influx (BD Biosciences, United States) and analyzed with BD FACS software.
We have followed the methodology described by Dreissig et al. (2015) to perform FACS-based purification of single nuclei and WGA. From the nuclei suspension, individualized nuclei were sorted into individual wells of a 384-microwell plate containing 2 μl lysis solution, which was composed of 0.5 μl lysis buffer, 0.5 μl ddH2O, and 1 μl sample buffer (Genomiphi V2, GE Healthcare). Illustra GenomiPhi V2 DNA Amplification Kit (GE Healthcare, United States) was used to WGA Prior to calculate the DNA concentration of WGA products, each sample was diluted with 500 μl of ddH2O, and subsequently was measured by fluorometric quantification (Qubit, Life Technologies).
Fifty-four pollen nuclei of “CSO” tangor and 44 pollen nuclei of “Eureka” lemon were sorted and whole genome amplified. Additionally, 10 single pollen nuclei of the analyzed genotypes were mixed in the same well to be used as a positive control against amplification errors (Dreissig et al., 2015).
Genotyping Using Simple Sequence Repeat and Single Nucleotide Polymorphism Markers
The hybrids of “RTSO” progeny and their parents together with the amplified pollen nuclei of “CSO” tangor were genotyped using 30 molecular markers (23 SSRs and seven SNPs) heterozygous for “CSO” and with polymorphism with “RTO.” These markers are distributed across all linkage groups (LGs) of the Clementine genetic map (Ollitrault et al., 2012a). Isolated pollen nuclei of “Eureka” lemon were genotyped using seven SSRs and five SNPs heterozygous markers. These markers are distributed across LG1 of the Clementine genetic map (Ollitrault et al., 2012a). Detailed information about all markers used is given in Tables 1, 2.
Table 1. Information about molecular markers used for genotyping “CSO” tangor pollen nuclei, “RTSO” progeny and its parents, indicating accession number in Gene Bank or Phytozome, position in the reference clementine genetic map, noted alleles, and reference.
Table 2. Information about molecular markers used for genotyping “Eureka” lemon pollen nuclei, indicating accession number in Gene Bank or Phytozome, position in the reference clementine genetic map, noted alleles and reference.
Polymerase chain reactions (PCRs) using SSR markers were performed using a Thermocycler rep gradient S (Eppendorf®) using the following protocol: reaction volume, 15 μl containing 0.5 μl of 1 U/μl of Taq DNA polymerase (Fermentas®), 3 μL of citrus template DNA (10 ng/μl), 1.5 μl of 2 mM welled (Sigma®) dye-labeled forward primer, 1.5 μl of 2 mM non-dye-labeled reverse primer, 0.2 mM of each dNTP, 1.5 μl of PCR reaction buffer 10×, and 0.45 μl of 50 mM MgCl2. The cycling program was set as follows: denaturation for 5 min at 94°C followed by 40 cycles of 30 s at 94°C, 30 s at 50°C or 55°C, 30 s at 72°C; and a final elongation step of 8 min at 72°C. Separation was carried out by capillary gel electrophoresis using a Genetic Analysis System 8000 (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). GenomeLabTM v.10.0 (Beckman Coulter Inc.) genetic analysis software was used for data collection.
For SNP markers genotyping we used KASParTM technology by LGC Genomics1. Primers were directly designed by LGC Genomics from each SNP locus flanking sequence, considering approximately 50 nt on each side of the SNP. The KASPar genotyping system is a competitive allele-specific dual Förster resonance energy transfer (FRET)-based assay for SNP genotyping, and detailed explanation of specific conditions and reagents can be found in Cuppen (2007).
Analysis of Segregation Distortion and Crossovers
The potential allelic SD was analyzed in the “RTSO” progeny as well as in the pollen nuclei population from “CSO” tangor using the Chi-square test (χ2), assuming an expected allelic segregation ratio of 1:1 (p < 0.05) for each analyzed marker.
The crossover events along the chromosome 1 were detected in “Eureka” pollen nuclei by identifying changes in alleles inherited from sour orange to citron, and vice versa. The number of crossovers was estimated for each arm of the chromosome 1 according to the following calculation:
where n (CO) is the number of crossovers and n (M) is the total number of observations. The pairs of markers in the analyzed gametes that showed missing data points were omitted from the analysis.
Analysis of Population Diversity
Population diversity organization between hybrids of the “RTSO” progeny at the male gamete level and the “CSO” pollen nuclei population was examined using DARwin6 software (Perrier and Jacquemound-Collet, 2018). Neighbor-joining analysis using the simple matching dissimilarity index (di-j) between pairs of markers (units) was performed:
where di–j is the dissimilarity between units i and j, L is the number of loci, ml is the number of matching alleles for locus l, and π is the ploidy level, which is one in this case. From the dissimilarity matrix obtained, a Weighted neighbor-joining tree (Saitou and Nei, 1987) was computed.
Results and Discussion
Fluorescence Activated Cell Sorting of Pollen Nuclei, Whole Genome-Amplification, and Genotyping
In a previous work, we have demonstrated that FACS technique coupled with WGA is an adequate methodology for multi-locus SSR and SNP genotyping of citrus haploid pollen nuclei (Garavello et al., 2019). In the present work, a total of 44 and 54 haploid pollen nuclei from “Eureka” lemon and “CSO” tangor were analyzed, respectively.
“Clementine × sweet orange hybrid” tangor pollen nuclei were genotyped with the same set of molecular markers as the “RTSO” population (Table 1). For further analysis, we selected those gametes with a minimum of 65% positive amplifications (Figure 1) for the set of markers used. Thirty-four “Eureka” lemon (77.3%) and 48 “CSO” pollen nuclei (88.9%) showed at least 65% positive PCR reactions.
Figure 1. Representative dispersion diagram of PCR products obtained with CIC5087-01 SNP marker of haploid pollen grains nuclei amplified with WGA kit from the “CSO” diploid genotype.
Out of 408 PCR reactions performed for all marker combinations in “Eureka” lemon (Table 2), 362 were positive (88.7%). On the other hand, the “CSO” tangor showed 1,393 positive PCR reactions (96.7%) of a total of 1,440.
The percentage of WGA positive PCR in “Eureka” lemon was slightly lower than that found by Dreissig et al. (2015) when analyzed barley pollen nuclei with SNP markers. However, this value was more similar for “CSO.” In the same way, these results are in agreement with what was found by Garavello et al. (2019), who conducted a classification and genotyping study of individual haploid citrus pollen nuclei.
The genotyping of the 86 “RTSO” diploid hybrids and their parents, performed with the same set of markers used for “CSO” pollen nuclei (Table 1) allowed the unequivocal allelic differentiation between both parents and assessing the allelic contribution of the male parent to form each diploid hybrid (Figure 2).
Figure 2. Electropherograms of a diploid hybrid “RTSO-06” recovered from hybridization between “RTO” diploid tangor as female parent and the “CSO” diploid tangor as male parent with CX06F23 SSR marker.
Analysis of Crossover Events in Haploid Pollen Nuclei of “Eureka” Lemon
We analyzed the number of crossovers for both arms on chromosome 1 (Table 3). The analysis of crossover rates revealed the presence of up to five crossovers events on one arm and four on the other arm. Comparatively, the average number of crossovers on chromosome 1 in the pollen nuclei population from “Eureka” lemon was 1.97.
Table 3. Number of observed crossover events on each arm of chromosome 1 based on analysis of haploid pollen nuclei from “Eureka” lemon using twelve SSR and SNP markers.
Interestingly, we found five pollen nuclei (14.7%) with no crossover events between analyzed markers (Table 4), resulting with the same configuration as one of the lemon parents, i.e., citron or sour orange. Three pollen nuclei showed the same configuration as citron (samples Eur-P4, Eur-P9, and Eur-P19) while two pollen nuclei showed the same configuration as sour orange (samples Eur-P12 and Eur-P26) although it is possible that the crossovers in these pollen nuclei escaped to our detection, since in none of the five pollen nuclei all markers were amplified. The lack of amplification of these alleles may be attributed to the WGA kit performance at these specific loci since in the lemon leaf DNA PCR amplifications always happened (Table 4). In addition, other explanation could be related with the presence of null alleles in the parent species of the “Eureka” lemon. The remaining 29 pollen nuclei (85.3%) showed at least one recombination event between both genomes. However, our results indicate a tendency that were similar to those of Rouiss et al. (2017b) who, using six SSR markers detected up to four crossovers on one chromosome arm and three on the other arm when analyzing 27 unreduced gametes recovered from “Eureka Frost” lemon pollinated either with C. ichangensis or “Fortune” mandarin. Similarly, Aleza et al. (2015) and Cuenca et al. (2011) identified up to four crossovers on one chromosome arm in unreduced gametes of C. clementina and “Fortune” mandarin, respectively. Furthermore, Ollitrault et al. (2012a) identified up to two recombination break points on three LGs in the sweet orange gamete that originated C. clementina, and three crossover events on one LG in the C. clementina gamete that originated the haploid Clementine used by the International Citrus Genomic Consortium (ICGC) to establish the reference citrus haploid whole genome sequence (Aleza et al., 2009a).
Table 4. Multilocus configuration of the haploid pollen nuclei population from “Eureka” lemon analyzed with twelve SSR and SNP markers located on both arms of chromosome 1.
In addition, in several organisms the frequency of meiotic crossover formation is determined genetically and also by the influence of environmental conditions. In plants, exposure to biotic and abiotic stresses can modify the overall rate of recombination (De Storme and Geelen, 2014). For example, a change in temperature from 22 to 30°C caused a reduction and alter distribution of meiotic crossover formation in barley (Hordeum vulgare) (Higgins et al., 2012), whereas a change in temperature from 18 to 28°C caused an increase of meiotic crossover formation in Arabidopsis (Francis et al., 2007). Similarly, in maize (Zea mays), both low temperature and water deficit significantly enhance the frequency of crossovers (De Storme and Geelen, 2014).
On the other hand, recombination frequency was determined along chromosome 1 by counting the number of crossovers at intervals of neighboring markers (Figure 3). The highest frequencies of recombination were found toward the distal region of the chromosome, while the lowest frequencies were found near the centromere, which is in agreement with what was observed in C. clementina where centromeric areas showed low recombination rates (<1.0 cM/Mb), although this reduction changed greatly between chromosomes (Aleza et al., 2015). Suppression of crossovers in centromeric and pericentromeric regions has been displayed in many plant species like wheat, barley, and tomato (Sherman and Stack, 1995; Dreissig et al., 2015; Darrier et al., 2017; Blary and Jenczewski, 2019).
Figure 3. Recombination frequency along chromosome 1 determined by haploid pollen nuclei genotyping of “Eureka” lemon. Red point indicated relative position of centromere according to Aleza et al. (2015).
This methodology is of great interest to study in citrus how environmental stress can alter the number and the distribution of crossover in haploid gametes and its implication in breeding programs based on sexual hybridization.
Genotyping of “CSO” Haploid Pollen Nuclei Population and One Progeny of Diploid Hybrids Recovered With the Same Genotype as Male Parent
Conventional strategy to identify deviations from Mendelian segregation are based in genetic mapping populations from reciprocal and mixed crosses (Fishman and McIntosh, 2019). However, this conventional strategy can be improved by FACS, WGA, and genotyping of haploid pollen nuclei because with this new approach we might identify different mechanisms underlying SD that cannot be observed with the traditional way. To perform a comparative analysis between allele segregation in the “CSO” pollen nuclei population and the progeny recovered through hand-made pollination between “RTO” tangor as female parent and the “CSO” tangor as male parent (“RTSO” population), the same set of molecular markers was used distributed in the nine LGs of the Clementine genetic map (Ollitrault et al., 2012a). The preferential transmission of one allele over another in a statistically significant deviation from the expected Mendelian segregation ratio 1:1 is known as SD (Bodénès et al., 2016; Dai et al., 2017; Dreissig et al., 2017). Table 5 shows the number and all the Chi-square values of the “CSO” alleles in the pollen nuclei and the “RTSO” populations that were analyzed (Supplementary Tables 1, 2).
Table 5. Analysis of Mendelian allelic segregation (Chi-square test) for “CSO” tangor pollen nuclei population and “RTSO” progeny.
The markers displayed allelic SD in the two populations were analyzed. The rate of SD observed in “CSO” pollen nuclei (13.3%) was slightly lower than that recovered in “RTSO” population (20%). In the pollen nuclei population, SD was found on LGs 2 and 5, while the “RTSO” population showed SD on LGs 2 and 7. Globally, a total of seven analyzed markers showed SD. Among these markers, 2P21022555, CX6F23, and CIC3712-01, all located on LG 2, were synchronously distorted in the two populations. However, MEST88 marker, located on LG 5, shows SD for the pollen nuclei population, and not in the “RTSO” progeny. In addition, MEST107, MEST202, and CIC3674-02 markers, all located on LG 7, showed SD only in the “RTSO” population.
Segregation distortion can be originated by different processes that include non-random segregation of gametes during meiosis, alteration of viability or functionality of gametes after meiosis, and differential zygotic fitness (Fishman and McIntosh, 2019; Seymour et al., 2019). Synchronic SD visualized in the two gamete populations on LG 2 could be related with a gamete selection during pollen meiotic process. In other species like maize, monkeyflowers (Mimulus sp.), cotton (Gossypium hirsutum) it has been observed gametal factors that influence male gametes viability (Fishman and Saunders, 2008; Xu et al., 2013; Dai et al., 2017). For example, in monkey flowers Fishman and Saunders (2008) observed that the locus D strongly affects pollen viability and DD homozygotes suffered a 20% pollen viability reduction contributing to male fitness variation.
In addition, most loci with SD tend to cluster in segregation distortion regions (SDRs) as we have observed in both populations of “CSO” male gametes on LG 2. One explanation for SDRs could be that specific loci in the genome are conducted to viability differentiation (Luo and Xu, 2003; Li et al., 2011). The selection of an allele at the locus would result in nearby markers that deviate from the expected ratio, consistent with the theory of genetic hitchhiking (Dai et al., 2017).
Likewise, a strong SD has been observed on LG 7 only in “RTSO” progeny rather than “CSO” pollen nuclei population (Table 5, Supplementary Table 2). The gametophytic incompatibility system (GIS) was identified in citrus as one of the pollen–pistil interaction mechanisms that cause SD, defined as “the inability of a fertile hermaphrodite seed plant to produce zygotes after self-pollination” (Soost, 1965; de Nettancourt, 1977). Recently, Liang et al. (2019) have demonstrated that in citrus operate the S-RNase-based gametophytic self-incompatibility system that hold a S-RNase linked to about nine S-locus F-box genes. In this system, incompatibility reaction arises from the cytotoxic activity of S-RNase, meanwhile compatible pollen tubes avoid S-RNase cytotoxicity and grows into the style reaching the ovary (McClure et al., 2011). This system has been described in many genotypes including ancient and cultivated citrus species like mandarins and its hybrids (Liang et al., 2019). The GIS could be a factor for male gametic selection and this may lead to a complete exclusion of one allele for the concerned locus as we have noted in the three markers located on this LG with a very high distortion (4/82, 86/0, and 85/1 for MEST107, MEST202, and CIC3674-02 markers, respectively). In this context and taking into account the S-RNase-based GIS (Liang et al., 2019), “RTO” and “CSO” tangors have C. sinensis as the same ancestor and could share a S-RNase haplotype, whereby segregate as could be expected for the GSI system. These results could also suggest that in this genomic region could be located alleles related with the GIS in citrus. In fact, Liang et al. (2019) locate the S locus on LG 7 which is the same one we have observed SD. In previous results (data not shown) we could not recover any hybrid using “CSO” tangor as male parent in hand-made pollinations using clementines as female parent, suggesting incompatibility of this tangor with other genotypes. Nevertheless we cannot rule out zygotic or post-zygotic mechanisms involved in SD (Sweigart and Willis, 2012).
On Figure 4 we display the genetic relationship of these two populations that was calculated by neighbor-joining analysis using the simple matching dissimilarity index, allowing the differentiation of gamete groups within each population and determining their genetic diversity. This figure shows a gamete cluster produced by “CSO” haploid pollen nuclei that were absent in the “RTSO” gamete population, allowing to conclude that in citrus some kind of selection occurs during the progamic phase, either in the pollen grain germination on the stigmatic surface, pollen tube growth into the style (GIS) or differences in the zygotic or post-zygotic viability, precluding the recovery of hybrids with specific allelic configurations.
Figure 4. Neighbor-joining tree obtained from allelic data of “CSO” pollen nuclei population (red), and the male gametes from “RTSO” progeny (green).
Each gametes population showed different patterns of marker segregation, demonstrating that the genetic mechanism of SD has a specificity in each gametes population resulting from complex genetic system (Dai et al., 2017). These variations in the SD agree with different widely documented studies in plants, obtained through interspecific crosses for genetic mapping (Xu et al., 1997, 2013; Shirasawa et al., 2010; Bodénès et al., 2016; Dai et al., 2017; Fishman and McIntosh, 2019; Seymour et al., 2019), where large differences in SD values were found according to the direction of the crosses. The first distortion studies associated with sex were carried out in Drosophila obscura (Gershenson, 1928) because sexual dimorphism is common, and the sex ratio deviated a lot from 1:1 (Fishman and McIntosh, 2019; Seymour et al., 2019). Since then, several articles revealed a variety of meiotic and post-meiotic processes responsible for SD (Soltis et al., 1993; Fishman and Saunders, 2008; McDermott and Noor, 2010; Xu et al., 2013; Lindholm et al., 2016; Lambing et al., 2017). In citrus, Bernet et al. (2010) reported a higher SD in the male parent rather than in the female in a reciprocal hybridization between “Chandler” pummelo and “Fortune” mandarin. Subsequently, SD was evidenced by Ollitrault et al. (2012a) in clementine, pummelo and sweet orange when they were used either as female or male parents. From a breeding point of view, the presence SDRs implies that frequencies of relevant genes located there could be inherited at frequencies different from expected, affecting the efficiency of citrus breeding programs based on sexual hybridizations (Bernet et al., 2010).
In citrus, papers published until now about SD have been carried out mainly studying reciprocal crosses, as we have described before. In this work, we display for the first time a multi-locus genotyping study of the “CSO” haploid pollen gametes with and without the interference of the female parent. This approach allows us to identify mechanisms that could be related with specific genomic regions associated to non-random segregation of gametes during meiosis (as we have shown on LG 2) or with male and female gametic interactions, or zygotic mechanisms (LG 7). In addition, this methodology may have an important advantage for the achievement of sequencing projects. High heterozygosity is a general characteristic of Citrus species, making it difficult to assemble large genome sequences. For this reason, the ICGC decided to establish a reference whole citrus genome sequence from a clementine homozygous plant recovered by in situ gynogenesis induced by irradiated pollen Aleza et al. (2009a). Therefore, obtaining haplotypes from FACS coupled with WGA of haploid pollen nuclei would allow the sequencing projects to be approached more easily, from any genotype regardless its heterozygosity, and without the need to recover haploid plants using in vitro techniques, which in many cases is a very difficult task with very low efficiency.
Fluorescence activated cell sorting coupled with WGA allows the analysis of individual nuclei with a large number of molecular markers, offering the opportunity to efficiently determine the frequency of crossovers and the SD in haploid pollen nuclei without the need to generate segregating populations.
Until now, no study was conducted regarding the recombination in citrus pollen nuclei. The analysis of the “Eureka” lemon pollen nuclei allowed the identification of recombination points through the use of SSR and SNP markers, showing a greater number of crossovers in centromere distal regions of chromosome 1. In addition, SD has been observed either on LG 2 of both populations and only on LG 7 we have observed that the SD in the pollen of “CSO” tangor and the plants of the population differ in frequency or position. Potential male gametic selection mechanisms were distinguished in pollen grains, while in the population, mechanisms of gametophytic selection and/or zygotic selection were observed. The methodology presented here represents a very useful tool to facilitate research focused on the reproductive biology of citrus and study the mechanisms that affect crossovers and SD.
Data Availability Statement
All datasets generated for this study are included in the article/Supplementary Material.
MG and PA conceived the study and were in charge of the direction and planning. JF isolated pollen nuclei and conducted flow-sorting. SD performed whole genome amplification from single pollen nuclei. MG, JC, and PA genotyped pollen nuclei, diploid hybrids, analyzed the data, and took the lead writing the manuscript with input and review of LN, SD, JF, and AH. All authors read and approved the final version of this manuscript.
This work was supported by the project RTA2015-00069-00-00 from the Ministerio de Economía y Competividad, co-financed by FEDER (Euopean Fund For Regional Development) and by the grant “Programa de Perfeccionamiento” resolution 1177/14 from Consejo Directivo del Instituto Nacional de Tecnología Agropecuaria de Argentina.
Conflict of Interest
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.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2020.00615/full#supplementary-material
Aleza, P., Cuenca, J., Hernández, M., Juárez, J., Navarro, L., and Ollitrault, P. (2015). Genetic mapping of centromeres in the nine Citrus clementina chromosomes using half-tetrad analysis and recombination patterns in unreduced and haploid gametes. BMC Plant Biol. 15:80. doi: 10.1186/s12870-015-0464-y
Aleza, P., Froelicher, Y., Schwarz, S., Agustí, M., Hernández, M., Juárez, J., et al. (2011). Tetraploidization events by chromosome doubling of nucellar cells are frequent in apomictic citrus and are dependent on genotype and environment. Ann. Bot. 108, 37–50. doi: 10.1093/aob/mcr099
Aleza, P., Juárez, J., Hernández, M., Pina, J. A., Ollitrault, P., and Navarro, L. (2009a). Recovery and characterization of a Citrus clementina Hort. ex Tan. “Clemenules” haploid plant selected to establish the reference whole Citrus genome sequence. BMC Plant Biol. 9:110. doi: 10.1186/1471-2229-9-110
Bernet, G. P., Fernandez-Ribacoba, J., Carbonell, E. A., and Asins, M. J. (2010). Comparative genome-wide segregation analysis and map construction using a reciprocal cross design to facilitate citrus germplasm utilization. Mol. Breed. 25, 659–673. doi: 10.1007/s11032-009-9363-y
Bodénès, C., Chancerel, E., Ehrenmann, F., Kremer, A., and Plomion, C. (2016). High-density linkage mapping and distribution of segregation distortion regions in the oak genome. DNA Res. 23, 115–124. doi: 10.1093/dnares/dsw001
Chen, C., Bowman, K. D., Choi, Y. A., Dang, P. M., Rao, M. N., Huang, S., et al. (2008). EST-SSR genetic maps for Citrus sinensis and Poncirus trifoliata. Tree Genet. Genomes 4, 1–10. doi: 10.1007/s11295-007-0083-3
Choi, K., Yelina, N. E., Serra, H., and Henderson, I. R. (2008). “Quantification and sequencing of crossover recombinant molecules from Arabidopsis pollen DNA,” in Encyclopedia of Genetics, Genomics, Proteomics and Informatics, ed. G. P. Rédei (Dordrecht: Springer), 842–842. doi: 10.1007/978-1-4020-6754-9_7341
Cuenca, J., Froelicher, Y., Aleza, P., Juárez, J., Navarro, L., and Ollitrault, P. (2011). Multilocus half-tetrad analysis and centromere mapping in citrus: evidence of SDR mechanism for 2n megagametophyte production and partial chiasma interference in mandarin cv “Fortune”. Heredity 107, 462–470. doi: 10.1038/hdy.2011.33
Curk, F., Ancillo, G., Ollitrault, F., Perrier, X., Jacquemoud-Collet, J.-P., Garcia-Lor, A., et al. (2015). Nuclear species-diagnostic SNP markers mined from 454 amplicon sequencing reveal admixture genomic structure of modern Citrus varieties. PLoS One 10:e0125628. doi: 10.1371/journal.pone.0125628
Curk, F., Ollitrault, F., Garcia-Lor, A., Luro, F., Navarro, L., and Ollitrault, P. (2016). Phylogenetic origin of limes and lemons revealed by cytoplasmic and nuclear markers. Ann. Bot. 117, 565–583. doi: 10.1093/aob/mcw005
Curtolo, M., Cristofani-Yaly, M., Gazaffi, R., Takita, M. A., Figueira, A., and Machado, M. A. (2017). QTL mapping for fruit quality in Citrus using DArTseq markers. BMC Genomics 18:289. doi: 10.1186/s12864-017-3629-2
Dai, B., Guo, H., Huang, C., Ahmed, M. M., and Lin, Z. (2017). Identification and characterization of segregation distortion loci on cotton chromosome 18. Front. Plant Sci. 7:2037. doi: 10.3389/fpls.2016.02037
Darrier, B., Rimbert, H., Balfourier, F., Pingault, L., Josselin, A.-A., Servin, B., et al. (2017). High-resolution mapping of crossover events in the hexaploid wheat genome suggests a universal recombination mechanism. Genetics 206, 1373–1388. doi: 10.1534/GENETICS.116.196014
De Storme, N., and Geelen, D. (2014). The impact of environmental stress on male reproductive development in plants: biological processes and molecular mechanisms. Plant Cell Environ. 37, 1–18. doi: 10.1111/pce.12142
Dreissig, S., Fuchs, J., Cápal, P., Kettles, N., Byrne, E., and Houben, A. (2015). Measuring meiotic crossovers via multi-locus genotyping of single pollen grains in barley. PLoS One 10:e0137677. doi: 10.1371/journal.pone.0137677
Dreissig, S., Fuchs, J., Himmelbach, A., Mascher, M., and Houben, A. (2017). Sequencing of single pollen nuclei reveals meiotic recombination events at megabase resolution and circumvents segregation distortion caused by postmeiotic processes. Front. Plant Sci. 8:1620. doi: 10.3389/fpls.2017.01620
Francis, K. E., Lam, S. Y., Harrison, B. D., Bey, A. L., Berchowitz, L. E., and Copenhaver, G. P. (2007). Pollen tetrad-based visual assay for meiotic recombination in Arabidopsis. Proc. Natt. Acad. Sci. 104, 3913–3918. doi: 10.1073/pnas.0608936104
Froelicher, Y., Dambier, D., Bassene, J. B., Costantino, G., Lotfy, S., Didout, C., et al. (2008). Characterization of microsatellite markers in mandarin orange (Citrus reticulata Blanco). Mol. Ecol. Resour. 8, 119–122. doi: 10.1111/j.1471-8286.2007.01893.x
Froelicher, Y., Mouhaya, W., Bassene, J.-B., Costantino, G., Kamiri, M., Luro, F., et al. (2011). New universal mitochondrial PCR markers reveal new information on maternal citrus phylogeny. Tree Genet. Genomes 7, 49–61. doi: 10.1007/s11295-010-0314-x
Garavello, M., Cuenca, J., Dreissig, S., Fuchs, J., Houben, A., and Aleza, P. (2019). Assessing ploidy level analysis and single pollen genotyping of diploid and euploid Citrus genotypes by fluorescence-activated cell sorting and whole-genome amplification. Front. Plant Sci. 10:1174. doi: 10.3389/fpls.2019.01174
Garavello, M., Cuenca, J., Garcia-Lor, A., Ortega, N., Navarro, L., Ollitrault, P., et al. (2020). Male and female inheritance patterns in tetraploid ‘Moncada’ mandarin. Plant Cell Rep. 39, 335–349. doi: 10.1007/s00299-019-02494-y
Garcia-Lor, A., Ancillo, G., Navarro, L., and Ollitrault, P. (2013a). Citrus (Rutaceae) SNP markers based on competitive allele-specific PCR; transferability across the aurantioideae subfamily. Appl. Plant Sci. 1:1200406. doi: 10.3732/apps.1200406
Garcia-Lor, A., Curk, F., Snoussi-Trifa, H., Morillon, R., Ancillo, G., Luro, F., et al. (2013b). A nuclear phylogenetic analysis: SNPs, indels and SSRs deliver new insights into the relationships in the ‘true citrus fruit trees’ group (Citrinae, Rutaceae) and the origin of cultivated species. Ann. Bot. 111, 1–19. doi: 10.1093/aob/mcs227
García-Lor, A., Luro, F., Navarro, L., and Ollitrault, P. (2012). Comparative use of InDel and SSR markers in deciphering the interspecific structure of cultivated citrus genetic diversity: a perspective for genetic association studies. Mol. Genet. Genomics 287, 77–94. doi: 10.1007/s00438-011-0658-4
Gu, C., Liu, Q.-Z., Yang, Y.-N., Zhang, S.-J., Khan, M. A., Wu, J., et al. (2013). Inheritance of hetero-diploid pollen S-haplotype in self-compatible tetraploid Chinese cherry (Prunus pseudocerasus Lindl). PLoS One 8:e61219. doi: 10.1371/journal.pone.0061219
Higgins, J. D., Perry, R. M., Barakate, A., Ramsay, L., Waugh, R., Halpin, C., et al. (2012). Spatiotemporal asymmetry of the meiotic program underlies the predominantly distal distribution of meiotic crossovers in barley. Plant Cell 24, 4096–4109. doi: 10.1105/tpc.112.102483
Kijas, J. M. H., Thomas, M. R., Fowler, J. C. S., and Roose, M. L. (1997). Integration of trinucleotide microsatellites into a linkage map of Citrus. Theor. Appl. Genet. 94, 701–706. doi: 10.1007/s001220050468
Li, F., De Storme, N., and Geelen, D. (2017). Dynamics of male meiotic recombination frequency during plant development using fluorescent tagged lines in Arabidopsis thaliana. Sci. Rep. 7:42535. doi: 10.1038/srep42535
Li, G., Serba, D. D., Saha, M. C., Bouton, J. H., Lanzatella, C. L., and Tobias, C. M. (2014). Genetic linkage mapping and transmission ratio distortion in a three-generation four-founder population of Panicum virgatum (L.). G3 (Bethesda). 4, 913–923. doi: 10.1534/g3.113.010165
Li, X., Wang, X., Wei, Y., and Brummer, E. C. (2011). Prevalence of segregation distortion in diploid alfalfa and its implications for genetics and breeding applications. Theor. Appl. Genet. 123, 667–679. doi: 10.1007/s00122-011-1617-5
Lindholm, A. K., Dyer, K. A., Firman, R. C., Fishman, L., Forstmeier, W., Holman, L., et al. (2016). The ecology and evolutionary dynamics of meiotic drive. Trends Ecol. Evol. 31, 315–326. doi: 10.1016/j.tree.2016.02.001
Luro, F. L., Costantino, G., Terol, J., Argout, X., Allario, T., Wincker, P., et al. (2008). Transferability of the EST-SSRs developed on Nules clementine (Citrus clementina Hort ex Tan) to other Citrus species and their effectiveness for genetic mapping. BMC Genomics 9:287. doi: 10.1186/1471-2164-9-287
Mase, N., Sawamura, Y., Yamamoto, T., Takada, N., Nishio, S., Saito, T., et al. (2014). Direct genotyping of single pollen grains of a self-compatible mutant of Japanese pear (Pyrus pyrifolia) revealed inheritance of a duplicated chromosomal segment containing a second S-haplotype. Euphytica 200, 297–304. doi: 10.1007/s10681-014-1168-3
Morgan, C. H., Zhang, H., and Bomblies, K. (2017). Are the effects of elevated temperature on meiotic recombination and thermotolerance linked via the axis and synaptonemal complex? Philos. Trans. R. Soc. B Biol. Sci. 372:20160470. doi: 10.1098/rstb.2016.0470
Navarro, L., Aleza, P., Cuenca, J., Juárez, J., Pina, J. A., Ortega, C., et al. (2015). The mandarin triploid breeding program in Spain. Acta Hortic. 1065, 389–395. doi: 10.17660/ActaHortic.2015.1065.48
Nicolosi, E., Deng, Z. N., Gentile, A., La Malfa, S., Continella, G., and Tribulato, E. (2000). Citrus phylogeny and genetic origin of important species as investigated by molecular markers. Theor. Appl. Genet. 100, 1155–1166. doi: 10.1007/s001220051419
Ollitrault, F., Terol, J., Pina, J. A., Navarro, L., Talon, M., and Ollitrault, P. (2010). Development of SSR markers from Citrus clementina (Rutaceae) BAC end sequences and interspecific transferability in Citrus. Am. J. Bot. 97, e124–e129. doi: 10.3732/ajb.1000280
Ollitrault, P., Terol, J., Chen, C., Federici, C., Lotfy, S., Hippolyte, I., et al. (2012a). A reference genetic map of C. clementina hort. ex Tan.; citrus evolution inferences from comparative mapping. BMC Genomics 13:593. doi: 10.1186/1471-2164-13-593
Ollitrault, P., Terol, J., Garcia-Lor, A., Bérard, A., Chauveau, A., Froelicher, Y., et al. (2012b). SNP mining in C. clementina BAC end sequences; transferability in the Citrus genus (Rutaceae), phylogenetic inferences and perspectives for genetic mapping. BMC Genomics 13:13. doi: 10.1186/1471-2164-13-13
Perrier, X., and Jacquemound-Collet, J. P. (2018). DARwin 6.0.18. Available online at: http://darwin.cirad.fr (accessed March 4, 2020).
Rouiss, H., Cuenca, J., Navarro, L., Ollitrault, P., and Aleza, P. (2017a). Tetraploid citrus progenies arising from FDR and SDR unreduced pollen in 4x X 2x hybridizations. Tree Genet. Genomes 13:10. doi: 10.1007/s11295-016-1094-8
Rouiss, H., Cuenca, J., Navarro, L., Ollitrault, P., and Aleza, P. (2017b). Unreduced megagametophyte production in lemon occurs via three meiotic mechanisms, predominantly second-division restitution. Front. Plant Sci. 8:1211. doi: 10.3389/fpls.2017.01211
Seymour, D. K., Chae, E., Arioz, B. I., Koenig, D., and Weigel, D. (2019). Transmission ratio distortion is frequent in Arabidopsis thaliana controlled crosses. Heredity 122, 294–304. doi: 10.1038/s41437-018-0107-9
Sherman, J. D., and Stack, S. M. (1995). Two-dimensional spreads of synaptonemal complexes from solanaceous plants. VI. High-resolution recombination nodule map for tomato (Lycopersicon esculentum). Genetics 141, 683–708.
Shirasawa, K., Isobe, S., Hirakawa, H., Asamizu, E., Fukuoka, H., Just, D., et al. (2010). SNP discovery and linkage map construction in cultivated tomato. DNA Res. 17, 381–391. doi: 10.1093/dnares/dsq024
Sun, H., Rowan, B. A., Flood, P. J., Brandt, R., Fuss, J., Hancock, A. M., et al. (2019). Linked-read sequencing of gametes allows efficient genome-wide analysis of meiotic recombination. bioRxiv[Preprint]. doi: 10.1101/484022
Xu, X., Li, L., Dong, X., Jin, W., Melchinger, A. E., and Chen, S. (2013). Gametophytic and zygotic selection leads to segregation distortion through in vivo induction of a maternal haploid in maize. J. Exp. Bot. 64, 1083–1096. doi: 10.1093/jxb/ers393
Xu, Y., Zhu, L., Xiao, J., Huang, N., and McCouch, S. R. (1997). Chromosomal regions associated with segregation distortion of molecular markers in F 2, backcross, doubled haploid, and recombinant inbred populations in rice (Oryza sativa L.). Mol. Gen. Genet. 253, 535–545. doi: 10.1007/s004380050355
Zhou, W., Tang, Z., Hou, J., Hu, N., and Yin, T. (2015). Genetic map construction and detection of genetic loci underlying segregation distortion in an intraspecific cross of populus deltoides. PLoS One 10:e0126077. doi: 10.1371/journal.pone.0126077
Keywords: fluorescence-activated cell sorting, whole genome amplification, SSR and SNP markers, mandarin, lemon
Citation: Garavello M, Cuenca J, Dreissig S, Fuchs J, Navarro L, Houben A and Aleza P (2020) Analysis of Crossover Events and Allele Segregation Distortion in Interspecific Citrus Hybrids by Single Pollen Genotyping. Front. Plant Sci. 11:615. doi: 10.3389/fpls.2020.00615
Received: 02 March 2020; Accepted: 21 April 2020;
Published: 25 May 2020.
Edited by:Jaime Prohens, Universitat Politècnica de València, Spain
Reviewed by:Gaetano Distefano, University of Catania, Italy
Marco Caruso, Council for Agricultural and Economics Research (CREA), Italy
Copyright © 2020 Garavello, Cuenca, Dreissig, Fuchs, Navarro, Houben and Aleza. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Pablo Aleza, email@example.com