The Tetratricopeptide repeat protein TaTPR-B1 regulates spike compactness in bread wheat.

Spike compactness (SC) is strongly associated with wheat (Triticum aestivum L.) grain yield. In this study, we conducted a quantitative trait locus (QTL) analysis using a doubled haploid (DH) population derived from a cross between two common wheat varieties with contrasting spike morphology, revealing 16 stable QTLs associated with SC. The effect of a major QTL, QSc.cau-6B.1, was validated in 231 F7 recombinant inbred lines (RILs) derived from the same cross as the DH population. Using two residual heterozygous lines (RHLs), we delimited QSc.cau-6B.1 to an approximately 0.5-Mbp physical interval containing four high-confidence genes. The tetratricopeptide repeat-TraesCS6B03G1214400 (TaTPR-B1) was the priority candidate gene according to sequence and expression variations between near-isogenic lines. Accordingly, TaTPR-B1 knockout in the common wheat variety 'CB037' significantly increased SC compared to the wild type (WT). Conversely, TaTPR-B1 overexpression in the common wheat variety 'Fielder' significantly decreased SC compared to the WT. Moreover, we developed a PCR-based marker targeting the 32-bp insertion/deletion (InDel) between the two TaTPR-B1 alleles, which could be practical and valuable in modern wheat breeding programs for diagnostic purposes. Collectively, these findings provide insight into the genetic basis of SC in common wheat and present a valuable target with a breeder-friendly diagnostic marker for gene pyramid breeding.

On the evolution and genetic diversity of the bread wheat D genome.

Bread wheat (Triticum aestivum) became a globally dominant crop after incorporating the D genome from the donor species Aegilops tauschii, but the evolutionary history that shaped the D genome during this process remains to be clarified. Here, we propose a renewed evolutionary model linking Ae. tauschii and the hexaploid wheat D genome by constructing an ancestral haplotype map covering 762 Ae. tauschii and hexaploid wheat accessions. We dissected the evolutionary trajectories of Ae. tauschii lineages and reported a few independent intermediate accessions, demonstrating that low-frequency inter-sublineage gene flow had enriched the diversity of Ae. tauschii. We discovered that the D genome of hexaploid wheat was inherited from a unified ancestral template, but with a mosaic composition that was highly mixed and derived mainly from three Ae. tauschii L2 sublineages located in the Caspian coastal region. This result suggests that early agricultural activities facilitated innovations in D-genome composition and finalized the success of hexaploidization. We found that the majority (51.4%) of genetic diversity was attributed to novel mutations absent in Ae. tauschii, and we identified large Ae. tauschii introgressions from various lineages, which expanded the diversity of the wheat D genome and introduced beneficial alleles. This work sheds light on the process of wheat hexaploidization and highlights the evolutionary significance of the multi-layered genetic diversity of the bread wheat D genome.

The transcription factor CAMTA2 interacts with the histone acetyltransferase GCN5 and regulates grain weight in wheat.

Grain weight and size are major traits targeted in breeding to improve wheat (Triticum aestivum L.) yield. Here, we find that the histone acetyltransferase GENERAL CONTROL NONDEREPRESSIBLE 5 (GCN5) physically interacts with the calmodulin-binding transcription factor CAMTA2 and regulates wheat grain size and weight. gcn5 mutant grains were smaller and contained less starch. GCN5 promoted the expression of the starch biosynthesis genes SUCROSE SYNTHASE 2 (Sus2) and STARCH-BRANCHING ENZYME Ic (SBEIc) by regulating H3K9ac and H3K14ac levels in their promoters. Moreover, immunoprecipitation followed by mass spectrometry (IP-MS) revealed that CAMTA2 physically interacts with GCN5. The CAMTA2-GCN5 complex activated Sus2 and SBEIc by directly binding to their promoters and depositing H3K9ac and H3K14ac marks during wheat endosperm development. camta2 knockout mutants exhibited similar phenotypes to gcn5 mutants, including smaller grains that contained less starch. In gcn5 mutants, transcripts of high molecular weight (HMW) Glutenin (Glu) genes were downregulated, leading to reduced HMW glutenin protein levels, gluten content, and sodium dodecyl sulfate (SDS) sedimentation volume. However, the association of GCN5 with Glu genes was independent of CAMTA2, since GCN5 enrichment on Glu promoters was unchanged in camta2 knockouts. Finally, we identified a CAMTA2-AH3 elite allele that corresponded with enhanced grain size and weight, serving as a candidate gene for breeding wheat varieties with improved grain weight.

The TaWAK2-TaNAL1-TaDST pathway regulates leaf width via cytokinin signaling in wheat.

Leaves play a crucial role in photosynthesis and respiration, ultimately affecting the final grain yield of crops, including wheat (Triticum aestivum L.); however, the molecular mechanisms underlying wheat leaf development remain largely unknown. Here, we isolated a narrow-leaf gene, TaWAK2-A, through a map-based cloning strategy. TaWAK2-A encodes a wall-associated kinase (WAK), for which a single Ala-to-Val amino acid substitution reduces the protein stability, leading to a narrow-leaf phenotype in wheat. Further investigation suggests that TaWAK2 directly interacts with and phosphorylates TaNAL1, a trypsin-like serine/cysteine protease. The phosphorylated TaNAL1 is then involved in the degradation of the zinc finger transcription factor TaDST, which acts as a repressor of leaf expansion by activating the expression of the cytokinin oxidase gene TaCKX9 and triggering in vivo cytokinin degradation. Therefore, our findings elucidate a signaling cascade involving TaWAK2-TaNAL1-TaDST that sheds light on the regulation of wheat leaf development.

Genome-wide association study identifies loci and candidate genes for RVA parameters in wheat (Triticum aestivum L.).

The gelatinization and retrogradation characteristics of wheat starch affect the eating quality of Chinese-style food. Rapid Visco Analyzer (RVA) parameters have been widely used as important indicators to evaluate and improve the quality of wheat starch. However, the genetic basis of RVA parameters remains to be further explored. In the present study, a natural population was genotyped using 90K single nucleotide polymorphism (SNP) arrays, and the RVA parameters of this population grown in five environments were evaluated. The results showed that 22,068 high-quality SNP markers were identified and distributed unequally on the chromosomes. According to the genetic distance, 214 wheat materials were divided into four groups. Except for the pasting temperature (PTT), six parameters followed a normal distribution. Based on the general linear model, 969 significant association SNPs were detected by genome-wide association studies (GWAS), and chromosomes 7A and 2B had the most associated SNPs. Breakdown viscosity (BV) was associated with the most SNPs (n = 238), followed by PTT (n = 186), peak viscosity (PV; n = 156), trough viscosity (TV; n = 127), and final viscosity (FV; n = 126). According to the average linkage disequilibrium (LD), 33 stable quantitative trait loci (QTLs) were identified for single parameters in multiple environments, of which 12 were associated with BV, followed by peak time (PT; n = 8) and PTT (n = 7). On the other hand, 67 pleiotropic QTLs were identified for multiple parameters. Three candidate genes-TasbeIIa, TasbeI, and TassIIa-were screened for phenotyping analysis. The grain width and the weight of the TasbeIIa and TaSSIIa knockout (KO) lines were significantly lower than those of the TasbeI KO lines and the control (CK). The KO lines had smaller endosperm cells, smaller A-type starch granules, and higher amylose content. The TasbeI KO lines showed normal RVA curves, while the TasbeIIa KO lines showed flat curves. However, the TaSSIIa lines failed to paste under the RVA temperatures. Conclusively, the SNPs/QTLs significantly associated with the RVA parameters and genetic resources with novel haplotypes could be used to improve the quality of wheat starch.

Grain albumin content is affected by 1BL/1RS translocation and alters the processing quality of wheat (Triticum aestivum L.).

Grain albumin is highly nutritious and closely related to the processing quality of wheat. However, few studies have explored the grain albumin content (GAC) in wheat. This study aims to uncover quantitative trait loci (QTLs) linked to wheat GAC by analyzing a doubled haploid (DH) population derived from common wheat cultivars ShanNong23 and ZhouMai17. We detected six QTLs controlling GAC on chromosomes 1B, 5A, and 6D, with individual QTL explaining 5.78% to 22.29% of the GAC variation. The effect of QGac.cau-1B.1 on GAC is attributed to the presence of the 1BL/1RS translocation, indicating that the 1BL/1RS translocation increase of GAC compared with the non-1BL/1RS translocation lines. The higher GAC observed in 1BL/1RS lines could be primarily attributed to the increased accumulation of omega-secalin, omega-gliadin, low molecular weight glutenin subunit and ribosomal protein content. Additionally, we also found that the SDS-sedimentation value of whole wheat flour was decreased by adding albumin solution. These results advance our understanding of the genetic basis of GAC and offer novel perspectives for enhancing wheat quality through genetic enhancements.

Fine-mapping of LrN3B on wheat chromosome arm 3BS, one of the two complementary genes for adult-plant leaf rust resistance.

The common wheat line 4N0461 showed adult-plant resistance to leaf rust. 4N0461 was crossed with susceptible cultivars Nongda4503 and Shi4185 to map the causal resistance gene(s). Segregation of leaf rust response in F2 populations from both crosses was 9 resistant:7 susceptible, indicative of two complementary dominant resistance genes. The genes were located on chromosome arms 3BS and 4BL and temporarily named LrN3B and LrN4B, respectively. Subpopulations from 4N0461 × Nongda4503 with LrN3B segregating as a single allele were used to fine-map LrN3B locus. LrN3B was delineated in a genetic interval of 0.07 cM, corresponding to 106 kb based on the Chinese Spring reference genome (IWGSC RefSeq v1.1). Four genes were annotated in this region, among which TraesCS3B02G014800 and TraesCS3B02G014900 differed between resistant and susceptible genotypes, and both were required for LrN3B resistance in virus-induced gene silencing experiments. Diagnostic markers developed for checking the polymorphism of each candidate gene, can be used for marker-assisted selection in wheat breeding programs.

Tagging large CNV blocks in wheat boosts digitalization of germplasm resources by ultra-low-coverage sequencing.

BACKGROUND: The massive structural variations and frequent introgression highly contribute to the genetic diversity of wheat, while the huge and complex genome of polyploid wheat hinders efficient genotyping of abundant varieties towards accurate identification, management, and exploitation of germplasm resources. RESULTS: We develop a novel workflow that identifies 1240 high-quality large copy number variation blocks (CNVb) in wheat at the pan-genome level, demonstrating that CNVb can serve as an ideal DNA fingerprinting marker for discriminating massive varieties, with the accuracy validated by PCR assay. We then construct a digitalized genotyping CNVb map across 1599 global wheat accessions. Key CNVb markers are linked with trait-associated introgressions, such as the 1RS·1BL translocation and 2NvS translocation, and the beneficial alleles, such as the end-use quality allele Glu-D1d (Dx5 + Dy10) and the semi-dwarf r-e-z allele. Furthermore, we demonstrate that these tagged CNVb markers promote a stable and cost-effective strategy for evaluating wheat germplasm resources with ultra-low-coverage sequencing data, competing with SNP array for applications such as evaluating new varieties, efficient management of collections in gene banks, and describing wheat germplasm resources in a digitalized manner. We also develop a user-friendly interactive platform, WheatCNVb ( http://wheat.cau.edu.cn/WheatCNVb/ ), for exploring the CNVb profiles over ever-increasing wheat accessions, and also propose a QR-code-like representation of individual digital CNVb fingerprint. This platform also allows uploading new CNVb profiles for comparison with stored varieties. CONCLUSIONS: The CNVb-based approach provides a low-cost and high-throughput genotyping strategy for enabling digitalized wheat germplasm management and modern breeding with precise and practical decision-making.

Ploidy variation induces butterfly effect on chromatin topology in wheat.

Polyploidy is a prominent driver of plant diversification, accompanied with dramatic chromosomal rearrangement and epigenetic changes that affect gene expression. How chromatin interactions within and between subgenomes adapt to ploidy transition remains poorly understood. We generate open chromatin interaction maps for natural hexaploid wheat (AABBDD), extracted tetraploid wheat (AABB), diploid wheat progenitor Aegilops tauschii (DD) and resynthesized hexaploid wheat (RHW, AABBDD). Thousands of intra- and interchromosomal loops are de novo established or disappeared in AB subgenomes after separation of D subgenome, in which 37-95% of novel loops are lost again in RHW after merger of D genome. Interestingly, more than half of novel loops are formed by cascade reactions that are triggered by disruption of chromatin interaction between AB and D subgenomes. The interaction repressed genes in RHW relative to DD are expression suppressed, resulting in more balanced expression of the three homoeologs in RHW. The interaction levels of cascade anchors are decreased step-by-step. Leading single nucleotide polymorphisms of yield- and plant architecture-related quantitative trait locus are significantly enriched in cascade anchors. The expression of 116 genes interacted with these anchors are significantly correlated with the corresponding traits. Our findings reveal trans-regulation of intrachromosomal loops by interchromosomal interactions during genome merger and separation in polyploid species.

Molecular characterization of QTL for grain zinc and iron concentrations in wheat landrace Chinese Spring.

KEY MESSAGE: Three stable QTL for grain zinc concentration were identified in wheat landrace Chinese Spring. Favorable alleles were more frequent in landraces than in modern wheat cultivars. Wheat is a major source of dietary energy for the growing world population. Developing cultivars with enriched zinc and iron can potentially alleviate human micronutrient deficiency. In this study, a recombinant inbred line (RIL) population with 245 lines derived from cross Zhou 8425B/Chinese Spring was used to detect quantitative trait loci (QTL) for grain zinc concentration (GZnC) and grain iron concentration (GFeC) across four environments. Three stable QTL for GZnC with all favorable alleles from Chinese Spring were identified on chromosomes 3BL, 5AL, and 5BL. These QTL explaining maxima of 8.7%, 5.8%, and 7.1% of phenotypic variances were validated in 125 resequenced wheat accessions encompassing both landraces and modern cultivars using six kompetitive allele specific PCR (KASP) assays. The frequencies of favorable alleles for QGZnCzc.caas-3BL, QGZnCzc.caas-5AL and QGZnCzc.caas-5BL were higher in landraces (90.4%, 68.0%, and 100.0%, respectively) compared to modern cultivars (45.9%, 35.4%, and 40.9%), suggesting they were not selected in breeding programs. Candidate gene association studies on GZnC in the cultivar panel further delimited the QTL into 8.5 Mb, 4.1 Mb, and 47.8 Mb regions containing 46, 4, and 199 candidate genes, respectively. The 5BL QTL located in a region where recombination was suppressed. Two stable and three less stable QTL for GFeC with favorable alleles also from Chinese Spring were identified on chromosomes 4BS (Rht-B1a), 4DS (Rht-D1a), 1DS, 3AS, and 6DS. This study sheds light on the genetic basis of GZnC and GFeC in Chinese Spring and provides useful molecular markers for wheat biofortification.

A k-mer-based pangenome approach for cataloging seed-storage-protein genes in wheat to facilitate genotype-to-phenotype prediction and improvement of end-use quality.

Wheat is a staple food for more than 35% of the world's population, with wheat flour used to make hundreds of baked goods. Superior end-use quality is a major breeding target; however, improving it is especially time-consuming and expensive. Furthermore, genes encoding seed-storage proteins (SSPs) form multi-gene families and are repetitive, with gaps commonplace in several genome assemblies. To overcome these barriers and efficiently identify superior wheat SSP alleles, we developed "PanSK" (Pan-SSP k-mer) for genotype-to-phenotype prediction based on an SSP-based pangenome resource. PanSK uses 29-mer sequences that represent each SSP gene at the pangenomic level to reveal untapped diversity across landraces and modern cultivars. Genome-wide association studies with k-mers identified 23 SSP genes associated with end-use quality that represent novel targets for improvement. We evaluated the effect of rye secalin genes on end-use quality and found that removal of ω-secalins from 1BL/1RS wheat translocation lines is associated with enhanced end-use quality. Finally, using machine-learning-based prediction inspired by PanSK, we predicted the quality phenotypes with high accuracy from genotypes alone. This study provides an effective approach for genome design based on SSP genes, enabling the breeding of wheat varieties with superior processing capabilities and improved end-use quality.

Fine-mapping and validation of the major quantitative trait locus QFlANG-4B for flag leaf angle in wheat.

KEY MESSAGE: This study precisely mapped and validated a quantitative trait locus (QTL) located on chromosome 4B for flag leaf angle in wheat. Flag leaf angle (FLANG) is closely related to crop architecture and yield. We previously identified the quantitative trait locus (QTL) QFLANG-4B for FLANG on chromosome 4B, located within a 14-cM interval flanked by the markers Xbarc20 and Xzyh357, using a mapping population of recombinant inbred lines (RILs) derived from a cross between Nongda3331 (ND3331) and Zang1817. In this study, we fine-mapped QFLANG-4B and validated its associated genetic effect. We developed a BC3F3 population using ND3331 as the recurrent parent through marker-assisted selection, as well as near-isogenic lines (NILs) by selfing BC3F3 plants carrying different heterozygous segments for the QFLANG-4B region. We obtained eight recombinant types for QFLANG-4B, narrowing its location down to a 5.3-Mb region. This region contained 76 predicted genes, 7 of which we considered to be likely candidate genes for QFLANG-4B. Marker and phenotypic analyses of individual plants from the secondary mapping populations and their progeny revealed that the FLANG of the ND3331 allele is significantly higher than that of the Zang1817 allele in multiple environments. These results not only provide a basis for the map-based cloning of QFLANG-4B, but also indicate that QFLANG-4B has great potential for marker-assisted selection in wheat breeding programs designed to improve plant architecture and yield.

Positional cloning and characterization reveal the role of TaSRN-3D and TaBSR1 in the regulation of seminal root number in wheat.

Seminal roots play a critical role in water and nutrient absorption, particularly in the early developmental stages of wheat. However, the genes responsible for controlling SRN in wheat remain largely unknown. Genetic mapping and functional analyses identified a candidate gene (TraesCS3D01G137200, TaSRN-3D) encoding a Ser/Thr kinase glycogen synthase kinase 3 (STKc_GSK3) that regulated SRN in wheat. Additionally, experiments involving hormone treatment, nitrate absorption and protein interaction were conducted to explore the regulatory mechanism of TaSRN-3D. Results showed that the TaSRN-3D4332 allele inhibited seminal roots initiation and development, while loss-of-function mutants showed significantly higher seminal root number (SRN). Exogenous application of epi-brassinolide could increase the SRN in a HS2-allelic background. Furthermore, chlorate sensitivity and 15N uptake assays revealed that a higher number of seminal roots promoted nitrate accumulation. TaBSR1 (BIN2-related SRN Regulator 1, orthologous to OsGRF4/GL2 in rice) acts as an interactor of TaSRN-3D and promotes TaBSR1 degradation to reduce SRN. This study provides valuable insights into understanding the genetic basis and regulatory network of SRN in wheat, highlighting their roles as potential targets for root-based improvement in wheat breeding.

The potassium transporter TaNHX2 interacts with TaGAD1 to promote drought tolerance via modulating stomatal aperture in wheat.

Drought is a major global challenge in agriculture that decreases crop production. γ-Aminobutyric acid (GABA) interfaces with drought stress in plants; however, a mechanistic understanding of the interaction between GABA accumulation and drought response remains to be established. Here we showed the potassium/proton exchanger TaNHX2 functions as a positive regulator in drought resistance in wheat by mediating cross-talk between the stomatal aperture and GABA accumulation. TaNHX2 interacted with glutamate decarboxylase TaGAD1, a key enzyme that synthesizes GABA from glutamate. Furthermore, TaNHX2 targeted the C-terminal auto-inhibitory domain of TaGAD1, enhanced its activity, and promoted GABA accumulation under drought stress. Consistent with this, the tanhx2 and tagad1 mutants showed reduced drought tolerance, and transgenic wheat with enhanced TaNHX2 expression had a yield advantage under water deficit without growth penalty. These results shed light on the plant stomatal movement mechanism under drought stress and the TaNHX2-TaGAD1 module may be harnessed for amelioration of negative environmental effects in wheat as well as other crops.

Innovative computational tools provide new insights into the polyploid wheat genome.

Bread wheat (Triticum aestivum) is an important crop and serves as a significant source of protein and calories for humans, worldwide. Nevertheless, its large and allopolyploid genome poses constraints on genetic improvement. The complex reticulate evolutionary history and the intricacy of genomic resources make the deciphering of the functional genome considerably more challenging. Recently, we have developed a comprehensive list of versatile computational tools with the integration of statistical models for dissecting the polyploid wheat genome. Here, we summarize the methodological innovations and applications of these tools and databases. A series of step-by-step examples illustrates how these tools can be utilized for dissecting wheat germplasm resources and unveiling functional genes associated with important agronomic traits. Furthermore, we outline future perspectives on new advanced tools and databases, taking into consideration the unique features of bread wheat, to accelerate genomic-assisted wheat breeding.

Genomic insights into the origin and evolution of spelt (Triticum spelta L.) as a valuable gene pool for modern wheat breeding.

Spelt (Triticum aestivum ssp. spelta) is an important wheat subspecies mainly cultivated in Europe before the 20th century that has contributed to modern wheat breeding as a valuable genetic resource. However, relatively little is known about the origins and maintenance of spelt populations. Here, using resequencing data from 416 worldwide wheat accessions, including representative spelt wheat, we demonstrate that European spelt emerged when primitive hexaploid wheat spread to the west and hybridized with pre-settled domesticated emmer, the putative maternal donor. Genomic introgression regions from domesticated emmer confer spelt's primitive morphological characters used for species taxonomy, such as tenacious glumes and later flowering. We propose a haplotype-based "spelt index" to identify spelt-type wheat varieties and to quantify utilization of the spelt gene pool in modern wheat cultivars. This study reveals the genetic basis for the establishment of the spelt wheat subspecies in a specific ecological niche and the vital role of the spelt gene pool as a unique germplasm resource in modern wheat breeding.

TaNAM-6A is essential for nitrogen remobilisation and regulates grain protein content in wheat (Triticum aestivum L.).

Grain protein content (GPC) is a crucial quality trait in bread wheat, which is influenced by the key transcription factor TaNAM. However, the regulatory mechanisms of TaNAM have remained largely elusive. In this study, a new role of TaNAM was unveiled in regulating nitrogen remobilisation which impacts GPC. The TaNAM knockout mutants generated by clustered regularly interspaced short palindromic repeats/Cas9 exhibited significantly delayed senescence and lower GPC, while overexpression of TaNAM-6A resulted in premature senility and much higher GPC. Further analysis revealed that TaNAM directly activates the genes TaNRT1.1 and TaNPF5.5s, which are involved in nitrogen remobilisation. This activity aids in the transfer of nitrogen from leaves to grains for protein synthesis. In addition, an elite allele of TaNAM-6A, associated with high GPC, was identified as a candidate gene for breeding high-quality wheat. Overall, our work not only elucidates the potential mechanism of TaNAM-6A affecting bread wheat GPC, but also highlights the significance of nitrogen remobilisation from senescent leaves to grains for protein accumulation. Moreover, our research provides a new target and approach for improving the quality traits of wheat, particularly the GPC.

Methyltransferase TaSAMT1 mediates wheat freezing tolerance by integrating brassinosteroid and salicylic acid signaling.

Cold injury is a major environmental stress affecting the growth and yield of crops. Brassinosteroids (BRs) and salicylic acid (SA) play important roles in plant cold tolerance. However, whether or how BR signaling interacts with the SA signaling pathway in response to cold stress is still unknown. Here, we identified an SA methyltransferase, TaSAMT1 that converts SA to methyl SA (MeSA) and confers freezing tolerance in wheat (Triticum aestivum). TaSAMT1 overexpression greatly enhanced wheat freezing tolerance, with plants accumulating more MeSA and less SA, whereas Tasamt1 knockout lines were sensitive to freezing stress and accumulated less MeSA and more SA. Spraying plants with MeSA conferred freezing tolerance to Tasamt1 mutants, but SA did not. We revealed that BRASSINAZOLE-RESISTANT 1 (TaBZR1) directly binds to the TaSAMT1 promoter and induces its transcription. Moreover, TaBZR1 interacts with the histone acetyltransferase TaHAG1, which potentiates TaSAMT1 expression via increased histone acetylation and modulates the SA pathway during freezing stress. Additionally, overexpression of TaBZR1 or TaHAG1 altered TaSAMT1 expression and improved freezing tolerance. Our results demonstrate a key regulatory node that connects the BR and SA pathways in the plant cold stress response. The regulatory factors or genes identified could be effective targets for the genetic improvement of freezing tolerance in crops.

Natural variation of STKc_GSK3 kinase TaSG-D1 contributes to heat stress tolerance in Indian dwarf wheat.

Heat stress threatens global wheat (Triticum aestivum) production, causing dramatic yield losses worldwide. Identifying heat tolerance genes and comprehending molecular mechanisms are essential. Here, we identify a heat tolerance gene, TaSG-D1E286K, in Indian dwarf wheat (Triticum sphaerococcum), which encodes an STKc_GSK3 kinase. TaSG-D1E286K improves heat tolerance compared to TaSG-D1 by enhancing phosphorylation and stability of downstream target TaPIF4 under heat stress condition. Additionally, we reveal evolutionary footprints of TaPIF4 during wheat selective breeding in China, that is, InDels predominantly occur in the TaPIF4 promoter of Chinese modern wheat cultivars and result in decreased expression level of TaPIF4 in response to heat stress. These sequence variations with negative effect on heat tolerance are mainly introduced from European germplasm. Our study provides insight into heat stress response mechanisms and proposes a potential strategy to improve wheat heat tolerance in future.