Skip to main content

Identification of consistent QTL for time to maturation in Virginia-type Peanut (Arachis hypogaea L.)



Time-to-maturation (TTM) is an important trait contributing to adaptability, yield and quality in peanut (Arachis hypogaea L). Virginia market-type peanut belongs to the late-maturing A. hypogaea subspecies with considerable variation in TTM within this market type. Consequently, planting and harvesting schedule of peanut cultivars, including Virginia market-type, need to be optimized to maximize yield and grade. Little is known regarding the genetic control of TTM in peanut due to the challenge of phenotyping and limited DNA polymorphism. Here, we investigated the genetic control of TTM within the Virginia market-type peanut using a SNP-based high-density genetic map. A recombinant inbred line (RIL) population, derived from a cross between two Virginia-type cultivars ‘Hanoch’ and ‘Harari’ with contrasting TTM (12–15 days on multi-years observations), was phenotyped in the field for 2 years following a randomized complete block design. TTM was estimated by maturity index (MI). Other agronomic traits like harvest index (HI), branching habit (BH) and shelling percentage (SP) were recorded as well.


MI was highly segregated in the population, with 13.3–70.9% and 28.4–80.2% in years 2018 and 2019. The constructed genetic map included 1833 SNP markers distributed on 24 linkage groups, covering a total map distance of 1773.5 cM corresponding to 20 chromosomes on the tetraploid peanut genome with 1.6 cM mean distance between the adjacent markers. Thirty QTL were identified for all measured traits. Among the four QTL regions for MI, two consistent QTL regions (qMIA04a,b and qMIB03a,b) were identified on chromosomes A04 (118680323–125,599,371; 6.9Mbp) and B03 (2839591–4,674,238; 1.8Mbp), with LOD values of 5.33–6.45 and 5–5.35 which explained phenotypic variation of 9.9–11.9% and 9.3–9.9%, respectively. QTL for HI were found to share the same loci as MI on chromosomes B03, B05, and B06, demonstrating the possible pleiotropic effect of HI on TTM. Significant but smaller effects on MI were detected for BH, pod yield and SP.


This study identified consistent QTL regions conditioning TTM for Virginia market-type peanut. The information and materials generated here can be used to further develop molecular markers to select peanut idiotypes suitable for diverse growth environments.

Peer Review reports


Peanut (Arachis hypogaea L.) is an important grain legume and oilseed source for human nutrition. It is grown in more than 100 countries and plays a significant role in global trade. In peanut, as in other legume crops, the growing period, or the time-to-maturation (TTM), is an essential characteristic for adaptation and yield. Although TTM in peanut is influenced by environmental conditions and agricultural practices [1, 2], it has a substantial genetic component reflected by the wide range of TTM among varieties. TTM was one of the crucial traits selected during the several thousand years of domestication and diversification. Cultivated peanut is classified into two subspecies, i.e., A. hypogaea ssp. hypogaea comprising Virginia and Runner market-types and A. hypogaea ssp. fastigiata including Spanish and Valencia market-types. The two subspecies diverge in maturity level, flowering pattern, shoot determination, and plant architecture [3]. The ssp. fastigiata is characterized by early fruit maturation, sequential flowering pattern, determinate shoot formation and erect growth habit, whereas ssp. hypogaea is late in fruit maturation and exhibits an alternative flowering pattern and indeterminate spreading or bunch habit [4]. Early-maturation (90 to 120 days post-planting) is necessary for drought avoidance in areas with a short rainy season. For this reason, Spanish market-type peanuts are predominately grown in West Africa and India, where the drought stress level is high [5]. On the other hand, high-yielding but late-maturing Runner- and Virginia- market-type peanuts are widely grown in the USA and the Middle East, where irrigation is available to most farming areas.

Variation in TTM has also been found within the ssp. hypogaea. Cultivars within subsp. hypogaea are classified as early-maturing (130–140 DPP), medium-maturing (140–150 DPP), and late-maturing (150–170 DPP) ( Developing early maturing Runner- and Virginia-type cultivars with improved yield and excellent agronomic characteristics has been an important objective of peanut breeding programs. Early maturity is essential in areas that suffer from limited water supply or by end-of-season cooler temperatures and early frosts [6] that might retard maturation, cause incomplete seed filling, lower the yield and grade or quality (including oleic to linoleic acid ratios) [7, 8].

Despite its evident impact, very little is known regarding the genetic control of TTM in peanuts. The maturity level of peanut was reported as a quantitative trait with low heritability [9, 10] and influenced by many genes and environmental factors [11, 12]. A few attempts to define QTL with a small effect for early-maturation were reported using breeding materials from fastigiata X hypogaea crosses [13,14,15]. The utility of low-density genetic maps in these studies was one of the limiting factors for QTL discovery. Domesticated peanut is a self-pollinated allotetraploid (AABB, 2n = 4x = 40), originating from two diploid progenitors, A. duranensis (AA, 2n = 2x = 20) and A. ipaensis (BB, 2n = 2x = 20) [16, 17]. Cultivated peanut has a narrow genetic base caused by the bottleneck of a single hybridization event that gave rise to this species and the crossing barrier between cultivated peanut and wild diploids due to ploidy differences [18]. With the advancement of SNP array technology [19, 20], the limitation of genetic map density was alleviated by a drastic increase in genetic markers for map construction. Close to 1000 SNP markers were placed on peanut linkage maps recently [21, 22]. In addition to the low polymorphism, another challenge is to perform phenotyping on TTM due to the unique underground formation of fruit, the indeterminate nature of pod formation, and the application of the commonly used but laborious and somewhat subjective hull-scrape method to determine pod maturity [23].

Identifying genetic mechanisms controlling TTM in peanuts has a significant practical and scientific impact. Translating major QTL controlling the trait to user-friendly marker platforms will enable marker-assisted selection (MAS) to accelerate breeding for early- or late-maturity. In addition, underlying genetic mechanisms can be further investigated based on QTL mapping discoveries.

In the present study, a SNP-based linkage map was constructed for a RIL population derived from Virginia market-type parents. QTL mapping with two-year field phenotyping data led to the discovery of two consistent QTL for TTM, for the first time reported for Virginia-type peanut.


Phenotyping of the parents and the RIL population

A RIL population was developed from a cross between Hanoch (late-maturing) and Harari (early-maturing) cultivars (Fig. 1a). The maturity index (MI), which is determined by the percentage of pods with black and brown mesocarp, was documented to indicate TTM. Data were collected from field experiments in two different environments (i.e., year, location, soil). MI values of parental lines were collected (Fig. 1b). A highly significant difference was found between the parental lines in MI (P = < 0.0001), with 30.9 ± 6.82 and 53.97 ± 7.63 for Hanoch and Harari, respectively. In addition, significant differences were found between the parental lines for all of the other measured traits on the 2 years mean data, including pod yield (PY), 50-pod weight (50PW), 50-seed weight (50SW), and shelling percentage (SP), except for harvest index (HI) (Fig. 1b).

Fig. 1

Phenotypic characterization of different traits among ‘Hanoch’ and ‘Harari’. a, MI morphology of ‘Hanoch’ and ‘Harari’. b, Comparisons between ‘Hanoch’ and ‘Harari′ in MI, PY, HI, 50PW, 50SW and SP. Data are shown as mean from two years (n = 9). The Student’s t-test was used to generate the P values. MI, maturity index (%); PY (g), pod yield; HI (g.g-1), harvest index; 50PW (g), 50 pod weight; 50SW (g), 50 seed weight; SP (%), shelling percentage

Normal or close to normal distribution was found in the RIL population data for all measured traits (Fig. 2; Table 1).

Fig. 2

Phenotypic distribution of MI and the other traits in two consecutive years. Y-axis corresponds to the RIL population density, and X-axis corresponds to the original metric trait value based on the average of three replicates. Arrows indicate the phenotypic values for Hanoch (yellow) and Harari (green). A normal distribution curve is indicated in red. MI, maturity index; PY, pod yield; HI, harvest index; 50PW, 50 pod weight; 50SW, 50 seed weight; SP, shelling percentage. _18, year 2018; _19, year 2019

Table 1 Summary statistics of MI and other traits among parents and RILs

Parental values of MI were within the range of the RILs. Some RILs exhibited a MI value beyond parental values at each end of the curve in both years, suggesting transgressive segregation of MI in this population. A significant effect was found for the blocks, RIL, year, and RIL X year interaction (Table 2) from ANOVA analysis. Therefore, QTL analysis was performed with data from each year separately. The broad-sense heritability for MI was 0.39, indicating a moderate-to-low but significant genetic component underlying this trait. The heritability estimates for other traits ranged from 0.07 (PY) to 0.36 (HI) (Table 2).

Table 2 Analysis of variance and heritability for MI and the other traits for the Hanoch X Harari RIL population across two years. Block [Year] indicates the nested effect of the Blocks within each year

Pearson correlation among the traits was calculated in each year (Fig. 3). The correlation of the two-year MI measurements was 0.59 (p < 0.0001), suggesting a relatively higher genetic heritability of this trait than estimated by the ANOVA. MI was significantly correlated with HI and PY in both years and was correlated with SP in 2018. A small but significant correlation was found between MI in 2019 and 50PW in 2018. MI showed no correlation with 50SW. Significant correlations were observed among the other traits, such as SP with HI, SP with PY, and 50PW with 50SW. The Branching Habit (BH) phenotype effect (spreading vs. bunch) on MI was inspected by a T-test (Additional file 1: Fig. S1). A significant but small effect was found for the branching habit (BH) phenotype on MI in both years, in which spreading type lines had higher MI values than bunch types.

Fig. 3

Pearson correlations for MI and the other traits evaluated in Hanoch X Harari RIL over two years. MI, maturity index; PY, pod yield; HI, harvest index; 50PW, 50 pod weight; 50SW, 50 seed weight; SP, shelling percentage. _18; year 2018; _19; year 2019. Significance of correlations: * P < 0.05 and ** P < 0.001

Construction of the genetic map

Genotyping of Hanoch x Harari RILs was performed with version 2 of the Axiom Arachis_SNP array consisting of 47 K SNP markers (Thermofisher Scientific). A set of 3283 polymorphic SNP markers between the two parental lines [24] was used in this study. After filtering and removing the missing data and heterozygous calls, 3074 SNPs were retained for the RIL population. Twenty-five RILs with greater than 10% missing data and greater than 20% heterozygous SNP calls were removed from further analysis. Subsequently, a genetic map was constructed with 235 RILs. Also, 773 SNPs that did not obey the chi-square test in the JoinMap tool, 457 SNPs that were identical to other loci, and 11 SNPs that generated final small non-significant linkage groups were excluded.Therefore, the genetic map contained 1833 markers distributed on 24 linkage groups covering a total of 1773.5 cM (Fig. 4; Table 3) (Additional file 2: Table S1).

Fig. 4

An overview of the genetic map and the QTLs identified for MI and the other traits. MI, maturity index; PY, pod yield; HI, harvest index; 50PW, 50 pod weight; 50SW, 50 seed weight; SP, shelling percentage; BH, branching habit

Table 3 Description of genetic linkage groups. Physical distance was determined by blasting markers against the A. hypogaea reference genome (

The 24 linkage groups ranged in size from 9.4 cM (B10_2) to 244.6 cM (A06). The average number of loci per linkage group ranged from 76, reaching up to 385 loci in LG A06. The average distance between the neighboring markers was 1.6 cM, ranging from 0.6 cM in LGs A07, A06, B06, and B03 to 4.9 cM in A09_2 (Table 3). Aligning the 1833 markers to the A. hypogaea pseudomolecules ( resulted in a total physical distance of 1952.6 Mbp and an average physical interval of 2.8 Mbp between loci (Additional file 1: Fig. S2; Table 3). The percentage of a pseudomolecule covered by linkage groups varied; eight groups covered more than 80% of a pseudomolecule, seven more than 90%, and one group (A06) was close to 100%. The average recombination rate was 0.7 cM/Mbp. A08 had the maximum recombination rate, while the groups B10_2, A09_1, B03, and B05_2 had the lowest recombination rates.

The linkage map quality was assessed by analyzing the SNPs’ collinearity to their physical positions (Mbp) in the A. hypogaea genome (Additional file 1: Fig. S2). As expected, the saturation of the markers in the arms was higher than in the pericentromeric regions. Some rearrangements were exhibited in a few linkage groups, such as apparent inversions in the middle of A07 and at the end of B06 (Additional file 1: Fig. S2).

QTL identification

QTL mapping of the MI and the other traits resulted in identifying 30 QTLs, with the LOD scores ranging from 3.03 to 81.2, explaining 5.8 to 79.6% of the phenotypic variance (PVE) (Fig. 4; Table 4). Nine linkage groups had at least one QTL, with a maximum of eight QTLs in B06 and 7 QTLs in A06. HI had the maximum number of QTLs, 4 QTLs each in 2018 and 2019. Major QTLs were found for BH and HI, explaining 79.6 and 20.6% PVE, respectively.

Table 4 QTL identified for MI and the other traits in the Hanoch X Harari RIL population

For MI trait, a total of six QTLs were identified, three QTLs each in 2018 and 2019, respectively (Fig. 4), explaining 9.3 to 11.9% PVE. Two consistent QTL regions were found in both years. One was observed on LG A04 between AX-176819644_A04 - AX-176815499_A04, spanning 6.9 Mbp, with PVE values of 9.9 and 11.9% for 2018 and 2019, respectively. The other consistent QTL region was observed on LG B03 within marker interval of AX-176807311_B03 - AX-176806413_B03, spanning 1.8 Mbp, explaining 9.3 and 9.9% PVE. The other two QTLs were identified on LG B06 (qMIB06) and LG B05_2 (qMIB05_2), which were significant only in 2018 and 2019, respectively. Alleles from the early-maturing Harari parent contribute to the high percentage of mature pods measured by MI for four QTLs, qMIA04a, qMIA04b, qMIB03a and qMIB03b (Table 4). The late-maturing Hanoch parent contributed to qMIB06 and qMIB05_2.

For PY, three QTLs were detected in 2018 (Fig. 4; Table 4), two on LG B06 (qPYB06.1 and qPYB06.2), explaining 9.2 and 6.8% PVE, respectively. The other, qPYB05_2, was identified on group B05_2, explaining 6.7% PVE. Two of these QTLs, qPYB06.1 and qPYB05_2 overlapped with qMIB06 and qMIB05_2 of MI, respectively.

A total of eight QTLs were detected for HI (Fig. 4), explaining 5.8–20.6% PVE. QHIB05_2 in 2018 shared a common region with qMIB05_2 in 2019 and qPYB05_2 in 2018. Another common QTL region shared between qHIB03 in 2019 and MI trait QTLs, qMIB03a in 2018 and qMIB03b in 2019 (Table 4). These overlapping QTLs detected between HI and MI traits aligned with the significant correlation between the traits.

For 50PW and 50SW, eight QTLs were identified (Fig. 4). Q50PWA06 in 2018, q50SWA06.1 and q50SWA06.2 in 2018 identified on A06 showed 11.8, 6.4 and 6.2% PVE, respectively. Similarly, q50PWA04 and q50SWA04 in 2019 were identified on A04, explaining 9 and 13.9% PVE, respectively. Significant QTLs, q50PWA09_1 and q50SWA09_1 were identified in 2018 on LG A09_1 spanning 1.7 Mbp (AX-176821658_A09 - AX-177644544_A09) explaining 9.2 and 10.5% of variation respectively (Table 4). As expected, there was a strong co-localization of QTLs between 50PW and 50SW. An overlapping QTL region was observed in 2019 between q50SWA04 and MI trait QTLs (qMIA04a and qMIA04b) spanning 1 Mbp around marker interval AX-176802283_A04 - AX-176819644_A04.

For SP, four QTLs were identified (Fig. 4), three in 2018 and one in 2019. QSPA06.1 and qSPA06.2 were identified on A06 in 2018, explaining 9.3 and 9.1% PVE, respectively. Additionally, qSPB02 in 2018 and qSPA09_1 in 2019 were observed on B02 and A09_1, explaining 14.1 and 5.8% PVE, respectively (Table 4). An overlapping QTL region was found among q50PWA09_1 and q50SWA09_1 in 2018 and qSPA09_1 in 2019.

For the BH trait, one very significant and consistent QTL, qBHB05_2 with marker interval AX-147251167_B05 - AX-147251374_B05 on LG B05_2 spanning 3.1 Mbp explaining 79.6% PVE was found. Since almost no differences were found within the same RILs in BH phenotype between 2018 and 2019, the same locus was denoted for both years (qBHB05_2). QBHB05_2 overlapped with the MI QTL, qMIB05_2 on LG B05_2 in 2019, suggesting a possible BH effect on MI (Fig. 4; Table 4).


Time-to-maturation (TTM) is one of the crucial traits for adaptability and yield in legumes. Late-maturation usually is associated with increased yield and prolonged pod-filling processes [25, 26]. In contrast, early-maturation is associated with better adaptation to terminal stresses and to avoid lodging in some legumes [27, 28]. Genetic factors and their interactions with the environment play crucial roles in the control of TTM. Two main developmental aspects control TTM in legumes, flowering time and plant architecture. Flowering time models in legumes include the vernalization-responsive long-day model and the warm-season short-day model [29]. Inflorescence architecture, the second factor that regulates TTM in legumes, derives from the final identity of the shoot apical meristem [30]. Most legume plants have evolved to a “complex” indeterminate architecture with compound inflorescences [30]. Still, some determinate varieties have been selected in several legumes, such as beans [31], soybeans [32], and peas [33], displaying a shorter flowering time and compact canopy to facilitate mechanized harvesting [34].

Peanut exhibits a unique genetic TTM system for legumes. It is basal to the phaseoloid clade, and therefore should be considered a short-day plant. However, studies clearly show that the time to first flower is minimally affected by photoperiod in peanut [35]. Peanut genotypes typically initiate flowering at about ~ 30 days post sowing regardless of the growing season. Inflorescence architecture systems are more relevant to peanut since the two main cultivated peanut subspecies, fastigiata and hypogaea, differ in both TTM and flowering patterns. Yet, the lateral shoots of both subspecies have indeterminate growing tips [36], contradicting the “classic” inflorescence architecture system as a possible contributor to differences in TTM between cultivars.

In the current study, we used a unique system to analyze the TTM genetics in peanuts. The segregating RIL population is created from two closely related parental lines that flower at the same time (~ 30 DPP) and have very strong indeterminate growth habits. Additionally, the lines do not differ in the flowering patterns (both have alternate flowering) and branch length, which were previously shown to associate with pod maturity in fastigiata peanuts [13]. Other traits that may influence TTM, such as pod number/plant and flowering rate [9, 37, 38], are also irrelevant for this specific population. Therefore, this system is interesting for revealing new genetic components that control TTM in the Virginia-type peanut. Indeed, due to the low genetic variation between Hanoch and Harari, the constructed genetic map seems to have some gaps and low coverage in few chromosomes. However, the two parental lines were part of the 20 genotypes used to construct the Axiom Arachis_SNP array, increasing the chance to exploit the best potential of this low polymorphic background.

A significant difference was found between the parental lines in MI in multiple environments. This result agrees with previous observations from commercial fields, showing Harari 12–15 days earlier than Hanoch. The range of MI data in the RIL population extended beyond both parents’ means, suggesting transgressive segregation of MI in this population. Some inconsistency was found between the years in the differences between the parental lines in other traits such as 50SW and PY, partly explained by the relatively small plot size in the experiments. The broad-sense heritability estimates for MI were moderate (~ 0.4) but somewhat higher than reported in other studies [9, 10], demonstrating the relatively strong genetic component and adequate phenotyping in the current system.

The most notable finding in this study was the relatively high phenotypic correlation between MI and harvest index (HI), which was significant in both years. Also, qMIB03a and qMIB03b - qHIB03, qMIB05_2 – qHIB05_2, and qMIB06 – qHIB06a and qHIB06b were shared between the traits, indicating a pleiotropic effect. Indeed, the ratio between pod yield biomass and the entire plant biomass at the later stages of development may influence TTM by increasing the plants’ sink capacity and promoting the crop termination. Since the total pod yield (PY) effect was not as strong as HI, we speculate that the canopy biomass had a stronger contribution to TTM than PY. This phenomenon is documented in several legume crops [39,40,41]. We suggest that sink strength is the most potent effector for early maturation in Virginia-type peanuts, particularly in closely related germplasm, as were used in this study. Interestingly, HI and PY are both quite similar between the parents, although they have distinct MI. In contrast, a relatively large variation in HI and PY was found between the RILs. This suggests that, in addition to MI, other factors determine HI and PY, leading to significant transgressive segregation of these traits in the population.

Another trait that can have some effect on MI in our system is branching habit. Branching habit (BH) is an important descriptive and agronomic character of peanut. In a previous study [42], we showed that the BH in the Hanoch X Harari genetic system is controlled by a single gene that was named Bunch1. Using a previous version on the Arachis SNP-array, Bunch1was located to a ~ 1.1 Mbp segment on the same locus as was found here on B05. So, it is not surprising to find such a high %PVE for BH (Table 4). In the current study, we found that BH has a small but significant effect on MI, particularly in 2019 (Additional file 1: Fig. S1). This is also reflected by the fact that MI and BH shared the same QTL in 2019. The spreading form of BH was associated with higher MI. Interestingly, the allele from the late maturing parental line (Hanoch; spreading) contributed to early maturity at this QTL region. We speculate that pods are formed closer to the ground in the spreading form, and therefore, have the opportunity to mature earlier.

The majority of the QTLs identified have small to moderate effects on MI (Table 4), two of them were consistently detected in both years. Together they explained ~ 20% of the total phenotypic variation for MI and ~ 50% of the genetic variation (taking into account that the heritability estimate was ~ 0.4), indicating that they were not spurious. These two genomic locations are different from QTLs found in other studies [13,14,15], involving hypogaea X fastigiata crosses. Thorough literature screening showed that none of the previously described QTL for the other traits in this study (HI, PY, 50SW, 50PW and SP) matched the QTL found here, indicating that they are unique to the Virginia-type background as well.

An indication for the possible influence of these two QTLs in peanut maturation came from a retrospective analysis of the genotypic situation in newly bred Israeli cultivars, ‘Orit’ and ‘Einat’, both originating from similar crosses of Hanoch X Harari, and are presumed to be early-maturing [43]. However, recent multi-location observations of commercial plots have indicated that only Einat is an early maturing variety, while Orit is a more medium-maturing variety. Indeed, inspecting these two varieties’ genotypic information showed that only Einat contains the two consistent QTLs (data not presented). This example illustrates the potential of implementing MAS to provide efficient and unbiased selection of traits with high phenotyping costs such as early- or late-maturity in peanuts.

In conclusion, this study demonstrates the use of SNP-array technology for constructing and applying a genetic map in a biparental population with very low polymorphism. A new genetic map with 1833 SNP markers was constructed for the Virginia peanut background. A total of six QTLs regulating TTM were identified across 2 years of field tests. The novel information and materials generated here can promote the selection of peanut idiotypes using genetic markers associated with the QTLs discovered in this study.

Materials and methods

Plant material and growing conditions

A recombinant inbred line (RIL) population was developed from a cross between cv. Hanoch and cv. Harari [44], two closely related Israeli Virginia-type cultivars differing in TTM (Fig. 1a). RILs were obtained by a single seed descend procedure, up to F7, and then were multiplied as bulks for additional two generations (F7:9). ‘Hanoch’ has been the leading Israeli in-shell peanut cultivar for over two decades. It is a late-maturing cultivar with long, smooth, and hard pods. These qualities make it well received in the EU “in-shell” market. ‘Harari’ is an early-maturing cultivar grown in Israel’s northern part. Its growing season is limited by late sowing time (due to the double-cropping system) and early harvest (due to autumn rainfalls). Harari has reticulated and soft pod walls and is targeted for the local shelled industry. Harari has a bunch-type growth habit while Hanoch is spreading. Common traits between Hanoch and Harray include pod size, flowering time, flowering pattern, lateral branch length, pod number/plant and flowering rate. Both parental lines were part of the Arachis_SNP-array development panel [45].

A total of 260 RILs were planted in two successive years, respectively. The first year was planted in April 2018 in the Hula Valley, Northern Israel (33°11′17.7″N 35°34′25.6″E), characterized by heavy black soil. The second was planted in April 2019 in Urim, Western Negev, Israel (31°20′27.4″N 34°29′46.1″E), characterized by fine sandy-loam. Besides the soil type, these two regions are significantly different in environmental conditions. Hulla Vally is a typical semi-arid climate with rainy winters and high summer humidity, while Nirim is located in the Negev Desert with low humidity. A similar experimental design of randomized complete blocks with three replications was implemented in both seasons. Each ‘plot’ (line X block) consisted of two rows on a bed, 4 m in length, rows spaced 90 cm apart, and seeding rates of 10 seeds/m2 (total of 20 plants/plot). Parental lines were grown as control plots with nine replications. Fields were maintained under full-irrigation conditions, and all recommended agronomic practices were carried out as previously described [44]. All plant material, including the parental and the RILs, was originated from Hovav laboratory and there are part of the ARO breeding program.

Phenotyping the maturity trait and post-harvest traits

TTM was evaluated at ~ 140–145 days post-planting (DPP). The exact sampling date was determined by testing the parental lines every few days, starting at 125 DPP, up to the point where Harari was ~ 60% mature on average. This timing was chosen in order to capture the widest variation in maturation among the RILs. The hull-scrape method [23] was used to measure the maturity level by randomly sampling 2–3 plants per plot and removing the exocarp from all pods using a PICO water pressure machine (Idromatic®, Italy), 14 MPa with 9 l/minute flow rate. Pressure-washed pods were separated into five categories based on mesocarp color: white, yellow, orange, brown, or black. Pod number in each category was documented, and the maturity index (MI) was calculated as the percentage of pods in the brown and black categories. In total, 729 and 780 MI measurements were taken for 2018 and 2019, respectively.

Other traits potentially associated with TTM were also recorded. Branching habit (BH) was documented at ~ 50 DAP as a spreading or a bunch. Negligible variations were found between 2018 and 2019 for BH data. Therefore, one value was used for BH in both years. After digging, plots were dried for 7–10 days, and the entire biomass of each plot (excluding the tap-roots) was weighed. Each plot was then threshed by an experimental thresher (Kincade, USA), and total pod yield/plot (PY) was measured. Harvest index (HI) was calculated as the ratio of pod yield/plot over biomass weight/plot. 50-pod weight (50PW), 50-seed weight (50SW), and shelling percentage (SP) were recorded by randomly sampling 100 pods from each plot.

Statistical analysis of phenotypic data

Statistical difference between the parents was determined by Student’s t-test. As for the RILs, Anderson-Darling test was performed to determine the normality of distribution. The ANOVA analysis model included the RIL, Year, Year X RIL and Block [Year] effects. All effects were defined as random to calculate the heritability rates. Broad sense heritability (H2) was estimated with the equation H2 = σg2/(σg2 + σge2 + σe2), by the ANOVA analysis with QTL IciMapping v4.2 ( [46]. The σg2, σe2 and σge2 denoted the variances of genotypes (G), environment (E) and interaction of genotypes and environments (G x E). Correlation coefficients were calculated among all the traits across years. One way ANOVA was performed to check the effect of the BH phenotype on MI. Distribution and correlation statistics, histograms, and boxplots were performed with JMP® Pro 15 (SAS Institute Inc., Cary, NC, 1989–2019).

Genotyping and genetic map construction

Genomic DNA was extracted using DNeasy® Plant Mini Kit (Qiagen; Hilden, Germany) from young leaflets from each RIL and the two parents. DNA quantification was performed with Qubit (Invitrogen; CA, USA). The samples were diluted to 40 ng/μL according to protocol guidelines and genotyped by using the Affymetrix Axiom_Arachis2 SNP array comprising 47,837 SNPs, divided into their AA and BB subgenomic origin [20, 47]. Genotyping data were analyzed by the Axiom analysis suite Software 3.1 [24]. The polymorphic homozygous SNPs (AA and BB) and polymorphic heterozygous SNPs (AA or BB and AB) were retained with 65–35% call-rate frequencies among the RILs. Out of 260 RILs, 25 RILs with greater than 10% missing data and greater than 20% heterozygous SNP calls were removed from further analysis. Subsequently, the genetic map was constructed with 235 RILs. The genetic linkage map was constructed using Joinmap v4.1 [48] maximum likelihood (ML) algorithm with a minimum LOD of 3.0 and the Haldane mapping function. The graphical representation of the linkage maps was generated through Mapchart v2.3 [49]. Confirmation of the loci positions was done as previously described [47] with few modifications (BLASTN (e value < 1 × 10− 18) and mismatch of less than 2). Linkage groups (LG) generated were assigned to the pseudo-molecules of the tetraploid A. hypogaea cv. Tifrunner [17] ( Due to the high sequence similarity between the two subgenomes of cultivated peanut [16], the position assignment of the 70-bp-long SNP markers from the array can be ambiguous. Linkage groups were assigned to the respective pseudomolecules (chromosomes) of the sequenced A. hypogaea genome cv. Tifrunner (Table 3). If a LG had more than 51% of the SNPs representing a particular chromosome, then this LG was assigned to that chromosome. To assess the quality of the genetic map, a collinearity analysis was performed using the genetic distances (cM) versus the physical positions (Mbp).

QTL analysis

QTL mapping was performed on 235 RILs using MapQTL v6 [50] on the mean phenotypic data collected each year (Additional file 2: Table S2). A LOD score of 3 with 1000 permutations was used to confirm the presence of a putative QTL at a 95% significance level. QTLs were manually assigned to the genetic maps. The naming of the QTL follows the terminology of “q” as QTL, followed by an abbreviation of the trait. The last digit represents the LG, and repetition of the QTL in alphabetical order if in both years, or numerical order if there are more than one, on the same LG. SNP markers flanking the QTL were used to obtain the physical position from the A. hypogaea genome.

Availability of data and materials

The data sets supporting the results of this study are included in the manuscript and additional supporting files.



50-pod weight


50-seed weight


Additive effect


Harvest index


Linkage group


Logarithm of the odds


Maturity index


Pod yield


Phenotypic variation explained


Quantitative trait locus


Recombinant inbred line


Shelling percentage


Single nucleotide polymorphism


Time to maturation


  1. 1.

    Canavar O, Kaynak MA. Growing degree day and sunshine radiation effects on peanut pod yield and growth. Afr J Biotech. 2010;9:2234–41.

    Google Scholar 

  2. 2.

    Inanaga S, Yoshida T, Hoshino T, Nishihara T. The effect of mineral elements on the maturity of peanut seed. Plant Soil. 1988;106(2):263–8.

    CAS  Article  Google Scholar 

  3. 3.

    Nigam SN, Aruna R. In: Janick J, editor. Improving breeding efficiency for early maturity in peanut. In: Jules Janick, editor. Plant Breeding Reviews, Vol 30. Wiley; 2008. p. 295–322. ISBN 978040171523.

  4. 4.

    Stalker HT, Simpson CE. Genetic resources in Arachis, pp. 14-53. In: Pattee HE, Stalker HT, editors. Advances in Peanut Science. Stillwater: American Peanut Research and Education Society; 1995.

  5. 5.

    Stalker HT, Tallury SP, Seijo GR, Leal-Bertioli SC. Chapter 2 - Biology, Speciation, and Utilization of Peanut Species. In: Stalker HT, F. Wilson R, editors. Peanuts: AOCS Press; 2016. p. 27-66.

  6. 6.

    Bell MJ, Roy RC, Tollenaar M, Michaels TE. Importance of variation in chilling tolerance for peanut genotypic adaptation to cool, short-season environments. Crop Sci. 1994;34(4):1030–9.

    Article  Google Scholar 

  7. 7.

    Sanders TH. Maturity distribution in commercially sized Florunner peanuts. Peanut Sci. 1989;16(2):91–5.

    Article  Google Scholar 

  8. 8.

    Upadhyaya HD, Reddy LJ, Gowda CLL, Singh S. Identification of diverse groundnut germplasm: sources of early maturity in a core collection. Field Crop Res. 2006;97(2-3):261–71.

    Article  Google Scholar 

  9. 9.

    Ali N, Wynne JC. Heritability estimates and correlation studies of early maturity and other agronomic traits in two crosses of peanuts (Arachis hypogaea L.). Pak J Bot. 1994;26:75–82.

    Google Scholar 

  10. 10.

    Chiow HY, Wynne JC. Heritabilities and genetic correlations for yield and quality traits of advanced generations in a cross of peanut. Peanut Sci. 1983;10(1):13–7.

    Article  Google Scholar 

  11. 11.

    Pattee HE, Singleton JA, Sanders TH. Composition changes of peanut fruit parts during maturation. Peanut Sci. 1974;1(2):57–62.

    CAS  Article  Google Scholar 

  12. 12.

    Rowland DL, Sorensen RB, Butts CL, Faircloth WH. Determination of maturity and degree day indices and their success in predicting peanut maturity. Peanut Sci. 2006;33(2):125–36.[125:DOMADD]2.0.CO;2.

    Article  Google Scholar 

  13. 13.

    Chirinos FV, Isleib TG, Milla-Lewis SR. Breeding for early maturity in peanuts (Arachis hypogaea L.) using traditional methods and marker assisted selection (MAS). 2011;PhD thesis, North Carolina State University, Raleigh, North Carolina.

  14. 14.

    Fonceka D, Tossim HA, Rivallan R, Vignes H, Faye I, Ndoye O, et al. Fostered and left behind alleles in peanut: interspecific QTL mapping reveals footprints of domestication and useful natural variation for breeding. BMC Plant Biol. 2012;12(1):26.

    Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Selvaraj MG, Narayana M, Schubert AM, Ayers JL, Baring MR, Burow MD. Identification of QTLs for pod and kernel traits in cultivated peanut by bulked segregant analysis. Electron J Biotechnol. 2009;12.

  16. 16.

    Bertioli DJ, Cannon SB, Froenicke L, Huang GD, Farmer AD, Cannon EKS, et al. The genome sequences of Arachis duranensis and Arachis ipaensis, the diploid ancestors of cultivated peanut. Nat Genet. 2016;48(4):438–44.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Bertioli DJ, Jenkins J, Clevenger J, Dudchenko O, Gao DY, Seijo G, et al. The genome sequence of segmental allotetraploid peanut Arachis hypogaea. Nature Genet. 2019;51(5):877–84.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Zhuang WJ, Chen H, Yang M, Wang JP, Pandey MK, Zhang C, et al. The genome of cultivated peanut provides insight into legume karyotypes, polyploid evolution and crop domestication. Nat Genet. 2019;51(5):865–76.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Clevenger J, Chu Y, Chavarro C, Agarwal G, Bertioli DJ, Leal-Bertioli SCM, et al. Genome-wide SNP genotyping resolves signatures of selection and Tetrasomic recombination in Peanut. Mol Plant. 2017;10(2):309–22.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Clevenger JP, Korani W, Ozias-Akins P, Jackson S. Haplotype-based genotyping in polyploids. Front Plant Sci. 2018;9:564.

    Article  PubMed  PubMed Central  Google Scholar 

  21. 21.

    Chu Y, Chee P, Culbreath A, Isleib TG, Holbrook CC, Ozias-Akins P. Major QTLs for resistance to early and late leaf spot diseases are identified on chromosomes 3 and 5 in peanut (Arachis hypogaea). Front Plant Sci. 2019;10:883.

    Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Chavarro C, Chu Y, Holbrook C, Isleib T, Bertioli D, Hovav R, et al. Pod and seed trait QTL identification to assist breeding for peanut market preferences. G3 Genes Genomes Genet. 2020;10(7):2297–315.

    CAS  Google Scholar 

  23. 23.

    Williams EJ, Drexler JS. A non-destructive method for determining peanut pod maturity. Peanut Sci. 1981;8(2):134–41.

    Article  Google Scholar 

  24. 24.

    Patil A, Popovsky S, Levy Y, Chu Y, Clevenger J, Ozias-Akins P, et al. Genetic insight and mapping of the pod constriction trait in Virginia-type peanut. BMC Genet. 2018;19(1):93.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Clements JC, Dracup M, Galwey N. Effect of genotype and environment on proportion of seed hull and pod wall in lupin. Aust J Agric Res. 2002;53(10):1147–54.

    Article  Google Scholar 

  26. 26.

    Gupta K, Hedvat I, Faigenboim-Doron A, Clevenger J, Ozias-Akins P, Hovav H. Transcriptome profiling of peanut developing seed with a focus on duplicate oil related pathways. Plant Sci. 2016;248:116–27.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Siddique KHM, Loss SP, Thomson BD. Cool season grain legumes in dryland Mediterranean environments of Western Australia: significance of early flowering. In: Saxena N, editor. Management of Agricultural Drought. Enfield, NH: Enfield Publishers; 2003. p. 151–61.

    Google Scholar 

  28. 28.

    Duc G, Agrama H, Bao SY, Berger J, Bourion V, De Ron AM, et al. Breeding annual grain legumes for sustainable agriculture: new methods to approach complex traits and target new cultivar ideotypes. Crit Rev Plant Sci 2015;34:381–411, 1-3, DOI:

  29. 29.

    Weller JL, Ortega R. Genetic control of flowering time in legumes. Front Plant Sci. 2015;6.

  30. 30.

    Jonsell B, Weberling F. Morphology of flowers and inflorescences. Nord J Bot. 1989;1991.

  31. 31.

    Kelly JD. Remaking bean plant architecture for efficient production. Adav Agron. 2001;71:109–43.

    Article  Google Scholar 

  32. 32.

    Boote KJ, Jones JW, Batchelor WD, Nafziger ED, Myers O. Genetic coefficients in the CROPGRO-soybean model. Agron J. 2003;95(1):32–51.

    Article  Google Scholar 

  33. 33.

    Foucher F, Morin J, Courtiade J, Cadioux S, Ellis N, Banfield MJ, et al. DETERMINATE and LATE FLOWERING are two TERMINAL FLOWER1/CENTRORADIALIS homologs that control two distinct phases of flowering initiation and development in pea. Plant Cell. 2003;15(11):2742–54.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Tian ZX, Wang XB, Lee R, Li YH, Specht JE, Nelson RL, et al. Artificial selection for determinate growth habit in soybean. P Natl Acad Sci USA. 2010;107(19):8563–8.

    Article  Google Scholar 

  35. 35.

    Bagnall DJ, King RW. Response of peanut (Arachis hypogaea) to temperature, photoperiod and irradiance .2. Effect on peg and pod development. Field Crop Res. 1991;26(3-4):279–93.

    Article  Google Scholar 

  36. 36.

    Reddy PR. Physiology. In: Reddy PR, editor. Groundnut. New Delhi: Indian Council of Agricultural Research; 1988. p. 77–119.

    Google Scholar 

  37. 37.

    Cahaner A, Ashri A. Vegetative and reproductive development of Virginia-type peanut varieties in different stand densities. Crop Sci. 1974;14(3):412–6.

    Article  Google Scholar 

  38. 38.

    Puangbut D, Jogloy S, Vorasoot N, Kesmala T, Holbrook CC, Patanothai A. Response of reproductive parts of peanut genotypic variation and their contributions to yield after pre-flowering drought. Aust J Crop Sci. 2013;7:1627–33.

    Google Scholar 

  39. 39.

    Angus IF, Hasegawa S, Hsiao TC, Liboon SP, Zandstra HG. The water balance of post-monsoonal dryland crops. J Agri Sci. 1983;101(3):699–710.

    Article  Google Scholar 

  40. 40.

    Ayaz S. Variability of harvest index in four grain legume species. Canterbury, New Zealand: Unpublished PhD thesis, Lincoln University; 2001.

    Google Scholar 

  41. 41.

    Jasil P, Akkaseang C, Kesmala T, Jogloy S. Heritability and correlation for maturity and pod yield in peanut. J App Sc Res. 2011:134–40.

  42. 42.

    Kayam G, Brand Y, Doron-faigenboim A, Patil A, Hedvat I, Hovav R. Fine-mapping the branching habit trait in cultivated peanut by combining bulked Segregant analysis and high-throughput sequencing. Front Plant Sci. 2017;8:467.

    Article  Google Scholar 

  43. 43.

    Hovav R, Hedvat I, Levy Y (2017) The development and registration of two early-maturing Virginia-type peanuts in Israel (Hebrew). Israel Groundnuts Production & Marketing Board annual scientific reports Pp 12–16.

  44. 44.

    Patil A, Hedvat I, Levy Y, Galili S, Hovav R. Genotype-by-environment effects on the performance of recombinant inbred lines of Virginia-type peanut. Euphytica. 2018;214(5):83.

    CAS  Article  Google Scholar 

  45. 45.

    Clevenger J, Chu Y, Chavarro C, Agarwal G, Bertioli DJ, Leal-Bertioli SCM, et al. Genome-wide SNP genotyping resolves signatures of selection and tetrasomic recombination in peanut. Mol Plant. 2016;10:309–22.

    Article  Google Scholar 

  46. 46.

    Zhang S, Hu X, Miao H, Chu Y, Cui F, Yang W, et al. QTL identification for seed weight and size based on a high-density SLAF-seq genetic map in peanut (Arachis hypogaea L.). BMC Plant Biol. 2019;19:537.

    CAS  Article  Google Scholar 

  47. 47.

    Van Ooijen JW. JoinMap® 4 software for the calculation of genetic linkage maps in experimental populations. JoinMap 2006.

  48. 48.

    Voorrips RE, Chen X, Liang X, Liu H, Zhou G, et al. MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered. 2002;93(1):77–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  49. 49.

    Korani W, Clevenger J, Chu Y, Ozias-Akins P. Machine learning as an effective method for identifying true SNPs in polyploid plants. Plant Genome. 2019;12(1):180023.

    CAS  Article  Google Scholar 

  50. 50.

    Ooijen JW Van. MapQTL 6. Genome. 2009.

Download references


The authors express heartfelt gratitude to Dr. Chandrasekhar for his support in the field experiment.


This study was funded by the Israeli Ministry of Agriculture (grant no. 0142-01-18) and an Israel-USA Bilateral Agriculture Research and Development (BARD) grant (grant no. IS-5020-17). The funders did not have any scientific part in the study.

Author information




SK is the Ph.D. student that is in charge of the study and drafted the manuscript. SA and YL helped with the molecular labor. IC is the field experiment manager. YC performed SNP-array analysis and revised the manuscript. RH and POA are the PI and the CO-PI of this project, respectively. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Ran Hovav.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Fig. S1.

Box plot analysis to study the effect of branching habit (X-axis) on maturity index (Y-axis) across 2018 and 2019. P values were generated through Student’s t-test. The color of the boxes indicates as follows, blue, bunch; red, spreading. Fig. S2. Correlation between the genetic distance (cM) (x-axis) of markers on each linkage group (LG) and the physical genome position (Mbp) (y-axis) based on the Tifrunner reference genome. Black dots represent markers mapped to the respective chromosome, red dots indicate markers mapped to the homeologous chromosome and black circles represent markers mapped to other chromosomes.

Additional file 2: Table S1.

Population genotype data presented by the order of linkage map. Physical positions of markers were based on the tetraploid peanut genome sequence [17]. Table S2. Mean phenotypic values of measured traits over two years of field tests.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Kunta, S., Agmon, S., Chedvat, I. et al. Identification of consistent QTL for time to maturation in Virginia-type Peanut (Arachis hypogaea L.). BMC Plant Biol 21, 186 (2021).

Download citation