The study was conducted using an F₂ population derived from the cross between the cultivated variety Fleur 11 and the synthetic allotetraploid ISATGR 278-18 (Arachis duranensis × Arachis batizocoi)^4x. Fleur 11 is a Spanish type with erect growth habit widely cultivated in West Africa. ISATGR 278-18 was developed and kindly provided by ICRISAT-India (Mallikarjuna et al., 2011). The synthetic allotetraploid combines the AA genome of A. duranensis (ICG 8138; 2n = 2x = 20), a close wild relative and one of the most probable ancestors of A. hypogaea, and the KK genome of A. batizocoï (ICG 13160; 2n = 2x = 20), a wild relative taxa that was reported to pair with the B genome of the cultivated species during meiosis (Burow et al., 2001; Mallikarjuna et al., 2011; Leal-Bertioli et al., 2014).

ISATGR 278-18 was reported to have a normal chromosome configuration with 20 bivalents (Mallikarjuna et al., 2011). However, several cycles of self-pollination were performed at CERAAS prior to hybridization with Fleur 11. Five plants of the synthetic allotetraploid were used as male to cross with five plants of Fleur 11 used as female. The F₁ plants were differentiated from plants derived from self-pollination of Fleur 11 using morphological traits (dark green leaves and procumbent growth habit). The F₂ progeny were produced from the self-pollination of 15 F₁ plants. All crosses were performed in plastic pots under greenhouse conditions at the Centre d’Etudes Régional pour l’Amélioration de l’Adaptation à la Sécheresse (CERAAS) in Senegal during 2011 and 2013.

Genomic DNA of both parents and F₂ progeny was extracted from young leaves according to the MATAB protocol as described by Foncéka et al. (2009). Polymorphisms were assessed in the parents using 602 primer pairs, mainly selected from a previous study (Foncéka et al., 2009). The parents and F₂ progeny were genotyped with polymorphic SSR markers. PCR was carried out in 96-wellplates in a total volume of 10 μl consisting of 5 μl of 5 ng/μl of the DNA template and 5 μl of a mixture of 0.1 μM of each SSR primer, 0.2 mM of each dNTP, 1X PCR buffer, 2.5 mM MgCl2, 0.1U/μl of Taq polymerase and 0.1 mM of IR700 or IR800-labeled M13 primer (MWG Germany) for fluorescence detection of SSR amplicons. A forward primer pair was labeled with a fluorescent (LI-COR Biosciences). Reactions were performed in an Eppendorf Mastercycler epgradient thermocycler. The PCR products were separated by electrophoresis run at a constant 95 W for 1-2 h in a DNA Sequencer (LI-COR 4300 DNA Analyzer, Lincoln, NE, USA).

Scoring of the SSR bands was performed visually on electrophoresis profiles using the application Jelly 2.017b (Rami, unpublished). For a codominant marker in an F₂ population, we denoted “A” as the genotype of cultivated parent, “B” as the genotype of wild parent, and “H” as the genotype of their heterozygous hybrids. When one marker was dominant, the two non-separated genotypes “H” and “A” were denoted “D” and the two non-separated genotypes “H” and “B” were denoted “C”. For each SSR, the sub-genomic origin (A or K) of the wild alleles was assigned by comparing them to the alleles of the wild diploid progenitors. The sub-genomic origin of the cultivated allele was inferred by analyzing its co-inheritance with the wild alleles in the F₂ progeny. Loci were suffixed by <A> or <B> for differentiating subgenomic origin. In all scenarios, missing data were scored as “x” and the unexpected null bands were scored as “N”.

Based on the genotyping data for each SSR marker, we analyzed the genotype of the F₂ progeny and deduced the allelic constitution of the gamete produced by the F₁ hybrid. The genome of the cultivated species was noted “A₁A₁BB” and that of the wild species “A₂A₂KK”. Chromosomes “A₁” and “A₂” and “B” and “K” are homologous, while chromosomes “A₁” and “B”, “A₁” and “K”, “A₂” and “B”, and “A₂” and “K” are homeologous.

When recombination occurred only between homologous genomes, for SSRs that amplified only one locus on a given sub-genome, two segregating bands were expected and the segregation ratio in the F₂ progeny was 1A₁A₁: 2A₁A₂: 1A₂A₂ or 1BB: 2BK: 1KK. For SSR that amplified the two homeologous genomes up to four segregating bands were expected. The segregation ratio in the F₂ progeny is shown in Table 1. In this case, the cultivated and wild parents produced each one type of gamete (A₁B and A₂K, respectively) and the F₁ hybrid (A₁A₂BK) would produce 4 types of gametes (A₁B, A₁K, A₂B and A₂K).

When recombination occurred between homeologous genomes, for SSRs that amplified only one locus on a given sub-genome, a fourth ‘null’ genotype is expected to appear in addition to the three expected genotypes in an F₂ population. For SSRs that amplified the two loci, one on each of the homeologous genomes, 11 phenotypic bands are expected in the gel (Table 1). The segregation ratio in the F₂ progeny depends on the genomic constitution of the F₁ hybrids. For F₁ with normal genomic constitution (without homeologous recombination during parental meiosis), the expected segregation ratio in the F₂ progeny is shown in Table 1. It is not straightforward to estimate the genotypic frequencies in F₂ progeny derived from “aberrant” F₁ (resulting from homeologous recombination events during parental meiosis). In contrast to the homologous pairings, two specific phenotypes, A₁A₂ and BK, are observed in homeologous pairings. In this study, we use the terminology “homeologous-recombinant genotypes” for the progeny that carried these two peculiar phenotypes.

In the electrophoresis profiles, the band intensities of the homeologous-recombinant progeny were analyzed to determine the number of copies of each allele. The band intensities were determined using ImageJ version 1.46 (Ferreira and Rasband, 2012), and the dosage ratio between bands was compared to the relationships expected between alleles in the hypothetical configurations. For example, for the homeologous-recombinant progeny with A₁A₂ band-phenotype on the gels, three genotypes are possible with ratios of allelic combinations of 3:1 (A₁A₁A₁A₂), 1:1 (A₁A₁A₂A₂) or 1:3 (A₁A₂A₂A₂). The genotype was then determined based on allelic dosage. The F₂ progeny with four different alleles (i.e., A₁A₂BK) were used as a reference for determination, as these progeny involve single-copy alleles.

The segregation of codominant and dominant loci was compared to the segregation ratios expected under disomic inheritance (1:2:1 for codominant loci and 3:1 for dominant loci), using the chi-square test of the software “Calculation for the chi-square test” (Preacher, 2001). Data scored as “N” (null phenotype) as well as missing data were not considered in the analysis. Loci that deviated significantly (P < 0.05) from the theoretically expected ratios were considered distorted and were represented by an asterisk on the genetic map.

Linkage analysis was performed using the MapDisto software (Lorieux, 2012) using only the co-dominant loci. The linkage map was constructed in several steps. In the first step, markers that showed distorted segregation (P < 0.05) and those with unexpected null data were excluded. The non-distorted loci were grouped into LGs using a minimum LOD of 3, a maximum recombination frequency r of 0.3 and the Kosambi mapping function (Kosambi, 1943). The order of markers within each linkage group (LG) was determined using the “order” and “ripple” commands. In a second step, the distorted loci and loci with unexpected null data were progressively added into the established LGs if their presence did not significantly affect the marker order. The position of the loci with unexpected null data was adjusted iteratively and the poorly mapped loci were removed. The synteny analysis between the homeologous genome was performed on the basis of the mapped homeologous loci (Foncéka et al., 2009). The graphical linkage maps were drawn using SpiderMap software (Rami, unpublished) and the graphical genotype was drawn using the GGT software (van Berloo, 2008).

Among the 602 SSR markers screened, 447 (74%) detected a polymorphism between the parental genomes. A total of 431 (71.6%) SSR loci were polymorphic between the A homologous genomes (A₁ and A₂) and 465 (77.2%) between the B/K homologous genomes. From these polymorphic SSR markers, 380 that provided accurate amplification profiles were analyzed in the F₂ progeny and generated 562 loci, out of which 378 (67%) segregated as codominant and 184 (33%) segregated as dominant. Among the 378 codominant loci, 194 (51%) were assigned to the A-genomes and 184 (49%) were assigned to the B/K-genomes. Among the 184 dominant loci, 70 (38%) were from the A-genomes, whereas 114 (62%) were from the B/K-genomes.

The segregation of the SSRs was assessed with regard to their informativity for distinguishing polysomic versus disomic inheritance. This allowed for the classification of the SSR into ten classes (Table 2). The SSR markers within class 1 appeared as the most informative ones because they marked the homeologous genomes in both parents, allowing the identification of the expected genotypes for each of the meiotic behaviors in the F₂ progeny. They were chosen to illustrate the type of inheritance.

The segregation analysis of 380 SSR markers revealed that 338 (88.9%) were inherited in an F₂ population, as expected under disomic behavior. Two compelling examples of disomic segregations are shown in Figure 1 for the SSR markers Seq8D09 and AC3C02. These primer pairs marked the homeologous genomes in both parents and amplified up to four segregating bands. Each band was assigned to a subgenomic allele (“A₁”, “B”, “A₂” and “K”). The homologous bands of A-genomes (“A₁” and “A₂”) and B/K-genomes (“B” and “K”) are inherited as two independent codominant loci in the F₂ progeny, Aa at one locus and Bb at the other (Figure 1).

Interestingly, 42 (11.1%) SSR markers exhibited several unexpected F₂ band-phenotypes, which are impossible to explain under disomic behavior. Null bands for one of the homologous genomes were observed in some progeny, although the primer pair amplified the two homeologous genomes. Two examples of these unexpected F₂ phenotypes are shown in Figure 1 for the SSR markers IPAHM108 and Ah3TC39B04. These markers amplified up to 4 segregating bands. However, the homologous bands of the A-genomes (“A1” and “A2”) are lacking for some F₂ progeny, whereas the homologous bands of the B/K-genomes (“B” and “K”) are observed, excluding the missing data assumption. Those F₂ progeny are indicated with red arrows on Figure 1. Reciprocally, for the Ah3TC39B04 marker, the homologous bands of the B/K-genomes are missing, while the homologous bands of the A-genome are observed in one F₂ progeny (marked with a green arrow on Figure 1). These A₁A₂- and BK-phenotypes are unexpected under disomic behavior, but are expected under polysomic behavior (Table 1). These phenotypes occurred if and only if the homologous alleles ended up in the same gamete during the parental and/or the F₁ meiosis.

Remarkably, the value of the intensity of the “B” and “K” bands varied significantly among the BK-phenotypes for the SSR markers that segregated in a polysomic-like pattern, suggesting differences in the number of alleles. For example, in some BK-phenotypes, the more intense bands (“K”) represent two and three doses, while the less intense bands (“B”) represent a single dose (Supplementary Tables S1-S10). The K/B intensity ratio ranged from 1 to 3, indicating that the genotype of the BK F₂ band-phenotypes was either BBKK (1:1) or BKKK (1:3). The BKKK genotype is unexpected in the F₂ progeny derived from the A₁A₂BK F₁ hybrid (Table 1). The presence of BKKK genotypes suggests the occurrence of homeologous recombination during the meiosis of the parents, particularly in the synthetic allotetraploid.

Moreover, for three polysomic SSR such as RN13D04, only one allele was observed in one F₂ progeny whereas two alleles were present in both parents with up to four alleles segregating in the population (data not shown).

Among the 378 codominant loci scored in the F₂ progeny, 336 showed disomic inheritance. Of these, 279 (83%) followed the 1:2:1 Mendelian segregation ratio expected in an F₂ population and the remaining 57 (16%) loci deviated significantly from it (P < 0.05). A total of 17 and 40 loci were distorted for the A and B/K-genomes, respectively (Table 3). Of the 17 distorted loci for the A-genome, 14 were skewed toward the cultivated genotypes, whereas one and two loci were distorted in favor of the wild and heterozygous genotypes, respectively. Inversely, among the 40 loci distorted for the B/K-genome, 26 and 14 were skewed toward the wild and heterozygous genotypes, respectively. These results indicated differences between genomes and genotypes for SD.

All the 42 codominant SSR markers that showed an unexpected genotype under disomic inheritance were distorted when tested under the 1:2:1 Mendelian segregation ratio. We were not able to test for the segregation ratio expected under polysomic behavior (1:34:1) since the genomic constitution of the F₁ hybrids was unknown and unexpected genotypes could result from “aberrant” F₁ hybrids. Nonetheless, as indicated on the map below, these codominant distorted markers were clustered in particular regions on different LG, suggesting that such distortion likely originated from homeologous pairing.

In our study, the genetic map was developed progressively due to the mixture in the dataset of codominant loci that segregated in a disomic or polysomic form, combined in some cases with the distortion of segregation. The map constructed without the distorted and the polysomic loci comprised 270 loci distributed on 16 LGs, while 34 loci remained unlinked. Distorted loci were then included in the map using LOD score values comprised between 4 and 7, r = 0.3. This allowed for the addition of two new LGs (LGB2 and LGB9) that were mainly formed by distorted loci (Figure 2). At this step, the map comprised 304 loci clustered in 18 LGs. In the third step, when adding the polysomic loci, two new LGs (LGA4 and LGB4) were formed and two initially distinct LGs (LGA3 and LGB3) clustered together. The clustered LGs were dissociated by removing three polysomic loci (Seq16F01, Ah3TC31H02, and Seq2H11). At this step, the resulted genetic map comprised 20 LGs, and 12 loci remained unlinked. The loci were then ordered within each LG. Finally, the estimated linkage map included 357 codominant loci distributed into 20 LGs and covering a total genetic distance of 1,728 cM, with an average interval of 5.1 cM between two adjacent markers.

We analyzed the distribution of the loci in the A and B/K-genomes. For the A-genome, 179 loci were mapped in 10 LGs with an average number of 18 markers per LG, ranging from 10 (LG A5) to 36 (LG A1). The length of the LGs ranged from 28.0 cM (LG A5) to 165.2 cM (LG A1), with an average of 104.3 cM. For the B/K-genome, 178 loci were mapped on 10 LGs with an average number of 18 markers per LGs ranging from 13 (LG B10) to 28 (LG B1). The length of the LGs ranged from 33.3 cM (LG B10) to 111.4 cM (LG B1), with an average of 68.5 cM (Table 4). The average interval between adjacent markers was 6.4 and 4.0 cM for the A- and B/K-genomes, respectively.

Ten homeologous LGs were clearly identified based on the common homeologous loci (Figure 2). The number of bridge markers within each pair of homeologous LGs ranged from 3 to 16 (Figure 2). A good collinearity was observed among seven pairs of homeologous LGs (A2/B2; A3/B3; A4/B4; A5/B5; A7/B7; A8/B8; and A10/B10). Except for one major inversion on LG6, two local rearrangements on LG9 and three on LG1, a good synteny of the markers along the LGs was observed (Figure 2).

Except for three loci, all distorted loci (P < 0.05) were clustered along the LGs. Only 17 distorted loci were mapped in the A-genome, whereas 40 were mapped for the B genome. The distorted loci mapped in the A-genome were skewed toward the cultivated genotype (A4, A7, and A9 LGs), whereas those distorted in the B-genome were mainly skewed toward the wild genotype (B2 and B9 LGs) or the heterozygous genotype (B3, B4, B6, B7, and B10 LGs) (Figure 2). Overall, with few exceptions, we found that SD occurred in the regions that displayed homeologous recombination uncovering that homeologous pairing plays an important role in shaping SD in the tetraploid peanut genome.

One of the most remarkable features of this map is the distribution of the disomic and polysomic loci along the LGs. Among the ten homeologous sets of LGs, six sets consisted of markers that segregated in a likely disomic pattern and three in a mixed pattern with the disomic and polysomic loci clustered on the same LG (Figure 2). Surprisingly, one set (LG4) consisted of markers that all segregated in a polysomic way. The percentage of polysomic loci along the LGs ranged from 13.3% (LG6) to 100% (LG4) with an average of 39.9%.

The co-inheritance of the polysomic loci allowed for the estimation of the portion of the genome that underwent homeologous recombination. These portions ranged from 14.9 cM for LG6 to 95.5 cM for LG4, with an average of 43.7 cM (Table 3). Our study clearly showed that for some progeny, the chromosomes were a mosaic of homologous and homeologous regions (Figure 3). Taking the physical distance covered by the LGs into account, a completed substitution of the A chromosome by its homeologous counterparts was observed in three F₂ genotypes (Figure 3).

The location of the homeologous recombination breakpoints along the chromosomes was assessed for LG4, which displayed a full homeologous pairing. Of the 17 regions in which the homeologous exchanges occurred, 10 regions of the B-genome were replaced by the A corresponding region and the seven others exhibited the reciprocal situation (Figure 3). The number of the homeologous recombination breakpoints along the LGs 4A and 4B ranged from 1 to 3.

In our study, we showed that a large number of SSR markers were inherited likely in a disomic way, but some others showed genetic inheritance as a polysomic mode leading to a complex meiotic behavior.

Inheritance patterns of molecular markers have been considered a powerful method for determining meiotic behavior in polyploidy species (Lerceteau-Köhler et al., 2003) and have generated interesting conclusions about the genome behavior of several species including bermugrass (Guo et al., 2015), mimulus (Modliszewski and Willis, 2014), chrysanthemum (Klie et al., 2014), kiwifruit (Wu et al., 2013), roses (Koning-Boucoiran et al., 2012), yam (Bousalem et al., 2006; Nemorin et al., 2012), citrus (Kamiri et al., 2011), swithgrass (Okada et al., 2010), Yellow Cress (Stift et al., 2008), tomato (Barone et al., 2002), birdsfoot trefoil (Fjellstrom et al., 2001), sugar cane (Hoarau et al., 2001), and alfalfa (Diwan et al., 2000). The molecular methods used to distinguish disomy and polysomy are usually based on the signal intensity of PCR products and comparison of the number of loci linked in coupling versus the repulsion phase. However, the interpretation of multiple dose markers is often difficult in polyploids and is impossible in some species with polysomic inheritance (Esselink et al., 2004; Landergott et al., 2006).

In our study, to partially overcome the problem of complex electrophoresis profiles, we undertook a direct interpretation of SSR bands, suitably assigned to the subgenome of both parents. Using this approach, we observed that some markers followed the disomic inheritance as usually reported in allotetraploid peanut (Burow et al., 2001; Foncéka et al., 2009; Shirasawa et al., 2013) but some others exhibited an F₂ genotype unexpected under disomic inheritance, but that fitted with a polysomic segregation. The exclusive presence of homologous alleles among some F₂ genotypes for SSRs that mark homeologous genomes indicates that the homologous alleles ended up in the same gamete during the parental and/or the F₁ meiosis. Our findings are consistent with the recent study published by Leal-Bertioli et al. (2015) that explained the missing data observed among genotypes by tetrasomic recombination.

In some cases, the segregation patterns of SSR loci were similar to that observed in case of double reduction. We observed a phenotypic class with only one allele in some F₂ progeny, whereas up to 4 alleles segregate in the population. Double reduction refers to the fact that for a specific locus, the sister alleles come together in the same gamete during meiosis. It has been reported that the rate of double reduction is expected to increase towards the telomeres (Koning-Boucoiran et al., 2012; Bourke et al., 2015). In our study, the loci that showed this peculiar pattern of segregation are located at the extremity of LG3 (IPAHM103 and IPAHM282) and LG7 (RN13D04). However, since the genotype of the F₁ plant that gave rise to the F₂ progeny is unknown one cannot exclude that this pattern of segregation arose from the fusion of male and female gametes with the same haplotype (i.e., KK × KK).

The construction of the linkage map allowed us to locate the regions where the genetic exchanges occurred between the homeologous genomes. The LGs 3, 4, 6, and 7, which underwent homeologous recombination in our study, are consistent with the findings of Leal-Bertioli et al. (2015). However, the percentage of markers involved in homeologous recombination was higher in our study compared to that mentioned in Leal-Bertioli et al. (2015) (11% vs. 3%). Moreover, in LG4, we report a complete substitution of the A-chromosome by its B/K homeologous genome in three progeny. In the study of Leal-Bertioli et al. (2015), the LG4 of the induced tetraploid parent was involved in an almost complete substitution of the B-chromosome by the A-chromosome. The similarity of the results between these two studies using different synthetic allotetraploids raised questions about the factors that drive homeologous recombination.

Pairing affinity between different sets of chromosomes was reported to be influenced by structural homology (Mason et al., 2010, 2014; Mandáková, Marhold and Lysak, 2014). Bertioli et al. (2016) reported a close similarity between the A and B wild species genomes based on sequence data comparison of A. duranensis and A. ipaensis. Moreover, in the present study and in many other genetic mapping studies in peanut (Burow et al., 2001; Foncéka et al., 2009; Shirasawa et al., 2013), a good collinearity was found between the homeologous genomes. These results are in favor of a homology-driven homeologous chromosome pairing.

However, we found some loci rearrangements between homeologous LGs that displayed a mixed inheritance. This was particularly the case for LG6, in which Bertioli et al. (2016) also reported a large inversion. Thus, there are probably others factors that drive the pairing between homeologous genomes in polyploid species. Those factors can be genetic (Cifuentes et al., 2010; Gaeta and Chris Pires, 2010), similar to that exerted by the Ph-1 locus in wheat (Moore, 2014) in rye (Lukaszewski and Kopecký, 2010) and Brassica (Nicolas et al., 2009). More studies will be needed to decipher the molecular forces that drive homeologous pairing in peanut.

The results of this study suggest that, at least in an interspecific context, the meiotic pairing in tetraploid peanut fits with the intermediate inheritance where pairing of chromosome sets ranged from strict disomic inheritance as in diploids to full polysomic as in autopolyploids (Stebbins, 1947; Sybenga, 1996; Catalán et al., 2006; Xu et al., 2015). Many important plants are polyploids and some of them have been clearly identified as intermediate genetic pattern such as chrysanthemum (Klie et al., 2014), banana (Jeridi et al., 2012), and strawberry (Lerceteau-Köhler et al., 2003). These findings also raised the question of chromosome pairing during meiosis in crosses involving only cultivated peanut varieties. Until now, inter-genomic exchanges have never been reported using genetic mapping approaches in cultivated peanut. Recently, thanks to tetraploid versus diploid sequence comparisons, genomic exchanges between A and B sub-genomes have been reported (Bertioli et al., 2016). However, one still has to puzzle out whether these exchanges resulted from ancient recombination events that occured at the infancy of cultivated peanut as a consequence of the genomic shock due to the first hybridization between the A and B subgenomes or were due to something more recent that is still occurring in peanut varieties.

The reconstruction of the genotype of the homeologous recombinants suggests that, homeologous pairings have arisen in the synthetic allotetraploid parent. Indeed, the genotype of some homeologous recombinant F₂ progeny was either BKKK or BBKK. If homeologous recombination occurred only during the meiosis of the F₁ plants and not in the meiosis of the parents, the genotype of the homeologous recombinant in the F₂ progeny would only be A₁A₁A₂A₂ or BBKK (Table 1). The BKKK genotypic composition was regarded as a consequence of homeologous recombination during the meiosis of the parents particularly in the synthetic allotetraploid. In some species, homeologous pairing during meiosis is suspected to occur in the first generations following polyploid formation (Ramsey and Schemske, 2002; Soltis et al., 2010; Szadkowski et al., 2010).

The pattern of recombination between chromosomes of related species is a key point to transfer genes between the species. The knowledge about the inheritance mode is essential information not only because it sheds light on homeologous chromosome pairing behavior but can also influence the breeding strategies that are used for cultivar development. In this study, homeologous recombination events have been located in some genomic regions that are common with the ones reported in the study of Leal-Bertioli et al. (2015). One can suppose that these regions are particularly sensitive to homeologous recombination especially when synthetic allotetraploid are used. The SSRs that are located in these specific regions, particularly those that marked the homeologous genome in both parents, have proven to be very efficient in revealing the homeologous recombination events. Thus, they can be used to detect and trace back those events in the synthetic allotetraploid parental lines before crossing, as well as in their related F₁ hybrids.

Mixed inheritance may have unwanted impacts on aberrant meiotic behavior, karyotype destabilization, and fertility reduction. However, it could also speed up the accumulation of rare but favorable alleles through homeologous recombination and marker-assisted introgression. The interspecific breeding population developed in this study would be an ideal genetic material to study the genetic effects of cumulative homologous and homeologous alleles and to transfer valuable alleles from wild species into cultigens. Despite conflicting results around the A. batizocoi taxon and its potential usefulness (Leal-Bertioli et al., 2014), its utilization for peanut breeding has been successful reported (Burow et al., 2014; Sukruth et al., 2015). Although erratic fertility was observed in some lines, the advanced backcross population developed from the same cross has a high phenotypic variation for many important agronomic traits, such as plant architecture, yield related traits, drought tolerance and resistance to leaf spot (Nguepjop et al., in preparation). We believe that the genetic variation among AB-QTL lines is increased when homeologous chromosomes pair. The mosaic compositions of the genome and the homeologous chromosome substitutions may speed up novel genetic combinations, opening new horizons for peanut breeding.