Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Plant Sci., 24 July 2025

Sec. Plant Breeding

Volume 16 - 2025 | https://doi.org/10.3389/fpls.2025.1635721

Genome-wide association analysis of wheat stem traits using 55K microarrays

Wei Wang&#x;Wei Wang1†Na Sun&#x;Na Sun2†Kai ZhaoKai Zhao1JiKun SongJiKun Song3Hui FangHui Fang4Guiqiang FanGuiqiang Fan4Yonghong GaoYonghong Gao4Tianrong HuangTianrong Huang4Yindeng Ding*Yindeng Ding4*
  • 1School of Computer Science and Information Engineering, Anyang Institute of Technology, Anyang, China
  • 2Yili Prefecture Institute of Agricultural Science, Yining, China
  • 3Cotton Research Institute, Chinese Academy of Agricultural Sciences, AnYang, China
  • 4Institute of Crops, Xinjiang Academy of Agricultural Sciences, Xinjiang, China

Wheat lodging poses a severe threat to yield and quality, with the morphological and structural traits of the basal second internode being critical determinants of stem strength and lodging resistance. We conducted phenotypic analysis on 239 wheat varieties (lines) collected from around the world. This analysis was complemented by genotyping using the wheat 55K SNP chip. Genome-wide association analysis (GWAS) was executed employing the MLM (Q+K) algorithm within the TASSLE software suite. The findings unveiled pronounced phenotypic variability in the stem diameter of the second internode across disparate temporal intervals, characterized by a coefficient of variation spanning from 11.31% to 13.95%, alongside robust inter-year correlations. Furthermore, the genome-wide linkage disequilibrium (LD) decay distance was ascertained to 3 Mb. Analyses of population structure, phylogenetic dendrograms, and principal components revealed that the intrinsic population structure of the 239 wheat varieties (lines) was markedly simplified, segregating into three discernible subgroups. GWAS analysis identified 118 SNPs significantly associated with the stem diameter of the second internode (P ≤ 0.001). Notably, among these loci, three SNPs (AX-111557672, AX-94584919, and AX-109819835) overlapped with previously reported associations, while the remaining 115 SNPs represented novel discoveries distributed across chromosomes 1A, 1B, 1D, 2A, 2B, 2D, 3A, 3B, 4A, 4B, 4D, 5A, 5B, 6A, 6B, 6D, 7A, 7B, and 7D. These newly identified loci exhibited substantial phenotypic variance explained, ranging from 8.09% to 29.14% for individual SNPs. Subsequent rigorous screening of loci showing significant phenotypic impacts and stability across diverse environmental contexts culminated in the identification of seven candidate genes implicated in the stem diameter of the second internode. This investigation provides new SNPs markers for enhancing lodging resistance in wheat, highlighting substantial practical implications.

Introduction

In the context of escalating global food demands and climate change, ensuring the sustainable production of staple crops is critical for global food security. Wheat (Triticum aestivum L.), as one of the most widely cultivated cereal crops, accounts for approximately 21% of human dietary calories and 20% of protein intake (Wang et al., 2022a, b, 2024). Its unique gluten-forming proteins endow wheat-based products with distinct textures, making it indispensable in global diets (Liu et al., 2017, 2023). However, the pursuit of high-yield wheat breeding has led to a significant trade-off between yield potential and resistance to abiotic stresses. Lodging, defined as the displacement of the upright plant structure from its normal position, has emerged as a major constraint to wheat productivity. When wheat plants encounter adverse environmental conditions, such as strong winds during the grain-filling stage or waterlogging-induced soil instability, lodging can cause a yield reduction ranging from 10% to 40% (Zhu et al., 2008; Sun et al., 2017). Beyond yield losses, lodging compromises grain quality by reducing kernel plumpness, increases disease susceptibility, and complicates mechanized harvesting (Tian et al., 2017; Kong et al., 2024). Among different types of lodging, stem lodging, primarily characterized by the bending or breaking of basal internodes, is particularly prevalent during the milk-ripening stage (Anderson et al., 2025). Despite its economic significance, the genetic architecture underlying stem-related traits, such as internode strength, height, and mechanical properties, remains incompletely understood. Therefore, the study of genetic mechanisms of wheat stalk traits and improvement of stalk characteristics can help to improve the resistance of wheat to failure.

The second internode from the base of the stem is a critical support point for the wheat plant. Its morphological and structural characteristics have a significant impact on the plant’s lodging resistance (Gong, 1982; Kelbert et al., 2004; Zhang et al., 2016; Crook and Ennos, 1994). Research findings indicate that stem diameter, stem wall thickness, and mechanical strength at the second internode of the basal stem exhibit a significant positive correlation with lodging resistance (Rasheed et al., 2014). Plants with greater stem diameter and thicker stem walls generally exhibit enhanced lodging resistance. These features enhance the stem’s mechanical strength, allowing the plant to better endure external forces and resist bending or breaking under unfavorable conditions (Cat, 2024; Chen et al., 2024). Moreover, the composition of the cell wall at the second internode of the basal stem, particularly the content of cellulose and lignin, significantly influences lodging resistance. These structural components play a pivotal role in enhancing the mechanical properties of the stem, thereby improving its ability to withstand external stresses such as wind or heavy rainfall (Zhu et al., 2008; Tolossa et al., 2024; Xu et al., 2024). The higher the cellulose content of wheat stalks, the stronger the compressive strength of the stalks, while varieties with lower cellulose content are prone to collapse because their stalks are more fragile and weaker (Mostafaie et al., 2024). In recent years, the rapid development of molecular biology techniques has provided new perspectives for studying the genetic mechanisms of wheat resistance traits to failure (Niu et al., 2022). Researchers localized genetic loci associated with stalk traits through different populations and marker techniques (Khan et al., 2024). A total of 23 QTLs were identified on chromosomes 1B, 3B, 4A, 4B, and 5A, associated with traits such as stem strength, diameter, and pith cavity diameter. The phenotypic contribution of each QTL ranged from 3.5% to 44.0% (Dhariwal et al., 2022). Further research identified two QTLs associated with stem strength on chromosomes 3B and 5A. (Liu et al., 2017) using the 90K SNP chip, detected nine QTLs on chromosome 4B that regulate stem strength and the wall thickness at the second internode from the base, with phenotypic variation explained ranging from 9.4% to 36.6%. These findings elucidate the polygenic regulation of stem-related traits and identify chromosomal hotspot regions. In natural wheat populations, 37 SNP loci significantly associated with stem strength have been identified, distributed across 15 chromosomes (1A, 1B, 2B, 2D, 3A, 3B, 4A, 4B, 5A, 5B, 5D, 6B, 7A, 7B, and 7D). The phenotypic contribution rate of individual markers ranged from 7.76% to 13.77%, highlighting the complex genetic architecture underlying stem mechanical properties (Esposito et al., 2022). Analysis of 105 materials by 90K SNP microarrays identified multiple stem strength-associated loci on chromosomes 2B, 3B, 4A, 6B, and 7B, further confirming the synergistic multi-chromosome regulation (Genievskaya et al., 2020). The genes involved in the lignin biosynthesis pathway play a central role in the regulation of stem strength. Comparative transcriptomic analysis between the dwarf mutant DC20 and the wild-type (WT) wheat revealed 2135 differentially expressed genes in the stems during the booting stage. These genes primarily participate in sugar metabolism, energy metabolism, and post-transcriptional modification processes. Notably, 47 genes were significantly implicated in the dynamic balance of gibberellin (GAS) and auxin (IAA), cell elongation, and signal transduction (Bainsla et al., 2020). The study revealed that the expression levels of most differential genes were downregulated, while the expression of inhibitory factors was upregulated, suggesting that the dwarfism in the DC20 mutant may be closely associated with the suppression of hormone signaling and dysregulation of cell cycle control (Saint Pierre et al., 2010). Through these studies, scientists have continued to explore in depth the genetic mechanism of wheat resistance traits, providing a theoretical basis for further variety improvement and molecular breeding.

Although numerous studies have explored the genetic mechanisms underlying wheat lodging resistance, the genetic basis of morphological and structural characteristics in the second internode at the stem base remains systematically understudied. This study used 239 wheat varieties (lines) from around the world as research material to conduct a comprehensive analysis of morphological and structural traits, including stem diameter, at the second internode from the base. GWAS was conducted using the 55K SNP genotyping chip data for wheat, which covers 16,649 SNP markers across the entire genome. The analysis was performed using the mixed linear model (MLM, Q+K model) in the TASSEL software. Through GWAS, we aimed to identify significant SNP loci associated with stem diameter at the second internode from the base and further explore potential candidate genes. These findings will help to reveal the genetic basis of stem mechanical properties and provide targets for molecular breeding of lodging-resistant wheat varieties, and offer genetic resources for the development of lodging-resistant wheat varieties, thereby accelerating the breeding process for such varieties.

Materials and methods

Materials

This study selected 239 wheat varieties as experimental materials, covering multiple regions both domestically and internationally. Among them, 89.1% (213 varieties) were from China, primarily from Xinjiang (115 varieties), Beijing (34 varieties), Hebei (19 varieties), and Henan (19 varieties), with the remaining varieties distributed across Anhui, Jiangsu, Shandong, Shanxi, Shaanxi, and Tianjin provinces. The remaining 10.90% (26 varieties) were from abroad, with 25 varieties from the United States and 1 variety from Ukraine. All materials used in the study were able to adapt to the environmental conditions of the experimental site and completed their normal growth cycle.

Experimental design

The study was conducted from September 2019 to July 2022 at two experimental sites of the Xinjiang Academy of Agricultural Sciences, namely the Zepu site (77°16’17.22”N, 38°11’21.65”E) and the Anningqu site (43°58’53.38”N, 87°30’17.72”E). Specifically, the Zepu site was designated as E1 in 2020, E2 in 2021, and E3 in 2022, while the Anningqu site was designated as E4 in 2020, E5 in 2021, and E6 in 2022. The six environmental conditions were tested over three years using a randomized block design, with two rows of each material planted per plot. The row length was 1 meter, the row spacing was 20 cm, and the sowing density was 350,000 seeds per acre. The planting direction was from north to south. Field management included control of light, water, fertilizer, and other routine practices. Each environmental condition was replicated three times to ensure consistency in cultivation and production conditions.

Measurement of stem diameter at the second internode from the base

Following the milky stage of wheat, the first internode counted from the ground upward was designated as the basal first internode, and the second internode as the basal second internode. Using vernier calipers, five plants with comparable growth vigor from each experimental plot were measured, and their average value was calculated and recorded as the stem diameter of the basal second internode for the corresponding plot.

Phenotyping methods

During the phenotypic analysis, variance analysis (ANOVA), descriptive statistics, significance tests for differences, and correlation analysis were conducted. All statistical procedures were performed using Python 3.8 in the Spyder IDE under Anaconda3, running on a Windows 10 workstation equipped with an Intel i7–6800 K 3.40 GHz CPU, 16 GB RAM, and a NVIDIA GeForce RTX 2080Ti GPU. The stem diameter traits of the basal second internode in wheat were analyzed using Pandas (v1.3.2), Matplotlib (v3.4.2), Scikit-learn (v0.24.2), and SPSS (v21.0) for statistical computation, correlation evaluation, and significance assessment.

Association analysis

Association analysis, also termed association mapping or linkage disequilibrium (LD) mapping, is a quantitative genetic analytical technique. It integrates phenotypic diversity of target traits with polymorphisms in genes or markers, and identifies loci or markers associated with target traits based on LD between alleles at different loci in natural populations.

In this study, association analyses were conducted on wheat populations under diverse growing conditions using TASSEL v5.0 software, with integration of phenotypic evaluation datasets and a 55K SNP array. The Mixed Linear Model (MLM) was implemented for populations association profiling. Computational analysis revealed SNP markers exhibiting a significant association with phenotypic traits at a threshold of P ≤ 0.001. Subsequent visualization generated Manhattan plots and QQ plots. QQ plots validated analytical robustness by filtering spurious associations, while Manhattan plots mapped SNP distributions across wheat’s 21 chromosomes (horizontal axis) against −log10(P-value) (vertical axis), pinpointing loci with trait correlations.

Candidate gene

The extended sequences of stable SNP markers were subjected to BLAST alignment in the Chinese Spring wheat genome database (https://urgi.versailles.inrae.fr/blast_iwgsc/). This was done to identify their locations and potential homologous sequences in the wheat genome. After obtaining the alignment results, gene function annotations were performed using the Wheat Omics 1.0 database (http://202.194.139.32/tools/intervalTools.html), which provided detailed information about the genes related to the SNP markers, including their functional roles, associated pathways, and other relevant biological information. This process helps to explore the genetic and functional significance of the identified SNP markers in wheat.

Results

Analysis of phenotypic variation in stem diameter at the second basal internode of wheat

Stem thickness of the 2nd internode at the base of wheat stalks was phenotyped at six environmental points, E1, E2, E3, E4, E5, and E6. The data were evaluated using four metrics: mean (μ), median, coefficient of variation (cv), and standard deviation (σ). As seen in Figure 1, the mean μ of stem thickness at the basal internode 2 of wheat stems in the E1 environment was 2.18, the median was 2.19, the coefficient of variation was 13.9%, and the standard deviation was 0.30; the mean μ of stem thickness at basal internode 2 of wheat stems in E2 environment was 2.29, median was 2.30, coefficient of variation was 13.7% and standard deviation was 0.31; the mean μ in the E3 environment was 2.25, the median was 2.27, the coefficient of variation was 11.3% and the standard deviation was 0.25; In the E4 environment, the mean (μ) was 2.26, with the median at 2.25, a coefficient of variation of 12.0%, and a standard deviation of 0.27. As for the E5 environment, it showed a mean (μ) of 2.31, a median of 2.33, a coefficient of variation of 11.2%, and a standard deviation of 0.26. The E6 environment, on the other hand, had a mean (μ) of 2.24, a median of 2.26, a coefficient of variation of 11.6%, and the standard deviation was 0.26. Overall, the data across different environments exhibit continuity and a normal distribution, which aligns with the typical characteristics of quantitative traits.

Figure 1
Bar chart comparing values for types E1 to E6. E6 has the highest value, while E1 has the lowest. Error bars suggest variability. Significance is indicated with asterisks, showing differing levels of statistical significance among the groups.

Figure 1. Trait distribution of the wheat stem diameter at the second internode from the base under different environmental conditions. E1: 2020 Zepu, Xinjiang; E2: 2021 Zepu, Xinjiang; E3: 2022 Zepu, Xinjiang; E4: 2020 Anningqu, Xinjiang; E5: 2021 Anningqu, Xinjiang; E6: 2022 Anningqu, Xinjiang. * p≤0.05; ** p≤0.01; *** p≤0.001; **** p≤0.0001.

In the present investigation, we undertake a rigorous analysis that amalgamates significance testing, as delineated in Figure 2, with variance decomposition within a Best Linear Unbiased Prediction (BLUP) paradigm, as tabulated in Table 1. Our analytical framework has unveiled nuanced variability patterns in the second internode stem diameter of wheat culms. Empirical data demonstrate that this trait exhibits a standard deviation ranging from 0.16 to 0.31, with coefficients of variation fluctuating between 6.98% and 13.95%, and a heritability estimate of 0.73. A thorough variance analysis elucidates that genotypic disparities are the predominant contributors to the observed variation in stem diameter at the second internode, with environmental factors exerting a secondary influence. Collectively, these findings emphatically underscore the central role of genetic effects in orchestrating the phenotypic differentiation of this critical trait.

Figure 2
Six histograms titled E1 to E6 display frequency distributions with orange dashed normal distribution curves overlaying blue bars. Each plot includes statistics: mean (μ), median, coefficient of variation (cv), and standard deviation (σ). Values for μ range from 2.18 to 2.31, median from 2.19 to 2.33, cv from 0.112 to 0.139, and σ from 0.25 to 0.31. Horizontal axes range from 1.5 to 3.0 or greater, and vertical axes range from 0 to 2.

Figure 2. Statistical significance chart for the wheat stem diameter at the second internode from the base under different environmental conditions.

Table 1
www.frontiersin.org

Table 1. Descriptive statistics for stem diameter trait at the second internode of wheat culm bases.

Correlation analysis of stem diameter traits at the second internode of wheat culm bases

In this study, the results of correlation analysis of stem thickness traits at the basal 2nd internode of wheat stalks at six experimental sites (E1-E6) are presented in Figure 3. The data showed that the correlation of this trait reached highly significant level (p<0.001) among all environmental points. The correlation coefficients for the stem thickness trait at the basal internode 2 of wheat stalks ranged from 0.50 to 0.57 in E1, 0.50 to 0.57 in E2, 0.43 to 0.54 in E3, 0.43 to 0.56 in E4, 0.52 to 0.60 in E5 and 0.49 to 0.60 in E6. Overall, the correlation coefficients of all environmental traits were distributed in the range of 0.43-0.60, indicating that this agronomic trait had stable and significant correlations under different environmental conditions.

Figure 3
Matrix scatter plot showing relationships among five variables (E1 to E5). Diagonal contains histograms with red density lines. Off-diagonal plots display scatter plots with red trend lines. Correlation coefficients range from 0.43 to 0.60. Significant correlations are marked with triple asterisks.

Figure 3. Correlation of stem thickness traits at the basal 2nd internode of wheat stalks in different environments.

Population structure and evolutionary tree analysis

A comprehensive analysis of the population structure of 239 wheat varieties was performed using Structure software. Through a careful integration of population structure analysis, phylogenetic tree reconstruction and principal component analysis (PCA), we systematically partitioned the population structure into distinct subgroups. Our results, which were consistently supported by all three analytical approaches, indicate that it is both scientifically sound and methodologically appropriate to divide the total population into three subgroups: Subgroup 1, comprising 95 varieties (or lines); Subgroup 2, comprising 89 varieties (or lines); and Subgroup 3, comprising 55 varieties (or lines). The distribution frequencies of these subgroups are 39.75%, 37.24% and 23.01% respectively. In addition, we performed rigorous calculations to determine the linkage disequilibrium (LD) decay distances for the 239 wheat varieties (or lines) across genomes A, B, D and the entire genome, yielding values of 3 Mb, 3 Mb, 2 Mb and 3 Mb, respectively. Based on the LD decay distance observed for the whole genome, we have established a criterion whereby loci located within a 3 Mb interval on the physical map are designated as candidate loci for further analysis and investigation (Ding et al., 2023).

GWAS analysis of stem thickness traits in the basal 2nd internode of wheat stalks

A genome-wide association study (GWAS) was conducted using TASSEL 5.0 software, combining the stem thickness at the second internode from the base of 239 wheat varieties (lines) with 16,649 high-quality SNP markers identified from the 55K SNP array. Based on the MLM (Q+K) model, markers exhibiting a significance threshold of P ≤ 0.001 were deemed significantly associated with the traits. Loci consistently detected across multiple environments were identified as stably inherited genetic regions (Figure 4; Table 2). GWAS Analysis of Wheat Stem Diameter at the Second Basal Internode. Based on genome-wide association study (GWAS) results, 139 single nucleotide polymorphism (SNP) markers were significantly associated with stem diameter at the second basal internode in wheat. These markers were distributed across 17 chromosomes: 1A, 1B, 1D, 2A, 2B, 2D, 3A, 3B, 4A, 4B, 4D, 5A, 5B, 6A, 6B, 6D, 7A, 7B, and 7D (Figure 1). Individual SNPs explained phenotypic variation rates ranging from 8.09% to 29.14%, highlighting their substantial contribution to trait heritability. AX-111472018 and AX-111589920, located on chromosome 1B, are separated by only 0.2 Mb, with chromosomal clusters identified, and the phenotypic variation rate ranging from 10.42% to 10.75%. AX-109449735, localized on chromosome 2D, exhibits chromosome clustering and is positioned closer to other chromosomes within a region of reduced linkage disequilibrium (LD). The phenotypic variation rate for this marker ranges from 13.47% to 22.26%. AX-11157475 and AX-111490213, located on chromosome 4D, are separated by 0.3 Mb, show chromosome clustering, and the phenotypic variation rate ranges from 11.85% to 17.81%. AX-111502310 and AX-110491224, mapped to chromosome 5A, are separated by 0.8 Mb, with chromosomal clustering observed, and the phenotypic variation rates range from 15.90% to 17.57%. The remaining markers were detected in only one environment. Among the 118 identified SNPs, 115 novel loci were discovered including AX-111089714 and AX-109453515, which have not been previously associated with stem thickness traits in the basal 2nd internode of wheat stalks.

Figure 4
The image displays several Manhattan plots and Q-Q plots across different environments and a BLUP analysis. Each Manhattan plot shows colorful data points representing -log10 p-values across chromosomes labeled 1B to 7D. The Q-Q plots compare observed and expected -log10 p-values, showing deviation from expected under the null hypothesis. Each section is labeled E1 to E6 and Blup, indicating different experimental conditions or analyses.

Figure 4. Manhattan plot of stem thickness traits in the basal internode 2 of wheat stalks in different environments.

Table 2
www.frontiersin.org

Table 2. Information on significantly associated loci for stem thickness traits in the basal 2nd internode of wheat stalks.

Functional prediction of candidate genes for stem thickness traits in the basal 2nd internode of wheat stalks

Significant SNP markers with substantial phenotypic effects and stable inheritance were retrieved from the Chinese Spring wheat genome database. Additionally, BLASTx sequence alignment was conducted against the NCBI database, leading to the identification of seven candidate genes most likely associated with the stem diameter at the second internode of the basal stem in wheat (Table 3). These candidate genes are primarily involved in the fibrillin family, calcium ion binding, disease resistance protein RPM1, and proteins containing the F-box domain. The gene TraesCS1A02G422100, located on chromosome 1A, encodes Asparagine synthase; TraesCS1A01G058600, also on chromosome 1A, is associated with Peroxisomal membrane 22 kDa (Mpv17/PMP22) family protein; TraesCS2D01G488200 on chromosome 2D encodes Plastid-lipid associated protein PAP/fibrillin family protein; TraesCS4D01G281500 on chromosome 4D encodes an F-box family protein; TraesCS5A01G073300 on chromosome 5A encodes Disease resistance protein RPM1; TraesCS6B02G050900 on chromosome 6B encodes a protein kinase catalytic domain; and TraesCS7A02G481800 on chromosome 7A encodes a SNARE motif protein.

Table 3
www.frontiersin.org

Table 3. Information of candidate genes for stem thickness trait in the 2nd internode at the base of wheat stalks.

Discussion

Association analysis of stem thickness traits in the basal second internode of wheat stalks

With the rapid development of biology and bioinformatics, GWAS (genome-wide association analysis) has become an important way to study quantitative traits in plants, and the mining of genes related to stem thickness traits in the 2nd internode at the base of wheat stalks has also been promoted to a greater extent (Li et al., 2024). Classical genetic studies have shown that stem diameter at the second internode of the basal stem in wheat is a complex trait controlled by both Mendelian and quantitative genes. This trait is closely associated with wheat lodging resistance, as the second internode at the base of the plant serves as a critical mechanical support, and its stem diameter directly affects the mechanical strength and bending resistance of the stem (Zhang et al., 2016; Song et al., 2021). This study conducted multi-environment genetic loci discovery across six environmental locations, identifying a total of 118 SNP markers distributed on chromosomes 1A, 1B, 1D, 2A, 2B, 2D, 3A, 3B, 4A, 4B, 4D, 5A, 5B, 6A, 6B, 6D, 7A, 7B, and 7D. The phenotypic variation explained by individual markers ranged from 8.99% to 85.48%. The widespread distribution of the 118 SNP markers suggests that the stem diameter trait in wheat is regulated by multiple genes, which may interact through a complex network to collectively influence the phenotype. The AX-110671322 SNP marker located on chromosome 1D may be linked to the known lodging resistance gene Rht-B1b, while other markers may regulate cell wall biosynthesis or hormone signaling pathways, further affecting stem diameter development. In the 1960s, the application of dwarfing genes such as Rht1 and Rht2 significantly enhanced wheat lodging resistance and increased yield, triggering the agricultural ‘Green Revolution.’ Current research has confirmed that the Rht-B1 and Rht-D1 genes located on chromosomes 4B and 4D have become core genetic resources for global wheat breeding (Xu et al., 2023). Genome-wide studies have revealed that plant height-associated loci are widely distributed across all chromosomes except for 3D, 4A, 5D, 7A, 7B, and 7D. Notably, the AX-109852602 locus on chromosome 7A overlaps with the AX-111557672 locus reported by Yu et al. (2024), while the AX-94584919 locus on chromosome 1A shows spatial overlap with the AX-109819835 locus identified by Wang (2022). These overlapping regions suggest that they may have conserved functions in stem diameter regulation, providing new insights for the analysis of lodging resistance genetic mechanisms and the development of molecular markers.

Functional analysis of candidate genes

Through the retrieval of stable genetic SNP markers in the Chinese Spring wheat genome database, seven candidate genes most likely associated with the second internode stem diameter at the base of the wheat stem were identified. These candidate genes are involved in several key biological processes, primarily related to the fibrillin family, calcium ion binding, the disease resistance protein RPM1, and F-box domain-containing proteins. The fibrillin family genes play a crucial role in cell wall construction and strengthening, while calcium ions are essential in cell signaling, regulating plant responses to environmental stresses. RPM1, as a resistance gene, is involved in wheat’s immune response to pathogens, and F-box proteins are important in plant protein degradation and signal transduction (Zhu et al., 2008; Xiong et al., 2021). The gene TraesCS1A02G422100, located on chromosome 1A, encodes Asparagine synthase, potentially regulating cytokinin precursor synthesis through amino acid metabolism, thereby affecting cell division and expansion. Another gene on chromosome 1A, TraesCS1A01G058600, is associated with Peroxisomal membrane 22 kDa (Mpv17/PMP22) family proteins. On chromosome 2D, TraesCS2D01G488200 encodes a Plastid-lipid associated protein PAP/fibrillin family protein, which plays an important role in cell wall construction and reinforcement. Fibrillin proteins affect the mechanical strength of the stem cell wall through cellulose microfibril alignment and deposition. This finding complements previous results from associative transcriptomics, which suggested that xylan acetylation regulates stem strength by influencing the interaction of cell wall polysaccharides, while fibrillin family proteins mediate the structural composition of the cell wall (Miller et al., 2016; Zhou et al., 2023). This synergistic interaction indicates that the mechanical properties of the cell wall are a complex trait regulated by a multi-gene network, involving polysaccharide modifications and structural protein assembly. TraesCS4D01G281500, located on chromosome 4D, encodes an F-box family protein. This protein, a key component of the ubiquitin-proteasome degradation pathway, targets auxin signaling inhibitors (such as AUX/IAA proteins), regulating the expression of auxin-responsive genes (Harper et al., 2016; Miller et al., 2016). This finding complements previous studies on the involvement of auxin response factors (ARF) and the SAUR gene cluster in plant height regulation, suggesting that stem diameter and plant height may share overlapping hormonal signaling pathways. Additionally, the identification of calcium ion-binding genes further supports the role of calcium signaling in responding to environmental stimuli, such as mechanical stress, and regulating cell expansion, offering insights into the adaptive mechanisms of wheat stems under lodging stress. On chromosome 5A, TraesCS5A01G073300 encodes the disease resistance protein RPM1. Its discovery suggests a potential link between disease resistance and stem development. RPM1 might indirectly modulate lignin deposition in the cell wall via the activation of the salicylic acid (SA) signaling pathway. It is well - documented that this pathway can induce the expression of genes implicated in lignin synthesis (Peng et al., 2014). This cross-regulation mechanism implies that wheat may enhance stem mechanical strength in response to pathogen infections by remodeling its stem structure, forming a synergistic defense and developmental strategy. Additionally, TraesCS6B02G050900 on chromosome 6B encodes a protein kinase catalytic domain, likely a downstream component of the RPM1 signaling pathway, which may regulate cell wall-synthesizing enzymes via phosphorylation. Notably, TraesCS7A02G481800 on chromosome 7A encodes a protein containing the SNARE motif, which is involved in vesicle transport and may mediate the transport of cell wall precursors to the plasma membrane, providing new insights into improving stem diameter by regulating material transport efficiency.

The identified genes are involved in cell wall synthesis, cellulose accumulation, and plant hormone regulation. Genes related to fibrillin, carbohydrate synthesis, and cell division and expansion play a crucial role in stem diameter development (Kelbert et al., 2004; Peng et al., 2014). Moreover, the regulation of stem diameter is closely related to plant hormone balance, particularly the regulation of gibberellin and cytokinin. These findings provide more effective genetic improvement strategies for enhancing crop resistance to lodging, increasing yield, and improving stem physiological traits.

Conclusion

In this study, a Q+K mixed linear model combined with a 55K SNP array was used to perform a genome-wide association analysis of stem diameter traits at the second internode of the basal stem in 239 wheat varieties (lines) from both domestic and international sources. A total of 118 SNP markers loci significantly associated with the trait (P ≤ 0.001) were successfully identified. Except for AX-111557672, AX-94584919, and AX-109819835, the other identified loci were not found in the previously reported articles and may be newly SNPs marker. These SNPs markers were distributed across multiple chromosomes of wheat, including 1A, 1B, 1D, 2A, 2B, 2D, 3A, 3B, 4A, 4B, 4D, 5A, 5B, 6A, 6B, 6D, 7A, 7B, and 7D. The phenotypic variance explained by individual SNP markers ranged from 8.09% to 29.14%. Furthermore, chromosomal clusters with large phenotypic effects were identified on chromosomes 2D, 4D, and 5A. By searching for stable, genetically significant SNP markers with large phenotypic effects in the Chinese Spring wheat genome database, seven candidate genes most likely related to the stem diameter trait at the second internode of the basal stem were ultimately selected. These findings provide essential theoretical support for further understanding the genetic mechanisms underlying stem diameter traits in wheat and offer potential molecular markers and target genes for wheat breeding improvements.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/supplementary material.

Author contributions

WW: Formal Analysis, Supervision, Project administration, Writing – original draft, Writing – review & editing, Resources, Conceptualization. NS: Writing – original draft, Supervision, Project administration, Resources, Writing – review & editing, Conceptualization. KZ: Supervision, Resources, Project administration, Conceptualization, Writing – review & editing, Writing – original draft. JS: Investigation, Data curation, Methodology, Writing – original draft. HF: Writing – original draft, Writing – review & editing, Resources, Project administration, Supervision, Conceptualization. GF: Supervision, Project administration, Writing – review & editing, Conceptualization, Writing – original draft, Resources, Formal Analysis. YG: Writing – original draft, Resources, Project administration, Conceptualization, Writing – review & editing, Supervision. TH: Writing – review & editing, Resources, Writing – original draft, Supervision, Project administration, Conceptualization. YD: Writing – original draft, Conceptualization, Resources, Project administration, Supervision, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research was funded by the Science and Technology Tackling in Henan Province(232102111132); Key Research Projects of Higher Education Institutions(23A210031, 23B520021); Department of Science and Technology Planning Project of Henan Province (242102110299);National Agricultural Biological Breeding Major Project Sub-theme ‘Breeding of New Ultra-high-yielding and Wide-adapted Wheat Varieties in Xinjiang Wheat Area’ (2023ZD040230310); National Agricultural Biological Breeding Major Project Sub-theme ‘Creation and Utilisation of New Materials of Ultra-high-yielding and Early-Maturing Wheat (2023ZD040230205); Xinjiang Uygur Autonomous Region Key Research and Development Programme ‘Research on the Mechanism of Important Trait Formation of Wheat in Xinjiang’ (2024B02007-2); Kashgar Regional Science and Technology Programme ‘Introduction of New Wheat Germplasm with Excellent Characteristics’. Introduction, Screening and Application of New Wheat Germplasm with Special Characteristics’ (KS2024044); Wheat Industry Technology System in Xinjiang Autonomous Region; Xinjiang Uygur Autonomous Region Key R&D Program Project “Research and Development of Key Technologies for Water Saving, Fertilizer Saving and Fertilizer Efficiency of Major Grain Crops in Yili River Valley” (2022B02015-2).

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.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

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.

References

Anderson, J. A., Wiersma, J. J., Reynolds, S. K., Conley, E. J., Stuart, N., Caspers, R., et al. (2025). Registration of ‘MN-Rothsay’ spring wheat with high grain yield and lodging resistance. J. Plant Registrat. 19, 13. doi: 10.1002/plr2.20400

Crossref Full Text | Google Scholar

Bainsla, N. K., Yadav, R., Singh, G. P., and Sharma, R. K. (2020). Additive genetic behavior of stem solidness in wheat (Triticum aestivum L.). Sci. Rep. 10, 9. doi: 10.1038/s41598-020-64470-x

PubMed Abstract | Crossref Full Text | Google Scholar

Cat, A. (2024). Evaluation of genetic variation and host resistance to wheat stem rust pathogen (Puccinia graminis f. sp. tritici) in bread wheat (Triticum aestivum L.) varieties grown in Türkiye. Peerj 12, 13. doi: 10.7717/peerj.17633

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, X., Liu, Q. Q., Zheng, B. Q., and Li, J. C. (2024). Optimized phosphorus application enhances wheat stem lodging resistance under spring low-temperature stress. Plants-Basel 13, 16. doi: 10.3390/plants13212980

PubMed Abstract | Crossref Full Text | Google Scholar

Crook, M. J. and Ennos, A. R. (1994). Stem and root characteristics associated with lodging resistance in four winter wheat cultivars. J. Agric. Sci. 123, 167–174. doi: 10.1017/s0021859600068428

Crossref Full Text | Google Scholar

Dhariwal, R., Hiebert, C. W., and Randhawa, H. S. (2022). QTL analysis identified two major all-internodes solidness loci from a completely solid-stemmed spring wheat line. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.1035620

PubMed Abstract | Crossref Full Text | Google Scholar

Ding, Y. D., Fang, H., Gao, Y. H., Fan, G. Q., Shi, X. L., Yu, S., et al. (2023). Genome-wide association analysis of time to heading and maturity in bread wheat using 55K microarrays. Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1296197

PubMed Abstract | Crossref Full Text | Google Scholar

Esposito, S., Taranto, F., Vitale, P., Ficco, D. B. M., Colecchia, S. A., Stevanato, P., et al. (2022). Unlocking the molecular basis of wheat straw composition and morphological traits through multi-locus GWAS. BMC Plant Biol. 22, 19. doi: 10.1186/s12870-022-03900-6

PubMed Abstract | Crossref Full Text | Google Scholar

Genievskaya, Y., Turuspekov, Y., Rsaliyev, A., and Abugalieva, S. (2020). Genome-wide association mapping for resistance to leaf, stem, and yellow rusts of common wheat under field conditions of South Kazakhstan. Peerj 8, 23. doi: 10.7717/peerj.9820

PubMed Abstract | Crossref Full Text | Google Scholar

Gong, Y.-R. (1982). Falls of wheat and preventive measures. Hubei. Agric. Sci. 09), 4–7.

Google Scholar

Harper, A. L., Trick, M., He, Z. S., Clissold, L., Fellgett, A., Griffiths, S., et al. (2016). Genome distribution of differential homoeologue contributions to leaf gene expression in bread wheat. Plant Biotechnol. J. 14, 1207–1214. doi: 10.1111/pbi.12486

PubMed Abstract | Crossref Full Text | Google Scholar

Kelbert, A. J., Spaner, D., Briggs, K. G., and King, J. R. (2004). The association of culm anatomy with lodging susceptibility in modern spring wheat genotypes. Euphytica 136, 211–221. doi: 10.1023/B:EUPH.0000030668.62653.0d

Crossref Full Text | Google Scholar

Khan, H., Krishnappa, G., Kumar, S., Devate, N. B., Rathan, N. D., Kumar, S., et al. (2024). Genome-wide association study identifies novel loci and candidate genes for rust resistance in wheat (Triticum aestivum L.). BMC Plant Biol. 24, 18. doi: 10.1186/s12870-024-05124-2

PubMed Abstract | Crossref Full Text | Google Scholar

Kong, X., Zhao, G. S., Sun, X. R., and Fu, Y. H. (2024). Coordinating lodging incidence and grain yield through wheat genetic diversity. Field Crops Res. 315, 13. doi: 10.1016/j.fcr.2024.109468

Crossref Full Text | Google Scholar

Li, J. G., Zhao, P. M., Zhao, L. Y., Chen, Q., Nong, S. K., Li, Q., et al. (2024). Integrated VIS/NIR spectrum and genome-wide association study for genetic dissection of cellulose crystallinity in wheat stems. Int. J. Mol. Sci. 25, 13. doi: 10.3390/ijms25053028

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, K., Harrison, M. T., Yan, H. L., Liu, D. L., Meinke, H., Hoogenboom, G., et al. (2023). Silver lining to a climate crisis in multiple prospects for alleviating crop waterlogging under future climates. Nat. Commun. 14, 13. doi: 10.1038/s41467-023-36129-4

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, Y. X., Lin, Y., Gao, S., Li, Z. Y., Ma, J., Deng, M., et al. (2017). A genome-wide association study of 23 agronomic traits in Chinese wheat landraces. Plant J. 91, 861–873. doi: 10.1111/tpj.13614

PubMed Abstract | Crossref Full Text | Google Scholar

Miller, C. N., Harper, A. L., Trick, M., Werner, P., Waldron, K., and Bancroft, I. (2016). Elucidation of the genetic basis of variation for stem strength characteristics in bread wheat by Associative Transcriptomics. BMC Genomics 17, 11. doi: 10.1186/s12864-016-2775-2

PubMed Abstract | Crossref Full Text | Google Scholar

Mostafaie, P., Afjeh, S. S., Ahmadi, A., and Oveisi, M. (2024). Exogenous spermine promotes grain filling in spring wheat by enhancing the accumulation and remobilization of stem reserves under drought stress conditions. J. Plant Growth Regul. 20, 2164–2183. doi: 10.1007/s00344-024-11535-0

Crossref Full Text | Google Scholar

Niu, Y. A., Chen, T. X., Zhao, C. C., Guo, C., and Zhou, M. X. (2022). Identification of QTL for stem traits in wheat (Triticum aestivum L.). Front. Plant Sci. 13. doi: 10.3389/fpls.2022.962253

PubMed Abstract | Crossref Full Text | Google Scholar

Peng, D. L., Chen, X. G., Yin, Y. P., Lu, K. L., Yang, W. B., Tang, Y. H., et al. (2014). Lodging resistance of winter wheat (Triticum aestivum L.): Lignin accumulation and its related enzymes activities due to the application of paclobutrazol or gibberellin acid. Field Crops Res. 157, 1–7. doi: 10.1016/j.fcr.2013.11.015

Crossref Full Text | Google Scholar

Rasheed, A., Xia, X., Ogbonnaya, F., Mahmood, T., Zhang, Z., Mujeeb-Kazi, A., et al. (2014). Genome-wide association for grain morphology in synthetic hexaploid wheats using digital imaging analysis. BMC Plant Biol. 14, 128. doi: 10.1186/1471-2229-14-128

PubMed Abstract | Crossref Full Text | Google Scholar

Saint Pierre, C., Trethowan, R., and Reynolds, M. (2010). Stem solidness and its relationship to water-soluble carbohydrates: association with wheat yield under water deficit. Funct. Plant Biol. 37, 166–174. doi: 10.1071/fp09174

Crossref Full Text | Google Scholar

Song, P. B., Wang, X., Wang, X. X., Zhou, F., Xu, X. W., Wu, B. J., et al. (2021). Application of 50K chip-based genetic map to QTL mapping of stem-related traits in wheat. Crop Pasture Sci. 72, 105–112. doi: 10.1071/cp20326

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

Tian, X. L., Wen, W. E., Xie, L., Fu, L. P., Xu, D. G., Fu, C., et al. (2017). Molecular mapping of reduced plant height gene rht24 in bread wheat. Front. Plant Sci. 8. doi: 10.3389/fpls.2017.01379

PubMed Abstract | Crossref Full Text | Google Scholar

Tolossa, M., Adugna, G., and Hundie, B. (2024). Identification of physiological race of wheat stem rust (Puccinia graminis f.sp. tritici) and reaction of bread wheat cultivars to emerging dominant races. Trop. Plant Pathol. 49, 786–795. doi: 10.1007/s40858-024-00670-4

Crossref Full Text | Google Scholar

Wang, Y. (2022). Genome-wide linkage analysis of traits related to resistance to heading in common wheat. J. Plant Genet. Resour. 23, 13.

Google Scholar

Wang, W., Cheng, Y. K., Ren, Y., Zhang, Z. H., and Geng, H. W. (2022a). Prediction of chlorophyll content in multi-temporal winter wheat based on multispectral and machine learning. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.896408

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, W., Gao, X., Cheng, Y. K., Ren, Y., Zhang, Z. H., Wang, R., et al. (2022b). QTL mapping of leaf area index and chlorophyll content based on UAV remote sensing in wheat. Agriculture-Basel 12, 19. doi: 10.3390/agriculture12050595

Crossref Full Text | Google Scholar

Wang, W., Sun, N., Bai, B., Wu, H., Cheng, Y. K., Geng, H. W., et al. (2024). Prediction of wheat SPAD using integrated multispectral and support vector machines. Front. Plant Sci. 15. doi: 10.3389/fpls.2024.1405068

PubMed Abstract | Crossref Full Text | Google Scholar

Xiong, H. C., Li, Y. T., Guo, H. J., Xie, Y. D., Zhao, L. S., Gu, J. Y., et al. (2021). Genetic mapping by integration of 55K SNP array and KASP markers reveals candidate genes for important agronomic traits in hexaploid wheat. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.628478

PubMed Abstract | Crossref Full Text | Google Scholar

Xu, D. A., Bian, Y. J., Luo, X. M., Jia, C. F., Hao, Q. L., Tian, X. L., et al. (2023). Dissecting pleiotropic functions of the wheat Green Revolution gene Rht-B1b in plant morphogenesis and yield formation. Development 150, 13. doi: 10.1242/dev.201601

PubMed Abstract | Crossref Full Text | Google Scholar

Xu, F., Shi, R. J., Liu, L. L., Li, S. F., Wang, J. M., Han, Z. H., et al. (2024). Fusarium pseudograminearum biomass and toxin accumulation in wheat tissues with and without Fusarium crown rot symptoms. Front. Plant Sci. 15. doi: 10.3389/fpls.2024.1356723

PubMed Abstract | Crossref Full Text | Google Scholar

Yu, H., Wu, W., Pei, X., Liu, X., Deng, H., Li, X., et al. (2024). Transcriptome sequencing and genome-wide association analysis of wheat stalk traits. J. Crop Sci. 50, 2187–2206. doi: 10.3724/sp.j.1006.2024.31076

Crossref Full Text | Google Scholar

Zhang, Y., Xu, W. G., Wang, H. W., Fang, Y. H., Dong, H. B., and Qi, X. L. (2016). Progress in improving stem lodging resistance of Chinese wheat cultivars. Euphytica 212, 275–286. doi: 10.1007/s10681-016-1768-1

Crossref Full Text | Google Scholar

Zhou, C., Xiong, H. C., Fu, M., Guo, H., Zhao, L. S., Xie, Y. D., et al. (2023). Genetic mapping and identification of rht8-b1 that regulates plant height in wheat. BMC Plant Biol. 23, 1–11. doi: 10.1186/s12870-023-04343-3

PubMed Abstract | Crossref Full Text | Google Scholar

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

Crossref Full Text | Google Scholar

Keywords: wheat, stem thickness, collapse, GWAS, candidate genes

Citation: Wang W, Sun N, Zhao K, Song J, Fang H, Fan G, Gao Y, Huang T and Ding Y (2025) Genome-wide association analysis of wheat stem traits using 55K microarrays. Front. Plant Sci. 16:1635721. doi: 10.3389/fpls.2025.1635721

Received: 26 May 2025; Accepted: 03 July 2025;
Published: 24 July 2025.

Edited by:

Linghe Zeng, United States Department of Agriculture (USDA), United States

Reviewed by:

Maria Khalid, National University of Sciences and Technology (NUST), Pakistan
Arpit Gaur, Montana State University, United States
Sukumar Taria, Central Agroforestry Research Institute (CAFRI), India

Copyright © 2025 Wang, Sun, Zhao, Song, Fang, Fan, Gao, Huang and Ding. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yindeng Ding, MTUyOTI0MzQ0MjVAMTYzLmNvbQ==

These authors have contributed equally to this work

Disclaimer: 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.