Identification and validation of QTLs for kernel number per spike and spike length in two founder genotypes of wheat
BMC Plant Biology volume 22, Article number: 146 (2022)
Kernel number per spike (KNS) and spike length (SL) are important spike-related traits in wheat variety improvement. Discovering genetic loci controlling these traits is necessary to elucidate the genetic basis of wheat yield traits and is very important for marker-assisted selection breeding.
In the present study, we used a recombinant inbred line population with 248 lines derived from the two founder genotypes of wheat, Bima4 and BainongAK58, to construct a high-density genetic map using wheat 55 K genotyping assay. The final genetic linkage map consists of 2356 bin markers (14,812 SNPs) representing all 21 wheat chromosomes, and the entire map spanned 4141.24 cM. A total of 7 and 18 QTLs were identified for KNS and SL, respectively, and they were distributed on 11 chromosomes. The allele effects of the flanking markers for 12 stable QTLs, including four QTLs for KNS and eight QTLs for SL, were estimated based on phenotyping data collected from 15 environments in a diverse wheat panel including 384 elite cultivars and breeding lines. The positive alleles at seven loci, namely, QKns.his-7D2–1, QKns.his-7D2–2, QSl.his-4A-1, QSl.his-5D1, QSl.his-4D2–2, QSl.his-5B and QSl.his-5A-2, significantly increased KNS or SL in the diverse panel, suggesting they are more universal in their effects and are valuable for gene pyramiding in breeding programs. The transmission of Bima4 allele indicated that the favorite alleles at five loci (QKns.his-7D2–1, QSl.his-5A-2, QSl.his-2D1–1, QSl.his-3A-2 and QSl.his-3B) showed a relatively high frequency or an upward trend following the continuity of generations, suggesting that they underwent rigorous selection during breeding. At two loci (QKns.his-7D2–1 and QSl.his-5A-2) that the positive effects of the Bima4 alleles have been validated in the diverse panel, two and one kompetitive allele-specific PCR (KASP) markers were further developed, respectively, and they are valuable for marker-assisted selection breeding.
Important chromosome regions controlling KNS and SL were identified in the founder parents. Our results are useful for knowing the molecular mechanisms of founder parents and future molecular breeding in wheat.
Wheat (Triticum aestivum L.) is a major cereal crop worldwide. The current yield trend in wheat is insufficient to meet the future demand of a growing world population, and wheat yield and total production must be increased continuously. The formation of wheat yield is a complex trait and is affected by three yield components, e.g., spike number per unit area, kernel weight, and kernel number per spike (KNS). Among them, kernel weight and KNS are closely related to spike morphology, which is primarily determined by spike length (SL), spikelet density, and fertile floret number. Previous studies showed that increasing KNS is an effective approach for wheat yield improvements compared to kernel weight [1, 2], and increasing SL without modification of the spikelet density can increase KNS and subsequently raise yield capacity . A positive correlation between SL and yield was also validated in some previous studies . Besides, the long spike is often associated with reduced severity of Fusarium head blight in wheat . Therefore, discovering genetic loci controlling KNS and SL is necessary to elucidate the genetic basis of wheat yield traits and is very important for marker-assisted selection (MAS) breeding.
Like other spike-related traits, KNS and SL are controlled by multiple genes and affected by environments. Quantitative trait loci (QTL) analysis using different genetic populations and diverse wheat panels provides an effective method to study the genes governing these traits. To date, numerous QTLs associated with KNS have been mapped on nearly all the 21 chromosomes in wheat, such as two QTLs identified on 2D and 4A , eight QTLs identified on 1A, 1B, 2B, 2D, 3B, 4B, 6A and 7B , one QTL identified on 3D , one QTL identified on 4A , four QTLs identified on 2A, 4B and 7A , six QTLs identified on 1D, 2A, 2D, 3A, 4D and 6D , three QTLs identified on 1D, 4D and 6B , 12 QTLs identified on 1A, 2D, 3B, 4A, 4B, 5A, 5B, 7A and 7B , and one QTL identified on 7A . Likewise, many previous studies have proven that almost all the 21 wheat chromosomes harbored factors affecting SL [3, 10, 15,16,17,18,19,20,21]. Yao et al.  reported that approximately 350 QTLs of SL have been identified currently, and some of them with relatively large effects were distributed on chromosomes 2D, 3A, 4A, 4B, 5A, 6A, 6B, 7A, 7B, and 7D. Briefly, because of the complexity of the wheat genome, although many QTL for KNS and SL have been reported, common QTLs across different mapping populations are limited, and few of them are used in practical wheat breeding.
Founder parents have played particularly crucial roles in the improvement of wheat worldwide. Many QTLs or chromosomal regions associated with important traits have been found in founder genotypes in wheat [23,24,25,26]. However, the knowledge of the molecular mechanisms for the formation of founder parents remains unclear. In China, Bima4 is one of the founder parents that played important roles in wheat breeding, used widely in the Yellow and Huai River Facultative Winter Wheat Region between 1950 and 1970 . It was obtained from the cross between another founder parent Mazhamai and Quality from the United States. More than 70 improved cultivars were developed from Bima4, and some of them such as Shijiazhuang54, Jinan2, Beijing8 and Taishan1, had annual maximum acreages over 667,000 ha and were grown for at least 12 years. Similarly, BainongAK58 is a famous cultivar released in 2003 by the Henan Institute of Science and Technology, and its maximum acreage was over 13,333,333 ha. It was also widely utilized as a crossing parent in wheat breeding, from which more than fifty improved cultivars were developed.
In the present study, we used a recombinant inbred line (RIL) population with 248 lines derived from the two founder genotypes of wheat, Bima4 and BainongAK58. The QTL analysis was conducted with a high-density genetic map by using the developed wheat 55 K genotyping assay to identify QTLs responsible for KNS and SL. These QTLs detected were further validated in a diverse wheat panel. Furthermore, we analyzed the transmission of Bima4 alleles to its derivative descendants, and two and one KASP markers for two important loci, QKns.his-7D2–1 and QSl.his-5A-2, were developed. This study is useful for knowing the molecular mechanisms of founder parents and future molecular breeding in wheat.
Linkage map construction in the RIL population
Out of 53,063 SNPs in the 55 k Infinium chip, 16,628 SNPs were polymorphic between the two parents and among the RIL population. These 16,628 markers were divided into 2488 bins. Only one marker was chosen to represent each bin for the genetic map construction. The final genetic linkage map consists of 2356 bin markers (14,812 SNPs) representing all 21 wheat chromosomes. Of them, 1147 bins include only one SNP marker, and the remaining comprises two or more SNP markers.
The 2356 bin markers were mapped on 28 linkage maps (Table 1 and Table S1). Each of the chromosomes 1A, 1D, 2D, 3D, 4D, 5D, and 7D was integrated by two linkage groups. The entire map spanned 4141.24 cM with six gaps (> 30 cM) distributed on chromosomes 2D, 6A, 4D, 5B, and 7D. The mean of genetic distance among adjacent bin markers across all chromosomes was 1.76 cM and varied among 28 linkage groups from 0.71 (1A2) to 5.98 (7D1). The bin markers mapped on the A genome (37.9%) were more than those on the B (34.2%) and D (28.0%) genome. Similarly, most of the mapped markers including bin and redundant markers were distributed on A (43.0%) and B genome (36.9%), and only 20.0% of the markers were mapped on D genome. The number of bin markers on 21 chromosomes ranged from 51 on 1D to 177 on 7D, however the number of the mapped markers ranged from 168 on 4D to 1479 on 2A.
Phenotypic analysis for KNS and SL in the RIL population
These two traits for the RIL populations and the two parents in the four environments are shown in Table S2. The SL and KNS showed inconsistency between the parental lines over environments, indicating strongly affected by the environment. In the RIL populations, the KNS and SL showed normal distributions in all the environments, suggesting the polygenic inheritance of these traits (Fig. 1). The transgressive inheritance was found in certain lines for SL and KNS (Fig. 2). The two traits showed strong correlations with each other in all environments. The correlation coefficients ranged from 0.86 to 0.96 for SL and from 0.50 to 0.86 for KNS. The SL had a strong positive correlation with KNS at 0.24 (P < 0.0001) (Table 2). The SL and KNS showed high broad-sense heritability at 0.95 and 0.85, respectively.
QTL detection for KNS and SL in the RIL population
A total of seven QTLs were detected for KNS on chromosomes 3A, 3D, 4A, 5A, and 7D (Table 3 and Fig. 3). A major stable QTL, Qkns.his-4A, was detected in all four environments and the average value and explained 9.78–24.24% of the phenotypic variance. Qkns.his-5A-2 was identified in three environments and the average value and explained 3.72–7.01% of phenotypic variation. The positive alleles of Qkns.his-4A and Qkns.his-5A-2 were contributed by Bima4. The two QTLs, Qkns.his-7D-1 and Qkns.his-7D-2, were detected in one environment and the average value, and the positive alleles were contributed by Bima4 and BainongAK58, respectively. The remaining three QTLs, Qkns.his-3A, Qkns.his-3D and Qkns.his-5A-1, were detected in a single environment, and they explain 5.29, 4.22, and 4.01% of the phenotypic variance, respectively.
A total of 18 QTLs were detected for SL on chromosomes 2D, 3A, 3B, 3D, 4A, 4D, 5A, 5B, 5D, 7B, and 7D with phenotypic variations ranging from 2.04 to 22.31% (Table 3 and Fig. 3). Among them, eight QTLs were detected in a single environment, explaining 2.93–4.85% of the phenotypic variance. Two stable QTLs, Qsl.his-2D1–1 and Qsl.his-5A-2, were detected in all four environments and the average values, and the positive alleles were from Bima4. Of these, the major QTL, Qsl.his-2D1–1, explained 11.03–22.31% of the phenotypic variance. Two QTLs, Qsl.his-4A-1 and Qsl.his-5D1, were identified in three environments and the average values. They accounted for 3.43–8.21% of the phenotypic variance and the positive alleles were contributed by BainongAK58. The two QTLs, Qsl.his-3A-2 and Qsl.his-3B, were detected in three environments and the average values, and they explained 2.18–3.65% of phenotypic variation and the positive alleles were from Bima4.
Validation of the QTL effects in the diverse wheat panel
The allele effects of the flanking markers for 12 stable QTLs, including four QTLs for KNS and eight QTLs for SL, were estimated based on phenotyping data in the diverse wheat panel. T-test analyses were used to compare the two different allele groups in the same locus (P < 0.05). For each QTL investigated, the QTL-associated SNP markers for which differences of phenotypic values showed significance in most environments in the diverse wheat panel were analyzed (Table 4). For KNS, three QTL-associated SNP markers of QKns.his-7D2–1 were analyzed. The marker AX-110196726 showed significance in 14 environments in the diverse wheat panel, while the other two markers (AX-109475040 and AX-111094913) showed a significant difference of KNS only in six environments, indicating that the former was closer to this QTL compared to the latter. The positive allele of AX-110196726 contributed by Bima4 had a more KNS than BainongAK58 allele in the diverse wheat panel. Similarly, two QTL-associated SNP markers (AX-108912162 and AX-111577597) of QKns.his-7D2–2 were analyzed. The positive alleles contributed by BainongAK58 showed significant across all 15 environments more KNS at 5.17 and 4.85 than the Bima4 alleles at these two loci in the diverse wheat panel, respectively. In addition, among six QTL-associated SNP markers of QKns.his-5A-2 where the favourable alleles were from Bima4, four showed significant differences of KNS in 15 environments, respectively. At the same time, the effects on KNS were inverse among them, i.e., the Bima4 alleles showed negative effects at AX-109980237 and AX-110121838, but positive effects at AX-111102726 and AX-109876198. In the same way, of four QTL-associated SNP markers of QKns.his-4A, two markers (AX-111508583 and AX-109332913) showed significant differences of KNS only in 6 environments respectively. At the same time, the effects on KNS were inverse between them.
For SL, six and two QTL-associated SNP markers were analyzed for QSl.his-4A-1 and QSl.his-5D1, respectively. Only AX-108955453 and AX-110409786 had significant differences of SL in 12 environments in the diverse wheat panel separately, and the positive alleles contributed by BainongAK58 increased 0.52 cm and 0.33 cm SL across these two loci, respectively. For QSl.his-4D2–2 and QSl.his-5B, two and four QTL-associated SNP markers were analyzed, respectively. Significant differences of SL were found at AX-109924587 and AX-108886889 in 15 and 10 environments in the diverse wheat panel, respectively, and the positive alleles from BainongAK58 increased 0.46 cm and 0.31 cm SL across these two loci, respectively. Four QTL-associated SNP markers of QSl.his-5A-2 were analyzed. Two markers (AX-109622137 and AX-110199675) showed significant differences of SL in 14 and 13 environments in the diverse wheat panel, respectively, and the positive alleles from Bima4 had a higher SL than BainongAK58 alleles. In addition, four QTL-associated SNP markers of QSl.his-3A-2 and QSl.his-3B were analyzed, respectively. Only AX-111618763 and AX-110931375 showed significant differences of SL in 8 and 12 environments separately, while the positive alleles obtained from Bima4 across the two loci were unfavorable in the diverse wheat panel. Likewise, of five QTL-associated SNP markers of QSl.his-2D1–1 analyzed, two had significant differences of SL in 15 environments, whereas the positive alleles from Bima4 across the two loci were unfavorable in the diverse wheat panel.
Tracking of Bima4 allele in its derivatives and development of KASP markers
Of the 12 stable QTLs analyzed above, the positive alleles for 7 QTLs including 3 QTLs of KNS and 4 QTLs of SL were contributed by Bima4 in the RIL population. The transmission of Bima4 alleles in the QTL-associated SNP markers was determined using its 70 descendants. Among the seven QTLs, the transmission of Bima4 alleles at five loci (QKns.his-7D2–1, QSl.his-5A-2, QSl.his-2D1–1, QSl.his-3A-2, and QSl.his-3B) to its derivative descendants showed a relatively high frequency or an upward trend following the continuity of generations. For example, the Bima4 allele at the QTL-associated SNP marker (AX-110196726) of QKns.his-7D2–1 showed an upward trend following the continuity of generations (Fig. 4A). The Bima4 alleles at two QTL-associated SNP markers (AX-109622137 and AX-110199675) of QSl.his-5A-2 also presented a relatively high frequency or an upward trend following the continuity of generations (Fig. 4A). For these two QTLs, the positive effects of Bima4 alleles have been validated in the diverse wheat panel. Furthermore, two flanking SNP (AX-110945813 and AX-111490337) of QKns.his-7D2–1, which were located in the same bin with AX-110196726 and the physical distances between these two markers and AX-110196726 were only 0.11 Mb and 0.20 Mb, were successfully converted to kompetitive allele-specific PCR (KASP) markers (Table 5). Likewise, a KASP marker was developed from the flanking SNP (AX-108964722) of QSl.his-5A-2, which were located in the same bin with AX-109622137 and the interval between them was only 1.57 Mb. Similarly, for the three SL QTLs (QSl.his-2D1–1, QSl.his-3A-2, and QSl.his-3B), the Bima4 alleles at four markers (AX-110332825, AX-108836084, AX-111618763, and AX-110931375) showed a high frequency in its four derivate generations, respectively (Fig. 4B). In addition, for QKns.his-4A, the Bima4 allele showed a relatively high frequency (100–80.0%) in its four derivate generations at Ax-109332913, but a low frequency (30.0%) across all derivatives at the other marker AX-111508583 (Fig. 4C). A similar result could be observed for QKns.his-5A-2, e.g., the Bima4 alleles showed a high frequency or an upward trend following the continuity of generations at two markers (AX-109980237 and AX-110121838) but had a low frequency (50.0%) across all derivatives at another two markers (AX-111102726 and AX-109876198) (Fig. 4D).
Comparison of the QTLs identified for KNS and SL with previous studies
Generally, the major QTLs consistent over environments may play a key role in modulating the agronomic traits of wheat cultivars and have great value for MAS in breeding programs. Based on genetic marker sequence flanking for KNS and SL QTLs and the genome sequence from Chinese Spring wheat (IWGSC V1.0) (http://www.wheatgenome.org/), physical positions of these stable QTLs detected in our study were compared with those reported previously. In the present study, four major QTLs, Qkns.his-4A, Qkns.his-5A-2, Qkns.his-7D2–1 and Qkns.his-7D2–2, were identified for KNS. Of these, the locus Qkns.his-4A in the interval 122.03–122.24 cM on 4A was identified in four environments and the average value, and it was located in the interval 642.37–672.88 Mb. Using a 660 K wheat SNP array, Cui et al.  identified a major stable QTL, qKnps-4A, for KNS in the interval 680.40–683.64 Mb. Gao et al.  detected a QTL, QKNS.caas-4AL, in the interval 626.32–660.99 Mb using a 90 K wheat SNP array. Kirigwi et al.  identified two simple sequence repeat (SSR) loci, Xwmc89 and Xwmc420, related to KNS at positions 515.85 Mb and 538.22 Mb, respectively. There were also some other reported QTLs for KNS on 4A in previous studies [6, 13, 30]. Nevertheless, further research is needed to identify whether these genes are identical. Qkns.his-5A-2, mapped in three environments and the average value in the interval 106.50–125.16 cM on 5A in the present study, was located in the interval 549.34–572.81 Mb, whereas only a minor QTL for KNS reported by Wang et al.  in a single environment in marker interval Xgwm126-Xgwm291 on 5A positioned in the interval 671.39–698.19 Mb, indicating that Qkns.his-5A-2 is likely to be a new KNS QTL. Likewise, Qkns.his-7D2–1 and Qkns.his-7D2–2, mapped in the 2018XI environment and the average values in the interval 3.00–6.00 cM and 159 cM on 7D in the present study separately, were located in the interval 4.39–7.45 Mb and 526.36–530.75 Mb, respectively. There were also some other reported QTLs for KNS on 7D using SSR or RFLP markers [6, 31], but these markers could not be obtained or precisely located in the reference genome. So, we cannot determine whether the loci were nearby or identical with our results or not.
Of 18 QTLs for SL identified in the present study, 10 were detected in at least two environments. Of these, Qsl.his-2D1–1 was identified in all four environments and the average value. The locus Qsl.his-2D1–1 at the interval 29.00–37.00 cM explained 11.03–22.31% of the phenotypic variance and was located on 2D in the interval 13.25–36.89 Mb. Wu et al.  identified an SL-associated gene, QSpl.nau-2D, near position 23.02 Mb. Chai et al.  identified two QTLs (QPht/Sl.cau-2D.1 and QPht/Sl.cau-2D.2) with pleiotropic effects on plant height and SL. QPht/Sl.cau-2D.1 is a novel QTL located between SNP makers BS00022234_51 and BobWhite_rep_c63957_1472 near position 20.77 Mb, whereas QPht/Sl.cau-2D.2 was located on the same genetic interval of Rht8. In addition, Sourdille et al. , Kumar et al.  and Suenaga et al.  identified one SSR locus, Xgwm261, associated with SL on 2D at position 19.6 Mb. The marker Xgwm261 is linked to the dwarf gene Rht8. Some previous studies [34,35,36] indicated that Rht8 does not affect SL, but contrasting with other recent studies showing that Rht8 introgression decreased SL with constant spikelet number . Therefore, further research is needed to identify whether Qsl.his-2D1–1 is nearby or identical with Rht8.
In the present study, Qsl.his-5A-2 was identified in all four environments and the average value. It was at the interval 94.00–96.00 cM explained 4.65–8.38% of the phenotypic variance and was located on 5A in the interval 506.65–524.73 Mb. Fan et al.  identified an SL-associated gene, qSl-5A.3, in the interval 478.65–541.29 Mb. Kumar et al.  identified an SSR locus, Xgwm186, related to SL on 5A at position 471.71 Mb. Cui et al.  detected two QTLs related to SL on 5A at positions 444.92 and 682.71 Mb, respectively. In addition, Liu et al.  found an SNP maker, IAAV8258, related to SL on 5A at position 572.84 Mb. In our study, the locus Qsl.his-5D1, which accounted for 3.97–8.21% of the phenotypic variance and was identified in three environments and the average value at the interval 151–156 cM, was located on 5D1 in the interval 446.80–475.31 Mb. At a similar location to Qsl.his-5D1, marker Xgwm182 (439.22 Mb) on 5D affecting SL was reported by Kumar et al. . Deng et al.  also reported a QTL, QSl.sdau-5D, linked to SL in marker interval Xbarc1097-Xcfd8 on 5D positioned in the interval 287.41–396.41 Mb.
QTL effects in the diverse wheat panel
In this study, we conducted the allelic analysis based on phenotyping data collected from 15 environments in the diverse wheat panel. As the result showed, the positive alleles of QKns.his-7D2–1, QKns.his-7D2–2, QSl.his-4A-1, QSl.his-5D1, QSl.his-4D2–2, QSl.his-5B, and QSl.his-5A-2 significantly increased KNS or SL in the diverse panel, suggesting that they are more universal in their effects. These important loci were very beneficial to pyramid breeding in wheat. On the other hand, the positive alleles of QSl.his-2D1–1, QSl.his-3A-2, and QSl.his-3B in the RIL population were unfavorable in the diverse wheat panel, indicating they may be population-specific QTL. In addition, for QKns.his-5A-2 where the favorable allele was obtained from Bima4 in the RIL population, the Bima4 alleles showed negative effects at two loci (AX-109980237 and AX-110121838), but positive effects at another two loci (AX-111102726 and AX-109876198) in the diverse wheat panel. QKns.his-5A-2 was mapped at the interval 106.00–125.00 cM and the physical distance between the marker AX-110121838 (549336395) and AX-111102726 (572237027) reached 22.90 Mb. These results indicated that there may be a great distance between these flanking markers and the peak markers for QKns.his-5A-2.
Transmission of Bima4 alleles to its derivative descendants
Bima4 possesses many superior agronomic traits, especially high resistance to stripe rust, and it has played a crucial role in Chinese wheat breeding and production. In this study, the transmission of Bima4 alleles which showed positive effects in the RIL population at five loci (QKns.his-7D2–1, QSl.his-5A-2, QSl.his-2D1–1, QSl.his-3A-2, and QSl.his-3B) to its derivative descendants showed a relatively high frequency or an upward trend following the continuity of generations, suggesting that they underwent rigorous selection during breeding. These important loci in Bima4 had a great effect on the improvement of wheat breeding and should be studied intensively. Our results also accorded with previous reports by Guo et al. , Li et al. , Russell et al. , Pestsova and Röder  and Sjakste et al. , who found that the alleles selected preferentially in progeny were associated with advantageous traits. More importantly, the positive effects of the Bima4 alleles at these two loci QKns.his-7D2–1 and QSl.his-5A-2 have been validated in the diverse panel. We further developed two and one KASP markers for these two loci, which are valuable for MAS breeding. Similarly, a few KASP markers were developed in some studies for yield-related traits such as thousand kernel weight , grain length , productive tiller and fertile spikelet numbers , and plant height, SL, and total spikelet number per spike . Compared with conventional molecular markers such as SSR, these KASP markers are more accurate and high-throughput, which can greatly improve the speed and efficiency of genomic selection for MAS breeding [44, 45].
A high-density genetic map, consisting 2356 bin markers (14,812 SNPs) and spanning 4141.24 cM, was constructed using the wheat 55 K genotyping assay in the RIL population with 248 lines derived from the two founder genotypes of wheat, Bima4 and BainongAK58. A total of seven and 18 QTLs were identified for KNS and SL, respectively, and they were distributed on 11 chromosomes. The allele effects of the flanking markers for 12 stable QTLs including four QTLs for KNS and eight QTLs for SL were estimated based on phenotyping data collected from 15 environments in a diverse wheat panel including 384 elite cultivars and breeding lines. The positive alleles at seven loci significantly increased KNS or SL in the diverse panel, suggesting that they are more universal in their effects and are valuable for gene pyramiding in breeding programs. The transmission of the Bima4 alleles indicated that the favorite alleles at five loci showed a relatively high frequency or an upward trend following the continuity of generations, suggesting that they underwent rigorous selection during breeding. The positive effects of the Bima4 alleles at two loci QKns.his-7D2–1 and QSl.his-5A-2 have been validated in the diverse panel, and two and one KASP markers were developed for these two loci. Our results are useful for knowing the molecular mechanisms of founder parents and future molecular breeding in wheat.
The QTL mapping population containing 248 RILs (F7) were derived from the F2 population of the cross BainongAK58 × Bima4 by the single seed descent method. Bima4 is both an important founder genotype and a widely grown cultivar with high yield potential and wide environmental adaptability. BainongAK58 has many important traits such as lodging resistance, disease resistance, and yield potential. A diverse wheat panel containing 384 elite cultivars and breeding lines was used for QTL validation in this study, and detailed information was described in Li et al. . Seventy cultivars derived from Bima4 were also included (Table S3), and there are 12, 35, 18, and 5 accessions in the first, second, third, and fourth generations of the derivatives, respectively. Seeds of all accessions were provided by the National Crop Gene Bank, Chinese Academy of Agricultural Sciences, Beijing.
Field trials and data analysis
Field experiments for the RIL population were performed at Xinxiang (117.17°E, 40.69°N) in 2017, 2018 and 2019 (2017XI, 2018XI and 2019XI) and Huixian (116.41°E, 39.91°N) in 2018 (2018HU) in Henan province in a randomized block design. Thirty seeds for each line were evenly planted in two rows of 2 m in length and 25 cm between rows. The main spikes of at least 6 plants in each plot were measured to investigate the SL and KNS when ripening. Broad-sense heritability across different environments was calculated based on the ANOVA model as described by Wu et al. .
The diverse wheat panel was planted in randomized complete blocks with two or three replicates in five major wheat ecological regions of China in the 2007, 2008, and 2009 planting seasons as described previously , including Yangling (108.08°E, 34.27°N) in Shaanxi Province, Tai′an (117.09°E, 36.21°N) in Shandong Province, Shijiazhuang (114.52°E, 38.05°N) in Hebei Province, Chengdu (104.08°E, 30.66°N) in Sichuan Province, and Yangzhou (119.42°E, 32.40°N) in Jiangsu Province. Two hundred seeds for each cultivar were evenly planted in five rows 2 m long and spaced 30 cm apart. The SL and KNS traits were assessed from 10 spikes randomly sampled from the centre of each plot before harvesting.
SNP genotyping, linkage map construction and QTL detection
The RIL lines and two parents were genotyped with the high-density Illumina Infinium iSelect 55 K SNP array by China Golden Marker (Beijing, China). The diverse wheat panel was also genotyped using the same SNP array . After excluding the monomorphic markers in the RIL population, markers retained were analyzed using the BIN function of IciMapping 4.2 (http://www.isbreeding.net) based on their segregation patterns with the parameters of “Missing Rates” and “Distortion Value” being set as 20 and 0.001, respectively. Only one marker with the least “Missing Rate” was chosen to represent each bin for constructing genetic maps and QTL mapping in this study. Linkage analysis was performed with IciMapping 4.2 using the default mapping function, and the resulting genetic map was displayed with MapChart v2.2 (http://www.biometris.nl/uk/Software/MapChart/). QTLs for SL and KNS in each environment and the average values across all environments were detected using the inclusive composite interval mapping (ICIM) function of IciMapping 4.2 and LOD score values ≥2.5.
QTL validation and development of KASP markers
For certain stable QTLs identified for SL and KNS in the RIL population, the QTL-associated flanking markers were validated using the diverse wheat panel. Furthermore, of the stable QTLs analyzed at which the positive alleles were contributed by Bima4 in the RIL population, the transmission of Bima4 alleles at the QTL-associated SNP markers were also determined using its 70 descendants. SNP markers highly associated with a specific QTL were selected and converted to KASP markers.
Availability of data and materials
The main datasets supporting the conclusions of this article are included within the article and its additional file. A small piece of data used for QTL validation in this study, including the diverse wheat panel and their phenotype and SNP genotyping data, is available in this reference which we published before (Li X, Xu X, Liu W, Li X, Yang X, Ru Z, Li L. Dissection of superior alleles for yield-related traits and their distribution in important cultivars of wheat by association mapping. Front Plant Sci. 2020;11: 175. DOI:https://doi.org/10.3389/fpls.2020.00175). All datasets are available from the corresponding author on reasonable request.
Kernel number per spike
Quantitative trait loci
Recombinant inbred line
Kompetitive allele-specific PCR
Würschum T, Leiser WL, Langer SM, Tucker MR, Longin CFH. Phenotypic and genetic analysis of spike and kernel characteristics in wheat reveals long-term genetic trends of grain yield components. Theor Appl Genet. 2018;131:2071–84.
Sadras V, Slafer G. Environmental modulation of yield components in cereals: heritabilities reveal a hierarchy of phenotypic plasticities. Field Crop Res. 2012;127:215–24.
Cui F, Ding A, Li J, Zhao C, Wang L, Wang X, et al. QTL detection of seven spike-related traits and their genetic correlations in wheat using two related RIL populations. Euphytica. 2012;186:177–92.
Donmez E, Sears R, Shroyer J, Paulsen G. Genetic gain in yield attributes of winter wheat in the Great Plains. Crop Sci. 2001;41:1412–9.
Jones S, Farooqi A, Foulkes J, Sparkes DL, Linforth R, Ray RV. Canopy and ear traits associated with avoidance of fusarium head blight in wheat. Front Plant Sci. 2018;9:1021.
Borner A, Schumann E, Furste A, Coster H, Leithold B, Röder S, et al. Mapping of quantitative trait loci determining agronomic important characters in hexaploid wheat (Triticum aestivum L.). Theor Appl Genet. 2002;105:921–36.
Marza F, Bai GH, Carver BF, Zhou WC. Quantitative trait loci for yield and related traits in the wheat population Ning7840 × Clark. Theor Appl Genet. 2006;112:688–98.
Narasimhamoorthy B, Gill BS, Fritz AK, Nelson JC, Brown-Guedira GL. Advanced backcross QTL analysis of a hard winter wheat × synthetic wheat population. Theor Appl Genet. 2006;112:787–96.
Kirigwi FM, Van Ginkel M, Brown-Guedira G, Gill BS, Paulsen GM, Fritz AK. Markers associated with a QTL for grain yield in wheat under drought. Mol Breed. 2007;20:401–13.
Kumar N, Kulwal PL, Balyan HS, Gupta PK. QTL mapping for yield and yield contributing traits in two mapping populations of bread wheat. Mol Breed. 2007;19:163–77.
Wang RX, Hai L, Zhang XY, You GX, Yan CS, Xiao SH. QTL mapping for grain filling rate and yield-related traits in RILs of the Chinese winter wheat population Heshangmai × Yu8679. Theor Appl Genet. 2009;118:313–25.
McIntyre CL, Mathews KL, Rattey A, Chapman SC, Drenth J, Ghaderi M, et al. Molecular detection of genomic regions associated with grain yield and yield-related components in an elite bread wheat cross evaluated under irrigated and rainfed conditions. Theor Appl Genet. 2010;120:527–41.
Wang J, Liu W, Wang H, Li L, Wu J, Yang X, et al. QTL mapping of yield-related traits in the wheat germplasm 3228. Euphytica. 2011;177:277–92.
Bennett D, Reynolds M, Mullan D, Izanloo A, Kuchel H, Langridge P, et al. Detection of two major grain yield QTL in bread wheat (Triticum aestivum L.) under heat, drought and high yield potential environments. Theor Appl Genet. 2012;125:1473–85.
Sourdille P, Cadalen T, Guyomarc'h H, Snape JW, Perretant MR, Charmet G, et al. An update of the Courtot × Chinese spring intervarietal molecular marker linkage map for the QTL detection of agronomic traits in wheat. Theor Appl Genet. 2003;106:530–8.
Suenaga K, Khairallah M, William HM, Hoisington DA. A new intervarietal linkage map and its application for quantitative trait locus analysis of "gigas" features in bread wheat. Genome. 2005;48:65–75.
Chu CG, Xu SS, Friesen TL, Faris JD. Whole genome mapping in a wheat doubled haploid population using SSRs and TRAPs and the identification of QTL for agronomic traits. Mol Breed. 2008;22:251–66.
Jantasuriyarat C, Vales M, Watson C, Riera-Lizarazu O. Identification and mapping of genetic loci affecting the free-threshing habit and spike compactness in wheat (Triticum aestivum L.). Theor Appl Genet. 2004;108:261–73.
Xie Q, Mayes S, Sparkes DL. Spelt as a genetic resource for yield component improvement in bread wheat. Crop Sci. 2015;55:2753–65.
Deng Z, Cui Y, Han Q, Fang W, Li J, Tian J. Discovery of consistent QTLs of wheat spike-related traits under nitrogen treatment at different development stages. Front Plant Sci. 2017;8:2120.
Liu K, Sun X, Ning T, Duan X, Wang Q, Liu T, et al. Genetic dissection of wheat panicle traits using linkage analysis and a genome-wide association study. Theor Appl Genet. 2018;131:1073–90.
Yao H, Xie Q, Xue S, Luo J, Lu J, Kong Z, et al. HL2 on chromosome 7D of wheat (Triticum aestivum L.) regulates both head length and spikelet number. Theor Appl Genet. 2019;132:1789–97.
Guo J, Shi W, Zhang Z, Cheng J, Sun D, Yu J, et al. Association of yield-related traits in founder genotypes and derivatives of common wheat (Triticum aestivum L.). BMC Plant Biol. 2018;18:38.
Pestsova E, Röder M. Microsatellite analysis of wheat chromosome 2D allows the reconstruction of chromosomal inheritance in pedigrees of breeding programmes. Theor Appl Genet. 2002;106:84–91.
Li XJ, Xu X, Yang XM, Li XQ, Liu WH, Gao AN, et al. Genetic diversity among a founder parent and widely grown wheat cultivars derived from the same origin based on morphological traits and microsatellite markers. Crop Pasture Sci. 2012;63:303–10.
Ma F, Xu Y, Ma Z, Li L, An D. Genome-wide association and validation of key loci for yield-related traits in wheat founder parent Xiaoyan 6. Mol Breed. 2018;38:91.
Zhuang QS. Chinese wheat improvement and pedigree analysis. Beijing: China Agricultural Press; 2003. (In Chinese)
Cui F, Zhang N, Fan XL, Zhang W, Zhao CH, Yang LJ, et al. Utilization of a Wheat660K SNP array-derived high-density genetic map for high-resolution mapping of a major QTL for kernel number. Sci Rep. 2017;7:3788.
Gao F, Wen W, Liu J, Rasheed A, Yin G, Xia X, et al. Genome-wide linkage mapping of QTL for yield components, plant height and yield-related physiological traits in the Chinese wheat cross Zhou 8425B/Chinese spring. Front Plant Sci. 2015;6:1099.
Kuang CH, Zhao XF, Yang K, Zhang ZP, Ding L, Pu ZE, et al. Mapping and characterization of major QTL for spike traits in common wheat. Physiol Mol Biol Plants. 2020;26:1295–307.
Huang XQ, Kempf H, Ganal MW, Röder MS. Advanced backcross QTL analysis in progenies derived from a cross between a German elite winter wheat variety and a synthetic wheat (Triticum aestivum L.). Theor Appl Genet. 2004;109:933–43.
Wu X, Cheng R, Xue S, Kong Z, Wan H, Li G, et al. Precise mapping of a quantitative trait locus interval for spike length and grain weight in bread wheat (Triticum aestivum L.). Mol Breed. 2014;33:129–38.
Chai L, Chen Z, Bian R, Zhai H, Cheng X, Peng H, et al. Dissection of two quantitative trait loci with pleiotropic effects on plant height and spike length linked in coupling phase on the short arm of chromosome 2D of common wheat (Triticum aestivum L.). Theor Appl Genet. 2018;131:2621–37.
Korzun V, Röder M, Ganal M, Worland A, Law C. Genetic analysis of the dwarfing gene (Rht8) in wheat. Part I. molecular mapping of Rht8 on the short arm of chromosome 2D of bread wheat (Triticum aestivum L.). Theor Appl Genet. 1998;96:1104–9.
Worland A, Sayers E, Korzun V. Allelic variation at the dwarfing gene Rht8 locus and its significance in international breeding programmes. Euphytica. 2001;119:155–9.
Gasperini D, Greenland A, Hedden P, Dreos R, Harwood W, Griffiths S. Genetic and physiological analysis of Rht8 in bread wheat: an alternative source of semi-dwarfism with a reduced sensitivity to brassinosteroids. J Exp Bot. 2012;63:4419–36.
Kowalski AM, Gooding M, Ferrante A, Slafer GA, Orford S, Gasperini D, et al. Agronomic assessment of the wheat semi-dwarfing gene Rht8 in contrasting nitrogen treatments and water regimes. Field Crop Res. 2016;191:150–60.
Fan X, Cui F, Ji J, Zhang W, Zhao X, Liu J, et al. Dissection of pleiotropic QTL regions controlling wheat spike characteristics under different nitrogen treatments using traditional and conditional QTL mapping. Front Plant Sci. 2019;10:187.
Russell J, Ellis R, Thomas W, Waugh R, Provan J, Booth A, et al. A retrospective analysis of spring barley germplasm development from ‘foundation genotypes’ to currently successful cultivars. Mol Breed. 2000;6:553–68.
Sjakste T, Rashal I, Röder M. Inheritance of microsatellite alleles in pedigrees of Latvian barley varieties and related European ancestors. Theor Appl Genet. 2003;106:539–49.
Rasheed A, Wen W, Gao F, Zhai S, Jin H, Liu J, et al. Development and validation of KASP assays for genes underpinning key economic traits in bread wheat. Theor Appl Genet. 2016;129:1843–60.
Cao P, Liang X, Zhao H, Feng B, Xu E, Wang L, et al. Identification of the quantitative trait loci controlling spike-related traits in hexaploid wheat (Triticum aestivum L.). Planta. 2019;250:1967–81.
Wang R, Liu Y, Isham K, Zhao W, Wheeler J, Klassen N, et al. QTL identification and KASP marker development for productive tiller and fertiles pikelet numbers in two high-yielding hard white spring wheat cultivars. Mol Breed. 2018;38:135.
Thomson M. High-throughput SNP genotyping to accelerate crop improvement. Plant Breed Biotechnol. 2014;2:195–212.
Rasheed A, Hao Y, Xia X, Khan A, Xu Y, Varshney R, et al. Crop breeding chips and genotyping platforms: progress, challenges, and perspectives. Mol Plant. 2017;10:1047–64.
Li X, Xu X, Liu W, Li X, Yang X, Ru Z, et al. Dissection of superior alleles for yield-related traits and their distribution in important cultivars of wheat by association mapping. Front Plant Sci. 2020;11:175.
Wu Q, Chen Y, Fu L, Zhou S, Chen J, Zhao X, et al. QTL mapping of flag leaf traits in common wheat using an integrated high-density SSR and SNP genetic linkage map. Euphytica. 2015;208:337–51.
The funders had no role in the experiment design, data analysis, decision to publish, or preparation of the manuscript. This work was supported by grants from the Natural Science Foundation of China (31571752), the Science and Technology Research Project in Henan Province (212102110267, 222102110206 and 212102110273), the Science and Technology Research Project in Xinxiang City (GG2020015), and the Project of Henan Provincial Youth Backbone Teachers (2019GGJS246).
Ethics approval and consent to participate
This study on plants complies with relevant institutional, national, and international guidelines and legislation, including the collection of plant material.
Consent for publication
The authors declare that there are no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Genetic map of wheat developed using the RIL population derived from the cross between BainongAK58 and Bima4. Table S2. Phenotypic variation for SL and KNS in the RIL population. Table S3. Seventy wheat cultivars derived from Bima4.
About this article
Cite this article
Xu, X., Li, X., Zhang, D. et al. Identification and validation of QTLs for kernel number per spike and spike length in two founder genotypes of wheat. BMC Plant Biol 22, 146 (2022). https://doi.org/10.1186/s12870-022-03544-6
- Kernel number per spike
- Spike length
- KASP markers