QTL Mapping of Yield-Related Traits in Tetraploid Wheat Based on Wheat55K SNP Array

To enhance the understanding of yield-related traits in tetraploid wheat, it is crucial to investigate and identify genes that govern superior yield characteristics. This study utilized the wheat55K single nucleotide polymorphism array to genotype a recombinant inbred line (RIL) population consisting of 120 lines developed through the crossbreeding of two tetraploid wheat varieties, Qin Hei-1 (QH-1) and Durum Wheat (DW). An investigation and analysis were conducted on 11 yield-related traits, including peduncle length (PL), neck length (NL), spike length (SL), flowering date (FD), heading date (HD), thousand-kernel weight (TKW), kernel area ratio (KAR), kernel circumference (KC), kernel length (KL), kernel width (KW), and kernel length–width ratio (KL-WR), over a period of three years in two locations, Yang Ling, Shaanxi, and Lin He, Inner Mongolia. The analysis identified nine stable loci among eight agronomic traits, named QSL.QD-1A.1, QNL.QD-4B.2, QPL.QD-4B.1, QFD.QD-2B, QHD.QD-2B.1, QHD.QD-4B, QKC.QD-4B.2, QKL-WR.QD-4B.6, and QKL.QD-4B.2. Among them, the additive effects of three QTLs, QSL.QD-1A.1, QNL.QD-4B.2, and QFD.QD-2B, were positive, indicating that the enhancing alleles at these loci were derived from the parent line QH-1. These three QTLs showed significant positive effects on the phenotypes of the population materials. Furthermore, potential functional genes were identified within the mapping intervals of QSL.QD-1A.1 and QNL.QD-4B.2, which regulate the development of spike length and neck length, respectively. These results provide potential QTLs and candidate genes, which broaden the genetic basis of agronomic traits related to yield, such as SL, NL, PL, and FD, and benefits for wheat breeding and improvement.


Introduction
Wheat yield, a complex trait influenced by numerous genes and the interaction of environmental and genetic factors throughout plant growth stages [1], encompasses various components, including plant morphology such as peduncle length (PL), neck length (NL), and spike length (SL).Additionally, yield traits are influenced by growth stage indicators like flowering date (FD) and heading date (HD), along with grain development-related characteristics such as thousand-kernel weight (TKW), kernel area ratio (KAR), kernel circumference (KC), kernel length (KL), kernel width (KW), and kernel length-width ratio (KL-WR).Peduncle length, a component of plant height representing the distance from the panicle base to the flag leaf petiole [2], plays a pivotal role in shaping overall plant Plants 2024, 13, 1285 3 of 17 genes in spikelet growth and their impact on spikelet number per spike by comparing the genetic positions of tetraploid Triticum turgidum and hexaploid Chinese Spring [27].Liu J et al. identified twenty common genes in the physical interval between flanking markers on chromosome 4B of "Chinese Spring" and wild emmer, which regulate productive tiller number [28].Considering the phenomenon of shared genetic variation in yield-related traits between tetraploid and hexaploid wheat, and the relative simplicity of the tetraploid wheat genome, research on the localization of yield-related traits in tetraploid wheat is beneficial for accelerating the breeding process of common wheat.
QH-1, characterized as a spring wheat variant, exhibits notable phenotypic traits and resistance attributes, including a compact plant stature of approximately 90 cm, increased tillering, slender grains, and robust disease resistance.In contrast, DW, a tetraploid wheat variety, holds the second position after common wheat and accounts for approximately 5% of the global annual wheat yield.DW is extensively cultivated across Europe, North Africa, and West Asia, where it is predominantly used in local culinary practices [29,30].The plants are shorter in height, with sturdy stems that are resistant to lodging, and have larger grains.Furthermore, they exhibit strong stress resistance and high grain hardness, rendering them exceptional parental lines for genetic breeding.This study utilized a RIL population derived from two tetraploid wheat materials, "QH-1" and "DW", as parents.The investigation focused on 11 yield-related traits over a period of three years in two environmental locations, Yang Ling, Shaanxi, and Lin He, Inner Mongolia.By combining the survey results of the parents and RIL population materials with genotyping using the wheat55K SNP chip strategy, additive QTL analysis was conducted in single and multiple environments, aiming to identify the major genes or loci influencing yield-related traits in tetraploid wheat.This study conducted gene mining and material characteristic evaluation on two levels to provide a theoretical basis for genetic improvement of yield-related traits in wheat.

Phenotypic Differences and Analysis
To assess variances among materials, an examination of 11 traits within the recombinant inbred line (RIL) population comprising parental and hybrid combinations was conducted.The findings indicated a marked elevation in the lengths of PL, NL, and SL linked to yield in QH-1 in contrast to DW and KW in DW, which significantly surpassed that of QH-1, with no observable change in KL.Consequently, it is evident that TKW, KC, and KAR in DW exhibit slight superiority over those in QH-1 (Figure 1 and Table 1).Through BLUE analysis of parental and RIL populations across five environmental time points (2021-2023YL and 2022-2023LH), followed by normal distribution assessment, it was observed that all 11 yield-associated traits displayed a continuous distribution and transgressive segregation, indicative of quantitative trait attributes governed by multiple genes (Table 1 and Figure S1).Examination of population phenotypes via H 2 analysis unveiled that FD, HD, and KAR are prone to environmental influences, whereas the remaining traits are predominantly determined by genetic factors (Table 1).Through BLUE analysis of parental and RIL populations across five environmental time points (2021-2023YL and 2022-2023LH), followed by normal distribution assessment, it was observed that all 11 yield-associated traits displayed a continuous distribution and transgressive segregation, indicative of quantitative trait attributes governed by multiple genes (Table 1 and Figure S1).Examination of population phenotypes via H 2 analysis unveiled that FD, HD, and KAR are prone to environmental influences, whereas the remaining traits are predominantly determined by genetic factors (Table 1).

Correlation Analyses among Different Traits
The phenotypic correlations among the 11 agronomic traits, as determined by BLUE values, are visually depicted in Figure 2. The analysis indicates notably strong positive correlations (p < 0.01) between PL and NL; HD and FD; TKW and SL, NL, PL; KC and SL, TKW; KL-WR and FD, HD, KC; KL and SL, TKW, KC, KL-WR; KW and SL, NL, PL, TKW, KC, KL.Notably, the correlation coefficient between KC and KL is the most pronounced at 0.98, followed by 0.94 between PL and NL.Furthermore, a significant negative correlation (p < 0.01) is observed between KL-WR and NL, PL, and KW, with the strongest correlation coefficient being −0.49.Collectively, the holistic evaluation of yield traits indicates a marked positive correlation among the majority of traits within the population, exerting a supportive role in the developmental trajectory.
Plants 2024, 13, 1285 6 of 17 KC, KL.Notably, the correlation coefficient between KC and KL is the most pronounced at 0.98, followed by 0.94 between PL and NL.Furthermore, a significant negative correlation (p < 0.01) is observed between KL-WR and NL, PL, and KW, with the strongest correlation coefficient being −0.49.Collectively, the holistic evaluation of yield traits indicates a marked positive correlation among the majority of traits within the population, exerting a supportive role in the developmental trajectory.

Genetic Map Construction
Genotyping of all 120 families of RILs was performed using the Wheat 55K SNP chip.Following thorough filtration through software and manual curation, a total of 11,704 pure and polymorphic SNP markers differing between the parental lines were identified after excluding SNPs with unsuccessful genotype calls.Subsequently, employing the BIN function of IciMapping V4.1, measures were taken to eliminate missing values, segregational distortions surpassing the predefined threshold, and redundant markers, resulting in a refined collection of 913 high-quality SNP markers.The constructed genetic map spans a total length of 2141.26 cM, showcasing an average marker density of 2.35 cM per SNP.Through a comparative analysis of the probe information in the SUM file of the mapping outcomes with the genomic data of Chinese Spring wheat, a subset of 848 informative probes residing on the A/B chromosome groups was identified, distributed across 26 linkage groups spanning 14 chromosomes (Table S1).

QTL Mapping Analysis
The ICIM-BIP approach was employed to identify additive QTLs associated with the 11 specified traits, leading to the detection of a total of 100 additive QTLs as presented in Table S2.These QTLs exhibited a range of LOD values from 2.50 to 26.71 and were located on 14 chromosomes, explaining 0.23% to 55.42% of the phenotypic variance (Table S2).Notably, within the identified additive QTLs, a subset of 41 were categorized as major QTLs,  2 and Table S2, Figure 3).

Combined QTL-Environment Interaction Analysis
The MET-Add method was used for cQTL analysis, and a total of 259 additive cQTLs for SL (22), NL (26), PL (21), FD (10), HD (16), TKW (28), KAR (18), KC (27), KL-WR (42), KL (28), and KW (21) were mapped (Table S4).Most of these cQTLs were minor QTLs.Moreover, all of the nine stable QTL detected using the ICIM-BIP method were similarly detected using the MET-Add method (Tables S2 and S3), showing that they might be stable in expression and less susceptible to environmental influences.Among the nine stable QTLs located on eight traits, the proportion of phenotypic variance explained (PVE) by the additive genetic effect (PVE (A)) is higher than PVE (A/E) for six QTLs associated with SL, PL, HD, KC, and KL traits.This indicates that the additive effects of these six stable major loci have a greater impact on the phenotypic variation compared to environmental effects, suggesting that genetic factors play a dominant role in phenotypic variation.Therefore, these loci can be utilized for genetic improvement to effectively enhance the quality traits.On the other hand, for NL, FD, and KL-WR traits, the PVE (A) is lower than PVE (A/E) for the stable major loci associated with these traits.This suggests that environmental effects have a greater influence on phenotypic variation compared to additive effects, indicating that environmental factors dominate the phenotypic variation for these traits.
In addition, the analysis results of the ICIM-BIP method and the MET-Add method both show that the additive effect values of QSL.QD-1A.1,QNL.QD-4B.2,QPL.QD-4B.1, and QFD.QD-2B in the nine stable QTLs are positive, indicating that the enhancing alleles at these four loci originate from the parent QH-1.The additive effect values of QHD.QD-2B.1,QHD.QD-4B, QKC.QD-4B.2,QKL-WR.QD-4B.6,and QKL.QD-4B.2 are negative, indicating that the enhancing alleles at these five loci originate from the parent DW.The underlined QTLs indicate that they are the stable QTLs and the bold font identifies that they are the major QTLs in this study.Env: Environments.

Combined QTL-Environment Interaction Analysis
The MET-Add method was used for cQTL analysis, and a total of 259 additive cQTLs for SL (22), NL (26), PL (21), FD (10), HD (16), TKW (28), KAR (18), KC (27), KL-WR (42), KL (28), and KW (21) were mapped (Table S4).Most of these cQTLs were minor QTLs.Moreover, all of the nine stable QTL detected using the ICIM-BIP method were similarly detected using the MET-Add method (Tables S2 and S3), showing that they might be stable in expression and less susceptible to environmental influences.Among the nine stable QTLs located on eight traits, the proportion of phenotypic variance explained (PVE) by the additive genetic effect (PVE (A)) is higher than PVE (A/E) for six QTLs associated with SL, PL, HD, KC, and KL traits.This indicates that the additive effects of these six stable major loci have a greater impact on the phenotypic variation compared to environmental effects, suggesting that genetic factors play a dominant role in phenotypic variation.Therefore, these loci can be utilized for genetic improvement to effectively enhance the quality traits.On the other hand, for NL, FD, and KL-WR traits, the PVE (A) is lower than PVE (A/E) for the stable major loci associated with these traits.This suggests that environmental effects have a greater influence on phenotypic variation compared to additive effects, indicating that environmental factors dominate the phenotypic variation for these traits.

Analysis of the Effects of Stable Major Loci
To assess the impact of different genotypes on the associated phenotypes, families with predicted homozygous genotypes at relevant loci were identified using genotyping data generated by the wheat 55K SNP chip and genotypic values predicted by IciMapping.Families exhibiting plausible homozygous genotypes for the respective loci were selected.These families were then categorized into the Qin Hei (QH) type and Durum Wheat (DW) type based on the sources of homozygous alleles at the corresponding QTL loci within the population.Subsequently, the phenotypic traits linked to stable loci in the two genotype groups were subjected to analysis.
As can be seen from Figure 4, the genotypic effects predicted by IciMapping for the eight traits characterized by stable loci show consistently higher values in the QH type lines compared to the DW type.This observation suggests that the allelic genotype of the parental line QH-1 exerts an additive influence on these eight traits.Specifically, in the yield-related trait SL, noteworthy variations between the two types were apparent at environmental points 2023LH and 2022YL, with the QH type exhibiting superior performance by 6.55% and 9.97% (p < 0.01), respectively.For the NL trait, substantial differences between the two types were observed at environmental points 2023LH and 2022YL, with the QH type surpassing the DW type by 54.7% and 27.8% (p < 0.01), respectively.Furthermore, in the FD trait, significant differences between the two types were detected at environmental points 2023YL and 2022YL, with the QH type outperforming the DW type by 7.11% and 11.77% (p < 0.01), respectively.These three genotypes with notable or highly notable differences exhibit positive additive effects attributed to QH-1, corroborating the notion that the enhancing alleles at these three loci stem from the parent QH-1, as evidenced by the positive additive effect values of stable QTLs QSL.QD-1A.1,QNL.QD-4B.2, and QFD.QD-2B in the aforementioned investigation.
performance by 6.55% and 9.97% (p < 0.01), respectively.For the NL trait, substantial differences between the two types were observed at environmental points 2023LH and 2022YL, with the QH type surpassing the DW type by 54.7% and 27.8% (p < 0.01), respectively.Furthermore, in the FD trait, significant differences between the two types were detected at environmental points 2023YL and 2022YL, with the QH type outperforming the DW type by 7.11% and 11.77% (p < 0.01), respectively.These three genotypes with notable or highly notable differences exhibit positive additive effects attributed to QH-1, corroborating the notion that the enhancing alleles at these three loci stem from the parent QH-1, as evidenced by the positive additive effect values of stable QTLs QSL.QD-1A.1,QNL.QD-4B.2, and QFD.QD-2B in the aforementioned investigation.

Analysis and Discussion of Enrichment of QTL Clusters in Agronomic Traits and Correlation among Agronomic Traits
In most crop localization studies, QTLs related to some agronomic traits are clustered on chromosomes, indicating that genes controlling multiple traits may be located in the same marker interval or adjacent regions.Researchers refer to this phenomenon as QTL clusters [31,32].In addition, researchers have found a direct relationship between the enrichment of agronomic traits in QTL clusters and the traits themselves [33,34].In this study, through the analysis of QTL positioning results in a single environment, it was found that a total of 20 QTL clusters were located on six chromosomes (Table S6).These QTL clusters collectively contain 51 QTLs, with 7 stable QTLs distributed in 5 QTL clusters, accounting for about 77.8% of the stable QTLs detected by the ICIM-BIP method.In-depth analysis of these five QTL clusters revealed that traits such as FD and HD in cluster C5 belong to phenological traits with a significantly high correlation coefficient between them, but their correlation with the plant type trait SL in the same QTL cluster is not high.Cluster C18 contains four QTLs for FD, KW, NL, and PL, respectively.Combined with the analysis of the correlation between the corresponding agronomic traits, it was found that FD and KW are significantly negatively correlated with a correlation coefficient

Analysis and Discussion of Enrichment of QTL Clusters in Agronomic Traits and Correlation among Agronomic Traits
In most crop localization studies, QTLs related to some agronomic traits are clustered on chromosomes, indicating that genes controlling multiple traits may be located in the same marker interval or adjacent regions.Researchers refer to this phenomenon as QTL clusters [31,32].In addition, researchers have found a direct relationship between the enrichment of agronomic traits in QTL clusters and the traits themselves [33,34].In this study, through the analysis of QTL positioning results in a single environment, it was found that a total of 20 QTL clusters were located on six chromosomes (Table S6).These QTL clusters collectively contain 51 QTLs, with 7 stable QTLs distributed in 5 QTL clusters, accounting for about 77.8% of the stable QTLs detected by the ICIM-BIP method.In-depth analysis of these five QTL clusters revealed that traits such as FD and HD in cluster C5 belong to phenological traits with a significantly high correlation coefficient between them, but their correlation with the plant type trait SL in the same QTL cluster is not high.Cluster C18 contains four QTLs for FD, KW, NL, and PL, respectively.Combined with the analysis of the correlation between the corresponding agronomic traits, it was found that FD and KW are significantly negatively correlated with a correlation coefficient of −0.19; FD has no significant correlation with NL or PL; NL, with PL and KW, shows a significant positive correlation, with correlation coefficients of 0.94 and 0.34, respectively.Similarly, through the comprehensive analysis of the other three QTL clusters, it was found that there are singlelocus multiple-trait effects within the QTL clusters, and the performance of agronomic traits in the field cannot be directly correlated.The magnitude of the correlation coefficient of agronomic traits and whether they are positively or negatively correlated can only be used as a reference but cannot directly determine whether the traits will appear in the same QTL cluster.

Exploration of the Application of Different Types of Arrays in Wheat Materials with Different Ploidies
The application of marker-assisted selection (MAS) at the whole-genome level has accelerated the breeding process of wheat.Subsequently, there have been advancements ranging from PCR-based molecular markers and extensive coverage provided by Genotypingby-Sequencing (GBS) offered by NGS, to flexible and diverse SNP arrays.Of course, within these methods, there are certain advantages and limitations.For example, PCR-based molecular markers exhibit lower throughput and density, and many wheat traits are often influenced by multiple genes, leading to poor detection of minor genes by functional markers.Although NGS can significantly improve coverage and reduce time and cost, the poor detection ability of GBS, especially in heterozygous or early populations, is particularly evident.In addition, GBS is a time-consuming and complex process in terms of library construction, data analysis, and storage.On the other hand, SNP arrays offer great flexibility in sample and data point customizations.There are currently 9 K, 15 K, 35 K, 55 K, 90 K, 820 K, and 660 K arrays available for wheat.In the determination of SNP density the chromosomes were divided into 1 MB windows for calculation.It was observed that almost all bins in the Wheat 660 K SNP array contained SNPs, whereas the other four arrays detected many bins without SNPs (especially the Wheat 35 K, 55 K, and 90 K SNP arrays).In comparison to the Wheat 820 K SNP array, more SNPs were found at the top or end of the chromosomes in the Wheat 660 K SNP array, indicating a higher frequency of recombination events in these regions compared to the centromeric regions [35].Thus, it is suggested that a lower coverage array such as the 55 K array could satisfy the purpose of mapping loci in a small-sized population [36].The high-throughput SNP arrays are more suited for larger-sized fne-mapping populations or for detecting as many polymorphism loci as possible for wheat-related species in Triticeae such as Agropyron [37].
Currently, the application of the 55 K and 660 K arrays is primarily focused on hexaploid wheat, with relatively limited application in tetraploid wheat.One of the reasons for this is the reduced number and density of applied markers, which subsequently impacts the length of the genetic map and the size of the localization interval.In a study by Jiajun Liu et al., utilizing the 55 K array in common wheat (Triticum aestivum L.), a genetic map of an RIL population with 2524 bin markers was mapped on 34 linkage maps for the 21 chromosomes of wheat.The genetic map length was found to be 3021.04cM, with an average of 1.20 cM per marker [28].In a study by Ziqiang Mo, utilizing a RILs population created from tetraploid wheat materials, 150 bin markers were mapped onto 15 linkage groups across the 14 chromosomes of tetraploid wheat, with a total length of 2411.8 cM and an average marker density of 2.10 cM per bin [27].In this study, 486 effective probes were obtained from a RIL population created from tetraploid materials using the 55 K array.The genetic map generated from these effective probes spanned 21 linkage groups, with a total genetic distance of 4397.21 cM.The average marker density between adjacent markers was 9.05 cM.Another factor contributing to the limited application on tetraploid wheat is the interference of D genome group probes during database alignment, which may lead to the inclusion of some unreliable localization intervals.These intervals often exhibit discrepancies between the information of neighboring probes, necessitating their exclusion.

Candidate Gene Prediction Analysis within the QTL Mapping Interval
The stable locus QSL.QD-1A.1 on the spike length trait is located on 543283549 to 566194763 bp in the Durum Wheat (cv.Svevo) genome.Mengistu et al. identified QTL related to yield in this interval [38].In this mapping interval, M. Graziani et al. identified QTLs regulating grain weight and test weight [39].Within this interval, there is a potential gene, TRITD1Av1G218480, corresponding to Chinese Spring v2.1 genome TraesCS1A03G1004700, which belongs to the AP2/B3 transcription factor family protein and plays a role as an ethylene-responsive transcription factor.Wang Yuange et al. identified DUO-B1, encoding an AP2/ERF transcription factor, using wheat as a model.This transcription factor regulates the spikelet structure of bread wheat.Furthermore, mutations in DUO-B1 result in slightly more spikelets and increased grain number per spike [40].In this interval, there is also a gene, TraesCS1A02G413800, which is functionally annotated as DELLA protein GAI.The function of DELLA protein in dwarfing genes is well understood, and the plant height also has a direct impact on spike length.NL is part of PL and the peduncle is part of plant height (PH); they are significantly influenced by the GA metabolism and signaling [3,7] and thus tightly linked to the PH genes.Recently, a PL QTL independent of PH was identified on 5AS [8], providing a chance to modulate the PL.In this study, QNL.QD-4B.2 and QPL.QD-4B.1 are located in the same interval.In this interval, we speculate that the genes TRITD4Bv1G012860 or TRITD4Bv1G012960 may play a regulatory role in PL and NL.These genes have a Vps4 function.Vps4 is a key component required for ESCRT-III disassembly and endosomal vesicle formation.The plant-specific ESCRT component, known as PROS, is a positive regulator of Vps4 ATPase activity.PROS interacts with Vps4 and the Vps4 positive regulator LIP5/VTA1.PROS is expressed in various tissues and cell types in Arabidopsis, and its silencing results in reduced cell proliferation and abnormal organ growth [41,42]. of flowering time with favorable conditions for seed reproduction plays a vital role in plant species survival and is crucial for enhancing crop yield and global food security.The developmental process of wheat is intricate and governed by multiple genes categorized based on their response to vernalization (Vrn) and photoperiod (Ppd), or their influence on earliness per se (Eps).These genes orchestrate flowering time by regulating variations in the ultimate leaf count, phyllochron, and the duration of the flag leaf-anthesis phase.Notably, this study identified a previously documented locus associated with flowering time within the mapping region of QFD.QD-2B [43,44].There is a corresponding region identified by previous researchers within the mapping interval of QHD.QD-2B.1 [45].Within the mapping interval of QHD.QD-4B, there is a corresponding region where previous researchers have identified a similar mapping interval [44].
KC and KL are located in the same interval (653.1 MB-654.8MB, Durum Wheat reference genome), where known QTLs regulating grain weight are present [46].The aspect ratio of the kernel is often directly related to the length of the kernel and the weight of the kernel.Within this range, there is a mapping interval for the aspect ratio of the kernel, as studied by Peleg and others [47,48].In previous research reports, the interaction module between miR1432 and acyl-CoA synthetase significantly increased kernel weight by improving the grain filling rate, thereby increasing total grain yield [49].The ZmCesA5 gene in corn encodes a glycosyltransferase protein.Overexpression of the ZmCesA5 gene results in heavier grains, while its deletion affects the metabolism of starch and sucrose, leading to a reduction in corn grain weight [50].Therefore, we speculate that within the QKC.QD-4B.2 and QKL.QD-4B.2loci, the functional annotation of TRITD4Bv1G199600 as fatty acyl-CoA reductase may play a role in regulating KL and KC.Additionally, within the QKL-WR.QD-4B.6 locus, the functional annotation of TRITD4Bv1G017750 as glycosyltransferase may influence the KL-WR.

Plant Materials and Field Trials
In this study, two tetraploid wheat cultivars, namely "Qin Hei-1" (QH-1) and "Durum Wheat T6" (DW), were selected as the parental lines for crossbreeding.A population of 120 genotypes of recombinant inbred lines (RILs) was generated through the single seed descent method and employed as the experimental subjects.
Throughout the wheat growing seasons in Yang Ling, Shaanxi (108 • 07 ′ E, 34 • 30 ′ N) during 2020-2021 (YL21), 2021-2022 (YL22), and 2022-2023 (YL23), as well as in Lin He, Inner Mongolia (107 • 44 ′ E, 44 • 17 ′ N) during the 2021-2022 (LH22) and 2022-2023 (LH23) seasons, RILs along with the parental lines were cultivated for this study.Before sowing, the experimental fields were plowed, leveled, and fertilized.Throughout the growth phase of the plants, standard field management practices such as weeding and winter irrigation were implemented.The population was randomly segregated into three distinct small plots situated in various sections of the experimental field, with each small plot containing one replicate.Each line was sown in a single row, with a row length of 1.5 m, plant spacing set at 10 cm, and row spacing at 25 cm.

Phenotypic Evaluation and Statistical Analysis
Eleven significant agronomic traits related to yield were assessed across parental and population materials from diverse perspectives.These traits encompassed three plant type-related agronomic characteristics that were manually assessed post-flowering, namely, peduncle length (PL), neck length (NL), and spike length (SL).Specifically, PL denotes the distance from the spike's base to the initial internode, NL represents the measurement from the spike's base to the flag leaf ear, while SL indicates the length of the spike itself.Additionally, two traits linked to the growth phase, namely, flowering date (FD) and heading date (HD), were examined.For each trait, five plant materials were selected in each replicate for subsequent statistical analysis.After harvest, six grain-related agronomic traits were quantified using an SC-G type automatic grain analyzer from Hangzhou Wanshen Technology Co., Ltd.,China (Hangzhou, China).These traits comprised thousandkernel weight (TKW), kernel area ratio (KAR), kernel circumference (KC), kernel length (KL), kernel width (KW), and kernel length-width ratio (KL-WR).Each line material in the population was measured five times per replicate for statistical analysis.For data processing, we referred to the method proposed by Ma et al. [22].

Genotyping and Genetic Map Construction
The genomic DNA (gDNA) of each RIL and the wheat 55K SNP array work were entrusted to Beijing CapitalBio Corporation (https://www.capitalbiotech.com/,accessed on 8 September 2022) for execution.
After genetic typing results were obtained, a comparison of SNP markers was conducted between the parental entities, leading to the elimination of failed markers and the identification of polymorphic markers.Additionally, SNP markers showing heterozygous genotypes and allele frequencies less than 0.3 or greater than 0.7 between the parents were further excluded.Considering that the RIL population was derived from the tetraploid wheat background and the reference was the 55K chip based on the hexaploid wheat background, the probe positions in the BIN template file were all marked as "0" in the Anchor table using the QTL IciMapping V4.1 (https://isbreeding.caas.cn/,accessed on 26 September 2023) mapping software.The redundant markers were then removed using the "By Missing" command to obtain BIN markers for constructing the genetic map.In the MAP module, a LOD value of 7.5 was selected, and the nnTwoOpt algorithm and SARF standard were used for marker sorting and rippling.The output sum file from MAP was used for mapping and labeling screening, with reference to the Chinese Spring wheat reference genome (https://urgi.versailles.inrae.fr/download/iwgsc/IWGSC_RefSeq_Assemblies/v2.1/,accessed on 12 October 2023).By comparing and filtering, erroneous values were removed or replaced, and the positioning map was ultimately limited to only the 14 chromosomes of the A and B chromosome groups.The genetic distances were converted from recombination rates to centimorgans (cM) using the Kosambi mapping function, and a genetic linkage map was drawn using MapChart V2.3.

Quantitative Trait Loci Mapping
QTL analysis was undertaken for the 11 agronomic traits using IciMappingV4.1 [22,23,27].Among the detected additive QTLs, those with LOD > 3 and that explained phenotypic variation > 10% were defined as major QTLs.Additive QTLs that can be detected in two or more individual environments were also considered stable QTLs.The multi-environment trials (MET-Add) module of IcimappingV4.1 was used to identify cQTL (combined quantitative trait loci), and epistatic analysis was performed using the IciMappingV4.1 MET-EPI module and eQTL (epistatic quantitative trait loci) detection [22].QTLs were named based on the International Rules of Genetic Nomenclature, where "QD" is the abbreviation of the parents QH-1 and DW.
For flanking sequences of key or consistent QTL data processing, we referred to the method proposed by Ma et al. [22].Subsequently, Uniprot was employed for indepth annotation and functional scrutiny of the candidate genes within the QTL regions (https://www.uniprot.org/,accessed on 25 October 2023) [22].

Figure 2 .
Figure 2. Correlation analysis of different traits based on BLUE.* Significance at the 0.05 probability level; ** Significance at the 0.01 probability level; *** Significance at the 0.001 probability level; **** Significance at the 0.0001 probability level.The colored scale on the right represents the values of correlation coefficients.

Figure 2 .
Figure 2. Correlation analysis of different traits based on BLUE.* Significance at the 0.05 probability level; ** Significance at the 0.01 probability level; *** Significance at the 0.001 probability level; **** Significance at the 0.0001 probability level.The colored scale on the right represents the values of correlation coefficients.

Figure 3 .
Figure 3. Chromosome positions of the major or stable additive QTLs for the 11 traits described.The underlined QTLs indicate that they are the stable QTLs identified in this study.Different slash types on the QTLs represent QTL clusters at different locations.The colored marker names and bars to the right of each linkage group indicate the location of the stable QTLs or QTL clusters.Different colors are used to highlight and distinguish different traits on the same chromosome.

Figure 4 .
Figure 4. Effects of different allele genotypes of stable QTLs on correlated traits.* Significance at the 0.05 probability level; ** Significance at the 0.01 probability level.

Table 1 .
Phenotypic variation and heritability of characters in different environments in parents and RILs.

Table 1 .
Phenotypic variation and heritability of characters in different environments in parents and RILs.

Table 2 .
Stable QTLs for the quality-related trait under single-environment analysis.