Wheat heat tolerance is impaired by heightened deletions in the distal end of 4AL chromosomal arm

A Bread Wheat Line with the Substituted Wild Emmer Chromosome 4A Results in Fragment Deletions of Chromosome 4B and Weak Plants

In response to the growing genetic uniformity within wheat populations, developing efficient wheat-alien translocation strategies has become critically important. We observed that several offspring of the common wheat (Triticum aestivum L.)-wild emmer (Triticum turgidum L. var. dicoccoides) chromosome arm substitution line (CASL4AL) exhibited stunted growth, including significantly reduced plant height, spike length, spikelet number, and stem width compared to normal plants. Integrative transcriptomic analyses (RNA-Seq and BSR-Seq) revealed a statistically significant depletion (p < 0.01) of single nucleotide polymorphisms (SNPs) on chromosome 4B in compromised plants. Chromosome association analysis of differentially expressed genes (DEGs, up- or downregulated) revealed that downregulated genes were predominantly located on chromosome 4B. The 1244 downregulated DEGs on Chr4B were employed for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, and RNA metabolic processes, DNA repair, and transport systems were significantly enriched by GO analysis; however, only the mRNA surveillance pathway was enriched by KEGG enrichment. Molecular marker profiling showed a complete absence of target amplification in the critical 0-155 Mb region of chromosome 4B in all weak plants. Pearson's correlation coefficients confirmed significant associations (p < 0.01) between 4B-specific amplification and weak phenotypes. These results demonstrate that 4B segmental deletions drive weak phenotypes in CASL4AL progeny, and provide experimental evidence for chromosome deletions induced in wild emmer chromosome substitution lines. This study highlights the potential of wild emmer as a valuable tool for generating chromosomal variations in wheat breeding programs.

Deciphering the Vulnerability of Pollen to Heat Stress for Securing Crop Yields in a Warming Climate

Climate change is leading to more frequent and severe extreme temperature events, negatively impacting agricultural productivity and threatening global food security. Plant reproduction, the process fundamental to crop yield, is highly susceptible to heatwaves, which disrupt pollen development and ultimately affect seed-set and crop yields. Recent research has increasingly focused on understanding how pollen grains from various crops react to heat stress at the molecular and cellular levels. This surge in interest over the last decade has been driven by advances in genomic technologies, such as single-cell RNA sequencing, which holds significant potential for revealing the underlying regulatory reprogramming triggered by heat stress throughout the various stages of pollen development. This review focuses on how heat stress affects gene regulatory networks, including the heat stress response, the unfolded protein response, and autophagy, and discusses the impact of these changes on various stages of pollen development. It highlights the potential of pollen selection as a key strategy for improving heat tolerance in crops by leveraging the genetic variability among pollen grains. Additionally, genome-wide association studies and population screenings have shed light on the genetic underpinnings of traits in major crops that respond to high temperatures during male reproductive stages. Gene-editing tools like CRISPR/Cas systems could facilitate precise genetic modifications to boost pollen heat resilience. The information covered in this review is valuable for selecting traits and employing molecular genetic approaches to develop heat-tolerant genotypes.

Fine-mapping of PmHHM, a broad-spectrum allele from a wheat landrace conferring both seedling and adult resistance to powdery mildew

Introduction

Common wheat is a leading global food crop that impacts food security. Wheat powdery mildew (PM), caused by Blumeria graminis f. sp. tritici (Bgt), poses a significant threat to grain yield and flour quality. The identification and utilization of broad-spectrum resistance genes against PM are essential for effective disease control.

Methods

The resistance spectrum test during the seedling stage and the identification of resistance during the adult stage were conducted to evaluate the wheat landrace Honghuamai (HHM). Five segregating populations were investigated to assess the inheritance of PM resistance in HHM. To map its PM resitance gene, bulked segregant analysis, molecular mapping and comparative genomic analysis were also used in the present study.

Results

HHM shows remarkable adult resistance in the field and is nearly immune to all 25 Bgt isolates used in seedling tests, making it an excellent source of PM resistance. PM resistance in HHM was determined by a single dominant gene, temporarily named PmHHM. It was then fine-mapped to an interval with a genetic distance of 0.0031 cM and a physical distance of 187.4 kb on chromosome 4AL of the Chinese Spring reference sequence v.2.1. Four genes were identified in the target region, three of which encode nucleotide-binding leucine-rich repeat (NLR) proteins. Comparative genomic analysis revealed presence/absence variations (PAVs) of the PmHHM locus among common wheat varieties.

Discussion

These closely linked molecular markers will not only benefit the cloning of the gene underlying PmHHM but also facilitate the efficient utilization of the gene in breeding programs.

Wheat2035: Integrating pan-omics and advanced biotechnology for future wheat design

Wheat (Triticum aestivum) production is vital for global food security, providing energy and protein to millions of people worldwide. Recent advancements in wheat research have led to significant increases in production, fueled by technological and scientific innovation. Here, we summarize the major advancements in wheat research, particularly the integration of biotechnologies and a deeper understanding of wheat biology. The shift from multi-omics to pan-omics approaches in wheat research has greatly enhanced our understanding of the complex genome, genomic variations, and regulatory networks to decode complex traits. We also outline key scientific questions, potential research directions, and technological strategies for improving wheat over the next decade. Since global wheat production is expected to increase by 60% in 2050, continued innovation and collaboration are crucial. Integrating biotechnologies and a deeper understanding of wheat biology will be essential for addressing future challenges in wheat production, ensuring sustainable practices and improved productivity.

Heat stress tolerance in wheat seedling: Clustering genotypes and identifying key traits using multivariate analysis

Elevated atmospheric heat is considered as one of the bottlenecks for global wheat production. Screening potential wheat genotypes against heat stress and selecting some suitable indicators to assist in understanding thermotolerance could be crucial for sustaining wheat cultivation. Accordingly, 80 diverse bread wheat genotypes were evaluated in controlled lab condition by imposing a week-long heat stress (35/25 °C D/N) at the seedling stage. The response of heat stress was evaluated using multivariate analysis techniques on 20 morpho-physiological traits. Results showed significant variations in the studied traits due to the imposition of heat stress. Eleven seedling traits that contributed significantly to the genotypic variability were identified using principal component analysis (PCA). A substantial correlation between most of the selected seedling attributes was observed. Hierarchical cluster analysis identified three distinct clusters among the tested wheat genotypes. Cluster 1, consisting of 33 genotypes, exhibited the highest tolerance to heat stress, followed by Cluster 2 (18 genotypes) with moderate tolerance and Cluster 3 (29 genotypes) showing susceptibility. Linear discriminant analysis (LDA) approved that nearly 93 % of the wheat genotypes were appropriately ascribed to each cluster. The squared distance analysis confirmed the distinct nature of the clusters. Using multi-trait genotype-ideotype distance index (MGIDI), all 12 identified tolerant genotypes (BG-30, BD-468, BG-24, BD-9908, BG-32, BD-476, BD-594, BD-553, BD-488, BG-33, BD-495, and AS-10627) originated from Cluster 1. Selection gain in MGIDI analysis, broad-sense heritability, and multiple linear regression analysis together identified shoot and root dry and fresh weights, chlorophyll contents (a and total), shoot tissue water content, root-shoot dry weight ratio, and efficiency of photosystem II (PS II) as the most vital discriminatory factors explaining heat stress tolerance of 80 wheat genotypes. The identified genotypes with superior thermotolerance would offer resourceful genetic tools for breeders to improve wheat yield in warmer regions. The traits found to have greater contribution in explaining heat stress tolerance will be equally important in prioritizing future research endeavors.

Genome-wide association study for seedling heat tolerance under two temperature conditions in bread wheat (Triticum aestivum L.)

BACKGROUND

As the greenhouse effect intensifies, global temperatures are steadily increasing, posing a challenge to bread wheat (Triticum aestivum L.) production. It is imperative to comprehend the mechanism of high temperature tolerance in wheat and implement breeding programs to identify and develop heat-tolerant wheat germplasm and cultivars.

RESULTS

To identify quantitative trait loci (QTL) related to heat stress tolerance (HST) at seedling stage in wheat, a panel of 253 wheat accessions which were re-sequenced used to conduct genome-wide association studies (GWAS) using the factored spectrally transformed linear mixed models (FaST-LMM). For most accessions, the growth of seedlings was found to be inhibited under heat stress. Analysis of the phenotypic data revealed that under heat stress conditions, the main root length, total root length, and shoot length of seedlings decreased by 47.46%, 49.29%, and 15.19%, respectively, compared to those in normal conditions. However, 17 varieties were identified as heat stress tolerant germplasm. Through GWAS analysis, a total of 115 QTLs were detected under both heat stress and normal conditions. Furthermore, 15 stable QTL-clusters associated with heat response were identified. By combining gene expression, haplotype analysis, and gene annotation information within the physical intervals of the 15 QTL-clusters, two novel candidate genes, TraesCS4B03G0152700/TaWRKY74-B and TraesCS4B03G0501400/TaSnRK3.15-B, were responsive to temperature and identified as potential regulators of HST in wheat at the seedling stage.

CONCLUSIONS

This study conducted a detailed genetic analysis and successfully identified two genes potentially associated with HST in wheat at the seedling stage, laying a foundation to further dissect the regulatory mechanism underlying HST in wheat under high temperature conditions. Our finding could serve as genomic landmarks for wheat breeding aimed at improving adaptation to heat stress in the face of climate change.

Fine mapping of powdery mildew resistance gene PmXNM in a Chinese wheat landrace Xiaonanmai

KEY MESSAGE

Powdery mildew resistance gene PmXNM, originated from the Chinese wheat landrace Xiaonanmai, was delimited to a 300.7-kb interval enriched with resistance genes. Powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is a globally devastating disease threatening the yield and quality of wheat worldwide. The use of broad-spectrum disease resistance genes from wheat landraces is an effective strategy to prevent this pathogen. Chinese wheat landrace Xiaonanmai (XNM) was immune to 23 tested Bgt isolates at the seedling stage. The F1, F2, and F2:4 progenies derived from the cross between XNM and Chinese Spring (CS) were used in this study. Genetic analysis revealed that powdery mildew resistance in XNM was controlled by a single dominant gene, temporarily designated PmXNM. Bulked segregant analysis and molecular mapping delimited PmXNM to the distal terminal region of chromosome 4AL flanked by markers caps213923 and kasp511718. The region carrying the PmXNM locus was approximately 300.7 kb and contained nine high-confidence genes according to the reference genome sequence of CS. Five of these genes, annotated as disease resistance RPP13-like proteins 1, were clustered in the target region. Haplotype analysis using the candidate gene-specific markers indicated that the majority of 267 common wheat accessions (75.3%) exhibited extensive gene losses at the PmXNM locus, as confirmed by aligning the targeted genome sequences of CS with those of other sequenced wheat cultivars. Seven candidate gene-specific markers have proven effective for marker-assisted introgression of PmXNM into modern elite cultivars.

Genome resources for the elite bread wheat cultivar Aikang 58 and mining of elite homeologous haplotypes for accelerating wheat improvement

Despite recent progress in crop genomics studies, the genomic changes brought about by modern breeding selection are still poorly understood, thus hampering genomics-assisted breeding, especially in polyploid crops with compound genomes such as common wheat (Triticum aestivum). In this work, we constructed genome resources for the modern elite common wheat variety Aikang 58 (AK58). Comparative genomics between AK58 and the landrace cultivar Chinese Spring (CS) shed light on genomic changes that occurred through recent varietal improvement. We also explored subgenome diploidization and divergence in common wheat and developed a homoeologous locus-based genome-wide association study (HGWAS) approach, which was more effective than single homoeolog-based GWAS in unraveling agronomic trait-associated loci. A total of 123 major HGWAS loci were detected using a genetic population derived from AK58 and CS. Elite homoeologous haplotypes (HHs), formed by combinations of subgenomic homoeologs of the associated loci, were found in both parents and progeny, and many could substantially improve wheat yield and related traits. We built a website where users can download genome assembly sequence and annotation data for AK58, perform blast analysis, and run JBrowse. Our work enriches genome resources for wheat, provides new insights into genomic changes during modern wheat improvement, and suggests that efficient mining of elite HHs can make a substantial contribution to genomics-assisted breeding in common wheat and other polyploid crops.

Identification of Single-Nucleotide Polymorphisms (SNPs) Associated with Heat Tolerance at the Reproductive Stage in Synthetic Hexaploid Wheats Using GWAS

The projected rise in global ambient temperature by 3-5 °C by the end of this century, along with unpredicted heat waves during critical crop growth stages, can drastically reduce grain yield and will pose a great food security challenge. It is therefore important to identify wheat genetic resources able to withstand high temperatures, discover genes underpinning resilience to higher temperatures, and deploy such genetic resources in wheat breeding to develop heat-tolerant cultivars. In this study, 180 accessions of synthetic hexaploid wheats (SHWs) were evaluated under normal and late wheat growing seasons (to expose them to higher temperatures) at three locations (Islamabad, Bahawalpur, and Tando Jam), and data were collected on 11 morphological and yield-related traits. The diversity panel was genotyped with a 50 K SNP array to conduct genome-wide association studies (GWASs) for heat tolerance in SHW. A known heat-tolerance locus, TaHST1, was profiled to identify different haplotypes of this locus in SHWs and their association with grain yield and related traits in SHWs. There was a 36% decrease in grain yield (GY), a 23% decrease in thousand-grain weight (TKW), and an 18% decrease in grains per spike (GpS) across three locations in the population due to the heat stress conditions. GWASs identified 143 quantitative trait nucleotides (QTNs) distributed over all 21 chromosomes in the SHWs. Out of these, 52 QTNs were associated with morphological and yield-related traits under heat stress, while 15 of them were pleiotropically associated with multiple traits. The heat shock protein (HSP) framework of the wheat genome was then aligned with the QTNs identified in this study. Seventeen QTNs were in proximity to HSPs on chr2B, chr3D, chr5A, chr5B, chr6D, and chr7D. It is likely that QTNs on the D genome and those in proximity to HSPs may carry novel alleles for heat-tolerance genes. The analysis of TaHST1 indicated that 15 haplotypes were present in the SHWs for this locus, while hap1 showed the highest frequency of 25% (33 SHWs). These haplotypes were significantly associated with yield-related traits in the SHWs. New alleles associated with yield-related traits in SHWs could be an excellent reservoir for breeding deployment.

Wheat EARLY FLOWERING 3 affects heading date without disrupting circadian oscillations

Plant breeders have indirectly selected for variation at circadian-associated loci in many of the world's major crops, when breeding to increase yield and improve crop performance. Using an eight-parent Multiparent Advanced Generation Inter-Cross (MAGIC) population, we investigated how variation in circadian clock-associated genes contributes to the regulation of heading date in UK and European winter wheat (Triticum aestivum) varieties. We identified homoeologues of EARLY FLOWERING 3 (ELF3) as candidates for the Earliness per se (Eps) D1 and B1 loci under field conditions. We then confirmed a single-nucleotide polymorphism within the coding region of TaELF3-B1 as a candidate polymorphism underlying the Eps-B1 locus. We found that a reported deletion at the Eps-D1 locus encompassing TaELF3-D1 is, instead, an allele that lies within an introgression region containing an inversion relative to the Chinese Spring D genome. Using Triticum turgidum cv. Kronos carrying loss-of-function alleles of TtELF3, we showed that ELF3 regulates heading, with loss of a single ELF3 homoeologue sufficient to alter heading date. These studies demonstrated that ELF3 forms part of the circadian oscillator; however, the loss of all homoeologues was required to affect circadian rhythms. Similarly, loss of functional LUX ARRHYTHMO (LUX) in T. aestivum, an orthologue of a protein partner of Arabidopsis (Arabidopsis thaliana) ELF3, severely disrupted circadian rhythms. ELF3 and LUX transcripts are not co-expressed at dusk, suggesting that the structure of the wheat circadian oscillator might differ from that of Arabidopsis. Our demonstration that alterations to ELF3 homoeologues can affect heading date separately from effects on the circadian oscillator suggests a role for ELF3 in cereal photoperiodic responses that could be selected for without pleiotropic deleterious alterations to circadian rhythms.

Transcriptome Analysis Reveals Potential Mechanism in Storage Protein Trafficking within Developing Grains of Common Wheat

Gluten proteins are the major storage protein fraction in the mature wheat grain. They are restricted to the starchy endosperm, which defines the viscoelastic properties of wheat dough. The synthesis of these storage proteins is controlled by the endoplasmic reticulum (ER) and is directed into the vacuole via the Golgi apparatus. In the present study, transcriptome analysis was used to explore the potential mechanism within critical stages of grain development of wheat cultivar "Shaannong 33" and its sister line used as the control (CK). Samples were collected at 10 DPA (days after anthesis), 14 DPA, 20 DPA, and 30 DPA for transcriptomic analysis. The comparative transcriptome analysis identified that a total of 18,875 genes were differentially expressed genes (DEGs) between grains of four groups "T10 vs. CK10, T14 vs. CK14, T20 vs. CK20, and T30 vs. CK30", including 2824 up-regulated and 5423 down-regulated genes in T30 vs. CK30. Further, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment highlighted the maximum number of genes regulating protein processing in the endoplasmic reticulum (ER) during grain enlargement stages (10-20 DPA). In addition, KEGG database analysis reported 1362 and 788 DEGs involved in translation, ribosomal structure, biogenesis, flavonoid biosynthesis pathway and intracellular trafficking, secretion, and vesicular transport through protein processing within ER pathway (ko04141). Notably, consistent with the higher expression of intercellular storage protein trafficking genes at the initial 10 DPA, there was relatively low expression at later stages. Expression levels of nine randomly selected genes were verified by qRT-PCR, which were consistent with the transcriptome data. These data suggested that the initial stages of "cell division" played a significant role in protein quality control within the ER, thus maintaining the protein quality characteristics at grain maturity. Furthermore, our data suggested that the protein synthesis, folding, and trafficking pathways directed by a different number of genes during the grain enlargement stage contributed to the observed high-quality characteristics of gluten protein in Shaannong 33 (Triticum aestivum L.).

Heat Stress Tolerance 2 confers basal heat stress tolerance in allohexaploid wheat (Triticum aestivum L.)

Heat stress substantially reduces the yield potential of wheat (Triticum aestivum L.), one of the most widely cultivated staple crops, and greatly threatens global food security in the context of global warming. However, few studies have explored the heat stress tolerance (HST)-related genetic resources in wheat. Here, we identified and fine-mapped a wheat HST locus, TaHST2, which is indispensable for HST in both the vegetative and reproductive stages of the wheat life cycle. The studied pair of near isogenic lines (NILs) exhibited diverse morphologies under heat stress, based on which we mapped TaHST2 to a 485 kb interval on chromosome arm 4DS. Under heat stress, TaHST2 confers a superior conversion rate from soluble sugars to starch in wheat grains, resulting in faster grain filling and a higher yield potential. A further exploration of genetic resources indicated that TaHST2 underwent strong artificial selection during wheat domestication, suggesting it is an essential locus for basal HST in wheat. Our findings provide deeper insights into the genetic basis of wheat HST and might be useful for global efforts to breed heat-stress-tolerant cultivars.

Fine mapping of the Hairy glume (Hg) gene in a chromosome variation region at the distal terminus of 1AS

Trichomes are differentiated epidermal cells and exist on above-ground organs of nearly all land plants with important roles in resistance to a wide range of biotic and abiotic stresses. We attempted to obtain candidate gene (s) for Hairy glume (Hg), responsible for the trichome on wheat glume, by using bulked segregant exome capture sequencing (BSE-Seq), while Hg was only mapped in 0.52-3.26 Mb of 1AS. To further fine map this gene and identify candidate genes in this region, a near isogenic line-derived population consisting of 2,050 F2 lines was generated in the present study. By analyzing this population, Hg was fine mapped into a 0.90 cM region covering a physical distance of ~825.03 Kb encompassing 6 high- and 23 low-confidence genes in the reference genome of Chinese Spring. A presence-absence variation was identified in the fine mapping region through analyses of sequence-tagged sites markers and genome sequences of the hairy glume parent of the near isogenic lines. The results presented here will be useful for further cloning Hg in wheat.

Fine mapping of Pm58 from Aegilops tauschii conferring powdery mildew resistance

The powdery mildew resistance gene Pm58 was traced to a 141.3-kb interval with the co-segregating marker Xkasp68500 in wheat breeding. Pm58 is a powdery mildew resistance gene identified in Aegilops tauschii accession TA1662 and effective in a common wheat background. To finely map Pm58, an F₂ population of 676 plants derived from the cross T093 × TA1662 was used for recombinant screening. We obtained 13 recombinants that occurred between the flanking markers Xhnu670 and Xhnu186. Genotyping and phenotyping these recombinant F_2:3 families delimited Pm58 to a 0.22-cM interval (Xsts20220-Xkasp61553) on chromosome arm 2DS. The region carrying the Pm58 locus was approximately 141.3-kb, which contained eight annotated genes according to the reference genome sequence of Ae. tauschii AL8/78. Haplotype analysis of 178 Ae. tauschii accessions using the candidate gene-specific markers identified a disease resistance gene AET2Gv20068500 as a candidate for Pm58. Comparative mapping of the Pm58-containing interval revealed two presence/absence variations (PAVs) between AL8/78 and common wheat Chinese Spring. PAV-1 resides in the 3'-end of AET2Gv20068500. The majority of 158 common wheat cultivars (84.8%) displayed the absence of a 14.1-kb fragment in the PAV-1 region, which was confirmed by aligning the targeted genome sequences of the other sequenced Ae. tauschii accessions and common wheat cultivars. A co-segregating marker Xkasp68500 developed from AET2Gv20068500 can distinguish TA1662 from all randomly selected common wheat cultivars and will be instrumental for tracking Pm58 in breeding programs.

Impacts, Tolerance, Adaptation, and Mitigation of Heat Stress on Wheat under Changing Climates

Heat stress (HS) is one of the major abiotic stresses affecting the production and quality of wheat. Rising temperatures are particularly threatening to wheat production. A detailed overview of morpho-physio-biochemical responses of wheat to HS is critical to identify various tolerance mechanisms and their use in identifying strategies to safeguard wheat production under changing climates. The development of thermotolerant wheat cultivars using conventional or molecular breeding and transgenic approaches is promising. Over the last decade, different omics approaches have revolutionized the way plant breeders and biotechnologists investigate underlying stress tolerance mechanisms and cellular homeostasis. Therefore, developing genomics, transcriptomics, proteomics, and metabolomics data sets and a deeper understanding of HS tolerance mechanisms of different wheat cultivars are needed. The most reliable method to improve plant resilience to HS must include agronomic management strategies, such as the adoption of climate-smart cultivation practices and use of osmoprotectants and cultured soil microbes. However, looking at the complex nature of HS, the adoption of a holistic approach integrating outcomes of breeding, physiological, agronomical, and biotechnological options is required. Our review aims to provide insights concerning morpho-physiological and molecular impacts, tolerance mechanisms, and adaptation strategies of HS in wheat. This review will help scientific communities in the identification, development, and promotion of thermotolerant wheat cultivars and management strategies to minimize negative impacts of HS.

The genetic and molecular basis for improving heat stress tolerance in wheat

Wheat production requires at least ~ 2.4% increase per year rate by 2050 globally to meet food demands. However, heat stress results in serious yield loss of wheat worldwide. Correspondingly, wheat has evolved genetic basis and molecular mechanisms to protect themselves from heat-induced damage. Thus, it is very urgent to understand the underlying genetic basis and molecular mechanisms responsive to elevated temperatures to provide important strategies for heat-tolerant varieties breeding. In this review, we focused on the impact of heat stress on morphology variation at adult stage in wheat breeding programs. We also summarize the recent studies of genetic and molecular factors regulating heat tolerance, including identification of heat stress tolerance related QTLs/genes, and the regulation pathway in response to heat stress. In addition, we discuss the potential ways to improve heat tolerance by developing new technologies such as genome editing. This review of wheat responses to heat stress may shed light on the understanding heat-responsive mechanisms, although the regulatory network of heat tolerance is still ambiguous in wheat.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42994-021-00064-z.

Genome-wide association for heat tolerance at seedling stage in historical spring wheat cultivars

Increasing global temperature has adverse effects on crop health and productivity at both seedling and reproductivity stages. It is paramount to develop heat tolerant wheat cultivars able to sustain under high and fluctuating temperature conditions. An experiment was conducted to characterize 194 historical wheat cultivars of Pakistan under high temperature at seedling stage to identify loci associated with heat tolerance using genome-wide association studies (GWAS). A quantitative trait locus, TaHST1, on chr4A was also characterized to identify the haplotypes at this locus associated with heat tolerance in wheat from Pakistan. Initially, the diversity panel was planted under control conditions (25°C/20°C day and night temperature) in a glass house. At three leaf stage, plants were subjected to heat stress (HS) by increasing temperature (40°C/35°C day and night), while one treatment was kept at control condition. After 7 days of HS, data were collected for seedling morphology. Heat stress reduced these traits by 25% (root weight) to 40% (shoot weight), and shoot biomass was largely affected by heat stress. A GWAS model, fixed and random model circulating probability unification (FarmCPU), identified 43 quantitative trait nucleotides (QTNs) on all chromosomes, except chr7B, were associated under both HS and control conditions. Thirteen QTNs were identified in control, while 30 QTNs were identified in HS condition. In total, 24 haplotypes were identified at TaHST1 locus, and most of the heat tolerant genotypes were assigned to Hap-20 and Hap-21. Eleven QTNs were identified within 0.3-3.1 Mb proximity of heat shock protein (HSP). Conclusively, this study provided a detailed genetic framework of heat tolerance in wheat at the seedling stage and identify potential genetic regions associated with heat tolerance which can be used for marker assisted selection (MAS) in breeding for heat stress tolerance.