J Genomics 2022; 10:8-15. doi:10.7150/jgen.67253 This volume

Research Paper

Genetic Dissection of the Mixing Properties of Wheat Flour (Triticum aestivum L.) Using Unconditional and Conditional QTL Mapping

Haixia Yu1#, Yuling An1#, Aiping Wang2, Xin Guan1, Jichun Tian1, Tangyuan Ning1, Kexin Fan1, Hao Li1, Qianqian Liu1, Dongxue Wang1, Jiansheng Chen1 Corresponding address

1. State Key Laboratory of Crop Biology/Key Laboratory of Crop Water Physiology and Drought-tolerance Germplasm Improvement, Ministry of Agriculture/Group of Wheat Quality Breeding, Shandong Agricultural University, Tai'an 271018, P.R. China.
2. Dezhou Agricultural Protection and Technological Extension Center, Dezhou, 253000, P.R. China.
# indicates the authors who contributed equally to this study

This is an open access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/). See http://ivyspring.com/terms for full terms and conditions.
Citation:
Yu H, An Y, Wang A, Guan X, Tian J, Ning T, Fan K, Li H, Liu Q, Wang D, Chen J. Genetic Dissection of the Mixing Properties of Wheat Flour (Triticum aestivum L.) Using Unconditional and Conditional QTL Mapping. J Genomics 2022; 10:8-15. doi:10.7150/jgen.67253. Available from https://www.jgenomics.com/v10p0008.htm

File import instruction

Abstract

Graphic abstract

Wheat (Triticum aestivum L.) flour mixing properties are essential quality parameters in the dough development process. Limited research on superior alleles for mixing properties has restricted their molecular improvement, and other factors related to the complex traits have been ignored. A molecular map of 9576 polymorphic markers in the RIL population (F8:9) (Shannong01-35/Gaocheng9411) was constructed to evaluate mixing property effects in three environments. The parents were selected with markedly distinct high-molecular-weight glutenin subunits (HMW-GS). This study not only evaluated mixing properties using conventional unconditional QTL mapping but also evaluated the relationships between protein-related traits using conditional QTL mapping. The analyses identified most additive QTLs for major mixing properties on chromosomes 1A, 1B, and 1D. Two major loci (1A.1-15 and 1D-1) associated with mixing properties have confirmed the important contributions of Glu-A1 and Glu-D1 to wheat quality at the QTL level, which were mainly affected by the gluten index. Another important locus, 1B.1-24 (associated with midline peak value and midline peak width, with high phenotypic variations explained), might represent a new variation distinct from Glu-B1. The favored alleles came from Gaocheng9411. Several mixing properties shared the same QTLs (1B.1-6 and 1A.1-15), indicating tight linkage or pleiotropism. Genotype-by-environment (G×E) interactions were also investigated in the present study. The QTL results in our study may improve our understanding of the genetic interrelationships between mixing properties and protein-related traits.

Keywords: flour mixing properties, unconditional QTL mapping, conditional QTL mapping, HMW-GS

Introduction

Wheat flour mixing properties have been widely applied to test the rheological characteristics of food end products, especially in the baking industry [1, 2]. They are typically quantitative traits controlled by multiple genes and susceptible to environmental factors. QTLs (quantitative trait loci) for mixing properties were identified on chromosomes 1B, 1D, 3B, 4D, and 5D using a doubled haploid (DH) population [3]. Nelson et al. (2006) found QTLs for mixograph properties with strong effects on chromosomes 1AS, 1BS and 6DS [4]. QTLs for nine mixing properties were identified on chromosomes 1B, 2A, 2B, 3A, 5B, 5D, 6B and 6D [5]. QTLs for mixing time and mixing 8 min width were identified on chromosomes 1B and 1D from 240 RIL lines (F5:6) of PH82-2/Neixiang 188 [6]. Additionally, QTLs for midline peak time, mixing peak value and midline peak width were identified on chromosomes 1A, 1B, 1D, 2A, 6D and 7D [7].

To date, little information on QTL mapping for mixing properties is available, and larger differences among previous reports were mostly caused by genetic populations or maps. Conventional QTL studies do not actually evaluate closely related traits and the genetic contribution to a single trait. Protein content and quality are strongly impacted by mixing properties. Zhu and Keightley (1995) proposed the conditional genetic analysis method [8]. It could study the contributions of single trait to complex traits and the genetic relationships among more than two traits. This method has been widely used to investigate the genetic basis of complex traits in many crops [9-12]. In wheat, conditional QTL analyses have also been conducted on plant height [13], kernel weight [14], and protein content [15]. At present, few conditional QTL analyses have evaluated the genetic relationships between mixing properties and protein-related traits.

In the present study, both unconditional and conditional QTL mapping were conducted to investigate the underlying genetic basis of mixing properties and to dissect the genetic relationships between these protein-related traits at the QTL level. The results could contribute to marker-assisted breeding for wheat quality, especially for mixing properties.

Materials and methods

Plant materials

The RIL population (F8:9) with 173 lines was generated by the cross of Shannong01-35 (SN01-35) and Gaocheng9411 (GC9411) [16]. The parents were selected with marked distinct high-molecular-weight glutenin subunits (HMW-GS).

Field trials

The RIL population was harvested at Tai'an (116°36E, 36°57N) in 2009 (E1), 2010 (E2) and Suzhou (116°58E, 33°38N) in 2010 (E3). The materials were planted in a randomized complete block design with two replicates in each environment. Every pot had three rows, and each row was two meters long with a row-to-row distance of 21 centimeters. Crop management was performed according to local practices. The average value of the mixing properties in the above environments was also calculated (AE).

Measurement of mixing properties and protein-related traits

Harvested seeds were grinded into powder by a Buhler experimental mill (Buhler, Buhler-Miag Co., Germany). The flours were then sealed and kept at 8-10℃ for analysis. Mixing properties were determined by Mixograph (National Mfg. Co., Lincoln, NE) using a 10 g sample of flour system according to AACC-54-40. The mixograph curves (two envelopes and one midline) were computed with Mixsmart software (National Mfg. Co.). Parameters were recorded as follows: MPV (midline peak value), MPW (midline peak width), MPT (midline peak time), EW (end width), MTxV (midline time-8 min value), IT (integral of curve tail), CTV (curve tail value), CTW (curve tail width), SL (left slope of peak) and SR.

Traits PC (protein content), SV (sedimentation volume), WG (wet gluten), DG (dry gluten), and GI (gluten index) were measured. The measurement methods and phenotypic values among the above traits could be viewed in a previous study [16].

Genotyping and linkage map construction

A total of 9576 polymorphic molecular markers, including 9072 SNPs [17, 18], 442 DArT markers (http://www.triticarte.com.au), 59 SSR markers, 2 biochemical markers and 1 CAPS marker, were used to genotype the RIL population. Genetic linkage groups were constructed using Mapmaker/Exp 3.0. The threshold for marker loci assignment was an LOD score of ≥3.0. The linkage map was ultimately generated using MapChart 2.1 software [19]. The genetic map spanned approximately 4048.5 cM, and the average interval of markers was 8.13 cM [16].

Statistical analysis and QTL mapping

Phenotypic values were analyzed using SPSS 17.0 software. QTL analyses were conducted by QTL Network 2.0 using the mixed-linear model approach. QTLs were designated as follows: Q, followed by the abbreviated trait name, then the wheat chromosome on which the QTL was identified; the last number indicated the marker interval on the chromosome. For example, QMPV1B.1-6 indicated a QTL for MPV identified on the sixth interval of chromosome 1B.1.

Results

Phenotypic analysis of mixing properties

The phenotypic values were listed in Table 1. The two parents had marked differences in MPT, CTV, MTxV, CTW and EW. GC9411 exhibited higher parameter values than SN01-35 in all environments. Strong transgressive segregation was detected. This phenomenon may be affected by differences in HMW-GS subunits between parents. Differences in the CV (coefficients of variation) were found for all the traits in all three environments. Among these, CTW possessed the highest CV of 80.32% and a mean of 41.86%, followed by EW, SL and SR, with CV ranges of 37.83-49.27%, 39.46-58.24% and 26.82-36.97%, respectively. The CVs of MPV and MTxV were less than 10%, indicating that these three parameters had little potential for improvement. According to the analysis, the traits followed a normal distribution and were suitable for QTL analysis.

QTLs for mixing properties

Additive QTLs and genotype-by-environment interactions (G×E) for mixing properties MPV, CTV, and MTxV

Four QTLs for MPV were detected (primarily on 1B) with net additive effect values of 0.83-2.17%, and the phenotypic variation explained (PVE) value of 6.9-29.6% (Table 2, Figure 1). The effect value of QMPV1B.1-4 was obviously changed when conditioned on SV, GI and PC, indicating that QMPV1B.1-4 was partially influenced by SV and GI and PC; however, it was not detected when conditioned on WG and DG, indicating that it was completely influenced by them. The PVE of QMPV1B.1-6 (PVE of 26.1%) was changed when conditioned on SV (PVE of 9%), but it showed no significant differences when conditioned on PC (PVE of 26.2%), GI (PVE of 25.8%) or DG (PVE of 25.6%), indicating that QMPV1B.1-6 was not influenced by PC, GI or DG but was partially influenced by SV. Notably, QMPV1B.1-24 was identified when conditioned on WG and DG in three environments, with PVEs of 8.9-27.1%. The positive alleles of the abovementioned QTL were derived from GC9411.

Six QTLs on chromosomes 1A, 1B, 3A and 5A were identified to associate with CTV, with PVEs of 1.1-21.4%. QCTV1A.1-15 was not identified when conditioned on SV, GI, WG or DG, indicating that QCTV1A.1-15 was completely influenced by the above related traits. QCTV1B.1-31 showed positive additive effects with a higher PVE of 17.3%, indicating that favorable alleles were contributed by GC9411. QCTV1B.1-26 was only identified when conditioned on PC, WG and DG, indicating that this QTL was detected when the influence of the mentioned traits was excluded. QCTV5A.1-5 was only identified when conditioned on GI, with a PVE of 5.2%, indicating that QCTV5A.1-5 was identified when the influence of GI was excluded. The favorable alleles of all the detected QTLs for CTV were contributed by GC9411.

Five QTLs associated with MTxV were identified on 1A, 1B and 1D; these QTLs were not related to DG, WG, or PC. Of these, QMTxV1B-31 exhibited the maximum PVE of 16.4%, which was strongly influenced by GI and SV. QMTxV1D-1 (Glu-D1-WPT-667287) was not identified when conditioned on SV, GI, DG or WG, with a PVE of 7.6-10.4%; this indicates that the QTLs associated with MTxV were completely influenced by SV, GI, DG and WG. QMTxV1A.1-15 (Glu-A1-WPT-672089) was identified with a PVE of 8.07%; this QTL was not related to PC, SV, WG or DG, but it was crucially affected by GI.

Genotype-by-environment (G×E) interactions were investigated in the present study. Of the QTLs identified in the G×E analysis, QMPV1B.1-4, QCTV1B.1-26 and QMTxV1B.1-26 showed A×E% effects, with PVEs of 0.1%, 0.4% and 0.1%, respectively. Although the PVEs of the G×E QTLs were relatively low, these QTLs warrant attention for their possible role in wheat quality improvement [19].

 Table 1 

Phenotypic values for mixographic parameters of the RIL population in different environments.

TraitEnvSN01-35GC9411MinMaxMeanCV%H2 (%)
MPTE11.83.51.54.32.522.6257.82
E21.55.51.56.33.232.16
E31.64.91.75.62.925.26
AE1.64.61.54.92.924.94
MPVE163.167.850.976.567.06.8765.23
E254.951.741.259.752.66.47
E356.457.448.064.955.96.11
AE58.259.047.970.558.85.98
CTVE120.728.320.461.951.79.9549.89
E235.248.330.948.743.07.12
E340.547.830.850.543.17.61
AE32.141.533.056.146.27.27
MTxVE148.457.739.263.853.68.1267.13
E236.350.232.950.144.56.82
E341.750.235.053.544.77.12
AE42.152.736.058.147.86.87
MPWE120.728.35.240.625.319.8962.19
E215.118.49.523.517.315.59
E321.118.112.230.919.418.48
AE19.021.69.928.920.913.70
CTWE13.911.13.251.86.780.3253.18
E22.911.92.513.96.036.01
E33.16.62.917.65.035.64
AE3.39.93.320.35.941.86
EWE14.214.23.321.67.249.2749.82
E23.113.73.014.06.837.54
E33.410.03.113.15.337.83
AE3.612.63.316.86.538.87
SLE1-8.3-0.615.77.358.2450.97
E2-1.20.29.94.546.13
E3-4.61.116.36.139.46
AE-4.71.116.36.141.76
SRE1-3.1-3.5-5.6-0.5-3.130.2053.85
E2-1.7-1.1-3.6-0.4-2.036.97
E3-2.1-2.1-4.1-0.4-2.526.82
AE-2.3-2.3-3.8-0.9-2.523.30
ITE1515.2585.9430.8661.7565.37.2562.10
E2405.3470.5356.4506.3458.75.69
E3448.7489.6403.5544.9472.35.25
AE456.4515.3398.7605.4500.95.75

In the environment column, E1, E2, E3 refer to 2009 Tai'an, 2010 Tai'an, and 2010 Suzhou; AE refer to the average data of three environments. MPT, midline peak time; MPV, midline peak value; CTV, curve tail value; MTxV, midline time-8 min value; MPW, midline peak width; EW, end width; CTW, curve tail width; SL, left slope of peak; IT, integral of curve tail; CV, coefficients of variation; H2, the broad-sense heritability.

 Table 2 

Unconditional and conditional additive QTLs for mixing parameters MPV, CTV, MTxV.

TraitQTLEFlanking markersPositionAaUnconditional PVE (A)%bPVE
(AE)%c
Conditional PVE (A)%d
PCSVGIWGDG
MPVQMPV1B.1-4E1/E2/AEBOBWHITE_C6347_307-WSNP_BE637864B_TA_1_19.2-0.8322.3/19.5/-0.129.6/30/28.915.9/12.6/14.66.9/15.9/21.1-/-/7.5-/25.5/21.4
QMPV1B.1-6E3BS00066444_51-IAAV498511.1-1.5826.1-26.29.025.822.725.6
QMPV 1B.1-24E1/E2/AEWPT-8971-WPT-556236.2/--2.17-----26.5/17.5/8.927.1/-/-
QMPV 7B.4-22E1RAC875_REP_C72877_159-WPT-546320.1-1.2313.0-11.2----
CTVQCTV1A.1-6AEBS00022677_51-TDURUM_CONTIG44888_8377.1-0.614.57------
QCTV1A.1-15AEGLU-A1-WPT-67208922-1.177.00------
QCTV 1B.1-26E2/AEWPT-8682-TDURUM_CONTIG98378_45238.2--0.420.621.420.9
QCTV 1B.1-31E2KUKRI_C8255_1324-TDURUM_CONTIG30113_21441.5-1.3617.3------
QCTV 3A.2-236E3EXCALIBUR_C31571_136-XGPW8072123.3-0.671.1------
QCTV 5A.1-5AEBS00001525_51-BOBWHITE_C31599_60414.4-0.6----5.2--
MTxVQMTxV 1A.1-15AE/E1GLU-A1-WPT-67208922-1.36/--8.07/--7.34/5.738.29/6.18-6.25/5.526.59/5.65
QMTxV1B.1-26E2WPT-8682-TDURUM_CONTIG98378_45239.2-1.98-0.120.5--21.421
QMTxV 1B.1-31E3KUKRI_C8255_1324-TDURUM_CONTIG30113_21441.5-1.2916.4-13.7--14.614.4
QMTxV 1D-1E2/E3/AEGLUD1-WPT-66728711-1.04/-1.112/-1.0210.4/11.5/
12.0
-7.6/-/-----
QMTxV 1D-4AEWPT-3743-RAC875_c48669_29215.4-0.98-----7.20.92

a Additive effect of the QTL. A positive number indicates that the Shannong01-35 allele was associated with larger trait values than the Gaocheng9411 allele; a negative number indicates that the Gaocheng9411 allele was associated with larger trait values than Shannong01-35; b Phenotypic variance explained by additive effects in unconditional analysis; c Phenotypic variance explained by additive-by-environment interaction effects; d Phenotypic variance explained by additive effects in conditional analysis. PC, Protein content; SV, sedimentation volume; GI, gluten index; WG, wet gluten content; DG, dry gluten content.

Additive QTLs and genotype-by-environment interactions (G×E) for mixing properties MPW, CTW and EW

Important QTL clusters associated with MPW were identified on chromosomes 1B, 4B, and 1D (Table 3, Figure 1). QMPW1B.1-4 was identified in unconditional analysis, with PVEs of 4.7-12.1% and net synergistic effects of 0.94-1.48%. The data indicated that cluster QMPW1B.1-4 was completely influenced by SV, WG, and DG and was not detected when conditioned on these traits. QMPW1B.1-24 was identified only in the conditional analysis, with PVEs of 6.1-15.9%, indicating that this QTL was identified when the influence of PC and WG was excluded. QMPW4B.5-112 was identified only when conditioned on WG and DG. QMPW1D-1, with the PVE of 11.2%, was identified only when conditioned on WG and DG.

Five QTLs were identified as associated with CTW. These QTLs, with PVEs of 3.8-11.2%, were identified only in the conditional analysis, indicating that they were influenced by the quality and quantity of protein. QCTW3B.3-31 was a stable QTL identified in multiple environments. QCTW1D-1 (PVE of 13.04%) was detected in two environments, with an increased PVE compared with the corresponding unconditional QTL (unconditional PVE of 7.5%). QCTW1D-1 was not identified when conditioned on SV, GI, DG or WG, indicating that QCTW1D-1 was completely influenced by the above traits and was partially dependent on PC.

QTLs associated with EW were located on chromosomes 1A, 1D and 3B. QEW1D-1 showed the greatest contribution to EW variation, with PVEs of 16.7%. The PVEs of QEW1D-1 when conditioned on protein content did not differ significantly, and its PVE was decreased when conditioned on DG (PVE of 11.9%) compared to the corresponding unconditional QTL (unconditional PVE of 15.6%). QEW1D-1 was not identified when conditioned on SV, GI or WG, demonstrating that this QTL was independent of PC, completely dependent on SV, GI and WG and partially influenced by DG. QEW1A.1-8 was identified in multiple environments and showed PVEs of 6.63-10.04% in both conditional and unconditional mapping. QEW3B.3-31 was the major QTL identified in multiple environments (E1 and AE), with PVEs of 2.4-11.2%; this QTL was independent of SV.

A G×E interaction of QEW3B.3-31 was identified with an AE of 1.2%.

 Table 3 

Unconditional and conditional additive QTLs for mixing parameters MPW, CTW, and EW.

TraitQTLEFlanking markersPositionAaUnconditional PVE (A)%bPVE
(AE)%c
Conditional PVE (A)%d
PCSVGIWGDG
MPWQMPW1B.1-4E1/E3/AEBOBWHITE_C6347_307-WSNP_BE637864B_TA_1_18.2-0.94/-1.48/-1.012.1/-/4.7--/19.9/8.3/-----
QMPW1B.1-5E3WSNP_BE637864B_TA_1_1-BS00066444_5111.1-1.6612.1-12.1--19.7-20.2-20.1-
QMPW 1B.1-24E2/AEWPT-8971-WPT-556232.2-1.01/---14.7/--/14.6/--/14.315.0/12.5-/15.9
QMPW1D.1-1E3GLU-D1-WPT-6672879.01.25--11.2----
QMPW 4B.5-59E2RAC875_C2545_1127-WSNP_EX_C12526_1995164013.91.013.80------
QMPW 4B.5-58E2RA_C41921_1056-RAC875_C2545_112713.91.18---20.3---
QMPW 4B.5-60E2WSNP_EX_C12526_19951640-TDURUM_CONTIG67477_13613.91.01----13.4--
QMPW4B.5-112E2TDURUM_CONTIG42107_2206-WSNP_EX_C16825_2538784165.80.78-----9.510.9
CTWQCTW1A.1-8E1/AEEXCALIBUR_C11258_1700-EX_C4206_5598.1-0.83--8.64--6.937.27
QCTW1 A.1-11E1/E2BS00026456_51-WPT_40299.7-0.77--9.36/6.68----
QCTW1D-1E2/AEGLU-D1-WPT-66728711.0/6-0.92/-0.86-7.5/2.5--/12.35----
QCTW 3B.3-31E1/AE/E2BS00064177_51-EXCALIBUR_C25566_42333.3-0.93/-0.57/-0.57--11.1/7.3/5.8--9.810.5/
3.8/11.2
QCTW6B.5-409AETDURUM_CONTIG30640_263-WSNP_CAP11_C1541_857160990.52/0.52--4.6/--/4.6--/3.8-/4.0
EWQEW1A.1-8E1/E2/AE/E3EXCALIBUR_C11258_1700-EX_C4206_5598.1-0.70/-7.56/--9.05/8.07/6.31/5.99-/10.04/
8.29/6.63
-/1.81-/7.49/
6.31/4.43
-/7.86/
6.63/5.13
QEW1D-1E1/E2/
AE
GLU-D1-WPT-6672877-1.60/-0.96/-1.16/-0.9315.6/11.4/
16.5/10.6
-16.5/13.0/16.7/----11.93/
11.2/-
QEW3B.3-31E1/AEBS00064177_51-EXCALIBUR_C25566_42333.3-/-0.83-/10.70-/1.211.2/
2.4
--9.6/-10.3/-

a b c d Same as Table 2.

 Figure 1 

Most significant additive QTLs for mixing properties mapped in the RIL population based on unconditional and conditional QTLs.

J Genomics Image

(View in new window)

 Table 4 

Unconditional and conditional additive QTLs for mixograph parameters MPT, SR and IT.

TraitQTLEFlanking markersPositionAaUnconditional
PVE (A)%b
PVE
(AE)%c
Conditional PVE (A)%d
PCSVGIWGDG
MPTQMPT1A.1-14AE/E1/E3WPT-8455-GLU-A116.1-0.35/-/-11.14/--12.69/14.36/13.8412.05/14.13/
15.76
--/-/9.04-/7.23/-
QMPT1A.1-15AEGLU-A1-WPT-67208937.2-0.17--5.2----
QMPT1A.1-24E2KUKRI_C13136_1363-BOBWHITE_C4384_26258.4-0.276.93--6.73---
QMPT1D-1E1/E2/E3/AEGLU-D1-WPT-66728711-0.35/-0.67
/-0.45/-0.47
33.4/27.1/33.7/
36.7
-37.9/33.0/
35.5/41.0
20.7/-/27.1/
26.7
-/-/-/13.214.9/-/22.2/
28.9
20.4/-/
-/31.25
SRQSR1A.1-14E2/AEWPT-8455-GLU-A116.1-0.29/-0.25/-7.41/8.10-/9.478.50/3.90--5.82/-
QSR1B.1-6E3BS00066444_51-IAAV498511.10.28.4-7.8/-10.8/3.2--
QSR1D-1E1/E2GLU-D1-WPT-66728710-0.31/-0.26-/11.59.5/15.5-/11.7---
QSR 2B.7-34E1/AEEXCALIBUR_C181_265-EXCALIBUR_REP_C107306_13032.70.31/0.21---11.3/12.1/----/3.2
ITQIT 1B.1-3AETPLB0025B13_372-BOBWHITE_C6347_3076.9-15.9230.6-40.9----
QIT 1B.1-4E1/E3/AEBOBWHITE_C6347_307-WSNP_BE637864B_TA_1_110.2-13.36/----15.2/-/16.725.2/-/23.0--/12.9/17.7
QIT 1B.1-6E3BS00066444_51-IAAV498511.1-11.621.9-13.8-11.518.812.9

a b c d Same as Table 2.

Additive QTLs and genotype-by-environment interactions (G×E) for mixing properties MPT, SR and IT

Four QTLs were detected as associated with MPT (Table 4, Figure 1). The PVE of QMPT1A.1-14 was changed when conditioned on PC and SV, indicating that this QTL was partially influenced by PC and SV and completely influenced by GI. QMPT1A.1-15 was not identified when conditioned on SV, GI, WG or DG. QMPT1D-1 explained a PVE of 13.2-41.0%. Compared to the corresponding unconditional QTL, the PVE of QMPT1D-1 was not significantly different when conditioned on PC, but it showed a high PVE of more than 20% when conditioned on SV, WG and DG. QMPT1D-1 was not identified when conditioned on GI. The above results indicate that QMPT1D-1 was completely influenced by GI, partially influenced by SV and gluten content, and not influenced by PC. QMPT1A.1-14 and QMPT1D-1 showed positive additive effects, indicating that favorable alleles were contributed by GC9411.

Four QTLs on chromosomes 1A, 1B, 1D and 2B were discovered to be related to SR. Among them, QSR1A.1-14 was partially influenced by PC, SV and DG, whereas it was completely dependent on GI and WG. QSR1B.1-6 and QSR2B.7-34 were identified in multiple environments with a high PVE of more than 10%, indicating that these QTLs were major loci with strong stability. QSR1D-1, which was identified by both conditional and unconditional QTL mapping, was influenced by PC. The favorable alleles of all the QTLs, except for QSR1B.1-6 and QSR2B.7-34, were contributed by GC9411.

One major QTL cluster on chromosome 1B was identified as associated with IT. The three QTLs were located close to each other, with a maximum PVE of 40.9% and a net additive effect of 15.9% min. QIT1B.1-4 was detected only when conditioned on GI, WG and DG, which should attract more attention in breeding programs.

The abovementioned QTLs for mixing properties were mainly distributed on chromosomes 1A, 3A, 5A, 1B, 2B, 7B and 1D. Among them, the favorable alleles were mostly contributed by GC9411.

Discussion

Although several QTLs for mixing properties have been identified in previous research, there is no consensus on the merits of mixing properties in the selection of optimal dough properties. Furthermore, most of these studies only focused on additive effects using unconditional QTL mapping methods without considering the epistatic or G×E interaction effects. In this study, additive and epistatic-interaction QTLs for mixing properties, as well as their environmental interaction effects, were analyzed. In addition, we combined conditional and unconditional QTL analyses to investigate the molecular mechanisms governing mixing properties.

Important QTLs associated with mixing properties

According to the present results, several QTL clusters were identified with high PVEs and strong stability, particularly those distributed on chromosomes 1A, 1B and 1D (Fig. 1). Locus 1A.1-14 within the genomic region WPT-8455-GLU-A1 was detected in multiple environments; this QTL was related to MPT and SR and explained a PVE of 7.41-12.69%. The effect was speculated to be closely related to Glu-A1. A comparison of map positions revealed that 1A.1-8 was 13.9 cM away from 1A.1-14 (WPT-8455-GLU-A1), which was predicted to be different from Glu-A1. Further investigations are needed to completely rule out any such possibility.

QTLs related to mixing parameters (MPV, CTV, MTxV, MPW, SR, and IT) were detected on chromosome 1B and have been reported in several previous studies [2-4]. Here, using SDS-PAGE analysis, we confirmed that GC9411 and SN01-35 have the same HMW glutenin subunits, 7+8, which are controlled by Glu-B1 on 1BL. Therefore, the QTLs on chromosome 1BS identified as controlling mixing properties in the present study are most likely not associated with the Glu-B1 locus in the population. Consensus results were obtained by Echeverry-Solarte [20] and Liu [16], who detected loci for increased IT located on 1BS. Further investigations are needed.

Furthermore, two major pleiotropic QTL clusters independent of Glu-B1 were detected on 1BS. One cluster, including 1B.1-24 and 1B.1-26, was associated with MPV and CTV. The other cluster, containing 1B.1-3 and 1B.1-4, was identified as associated with MPV and MPW and IT in multiple environments. In addition, 1B.1-6 and 1B.1-26 were detected in a single environment, of which 1B.1-6 and 1B.1-26 explained higher phenotypic variation. One major cluster that included three major QTLs (QIT1B.1-3, QIT1B.1-4 and QIT1B.1-6) for IT was located on chromosome 1B. Its maximum PVE was 40.9%, with a net additive effect of 40.9%, and its favorable alleles were contributed by GC9411. This cluster warrants greater attention in breeding programs.

The genomic region most consistently associated with MPT, MTxV, MPW, CTW, EW and SR were identified on chromosome 1D, especially 1D.1-1 (GLU-D1-WPT-667287), which was detected in all environments, with PVEs of 14.9%-20.7%. The high PVEs of these QTLs and their stability across environments further confirm the importance of HMW-GS encoded by Glu-D1. The closest region to these QTLs was previously identified as associated with mixing properties and other quality traits [3, 4, 21, 22, 23].

The high-molecular-weight glutenin subunit (HMW-GS) composition of Shannong01-35 was null, 7+8, 2+12; that of Gaocheng9411 was 1, 7+8, 5+10. The effect of Glu-D1 on quality was mainly due to the difference between HMW-GS 2+12 and 5+10. Lefebvre et al. (2000) showed that large gluten compositions dominated dough rheology [24], and glutenin subunit 5+10 was larger than 2+12 [25]. Many researchers have pointed out that subunit 5+10 had better quality characteristics of strong dough strength and food quality [26, 27]. Not surprisingly, GC9411 exhibited higher mixing parameter values than SN01-35 in all environments.

Conditional QTLs for mixing properties

Conditional analysis has the ability to discern the contribution of each component trait to a complex trait. Based on the QTL analysis of mixing properties, locus 1A.1-14, which controlled MPT and SR, was identified as partially influenced by PC, SV, and DG and completely influenced by GI. An important QTL, 1D-1, that was associated with MPT, EW, CTW and MTxV was completely influenced by GI, partially influenced by SV, WG and DG, and slightly influenced by PC. A QTL that controlled MPV and MPW, 1B.1-4, was partially influenced by PC, SV, and GI and was crucially affected by WG and DG. It is worth noting that 1B.1-24, an additional locus controlling MPV and MPW, was identified in multiple environments when the influence of the examined traits was excluded, but it was not identified in the unconditional analysis. QCTV1A.1-15, which was completely influenced by protein-related traits, was detected only in unconditional mapping. QCTV5A.1-5 was identified only when the influence of GI was excluded.

QMPW1B.1-4 was partially influenced by SV, PC, GI, and WG. QTL clusters 4B.5-59, 4B.5-58, and 4B.5-60 were partially influenced by all the examined traits. QCTW3B.3-31 and QCTW6B.5-409 were identified only in the conditional analysis. QEW1D-1 was completely influenced by SV and GI. Only two QTLs (QSR1B.1-6 and QSR2B.7-34) were found to be partially influenced by protein-related traits.

Conclusion

1A.1-15 and 1D.1-1 were major loci identified as associated with rheological properties influenced only by SV and GI. This finding demonstrates that these two QTL loci provide an important contribution to wheat quality at the genetic level. The combination of conditional and unconditional QTL mapping in the present study provided an opportunity for the detection of the mechanisms underlying the mixing properties of wheat.

Acknowledgements

This work was supported by the National Natural Science Foundation of China (31971936) and the Science and Technology of Shandong projects 2019YQ028, 2019LZGC001-3 and ZR2020MC093. SNP analysis and the construction of genetic maps were kindly conducted by Dr. Mingcheng Luo from the University of California, Davis, and by Dr. Jirui Wang of Sichuan Agricultural University. Prof. Wujun Ma from Murdoch University and Wolfgang Friedt from University Giessen provided valuable revision suggestions.

Competing Interests

The authors have declared that no competing interest exists.

References

1. Goesaert H, Brijs K, Veraverbeke WS. et al. Wheat flour constituents: how they impact bread quality, and how to impact their functionality. Trends Food Sci. Technol. 2005;16:12-30

2. Sudha ML, Srivastava AK, Vetrimani R. et al. Fat replacement in soft dough biscuits: its implications on dough rheology and biscuit quality. J Food Eng. 2007;80:922-930

3. Huang XQ, Cloutier S, Lycar L. et al. 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

4. Nelson JC, Andreescu C, Breseghello F. et al. Quantitative trait locus analysis of wheat quality traits. Euphytica. 2006;149:145-159

5. Martinant JP, Nicolas Y, Bouguennec A. et al. Relationships between mixograph parameters and indices of wheat grain quality. J Cereal Sci. 1998;27:179-189

6. Wu YP, Zhang YL, Xiao YG. et al. QTL Mapping for Important Quality Traits in Common Wheat, Sci Agr Sinica. 2008; 41:331-339 (In Chinese with English abstract).

7. Tsilo TJ, Simsek S, Ohm JB. et al. Quantitative trait loci influencing endosperm texture, dough-mixing strength, and bread-making properties of the hard red spring wheat breeding lines. Genome. 2011;54:460C

8. Zhu J. keightley PD. Analysis of conditional genetic effects and variance components in developmental Genetics. Genetics. 1995;141:1633-1639

9. Wen YX, Zhu J. Multivariable conditional analysis for complex trait and its components, Acta Genet Sin. 2005; 32: 289-296.

10. Zhao J, Becker HC, Zhang D. Conditional QTL mapping of oil content in rapeseed with respect to protein content and traits related to plant development and grain yield, Theor Appl Genet. 2006; 113: 33-38.

11. Liu GF, Zhu HT, Liu SW. et al. Unconditional and conditional QTL mapping for the developmental behavior of tiller number in rice (Oryza sativa L.). Genetica. 2010;138:885-893

12. Qu YY, Pu M, Zhang HL. et al. Mapping QTLs of root morphological traits at different growth stages in rice. Genetica. 2008;133:187-200

13. Cui F, Li J, Ding AM. et al. Conditional QTL mapping for plant height with respect to the length of the spike and internode in two mapping populations of wheat. Theor Appl Genet. 2011;122:1517-1536

14. Li QF, Zhang Y, Liu TT, et al. Genetic analysis of kernel weight and kernel size in wheat (Triticum aestivum, L.) using unconditional and conditional QTL mapping, Mol Breeding. 2015; 35: 1-15

15. Wang L, Cui FA, Wang JP. et al. Conditional QTL mapping of protein content in wheat with respect to grain yield and its components. J Genetics. 2012;91:303-312

16. Liu T, Liu K, Wang F, et al. Conditional and unconditional QTLs mapping of gluten strength in common wheat (Triticum asetivum L.), J Integr Agri. 2017; 16: 2145-2155

17. Chao SM, Zhang WJ, Akhunov E. et al. Analysis of gene-derived SNP marker polymorphism in US wheat (Triticum aestivum L.) cultivars. Mol Breeding. 2009;23:23-33

18. Wang S, Wong D, Forrest K. 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.

19. Voorrips RE. Mapchart: software for the graphical presentation of linkage maps and QTLs. J Heredity. 2002;93:77-78

20. Echeverry-Solarte M, Kumar A, Kianian S. et al. New QTL alleles for quality-related traits in spring wheat revealed by RIL population derived from supernumerary × non-supernumerary spikelet genotypes, Theor Appl Genet. 2015; 128: 893-912.

21. Campbell KG, Bergman CJ, Gualberto DG. et al. Quantitative trait loci associated with kernel traits in a soft x hard wheat cross, Crop Sci. 2001; 41: 1275-1285.

22. Mann G, Diffey S, Cullis B. et al. Genetic control of wheat quality: interactions between chromosomal regions determining protein content and composition, dough rheology, and sponge and dough baking properties, Theor Appl Genet. 2009; 118: 1519-1537.

23. Simons K, Anderson JA, Mergoum M. et al. Genetic mapping analysis of bread-making quality traits in spring wheat, Crop Sci. 2012; 52: 2182-2197.

24. Lefebvre J, PopineauY, Deshayes G. et al. Temperature-induced changes in the dynamic rheological behaviour and size distribution of polymeric proteins for glutens from wheat near isogenic lines differing in HMW-GS composition. Cereal Chemistry. 2000;77:193-201

25. Don C, Lichtendonk WJ, Plijter JJ. et al. The effect of mixing on glutenin particle properties: aggregation factors that affect gluten function in dough Journal of Cereal Science. 2005; 41:69-83.

26. Shewry PR, Halford NG, Tatham AS. et al. The high molecular weight subunits of wheat glutenin and their role in determining wheat processing properties. Advances in Food and Nutrition Research. 2003;45:219-302

27. Maucher T, Figueroa JDC, Reule W. et al. Allelic variations of functional markers for high molecular weight glutenin genes in Indian wheat (Triticum aestivum L.) cultivars and their correlation with bread loaf volume. Plant Physiology. 2015;20(1):97-102

Author contact

Corresponding address Corresponding author: Jiansheng Chen, E-mail: jshchedu.cn


Received 2021-9-18
Accepted 2021-11-17
Published 2022-1-1