# Preferential Homologous Chromosome Pairing in a Tetraploid Intergeneric Somatic Hybrid (Citrus reticulata plus Poncirus trifoliata) Revealed by Molecular Marker Inheritance

Mourad Kamiri, Marc Stift, Gilles Costantino, Dominique Dambier, Tariq Kabbage, Patrick Ollitrault, Yann Froelicher

<figure><img src="image_1.png" /></figure>

► To cite this version:

Mourad Kamiri, Marc Stift, Gilles Costantino, Dominique Dambier, Tariq Kabbage, et al.. Preferential Homologous Chromosome Pairing in a Tetraploid Intergeneric Somatic Hybrid (Citrus reticulata plus Poncirus trifoliata) Revealed by Molecular Marker Inheritance. Frontiers in Plant Science, 2018, 9, 10.3389/fpls.2018.01557. [hal-02623132](https://hal.inrae.fr/hal-02623132)

HAL Id: hal-02623132

https://hal.inrae.fr/hal-02623132

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Check for
updates

# Preferential Homologous Chromosome Pairing in a Tetraploid Intergeneric Somatic Hybrid (Citrus reticulata + Poncirus trifoliata) Revealed by Molecular Marker Inheritance

Mourad Kamiri<fnref n="1" />, Marc Stift², Gilles Costantino<fnref n="3" />, Dominique Dambier⁴, Tariq Kabbage⁵, Patrick Ollitrault<fnref n="1" />* and Yann Froelicher<fnref n="1" />*

**OPEN ACCESS**

**Edited by:**

Sergio J. Ochatt,
INRA UMR1347 Agroécologie, France

**Reviewed by:**

Rosario Muleo,
Università degli Studi della Tuscia, Italy
Wen-Wu Guo,
Huazhong Agricultural University,
China

*Correspondence:

Patrick Ollitrault
patrick.ollitrault@cirad.fr

Yann Froelicher
yann.froelicher@cirad.fr

**Specialty section:**

This article was submitted to
Plant Breeding,
a section of the journal
Frontiers in Plant Science

Received: 13 July 2018

Accepted: 04 October 2018

Published: 02 November 2018

**Citation:**

Kamiri M, Stift M, Costantino G, Dambier D, Kabbage T, Ollitrault P and Froelicher Y (2018) Preferential Homologous Chromosome Pairing in a Tetraploid Intergeneric Somatic Hybrid (Citrus reticulata + Poncirus trifoliata) Revealed by Molecular Marker Inheritance.

Front. Plant Sci. 9:1557.
doi: 10.3389/fpls.2018.01557

The creation of intergeneric somatic hybrids between *Citrus* and *Poncirus* is an efficient approach for citrus rootstock breeding, offering the possibility of combining beneficial traits from both genera into novel rootstock lineages. These somatic hybrids are also used as parents for further tetraploid sexual breeding. In order to optimize these latter breeding schemes, it is essential to develop knowledge on the mode of inheritance in the intergeneric tetraploid hybrids. We assessed the meiotic behavior of an intergeneric tetraploid somatic hybrid resulting from symmetric protoplast fusion of diploid *Citrus reticulata* and diploid *Poncirus trifoliata*. The analysis was based on the segregation patterns of 16 SSR markers and 9 newly developed centromeric/pericentromeric SNP markers, representing all nine linkage groups of the *Citrus* genetic map. We found strong but incomplete preferential pairing between homologues of the same ancestral genome. The proportion of gametes that can be explained by random meiotic chromosome associations (τ) varied significantly between chromosomes, from 0.09 ± 0.02 to 0.47 ± 0.09, respectively, in chromosome 2 and 1. This intermediate inheritance between strict disomy and tetrasomy, with global preferential disomic tendency, resulted in a high level of intergeneric heterozygosity of the diploid gametes. Although limited, intergeneric recombinations occurred, whose observed rates, ranging from 0.09 to 0.29, respectively, in chromosome 2 and 1, were significantly correlated with τ. Such inheritance is of particular interest for rootstock breeding because a large part of the multi-trait value selected at the teraploid parent level is transmitted to the progeny, while the potential for some intergeneric recombination offers opportunities for generating plants with novel allelic combinations that can be targeted by selection.

**Keywords:** Citrus, somatic hybrid, tetraploid, disomic, tetrasomic, intermediate inheritance, SSR markers, SNP markers

<footnote n="1">¹ UMR AGAP, CIRAD, San Giuliano, France, ² Ecology, Department of Biology, University of Konstanz, Konstanz, Germany,</footnote>

<footnote n="3">³ UMR AGAP, INRA, San Giuliano, France, ⁴ UMR AGAP, CIRAD, Montpellier, France, ⁵ Domaines Abbes Kabbage, Agadir, Morocco</footnote>
INTRODUCTION

Polyploidization is a major evolutionary pathway (Ramsey
and Schemske, 2002; Adams and Wendel, 2005; Soltis et al.,
2009) with important implications for breeding. In plants,
polyploid lineages arise naturally through the union of
unreduced gametes (Harlan and deWet, 1975; Ramsey and
Schemske, 1998) and chromosome doubling of somatic cells
(Aleza et al., 2011). They can also be artificially produced by
treatment of cell tissue with colchicine (Aleza et al., 2009)
or somatic hybridization (Grosser et al., 2000). Natural
polyploidization is often associated with hybridization
between species, resulting in allopolyploidy. Allopolyploids
typically show strict preferential pairing in meiosis, resulting
in disomic inheritance, without an opportunity for interspecific
recombination (Soltis and Soltis, 2000; Pairon and Jacquemart,
2005). In allotetraploids, disomic inheritance leads diploid
gametes to display interspecific heterozygosity throughout
the genome. However, it has been recognized that even
allopolyploids that combine strongly diverged parental genomes
may have occasional nonhomologous chromosome pairing,
leading to nonstrict disomic inheritance and intergenomic
recombination (Stebbins, 1947; Sybenga, 1994, 1996; Udall
and Wendel, 2006; Stift et al., 2008; Jeridi et al., 2012).
Differential chromosome pairing affinities is therefore an
essential component determining species evolution in natural
polyploid populations. It is also a key parameter to design
efficient polyploid breeding strategies, as it strongly impacts
the inheritance of chromosome fragments and agronomical
traits. Molecular marker analysis is a powerful approach to study
the genetic structures of gametes produced by allopolyploid
organisms, and in turn to estimate the preferential pairing
pattern and its impacts on genome fragment inheritance and
recombination (Allendorf and Danzmann, 1997; Barone et al.,
2002; Jannoo et al., 2004; Bousalem et al., 2006; Landergott et al.,
2006).

In citrus, polyploidy has been described to improve adaptation
to different stresses and resilience (Mouhaya et al., 2010; Podda
et al., 2013; Allario et al., 2013; Tan et al., 2015; Tan et al.,
2017; Dutra de Souza et al., 2017; Oliveira et al., 2017; Oustric
et al., 2017, 2018), and several teams worldwide have been
developing rootstock breeding at the tetraploid level (Grosser
and Chandler, 2000; Grosser et al., 2000, 2015; Ollitrault et al.,
2000, 2007; Guo et al., 2007; Grosser and Gmitter, 2010; Dambier
et al., 2011; Guerra et al., 2016). Many of these programs
focus on combinations of *Citrus* species with *Poncirus*, a related
genus. *Poncirus* and *Citrus* species are sexually compatible
(Spiegel-Roy and Goldschmidt, 1996). However, molecular
phylogenetic studies based on whole genome resequencing data
have demonstrated substantial genetic differentiation (Carbonell-
Caballero et al., 2015; Wu et al., 2018). *Poncirus* has agriculturally
useful traits such as cold adaptation, tolerances to *Phytophthora*
species and nematodes, and resistance to the citrus tristeza
virus (CTV; Yang et al., 2001). It has also been described
to provide some tolerance to Huanglongbing, a devastating
citrus disease caused by the phloem bacterium *Candidatus
Liberibacter* sp. (Stover et al., 2010). However, *Poncirus* suffers

iron chlorosis on alkaline soils and it is susceptible to
salinity, which limits its use in some areas, particularly in
the Mediterranean Basin (Spiegel-Roy and Goldschmidt, 1996).
Rootstock breeding programs attempt to combine beneficial
Poncirus traits with abiotic stress tolerance traits present in
Citrus species. Although some interesting diploid intergeneric
sexual hybrids [e.g., citrange (P. trifoliata × C. sinensis) and
citrumelo (P. trifoliata × C. paradisi)] have been produced,
the long generation times, partial apomixes, and segregation of
beneficial allelic combinations due to the high heterozygosity
of parental genotypes (Herrero et al., 1996; Ollitrault et al.,
2003; Barkley et al., 2006; Wu et al., 2018) limit the efficiency
of conventional diploid breeding. Somatic hybridization is an
efficient alternative for citrus rootstock breeding as it allows
breeders to combine favorable parental genes regardless of their
heterozygosity level (Ollitrault et al., 1998; Grosser and Gmitter,
2010; Dambier et al., 2011; Ruiz et al., 2018). In citrus, somatic
hybridization has also been successfully applied to diversify the
tetraploid gene pool used as parent for triploid breeding (Grosser
et al., 2000; Ollitrault et al., 2007, 2008) and to generate hybrids
(Aleza et al., 2016b). The “Tetrazyg” approach (Grosser and
Gmitter, 2010) was introduced to breed new tetraploid rootstock
by sexual hybridization using selected allotetraploid somatic
hybrid rootstock as parent material. In order to optimize the
efficiency of the “Tetrazyg” strategy, it is essential to develop
knowledge on the mode of inheritance in the allotetraploid
hybrids.

Previous studies on tetraploid citrus revealed different meiotic
behavior. Weak or no preferential pairing of homologous
chromosomes (i.e., mostly tetrasomic inheritance) with
bivalents and quadrivalent meiotic configurations was
observed for a mandarin + lemon somatic hybrid (Kamiri
et al., 2011), for a somatic hybrid between tangelo and
pummelo (Xie et al., 2015), and doubled-diploid clementine
(Aleza et al., 2016a). By contrast, the doubled-diploid
“Mexican” lime had predominantly disomic segregation
(Rouiss et al., 2018). However, no data were available on the
inheritance mode of tetraploid *Citrus-Poncirus* intergeneric
hybrids.

The aim of the present study was to investigate the mode
of inheritance in a tetraploid intergeneric somatic hybrid called
Flhorag1 obtained through protoplast fusion between *Citrus
reticulata* cv "Willowleaf" mandarin and *Poncirus trifoliata* cv
"Pomeroy" (Ollitrault et al., 2000). This somatic hybrid provided
improved agronomic traits when used as rootstock with sweet
orange (Dambier et al., 2011). For this purpose, a triploid
progeny population (2n = 3x = 27) resulting from "Chandler"
pummelo × Flhorag1 sexual hybridization was genotyped at
19 simple sequence repeat (SSR) and 9 single nucleotide
polymorphism (SNP) loci to infer the allelic constitution of
gametes produced by the somatic hybrid. The likelihood-based
approaches proposed by Stift et al. (2008) for multi allelic
loci and Aleza et al. (2016a) for di-allelic loci in duplex
tetraploid were applied to analyze the meiotic behavior of
Flhorag1. The implications for citrus rootstock breeding are
discussed with a special focus on intergeneric heterozygosity
restitution.
# MATERIALS AND METHODS

## Plant Materials

An intergeneric somatic hybrid between a diploid (2n = 2x = 18) *Citrus reticulata* Blanco (“Willowleaf” mandarin SRA 133, hereafter called WLM) and a diploid (2n = 2x = 18) *Poncirus trifoliata* L. (“Pomeroy” trifoliate orange SRA 1074, hereafter called PON) was previously obtained by protoplast electrofusion (Ollitrault et al., 2000). To assess the inheritance of the tetraploid somatic hybrid (2n = 4x = 36; hereafter called Flhorag1), we performed a cross with diploid *Citrus maxima* (Burm.) Merr (2n = 2x = 18) (“Chandler” pummelo SRA 608, hereafter called CHA), using Flhorag1 as pollen donor. The pollen was collected on the mother tree of Flhorag1 (the plant directly regenerated from protoplast fusion).

CHA was chosen because it is self-incompatible, not apomictic, and it is genetically well differentiated from both WLM and PON. Cross was performed at the San Giuliano Research Station (Corsica, France). Recovered mature seeds were germinated *in vitro* in MT medium (Murashige and Tucker, 1969) supplemented with 30 g l⁻¹ sucrose and 1 mg l⁻¹ GA₃ (Ollitrault et al., 1996). The obtained plants were grafted on “Volkamer” lemon (*C. limonia* Osbeck) and further grown in a growth chamber.

## Flow Cytometry and Cytogenetic Analyses

The ploidy of the progeny was determined by flow cytometry and confirmed by chromosome counts. For flow cytometry, approximately 0.5 cm² of plantlet leaf and a similar amount of leaf tissue from a diploid reference *C. madurensis* (2n = 2x = 18) were chopped with a sharp razor blade in 250 µl of extraction buffer (Partec Cystain UV PreciseP) to isolate intact nuclei. After filtering the resulting suspension (30 µm pore size), 800 µl of DAPI (4-6-diamine-2-phenylindol) staining buffer was added (Partec, Cystain UV Precise P Staining Buffer). Samples were analyzed with a PA-I flow cytometer (Partec, Germany). We followed the protocol of D'Hont et al. (1996) for chromosome preparations. Briefly, fresh young leaf tissues were treated in 0.04% hydroxyquinoline for 4 h at room temperature and fixed for 48 h in 3:1 ethanol:acetic acid and stored at 4°C in 70% ethanol. The preparations were then treated for 20 min in 5N HCl and washed with distilled water. Finally, the tissue was deposited on microscope slides, stained with a drop of DAPI staining buffer (Partec, Cystain UV Precise P Staining Buffer), and squashed. Chromosomes were counted under an Eclipse 80i fluorescence microscope (Nikon Instruments, France).

## DNA Extraction

DNA was extracted from leaves using a modified mixed alkyl trimethyl ammonium bromide (MATAB) procedure (Gawel and Jarret, 1991). The DNA concentration was determined using the Hoechst 33258 (Sigma Chemical Co., MO, United States)

protocol (Sambrook and Russell, 2006). Samples were diluted with MQ sterile water and stored at -20°C until use.

## SSR Amplification

From previously identified, characterized, and mapped SSR loci (Froelicher et al., 2008; Luro et al., 2008; Ollitrault et al., 2010, 2012), we selected 19 loci (Table 1) that (i) were polymorphic between WLM and PON (diploid parents of the tetraploid somatic hybrid) and the seed parent CHA, and (ii) represented each of the nine linkage groups of the clementine genetic map (Ollitrault et al., 2012). Primers were labeled with WELLRED fluorochrome PA-2(dye2), PA-3(dye3), or PA-4(dye4) (Beckman-Coulter, CA, United States) and synthesized by Sigma-Aldrich (France). Among the selected markers, the mCrCIR02F12 locus was known to have a frequent null allele in *Poncirus* accessions (unpublished data), and the PON parent was suspected to be heterozygous (A0).

PCR reactions were performed in 20 µl with 1x Taq buffer, 1.5 mM MgCl₂, 0.8 U Taq DNA polymerase, 0.5 ng/µl template DNA, 0.2 mM dNTPs, 0.4 µM forward primer, and 0.4 µM reverse primer in an AG Primus 96 plus thermocycler (MWG, Germany). The PCR program consisted of 5 min initial denaturation (94°C), followed by 40 cycles of 30 s denaturation (94°C), 1 min primer annealing (50 or 55°C depending on the primers, Table 1), and 45 s extension (72°C), and a final extension of 4 min (72°C). Fragment analysis and allele calling were done by capillary electrophoresis using the CEQ 8800 genetic analyzer and software (Beckman Coulter, CA, United States). As CHA is totally differentiated from the WLM + PON tetraploid parent, the diploid gamete genotypes were directly inferred from the triploid hybrid genotypes by removing the specific CHA alleles.

## Centromeric SNP Development and Analysis

We have developed new SNP markers in centromeric region of the nine citrus chromosomes. “Pomeroy” trifoliate, “Willow leaf” mandarin and “Chandler” pummelo whole genome resequencing data (available, respectively, with SRX2442480, SRX372685, and SRX372688 SRA number in NCBI database) were mapped on the haploid clementine reference genome (Wu et al., 2014) using “BWA-MEM, v0.7.12-r1039” (Li and Durbin, 2010) and variant calling was performed with GATK (McKenna et al., 2010). To be able to infer the intergeneric gamete structure from the genotyping of the triploid issued from diploid Chandler diploid X (Willow Leaf + Pomeroy) tetraploid hybridizations, we selected SNPs homozygous for “Willow leaf” (AA), “Chandler” (AA), and “Pomeroy” (BB) in gene sequences of the identified centromeric/pericentromeric regions (Wu et al., 2014; Aleza et al., 2015). The genetic distances to the centromeres were inferred from available genetic mapping data (Ollitrault et al., 2012; Aleza et al., 2015) considering the flanking markers in the physical sequence (Wu et al., 2014). We were able to select efficient markers located at less than 1 cM for seven chromosomes, at less than 2 cM
<table><caption>**TABLE 1** | Characteristics of selected SSR markers used for the Willowleaf mandarin + Pomeroy Poncirus somatic hybrid (Flhorag1) allelic inheritance.</caption><thead><tr><th>Locus</th><th>EMBL</th><th>Linkage Group</th><th>Primer Sequence (5'-3')</th><th>Tm (C°)</th></tr></thead><tbody><tr><td>CiBE5055</td><td>ET111355</td><td>1</td><td>F:AACAGTGGTTCTGGGAAATAG<br>R:GGTGGTCTCAAAGTCATCATC</td><td>55</td></tr><tr><td>MEST431</td><td>FC883898</td><td>1</td><td>F: GAGCTCAAAACAATAGCCGC<br>R: CATACCTCCCCGTCCATCTA</td><td>55</td></tr><tr><td>mCrCIR02D09</td><td>FR677569</td><td>2</td><td>F:AATGATGAGGGTAAAGATG<br>R:ACCCATCACAAAACAGA</td><td>55</td></tr><tr><td>MEST46</td><td>FC901824</td><td>2</td><td>F: AACCAGAATCAGAACCCGA<br>R: GGTGAGCATCTGGACGACTT</td><td>55</td></tr><tr><td>mCrCIR02G12</td><td>FR677575</td><td>3</td><td>F:AAACCGAAATACAAGAGTG<br>R:TCCACAAACAATACAACG</td><td>55</td></tr><tr><td>mCrCIR02D04b</td><td>FR677564</td><td>4</td><td>F:CTCTCTTTCCCCATTAGA<br>R:AGCAAACCCCACAAC</td><td>50</td></tr><tr><td>mCrCIR03D12a</td><td>FR677577</td><td>4</td><td>F:GCCATAAGCCCTTTCT<br>R:CCCACAACCATCACC</td><td>50</td></tr><tr><td>mCrCIR07D06</td><td>FR677581</td><td>4</td><td>F:CCTTTTCACAGTTTGCTAT<br>R:TCAATTCCTCTAGTGTGTGT</td><td>55</td></tr><tr><td>mCrCIR03F05</td><td>FR692364</td><td>4</td><td>F:CTAAGGAAGAGTAGAGAGCA<br>R:TAAAATCCAAGGTTCCA</td><td>50</td></tr><tr><td>mCrCIR01F08a</td><td>AM489737</td><td>5</td><td>F:ATGAGCTAAAGAGAAGAGG<br>R:GGACTCAACACAACACAA</td><td>50</td></tr><tr><td>mCrCIR07E12</td><td>AM489750</td><td>5</td><td>F:TGTAGTCAAAAGCATCAC<br>R:TCTATGATTCCTGACTTTA</td><td>50</td></tr><tr><td>mCrCIR02F12</td><td>FR677570</td><td>6</td><td>F:GGCCATTTCTCTGATG<br>R:TAACTGAGGGATTGGGTTT</td><td>55</td></tr><tr><td>mCrCIR02D03</td><td>FR692360</td><td>7</td><td>F:CAGACAACAGAAAACCAA<br>R:GACCATTTTCCACTCAA</td><td>55</td></tr><tr><td>mCrCIR07E05</td><td>AM489749</td><td>7</td><td>F:GGAGAACAAAACACAATG<br>R:ATCTTTCGGACAATCTT</td><td>50</td></tr><tr><td>mCrCIR02A09</td><td>FR677568</td><td>8</td><td>F:ACAGAAGGTAGTATTTTAGGG<br>R:TTGTTTGGATGGGAAG</td><td>50</td></tr><tr><td>mCrCIR07B05</td><td>AM489747</td><td>8</td><td>F:TTTGTTCTTTTTGGTCTTTTT<br>R:CTTTTCTTTTCCTAGTTTCCC</td><td>50</td></tr><tr><td>mCrCIR02G02</td><td>FR677572</td><td>8</td><td>F:CAATAAGAAAACGCAGG<br>R:TGGTAGAGAAACAGAGGTG</td><td>55</td></tr><tr><td>mCrCIR07C09</td><td>AJ567410</td><td>9</td><td>F:GACCCTGCCTCCAAAGTATC<br>R:GTGGCTGTTGAGGGGTTG</td><td>55</td></tr><tr><td>mCrCIR02B07</td><td>AJ567403</td><td>9</td><td>F:CAGCTCAACATGAAAGG<br>R:TTGGAGAACAGGATGG</td><td>50</td></tr></tbody></table>

EMBL, accession number in the EMBL-EBI database (https://www.ebi.ac.uk/).

for chromosome 5 and less than 5cM for chromosome 8 (Table 2). The selected SNPs were analyzed with KASParTM Genotyping System (a competitive, allele-specific dual Förster resonance energy transfer-based assay). Primers were designed by LGC Genomics from the SNP locus flanking sequence (Supplementary Material 1). Detailed explanation on the specific conditions and reagents used can be found in Cuppen (2007). Identification of allele dosages in heterozygous triploid hybrids was carried out on the basis of relative allele signals, as described by Cuenca et al. (2013). For the nine SNP markers, the allelic configuration was as follow: “Chandler” (AA) × “Flhorag1” (AABB) producing AAA, AAB, and ABB triploid progenies with direct inference of the “Flhorag1” diploid gametes (AA, AB, and BB, respectively).

## Data Analysis

Pluri-allelic SSR markers: first we tested if the SSR allele frequencies observed in gametes deviated from those expected based on the parental genotypes (chi-square goodness-of-fit test), which would indicate the presence of null alleles or selection against particular alleles. Subsequently, we applied a likelihood based approach that simultaneously estimates two parameters: (1) $\tau$ – the proportion of gamete allelic constitutions that can be explained by random tetrasomic segregation and (2) $\beta$ – the double reduction frequency relative to $\tau$ (Stift et al., 2008). This was done for each preferential pairing scenario (e.g., for a marker displaying the ABCD genotype for WLM + PON, we evaluated three pairing scenarios: A with B, C with D; A with C, B with D; A with D, B
<table><caption>**TABLE 2** | Location of nine centromeric/pericentromeric SNP loci fully distinguishing “Pomeroy” poncirus from “willow leaf” mandarin and “Chandler” pummel.</caption><thead><tr><th>Marker</th><th>Scaffold</th><th>Position</th><th>SNP</th><th>Gene ID</th><th>D/Cent (cM)</th></tr></thead><tbody><tr><td>P1_16582061</td><td>1</td><td>16582061</td><td>[C/T]</td><td>Ciclev10007229m.g</td><td>&lt;1</td></tr><tr><td>P2_19903193</td><td>2</td><td>19903193</td><td>[T/C]</td><td>Ciclev10014147m.g</td><td>&lt;1</td></tr><tr><td>P3_16287238</td><td>3</td><td>16287238</td><td>[T/C]</td><td>Ciclev10024301m.g</td><td>&lt;1</td></tr><tr><td>P4_9505439</td><td>4</td><td>9505439</td><td>[A/G]</td><td>Ciclev10031824m.g</td><td>&lt;1</td></tr><tr><td>P5_20142568</td><td>5</td><td>20142568</td><td>[T/A]</td><td>Ciclev10000575m.g</td><td>&lt;2</td></tr><tr><td>P6_3130249</td><td>6</td><td>3130249</td><td>[G/A]</td><td>Ciclev10011214m.g</td><td>&lt;1</td></tr><tr><td>P7_15458045</td><td>7</td><td>15458045</td><td>[A/G]</td><td>Ciclev10025280m.g</td><td>&lt;1</td></tr><tr><td>P8_16631925</td><td>8</td><td>16631925</td><td>[A/T]</td><td>Ciclev10030352m.g</td><td>&lt;5</td></tr><tr><td>P9_12062066</td><td>9</td><td>12062066</td><td>[C/T]</td><td>Ciclev10006749m.g</td><td>&lt;1</td></tr></tbody></table>

SNP positions, gene ID, and flanking position came from the haploid clementine reference sequence (Wu et al., 2014; https://phytozome.jgi.doe.gov/pz/portal.htm#!info?alias=Org_Cclementina).

with C). For parameter estimation, we used the constrained non-linear regression (CNLR) function implemented in SPSS 15.0 (SPSS syntax file<fnref n="1" />). We then used a likelihood ratio test (LRT) evaluated against a compound distribution of 1/2χ₀² + 1/2χ₁² (Self and Liang, 1987) to determine if the model with the estimated τ explained the data significantly better than a random null model (i.e., a model with strict tetrasomic segregation, τ = 1). For loci for which “Flhorag1” had less than four alleles, we considered the possibility of null alleles using a G-test for contingency tables to evaluate the most probable constitution. In case of three alleles, we compared abc0 and abcc (for parents of the somatic hybrid genotyped ab and cc), and in the case of two allele aabb, aab0, abb0, or ab00 (parents of the somatic hybrid aa and bb).

Di-allelic centromeric/pericentromeric SNP markers: for centromeric diallelic markers the segregation model for duplex tetraploid parents (AABB configuration) is greatly simplified because the double reduction effect can be missed and the segregation pattern is a direct function of τ. Therefore, for these markers τ was estimated by the maximum likelihood approach proposed by Aleza et al. (2016a).

Genetic dissimilarities between diploid gametes were estimated with the DARwin 6.0 software (Perrier and Jacquemoud-Collet, 2006) using the “simple matching” dissimilarity index using the pluriallelic SSRs markers to integrate intergeneric and intraspecific polymorphisms.

# RESULTS

## Ploidy Level Determination

The cross between the diploid maternal parent (CHA) and the tetraploid somatic hybrid (“Flhorag1”) resulted in 63 fully developed seeds collected from 41 fruits, from which 59 plantlets germinated. Flow cytometry indicated that the plantlets were triploid, with one tetraploid exception (data not shown).

Chromosome counts (Supplementary Material 2) confirmed the flow cytometry results. The segregation analyses were performed for the 58 triploid hybrids.

## SSR Genotyping

For 15 out of the 19 SSRs markers, allele frequencies observed in the gametes did not significantly differ from expected frequencies based on the “Flhorag1” genotype, that is, there was no segregation distortion. This confirmed the correct determination of allele doses in the somatic hybrid and the absence of null alleles. For one marker (mCrCIR02F12), a scenario assuming the previously suspected presence of a null allele in the PON genome could explain the observed allele frequencies. Moreover, the “microsatellite DNA allele counting—peak ratios” (MAC-PR) method based on relative areas of the peaks of the different alleles provided additional evidence of a null allele at this locus.

At the three remaining loci (mCrCIR02D03, mCrCIR07C09, and mCrCIR03F05), the progeny allele frequencies were significantly distorted (P < 0.05) and remained so when scenarios were considered involving null alleles in either WLM or PON (Table 3). In addition, MAC-PR provided no evidence for null alleles at these loci. They were excluded from further analyses, so the segregation analyses were performed with 16 loci with Mendelian segregation. For these 16 SSR markers, the MAC-PR analysis of “Flhorag1” confirmed a full addition of parental (WLM and PON) alleles.

## SSR Inheritance in the Flhorag1 Tetraploid Somatic Hybrid

For all 16 analyzed loci, the estimated value of τ (proportion of gametes that can be explained by random meiotic chromosome associations) ranged from 0.07 to 0.58 (Figure 1 and Table 4). For each marker, fits between full tetrasomic inheritance model (τ = 1) and the best fitting intermediate model were compared. The LRT values were all highly significant, with *p* values ranging from 7.56E-09 to 1.78E-02 (Table 4). Preferential pairing was always between homologous chromosomes (between chromosomes derived from the same parental genome; i.e., mandarin or Poncirus).

The double reduction rate could only be estimated at 10 loci. Indeed, for MEST046, MEST431, mCrCIR01F08a, mCrCIR07B05, and mCrCIR07E05 markers, WLM and PON were both homozygous with a double homoduplex conformation (AA, BB). For the mCrCIR02F12 locus, null alleles prevented double reduction gamete identification. Double reduction was detected for five loci (mCrCIR02D09, mCrCIR02G12, mCrCIR03D12a, mCrCIR02A09, and mCrCIR 02B07) with *β* estimates ranging from 0.03 to 0.07 (Table 4).

## SNP Inheritance in the Flhorag1 Tetraploid Somatic Hybrid

For each SNP locus, “Flhorag1” was heterozygous with equilibrated allelic doses (AABB). No allelic distortion was

<footnote n="1">¹http://www.genetics.org/cgi/content/full/genetics.107.085027/DC1</footnote>
<table><caption>TABLE 3 | Distribution of inferred diploid gamete genotypes that produced CHA × (WLM + PON) progeny at the 19 studied loci.</caption><thead><tr><th rowspan="2">Locus</th><th rowspan="2">Linkage Group</th><th colspan="2">Parental Genotypes</th><th colspan="5">Inferred Genotypes of Diploid Gametes</th><th rowspan="2">Allelic Distortion p Value</th></tr><tr><th>WLM</th><th>PON</th><th>AA</th><th>AB</th><th>AC</th><th>BB</th><th>CC</th><th>BC</th></tr></thead><tbody><tr><td colspan="10"><strong>Segregation with two or three different alleles</strong></td></tr><tr><td>CiBE5055</td><td>1</td><td>BC</td><td>AA</td><td>5</td><td>22</td><td>28</td><td>0</td><td>0</td><td>3</td><td>0.68</td></tr><tr><td>MEST431</td><td>1</td><td>BB</td><td>AA</td><td>5</td><td>46</td><td>n.a.</td><td>7</td><td>n.a.</td><td>n.a.</td><td>0.68</td></tr><tr><td>mCrCIR02D09</td><td>2</td><td>AB</td><td>CC</td><td>1</td><td>0</td><td>25</td><td>1</td><td>1</td><td>30</td><td>0.79</td></tr><tr><td>MEST46</td><td>2</td><td>BB</td><td>AA</td><td>2</td><td>56</td><td>n.a.</td><td>0</td><td>n.a.</td><td>n.a.</td><td>0.71</td></tr><tr><td>mCrCIR02G12</td><td>3</td><td>CC</td><td>AB</td><td>2</td><td>3</td><td>25</td><td>1</td><td>0</td><td>27</td><td>0.54</td></tr><tr><td>mCrCIR02D04b</td><td>4</td><td>AA</td><td>BC</td><td>5</td><td>28</td><td>24</td><td>0</td><td>0</td><td>1</td><td>0.66</td></tr><tr><td>mCrCIR03D12a</td><td>4</td><td>AA</td><td>BC</td><td>1</td><td>28</td><td>24</td><td>1</td><td>1</td><td>3</td><td>0.66</td></tr><tr><td>mCrCIR07D06</td><td>4</td><td>AB</td><td>CC</td><td>0</td><td>0</td><td>26</td><td>0</td><td>6</td><td>26</td><td>0.54</td></tr><tr><td>mCrCIR03F05</td><td>4</td><td>BB</td><td>AC</td><td>2</td><td>27</td><td>8</td><td>0</td><td>4</td><td>17</td><td>0.02*</td></tr><tr><td>mCrCIR01F08a</td><td>5</td><td>BB</td><td>AA</td><td>3</td><td>50</td><td>n.a.</td><td>5</td><td>n.a.</td><td>n.a.</td><td>0.71</td></tr><tr><td>mCrCIR07E12</td><td>5</td><td>AC</td><td>BB</td><td>0</td><td>26</td><td>5</td><td>2</td><td>0</td><td>25</td><td>0.79</td></tr><tr><td>mCrCIR02D03</td><td>7</td><td>AB</td><td>AC</td><td>1</td><td>5</td><td>28</td><td>1</td><td>1</td><td>22</td><td>1.14E-06*</td></tr><tr><td>mCrCIR07E05</td><td>7</td><td>BB</td><td>AA</td><td>4</td><td>51</td><td>n.a.</td><td>3</td><td>n.a.</td><td>n.a.</td><td>0.85</td></tr><tr><td>mCrCIR02A09</td><td>8</td><td>CC</td><td>AB</td><td>1</td><td>1</td><td>28</td><td>0</td><td>2</td><td>26</td><td>0.87</td></tr><tr><td>mCrCIR07B05</td><td>8</td><td>BB</td><td>AA</td><td>1</td><td>55</td><td>n.a.</td><td>2</td><td>n.a.</td><td>n.a.</td><td>0.85</td></tr><tr><td>mCrCIR07C09</td><td>9</td><td>AC</td><td>BC</td><td>1</td><td>20</td><td>7</td><td>1</td><td>11</td><td>18</td><td>0.04*</td></tr><tr><td>mCrCIR02B07</td><td>9</td><td>AC</td><td>BB</td><td>0</td><td>25</td><td>3</td><td>5</td><td>2</td><td>23</td><td>0.97</td></tr><tr><td colspan="10"><strong>Segregation with four different alleles</strong></td></tr><tr><td>mCrCIR02G02</td><td>8</td><td>AB</td><td>CD</td><td>0</td><td>2</td><td>23</td><td>7</td><td>0</td><td>13</td><td>0.14</td></tr><tr><td colspan="10"><strong>Segregation including a null allele</strong></td></tr><tr><td>mCrCIR02F12</td><td>6</td><td>AB</td><td>C0</td><td>0</td><td>3</td><td>14</td><td>10</td><td>0</td><td>21</td><td>0.16</td></tr></tbody></table>

WLM, *Citrus reticulata cv Willowleaf mandarin*; PON, *Poncirus trifoliata cv Pomeroy*; *Fhorag1*, *Citrus reticulata cv Willowleaf mandarin + Poncirus trifoliata cv Pomeroy*; A. B. C. D, allele length (A represents the longest allele); 0, null allele. For three of these loci, allelic distortion was observed (*). Corresponding P values for allelic distortion are given.

observed for the nine diallelic SNPs. They revealed similar intergeneric heterozygosity restitution and τ values than the ones estimated on the same LGs with SSRs (Table 5). Intergeneric heterozygosity varied between 86% in LGs 1 and 9 to 97% in LG 2. In accordance, the lower τ value (0.11) was inferred for LG2 and the higher (0.42) for LG1 and LG9.

## Distribution of τ Variations of SSR and SNP Loci Among the Different LGs

Loci belonging to the same linkage group (LG) tended to have similar τ values with significant differences between LGs. LG1 (CiBE5055, MEST431, P1_16582061; τ = 0.47 ± 0.09); LG2 (mCrCIR02D09, MEST046, P2_19903193; τ = 0.09 ± 0.02); LG3 (mCrCIR02G12, P3_16287238; τ = 0.21 ± 0.00); LG4 (mCrCIR02D04b, mCrCIR03D12a, mCrCIR07D06, P4_9505439; τ = 0.29 ± 0.03); LG5 (mCrCIR01F08a, mCrCIR07E12, P5_20142568; τ = 0.36 ± 0.05); LG6 (mCrCIR02F12; P6_3130249; τ = 0.18 ± 0.04); LG7 (mCrCIR07E05; P7_15458045; τ = 0.34 ± 0.03); LG8 (mCrCIR02A09, mCrCIR07B05, mCrCIR02G02, P8_16631925; τ = 0.15 ± 0.01); LG9 (mCrCIR02B07, P9_12062066; τ = 0.40 ± 0.02). Therefore, we can consider that LG1 has an intermediary inheritance while the other LGs had a preferential disomic tendency highly marked for LG2 and LG8.

## Heterozygosity Transmission and Gamete Diversity

Diploid gametes from this tetraploid intergeneric somatic hybrid transmitted both intraspecific and intergeneric heterozygosity. The intergeneric heterozygosity transmission rate was directly linked with τ, and decreased linearly from 97 to 79% with increasing τ values (Figure 2). When analyzing the multilocus gametic structure over the 25 loci with complete allelic differentiation between *Citrus* and *Poncirus*, we found five diploid gametes (among the 58 analyzed) with complete intergeneric heterozygosity (Supplementary Material 3). The average intergeneric and intragenic heterozygosity transmission were, respectively, 90.3 and 4.1%. Intergeneric recombinations were detected from multiloci pattern within each LG (Supplementary Material 3), whose frequencies (0.29, 0.09, 0.14, 0.24, 0.26, 0.10, 0.16, 0.14, 0.26 in LGs 1 to 9, respectively) were significantly related ($R^2$ = 0.82) with the τ average values by LGs. We compared the distribution of dissimilarities between gametes considering all SSR alleles (intraspecific and intergeneric diversity) or solely the intergeneric origin of alleles (i.e., only two alleles considered: a poncirus allele P and a mandarin allele M and therefore segregation between PP, MM, and PM genotypes). Preferential pairing and consecutive predominant intergeneric heterozygosity transmission resulted in a relatively low intergeneric diversity
<figure><img src="image_4.png" /><figcaption>**FIGURE 1** | Example of deviance (G) for six observed SSR loci segregations in the "Flhorag1" somatic hybrid gametes with inheritance models ranging from τ = 0 (full disomic) to τ = 1 (full tetrasomic). Among the 16 analyzed SSRs markers Mest431 and mCrCIR02D09 display, respectively, the higher (0.58) and lower (0.07) τ values.</figcaption></figure>

contribution and gamete genetic diversity seemed to be mainly due to parental intraspecific diversity segregation (Figure 3).

DISCUSSION

Intermediate Inheritance With a
Preferential Disomic Trend Detected in
the *C. reticulata* + *P. trifoliata* (*Flhorag1*)
Intergeneric Tetraploid Somatic Hybrid

All 16 SSR and 9 SNP markers studied displayed intermediate
inheritance with a preferential disomic tendency. The gamete
proportions explained by random meiotic chromosome
association (τ) ranged from 7 to 58%, and for all markers
preferential pairing occurred between chromosomes derived
from the same parent (i.e., WLM or PON). Such intermediate
inheritance with a disomic tendency was also described
for a doubled-diploid "Mexican lime" (Rouiss et al., 2018)
belonging to *C. x aurantiifolia* (Christm.) Swingle, a species of
interspecific origin (*C. micrantha* Wester *x* *C. medica* L.; Nicolosi
et al., 2000; Curk et al., 2016; Wu et al., 2018). The levels of
preferential pairing revealed for "Flhorag1" and the doubled-
diploid "Mexican lime" are much higher than that observed
for (i) an interspecific somatic hybrid (*C. reticulata* + *C.
lemon*) within the *Citrus* genus, where random tetrasomic
inheritance (τ = 1) could only be rejected for 8 out of 17
markers, with τ ranging from 0.24 to 0.95 (Kamiri, et al. 2011),
and (ii) a doubled-diploid clementine displaying tetrasomic
segregation for five chromosomes, intermediate segregation with
a tetrasomic tendency for three chromosomes and intermediate
segregation with a disomic tendency for only one chromosome

(Aleza et al., 2016a). Preferential pairing variations between
tetraploid citrus reveal different levels of divergence between
their constitutive genomes.

Cytogenetic study of meiosis behavior is another approach to
analyze chromosome affinity in polyploids. Indeed, tetravalent
formation testifies for chromosome homology between parents
while they should be precluded in wide allotetraploids (assuming
the absence of large structural variation). Few cytogenetic
studies were performed in tetraploid somatic hybrids and most
concerned interspecific hybrids within the *Citrus* genus. Del
Bosco et al. (1999) were the first to report the microsporogenesis
in a Citrus interspecific tetraploid somatic hybrid between
“Valencia” sweet orange and “Femminello” lemon. They observed
frequent tetravalents. One of these tetravalents by meiotic cell
was related to a reciprocal translocation in sweet-orange but
the additional tetravalents were considered as a consequence
of intergenomic pairing. Similar observation for frequent
tetravalents were reported for “Hamlin” sweet orange + “Rough”
Lemon and “Key” lime + “Valencia” sweet orange (Chen et al.,
2004) and “Willow leaf” mandarin + “Eureka” lemon (Kamiri
et al., 2011). For the last somatic hybrid, this cytogenetic
observation was associated with preferential tetrasomic
inheritance of molecular markers. However, as previously
mentioned, reciprocal or inverted translocations can result in
tetravalent formation even in diploid parents as observed in
“Valencia” sweet orange (Del Bosco et al., 1999) or Mexican
lime (Rouiss et al., 2018). Therefore, cytogenetic observations
should be associated with marker segregation analysis for a
full understanding of tetraploid meiosis behavior. In the case
of mandarin + poncirus somatic hybrids, the predominant
preferential pairing between homologous chromosomes inferred
in “Flhorag1” from molecular marker segregation data is in
agreement with the previous cytogenetic study conducted in an
intergeneric citrus somatic hybrid between “Cleopatra” mandarin
and “Argentine” poncirus (unpublished results mentioned in
Chen et al., 2004). Indeed, they observed a high percentage of
bivalents, suggesting low chromosome homology between the
fusion parents.

Recent whole nuclear genome resequencing data revealed an
important differentiation of *Poncirus trifoliata* from the *Citrus*
species clade (Wu et al., 2018). This genomic differentiation
level is in line with the preferential disomic inheritance observed
for our *C. reticulata* + *P. trifoliata* intergeneric somatic hybrid.
Interestingly, our finding that preferential pairing was not
complete (i.e., inheritance is not fully disomic and intergeneric
recombinations are observed in our study), indicated that the
remaining homology is still enough to allow *Citrus* and *Poncirus*
chromosome pairing, and helps understand the long-known
intergeneric sexual compatibility at the diploid level (Cameron
and Garber, 1968). Moreover, comparative studies of *Citrus* and
*Poncirus* genetic maps revealed a high level of synteny and
collinearity of markers. From unsaturated maps, Chen et al.
(2008) observed only a few inversions between shared loci.
Bernet et al. (2010) reported high collinearity between Fortune
“mandarin” and *Poncirus trifoliata*. More recently, the availability
of reference whole genome sequence allowed to compare genetic
and physical maps and globally confirmed the good conservation
of marker order between *P. trifoliata* and “Sunki” mandarin (Curtolo et al., 2018) and *P. trifoliata* and sweet orange (Huang et al., 2018).

Although inheritance was consistently intermediate between disomic and tetrasomic for all markers, there was considerable variation in the degree of preferential pairing. Such variation among loci is common, and has for example been observed in autotetraploid Pacific oysters (Curole and Hedgecock, 2005) and sugarcane (Jannoo et al., 2004). Variation among markers could merely reflect stochasticity, but also true differences in homology for different parts of the genome. Significant differences of the τ values were observed between LGs. LG1 displayed the higher value of random association (close to 50%) while the other LGs had a preferential disomic tendency. Disomy was very high for LG2 and LG8 with more than 85% of preferential chromosome pairing. Although an analysis of more markers would be needed

to confirm the pattern, our results suggest that the divergence between “Willow leaf” mandarin and “Pomeroy” poncirus varies between chromosomes.

Double Reduction Rates and Heterozygosity Transmission

Gametes resulting from double reduction have been detected in five loci. Double reduction implies multivalent formation and a crossover between the considered locus and its centromere with further adjacent segregation (Haynes and Douches, 1993). It results in increased homozygosity and its maximum frequency (1/6) is reached in case of systematic quadrivalent formation at meiosis (Stift et al., 2008). Our progeny sample size was not enough for accurate estimation of the double reduction rate but highlighted the possibility of tetravalent formation in the intergeneric somatic hybrid.

<table><caption>TABLE 4 | Fitting the inheritance model and intraparental and interparental heterozygosity rates on the segregation of 16 SSR loci in the "Flhorag1" diploid gametes.</caption><thead><tr><th rowspan="2">SSR</th><th rowspan="2">LG</th><th rowspan="2">WLM</th><th rowspan="2">PON</th><th rowspan="2">Pref Pairing</th><th colspan="4">Best model</th><th colspan="2">% Heterozygosity</th></tr><tr><th>β</th><th>τ</th><th>LRT</th><th>p value</th><th>Intra-specific</th><th>Inter-generic</th></tr></thead><tbody><tr><td>CIBE5055</td><td>1</td><td>BC</td><td>AA</td><td>BC/AA</td><td>–</td><td>0.41</td><td>5.79</td><td>8.04E-03</td><td>5%</td><td>86%</td></tr><tr><td>MEST431</td><td>1</td><td>AA</td><td>BB</td><td>AA/BB</td><td>–</td><td>0.58</td><td>4.41</td><td>1.78E-02</td><td>0%</td><td>79%</td></tr><tr><td>mCrCIR02D09</td><td>2</td><td>AB</td><td>CC</td><td>AB/CC</td><td>0.07</td><td>0.07</td><td>20.03</td><td>4.00E-06</td><td>2%</td><td>95%</td></tr><tr><td>MEST46</td><td>2</td><td>BB</td><td>AA</td><td>AA/BB</td><td>–</td><td>0.1</td><td>32.04</td><td>7.56E-09</td><td>0%</td><td>97%</td></tr><tr><td>mCrCIR02G12</td><td>3</td><td>CC</td><td>AB</td><td>CC/AB</td><td>0.05</td><td>0.21</td><td>12.74</td><td>1.79E-04</td><td>5%</td><td>90%</td></tr><tr><td>mCrCIR02D04b</td><td>4</td><td>AA</td><td>BC</td><td>AA/BC</td><td>–</td><td>0.31</td><td>8.39</td><td>1.89E-03</td><td>2%</td><td>90%</td></tr><tr><td>mCrCIR03D12a</td><td>4</td><td>AA</td><td>BC</td><td>AA/BC</td><td>0.04</td><td>0.23</td><td>11.42</td><td>3.64E-04</td><td>5%</td><td>90%</td></tr><tr><td>mCrCIR07D06</td><td>4</td><td>AB</td><td>CC</td><td>AB/CC</td><td>–</td><td>0.31</td><td>8.39</td><td>1.89E-03</td><td>0%</td><td>90%</td></tr><tr><td>mCrCIR01F08a</td><td>5</td><td>BB</td><td>AA</td><td>AA/BB</td><td>–</td><td>0.41</td><td>11.45</td><td>3.57E-04</td><td>0%</td><td>86%</td></tr><tr><td>mCrCIR07E12</td><td>5</td><td>AC</td><td>BB</td><td>AC/BB</td><td>–</td><td>0.36</td><td>7.01</td><td>4.06E-03</td><td>9%</td><td>88%</td></tr><tr><td>mCrCIR02F12</td><td>6</td><td>AB</td><td>C0</td><td>AB/C0</td><td>–</td><td>0.21</td><td>11.73</td><td>3.07E-04</td><td>7%</td><td>93%</td></tr><tr><td>mCrCIR07E05</td><td>7</td><td>BB</td><td>AA</td><td>AA/BB</td><td>–</td><td>0.36</td><td>13.86</td><td>9.90E-05</td><td>0%</td><td>88%</td></tr><tr><td>mCrCIR02A09</td><td>8</td><td>CC</td><td>AB</td><td>CC/AB</td><td>0.03</td><td>0.14</td><td>14.19</td><td>8.20E-05</td><td>2%</td><td>93%</td></tr><tr><td>mCrCIR07B05</td><td>8</td><td>BB</td><td>AA</td><td>AA/BB</td><td>–</td><td>0.16</td><td>27.26</td><td>8.88E-08</td><td>0%</td><td>95%</td></tr><tr><td>mCrCIR02G02</td><td>8</td><td>AB</td><td>CD</td><td>AB/CD</td><td>–</td><td>0.16</td><td>13.79</td><td>1.02E-04</td><td>5%</td><td>95%</td></tr><tr><td>mCrCIR02B07</td><td>9</td><td>AC</td><td>BB</td><td>AC/BB</td><td>0.07</td><td>0.38</td><td>6.34</td><td>5.89E-03</td><td>5%</td><td>83%</td></tr></tbody></table>

Comparisons of fit were performed between the tetrasomic null model and the best fitting intermediate model. LRT values were evaluated as described in the Materials and methods and were significant at the indicated level. WLM, Citrus reticulata cv Willowleaf mandarin; PON, Poncirus trifoliata cv Pomeroy; A. B. C. D, allele length (A represents the longest allele); 0, null allele.

<table><caption>TABLE 5 | Fitting the inheritance model and intergeneric heterozygosity rates on the segregation of nine centromeric SNP loci in the "Flhorag1" diploid gametes.</caption><thead><tr><th rowspan="2"></th><th colspan="4">Parental Genotypes</th><th colspan="3">Diploid Gametes<br>Genotypes</th><th rowspan="2">Allelic Distortion p value</th><th rowspan="2">τ</th><th rowspan="2">Intergeneric Heterozygosity</th></tr><tr><th>LG</th><th>WLM</th><th>PON</th><th>Flhorag1</th><th>AA</th><th>AB</th><th>BB</th></tr></thead><tbody><tr><td>P1_16582061</td><td>1</td><td>AA</td><td>BB</td><td>AABB</td><td>5</td><td>50</td><td>3</td><td>0.71</td><td>0.42</td><td>86%</td></tr><tr><td>P2_19903193</td><td>2</td><td>AA</td><td>BB</td><td>AABB</td><td>1</td><td>56</td><td>1</td><td>1.00</td><td>0.11</td><td>97%</td></tr><tr><td>P3_16287238</td><td>3</td><td>AA</td><td>BB</td><td>AABB</td><td>3</td><td>54</td><td>1</td><td>0.71</td><td>0.21</td><td>93%</td></tr><tr><td>P4_9505439</td><td>4</td><td>AA</td><td>BB</td><td>AABB</td><td>4</td><td>52</td><td>2</td><td>0.71</td><td>0.31</td><td>90%</td></tr><tr><td>P5_20142568</td><td>5</td><td>AA</td><td>BB</td><td>AABB</td><td>2</td><td>52</td><td>4</td><td>0.71</td><td>0.31</td><td>90%</td></tr><tr><td>P6_3130249</td><td>6</td><td>AA</td><td>BB</td><td>AABB</td><td>2</td><td>55</td><td>1</td><td>0.85</td><td>0.16</td><td>95%</td></tr><tr><td>P7_15458045</td><td>7</td><td>AA</td><td>BB</td><td>AABB</td><td>2</td><td>52</td><td>4</td><td>0.71</td><td>0.31</td><td>90%</td></tr><tr><td>P8_16631925</td><td>8</td><td>AA</td><td>BB</td><td>AABB</td><td>1</td><td>55</td><td>2</td><td>0.85</td><td>0.16</td><td>95%</td></tr><tr><td>P9_12062066</td><td>9</td><td>AA</td><td>BB</td><td>AABB</td><td>5</td><td>50</td><td>3</td><td>0.71</td><td>0.42</td><td>86%</td></tr></tbody></table>
Inter-generic heterozygosity

<figure><img src="image_5.png" /><figcaption>**FIGURE 2** | Correlation between the intergeneric heterozygosity transmission rate and the τ value for the 16 SSR and 9 SNP markers.</figcaption></figure>

Frequency

<figure><img src="image_6.png" /><figcaption>**FIGURE 3** | Distribution of genetic dissimilarities between diploid gametes of the Flhorag1 somatic hybrid based on intergeneric segregation (blue) and all allele segregation (red).</figcaption></figure>

The diploid gametes produced by tetraploid somatic hybrids are highly heterozygous. Under tetrasomic inheritance, they should transmit either intraparental or interparental heterozygosity (Ollitrault et al., 2008), in contrast to strict disomic inheritance, which leads to exclusive transmission of interparental heterozygosity for wide interspecific somatic hybrids. When inheritance is intermediate, as found in the intergeneric *Citrus-Poncirus* Flhorag1 somatic hybrid, transmission of intraspecific and intergeneric heterozygosity depends on the degree of preferential chromosome pairing and the double reduction rates in case of tetravalent formation. As no double reduction is possible in centromeric areas, they should display higher intergeneric restitution values.

Assuming preferential pairing between parental homologous chromosomes, interparental heterozygosity transmission decreases with increasing τ values. Therefore, it is logical that the interparental heterozygosity transmission (in this case intergeneric heterozygosity) observed in the intergeneric *Flhorag1* somatic hybrid (90%) was much higher than the 64.1% observed for a *C. reticulata* + *C. lemon* interspecific somatic hybrid (Kamiri et al., 2011). As a consequence of this high level of intergeneric heterozygosity inheritance, the genetic diversity of the gamete population results mainly from parental intraspecific diversity segregation. However, these conclusions for genetic markers should not be extrapolated for phenotypic traits because phenotypic differentiation at the intergeneric level is very high compared with the intra-poncirus or mandarin variability.

## Implications for Breeding

Citrus breeding is hampered by complex genome structures, features of reproductive biology and the long juvenile phase, but can take advantage of vegetative propagation including apomictic seeds for rootstock multiplication, allowing clonal propagation of elite genotypes whatever its genome complexity (Ollitrault and Navarro, 2012). Therefore, breeding strategies are generally based on one (or very few) cycle of variability induction followed by direct selection of cultivars or rootstock. In such context, it is essential to optimize the transfer to the progenies, of the genetic gains obtained by phenotypic selection at the parental level. In case of *Citrus-Poncirus* intergeneric polyploid breeding, most of the genetic value comes from the combination of favorable traits from *Citrus* (tolerance to abiotic stresses such as salinity, water deficit, calcareous soils) and *Poncirus* (resistance/tolerance to diseases and pests such as tristeza virus, *Phytophthora*, nematodes; cold tolerance). It is based on dominant inheritance in highly intergeneric heterozygous structures. For further breeding using these intergeneric tetraploids as sexual parents (tetrazyg strategy; Grosser and Gmitter, 2010), it is therefore essential to transmit a large part of the parental intergeneric heterozygosity to the progenies in order to prevent overall breakage of the favorable complex multilocus genotypic structure selected at the allotetraploid parental level. From this work, it appears that the differentiation between *C. reticulata* and *Poncirus trifoliata* genomes results in preferential homologous pairing and predominantly intergeneric heterozygosity transmission. Most of the value of the somatic hybrid should thus be transmitted to its progeny. Moreover, the infrequent occurrence of non-homologous chromosome pairing offers an opportunity for intergeneric recombination, and the generation of novel allelic combinations. Several *Citrus* × *Poncirus* diploid intergeneric hybrids such as citrumello (*C. paradisi* × *P. trifoliata*), citrange (*C. sinensis* × *P. trifoliata*), and citrandarin (*C. reticulata* × *P. trifoliata*) proved their interest as rootstock and are widely used worldwide. Taking advantage of spontaneous chromosome doubling in nuclear cells, doubled diploid lines were selected for most of them (Aleza et al., 2011). If they display similar preferential disomic tendency than
the “Flhorag1” somatic hybrid, they should produce gametes in which a very high proportion of the genome of the initial diploid intergeneric hybrid has been transferred. Thus, the high phenotypic value deeply selected at the diploid intergeneric parent level should have a very marked positive impact on the products of the “Tetrazyg” strategy when using these doubled-diploid parents.

## AUTHOR CONTRIBUTIONS

PO, YF, and TK designed the experiment. DD created and provided the parental somatic hybrids the parental somatic hybrid. MK and GC performed the SSRs analysis. MK and MS analyzed the results. MK, MS, YF, and PO wrote the paper.

## REFERENCES

Adams, K. L., and Wendel, J. F. (2005). Polyploidy and genome evolution in plants. *Curr. Opin. Plant Biol.* 8, 135-141. doi: 10.1016/j.pbi.2005.01.001

Aleza, P., Cuenca, J., Hernandez, M., Juarez, 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., Cuenca, J., Juárez, J., Navarro, L., and Ollitrault, P. (2016a). Inheritance in doubled-diploid clementine and comparative study with SDR unreduced gametes of diploid clementine. *Plant Cell Rep.* 35, 1573-1586. doi: 10.1007/s00299-016-1972-4

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., Garcia-Lor, A., and Juárez, J. (2016b). Recovery of citrus cybrid plants with diverse mitochondrial and chloroplastic genome combinations by protoplast fusion followed by in vitro shoot, root, or embryo micrografting. *Plant Cell Tissue Organ Cult.* 126, 205-217. doi: 10.1007/s11240-016-0991-8

Aleza, P., Juárez, J., Ollitrault, P., and Navarro, L. (2009). Production of tetraploid plants of non-apomictic citrus genotypes. *Plant Cell Rep.* 28, 1837-1846. doi: 10.1007/s00299-009-0783-2

Allario, T., Brumos, J., Colmenero-Flores, J. M., Iglesias, D. J., Pina, J. A., Navarro, L., et al. (2013). Tetraploid Rangpur lime rootstock increases drought tolerance via enhanced constitutive root abscisic acid production. *Plant Cell Environ.* 36, 856-868. doi: 10.1111/pce.12021

Allendorf, F. W., and Danzmann, R. G. (1997). Secondary tetrasomic segregation of MDH-B and preferential pairing of homeologues in rainbow trout. *Genetics* 145, 1083-1092.

Barkley, N. A., Roosevelt, M. L., Krueger, R. R., and Federici, C. T. (2006). Assessing genetic diversity and population structure in a citrus germplasm collection utilizing simple sequence repeat markers (SSRs). *Theor. Appl. Genet.* 112, 1519-1531. doi: 10.1007/s00122-006-0255-9

Barone, A., Li, J., Sebastiano, A., Cardi, T., and Frusciante, L. (2002). Evidence for tetrasomic inheritance in a tetraploid *Solanum commersonii* (+) *S. tuberosum* somatic hybrid through the use of molecular markers. *Theor. Appl. Genet.* 104, 539-546. doi: 10.1007/s00122-001-0792-1

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

Bousalem, M., Arnau, G., Hochu, I., Arnolin, R., Viader, V., Santoni, S., et al. (2006). Microsatellite segregation analysis and cytogenetic evidence for tetrasomic inheritance in the American yam *Dioscorea trifida* and a new basic chromosome number in the *Dioscoreae*. *Theor. Appl. Genet.* 113, 439-451. doi: 10.1007/s00122-006-0309-z

## FUNDING

The Domaines Abbes Kabbage, Agadir (Morocco), and the Comité Mixte Inter Universitaire Franco-Marocain (Program Volubils, MA/08/196) funded a grant awarded to MK. This study was also supported by the Collectivité Territoriale de Corse (CTC), “Sélection de Porte-Greffe pour le Développement d’une Agrumiculture Durable” ARR-15-05332 and ARR-17001001.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.01557/full#supplementary-material

Cameron, J. W., and Garber, M. J. (1968). Identical-twin hybrids of *Citrus x Poncirus* from strictly sexual seed parents. *Am. J. Bot.* 55, 199-205. doi: 10.1002/j.1537-2197.1968.tb06962.x

Carbonell-Caballero, J., Alonso, R., Ibañez, V., Terol, J., Talon, M., and Dopazo, J. A. (2015). Phylogenetic analysis of 34 chloroplast genomes elucidates the relationships between wild and domestic species within the genus *Citrus*. *Mol. Biol. Evol.* 32, 2015-2035. doi: 10.1093/molbev/msv082

Chen, C. L., Guo, W. W., Yi, H. L., and Deng, X. X. (2004). Cytogenetic analysis of two interspecific citrus allotetraploid somatic hybrids and their diploid fusion parents. *Plant Breed.* 123, 332-337. doi: 10.1111/j.1439-0523.2004.00972.x

Chen, X., 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

Cuenca, J., Aleza, P., Navarro, L., and Ollitrault, P. (2013). Assignment of SNP allelic configuration in polyploids using competitive allele-specific PCR: application to citrus triploid progeny. *Ann. Bot.* 111, 731-742. doi: 10.1093/aob/mct032

Cuppen, E. (2007). Genotyping by allele-specific amplification (KASPar). *CSH Protoc.* 2007:4841. doi: 10.1101/pdb.prot4841

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

Curole, J. P., and Hedgecock, D. (2005). Estimation of preferential pairing rates in second-generation autotetraploid pacific oysters (*Crassostrea gigas*). *Genetics* 171, 855-859. doi: 10.1534/genetics.105.043042

Curtolo, M., Soratto, T. A. T., Gazaffi, R., Takita, M. A., Machado, M. A., and Cristofani-Yaly, M. (2018). High-density linkage maps for *Citrus sunki* and *Poncirus trifoliata* using DArTseq markers. *Tree Genet. Genomes* 14:5. doi: 10.1007/s11295-017-1218-9

Dambier, D., Benyahia, H., Pensabene-Bellavia, G., Aka, K., Froelicher, Y., Belfalah, Z., et al. (2011). Somatic hybridization for Citrus rootstock breeding: an effective tool to solve some important issues of the Mediterranean citrus industry. *Plant Cell Rep.* 30, 883-900. doi: 10.1007/s00299-010-1000-z

Del Bosco, S. F., Tusa, N., and Conicella, C. (1999). Microsporogenesis in a Citrus interspecific tetraploid somatic hybrid and its fusion parents. *Heredity* 83, 373-377. doi: 10.1038/sj.hdy.6885770

D'Hont, A., Grivet, L., Feldmann, P., Rao, S., Berding, N., and Glaszmann, J. C. (1996). Characterisation of the double genome structure of modern sugarcane cultivars (*Saccharum spp.*) by molecular cytogenetics. *Mol. Gen. Genet.* 250, 405-413. doi: 10.1007/BF02174028

Dutra de Souza, J., de Andrade Silva, E. M., Coelho Filho, M. A., Morillon, R., Bonatto, D., Micheli, F., et al. (2017). Different adaptation strategies of two citrus scion/rootstock combinations in response to drought stress. *PLoS One* 12:e0177993. doi: 10.1371/journal.pone.0177993

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
Gawel, N. J., and Jarret, R. L. (1991). A modified CTAB DNA extraction procedure for *Musa* and *Ipomoea*. *Plant Mol. Biol. Rep.* 9, 374-380. doi: 10.1007/BF02672076

Grosser, J., and Gmitter, F. (2010). Protoplast fusion for production of tetraploids and triploids: applications for scion and rootstock breeding in citrus. *Plant Cell Tissue Organ Cult.* 104, 343-357. doi: 10.1007/s11240-010-9823-4

Grosser, J. W., Barthe, G. A., Castle, B., Gmitter, F. G. Jr., and Lee, O. (2015). The development of improved tetraploid citrus rootstocks to facilitate advanced production systems and sustainable citriculture in Florida. *Acta Hortic.* 1065, 319-327. doi: 10.17660/ActaHortic.2015.1065.38

Grosser, J. W., and Chandler, J. L. (2000). Somatic hybridization of high yield, cold-hardy and disease resistant parents for citrus rootstock improvement. *J. Hortic. Sci. Biotechnol.* 75, 641-644. doi: 10.1080/14620316.2000.11511300

Grosser, J. W., Ollitrault, P., and Olivares-Fuster, O. (2000). Somatic hybridization in citrus: an effective tool to facilitate variety improvement. *In Vitro Cell. Dev. Biol. Plant* 36, 434-439. doi: 10.1007/s11627-000-0080-9

Guerra, D., Schifino-Wittmann, M. T., Schwarz, S. F., Weiler, R. L., Dahmer, N., and Dutra Souza, P. V. (2016). Tetraploidization in citrus rootstocks: effect of genetic constitution and environment in chromosome duplication. *Crop Breed. Appl. Biotechnol.* 16, 35-41. doi: 10.1590/1984-70332016v16n1a6

Guo, W. W., Wu, R. C., Cheng, Y. J., and Deng, X. X. (2007). Production and molecular characterization of *Citrus* intergeneric somatic hybrids between red tangerine and citrange. *Plant Breed.* 126, 72-76. doi: 10.1111/j.1439-0523.2006.01315.x

Harlan, J. R., and deWet, J. M. (1975). On O. Winge and a prayer: the origins of polyploidy. *Bot. Rev.* 41, 361-390. doi: 10.1007/BF02860830

Haynes, K. G., and Douches, D. S. (1993). Estimation of the coefficient of double reduction in the cultivated tetraploid potato. *Theor. Appl. Genet.* 85, 857-862. doi: 10.1007/BF00225029

Herrero, R., Asins, M. J., Pina, J. A., Carbonell, E. A., and Navarro, L. (1996). Genetic diversity in the orange subfamily Aurantoideae. II. Genetic relationships among genera and species. *Theor. Appl. Genet.* 93, 1327-1334. doi: 10.1007/BF00223466

Huang, M., Moose, M., Yu, Q., Du, D., Zhang, Y., Deng, Z., et al. (2018). Construction of high-density genetic maps and detection of QTLs associated with Huangongbl eing infection in citrus. *bioRxiv* [Preprint]. doi: 10.1101/330753

Jannoo, N., Grivet, L., David, J., D'Hont, A., and Glaszmann, J. C. (2004). Differential chromosome pairing affinities at meiosis in polyploid sugarcane revealed by molecular markers. *Heredity* 93, 460-467. doi: 10.1038/sj.hdy.6800524

Jeridi, M., Perrier, X., Rodier-Goud, M., Ferchichi, A., D'Hont, A., and Bakry, F. (2012). Cytogenetic evidence of mixed disomic and polysomic inheritance in an allotetraploid *AABB* *Musa* genotype. *Ann. Bot.* 1108, 1593-1606. doi: 10.1093/aob/mcs220

Kamiri, M., Stift, M., Srairi, I., Costantino, G., El Moussadik, A., Hmyene, A., et al. (2011). Evidence for non-disomic inheritance in a Citrus interspecific tetraploid somatic between *C. reticulata* and *C. lemon* hybrid using SSR markers and cytogenetic analysis. *Plant Cell Rep.* 30, 1415-1425. doi: 10.1007/s00299-011-1050-x

Landergott, U., Naciri, Y., Schneller, J. J., and Holderegger, R. (2006). Allelic configuration and polysomic inheritance of highly variable microsatellites in tetraploid gynodioecious *Thymus praecox* agg. *Theor. Appl. Genet.* 113, 453-465. doi: 10.1007/s00122-006-0310-6

Li, H., and Durbin, R. (2010). Fast and accurate long-read alignment with Burrows-Wheeler transform. *Bioinformatics* 26, 589-595. doi: 10.1093/bioinformatics/btp698

Luro, F., 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

McKenna, A., Hanna, M., Banks, E., Sivachenko, A., Cibulskis, K., Kernytsky, A., et al. (2010). The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. *Genome Res.* 20, 1297-1303. doi: 10.1101/gr.107524.110

Mouhaya, W., Allario, T., Brumos, J., Andrés, F., Froelicher, Y., Luro, F., et al. (2010). Sensitivity to high salinity in tetraploid citrus seedlings increases with water availability and correlates with expression of candidate genes. *Funct. Plant Biol.* 37, 674-685. doi: 10.1071/FP10035

Murashige, T., and Tucker, D. P. H. (1969). Growth factor requirements of tissue culture. *Proc. First Int. Citrus Symp.* 3, 1155-1161.

Nicolosi, E., Deng, Z. N., Gentile, A., Malfa La, 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

Oliveira, T. M., Ben Yahmed, J., Dutra, J., Maserti, B. E., Talon, M., Navarro, L., et al. (2017). Better tolerance to water deficit in doubled diploid 'Carrizo citrange' compared to diploid seedlings is associated with more limited water consumption and better H₂O₂ scavenging. *Acta Physiol. Plant.* 39:e204. doi: 10.1007/s11738-017-2497-3

Ollitrault, F., Terol, J., Pina, J. A., Navarro, L., Talon, M., and Ollitrault, P. (2010). Development of SSRs markers from *Citrus clementina* (Rutaceae) BAC end sequences and interspecific transferability in *Citrus*. *Am. J. Bot.* 97, 124-129. doi: 10.3732/ajb.1000280

Ollitrault, P., Dambier, D., Allent, V., Luro, F., and Jacquemond, C. (1996). "In vitro rescue and selection of spontaneous triploids by flow cytometry for easy peeler Citrus breeding," in *Proceedings of the Eighth International Citrus Congress International Society of Citriculture*, Vol. 2, Sun City, 913-917.

Ollitrault, P., Dambier, D., Froelicher, Y., Bakry, F., and Aubert, B. (1998). Rootstock breeding strategies for the Mediterranean citrus industry: the somatic hybridization potential. *Fruits* 53, 335-344.

Ollitrault, P., Dambier, D., Luro, F., and Froelicher, Y. (2008). "Ploidy manipulation for breeding seedless triploid citrus," in *Plant Breeding Reviews*, Vol. 30, ed. J. Janick (Hoboken, NJ: J. Wiley), 323-352.

Ollitrault, P., Froelicher, Y., Dambier, D., and Seker, M. (2000). "Rootstock breeding by somatic hybridization for the Mediterranean citrus industry," in *Proceedings of the First International Symposium on Citrus Biotechnology*, Eilat, 157-162.

Ollitrault, P., Guo, W., and Grosser, J. W. (2007). "Somatic hybridization," in *Citrus Genetics, Breeding and Biotechnology*, ed. I. A. Khan (Wallingford: CAB International), 235-260. doi: 10.1079/9780851990194.0235

Ollitrault, P., Jacquemond, C., Dubois, C., and Luro, F. (2003). "Citrus," in *Genetic Diversity of Cultivated Tropical Plants*, eds P. Hamon, M. Seguin, X. Perrier, and J.-C. Glaszmann (Montpellier: CIRAD), 193-217.

Ollitrault, P., and Navarro, L. (2012). "Citrus," in *Fruit Breeding*, eds M. L. Badenes and D. H. Byrne (New York, NY: Springer), 623-662. doi: 10.1007/978-1-4419-0763-9_16

Ollitrault, P., Terol, J. F., Chen, C., Federici, C. T., Lotfy, S., Hippolyte, I., et al. (2012). 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

Ouatric, J., Morillon, R., Luro, F., Herbette, S., Lourkistia, R., Giannettini, J., et al. (2017). Tetraploid Carrizo citrange rootstock (*Citrus sinensis* Osb. × *Poncirus trifoliata* L. Raf.) enhances natural chilling stress tolerance of common clementine (*Citrus clementina* Hort. ex Tan). *J. Plant Physiol.* 214, 108-115. doi: 10.3389/fpls.2016.00701

Ouatric, J., Morillon, R., Ollitrault, P., Herbette, S., Luro, F., Froelicher, Y., et al. (2018). Somatic hybridization between diploid *Poncirus* and *Citrus* improves natural chilling and light stress tolerances compared with equivalent doubled-diploid genotypes. *Trees* 32, 883-895. doi: 10.1007/s00468-018-1682-3

Pairon, M. C., and Jacquemart, A. L. (2005). Disomic segregation of microsatellites in the tetraploid *Prunus serotina* Ehrh. (Rosaceae). *J. Am. Soc. Hortic. Sci.* 130, 729-734.

Perrier, X., and Jacquemoud-Collet, J. P. (2006). *DARwin Software*. Available at: http://darwin.cirad.fr/darwin

Podda, A., Checcucci, G., Mouhaya, W., Centeno, D., Rofidal, V., Del Carratore, R., et al. (2013). Salt-stress induced changes in the leaf proteome of diploid and tetraploid mandarins with contrasting Na(+) and Cl(-) accumulation behaviour. *J. Plant Physiol.* 170, 1101-1112. doi: 10.1016/j.jplph.2013.03.006

Ramsey, J., and Schemske, D. W. (1998). Pathways, mechanisms, and rates of polyploid formation in flowering plants. *Annu. Rev. Ecol. Syst.* 29, 467-501. doi: 10.1104/pp.16.00176

Ramsey, J., and Schemske, D. W. (2002). Neopolyploidy in flowering plants. *Annu. Rev. Ecol. Syst.* 33, 589-639. doi: 10.1146/annurev.ecolsys.33.010802.150437

Rouiss, H., Bakry, F., Froelicher, Y., Navarro, L., Aleza, P., and Ollitrault, P. (2018). Origin of *C. latifolia* and *C. aurantiifolia* triploid limes: the preferential disomic
inheritance of doubled-diploid 'Mexican' lime is consistent with an interploid
hybridization hypothesis. *Ann. Bot.* 121, 571–585. doi: 10.1093/aob/mcx179

Ruiz, M., Pensabene-Bellavia, G., Quinones, A., Garcia-Lor, A., Morillon, R., Ollitrault, P., et al. (2018). Molecular characterization and stress tolerance evaluation of new allotetraploid somatic hybrids between carrizo citrange and *Citrus macrophylla* W. rootstocks. *Front. Plant Sci.* 9:901. doi: 10.3389/fpls.2018.00901

Sambrook, J., and Russell, D. W. (eds). (2006). "Fluorometric quantitation of DNA using hoechst 33258," in *The Condensed Protocols from Molecular Cloning: a Laboratory Manual*, (Cold Spring Harbor, NY: Cold Spring Harbor Laboratory Press), 742–745.

Self, S. G., and Liang, K. Y. (1987). Asymptotic properties of maximum likelihood estimators and likelihood ratio tests under nonstandard conditions. *J. Am. Stat. Assoc.* 82, 605–610. doi: 10.1080/01621459.1987.10478472

Soltis, D. E., Albert, V. A., Leebens-Mack, J., Bell, C. D., Paterson, A. H., Zheng, C., et al. (2009). Polyploidy and angiosperm diversification. *J. Bot.* 96, 336–348. doi: 10.3732/ajb.0800079

Soltis, P. S., and Soltis, D. E. (2000). The role of genetic and genomic attributes in the success of polyploids. *Proc. Natl. Acad. Sci. U.S.A.* 97, 7051–7057. doi: 10.1073/pnas.97.13.7051

Spiegel-Roy, P., and Goldschmidt, E. E. (1996). *Biology of Citrus*. Cambridge: Cambridge University Press. doi: 10.1017/CBO9780511600548

Stebbins, G. (1947). Types of polyploids: their classification and significance. *Adv. Genet.* 1, 403–429. doi: 10.1016/S0065-2660(08)60490-3

Stift, M., Berenos, C., Kuperus, P., and Van Tienderen, P. (2008). Segregation models for disomic, tetrasomic and intermediate inheritance in tetraploids: a general procedure applied to Rorippa (yellow cress) microsatellite data. *Genetics* 179, 2113–2123. doi: 10.1534/genetics.107.085027

Stover, E., Shatters, R., McCollum, G., Hall, D. G., and Duan, Y. (2010). Evaluation of *Candidatus* Liberibacter asiaticus titer in field-infected trifoliate cultivars: preliminary evidence for HLB resistance. *Proc. Fla. State Hortic. Soc.* 123, 115–117.

Sybenga, J. (1994). Preferential pairing estimates from multivalent frequencies in tetraploids. *Genome* 37, 1045–1055. doi: 10.1139/g94-149

Sybenga, J. (1996). Chromosome pairing affinity and quadrivalent formation in polyploids: do segmental allopolyploids exist? *Genome* 39, 1176–1184.

Tan, F. Q., Tu, H., Liang, W. J., Long, J. M., Wu, X. M., Zhang, H. Y., et al. (2015). Comparative metabolic and transcriptional analysis of a doubled diploid

and its diploid citrus rootstock (C. junos cv. Ziyang xiangcheng) suggests its
potential value for stress resistance improvement. BMC Plant Biol. 15:89. doi:
10.1186/s12870-015-0450-4

Tan, F. Q., Tu, H., Wang, R., Wu, X. M., Xie, K. D., Chen, J. J., et al. (2017). Metabolic adaptation following genome doubling in citrus doubled diploids revealed by non-targeted metabolomics. *Metabolomics* 13:143. doi: 10.1007/s11306-017-1276-x

Udall, J. A., and Wendel, J. F. (2006). Polyploidy and crop improvement. *Crop Sci.* 46, 3–14. doi: 10.2135/cropsci2006.07.0489tpg

Wu, G. A., Prochnik, S., Jenkins, J., Salse, J., Hellsten, U., Murat, F., et al. (2014). Sequencing of diverse mandarin, pummelo and orange genomes reveals complex history of admixture during citrus domestication. *Nat. Biotechnol.* 2014, 656–662. doi: 10.1038/nbt.2906

Wu, G. A., Terol, J. F., Ibáñez, V., López-García, A., Pérez-Román, E., Borredá, C., et al. (2018). Genomics of the origin and evolution of Citrus. *Nature* 554, 311–316. doi: 10.1038/nature25447

Xie, K. D., Xia, Q. M., Wang, X. P., Liang, W. J., Wu, X. M., Grosser, J. W., et al. (2015). Cytogenetic and SSR-marker evidence of mixed disomic, tetrasomic, and intermediate inheritance in a citrus allotetraploid somatic hybrid between 'Nova' tangelo and 'HB' pummelo. *Tree Genet. Genomes* 11:112. doi: 10.1007/s11295-015-0940-4

Yang, Z. N., Ye, X. R., Choi, S., Molina, J., Moonan, F., Wing, R. A., et al. (2001). Construction of a 1.2-Mb contig including the citrus tristeza virus resistance gene locus using a bacterial artificial chromosome library of *Poncirus trifoliata* (L.) Raf. *Genome* 44, 382–393. doi: 10.1139/g01-021

**Conflict of Interest Statement:** 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.

Copyright © 2018 Kamiri, Stift, Costantino, Dambier, Kabbage, Ollitrault and Froelicher. 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.