Map construction
All markers evaluated in this study were amplified by the use of heterologous primers. With the use of minimum LOD score of 6.0 and a maximum recombination fraction (θ) of 0.35, and after the exclusion of co-segregating markers, we found that 399 markers mapped onto 10 LGs, spanning a total map distance of 678.2 cM. These markers included 378 microsatellites and 21 transposon (miniature inverted-repeat transposable element) markers. LGs were numbered and oriented essentially according to Moretzsohn et al. (2009), but with B05 and B08 reoriented (“flipped”) to ensure forward compatibility with pseudomolecules produced by the Peanut Genome Sequencing Consortium (www.peanutbase.org; Genbank GCA_000816755.1). LGs ranged from 41.5 cM (with 35 markers) to 139.2 cM (55 markers), with an average distance of 1.7 cM between adjacent markers. A total of 91 (22.8%) of the 399 mapped markers deviated from the expected 1:1 ratio at P < 0.05 level. Of these, 79 markers were skewed toward A. ipaënsis and only 12 markers toward A. magna. Most LGs have few distorted markers, and LG B05 and B08 had no distorted markers. In contrast, LG B02 and the upper portion of LG B04 were almost entirely composed of distorted markers. These two LGs and LG B06 were distorted toward A. ipaënsis alleles. Distorted markers at P < 0.05 were identified by # (Figure 3).
Figure 3  A genetic linkage map of the B-genome of Arachis obtained through the analysis of 94 F6 plants, generated from a cross between A. ipaënsis K 30076 and A. magna K 30097. Numbers on the left of each group are Kosambi map distances (cM). Markers that amplified more than one loci have numbers _1 and _2 after the marker name. Quantitative trait loci (QTL) are indicated as colored bars running alongside linkage groups. Colors/textures are according to categories: black, rust resistance [total number of lesions/leaf area (TL/LA), number of sporulated lesions/leaf area (SL/LA), susceptibility index (SI), and Incubation Period (IncPer)]; white, seed characteristics [seed number (SN) and 10-seed weight (10_SW)]; textured, plant architecture [main stem height (MSH)]; and domestication traits [peg length (PL) and pod constriction (PC)]. Distorted alleles (P < 0.05) are indicated by #.

QTL identification:
The framework map, containing 399 markers, was used for QTL analysis. LOD significance threshold estimated for each trait ranged from 2.5 to 3.2, and only QTL with LOD values exceeding these values were included. At least one QTL was detected for each of the 15 traits analyzed, with a total of 38 QTL mapped. A summary of QTL is provided in Table 1. Detailed information about markers and QTL are presented in File S1.

QTL for rust resistance:
For rust resistance, 13 QTLs were identified on the two bioassays. A major QTL for the four components of rust resistance (SI, TL/LA, SL/LA, and IncPer) was consistently identified in both bioassays and mapped at the same marker interval on map position 35.1−42.9 cM on LG B08, with LOD scores between 2.9 and 8.2. Its closest marker is Ah-280. This QTL explained 5.8–59.3% of the phenotypic variance of the four different traits (SI, TL/LA, SL/LA, and IncPer). Another QTL for SI and TL/LA evaluated in 2012 and for IncPer_2013 was found at 25.4−33.1 cM on the same LG B08 (Table 1, Figure 3), explaining 13.2–34.8% of the total variance in 2012 and 2013, respectively. The closest markers are AHGS1350 and AHGS2541. In addition to these two QTL, two minor QTL were identified in LGs B04 and B07 for SI_2012. For all QTL, resistance was derived from A. magna. The contribution of A. magna alleles for the rust resistance-related traits was evaluated by calculating the average phenotypes of the homozygote plants for the A. magna allele vs. the average phenotype of those homozygote for the A. ipaënsis allele (Figure 4). Arachis magna alleles contributed significantly for the reduction of SI, TL/LA, and SL/LA. The effect of A. magna alleles was less pronounced for IncPer (Figure 4).
Figure 4  Bar graph of contribution of A. ipaënsis K 30076 (blue bars) vs. A. magna K 30097 (orange bars) alleles to the rust resistance-related traits: susceptibility index (SI), number of sporulated lesions/leaf area (SL/LA), total number of lesions/leaf area (TL/LA). and incubation period (days/20) (IncPer/20).

QTL for other agronomic/domestication traits:
Agronomic and domestication traits were evaluated in three different years. QTL were consistent between years; therefore, for final QTL analyses, data from different years were pooled and considered different replicates of the same experiment. The exception was MSH, which was not pooled, because measures were performed in different stages of plant growth in each year. In total, 25 QTL were identified. Interestingly, the QTL for two domestication traits (PC and PL) were placed in the same map positions: 37.1−40.8 cM in LG B01 and 64.3−64.9 cM in LG B04 (Table 1, Figure 3). LG B04 also harbors a cluster of QTL for MSH at 64.3−64.9 cM, explaining 10.7–30.4% of the total phenotypic variance. For all these traits, the main markers associated were AHGS1917 and AHGS2155. Alleles increasing values of domestication-related traits originated from both parents. However, alleles increasing seed number are all derived from A. ipaënsis. No QTL for domestication/productivity traits co-localized with the rust resistance QTL.