Genetic analysis of morphological traits in spring wheat from the Northeast of China by a genome-wide association study

Identification of the gene for agronomic traits is important for the wheat marker-assisted selection (MAS) breeding. To identify the new and stable loci for agronomic traits, including flag leaf length (FLL), flag leaf width (FLW), uppermost internode length (UIL), and plant morphology (PM, including prostrate, semi-prostrate, and erect). A total of 251 spring wheat accessions collected from the Northeast of China were used to conduct genome-wide association study (GWAS) by 55K SNP arrays. A total of 30 loci for morphological traits were detected, and each explained 4.8–17.9% of the phenotypic variations. Of these, 13 loci have been reported by previous studies, and the other 17 are novel. We have identified seven genes involved in the signal transduction, cell-cycle progression, and plant development pathway as candidate genes. This study provides new insights into the genetic basis of morphological traits. The associated SNPs and accessions with more of favorable alleles identified in this study could be used to promote the wheat breeding progresses.

Background

Common wheat (Triticum aestivum L.) is an essential cereal crop and provides nearly 20% of the total caloric input to the global population (He et al., 2010; Wang et al., 2019). Due to the complicated genetic architecture (Gao et al., 2017; Li et al., 2021), the progresses for the improvement of morphological traits are difficult. Heilongjiang and Jilin provinces located at the Northeast of China are the major spring wheat-producing regions (Li et al., 2021). Although the wheat production has been improved largely in the past decades, it is still not enough to meet the needs of the people (He et al., 2010; Xu et al., 2021).

Morphological traits include various traits, such as flag leaf length (FLL), flag leaf with (FLW), uppermost internode length (UIL), and plant morphology (PM, including prostrate, semi-prostrate, and erect) (Li et al., 2018; Li et al., 2021). Both the genetic and environmental factors influenced morphological traits. Marker-assisted selection (MAS) based on molecular marker is an available and effective approach for the further improvement of morphological traits of wheat (He et al., 2010; Rasheed et al., 2016; Wang et al., 2020). However, the effective and reliability of MAS depend on the number and quality of available genes/QTLs and associated markers for target traits. Until now, over 30 genes associated with morphological traits were cloned by homologous cloning or map-based cloning, and over 50 functional markers or kompetitive allele-specific PCR (KASP) markers were developed (Cui et al., 2014; Rasheed et al., 2016; Nadolska-Orczyk et al., 2017; Lozada et al., 2019). Besides, over 120 loci for morphological traits were detected by linkage mapping and association analysis (Gao et al., 2017; Würschum et al., 2017; Rahimi et al., 2019; Pang et al., 2020; Li et al., 2021). However, the loci/gene is still not enough for the improvement of wheat morphological traits. Identifying the novel genes or loci for morphological traits is urgent.

Single nucleotide polymorphism (SNP), insertion and deletion (InDel), simple sequence repeats (SSRs), and diversity array technology (DArT) markers are the mostly used molecular marker for genotyping (Wang et al., 2014; Rasheed et al., 2016). Compared with SSR and DArT markers, SNPs have more abundant and higher coverage. With the development of the gene chip and NGS technology, getting SNP quickly becomes feasible and provides an effective way to identifying genes/QTLs for complex traits (Zhu et al., 2008; Liu et al., 2017; Wang et al., 2020). Now, the 55, 90, and 660K wheat SNP arrays are gradually replacing SSR and DArT markers in genetic analysis and have been widely used in the genetic analysis for yield (Gao et al., 2015; Sun et al., 2017; Beyer et al., 2019; Li et al., 2019), disease resistance (stripe rust, leaf rust, or powdery mildew), end-use quality, procession quality, and biotic or abiotic stress tolerance (drought or flood)-related traits (Liu et al., 2016a; Liu et al., 2017; Valluru et al., 2017; Liu et al., 2019; Li et al., 2021; Quan et al., 2021).

Linkage analysis and association mapping are the two main ways to uncover the genetic mechanism of complex traits (Zhu et al., 2008; Liu et al., 2016b; Liu et al., 2017). Compared with the bi-parental linkage analysis, association mapping is based on natural population (including wild types, landraces, released cultivars, and improved accessions) and offers an effective and reliable approach to uncover the genetic architecture of complex traits (Liu et al., 2017; Wang et al., 2020). Besides, linkage analysis focused on specific traits, whereas the GWAS could be used to analyze all traits based on the same set of genotype data (Zhu et al., 2008). Nowadays, association mapping has been widely used in the genetic analysis of complex traits in wheat, including grain yield-related traits, disease-related traits, and biotic and abiotic stresses (Cui et al., 2014; Zuo et al., 2020; Quan et al., 2021). The rapid, efficient, and accurate genotyping is the basis to conduct GWAS. Thus, the SNP chip provides an effective and feasible way for association mapping.

In this study, 251 spring wheat accessions mainly collected from the Northeast of China (Heilongjiang and Jilin province) were used to 1) identify loci for morphological traits in spring wheat, 2) get new insights into the genetic architecture of target traits, and 3) search for candidate genes for further study.

Materials and methods

Plant materials and field trials

A total of 251 spring wheat accessions from the Northeast of China (Heilongjiang or Jilin province) were collected for the GWAS of morphological traits (Supplementary Table S1). The diverse panel was grown at the Harbin and Keshan experimental station of the Heilongjiang academy of agricultural science in Heilongjiang province during the 2018–2019 and 2019–2020 cropping seasons. A complete randomized block design with three replicates was employed in field trials. For both Harbin and Keshan experimental stations, each plot comprised four 2.0 m rows spaced 20 cm apart, with 40 seeds in each row. Agronomic management was performed according to local practices.

Phenotyping and statistical analysis for flag leaf length, flag leaf width, uppermost internode length, and plant morphology

Four traits related to morphological were conducted in all four environments, including FLL, FLW, UIL, and PM. For FLL and FLW, 10 random flag leaves in each plot at the mid-grain-fill stage were used to measure FLL and FLW, represented by the distance between the base and the tip, and width at the widest point, respectively. UIL is the mean distance between the stem base and the top of spikes excluding awns and the mean length of the uppermost internode. Ten single plants in each plot were randomly selected at physiological maturity for measuring UIL. BLUE for four traits among four environments was calculated by a one-stage approach using the R package sommer (https://cran.r-project.org/web/packages/sommer/index.html) as $\mathit{y}=1_{\mathit{n}}\mu +\mathit{Xg}+\mathit{Ze}+\mathit{v}+\mathit{\epsilon }$, where $\mathit{y}$ is the n-dimensional multi-environment phenotypic records, n is the number of multi-environment phenotypic records, $1_{\mathit{n}}$ is an n-dimensional vector of ones, $\mu$ is the intercept, $\mathit{g}$ and $\mathit{e}$ are the vector of genetic and environment effects, respectively, $\mathit{v}$ is the vector of genotype-by-environment interaction effects, $\mathit{X}$ and $\mathit{Z}$ are the design matrices of $\mathit{g}$ and $\mathit{e}$, respectively, and$\mathit{\epsilon }$ is the random residuals. $\mathit{g}$, $\mathit{e}$, $\mathit{v}$, and $\mathit{\epsilon }$ were all assumed as random effects following normal distributions.

Genotyping, population structure, and linkage disequilibrium

A total of 3000 polymorphic and evenly distributed markers on 21 chromosomes were used to conduct population structure analysis by Structure v2.3.4 (Pritchard et al., 2000) (http://pritchardlab.stanford.edu/structure.html). Besides, principal component analysis (PCA) and neighbor-jointing (NJ) trees were also estimated using the Tassel v5.0 to validate the results of population structure analysis (Breseghello and Sorrells, 2006). The LD decay analysis was calculated for the whole genomes using the full matrix and sliding window options in Tassel v5.0 (Breseghello and Sorrells, 2006). The results of LD decay, population analysis, PCA, and NJ-tree analysis for the 251 spring wheat accessions have been reported by Li et al. (2021). The multi-environment phenotypic data were analyzed by a one-stage approach as $\mathit{y}=1_{\mathit{n}}\mu +\mathit{Xg}+\mathit{Ze}+\mathit{v}+\mathit{\epsilon }$, where $\mathit{y}$ is the n-dimensional multi-environment phenotypic records, n is the number of multi-environment phenotypic records, $1_{\mathit{n}}$ is an n-dimensional vector of ones, $\mu$ is the intercept, $\mathit{g}$ and $\mathit{e}$ are the vector of genetic and environment effects, respectively, $\mathit{v}$ is the vector of genotype-by-environment interaction effects, $\mathit{X}$ and $\mathit{Z}$ are the design matrices of $\mathit{g}$ and $\mathit{e}$, respectively, and $\mathit{\epsilon }$ is the random residuals. $\mathit{g}$, $\mathit{e}$, $\mathit{v}$, and $\mathit{\epsilon }$ were all assumed as random effects following normal distributions. The heritability was estimated using the entry-mean basis formula $h^2=\frac{{\sigma }_g^2}{{\sigma }_g^2+\frac{{\sigma }_v^2}{n_e}+\frac{{\sigma }_{\epsilon }^2}{n_e\times n_r}}$, where ${\sigma }_g^2$, ${\sigma }_v^2$, and ${\sigma }_{\epsilon }^2$ are the variance components of genetic effect, environment effect, and residual, $n_e$ and $n_e$ are the number of environments and number of replicates per accession in each environment, respectively. The one-stage phenotypic analysis was realized using the R package sommer (https://cran.r-project.org/web/packages/sommer/index.html) in statistical software R. We have added the statistical method in the M&M section. The SNP-based heritability was estimated by the GREML-LDMS method based on GCTA (Yang et al., 2011) (https://yanglab.westlake.edu.cn/software/gcta/) software.

Genome-wide association study and the identification of candidate genes

To eliminate the spurious marker-trait associations (MTAs) caused by environment variation, the mixed linear model (MLM, PCA + K) in Tassel v5.0 (Breseghello and Sorrells, 2006) was used as follows: y = μ + xβ + u + e. Of these, y is the vector of phenotype; µ is the mean; x represents the genotype; β is the effect of the SNP; and u is the random effects. The kinship matrix was treated as a random-effect factor, whereas the PCA was considered as a fixed-effect factor. Both the kinship and PCA matrix were calculated by the software Tassel v5.0. In this study, markers with a -log₁₀ (p-value) ≥ 3.0 were regarded as MTAs. Manhattan plots and Q–Q plots were drawn by the CMplot package (https://cran.r-project.org/web/packages/CMplot/index.html) based on R 3.6.5.

The flanking sequences corresponding to the SNP markers (the SNPs in the LD decay interval) significantly associated with morphological traits were used in BLASTn and BLASTx searches against NCBI (http://www.ncbi.nlm.nih.gov/) databases. Besides, the annotation for IWGSC v2.1 was also used to identify candidate genes. Potential candidate genes were then selected with significant MTAs corresponding to non-synonymous SNPs in the coding region of the genes.

Quantitative real-time PCR (qRT-PCR) was conducted to test expression differences of the candidate genes in the accessions with extreme traits (Supplementary Table S2). The flag leaf used for test FLL and FLW was sampled for RNA extraction during the flag leaf fully drawn out; the stem used for test PM was sampled for RNA extraction during the erecting period; the uppermost internode used for test UIL was sampled for RNA extraction during the heading stage. Total RNA was extracted according to the Trizol method, whereas cDNA was synthesized with the HiScript II 1st Strand cDNA Synthesis Kit (Vazyme, Nanjing, China). The primers were designed with Primer Premier 5.0 software9 (Supplementary Table S3). PCR procedure was conducted in a volume of 20 μl, containing 2 μl cDNA, 0.4 μl of each primer, and a 10 μl ChamQ Universal SYBR qPCR Master Mix. The reaction was conducted in the ABI StepOnePlus Real-Time PCR System with Tower. The gene expression level was analyzed with 2^–ΔΔCT method. Actin1 was used as an internal control to normalize the expression levels of different samples. All assays were performed in two independent experiments with three repetitions.

Results

Phenotypic evaluation

Significant and continuous variations of adaptive traits were exhibited in the diverse panel. The BLUE values of FLL, FLW, UIL, and PM (0 means prostrate, one means semi-prostrate, two means erect) were 25.0 cm (19.4–29.0 cm), 14.7 mm (10.8–19.1 mm), 22.3 cm (14.1–33.3 cm), and 2.70 (1.3–3.0). The standard deviation and coefficient of variation of FLL, FLW, UIL, and PM were 7.25 cm (0.29), 2.94 mm (0.20), 5.58 cm (0.25), and 0.59 (0.22). The UIL was negatively correlated with FLW (−0.557, p < 0.001). FLL with FLW, PM and UIL, FLW with PM, and UIL with PM show no significant correlation (Supplementary Table S4). ANOVA showed highly significant effects (p < 0.01) of genotypes, environments, and genotype × environment interactions on all traits (Table 1). The h² estimates by a one-stage approach for FLL, FLW, UIL, and PM were 0.78, 0.50, 0.93, and 0.79, whereas the SNP-based heritabilities for FLL, FLW, UIL, and PM were 0.77, 0.50, 0.91, and 0.78, respectively, indicating that morphological traits were determined by genetic factors and affected by environment.

TABLE 1

TABLE 1. ANOVA analysis for the morphological traits in 251 spring wheat accessions from the Northeast of China.

Genotyping, population structure, and linkage disequilibrium decay analysis

Totally, 52,503 polymorphic SNPs after filtering (MAF <0.05, missing rate >0.1) were conducted for further GWAS (Li et al., 2021). Among the SNP markers, 18323, 18691, and 15489 were from A, B, and D genomes, respectively, with an average marker density of 0.273 Mb per marker. All the 251 spring wheat accessions could be divided into three subgroups, subgroups I, II, and III. Of these, subgroup I contained 126 varieties mainly from Heilongjiang ranging from 1950s to 1980s; subgroup II had 75 varieties mainly from Heilongjiang ranging from 1990s to 2010s; whereas subgroup III comprised 50 varieties mainly from Jilin and foreign counties, including United States, Canada, and Japan. Furthermore, the average LD for the whole genome was 8 Mb according to the LOESS curve (Li et al., 2021).

Genome-wide association study

Thirty loci were detected associated with morphological traits (Table 2; Supplementary Table S5; Figure 1; Supplementary Figure S1). Among these, nine loci for FLL were identified on chromosomes 1, 2, 2A, 2B, 3B, 3D, 5D, 6A, and 7D, and each explained 5.7–9.5% of the phenotypic variances, respectively; seven for FLW were identified on chromosomes 1A (2), 5A, 5A, 5B, 5D, and 7A, and explained 5.1–17.9% of the phenotypic variances, respectively; seven for FLL were identified on chromosomes 2B, 5A, 6A, 6D, 7A, 7B, and 7D and explained 4.8–8.8% of the phenotypic variances, respectively; seven for PM were identified on chromosomes 2B, 2B, 3D, 4A, 5A, 5B, and 6A, and each explained 5.3–12.5% of the phenotypic variances, respectively.

TABLE 2

TABLE 2. Loci for morphological traits in 251 spring wheat accessions from the Northeast of China by association analysis.

FIGURE 1

FIGURE 1. Manhattan and Q–Q plot for morphological traits in 251 spring wheat accessions from the Northeast of China analyzed by the mixed linear model (MLM) in Tassel v5.0. FLL, Flag leaf length; FLW, Flag leaf width; UIL, Uppermost internode length; PM, Plant morphology. The 1, 2, 3, 4, and 5 indicated the Harbin 2018, Harbin 2019, Keshan 2018, Keshan 2019, and the best linear unbiased estimation (BLUE), respectively.

Of these, the locus on chromosome 1A (556.0–577.7 Mb) is an identical locus, which showed effects on FLW and FLL; a locus on chromosome 2B (652.7–664.2 Mb) controlled both the FLL and ST; a locus on chromosome 5D (551.2–562.7 Mb) is the same loci for FLL and FLW. Besides, FLL and UIL have the same locus on chromosome 7D (5.4–8.3 Mb).

Candidate genes

Totally, seven candidate genes for morphological traits were identified (Table 3). Two candidate genes encoded for the E3 ubiquitin-protein ligase-like protein (TraesCS1A01G164400 and TraesCS2B01G466600) were identified in the LD decay of the loci on 1A (296.9–297.7 Mb) and 2B (652.7–656.0 Mb) chromosomes. Another gene encoding a cytokinin riboside 5′-monophosphate phosphoribohydrolase (TraesCS1A01G362500) was identified in the LD decay of the loci on chromosome 1A (556.0–587.4 Mb). For the loci on chromosome 5A (383.4–385.0 Mb and 595.4–597.7 Mb), candidate genes TraesCS5A01G182500 and TraesCS5A01G406800 were identified, which encode the F-box family protein. Besides, the genes TraesCS2B01G123800 and TraesCS4A01G055000 encoding the C2H2 zinc finger and leucine-rich repeat receptor protein kinase were identified as the candidate gene for the loci 2B (91.0–92.6 Mb) and 4A (41.4–46.2 Mb). The qRT-PCR of seven candidate genes showed that TraesCS2B01G466600 and TraesCS5A01G406800 showed no significant differences between the extreme accessions; TraesCS1A01G164400, TraesCS4A01G055000, and TraesCS1A01G362500 showed 1.6–3.2-fold higher expressions among the extreme accessions, whereas TraesCS5A01G182500 and TraesCS2B01G123800 showed more than 1.9–7.6-fold lower expressions among extreme accessions (Supplementary Figure S1).

TABLE 3

TABLE 3. The details for the candidate genes of morphological traits.

Discussion

All the 251 spring accessions could be divided into three subgroups (Li et al., 2021), and the characterization of the subgroups was largely consistent with geographic origins, released years, and pedigrees. Most of the cultivars from Heilongjiang ranging from 1950s to 1980s belonged to subgroup-1; the accessions from Heilongjiang ranging from 1990s to 2010s belonged to subgroup-2, whereas subgroup-3 mainly including the accessions from the Jilin and foreign counties (United States, Canada). Population analysis (PCA, NJ-tree, and Kinship) indicated that a significant population structure existed in the 251 spring wheat accessions. The lack of appropriate correction for population structure can lead to spurious MTAs (Zhu et al., 2008). Thus, to eliminate spurious MTAs, an MLM model with subpopulation data (Q) (fixed-effect factors) and kinship matrix (random-effect factor) was conducted. LD decay was influenced by population structure, allele frequency, recombination rate, and selection progresses and significantly affected the precision of association mapping. The LD decay for the whole genome was about 8 Mb, consistent with previous studies (Liu et al., 2017). The number of markers is enough for the subsequent association analysis.

Comparison with the QTL or gene in previous studies

The genes or loci associated with morphological traits were extensively reported previously (Wu et al., 2016; Li et al., 2018; Li et al., 2021). In this study, the association of morphological traits was performed and 30 significant loci were detected.

Several studies for flag leaf-related traits in common wheat have reported. We have identified nine loci for FLL in eight different chromosomes. Li et al. (2020) have identified eight loci for FLL on chromosomes 1A, 2A (2), 2B, 3A, 5A, 6B, and 6D, and each explained 6.9–19.6% of the phenotypic variances. Of these, the loci identified in the 1A and 2A were overlapped with the regions 1A (572.9–577.7 Mb) and 2A (554.6–560.2 Mb) identified in this study. Wu et al. (2016) mapped two FLL QTL on chromosomes 3B and 7D, which is overlapped with the FLL locus (3B: 121.4–129.6 Mb and 7D: 5.4–10.5 Mb) in the present study. Another FLL QTL was previously identified on chromosome 2B linked with the SSR marker barc318 (Liu et al., 2018); it is coincided with the present FLL locus (2B: 652.7–656.0 Mb) based on the consensus linkage map (Maccaferri et al., 2015). These loci on chromosomes 2A (66.3–68.2 Mb), 3D (82.8–94.5 Mb), 5D (551.2–556.2 Mb), and 6A (5.1–17.9 Mb) were potential novel loci. Li et al. (2020) have identified five loci for FLW on chromosomes 1A, 3B, 5B (2), and 6B, accounting for 6.9–11.4% of the phenotypic variances. The locus on chromosome 5B (AX_109519234) was coincided with the locus (5B: 615.7–617.9 Mb) detected in this study. Besides, the locus on chromosome 1A (296.9–297.7 Mb) was also overlapped with the locus on 1A (AX_111540798) identified by Li et al. (2018). Wu et al. (2016) reported a locus for both FLW and flag leaf angle at the position (IWB4576) on chromosome 5A, showing a pleiotropic effect, which coincided with the loci identified in this study (5A: 290.1–299.5 Mb). In conclusion, the loci identified in 1A (556.0–587.4 Mb), 5A (595.4–597.7 Mb), 5D (562.0–562.7 Mb), and 7A (18.8–25.9 Mb) were potentially new.

UIL is important for the construction of plant architecture and influences the yield. Several studies have focused on the genetic basis of UIL. Wu et al. (2016) have identified seven loci for UIL. Li et al. (2020) have identified 12 loci for UIL on chromosomes 1A (2), 1B (2), 3A, 5A, 6B (3), 6D (2), and 7B, with a single locus explaining 6.7–16.4% of the phenotypic variances. Of these, the loci 1A (AX_109449226), 6D (AX_109331000), and 7B (AX_186165710) coincided with 5A (383.4–385.0 Mb), 6D (388.2–407.5 Mb), and 7B (701.1 Mb) identified in this study. Besides, the UIL locus (6D: 388.2–407.5 Mb) is about one LD from the QTL associated with both PH and the third internode length (Cui et al., 2011); they are likely to be the same. Except for the three loci talked above, the remaining four loci (2B: 409.9–439.5 Mb; 6A: 606.6–611.7 Mb; 7A: 496.3–512.9 Mb, and 7D: 8.3 Mb) are all likely to be novel.

Until today, a few studies have been conducted on plant morphology for common wheat. Thus, it is difficult to compare with the present results. Li et al. (2002) centered to the Gli-A2 gliadin locus and associated with a QTL affecting prostrate growth trait on chromosome 6A and nearly with the locus identified in this study in 6A (23.8–24.4 Mb) according to Maccaferri et al., 2015. Liu et al. (2016) have identified a region associated with plant morphology on the 4A chromosome, which maybe coincided with the locus identified in our study (4A: 41.4–46.2 Mb). Although some studies have been focused on the plant habit growth in bread wheat, no marker information was available in order to confirm our results. Thus, we think the loci 2B (91–92.6 Mb, 663.1–664.2 Mb), 3D (602.8–606.7 Mb), 5A (688.9 Mb), and 5B (696.7–697 Mb) were potential be novel.

Among the 30 loci for morphological traits, 13 loci talked above (1A: 572.9–577.7 Mb, 2A: 554.6–560.2 Mb, 2B: 652.7–656.0 Mb, 3B: 121.4–129.6 Mb, 7D: 5.4–10.5 Mb, 1A: 296.9–297.7 Mb, 5A: 290.1–299.5 Mb, 5B: 615.7–617.9 Mb, 5A: 383.4–385.0 Mb, 6D: 388.2–407.5 Mb, 7B: 701.1–701.1 Mb, 4A: 41.4–46.2 Mb, and 6A: 23.8–24.4 Mb) should be the same as the QTL reported in previous studies, whereas the remaining 17 are likely to be new. The stable loci validated by our studies and previous studies indicated that they are widespread and maybe stable in various varieties.

Candidate gene analysis

To identify candidate genes for morphological traits, the flanking sequences of SNP markers (in the LD decay interval and corresponding non-synonymous SNPs in the coding region of the genes) significantly associated with morphological traits were used as queries to BLAST against the NCBI. Totally, seven candidate genes were identified. For loci 1A (296.9–297.7 Mb) and 2B (652.7–656.0 Mb), the candidate genes for the E3 ubiquitin-protein ligase-like protein (TraesCS1A01G164400 and TraesCS2B01G466600) were identified in the LD decay distance. E3 ubiquitin-protein ligase-like protein plays a crucial role in plant growth and development (Karki et al., 2021; Zhang et al., 2022). For loci 4A (41.4–46.2 Mb), candidate gene TraesCS4A01G055000 encoded leucine-rich repeat receptor-like protein kinase family protein, which may trigger multiple physiological pathways (Coleman et al., 2021). For the loci on chromosome 5A (383.4–385.0 Mb and 595.4–597.7 Mb), candidate genes (TraesCS5A01G182500 and TraesCS5A01G406800) for F-box proteins were identified. F-box proteins play crucial roles in cell-cycle progression, transcriptional regulation, flower formation, signal transduction, and many other cellular processes (El-Sharkawy et al., 2021; Guérin et al., 2021). Of these, a cytokinin riboside 5′-monophosphate phosphoribohydrolase (TraesCS1A01G362500) was identified in the LD decay of the loci on chromosome 1A (556.0–587.4 Mb). The cytokinin is a positive regulator of shoot growth (Chen et al., 2021) and response to biotic and abiotic stressors (Ellis et al., 2005; Chen et al., 2022). Besides, the gene TraesCS2B01G123800 was identified as the candidate gene for the loci 2B (91.0–92.6 Mb) and encodes the C2H2 zinc finger, which is important for determining the prostrate/erect plant morphology (Li et al., 2021). The expressions of seven candidate genes showed that TraesCS2B01G466600 and TraesCS5A01G406800 showed no significant differences between the extreme accessions; TraesCS1A01G164400, TraesCS4A01G055000, and TraesCS1A01G362500 showed 1.6–3.2-fold higher expressions among the extreme accessions; whereas TraesCS5A01G182500 and TraesCS2B01G123800 showed more than 1.9- to 7.6-fold lower expressions among extreme accessions. Thus, TraesCS1A01G164400, TraesCS4A01G055000, TraesCS1A01G362500, TraesCS5A01G182500, and TraesCS2B01G123800 are the candidate genes in our study.

Potential implications in wheat breeding

Although conventional breeding has led to improved morphological traits on wheat, selection is time-consuming and not very efficient (He et al., 2010; Rasheed et al., 2016). The SNPs associated with morphological traits identified in this study should facilitate the progresses of MAS and pyramiding favorable alleles will improve morphological traits (Li et al., 2018). Accessions with superior morphological traits and high numbers of favorable alleles (such as Hechun 5, Jichun 158, LongFu 18-171, LongFu 17-5277-12-24 and Longchun 1 with longer FLL; Xinshuguang 4, Dongnong 156597, Kechun 151350, Kenda 163672 and Gang 09-558 with wider FLW; Longmai 17, Kehua, Xiaomaixiaobing 4, Xiaomaixiaobing seven and Yongjie with longer UIL; Kehan 1, Shen 68-71, Mailiduo, Jichun101 and Jichun1 with erect, whereas Beimai6, Lianfeng, Kechun 140865, Beimai one and LongFumai8 with prostrate) could be used as parental lines for the improvement of morphological traits in wheat breeding. The loci with pleiotropic and consistent effects across each environment in this study should be amenable to MAS.

Conclusion

We have identified 30 loci for morphological traits in spring wheat accessions by GWAS. Of these, 17 loci were likely to be new. Besides, five candidate genes were identified for morphological traits. The associated markers and varieties with favorable alleles could be used to accelerate the progresses of wheat MAS breeding.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.

Author contributions

WL and YL designed the research. YL and YS analyzed the physiology data. WL and YL drafted the article. WL, JT, JC, SY, XW, and RZ performed the experiment. HZ and YL revised the article. All authors have read, edited, and approved the current version of the article.

Funding

This work was funded by the “Crop Varietal Improvement and Insect Pests Control by Nuclear Radiation.”

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.

Supplementary material

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

SUPPLEMENTARY FIGURE S1 | The qRT-PCR results of the seven candidate genes for the extreme accessions of the 25 spring wheat accessions from the Northeast of China. The extreme accessions of the 25 spring wheat accessions from the Northeast of China were according to Supplementary Table S2.

Abbreviations

BLUE, Best linear unbiased estimation; FLL, Flag leaf length; FLW, Flag leaf width; PM, Plant morphology; GWAS, Genome-wide association study; hb2, Broad-sense heritability; KASP, Kompetitive allele-specific PCR; LD, Linkage disequilibrium; MAS, Marker-assisted selection; QTL, Quantitative trait loci; R2, Phenotypic variance explained; RIL, Recombinant inbred line; SNP, Single nucleotide polymorphism; UIL, Uppermost internode length.

References

Beyer, S., Daba, S., Tyagi, P., Bockelman, H., Brown-Guedira, G., and Mohammadi, M. (2019). Loci and candidate genes controlling root traits in wheat seedlings-a wheat root GWAS. Funct. Integr. Genomics 19, 91–107. doi:10.1007/s10142-018-0630-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Breseghello, F., and Sorrells, M. E. (2006). Association mapping of kernel size and milling quality in wheat (Triticum aestivum L.) cultivars. Genetics 172, 1165–1177. doi:10.1534/genetics.105.044586

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, L., Jameson, G. B., Guo, Y., Song, J., and Jameson, P. E. (2022). The LONELY GUY gene family: From mosses to wheat, the key to the formation of active cytokinins in plants. Plant Biotechnol. J. 20, 625–645. doi:10.1111/pbi.13783

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, L., Zhao, J., Song, J., and Jameson, P. E. (2021). Cytokinin glucosyl transferases, key regulators of cytokinin homeostasis, have potential value for wheat improvement. Plant Biotechnol. J. 19, 878–896. doi:10.1111/pbi.13595

PubMed Abstract | CrossRef Full Text | Google Scholar

Coleman, A. D., Maroschek, J., Raasch, L., Takken, F. L., Ranf, S., Hückelhoven, R., et al. (2021). The Arabidopsis leucine‐rich repeat receptor‐like kinase MIK2 is a crucial component of early immune responses to a fungal‐derived elicitor. New Phytol. 229, 3453–3466. doi:10.1111/nph.17122

PubMed Abstract | CrossRef Full Text | Google Scholar

Cui, F., Li, J., Ding, A., Zhao, C., Wang, L., Wang, X., et al. (2011). Conditional QTL mapping for plant height with respect to the length of the spike and internode in two mapping populations of wheat. Theor. Appl. Genet. 122, 1517–1536. doi:10.1007/s00122-011-1551-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Cui, F., Zhao, C., Ding, A., Li, J., Wang, L., Li, X., et al. (2014). Construction of an integrative linkage map and QTL mapping of grain yield-related traits using three related wheat RIL populations. Theor. Appl. Genet. 127, 659–675. doi:10.1007/s00122-013-2249-8

PubMed Abstract | CrossRef Full Text | Google Scholar

El-Sharkawy, I., Ismail, A., Darwish, A., El Kayal, W., Subramanian, J., Sherif, S. M., et al. (2021). Functional characterization of a gibberellin F-box protein, PslSLY1, during plum fruit development. J. Exp. Bot. 72, 371–384. doi:10.1093/jxb/eraa438

PubMed Abstract | CrossRef Full Text | Google Scholar

Ellis, M. H., Rebetzke, G. J., Azanza, F., Richards, R. A., and Spielmeyer, W. (2005). Molecular mapping of gibberellin-responsive dwarfing genes in bread wheat. Theor. Appl. Genet. 111, 423–430. doi:10.1007/s00122-005-2008-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, F., Ma, D., Yin, G., Rasheed, A., Dong, Y., Xiao, Y., et al. (2017). Genetic progress in grain yield and physiological traits in Chinese wheat cultivars of Southern Yellow and Huai Valley since 1950. Crop Sci. 57, 760–773. doi:10.2135/cropsci2016.05.0362

CrossRef Full Text | Google Scholar

Gao, F., Wen, W., Liu, J., Rasheed, A., Yin, G., Xia, X., et al. (2015). Genome-wide linkage mapping of QTL for yield components, plant height and yield-related physiological traits in the Chinese wheat cross Zhou 8425B/Chinese Spring. Front. Plant Sci. 6, 1099. doi:10.3389/fpls.2015.01099

PubMed Abstract | CrossRef Full Text | Google Scholar

Guérin, C., Mouzeyar, S., and Roche, J. (2021). The landscape of the genomic distribution and the expression of the F-Box genes unveil genome plasticity in hexaploid wheat during grain development and in response to heat and drought stress. Int. J. Mol. Sci. 22, 3111. doi:10.3390/ijms22063111

PubMed Abstract | CrossRef Full Text | Google Scholar

He, Z. H., Xia, X. C., and Bonjean, A. P. A. (2010). Wheat improvement in China. Mexico: CIMMYT, 51–68.

Google Scholar

Karki, S. J., Reilly, A., Zhou, B., Mascarello, M., Burke, J., Doohan, F., et al. (2021). A small secreted protein from Zymoseptoria tritici interacts with a wheat E3 ubiquitin ligase to promote disease. J. Exp. Bot. 72, 733–746. doi:10.1093/jxb/eraa489

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, F., Wen, W., He, Z., Liu, J., Jin, H., Cao, S., et al. (2018). Genome-wide linkage mapping of yield related traits in three Chinese bread wheat populations using high-density SNP markers. Theor. Appl. Genet. 131, 1903–1924. doi:10.1007/s00122-018-3122-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, F., Wen, W., Liu, J., Zhang, Y., Cao, S., He, Z., et al. (2019). Genetic architecture of grain yield in bread wheat based on genome-wide association studies. BMC Plant Biol. 19, 168. doi:10.1186/s12870-019-1781-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, W. L., Nelson, J. C., Chu, C. Y., Shi, L. H., Huang, S. H., Liu, D. J., et al. (2002). Chromosomal locations and genetic relationships of tiller and spike characters in wheat. Euphytica 125, 357–366. doi:10.1023/a:1016069809977

CrossRef Full Text | Google Scholar

Li, Y., Sun, A., Wu, Q., Zou, X., Chen, F., Cai, R., et al. (2021). Comprehensive genomic survey, structural classification and expression analysis of C2H2-type zinc finger factor in wheat (Triticum aestivum L.) BMC Plant Biol. 21, 380. doi:10.1186/s12870-021-03016-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J., He, Z., Wu, L., Bai, B., Wen, W., Xie, C., et al. (2016a). Genome-wide linkage mapping of QTL for black point reaction in bread wheat (Triticum aestivum L.) Theor. Appl. Genet. 129, 2179–2190. doi:10.1007/s00122-016-2766-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, X., Huang, M., Fan, B., Buckler, E. S., and Zhang, Z. (2016b). Iterative usage of fixed and random effect models for powerful and efficient genome-wide association studies. PLoS Genet. 12, e1005767. doi:10.1371/journal.pgen.1005767

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J., He, Z., Rasheed, A., Wen, W., Yan, J., Zhang, P., et al. (2017). Genome-Wide Association Mapping of Black Point Reaction in Common Wheat (Triticum aestivum L.. BMC Plant Bio. 17, 1–12.

CrossRef Full Text | Google Scholar

Liu, J., Rasheed, A., He, Z., Imtiaz, M., Arif, A., Mahmood, T., et al. (2019). Genome-Wide Variation Patterns Between Landraces and Cultivars Uncover Divergent Selection During Modern Wheat Breeding. Theor. Appl. Genet. 132, 2509–2523.

PubMed Abstract | CrossRef Full Text | Google Scholar

Lozada, D. N., Mason, R. E., Sarinelli, J. M., and Brown-Guedira, G. (2019). Accuracy of genomic selection for grain yield and agronomic traits in soft red winter wheat. BMC Genet. 20, 82. doi:10.1186/s12863-019-0785-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Nadolska-Orczyk, A., Rajchel, I. K., Orczyk, W., and Gasparis, S. (2017). Major Genes Determining Yield-Related Traits in Wheat and Barley. Theor. Appl. Genet. 130, 1081–1098.

PubMed Abstract | CrossRef Full Text | Google Scholar

Pang, Y., Liu, C., Wang, D., Amand, P. S., Bernardo, A., Li, W., et al. (2020). High-resolution genome-wide association study identifies genomic regions and candidate genes for important agronomic traits in wheat. Mol. Plant 13, 1311–1327. doi:10.1016/j.molp.2020.07.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Pritchard, J. K., Stephens, M., Rosenberg, N. A., and Donnelly, P. (2000). Association mapping in structured populations. Am. J. Hum. Genet. 67, 170–181. doi:10.1086/302959

PubMed Abstract | CrossRef Full Text | Google Scholar

Quan, X., Dong, L. J., Zhang, N., Xie, C., Li, H., Xia, X., et al. (2021). Genome-wide association study uncover the genetic architecture of salt tolerance-related traits in common wheat (Triticum aestivum L.) Front. Genet. 12, 663941. doi:10.3389/fgene.2021.663941

PubMed Abstract | CrossRef Full Text | Google Scholar

Rahimi, Y., Bihamta, M. R., Taleei, A., Alipour, H., and Ingvarsson, P. K. (2019). Genome-wide association study of agronomic traits in bread wheat reveals novel putative alleles for future breeding programs. BMC Plant Biol. 19 (1), 541. doi:10.1186/s12870-019-2165-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Rasheed, A., Wen, W., Gao, F., Zhai, S., Jin, H., Liu, J., et al. (2016). Development and validation of KASP assays for genes underpinning key economic traits in bread wheat. Theor. Appl. Genet. 129, 1843–1860. doi:10.1007/s00122-016-2743-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, C. W., Zhang, F. Y., Yan, X. F., Zhang, X. F., Dong, Z. D., Cui, D. Q., et al. (2017). Genome-wide association study for 13 agronomic traits reveals distribution of superior alleles in bread wheat from the Yellow and Huai Valley of China. Plant Biotechnol. J. 15, 953–969. doi:10.1111/pbi.12690

PubMed Abstract | CrossRef Full Text | Google Scholar

Valluru, R., Reynolds, M. P., Davies, W. J., and Sukumaran, S. (2017). Phenotypic and genome-wide association analysis of spike ethylene in diverse wheat genotypes under heat stress. New Phytol. 214, 271–283. doi:10.1111/nph.14367

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Wang, R., Mao, X., Zhang, J., Liu, Y., Xie, Q., et al. (2020). RING finger ubiquitin E3 ligase gene TaSDIR1-4A contributes to determination of grain size in common wheat. J. Exp. Bot. 71, 5377–5388. doi:10.1093/jxb/eraa271

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S. C., Wong, D., Forrest, K., Allen, A., Chao, S., Huang, B. E., et al. (2014). Characterization of polyploid wheat genomic diversity using a high-density 90000 single nucleotide polymorphism array. Plant Biotechnol. J. 12, 787–796. doi:10.1111/pbi.12183

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Hou, J., Liu, H., Li, T., Wang, K., Hao, C., et al. (2019). TaBT1, affecting starch synthesis and thousand kernel weight, underwent strong selection during wheat improvement. J. Exp. Bot. 70, 1497–1511. doi:10.1093/jxb/erz032

PubMed Abstract | CrossRef Full Text | Google Scholar

Würschum, T., Leiser, W. L., Langer, S. M., Tucker, M. R., and Longin, C. F. H. (2018). Phenotypic and Genetic Analysis of Spike and Kernel Characteristics in Wheat Reveals Long-Term Genetic Trends of Grain Yield Components. Theor. Appl. Genet. 131, 2071–2084.

PubMed Abstract | Google Scholar

Xu, X. P., Ping, H., Chuan, L. M., Liu, X. Y., Liu, Y. X., Zhang, J. J., et al. (2021). Regional distribution of wheat yield and chemical fertilizer requirements in China. J. Integr. Agric. 20, 2772–2780. doi:10.1016/s2095-3119(20)63338-x

CrossRef Full Text | Google Scholar

Yang, J. A., Lee, S. H., Goddard, M. E., and Visscher, P. M. (2011). Gcta: A tool for genome-wide complex trait analysis. Am. J. Hum. Genet. 88, 76–82. doi:10.1016/j.ajhg.2010.11.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, D., Zhang, X., Xu, W., Hu, T., Ma, J., Zhang, Y., et al. (2022). TaGW2L, a GW2-like RING finger E3 ligase, positively regulates heading date in common wheat (Triticum aestivum L.) Crop J. doi.org/doi:10.1016/j.cj.2021.12.002

CrossRef Full Text | Google Scholar

Zhu, C., Gore, M., Buckler, E. S., and Yu, J. (2008). Status and prospects of association mapping in plants. Plant Genome 1, plantgenome2008.02.0089. doi:10.3835/plantgenome2008.02.0089

CrossRef Full Text | Google Scholar