Identification and Validation of Stable Quantitative Trait Loci for SDS-Sedimentation Volume in Common Wheat (Triticum aestivum L.)

Sodium dodecyl sulfate-sedimentation volume is an important index to evaluate the gluten strength of common wheat and is closely related to baking quality. In this study, a total of 15 quantitative trait locus (QTL) for sodium dodecyl sulfate (SDS)-sedimentation volume (SSV) were identified by using a high-density genetic map including 2,474 single-nucleotide polymorphism (SNP) markers, which was constructed with a doubled haploid (DH) population derived from the cross between Non-gda3753 (ND3753) and Liangxing99 (LX99). Importantly, four environmentally stable QTLs were detected on chromosomes 1A, 2D, and 5D, respectively. Among them, the one with the largest effect was identified on chromosome 1A (designated as QSsv.cau-1A.1) explaining up to 39.67% of the phenotypic variance. Subsequently, QSsv.cau-1A.1 was dissected into two QTLs named as QSsv.cau-1A.1.1 and QSsv.cau-1A.1.2 by saturating the genetic linkage map of the chromosome 1A. Interestedly, favorable alleles of these two loci were from different parents. Due to the favorable allele of QSsv.cau-1A.1.1 was from the high-value parents ND3753 and revealed higher genetic effect, which explained 25.07% of the phenotypic variation, mapping of this locus was conducted by using BC₃F₁ and BC₃F₂ populations. By comparing the CS reference sequence, the physical interval of QSsv.cau-1A.1.1 was delimited into 14.9 Mb, with 89 putative high-confidence annotated genes. SSVs of different recombinants between QSsv.cau-1A.1.1 and QSsv.cau-1A.1 detected from DH and BC₃F₂ populations showed that these two loci had an obvious additive effect, of which the combination of two favorable loci had the high SSV, whereas recombinants with unfavorable loci had the lowest. These results provide further insight into the genetic basis of SSV and QSsv.cau-1A.1.1 will be an ideal target for positional cloning and wheat breeding programs.

Introduction

Common wheat is one of the most widely cultivated food crops in the world and one of the important sources of carbohydrates and proteins for human beings (Osinowo, 2011). For a long time, breeders of common wheat (Triticum aestivum L.) have mainly committed to improving wheat yield and disease resistance (Curtis and Halford, 2014). However, high-quality wheat has been demanded by consumers and industries, and wheat quality improvement has attracted increasing attention among breeders (Guzman et al., 2016). Gluten strength has a considerable influence on the end-use quality of wheat (Rubenthaler et al., 1990; Addo et al., 1991; Slade and Levine, 1994; Kweon et al., 2011; Liu et al., 2017b), which can be measured by various tests such as sodium dodecyl sulfate (SDS)-sedimentation volume (SSV), extensograph, farinograph, alveograph, and gluten index (Huang et al., 2006; Elangovan et al., 2008; Li et al., 2009; Kerfal et al., 2010; Tsilo et al., 2011). The SSV test is well correlated with gluten strength and bread-making quality of wheat (Axford et al., 1979; He et al., 2004; Ozturk et al., 2008), and exhibits advantages such as simplicity, low cost, small sample size requirement, and high efficiency. Therefore, it has been widely used for evaluating the content and quality of gluten protein and for fast screening desired cultivars in wheat breeding programs (Clarke et al., 2000).

Some recent studies have found that the quantitative nature of SSV is closely correlated to multiple genes encoding glutenins and gliadins, such as Glu-1, Glu-A3, Glu-B3, and Gli-B1 (Payne and Lawrence, 1983; Payne et al., 1984; Shewry et al., 2003; Maucher et al., 2009; Reif et al., 2011; Deng et al., 2015; Guo et al., 2020). Glutenins and gliadins are not only the most important storage proteins of wheat but also the main components of gluten protein (Gianibelli et al., 2001; Kerfal et al., 2010). Glutenins are related to the extensibility of gluten, while gliadins are associated with the elasticity of gluten (MacRitchie, 1995; Veraverbeke and Delcour, 2002; Van Der Borght et al., 2005; Rasheed et al., 2014). The content and ratio of glutenins and gliadins are the main factors that determine the wheat processing quality (Yang et al., 2014). In addition, puroindolines are a component of wheat grain protein and are closely related to grain hardness. Puroindoline b (Pinb-D1) gene was found to be related to the variation of SSV in recent studies (Park et al., 2010, 2012; Ahn et al., 2014; Würschum et al., 2016).

Sodium dodecyl sulfate (SDS)-sedimentation volume (SSV) is a complex quantitative trait affected by both environmental and genetic factors. Quantitative trait locus (QTL) analysis is an effective approach for examining the genetic basis of quantitative traits (Doerge, 2002). Many studies have analyzed QTLs for SSV. To date, QTLs for SSV have been detected on almost all chromosomes, explained 2.2–41.4% of the phenotypic variation (Li et al., 2009; Kerfal et al., 2010; Reif et al., 2011; Deng et al., 2015; Würschum et al., 2016; Liu et al., 2017a; Mir Drikvand et al., 2018; Goel et al., 2019; Guo et al., 2020; Yang et al., 2020). However, most of the previously identified QTLs were detected only in one or two environments and could not be detected in multiple genetic backgrounds, which were not ideal targets for fine mapping and map-based cloning.

Here, to understand the genetic basis underlying SSV and provide molecular markers linked to QTL for wheat quality breeding, a doubled haploid (DH) population derived from a cross between Non-gda3753 (ND3753) and Liangxing99 (LX99) was employed to detect the QTLs associated with SSV variation. The genetic effect of two major SSV QTLs was validated.

Materials and Methods

Plant Materials

The DH population consisting of 123 individuals was developed through in vitro anther culture (De Buyser and Henry, 1980) of the F₁ hybrids from a cross between ND3753 and LX99. The DH population and two parents were used for genome-wide identification of QTLs related to SSV.

For mapping of the possible QTL, ND3753 that carried the positive allele of the QTL in the confidence interval was crossed with LX99. F₁ plants were backcrossed with LX99 for three generations with insertion-deletion (InDel) marker-assisted selection to generate a BC₃F₁ population containing 418 plants. Subsequently, 126 heterozygotes lines at the QTL-anchored region were self-pollinated to BC₃F₂ containing 1,081 plants. This population with LX99 background is presented as BC₃F₂-L in the present paper.

In addition, in order to evaluate the effects of two QTL, LX99 was crossed with recurrent parent ND3753 and 64 BC₃F₁ heterozygotes lines at the QTL-anchored region were self-pollinated to construct another BC₃F₂ population containing 387 plants with marker-assisted selection. This population with the background of ND3753 is presented as BC₃F₂-N.

Field Trials

The DH population and two parents were planted in seven environments during the wheat-growing seasons of 2016, 2017, and 2018 in Beijing (BJ) (40°08′N, 116°10′E), Linfen (LF) (36°04′N, 111°31′ E), Xi’an (XA) (34°16′N, 108°55′E), and Cangzhou (CZ) (38°18′N, 116°49′E), China. The seven environments, namely 2016BJ, 2017BJ, 2017LF, 2017XA, 2018BJ, 2018LF, and 2018CZ were presented in this study as E1, E2, E3, E4, E5, E6, and E7, respectively. The field trials were conducted following a complete random block design with three biological replicates. However, only one biological duplication was harvested in E7 due to an accident of field management. Then 60 seeds for each of the lines and two parents were planted in two rows of 1.5 long and the row space was 20 cm.

The BC₃F₁, BC₃F₂-L, and BC₃F₂-N populations were all planted in Beijing. The BC₃F₁ population was planted in the wheat-growing seasons of 2018, while the BC₃F₂-L and BC₃F₂-N populations were planted in 2019. All of these backcross populations were sown in rows of 1.5 m long and 30 cm row space with a sowing density of 20 seeds per row. The BC populations were all planted in one trial and designed as a single replicate. During the whole growing season, the local standard field management methodologies were adopted for plant cultivation.

Evaluation of Traits

For DHs, 80 plants were harvested in each line of the seven environments and wheat flour of each line planted in E1 was obtained with a CD1 Quadrumat Junior laboratory mill (Chopin Technology, Paris, France), while the whole wheat flour of each line in the other six environments was produced by an XF-98B experimental mill (Zhenxing Electromechanical Instrument Factory, Cangzhou, China). SSV was determined according to a modified protocol of Axford et al. (1979) and Preston et al. (1982) using 2 g of samples. The specific procedure of the SSV test was similar to that described by Li et al. (2009). In particular, SSVs of each DH line with only one biological duplication harvesting from E7 were measured. SSV of the BC₃F₁, BC₃F₂-L, and BC₃F₂-N populations was measured with whole wheat flour from a single plant.

In addition, 300 g of grains of 30 randomly selected DH lines based on the minimum sample required for the Pearson’s correlation coefficient calculated according to the formula provided by Mangard et al. (2007) and Chen et al. (2011) in E5 were ground into flour with a flour yield of approximately 60% in all samples. Their farinograph parameters (GB/ICC) were recorded by a Farinograph (DongFu JiuHeng, Beijing, China) to evaluate the correlation with SSV (Chicago, IL, United States) (ICC, 1996; Luo et al., 2018).

Genetic Map Construction

Deoxyribonucleic acid (DNA) was extracted from fresh leaves of individual DH lines and two parents using the hexadecyltrimethy ammonium bromide (CTAB) method (Allen et al., 2006). The 15 K Axiom^® Wheat Breeder single-nucleotide polymorphism (SNP) Genotyping Array (China Golden Marker Co., Beijing, China) containing 13,947 SNP markers was used to genotype the DH population and parents. SNP markers with a missing data rate > 20% were removed, and the remaining polymorphic markers were used to construct a wheat genetic map based on the inclusive composite interval mapping (ICIM) method using IciMapping v4.1 (Chinese Academy of Agricultural Sciences, China) and MapChart v2.32 (Plant Research International, P.O. Box 16, 6700 AA Wageningen) (Voorrips, 2002). The physical locations of unique SNP markers were obtained from the International Wheat Genome Sequencing Consortium (Appels et al., 2018).

Quantitative Trait Locus Mapping

The average value of SSV in each environment and the BLUP were employed for QTL analysis using inclusive composite interval mapping (ICIM¹) method in software IciMapping v4.1 (Meng et al., 2015). A QTL with LOD ≥ 2.5 was defined as a significant QTL. The confidence intervals (±2 LOD away from the peaks of likelihood ratios) of several QTLs were coincident, which were preliminarily considered as the same QTL. In this study, the QTL that can be detected in three or more environments is defined as an environmentally stable QTL.

Re-sequencing and InDel Markers Development

High-quality genomic DNA of ND3753 and LX99 was extracted to construct paired-end-sequencing libraries. According to the procedures described by Li et al. (2020), the parents were re-sequenced with an average sequencing depth of 6 × and paired-end reads of length 150 bp for two parents using the Illumina HiSeq X Ten platform (Illumina, California, United States), and the re-sequencing data were processed. The InDels were identified using the HaplotypeCaller module of the Genome Analysis Toolkit (GATK). The InDel markers were developed based on the sequence difference between the parents around the target region. Primer3 version 0.4.0² was used to design the sequences of InDel primers.

Deoxyribonucleic acid (DNA) amplification was programmed for an initial 5 min at 94°C, then followed by 35 cycles of 30 s at 94°C, 30 s at 56°C, and 30 s at 72°C, and finally 5 min at 72°C. A 10 μL PCR reaction system was used, containing 5 μL of 2 × Taq PCR StarMix (GenStar, Beijing, China) (for PAGE), 1.5 μL of DNA template (about 50–100 ng), 1.5 μL of each InDel primer, and double-distilled H₂O. The PCR products were analyzed on 8% non-denaturing polyacrylamide gels with silver staining.

Validation and Mapping and Annotation of Putative Genes

Insertion-deletion (InDel) markers tightly linked to QSsv.cau-1A.1.1 and QSsv.cau-1A.1.2 were used to genotype DH, BC₃F₁, BC₃F₂-L, and BC₃F₂-N populations. The putatively annotated high confidence (HC) genes located between the flanking markers of QSsv.cau-1A.1.1 were acquired based on the reference genome of T. aestivum cv. Chinese Spring³.

Statistical Analysis

To conduct phenotypic statistical analysis and calculate correlation coefficients between SSV and farinograph parameters, IBM SPSS Statistics 21.0 (SPSS, Chicago, United States) was used. R software v3.6.2⁴ was used to perform the Shapiro-Wilk test across seven environments and the best linear unbiased prediction (BLUP), as well as to estimate the broad-sense heritability (h2B) following the formula: $h_B^2={\mathrm{\sigma }}_g^2/({\mathrm{\sigma }}_g^2+{\mathrm{\sigma }}_{ge}^2/n+{\mathrm{\sigma }}^2/nr)$, as described by Liu et al. (2014).

Results

Phenotypic Analysis

Descriptive statistics for SSV of two parents and the DH population in the seven environments are shown in Table 1. SSV of ND3753 was significantly higher than that of LX99 in all environments. SSV showed bi-directional transgressive segregation, suggesting that both parents have increasing alleles for SSV. The broad-sense heritability values in all environments were greater than 0.8, indicating that SSV was mainly controlled by genetic factors. The result of the Shapiro-Wilk test displayed that SSV exhibited normal distribution under six environments and BLUP value, indicating SSV was determined by many genes (Figure 1). Pearson’s correlation coefficient analysis was carried out between farinograph parameters and SSVs of 30 DH lines planted in E5 which indicated a significant positive correlation between the two (Table 2).

TABLE 1

Table 1. Descriptive statistics of two parents and DH population for SDS-sedimentation volume (SSV) under seven environments.

FIGURE 1

Figure 1. Histograms of the ND3753/LX99 DH population for SDS-sedimentation volume (SSV) under BLUP data. The Y-axis represents the density (the ratio of frequency to group distance) of each trait and the X-axis represents the phenotypic data.

TABLE 2

Table 2. Pearson’s correlation analysis among SDS-sedimentation volume (SSV), dough stability time (DST), dough developing time (DDT), and water absorption (WA) in E5/2018–2019 (Beijing).

Linkage Map Construction

A total of 2,523 SNP markers showed polymorphisms between the two parents. Finally, 2,474 SNP markers participated in the map construction and were mapped to 21 linkage groups, covering the 21 chromosomes of common wheat (Supplementary Table 1). The total length of the map was 7,349.01 cm, and the average interval distance between two adjacent markers was 7.24 cm (Supplementary Table 1). The A genome contained the most SNP markers (963), followed by the B genome (902), while the D genome had the least (609) (Supplementary Table 1). The total length of chromosome 7D was the largest (619.02 cm), while that of chromosome 4B was the smallest (180.86 cm) (Supplementary Table 1). Chromosome 5A harbored the most SNP markers (236), while chromosome 6D contained the least (42) (Supplementary Table 1).

Quantitative Trait Locus Analysis

A total of 15 QTLs were detected on 11 chromosomes (1A, 1B, 1D, 2A, 2D, 4B, 4D, 5A, 5D, 6B, and 6D) in the seven environments (Table 3 and Supplementary Table 2). Four environmentally stable QTLs (QSsv.cau-1A.1, QSsv.cau-1A.2, QSsv.cau-2D and QSsv.cau-5D.1) were identified on chromosomes 1A, 1A, 2D, and 5D, respectively (Table 3). The favorable allele of QSsv.cau-2D came from LX99, while the superior alleles of the other three QTLs were contributed by ND3753. The major QTL QSsv.cau-1A.1 was repeatedly detected in five environments and the BLUP data, explaining 39.67% of the phenotypic variation in the BLUP analysis (Table 3). QSsv.cau-2D and QSsv.cau-5D.1 contributed 3.17 and 4.82% of the phenotypic variation in the BLUP analysis, respectively. QSsv.cau-1A.2 explained 8.17–18.62% of the phenotypic variation. The remaining 11 were putative QTLs (Supplementary Table 2).

TABLE 3

Table 3. The QTL regions harboring environmentally stable QTLs for SSV in the ND3753/LX99 DH population.

Verification of QSsv.cau-1A.1

Considering its stability and the genetic effect, QSsv.cau-1A.1 was chosen to saturate the positioning interval. Then 15 InDel markers near this region were developed according to the re-sequencing results of two parents (Table 3 and Supplementary Table 3) and a new genetic linkage map of 1A long arm was constructed. QTL remapping detected an additional QTL located next to the original interval. These two QTLs were named QSsv.cau-1A.1.1 and QSsv.cau-1A.1.2 (Figure 2).

FIGURE 2

Figure 2. Genetic and physical locations of QSsv.cau-1A.1.1 and QSsv.cau-1A.1.2. (A) QTL mapping for SSV in seven individual environments (E1–E7) and BLUP using a saturated genetic map of chromosome arm 1AL. (B) Saturated genetic map of chromosome arm 1AL with newly integrated InDel markers in the DH population. The black and red rectangles indicate QTLs with positive alleles from the parent ND3753 and the parent LX99, respectively. The black horizontal lines and bars above the genetic map represent the confidence interval of the two QTLs. (C) Corresponding physical positions according to the Chinese Spring IWGSC RefSeq v1.0 sequence.

QSsv.cau-1A.1.1 was detected in six environments and BLUP, explaining 17.21–26.47% of the phenotypic variation, and the favorable allele was from ND3753. The confidence interval was between the markers 1A248 and 1A332 corresponding CS physical position of 386,222,224–463,923,653 bp (Appels et al., 2018; Figure 2 and Supplementary Table 5). QSsv.cau-1A.1.2 was repeatedly detected in three environments as well as BLUP data, contributing 7.02–12.13% of the phenotypic variation and LX99 contributed the favorable allele. The physical position of QSsv.cau-1A.1.2 located on 462,634,655–492,004,197 bp by comparing flanking markers 1A36 and 1A215 to CS RefSeqv1.0 (Appels et al., 2018; Figure 2 and Supplementary Table 5).

Effects of QSsv-cau-1A.1.1 and QSsv-cau-1A.1.2 in Different Genetic Backgrounds

The flanking markers 1A248, 1A332 and 1A36, 1A215 delimiting confidence intervals of QSsv.cau-1A.1.1 and QSsv.cau-1A.1.2 (Figure 2) separately were used to detect genotypes in the DH and BC₃F₂-N populations. AA and aa represented genotypes with homozygous favorable and unfavorable alleles of QSsv.cau-1A.1.1 from ND3753, respectively, whereas BB and bb were symbols of that of QSsv.cau-1A.1.2 from LX99.

DH lines could be grouped into three genotypes which included two parental genotypes AAbb, aaBB, and one recombined genotype aabb, and each contained 51, 59, and 13 lines. The reason for the absence of genotype AABB remained unknown. The average SSV values in BLUP of genotype AAbb, aaBB, and aabb were 21.2, 19.9, and 17.8 mL, respectively, of which AAbb was significantly higher than aaBB and both were significantly higher than that of aabb (Figure 3A). This suggested that QSsv.cau-1A.1.1 had a stronger effect on SSV than QSsv.cau-1A.1.2, which was consistent with their contribution rates of phenotypic variation and additive effects in QTL analysis.

FIGURE 3

Figure 3. The phenotypic effect of QSsv.cau-1A.1.1 and QSsv.cau-1A.1.2 in the DH and BC₃F₂-N populations according to the BLUP value for SSV and the means of SSVs of different types, respectively. (A) DH population; (B) BC₃F₂-N population; *, **, and *** indicate significant differences at the.05, 0.01, and.001 levels, respectively (Student’s t-test). aabb: QSsv.cau-1A.1.1 carrying homozygous alleles from LX99 and QSsv.cau-1A.1.2 carrying homozygous alleles from ND3753; aaBB: QSsv.cau-1A.1.1 carrying homozygous alleles from LX99 and QSsv.cau-1A.1.2 carrying homozygous alleles from LX99; AAbb: QSsv.cau-1A.1.1 carrying homozygous alleles from ND3753 and QSsv.cau-1A.1.2 carrying homozygous alleles from ND3753; AaBB: QSsv.cau-1A.1.1 carrying heterozygous alleles and QSsv.cau-1A.1.2 carrying homozygous alleles from LX99. The numbers in parentheses indicate sample size.

In BC₃F₂-N population four allele combinations AaBB, AAbb, aaBB, and aabb were identified, which had SSV average values of 20.4, 17.5, 17.9, and 15.8 mL from 2, 23, 20, and 19 plants, respectively. Similarly, the SSV values of genotypes AaBB, AAbb, and aaBB were significantly higher than that of aabb (Figure 3B). However, although genotype AaBB had the distinct highest value, three genotypes with favorable alleles had no significant difference. This was possibly due to the too-small sample number of genotypes AaBB. Nevertheless, all results above could still prove that the favorable allele has positive effects.

This study did not find the combination type of AABB but found the type AaBB in the BC₃F₂-N population. In summary, the combination of two favorable loci had the high SSV, whereas recombinants with unfavorable loci had the lowest. In future research, the homozygous lines of the BC₃F₂-N population will be extracted to develop the corresponding NIL pairs to further verify their effects and examine their genetic effect on some wheat qualities, such as gluten content, extensograph, farinograph, alveograph, and gluten index.

Mapping of QSsv-cau-1A.1.1

In order to further verify and narrow down the confidence interval of QSsv.cau-1A.1.1, the BC₃F₁ population in the background of LX99 was genotyped using six InDel markers, and four recombinant types were obtained (Figure 4B). The SSV values of types 3 and 4 were similar and were significantly lower than that of types 1 and 2, indicating that QSsv.cau-1A.1.1 was delimited to the interval between markers 1A1 and 1A366.

FIGURE 4

Figure 4. Mapping of QSsv-cau-1A.1.1. (A) Genetic location of the region of interest on chromosome arm 1AL. Graphical illustration of recombinant genotypes from (B) the BC₃F₁ population and (C) the BC₃F₂ population in the LX99 background. (D) SSV values (mean ± SD). Black, gray, and white bars represent the ND3753 genotype, heterozygous genotype, and the LX99 genotype, respectively. The arrow indicates the 14.9-Mb mapping interval. Significant differences by comparing with the ND3753 genotype are indicated by * (p < 0.05), and ** (p < 0.05) (Student’s t-test). The numbers in and outside the parentheses indicate sample size and recombinant types, respectively.

The BC₃F₁ individuals with heterozygous genotypes between markers 1A1 and 1A39 were selected to generate the BC₃F₂-L population. The six markers between 1A1 and 1A39 were used to genotype the BC₃F₂-L population, and six recombinant types were obtained (Figure 4C). The SSV value of type 6 was similar to that of types 4 and 5 but was significantly lower than that of the other types (Figure 4D). These results once again proved the effectiveness of QSsv.cau-1A.1.1 and further confirmed that its location was between markers 1A1 and 1A366 (Figure 4C). There was no significant phenotypic difference between the heterozygous genotype and the ND3753 genotype, which implied that the ND3753 genotype was dominant. The corresponding physical interval of CS between these two markers was 14.9 Mb, with 89 putative high-confidence annotated genes (Supplementary Table 6). In addition, further fine mapping of QSsv-cau-1A.1.1 is under research.

Discussion

Correlation of Sodium Dodecyl Sulfate-Sedimentation Volume and Farinograph Parameters

Sodium dodecyl sulfate (SDS)-sedimentation volume (SSV) is a comprehensive indicator for indirectly testing wheat quality and one of the important tests to evaluate the gluten strength of flour and is closely related to the processing and baking quality of flour (Axford et al., 1978; Peña-Bautista, 2002; He et al., 2004). SSV is well correlated with other quality traits, such as grain protein content, gluten index, wet gluten content, bread volume, and farinograph parameters (Cubadda et al., 1992).

Our results showed that SSV was significantly positively correlated with stability time, formation time, and water absorption of dough, which is consistent with previous studies (Table 2; Cubadda et al., 1992). This confirms that SSV can be used as a suitable substitute for farinograph indicators that require a great number of samples to evaluate the rheological properties of wheat dough, thereby indirectly measuring the baking and processing quality of wheat flour.

Novel Quantitative Trait Locus for Sodium Dodecyl Sulfate-Sedimentation Volume

We compared the physical locations of QTLs for SSV reported in previous studies and those revealed in this study (Li et al., 2009; Reif et al., 2011; Deng et al., 2015; Würschum et al., 2016; Yang et al., 2020). The physical position of QSsv.cau-1A.2 was agreed with that of a previously reported QTL controlling SSV (Yang et al., 2020). SSV was found to be affected by allelic variations at Glu-A1 (508,726,618–508,725,448 bp, RefSeqv1.0) and Glu-A3 (4,203,001–4,202,275 bp, RefSeqv1.0) loci in several previous studies (Li et al., 2009; Reif et al., 2011; Deng et al., 2015; Würschum et al., 2016). Some QTLs associated with SSV was reported on chromosome 1A. For instance, Yang et al. (2020) identified a QTL (540,660,000–544,610,000 bp, RefSeqv1.0) for SSV that is located on chromosome 1A by genome-wide association study (GWAS). However, the physical position of QSsv.cau-1A.1.1 (371,573,909–386,426,688 bp, RefSeqv1.0) and QSsv.cau-1A.1.2 (419,490,584–492,004,197 bp, RefSeqv1.0) did not overlap with those of the above-mentioned QTLs/genes, suggesting these two QTLs may be novel. Li et al. (2009) identified a QTL (470,230,000–570,420,000 bp, RefSeqv1.0) for SSV on chromosome 2D using a recombinant inbred line population. Four QTLs (16,340,000, 59,102,000, 615,470,000, 646,600,000 bp, RefSeqv1.0) controlling SSV were reported to be located on chromosome 2D by multi-locus GWAS (Yang et al., 2020). However, the physical position of these QTLs and that of QSsv.cau-2D (140,759,212–467,689,413 bp, RefSeqv1.0) were not consistent, indicating that QSsv.cau-2D may also be a novel QTL. Li et al. (2009) identified a QTL for SSV on chromosome 5DS. SSV was found to be affected by the allelic variation at the Pinb-D1 locus on chromosome 5DS in some previous studies (Li et al., 2009; Reif et al., 2011; Deng et al., 2015; Würschum et al., 2016). However, the physical locations of these QTLs/genes and QSsv.cau-5D.1 does not match, implying that QSsv.cau-5D.1 on chromosome 5DL may be a novel QTL.

Sodium dodecyl sulfate (SDS)-sedimentation volume (SSV) is a quantitative trait affected by both environmental and genetic factors; thus, some QTLs can only be detected in specific environments (Supplementary Table 2). We found 11 such QTLs located on chromosomes 1B, 1D, 2A, 4B, 4D, 5A, 5D, 6B, and 6D, which is consistent with previous results (Li et al., 2009; Kerfal et al., 2010; Reif et al., 2011; Deng et al., 2015; Würschum et al., 2016; Liu et al., 2017a; Mir Drikvand et al., 2018; Goel et al., 2019; Guo et al., 2020; Yang et al., 2020). The Glu-D1 gene, which is located in the interval of QSsv.cau-1D, may be a candidate gene for QSsv.cau-1D. QSsv.cau-4B.1 and QSsv.cau-4D was located next to the dwarf genes Rht-B1 and Rht-D1, respectively. Previous studies have also revealed QTL-enrichment areas near Rht-B1 and Rht-D1, which are associated with kernel size, kernel hardness, kernel protein, pasting properties, and mixing properties (Shanhong et al., 2001; Li et al., 2006; Wang et al., 2012, 2017; Patil et al., 2013; Zhang et al., 2013; Jin et al., 2016; Liu et al., 2017b).

However, there is not enough evidence to support the correlation between Rht-B1/D1 and SSV, and the gene that controls SSV near Rht-B1 and Rht-D1 has not been cloned. We hypothesized: (1) there may be other genes affecting quality traits near Rht-B1 and Rht-D1; (2) allelic variations between Rht-B1 and Rht-D1 may also regulate certain quality traits, such as SSV. However, these hypotheses need to be further tested. The co-localization of dwarf genes and QTLs related to quality traits may also remind breeders to consider the selection of plant height and grain quality in the wheat breeding program.

Genetic Effects and Putative Annotated Genes of the Major Quantitative Trait Locus

In recent years, a large number of QTLs for SSV have been identified and characterized through GWAS and linkage analysis, and some of these QTLs are related to allelic variants of Glu-1, Glu-A3, Glu-B3, Gli-B1, and Pina-D1 (Ahn et al., 2014; Deng et al., 2015; Guo et al., 2020). However, most of the other QTLs have not been further verified or fine mapped. SSV is a typical quantitative trait with a complex genetic mechanism. The lack of information on the authenticity and genetic effects of these QTLs for SSV not only hinders the exploration of their genetic and molecular mechanisms but also fails to provide breeders with sufficient new high-quality genetic resources for wheat quality improvement. In this study, we verified the effects of QSsv.cau-1A.1.1 and QSsv.cau-1A.1.2 on SSV and the interaction between the two QTLs in the DH and BC3F2-N populations. Therefore, the InDel markers are closely linked to QSsv.cau-1A.1.1 and QSsv.cau-1A.1.2 developed in this study can be used by breeders to aggregate high-quality genes for wheat quality improvement.

In particular, QSsv.cau-1A.1.1 was delimited to an approximate 14.9 Mb between markers 1A1 and 1A366 (Figure 4). We conduct an orthologous analysis for the candidate region to predict HC genes in QSsv.cau-1A.1.1, but no ones are associated with SSV in Oryza sativa and Arabidopsis thaliana (Supplementary Table 6). This may be because Oryza sativa and Arabidopsis thaliana do not have gluten, and SSV is related to gluten strength in the common wheat.

Data Availability Statement

The raw data of Nongda3753 and Liangxing99 presented in the study are deposited in the NCBI Sequence Read Archive repository, accession number PRJNA722149.

Author Contributions

MY conceived the project. ST performed the research, constructed the linkage map, developed an InDel marker of the QTL region of interest, and developed the BC₃F₁ and BC₃F₂ populations. JL, MZ, and SW participated in the field trials. ST, CB, HZ, CW, ZC, JY, and XS performed the phenotypic analysis. YZ, CX, BL, RL, and QS assisted in revising the manuscript. ST and YZ analyzed the experimental results. ST and MY wrote the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was financially supported by the National Key Research and Development Program of China (Grant No. 2016YFD0101602).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors would like to thank Mingyi Zhang (Shanxi Academy of Agricultural Sciences) for help with field trials.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.747775/full#supplementary-material

Footnotes

1. ^ http://www.isbreeding.net/software/?type=detail&id%20=18
2. ^ http://bioinfo.ut.ee/primer3-0.4.0/
3. ^ https://urgi.versailles.inra.fr/download/iwgsc/IWGSC_RefSeq_Annotations/v1.0/
4. ^ https://www.r-project.org/

References

Addo, K., Pomeranz, Y., Huang, M. L., Rubenthaler, G. L., and Jeffers, H. C. (1991). Steamed Bread. II. Role of protein content and strength. Cereal Chem. 68, 39–42.

Google Scholar

Ahn, J. H., Kang, C. S., Jeung, J. U., Baik, B. K., and Park, C. S. (2014). Effect of allelic variations at the Glu-D1, Glu-A3, Glu-B3 and Pinb-D1 loci on flour characteristics and bread loaf volume. Int. Food Res. J. 21, 1141–1149.

Google Scholar

Allen, G. C., Flores-Vergara, M. A., Krasynanski, S., Kumar, S., and Thompson, W. F. (2006). A modified protocol for rapid DNA isolation from plant tissues using cetyltrimethylammonium bromide. Nat. Protoc. 1, 2320–2325. doi: 10.1038/nprot.2006.384

PubMed Abstract | CrossRef Full Text | Google Scholar

Appels, R., Eversole, K., Feuillet, C., Keller, B., Rogers, J., Stein, N., et al. (2018). Shifting the limits in wheat research and breeding using a fully annotated reference genome. Science 361:eaar7191. doi: 10.1126/science.aar7191

PubMed Abstract | CrossRef Full Text | Google Scholar

Axford, D., Mcdermott, E. E., and Redman, D. G. (1979). Note on sodium dodecyl sulfate test of breadmaking quality; Comparison with Pelshenke and Zeleny test. Cereal Chem. 56, 582–584.

Google Scholar

Axford, D. W. E., McDermott, E. E., and Redman, D. G. (1978). Small scale tests of bread making quality. Milling Feed Fert. 13, 18–20.

Google Scholar

Chen, K., Zhang, P., Deng, G., and Zou, C. (2011). Electromagnetic correlation attacks and the minimum number of attack traces analyzing. J. Huazhong Univ. Sci. Technol. 39:1009.

Google Scholar

Clarke, J. M., Clarke, F. R., Ames, N. P., McCaig, T. N., and Knox, R. E. (2000). Evaluation of predictors of quality for use in early generation selection. Options Méditer. Sér. A Sémin. Méditerranéens 40, 439–446.

Google Scholar

Cubadda, R., Carcea, M., and Pasqui, L. A. (1992). Suitability of the Gluten Index method for assessing gluten strength in durum-wheat and semolina. Cereal Foods World 37, 866–869.

Google Scholar

Curtis, T., and Halford, N. G. (2014). Food security: the challenge of increasing wheat yield and the importance of not compromising food safety. Ann. Appl. Biol. 164, 354–372. doi: 10.1111/aab.12108

PubMed Abstract | CrossRef Full Text | Google Scholar

De Buyser, J., and Henry, Y. (1980). Induction of haploid and diploid plants though in vitro anther culture of haploid wheat (n=3x=21). Theor. Appl. Genet. 57, 57–58. doi: 10.1007/BF00745029

CrossRef Full Text | Google Scholar

Deng, Z., Tian, J., Chen, F., Li, W., Zheng, F., Chen, J., et al. (2015). Genetic dissection on wheat flour quality traits in two related populations. Euphytica 203, 221–235. doi: 10.1007/s10681-014-1318-7

CrossRef Full Text | Google Scholar

Doerge, R. W. (2002). Mapping and analysis of quantitative trait loci in experimental populations. Nat. Rev. Genet. 3, 43–52. doi: 10.1038/nrg703

PubMed Abstract | CrossRef Full Text | Google Scholar

Elangovan, M., Rai, R., Dholakia, B. B., Lagu, M. D., Tiwari, R., Gupta, R. K., et al. (2008). Molecular genetic mapping of quantitative trait loci associated with loaf volume in hexaploid wheat (Triticum aestivum). J. Cereal Sci. 47, 587–598. doi: 10.1016/j.jcs.2007.07.003

CrossRef Full Text | Google Scholar

Gianibelli, M. C., Larroque, O. R., MacRitchie, F., and Wrigley, C. W. (2001). Biochemical, genetic, and molecular characterization of wheat glutenin and its component subunits. Cereal Chem. 78, 635–646. doi: 10.1094/CCHEM.2001.78.6.635

CrossRef Full Text | Google Scholar

Goel, S., Singh, K., Singh, B., Grewal, S., Dwivedi, N., Alqarawi, A. A., et al. (2019). Analysis of genetic control and QTL mapping of essential wheat grain quality traits in a recombinant inbred population. PLoS One 14:e0200669. doi: 10.1371/journal.pone.0200669

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, Y., Zhang, G., Guo, B., Qu, C., Zhang, M., Kong, F., et al. (2020). QTL mapping for quality traits using a highdensity genetic map of wheat. PLoS One 15:e0230601. doi: 10.1371/journal.pone.0230601

PubMed Abstract | CrossRef Full Text | Google Scholar

Guzman, C., Peña, R. J., Singh, R., Autrique, E., Dreisigacker, S., Crossa, J., et al. (2016). Wheat quality improvement at CIMMYT and the use of genomic selection on it. Appl. Transl. Genom. 11, 3–8. doi: 10.1016/j.atg.2016.10.004

PubMed Abstract | CrossRef Full Text | Google Scholar

He, Z. H., Yang, J., Zhang, Y., Quail, K. J., and Peña, R. J. (2004). Pan bread and dry white Chinese noodle quality in Chinese winter wheats. Euphytica 139, 257–267. doi: 10.1007/s10681-004-3283-z

CrossRef Full Text | Google Scholar

Huang, X. Q., Cloutier, S., Lycar, L., Radovanovic, N., Humphreys, D. G., Noll, J. S., et al. (2006). 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. 113, 753–766. doi: 10.1007/s00122-006-0346-7

PubMed Abstract | CrossRef Full Text | Google Scholar

ICC (1996). Standard Methods of The International Association for Cereal Science and Technology. Vienna: ICC.

Google Scholar

Jin, H., Wen, W., Liu, J., Zhai, S., Zhang, Y., Yan, J., et al. (2016). Genome-wide QTL mapping for wheat processing quality parameters in a Gaocheng 8901/Zhoumai 16 recombinant inbred line population. Front. Plant Sci. 7:1032. doi: 10.3389/fpls.2016.01032

PubMed Abstract | CrossRef Full Text | Google Scholar

Kerfal, S., Giraldo, P., Rodriguez-Quijano, M., Vázquez, J. F., Adams, K., Lukow, O. M., et al. (2010). Mapping quantitative trait loci (QTLs) associated with dough quality in a soft$\times$ hard bread wheat progeny. J. Cereal Sci. 52, 46–52. doi: 10.1016/j.jcs.2010.03.001

CrossRef Full Text | Google Scholar

Kweon, M., Slade, L., and Levine, H. (2011). Solvent retention capacity (SRC) testing of wheat flour: principles and value in predicting flour functionality in different wheat-based food processes and in wheat breeding-A review. Cereal Chem. 88, 537–552. doi: 10.1094/CCHEM-07-11-0092

CrossRef Full Text | Google Scholar

Li, J., Wen, S., Fan, C., Zhang, M., Tian, S., Kang, W., et al. (2020). Characterization of a major quantitative trait locus on the short arm of chromosome 4B for spike number per unit area in common wheat (Triticum aestivum L.). Theor. Appl. Genet. 133, 2259–2269. doi: 10.1007/s00122-020-03595-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X. P., Lan, S. Q., Liu, Y. P., Gale, M. D., and Worland, T. J. (2006). Effects of different Rht-B1b, Rht-D1b and Rht-B1c dwarfing genes on agronomic characteristics in wheat. Cereal Res. Commun. 34, 919–924. doi: 10.1556/CRC.34.2006.2-3.220

CrossRef Full Text | Google Scholar

Li, Y., Song, Y., Zhou, R., Branlard, G., and Jia, J. (2009). Detection of QTLs for bread-making quality in wheat using a recombinant inbred line population. Plant Breed. 128, 235–243. doi: 10.1111/j.1439-0523.2008.01578.x

CrossRef Full Text | Google Scholar

Liu, G., Jia, L., Lu, L., Qin, D., Zhang, J., Guan, P., et al. (2014). Mapping QTLs of yield-related traits using RIL population derived from common wheat and Tibetan semi-wild wheat. Theor. Appl. Genet. 127, 2415–2432. doi: 10.1007/s00122-014-2387-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, T., An, Y., Liu, K., Wang, F., Xie, C., Zhang, Y., et al. (2017a). A genetic analysis of the quality of northern-style Chinese steamed bread. Mol. Breed. 37:41. doi: 10.1007/s11032-016-0593-5

CrossRef Full Text | Google Scholar

Liu, T., Liu, K., Wang, F., Zhang, Y., Li, Q., Zhang, K. R., et al. (2017b). Conditional and unconditional QTLs mapping of gluten strength in common wheat (Triticum aestivum L.). J. Integr. Agric. 16, 2145–2155. doi: 10.1016/S2095-3119(16)61564-2

CrossRef Full Text | Google Scholar

Luo, D., Kou, X., Zhang, T., Nie, Y., Xu, B., Li, P., et al. (2018). Effect of inulin on rheological properties of soft and strong wheat dough. Int. J. Food Sci. Technol. 53, 1648–1656. doi: 10.1111/ijfs.13748

CrossRef Full Text | Google Scholar

MacRitchie, F. (1995). “Role of polymeric proteins in flour functionality,” in Proceedings of the International Symposium, Wheat Kernel Proteins: Molecular and Functional Aspects, Viterbo, 145–150.

Google Scholar

Mangard, S., Oswald, E., and Popp, T. (2007). Power Analysis Attacks: Revealing the Secrets of Smart Cards. Berlin: Springer.

Google Scholar

Maucher, T., Figueroa, J. D. C., Reule, W., and Pen,  a, R. J. (2009). Influence of low molecular weight glutenins on viscoelastic properties of intact wheat kernels and their relation to functional properties of Wheat Dough. Cereal Chem. 86, 372–375. doi: 10.1094/CCHEM-86-4-0372

CrossRef Full Text | Google Scholar

Meng, L., Li, H., Zhang, L., and Wang, J. (2015). QTL IciMapping: integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental populations. Crop J. 3, 269–283. doi: 10.1016/j.cj.2015.01.001

CrossRef Full Text | Google Scholar

Mir Drikvand, R., Najafian, G., Bihamta, M. R., and Ebrahimi, A. (2018). Mapping some seed quality traits in bread wheat (Triticum aestivum L.) by association mapping using ssr markers. J. Appl. Biotechnol. Rep. 5, 92–99. doi: 10.29252/JABR.05.03.02

CrossRef Full Text | Google Scholar

Osinowo, F. (2011). Food and Agricultural Organization of the United Nations, State of the World’s Forests.

Google Scholar

Ozturk, S., Kahraman, K., Tiftik, B., and Koksel, H. (2008). Predicting the cookie quality of flours by using Mixolab^®. Eur. Food Res. Technol. 227, 1549–1554. doi: 10.1007/s00217-008-0879-x

CrossRef Full Text | Google Scholar

Park, C. S., Kang, C. S., Cheong, Y. K., Jung, W., and Woo, S. H. (2010). Influence of puroindoline genotypes on grain characteristics, physico-chemical properties of flour and end-use quality of Korean wheats. Breed. Sci. 60, 233–242. doi: 10.1270/jsbbs.60.233

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, C.-S., Kim, H.-S., Kim, D.-H., Hyun, J.-N., and Kang, C.-S. (2012). Environmental impacts of Korean and CIMMYT wheat lines on protein characteristics and bread making quality. Korean J. Crop Sci. 57, 60–70. doi: 10.7740/kjcs.2012.57.1.060

CrossRef Full Text | Google Scholar

Patil, R. M., Tamhankar, S. A., Oak, M. D., Raut, A. L., Honrao, B. K., Rao, V. S., et al. (2013). Mapping of QTL for agronomic traits and kernel characters in durum wheat (Triticum durum Desf.). Euphytica 190, 117–129. doi: 10.1007/s10681-012-0785-y

CrossRef Full Text | Google Scholar

Payne, P. I., Holt, L. M., Jackson, E. A., Law, C. N., and Damania, A. B. (1984). Wheat storage proteins: their genetics, and their potential for manipulation by plant breeding. Philos. Trans. R. Soc. Lond. B Biol. Sci. 304, 359–371. doi: 10.1098/rstb.1984.0031

CrossRef Full Text | Google Scholar

Payne, P. I., and Lawrence, G. J. (1983). Catalogue of alleles for the complex gene loci, Glu-A1, Glu-B1, and Glu-D1 which code for high-molecular-weight subunits of glutenin in hexaploid wheat. Cereal Res. Commun. 11, 29–35. doi: 10.1021/jf990151p

PubMed Abstract | CrossRef Full Text | Google Scholar

Peña-Bautista, R. J. (2002). “Wheat for bread and other foods,” in Bread Wheat: Improvement and Production, eds B. C. Curtis, S. Rajaram, and H. G. Macpherson (Rome: FAO Plant Production and Protection Series).

Google Scholar

Preston, K. R., March, P. R., and Tipples, K. H. (1982). An assessment of the sds-sedimentation test for the prediction of canadian bread wheat quality. Can. J. Plant Sci. 62, 545–553. doi: 10.4141/cjps82-083

CrossRef Full Text | Google Scholar

Rasheed, A., Xia, X., Yan, Y., Appels, R., Mahmood, T., and He, Z. (2014). Wheat seed storage proteins: advances in molecular genetics, diversity and breeding applications. J. Cereal Sci. 60, 11–24. doi: 10.1016/j.jcs.2014.01.020

CrossRef Full Text | Google Scholar

Reif, J. C., Gowda, M., Maurer, H. P., Longin, C. F. H., Korzun, V., Ebmeyer, E., et al. (2011). Association mapping for quality traits in soft winter wheat. Theor. Appl. Genet. 122, 961–970. doi: 10.1007/s00122-010-1502-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Rubenthaler, G. L., Huang, M. L., and Pomeranz, Y. (1990). Steamed bread. I. Chinese steamed bread formulation and interactions. Cereal Chem. 67, 471–475.

Google Scholar

Shanhong, W., Fanhua, M., Li, Y., and Binghua, L. (2001). The effect of dwarf genes on agronomy characters in wheat. J. Triticeae Crop. 21, 5–9.

Google Scholar

Shewry, P. R., Halford, N. G., and Lafiandra, D. (2003). Genetics of wheat gluten proteins. Adv. Genet. 49, 111–184. doi: 10.1016/S0065-2660(03)01003-4

CrossRef Full Text | Google Scholar

Slade, L., and Levine, H. (1994). Structure-function relationships of cookie and cracker ingredients. Sci. Cookie Cracker Prod. 9, 23–141.

Google Scholar

Tsilo, T. J., Simsek, S., Ohm, J. B., Hareland, G. A., Chao, S., and Anderson, J. A. (2011). Quantitative trait loci influencing endosperm texture, dough-mixing strength, and bread-making properties of the hard red spring wheat breeding lines. Genome 54:460470. doi: 10.1139/g11-012

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Der Borght, A., Goesaert, H., Veraverbeke, W. S., and Delcour, J. A. (2005). Fractionation of wheat and wheat flour into starch and gluten: overview of the main processes and the factors involved. J. Cereal Sci. 41, 221–237. doi: 10.1016/j.jcs.2004.09.008

CrossRef Full Text | Google Scholar

Veraverbeke, W. S., and Delcour, J. A. (2002). Wheat protein composition and properties of wheat glutenin in relation to breadmaking functionality. Crit. Rev. Food Sci. Nutr. 42, 179–208. doi: 10.1080/10408690290825510

PubMed Abstract | CrossRef Full Text | Google Scholar

Voorrips, R. E. (2002). MapChart: software for the graphical presentation of linkage maps and QTLs. J. Hered. 93, 77–78. doi: 10.1093/jhered/93.1.77

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, F., Liu, T., Li, Q., An, Y., Xie, C., Sun, X., et al. (2017). QTL mapping of the pasting properties of wheat flour treated by papain digestion. Starch Stärke 69:1600077. doi: 10.1002/star.201600077

CrossRef Full Text | Google Scholar

Wang, L., Cui, F., Wang, J., Jun, L., Ding, A., Zhao, C., et al. (2012). Conditional QTL mapping of protein content in wheat with respect to grain yield and its components. J. Genet. 91, 303–312. doi: 10.1007/s12041-012-0190-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Würschum, T., Leiser, W. L., Kazman, E., and Longin, C. F. H. (2016). Genetic control of protein content and sedimentation volume in European winter wheat cultivars. Theor. Appl. Genet. 129, 1685–1696. doi: 10.1007/s00122-016-2732-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y., Chai, Y., Zhang, X., Lu, S., Zhao, Z., Wei, D., et al. (2020). Multi-Locus GWAS of quality traits in bread wheat: mining more candidate genes and possible regulatory network. Front. Plant Sci. 11:1091. doi: 10.3389/fpls.2020.01091

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y., Li, S., Zhang, K., Dong, Z., Li, Y., An, X., et al. (2014). Efficient isolation of ion beam-induced mutants for homoeologous loci in common wheat and comparison of the contributions of Glu-1 loci to gluten functionality. Theor. Appl. Genet. 127, 359–372. doi: 10.1007/s00122-013-2224-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Dell, B., Biddulph, B., Drake-Brockman, F., Walker, E., Khan, N., et al. (2013). 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. 32, 771–783. doi: 10.1007/s11032-013-9905-1

CrossRef Full Text | Google Scholar