Next Article in Journal
Essential Oil and Major Non-Volatile Secondary Metabolites from the Leaves of Amazonian Piper subscutatum
Next Article in Special Issue
The Soybean High Density ‘Forrest’ by ‘Williams 82’ SNP-Based Genetic Linkage Map Identifies QTL and Candidate Genes for Seed Isoflavone Content
Previous Article in Journal
First Report of Sesame Mutants Tolerant to Severe Drought Stress during Germination and Early Seedling Growth Stages
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Utilization of a Wheat50K SNP Microarray-Derived High-Density Genetic Map for QTL Mapping of Plant Height and Grain Traits in Wheat

1
College of Agronomy, Northwest A&F University, Xianyang 712100, China
2
Zhumadian Academy of Agricultural Sciences, Zhumadian 463000, China
3
Nanyang Academy of Agricultural Sciences, Nanyang 473000, China
*
Authors to whom correspondence should be addressed.
These authors have equal contribution.
Plants 2021, 10(6), 1167; https://doi.org/10.3390/plants10061167
Submission received: 17 April 2021 / Revised: 18 May 2021 / Accepted: 26 May 2021 / Published: 8 June 2021
(This article belongs to the Special Issue QTL Mapping of Seed Quality Traits in Crops)

Abstract

:
Plant height is significantly correlated with grain traits, which is a component of wheat yield. The purpose of this study is to investigate the main quantitative trait loci (QTLs) that control plant height and grain-related traits in multiple environments. In this study, we constructed a high-density genetic linkage map using the Wheat50K SNP Array to map QTLs for these traits in 198 recombinant inbred lines (RILs). The two ends of the chromosome were identified as recombination-rich areas in all chromosomes except chromosome 1B. Both the genetic map and the physical map showed a significant correlation, with a correlation coefficient between 0.63 and 0.99. However, there was almost no recombination between 1RS and 1BS. In terms of plant height, 1RS contributed to the reduction of plant height by 3.43 cm. In terms of grain length, 1RS contributed to the elongation of grain by 0.11 mm. A total of 43 QTLs were identified, including eight QTLs for plant height (PH), 11 QTLs for thousand grain weight (TGW), 15 QTLs for grain length (GL), and nine QTLs for grain width (GW), which explained 1.36–33.08% of the phenotypic variation. Seven were environment-stable QTLs, including two loci (Qph.nwafu-4B and Qph.nwafu-4D) that determined plant height. The explanation rates of phenotypic variation were 7.39–12.26% and 20.11–27.08%, respectively. One QTL, Qtgw.nwafu-4B, which influenced TGW, showed an explanation rate of 3.43–6.85% for phenotypic variation. Two co-segregating KASP markers were developed, and the physical locations corresponding to KASP_AX-109316968 and KASP_AX-109519968 were 25.888344 MB and 25.847691 MB, respectively. Qph.nwafu-4B, controlling plant height, and Qtgw.nwafu-4B, controlling TGW, had an obvious linkage relationship, with a distance of 7–8 cM. Breeding is based on molecular markers that control plant height and thousand-grain weight by selecting strains with low plant height and large grain weight. Another QTL, Qgw.nwafu-4D, which determined grain width, had an explanation rate of 3.43–6.85%. Three loci that affected grain length were Qgl.nwafu-5A, Qgl.nwafu-5D.2, and Qgl.nwafu-6B, illustrating the explanation rates of phenotypic variation as 6.72–9.59%, 5.62–7.75%, and 6.68–10.73%, respectively. Two QTL clusters were identified on chromosomes 4B and 4D.

1. Introduction

Wheat (Triticum aestivum L.) is a major food crop globally, providing carbohydrates and protein for 35% of the global population. It is estimated that wheat production will increase by more than 70% in the next 30 years to meet the needs of the growing population [1]. To ensure global food security, genetic improvement of food production will be one of the main goals of wheat breeding programs [2,3,4].
Both 1000-grain weight (TGW) and the genetic improvement of related traits, which play a vital role in wheat yield, are applicable to increasing wheat yield. TGW is mainly affected by grain morphological parameters, such as grain length and grain width [4,5,6]. TGW-related genes, including sucrose synthase genes, encode cell wall invertase and cytokinin oxidase/dehydrogenase. The sucrose synthase genes TaSus1-7A, -7B and TaSus2-2A, -2B determine TGW and grain size [7,8], TaGW2-6A, -6B the grain width [9,10], and TaGS-D1 the grain size [11]. TaCwi-A1 encodes cell wall invertase [12], TaCKX6-D1 encodes cytokinin oxidase/dehydrogenase [13], and TaGASR-A1 is a putative Snakin/GASA protein associated with grain length (GL) (Dong et al., 2014). The inheritance of grain traits is relatively stable, forming a higher heritability than overall yield [14]. The method is suitable for QTL analysis of wheat samples planted and collected from different places and years, and a stable QTL can be retrieved and detected. Over the past 20 years, more than 150 QTLs related to TGW, grain length, and grain width have been identified, which are distributed on 21 chromosomes of wheat [5,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45]. Some studies have shown that there is a significant positive correlation between plant height and TGW [19,32,33,39,46,47]. The application of Rht1 (RHT-B1b) and Rht2 (RHT-D1b) in the 1960s set off a green revolution in wheat breeding. So far, 25 Rht genes have been identified in wheat [48,49]. Amongst these 25 genes, Rht1 and Rht2 are dwarfing genes that show insensitivity to gibberellins located on chromosomes 4BS and 4DS, respectively [13]. The wild alleles Rht-B1a and Rht-D1a also have a significant positive correlation with TGW [32,50]. Another gene, called Rht8, is sensitive to gibberellins for reducing plant height and is located on the 2DS chromosome. Rht8 is another widely applied dwarfing gene that has no obvious negative effect on TGW, but affects panicle length. Thus, Rht8 is a typical pleiotropism gene [6,51]. The genetic relationship can be investigated by targeting gene loci related to TGW and plant height, obtained by QTL mapping [30,52].
QTL genetic mapping is a crucial means to analyze functional loci [28]. Constructing a saturated genetic map is the key to QTL mapping, and molecular markers are the genetic map carrier. Triticum aestivum L. is a typical allohexaploid (AABBDD) composed of three subunits, and it represents the largest crop genome. Moreover, it is also the genome with the highest proportion of repetitive sequences such as transposable elements (84.7%) (IWGSC2018). Multitudes of SNP markers bear abundant polymorphism [53], and mapping results are quite advantageous in terms of accuracy and precision, especially for QTL mapping of quantitative traits [53,54]. By constructing a high-density genetic map to target the SNPs’ genetic and physical loction, collinearity analysis is performed, and then the recombination rate in different regions of the chromosome can be judged. After comparing the genetic and physical distances between adjacent markers, the relative changes of recombination rates in each chromosome can be further investigated and analyzed. The range of the mating population required for a recombination event in a specific region can be estimated. Scientific and accurate estimation for breaking the chain of specific target areas can be provided, and accurate judgments for evaluating genetic linkage drag, together with guidance for improving breeding efficiency, can be achieved [54,55].
Until now, couples of common wheat SNP microarrays, including Wheat9K [56], Wheat90K [37,57,58], Wheat820K [59], Wheat660K (http://bioservices.capitalbio.com/index.shtml) [37,57,58], and the Wheat55K SNP array, have been developed based on the 660K SNP array [60,61,62,63,64,65].The Wheat50K SNP array is a high-efficiency genotyping technology completed by the Institute of Crop Science of the Chinese Academy of Agricultural Sciences and Affymetrix. The technology was developed using high-quality SNP markers selected from Wheat90K SNP arrays, 660K SNP arrays, and 35K SNP arrays. In the 50K SNP array, there are 135 functional markers and 700 SNP markers closely linked to known QTLs [66]. The functional markers covering ten TGW-related genes and two plant height-related genes are shown in Table S1.
In this study, a Kompetitive allele-specific PCR (KASP) marker was used, which is a polymerase chain reaction-based (PCR) technology using fluorescence for single nucleotide polymorphism (SNP) and small insertion and deletion (InDel). KASP markers have the advantage of a low error rate and a relatively low cost compared to other SNP genotyping platforms such as TaqMan systems. According to the method of Ma et al. [63], SNPs located in the main QTL interval were selected to develop KASP markers.
This project aims to determine the chromosome recombination rates in different regions using collinearity analysis of the genetic positions and physical locations of the SNP markers. By mapping the environment-stable QTL region of grain-related traits, whether corresponding loci are located in the recombination-rich or recombination-barren area can be confirmed, and a reasonable judgment for further fine mapping can be fulfilled. By traits and linkage analysis of the relationship between plant height and grain traits, useful insights for the next steps of molecular breeding can also be provided.

2. Results

2.1. Agronomic Traits Analysis

As was shown in Table 1, significant differences when p = 0.01 in the four environments appeared in relation to the plant height, TGW, grain length, and grain width of the two-parent materials. In Table 1 and Figure S1, we can see that fluctuations occurred in the same traits in different environments, indicating that these four traits were easily affected by the environment. The agronomic traits failed to accord with a strictly normal distribution (p < 0.05). The trait heritability values of plant height, TGW, grain length, and grain width were 0.73, 0.62, 0.61, and 0.72, respectively. As can be seen, those of plant height and grain width were relatively high.
As was shown in Figure S1, there was a significant positive correlation between the same traits and different environments when p = 0.001. The correlation between different years in the same place was higher than that in other combinations, indicating that a high degree of environmental similarity was present in the same place but in different years. The correlation between plant height and grain length was negative, but there was a significant positive correlation between TGW and grain width. TGW had a significant positive correlation with the other three traits, and a higher correlation with grain width than that with other traits. The correlation between grain length and grain width was different in different environments.

2.2. Construction of a Genetic Map

2.2.1. Description and Illustration of a Genetic Map

66,832 markers were subject to polymorphism analysis of population genotype by 50K gene microarray. A total of 19,601 SNP markers with differences were screened in the derived RIL populations of Xinong1376 and Xiaoyan81, while the remaining 15,822 markers were filtered by Chi-square test, and redundant markers were eliminated using the bin function of IciMapping. A total of 3136 bin markers, including 15,576 SNP markers, were eventually anchored to the genetic map. In addition, the genotyping, polymorphism marker, data filtering, physical map, genetic map, and bin map are all shown and illustrated in Table S2. Based on the 660K chip labeling, the SNP markers that differed between the two parents were detected and stored in Figure S2. The total length of the linkage map was 4512.79 cM, the average map distance was 1.44 cM, and the maximum gap was 26.86 cM, which covered 21 wheat chromosomes. According to linkage lengths in the homologous groups, their sequence in descending order was the fifth, the seventh, the third, the second, the fourth, the sixth, and the first. The linkage lengths were 813.14 cM, 794.35 cM, 703.96 cM, 631.98 cM, 563.99 cM, 537.27 cM, and 468.12 cM, and the numbers of bin markers were 621 (2947 SNP markers included), 549 (2193 SNP markers included), 524 (2846 SNP markers included), 327 (1865 SNP markers included), 393 (2002 SNP markers included), 372 (1865 SNP markers included), and 272 (2016 SNP markers included), respectively.
The numbers of bin markers located in wheat A, B, and D chromosome groups were 1231, 1197, and 708, respectively. The linkage lengths were 1703.69 cM, 1298.23 cM and 1510.87 cM, and the average map distances were 1.38 cM, 1.08 cM, and 2.13 cM, respectively. Molecular markers in the D genome were no more than those in the other two subgroups. In addition, the longest linkage group corresponding to chromosome 3A was 312.11 cM, and the shortest corresponding to chromosome 1D was 130.85 cM. The maps of each linkage group were shown and illustrated in Table 2 and Figure S2.

2.2.2. Collinearity Analysis of the Genetic Map

In this research, the genetic map and the collinearity map of the reference genome were analyzed as follows: The whole chromosome was included in the genetic map, the genetic map and the physical map were collinear, and the linkage map and the physical map were not linearly related. The recombination exchange on chromosomes was unbalanced, and the collinear diagrams of other chromosomes except for chromosome 1B appeared by and large S-shaped. The genetic positions of chromosomes increased linearly with the increase in physical locations, and the rest of the genetic positions aligned constantly with the increase in physical locations. This indicated that the two ends of the chromosome were recombination-rich areas and that the middle region was a recombination-barren area. A significant correlation of the genetic map and the physical one appeared when p = 0.001, the correlation coefficient ranged from 0.63 to 0.99, and the correlation coefficient of chromosome 1B was 0.63. The distribution presentation of bin markers on the reference genome showed that the number of bin markers on both ends of the chromosome was significantly higher than that of the middle region. The recombination rate of the two sides with a U-shaped distribution was significantly higher than that of the middle region, which confirmed that the ends of the chromosome were recombination-rich areas and the middle was the recombination-barren area. The reason for these findings was the inhibitory effect of centromere recombination.
No markers could be detected in the middle regions (more than 200 MB) of chromosomes 1D, 5A, and 6A. However, the linkage group was not divided into two parts in these regions, which were supposed to be recombination-barren regions. For nine chromosomes (2D, 3D, 5A, 5B, 5D, 6A, 6D, 7A and 7D), each chromosome included two linkage groups. For different linkage groups corresponding to the same chromosome, the grouping regions all appeared at both ends of the chromosome as the recombination-rich area, and the physical distance between the markers was less than 30 MB.
The collinearity map of chromosome 1B from 0 to 480 MB presented as an L-type curve. Although the gradual numerical values of physical location increased, the genetic distances were almost unchanged, and thus homologous recombination hardly occurred in the region. Xinong1376 belonged to 1BL/1RS translocation line, 1RS and 1BS hardly recombined, and the centromere’s inhibition of recombination happened in the middle region, making the collinearity map L-shaped.

2.2.3. Effects of 1B/1R on Traits Related to Plant Height and TGW

1RS specific marker was used to detect the population, the strains containing 1RS and 1BS were 51 and 147, respectively, and the p value of the chi-square test was 8.95 × 10−12, which proved to be a severely segregated marker that couldn’t be linked to the linkage group. According to the typing of the specific markers, the unpaired data T test was performed on the traits related to plant height and TGW, and there was no significant difference between 1RS and 1BS. According to the typing of specific markers, a two-factor analysis of variance was performed on the agronomic traits, and the TGW and grain width were not affected by the genotype. According to the results of the variance analysis, Duncan’s new multiple range test comparison of plant height and grain length was conducted. In terms of plant height, 1RS contributed to the reduction of plant height by 3.43cm. In terms of grain length, 1RS contributed to the elongation of grain by 0.11mm (shown in Table S3 and Figure S3).

2.3. QTL Mapping Analysis

A total of 43 QTLs for PH, TGW, GL, and GW were identified by QTL mapping analysis (Table 3 and Figure S4). These QTLs with LOD values ranging from 2.51 to 53.34 were distributed on 15 chromosomes and explained 1.36–33.08% of the phenotypic variation (Table 3 and Figure S4). There were 8, 11, 15, and 9 QTLs detected for PH, TGW, GL, and GW, respectively (Table 3 and Figure S4).
Inclusive composite interval mapping (ICIM) for PH identified a total of eight QTLs, which were located on six different chromosomes (Table 3 and Figure S4): 2D(2), 4B, 4D, 5B, 5D, and 6B(2). The QTL on 4B, Qph.nwafu-4B, was detected in four environments. Qph.nwafu-4B was thus treated as a major QTL, which explained 9.32–13.76% of phenotypic variance with LOD values ranging from 7.93 to 26.85. As was expected, the positive allele of Qph.nwafu-4B was contributed by Xiaoyan81 (Table 3 and Figure S4). The QTL on 4D, Qph.nwafu-4D, was detected in each of four environments. Qph.nwafu-4D was thus treated as a major QTL, which explained 20.11–27.09% of phenotypic variance with LOD values ranging from 16.78 to 42.21. As we expected, the positive allele of Qph.nwafu-4D was contributed by Xinong1376 (Table 3 and Figure S4).
One QTL, Qph.nwafu-2D.1, for PH was detected in two environments, which explained 3.3–3.73% of phenotypic variance. The remaining QTLs were detected only in a single environment (Table 3 and Figure S4).
ICIM for TGW identified a total of eleven QTLs, which were located on eight different chromosomes (shown in Table 3 and Figure S4): 2A, 2B, 3A, 4B, 4D(2), 5A, 5D(3), and 6A. The QTL on 4B, Qtgw.nwafu-4B, was detected in three environments. Qtgw.nwafu-4B was thus treated as a stable QTL, which explained 3.43–6.85% of phenotypic variance with LOD values ranging from 2.85 to 4.37. As was expected, the positive allele of Qtgw.nwafu-4B was contributed by Xinong1376 (shown in Table 3 and Figure S4). Based on the initial QTL mapping results, we developed two KASP markers, KASP_AX-109316968 and KASP_AX-109333198, and integrated them into the genetic map. When remapping with this integrated KASP marker, it was indicated that Qtgw.nwafu-4B was located in a 5 cM interval on chromosome arm 4BS, between the markers of AX-111494900 and AX-94438527, containing the newly developed KASP markers, including KASP_AX-109316968 and KASP_AX-109333198 (Figure S5 and Table S3).Three QTLs, Qtgw.nwafu-4D.1, Qtgw.nwafu-5A, and Qtgw.nwafu-5D.1, for TGW were detected in each of two environments, which explained 2.85–14.79% of phenotypic variance. The remaining QTLs were detected only in a single environment (Table 3).
ICIM for GL identified a total of fifteen QTLs, which were located on ten different chromosomes (Table 3 and Figure S4): 1A, 1B(2), 3A, 4A, 4B(2), 4D, 5A, 5B, 5D(4), and 6B. The QTL on 6B, Qgl.nwafu-6B, was detected in four environments. Qgl.nwafu-6B was thus treated as a major QTL, which explained 6.68–10.73% of phenotypic variance with LOD values ranging from 3.35 to 8.79. As was expected, the positive allele of Qgl.nwafu-6B was contributed by Xinong1376 (Table 3). The QTL on 5A, Qgl.nwafu-5A, was detected in in three environments. Qgl.nwafu-5A was thus treated as a stable QTL, which explained 6.72–9.59% of phenotypic variance with LOD values ranging from 5.8 to 6.93. As we expected, the positive allele of Qgl.nwafu-5A was contributed by Xinong1376 (Table 3 and Figure S4). The QTL on 5D, Qgl.nwafu-5D.2, was detected in three environments. Qgl.nwafu-5D.2 was thus treated as a stable QTL, which explained 5.62–7.75% of phenotypic variance with LOD values ranging from 4.36 to 6.13. As was expected, the positive allele of Qgl.nwafu-5D.2 was contributed by Xinong1376 (Table 3 and Figure S4). Four QTLs, Qgl.nwafu-1B.2, Qgl.nwafu-3A, Qgl.nwafu-4A, and Qgl.nwafu-4B.2, for GL were detected in two environments, explaining 3.51–6.13% of phenotypic variance. The remaining QTLs were detected only in a single environment (Table 3 and Figure S4).
ICIM for GW identified a total of nine QTLs, which were located on seven different chromosomes (Table 3, Figure S4): 2B, 2D, 3A, 4B(3), 4D, 5D, and 6D. The QTL on 4D, Qgw.nwafu-4D, was detected in each of the four environments. Qgw.nwafu-4D was thus treated as a major QTL, which explained 6.32–12.12% of phenotypic variance with LOD values ranging from 3.2 to 7.93. As we expected, the positive allele of Qgw.nwafu-4D was contributed by Xinong1376 (Table 3). One QTL, Qgw.nwafu-4B.1, for GW was detected in two environments, which explained 6.85–6.95% of phenotypic variance. The remaining QTLs were detected only in a single environment (shown in Table 3 and Figure S4).
Two QTL clusters were identified on chromosomes 4B and 4D (Table 3 and Figure S4). For the QTL cluster on chromosome 4B, Qtgw.nwafu-4B for TGW was co-localized with Qgl.nwafu-4B.1 for GL, and Qph.nwafu-4B and Qgl.nwafu-4B.2 for GL were co-localized with Qgl.nwafu-4B.2 and Qgl.nwafu-4B.3 for GL in a region ranging from 51 cM to 77 cM. On chromosome 4D, Qph.nwafu-4D for PH was clustered with Qtgw.nwafu-4D.1 for TGW, and Qgw.nwafu-4D for GW was clustered with with the alleles from Xiaoyan81 increasing PH, TGW and GW.

3. Discussion

3.1. The Impact of Linkage Map on QTL Mapping

In this research, a linkage map, based on 50K microarray markers, was constructed from 198 F8 RIL lines derived from the combination of two parents, Xinong1376 and Xiaoyan81. The linkage map had a total length of 4512.79 cM, covering 21 chromosomes of wheat. The reason why no marks could be targeted in the regions of more than 200 MB in the middle of the four chromosomes 1D, 5A, and 6A was that a recombination-barren area near the centromere appeared in the above regions, as was shown in Figure 1. Both parents were derived from the backbone parent Xiaoyan6, and a region with the same haplotype was formed rapidly [68], so that the two parents had no markers with polymorphic differences in the above regions. There was a long, excellent haplotype segment on chromosome 6A [60,69].
In this study, 43 QTLs were located. The genetic distance confidence interval was 0.5–12.5 cM, and the physical distance of the markers on both sides was 0.0201 MB–414.88328 MB. As was shown in Table 2, the genetic distance confidence interval was not proportional to the physical distance, which reflected the imbalance of the recombination exchange on the chromosomes.
By combining Figure S4 and Figure 1, it appeared that there were 5 QTLs located in the recombination-barren region of the reference genome, and more than 20 MB QTLs were distributed in this candidate region. The linkage interval of Qgl.nwafu-1B.1 was 0–0.5 cM, while the physical interval was 59.47117 MB–94.978091 MB and the interval physical distance was 35.506914 MB. The reason was that Xinong1376 belonged to the 1BL/1RS translocation line, and there was almost no recombination or recombination disorder between 1RS and 1BS [6,38,70,71]. Although the genetic distance of the confidence interval was short, the corresponding physical distance of it was far. As was shown in Figure S4, the linkage region of Qtgw.nwafu-2B was 95.5 cM–106.5 cM, and no marks could be targeted in this region. This area belongs to the reorganization cold spot area, and the corresponding physical distance was 153.585606 MB–568.468886 MB. The linkage regions of Qtgw/gl.nwafu-3A, Qgl.nwafu-4A, and Qgw.nwafu-4B.2 were 132.5 cM–134.5 cM, 46.5 cM–53.5 cM, and 67.5 cM–69.5 cM, respectively, and the corresponding regions were 457.796943 MB–431.074614 MB, 407.389107 MB–129.089816 MB, and 114.952789 MB–161.548436 MB, respectively. As was shown in Table 2, the above three QTLs all fell in the recombination-barren region of linkage groups with a large physical interval. The confidence interval of Qgw.nwafu-6D, which was the largest, was 0 cM–12.5 cM, but the corresponding physical region was 12.650045 MB–8.255713 MB, and the interval was only 4.4 MB. Qgw.nwafu-6D was located at the top of the chromosome, and belonged to the recombination-rich region, with a big genetic distance but a short corresponding physical distance.

3.2. Comparison with Previous Research Results

Two loci as environment-stable QTLs, targeted in three or four kinds of environments, were Qph.nwafu-4B and Qph.nwafu-4D, which control plant height. In the confidence interval, the function markers including Rht-1 and Rht-2 were AX-179477460 and AX-86170701, respectively. According to the additive effect, the effect of the Qph.nwafu-4D mutant in lowering plant height was stronger than that of the Qph.nwafu-4B mutant, which was consistent with the results of Zhai et al. [6] The locus, Qgl.nwafu-5A, which controlled the grain length, corresponded to the physical location of 698.508129 MB–700.34701 MB, which was located at the end of the chromosome. Compared with the results of previous studies [23,29,30,31,32,33,34,35,36,37,38,39,40,41,42], Qgl.nwafu-5A was a new QTL. The location of Qgl.nwafu-5D.2 which controlled the length of the grain corresponded to the physical location of 370.135626 MB-386.126855 MB. Based on previous research [22,24,35,42,43], Qgl.nwafu-5D.2 was defined as a new QTL as well. The location of Qgl.nwafu-6B, which controls grain length, corresponded to the physical location of 704.884934 MB–718.376276 MB. Compared with the results of previous studies [35], the physical location marked by IWB2746 was 701.387367 MB. As was shown in Figure S4, the collinearity between the linkage group and the physical position was relatively disordered at the end of chromosome 6B, and it was not clear whether they were the same QTL.
Qph.nwafu-4B (controlling plant height) and Qtgw.nwafu-4B (controlling TGW) had an obvious linkage relationship, with a distance of 7–8 cM. The physical location corresponding to this location of Qph.nwafu-4B was 30.805339 MB–32.961929 MB, and the physical position corresponding to the location of Qtgw.nwafu-4B was 25.847125 MB–26.491497 MB. Guan’s QTL mapping results were marked as BS00084904_51 and BS00011338_51 on both sides, and the physical location was 28.954526 MB–66.811785 MB [30]. Cui Fa’s QTL mapping results were marked as Rht-B1 and Xmag2055 on both sides, and the physical location was 30.860778 MB–20.741542 MB [70]. Quarrie’s QTL mapping results were marked as Rht-B1 and gwm165.1 on both sides, and the physical location was 30.860778 MB–269.948831 MB [42] (The results of previous studies on chromosome 4B and the specific QTL information related to TGW are shown and illustrated in Table S4). From the QTL mapping results in this study and the above three research results, it was suggested that the confidence interval had this overlap while the confidence interval of this study was the shortest. Based on heredity Doumai/Shi 41875, Li mapped the plant height and TGW. The physical location on chromosome 4B was 46.621203 MB [35], which was not the same QTL. The confidence intervals of Qph.nwafu-4D, Qtgw.nwafu-4D.1, and Qgw.nwafu-4D had clear overlaps and were stably expressed in multiple environments. The mutant at this locus lowered plant height while also decreasing TGW and grain width. Rht2 had a significant effect on TGW, as previously shown by Mohler et al. [32]. There was a significant overlap in the confidence interval of Qph.nwafu-5D controlling plant height and Qgl.nwafu-5D.3 controlling grain length, with a typical pleiotropism. This locus’s physical position was 466.230408 MB–469.357817 MB, and its additive effect was opposite, so physiological antagonism occurred. The location of wmc215 targeted by Hai et al. was 472.369175 MB, and that of gwm212 targeted by Quarrie was 472.630187 MB, which was in line with previous localization results [42,43]. The difference in physical location was 3 MB. Since subgroup D had a large linkage disequilibrium [72], it was impossible to determine whether these loci were the same one. Qtgw.nwafu-5D.1 controlling TGW and Qgl.nwafu-5D.1 controlling grain length were located in the region from 38.070293 MB–41.294446 MB, neither of which belonged to the same region of the 5D chromosome, compared with the results of previous studies [35,42,43,73].

3.3. Qtgw.Nwafu-4B Molecular Marker Development

Based on the confidence interval of the parental 660K chip marker, two co-segregating KASP markers were developed. Two KASP molecular markers were inserted into the original genetic map, and the genetic map of chromosome 4B maintained a high degree of collinearity. Two KASP molecular markers were inserted into the original genetic map, and the genetic map of chromosome 4B maintained a high degree of collinearity. The primer sequences and typing information of the two molecular markers of KASP_AX-109316968 and KASP_AX-109333198 are shown in Figure S5 and Table S5. Qph.nwafu-4B (controlling plant height) and Qtgw.nwafu-4B (controlling TGW) had an obvious linkage relationship, with a distance of 7–8 cM. Breeding is based on molecular markers that control plant height and thousand-grain weight to select strains with low plant height and large grain weight.

4. Materials and Methods

4.1. Plant Materials, Experimental Design, and Investigation of Agronomic Traits

Xinong 1376 is the female parent and Xiaoyan 81 is the male parent. Based on the single-grain transmission method, 198 RIL lines were generated. There were planted in Yangling, Shaanxi province and Nanyang, Henan province, from October 2018 to June 2019 and from October 2019 to June 2020, respectively. A randomized block design (repeated five times, with two rows of districts, 2 m row length, 70 plants per row, and 0.3 m row spacing) was adopted in each experimental site. The other field managements were subject to the same treatment as the local. During the wax maturity period of wheat, five individual plants were sampled in sequence from the fifth plant of each family. Plant height, TGW, grain length and grain width were also measured. By R/lme4 [73], each environment’s agronomic traits were obtained for W-test, and then multiple comparisons of parental traits and calculation of heritability were completed. The heritability of the two traits was calculated by using the formula as follows:
H2 = VG∕(VG + VGY/y + VGE/e + VE/nr) × 100%
where y is the number of years, e is the number of environments, and n is the number of repetitions.
The pedigrees of Xinong1376 and Xiaoyan81 are illustrated in Figure S2.

4.2. Construction and Evaluation of Genetic Maps

The wheat genomic DNA, with tender wheat leaves as the plant material, was extracted by CTAB, and the quality and quantity of DNA were detected and confirmed. Meanwhile, the DNA of each line was hybridized on the wheat 50K SNP array containing 66,832 markers using Burdock Biotechnology (Beijing, China).
The course of constructing the map was conducted as follows: The BIN function of IciMapping 4.1 [70] was utilized to analyze the markers, and the markers with partial separation rate (p < 0.001) and missing rate (>15%) were removed. The Kosambi function with LOD ≥ 5 was applied to group the combined marker groups in JoinMap 4.0; Kosambi mapping of MSTmap [74], according to the clustering results, was used in the markers’ ordination. The flanking sequences of SNPs were BLAST aligned with the genome of IWGSC RefSeq v1.0 (http://www.wheatgenome.org/News/Latest-news/All-IWGSC-data-related-to-the-reference-sequence-of-bread-wheat-IWGSC-RefSeq-v1.0-publicly-available-at-URGI) to obtain their physical locations. The version of BLAST used was 2.2.31 –outfmt 3–num_alignments 5.

4.3. Identification of 1BL/1RS Translocation

1RS, applied to identify parents and populations as x-sec-p1/x-sec-p2, respectively, was a specific marker [75]. Xinong1376 was identified as a 1BS/1RS translocation line. 1B/1R genotyping and traits data were stored in Table S4 and Figure S3. Analysis of variance and Duncan’s new multiple range test comparisons based on genotype and trait were conducted.

4.4. Detection of Quantitative Loci

IciMapping 4.2 based on the biparental population (BIP) module with the inclusive composite interval mapping (ICIM, http://www.isbreeding.net/software/?type=detail&id=28) was used for QTL mapping on data obtained from different environments. QTL mapping of the phenotypic values in the four environments was carried out. The LOD value was determined in 1000 permutation tests with a = 0.05 (Type I Error) as the parameter, and the background was set and controlled by the positive and negative stepwise regression, with the step width set to 1cM. QTLs were named based on the International Rules of Genetic Nomenclature (http://whea.pw.usda.gov/ggpages/wgc/98/Intro.htm). Mapchart2.3 (https://www.wur.nl/en/Research-Results/Research-Institutes/plant-research/Biometris-1/SoftwareService/Download-MapChart.htm) was used for the drawing of the genetic and QTL mapping. The collinearity drawing of genetic and physical maps, and the calculation of correlation coefficient were conducted by package plotrix (https://cran.r-project.org/src/contrib/Archive/plotrix/) and package (https://github.com/braverock/PerformanceAnalytics) of R software.

4.5. Breeding Molecular Marker Development

After obtaining the preliminary QTL mapping results, we anchored the flanking markers to the physical map. In order to develop a competitive allele-specific PCR (KASP) marker that can be used to track stable TGW QTLs, we used the Wheat660K SNP array to further genotype the parents of the Xinong1376/Xiaoyan81 population [63,71]. According to the method of Ma et al. [63], SNPs located in the main QTL interval were selected to develop KASP markers. The developed integrated genetic map of KASP markers was applied to relocate the target QTL.

5. Conclusions

In this research, a genetic map covering the entire wheat genome was constructed, with a total of 3136 bin markers, including 15576 SNP markers, and the total length of the linkage map was 4512.79 cM. Except for chromosome 1B, the ends of chromosomes were identified as recombination-rich areas, while the middle areas were recombination-barren. Both the genetic map and the physical map showed a significant correlation when p = 0.001. The correlation coefficient ranged from 0.63 to 0.99. There was almost no recombination between 1RS and 1BS. Among 43 QTLs indirectly compared by reference genome, only 13 QTLs were consistent with previous mapping results, and 30 QTLs were defined as new QTLs. Seven environment-stable QTLs were detected in this population, including Qph.nwafu-4B, Qtgw.nwafu-4B, Qgw.nwafu-4D, Qph.nwafu-4D, Qgl.nwafu-5A, Qgl.nwafu-5D.2, and Qgl.nwafu-6B. Qtgw.nwafu-4B, which influenced TGW, showed an explanation rate of 3.43–6.85% for phenotypic variation, with two co-segregating KASP markers developed, and the physical locations corresponding to KASP_AX-109316968 and KASP_AX-109519968 were 25.888344 MB and 25.847691 MB, respectively, for details, see Figure 2. Qph.nwafu-4B (controlling plant height) and Qtgw.nwafu-4B (controlling TGW) had an obvious linkage relationship, with a distance of 7–8 cM. The physical location corresponding to this location of Qph.nwafu-4B was 30.805339 MB–32.961929 MB, and the physical position corresponding to this location of Qtgw.nwafu-4B was 25.847125 MB–26.491497 MB. There is a functional marker (AX-179477460) for the control value of plant height in the Qph.nwafu-4B confidence interval, and this locus can be determined to be Rht-B1. The physical locations of Qph.nwafu-4B, Qph.nwafu-4D, and Qgw.nwafu-4D were consistent with previous mapping results. For Qgl.nwafu-6B, it couldn’t be accurately determined whether it was a new QTL or not. Two QTL clusters were identified on chromosomes 4B and 4D (Table 3 and Figure S4).

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/plants10061167/s1, Table S1. Statistical analysis of parent and RIL lines for traits, Table S2. Single-nucleotide polymorphism (SNP) marker statistics about distribution and density on 21 wheat chromosomes derived from crossing between Xinong1376 and Xiaoyan81, Table S3. Full genomic QTL mapping results of plant height and grain related traits in the F8 RIL lines between Xinong1376 and Xiaoyan81, Figure S1. Collinearity analysis of genetic map and reference genome, Figure S2. The linkage group corresponding to chromosome 4B.

Author Contributions

D.S., D.H. and J.W. designed the research. C.Z. and D.L. conducted genotyping of the population. C.Z., D.L., R.Y., J.Y., X.S., J.J., P.S. and Z.Z. conducted phenotyping of the population. C.Z. analyzed all data. D.L. wrote the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by grants from the National Key Research and Development Program of China (2016YFD0101802).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Please refer to suggested Data Availability Statements in section “MDPI Research Data Policies” at https://www.mdpi.com/ethics.

Conflicts of Interest

This study did not have any conflict of interest.

References

  1. Bailey-Serres, J.; Parker, J.E.; Ainsworth, E.A.; Oldroyd, G.E.D.; Schroeder, J.I. Genetic strategies for improving crop yields. Nature 2019, 575, 109–118. [Google Scholar] [CrossRef] [Green Version]
  2. Hanif, M.; Gao, F.; Liu, J.; Wen, W.; Zhang, Y.; Rasheed, A.; Xia, X.; He, Z.; Cao, S. TaTGW6-A1, an ortholog of rice TGW6, is associated with grain weight and yield in bread wheat. Mol. Breed. 2016, 36, 1. [Google Scholar] [CrossRef]
  3. Kumari, S.; Jaiswal, V.; Mishra, V.K.; Paliwal, R.; Balyan, H.S.; Gupta, P.K. QTL mapping for some grain traits in bread wheat (Triticum aestivum L.). Physiol. Mol. Biol. Plants 2018, 24, 909–920. [Google Scholar] [CrossRef]
  4. Xie, Q.; Mayes, S.; Sparkes, D.L. Carpel size, grain filling, and morphology determine individual grain weight in wheat. J. Exp. Bot. 2015, 66, 6715–6730. [Google Scholar] [CrossRef] [Green Version]
  5. Guan, P.; Di, N.; Mu, Q.; Shen, X.; Wang, Y.; Wang, X.; Yu, K.; Song, W.; Chen, Y.; Xin, M.; et al. Use of near-isogenic lines to precisely map and validate a major QTL for grain weight on chromosome 4AL in bread wheat (Triticum aestivum L.). Theor. Appl. Genet. 2019, 132, 2367–2379. [Google Scholar] [CrossRef]
  6. Zhai, H.; Feng, Z.; Li, J.; Liu, X.; Xiao, S.; Ni, Z.; Sun, Q. QTL Analysis of Spike Morphological Traits and Plant Height in Winter Wheat (Triticum aestivum L.) Using a High-Density SNP and SSR-Based Linkage Map. Front. Plant Sci. 2016, 7, 1617. [Google Scholar] [CrossRef]
  7. Jiang, Q.; Hou, J.; Hao, C.; Wang, L.; Ge, H.; Dong, Y.; Zhang, X. The wheat (T. aestivum) sucrose synthase 2 gene (TaSus2) active in endosperm development is associated with yield traits. Funct. Integr. Genom. 2011, 11, 49–61. [Google Scholar] [CrossRef]
  8. Hou, J.; Jiang, Q.; Hao, C.; Wang, Y.; Zhang, H.; Zhang, X. Global selection on sucrose synthase haplotypes during a century of wheat breeding. Plant Physiol. 2014, 164, 1918–1929. [Google Scholar] [CrossRef] [Green Version]
  9. Su, Z.; Hao, C.; Wang, L.; Dong, Y.; Zhang, X. Identification and development of a functional marker of TaGW2 associated with grain weight in bread wheat (Triticum aestivum L.). Theor. Appl. Genet. 2011, 122, 211–223. [Google Scholar] [CrossRef]
  10. Yang, Z.; Bai, Z.; Li, X.; Wang, P.; Wu, Q.; Yang, L.; Li, L.; Li, X. SNP identification and allelic-specific PCR markers development for TaGW2, a gene linked to wheat kernel weight. Theor. Appl. Genet. 2012, 125, 1057–1068. [Google Scholar] [CrossRef]
  11. Ma, L.; Li, T.; Hao, C.; Wang, Y.; Chen, X.; Zhang, X. TaGS5–3A, a grain size gene selected during wheat improvement for larger kernel and yield. Plant Biotechnol. J. 2016, 14, 1269–1280. [Google Scholar] [CrossRef]
  12. Ma, D.; Yan, J.; He, Z.; Wu, L.; Xia, X. Characterization of a cell wall invertase gene TaCwi-A1 on common wheat chromosome 2A and development of functional markers. Mol. Breed. 2010, 29, 43–52. [Google Scholar] [CrossRef]
  13. Zhang, L.; Zhao, Y.L.; Gao, L.F.; Zhao, G.Y.; Zhou, R.H.; Zhang, B.S.; Jia, J.Z. TaCKX6-D1, the ortholog of rice OsCKX2, is associated with grain weight in hexaploid wheat. New Phytol. 2012, 195, 574–584. [Google Scholar] [CrossRef] [PubMed]
  14. Gupta, P.K.; Rustgi, S.; Kumar, N. Genetic and molecular basis of grain size and grain number and its relevance to grain productivity in higher plants. Genome 2006, 49, 565–571. [Google Scholar] [CrossRef] [Green Version]
  15. Varshney, R.K.; Prasad, M.; Roy, J.K.; Kumar, N.; Harjit, S.; Dhaliwal, H.S.; Balyan, H.S.; Gupta, P.K. Identification of eight chromosomes and a microsatellite marker on 1AS associated with QTL for grain weight in bread wheat. Theor. Appl. Genet. 2000, 100, 1290–1294. [Google Scholar] [CrossRef]
  16. Cao, P.; Liang, X.; Zhao, H.; Feng, B.; Xu, E.; Wang, L.; Hu, Y. Identification of the quantitative trait loci controlling spike-related traits in hexaploid wheat (Triticum aestivum L.). Planta 2019, 250, 1967–1981. [Google Scholar] [CrossRef]
  17. Campbell, B.T.; Baenziger, P.S.; Gill, K.S.; Eskridge, K.M.; Budak, H.; Erayman, M.; Dweikat, I.; Yen, Y. Identification of QTLs and Environmental Interactions Associated with Agronomic Traits on Chromosome 3A of Wheat. Crop Sci. 2003, 43, 1493–1505. [Google Scholar] [CrossRef]
  18. Borner, A.; Schumann, E.; Furste, A.; Coster, H.; Leithold, B.; Roder, S.; Weber, E. Mapping of quantitative trait loci determining agronomic important characters in hexaploid wheat (Triticum aestivum L.). Theor. Appl. Genet. 2002, 105, 921–936. [Google Scholar] [CrossRef]
  19. Huang, X.Q.; Coster, H.; Ganal, M.W.; Roder, M.S. Advanced backcross QTL analysis for the identification of quantitative trait loci alleles from wild relatives of wheat (Triticum aestivum L.). Theor. Appl. Genet. 2003, 106, 1379–1389. [Google Scholar] [CrossRef] [PubMed]
  20. Huang, X.Q.; Cloutier, S.; Lycar, L.; Radovanovic, N.; Humphreys, D.G.; Noll, J.S.; Somers, D.J.; Brown, P.D. Molecular detection of QTLs for agronomic and quality traits in a doubled haploid population derived from two Canadian wheats (Triticum aestivum L.). Theor. Appl. Genet. 2006, 113, 753–766. [Google Scholar] [CrossRef] [PubMed]
  21. Narasimhamoorthy, B.; Gill, B.S.; Fritz, A.K.; Nelson, J.C.; Brown-Guedira, G.L. Advanced backcross QTL analysis of a hard winter wheat × synthetic wheat population. Theor. Appl. Genet. 2006, 112, 787–796. [Google Scholar] [CrossRef] [PubMed]
  22. Breseghello, F.; Sorrells, M.E. Association mapping of kernel size and milling quality in wheat (Triticum aestivum L.) cultivars. Genetics 2006, 172, 1165–1177. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Breseghello, F.; Sorrells, M.E. QTL analysis of kernel size and shape in two hexaploid wheat mapping populations. Field Crops Res. 2007, 101, 172–179. [Google Scholar] [CrossRef]
  24. Ramya, P.; Chaubal, A.; Kulkarni, K.; Gupta, L.; Kadoo, N.; Dhaliwal, H.S.; Chhuneja, P.; Lagu, M.; Gupta, V. QTL mapping of 1000-kernel weight, kernel length, and kernel width in bread wheat (Triticum aestivum L.). J. Appl. Genet. 2010, 51, 421–429. [Google Scholar] [CrossRef] [PubMed]
  25. Mir, R.R.; Kumar, N.; Jaiswal, V.; Girdharwal, N.; Prasad, M.; Balyan, H.S.; Gupta, P.K. Genetic dissection of grain weight in bread wheat through quantitative trait locus interval and association mapping. Mol. Breed. 2012, 29, 963–972. [Google Scholar] [CrossRef] [Green Version]
  26. Rasheed, A.; Xia, X.; Ogbonnaya, F.; Mahmood, T.; Zhang, Z.; Mujeeb-Kazi, A.; He, Z. Genome-wide association for grain morphology in synthetic hexaploid wheats using digital imaging analysis. BMC Plant Biol. 2014, 14, 128. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Williams, K.; Munkvold, J.; Sorrells, M. Comparison of digital image analysis using elliptic Fourier descriptors and major dimensions to traits seed shape in hexaploid wheat (Triticum aestivum L.). Euphytica 2012, 190, 99–116. [Google Scholar] [CrossRef]
  28. Zhang, Y.; Liu, J.; Xia, X.; He, Z. TaGS-D1, an ortholog of rice OsGS3, is associated with grain weight and grain length in common wheat. Mol. Breed. 2014, 34, 1097–1107. [Google Scholar] [CrossRef]
  29. Jia, H.; Wan, H.; Yang, S.; Zhang, Z.; Kong, Z.; Xue, S.; Zhang, L.; Ma, Z. Genetic dissection of yield-related traits in a recombinant inbred line population created using a key breeding parent in China’s wheat breeding. Theor. Appl. Genet. 2013, 126, 2123–2139. [Google Scholar] [CrossRef] [PubMed]
  30. Guan, P.; Lu, L.; Jia, L.; Kabir, M.R.; Zhang, J.; Lan, T.; Zhao, Y.; Xin, M.; Hu, Z.; Yao, Y.; et al. Global QTL Analysis Identifies Genomic Regions on Chromosomes 4A and 4B Harboring Stable Loci for Yield-Related Traits Across Different Environments in Wheat (Triticum aestivum L.). Front. Plant Sci. 2018, 9, 529. [Google Scholar] [CrossRef] [Green Version]
  31. Cuthbert, J.L.; Somers, D.J.; Brule-Babel, A.L.; Brown, P.D.; Crow, G.H. Molecular mapping of quantitative trait loci for yield and yield components in spring wheat (Triticum aestivum L.). Theor. Appl. Genet. 2008, 117, 595–608. [Google Scholar] [CrossRef] [PubMed]
  32. Mohler, V.; Albrecht, T.; Castell, A.; Diethelm, M.; Schweizer, G.; Hartl, L. Considering causal genes in the genetic dissection of kernel traits in common wheat. J. Appl. Genet. 2016, 57, 467–476. [Google Scholar] [CrossRef]
  33. Brinton, J.; Simmonds, J.; Minter, F.; Leverington-Waite, M.; Snape, J.; Uauy, C. Increased pericarp cell length underlies a major quantitative trait locus for grain weight in hexaploid wheat. New Phytol. 2017, 215, 1026–1038. [Google Scholar] [CrossRef] [Green Version]
  34. Onyemaobi, I.; Ayalew, H.; Liu, H.; Siddique, K.H.M.; Yan, G. Identification and validation of a major chromosome region for high grain number per spike under meiotic stage water stress in wheat (Triticum aestivum L.). PLoS ONE 2018, 13, e0194075. [Google Scholar] [CrossRef] [Green Version]
  35. Li, F.; Wen, W.; He, Z.; Liu, J.; Jin, H.; Cao, S.; Geng, H.; Yan, J.; Zhang, P.; Wan, Y.; et al. Genome-wide linkage mapping of yield-related traits in three Chinese bread wheat populations using high-density SNP markers. Theor. Appl. Genet. 2018, 131, 1903–1924. [Google Scholar] [CrossRef]
  36. Kumar, A.; Mantovani, E.E.; Seetan, R.; Soltani, A.; Echeverry-Solarte, M.; Jain, S.; Simsek, S.; Doehlert, D.; Alamri, M.S.; Elias, E.M.; et al. Dissection of Genetic Factors underlying Wheat Kernel Shape and Size in an Elite x Nonadapted Cross using a High Density SNP Linkage Map. Plant. Genome 2016, 9, 1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Wu, Q.H.; Chen, Y.X.; Zhou, S.H.; Fu, L.; Chen, J.J.; Xiao, Y.; Zhang, D.; Ouyang, S.H.; Zhao, X.J.; Cui, Y.; et al. High-density genetic linkage map construction and QTL mapping of grain shape and size in the wheat population Yanda1817 x Beinong6. PLoS ONE 2015, 10, e0118144. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Cui, F.; Zhao, C.; Ding, A.; Li, J.; Wang, L.; Li, X.; Bao, Y.; Li, J.; Wang, H. Construction of an integrative linkage map and QTL mapping of grain yield-related traits using three related wheat RIL populations. Theor. Appl. Genet. 2014, 127, 659–675. [Google Scholar] [CrossRef]
  39. Fan, X.; Cui, F.; Ji, J.; Zhang, W.; Zhao, X.; Liu, J.; Meng, D.; Tong, Y.; Wang, T.; Li, J. 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. [Google Scholar] [CrossRef]
  40. Liu, G.; Jia, L.; Lu, L.; Qin, D.; Zhang, J.; Guan, P.; Ni, Z.; Yao, Y.; Sun, Q.; Peng, H. Mapping QTLs of yield-related traits using RIL population derived from common wheat and Tibetan semi-wild wheat. Theor. Appl. Genet. 2014, 127, 2415–2432. [Google Scholar] [CrossRef]
  41. Liu, J.; Wu, B.; Singh, R.P.; Velu, G. QTL mapping for micronutrients concentration and yield component traits in a hexaploid wheat mapping population. J. Cereal Sci. 2019, 88, 57–64. [Google Scholar] [CrossRef] [PubMed]
  42. Quarrie, S.A.; Steed, A.; Calestani, C.; Semikhodskii, A.; Lebreton, C.; Chinoy, C.; Steele, N.; Pljevljakusic, D.; Waterman, E.; Weyen, J.; et al. A high-density genetic map of hexaploid wheat (Triticum aestivum L.) from the cross Chinese Spring x SQ1 and its use to compare QTLs for grain yield across a range of environments. Theor. Appl. Genet. 2005, 110, 865–880. [Google Scholar] [CrossRef]
  43. Hai, L.; Guo, H.; Wagner, C.; Xiao, S.; Friedt, W. Genomic regions for yield and yield parameters in Chinese winter wheat (Triticum aestivum L.) genotypes tested under varying environments correspond to QTL in widely different wheat materials. Plant Sci. 2008, 175, 226–232. [Google Scholar] [CrossRef]
  44. Wang, J.; Liu, W.; Wang, H.; Li, L.; Wu, J.; Yang, X.; Li, X.; Gao, A. QTL mapping of yield-related traits in the wheat germplasm 3228. Euphytica 2010, 177, 277–292. [Google Scholar] [CrossRef]
  45. Lee, H.S.; Jung, J.-U.; Kang, C.-S.; Heo, H.-Y.; Park, C.S. Mapping of QTL for yield and its related traits in a doubled haploid population of Korean wheat. Plant Biotechnol. Rep. 2014, 8, 443–454. [Google Scholar] [CrossRef]
  46. Roder, M.S.; Huang, X.Q.; Borner, A. Fine mapping of the region on wheat chromosome 7D controlling grain weight. Funct. Integr. Genom. 2008, 8, 79–86. [Google Scholar] [CrossRef]
  47. Tian, X.; Wen, W.; Xie, L.; Fu, L.; Xu, D.; Fu, C.; Wang, D.; Chen, X.; Xia, X.; Chen, Q.; et al. Molecular Mapping of Reduced Plant Height Gene Rht24 in Bread Wheat. Front. Plant Sci. 2017, 8, 1379. [Google Scholar] [CrossRef] [Green Version]
  48. Mo, Y.; Vanzetti, L.S.; Hale, I.; Spagnolo, E.J.; Guidobaldi, F.; Al-Oboudi, J.; Odle, N.; Pearce, S.; Helguera, M.; Dubcovsky, J. Identification and characterization of Rht25, a locus on chromosome arm 6AS affecting wheat plant height, heading time, and spike development. Theor. Appl. Genet. 2018, 131, 2021–2035. [Google Scholar] [CrossRef] [Green Version]
  49. Agarwal, P.; Balyan, H.S.; Gupta, P.K. Identification of modifiers of the plant height in wheat using an induced dwarf mutant controlled by RhtB4c allele. Physiol. Mol. Biol. Plants 2020, 26, 2283–2289. [Google Scholar] [CrossRef] [PubMed]
  50. Zhang, J.; Dell, B.; Biddulph, B.; Drake-Brockman, F.; Walker, E.; Khan, N.; Wong, D.; Hayden, M.; Appels, R. Wild-type alleles of Rht-B1 and Rht-D1 as independent determinants of thousand-grain weight and kernel number per spike in wheat. Mol. Breed. 2013, 32, 771–783. [Google Scholar] [CrossRef]
  51. Chai, L.; Chen, Z.; Bian, R.; Zhai, H.; Cheng, X.; Peng, H.; Yao, Y.; Hu, Z.; Xin, M.; Guo, W.; 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. 2019, 132, 1815–1831. [Google Scholar] [CrossRef] [Green Version]
  52. Chen, Z.; Cheng, X.; Chai, L.; Wang, Z.; Bian, R.; Li, J.; Zhao, A.; Xin, M.; Guo, W.; Hu, Z.; et al. Dissection of genetic factors underlying grain size and fine mapping of QTgw.cau-7D in common wheat (Triticum aestivum L.). Theor. Appl. Genet. 2020, 133, 149–162. [Google Scholar] [CrossRef]
  53. Sun, C.; Dong, Z.; Zhao, L.; Ren, Y.; Zhang, N.; Chen, F. The Wheat 660K SNP array demonstrates great potential for marker-assisted selection in polyploid wheat. Plant Biotechnol. J. 2020, 18, 1354–1360. [Google Scholar] [CrossRef] [PubMed]
  54. Wang, S.; Chen, J.; Zhang, W.; Hu, Y.; Chang, L.; Fang, L.; Wang, Q.; Lv, F.; Wu, H.; Si, Z.; et al. Sequence-based ultra-dense genetic and physical maps reveal structural variations of allopolyploid cotton genomes. Genome Biol. 2015, 16, 108. [Google Scholar] [CrossRef] [Green Version]
  55. Schwarzkopf, E.J.; Motamayor, J.C.; Cornejo, O.E. Genetic differentiation and intrinsic genomic features explain variation in recombination hotspots among cocoa tree populations. BMC Genom. 2020, 21, 332. [Google Scholar] [CrossRef]
  56. Cavanagh, C.R.; Chao, S.; Wang, S.; Huang, B.E.; Stephen, S.; Kiani, S.; Forrest, K.; Saintenac, C.; Brown-Guedira, G.L.; Akhunova, A.; et al. Genome-wide comparative diversity uncovers multiple targets of selection for improvement in hexaploid wheat landraces and cultivars. Proc. Natl. Acad. Sci. USA 2013, 110, 8057–8062. [Google Scholar] [CrossRef] [Green Version]
  57. Wang, S.; Wong, D.; Forrest, K.; Allen, A.; Chao, S.; Huang, B.E.; Maccaferri, M.; Salvi, S.; Milner, S.G.; Cattivelli, L.; et al. Characterization of polyploid wheat genomic diversity using a high-density 90,000 single nucleotide polymorphism array. Plant Biotechnol. J. 2014, 12, 787–796. [Google Scholar] [CrossRef] [Green Version]
  58. Sun, C.; Zhang, F.; Yan, X.; Zhang, X.; Dong, Z.; Cui, D.; Chen, F. Genome-wide association study for 13 agronomic traits reveals distribution of superior alleles in bread wheat from the Yellow and Huai Valley of China. Plant Biotechnol. J. 2017, 15, 953–969. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Winfield, M.O.; Allen, A.M.; Burridge, A.J.; Barker, G.L.; Benbow, H.R.; Wilkinson, P.A.; Coghill, J.; Waterfall, C.; Davassi, A.; Scopes, G.; et al. High-density SNP genotyping array for hexaploid wheat and its secondary and tertiary gene pool. Plant Biotechnol. J. 2016, 14, 1195–1206. [Google Scholar] [CrossRef] [PubMed]
  60. Liu, S.; Huang, S.; Zeng, Q.; Wang, X.; Yu, R.; Wang, Q.; Singh, R.P.; Bhavani, S.; Kang, Z.; Wu, J.; et al. Refined mapping of stripe rust resistance gene YrP10090 within a desirable haplotype for wheat improvement on chromosome 6A. Theor. Appl. Genet. 2021, 3, 1–17. [Google Scholar]
  61. Liu, J.; Luo, W.; Qin, N.; Ding, P.; Zhang, H.; Yang, C.; Mu, Y.; Tang, H.; Liu, Y.; Li, W.; et al. A 55 K SNP array-based genetic map and its utilization in QTL mapping for productive tiller number in common wheat. Theor. Appl. Genet. 2018, 131, 2439–2450. [Google Scholar] [CrossRef]
  62. Ma, J.; Zhang, H.; Li, S.; Zou, Y.; Li, T.; Liu, J.; Ding, P.; Mu, Y.; Tang, H.; Deng, M.; et al. Identification of quantitative trait loci for kernel traits in a wheat cultivar Chuannong16. BMC Genet. 2019, 20, 77. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Ma, J.; Ding, P.; Liu, J.; Li, T.; Zou, Y.; Habib, A.; Mu, Y.; Tang, H.; Jiang, Q.; Liu, Y.; et al. Identification and validation of a major and stably expressed QTL for spikelet number per spike in bread wheat. Theor. Appl. Genet. 2019, 132, 3155–3167. [Google Scholar] [CrossRef]
  64. Ren, T.; Hu, Y.; Tang, Y.; Li, C.; Yan, B.; Ren, Z.; Tan, F.; Tang, Z.; Fu, S.; Li, Z. Utilization of a Wheat55K SNP Array for Mapping of Major QTL for Temporal Expression of the Tiller Number. Front. Plant Sci. 2018, 9, 333. [Google Scholar] [CrossRef] [PubMed]
  65. Huang, S.; Wu, J.; Wang, X.; Mu, J.; Xu, Z.; Zeng, Q.; Liu, S.; Wang, Q.; Kang, Z.; Han, D. Utilization of the Genomewide Wheat 55K SNP Array for Genetic Analysis of Stripe Rust Resistance in Common Wheat Line P9936. Phytopathology 2019, 109, 819–827. [Google Scholar] [CrossRef] [PubMed]
  66. Rasheed, A.; Xia, X. From markers to genome-based breeding in wheat. Theor. Appl. Genet. 2019, 132, 767–784. [Google Scholar] [CrossRef]
  67. Cui, F.; Fan, X.; Chen, M.; Zhang, N.; Zhao, C.; Zhang, W.; Han, J.; Ji, J.; Zhao, X.; Yang, L.; et al. QTL detection for wheat kernel size and quality and the responses of these traits to low nitrogen stress. Theor. Appl. Genet. 2016, 129, 469–484. [Google Scholar] [CrossRef]
  68. Hao, C.; Jiao, C.; Hou, J.; Li, T.; Liu, H.; Wang, Y.; Zheng, J.; Liu, H.; Bi, Z.; Xu, F.; et al. Resequencing of 145 Landmark Cultivars Reveals Asymmetric Sub-genome Selection and Strong Founder Genotype Effects on Wheat Breeding in China. Mol. Plant 2020, 12, 1733–1751. [Google Scholar] [CrossRef]
  69. Brinton, J.; Ramirez-Gonzalez, R.H.; Simmonds, J.; Wingen, L.; Orford, S.; Griffiths, S.; Wheat Genome, P.; Haberer, G.; Spannagl, M.; Walkowiak, S.; et al. A haplotype-led approach to increase the precision of wheat breeding. Commun. Biol. 2020, 3, 712. [Google Scholar] [CrossRef]
  70. Meng, L.; Li, H.; Zhang, L.; Wang, J. QTL IciMapping: Integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental populations. Crop J. 2015, 3, 269–283. [Google Scholar] [CrossRef] [Green Version]
  71. Cui, F.; Zhang, N.; Fan, X.L.; Zhang, W.; Zhao, C.H.; Yang, L.J.; Pan, R.Q.; Chen, M.; Han, J.; Zhao, X.Q.; 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. [Google Scholar] [CrossRef] [PubMed]
  72. Wu, J.; Yu, R.; Wang, H.; Zhou, C.; Huang, S.; Jiao, H.; Yu, S.; Nie, X.; Wang, Q.; Liu, S.; et al. A large-scale genomic association analysis identifies the candidate causal genes conferring stripe rust resistance under multiple field environments. Plant Biotechnol. J. 2021, 19, 177–191. [Google Scholar] [CrossRef] [PubMed]
  73. Banta, J.A.; Stevens, M.H.H.; Pigliucci, M. A comprehensive test of the ‘limiting resources’ framework applied to plant tolerance to apical meristem damage. Oikos 2010, 119, 359–369. [Google Scholar] [CrossRef]
  74. Wu, Y.; Bhat, P.R.; Close, T.J.; Lonardi, S. Efficient and accurate construction of genetic linkage maps from the minimum spanning tree of a graph. PLoS Genet. 2008, 4, e1000212. [Google Scholar] [CrossRef] [PubMed]
  75. Yu, L.; He, F.; Chen, G.L.; Cui, F.; Li, X.F. Identification of 1BL·1RS Wheat-Rye Chromosome Translocations via 1RS Specific Molecular Markers and Genomic in situ Hybridization. Acta Agron. Sin. 2011, 37, 563–569. [Google Scholar] [CrossRef]
Figure 1. Collinearity analysis of genetic map and reference genome. NOTE: The genetic distances of the linkage group are shown as the left Y-axis, the recombination rate of bin markers as the right Y-axis, the physical location of the markers as the x-axis, the collinearity as the red scatter dots, and the recombination rate of bin markers on the reference genome as the black histogram. A, B, D are the three subgroups of common wheat.
Figure 1. Collinearity analysis of genetic map and reference genome. NOTE: The genetic distances of the linkage group are shown as the left Y-axis, the recombination rate of bin markers as the right Y-axis, the physical location of the markers as the x-axis, the collinearity as the red scatter dots, and the recombination rate of bin markers on the reference genome as the black histogram. A, B, D are the three subgroups of common wheat.
Plants 10 01167 g001
Figure 2. The linkage group corresponding to chromosome 4B. Note: (A) represents the original linkage group, and (B) the linkage group after the addition of the KASP marker.
Figure 2. The linkage group corresponding to chromosome 4B. Note: (A) represents the original linkage group, and (B) the linkage group after the addition of the KASP marker.
Plants 10 01167 g002
Table 1. Statistical analysis of parent and RIL lines for traits.
Table 1. Statistical analysis of parent and RIL lines for traits.
TraitsEnvironmentXinong1376Xiaoyan81Mean ± SDMinimumMaximump-ValueHeritability
Plant height19NY65.2577.75 **67.08 ± 13.7832.296.82.19 × 10−30.73
20NY68.2481.22 **80.03 ± 14.4340.2109.85.66 × 10−6
19YL68.3678.23 **65.78 ± 12.7834.690.96.14 × 10−4
20YL72.3383.25 **72.24 ± 15.0838.3109.23.82 × 10−2
Thousand Grain Weight19NY41.35 **36.2340.72 ± 4.3727.8152.191.12 × 10−10.62
20NY42.13 **39.4842.62 ± 4.5126.2851.761.24 × 10−3
19YL44.32 **41.7545.32 ± 4.4134.2255.052.80 × 10−2
20YL46.23 **42.3245.21 ± 4.4029.554.833.68 × 10−1
Grain length19NY7.12 **6.877.23 ± 0.376.278.043.93 × 10−20.61
20NY7.32 **6.757.14 ± 0.356.348.036.79 × 10−2
19YL7.51 **7.247.44 ± 0.346.688.231.81 × 10−2
20YL7.36 **7.147.51 ± 0.386.678.511.96 × 10−1
Grain width19NY3.313.213.37 ± 0.152.883.691.45 × 10−30.72
20NY3.88 **3.623.45 ± 0.182.813.832.18 × 10−3
19YL3.51 **3.283.60 ± 0.163.113.91.38 × 10−2
20YL3.66 **3.423.60 ± 0.163.163.953.02 × 10−3
Note: ** represents a significant difference between the two parents when p = 0.01.
Table 2. Single-nucleotide polymorphism (SNP) marker statistics about distribution and density on 21 wheat chromosomes derived from crossing between Xinong1376 and Xiaoyan81.
Table 2. Single-nucleotide polymorphism (SNP) marker statistics about distribution and density on 21 wheat chromosomes derived from crossing between Xinong1376 and Xiaoyan81.
ChromosomeLinkage GroupLength(cM)Maker NumbersBin NumberInsinuation MarkersMaximum ClearanceAverage BinBin Density
1ALG1A192.661064112104525.681.720.58
1BLG1B144.6155811844726.861.230.82
1DLG1D130.853944233618.013.120.32
2ALG2A215.9795114094023.461.540.65
2BLG2B244.4367617359725.441.410.71
2DLG2D.1132.891614815425.422.770.36
LG2D.238.6977117510.063.520.28
3ALG3A311.231322285130116.81.090.92
3BLG3B160.6148714445812.591.120.9
3DLG3D.117.463883613.712.180.46
LG3D.2214.6699987102622.842.470.41
4ALG4A228.4261412359224.851.860.54
4BLG4B169.56118519311568.570.881.14
4DLG4D166.012037719916.522.160.46
5ALG5A.1234.1896916996316.261.390.72
LG5A.252.94139391349.871.360.74
5BLG5B.168.44682886758.090.781.29
LG5B.2172.453816452915.431.050.95
5DLG5D.1223.5819211917113.691.880.53
LG5D.261.6427424158.031.470.68
6ALG6A.1112.711545013720.912.250.44
LG6A.254.951613615117.461.530.66
6BLG6B167.658521887837.70.891.12
6DLG6D.131.083473410.44.440.23
LG6D.2170.8850612449712.411.380.73
7ALG7A.175.251947617614.910.991.01
LG7A.2225.3864720163318.91.120.89
7BLG7B170.5488212984518.191.320.76
7DLG7D.1237.845313044615.511.830.55
LG7D.285.3817131624.966.570.15
1st homologous3468.122016272182826.861.720.58
2nd homologous4631.981865372176625.441.70.59
3rd homologous4703.962846524282122.841.340.74
4th homologous3563.992002393194724.851.440.7
5th homologous6813.142947621288715.431.310.76
6th homologous4537.271707405160220.911.330.75
7th homologous5794.352193549211624.961.450.69
A genome101703.6962151231607225.681.380.72
B genome81298.2358601197549026.861.080.92
D genome121510.873501708340525.442.130.47
TOTAL304512.791557631361496726.441.440.69
Table 3. Full genomic QTL mapping results of plant height and grain related traits in the F8 RIL lines between Xinong1376 and Xiaoyan81.
Table 3. Full genomic QTL mapping results of plant height and grain related traits in the F8 RIL lines between Xinong1376 and Xiaoyan81.
TraitQTLs NameEnvironmentPositionLODPVE (%)AddLeft and Right MarkerIntervalPhysical IntervalReference
PHQph.nwafu-2D.119YL1710.73.73−3.77AX-111561744/AX-17955774816.5–20.523.416254/28.417456(Zhai et al., 2016)
PH 20YL178.823.3−3.75AX-111561744/AX-17955774816.5–20.523.416254/28.417456
PHQph.nwafu-2D.219YL10353.3433.0811.29AX-94570302/AX-109998182102.5–103.5413.778968/425.474614
PHQph.nwafu-4B19NY597.3910.234.93AX-179477460/AX-11098406558.5–59.530.805339/32.961929(Mohler et al., 2016)
PH 20NY599.3613.765.9AX-179477460/AX-11098406558.5–59.530.805339/32.961929
PH 19YL5923.339.326.28AX-179477460/AX-11098406558.5–59.530.805339/32.961929
PH 20YL5926.8512.267.62AX-179477460/AX-11098406558.5–59.530.805339/32.961929
PHQph.nwafu-4D19NY6217.2727.08−7.69AX-86170701/AX-8944520161.5–62.518.781207/19.459614(Zhang et al., 2013)
PH 20NY6216.7827.09−7.93AX-86170701/AX-8944520161.5–62.518.781207/19.459614
PH 19YL6240.1720.11−8.84AX-86170701/AX-8944520161.5–62.518.781207/19.459614
PH 20YL6242.2123.71−10.16AX-86170701/AX-8944520161.5–62.518.781207/19.459614
PHQph.nwafu-5B20YL554.321.492.52AX-109908739/AX-8617461254.5–55.5422.122099/425.671678
PHQph.nwafu-5D20YL19046.9229.23−11.18AX-94390434/AX-110033637189.5–190.5466.230408/469.357817(Quarrie et al., 2005; Hai et al., 2008)
PHQph.nwafu-6B.120YL1396.192.193.07AX-109987590/AX-86162252137.5–139.5687.177084/688.20385
PHQph.nwafu-6B.219YL1604.151.362.27AX-110632551/AX-109509377159.5–160.5712.125253/711.370298
TGWQtgw.nwafu-2A20YL1862.65.080.95AX-95103231/AX-94508212185.5–186.5733.854404/734.347961(Cui et al., 2014)
TGWQtgw.nwafu-2B20YL1014.2410.141.34AX-108905289/AX-9523562695.5–106.5153.585606/568.468886(Li et al., 2018)
TGWQtgw.nwafu-3A19YL1332.74.14−0.89AX-179477407/AX-94457296132.5–134.5457.796943/431.074614
TGWQtgw.nwafu-4B20NY514.183.431.23AX-111494900/AX-9443852748.5–53.525.847125/26.491497[67]
TGW 19YL514.376.851.18AX-111494900/AX-9443852749.5–53.525.847125/26.491497
TGW 19NY522.855.021.06AX-94438527/AX-11038363448.5–55.526.491497/28.71668
TGWQtgw.nwafu-4D.119NY605.879.73−1.44AX-89703298/AX-8617070156.5–60.516.926631/18.781207(Mohler et al., 2016)
TGW 20NY606.25.25−1.48AX-89703298/AX-8617070156.5–60.516.926631/18.781207
TGWQtgw.nwafu-4D.219YL1113.545.55−1.03AX-111926032/AX-94818797107.5–112.5476.884228/477.371597
TGWQtgw.nwafu-5A19YL447.1811.941.51AX-95510385/AX-9511718843.5–45.5698.508129/702.466804
TGW 20YL443.516.991.11AX-95510385/AX-9511718843.5–45.5698.508129/702.466804
TGWQtgw.nwafu-5D.119NY375.519.321.4AX-111543112/AX-11057607434.5–38.538.070293/41.294446
TGW 20NY3714.7914.242.43AX-111543112/AX-11057607435.5–38.538.070293/41.294446
TGWQtgw.nwafu-5D.220NY466.465.56−1.51AX-111019963/AX-11008549944.5–49.542.928674/44.192407
TGWQtgw.nwafu-5D.319YL813.665.831.05AX-110867187/AX-10882729779.5–81.5369.202139/370.064947
TGWQtgw.nwafu-6A20YL293.066.26−1.06AX-109431286/AX-10935866727.5–30.5606.979733/608.046298(Cui et al., 2014)
GLQgl.nwafu-1A19YL1503.33.11−0.06AX-95682344/AX-108726119148.5–150.5572.350803/572.658176(Mir et al., 2012)
GLQgl.nwafu-1B.120YL03.624.510.09AX-94835306/AX-1794762790–0.559.471177/94.978091
GLQgl.nwafu-1B.219NY653.514.490.08AX-94650293/AX-11228850164.5–66.5640.845515/641.632325
GL 19YL655.235.140.08AX-94650293/AX-11228850164.5–65.5640.845515/641.632325
GLQgl.nwafu-3A20YL1355.717.61−0.1AX-94426283/AX-110122062134.5–136.5511.755031/510.853056
GL 20NY1373.786.88−0.09AX-179557644/AX-94387510136.5–137.5541.482465/540.048345
GLQgl.nwafu-4A20NY494.419.110.1AX-111251110/AX-17947667346.5–53.5407.389107/129.089816
GL 19NY504.576.070.1AX-111251110/AX-17947667347.5–52.5407.389107/129.089816
GLQgl.nwafu-4B.119YL514.284.170.07AX-179476673/AX-11017314047.5–52.5129.089816/140.310606
GLQgl.nwafu-4B.219YL6854.94−0.08AX-109507847/AX-10942790067.5–69.5114.952789/161.548436(Wang et al., 2010)
GL 20YL684.035.11−0.08AX-109507847/AX-10942790067.5–69.5114.952789/161.548436
GLQgl.nwafu-4D19YL162.672.550.05AX-108892806/AX-10944799715.5–18.56.598631/7.048661
GLQgl.nwafu-5A19NY446.939.590.12AX-95510385/AX-9511718843.5–44.5698.508129/698.508129
GL 19YL446.476.720.09AX-95510385/AX-9511718843.5–44.5698.508129/700.34701
GL 20YL445.87.730.1AX-95510385/AX-9511718843.5–45.5698.508129/700.34701
GLQgl.nwafu-5B19YL62.752.56−0.05AX-112288130/AX-956315255.5–6.56.654131/8.917454
GLQgl.nwafu-5D.120NY373.135.760.08AX-111543112/AX-11057607433.5–39.538.070293/41.294446
GLQgl.nwafu-5D.219YL826.135.990.08AX-111496494/AX-10970791381.5–84.5370.135626/379.028214
GL 20YL8267.750.1AX-111496494/AX-10970791381.5–84.5370.135626/379.028214
GL 19NY894.365.620.09AX-110558491/AX-11190391788.5–91.5385.893875/386.126855
GLQgl.nwafu-5D.320YL1917.459.830.11AX-110033637/AX-110830424190.5–191.5469.357817/469.523881
GLQgl.nwafu-5D.419YL2186.857.840.09AX-110777538/AX-111512534215.5–221.5485.909071/491.01105
GLQgl.nwafu-6B19YL1628.798.950.1AX-110287286/AX-111572797161.5–162.5712.125253/712.245125(Li et al., 2018)
GL 19NY1677.8810.730.13AX-89379712/AX-94499484166.5–167704.884934/718.376276
GL 20NY1673.556.680.09AX-89379712/AX-94499484166.5–167704.884934/718.376276
GL 20YL1677.7310.330.12AX-89379712/AX-94499484166.5–167704.884934/718.376276
GWQgw.nwafu-2B20YL943.866.130.04AX-109423066/AX-10899083293.5–94.5152.611396/153.128588
GWQgw.nwafu-2D19NY122.554.20.04AX-179477408/AX-11136773811.5–12.520.768547/21.405473(Huang et al., 2006; Guan et al., 2018; Wu et al., 2015)
GWQgw.nwafu-3A19YL3113.184.660.04AX-110915909/AX-110475339308.5–311746.360221/749.849798(Lee et al., 2014)
GWQgw.nwafu-4B.119NY514.156.950.05AX-111494900/AX-9443852749.5–54.525.847125/26.491497
GW 20NY513.56.850.05AX-111494900/AX-9443852748.5–54.525.847125/26.491497
GWQgw.nwafu-4B.220YL683.916.230.04AX-109507847/AX-10942790067.5–68.5114.952789/161.548436(Wang et al., 2010)
GWQgw.nwafu-4B.319YL779.2315.240.07AX-179559104/AX-9565879876.5–77.5520.214474/523.447693
GWQgw.nwafu-4D20NY603.26.32−0.05AX-89703298/AX-8617070159.5–61.516.926631/18.781207
GW 19NY614.157.22−0.05AX-86170701/AX-11057200659.5–61.518.781207/19.179341
GW 19YL637.9312.12−0.06AX-86170701/AX-8944520161.5–64.518.781207/19.459614
GW 20YL636.3710.37−0.05AX-86170701/AX-8944520161.5–63.518.781207/19.459614
GWQgw.nwafu-5D19NY1633.56.010.04AX-109317498/AX-109855976159.5–166.5448.686533/449.292436
GWQgw.nwafu-6D20NY42.515.80.05AX-111594857/AX-1094060810–12.512.650045/8.255713
Note: PH, TGW, GL, and GW represent plant height, thousand-grain weight, grain length, and grain width, respectively. Reference represents that the confidence interval of this study overlaps with that of previous studies.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lv, D.; Zhang, C.; Yv, R.; Yao, J.; Wu, J.; Song, X.; Jian, J.; Song, P.; Zhang, Z.; Han, D.; et al. Utilization of a Wheat50K SNP Microarray-Derived High-Density Genetic Map for QTL Mapping of Plant Height and Grain Traits in Wheat. Plants 2021, 10, 1167. https://doi.org/10.3390/plants10061167

AMA Style

Lv D, Zhang C, Yv R, Yao J, Wu J, Song X, Jian J, Song P, Zhang Z, Han D, et al. Utilization of a Wheat50K SNP Microarray-Derived High-Density Genetic Map for QTL Mapping of Plant Height and Grain Traits in Wheat. Plants. 2021; 10(6):1167. https://doi.org/10.3390/plants10061167

Chicago/Turabian Style

Lv, Dongyun, Chuanliang Zhang, Rui Yv, Jianxin Yao, Jianhui Wu, Xiaopeng Song, Juntao Jian, Pengbo Song, Zeyuan Zhang, Dejun Han, and et al. 2021. "Utilization of a Wheat50K SNP Microarray-Derived High-Density Genetic Map for QTL Mapping of Plant Height and Grain Traits in Wheat" Plants 10, no. 6: 1167. https://doi.org/10.3390/plants10061167

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop