Next Article in Journal
Expanding the Natural History of SNORD118-Related Ribosomopathy: Hints from an Early-Diagnosed Patient with Leukoencephalopathy with Calcifications and Cysts and Overview of the Literature
Previous Article in Journal
Association Mapping of Candidate Genes Associated with Iron and Zinc Content in Rice (Oryza sativa L.) Grains
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Association Mapping of Processing Quality Traits in Common Wheat (Triticum aestivum L.)

1
Institute of Forage and Grassland Sciences, Heilongjiang Academy of Agricultural Sciences, Harbin 150086, China
2
National Wheat Improvement Center, Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100000, China
3
Horticultural Branch of Heilongjiang Academy of Agricultural Sciences, Harbin 150086, China
*
Authors to whom correspondence should be addressed.
Genes 2023, 14(9), 1816; https://doi.org/10.3390/genes14091816
Submission received: 4 August 2023 / Revised: 13 September 2023 / Accepted: 14 September 2023 / Published: 18 September 2023
(This article belongs to the Special Issue Genetics and Breeding of Grains)

Abstract

:
Processing quality is an important economic wheat trait. The marker-assisted selection (MAS) method plays a vital role in accelerating genetic improvement of processing quality. In the present study, processing quality in a panel of 165 cultivars grown in four environments was evaluated by mixograph. An association mapping analysis using 90 K and 660 K single nucleotide polymorphism (SNP) arrays identified 24 loci in chromosomes 1A, 1B (4), 1D, 2A, 2B (2), 3A, 3B, 3D (2), 4A (3), 4B, 5D (2), 6A, 7B (2) and 7D (2), explaining 10.2–42.5% of the phenotypic variances. Totally, 15 loci were stably detected in two or more environments. Nine loci coincided with known genes or QTL, whereas the other fifteen were novel loci. Seven candidate genes encoded 3-ketoacyl-CoA synthase, lipoxygenase, pyridoxal phosphate-dependent decarboxylase, sucrose synthase 3 and a plant lipid transfer protein/Par allergen. SNPs significantly associated with processing quality and accessions with more favorable alleles can be used for marker-assisted selection.

1. Introduction

Dough rheological properties with significant effects on end-use products can be evaluated by mixograph [1,2]. Midline peak time (MPT), midline 8 min band width (MTxW), mid-line peak width (MPW) and midline peak value (MPV) are related to processing quality [3,4]. Dough properties are highly variable among wheat cultivars [5,6] and processing quality is a major breeding objective in wheat breeding programs [7,8].
Dough rheological properties are quantitatively inherited and largely controlled by multiple minor genes [9,10,11,12]. It is difficult to evaluate dough rheological properties in traditional breeding because they cannot be measured in the early segregating generations of crosses due to limited quantities of seed. Moreover, measurement of mixograph-related traits requires a professional grain chemistry laboratory. Marker-assisted selection (MAS) could be a useful approach to improve wheat processing quality [7]. Dough strength is a typical quantitative trait controlled by multiple minor genes. A number of studies have been conducted which focus on identifying QTLs for mixographs and a series of QTLs for mixograph-related traits [2,13,14,15,16,17] were reported mainly by bi-parental linkage mapping. Barakat et al. [2] reported 108 loci for farinograph- and mixograph-related traits on all 21 wheat chromosomes in two double haploid (DH) populations. However, the short-coming of the method is that only two alleles at any single locus can be evaluated in each cross. Echeverry-Solarte et al. [17] reported 31 loci for mixograph-related traits by an RIL population developed from a cross of an elite wheat line (WCB414) and an exotic genotype with supernumerary spikelets, and each explained 3.2–41.2% phenotypic variations.
The results of those studies were mainly based on relatively low numbers of simple sequence repeat (SSR) or diversity array technology (DArT) markers [13,14] and were difficult to apply in gene cloning and MAS. The development of the wheat 90 K [18], 660 K [19] and 55 K SNP genotyping assays [19,20] has made it possible to genotype large populations with high-density SNPs. As a result, genome-wide association studies (GWAS) have been extensively conducted to explore the extant allelic diversity concerning numerous agronomic traits [21,22]. GWAS utilizes large amounts of markers distributed across the entire chromosomes of a species genome to find marker–trait relationships using LD as its basis, thereby uncovering significant positions. LD represents a naturally occurring phenomenon within a group during selection and evolution processes whereby nonrandom correlations occur among genes located at various sites within the same individual. If the probability of one specific allele existing at two distinct locations exceeds random expectations, then these two places exhibit LD. The process from unlinked imbalance to balanced linkage occurs throughout LD attenuation. Species differences exist regarding their distance of LD decay, with interspecific species such as maize and rapeseed exhibiting much farther distances than intraspecific ones such as wheat and rice. The level of LD also depends on factors including recombination, mutation, population structure, sample size, selective pressure, genetic drift, founder effects, admixture rates, etc. Among them, recombination and mutations play crucial roles in influencing LD levels [21,22,23].
To avoid false positives resulting due to any type of familiarity found amidst members comprising studied collectives alongside similar concerns about establishing appropriate mathematical formulas proves useful. General Linear Model (GLM), which was initially introduced, only considers impacts stemming from social arrangements, henceforth causing some degree of inaccurate outcomes because GLM does not take enough measures aimed towards controlling consequences originating via interactions caused jointly by kin and societal constructs. In contrast, MLM has been developed more recently since this methodology incorporates Q (population structure)-based assessments too—K (kinship); therefore, it better manages situations arising out of overlapping ramifications exerted either way by those concepts contributing toward finalized conclusions drawn downstream. At current times, researchers have implemented applicable usage of all sorts of techniques encompassing crop quantity features regularly utilizing mixed linier regression versions (MLM) [21,22,23].
Association analysis is a quantitative genetic analysis method based on LD among allelic variants at the same locus, using natural populations as research materials. By investigating the association between group genotype data and phenotypic data, target trait genes can be discovered. The advantage of this approach includes, firstly, that it does not require constructing biparental populations like linkage analysis but uses existing natural populations, high-generation crosses, local varieties and wild species instead to significantly shorten study cycles and improve work efficiency [21]. Secondly, its use of diverse sources for inheritance variation allows simultaneous detection of multiple alleles at the same position with increased potential applicability across different breeding backgrounds. Thirdly, because complexity exists in interspecific comparisons due to differences in trait expression patterns within various breeds or germplasm pools, one sample may serve several purposes during an assay. Fourthly, higher resolution occurs when utilizing naturally occurring recombination events that are more informative than artificial ones used by breeders. However, there also exist limitations such as low precision caused by limited diversity levels, which leads to reduced accuracy compared to other methods involving marker-assisted selection. Additionally, challenges include complicated sampling histories from variable geographic locations along with concerns about confounding effects related to family structure and environmentally driven influences upon gene frequencies. Furthermore, rare variant discovery suffers lower sensitivity rates, resulting in possible loss of important heritable variations. Combining both methods enables effective utilization of their respective strengths while addressing deficiencies encountered throughout each process, ultimately improving overall effectiveness regarding efficient identification of key hereditary factors underlying multifactorial characteristics [22,23]. Moreover, GWAS can be performed much faster and at lower cost because it bypasses the time of developing biparental populations [24,25]. GWAS has been used to conduct genetic analysis of a wide range of agronomic traits and resistance to diseases [25,26], grain processing and end-use quality [27,28], tolerance to abiotic stress [29], and yield-related traits [30,31].
The Yellow and Huai River Valleys Facultative Wheat Region (YHVFWR) is the largest wheat production region in China. Breeding cultivars with superior processing quality could be greatly enhanced using markers developed from single-nucleotide polymorphisms (SNPs). In this study, SNPs associated with processing quality in a panel of 165 elite wheat accessions mainly from the YHVFWR were used to (1) dissect the genetic architecture of mixograph-related traits, (2) identify SNPs significantly associated with mixograp-related traits and (3) search for candidate mixograph-related traits genes for further study.

2. Materials and Methods

2.1. Plant Materials and Field Trials

In the current investigation, a comprehensive collection consisting of 165 diversified varieties was utilized; specifically, they were comprised of 143 germplasms sourced from both the Yellow and Huai River Valley Facultative Wheat Region located across mainland China while additionally integrating another 22 samples originating from different nations, such as Italy (9); Argentina (7); Japan (4); Australia (1); along with one sample collected each from Turkish territory. These hexaploid wheats represented all types found worldwide based upon their geographical origin or agronomic performance. Furthermore, every single variety had already received permission before being stored into the national seed bank maintained by the Institute of Plant Industry affiliated with the Crop Science Society of China.
The 165 accessions used in GWAS to identify the loci for quality traits assessed by mixograph were grown at Suixi in Anhui province and Anyang in Henan in the 2012–2013 and 2013–2014 cropping seasons. Field trials were conducted in randomized complete blocks design (RCBD). Agronomic management followed local practices. Each plot contained three 2 m rows spaced 20 cm apart, and there were 3 replications.

2.2. Mixograph-Related Traits Evaluated and Statistical Analyses

Clean samples of at least 300 grains were tempered overnight to the 14%, 15% and 16% moisture contents normally used to mill soft, medium and hard types, respectively. All samples were milled at 60% flour extraction using a Brabender Quadrumat Junior Mill (Brabender Inc., Duisberg, Germany). Mixographs have several advantages, including small sample amounts (usually 10 g, minimum 2 g), high daily processing samples, simple operation and rich curve information. They are increasingly widely used in research on the rheological properties of dough due to their excellent performance characteristics. Mixograph measurements were taken with 10 g of flour per sample on a 14% moisture basis using the National Manufacturing Mixograph (National Manufacturing, TMCO Division, Lincoln, NE, USA), according to the AACC (2000) method 54-40A. BLUP across environments was analyzed using the PROCMIXED function in SAS v9.3.

2.3. Genotyping and Population Structure

Cultivars were genotyped by the 90 K SNP and 660 K SNP arrays by CapitalBio (Beijing, China). The SNP chip genotyping procedure involves six steps: (1) sample whole-genome amplification at approximately 1000 folds; (2) fragmentation treatment on the diluted PCR product; (3) hybridization between segmented DNAs and chips; (4) single base extension reactions at specific loci; (5) staining; and (6) scanning imaging. Due to common wheat’s hexaploid nature, we use combinations of GeneStudio v2011.1 and GeneStudio Polyploid Clustering V1.0 for genetic typing. First, raw image scans from gene expressions were read out through GeneStudio v2011.1; then, clustering based on ploidy level was performed via GeneStudio Polyploid Clustering V1.0. Filter criteria followed four standards: (1) eliminating markers without differences among parents; (2) homozygous alleles assumed missing data; (3) removing marker datasets where more than 10 percentage points of values have been removed; and (4) one particular variant site should meet either condition that either one or two types account for no greater than 0.7 or equal to or larger than 0.3, respectively.
The Chinese Spring (IWGSC v1.0) reference genome was used for GWAS. Population structure, principal components analysis (PCA), NJ-tree and LD decay analysis were reported in a previous study [25]. Population structure was assessed utilizing 2000 polymorphic SNP markers sourced from the 660 K SNP arrays. The analysis was conducted using Structure v2.3.4 (http://pritchardlab.stanford.edu/struc-ture.html) (accessed on 5 July 2022). For each K value ranging from 2 to 12, five independent runs were executed employing an admixture model. Each run consisted of 100,000 Markov Chain iterations that were recorded, preceded by 10,000 burn-in periods. To anticipate the actual count of subpopulations, an ad hoc quantity statistic denoted as ΔK, reliant on the rate of logarithmic probability alteration between consecutive K values, was employed [25]. Broad-sense heritability (hb2) of mixograph-related traits was calculated as hb2 = σg2/(σg2 + σge2/r +σε2/re), where σge2, σε2 and σg2 mean the genotype × environment interaction, residual error variances and genotype, respectively. Of these, e and r were the number of environments and the number of replicates per environment, respectively.

2.4. Association Mapping and the Identification of Candidate Genes

The mixed linear model (MLM, PCA + K) was used to avoid the spurious marker–trait associations (MTAs) by Tassel v5.0 [27]. Both the kinship matrix and PCA were estimated by Tassel v5.0. Markers with an adjusted −log10 (p-value) ≥ 3.0 were regarded as MTAs because Bonferroni–Holm correction was too conservative. Manhattan and Q-Q plots were drawn by the CMplot (R 3.6.5).
Candidate genes associated with loci consistently identified across two or more environments were pinpointed. The ensuing procedures were undertaken to ascertain candidate genes for noteworthy or steadfast quantitative trait loci (QTL). Firstly, a thorough search was conducted to retrieve all genes located within the linkage disequilibrium (LD) block vicinity surrounding the peak single-nucleotide polymorphism (SNP) (within a ± 3.0 Mb range, based on prior LD decay analysis) of each significant QTL from the IWGSC V1.0 dataset. Subsequently, all accessible SNPs located within these genes were scrutinized. Genes (excluding those encoding hypothetical proteins, transposon proteins and retrotransposon proteins) harboring SNPs within coding regions, with the potential to induce missense mutations, were designated as candidate genes. Given the substantial regulation of processing quality traits by diverse phytohormones, as well as factors such as glycolysis, signal transduction and cell growth, genes participating in these pathways were classified as high-confidence candidate genes for processing quality traits. Flanking sequences of significantly associated SNPs (including the LD decay interval of peak markers around 3.0 Mb) were used in BLASTx against the NCBI database and reference genome annotations from IWGSC v1.0 was used to predict candidate genes.

3. Results

3.1. Genotyping and Population Structure Analysis

After filtering, 259,922 SNPs were used in GWAS of mixograph-related traits [25]. Population structure, neighbor-joining (NJ) tree and PCA analysis identified three subgroups of accessions [25]; Subgroup I mainly originated from Shandong; Subgroup II included cultivars mainly from Henan, Anhui and Shaanxi; and most Subgroup III cultivars were from Henan (Figure S2) [25]. LD decay for the whole genome was about 8 Mb, with the D genome 11 Mb, the A genome 6 Mb and the B genome 4 Mb (Figure S3) [25]. The SNP density across the whole genome was about 18.5 SNPs/Mb [25].

3.2. Phenotypic Evaluation

All mixograph-related traits showed continuous variation (Table S1; Figure S1). The mean values of MPT, MPV, MPW and MTxW were 3.16 (1.52–6.49), 47.84 (36.3–63.76), 18.32 (10.19–26.75) and 6.29 (2.5–17.67), respectively. ANOVA revealed significant effects (p <0.01) of genotypes, environments and genotype × environment interactions on each processing-related quality trait (Table 1). Broad-sense heritability (h2) estimated for MPW, MPV, MPT and MTxW was 0.75, 0.72, 0.69 and 0.71, respectively. MPT and MPW showed significant (p < 0.01) and positive correlations with MTxW (r = 0.885 and 0.448), and the MPV showed a significant (p < 0.01) and positive correlation with MPW (r = 0.867).

3.3. Genome-Wide Association Studies

Twenty-four loci associated with mixograph-related traits were detected. There were nine loci for MPT distributed across chromosomes 1A, 1B, 1D, 2B, 3A, 3D, 4A, 5D and 7D explaining 10.6–42.5% of the phenotypic variances. A single locus affecting MPV was located on chromosome 5D, explaining 13.4–16.8% of the phenotypic variance. Eight loci for MPW on chromosomes 1B (2), 2A, 3B, 4A, 4B, 6A and 7B explained 10.2–15.8% of the phenotypic variance. Twelve loci for MTxW were identified on chromosomes 1A, 1B (2), 1D, 2B, 3D (2), 4A (2), 7B and 7D (2), accounting for 10.5–27.3% of the phenotypic variances. The 1A (506.9 Mb), 1B (553.6 Mb), 1D (407.9–416.5 Mb), 3D (191.1 Mb) and 7D (321.0 Mb) loci showed pleiotropic effects on MPT and MTxW, whereas the 4A locus (610.1–621.6 Mb) controlled both MPW and MPT (Table 2 and Table S2, Figure 1 and Figure S4).

3.4. Candidate Genes

Seven candidate genes were identified for wheat progressing quality-related traits. Three 3-ketoacyl-CoA synthase genes (TraesCS2B01G535700, TraesCS3A01G451100 and TraesCS3D01G444000) were identified on chromosomes 2B (731.5 Mb), 3A (689.5 Mb) and 3D (553.2 Mb), and another gene encoding lipoxygenase (TraesCS4A01G359800) was detected in the QTL (632.9 Mb) on chromosome 4A. Candidate genes for plant lipid transfer protein/Par allergen (TraesCS4A01G021300), pyridoxal phosphate-dependent decarboxylase (TraesCS1A01G329500) and sucrose synthase 3 (TraesCS3D01G184500) were identified on chromosomes 4A (14.8 Mb), 1A (518.5 Mb) and 3D (170.2 Mb), respectively (Table 3 and Table S3).

4. Discussion

As a vital staple crop and an engine for advancing agricultural high-quality development, particularly with regard to enhancing wheat quality, there exists an increased urgency towards such efforts. Over the past several years, advancements have been achieved by means of genetic modifications within the context of complicated wheat quality characteristics. Nonetheless, owing to the challenges posed by the lack of accurate trait diagnosis methods when it comes to wheat quality, studies on its underlying genetics are moving at a snail’s pace, further impeding the overall wheat quality breeding processes. To address these issues, researchers should explore the highly correlated marker systems linked to various aspects of wheat quality via genealogical analyses so as to determine relevant regions that could be employed later during wheat quality breeding—ultimately speeding up the whole selection procedure.
The rheological properties of dough are an essential quality trait for wheat, determining not only its own processing qualities but also those of breads and cooked products made using it. Currently, various instruments such as the farinograph, rheometer and mixograph are primarily used to assess these traits. However, due to the relatively high quantity of flour needed (normally exceeding 100 g) and the lengthy testing times involved, these devices prove difficult to apply when confronted with numerous samples during breeding generations. In contrast, the mixograph calls for less material and features shorter test times than the aforementioned equipment. Furthermore, it displays highly significant correlations with key parameters measured by means of the farinograph and rheometer, while simultaneously demonstrating a strong connection with actual kneading times during baking. Consequently, this tool is presently extensively implemented in both domestic and international wheat breeding initiatives [13,14].
In this study, most of the 165 accessions mainly from Shandong and foreign cultivars were classified into subgroup 1; subgroup 2 consisted of 54 accessions, mainly from Henan, Anhui and Shaanxi provinces, whereas subgroup 3 mainly included the accessions from Henan province. Thus, the MLM model with population structure and kinship matrix settings was performed in this study to avoid spurious results [22]. The LD decay distance for the whole genome of about 8 Mb indicated that the marker density was adequate for the further association analysis [25].

4.1. Comparison with QTLs or Genes in Previous Studies

Dough strength constitutes a prototypical quantitative trait under the influence of multiple minor genes. Numerous QTLs associated with mixograph-related traits have been previously documented [13,14,15,16,17]. However, those studies were mainly based on simple sequence repeat (SSR) markers and not tightly associated with QTLs. In this study, association analysis of mixograph-related traits was performed using high-density SNP arrays. Barakat et al. [2] reported 108 loci for farinograph- and mixograph-related traits on all 21 wheat chromosomes in two double haploid (DH) populations. Three loci on chromosomes 1B, 2B and 3B overlapped present loci on chromosomes 1B (13.5–14.5 Mb) for MTxW, 2B (738.2–752.8 Mb) for MPT and 3B (561.0–579.3 Mb) for MPW. Zhang-Biehn et al. [28] have reported 11 loci for mixograph-related traits on chromosomes 1A, 1B, 1D, 2A, 5B, 6A and 6D in a panel of 462 advanced breeding lines; three overlapped SNPs were present in 1B locus (553.6 Mb) for MPT and MTxW, 1D (407.9–416.5 Mb) for MPT and MTxW, and 6A (23.4–26.9 Mb) for MPW. The locus for MPT and MTxW on chromosome 1D (407.9–416.5 Mb) was in the same position as Qmlt.tamu.1d identified by Yu et al. [32], who also reported Qgcd.tamu.7B, which overlapped with the present locus for MPW on chromosome 7B (551.9–557.4 Mb). Zhou et al. [33] reported two loci for mixograph-related traits on chromosomes 4A and 7B that coincided with the present 4A (621.2–667.7 Mb) for MPT and 7B (321.0 Mb) loci for MPT and MTxW. However, due to the fact that most genetic analyses of wheat quality traits have been conducted using traditional SSR and DArT markers, which lack specific physical locations and cannot be directly compared with the SNP information employed in this study, we are limited to the following inferences: among the 24 loci for mixograph-related traits, nine were reported previously; the remaining loci may be novel.

4.2. Candidate Gene Analysis

Candidate gene TraesCS1A01G329500 on chromosome 1A encoding a pyridoxal phosphate-dependent decarboxylase involved in nitrogen metabolism affects the amount and composition of dough development and gluten matrix [34]. TraesCS4A01G021300 on chromosome 4A in an LD decay with a nonspecific lipid transfer protein (LTP) gene promotes the inter-membrane transfer of lipids that plays an important role in dough mixing. Lipids not only influence the dough capacity, but also interact with glutenin and gliadin proteins [35]. TraesCS3D01G184500 on chromosome 3D encodes sucrose phosphate synthase III (SPS III), which plays an important role in the sucrose metabolic pathway in plants [36,37]. SPS III is an important enzyme involved in the metabolism of sugars in plant cells. Its main function is to convert glucose and fructose into sucrose and, at the same time, produce a phosphate group. Fructose 6-phosphate and UDP glucose are catalyzed by sucrose phosphate synthase to form sucrose 6-phosphate, which is catalyzed by sucrose phosphatase to form sucrose. Within the LD decay of a locus on chromosome 4A (632.9 Mb), there was (TraesCS4A01G359800, LOX 2) encoding lipoxygenase 2 that catalyzes bound fatty acids [38,39]. Previous studies showed that lipids form complexes with amylose and LOX 2 acts on the complex formed by starch and lipids to increase the solubility of starch [40]. TraesCS2B01G535700 on chromosome 2B (731.5 Mb), TraesCS3A01G451100 on chromosome 3A (689.5 Mb) and TraesCS3D01G444000 on chromosome 3D (553.2 Mb) encode 3-ketoyl COA synthases involved in the fatty acid biosynthesis pathway [41]. Previous studies indicated that interaction between lipids and proteins has a significant impact on the rheological properties of dough. Further experiments are needed to verify the functions of the candidate genes due to the complexity of metabolic pathways and genetic backgrounds.
While conventional wheat breeding practices have contributed to the enhancement of processing quality, the efficacy of early-generation selection remains limited. Additive effects have been discerned between traits associated with processing quality and advantageous alleles recognized through genome sequence variations. The accumulation of multiple favorable alleles, such as those pinpointed in this investigation, holds the potential to enhance processing quality-linked traits. Loci with stable effects across environments should be more appliable in MAS, such as those on chromosomes1A (506.9 Mb), 1B (553.6 Mb), 1D (407.9–416.5 Mb) and 7D (321.0 Mb), showing pleiotropic effects on MPT and MTxW, along with the 4A (610.1–621.6 Mb) locus with pleiotropic effects on MPW and MPT. The loci that have been confirmed through conventional linkage mapping or GWAS in prior research could serve as subjects for subsequent investigations. Finally, accessions with superior processing quality-related traits and larger numbers of favorable alleles, such as Aca 601, Klein Jaba l1, Shanyou 225, Jishi 02-1, Sagittario, Klein Flecha, Wanmai 33, Sunong6, Libero Mantol, ProINTA Colibr 1 and Nidera Baguette 20, are recommended as parental lines for breeding cultivars with superior quality-related traits by MAS (Table 4).

5. Conclusions

A GWAS for mixograph-related traits was conducted on a panel of 165 varieties using the wheat 90K and 660K SNP arrays. In total, 15 of the 24 loci identified were novel. Seven candidate genes for mixograph-related traits were predicted. SNPs in genes associated with favorable mixograph traits can be converted to selection markers that can be used in the early generations of breeding.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes14091816/s1, Figure S1. The histogram for the progressing quality-related traits in 165 wheat accessions. MPT: mixograph midline peak time; MPV: mixograph midline peak value; MPW: mixograph midline peak width; MTxW: mixograph midline 8 min band width. 13SX, 13AY, 14SX, 14AY and BLUP indicate Suixi 2013, Anyang 2013, Suixi 2014, Anyang 2014 and the best linear unbiased prediction (BLUP), respectively. Figure S2. Population structure for the 165 wheat accessions [25]. a Delta K for structure analysis; b population structure analysis; c neighbor-joining (NJ) tree; d principal components analysis (PCA) plots. Figure S3. LD decay for the 165 wheat accessions LD [25]. Figure S4. Q-Q plot for the progressing quality-related traits in 165 wheat accessions. MPT: mixograph midline peak time; MPV: mixograph midline peak value; MPW: mixograph midline peak width; MTxW: mixograph midline 8 min band width. 13SX, 13AY, 14SX, 14AY and BLUP indicate Suixi 2013, Anyang 2013, Suixi 2014, Anyang 2014 and the best linear unbiased prediction (BLUP), respectively. Table S1. The details for mixograph-related traits in 165 wheat accessions. MPT: mixograph midline peak time; MPV: mixograph midline peak value; MPW: mixograph midline peak width; MTxW: mixograph midline 8 min band width. E1, E2, E3, E4 and E5 indicate Suixi 2013, Anyang 2013, Suixi 2014, Anyang 2014 and the best linear unbiased prediction (BLUP), respectively. Table S2. Marker–trait associations for mixograph-related traits in 165 wheat accessions. Table S3. GO annotation for the candidate genes.

Author Contributions

H.J. carried out the experiment and wrote the paper. Y.T., J.L., Y.Z., H.Z. and R.Z. participated in field trials. X.Y., X.S., Y.D., Y.W. and Q.G. contributed to flour milling. J.Z. and Z.H. designed the experiment and assisted in writing the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Natural Science Foundation of Heilongjiang (YQ2022C030), Heilongjiang Academy of Agricultural Sciences Project (2020FJZX011), the Heilongjiang Scientific Research Institution Foundation (CZKYF2021C001) and Establishment of Joint Chinese–Bulgarian Laboratory for Molecular Biology of Crop Germplasm Resources (KY201901009).

Institutional Review Board Statement

We declare that these experiments complied with the ethical standards in China.

Informed Consent Statement

Not applicable.

Data Availability Statement

All datasets generated for this study are included in the article/Supplementary Material; further inquiries can be directed to the first author.

Conflicts of Interest

The authors declare that they have no competing interests.

References

  1. Jin, H.; Zhang, Y.; Li, G.; Mu, P.; Fan, Z.; Xia, X.; He, Z. Effects of allelic variation of HMW-GS and LMW-GS on mixograph properties and Chinese noodle and steamed bread qualities in a set of Aroona near-isogenic wheat lines. J. Cereal Sci. 2013, 57, 146–152. [Google Scholar] [CrossRef]
  2. Barakat, M.; Al-Doss, A.; Moustafa, K.; Motawei, M.; Alamri, M.; Mergoum, M.; Sallam, M.; Al-Ashkar, I. QTL analysis of farinograph and mixograph related traits in spring wheat under heat stress conditions. Mol. Biol. Rep. 2020, 47, 5477–5486. [Google Scholar] [CrossRef] [PubMed]
  3. Kaur, A.; Singh, N.; Kaur, S.; Ahlawat, A.K.; Singh, A.M. Relationships of flour solvent retention capacity, secondary structure and rheological properties with the cookie making characteristics of wheat cultivars. Food Chem. 2014, 158, 48–55. [Google Scholar] [CrossRef] [PubMed]
  4. Boehm, J.D.; Ibba, M.I.; Kiszonas, A.M.; Morris, C.F. End-use quality of CIMMYT-derived soft-kernel durum wheat germplasm: II. Dough strength and pan bread quality. Crop Sci. 2017, 57, 1485–1494. [Google Scholar] [CrossRef]
  5. Pu, H.; Wei, J.; Wang, L.; Huang, J.; Chen, X.; Luo, C.; Liu, S.; Zhang, H. Effects of potato/wheat flours ratio on mixing properties of dough and quality of noodles. J. Cereal Sci. 2017, 76, 236–242. [Google Scholar] [CrossRef]
  6. Liu, Y.; He, Z.; Appels, R.; Xia, X. Functional markers in wheat: Current status and future prospects. Theor. Appl. Genet. 2012, 125, 1–10. [Google Scholar] [CrossRef]
  7. He, Z.H.; Zhuang, Q.S.; Cheng, S.H.; Yu, Z.W.; Zhao, Z.D.; Liu, X. Wheat production and technology improvement in China. J. Agric. 2018, 8, 107, (In Chinese with English abstract). [Google Scholar]
  8. Rasheed, A.; Hao, Y.; Xia, X.; Khan, A.; Xu, Y.; Varshney, R.K.; He, Z. Crop breeding chips and genotyping platforms: Progress, challenges, and perspectives. Mol. Plant 2017, 10, 1047–1064. [Google Scholar] [CrossRef]
  9. Prashant, R.; Mani, E.; Rai, R.; Gupta, R.K.; Tiwari, R.; Dholakia, B.; Oak, M.; Röder, M.; Kadoo, N.; Gupta, V. Genotype × environment interactions and QTL clusters underlying dough rheology traits in Triticum aestivum L. J. Cereal Sci. 2015, 64, 82–91. [Google Scholar] [CrossRef]
  10. Sakr, N.; Rhazi, L.; Aussenac, T. Bread wheat quality under limiting environmental conditions: II—Rheological properties of Lebanese wheat genotypes. J. Saudi Soc. Agric. Sci. 2021, 20, 235–242. [Google Scholar] [CrossRef]
  11. Williams, R.M.; O’Brien, L.; Eagles, H.A.; Solah, V.A.; Jayasena, V. The influences of genotype, environment, and genotype × environment interaction on wheat quality. Aust. J. Agr. Res. 2008, 59, 95–111. [Google Scholar] [CrossRef]
  12. Cappelli, A.; Oliva, N.; Cini, E. Stone milling versus roller milling: A systematic review of the effects on wheat flour quality, dough rheology, and bread characteristics. Trends Food Sci. Technol. 2020, 97, 147–155. [Google Scholar] [CrossRef]
  13. McCartney, C.A.; Somers, D.J.; Lukow, O.; Ames, N.; Noll, J.; Cloutier, S.; Humphreys, D.G.; McCallum, B.D. QTL analysis of quality traits in the spring wheat cross RL4452 × ‘AC Domain’. Plant Breed. 2006, 125, 565–575. [Google Scholar] [CrossRef]
  14. Zhang, Y.; Wu, Y.; Xiao, Y.; Yan, J.; Zhang, Y.; Zhang, Y.; Ma, C.; Xia, X.; He, Z. QTL mapping for milling, gluten quality, and flour pasting properties in a recombinant inbred line population derived from a Chinese soft × hard wheat cross. Crop Pasture Sci. 2009, 60, 587–597. [Google Scholar] [CrossRef]
  15. Kerfal, S.; Giraldo, P.; Rodríguez-Quijano, M.; Vázquez, J.F.; Adams, K.; Lukow, O.M.; Röder, M.S.; Somers, D.J.; Carrillo, J.M. Mapping quantitative trait loci (QTLs) associated with dough quality in a soft × hard bread wheat progeny. J. Cereal Sci. 2010, 52, 46–52. [Google Scholar] [CrossRef]
  16. Zheng, F.F.; Deng, Z.Y.; Shi, C.L.; Zhang, X.Y.; Tian, J.C. QTL mapping for dough mixing characteristics in a recombinant inbred population derived from a waxy × strong gluten wheat (Triticum aestivum L.). J. Integr. Agr. 2013, 12, 951–961. [Google Scholar] [CrossRef]
  17. Echeverry-Solarte, M.; Kumar, A.; Kianian, S.; Simsek, S.; Alamri, M.S.; Mantovani, E.E.; McClean, P.E.; Deckard, E.L.; Elias, E.; Schatz, B.; 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. [Google Scholar] [CrossRef]
  18. Wang, S.; Wong, D.; Forrest, K.; Allen, A.; Chao, S.; Huang, B.E.; Maccaferri, M.; Salvi, S.; Milner, S.G.; Cattivelli, L.; et al. Characterization of polyploid wheat genomic diversity using a high-density 90,000 single nucleotide polymorphism array. Plant Biotechnol. J. 2014, 12, 787–796. [Google Scholar] [CrossRef]
  19. Sun, C.; Dong, Z.; Zhao, L.; Ren, Y.; Zhang, N.; Chen, F. The Wheat 660K SNP array demonstrates great potential for marker-assisted selection in polyploid wheat. Plant Biotechnol. J. 2020, 18, 1354–1360. [Google Scholar] [CrossRef]
  20. Huang, S.; Wu, J.; Wang, X.; Mu, J.; Xu, Z.; Zeng, Q.; Liu, S.; Wang, Q.; Kang, Z.; Han, D. Utilization of the genomewide wheat 55K SNP array for genetic analysis of stripe rust resistance in common wheat line P9936. Phytopathology 2019, 109, 819–827. [Google Scholar] [CrossRef]
  21. Xu, Y.; Crouch, J.H. Marker-assisted selection in plant breeding: From publications to practice. Crop Sci. 2008, 48, 391–407. [Google Scholar] [CrossRef]
  22. Zhu, C.; Gore, M.; Buckler, E.S.; Yu, J. Status and prospects of association mapping in plants. Plant Genome 2008, 1, 5–20. [Google Scholar] [CrossRef]
  23. Flint-Garcia, S.A.; Thornsberry, J.M.; Buckler, E.S. Structure of linkage disequilibrium in plants. Annu. Rev. Plant Biol. 2003, 54, 357–374. [Google Scholar] [CrossRef] [PubMed]
  24. Breseghello, F.; Sorrells, M.E. Association mapping of kernel size and milling quality in wheat (Triticum aestivum L.) cultivars. Genetics 2006, 172, 1165–1177. [Google Scholar] [CrossRef] [PubMed]
  25. Liu, J.; He, Z.; Rasheed, A.; Wen, W.; Yan, J.; Zhang, P.; Wan, Y.; Zhang, Y.; Xie, C.; Xia, X. Genome-wide association mapping of black point reaction in common wheat (Triticum aestivum L.). BMC Plant Biol. 2017, 17, 220. [Google Scholar] [CrossRef] [PubMed]
  26. Du, X.; Xu, W.; Peng, C.; Li, C.; Zhang, Y.; Hu, L. Identification and validation of a novel locus, Qpm-3BL, for adult plant resistance to powdery mildew in wheat using multilocus GWAS. BMC Plant Biol. 2021, 21, 357. [Google Scholar] [CrossRef] [PubMed]
  27. Battenfield, S.D.; Sheridan, J.L.; Silva, L.D.C.E.; Miclaus, K.J.; Dreisigacker, S.; Wolfinger, R.D.; Peña, R.J.; Singh, R.P.; Jackson, E.W.; Fritz, A.K.; et al. Breeding-assisted genomics: Applying meta-GWAS for milling and baking quality in CIMMYT wheat breeding program. PLoS ONE 2018, 13, e0204757. [Google Scholar] [CrossRef]
  28. Zhang-Biehn, S.; Fritz, A.K.; Zhang, G.; Evers, B.; Regan, R.; Poland, J. Accelerating wheat breeding for end-use quality through association mapping and multivariate genomic prediction. Plant Genome 2021, 14, e20164. [Google Scholar] [CrossRef]
  29. Thoen, M.P.M.; Olivas, N.H.D.; Kloth, K.J.; Coolen, S.; Huang, P.; Aarts, M.G.M.; Bac-Molenaar, J.A.; Bakker, J.; Bouwmeester, H.J.; Broekgaarden, C.; et al. Genetic architecture of plant stress resistance: Multi-trait genome-wide association mapping. New Phytol. 2017, 213, 1346–1362. [Google Scholar] [CrossRef]
  30. Ma, F.; Xu, Y.; Ma, Z.; Li, L.; An, D. Genome-wide association and validation of key loci for yield-related traits in wheat founder parent Xiaoyan 6. Mol. Breed. 2018, 38, 1–15. [Google Scholar] [CrossRef]
  31. Yang, Y.; Amo, A.; Wei, D.; Chai, Y.; Zheng, J.; Qiao, P.; Cui, C.; Lu, S.; Chen, L.; Hu, Y.G.; et al. Large-scale integration of meta-QTL and genome-wide association study discovers the genomic regions and candidate genes for yield and yield-related traits in bread wheat. Theor. Appl. Genet. 2021, 134, 3083–3109. [Google Scholar] [CrossRef] [PubMed]
  32. Yu, S.; Assanga, S.O.; Awika, J.M.; Ibrahim, A.M.H.; Rudd, J.C.; Xue, Q.; Guttieri, M.J.; Zhang, G.; Baker, J.A.; Jessup, K.E.; et al. Genetic mapping of quantitative trait loci for end-use quality and grain minerals in hard red winter wheat. Agronomy 2021, 11, 2519. [Google Scholar] [CrossRef]
  33. Zhou, Z.; Zhang, Z.; Jia, L.; Qiu, H.; Guan, H.; Liu, C.; Qin, M.; Wang, Y.; Li, W.; Yao, W.; et al. Genetic basis of gluten aggregation properties in wheat (Triticum aestivum L.) dissected by QTL mapping of glutopeak parameters. Front. Plant Sci. 2021, 11, 611605. [Google Scholar] [CrossRef] [PubMed]
  34. Kumar, R. Evolutionary trails of plant group II pyridoxal phosphate-dependent decarboxylase genes. Front. Plant Sci. 2016, 7, 1268. [Google Scholar] [CrossRef] [PubMed]
  35. Douliez, J.P.; Michon, T.; Elmorjani, K.; Marion, D. Mini review: Structure, biological and technological functions of lipid transfer proteins and indolines, the major lipid binding proteins from cereal kernels. J. Cereal Sci. 2000, 3, 1–20. [Google Scholar] [CrossRef]
  36. Castleden, C.K.; Aoki, N.; Gillespie, V.J.; MacRae, E.A.; Quick, W.P.; Buchner, P.; Foyer, C.H.; Furbank, R.T.; Lunn, J.E. Evolution and function of the sucrose-phosphate synthase gene families in wheat and other grasses. Plant Physiol. 2004, 135, 1753–1764. [Google Scholar] [CrossRef] [PubMed]
  37. Wang, L.; Cui, N.; Zhang, K.Y.; Fan, H.Y.; Li, T.L. Research advance of sucrose phosphate synthase (SPS) in higher plants. Int. J. Agric. Biol. 2013, 15, 1221–1226. [Google Scholar]
  38. Borrelli, G.M.; Troccoli, A.; Di Fonzo, N.; Fares, C. Durum wheat lipoxygenase activity and other quality parameters that affect pasta color. Cereal Chem. 1999, 76, 335–340. [Google Scholar] [CrossRef]
  39. Mousavi, B.; Kadivar, M. Effect of brine solution as a wheat conditioner, on lipase, amylase, and lipoxygenase activities in flour and its corresponding dough rheological properties. J. Food Process Pres. 2018, 42, e13631. [Google Scholar] [CrossRef]
  40. Bahal, G.; Sudha, M.L.; Ramasarma, P.R. Wheat germ lipoxygenase: Its effect on dough rheology, microstructure, and bread making quality. Int. J. Food Prop. 2013, 16, 1730–1739. [Google Scholar] [CrossRef]
  41. Jing, Y.F.; Wei, H.X.; Liu, F.F.; Liu, Y.F.; Zhou, L.; Liu, J.F.; Yang, S.Z.; Zhang, H.Z.; Mu, B.Z. Genetic engineering of the branched-chain fatty acid biosynthesis pathway to enhance surfactin production from Bacillus subtilis. Biotechnol. Appl. Biochem. 2023, 70, 238–248. [Google Scholar] [CrossRef]
Figure 1. Manhattan plots for mixograph-related traits in 165 wheat accessions by the mixed linear model (MLM) in Tassel v5.0. MPT: mixograph midline peak time; MPV: mixograph midline peak value; MPW: mixograph midline peak width; MTxW: mixograph midline 8 min band width.13SX, 13AY, 14SX, 14AY and BLUP indicate Suixi 2013, Anyang 2013, Suixi 2014, Anyang 2014 and the best linear unbiased prediction (BLUP), respectively.
Figure 1. Manhattan plots for mixograph-related traits in 165 wheat accessions by the mixed linear model (MLM) in Tassel v5.0. MPT: mixograph midline peak time; MPV: mixograph midline peak value; MPW: mixograph midline peak width; MTxW: mixograph midline 8 min band width.13SX, 13AY, 14SX, 14AY and BLUP indicate Suixi 2013, Anyang 2013, Suixi 2014, Anyang 2014 and the best linear unbiased prediction (BLUP), respectively.
Genes 14 01816 g001
Table 1. ANOVA for mixograph-related traits in a panel of 165 wheat accessions.
Table 1. ANOVA for mixograph-related traits in a panel of 165 wheat accessions.
Source of VariationdfMS
MPWMPVMPTMTxW
Genotypes1641.88 **46.3 **32.5 **20.3 **
Environments30.42 **295.0 **896.3 **56.4 **
Replicates (nested in environments)20.15 **6.1 **7.1 **3.2 **
Genotype*Environment9830.13 **3.8 **5.9 **2.1 **
Error1425
** significant at p = 0.01. MPW: mixograph midline peak width; MPV: mixograph midline peak value; MPT: mixograph midline peak time; MTxW: mixograph midline 8 min band width.
Table 2. Mixograph-related traits identified in the wheat accession panel by association analysis.
Table 2. Mixograph-related traits identified in the wheat accession panel by association analysis.
LociTraitChr.Start (Mb)R2p-ValueEnvironmentsFavorable AlleleReference
MinMaxMinMax
qM1MPT1A506.914.90%23.70%6.70 × 10−81.90 × 10−5E1; E3; E4; E5C
MTxW1A506.913.80%15.90%9.40 × 10−63.90 × 10−5E1; E3; E5T
qM2MTxW1B13.5–14.510.80%17.00%7.20 × 10−76.40 × 10−5E1; E4; E5T
qM3MPW1B130.610.30%12.10%1.70 × 10−57.40 × 10−5E1; E4; E5A[2]
MPT1B553.631.50%37.70%9.10 × 10−124.40 × 10−10E1; E3; E4; E5G
qM4MTxW1B553.614.80%27.10%1.90 × 10−93.70 × 10−6E1; E2; E3; E4; E5G
qM5MPW1B673.4–674.512.70%13.90%3.70 × 10−58.10 × 10−5E2A[28]
qM6MPT1D407.9–416.512.60%42.50%1.40 × 10−139.80 × 10−5E1; E2; E3; E4; E5G[32]
MTxW1D407.9–416.510.80%27.30%8.30 × 10−109.40 × 10−5E1; E2; E3; E4; E5G
qM7MPW2A191.9–199.612.50%13.80%3.50 × 10−59.00 × 10−5E2A[28]
qM8MTxW2B4.610.60%11.20%6.30 × 10−59.40 × 10−5E2; E3; E5G[32,33]
qM9MPT2B738.2–752.814.90%18.20%4.70 × 10−62.70 × 10−5E4G[2]
qM10MPT3A709.6–710.712.90%11.60%4.20 × 10−59.70 × 10−5E4G
qM11MPW3B561.0–579.313.00%13.60%4.00 × 10−57.80 × 10−5E2A
qM12MPT3D191.122.70%28.90%3.60 × 10−101.70 × 10−8E1; E3; E4; E5A[33]
MTxW3D191.114.70%17.30%5.60 × 10−73.70 × 10−6E1; E2; E3; E5A
qM13MTxW3D578.410.90%11.70%3.10 × 10−56.90 × 10−5E1; E4; E5A
qM14MTxW4A12.4–12.612.60%13.50%8.20 × 10−69.50 × 10−5E1; E2; E3; E5A
qM15MTxW4A89.9–90.410.50%14.40%4.60 × 10−68.10 × 10−5E1; E5A
qM16MPW4A610.1–621.612.80%14.40%2.60 × 10−59.00 × 10−5E2G[2]
MPT4A621.2–667.714.60%17.60%4.30 × 10−64.90 × 10−5E1; E4; E5G
qM17MPW4B12.9–25.810.20%15.80%1.30 × 10−67.70 × 10−5E2; E3; E5G
qM18MPV5D3.613.40%16.80%8.20 × 10−79.50 × 10−6E1; E3; E4; E5C
qM19MPT5D454.110.60%12.30%2.30 × 10−58.50 × 10−5E1; E3; E5G
qM20MPW6A23.4–26.912.70%14.10%3.70 × 10−59.60 × 10−5E2G
qM21MTxW7B126.614.80%18.20%3.20 × 10−73.20 × 10−6E1; E3; E4; E5C
qM22MPW7B551.9–557.412.50%15.20%2.00 × 10−59.90 × 10−5E2G
qM23MTxW7D57.513.90%22.30%3.40 × 10−76.00 × 10−5E2; E4; E5C
qM24MPT7D32120.80%25.10%3.70 × 10−95.90 × 10−8E1; E3; E4; E5G
MTxW7D32112.40%17.80%4.00 × 10−71.90 × 10−5E1; E2; E3; E5C[33]
MPW: mixograph midline peak width; MPT: mixograph midline peak time; MPV: mixograph midline peak value; MTxW: mixograph midline 8 min band width. E1, E2, E3, E4 and E5 indicate Suixi 2013, Anyang 2013, Suixi 2014, Anyang 2014 and the best linear unbiased prediction (BLUP), respectively.
Table 3. Candidate genes for mixograph-related traits.
Table 3. Candidate genes for mixograph-related traits.
Candidate GeneChromosomePosition (Mb)Annotation
TraesCS1A01G3295001A518.5Pyridoxal phosphate-dependent decarboxylase
TraesCS2B01G5357002B731.5 3-ketoacyl-CoA synthase
TraesCS3A01G4511003A689.5 3-ketoacyl-CoA synthase
TraesCS3D01G4440003D553.2 3-ketoacyl-CoA synthase
TraesCS3D01G1845003D170.2 sucrose synthase 3
TraesCS4A01G0213004A14.8Plant lipid transfer protein/Par allergen
TraesCS4A01G3598004A632.9 Lipoxygenase
Table 4. Accessions could be used for wheat processing quality improvement.
Table 4. Accessions could be used for wheat processing quality improvement.
CultivarBLUP-MPTBLUP-MPVBLUP-MPWBLUP-MTxW
Aca 6015.6 50.0 26.7 17.7
Klein Jabal 16.5 50.4 20.6 17.2
Shanyou 2254.6 54.0 24.4 16.6
Jishi 02-16.0 52.8 22.5 16.0
Sagittario5.1 50.6 23.4 15.6
Klein Flecha5.6 49.5 23.4 15.3
Wanmai 336.3 51.9 21.5 15.1
Sunong 66.2 48.4 18.9 14.8
Libero5.6 44.2 18.4 14.8
Mantol5.5 47.3 19.5 14.5
ProINTA Colibr 15.2 51.4 21.2 14.4
Nidera Baguette 205.0 47.2 21.6 13.4
Barra6.3 44.0 15.1 13.2
Aca 8015.4 48.3 23.8 13.0
Jimai 204.0 50.6 23.5 12.8
Nidera Baguette 105.8 42.4 15.3 12.5
Shanmai 944.6 52.9 22.8 11.8
Kitanokaori4.0 49.4 17.6 11.5
Xinong 979-0054.3 52.5 25.2 11.4
Norin 673.0 54.6 23.5 11.1
Shanmai 5096.1 44.0 16.4 11.0
Zhoumai 265.6 41.6 14.5 10.8
Jinan 173.5 52.6 25.7 10.6
Sunstate4.4 45.3 16.7 10.6
Zhoumai 194.3 49.5 19.5 10.3
Gaocheng 89014.9 52.4 19.2 10.3
Shannong 9813.2 63.8 30.4 10.2
Genio3.5 48.4 23.5 10.1
Shiyou 174.5 47.2 17.7 9.9
Jining 164.6 49.1 18.8 9.7
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jin, H.; Tian, Y.; Zhang, Y.; Zhang, R.; Zhao, H.; Yang, X.; Song, X.; Dimitrov, Y.; Wu, Y.-e.; Gao, Q.; et al. Genome-Wide Association Mapping of Processing Quality Traits in Common Wheat (Triticum aestivum L.). Genes 2023, 14, 1816. https://doi.org/10.3390/genes14091816

AMA Style

Jin H, Tian Y, Zhang Y, Zhang R, Zhao H, Yang X, Song X, Dimitrov Y, Wu Y-e, Gao Q, et al. Genome-Wide Association Mapping of Processing Quality Traits in Common Wheat (Triticum aestivum L.). Genes. 2023; 14(9):1816. https://doi.org/10.3390/genes14091816

Chicago/Turabian Style

Jin, Hui, Yuanyuan Tian, Yan Zhang, Rui Zhang, Haibin Zhao, Xue Yang, Xizhang Song, Yordan Dimitrov, Yu-e Wu, Qiang Gao, and et al. 2023. "Genome-Wide Association Mapping of Processing Quality Traits in Common Wheat (Triticum aestivum L.)" Genes 14, no. 9: 1816. https://doi.org/10.3390/genes14091816

APA Style

Jin, H., Tian, Y., Zhang, Y., Zhang, R., Zhao, H., Yang, X., Song, X., Dimitrov, Y., Wu, Y. -e., Gao, Q., Liu, J., Zhang, J., & He, Z. (2023). Genome-Wide Association Mapping of Processing Quality Traits in Common Wheat (Triticum aestivum L.). Genes, 14(9), 1816. https://doi.org/10.3390/genes14091816

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

Article Metrics

Back to TopTop