# Meiotic Behaviors of Allotetraploid Citrus Drive the Interspecific Recombination Landscape, the Genetic Structures, and Traits Inheritance in Tetrazyg Progenies Aiming to Select New Rootstocks

Lény Calvez, Alexis Dereeper, Aude Perdereau, Pierre Mournet, Maëva

Miranda, Saturnin Bruyère, Barbara Hufnagel, Yann Froelicher, Arnaud

Lemainque, Raphaël Morillon, et al.

► To cite this version:

Lény Calvez, Alexis Dereeper, Aude Perdereau, Pierre Mournet, Maëva Miranda, et al. Meiotic Behaviors of Allotetraploid Citrus Drive the Interspecific Recombination Landscape, the Genetic Structures, and Traits Inheritance in Tetrazyg Progenies Aiming to Select New Rootstocks. Plants, 2023, 12 (8), pp.1630. 10.3390/plants12081630. [https://doi.org/10.3390/plants12081630](https://doi.org/10.3390/plants12081630)

HAL Id: hal-04109398

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

Submitted on 30 May 2023

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_1.png" /></figure>

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

Article

# Meiotic Behaviors of Allotetraploid Citrus Drive the Interspecific Recombination Landscape, the Genetic Structures, and Traits Inheritance in Tetrazyg Progenies Aiming to Select New Rootstocks

Lény Calvez<fnref n="1" />,<fnref n="2" /> ID, Alexis Dereeper <fnref n="1" />,<fnref n="2" />, Aude Perdereau <fnref n="3" /> ID, Pierre Mournet <fnref n="2" />,<fnref n="4" /> ID, Maëva Miranda <fnref n="2" />,<fnref n="4" />, Saturnin Bruyère <fnref n="1" />,<fnref n="2" />, Barbara Hufnagel <fnref n="1" />,<fnref n="2" />, Yann Froelicher <fnref n="2" />,<fnref n="5" />, Arnaud Lemainque <fnref n="3" />, Raphaël Morillon <fnref n="2" />,<fnref n="4" /> and Patrick Ollitrault <fnref n="2" />,<fnref n="4" />,* ID

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

check for
updates

**Citation:** Calvez, L.; Dereeper, A.; Perdereau, A.; Mournet, P.; Miranda, M.; Bruyère, S.; Hufnagel, B.; Froelicher, Y.; Lemainque, A.; Morillon, R.; et al. Meiotic Behaviors of Allotetraploid Citrus Drive the Interspecific Recombination Landscape, the Genetic Structures, and Traits Inheritance in Tetrazyg Progenies Aiming to Select New Rootstocks. *Plants* 2023, *12*, 1630.
https://doi.org/10.3390/plants12081630

Academic Editor: Sonja Siljak-Yakovlev

Received: 25 February 2023

Revised: 29 March 2023

Accepted: 6 April 2023

Published: 12 April 2023

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

**Copyright:** © 2023 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** Sexual breeding at the tetraploid level is a promising strategy for rootstock breeding in citrus. Due to the interspecific origin of most of the conventional diploid citrus rootstocks that produced the tetraploid germplasm, the optimization of this strategy requires better knowledge of the meiotic behavior of the tetraploid parents. This work used Genotyping By Sequencing (GBS) data from 103 tetraploid hybrids to study the meiotic behavior and generate a high-density recombination landscape for their tetraploid intergenic Swingle citrumelo and interspecific Volkamer lemon progenitors. A genetic association study was performed with root architecture traits. For citrumelo, high preferential chromosome pairing was revealed and led to an intermediate inheritance with a disomic tendency. Meiosis in Volkamer lemon was more complex than that of citrumelo, with mixed segregation patterns from disomy to tetrasomy. The preferential pairing resulted in low interspecific recombination levels and high interspecific heterozygosity transmission by the diploid gametes. This meiotic behavior affected the efficiency of Quantitative Trait Loci (QTL) detection. Nevertheless, it enabled a high transmission of disease and pest resistance candidate genes from *P. trifoliata* that are heterozygous in the citrumelo progenitor. The tetrazyg strategy, using doubled diploids of interspecific origin as parents, appears to be efficient in transferring the dominant traits selected at the parental level to the tetraploid progenies.

**Keywords:** polyploidy; rootstock breeding; meiosis; preferential pairing; recombination landscape; quantitative traits loci; citrus; *P. trifoliata*; *C. medica*; inter-specificity

## 1. Introduction

Citrus cultivation is one of the major fruit productions in the world. However, this production is confronted with numerous biotic and abiotic constraints, including drought, extreme temperature, salinity, and pest and disease. These stresses devastate orchards by severely influencing tree growth and development, fruit production, and quality. Rootstocks are widely used for citrus cultivation because they improve trees' performance and aim at tolerance or even resistance to various stresses. Citrus breeding programs are essential to create new rootstock varieties better adapted to current global issues, such as

<footnote n="1">¹ UMR AGAP, CIRAD, F-97170 Petit-Bourg, France; leny.calvez@cirad.fr (L.C.); alexis.dereeper@ird.fr (A.D.); saturnin.bruyere@cirad.fr (S.B.); barbara.hufnagel@cirad.fr (B.H.)</footnote>

<footnote n="2">² UMR AGAP, Institut Agro, CIRAD, INRAE, University of Montpellier, F-34060 Montpellier, France; pierre.mournet@cirad.fr (P.M.); maeva.miranda@cirad.fr (M.M.); yann.froelicher@cirad.fr (Y.F.); raphael.morillon@cirad.fr (R.M.)</footnote>

<footnote n="3">³ Genoscope, Institut de Biologie François-Jacob, Commissariat à l'Energie Atomique (CEA), Université Paris-Saclay, F-91000 Evry, France; aperdere@genoscope.cns.fr (A.P.)</footnote>

<footnote n="4">⁴ UMR AGAP, CIRAD, F-34398 Montpellier, France</footnote>

<footnote n="5">⁵ UMR AGAP, CIRAD, F-20230 San Giuliano, France</footnote>

<footnote>* Correspondence: patrick.ollitrault@cirad.fr</footnote>
environmental pressures exacerbated by global warming [1] and the emergence of certain
diseases such as Huanglongbing (HLB). Due to a sizeable sexual compatibility within the
"true citrus" gene pool as defined by [2], many natural or man-driven interspecific and
even intergeneric hybridizations resulted in interesting rootstocks [3]. Volkamer lemon
(Citrus limonia Osb.) is one of the parental rootstocks used in the present work. It is an
example of a natural interspecific F1 hybrid between a mandarin (C. reticulata Blanco) and
a citron (C. medica L.) [4,5]. Used as rootstock, Volkamer lemon confers strong growth
and good fruit yield. It is adapted to dry, calcareous, and saline soil and presents good
tolerances against Citrus tristeza virus (CTV), exocortis viroid, and Phytophthora spp. [6].

Conventional citrus breeding by sexual hybridizations has been practiced since the end of the 19th century and is still commonly used today. Among the citrus relatively sexually compatible with the *Citrus* genus [7], *Poncirus trifoliata* (L.) Raf. has good agronomic qualities that make it a crucial genetic resource in citrus rootstock improvement programs. Indeed, *P. trifoliata* is resistant to nematodes [8] and CTV, which is the most critical citrus viral disease [9]. *P. trifoliata* is highly tolerant to *Phytophthora* spp. [10] and cold temperatures [11]. In addition, it has been described as tolerant to HLB [12–17], but this observation is still controversial [13]. However, *P. trifoliata* is susceptible to alkaline soils and salinity, which limits its use in many citrus-growing regions [7,18]. As a seedling, *P. trifoliata* genotypes are quite tolerant to water deficit [19]. Intergeneric hybridizations between *Poncirus* and *Citrus* genera allow the combination of several traits. Intergeneric hybridization (*Poncirus* × *Citrus*) seems to be an efficient way to develop new rootstock genotypes responding to the new challenges of the citrus industry. New intergeneric citrus hybrids such as Swingle citrumelo (C. paradisi Macf. × *P. trifoliata* (L.) Raf.), citrandarin (C. reticulata × *P. trifoliata*) Carrizo, Troyer, and C35 citrange (C. sinensis (L.) Osb. × *P. trifoliata*) have good agronomic properties and are widely used around the world [6]. The Swingle citrumelo rootstock, one of the genotypes used in this work, is resistant to nematodes, *phytophthora*, and CTV and also confers a good cold tolerance [6].

Most of the citrus germplasm and worldwide-used citrus rootstocks are diploid.
However, during the last decade, several studies have shown that tetraploid citrus rootstock,
compared to their diploid counterpart, showed increased tolerance to different pedoclimatic
constraints such as salt stress [20–23], low temperature [24], nutrient deficiency [25,26],
and water deficit [27–29]. Sivager et al. (2022) showed that tetraploid Swingle citrumelo
rootstocks improve HLB tolerance in diploid and triploid grafted scions. Interestingly, the
better tolerance in polyploid root and the scion is related to maintaining the phloem sap
flux by limiting the callose deposition [30] and a capacity to cope more effectively with the
oxidative stress induced by the disease at the leaf level [30] and the root level [31]. As a
consequence of this finding, ploidy manipulations have become an integral aspect of citrus
breeding programs [32].

Polyploids can result from sexual or somatic polyploidization events. Citrus sexual
polyploidization mainly results in the formation of triploids through the union of unre-
duced gametes with normal haploid gametes [33–35]. Due to the seedlessness of most
triploid hybrids, sexual polyploidization is actively exploited for seedless variety breed-
ing [36] but not for rootstock breeding. Somatic polyploidization is frequently observed
in apomictic citrus species, with a chromosomal doubling of nucellar tissues resulting in
spontaneous doubled diploid genotypes [37,38]. Doubled diploids can also be artificially
induced by chemical treatment with colchicine or oryzalin [39–41]. Several teams have
widely used somatic hybridization to create tetraploid rootstocks combining characters of
interest from different genetic resources. This method allows researchers to add genomes
of two diploid genotypes without sexual recombination [42–49] and, therefore, to combine
favorable dominant traits of the two parents. Recent genomic studies revealed chromosome
instability in somatic hybrids, which can lead to substantial loss of genomic regions of
parental lines [42].
Another way to combine the diversity of two parents for tetraploid rootstock breeding is sexual hybridization between two tetraploid hybrids. This "tetrazyg" breeding strategy was successfully used to select some new tetraploid hybrid rootstocks [49–52].

It is essential to develop knowledge of the mode of inheritance in tetraploid hybrids to optimize the efficiency of the "tetrazyg" strategy. Tetraploids can be distinguished as autotetraploids or allotetraploids based on the composition of their genomes [53,54]. Autotetraploids are derived from genomic material from a single species, while allotetraploids result from combining genomes from different species. The meiotic behavior of tetraploids, particularly the transmission of parental heterozygosity, depends significantly on the phylogenomic structure of the polyploid. Thus, two extreme models are considered: disomy in allotetraploids and tetrasomy in autotetraploids. In allotetraploids, there are two diploid sets of homologous chromosomes. During meiosis, homologous chromosomes are preferentially paired, and only bivalents are formed [53]. This results in a disomic inheritance that, in the extreme case, allows a total restitution of the interspecific parental heterozygosity [55]. In autotetraploids, there are four homologous versions of each chromosome. During meiosis, there is a random pairing with any of its homologs, forming bivalents or quadrivalents. This configuration leads to tetrasomic inheritance, where all allelic combinations are produced [55]. Heterozygosity restitution ranges from 66% for the extreme tetrasomic pattern, in the absence of double reduction, to 55% in the presence of double reduction [56]. In cases where parents are divergent but have retained enough homology to prevent exclusive preferential pairing, chromosome pairing is neither exclusively preferential nor entirely random, and intermediate patterns of inheritance between di- and tetrasomic are expected [5, 18, 53, 55, 57]. Such intermediate inheritance enables the diploid gamete to transmit an important part of the interspecific parental heterozygosity and provides opportunities for interspecific recombination [3]. Intermediate inheritance with disomic tendency has been described in several tetraploid hybrids resulting from somatic hybridization between diploid *C. reticulata* and *P. trifoliata* [18] or chromosome doubling of diploid *Citrus × Poncirus* intergeneric hybrids [3,18]. Intermediate inheritance was also observed for tetraploid hybrids resulting from chromosome doubling of an interspecific *C. micrantha* Wester × *C. medica* L. hybrid [57].

The present study focuses on a tetraploid hybrid population obtained by sexual crossing between a tetraploid Swingle citrumelo and a tetraploid Volkamer lemon. Both tetraploid parents result from chromosome doubling of the corresponding diploid germplasm of intergeneric (*C. paradisi* × *P. trifoliata*) and interspecific origins (*C. reticulata* × *C. medica*). Therefore, both should be considered allotetraploid. According to the agronomic behavior of the diploid rootstock lines described above and their complementarity, it is expected that the hybridization between these two tetraploid parents could help create new rootstocks corresponding to the demand of citrus growers.

The first study of the meiotic behavior of the tetraploid Swingle citrumelo was performed in this population by analyzing the segregation of 159 SNP markers that fully distinguished the genus *Poncirus* from the genus *Citrus* [3]. It revealed an intermediate inheritance with disomic tendency, allowing intergeneric recombination and opening the opportunity for phenotype-genetic association studies.

In the present work, we used a GBS approach anchoring on the *P. trifoliata* reference genome (https://phytozome.jgi.doe.gov/; accessed on 3 March 2022) for simultaneous SNP identification and high-throughput genotyping of the tetraploid population. Diagnostic SNPs for the four ancestral species constitutive of the parents (*C. reticulata*, *C. maxima*, *C. medica*, and *P. trifoliata*) were identified, and the TraceAncestor software [5] was applied to analyze the phylogenomic structure of each hybrid. This analysis permitted the identifications and locations of the interspecific (*C. reticulata*/*C. medica*) and intergeneric (*Citrus*/*Poncirus*) recombination points for the Swingle citrumelo and Volkamer lemon, respectively. Linkage maps and recombination landscapes were established, as well as the preferential pairing and the parental heterozygosity restitution along the genome for the two parents.
Size, density, length, proliferation, expansion, and growth rate of roots are essential characteristics for a good adaptation to drought and an important component for less impact of HLB disease during the dry season. We, therefore, performed phenotyping to evaluate the root development of each hybrid to characterize and select the genotypes with the best root system conformation. These phenotypic data were coupled with GBS genotyping data to search for QTL responsible for root elongation. Finally, we studied the inheritance in each hybrid of some genes and QTL of interest previously located in the *P. trifoliata* genome by [58].

## 2. Results

### 2.1. DSNP Mining

DSNPs were identified from GBS data of eight accessions representative of the basic *Citrus* and *Poncirus* diversity. After filtering with less than 20% of missing data, the VCFHunter pipeline delivered 144,871 dialelic SNPs. With a GST value of 1, we identified a total of 72,452 DSNPs along the 9 chromosomes, including 9444 DSNPs for *C. maxima*, representing 13.03% of the total number of DSNPs. A total of 13,941 DSNPs were found for *C. reticulata* (19.24%), 21,230 for *C. medica* (29.30%), and 27,837 for *P. trifoliata* (38.42%). All the detailed information on DSNPs is provided in Table 1.

<table><caption>**Table 1.** Distribution of the 72,452 diagnostic SNPs (DSNPs) per taxon and per chromosome.</caption><thead><tr><th>Chromosome</th><th>C. maxima</th><th>C. medica</th><th>C. reticulata</th><th>P. trifoliata</th><th>Total</th></tr></thead><tbody><tr><td>C1</td><td>1101</td><td>2220</td><td>1338</td><td>3028</td><td>7687</td></tr><tr><td>C2</td><td>1153</td><td>2654</td><td>1923</td><td>3724</td><td>9454</td></tr><tr><td>C3</td><td>1579</td><td>3938</td><td>2420</td><td>5024</td><td>12,961</td></tr><tr><td>C4</td><td>977</td><td>2221</td><td>1577</td><td>2962</td><td>7737</td></tr><tr><td>C5</td><td>1118</td><td>2193</td><td>1471</td><td>2970</td><td>7752</td></tr><tr><td>C6</td><td>885</td><td>1924</td><td>1325</td><td>2485</td><td>6619</td></tr><tr><td>C7</td><td>880</td><td>1975</td><td>1320</td><td>2626</td><td>6801</td></tr><tr><td>C8</td><td>935</td><td>2018</td><td>1319</td><td>2570</td><td>6842</td></tr><tr><td>C9</td><td>816</td><td>2087</td><td>1248</td><td>2448</td><td>6599</td></tr><tr><td>Total</td><td>9444</td><td>21,230</td><td>13,941</td><td>27,837</td><td>72,452</td></tr><tr><td>%</td><td>13.03</td><td>29.3</td><td>19.24</td><td>38.42</td><td>100</td></tr></tbody></table>

C1–C9: the nine chromosomes of the reference *P. trifoliata* genome (Peng et al. 2020; [58]); %: percentage of DSNPs for the taxon.

### 2.2. Genetic Linkage Maps of Tetraploid Swingle Citrumelo and Volkamer Lemon

GBS data from tetraploid rootstock hybrids were filtered with less than 20% missing data, and VCFHunter software identified 212,072 polymorphisms. For genetic mapping, data were filtered with less than 15% missing data for genomic windows (of 20 DSNPs) and less than 20% for individuals. The method we used to estimate the doses is based on assuming the physical ordering is correct. To generate the genetic map, we used a matrix of 1318 genomics windows and 95 hybrids for the tetraploid Swingle citrumelo and 1016 genomics windows and 96 individuals for Volkamer lemon. Marker-pairwise recombination frequencies were calculated, and nine linkage groups were identified during the linkage grouping constitution. Without preconceptions, we found that all the windows of the same chromosome were localized in the same linkage group. Thus, it confirms the total symmetry between linkage groups and physical assembly (Supplementary Figures S1 and S2). In the Swingle citrumelo genetic map, the number of mapped windows varied from 116 (LG6 and LG9) to 236 (LG3). The genetic map's total size was 435.64 cM, ranging from 28.20 cM (LG8) to 98.46 cM (LG3) (Table 2). In the Volkamer lemon genetic map, the number of mapped windows varied from 94 to 187 for LG6 and LG3, respectively. The total length of the map was 425 cM, with the smallest linkage group of 20.10 cM (LG8) and the largest one of 105.41 cM (LG3) (Table 2).
<table><caption>**Table 2.** Summary of tetraploid Swingle citrumelo and tetraploid Volkamer lemon genetic maps.</caption><thead><tr><th rowspan="2">LG</th><th colspan="2">Swingle Citrumelo</th><th colspan="2">Volkamer Lemon</th></tr><tr><th>Mks/LG</th><th>LG Size cM</th><th>Mks/LG</th><th>LG Size cM</th></tr></thead><tbody><tr><td>1</td><td>149</td><td>36.47</td><td>106</td><td>49.16</td></tr><tr><td>2</td><td>172</td><td>68.94</td><td>126</td><td>27.10</td></tr><tr><td>3</td><td>236</td><td>98.46</td><td>187</td><td>105.41</td></tr><tr><td>4</td><td>143</td><td>37.09</td><td>108</td><td>26.40</td></tr><tr><td>5</td><td>139</td><td>43.34</td><td>105</td><td>79.72</td></tr><tr><td>6</td><td>116</td><td>46.59</td><td>94</td><td>30.25</td></tr><tr><td>7</td><td>128</td><td>34.59</td><td>96</td><td>30.09</td></tr><tr><td>8</td><td>119</td><td>28.20</td><td>97</td><td>20.10</td></tr><tr><td>9</td><td>116</td><td>38.97</td><td>97</td><td>56.78</td></tr></tbody><tfoot><tr><td>Total</td><td>1318</td><td>432.64</td><td>1016</td><td>425</td></tr></tfoot></table>

LG: linkage group; Mks: number of markers; cM: centi morgan.

## 2.3. Recombination Rate

The recombination rates for tetraploid Swingle citrumelo and tetraploid Volkamer lemon were estimated and visualized in Circos (Figure 1).

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Circos plot of the recombination rate along the nine chromosomes. The red lines represent the recombination rate of the Swingle citrumelo and the blue lines Volkamer lemon (Min = 0; Max = 11.72 cM/Mb). The grey histogram represents the gene density in the *P. trifoliata* reference genome (Min = 0; Max = 0.57).</figcaption></figure>

For all chromosomes, we observed a close link between the gene density in the *P. trifoliata* reference genome and the recombination rate of Swingle citrumelo and Volkamer lemon gametes with a decreasing recombination rate around the centromeric regions. Comparison between recombination rates of Swingle citrumelo and Volkamer lemon gametes showed, for chromosome 1, predominantly greater recombination values for Volkamer lemon gametes, except at the chromosome end, where the recombination rate was greater in Swingle citrumelo. Chromosomes 2, 4, 6, and 8 exhibited higher recombination values in Swingle citrumelo. However, the recombination rates over chromosomes 5 and 9 were lower in Swingle citrumelo than Volkamer lemon. For chromosome 3, higher recombination
values were observed for Swingle citrumelo between the beginning of the chromosome and the centromere, while, for the other chromatid of the chromosome, the recombination rate was higher in Volkamer gametes. For chromosome 7, a high recombination rate was recorded at the beginning and a low value at the end of the chromosomes, with a similar profile between both genetic maps. For both parents, we observed a direct association between low combination rates and centromeric and pericentromeric genomic regions that also exhibit low gene sequence density on the *P. trifoliata* genome [58].

The variations of the estimated recombination rates agree with the Marey maps presented in Figure 2.

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

**Figure 2.** Marey map plot of the nine linkage groups of the Swingle citrumelo genetic map (red dot), Volkamer lemon genetic map (blue dot), and trifoliate orange reference genetic map generated in Huang et al. (2018; [59]) (grey dot). The x-axis represents the physical positions (in Mb) on the *P. trifoliata* reference genome (Peng et al. 2020; [58]), and the y-axis represents the position on each genetic map (in cM).

The genetic positions anchored on the physical position of the *P. trifoliata* reference genome [58] were represented for the Swingle citrumelo, Volkamer lemon, and trifoliate orange reference [59] genetic maps. A high recombination rate results in larger genetic distances by physical unit. We observe that the *P. trifoliata* reference genetic map is much larger for all chromosomes. Plateaus corresponding to very low recombination rates were observed in the three maps "for analogous regions of the *P. trifoliata* genome."

## 2.4. *Parental Heterozygosity Restitution and Meiotic Inheritance Analysis*

The Parental Heterozygosity Restitution (PHR) and Preferential chromosome Pairing (PP) were analyzed for both Swingle citrumelo and Volkamer lemon parents along the genome. We used the same matrix we used for genetic mapping to estimate PHR and PP. PHR was determined for each genomic window and chromosome, and PP was calculated from three markers located around the centromere, according to [56,60].

For the diploid Swingle citrumelo gamete, the average PHR over the 1318 genomic windows and the 95 gametes was 85.9%. PHR varied between chromosomes from 0.790 ± 0.007 for chromosome 3 to 0.933 ± 0.006 for chromosome 8. Significant differences in PHR between chromosomes were revealed by a Kruskal-Wallis test (*p*-value < 2.2 × 10⁻¹⁶) followed by a Wilcoxon test (*p*-value < 0.05) (Figure 3a).
<figure><img src="image_7.png" /><figcaption>**Figure 3.** Circos plot of meiotic parameter in (a) Swingle citrumelo and (b) Volkamer lemon gametes along the nine chromosomes. In the first section: violet dots in the scatter plot represent the Parental Heterozygosity Restitution; grey histograms represent the gene sequence density in the *P. trifoliata* reference genome (Peng et al. 2020; [58]), and the blue lines represent the marker density. In the inter section: green dots represent the q-value for the deviation to the tetrasomic gametic segregation model and the red lines represent the limit of significant value.</figcaption></figure>

The average PP was 76.3%. Significant PP variations between chromosomes were confirmed by a one-way ANOVA (p-value < 4.12 × 10⁻¹⁵) followed by a Newman-Keuls test (p-value < 0.05). Chromosomes 6 and 8 had the highest PP values, with 0.925 ± 0.086 and 0.923 ± 0.050, respectively. PP rates are significantly similar for chromosomes 4 and 5, with 0.850 ± 0.043 and 0.848 ± 0.047, and chromosomes 9 and 1, with 0.760 ± 0.045 and 0.742 ± 0.007, respectively. Chromosome 3 had a PP value of 0.690 ± 0.054, followed by chromosome 2 with 0.625 ± 0.033, and chromosome 7 had the lowest PP rate of 0.502 ± 0.029 (Table 3).

<table><caption>**Table 3.** Parameters of the meiotic behaviors of the tetraploid Swingle citrumelo and tetraploid Volkamer lemon.</caption><thead><tr><th rowspan="2">Chromosome</th><th colspan="2">Swingle Citrumelo</th><th colspan="2">Volkamer Lemon</th></tr><tr><th>PHR (Rank)</th><th>PP (Rank)</th><th>PHR (Rank)</th><th>PP (Rank)</th></tr></thead><tbody><tr><td>Chr1</td><td>0.889 ± 0.004 (a)</td><td>0.742 ± 0.007 (d)</td><td>0.826 ± 0.008 (a)</td><td>0.583 ± 0.019 (d)</td></tr><tr><td>Chr2</td><td>0.802 ± 0.006 (b)</td><td>0.625 ± 0.033 (b)</td><td>0.911 ± 0.008 (b)</td><td>0.81 ± 0 (f)</td></tr><tr><td>Chr3</td><td>0.790 ± 0.007 (c)</td><td>0.69 ± 0.054 (c)</td><td>0.719 ± 0.008 (c)</td><td>0.317 ± 0.061 (c)</td></tr><tr><td>Chr4</td><td>0.919 ± 0.004 (d)</td><td>0.85 ± 0.043 (e)</td><td>0.872 ± 0.005 (d)</td><td>0.548 ± 0.061(d)</td></tr><tr><td>Chr5</td><td>0.849 ± 0.011 (e)</td><td>0.848 ± 0.040 (e)</td><td>0.692 ± 0.008 (e)</td><td>0 ± 0 (a)</td></tr><tr><td>Chr6</td><td>0.890 ± 0.012 (adf)</td><td>0.925 (f)</td><td>0.869 ± 0.006 (d)</td><td>0.675 ± 0.012 (e)</td></tr><tr><td>Chr7</td><td>0.817 ± 0.003 (g)</td><td>0.502 ± 0.029 (a)</td><td>0.874 ± 0.007 (d)</td><td>0.663 ± 0.036 (e)</td></tr><tr><td>Chr8</td><td>0.933 ± 0.006 (h)</td><td>0.923 ± 0.050 (f)</td><td>0.897 ± 0.006 (f)</td><td>0.672 ± 0.057 (e)</td></tr><tr><td>Chr9</td><td>0.917 ± 0.003 (f)</td><td>0.76 ± 0.045 (d)</td><td>0.764 ± 0.008 (g)</td><td>0.25 ± 0.054 (b)</td></tr><tr><td>Whole genome</td><td>0.859 ± 0.004</td><td>0.763 ± 0.054</td><td>0.817 ± 0.005</td><td>0.502 ± 0.099</td></tr></tbody></table>

PHR: parental heterozygosity restitution; PP: preferential pairing. PHR rank based on Kruskal-Wallis and Wilcoxon test and PP rank based on ANOVA followed by a Newman-Keuls test. Values with different letters in brackets present significant differences.

The results for the deviation to the tetrasomic gametic segregation model at the individual window level based on q-value [61] are in accordance with previous analysis
(Figure 3a). Among the 1318 windows, 1109 (84.14%) showed significant segregation distortion. Low frequency of significant distortions to the gametic tetrasomic segregation were observed in chromosome 3 (45.76%). The other chromosomes presented a high value of significant gametic segregation distortion, ranging between 72.09% (chromosome 2) and 100% (chromosomes 1, 4, 8, and 9). The lowest average values were observed for chromosomes 3, 2, and 7, with 1.634 ± 0.171, 1.874 ± 0.129, and 2.084 ± 0.102. The highest average value was found in chromosomes 4, 8, and 9 (−log(q-value) > 5) (Table 4a). On the whole, Swingle citrumelo showed an intermediate inheritance with a disomic tendency for all chromosomes and particularly high for chromosomes 6 and 8.

<table><caption>**Table 4.** Deviation to the tetrasomic gametic segregation model in (a) Swingle citrumelo and (b) Volkamer lemon gametes.</caption><thead><tr><th>(a) Chromosome</th><th>N</th><th>-log(q-Value)</th><th>N &lt; -log(0.05)</th><th>%</th></tr></thead><tbody><tr><td>Chr1</td><td>149</td><td>4.276 ± 0.149 (c)</td><td>149</td><td>100%</td></tr><tr><td>Chr2</td><td>172</td><td>1.874 ± 0.129 (ef)</td><td>124</td><td>72.09%</td></tr><tr><td>Chr3</td><td>236</td><td>1.634 ± 0.171 (f)</td><td>108</td><td>45.76%</td></tr><tr><td>Chr4</td><td>143</td><td>5.300 ± 0.115 (b)</td><td>143</td><td>100%</td></tr><tr><td>Chr5</td><td>139</td><td>3.345 ± 0.275 (d)</td><td>122</td><td>87.77%</td></tr><tr><td>Chr6</td><td>116</td><td>4.368 ± 0.378 (c)</td><td>108</td><td>93.10</td></tr><tr><td>Chr7</td><td>128</td><td>2.084 ± 0.102 (e)</td><td>120</td><td>93.75%</td></tr><tr><td>Chr8</td><td>119</td><td>5.628 ± 0.201 (a)</td><td>119</td><td>100%</td></tr><tr><td>Chr9</td><td>116</td><td>5.097 ± 0.111 (b)</td><td>116</td><td>100%</td></tr><tr><td>Total</td><td>1318</td><td>3.492 ± 0.104</td><td>1109</td><td>84.14%</td></tr></tbody></table><table><thead><tr><th>(b) Chromosome</th><th>N</th><th>-log(q-value)</th><th>N &gt; -log(0.05)</th><th>%</th></tr></thead><tbody><tr><td>Chr1</td><td>106</td><td>2.565 ± 0.178 (d)</td><td>94</td><td>88.68%</td></tr><tr><td>Chr2</td><td>126</td><td>4.700 ± 0.222 (a)</td><td>126</td><td>100%</td></tr><tr><td>Chr3</td><td>187</td><td>0.549 ± 0.087 (f)</td><td>18</td><td>9.63%</td></tr><tr><td>Chr4</td><td>108</td><td>3.738 ± 0.113 (c)</td><td>108</td><td>100%</td></tr><tr><td>Chr5</td><td>105</td><td>0.407 ± 0.098 (f)</td><td>7</td><td>6.67%</td></tr><tr><td>Chr6</td><td>94</td><td>3.594 ± 0.153 (c)</td><td>94</td><td>100%</td></tr><tr><td>Chr7</td><td>96</td><td>3.578 ± 0.210 (c)</td><td>96</td><td>100%</td></tr><tr><td>Chr8</td><td>97</td><td>4.346 ± 0.190 (b)</td><td>97</td><td>100%</td></tr><tr><td>Chr9</td><td>97</td><td>0.959 ± 0.157 (e)</td><td>24</td><td>24.74%</td></tr><tr><td>Total</td><td>1016</td><td>2.568 ± 0.113</td><td>664</td><td>65.35%</td></tr></tbody></table>

N: number of genomic windows; -log(q-value): segregation distortion to tetrasomic model estimation with 95% confidence intervals; N > -log(0.05): number of genomic windows with a segregation distortion to tetrasomic model above the significance threshold; %: percentage of significant value.

For the Volkamer lemon diploid gamete, the average PHR over the 1016 genomic windows and the 96 gametes was 81.7%. PHR varied between markers from 0.692 ± 0.008 to 0.911 ± 0.008, and significant differences in PHR between chromosomes were revealed by a Kruskal-Wallis test (p-value < 2.2 × 10⁻¹⁶) followed by a Wilcoxon test (p-value < 0.05). The average PP for Volkamer lemon was 50.2%. Significant PP variations between chromosomes were confirmed by a one-way ANOVA (p-value < 2.2 × 10⁻¹⁶) followed by a Newman-Keuls test (p-value < 0.05). Chromosome 2 had the highest PP value of 0.810. Chromosomes 6, 8, and 7 had PP values significantly similar, from 0.675 ± 0.012 to 0.663 ± 0.036. Chromosomes 1 and 4 showed similar PP rates of 0.583 ± 0.019 and 0.548 ± 0.061. Chromosome 3 had a PP value of 0.317 ± 0.061, followed by chromosome 9 with 0.250 ± 0.054, and chromosome 5 with the lowest PP rate of 0 (Table 3). From the
1016 windows, 664 (65.35%) showed significant deviation of gametic tetrasomic segregation. A low percentage of significant tetrasomic inheritance deviation was observed for chromosomes 5 (6.67%), 3 (9.63%), and 9 (24.74%). Chromosome 1 displayed 88.68% of windows with such significant deviation. All the windows presented significant deviation for all the other chromosomes (chromosomes 2, 4, 6, 7, and 8). The lowest average q-values were observed for chromosomes 5, 3, and 9, with 0.407 ± 0.171, 0.549 ± 0.129, and 0.959 ± 0.157. The highest average values were found in chromosomes 2 and 8, with a q-value over 4 (Table 4b). Volkamer lemon presented intermediate segregation throughout, with critical inter-chromosomal differences (Figure 3b). Chromosomes 1, 2, 4, 6, 7, and 8 showed a disomic tendency, while chromosomes 3 and 9 displayed a tetrasomic tendency. Chromosome 5 exhibited a fully tetrasomic inheritance behavior.

A two-way ANOVA (p-value < 2.2 × 10⁻¹⁶) followed by a Newman-Keuls test (p-value < 0.05) were conducted to compare PP values between Swingle citrumelo and Volkamer lemon gamete populations. Chromosomes 1, 3, 4, 5, 6, 8, and 9 had significantly higher PP in Swingle citrumelo, while chromosomes 2 and 7 had higher PP in Volkamer lemon (Supplementary Figure S3).

## 2.5. *Impact of Preferential Pairing on Tetraploid Hybrid Diversity Structure*

A factorial analysis on the dissimilarity matrix was performed to investigate the genetic distance of the hybrid population from parentage and other genotypes representative of *Citrus* and *Poncirus* diversity. Our goal was to analyze the meiotic behavior and its implication on hybrid population diversity. As we have previously identified interesting preferential pairing differences between chromosomes 2 and 5, we focused our study on these two.

For chromosome 2, factorial analysis (Figure 4a) was represented on three axes (first axis: 25.38%, second axis: 19.45%, and third axis: 17.14%).

<figure><img src="image_8.png" /><figcaption>**Figure 4.** Factorial analysis based on a dissimilarity matrix established from ancestral doses along the genome for (a) chromosome 2 and (b) chromosome 5. Color dots represent parentage and grey dots represent the Volkamer lemon × Swingle citrumelo tetraploid hybrids.</figcaption></figure>

The first axis opposes the *Citrus* species and the *Poncirus* species. Tetraploid Swingle citrumelo and Volkamer lemon were differentiated by the second axis, with the intergeneric hybrids in intermediate positions. Hybrids were extended with a large diversity between *P. trifoliata* and ‘Marsh’ pummelo, which corresponds to the recombination of Swingle citrumelo. For this chromosome, Volkamer lemon displayed a high value of preferential pairing, which limits recombination events.
For chromosome 5, factorial analysis (Figure 4b) showed a first axis (27.72%) as opposed to *C. reticulata* and *C. maxima* against *C. medica*. The second axis primarily differentiated *Poncirus* and *Citrus* genotypes (20.91%). Hybrid orientation was done on the first axis with important diversity between *C. reticulata* and *C. medica*, which correspond to the recombination of Volkamer lemon. The null preferential paring estimated in Volkamer lemon leads to increased *C. medica/C. reticulata* recombination rates and Volkamer gamete diversity.

## 2.6. Inheritance in the Tetraploid Hybrids of Candidate Genes for Pest and Disease Resistance

The inheritance in the hybrid population was analyzed for the candidate genes associated with HLB tolerance, CTV resistance, and nematode resistance identified by Peng et al. (2020) [58]. The localization of these genes on the *Poncirus* genome was determined, and we analyzed the phylogenomic structure of hybrids at these positions. As the tolerance traits come from the *Poncirus* ancestor, hybrids carrying the *Poncirus* haplotype have potentially inherited the favorable trait carried by the candidate gene.

According to Huang et al. (2018) [59], HLB tolerance is associated with 11 genes localized in scaffolds 6 and 9. In scaffold 6, 4.21% of the hybrids did not inherit the *Poncirus* genome on the positions of the candidate gene, 83.16% had one dose of *Poncirus*, and 9.47% of the hybrids inherited two doses of *Poncirus*. *Poncirus* dose remained undetermined for 3.16% of the hybrids (Table 5).

<table><caption>**Table 5.** Inheritance in the hybrids population of candidate genes (according to Peng et al., 2020, [58]) located in the *P. trifoliata* reference genome.</caption><thead><tr><th rowspan="2">Gene ID</th><th rowspan="2">Trait</th><th colspan="3">Location on <i>P. trifoliata</i> Reference Genome</th><th colspan="4">Hybrids % with n Dose of <i>Poncirus</i></th></tr><tr><th>Scaffold</th><th>Start</th><th>End</th><th>0</th><th>1</th><th>2</th><th>N</th></tr></thead><tbody><tr><td>Ptrif.0006s1501.1</td><td rowspan="10">Huanglongbing tolerance</td><td>Scaffold 6</td><td>17843557</td><td>17847864</td><td rowspan="5">4.21</td><td rowspan="5">83.16</td><td rowspan="5">9.47</td><td rowspan="5">3.16</td></tr><tr><td>Ptrif.0009s1449.1</td><td>Scaffold 9</td><td>16743738</td><td>16750183</td></tr><tr><td>Ptrif.0009s1451.1</td><td>Scaffold 9</td><td>16761424</td><td>16771587</td></tr><tr><td>Ptrif.0009s1453.1</td><td>Scaffold 9</td><td>16803290</td><td>16817272</td></tr><tr><td>Ptrif.0009s1458.2</td><td>Scaffold 9</td><td>16966396</td><td>16974146</td></tr><tr><td>Ptrif.0009s2550.1</td><td>Scaffold 9</td><td>16873308</td><td>16878288</td></tr><tr><td>Ptrif.0009s2650.1</td><td>Scaffold 9</td><td>16878879</td><td>16881473</td></tr><tr><td>Ptrif.0009s1595.1</td><td>Scaffold 9</td><td>19302619</td><td>19305258</td><td rowspan="4">4.21</td><td rowspan="4">88.42</td><td rowspan="4">2.11</td><td rowspan="4">5.26</td></tr><tr><td>Ptrif.0009s1596.1</td><td>Scaffold 9</td><td>19351433</td><td>19355794</td></tr><tr><td>Ptrif.0009s1599.1</td><td>Scaffold 9</td><td>19383517</td><td>19386756</td></tr><tr><td>Ptrif.0009s1600.1</td><td>Scaffold 9</td><td>19400932</td><td>19403616</td></tr><tr><td>Ptrif.0007s1586.1</td><td rowspan="4">Citrus tristeza virus disease resistance</td><td>Scaffold 7</td><td>11780276</td><td>11787445</td><td rowspan="4">10.53</td><td rowspan="4">77.89</td><td rowspan="4">9.47</td><td rowspan="4">2.11</td></tr><tr><td>Ptrif.0007s1587.1</td><td>Scaffold 7</td><td>11823916</td><td>11826636</td></tr><tr><td>Ptrif.0007s1590.1</td><td>Scaffold 7</td><td>11860949</td><td>11863630</td></tr><tr><td>Ptrif.0007s1595.1</td><td>Scaffold 7</td><td>11908817</td><td>11912236</td></tr><tr><td>Ptrif.0007s1378.1</td><td rowspan="12">Nematode resistance</td><td>Scaffold 7</td><td>9840644</td><td>9843519</td><td rowspan="5">11.58</td><td rowspan="5">78.95</td><td rowspan="5">7.37</td><td rowspan="5">2.11</td></tr><tr><td>Ptrif.0007s1394.1</td><td>Scaffold 7</td><td>9974384</td><td>9977152</td></tr><tr><td>Ptrif.0007s1395.1</td><td>Scaffold 7</td><td>9985962</td><td>9988577</td></tr><tr><td>Ptrif.0007s1396.1</td><td>Scaffold 7</td><td>9990366</td><td>9991216</td></tr><tr><td>Ptrif.0007s1398.1</td><td>Scaffold 7</td><td>10007442</td><td>10010437</td></tr><tr><td>Ptrif.0007s1402.1</td><td>Scaffold 7</td><td>10043069</td><td>10045927</td><td rowspan="5">11.58</td><td rowspan="5">76.84</td><td rowspan="5">9.47</td><td rowspan="5">2.11</td></tr><tr><td>Ptrif.0007s1404.1</td><td>Scaffold 7</td><td>10056890</td><td>10060283</td></tr><tr><td>Ptrif.0007s1406.2</td><td>Scaffold 7</td><td>10065703</td><td>10067962</td></tr><tr><td>Ptrif.0007s1411.2</td><td>Scaffold 7</td><td>10097497</td><td>10100154</td></tr><tr><td>Ptrif.0007s1415.1</td><td>Scaffold 7</td><td>10123503</td><td>10126433</td></tr><tr><td>Ptrif.0007s2703.1</td><td>Scaffold 7</td><td>10139209</td><td>10144561</td><td rowspan="4">11.58</td><td rowspan="4">78.95</td><td rowspan="4">9.47</td><td rowspan="4">0</td></tr><tr><td>Ptrif.0007s1481.1</td><td>Scaffold 7</td><td>10687572</td><td>10691702</td></tr><tr><td>Ptrif.0007s1484.1</td><td>Scaffold 7</td><td>10729463</td><td>10732201</td></tr><tr><td>Ptrif.0007s1501.1</td><td>Scaffold 7</td><td>10812404</td><td>10852651</td></tr><tr><td>Ptrif.0007s1502.1</td><td>Scaffold 7</td><td>10860634</td><td>10863348</td><td></td><td></td><td></td><td></td></tr></tbody></table>

Hybrids % with n dose of *P. trifoliata*: percentage of hybrids containing 0, 1, or 2 doses of *P. trifoliata* at the corresponding genes position. N: percentage of undetermined origin.

In scaffold 9, 10 candidate genes are localized in two regions where 5.26% and 4.21% of the hybrids did not inherit the *Poncirus* genome on the desired positions, 90.53% and 88.42% had one dose of *Poncirus*, and 1.05 and 2.11% of the hybrids inherited two doses of *Poncirus*. *Poncirus* doses remained undetermined for 3.16% and 5.26% of the hybrids.
Four CTV resistance candidate genes are localized in scaffold 7 (11780276–11912236). For the positions of interest, 10.53% of the hybrids did not inherit *Poncirus*, 77.89% inherited one dose, and 9.47% two doses of *Poncirus*. In addition, 2.11% remained undetermined.

Fifteen candidate genes for nematode resistance are localized in a 1 Mb region of scaffold 7 of the *Poncirus* genome, and 11.58% of the hybrids did not inherit the *Poncirus* genome. According to the genes, between 76.84% and 78.95% of the hybrids inherited one dose and 7.37% to 9.47% inherited two doses of *Poncirus*. Indetermination concerned between 0 and 2.11% of the hybrids.

A total of 77.89% of hybrids have inherited at least one dose of *Poncirus* for all genomic regions considered.

## 2.7. Root Phenotyping and Genomic Association

The root length, root area, number of secondary roots, and root diameter were evaluated in the population. Unfortunately, we were not able to propagate cuttings of Swingle citrumelo 4×. The measurements were highly heterogeneous between hybrids, and we observed a large diversity in root system morphology (Supplementary Figure S4). Statistical analyses using linear-mixed effects models showed a significant effect of genotype (*p*-value < 0.001) on the different measured traits. The hybrid values were distributed mainly outside the Swingle citrumelo 2× and Volkamer lemon 2×/4× range (Figure 5).

<figure><img src="image_9.png" /><figcaption>**Figure 5.** Distribution of the four root phenotypic traits in the 4× Volkamer lemon × 4× Swingle citrumelo progenies. (a) Number of secondary roots; (b) root surface; (c) root length; (d) root diameter.</figcaption></figure>

The tetraploid hybrids displayed mostly lower values than the diploid and tetraploid parents regarding the number of secondary roots, root surface, and root length, while showing greater root diameters. A PCA analysis confirmed the relationships between the different measured variables (Supplementary Figure S5). The first axis, explaining 46.67% of the variables, is characterized by root area, number of secondary roots, and root length. The arrows related to these variables are clustered around the first axis, so these traits are highly correlated. Root diameter drives the second principal axis, explaining 21.24% of the root variables. However, on axis 1, this variable also seems opposed to root length.

With these phenotyping data, we tried to identify some genetic determinants of the root system architecture in the hybrids through genetic association analyses. We used ancestral doses of *C. medica* and *P. trifoliata* in the diploid gamete populations of Volkamer lemon and Swingle citrumelo that generated our tetraploid hybrids. A significant genome
association signal (*p*-value > 0.01) was observed for the two gamete populations at the end of chromosome 1 (Figure 6) for the number of secondary roots and root surfaces.

Another significant signal was noticed for the Volkamer lemon gamete at the beginning of chromosome 9 for root diameter. In addition, a weaker significant signal was detected for the Swingle citrumelo gamete at the end of chromosome 5 for the root surface trait.

Considering the marker with the lower *p*-value, the boxplot analyses of the hybrids (Supplementary Figure S6) reveal a positive effect of the *P. trifoliata* haplotype and a negative effect of *C. medica* for the three traits in all identified regions. However, the distribution of the different classes of ancestral doses (zero, one, or two) in the diploid gametes is highly unbalanced, with a strong excess of one dose (heterozygous gametes) resulting from the intermediate segregation of the two allotetraploid parents.

<figure><img src="image_10.png" /><figcaption>**Figure 6.** Evaluation of the association between ancestors doses in diploid gametes and root traits. Manhattan plots from GLM procedure: X axis: position on the nine *P. trifoliata* chromosome; Y axis: -log₁₀(*p*-value); (a,e): root length, (b,f): root diameter, (c,g): number of secondary roots, (d,g): root surface. (a–d): *C. medica* doses in Volkamer lemon gametes. (e–h): *P. trifoliata* allelic doses in Swingle citrumelo gametes.</figcaption></figure>

**3. Discussion**

3.1. *GBS Coupled with TraceAncestor Was Efficient for Estimating Haplotype Ancestral Doses All along the Genome of the Tetraploid Hybrids*

The genetics of polyploid species are much more complex than diploid ones. In partic-
ular, the accurate study of polyploid populations requires the estimation of the allelic copy
number for heterozygous genotypes. Allele dosage can affect gene expression and, there-
fore, the phenotype. Their determination is, therefore, essential for genotype/phenotype
association studies [62]. Various methods for estimating allelic dosage in citrus polyploid
plants have been developed using co-dominant markers. For Simple Sequence Repeat
(SSR) markers, the MAC-PR method [63], proposed to deal with differential amplification
intensities between alleles, has been successfully used in citrus [34,64]. The KASPar geno-
typing technique, based on competitive allele-to-allele PCR, is also suitable for studying
SNP marker allelic patterns within polyploid populations [65]. The KASPar genotyping
technique has been successfully applied for the analysis of the origin of 2n gametes and
their genetic structures [33,35,60,66] as well as for the analysis of the meiotic behavior of
tetraploid progenitors [3,18,56,57,67,68].

However, the two previous methods are adapted for studies requiring a limited num-
ber of markers, but not for high throughput analysis. Taking advantage of Next Generation
Sequencing (NGS), the authors of [5] developed a likelihood method based on the multilo-
cus number of ancestral reads of successive phased markers to estimate the contribution of ancestral haplotypes in diploid and polyploid admixed plants. Our tetraploid population is the result of crossing tetraploid Swingle citrumelo (C. paradisi × P. trifoliata) and Volkamer lemon (C. medica × C. reticulata) rootstocks [3]. According to recent phylogenomic studies, the grapefruit parent (C. paradisi Macf.) of Swingle citrumelo results from hybridization between a pummelo (C. maxima) and a sweet orange (C. sinensis), sweet orange being itself an admixture between C. reticulata and C. maxima [5,69,70]. Therefore, our Volkamer lemon × citrumelo hybrids have a complex genetic structure with inheritance from four ancestral parents, including C. medica, C. reticulata, C. maxima, and P. trifoliata. From our GBS study, we identified 72,452 DSNPs of these 4 ancestral species along the 9 chromosomes: 9444 DSNPs for the C. maxima ancestor, 13,941 DSNPs for C. reticulata, 21,230 for C. medica, and 27,837 for P. trifoliata. Using the TraceAncestor tool [5], the estimation of the allelic doses of the C. medica and P. trifoliata ancestors in each tetraploid hybrid allowed us to deduce the phylogenomic structure of the corresponding Volkamer lemon ovule and Swingle citrumelo pollen that generate each hybrid and, therefore, to identify and locate the interspecific (C. reticulata/C. medica) and intergeneric (Citrus/Poncirus) recombination events, respectively, for the Volkamer lemon and Swingle citrumelo parent.

## 3.2. *The Tetraploid Volkamer Lemon and Swingle Citrumelo Parents Display Intermediate Inheritance with a Disomic Tendency and Variability between Parents and Chromosomes for Preferential Chromosome Pairing*

According to [56], we analyzed preferential pairing in the Swingle citrumelo and Volkamer lemon using three centromeric and pericentromeric genetic windows with the assumption of no double reduction. In addition, Parental Heterozygosity Restitution (PHR) of the diploid gamete was computed for each marker. For most of the genome, intermediate segregation with a clear disomic tendency was identified in Swingle citrumelo (PHR: 0.859 ± 0.004; PP: 0.763 ± 0.054). Volkamer lemon displayed an intermediate inheritance with slight disomic behavior (PHR: 0.817 ± 0.005; PP: 0.502 ± 0.099). Such intermediate inheritance with a disomic tendency has already been described in the 'Giant key' doubled-diploid Mexican lime [57,67] and in Flhorag1 allotetraploid rootstock [18]. Inter-chromosomal comparisons revealed significant differences in PP and PHR, resulting in segregation pattern differences between chromosomes.

For tetraploid Swingle citrumelo, LG7 displayed a clear intermediate inheritance (PP = 0.502). Furthermore, other LGs (LG1, LG2, LG3, LG4, LG5, LG6, LG8, and LG9) displayed an intermediate inheritance with a disomic inheritance tendency (PP comprised between 0.625 and 0.925). Disomy was particularly marked for LG6 and LG8, with more than 92% preferential chromosome pairing.

For tetraploid Volkamer lemon, a complete tetrasomic inheritance was identified for LG5 with an absence of PP. LG3 and LG9 showed an intermediate inheritance with a tendency towards a tetrasomic inheritance (PP [LG9] = 25%; PP [LG3] = 32%). LG1 and LG4 presented a clear intermediate inheritance with PPs of 58% and 55%, respectively. LG2, LG6, LG7, and LG8 exhibited an intermediate inheritance with a tendency towards a disomic inheritance with PPs between 66% and 81%.

Such inter-chromosomal differences in segregation pattern have already been observed in a doubled-diploid clementine displaying tetrasomic segregation, intermediate segregation with a tetrasomic tendency, and intermediate segregation with a disomic tendency [56] and in a citrus allotetraploid somatic hybrid ('Nova' tangelo + 'HB' pummelo) that displayed a mixed disomic, tetrasomic, and intermediate inheritance [71]. In addition, this kind of segregation was observed in other species, such as in colchicine-induced allotetraploid *Musa* [72].

## 3.3. *Apparent Interspecific and Intergeneric Recombination Rates Are Limited by Preferential Pairing and Are Very Low in Centromeric-Pericentromeric Areas*

Despite predominant disomic inheritance, we observed interspecific recombination in the nine citrus chromosomes for the two allotetraploid parents, which allowed us to
establish the two parental genetic maps. The construction of high-density genetic maps is valuable for breeders because it will potentially facilitate genomic studies and identify genomic regions with interesting agronomic traits in citrus. Our genetic maps have been anchored on the *P. trifoliata* reference genome [58]. During the elaboration of the genetic maps, we observed that all the windows of the same chromosome were in the same linkage group. This observation confirms the total synteny of the genetic linkage maps with the *P. trifoliata* reference genome. Due to the quality of the *P. trifoliata* assembly and the collinearity that has already been demonstrated between the genus *Citrus* and *Poncirus* [3,59,73–75], we assumed that the physical ordering of *P. trifoliata* was the correct ordering. For Swingle citrumelo, the genetic map spanned 432.64 cM and consisted of 1318 windows of 20 DSNPs. For Volkamer lemon, the map consisted of 1016 markers spanning 425 cM. Both tetraploid Volkamer lemon and Swingle citrumelo genetics maps are more than two times smaller than the diploid trifoliate orange genetic map and most of the saturated genetic maps published for *Citrus* [59,74,76–78]. Our tetraploid maps are also half the size of the genetic map of the tetraploid clementine, which displayed a mostly tetrasomic inheritance with very low levels of preferential pairing [56].

Ahmed et al. (2020) stated that the high level of chromosome pairing in the allotetraploid giant key lime was the main factor resulting in the low apparent interspecific recombination between *C. medica* and *C. micrantha*, the two progenitors. Indeed, high preferential pairing strongly limits the proportion of gametes that can undergo interspecific chiasmata and, consequently, interspecific recombination. This conclusion is supported by our observation of chromosomes with marked differences for preferential pairing between Volkamer lemon and Swingle citrumelo. By comparing the genetic map and chromosomic preferential pairing between Volkamer lemon and Swingle citrumelo, a clear negative relation is found between genetic distance and preferential pairing. In chromosome 2, Volkamer lemon displayed a higher value of preferential pairing; thus, we observed a lower recombination rate and smaller genetic distance. In contrast, in chromosome 5, a null value of PP was identified in Volkamer lemon coupled with a higher recombination rate and larger genetic distance. For this chromosome, the apparent intergeneric recombination between *Citrus* and *Poncirus* was severely limited in citrumelo by preferential chromosome pairing.

Recombination rates in interspecific hybrids [79] may also be affected by the divergence between parental genomes, whatever the ploidy level. Decrease of recombination linked with genome sequence divergence has been described at the diploid level for several species [80–84] and suspected in *Citrus* from a comparative genetic mapping study between clementine and sweet orange [78]. The authors suggested that higher recombination rates in some genomic regions of sweet orange were associated with the interspecific heterozygosity of *C. reticulata/C. maxima*, which is higher in sweet orange than in clementine.

We cannot discard that, in addition to the predominant preferential pairing, the divergence between *Citrus* and *Poncirus* or between *C. reticulata* and *C. medica* contributed to limiting the recombination rates for the tetraploid Swingle citrumelo and Volkamer lemon, respectively, even for gametes resulting from the homologous pairing. Comparative genetic mapping between diploid and tetraploid Swingle citrumelo and Volkamer lemon should allow this hypothesis to be tested. Some works revealed increased recombination in polyploid compared to diploid. Recombination between homologous chromosomes was more frequent in autotetraploid and allotetraploid *Arabidopsis* than in diploids [85]. Similarly, *Brassica* allotriploid hybrids (AAC) showed increased recombination between homologous A chromosomes compared to AA diploids [86]. In these two works, the increased recombination concerned homologous chromosomes. In contrast, in our work with interspecific doubled diploid, we can only see the interspecific recombination between heterologous chromosomes, the homologous chromosome being fully identical. Our conclusion of reduced recombination in allotetraploid concerns, therefore, only the interspecific recombination. In *Brassica*, while the diploid AA displays a classical marked decrease of recombination in pericentromeric areas, a homogeneous distribution of recombination all
along the genome between the homologous A chromosome was observed in the AAC triploid [86]. In the case of our two interspecific doubled diploids, the recombination pattern displayed similar strong limitations in centromeric and pericentromeric regions, as already described in diploid citrus mapping studies [56,67,78].

### 3.4. Preferential Chromosome Pairing in Tetraploid Rootstock Resulting from Chromosome Doubling of Interspecific Hybrids Is Unfavorable for QTL Analysis

Genome-Wide Association Studies (GWAS) were performed by combining phenotyping and genotyping data at the gamete population level to identify and characterize genes involved in root growth and development. Preliminary results of GWAS analyses showed mostly weak signals. However, a dose effect of *Poncirus* and *C. medica* ancestors was identified at the end of chromosome 1 for the number of secondary roots and root surface. In Volkamer lemon gametes, a significant signal was detected for the root diameter trait at the beginning of chromosome 9. Furthermore, a weaker signal was also observed in the Swingle citrumelo gamete population for roots surface at the end of chromosome 5. *C. medica* haplotypes appeared unfavorable, while *P. trifoliata* haplotypes were favorable for all identified QTLs. Strong preferential chromosome pairing during meiosis of interspecific doubled diploid complicates QTLs discovery due to reduced recombination limiting genetic disparities between hybrids and resulting in a very important excess of heterozygous gametes compared to the two classes of homozygotes. More robust results could be obtained in such allotetraploid populations by increasing the number of individuals. Another option, even if it will not allow for analysis of the ploidy and allele doses effects, should be to analyze the architecture of targeted traits using diploid progenies resulting from hybridization between diploid Volkamer lemon and Swingle citrumelo. After validation at the tetraploid level, genetic markers linked to the identified QTL may be used for marker-assisted selection on the tetraploid progenies.

### 3.5. Implication for Rootstock Breeding Programs

The 'tetrazyg' strategy has been used previously for rootstock breeding [49–51] and has been shown to be effective in generating interesting recombinant tetraploid progenies. Tetrazyg breeding can be performed using somatic hybrids or doubled diploids of preselected diploid rootstock arising from chromosome doubling in nucellar cells, such as the two parents used in this work. We took advantage of the very reduced rate of the apomixis of the tetraploid Volkamer lemon compared to other tetraploid rootstocks [3] to obtain numerous tetrazyg hybrids in the 4× Volkamer lemon × 4× Swingle citrumelo hybridization. Intergeneric hybridization between the genus *Poncirus* and *Citrus* is commonly performed for rootstock improvement because it allows for the combination of interesting characteristics of both genera. When using intergeneric hybrids such as Swingle citrumelo, citrange, and citrandarin for rootstock breeding, it is, therefore, essential to transmit the complementary favorable genes of the *Citrus* and *Poncirus* grandparents to the progenies.

Chromosome preferential pairing has an effect on the transmission of parental heterozygosity. It is also strongly involved in reducing effective interspecific recombination and, thus, the genetic and phenotypic diversity of the hybrid population. In this study, we identified the high value of preferential pairing, leading mainly to an intermediate inheritance for both parents. This is an interesting feature in the case of rootstock breeding. Indeed, PP strongly limits the recombination rate, so a great part of the genetic value of the parents is transmitted to the tetraploid sexual progeny. The highly marked disomy observed in 4× Swingle citrumelo will allow retention of many genes of interest from the original Swingle citrumelo that required a long selection effort.

In a recent study [58] of the *Poncirus* reference genome, various loci of interest were involved in CTV resistance [9], nematode resistance [8], and HLB tolerance [59]. Thirty NBS-encoding genes involved in these traits of interest were located in the *P. trifoliata* genome on chromosomes 6, 7, and 9. Thus, these chromosomes would play a very important role in transmitting resistance or tolerance to HLB, CTV, and nematodes [58]. The inheritances of
these genes were studied within the tetraploid rootstock population. A total of 77.89% of the hybrids inherited at least one copy of the *P. trifoliata* genome for all regions of interest. These tetraploid hybrids may have inherited the trait of interest and will be submitted to further wide evaluation for resistance to abiotic and biotic stresses and agronomic behavior.

A vigorous root system is one of the main selection criteria in rootstock improvement programs. Previous studies suggested that, in the field, doubled-diploid rootstocks presented a superficial ramified root system and had less ability to capture soil resources during drought stress. We evaluated the root system architecture in our hybrids population in the present study. Significant phenotypic variability between hybrids was observed for all measured traits (number of secondary roots, root length, root surface area, and root diameter). Even though most of the hybrids presented a limited root development compared with diploid Citrumelo and Volkamer lemon, suggesting a limited capacity to go deep in the soil to absorb water, a limited number of hybrids presented root development higher or similar to Swingle citrumelo 2×, making these genotypes of interest for further evaluation.

# 4. Materials and Methods

## 4.1. Plant Material

Sexual hybridization between the tetraploid Volkamer lemon, SRA 1122 (C. reticulata × C. medica), as the female parent and tetraploid intergeneric Swingle citrumelo, SRA 1112 (C. paradisi × P. trifoliata), as a male parent was performed to obtain 103 tetraploid rootstocks hybrids. The two parents were selected several years ago among their diploid corresponding rootstock seedlings [37]. As described by [3], zygotic seedlings were selected by SNP marker genotyping, and tetraploidy of the 103 hybrids was confirmed by flow cytometry analysis of young leaves.

Other genotypes representative of the *Citrus* and *Poncirus* diversity were added to the study for diversity analysis, including three accessions representative of the trifoliate orange horticultural group ('Flying Dragon', 'Pomeroy', and 'Kryder' trifoliate oranges), one accession representative of the pummelo horticultural group ('Seedless' pummelo), two accessions representative of the citron horticultural group ('Humpang' and 'Poncire commun' citron), and two accessions representative of the mandarin horticultural group ('Cleopatra' and 'Sunki' mandarin) and the 'Marsh' grapefruit as a representative of *C. paradisi* (the second parent of the Swingle citrumelo). The Inrae-Cirad Citrus Biological Resources Center provided all these genotypes in San Giuliano, Corsica, France.

## 4.2. DNA Extraction and Genotyping by Sequencing (GBS)

Following the protocol of [69], genomic DNA was isolated using the Plant DNAeasy® kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions. The concentration of genomic DNA was adjusted to 20 ng/µL, and ApeKI GBS libraries were prepared following the protocol described by Elshire et al. [87] and adapted to citrus by Ahmed et al. [67] with the use of a selective base (A) during the PCR step to reduce complexity as recommended by Sonah et al. [88]. Single-end and pair-end sequencing were performed on a double and single lane of an Illumina HiSeq4000.

## 4.3. SNP Calling

Raw sequencing data were demultiplexed using GBSX v1.1.2 [89] for paired-end reads and using Sabre (https://github.com/najoshi/sabre; accessed on 3 May 2022) for single-end reads. Cutadapt v1.8 [90] was used to remove the adapter sequences and low-quality bases (-q 20,20). Processed reads shorter than 30 were discarded. SNP calling and filtering were performed with the VcfHunter package described in [91] (https://github.com/SouthGreenPlatform/VcfHunter/; accessed on 3 May 2022) using VcfPreFilter.1.0.py and vcfFilter.1.0.py utilities. This procedure included mapping reads against the *P. trifoliata* reference genome (v1.3.1, https://phytozome.jgi.doe.gov/; accessed on 3 May 2022) using the BWA software [92]. For SNP calling, a minimal coverage by accession to keep genotype calling was set to 10, and a Minimal Allele Frequency (MAF) to keep genotype calling
was set to 0.05. We only considered diallelic polymorphic positions. SNPs and genotype information (VCF file) were then imported into the SNIPlay3 application [93] for subsequent SNP exploration and analysis. The Illumina Hiseq sequencing raw data are available in the NCBI SRA (Sequence Read Archive) under the accession number PRJNA898312.

### 4.4. DSNPs Identification

Ancestral dose analysis along the genome was performed based on the number of allelic reads of ancestral diagnostic SNPs (DSNPs). According to the known phylogenomic structures of the parents, four ancestral species were considered: *P. trifoliata*, *C. medica*, *C. maxima*, and *C. reticulata*. Diagnostic SNPs for the ancestral species were selected using GBS data from eight accessions representative of the considered taxa: three for *P. trifoliata*, one for *C. maxima*, two for *C. medica*, and two for *C. reticulata*.

The search for diagnostic SNPs for the three *Citrus* ancestral species and *P. trifoliata* was based on the estimation of the inter-population differentiation parameter (GST) defined by [94]. For each taxon, GST estimation was performed considering two sub-populations: (1) the taxon concerned ($T_i$) and (2) a theoretical population of the three other taxa ($T_{-i}$). Analyses were performed from the estimated allele frequency of each taxon considering the same population size for each taxon to estimate the frequency of the two sub-populations ($T_i$ and $T_{-i}$) and the frequency of the whole population (Tot):

$$GST = \frac{He_{Tot} - \frac{(He_{T_i} + He_{T_{-i}})}{2}}{He_{Tot}} \quad (1)$$

where He is the expected proportion of heterozygous loci per individual under Hardy-Weinberg equilibrium ($He = 1 - \sum p_i^2$, $p_i$ is the frequency of a given allele in the population or sub-population considered). GST values ranged from zero to one. GST = 1 means that the two taxa were fully differentiated and were the criterion for DSNP identification.

Neighbor-joining and factorial analysis were computed using DARwin software version 6.0 [95]. Both were based on the Manhattan index:

$$D_{i-j} = 1/K \sum_{1}^{K} |x_{ik} - x_{jk}| \quad (2)$$

where i and j are the two considered individuals, k is the considered locus, K is the total number of loci, and $x_{ik}$ is the frequency of the alternative allele at locus k for the individual i. Three-dimensional (3D) visualization was done using the R package {rgl} [96].

### 4.5. Matrix Preparation

A previous study [5] showed that polyploid allelic doses could not be estimated accurately at a single SNP locus from the GBS data generated with our protocol. Therefore, we developed an efficient approach based on the frequency of reads of phased successive markers to estimate the haplotypic doses in successive genomic windows of n markers. The published tool, called TraceAncestor, was successfully applied to analyze the phylogenomic structure of *Citrus* polyploid germplasm [5]. Our Volkamer lemon × Swingle citrumelo progeny has an adequate phylogenomic structure to infer allelic doses at the gamete level using the TraceAncestor tool and ancestral DSNP information. Indeed, Volkamer lemon is in *C. reticulata*/*C. medica* heterozygosity all along the genome, while the Swingle citrumelo genome alternates between *P. trifoliata*/*C. reticulata* and *P. trifoliata*/*C. maxima* heterozygosity. Therefore, estimating the allelic doses of *C. medica* and *P. trifoliata* in each tetraploid hybrid allows us to infer without ambiguity the phylogenomic structure of the corresponding Volkamer lemon ovule and Swingle citrumelo pollen (Figure 7).

The GBS data, the DSNP matrix, and TraceAncestor software were used to analyze the contribution of the four ancestral species along the genome of the intergeneric hybrids, following the method described by [5] using genomic windows of 20 DSNPs. We generated
two matrixes from the hybrid mosaic structure: one for the Volkamer lemon parent with the dose estimation of *C. medica* for each window and another for the Swingle citrumelo parent with the dose estimation of *P. trifoliata*. Then, singletons were identified and replaced with missing data, as recommended by Van Os et al. [97], using a homemade Excel page routine. At the same time, a few individuals displaying an aberrant number of recombination, set by examining the global recombination distribution, were removed. These two matrixes were used for genetic mapping and the analysis of the meiotic behavior of each parent.

<figure><img src="image_11.png" /><figcaption>**Figure 7.** Phylogenomic structure of tetraploid parents and inference of phylogenomic structure of tetrazyg hybrids. (a) Origin and phylogenetic karyotypes of the tetraploid parents used in this paper. (b) Examples of Volkamer lemon ovule and Swingle citrumelo pollen phylogenetic structure for chromosome 1 and (c) examples of phylogenetic structure of chromosome 1 for four intergeneric Volkamer lemon × Swingle citrumelo tetraploid hybrids. C.D.: chromosome doubling of nucellar cells. Ancestral genome contribution: *C. reticulata* (red), *C. maxima* (blue), *C. medica* (yellow), *P. trifoliata* (purple), indeterminate area (grey).</figcaption></figure>

## 4.6. *Mapping Analysis*

The genetic map of the tetraploid Swingle citrumelo and Volkamer lemon was built using the {pergola} R package [98]. Genotyping data were filtered to less than 15% of missing data for markers and 20% for hybrids.
The good quality assembly of *Poncirus*, the strong collinearity already demonstrated between the genus *Poncirus* and *Citrus* [3,59,73–75], and our choice to work on genomic windows of 20 successive markers led us to work based on physical ordering. In the present mapping analysis, we did not try to establish the order of the mapped windows, but we wanted to study the recombination rates along the genome.

The map was created using the multidimensional scaling algorithm and the Kosambi mapping function, allowing incomplete interference among the recombination events. The synteny and collinearity of both Swingle citrumelo and Volkamer lemon genetic maps with the reference *Poncirus* genome were visualized using Circos [99]; http://circos.ca (accessed on 25 February 2023), in Galaxy [100]. In addition, Marey maps were drawn using Microsoft Excel 2016 to visualize the collinearity between the genetic and physical position of the markers on the *Poncirus* reference genome [58].

## 4.7. Recombination Rate

Changes in the meiotic recombination rate along the genome are estimated with the online software « MareyMap Online » https://lbbe-shiny.univ-lyon1.fr/MareyMapOnline/ [101] accessed on 9 June 2022, with the following parameters: Loess estimation method with Span = 0.15. Recombination rates were estimated and visualized in Circos for the Swingle citrumelo and Volkamer lemon genetic map and the trifoliate orange genetic map generated by Huang et al. (2018; [59]).

## 4.8. Parental Heterozygosity Restitution (PHR)

The Parental Heterozygosity Restitution (PHR) was calculated at each locus as the percentage of individuals with the heterozygous allelic configuration. The estimation of PHR for each chromosome is the average of the values for all the markers on the chromosome.

## 4.9. Estimation of Preferential Pairing (PP)

Preferential Pairing (PP) defines the proportion of gametes resulting from the exclusive pairing of homologous chromosomes and ranges between 0 (full tetrasomic inheritance characteristic of autotetraploids) and 1 (full disomic inheritance in strict allotetraploids).

For a doubled diploid, under intermediate inheritance, marker segregation depends directly on the rates of chromosome Preferential Pairing (PP) and double reduction for the marker considered [56]. In the absence of recombination in the centromeric area and, hence, a null value for double reduction, there is a direct relationship between PP and PHR for the centromeric markers [56]. The PP values for three centromeric and peri-centromeric markers were estimated for each chromosome using the maximum likelihood approach proposed for centromeric loci by [56], and the average value was retained. Centromeric positions were estimated for areas with low recombination rates and gene density. For each locus, the estimation of PP was performed as follows:

Considering the tetrasomic parameter ($\tau = 1 - PP$) that defines the proportion of gametes resulting from random meiotic chromosomal pairing used in several previous studies of meiotic behavior of tetraploid citrus [18,55–57,67,102] and according to [56], for centromeric loci, the probability of the three possible diploid gametes produced for an A1A1A2A2 duplex marker can be calculated using the following formulae:

$$P(A_1A_1) = \tau/6 \qquad (3)$$

$$P(A_2A_2) = \tau/6 \qquad (4)$$

$$P(A_1A_2) = 1 - \tau/3 \qquad (5)$$

The probability of obtaining the observed gamete population as a function of $\tau$ is:

$$L(\tau) = C \times (P(A_1A_1)^{x_1}) \times (P(A_2A_2)^{x_2}) \times (P(A_1A_2)^{x_3}) = C \times (1/6\tau)^{x_1+x_2} \times (1-1/3\tau)^{x_3} \qquad (6)$$
where C is a combinatory coefficient constant for the observed data and x1, x2, and x3 are the number of observed A₁A₁, A₂A₂, and A₁A₂ gametes. To estimate τ by a maximum likelihood approach, for each marker, we analyzed the L(τ) functions with τ values varying from 0 to 1 with a 0.005 interval. The estimated τ value was the one maximizing L(τ). PP was then calculated as 1 − τ.

R Studio V3.5.3 environment (R Studio Team, 2016, R Core Team, 2018, Vienna, Austria) was used for statistical analysis for PP and PHR. QQplot and a Shapiro-Wilk test were performed to test the normal distributions of PHR and PP. Normal distribution was confirmed for PP but rejected for PHR (*p*-value < 0.05). Analyses of variance (ANOVA) were performed to determine significant PP differences between chromosomes using a Newman-Keuls test for mean comparisons. A non-parametric Kruskal-Wallis test was used, followed by a Wilcoxon test, in order to study significant PHR differences between chromosomes.

4.10. *Analysis of the Deviation from Expected Gametic Segregation under a Tetrasomic Model*

The same matrix used for genetic mapping was used to study the deviation from
expected gametic segregation under a tetrasomic model without double reduction all along
the genome for each parent. The *p*-values for the Chi2 test according to the tetrasomic
theoretical frequency for each possible gamete (for a duplex locus A₁A₁A₂A₂: A₁A₁ = 1/6;
A₂A₂ = 1/6; and A₁A₂ = 4/6) were computed with Excel, and we used the approach
proposed by [103] to limit the False Discovery Rate (FDR) in multiple testing; the procedure
was performed according to the method of [61] with a *q*-value threshold of 0.05. The results
were visualized in a Circos plot.

4.11. *Inheritance of Candidate Genes in the Hybrid Population*

Loci associated with *Citrus trizteza virus* (CTV) resistance [104], nematode resistance (Tyr1) [8], and Huanglongbing (HLB) tolerance [59] were previously localized in *P. trifoliata*. In their study, [58] identified candidate genes in the *Poncirus* reference genome associated with these traits. The inheritance of these candidate genes was investigated in the tetraploid hybrid rootstock population. For this purpose, the positions of the genes were located on the genome, and the contribution of the *Poncirus* ancestor to this position was defined (zero, one, or two doses of *Poncirus* were administered to each hybrid).

4.12. *Root System Architecture Phenotyping*

A phenotypic analysis of the root architecture of the hybrids was performed on four-month-old cuttings. The analyses were performed on 100 genotypes with 6 replicates per genotype. Photographs of the root system of each individual were taken and analyzed using ImageJ software [105] to determine the total root area, length of the longest root, number of secondary roots, and root diameter. Statistical analyses were then performed in R using the R Studio V3.5.3 environment (R Studio Team, 2016, R Core Team, 2018, Vienna, Austria).

4.13. *Quantitative Trait Locus Analysis*

We set up a strategy to evaluate the association between the allelic dose of the ances-
tors and the phenotypic information collected for the root system architecture. We took
advantage of the fact that Volkamer lemon and citrumelo were fully heterozygous all along
the genome for *C. medica* and *P. trifoliata*, respectively, to perform an association study
between the doses (0, 1, or 2) of these two ancestors in the two gamete populations and
the phenotypic traits of the tetraploid hybrids (root surface, root length, root diameter, and
the number of secondary roots). We used the ancestor doses previously estimated with
TraceAncestor in the successive windows of 20 ancestor diagnostic SNPs. We considered
these window markers (positioned at the average value of the positions of the consid-
ered 20 DSNPs) to perform the association study using the GLM (General Linear Model)
procedure under the default settings in TASSEL 5.0 [106].
**5. Conclusions**

In this study, we analyzed GBS data of a tetraploid hybrid rootstock population of 103 progenies from the cross between Swingle citrumelo and Volkamer lemon. Estimation of allelic doses of *C. medica* and *P. trifoliata* ancestors in each tetraploid hybrid allowed the deduction of the phylogenomic structure of gametes of both Volkamer lemon and Swingle citrumelo, respectively. Two high-density recombination landscapes and a meiotic behavior analysis of tetraploid Swingle citrumelo and Volkamer lemon were established. Concerning the *Poncirus* diploid genetic map previously published, the recombination rate was lower compared to tetraploid Swingle citrumelo and Volkamer lemon. This can be mainly explained by the high preferential chromosome pairing detected. The measured meiotic parameter is defined primarily as an intermediate inheritance with a high disomic tendency for Swingle citrumelo. In Volkamer lemon, a more complicated meiotic inheritance was highlighted with very different patterns between chromosomes, including tetrasomic segregation, intermediate inheritance with tetrasomic tendency, intermediate inheritance, and intermediate inheritance with disomic tendency. Significant variability among progenies was found for root architecture traits, but the substantial reduction of interspecific recombination by the preferential chromosome pairing affected the ability to detect QTL for root architecture. Conversely, the segregation pattern with high preferential pairing appears beneficial for rootstock breeding at the tetraploid level. Indeed, there is a high level of parental interspecific heterozygosity restitution, and the lower apparent recombination rate limits allele shuffling. Thus, a large part of the genetic value of the parents is transmitted to the descendants. In that way, we found that genome regions present in the *Poncirus* genome involved in CTV and nematode resistance and Huanglongbing tolerance had a strong level of inheritance in the offspring. Therefore, we can conclude that the “tetrazyg” strategy with interspecific doubled-diploid parents enhances the transfer of the genetic value selected at the level of initial diploid rootstock. The tetrazyg method also allows the generation of some tetraploid hybrids with vigorous root systems contrary to doubled-diploid rootstocks. It is, therefore, highly promising for rootstock breeding. Our tetrazyg progenies were recently propagated by cutting with the perspective to perform field evaluation of agronomic performances and evaluation for HLB tolerance.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants12081630/s1, Figure S1: Links between the position of markers on the (A) tetraploid Swingle citrumelo and (B) tetraploid Volkamer lemon genetic map (LGs) and on the chromosome assembly (Chr) of *P. trifoliata* genome (Peng et al. 2020, [58]); Figure S2: Marey map plot of the nine linkage groups of (A) Swingle citrumelo and (B) Volkamer lemon compared with the *P. trifoliata* reference genome (Peng et al. 2020, [58]); Figure S3: Graphical representations showing the meiotic behavior of the tetraploid Swingle citrumelo and tetraploid Volkamer lemon. PHR: parental heterozygosity restitution; PP: preferential pairing. Figure S4: Photographs showing root system diversity in three 4× Volkamer lemon × 4× Swingle citrumelo hybrids (a, b, c) and in 4× Volkamer lemon (d). Figure S5: PCA showing correlation between measured root traits in the 4× Volkamer lemon × 4× Swingle citrumelo population; Figure S6: Boxplots showing the *C. medica* and *P. trifoliata* ancestors dose effect on two root traits (number of secondary roots and root surface) for four significant markers of chromosome 1.

**Author Contributions:** Conceptualization, P.O.; methodology and GBS data production, P.M., A.P., A.L. and M.M.; investigation, L.C., A.D., B.H. and P.O.; tetrazyg progeny production, Y.F.; plant management and phenotyping, S.B. and L.C.; data curation, A.D. and L.C.; formal analysis, L.C., A.D., and P.O.; writing—original draft preparation, L.C. and P.O.; writing—review and editing, L.C., P.O., B.H. and R.M.; visualization, L.C. and A.D.; supervision, P.O.; funding acquisition, R.M., P.O. and A.L.; project administration, A.L., R.M. and P.O. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Feder-Guadeloupe region “Cavalbio” project, the FRANCE GENOMIQUE “Dynamo” project, and the Région Guadeloupe for Leny Calvez’s doctoral fellowship.
**Data Availability Statement:** The Illumina Hiseq sequencing raw data from the GBS library are available in the NCBI SRA (Sequence Read Archive) under the accession number PRJNA898312.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

## References

1. Shafqat, W.; Naqvi, S.A.; Maqbool, R.; Haider, M.S.; Jaskani, M.J.; Khan, I.A. Climate Change and Citrus. *Citrus* **2021**, 147. [CrossRef]

2. Swingle, W.T.; Reece, P.C. The Botany of Citrus and Its Wild Relatives. In *The Citrus, Industry*; Reuther, W., Webber, H.J., Batchelor, L.D., Eds.; University of California Press: Berkeley, CA, USA, 1967.

3. Calvez, L.; Dereeper, A.; Mournet, P.; Froelicher, Y.; Bruyère, S.; Morillon, R.; Ollitrault, P. Intermediate Inheritance with Disomic Tendency in Tetraploid Intergeneric *Citrus × Poncirus* Hybrids Enhances the Efficiency of Citrus Rootstock Breeding. *Agronomy* **2020**, *10*, 1961. [CrossRef]

4. Curk, F.; Ollitrault, F.; Garcia-Lor, A.; Luro, F.; Navarro, L.; Ollitrault, P. Phylogenetic Origin of Limes and Lemons Revealed by Cytoplasmic and Nuclear Markers. *Ann. Bot.* **2016**, *117*, 565–583. [CrossRef] [PubMed]

5. Ahmed, D.; Comte, A.; Curk, F.; Costantino, G.; Luro, F.; Dereeper, A.; Mournet, P.; Froelicher, Y.; Ollitrault, P. Genotyping by Sequencing Can Reveal the Complex Mosaic Genomes in Gene Pools Resulting from Reticulate Evolution: A Case Study in Diploid and Polyploid Citrus. *Ann. Bot.* **2019**, *123*, 1231–1251. [CrossRef] [PubMed]

6. Bowman, K.D.; Joubert, J. Chapter 6—Citrus Rootstocks. In *The Genus Citrus*; Talon, M., Caruso, M., Gmitter, F.G., Eds.; Woodhead Publishing: Duxford, UK, 2020; pp. 105–127, ISBN 978-0-12-812163-4.

7. Spiegel-Roy, P.; Goldschmidt, E.E. The Biology of Horticultural Crops. In *The Biology of Citrus*; Cambridge University Press: Cambridge, UK, 1996; ISBN 978-0-521-33321-4.

8. Ling, P.; Duncan, L.W.; Deng, Z.; Dunn, D.; Hu, X.; Huang, S.; Gmitter Jr, F.G. Inheritance of Citrus Nematode Resistance and Its Linkage with Molecular Markers. *Theor. Appl. Genet.* **2000**, *100*, 1010–1017. [CrossRef]

9. Gmitter, F.G.; Xiao, S.Y.; Huang, S.; Hu, X.L.; Garnsey, S.M.; Deng, Z. A Localized Linkage Map of the *Citrus Tristeza Virus* Resistance Gene Region. *TAG Theor. Appl. Genet. Theor. Angew. Genet.* **1996**, *92*, 688–695. [CrossRef]

10. Li, T.; Wu, Q.-S.; Kuča, K.; Rahman, M.M. Responses of Four Citrus Plants to Phytophthora-Induced Root Rot. *Sains Malays.* **2018**, *47*, 1693–1700. [CrossRef]

11. Zhang, C.; Lang, P.; Dane, F.; Ebel, R.C.; Singh, N.K.; Locy, R.D.; Dozier, W.A. Cold Acclimation Induced Genes of Trifoliate Orange (*Poncirus Trifoliata*). *Plant Cell Rep.* **2005**, *23*, 764–769. [CrossRef]

12. Folimonova, S.Y.; Robertson, C.J.; Garnsey, S.M.; Gowda, S.; Dawson, W.O. Examination of the Responses of Different Genotypes of Citrus to Huanglongbing (Citrus Greening) under Different Conditions. *Phytopathology* **2009**, *99*, 1346–1354. [CrossRef]

13. Alves, M.N.; Lopes, S.A.; Raiol-Junior, L.L.; Wulff, N.A.; Girardi, E.A.; Ollitrault, P.; Pena, L. Resistance to "*Candidatus* Liberibacter Asiaticus," the Huanglongbing Associated Bacterium, in Sexually and/or Graft-Compatible Citrus Relatives. *Front. Plant Sci.* **2021**, *11*, 617664. [CrossRef]

14. Soratto, T.A.T.; Curtolo, M.; Marengo, S.; Dezotti, A.L.; Lima, R.P.M.; Gazaffi, R.; Machado, M.A.; Cristofani-Yaly, M. QTL and EQTL Mapping Associated with Host Response to *Candidatus* Liberibacter Asiaticus in Citrandarins. *Trop. Plant Pathol.* **2020**, *45*, 626–645. [CrossRef]

15. Albrecht, U.; McCollum, G.; Bowman, K.D. Influence of Rootstock Variety on Huanglongbing Disease Development in Field-Grown Sweet Orange (*Citrus Sinensis* [L.] Osbeck) Trees. *Sci. Hortic.* **2012**, *138*, 210–220. [CrossRef]

16. Stover, E.; McCollum, G. Incidence and Severity of Huanglongbing and *Candidatus* Liberibacter Asiaticus Titer among Field-Infected Citrus Cultivars. *HortScience* **2011**, *46*, 1344–1348. [CrossRef]

17. George, J.; Lapointe, S.L. Host-Plant Resistance Associated with *Poncirus Trifoliata* Influence Oviposition, Development and Adult Emergence of *Diaphorina Citri* (Hemiptera: Liviidae). *Pest Manag. Sci.* **2019**, *75*, 279–285. [CrossRef]

18. Kamiri, M.; Stift, M.; Costantino, G.; Dambier, D.; Kabbage, T.; Ollitrault, P.; Froelicher, Y. Preferential Homologous Chromosome Pairing in a Tetraploid Intergeneric Somatic Hybrid (*Citrus Reticulata* + *Poncirus Trifoliata*) Revealed by Molecular Marker Inheritance. *Front. Plant Sci.* **2018**, *9*, 1557. [CrossRef]

19. Ben Yahmed, J.; Costantino, G.; Amiel, P.; Talon, M.; Ollitrault, P.; Morillon, R.; Luro, F. Diversity in the Trifoliate Orange Taxon Reveals Two Main Genetic Groups Marked by Specific Morphological Traits and Water Deficit Tolerance Properties. *J. Agric. Sci.* **2016**, *154*, 495–514. [CrossRef]

20. Khalid, M.F.; Hussain, S.; Anjum, M.A.; Ahmad, S.; Ali, M.A.; Ejaz, S.; Morillon, R. Better Salinity Tolerance in Tetraploid vs Diploid Volkamer Lemon Seedlings Is Associated with Robust Antioxidant and Osmotic Adjustment Mechanisms. *J. Plant Physiol.* **2020**, *244*, 153071. [CrossRef]

21. Khalid, M.F.; Morillon, R.; Anjum, M.A.; Ejaz, S.; Rao, M.J.; Ahmad, S.; Hussain, S. Volkamer Lemon Tetraploid Rootstock Transmits the Salt Tolerance When Grafted with Diploid Kinnow Mandarin by Strong Antioxidant Defense Mechanism and Efficient Osmotic Adjustment. *J. Plant Growth Regul.* **2022**, *41*, 1125–1137. [CrossRef]
22. Mouhaya, W.; Allario, T.; Costantino, G.; Andres, F.; Talon, M.; Luro, F.; Ollitrault, P.; Morillon, R. Citrus Tetraploid Rootstocks Are More Tolerant to Salt Stress than Diploid. Available online: https://agritrop.cirad.fr/547588/ (accessed on 13 September 2021).

23. Saleh, B.; Allario, T.; Dambier, D.; Ollitrault, P.; Morillon, R. Tetraploid Citrus Rootstocks Are More Tolerant to Salt Stress than Diploid. *Comptes. Rendus Biol.* **2008**, *331*, 703–710. [CrossRef]

24. Oustric, J.; Morillon, R.; Luro, F.; Herbette, S.; Lourkisti, R.; Giannettini, J.; Berti, L.; Santini, J. 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.* **2017**, *214*, 108–115. [CrossRef]

25. Oustric, J.; Quilichini, Y.; Morillon, R.; Herbette, S.; Luro, F.; Giannettini, J.; Berti, L.; Santini, J. Tetraploid Citrus Seedlings Subjected to Long-Term Nutrient Deficiency Are Less Affected at the Ultrastructural, Physiological and Biochemical Levels than Diploid Ones. *Plant Physiol. Biochem.* **2019**, *135*, 372–384. [CrossRef] [PubMed]

26. Oustric, J.; Herbette, S.; Morillon, R.; Giannettini, J.; Berti, L.; Santini, J. Influence of Rootstock Genotype and Ploidy Level on Common Clementine (*Citrus clementina* Hort. Ex Tan) Tolerance to Nutrient Deficiency. *Front. Plant Sci.* **2021**, *12*, 403. [CrossRef] [PubMed]

27. Allario, T.; Brumos, J.; Colmenero-Flores, J.M.; Iglesias, D.J.; Pina, J.A.; Navarro, L.; Talon, M.; Ollitrault, P.; Morillon, R. Tetraploid Rangpur Lime Rootstock Increases Drought Tolerance via Enhanced Constitutive Root Abscisic Acid Production. *Plant Cell Environ.* **2013**, *36*, 856–868. [CrossRef] [PubMed]

28. Oliveira, T.M.; Yahmed, J.B.; Dutra, J.; Maserti, B.E.; Talon, M.; Navarro, L.; Ollitraut, P.; da S. Gesteira, A.; Morillon, R. Better Tolerance to Water Deficit in Doubled Diploid 'Carrizo Citrange' Compared to Diploid Seedlings Is Associated with More Limited Water Consumption. *Acta Physiol. Plant.* **2017**, *39*, 204. [CrossRef]

29. Wei, T.; Wang, Y.; Xie, Z.; Guo, D.; Chen, C.; Fan, Q.; Deng, X.; Liu, J. Enhanced ROS Scavenging and Sugar Accumulation Contribute to Drought Tolerance of Naturally Occurring Autotetraploids in *Poncirus Trifoliata*. *Plant Biotechnol. J.* **2019**, *17*, 1394–1407. [CrossRef]

30. Sivager, G.; Calvez, L.; Bruyere, S.; Boisne-Noc, R.; Brat, P.; Gros, O.; Ollitrault, P.; Morillon, R. Specific Physiological and Anatomical Traits Associated With Polyploidy and Better Detoxification Processes Contribute to Improved Huanglongbing Tolerance of the Persian Lime Compared With the Mexican Lime. *Front. Plant Sci.* **2021**, *12*, 1343. [CrossRef]

31. Sivager, G.; Calvez, L.; Bruyere, S.; Boisne-Noc, R.; Hufnagel, B.; Cebrian-Torrejon, G.; Doménech-Carbó, A.; Gros, O.; Ollitrault, P.; Morillon, R. Better Tolerance to Huanglongbing Is Conferred by Tetraploid Swingle Citrumelo Rootstock and Is Influenced by the Ploidy of the Scion. *Front. Plant Sci.* **2022**, *13*, 14. [CrossRef]

32. Ollitrault, P.; Germanà, M.A.; Froelicher, Y.; Cuenca, J.; Aleza, P.; Morillon, R.; Grosser, J.W.; Guo, W. Ploidy Manipulation for Citrus Breeding, Genetics, and Genomics. In *The Citrus Genome*; Gentile, A., La Malfa, S., Deng, Z., Eds.; Compendium of Plant Genomes; Springer International Publishing: Cham, The Switzerland, 2020; pp. 75–105, ISBN 978-3-030-15308-3.

33. Rouiss, H.; Cuenca, J.; Navarro, L.; Ollitrault, P.; Aleza, P. Tetraploid Citrus Progenies Arising from FDR and SDR Unreduced Pollen in 4x X 2x Hybridizations. *Tree Genet. Genomes* **2017**, *13*, 10. [CrossRef]

34. Cuenca, J.; Froelicher, Y.; Aleza, P.; Juárez, J.; Navarro, L.; Ollitrault, P. 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* **2011**, *107*, 462–470. [CrossRef]

35. Cuenca, J.; Aleza, P.; Juárez, J.; García-Lor, A.; Froelicher, Y.; Navarro, L.; Ollitrault, P. Maximum-Likelihood Method Identifies Meiotic Restitution Mechanism from Heterozygosity Transmission of Centromeric Loci: Application in Citrus. *Sci. Rep.* **2015**, *5*, 9897. [CrossRef]

36. Aleza, P.; Juárez, J.; Cuenca, J.; Ollitrault, P.; Navarro, L. Recovery of Citrus Triploid Hybrids by Embryo Rescue and Flow Cytometry from 2x x 2x Sexual Hybridisation and Its Application to Extensive Breeding Programs. *Plant Cell Rep.* **2010**, *29*, 1023–1034. [CrossRef]

37. Aleza, P.; Froelicher, Y.; Schwarz, S.; Agustí, M.; Hernández, M.; Juárez, J.; Luro, F.; Morillon, R.; Navarro, L.; Ollitrault, P. Tetraploidization Events by Chromosome Doubling of Nucellar Cells Are Frequent in Apomictic Citrus and Are Dependent on Genotype and Environment. *Ann. Bot.* **2011**, *108*, 37–50. [CrossRef]

38. Cameron, J.W.; Frost, H.B. Genetics, Breeding and Nucellar Embryony. *Citrus Ind.* **1968**, *1*, 325–370.

39. Barrett, H.C. Colchicine-Induced Polyploidy in Citrus. *Bot. Gaz.* **1974**, *135*, 29–41. [CrossRef]

40. Aleza, P.; Juárez, J.; Ollitrault, P.; Navarro, L. Production of Tetraploid Plants of Non Apomictic Citrus Genotypes. *Plant Cell Rep.* **2009**, *28*, 1837–1846. [CrossRef]

41. Grosser, J.W.; Kainth, D.; Dutt, M. Production of Colchicine-Induced Autotetraploids in Pummelo (*Citrus grandis* Osbeck) through Indirect Organogenesis. *HortScience* **2014**, *49*, 944–948. [CrossRef]

42. Dambier, D.; Barantin, P.; Boulard, G.; Costantino, G.; Mournet, P.; Perdereau, A.; Morillon, R.; Ollitrault, P. Genomic Instability in Somatic Hybridization between Poncirus and Citrus Species Aiming to Create New Rootstocks. *Agriculture* **2022**, *12*, 134. [CrossRef]

43. Dambier, D.; Benyahia, H.; Pensabene-Bellavia, G.; Aka Kaçar, Y.; Froelicher, Y.; Belfalah, Z.; Lhou, B.; Handaji, N.; Printz, B.; Morillon, R.; et al. Somatic Hybridization for Citrus Rootstock Breeding: An Effective Tool to Solve Some Important Issues of the Mediterranean Citrus Industry. *Plant Cell Rep.* **2011**, *30*, 883–900. [CrossRef]

44. Grosser, J.W.; Ollitrault, P.; Olivares-Fuster, O. Somatic Hybridization in Citrus: An Effective Tool to Facilitate Variety Improvement. *Vitr. Cell. Dev. Biol.*—*Plant* **2000**, *36*, 434–449. [CrossRef]
45. Grosser, J.W.; Chandler, J.L. Somatic Hybridization of High Yield, Cold-Hardy and Disease Resistant Parents for Citrus Rootstock Improvement. *J. Hortic. Sci. Biotechnol.* **2000**, *75*, 641–644. [CrossRef]

46. Ollitrault, P.; Vanel, F.; Froelicher, Y.; Dambier, D. Creation of Triploid Citrus Hybrids by Electrofusion of Haploid and Diploid Protoplasts. Available online: https://agritrop.cirad.fr/477415/ (accessed on 13 September 2021).

47. Ollitrault, P.; Guo WenWu, G.W.; Grosser, J.W. Somatic Hybridization. In *Citrus Genetics, Breeding and Biotechnology*; Khan, I.A., Ed.; CABI: Wallingford, UK, 2007; pp. 235–260, ISBN 978-0-85199-019-4.

48. Guo, W.W.; Wu, R.C.; Cheng, Y.J.; Deng, X.X. Production and Molecular Characterization of Citrus Intergeneric Somatic Hybrids between Red Tangerine and Citrange. *Plant Breed.* **2007**, *126*, 72–76. [CrossRef]

49. Grosser, J.W.; Gmitter, F.G. Protoplast Fusion for Production of Tetraploids and Triploids: Applications for Scion and Rootstock Breeding in Citrus. *Plant Cell Tissue Organ Cult.* **(PCTOC)** **2011**, *104*, 343–357. [CrossRef]

50. Grosser, J.W.; Graham, J.H.; McCoy, C.W.; Hoyte, A.; Rubio, H.M.; Bright, D.B.; Chandler, J.L. Development of "Tetrazyg" Rootstocks Tolerant of the Diaprepes/Phytophthora Complex under Greenhouse Conditions. *Proc. Fla. State Hortic. Soc.* **2003**, *116*, 263–267.

51. Grosser, J.W.; Barthe, G.A.; Castle, B.; Gmitter, F.G.J.; Lee, O. The Development of Improved Tetraploid Citrus Rootstocks to Facilitate Advanced Production Systems and Sustainable Citriculture in Florida. *Acta Hortic.* **2015**, *1065*, 319–327. [CrossRef]

52. Grosser, J.W.; Graham, J.H.; Hoyte, A.; Rubio, H.M.; Bright, D.B.; Gmitter, J.; Chen, C.X.; Gmitter, F.G. Continued Development of Rootstocks Tolerant of the Phytophthora-Diaprepes Complex via Greenhouse Screening. *Proc. Fla. State Hortic. Soc.* **2007**, *120*, 103–109.

53. Stebbins, G.L. Types of Polyploids: Their Classification and Significance. In *Advances in Genetics*; Demerec, M., Ed.; Academic Press: Cambridge, MA, USA, 1947; Volume 1, pp. 403–429.

54. Chen, Z.J. Genetic and Epigenetic Mechanisms for Gene Expression and Phenotypic Variation in Plant Polyploids. *Annu. Rev. Plant Biol.* **2007**, *58*, 377–406. [CrossRef]

55. Stift, M.; Berenos, C.; Kuperus, P.; van Tienderen, P.H. Segregation Models for Disomic, Tetrasomic and Intermediate Inheritance in Tetraploids: A General Procedure Applied to Rorippa (Yellow Cress) Microsatellite Data. *Genetics* **2008**, *179*, 2113–2123. [CrossRef]

56. Aleza, P.; Cuenca, J.; Juárez, J.; Navarro, L.; Ollitrault, P. Inheritance in Doubled-Diploid Clementine and Comparative Study with SDR Unreduced Gametes of Diploid Clementine. *Plant Cell Rep.* **2016**, *35*, 1573–1586. [CrossRef]

57. Rouiss, H.; Bakry, F.; Froelicher, Y.; Navarro, L.; Aleza, P.; Ollitrault, P. 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.* **2018**, *121*, 571–585. [CrossRef]

58. Peng, Z.; Bredeson, J.V.; Wu, G.A.; Shu, S.; Rawat, N.; Du, D.; Parajuli, S.; Yu, Q.; You, Q.; Rokhsar, D.S.; et al. A Chromosome-Scale Reference Genome of Trifoliate Orange (*Poncirus Trifoliata*) Provides Insights into Disease Resistance, Cold Tolerance and Genome Evolution in Citrus. *Plant J.* **2020**, *104*, 1215–1232. [CrossRef]

59. Huang, M.; Roosevelt, M.L.; Yu, Q.; Du, D.; Yu, Y.; Zhang, Y.; Deng, Z.; Stover, E.; Gmitter, F.G.J. Construction of High-Density Genetic Maps and Detection of QTLs Associated With Huanglongbing Tolerance in Citrus. *Front. Plant Sci.* **2018**, *9*, 1694. [CrossRef]

60. Aleza, P.; Cuenca, J.; Hernández, M.; Juárez, J.; Navarro, L.; Ollitrault, P. 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.* **2015**, *15*, 80. [CrossRef]

61. Storey, J.D. A Direct Approach to False Discovery Rates. *J. R. Stat. Soc. Ser. B Stat. Methodol.* **2002**, *64*, 479–498. [CrossRef]

62. De Jong, W.S.; De Jong, D.M.; Bodis, M. A Fluorogenic 5' Nuclease (TaqMan) Assay to Assess Dosage of a Marker Tightly Linked to Red Skin Color in Autotetraploid Potato. *TAG Theor. Appl. Genet. Theor. Angew. Genet.* **2003**, *107*, 1384–1390. [CrossRef]

63. Esselink, G.D.; Nybom, H.; Vosman, B. Assignment of Allelic Configuration in Polyploids Using the MAC-PR (Microsatellite DNA Allele Counting-Peak Ratios) Method. *TAG Theor. Appl. Genet. Theor. Angew. Genet.* **2004**, *109*, 402–408. [CrossRef]

64. Ferrante, S.P.; Lucretti, S.; Reale, S.; De Patrizio, A.; Abbate, L.; Tusa, N.; Scarano, M.-T. Assessment of the Origin of New Citrus Tetraploid Hybrids (2n = 4x) by Means of SSR Markers and PCR Based Dosage Effects. *Euphytica* **2010**, *173*, 223–233. [CrossRef]

65. Cuenca, J.; Aleza, P.; Navarro, L.; Ollitrault, P. Assignment of SNP Allelic Configuration in Polyploids Using Competitive Allele-Specific PCR: Application to Citrus Triploid Progeny. *Ann. Bot.* **2013**, *111*, 731–742. [CrossRef]

66. Rouiss, H.; Cuenca, J.; Navarro, L.; Ollitrault, P.; Aleza, P. Unreduced Megagametophyte Production in Lemon Occurs via Three Meiotic Mechanisms, Predominantly Second-Division Restitution. *Front. Plant Sci.* **2017**, *8*, 1211. [CrossRef]

67. Ahmed, D.; Curk, F.; Evrard, J.C.; Froelicher, Y.; Ollitrault, P. Preferential Disomic Segregation and *C. micrantha/C. medica* Interspecific Recombination in Tetraploid 'Giant Key' Lime; Outlook for Triploid Lime Breeding. *Front. Plant Sci.* **2020**, *11*, 939. [CrossRef]

68. Garavello, M.; Cuenca, J.; Garcia-Lor, A.; Ortega, N.; Navarro, L.; Ollitrault, P.; Aleza, P. Male and Female Inheritance Patterns in Tetraploid 'Moncada' Mandarin. *Plant Cell Rep.* **2020**, *39*, 335–349. [CrossRef]

69. Oueslati, A.; Salhi-Hannachi, A.; Luro, F.; Vignes, H.; Mournet, P.; Ollitrault, P. Genotyping by Sequencing Reveals the Interspecific *C. maxima/C. reticulata* Admixture along the Genomes of Modern Citrus Varieties of Mandarins, Tangors, Tangelos, Orangelos and Grapefruits. *PLoS ONE* **2017**, *12*, e0185618. [CrossRef] [PubMed]

70. Wu, G.A.; Terol, J.; Ibanez, V.; López-García, A.; Pérez-Román, E.; Borredá, C.; Domingo, C.; Tadeo, F.R.; Carbonell-Caballero, J.; Alonso, R.; et al. Genomics of the Origin and Evolution of Citrus. *Nature* **2018**, *554*, 311–316. [CrossRef] [PubMed]
71. Xie, K.-D.; Xia, Q.-M.; Wang, X.-P.; Liang, W.-J.; Wu, X.-M.; Grosser, J.W.; Guo, W.-W. 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* **2015**, *11*, 112. [CrossRef]

72. Jeridi, M.; Perrier, X.; Rodier-Goud, M.; Ferchichi, A.; D'Hont, A.; Bakry, F. Cytogenetic Evidence of Mixed Disomic and Polysomic Inheritance in an Allotetraploid (AABB) Musa Genotype. *Ann. Bot.* **2012**, *110*, 1593–1606. [CrossRef]

73. Chen, C.; Bowman, K.D.; Choi, Y.A.; Dang, P.M.; Rao, M.N.; Huang, S.; Soneji, J.R.; McCollum, T.G.; Gmitter, F.G. EST-SSR Genetic Maps for *Citrus sinensis* and *Poncirus trifoliata*. *Tree Genet. Genomes* **2008**, *4*, 1–10. [CrossRef]

74. Curtolo, M.; Soratto, T.A.T.; Gazaffi, R.; Takita, M.A.; Machado, M.A.; Cristofani-Yaly, M. High-Density Linkage Maps for *Citrus sunki* and *Poncirus trifoliata* Using DArTseq Markers. *Tree Genet. Genomes* **2017**, *14*, 5. [CrossRef]

75. Bernet, G.P.; Fernandez-Ribacoba, J.; Carbonell, E.A.; Asins, M.J. Comparative Genome-Wide Segregation Analysis and Map Construction Using a Reciprocal Cross Design to Facilitate Citrus Germplasm Utilization. *Mol. Breed.* **2010**, *25*, 659–673. [CrossRef]

76. Xu, Y.-Y.; Liu, S.-R.; Gan, Z.-M.; Zeng, R.-F.; Zhang, J.-Z.; Hu, C.-G. High-Density Genetic Map Construction and Identification of QTLs Controlling Leaf Abscission Trait in Poncirus Trifoliata. *Int. J. Mol. Sci.* **2021**, *22*, 5723. [CrossRef]

77. Ollitrault, P.; Ahmed, D.; Costantino, G.; Evrard, J.-C.; Cardi, C.; Mournet, P.; Perdereau, A.; Froelicher, Y. Segregation Distortion for Male Parents in High Density Genetic Maps from Reciprocal Crosses between Two Self-Incompatible Cultivars Confirms a Gametophytic System for Self-Incompatibility in Citrus. *Agriculture* **2021**, *11*, 379. [CrossRef]

78. Ollitrault, P.; Terol, J.; Chen, C.; Federici, C.T.; Lotfy, S.; Hippolyte, I.; Ollitrault, F.; Bérard, A.; Chauveau, A.; Cuenca, J.; et al. A Reference Genetic Map of *C. clementina* Hort. Ex Tan.; Citrus Evolution Inferences from Comparative Mapping. *BMC Genom.* **2012**, *13*, 593. [CrossRef]

79. Parker, J.S.; Palmer, R.W.; Whitehorn, M.A.F.; Edgar, L.A. Chiasma Frequency Effects of Structural Chromosome Change. *Chromosoma* **1982**, *85*, 673–686. [CrossRef]

80. Chambers, S.R.; Hunter, N.; Louis, E.J.; Borts, R.H. The Mismatch Repair System Reduces Meiotic Homeologous Recombination and Stimulates Recombination-Dependent Chromosome Loss. *Mol. Cell. Biol.* **1996**, *16*, 6110–6120. [CrossRef]

81. Liharska, T.; Wordragen, M.; Kammen, A.; Zabel, P.; Koornneef, M. Tomato Chromosome 6: Effect of Alien Chromosomal Segments on Recombinant Frequencies. *Genome* **1996**, *39*, 485–491. [CrossRef]

82. Chetelat, R.T.; Meglic, V.; Cisneros, P. A Genetic Map of Tomato Based on BC(1) *Lycopersicon Esculentum* x *Solanum Lycopersicoides* Reveals Overall Synteny but Suppressed Recombination between These Homeologous Genomes. *Genetics* **2000**, *154*, 857–867. [CrossRef]

83. Opperman, R.; Emmanuel, E.; Levy, A.A. The Effect of Sequence Divergence on Recombination between Direct Repeats in Arabidopsis. *Genetics* **2004**, *168*, 2207–2215. [CrossRef]

84. Li, L.; Jean, M.; Belzile, F. The Impact of Sequence Divergence and DNA Mismatch Repair on Homeologous Recombination in Arabidopsis. *Plant J. Cell Mol. Biol.* **2006**, *45*, 908–916. [CrossRef]

85. Pecinka, A.; Fang, W.; Rehmsmeier, M.; Levy, A.A.; Mittelsten Scheid, O. Polyploidization Increases Meiotic Recombination Frequency in Arabidopsis. *BMC Biol.* **2011**, *9*, 1–7. [CrossRef]

86. Pelé, A.; Falque, M.; Trotoux, G.; Eber, F.; Nègre, S.; Gilet, M.; Huteau, V.; Lodé, M.; Jousseaume, T.; Dechaumet, S.; et al. Amplifying Recombination Genome-Wide and Reshaping Crossover Landscapes in Brassicas. *PLoS Genet.* **2017**, *13*, e1006794. [CrossRef]

87. Elshire, R.J.; Glaubitz, J.C.; Sun, Q.; Poland, J.A.; Kawamoto, K.; Buckler, E.S.; Mitchell, S.E. A Robust, Simple Genotyping-by-Sequencing (GBS) Approach for High Diversity Species. *PLoS ONE* **2011**, *6*, e19379. [CrossRef]

88. Sonah, H.; Bastien, M.; Iquira, E.; Tardivel, A.; Légaré, G.; Boyle, B.; Normandeau, É.; Laroche, J.; Larose, S.; Jean, M.; et al. An Improved Genotyping by Sequencing (GBS) Approach Offering Increased Versatility and Efficiency of SNP Discovery and Genotyping. *PLoS ONE* **2013**, *8*, e54603. [CrossRef]

89. Herten, K.; Hestand, M.S.; Vermeesch, J.R.; Van Houdt, J.K.J. GBSX: A Toolkit for Experimental Design and Demultiplexing Genotyping by Sequencing Experiments. *BMC Bioinform.* **2015**, *16*, 73. [CrossRef] [PubMed]

90. Martin, M. Cutadapt Removes Adapter Sequences from High-Throughput Sequencing Reads. *EMBnet. J.* **2011**, *17*, 10–12. [CrossRef]

91. Garsmeur, O.; Droc, G.; Antonise, R.; Grimwood, J.; Potier, B.; Aitken, K.S.; Jenkins, J.; Martin, G.; Charron, C.; Hervouet, C.; et al. A Mosaic Monoploid Reference Sequence for the Highly Complex Genome of Sugarcane. *Nat. Commun.* **2018**, *9*, 2638. [CrossRef] [PubMed]

92. Li, H.; Durbin, R. Fast and Accurate Short Read Alignment with Burrows-Wheeler Transform. *Bioinformatics* **2009**, *25*, 1754–1760. [CrossRef]

93. Dereeper, A.; Homa, F.; Andres, G.; Sempere, G.; Sarah, G.; Hueber, Y.; Dufayard, J.-F.; Ruiz, M. SNiPlay3: A Web-Based Application for Exploration and Large Scale Analyses of Genomic Variations. *Nucleic Acids Res.* **2015**, *43*, W295–W300. [CrossRef]

94. Nei, M. Analysis of Gene Diversity in Subdivided Populations. *Proc. Natl. Acad. Sci. USA* **1973**, *70*, 3321–3323. [CrossRef]

95. Perrier, X.; Jacquemoud-Collet, J.P. DARwin Software. 2006. Available online: http://darwin.cirad.fr/Darwin (accessed on 6 September 2022).

96. Adler, D.; Nenadí, O.; Zucchini, W. Rgl: A r-Library for 3d Visualization with Opengl. In Proceedings of the 35th Symposium of the Interface: Computing Science and Statistics, Salt Lake City, UT, USA, 12–15 March 2003; Volume 35, pp. 419–429.
97. Van Os, H.; Stam, P.; Visser, R.G.F.; van Eck, H.J. SMOOTH: A Statistical Method for Successful Removal of Genotyping Errors from High-Density Genetic Linkage Data. *TAG Theor. Appl. Genet. Theor. Angew. Genet.* **2005**, *112*, 187–194. [CrossRef]

98. Grandke, F.; Ranganathan, S.; van Bers, N.; de Haan, J.R.; Metzler, D. PERGOLA: Fast and Deterministic Linkage Mapping of Polyploids. *BMC Bioinform.* **2017**, *18*, 12. [CrossRef]

99. Krzywinski, M.; Schein, J.; Birol, İ.; Connors, J.; Gascoyne, R.; Horsman, D.; Jones, S.J.; Marra, M.A. Circos: An Information Aesthetic for Comparative Genomics. *Genome Res.* **2009**, *19*, 1639–1645. [CrossRef]

100. Rasche, H.; Hiltemann, S. Galactic Circos: User-Friendly Circos Plots within the Galaxy Platform. *GigaScience* **2020**, *9*, giaa065. [CrossRef]

101. Siberchicot, A.; Bessy, A.; Guéguen, L.; Marais, G.A. MareyMap Online: A User-Friendly Web Application and Database Service for Estimating Recombination Rates Using Physical and Genetic Maps. *Genome Biol. Evol.* **2017**, *9*, 2506–2509. [CrossRef]

102. Kamiri, M.; Stift, M.; Srairi, I.; Costantino, G.; Moussadik, A.E.; Hmyene, A.; Bakry, F.; Ollitrault, P.; Froelicher, Y. Evidence for Non-Disomic Inheritance in a Citrus Interspecific Tetraploid Somatic Hybrid between C. *Reticulata* and C. *Limon* Using SSR Markers and Cytogenetic Analysis. *Plant Cell Rep.* **2011**, *30*, 1415–1425. [CrossRef]

103. Benjamini, Y.; Hochberg, Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. *J. R. Stat. Soc. Ser. B Methodol.* **1995**, *57*, 289–300. [CrossRef]

104. Yang, Z.-N.; Ye, X.-R.; Molina, J.; Roosevelt, M.L.; Mirkov, T.E. Sequence Analysis of a 282-Kilobase Region Surrounding the Citrus Tristeza Virus Resistance Gene (Ctv) Locus in *Poncirus Trifoliata* L. *Raf. Plant Physiol.* **2003**, *131*, 482–492. [CrossRef]

105. Schneider, C.A.; Rasband, W.S.; Eliceiri, K.W. NIH Image to ImageJ: 25 Years of Image Analysis. *Nat. Methods* **2012**, *9*, 671–675. [CrossRef]

106. Bradbury, P.J.; Zhang, Z.; Kroon, D.E.; Casstevens, T.M.; Ramdoss, Y.; Buckler, E.S. TASSEL: Software for Association Mapping of Complex Traits in Diverse Samples. *Bioinformatics* **2007**, *23*, 2633–2635. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.