ORIGINAL RESEARCH article

Front. Plant Sci., 17 January 2025

Sec. Plant Breeding

Volume 15 - 2024 | https://doi.org/10.3389/fpls.2024.1471636

Performing whole-genome association analysis of winter wheat plant height using the 55K chip

  • 1. Institute of Grain Crops, Xinjiang Academy of Agricultural Sciences, Xinjiang, China

  • 2. College of Agriculture, Xinjiang Agricultural University, Xinjiang, China

  • 3. Institute of Crop Variety Resources, Xinjiang Academy of Agricultural Sciences, Xinjiang, China

  • 4. Department of Computer Science and Information Engineering, Anyang Institute of Technology, Anyang, China

  • 5. Cotton Research Institute, Chinese Academy of Agricultural Sciences, Anyang, China

  • 6. Yili Prefecture Institute of Agricultural Sciences/Yili Prefecture Key Laboratory of Crop Breeding and Quality Detection, Yining, China

Abstract

Plant height is a critical agronomic that affects both plant architecture and yield. To decipher the genetic mechanisms underlying winter wheat plant height and identify candidate genes associated with this trait, 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. A Wholegenome association analysis (GWAS) of wheat plant height was conducted utilizing the MLM (Q+K) model within TASSLE software. The results revealed significant phenotypic variation in wheat plant height across different years, with coefficients of variation ranging from 0.96% to 1.97%. Additionally, there was a strong correlation in plant height measurements between different years. GWAS identified 44 SNP markers significantly associated with wheat plant height across various environments (P ≤ 0.00001), predominantly distributed on chromosomes 1B, 1D, 2A, 2B, 2D, 3B, 3D, 4A, 4B, 6B, 6D, and 7D, explaining individual phenotypic variance rates ranging from 5.00% to 11.11%. Further, by mining association loci with substantial phenotypic effects and stability across multiple environments, seven candidate genes related to wheat plant height have been identified. This study provides new genetic markers and resources for improving wheat plant height.

Introduction

As a highly adaptable and widely grown food crop in the world, wheat provides about 21 percent of food calories and 20 percent of protein for humans (Maccaferri et al., 2014). Its high nutritional value and unique gluten properties permit the preparation of a wide variety of food products, therefore it has become one of the most important food crops (Wang et al., 2022; Liu et al., 2023). The key goal of their breeding has been to increase yield (Shiferaw et al., 2013; Jin et al., 2019; Wang et al., 2024). Wheat plant height is a key agronomic trait affecting plant conformation and yield. As an aboveground support organ, the wheat stem connects various physiological activities between the leaves and the grain. If wheat plant height is too high, it can lead to plant collapse, which directly reduces yield, while too short a plant can lead to overcrowding of canopy leaves, which affects photosynthetic rate, thus making biomass unable to meet the demand, which in turn reduces yield (Wang and Li, 2008; Godin, 2000). Plant height is simultaneously regulated by both genes and the environment (Li et al., 2020). Dwarf wheat has dramatically increased wheat yields in practice due to its strong resistance to felling (Miao et al., 2022). About 30 dwarf genes have been reported, mainly on chromosomes 2A, 2B, 2D, 3B, 4B, 4D, 5A, 6A, 7A, and 7B, of which only the dwarf genes Rht1, Rht2, Rht8, Rht11, Rht13, and Rht24 have been applied in wheat breeding both at home and abroad, and have affected agronomic traits to varying extents to improve yields (Sun et al., 2017). Throughout the world, 70% of wheat varieties contain Rhtl or Rht2. In the Huanghuai wheat region of China, in addition to Rhtl and Rht2, a few short-stalked genes, such as Rht8 and Rht24, have no significant negative effect on yield traits while reducing plant height, and have been selected for use in breeding (Li et al., 2018; Liu et al., 2017). The discovery and utilization of new excellent dwarf genes not only can enrich the genetic resources of existing dwarf genes, but also is an important basis for the sustainable development of wheat highyield breeding. It is of great significance to analyse the genetic effect of plant height and to explore the functional genes related to plant height for the improvement of wheat varieties.

To the present, many scholars at home and abroad have conducted QTL localization studies on wheat plant height using linkage analysis and association analysis methods, and identified a large number of primary QTL loci for controlling plant height. Through QTL chain localization, Gao et al. (2015) used a population of recombinant inbred lines (RILs) constructed with ‘Zhou 8425’ and ‘China Spring’ as the parents as the material to localize six stable QTLs related to plant height, and a single QTL could explain 2.3%-33.2% of the A single QTL could explain 2.3%-33.2% of the phenotypic variation, among which two QTLs were colocalized with the short culm genes Rhtl and Rht2, respectively. Xiong et al. (2021) utilized a population of RILs constructed from an early tasselling mutant (eh1) and ‘Rotary Selection 987’ to identify seven QTLs located on chromosomes 2A, 4A, 4B, and 6B for controlling the height of wheat plants, with phenotypic variance contributions ranging from 5.5% to 34.4%. Li et al. (2020) obtained four QTLs controlling plant height and predicted four candidate genes encoding IAAamino acid hydrolase ILR1, growth regulator, flowering promoter, and phosphosucrose synthase in a population of RILs constructed from ‘Semimanzan’ and ‘Jeimai 22’. Li et al. (2020) identified three stable QTLs tightly linked to plant height using two recombinant inbred line (RIL) populations. They developed competitive allelespecific PCR (kompetitive allelespecific PCR, KASP) markers for OSE-Liscau-2CN-5A within this interval. Three candidate genes, all encoding cytochrome P450 and highly expressed in wheat stems, were discovered in this region. Qu et al. (2021) identified 21 QTLs closely linked to plant height in three populations of RILs, and identified a gene encoding a gibberellinregulated protein in the qPH-4D-1 interval of chromosome 4D, which is related to the regulation of plant height in wheat. The association analysis has been widely used in the genetic analysis of wheat plant height. Sukumaran et al. (2015) conducted a genome-wide association study (GWAS) on 287 highyielding wheat varieties, identifying two significant SNP loci associated with plant height. These loci were located on chromosomes 1A and 6A, contributing 6% and 7% of the phenotypic variance, respectively. Sun et al. (2017) conducted a genome-wide association study (GWAS) using SNP markers and 163 common wheat varieties from the Huanghuai region, identifying a total of 143 significant loci associated with wheat plant height. These loci were distributed across chromosomes 1B, 2A, 2B, 2D, 4B, 4D, and 6D, explaining 35.8% to 65.6% of the phenotypic variance. Rabbi et al. (2021), on the other hand, performed GWAS analysis on 361 US wheat varieties, detecting 20 significant SNP loci associated with plant height. These loci explained 4.54% to 48.01% of the phenotypic variance and predicted four candidate genes related to plant height on chromosomes 2A, 3A, and 6A.

Although some progress has been made in the basic research on plant height genetics, due to the large wheat genome and the complexity of the quantitative genetic mechanism of plant height, most of the genetic variation cannot be directly applied to molecular breeding, and its genetic analysis needs to be further expanded and deepened. Therefore, in this study, 239 wheat varieties (lines) were used as materials, and 16649 SNP markers covering the whole genome were used for genomewide association analysis of plant height in four environments, to detect the significant association loci and discover the relevant candidate genes, with a view to providing theoretical basis for the molecular breeding of plant height in wheat, and the excavation and utilisation of the candidate genes.

Materials and methods

Materials

A total of 239 natural populations consisting of wheat breeding varieties and foreign introduced varieties promoted in winter wheat areas of China were used as test materials. Of these, 213 were from China, including 6 from Anhui, 34 from Beijing, 19 from Hebei, 19 from Henan, 5 from Jiangsu, 11 from Shandong, 1 from Shanxi, 3 from Shaanxi, 1 from Tianjin and 115 from Xinjiang. There were 26 foreign materials, including 25 from the USA and 1 from Ukraine. All these materials could grow and develop normally in the test site.

Multi-environment trials

The experiments were conducted from September 2019 to July 2022 at the ZePu Breeding base of Xinjiang Academy of Agricultural Sciences (77°16’17.22”N, 38°11’21.65”E), and in June 2021 at the Anningqu base of Xinjiang Academy of Agricultural Sciences (43°58’53.38”N, 87°30’17.72”E). The ZePu region exhibits a continental warm-temperate arid climate characterized by low annual precipitation. Overall, rainfall levels remain minimal, consistent with the arid nature of the Xinjiang region. The area experiences distinct four seasons, marked by significant diurnal temperature variations, with summers being notably warm and winters relatively cold. The Anningqu area is characterized by a temperate continental climate, with perennial low precipitation levels primarily concentrated during the summer months of June to August. Spring (March to May) experiences a gradual rebound in temperatures, although the rate of warming is relatively slow. In summer (July), the region is influenced by subtropical high-pressure systems, resulting in hot weather. Autumn (September to October) sees a decrease in rainfall alongside a gradual decline in temperatures. Winter (November to March) is dominated by prevailing northwesterly winds, leading to cold conditions, with snowfall being a prominent feature. During this season, temperatures can drop significantly, reaching as low as -20°C. In these trials, ZePu in Xinjiang served as E1 in 2020, E2 in 2021, and E4 in 2022, while Anningqu in Xinjiang was E3 in 2021. The best linear unbiased prediction (BLUP) of value were obtained using R 4.0 (the Matrix and lme4 packages). Each environment point utilized a randomized complete block design over three years. Each plot consisted of two rows, each 1 meter long and spaced 0.2 meters apart, with a sowing density of 350,000 seeds per hectare, oriented in a north-south direction. Standard field management practices including light, water, fertilizer, and other agronomic practices were followed consistently across the three replicates. The cultivation and production conditions were identical for each replicate.

Measurement of plant height

After the milky stage of wheat, a straightedge was purchased to measure the length from the ground to the top of the spike (excluding the awns) in cm. 5 plants with similar growth were measured in each plot, and the average was taken as the height of the plant in the plot (Supplementary Table S2).

Statistical analysis

The process of phenotyping was mainly ANOVA as well as descriptive statistics of phenotypes, significance test of difference and correlation analysis, in addition all statistical analyses of data were implemented in IDE Spyder under Anacondas3, using Python 3.8.8 for data processing, on a computer with an Intel i7-6800 K 3.40 GHz CPU, 16 GB of RAM, and Nvidia GeForce GTX 2080Ti on a graphics workstation running the Win10 operating system. Statistical analyses, correlation analyses and significance of difference tests for wheat plant height were performed using the application software Pandas 1.3.2, Matplotlib 3.4.2, Scikit-Learn 0.24.2 and SPSS 21.0.

GWAS analysis

A total of 16649 SNP markers were identified and used for subsequent analysis as previously described by Ding et al. (2023) (Supplementary Table S3). The plant height (PH) in the four different environments were used for GWAS analysis. A previous study analysed the population structure and linkage disequilibrium (LD). The mixed linear model (MLM) in TASSEL 5.0 was used for GWAS, and the population structure was considered as a fixed effect (Ding et al., 2023). The Bonferroni threshold of SNP significance was P ≤ 0.00001. Manhattan plots were generated by using the “CMplot” package for R (Zhu et al., 2008).

Results

Analysis of phenotypic variation of plant height

Wheat plant height was phenotyped at four environmental sites, E1, E2, E3, and E4. The data were evaluated through four dimensions: mean (denoted by μ), median (denoted by median), coefficient of variation (denoted by cv), and standard deviation (denoted by σ). As seen in Figure 1, the mean μ of plant height in the E1 environment was 77.18, the median was 75.12, the coefficient of variation was 15.0%, and the standard deviation was 11.56; the mean μ of plant height in the E2 environment was 83.94, the median was 82.38, and the coefficient of variation was 11.3%, and the standard deviation was 9.49; and the mean μ of the germination-twigging stage in the E3 environment was 69.43, the median was 68.25, coefficient of variation was 10.80%, and standard deviation was 7.48; the mean μ of plant height under E4 environment was 72.72, median was 71.90, coefficient of variation was 12.50%, and standard deviation was 9.08 (Table 1). The ANOVA indicated that the environment and genotype interaction had significant effects on PH (P < 0.01) (Supplementary Table S1). Overall, the data under different environments showed continuous and normal distribution, which was in line with the characteristics of typical quantitative traits. These results provide reliable phenotypic data for the genetic analysis of PH.

Figure 1

Table 1

EnvsMin
(cm)
Max
(cm)
Average (cm)Standard DevitionCoefficient of Variation
E150.41146.3877.1811.5915.00%
E255.13115.0683.949.5111.30%
E352.2098.1069.437.5010.80%
E453.81107.0272.729.0812.50%

Descriptive statistics of wheat plant height.

Correlation analysis of plant height

The results of correlation analyses of plant height at four environmental points, E1, E2, E3 and E4, are shown in Figure 2. The figure compares the correlations of plant height in winter wheat under different environments, and the results show that the correlations were significant (p<0.001) in most cases. Specifically, the correlations of plant height ranged from -0.59 to -0.65 in E1, 0.65 to 0.79 in E2, 0.59 to 0.72 in E3, and 0.65 to 0.79 in E4. Taken together, the range of correlations of plant height in different environments was from 0.59 to 0.79, which were all significant (p<0.001).

Figure 2

GWAS analysis of plant height

Using TASSEL 5.0 software, a genome-wide association analysis was conducted on 239 wheat varieties (lines) using 16,649 high-quality SNP markers identified from a 55K SNP chip and combined with SPAD content measurements. The MLM (Q+K) mod-el was employed, with markers considered significantly associated with traits when (P ≤ 0.00001). Markers detected across multiple environments were considered stable genetic loci (Figure 3, Table 2). The GWAS results revealed a total of 44 markers associated with plant height, distributed across chromosomes 1B, 1D, 2A, 2B, 2D, 3B, 3D, 4A, 4B, 6B, 6D, and 7D. Individual markers explained phenotypic variance rates ranging from 5.00% to 11.11%. Notably, AX-110425403 located on chromosome 1B was detected in environments E1, E2, E3 and BLUP, with phenotypic variance rates ranging from 5.00% to 9.68%. AX-108850618 located on chromosome 1D was detected in environments E2 BLUP, with phenotypic variance rates ranging from 5.09% to 5.37%. AX-94748732 on chromosome 3D was detected in environments E1,E4 and BLUP, with phenotypic variance rates of 4.04% and 5.39%, respectively. AX-109306628 on chromosome 4A was detected in environments E2 and BLUP, with phenotypic variance rates of 5.41% and 5.47%, respectively. AX-109392949 on chromosome 4B was detected in environments E4 and BLUP, with phenotypic variance rates of 6.29% and 6.95%, respectively. AX-111046577 on chromosome 4B was detected in environments E3 and BLUP, with phenotypic variance rates of 5.09% and 5.24%, respectively. AX-111230804 on chromosome 4B was detected in environments E3 and BLUP, with phenotypic variance rates of 5.09% and 5.24%, respectively. AX-89577308 on chromosome 4B was detected in environments E4 and BLUP, with phenotypic variance rates of 4.61% and 5.21%, respectively. AX-109934491 on chromosome 5A was detected in environments E2 and BLUP, with phenotypic variance rates of 5.29% and 5.38%, respectively.AX-110644627 on chromosome 5A was detected in environments E4 and BLUP, with phenotypic variance rates of 6.03% and 6.74%, respectively. AX-110616775 on chromosome 6B was detected in environments E4 and BLUP, with phenotypic variance rates of 4.70% and 5.04%, respectively. AX-111593090 on chromosome 6B was detected in environments E2 and BLUP, with phenotypic variance rates of 4.96% and 5.07%, respectively. Other markers were detected in single environments.

Figure 3

Table 2

MarkerChr.Position (Mb)P-valueR2 (%)Environment
AX-1104254031B581.1666589.9313E-065.00E1
AX-1104254031B581.1666583.6708E-065.44E2
AX-1104254031B581.1666582.0747E-109.68E4
AX-1104254031B581.1666581.86E-087.73BLUP
AX-1112284201B512.311472.0967E-065.68E1
AX-1088506181D42.6122938.1813E-065.09E2
AX-1088506181D42.6122934.22E-065.37BLUP
AX-1094727122A767.9399961.851E-1110.73E1
AX-1104727532A227.2669113.0102E-076.52E1
AX-1098526022B778.0142680.0000051945.28E4
AX-1087644062D641.1223384.1945E-087.38E1
AX-1094628412D71.7687886.1277E-065.21E4
AX-1094858642D643.0491517.8226E-1211.11E1
AX-1114615892D645.368917.8226E-1211.11E1
AX-1119560722D34.4288381.1099E-087.95E4
AX-896447842D73.1111037.0528E-065.15E4
AX-945762662D646.03238.4108E-076.08E1
AX-1099119143B701.4024178.3294E-065.08E1
AX-1088912933D418.7882682.2075E-065.66E4
AX-944651903D603.1162861.6721E-065.78E1
AX-947487323D570.8051137.2728E-065.14E1
AX-945362913D571.369.09E-054.04BLUP
AX-949422513D569.5242324.0844E-065.39E4
AX-1088580244A237.6877347.8989E-076.10E1
AX-1093066284A88.0194543.8606E-065.41E2
AX-1093066284A88.0194543.39E-065.47BLUP
AX-1098441864A260.1864777.4415E-1211.13E1
AX-1099069214A175.8503594.4213E-087.35E1
AX-1104534584A110.2809624.4504E-1110.35E1
AX-1092700004B604.8668517.7292E-065.11E4
AX-1092995134B495.5951814.1025E-1110.39E1
AX-1093846914B604.9248327.7292E-065.11E4
AX-1093929494B481.5075231.20E-076.92E4
AX-1093929494B481.5075235.45E-076.26BLUP
AX-1101446774B605.1565726.7405E-065.17E4
AX-1110465774B522.8013285.7912E-065.24E3
AX-1110465774B522.8013288.09E-065.09BLUP
AX-1112308044B526.1591515.7912E-065.24E3
AX-1112308044B526.1591518.09E-065.09BLUP
AX-1116496984B609.4152968.586E-065.07E4
AX-895773084B559.6478846.1248E-065.21E4
AX-895773084B559.6478842.43E-054.61BLUP
AX-1099166895A10.3011123.4886E-087.46E1
AX-1099344915A618.8456095.1763E-065.29E2
AX-1099344915A618.8456094.20E-065.38BLUP
AX-1106446275A3.5559.4213E-076.03E4
AX-1106446275A3.5551.80E-076.74BLUP
AX-1088035845B707.2616361.1583E-076.94E1
AX-1088298356B151.1965920.0000072975.14E4
AX-1106167756B708.6595419.1494E-065.04E4
AX-1106167756B708.6595411.99E-054.70BLUP
AX-1115930906B526.7409388.4337E-065.07E2
AX-1115930906B526.7409381.11E-054.96BLUP
AX-1110289676D29.4671677.8226E-1211.11E1
AX-1119311366D222.5729885.7479E-065.24E1
AX-897428387D199.4138384.2648E-087.37E4

Information of significantly associated loci for wheat plant height.

Functional prediction of candidate genes for strain height

SNP markers with large phenotypic effect values that could be stably inherited were searched in the Chinese Spring Genome Database of common wheat and BLASTx sequence comparison was performed in the NCBI database, and a total of seven candidate genes most likely to be associated with plant height were mined (Table 3). The candidate genes for plant height were mainly related to photosynthesis, potassium transporters, zinc finger family proteins and expressed proteins containing F-box structural domains in crops. Among them, the gene TraesCS1B01G350500 localised on chromosome 1B, encodes a transcription repressor; TraesCS2A01G570600 localised on chromosome 2A, is associated with Agmatine coumaroyltransferase-2; and TraesCS2D01G570600 localised on chromosome 2D, is associated with Agmatine coumaroyltransferase-2. TraesCS2D01G581800 on chro-mosome 2D encodes GRF zinc finger family protein; TraesCS3D01G466000 on chromosome 3D encodes F-box domain containing protein; TraesCS4B01G466000 on chromosome 4B encodes F-box domain containing protein; and TraesCS6D01G061800 on chromosome 6D encodes Potassium transporter.

Table 3

MarkerChr.Position (Mb)GenesGene annotation
AX-1104254031B581.17TraesCS1B01G350500transcription repressor
AX-1094727122A767.94TraesCS2A01G570600Agmatine coumaroyltransferase-2
AX-1114615892D645.37TraesCS2D01G581800GRF zinc finger family protein
AX-947487323D570.81TraesCS3D01G466000F-box domain containing protein
AX-1092995134B495.6TraesCS4B01G238200Photosystem I reaction center subunit IV
AX-1106446275A3.56TraesCS5A01G013800F-box domain containing protein
AX-1110289676D29.47TraesCS6D01G061800Potassium transporter

Information of candidate genes for wheat plant height.

Discussion

Association analysis of plant height

With the rapid development of biology and bioinformatics, GWAS analysis has become an important way to study quantitative traits in plants, and the mining of genes related to wheat plant height has been promoted to a greater extent. Classical genetic studies have shown that plant height is a complex trait jointly controlled by Mendelian and quantitative genes. In this study, multi-environmental genetic locus mining was carried out through four environmental sites, and a total of 44 SNP markers were mined, which were distributed on chromosomes 1B, 1D, 2A, 2B, 2D, 3B, 3D, 4A, 4B, 4A, 6B, 6D, and 7D, with a rate of variation of 5.00-11.11% for the individual inter-pretable phenotypes. In the 1960s, the utilization of dwarf genes such as Rht1 (Reduced Height 1) and Rht2 significantly enhanced wheat’s resistance to stunting and increased yield per unit area. This advancement played a pivotal role in the development of agricultural biotechnology and contributed to the onset of the ‘Green Revolution’ in agricultural production. Among the currently available genes related to the study of wheat plant height, the Rht-B1 and Rht-D1 are located on chromosomes 4B and 4D, respectively, and have been successfully used in wheat breeding worldwide (Xu et al., 2023). In addition to Rht-B1 and Rht-D1 genes (Bazhenov et al., 2015), other PH-related genes such as Rht 8, Rht 13, Rht 22, Rht 24, etc. have also been reported (Divashuk et al., 2012). Many loci were also identified as being located on all chromosomes except chromosomes 3D, 4A, 5D, 7A, 7B and 7D. In addition, the locus AX-109852602 on chromosome 2B mined in this study overlapped with the locus on plant height mined by Li et al. (2022) and the locus AX-111956072 on chromosome 2D overlapped with the locus on plant height mined by Guo et al. (2010).

Functional analysis of candidate genes

By querying stable inherited SNP markers against the common wheat Chinese Spring genome database, seven candidate genes most likely associated with plant height were identified. These candidate genes are primarily linked to photosynthesis, potassium transporters, zinc finger proteins, and expressed proteins containing the F-box domain. The gene TraesCS1B01G350500 located on chromosome 1B encodes a transcription repressor. Throughout the long evolutionary process, plants have developed intricate gene regulatory networks, where reversible transcription repression and relief responses play a crucial role in plant adaptation to the environment and meeting their own growth and development needs. Genes TraesCS2A01G570600 and Agmatine cou-maroyltransferase-2 associated with chromosome 2A; and TraesCS2D01G581800 encoding GRF zinc finger family protein on chromosome 2D. This protein participates in the physiological and biochemical regulatory mechanisms of plant growth and development. TraesCS3D01G466000 located on chromosome 3D and TraesCS5A01G013800 on chromosome 5A encode F-box domain containing proteins, which are expressed. The F-box protein family plays crucial roles in plant hormone signalling, light signal transduction, and physiological processes such as floral organ development (Kipreos and Pagano, 2000). TraesCS4B01G238200 on chromosome 4B encodes Photosystem I reaction centre subunit IV. Additionally, TraesCS6D01G061800 on chromosome 6D encodes a Potassium transporter. Potassium ions are essential monovalent cations for plant growth and development, being among the most abundant cations in plant cells (Himi et al., 2011). Potassium nutrient uptake is critically important for normal plant growth and development. Potassium ions participate in various processes including enzyme activation, regulation of cell turgor pressure, enhancement of substance synthesis and transport within plants, promotion of photosynthesis, and bolstering crop resistance to diseases and stresses.

Conclusion

This study utilized a Q+K mixed linear model and employed a 55K SNP chip to conduct a genome-wide association analysis of heading date and maturity in 239 domestic and international wheat varieties. A total of 44 SNP markers significantly associated with wheat plant height were identified (P ≤ 0.00001). These markers were distributed across chromosomes 1B, 1D, 2A, 2B, 2D, 3B, 3D, 4A, 4B, 6B, 6D, and 7D, individually explaining phenotypic variance rates ranging from 5.00% to 11.11%. Notably, in multiple environments, AX-110425403 on chromosome 1B, AX-94748732 on chromosome 3D, and AX-110644627 on chromosome 5A exhibited phenotypic variance rates of 5.00% to 9.68%. Subsequently, SNPs showing substantial phenotypic effects and stable heritability were queried against the common wheat Chinese Spring genomic database, leading to the discovery of 7 promising candidate genes most likely associated with wheat plant height.

Statements

Data availability statement

The datasets generated during and/or analyzed during the current study are available from the corresponding author upon reasonable request.

Author contributions

YD: Resources, Writing – original draft, Writing – review & editing. GF: Writing – review & editing. YG: Conceptualization, Writing – review & editing. TH: Data curation, Writing – review & editing. AZ: Writing – review & editing. SY: Writing – review & editing. LZ: Formal analysis, Writing – review & editing. XS: Writing – review & editing. SD: Writing – review & editing. JH: Writing – review & editing. WW: Writing – review & editing. JS: Writing – review & editing. NS: Writing – review & editing. HF: Writing – original draft, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by the Xinjiang Uygur Autonomous Region Major Science and Technology Special Project (2022A03004-3); Xinjiang Academy of Agricultural Sciences Youth Scientific and Technological Backbone Innovation Ability Cultivation Project (xjnkq-2021006); Xinjiang Uygur Autonomous Region Agricultural and Rural Affairs Department’s ‘Unveiling and Hanging the Marshal’ Tackling Project in 2022 (Early Maturity of Southern Xinjiang, High-yielding, high-quality winter wheat varieties and core technologies); Major Science and Technology Special Project of Xinjiang Uygur Autonomous Region (2021A02001-1); Key Research and Development Special Project of Xinjiang Uygur Autonomous Region (2022B02015-2); Key Disciplines of Anyang Institute of Technology; the Natural Science Foundation of Henan (242300420459); Autonomous Region Young Scientific and Technological Talents Project - ‘Tianshan Talents’ Training Programme (2022TSYCJC0051).

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/fpls.2024.1471636/full#supplementary-material

References

  • 1

    BazhenovM. S.DivashukM. G.AmagaiY.WatanabeN.KarlovG. I. (2015). Isolation of the dwarfing Rht-B1p (Rht17) gene from wheat and the development of an allele-specific PCR marker. Mol. Breed.35, 18. doi: 10.1007/s11032-015-0407-1

  • 2

    DingY.FangH.GaoY.FanG.ShiX.YuS.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

  • 3

    DivashukM. G.VasilyevA. V.BespalovaL. A.KarlovG. I. (2012). Identity of the Rht-11 and Rht-B1e reduced plant height genes. Russian J. Genet.48, 761763. doi: 10.1134/S1022795412050055

  • 4

    GaoF.WenW.LiuJ.RasheedA.YinG.XiaX.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. doi: 10.3389/fpls.2015.01099

  • 5

    GodinC. (2000). Representing and encoding plant architecture: a review. Ann. For. Sci.57, 413438. doi: 10.1051/forest:2000132

  • 6

    GuoZ.SongY.ZhouR.RenZ.JiaJ. (2010). Discovery, evaluation and distribution of haplotypes of the wheat. Ppd-D1 gene. New Phytol.185, 841851. doi: 10.1111/j.1469-8137.2009.03099.x

  • 7

    HimiE.MaekawaM.MiuraH.NodaK. (2011). Development of PCR markers for Tamyb10 related to R-1, red grain color gene in wheat. Theor. Appl. Genet.122, 15611576. doi: 10.1007/s00122-011-1555-2

  • 8

    JinY.TangJ.SablaniS. S. (2019). Food component influence on water activity of low-moisture powders at elevated temperatures in connection with pathogen control. Lwt112, 108257. doi: 10.1016/j.lwt.2019.108257

  • 9

    KipreosE. T.PaganoM. (2000). The F-box protein family. Genome Biol.1, 17. doi: 10.1186/gb-2000-1-5-reviews3002

  • 10

    LiA.HaoC.WangZ.GengS.JiaM.WangF.et al. (2022). Wheat breeding history reveals synergistic selection of pleiotropic genomic sites for plant architecture and grain yield. Mol. Plant15, 504519. doi: 10.1016/j.molp.2022.01.004

  • 11

    LiC.TangH.LuoW.ZhangX.MuY.DengM.et al. (2020). A novel, validated, and plant height-independent QTL for spike extension length is associated with yield-related traits in wheat. Theor. Appl. Genet.133, 33813393. doi: 10.1007/s00122-020-03675-0

  • 12

    LiS.TianY.WuK.YeY.YuJ.ZhangJ.et al. (2018). Modulating plant growth-metabolism coordination for sustainable agriculture. Nature560, 595600. doi: 10.1038/s41586-018-0415-5

  • 13

    LiuK.HarrisonM. T.YanH.LiuL.MeinkeH.HoogenboomG.et al. (2023). Silver lining to a climate crisis in multiple prospects for alleviating crop waterlogging under future climates. Nat. Commun.14, 765. doi: 10.1038/s41467-023-36129-4

  • 14

    LiuY.LinY.GaoS.LiZ.MaJ.DengM.et al. (2017). A genome-wide association study of 23 agronomic traits in Chinese wheat landraces. Plant J.91, 861873. doi: 10.1111/tpj.13614

  • 15

    MaccaferriM.Cane’M. A.SanguinetiM. C.SalviS.ColalongoM. C.MassiA.et al. (2014). A consensus framework map of durum wheat (Triticum durum Desf.) suitable for linkage disequilibrium analysis and genome-wide association mapping. BMC Genomics15, 873. doi: 10.1186/1471-2164-15-873

  • 16

    MiaoY.JingF.MaJ.LiuY.ZhangP.ChenT.et al. (2022). Major genomic regions for wheat grain weight as revealed by QTL linkage mapping and meta-analysis. Front. Plant Sci.13. doi: 10.3389/fpls.2022.802310

  • 17

    QuP.WangJ.WenW.GaoF.LiuJ.XiaX.et al. (2021). Construction of consensus genetic map with applications in gene mapping of wheat (Triticum aestivum L.) using 90K SNP array. Front. Plant Sci.12. doi: 10.3389/fpls.2021.727077

  • 18

    RabbiS. M. H. A.KumarA.Mohajeri NaraghiS.SimsekS.SapkotaS.SolankiS.et al. (2021). Genome-wide association mapping for yield and related traits under drought stressed and non-stressed environments in wheat. Front. Genet.12. doi: 10.3389/fgene.2021.649988

  • 19

    ShiferawB.SmaleM.BraunH. J.DuveillerE.ReynoldsM.MurichoG. (2013). Crops that feed the world 10. Past successes and future challenges to the role played by wheat in global food security. Food Secur.5, 291317. doi: 10.1007/s12571-013-0263-y

  • 20

    SukumaranS.DreisigackerS.LopesM.ChavezP.ReynoldsM. P. (2015). Genome-wide association study for grain yield and related traits in an elite spring wheat population grown in temperate irrigated environments. Theor. Appl. Genet.128, 353363. doi: 10.1007/s00122-014-2435-3

  • 21

    SunC.ZhangF.YanX.ZhangX.DongZ.CuiD.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, 953969. doi: 10.1111/pbi.12690

  • 22

    WangW.GaoX.ChengY.RenY.ZhangZ.WangR.et al. (2022). QTL mapping of leaf area index and chlorophyll content based on UAV remote sensing in wheat. Agriculture12, 595. doi: 10.3390/agriculture12050595

  • 23

    WangW.SunN.BaiB.WuH.ChengY.GengH.et al. (2024). Prediction of wheat SPAD using integrated multispectral and support vector machines. Front. Plant Sci.15. doi: 10.3389/fpls.2024.1405068

  • 24

    WangY.LiJ. (2008). Molecular basis of plant architecture. Annu. Rev. Plant Biol.59, 253279. doi: 10.1146/annurev.arplant.59.032607.092902

  • 25

    XiongH.LiY.GuoH.XieY.ZhaoL.GuJ.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

  • 26

    XuD.BianY.LuoX.JiaC.HaoQ.TianX.et al. (2023). Dissecting pleiotropic functions of the wheat green revolution gene Rht-B1b in plant morphogenesis and yield formation. Development150, dev201601. doi: 10.1242/dev.201601

  • 27

    ZhuC.GoreM.BucklerE. S.YuJ. (2008). Status and prospects of association mapping in plants. Plant Genome1, 520. doi: 10.3835/plantgenome2008.02.0089

Summary

Keywords

winter wheat, plant height, 55K SNP array, GWAS, candidate genes

Citation

Ding Y, Fan G, Gao Y, Huang T, Zhou A, Yu S, Zhao L, Shi X, Ding S, Hao J, Wang W, Song J, Sun N and Fang H (2025) Performing whole-genome association analysis of winter wheat plant height using the 55K chip. Front. Plant Sci. 15:1471636. doi: 10.3389/fpls.2024.1471636

Received

28 July 2024

Accepted

19 December 2024

Published

17 January 2025

Volume

15 - 2024

Edited by

Andrés J. Cortés, Colombian Corporation for Agricultural Research (AGROSAVIA), Colombia

Reviewed by

Jindong Liu, Chinese Academy of Agricultural Sciences, China

Felipe López-Hernández, Colombian Agricultural Research Corporation (CORPOICA), Colombia

Updates

Copyright

*Correspondence: Hui Fang,

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics