Mol Breeding (2017) 37: 17 DOI 10.1007/s11032-016-0587-3

Genetic mapping of yield traits using RIL population derived from Fuchuan Dahuasheng and ICG6375 of peanut (Arachis hypogaea L.) Yuning Chen & Xiaoping Ren & Yanli Zheng & Xiaojing Zhou & Li Huang & Liying Yan & Yongqing Jiao & Weigang Chen & Shunmou Huang & Liyun Wan & Yong Lei & Boshou Liao & Dongxin Huai & Wenhui Wei & Huifang Jiang

Received: 1 March 2016 / Accepted: 1 November 2016 / Published online: 30 January 2017 # The Author(s) 2017. This article is published with open access at Springerlink.com

Abstract The genetic architecture determinants of yield traits in peanut (Arachis hypogaea L.) are poorly understood. In the present study, an effort was made to map quantitative trait loci (QTLs) for yield traits using recombinant inbred lines (RIL). A genetic linkage map was constructed containing 609 loci, covering a total of

1557.48 cM with an average distance of 2.56 cM between adjacent markers. The present map exhibited good collinearity with the physical map of diploid species of Arachis. Ninety-two repeatable QTLs were identified for 11 traits including height of main stem, total branching number, and nine pod- and seed-related traits.

Electronic supplementary material The online version of this article (doi:10.1007/s11032-016-0587-3) contains supplementary material, which is available to authorized users. Y. Chen : X. Ren : Y. Zheng : X. Zhou : L. Huang : L. Yan : Y. Jiao : W. Chen : S. Huang : L. Wan : Y. Lei : B. Liao : D. Huai : W. Wei : H. Jiang (*) Oil Crop Research Institute, Chinese Academy of Agricultural Sciences/Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agricultural, Wuhan 430062, People’s Republic of Chinae-mail: [email protected]

Y. Jiao e-mail: [email protected]

Y. Chen e-mail: [email protected]

S. Huang e-mail: [email protected]

X. Ren e-mail: [email protected]

L. Wan e-mail: [email protected]

Y. Zheng e-mail: [email protected]

Y. Lei e-mail: [email protected]

X. Zhou e-mail: [email protected]

B. Liao e-mail: [email protected]

L. Huang e-mail: [email protected]

D. Huai e-mail: [email protected]

L. Yan e-mail: [email protected]

W. Wei e-mail: [email protected]

W. Chen e-mail: [email protected]

17 Page 2 of 14

Of the 92 QTLs, 15 QTLs were expressed across three environments and 65 QTLs were newly identified. Twelve QTLs for the height of main stem and the podand seed-related traits explaining more than 10 % of phenotypic variation showed a great potential for marker-assisted selection in improving these traits. The trait-by-trait meta-analysis revealed 33 consensus QTLs. The consensus QTLs and other QTLs were further integrated into 29 pleiotropic unique QTLs with the confidence interval of 1.86 cM on average. The significant co-localization of QTLs was consistent with the significant phenotypic correlations among these traits. The complexity of the genetic architecture of yield traits was demonstrated. The present QTLs for pod- and seedrelated traits could be the most fundamental genetic factors contributing to the yield traits in peanut. The results provide a good foundation for fine mapping, cloning and designing molecular breeding of favorable genes in peanut. Keywords Peanut (Arachis hypogaea L) . QTL analysis . Yield traits . Seed length . Seed width . Pod weight . Seed weight

Introduction The cultivated peanut (Arachis hypogaea L.) is a major oil and food crop in most tropical and subtropical areas of the world. The major objectives are to increase grain yield and improve the selection efficiency in peanut breeding. The yield traits, such as the height of main stem (HMS), total branching number (TBN), and the pod and seed traits directly contributed to grain yield in peanut (Holbrook and Stalker 2003; Shirasawa et al. 2012). The 100 pod weight (100PW), 100 seed weight (100SW), and shelling percentage (SP) are important index for grain yield of peanut. They were determined by pod and seed size, which were quantified by the pod length (PL), pod width (PW), pod thickness, seed length (SL), seed width (SW), and seed thickness (Shirasawa et al. 2012). The pod and seed traits are the most directly selected during long-term domestication and breeding. Research on them is also helpful to understand the yield component and evolution of crop species (Moles et al. 2005; Meyer and Purugganan 2013).

Mol Breeding (2017) 37: 17

Marker-assisted selection (MAS) is a powerful tool for enhancing selection efficiency. The basis of MAS is the construction of genetic map and the identification of the major genes or markers which are directly associated with the objective traits (Knapp 1998; Francia et al. 2005). Enough simple sequence repeat (SSR) markers have been developed in peanut research community for its virtues of transferable nature and the practical handle, and the construction of integrated consensus map also allowed a transverse comparison between populations (Gautami et al. 2012; Qin et al. 2012; Shirasawa et al. 2013; Huang et al. 2016; Zhou et al. 2016). The remarkable progress had been made on construction of genetic map and quantitative traits loci (QTLs) detection for important traits in peanut (Burow et al. 2001; Garcia et al. 2005; Gomez Selvaraj et al. 2009; Hong et al. 2010; Gautami et al. 2012; Qin et al. 2012; Shirasawa et al. 2012; Wang et al. 2012a, b, 2015a; Jiang et al. 2014; Pandey et al. 2014; Zhou et al. 2014; Huang et al. 2015, 2016). MAS had been used successfully in quality improvement, nematode, and rust resistance breeding in peanut (Chu et al. 2009, 2011; Sujay et al. 2012; Varshney et al. 2014; Wang et al. 2015a; Janila et al. 2016). Markers had also been developed for drought tolerance in peanut (Gautami et al. 2012). Up to date, the MAS for peanut yield breeding is still in progress. For the HMS and TBN, a number of QTLs were identified (Fonceka et al. 2012; Shirasawa et al. 2012; Jiang et al. 2014; Huang et al. 2015). QTLs for the pod and seed traits were also identified using linkage and association mapping. Eighteen QTLs for PL and 24 QTLs for PW were identified, which explained 1.24–28.2 and 1.13–22.3 % of phenotypic variation, respectively (Gomez Selvaraj et al. 2009; Fonceka et al. 2012; Shirasawa et al. 2012; Jiang et al. 2014; Huang et al. 2015; Chen et al. 2016a, b). Eighteen QTLs for SL and 15 QTLs for SW were located on different linkage groups, which explained 9.86–16.3 and 6.39–23.7 % of phenotypic variation, respectively (Gomez Selvaraj et al. 2009; Huang et al. 2015; Chen et al. 2016a, b). There were also 13 and 5 QTLs for SL and SW, respectively, which were identified using association analysis (Jiang et al. 2014; Pandey et al. 2014). Up to date, as an important index of exterior quality, the study on the ratio of PL to PW and SL to SW was not reported. For 100PW, six QTLs explaining 1.69–20.6 % of phenotypic variation were located on linkage A2, A8, B2, B3, and B5 linkage groups (Fonceka et al. 2012;

Mol Breeding (2017) 37: 17

Huang et al. 2015). The minor QTLs with 1.22– 2.26 % of phenotypic variation explained were also detected using association analysis (Jiang et al. 2014). For 100SW, six QTLs with 8.02–19.1 % of phenotypic variation explained were located on five linkage groups (Gomez Selvaraj et al. 2009; Fonceka et al. 2012; Shirasawa et al. 2012; Huang et al. 2015), and five QTLs with 12.73–26.08 % of phenotypic variation explained were identified using association analysis (Pandey et al. 2014). The genetic basis of seed size and weight were also characterized well in other crops and the pleiotropy for these traits was observed. In soybean, the QTLs for seed size and weight were distributed over 16 chromosomes and many of them explained minor or moderate phenotypic variation (Salas et al. 2006; Xu et al. 2011; Niu et al. 2013; Kato et al. 2014). The QTL X40, X53, X83– X85, and X92 simultaneously controlled two of the SL, SW, and 100SW (Xie et al. 2014). In rice, the genes controlling seed-related traits acted in independent genetic pathways. The GS3, GS5, qGL3/qGL3.1, GW5/ qSW5, and TGW6 simultaneously affected two or more seed size traits and grain quality, indicating a complex genetic basis of these traits (Zuo and Li 2014; Peng et al. 2015; Wang et al. 2015b). The tight linkage and pleiotropy were observed on the silique length, seed size, and weight in Brassica napus. The gene ARF18 simultaneously affected seed weight and silique length (Liu et al. 2015). Three major QTLs simultaneously controlled seed weight and silique length (Yang et al. 2012; Li et al. 2015). However, compared to the comprehensive QTLs research on rice, soybean, and oilseed rape, it is still far from comprehensively understanding the genetic factors controlling these traits in peanut (Varshney et al. 2010, 2013; Pandey et al. 2014). Recently, the genome sequences of the two diploid species of Arachis were published (Bertioli et al. 2016; Chen et al. 2016a, b) which are the common ancestors of cultivated peanut. It would be greatly facilitated for genetic map validation and QTLs mapping. The main objectives of the present study were (1) construction of a genetic map using recombinant inbred lines (RIL) population; (2) mapping of QTLs for yield traits including the HMS, TBN, pod- and seed-related traits; and (3) comparative mapping of QTLs between the physical and the genetic maps. These results would provide comprehensive information for QTL comparison among similar studies of peanut, MAS of yield breeding, and mapbased cloning of the candidate genes in peanut.

Page 3 of 14 17

Materials and methods Plant material A RIL population including 188 lines was developed by single seed descent from a cross between female parent Fuchuan Dahuasheng (A. hypogaea L. subsp. hypogaea L . v a r. h i r s u t a K o h l e ) a n d m a l e p a r e n t ICG6375(A. hypogaea L. subsp. fastigata Waldron var. vulgaris Harz) (Fig. S1). Fuchuan Dahuasheng is a local cultivar in Guangxi Province of southwest China, ICG6375 is introduced from International Crops Research Institute for the Semi-arid Tropics (Fig. S1). This combination was chosen because of their big differences in major agronomic traits including the pod’s shape and size; seeds’ shape, size, and coat color; flowering time, number of seeds per pod (50 % pods of Fuchuan Dahuasheng were three-seeded); and total branching number. The phenotypic value of traits in the present study was listed in Table S1. DNA extraction The young unopened leaves were collected from each accession for genomic DNA extraction using a modified CTAB method. Total DNA was digested with RNase and then quantified with a Beckman DU-650 spectrophotometer. Agarose gel electrophoresis was employed for evaluation of quality and integrity. The PCR product were visualized by 6 % denaturing polyacrylamide gel (PAGE) followed by silver staining as described by Chen (2008). Field trials and trait phenotyping The population comprising 188 lines was grown in 2013 (F5) and 2014 (F6) seasons. Ten individuals of each line were evaluated for the HMS, TBN, PL, PW, PL/PW, SL, SW, SL/SW, 100PW, 100SW, and SP. Phenotyping of F5 RIL population was carried in Wuchang (environment 1, E 114° 34′/N 30° 59′) and Yangluo (environment 2, E 114° 52′/N 30° 59′) in 2013, and phenotyping of F6 RIL population was carried in Wuchang in 2014 (Environment 3, E 114° 34′/N 30° 59′) for two trials. The two-seeded pods and seeds were measured by their length and width using a parallel rule. The pods and seeds weight were measured by an electrical scale. The PL/PWand SL/SW were calculated using PL and PW and SL and SW. The SP was calculated as 100PW divided by the seed weight of 100 pods.

17 Page 4 of 14

Statistical analysis Phenotypic data were tested for normality using the PROC UNIVARIATE procedure of SAS 9.3 (SAS Institute, Cary, NY, USA). Correlation coefficients among 11 traits were estimated using the PROC CORR procedure of SAS. The broad-sense heritability for each trait was calculated by a method described by Wu et al. (2009). SSR markers analysis A total of 8456 SSR primer pairs were employed to detect their polymorphism between parents of RIL population. The selected polymorphic markers were assigned to genotype the RIL mapping population. PCR was performed in a total volume of 10 μl of reaction mixture consisted of 0.4 mM dNTPs, 2.0 mM Mg2+, 1× Taq buffer, 0.5 U Taq DNA polymerase (Tiangen), 0.3 μM primers, and 25 ng genomic DNA. Amplification was carried out on a T100™ Thermo cycler (BIO-RAD). The parameters of thermal cycle were 95 °C for 3 min, followed by 10 cycles of 94 °C for 30 s, 65 °C for 30 s (decreasing in decrements of 0.5 °C per cycle), and 72 °C extension for 60 s; then followed by another 30 cycles of 94 °C for 30 s, 55 °C for 35 s, 72 °C for 60 s, and final 72 °C for 5 min. Construction of genetic map and comparison between the genetic and physical maps Genotyping data obtained from the RIL population was used for linkage analysis using JoinMap 4.0 (Van Oojen and Voorips 2006). Segregation distortion of each marker was examined at an expected 1:1 using chi-squared (χ2 ) testing under BLocus genotype frequency^ function. The BKosambi mapping^ function (Kosambi 1944) was used to transform recombination frequencies into centiMorgan (cM) map distances with a minimum log-of-odds (LOD) threshold of 4.0 and recombination frequency of 0.35. Marker order in groups was corrected under the BCalculate Map^ command. The comparison between the present genetic map and the physical map was based on the markers assigned to the chromosomes of A genome of Arachis duranensis and B genome of Arachis ipaensis (Bertioli et al. 2016, http://peanutbase.org/download). The corresponding genomic position of markers on the present genetic map was determined using blast/e-PCR

Mol Breeding (2017) 37: 17

with the primer sequence against the physical map of A. duranensis and A. ipaensis. The comparison between the present genetic map and the integrated consensus map (the INT map, Shirasawa et al. 2013) was based on the position of common markers. The compared map was generated using MapChart for Windows (Voorrips 2002). QTLs mapping The respective data of single environment and the averaged data of each 2013 and 2014 were combined with marker genotyping data to identify QTLs, respectively. The average value of phenotyping data across three environments was also employed to detect the QTLs which were expressed across the three environments. QTL analysis was conducted by composite interval mapping (CIM) using Windows QTL Cartographer 2.5 (Wang et al. 2012a, b). The LOD value was chosen to be 2.9 to declare a QTL significant based on a permutation test (Churchill et al. 1994) with 1000 runs to determine the P = 0.05 genome-wide significance level. To evaluate more potential QTLs, the QTLs with a LOD value more than 2.5 were also referenced. Based on the QTL detection above using Windows QTL Cartographer 2.5, algorithms for QTL meta-analysis were used to estimate the number and position of QTLs which were located at same region on the linkage group (Goffinet and Gerber 2000). This approach provides a modified Akaike criterion that can be used to determine the number of meta-QTL that best fitted the results on a given linkage group. The two-round strategy of QTL meta-analysis was performed to integrate the QTLs using BioMeractor2.1 software (Arcade et al. 2004; Chardon et al. 2004). The BioMeractor2.1 set a 95 % confidence intervals 3:92 ffi pffiffiffiffiffiffiffiffiffiffiffiffiffiffi i¼k 2 , for the meta-QTL using the formula C:I:¼ ∑i¼1 ð1=S i Þ

S 2i

presents the variance of the position of where QTLi, and k presents the total number of QTL integrated into the meta-QTLs. Firstly, the additive QTLs expressed in different environments for the same trait were integrated into the consensus QTL. Secondly, the overlapped consensus QTLs and the other QTLs for different traits were further integrated into the unique QTL. The QTL nomenclature was referenced to the research reported by Shi et al. (2009). For the consensus QTL, the QTL nomenclature corresponded to q (abbreviation of QTL), followed by the abbreviation of

Mol Breeding (2017) 37: 17

Page 5 of 14 17

values of all traits across three environments were used in subsequent variance analysis. All traits showed continuous normal or near-normal distribution, indicating a polygenic inheritance of them and an ideal model for QTL analysis (Fig. 1). All traits showed high broad sense heritability ranged from 0.75 to 0.95 except HMS (0.40) and SP (0.59), and SW showed the highest heritability. The correlation coefficient analysis discovered a significant correlation among most of the traits (Table 1). The PL, PW, SL, and SW showed significant positive correlation (P < 0.01) with 100PW and 100SW, and PL showed significant negative correlation (P < 0.01) with SP. The HMS and TBN were relatively independent of other traits.

the trait (e.g., HMS); then, the linkage group number; and finally, the serial number. For the unique QTL, the nomenclature corresponded to uq, followed by the linkage group number, and then the serial number.

Results Phenotypes of mapping population The population was phenotyped extensively for 11 traits. For all of the traits except SW and SP, significant differences were observed between Fuchuan Dahuasheng and ICG6375 (Table S1). The average HMS 50

60

TBN

60

PW

PL 50

50

50

40

40

30

30

20

20

10

10

10

10

0

0

0

0

40

40

30 20

PL/PW 50

70

SL

70

30 20

SL/SW

SW 60

60

60

40

50

50

30

40

40

30

30

20

20

10

10

10

0

0

0

20 10 0

70

SP

50

40 30

50

40

40

30

30

20

20

10

10

0

0

20 10 0

30 20

100SW

100PW 50

60 50

40

Fig. 1 Phenotypic distribution of 11 yield traits in the RIL population

17 Page 6 of 14

Mol Breeding (2017) 37: 17

Table 1 Correlation analysis of the 11 traits in present population across three environments Correlation

HMS

TBN

HMS

PL

*

PW

PL/PW

SL

SW

SL/SW

100PW

100SW

SP

ns

ns

ns

ns

ns

ns

ns

ns

ns

ns

ns

ns

ns

ns

ns

ns

ns

ns

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

**

ns

ns

ns

**

**

**

**

ns

**

**

**

ns

**

**

ns

**

**

TBN

0.180

PL

0.004

−0.088

PW

0.046

−0.130

0.725

PL/PW

−0.058

0.109

0.463

−0.267

SL

−0.011

−0.085

0.856

0.733

0.247

SW

−0.053

−0.081

0.482

0.820

−0.376

0.604

SL/SW

0.043

0.065

0.559

0.084

0.670

0.620

−0.246

100PW

0.002

−0.057

0.828

0.881

0.017

0.837

0.786

0.234

100SW

−0.022

−0.105

0.721

0.863

−0.101

0.851

0.857

0.329

0.914

SP

−0.106

−0.104

−0.259

−0.224

−0.063

−0.136

−0.019

−0.046

−0.326

ns −0.087

ns non-significant at P < 0.05 *Significant at P < 0.05;**Significant at P < 0.01

Characteristics of genetic map

mapped onto the 26 linkage groups with an average of 23.4 loci per group (Fig. 2, File S1). The map covered a total of 1557.48-cM genetic distance ranging from 14.25 to 95.24 cM with averaged 59.90 cM for each group and 2.56 cM between adjacent markers (Fig. 2). A

A total of 718 of 8456 SSR markers were polymorphic between Fuchuan Dahuasheng and ICG6375. Under the threshold of more than 6 loci for grouping, 609 loci were 0.0 4.5 7.1 9.7 10.5 11.6 13.7 14.6 15.3 15.9 16.3 16.9 17.6 17.9 18.1 18.4 18.9 19.7 20.1 20.8 21.5 22.2 22.7 23.5 24.1 25.2 26.1 28.2 29.8

A1

AGGS874 AHGS1846 AGGS633 AHGA164454 AGGS940 AHGA284837 AHGS1974 AHGA329002 AHGA164439 ARS789 ARS731 AHGA164444 AHGA164437 AGGS1572 A01A246 AGGS60 AHGS2535 AHGS1646 AGGS1213 AHGS1668 AHGA367823 TC4H02 GA110 AHGA368803 AHGA374228 AHGA374249 TC4G05 GM2607 AGGS331

A2 0.0

AHGA23452

3.3

AGGS1413

6.5

ARS741

13.2

pPGSseq11G3

16.6

AGGS1404

21.3 24.1 25.6

AGGS679 AGGS1085-2 AHGS2026

32.4

AGGS1168

35.9 38.1

AGGS1472 AHGA362472

45.2 47.0

AGGS638-2 AHGA362499

55.2

pPGSseq11H1-1

62.0

0.0 6.8 8.9 9.9 14.6 16.8 17.7 18.9 19.8 20.5 21.2 22.0 22.8 23.6 24.7 25.3 26.8 30.4 32.5 33.9 34.9 36.4 38.0 43.6 46.6 49.2 53.2 55.3 57.4 62.7

A3

AHGS2019 TC27H12-2 AHGS1813 AGGS875 AGGS21 TC3H02 AHGA335258 AGGS1578 pPGPseq8E12 AGGS925-2 TC23C08 AGGS925-1 AGGS947 AHGA364951 AHGS3715 AHGS1130 AGGS188 AHGA364936 AHGA11190 TC27H12-1 TC1A08 AGGS1376 AHGA24894 AHGS2374 ARS545 AHGS3627 AHGS1829 AGGS707 AGGS1287 AHGS2579 AGGS408 AGGS1612-1 AGGS1110 AHGS1462 AHGA25786 AGGS1056 AC2C08 ARS627-2 AGGS1363 AHGS1524 AGGS927 AGGS305 AGGS1395 AGGS1122 AGGS992 AGGS505 AHGA68610

A4

AHGA170478 A04A168 AGGS383 AGGS477 AHGS1802 AHGS1679 GNB533-1 A04A301 AGGS1547 AGGS1326 AGGS1520 AHGA354327 AHGA220404 AHGA38598 AHGA36571 AGGS580 AGGS570 POCR536 AGGS1414 AGGS1458 pPGSseq18C5 AGGS1398 AHGA331142 AHGA331177 AGGS1612-2

68.6

A04A552 AGGS1211 GM2637

84.3

GM2673

61.4 64.8

pPGSseq9H8

76.0

B1

0.0 1.0 2.7 4.0 5.1 8.4 9.6 11.2 12.0 13.0 14.1 14.8 15.7 17.0 17.8 19.0 20.3 22.9 25.4 26.9 30.9 34.1 35.5 37.6 44.2

AGGS1608

AHGS2429

39.0

0.0 2.6 5.5 8.1 11.2 14.7 16.0 18.0 19.8 27.7 33.2 45.8 53.6 55.8 57.9 59.0 60.3 60.8 62.0 65.1 66.7 67.7 68.5 69.9 70.6 72.0 72.5 73.3 73.6 74.7 75.1 75.6 75.9 76.2 76.7 78.1 79.0 80.2 81.4 82.0 83.4 85.3 86.1 88.7 90.4 92.1 93.2

AGGS981 AGGS1270 AHGS1672 AGGS1257 AGGS1132 AGGS435 AHGS1792 AGGS933-2 AGGS1540 AHGS3647 AHGA40118 AGGS1374 ARS761 AHGS1106 AHGA40112 AHGA62112 GNB167 A03A375 GNB675 A03A335 AHGS1847 AGGS1483 AHGS1674 AGGS986 AHGS0132 AGGS1149 PM238-2 AGGS1139 AHGS2025 TC1E06

B2 0.0 1.7 8.8 13.2 13.7 14.8 15.7 17.7 18.9 19.9 20.5 21.8 22.3 23.1 23.5 24.1 25.0 26.0 26.8 27.7 28.5 30.2 31.5 33.2 34.4 35.8 38.8 40.5 41.0 43.2 47.6 48.5 53.6 58.9 67.8 75.1

AGGS1409

0.0

AHGA96496

4.9

AGGS1199

9.8

AGGS1119

3.1 5.7 6.4 8.4

AHGA96466 AHGA96458 AHGA96464 AGGS283

AHGS1429 AGGS933-1 11.0 AHGA47958 14.3 AGGS416 AGGS57 AHGS1722 0.0 AGGS1289 AHGS1140 4.8 AGGS1407 8.2 PM238-1 AGGS270 AGGS325 16.1 AGGS717 GA27 AGGS954 23.9 AGGS1325 AHGA300244 pPGPseq3F5 AHGA362520 39.9 AHGA300241 AGGS1083 AHGS2526 AGGS1010 0.0 ARS534 2.3 pPGPseq3E10 2.6 AGGS1161

A6

AGGS1622-2 ARS152-2

B32

AHGS1561 AGGS1276

0.0

9.8

AHGA26512 AGGS351 ARS777 ARS202 ARS715 ARS760 AGGS1187 AGGS671 AGGS438 AGGS738 AHGA44674-1 GNB533-2 AGGS1457 AHGA9104 AHGS0599 AGGS984 AHGS1615 A05A295 PM45 AGGS1525 TC2B09 ARS702 AGGS1521 AGGS1335 AHGA44691 AGGS451

TC11H06

AGGS706

14.1

A04A55

18.0

AHGS0288

20.9

AHGA161495

25.2

AGGS461

34.1 35.9 37.2 38.7 42.1 43.5 45.5

AGGS186 AGGS1101 AHGS1382 pPGSseq14D11 AHGS1579 AHGS2550 AGGS1601

AGGS1369 AHGA227117 AHGA105549

AGGS1016

B33

AGGS1622-1 ARS152-1 ARS152

13.6

IPAHM103

17.7

AHGS1788

21.9

AGGS824

27.5

AGGS315

33.3

AHGS0626

0.0

AHGS1507

B4

0.0

14.5 17.8 21.5 24.4 29.4 30.7 34.6 37.3 39.8 41.6 43.7 45.9 49.0 50.1 51.3 52.4 57.2 58.7 60.2 61.7 63.6 69.8 71.2 76.1 79.6 84.6

13.7 14.7 15.8 19.7 23.4 25.4 26.1 27.3 28.4 30.0 31.0 32.1 33.2 34.6 35.8 37.3 38.2 38.6 39.0 39.7 40.7 43.2 45.8 46.8 47.5 52.2

B31

B3 AGGS1468 pPGSseq4E8 AGGS1085-1 AGGS78 AGGS1356 AHGS1853 AGGS399 AHGS2557 AGGS1484 AGGS1238 AGGS1175 AHGS1530 AHGS1393 AGGS1461 AHGS1928 AHGS2736 AHGS2200 AC2C05 AHM022 AGGS1532 AGGS1445 AHGS2088 AHGS1845 AHGS1419 AGGS1399 AHGS0357 AGGS1400 AGGS638-1 pPGPseq11H1-2 AHGS3626 AHGS1940 AGGS231 AHGA327396 AHGS2344 GM2808 AGGS1092

A5 0.0

53.2

ARS542

56.5

AGGS81

62.2

AGGS890

16.4

AGGS705

33.6

AGGS1479

37.1 41.7 43.8 45.8 47.3 49.1 50.1 51.6 53.2 56.9 59.3 61.6 63.8

AGGS750 AGGS1450-2 AGGS997 ARS816 AHGA193640 AGGS500 A06A285 GM623 AHGA193642 AGGS1209 AGGS1314 AHGA141434 AGGS1405

72.2

0.0 12.1 16.2 19.7 24.0 25.9 27.5 30.3 31.8 33.0 34.3 35.4 36.5 37.6 38.9 40.4 41.6 43.8 44.5 46.2 48.3 49.1 50.3 53.1 55.0 56.8 66.2 70.3 76.0 77.4 78.8 81.5 83.4 85.4 89.1 95.2

B5

0.0 2.9 4.5 8.3 11.7 12.4 13.8 15.8 16.3 16.9 17.8 18.8 19.3 19.8 20.4 20.7 21.1 21.7 22.1 22.6 23.0 23.6 24.3 25.7 26.6 27.1 29.0 29.7 32.0 32.3 34.2 35.5 37.5 38.9 42.5 48.7 58.3

AGGS1143

PM377

AHGA152207 AHGS1391 AGGS973 AHM082 AGGS1481 AGGS243 AGGS1299 AHGS2141 AHGS1925 AHGS1914 AGGS1262 AGGS244 AGGS1464 Aill9F10 IPAHM282 AHGS2395 PM36 AHGS1682 TC6E01 AGGS979 AHGA221480 AHGA221476 AGGS1129 AGGS780 AHGA68622 ARS627-1 AGGS670 AGGS931 AHGA363495 AHGA363491 IPAHM73 EM-106 ARS173 AGGS729 AGGS1422 AGGS955

A7

B6 0.0 5.9 8.4 11.0 16.0 18.1 21.1 22.7 24.5 26.1 27.9 29.0 30.7 32.2 33.3 35.1 36.2 37.4 39.5 40.6 41.6 43.1 44.4 48.1 48.7 50.1 51.4 53.1 55.6 62.6 63.8

A71

AHGS1912 GM1494 AHGA102053 AGGS1365 AHGS1954 AHGS1475 AHGS1436 PM204 AHGS1692 pPGSseq3A1 AHGS2478 ARS127 AHGA65329 AHGA65328 AHGA65333 AGGS987 AHGS1980 PM358 TC9C08 ARS753 AGGS410 AHGA363105 AGGS1123 AGGS1389 AHGS2153 AGGS1392 AGGS1162 AHGS2413 AHGS0546 AHGSO346 AGGS1585 AHGA373382 AGGS1433 AHGS2298 AGGS1577 AGGS1237 AHGS0228

EE-95

7.2

AGGS847

10.4 12.0 13.3

AGGS846 AGGS831 AHS0035

16.3

AGGS848

22.8

AGGS836

30.9

AHS0001

34.6

ARS734

46.1

AGGS1272

B7 AGGS1455 PM210 AGGS2482 AGGS1091 AGGS50 AGGS138 AHGS1900 AGGS1437 AGGS1517 AGGS1267 AGGS1105 AHGS2093 AHGS1930 AHGS1361 AGGS1450-1 AGGS1467 TC24B05 AGGS1093 AHGS1742 AHGA105542 TC4F10-2 AHGS1960 AGGS1541 AHGS2553 AGGS237 AHGS1756 AHGA367938 AGGS1539 AGGS1397 GM2444 AHGS1511

A81

A8

0.0

0.0

AHGS1444

8.9

AGGS1077

15.9

AGGS1197

20.8

AHGA109663

26.7

A09A39

34.1

AHGS1383

39.1

pPGPseq3A8

44.6

AGGS704

48.1

AHGS1434

52.2 53.6

ARS120 A08A119

57.4

AGGS950

65.0

AGGS1495

68.5

AHGS2319

B8

0.0

TC2G05

8.0

AHGA59797

14.1 15.9 17.3 18.4 20.4 21.7 23.1 24.8 26.5 28.3 30.7 33.6 37.8 42.1

AHGS1915 AGGS522 AHGS1952 TC9H09 AHGA109692 AGGS1021 AHGA115756 AHGS2471 AGGS1631 AHGS2419 STS20-1 Ahm054 AGGS662 AHGS1864

47.2

TC3B04

56.1

TC7A02

62.5

AGGS965

67.3

AGGS967

0.0

AGGS1605

4.4 9.6 11.7 12.3 13.7 14.6 15.7 16.6 17.1 18.5 19.5 20.4 21.3 22.1 23.3 24.5 25.2 26.7 28.3 29.8 32.7 34.8 36.5 37.7 42.9

AHGA65348 TC6H03 AHGS1539 AHGS2147 AHGS2751 AHGS2069 AHGS1385 AGGS233 AGGS249 AGGS1452 GA80 AHGS1510 AGGS1215 AGGS1342 AHGS2111 AHGA7426 AHGA7429 AGGS1576 AHGA7413 AGGS1212 AGGS189 AGGS1323 AGGS443 AHGS1141 AGGS949

47.4

AGGS930

A9

0.0

AGGS953

7.8

AHGS1578

11.3

AGGS400

19.0

AGGS184

29.3

AGGS1383

38.8

AGGS1428

48.3

AHGA301840

0.0 7.7 13.7 14.9 17.1 18.0 19.2 21.0 21.6 22.4 23.4 24.3 24.8 25.7 26.8 27.6 28.3 28.6 29.6 30.5 31.0 31.9 32.6 33.4 34.3 35.0 36.0 36.9 37.6 38.8 39.4 41.1 43.0 47.2 48.4

B9

A10 AGGS1463

0.0

AHGA364794

0.0

10.4 11.2 13.0 16.0 16.9 19.9 21.3 22.7 23.5 24.4 26.4 27.5 28.9 29.9 30.8 32.2 33.0 35.2 36.6 37.9 40.6 42.0 47.8 50.9 54.1 58.1

AHGA72170 AGGS694-2 AHGS0695 AHGS2635 GM2839 AGGS1317 AGGS1193 AGGS1248 AGGS385 ARS730 AGGS1498-2 AGGS449 GM2792 AHGA98576 EM-87 AGGS1128 ARS742 ARS766 AHGS1478 AGGS957 AHGA57727 AHGA22895 ARS205 AGGS389 AHGS1683 TC1D02

6.7 11.6 14.2 15.0 16.1 17.6 19.2 20.3 21.3 23.0 23.5 24.3 24.9 25.4 26.3 26.7 27.4 28.6 29.5 31.1 31.6 31.9 33.7 34.6 35.8 38.7

TC1G04 AGGS1301 AHGA364906 pPGPseq4H11 AHGS1390 AGGS299 AGGS805 AHGS1777 AHGS2024 AGGS1015 AHGS1606 GM2084 AGGS405 AHGS0648 GM692 AGGS357 AGGS1147 AGGS1345 AHGS2297 ARS710 AHGA195553 TC31G11 AGGS1185 AGGS1203 AGGS751 AHGA195528

68.1

AGGS285

45.6

AGGS1232

1

AGGS1598 AHGS1969 AGGS1380 AGGS308 AHGS1156 AGGS93 AGGS428 AHGS2325 POCR164 AHGS1998 AGGS37 AGGS1239 AGGS347 GM2553 AHGS1634 TC4F10-1 AHGS1790 AGGS1308 AGGS1300 AHGS1754 AHGS2529 AGGS358 pPGSseq14H6 AGGS1498-1 AGGS1384 AHGS1394 AGGS1095 PM436 Gi-1107 IPAHM166 AGGS309 AGGS1253 AHGS2324 AHGS1804 AHGS1715

0.0

8.4

B9

AGGS1434

AGGS618

20.4

AHGA352265

24.7

AHGA352202

28.5

AGGS694-1

31.7

AHGA352262

41.0

AHGA352263

0.0

B10

EE-16

7.0

AHGS2543

18.0 24.4 33.5 35.1 39.3 44.0 45.2 46.3 49.5 50.5 52.3 52.8 53.3 53.8 54.1 55.0 55.4 55.8 56.5 56.7 57.6 58.2 60.0 61.2 63.2 64.7 65.8 71.1

AGGS1194 AGGS1355 AHGS1395 A01A157 AGGS1160 AGGS445 AGGS1420 AGGS1442 AHGS1564 AGGS1359 AGGS1118 AGGS675 AGGS1176 AHGS2096 AGGS1273 AHM062 IPAHM23 TC3E05 AHGA75538 AHGA75537 AHGS1696 A10A280 AGGS969 AGGS58 AGGS1393 TC7H11 GM690 AHGA61572

Fig. 2 Genetic linkage map and QTL detection of the 11 traits in an RIL population derived from Fuchuan Dahuasheng × ICG6375

Mol Breeding (2017) 37: 17

Page 7 of 14 17

total of 321/609 (52.71 %) loci showed distorted segregation ration (P < 0.01), 238/321 (74.14 %), and 83/321 (25.86 %) loci skewed toward female and male parents, respectively. The groups were designated as A1–A10 and B1–B10 for A and B genome by assigning common markers to the INT map (Table S2, File S2, Figs. S2 and 4). Especially, more than one group shared the same designation and were marked using superscript number when they were assigned to the same one group on INT

map. As a result, the short group A71 was attached to A7 group, B31, B32, and B33 were attached to B3 group. Homologous relationships of FI map to the physical map and the INT map The availability of the genome sequence of diploid species of Arachis is helpful to compare the present genetic map and the physical map in A. hypogaea

Table 2 Summary statistical of the present linkage map Linkage group Markers Genetic distance Marker Anchored Recombination Coverage rationd Number of marker distorted to frequencyc number densitya markers numberb Fuchuang ICG6375 Total A1

30

39.02

1.3

24

2.83

0.79

2

16

18

A2

15

62.02

4.43

11

1.37

0.85

3

3

6

A3

31

75.96

2.53

14

1.80

0.83

4

5

9

A4

29

84.3

3.01

16

2.18

0.66

8

1

9

A5

28

52.19

1.93

18

1.44

0.68

0

12

12

A6

17

72.2

4.51

11

1.47

0.80

4

6

10

A7

38

58.25

1.57

23

1.17

0.85

16

2

18

A71

10

46.14

5.13

9

0.06

0.04

3

0

3

A8

14

68.51

5.27

10

0.88

0.37

7

1

8

A81

7

48.27

8.04

4

0.16

0.16

1

1

2

A9

28

68.13

2.52

18

1.53

0.87

14

0

14

A10

29

45.63

1.63

17

2.68

0.85

13

0

13

B1

47

93.2

1.98

30

1.44

0.95

14

5

19

B2

36

75.12

2.15

28

1.32

0.91

16

5

21

B3

30

84.61

2.92

25

1.29

0.80

26

0

26

1

B3

7

14.25

2.38

5

0.29

0.02

1

2

3

B32

6

39.92

7.98

5

0.10

0.03

0

0

0

B33

8

33.34

4.76

6

0.09

0.02

4

0

4

B4

16

62.21

12.44

11

1.90

0.88

2

4

6

B5

36

95.24

2.72

28

1.38

0.88

20

0

20

B6

31

63.82

2.13

24

1.84

0.84

0

12

12

B7

20

67.32

3.54

24

1.60

0.85

12

2

14

B8

27

47.37

1.82

17

3.14

0.88

5

6

11

B9

35

48.4

1.42

29

2.71

0.89

33

0

33

B91

7

40.96

6.83

6

0.06

0.01

1

0

1

B10

30

71.1

2.45

20

1.69

0.88

29

0

29

Total

612

93.2

2.13

437

1.40

0.65

238

83

321

a

Marker density is the ration of markers number to genetic distance

b

Anchor markers number represents the number of markers which was anchored to the sequenced A genome of A.duranensis and B genome of A.ipaensis c d

Recombination frequency is the ratio of physical distance to the genetic distance

Coverage ration is calculated as the covered physical distance (Mb) of linkage group divided by the length (Mb) of the corresponding chromosome

17 Page 8 of 14

Mol Breeding (2017) 37: 17

Fig. 3 Comparison of the present map and the physical map of A and B genomes of diploid species A. duranensis and A. ipaensis. Collinear SSR loci between linkage groups of present map and the counterpart chromosomes of A and B genome were indicated by lines

(Bertioli et al. 2016, http://peanutbase.org/download). A total of 428 loci presenting 1432.86 cM were successfully aligned to the chromosomes and covered 2007.86 Mb with an average of 1.4 Mb/cM. There were 175 loci representing 626.43 cM and 253 loci representing 806.43 cM were assigned to A and B genomes, and covered 818.0 and 1189.86 Mb, respectively, (Table 2, File S3 and 4). The high collinearity was also observed between the present and the physical maps (Fig. 3, Fig. S3, File S3 and 4). Generally, the good collinearity was observed between A3 and A03 , A8 and A08, A9 and A09, B5 and B05, B8 and B08, B9 and B09, and B10 and B10 (Fig. 3, Fig. S4). For a huge and poorly characterized genome of allotetraploid A. hypogaea, it would be helpful to well understand the genetic make-up of tetraploid cultivars of peanut by comparison between the present map and the INT map which showed the highest marker density and the largest genetic distance (Shirasawa et al. 2013). The good agreement of collinear was observed between the present map and the INT map (Fig. S1, File S2). Twenty-four groups including the attached groups of A7, A8, and B3 were homologous to the 20 groups on the INT map. For the A2 group, only one marker was assigned to A02 group of the INT map. A total of 167 loci covered 904.63/714.96 cM in the present map and the INT map, and the collinear region was narrowed in the INT map. Ten groups of 388.0 cM with 72 loci were assigned to A genome, and 12 groups of 516.63 cM with 95 loci were assigned to B genome (Fig. S2 and 4, Table S2, File S2). The fine collinearity was observed between A1

and A01, A7 and A07, A10 and A10, B2 and B02, and B6 and B06 (Figs. S2 and 4). QTLs detection and meta-analysis For all 11 traits, a total of 175 QTLs were detected (File S5). And 92 of them were defined as repeatable QTLs for they were detected in more than two environments or using average data or their overlapped confidence interval for different traits (File S6). The identified QTLs explained 5.1–13.2 % of the phenotypic variation. Eleven QTLs were identified as major QTLs with more than 10 % of phenotypic variation explained, and 15/92 QTLs were also meanwhile detected across three environments (File S6). Where the confidence interval of the QTLs for the same trait in different environments overlapped, meta-analysis was employed to integrate these QTLs into consensus QTL. As a result of the first round of metaanalysis, the 70 QTLs were integrated into 33 consensus QTLs (Table. S3, File S7). For example, eight QTLs for SL were integrated into three consensus QTLs on A7 group (Fig. S5). In the second round of meta-analysis, besides the consensus QTLs, a total of 53 single QTLs for different traits detected were also employed because of their contribution to the correlated traits and the same or very close position on the linkage group. Finally, the 30/ 33 consensus QTLs and the 53 QTLs were further integrated into 29 unique QTLs (Table S3, File S8). For example of the integration of unique QTLs uqA7–1, the five QTLs including three single QTLs qPLA7.1a, qPL/ PWA7.1a, and q100SWA7.1a, and two consensus QTLs qSLA7.1a and q100PWA7.1a, were further integrated

Mol Breeding (2017) 37: 17

into the uqA7–1 (Fig. S5, Table S3). A total of ten unique QTLs were integrated by the QTLs detected only once, and the remaining 19 unique QTLs were integrated by both the consensus and the single QTLs (Table S3). The average confidence interval of all unique QTLs was reduced from 3.03 to 1.86 cM, and the average genetic distance of linked marker to the peak position was 0.45 cM. A total of 29 markers linked to the unique QTLs were also identified (File S9). The QTLs were mainly located on A2, A3, A5, A6, A7, A10, B1, B2, B4, B6, B7, and B8 groups (Fig. 2). As a result of the two rounds of meta-analysis, the significant co-localization of QTLs was observed especially for pod and seed traits (Fig. 2, Table S3). A total of 83 QTLs were co-localized at 29 pleiotropic loci which were distributed on the 11 linkage groups. And 13 of 29 unique QTLs controlled at least three different traits. For example, the uqA5.1a, uqA5.1b, uqA7.1a, and uqB1.1a represented QTL for five traits. The co-localization of these QTLs was in accordance with the signs of the genetic-correlation coefficients of pod- and seedrelated traits. The HMS and TBN are the most important plant architecture traits. For the HMS, a total of six QTLs were detected repeatedly. And four consensus QTLs were detected and explained 5.95–9.14 % of phenotypic variation. Among them, one was on A6 linkage group and three were on B6 linkage group (Fig. 2, Table S4). QTL qHMSA6.1a, qHMSA6.1b, qHMSB6.1a, and qHMSB6.1b were detected across three environments. All four QTLs showed negative additive effect (Table S4). For the TBN, three QTLs, qTBNB1.1a, qTBNB6.1b, and qTBNB9.1a were detected using pooled data of E1 and E2, explaining 5.08–6.03 % of phenotypic variation (Fig. 2). The qTBNB1.1a and qTBNB9.1a showed positive additive effects of 2.07 and 2.50, while the qTBNB6.1b showed negative additive effect of 2.20 (File S5). In the present study, 83 of 92 QTLs were detected repeatedly for pod- and seed-related traits, and 11 of them were also detected across three environments. Finally, 30/33 consensus QTLs and all 29 unique QTLs were all accounted to these traits (Table S3, Table S4). Almost all the QTLs for pod- and seed-related traits were involved into the generation of unique QTLs, except the qPWA6.1a, qSWA8.1a, and qSPA10.1a (Fig. 2, Table S3). For the PL, 12 QTLs explained 5.6–13.2 % of phenotypic variation, the qPLA7.1a was detected across three environments, and four

Page 9 of 14 17

consensus QTLs were located on A5 and A7 groups. Seven QTLs explained 5.6–8.6 % of phenotypic variation of PW, the qPWB8.1a was detected across the three environments, and two consensus QTLs were located on B8 group. Thirteen QTLs explained 5.4–12.1 % of phenotypic variation of PL/PW, only one consensus QTLs qPL/PWB7.1a was located on the B7 group. Twelve QTLs explained 5.44–13.2 % of phenotypic variation of SL. And as many as seven consensus QTLs for SL were generated, and the consensus QTL qSLA7.1b was integrated by four QTLs. Ten QTLs explained 6.0–12.8 % of phenotypic variation of SW, three of which were consensus QTL and located on A10 and B2 groups. For the SL/SW, nine QTLs explained 5.5–11.8 % of phenotypic variation. Of them, the qSL/ SWA2.1a and qSL/SWA3.1a were expressed across three environments, and the qSL/SWA7.1a and qSL/SWA7.1b were consensus QTL. Interestingly, unlike the qPLB1.1a, qPWB1.1a, and qPL/PWB1.1a were detected independently at the same one locus, although the qPLB8.1a and qPWB8.1a were detected at the same one locus on B8 group; however, the QTL for PL/PW was not detected on B8 group. The 100PW, 100SW, and SP are the most direct yield traits. For the 100PW, eight QTLs explained 5.3–11.2 % of phenotypic variation. QTLs q100PWA7.1a and q100PWB7.1a were detected across three environments, the q100PWA5.1a, q100PWA7.1a, and q100PWB8.1a were consensus QTL. For the 100SW, nine QTLs explained 5.2–10.8 % of phenotypic variation, three of which were consensus QTL. The QTLs q100SWA7.1a, q100PWA7.1b, q100PWB6.1a, and q100PWB8.1a were detected across three environments. Three QTLs explained 5.6–10.9 % of phenotypic variation of SP, and the qSPA10.1a and qSPB6.1a were consensus QTL.

Map-based comparison of the present and previous QTLs for yield traits The availability of the common markers enabled the comparison of the present and previous QTLs for yield traits. By common marker and the position of the QTL confidence interval, the present QTLs qSWA3.1a, qPWA5.1b q100PWA7.1a, and q100PWA7.1b were identical to the qSWA3, qPWA5, qHPWA7.1 and qHPWA7.2 (Huang et al. 2015). One QTL for seed weight identified previously (Pandey et al. 2014) was identical to the present qPWA2.1a and q100SWA2.1a.

17 Page 10 of 14

The QTL qPLA5.1a was identical to the QTL for pod maturity (Gomez Selvaraj et al. 2009). One QTL for the HMS on A6 identified by Shirasawa et al. (2013) was also possible identical to the one of present QTLs for the HMS. And seven QTLs identified by Fonceka et al. (2012) were possibly identical to the present QTLs on A7, A10, and B6 groups. A total of the present 12 QTLs were also detected using F2:3 populations of the present RIL population (Chen et al. 2016a, b). Totally, there were 27 previous QTLs were reproduced in the present study, and 65 QTLs were newly identified in the present study, and 56 of them were for pod and seed traits which directly contributed to peanut grain yield.

The comparative mapping of QTLs on the physical map Due to the good collinearity between the present map and their counterpart chromosomes of the physical map of diploid species of Arachis, an attempt was made to evaluate the physical distance of QTLs based on the e-PCR with the present primer sequence against the physical map of A. duranensis and A. ipaensis (Bertioli et al. 2016). The physical distance of the interval of flanking markers and the confidence interval of QTLs were evaluated based on the position of markers linked to the QTLs on the physical map. A total of 11 unique QTLs with flanking markers were directly assigned to the physical map (File S9). The physical distance ranged from 0.71 to 28.2 Mb with an averaged distance of 9.82 Mb. Interestingly, the uqA7–1 and uqA7–2 were located in the region of 1.29 and 0.71 Mb. For tracking the result in further study, the physical position of the markers linked to the peak position was also listed (File S9).

Mol Breeding (2017) 37: 17

structure of peanut. The good collinearity between the present map, the physical and INT map indicated the high quality of the present linkage map (Fig. 3, Figs. S2, 3 and 4). It still needs further work to make it clear whether the insufficient number of markers or the low frequency of recombination caused the separation between A7and A71, B3 to B31, B32 and B33. The good collinearity between the present and the physical map revealed an evolutionarily highly conserved genome sequences between diploid and tetraploid species of Arachis (Fig. 3). It will be significant for guiding the fine mapping and map-based cloning of QTLs by comparative mapping between the present and the physical map. The 9.82 Mb of average physical distance of the flanking markers of unique QTLs also indicated a great potential and challenge to carry out comparative mapping by using the physical map of diploid species of Arachis (Bertioli et al. 2016). The INT map was constructed using 16 populations and involved four diploid species A. duranensis, Arachis stenosperma, A. ipaensis, and Arachis magna and tetraploid species A. hypogaea (Shirasawa et al. 2013). It presented almost all the possible genetic recombination reported by all past researchers with the highest resolution of 0.7 cM/locus and the largest genetic distance of 2651 cM with 3693 loci. The high collinearity between the present map and the INT map validated the reliability of the present map and indicated a genome-wide conserved genetic make-up of tetraploid species of Arachis (Fig. S2). The frequent chromosome recombination between the present and the INT map might indicate an intrinsic characteristic of the present RIL populations. The comparative mapping significantly narrowed the collinear region on INT map in several groups and should be helpful to more precise genetic mapping of the present QTLs in future study.

Discussion

QTLs for yield traits

It is essential to employ desired materials for genetic mapping of QTLs. The present RIL showed the most abundant variations in all involved traits and thus were ideal for genetic map construction and QTL detection (Fig. 1). The present map was comprised of 26 groups and 609 loci. These valuable characteristics will be of benefit to marker transfer and the integration/ comparison among different maps in peanut (Fig. 2). The genetic separation between A7 and A71, B3 and B31, B32, and B33 gave more insight into the genetic

The present study revealed 92 QTLs for 11 traits. Among them, 15 QTLs were detected across three environments. Due to the poor knowledge of QTLs for peanut yield traits (Gomez Selvaraj et al. 2009; Varshney et al. 2010; Fonceka et al. 2012; Pandey et al. 2014; Huang et al. 2015; Chen et al. 2016a, b), up to 65 QTLs were newly identified. After two rounds of meta-analysis, the 70 overlapping QTLs were integrated into 33 consensus QTLs (Table S4); finally, a total of 29 pleiotropic QTLs were identified with significantly narrowed confidence

Mol Breeding (2017) 37: 17

interval (1.86 cM) of QTLs. The employment of metaanalysis improved the location accuracy of the QTLs and facilitated the dissection of the genetic basis of the yield traits. The present study also gave a more precise location for QTL analysis than previous researches. For the two plant architecture traits, HMS and TBN, QTLs for both traits were located on A2, A6, B1, B6, and B9 groups with moderate phenotypic variation. The previous QTLs for the HMS were located on A3, A4, A7, B4, B7, B8, and B10 groups, and the QTLs for the TBN were located on A1 and A8 with moderate phenotypic variation (Fonceka et al. 2012; Huang et al. 2015). It indicated that the genetics basis of the HMS was dependent on peanut genotypes. It was consistent to the results in other crops that the HMS was controlled by multiple QTLs (Peiffer et al. 2014; Zhang et al. 2015). It seemed that the present QTLs for the HMS were controlled mainly by genotype because of their expression across three environments. In the present study, nine pod and seed traits, the PL, PW, PL/PW, SL, SW, SL/SW, 100SW, 100PW, and SP were characterized. Most of them were newly identified, and the general correlations among them were observed (Table 1). The previously reported QTLs for these traits were located on A2, A3, A5, A7, A8, A9, A10, B2, B3, B5, B6, and B9 groups (Fonceka et al. 2012; Shirasawa et al. 2012; Jiang et al. 2014; Pandey et al. 2014; Huang et al. 2015, 2016; Chen et al. 2016a, b). In the present study, a total of 83 QTLs for pod and seed traits were identified, and 56 of them were newly identified; most of the identified QTLs could only explain moderate phenotypic variation. Besides the previously reported A2, A3, A5, A7, A10, and B6 groups, the QTLs were observed on B1, B4, B7, and B8 for the first time. It seemed that the loci controlling the yield-determining traits distributed all over the whole genome, which also indicated a complex genetic basis of the yield-determining traits. It also was practicable to clone the QTLs with moderate phenotypic variation for seed size and weight in peanut as well as in soybean (Xie et al. 2014). In the present study, the 100PW and 100SW was more than twofold for Fuchuan Dahuasheng than those for ICG6375, and the measurement of 64.3 g for 100PW and 25.5 g for 100SW of ICG6375 were almost the lowest level of these two traits in cultivated peanut (Table S1). The identified QTLs for yield traits were contributed by Fuchuan Dahuasheng with additive effect and moderate phenotypic variation explained, no QTLs for 100PW and 100SW was observed from ICG6375; the

Page 11 of 14 17

question was whether the genetic factors of these two traits could not be further dissected in ICG6375. It could be deduced that the present QTLs for pod- and seedrelated traits might be the part of the most fundamental factors contributing to the yield traits in peanut. The QTL qPLA7.1b, qPLB1.1a, qSLB1.1a, qSWB1.1a, and qSL/SWA2.1a should be candidate for MAS and mapbased cloning in further study for their more than 10 % of phenotypic variation. For the PL, only the qPLA7.1c was detected across all environments, it seemed that the PL was controlled by intensive interaction between genotype and environments. Interestingly, the qPLA5.1a was overlapped with the QTL for maturity by common marker PM45 associated with pod maturity, and the correlation was not observed between the PL and the maturity (Gomez Selvaraj et al. 2009). It was a common sense that the completion of pod expanding often symbolized the maturity of pod during normal peanut ontogenesis. Further study is still needed to clarify whether the same one QTL contributed to both the maturity and the PL in present RIL population. In B. napus, a similar trait, the silique length was also controlled by multiple QTLs (Shi et al. 2009; Yang et al. 2012; Li et al. 2014). It was expected that more QTLs of PL should be identified using different peanut genotypes. The QTLs of seed size were well characterized in soybean. The seed sizerelated QTLs were distributed over 16 chromosomes. Many of them accounted minor or moderate phenotypic variation, indicating a comprehensive genetic basis of seed size (Salas et al. 2006; Xu et al. 2011; Niu et al. 2013). Because of the close genetic relationship between soybean and peanut (Shirasawa et al. 2013), it is possible that multiple QTLs actually contributed to pod and seed size in peanut. And up to date, more than ten genes associated with grain shape and size have been characterized well in rice, and these genes acted in independent genetic pathways (Zuo and Li 2014; Wang et al. 2015b). More interestingly, the QTLs GS5, qSW/GW and qGW7 controlled seed size and affected the grain quality in rice (Peng et al. 2015; Wang et al. 2015b). It still needs further study on whether these QTLs for seed traits affected quality traits in peanut as well as GS5 and qSW/GW in rice. For the two yield-directly related traits 100PW and 100SW, all QTLs were newly identified. The QTLs q100SWA7.1c and q100PWB8.1a could be further studied for their more than 10 % of phenotypic variation. And the present 100PW and 100SW were possibly controlled mainly by genotype because the q100SWA7.1a, q100SWB6.1a, q100SWb8.1a, q100PWA7.1a, q100SWA7.1a, and q100SWB6.1b were expressed across

17 Page 12 of 14

the three environments. The 100PW and 100SW were controlled by multiple QTLs in soybean (Han et al. 2012; Kato et al. 2014; Xie et al. 2014) and B. napus (Li et al. 2014; Liu et al. 2015). In the present study, the significant correlation was observed between 100PW and 100SW, and the QTLs controlling both of them were often co-localized on A7 and B8 groups. However, the QTLs for 100PW were also located on A5, A6, and B7 groups; meanwhile, the QTLs for 100SW were located on A2 and B6 groups. The unanswered question was whether the shell contributed to the different locations of QTLs for 100PW and 100SW? The co-localization of the QTLs for yield-related traits As the results of meat analysis, the significant colocalization of QTLs was observed, and 29 pleiotropic QTLs were identified. The present pleiotropic QTLs should contribute to better understanding of yield components and the linked markers will facilitate MAS breeding in peanut. The genetic bases of yield component had been studied, and the pod- and seed-related traits were highly valued in peanut (Gomez Selvaraj et al. 2009; Shirasawa et al. 2012; Jiang et al. 2014; Pandey et al. 2014; Huang et al. 2015), soybean(Xu et al. 2011; Niu et al. 2013; Kato et al. 2014; Nemli et al. 2014; Xie et al. 2014), rapeseed (Chen et al. 2007; Li et al. 2014, 2015; Liu et al. 2015; Shi et al. 2009, 2015), and rice (Zuo and Li 2014; Wang et al. 2015b). The significant colocalizations indicated that the yield traits were dependent on each other; the pleiotropic QTLs were the important genetic factors which contributed to the yield traits in the present RIL population. It means that the locus of pleiotropic QTL containing multiple, tightly linked, traitspecific genes or the genes that affect multiple traits (Hall et al. 2006). It is often observed that the QTLs and genes acted as affecting yield exhibiting pleiotropic effects on more than one trait (Shi et al. 2009). In soybean, it was the pleiotropy contributed to the seed size and weight, the QTL X40, X53, X83–X85, and X92 simultaneously controlled two of the seed length, seed width, and 100-seed weight (Xie et al. 2014). In Brassina napus, the pleiotropic QTLs were also identified. For example, the QTLs uq. A09–1 and uq. A09–3 contributed to both the seed weight and the silique length (Li et al. 2014), the gene ARF18 simultaneously affected the seed weight and silique length (Liu et al. 2015). And in rice, the QTLs contributed to the seed size and weight were reported (Zuo and Li 2014; Wang et al. 2015b). It could

Mol Breeding (2017) 37: 17

be reasoned that the significant pleiotropic QTLs might be resulted from the artificial selection in long term of peanut breeding, because the yield-determining QTLs were all from Fuchuan Dahuasheng. The QTL cluster was also a mechanism for rapid multiple traits evolution (Yoshizawa et al. 2013). Acknowledgments The research was funded by the National Key Basic Research Special Foundation of China B973^ Project (2011CB109304), the National Natural Science Foundation of China (31000724, 31271764), and the National Science and Technology Infrastructure (NICGR2015-036). Authors’ contribution YNC and HFJ designed the experiment. YNC and XPR developed the mapping population. YLZ and WGC performed SSR genotyping. LH and WGC performed the part of phenotypic analysis. XJZ, LYY, and YL performed the agronomic traits measurements of the RIL population. YNC constructed the genetic linkage mapping, performed the QTL analysis and wrote the manuscript. SMH supported the bioinformatics analysis. DXH, YQJ, BSL, and WHW revised the manuscript. All of the authors read and approved the final manuscript. Compliance with ethical standards Conflict of interest The authors declared that they have no conflicts of interest.

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

References Arcade A, Labourdette A, Falque M, Mangin B, Chardon F, Charcosset A, Joets J (2004) BioMercator: integrating genetic maps and QTL towards discovery of candidate genes. Bioinformatics 20:2324–2326 Bertioli DJ, Cannon SB, Froenicke L et al (2016) The genome sequences of Arachis duranensis and Arachis ipaensis, the diploid ancestors of cultivated peanut. Nat Genet 48:438–446 Branch WD (2008) Inheritance of a one-seeded pod trait in peanut. J Hered 99:221–222 Burow MD, Simpson CE, Starr JL, Paterson AH (2001) Transmission genetics of chromatin from a synthetic amphidiploid to cultivated peanut (Arachis hypogaea L.): broadening the gene pool of a monophyletic polyploidy species. Genetics 159:823–837 Chardon F, Virlon B, Moreau L, Falque M, Joets J, Decousset L, Murigneux A, Charcosset A (2004) Genetic architecture of flowering time in maize as inferred from quantitative trait loci

Mol Breeding (2017) 37: 17 meta-analysis and synteny conservation with the rice genome. Genetics 168:2169–2185 Chen W, Zhang Y, Liu XP, Chen BY, Tu JX, Fu TD (2007) Detection of QTL for six yield-related traits in oilseed rape (Brassica napus) using DH and immortalized F2 populations. Theor Appl Genet 115:849–858 Chen BY, Jiang HF, Liao BS, Ren XP, Huang JQ, Lei Y, Wang SY (2008) Genetic diversity analysis of Arachis gerplasm by SSR. J Trop Subtrop Bot 16:296–303 Chen WG, Jiao YQ, Cheng LQ, Huang L, Liao BS, Tang M, Ren XP, Zhou XJ, Chen YN, Jiang HF (2016a) Quantitative trait locus analysis for pod- and kernel-related traits in the cultivated peanut (Arachis hypogaea L.). BMC Genet 17:25 Chen XP, Li HJ, Pandey MK et al (2016b) Draft genome of the peanut A-genome progenitor (Arachis duranensis) provides insights into geocarpy, oil biosynthesis, and allergens. PNAS 113:6785–6790 Chu Y, Holbrook C, Ozias-Akins P (2009) Two alleles of ahFAD2B control the high oleic acid trait in cultivated peanut. Crop Sci 49:2029–2036 Chu Y, Wu CL, Holbrook CC, Tillman BL, Person G, Ozias-Akins P (2011) Marker-assisted selection to pyramid nematode resistance and the high oleic trait in peanut. Plant Genome 4:110–117 Churchill GA, Doerge RW (1994) Empirical threshold values for quantitative trait mapping. Genetics 13:963–971 Fonceka D, Tossim HA, Ronan R, Vignes H, Faye I, Ndoye O, Moretzsohn MC, Bertioli DJ, Glaszmann JC, Brigitte C (2012) Fostered and left behind alleles in peanut: interspecific QTL mapping reveals footprints of domestication and useful natural variation for breeding. BMC Plant Biol 12:26 Francia E, Tacconi G, Crosatti C, Barabaschi D, Bulgarelli D, Dallaglio E, Valè G (2005) Marker assisted selection in crop plants. Plant Cell Tissue Org 82:317–342 Garcia GM, Stalker HT, Schroeder E, Lyerly JH, Kochert G (2005) A RAPD-based linkage map of peanut based on a backcross population between the two diploid species Arachis stenosperma and A. cardenasii. Peanut Sci 32:1–8 Gautami B, Foncéka D, Pandey MK, Moretzsohn MC, Sujay V, Qin H, Hong Y, Faye I, Chen X, BhanuPrakash A (2012) An international reference consensus genetic map with 897 marker loci based on 11 mapping populations for tetraploid groundnut (Arachis hypogaea L.). PLoS One 7:e41213 Goffinet B, Gerber S (2000) Quantitative trait loci: a meta-analysis. Genetics155:463-473 Gomez Selvaraj M, Narayana M, Schubert AM, Ayers JL, Baring MR, Burow MD (2009) Identification of QTLs for pod and kernel traits in cultivated peanut by bulked segregant analysis. Electron J Biotech 12:3–4 Hall MC, Basten CJ, Willis JH (2006) Pleiotropic quantitative trait loci contribute to population divergence in traits associated with life-history variation in Mimulus guttatus. Genetics 172: 1829–1844 Halward T, Stalker H, Kochert G (1993) Development of an RFLP linkage map in diploid peanut species. Theor Appl Genet 87: 379–384 Han YP, Li DM, Zhu D, Li HY, Li XP, Teng WL, Li WB (2012) QTL analysis of soybean seed weight across multi-genetic backgrounds and environments. Theor Appl Genet 125:671– 683

Page 13 of 14 17 Holbrook CC, Stalker HT (2003) In peanut breeding and genetic resources. In: Janick J (ed) Plant breeding reviews, vol Volume 22. John Willey and Sons, New York, pp. 297–356 Hong YB, Chen XP, Liang XQ, Liu HY, Zhou GY, Li SX, Wen SJ, Holbrook CC, Guo BZ (2010) A SSR-based composite genetic linkage map for the cultivated peanut (Arachis hypogaea L.) genome. BMC Plant Biol 10:17 Huang L, He HY, Chen WG, Ren XP, Chen YN, Zhou XJ, Xia YL, Wang XL, Jiang XG, Liao BS, Jiang HF (2015) Quantitative trait locus analysis of agronomic and qualityrelated traits in cultivated peanut (Arachis hypogaea L.). Theor Appl Genet 128:1103–1115 Huang L, Wu B, Zhao JJ, Li HT, Chen WG, Zheng YL, Ren XP, Chen YN, Zhou XJ, Lei Y, Liao BS, Jiang HF (2016) Characterization and transferable utility of microsatellite markers in the wild and cultivated Arachis species. PLoS One 11:e0156633 Janila P, Pandey MK, Shasidhar Y et al (2016) Molecular breeding for introgression of fatty acid desaturase mutant alleles (ahFAD2A and ahFAD2B) enhances oil quality in high and low oil containing peanut genotypes. Plant Sci 242:203–213 Jiang HF, Huang L, Ren XP, Chen YN, Zhou XJ, Xia YL, Huang JQ, Lei Y, Yan LY, Wan LY, Liao BS (2014) Diversity characterization and association analysis of agronomic traits in a Chinese peanut (Arachis hypogaea L.) mini-core collection. J Integr Plant Biol 56:159–169 Kato S, Sayama T, Fujii K, Yumoto S, Kono Y, Hwang TY, Kikuchi A, Takada Y, Tanaka Y, Shiraiwa T et al (2014) A major and stable QTL associated with seed weight in soybean across multiple environments and genetic backgrounds. Theor Appl Genet 127:1365–1374 Knapp SJ (1998) Marker-assisted selection as a strategy for increasing the probability of selecting superior genotypes. Crop Sci 38:1164–1174 Kosambi D (1944) The estimation of map distances from recombination values. Ann Eugenics 12:172–175 Li N, Shi JQ, Wang XF, Liu GH, Wang HZ (2014) A combined linkage and regional association mapping validation and fine mapping of two major pleiotropic QTLs for seed weight and silique length in rapeseed (Brassica napus L.). BMC Plant Biol 14:114 Li SP, Chen L, Zhang LW, Li X, Liu Y, Wu ZK, Dong FM, Wan LL, Liu KD, Hong DF, Yang GS (2015) BnaC9.SMG7b functions as a positive regulator of number of seeds per silique in rapeseed (Brassica napus L.) by regulating the formation of functional female gametophytes. Plant Physiol 169:2744–2760 Liu J, Hua W, Hu ZY, Yang HL, Zhang L, Li RJ, Deng LB, Sun XC, Wang XF, Wang HZ (2015) Natural variation in ARF18 gene simultaneously affects seed weight and silique length in polyploid rapeseed. Proc Natl Acad Sci 112:E5123–E5132 Meyer RS, Purugganan MD (2013) Evolution of crop species: genetics of domestication and diversification. Nat Rev Genet 14:840–852 Moles AT, Ackerly DD, Webb CO, Tweddle JC, Dickie JB, Westoby M (2005) A brief history of seed size. Science 307:576–580 Nemli S, Asciogul TK, Kaya HB, Kahraman A, Eşiyok D, Tanyolac B (2014) Association mapping for five agronomic traits in the common bean (L.). Journal of the Science of Food and Agriculture 94 (15):3141–3151

17 Page 14 of 14 Niu Y, Xu Y, Liu XF, Yang SX, Wei SP, Xie FT, Zhang YM (2013) Association mapping for seed size and shape traits in soybean cultivars. Mol Breeding 31:785–794 Pandey MK, Upadhyaya HD, Rathore A et al (2014) Genomewide association studies for 50 agronomic traits in peanut using the ‘reference set’ comprising 300 genotypes from 48 countries of the semi-arid tropics of the world. PLoS One 9:e105228 Peiffer JA, Romay MC, Gore MA et al (2014) The genetic architecture of maize height. Genetics 196:1337–1356 Peng Y, Hu YY, Mao BG, Xiang HT, Shao Y, Pan YL, Sheng XB, Li YK, Ni XM (2015) Genetic analysis for rice grain quality traits in the YVB stable variant line using RAD-seq. Mol Gen Genomics 291:297–307 Qin HD, Feng SP, Chen C, Guo YF, Knapp S, Culbreath A, He GH, Wang ML, Zhang XY, Holbrook CC (2012) An integrated genetic linkage map of cultivated peanut (Arachis hypogaea L.) constructed from two RIL populations. Theor Appl Genet 124:653–664 Salas P, Oyarzo-Llaipen JC, Wang D, Chase K, Mansur L (2006) Genetic mapping of seed shape in three populations of recombinant inbred lines of soybean (Glycine max L. Merr.). Theor Appl Genet 113:1459–1466 Shi J, Zhan J, Yang Y, Ye J, Huang S, Li R, Wang X, Liu G, Wang H (2015) Linkage and regional association analysis reveal two new tightly-linked major-QTLs for pod number and seed number per pod in rapeseed (Brassica napus L.). Scientific Reports 5:14481 Shi JQ, Li RY, Qiu D, Jiang CC, Long Y, Morgan C, Bancroft I, Zhao JY, Meng JL (2009) Unraveling the complex trait of crop yield with quantitative trait loci mapping in Brassica napus. Genetics 182:851–861 Shirasawa K, Koilkonda P, Aoki K, Hirakawa H, Tabata S, Watanabe M, Hasegawa M, Kiyoshima H, Suzuki S, Kuwata C (2012) In silico polymorphism analysis for the development of simple sequence repeat and transposon markers and construction of linkage map in cultivated peanut. BMC Plant Biol 12:80 Shirasawa K, Bertioli DJ, Varshney RK, Moretzsohn MC, LealBertioli SC, Thudi M, Pandey MK, Rami J-F, Foncéka D, Gowda MV (2013) Integrated consensus map of cultivated peanut and wild relatives reveals structures of the A and B genomes of Arachis and divergence of the legume genomes. DNA Res 20:173–184 Sujay V, Gowda M, Pandey M, Bhat R, Khedikar Y, Nadaf H, Gautami B, Sarvamangala C, Lingaraju S, Radhakrishan T (2012) Quantitative trait locus analysis and construction of consensus genetic map for foliar disease resistance based on two recombinant inbred line populations in cultivated groundnut (Arachis hypogaea L.). Mol Breeding 30:773–788 Udall JA, Quijada PA, Lambert B, Osborn TC (2006) Quantitative trait analysis of seed yield and other complex traits in hybrid spring rapeseed (Brassica napus L.): 2. Identification of alleles from unadapted germplasm. Theor Appl Genet 113: 597–609 Van Ooijen J, Voorrips R (2006) Joinmap 4.0. Software for the calculation of genetic linkage map. Kyazma BV, Wageningen Varshney RK, Glaszmann JC, Leung H, Ribaut JM (2010) More genomic resources for less-studied crops. Trends Biotechnol 28:452–460

Mol Breeding (2017) 37: 17 Varshney RK, Mohan SM, Gaur PM et al (2013) Achievements and prospects of genomics-assisted breeding in three legume crops of the semi-arid tropics. Biotechnol Adv 31:1120–1134 Varshney RK, Pandey MK, Janila P, Nigam SN, Sudini H, Gowda MV, Sriswathi M, Radhakrishnan T, Manohar SS, Nagesh P (2014) Marker-assisted introgression of a QTL region to improve rust resistance in three elite and popular varieties of peanut (Arachis hypogaea L.). Theor Appl Genet 127: 1771–1781 Voorrips RE (2002) MapChart: software for the graphical presentation of linkage map and QTLs. J Hered 93:77–78 Wang H, Penmetsa RV, Yuan M, Gong LM, Zhao YL, Guo BZ, Farmer AD, Rosen BD, Gao JL, Isobe S (2012a) Development and characterization of BAC-end sequence derived SSRs, and their incorporation into a new higher density genetic map for cultivated peanut (Arachis hypogaea L.). BMC Plant Biol 12:–10 Wang S, Basten CJ, Zeng Z (2012b) Windows QTL Cartographer 2.5. Department of statistics. North Carolina State University, Raleigh, NC Wang ML, Khera P, Pandey MK et al (2015a) Genetic mapping of QTLs controlling fatty acids provided insights into the genetic control of fatty acid synthesis pathway in peanut (Arachis hypogaea L.). PLoS One 10:e0119454 Wang SK, Li S, Liu Q, Zhang JQ, Wang SS, Wang Y, Chen XB, Zhang Y, Gao CX (2015b) The OsSPL16-GW7 regulatory module determines grain shape and simultaneously improves rice yield and grain quality. Nat Genet 47:949–954 Wu XL, Blake S, Sleper DA, Shannon JG, Cregan P, Nguyen HT (2009) QTL, additive and epistatic effects for SCN resistance in PI 437654. Theor Applied Genet 118:1093–1105 Xie FT, Niu Y, Zhang J, Bu SH, Zhang HZ, Geng QC, Feng JY, Zhang YM (2014) Fine mapping of quantitative trait loci for seed size traits in soybean. Mol Breeding 34:2165–2178 Xu Y, Li HN, Li GJ, Wang X, Cheng LG, Zhang YM (2011) Mapping quantitative trait loci for seed size traits in soybean (Glycine max L. Merr.). Theor Appl Genet 122:581–594 Yang P, Shu C, Chen L, Xu JS, Wu JS, Liu KD (2012) Identification of a major QTL for silique length and seed weight in oilseed rape (Brassica napus L.). Theor Appl Genet 125:285–296 Yoshizawa M, O'Quin KE, Jeffery WR (2013) QTL clustering as a mechanism for rapid multi-trait evolution. Commun Integr Biol 6:e24548 Zhang J, Song Q, Cregan PB, Nelson RL, Wang X, Wu J, Jiang GL (2015) Genome-wide association study for flowering time, maturity dates and plant height in early maturing soybean (Glycine max) germplasm. BMC Genomics 16:217 Zhou XJ, Xia YL, Ren XP, Chen YN, Huang L, Huang SM, Liao BS, Lei Y, Yan LY, Jiang HF (2014) Construction of a SNPbased genetic linkage map in cultivated peanut based on large scale marker development using next-generation double-digest restriction-site-associated DNA sequencing (ddRADseq). BMC Genomics 15:351 Zhou XJ, Dong Y, Zhao JJ, Huang L, Ren XP, Chen YN, Huang SM, Liao BS, Lei Y, Yan LY, Jiang HF (2016) Genomic survey sequencing for development and validation of single-locus SSR markers in peanut (Arachis hypogaea L.). BMC Genomics 17:420 Zuo JR, Li JY (2014) Molecular genetic dissection of quantitative trait loci regulating rice grain size. Annu Rev Genet 48:99–118