A new source of root-knot nematode resistance from Arachis stenosperma incorporated into allotetraploid peanut (Arachis hypogaea)

Root-knot nematode is a very destructive pathogen, to which most peanut cultivars are highly susceptible. Strong resistance is present in the wild diploid peanut relatives. Previously, QTLs controlling nematode resistance were identified on chromosomes A02, A04 and A09 of Arachis stenosperma. Here, to study the inheritance of these resistance alleles within the genetic background of tetraploid peanut, an F2 population was developed from a cross between peanut and an induced allotetraploid that incorporated A. stenosperma, [Arachis batizocoi x A. stenosperma]4×. This population was genotyped using a SNP array and phenotyped for nematode resistance. QTL analysis allowed us to verify the major-effect QTL on chromosome A02 and a secondary QTL on A09, each contributing to a percentage reduction in nematode multiplication up to 98.2%. These were validated in selected F2:3 lines. The genome location of the large-effect QTL on A02 is rich in genes encoding TIR-NBS-LRR protein domains that are involved in plant defenses. We conclude that the strong resistance to RKN, derived from the diploid A. stenosperma, is transferrable and expressed in tetraploid peanut. Currently it is being used in breeding programs for introgressing a new source of nematode resistance and to widen the genetic basis of agronomically adapted peanut lines.

RKN resistance into peanut dates from 2001, with the release of the cultivar COAN 16 . The resistance was based on a chromosome A09 segment from the wild species Arachis cardenasii Krapov. & W.C. Greg. 17 . The segment was introgressed through a backcrossing scheme involving interspecific hybrids 18,19 . Since then, several other cultivars have been released in the U.S. using the same source of resistance: NemaTAM, Tifguard, Webb, Georgia-14N and TifNV-High O/L [20][21][22][23][24] .
Additional sources of resistance are important for the development of new high-yielding and nematode-resistant peanut cultivars, and to reduce the risk of resistance breakdown in the varieties currently used 8 . Previously, we studied the peanut wild relative Arachis stenosperma Krapov. & W. C. Greg., which harbors resistance to a number of pests, including the RKN Meloidogyne spp. and foliar diseases, such as late leaf spot and rust [25][26][27][28] . Subsequently, using diploid mapping populations, we identified a large-effect QTL controlling RKN resistance on chromosome A02 together with minor effect QTLs on A04 and A09 of A. stenosperma (accession V10309) 29 . Strong resistance to RKN also has been reported in the wild diploid species A. batizocoi Krapov. & W. C. Greg. (accession K9484) 30 .
In this study, an F 2 population derived from the cross of A. hypogaea with an induced allotetraploid (A. batizocoi x A. stenosperma) 4× 31 was used to identify genome regions that confer RKN resistance within a tetraploid genetic background; to develop reliable molecular markers tightly linked to the resistance loci for selection in breeding programs; and finally, to characterize the genetic behavior of wild-cultivated crosses. This research will contribute to the production of advanced peanut lines that incorporate wild-derived chromosome segments that confer a new source of resistance to RKN.

Results
nematode screening. Resistance to RKN was evaluated over three years by measuring three traits: EGR (eggs/gram of root), RF (reproductive factor), and GI (galling index). Arachis stenosperma V10309, A. batizocoi K9484 and the induced allotetraploid BatSten1 ([A. batizocoi x A. stenosperma] (2n=4×=40) were found to be resistant to RKN, with no or small galls and low egg production. Cultivar Runner-886 was susceptible in all the assays and presented the highest values for all traits (Fig. 1, Table S1 and Tab 'Traits for F 2 lines' in File S1). Individuals from the RBS-F 2 population (derived from a cross between A. hypogaea cv. Runner IAC-886 and the synthetic allotetraploid BatSten1) showed varying degrees of resistance to RKN. The frequency distributions of EGR, RF and GI traits were distinctly non-normal, being skewed towards resistance (Fig. 1). EGR, RF and GI values for the RBS-F 2 were significantly different at P < 0.05 for 2014 and 2016 assays, but not for 2015. Values for the two parents (Runner-886 and BatSten1) were significantly different (P < 0.05) for 2014 and 2016 assays (in 2015 data for BatSten1 was not available). Transgressive segregation was also observed: across years, on average 49.3%, 48.0% and 59.0% segregating individuals were as, or more, resistant than A. stenosperma for EGR, RF and GI, respectively. For the three measured traits, few lines were more susceptible than Runner-886; on average 0.7%, 0.3% and 0.3% individuals were as, or more, susceptible than Runner-886 for EGR, RF and GI, respectively. Several individuals that produced galls did not support the production of nematode eggs. These traits were highly correlated: Pearson correlation was calculated in Minitab v.15.1.0.0 between EGR, RF and GI (Tab 'Pearson Correlation' in File S1) and significant values (P ≤ 0.01) ranged between 0.281 and 0.887.
Genotyping and genetic mapping. Using the Affymetrix ' Axiom_Arachis' SNP array 32,33 196 F 2 progeny and controls were genotyped. Two individuals were eliminated since they had too many missing data points. For the rest of the analyses, data of 194 progenies were used. A total of 1587 polymorphic SNP markers were identified from our filtering strategy, 911 assigned to A-subgenome (A. stenosperma-specific markers) and 676 to B/K-subgenome (A. batizocoi-specific markers); additionally, 9696 A. hypogaea-specific markers were identified. Only markers derived from A. stenosperma and A. batizocoi were utilized for genetic mapping and QTL identification. The physical positions of the A-genome markers were determined according to the position of their homologues in the A. duranensis pseudomolecules and the B-genome markers based on the A. ipaensis pseudomolecules (www.peanutbase.org) 2 . After removing low-quality and unlinked markers, 1499 SNP markers were ordered into 20 linkage groups (LGs) that ranged in size from 100.7 cM (LG A05) to 359.5 cM (LG A02), spanning a total map distance of 3984.9. The average distance between adjacent markers ranged from 1.35 cM (LG A01) to 3.8 cM (LG A06) and the largest distance was 29.48 cM (LG A06). The number of markers in each linkage group varied from 36 (LG B04) to 217 (LG A01) (Table S2, Figure S1 and Tab 'Framework Map' in File S1). There was a strong relationship between genetic and physical positions of their homologues. As expected, higher recombination frequencies were observed in the distal parts of the chromosomes. LG A07 presents an inversion on the lower end relative to its diploid ancestor, but consistent with the genome sequence of A. hypogaea ( Figure S2) 1 . The heatmap of the marker-pairwise estimated recombination fractions versus LOD scores indicated that there is only linkage within LGs, with the exception of LGs A05 and B05, where markers seem to be tightly associated with each other, especially at the beginning and end of the linkage groups ( Figure S3). This is consistent with tetrasomic genome composition and recombination between alleles between these homeologous chromosomes 1,2 .

QTL identification.
A large-effect QTL for all three measurements of resistance to RKN was detected at the bottom of LG A02 with LOD scores above 5.9 and with a peak at marker A02_89159922_Sten (320.7 cM) at genome-wise α = 0.05 and 0.01 thresholds ( Table 1, Fig. 2, Figure S4A-C and Tab 'QTL' in File S1). Using EGR and RF data, the QTL was identified at 1% level of significance ( Figure S4A,B and Tab 'QTL' in File S1) and for GI data at 5% ( Figure S4C and Tab 'QTL' in File S1). A second QTL for EGR and RF was observed in LG A04 below the α = 0.05 genome-wise threshold. This QTL seems to comprise two separate loci with LOD scores above 3.2. The highest peak was near A04_111013470_Sten SNP marker (160.6 cM), while the lower peak was near A04_110684871_Sten (139.6 cM) ( Table 1, Fig. 2, Figure S5A,B and Tab 'QTL' in File S1). A third QTL, for RF, was present in LG A09 with LOD score of 3.7 at marker A09_116335836_Sten (139.6) ( Table 1, Fig. 2, Figure S5C and Tab 'QTL' in File S1).
All QTLs derived from A. stenosperma are colocalized with the ones found in a previous study on a A. stenosperma-derived diploid population ( Fig. 2 and Tab 'Colocalization' in File S1) 29 . The QTL on LG A02 contributed to a percentage reduction up to 98.2% in nematode multiplication. Although putative QTL on A09 did not pass the estimated threshold, it reduced the nematode development up to 97.8%. Putative QTL on A04 did not contribute to improve nematode resistance (Table 1). Although A. batizocoi was also found to be resistant to RKN 30 , A. batizocoi-derived QTLs associated with nematode resistance were not identified.
The analysis of phenotypic effects of markers tightly linked to QTLs contributing to nematode resistance reveals that most F 2 lines with A. stenosperma alleles at QTLs on LG A02 and A09 showed the lowest phenotypic mean scores for all the measurements of resistance (EGR, RF and GI); whereas, F 2 plants carrying A. hypogaea alleles had a higher incidence of the disease (Fig. 3A,C), Graphs for RF and GI are not shown, but have the same trend. This was also supported with the positive additive effect values ( Table 1), suggesting that these genome segments contribute to reduced egg production, nematode reproduction and gall formation, as previously described 29 . This tendency was also observed when analyzing combined effect of two QTLs. F 2 lines that were homozygous for A. stenosperma alleles on A02 and A09 together had low disease scores (Fig. 3D), on average presenting reduction of EGR, RP and GI, of 87.3%, 95.2% and 67.5%, respectively. Conversely, the putative QTL on LG A04 does not produce an improved resistance when the F 2 lines were homozygous or heterozygous for the A. stenosperma alleles (Fig. 3B,E,F), which is supported with the negative additive effect we found for this segment (Table 1). It is possible that this trend on the A04 putative QTL reveals underdominance, where heterozygous F 2 lines have an inferior performance for nematode infection than either the cultivated or wild homozygous genotypes (Fig. 3B). Alternatively, the lower resistance could have arisen from unfavorable gene interactions 34 . Validation of resistance. Progeny of F 2 lines produced on average 10 seeds per plant (median = 4). Six lines that produced between 12 and 36 seeds were selected to be planted in the field for yield and disease resistance evaluation. Two lines that carried markers associated to QTLs in LG02 and LG09 were selected and advanced to the F 2:3 generation. To confirm the RKN resistance on the F 2:3 lines, cuttings where phenotyped for RKN resistance using the same screening approach described in Methods but using eight replicates per genotype and inoculation with 5.1K second-stage juveniles (J2) per pot (Tab 'Traits for F 2 -derived F 3 lines' in File S1). According to Kruskal-Wallis test (P < 0.05), there was a significant difference in RKN resistance between genotypes. Further analysis included the Wilcoxon signed rank test for pairwise comparisons using FDR (false discovery rate) correction, in order to group the samples by significant similarity (P > 0.05). According to the grouping, for EGR, RF and GI all the F 2 -derived F 3 lines were significantly different from the TifGp-2 susceptible genotype ( Table 2). These lines were also genotyped using the ' Axiom_Arachis v01' 58 K high-density SNP array 32,33 and data processed similarly to the F 2 genotypic data. Only genotyping calls for the A. stenosperma and A. batizocoi-derived markers were taken into account, the data was then filtered as Batsten1 ≠ A. hypogaea (AppendixS2) and then visually inspected for the A. stenosperma introgressions along the three chromosome segments associated with nematode resistance in LG A02, A04 and A09. F 2:3 -7 and F 2:3 -34 were found to carry the resistance segments in heterozygous or homozygous states for A. stenosperma alleles ( Table 2). inheritance patterns and segregation distortion. In an F 2 population with disomic recombination, three main genotypic classes are expected: the two parental and the hybrid types (Fig. 4A). A set of 1156 SNPs (77.1%) were inherited as expected under disomic inheritance: 694 (46.3%) markers followed the Mendelian segregation ratio expected for an F 2 population (1:2:1) and 462 (30.8%) showed significant deviation (P < 0.01) ( Table 3, Tab 'Tetrasomic Recombination' in File S1). From the 234 (15.6%) distorted loci located in the A-genome, distorted blocks were skewed toward the cultivated genotype specially on LGs A04 and A05. Among the 228 (15.2%) loci in B/K genome, we found a biased segregation in favor of wild in B06 and an excess of cultivated alleles in B07 ( Figure S6 and Tab 'Tetrasomic Recombination' in File S1).
During marker analyses, it was also noted that several markers exhibited unexpected genotypic classes, in other words additional clusters to the three genotypic classes expected with disomic recombination. To investigate the nature of these markers, the clustering of 1499 polymorphic SNPs for the 194 F 2 lines were visually inspected and confirmed. The percentage of markers in each LG that showed tetrasomic recombination in at least one genotype, ranged from 4.3% (LG B03) to 75% (LG A04) with an average of 22.9% (343 SNPs) (Table S2). These informative SNPs showed unexpected clustering patterns (Fig. 4B-D), which was explained by recombination between homeologous chromosomes (Table 3, Tab 'Tetrasomic Recombination' in File S1). A total of 565 individual data points across the whole data set (1499 SNPs x 194 individuals) (0.19%) were informative to identify homeologous recombination. Lines that are triplex (single allele replaced by its homeolog) or quadriplex (both alleles replaced by homeolog) are expected only when individuals have undergone polysomic recombination; Fig. 4C,D shows examples of assays that can distinguish two products of homeologous recombination: triplex (green) and quadriplex (gray). In the example in Fig. 4B. It is unclear whether or not the individual that has undergone homeologous recombination for this locus is triplex or quadriplex, since there is no information about the allele dosage in the B-subgenome. Segregation ratio and distortion was not possible to describe for the markers with tetrasomic behavior, since the information about the F 1 hybrid and its gametic allelic constitution is unknown. Additionally, we observed blocks of genome substitution especially on LG A02/B02 (Fig. 5), A03, A04, A08 and A10/B10. A04 was the LG with highest number of markers exhibiting homeologous recombination (Tab 'Tetrasomic Recombination' in File S1). All of above is consistent with previous reports for F 2

Discussion
The peanut wild relative A. stenosperma has been shown to be resistant to multiple pests and pathogens including M. arenaria Race 1. This resistance is manifest at multiple stages of the infection cycle: both much lower rates of penetration of nematodes, and hypersensitive response to the few nematodes that do penetrate are observed 27 .

Trait Symbol
LG a www.nature.com/scientificreports www.nature.com/scientificreports/ Genes involved in hypersensitive response and secondary metabolite production for defense against nematode infection have been identified in differential analyses of gene expression and histology 27,28,37 . Prior to the present work, RKN resistance segments on A02, A04 and A09 of the wild A. stenosperma were discovered using the simplified genetic context of a diploid mapping population derived from a cross of A. duranensis x A. stenosperma 29 . Here, we used an F 2 population from the cross of the cultivar Runner-886 and the wild-derived allotetraploid Batsten1 31 to analyze the genome regions conferring resistance in a tetraploid background. The aim of this work was to provide a framework of knowledge to incorporate this new source of resistance into elite peanut cultivars through marker-assisted backcrossing schemes (MABC) 38 , work that is now in progress. Using the high-density genetic map and RKN phenotypic data measured on the F 2 progeny, allowed us to validate the previously described diploid QTLs, in the tetraploid context of cultivated peanut. We identified three A. stenosperma QTLs associated with nematode resistance: on the bottom of LG 02, middle of LG 04 and bottom of LG 09.
The QTL on LG A02 was consistently associated with lower EGR, RF and GI (measures of resistance) ( Table 1, Fig. 2, Figure S4, and Tab 'QTL' in File S1). The QTL located on LG A02 was found in the same region as the previously described diploid QTL (Fig. 2 and Tab 'Colocalization' in File S1). The diploid QTL peak was mapped at 87.4Mbp with an interval from 83.6Mbp to 92.5Mbp 2,29 . Here in the tetraploid context we located this resistance segment at 89.2Mbp with an interval between 84.3Mbp and 92.1Mbp. The second tetraploid QTL on LG A04 was near the diploid QTL at 111.0 Mbp. The QTL we identified in A09 was located at 116.3 Mbp, a little further from the previously reported diploid position (112.8 Mbp) (Fig. 2 and Tab 'Colocalization' in File S1) 2 . With this information, we could demonstrate that RKN resistance from A. stenosperma is transferable and stable in tetraploid genotypes.
The chromosome segment conferring resistance on LG A02 is homologous to the region of the reference genome of A. duranensis, that harbors multiple resistance-genes (R-genes) 2 . A gene encoding a toll/interleukin-1 receptor (TIR)-like-nucleotide binding (NBS)-leucine-rich repeat (LRR), associated with plant immune defenses 39-41 is close to the A02-QTL peak (89.2Mbp) (Fig. 2 and Tab 'Colocalization' in File S1). Recently, at www.nature.com/scientificreports www.nature.com/scientificreports/ a different genome location, a different TIR-NBS-LRR gene was observed to be constitutively expressed in the cultivar Tifguard, which carries RKN resistance from the wild species A. cardenasii and absent in the susceptible cultivar Gregory 42 .
The effects of the A02 and A09-QTL on nematode resistance was confirmed in selected F 3 lines (F 2:3 -7 and F 2:3 -34) carrying resistance loci in a heterozygous ( Table 2, "−/+") or homozygous ( Table 2, "+/+") state. The presence of A. stenosperma alleles at the QTL positions significantly reduced egg production (EGR), inhibited the nematode reproduction (RF) and decreased gall formation in comparison with susceptible genotype TifGP-2. We can infer that the presence of resistance segments is important to halt the completion of different steps of nematode life cycle. In the future we intend to carry out histological work on advanced backcross lines carrying different combinations of wild species chromosome segments, to provide insights into the timing and nature of the resistances conferred by these chromosome segments as in Proite et al. 27 . Currently lines F 2:3 -7 and F 2:3 -34 are being crossed and backcrossed with several agronomically elite peanut lines.
Although the main focus of this study was to discover, introgress and validate QTLs associated with nematode resistance, interesting non-disomic inheritance of markers was also detected in the F 2 genotyping data. In a tetraploid hybrid context with cultivated peanut, the A. stenosperma genome is expected to recombine predominantly with A. hypogaea A-subgenome and the K genome (B genome sensu lato) of A. batizocoi with the A. hypogaea B-subgenome. Previously, genetic mapping and QTL identification studies assumed disomic inheritance in Arachis tetraploid species [43][44][45][46][47] . However, Leal-Bertioli et al. 36 provided the first molecular evidence of non-homologous alleles recombination in peanut, when unexpected genotyping patterns where detected for some loci in cultivated x artificially induced allotetraploid RIL lines as well as the parent. Later, this phenomenon was reported in an F 2 progeny, also derived from a peanut by synthetic tetraploid cross 35 . More recently tetrasomic recombination has been reported in pure A. hypogaea crosses 1,2,32 .
In this study, using Affymetrix genotyping, most SNP markers were inherited as expected for disomic segregation, but others appeared to have undergone homeologous recombination. The percentage of markers showing www.nature.com/scientificreports www.nature.com/scientificreports/ tetrasomic recombination found in this study (22.9%) was higher than previously found by other groups 35 for an F 2 progeny (11.05%).
LG A04 had the highest number of markers showing tetrasomic behavior (Table S2); in one F 2 plant we observed almost a complete substitution of the B/K04 and B/K10 alleles by their homeologs on A04 and A10 (Tab 'Tetrasomic Recombination' in File S1). This recombination between homeologous chromosomes supports the idea that cultivated peanut has segmental genetic inheritance, where behavior is mostly disomic but partially polysomic. Since tetrasomic recombination occurs mostly in gene-rich regions (distal parts of chromosomes), understanding this phenomenon is not only important to avoid disregarding these genomic segments during genetic mapping 32 , but also to address peanut breeding strategies to transfer genes between species and to accelerate the accumulation of favorable alleles through marker-assisted introgression 35 .
From the analysis of segregation distortion (Chi-square test, P < 0.01), we observed that markers with significant deviation from expected Mendelian segregation ratio (1:2:1) were distributed as clusters, possibly located within segregation distortion regions (SDRs) 48 . The SDRs have been shown to be present in interspecific or wide crosses in plants [48][49][50][51] and also described in A. hypogaea intra and interspecific mapping populations 35,44,46,52 . Here, distorted markers on LG A04, A05 and B07 showed a bias toward A. hypogaea, indicating strong selection against the A. stenosperma alleles. In chromosome B06, we found a biased segregation in favor of wild A. batizocoi (FigureS 6 and Tab 'Tetrasomic Recombination' in File S1). In general, more distorted loci were found in favor of the A. hypogaea parent.
LG Marker  www.nature.com/scientificreports www.nature.com/scientificreports/ Materials and Methods plant materials. Wild Arachis accessions were obtained from the USDA-GRIN system (https://www. ars-grin.gov/). To introgress the diploid nematode resistance into tetraploid peanut, a synthetic allotetraploid (BatSten1) was created using the peanut wild relatives A. stenosperma PI666100 (original collection voucher V10309) and A. batizocoi PI298639 (original collection voucher K9484), as described in Leal-Bertioli et al. 36 . This wild-derived allotetraploid combines the A genome of A. stenosperma 6,53 and a K genome (B genome sensu lato) of A. batizocoi 6,54 . An F 2 segregating population was created by selfing an F 1 derived from a cross between A. hypogaea cv. Runner IAC-886 (herein called Runner-886) and BatSten1. This population was named RBS-F 2 and had initially, 196 individuals. To be able to perform the bioassays in different years, the RBS-F 2 population was maintained in the greenhouse by vegetative propagation. The majority of the individuals, but not all, were maintained for the duration of this work.
Root-knot nematode resistance evaluation. The RBS-F 2 population was evaluated for resistance to root-knot nematode (RKN) M. arenaria Chitwood race 1 under greenhouse conditions 29,55 in a randomized complete block design with five replicates per genotype. The tetraploid parents (BatSten1 and Runner-886) and the diploid wild species A. batizocoi and A. stenosperma were used as controls. RKN populations were maintained and extracted from eggplant (Solanum melongena cv. Black Beauty). Arachis seeds were planted in nursery pots   www.nature.com/scientificreports www.nature.com/scientificreports/ (15 cm in diameter 10 cm in height) filled with Promix (Premier Horticulture, Quakertown, PA) and maintained in a greenhouse. Two-month-old cuttings from each F 2 line were established in steam-sterilized sandy soil and inoculated with 10,000 second-stage RKN juveniles (J2) by distributing the inoculum in two 2-cm deep holes at the base of the plant. Eight weeks later, plants were uprooted, rinsed, and the roots weighed after removing excess water with a paper towel. Roots were stained with 0.05% phloxin B solution for 3 to 5 min. Nematode eggs were extracted from roots using 0.5% NaOCl 56,57 . Assays were conducted over three years, with 155, 105 and 99 segregating F 2 s, respectively, as a few individuals died with time. Resistance was assessed using three different traits: 1) Numbers of eggs per gram of root (EGR); 2) Nematode reproductive factor (RF = Pf/Pi; where Pf is the final egg population and Pi the initial J2 population 58 ; and 3) Root-galling index or egg masses: 0 (no galling or no egg masses), 1 (1-2 galls), 2 (3-10 galls), 3 (11-30 galls), 4 (31-100 galls) and 5 (more than 100 galls or egg masses per root system) 57 .

Statistical analysis. The Shapiro-Wilk test was used to test normality of phenotypic data. Non-parametric
Kruskal-Wallis one-way analysis of variance 59 was used to assess the global differences in phenotypic traits at a 5% level of significance (P < 0.05) among RBS-F 2 lines and controls for each year using the Statistical package R. Non-normal phenotype data were transformed to Log 10 and Log 10 (x + 1) for QTL identification. Additionally, the Best Linear Unbiased Predictors (BLUPs) of random effects were calculated for each trait using the ranef function in R. BLUPs were calculated in order to control for missing phenotypic data and transformed to Log 10 and Log 10 (x + 1).  www.nature.com/scientificreports www.nature.com/scientificreports/ with PicoGreen and samples were submitted for genotyping with the ' Axiom_Arachis v01' 58 K high-density SNP array 32,33 . The genotypic data were extracted, processed and analyzed using the Axiom Analysis Suite 2.0 software (http://www.affymetrix.com). Output was analyzed using Unix scripts (AppendixS1) and data were visualized as a color map in Microsoft Excel (Tab 'Framework Map' in File S1). The strategy to identify polymorphic SNP markers included three different steps:

SNP genotyping, analysis
Firstly, informative SNP assay results were extracted from SNP calling using a panel of diploid species plus a single tetraploid genotype (A. hypogaea Runner-886). This set of markers was filtered to reveal SNP markers specific to each of the three parental species in the pedigree of the F 2 population as follows: Secondly, SNP assay results were extracted from SNP calling of tetraploid genotypes only (Runner-886, BatSten1 and RBS-F 2 population). Finally, the three sets of informative SNP markers identified in the first step were retrieved using the panel of tetraploid genotypes, followed by merging and filtering as Batsten1 ≠ A. hypogaea.
Genetic Mapping and QtL discovery. Genetic maps for A and B subgenomes were constructed using Kosambi's genetic map function 60 and maximum likelihood algorithm in JoinMap v.4.0 61,62 . The goodness of fit Chi-square test was performed to evaluate the expected 1:2:1 segregation ratio for the F 2 population for each locus (P < 0.01). The genetic map was visualized by calculating pairwise logarithm of the odds (LOD) scores and recombination fractions using the plot.rf function in R/QTL 63,64 .
The genetic map in combination with transformed RKN measurements of resistance for three years were used for QTL identification using R/QTL software following the procedure described in "A guide to QTL mapping with R/qtl" 63 . Due to the spike observed in the phenotype data distribution at zero (null phenotype, Fig. 1), we employed a two-part binary plus normal analysis method using the scanone function [63][64][65] . The two-part model is suitable for data with non-normal distribution prior and after transformations, as we observed in our case (Fig. 1). This model performs two different analyses. First, the phenotype is treated as a binary trait (0 or >0), and then as a quantitative trait, for those individuals with phenotypic values above zero 65,66 . The two-part model calculates the LOD scores for each tested genome position to assess the significant association with the trait of interest [63][64][65] . 1000 permutations were used to identify genome-wide LOD significance thresholds for QTL identification at 1% and 5% level of significance [67][68][69] . 95% Bayesian credible interval was calculated with the bayesint function and LOD support interval with lodint function in R/qtl. The percentage of phenotypic variability explained by a QTL (R 2 ) and the estimated effect was assessed using the fitqtl function in R/qtl 63 . Physical positions for each marker on the A and B subgenomes were determined, respectively, based on the A. duranensis and A. ipaensis pseudomolecules (www.peanutbase.org) 2 .

Meiotic behavior analysis and tetrasomic inheritance.
Given that the genome of tetraploid peanut harbors regions where tetrasomic recombination can occur 1,2,32,35,36 , we examined our data for evidence of such chromosomal behavior. When assuming recombination only between homologous chromosomes, the expected segregation ratio in the F 2 population should be 1(A 1 A 1 ):2(A 1 A 2 ):1(A 2 A 2 ) for A-subgenome and 1(BB):2(BK):1(KK) for the B/K-subgenome, but when homeologous chromosomes recombine during the parental and/or F 1 meiosis, the expected segregation ratio in the F 2 progeny changes, as described by Nguepjop et al. 35 .
In order to determine lines exhibiting tetrasomic recombination and the markers associated with these regions, all 1499 polymorphic markers were visually and manually inspected. The criterion was to analyze data points with "No Call" data or forming unexpected genetic clusters, similar to the rationale described in Leal-Bertioli et al. 36 but applied to an F 2 population. Genotypes and markers were scored as "tetrasomic" on the color map (Tab 'Tetrasomic Recombination' in File S1, red on A-subgenome and blue in B-subgenome). Additionally, segregation distortion from 1:2:1 ratio (P < 0.01) was analyzed for markers segregating in a disomic manner.