Exploiting Agronomic and Biochemical Traits to Develop Heat Resilient Cotton Cultivars under Climate Change Scenarios

: The development of high-yielding heat-tolerant cotton cultivars harboring plastic phenotypes across warming climatic regions is prime objectives of today’s cotton breeding programs. We evaluated eight parents and 15 F 1 hybrids under normal and heat stress conditions. Agronomic and biochemical characters were analyzed using standard least square, correlation, principal component analysis (PCA), and hierarchical clustering. The results explained a signiﬁcant reduction in all traits except hydrogen peroxide contents, catalase, and peroxidase activities with a prominent increase under heat stress. A signiﬁcant positive correlation was observed among all agronomic and biochemical traits. POD was found to have a maximum positive correlation with CAT (0.947) and minimum with boll weight (0.050). PCA showed ﬁrst two components accounting for 78.64% of the total variation, with 55.83% and 22.80% of the total variation, respectively. Based on multivariate analyses methods 23 genotypes have been placed in 3 groups: tolerant (cluster-3), moderately tolerant (cluster-2), and susceptible (cluster-1). In a general perspective hybrids have better performance across normal and heat stress supports the idea of hybrid adaptability across stress environments. In speciﬁc FH-458 × FH-313 cross performed best across both conditions for yield and physiological traits. Hence, the generated information from the present study would support breeders in developing heat-resilient cultivars to endure the prevailing extreme environmental conditions.


Introduction
Cotton (Gossypium hirsutum L.) is a major source of oil and natural fiber [1]. Heat stress is the primary restraining factor affecting cotton growth and production [2]. Cotton is mainly cultivated across tropical and sub-tropical regions with temperatures ranging from 40 to 45 • C. Heat stress resulting from elevated temperature is considered a severe agricultural issue across many areas worldwide. The occurrence of elevated temperature causes an array of biological changes within plants, ultimately affecting plant growth and development, leading to a severe decline in productivity [3]. It is assumed that plant development and biomass accumulation are highly temperature-dependent across the growing season [4,5]. Plant response to heat stress mainly considers the peak temperature environments, including correlation and regression analyses and factor and cluster analyses to assess yield contributing characters for breeding programs [35].
Phenotypic flexibility that implies to traits is not constant across changing climatic conditions, especially rapidly warming. In recent years, it common, as one genotype may exhibit different phenotypes due to the prevailing environmental conditions [30,36]. In agreement with this phenomenon, assessment of the various morphological, biochemical traits may play a major role in maintaining germplasm resources, especially in terms of capturing resources as well as its conversion towards sustainable economic yield. It can provide a strong basis for the possibility and successful future crop improvement programs [37]. It can be achieved through analyzing these traits for the assessment of existing variability between and within populations through using mathematical models including linear regression models on the one hand whereas various multivariate techniques viz. correlation, principal component analysis (PCA) and cluster analyses on the other hand. These statistical techniques have been demonstrated as the most suitable in assessing variations and can efficiently discriminate multiple individual genotypes as well as different populations even though the measured variables may or may not be considered in tandem [35,[38][39][40].
To this end, the current study has been designed to identify suites of phenotypic and multivariate characteristics of available cotton genetic resources that may assist in developing new cultivars and contribute to the theoretical determination of the agronomical value of these genetic assets. A set of advance statistical procedures including least square fit and multivariate analysis techniques including correlation, principal component analysis (PCA) and cluster analyses have been utilized for the discrimination and grouping of heat tolerant and susceptible accessions for further selection. This kind of base information undoubtedly can outline the possibility of contribution in the existing pool of natural variability in economically important traits in cotton, which might be used to develop highyielding cultivars suitable for higher temperature areas and come across with significant evolutionary benefits crop can owe.

Materials and Methods
Eight cotton genotypes have been selected from a pre-screening experiment of cotton germplasm including 50 accessions for heat tolerance and already designated as tolerant and susceptible genotypes. Among them five heat-tolerant cotton genotypes (lines) and three heat-sensitive genotypes (testers) were crossed in line × tester mating design in 2018 at the experimental farm, Department of Plant Breeding and Genetics, University of Agriculture Faisalabad, Pakistan. During the next cotton growing season in June 2019, F 1 seeds of 15 crosses and their parents were sown in field conditions under normal and heat stress conditions with three replications under randomized complete block design (RCBD) following by split-plot arrangement. All agronomic practices were kept constant throughout the season. The list of genotypes is given in (Table S1).
At 50% flowering, heat stress was applied in September for 12 days. By constructing a polythene tunnel, the considerable temperature was elevated in the daytime and uncovered at night [41]. During heat stress, the temperature of control and tunnel was recorded twice a day [42]. The maximum and minimum temperature was measured inside the tunnel with the help of a digital thermometer ( Figure 1).

Biochemical and Yield-Related Parameters
After heat stress of 12 days, data were recorded for biochemical characters from five selected plants under normal and stressed conditions from each replication. Young leaves were collected from the top of the plant for biochemical analysis from each genotype. Enzyme extraction was performed with 0.5 g of cotton leaf sample in chilled condition with 1-2 mL of potassium phosphate buffer (pH 7.8). The prepared mixture was then centrifuged for 5 min at 1400 rpm. The supernatant was used for biochemical traits quantification by using a UV spectrometer at different wavelengths. The H 2 O 2 was determined by following the protocol proposed by [43], catalase and peroxidase by [44], and chlorophyll contents by [45]. Data regarding agronomic traits plant height (cm), number of bolls, boll weight (g), seed cotton yield (g) and seed index (g) were recorded at maturity from five randomly selected guarded plants of each accession.

Biochemical and Yield-Related Parameters
After heat stress of 12 days, data were recorded for biochemical characters from five selected plants under normal and stressed conditions from each replication. Young leaves were collected from the top of the plant for biochemical analysis from each genotype. Enzyme extraction was performed with 0.5 g of cotton leaf sample in chilled condition with 1-2 mL of potassium phosphate buffer (pH 7.8). The prepared mixture was then centrifuged for 5 min at 1400 rpm. The supernatant was used for biochemical traits quantification by using a UV spectrometer at different wavelengths. The H2O2 was determined by following the protocol proposed by [43], catalase and peroxidase by [44], and chlorophyll contents by [45]. Data regarding agronomic traits plant height (cm), number of bolls, boll weight (g), seed cotton yield (g) and seed index (g) were recorded at maturity from five randomly selected guarded plants of each accession.

Statistical Analysis
Initially, the data collected for all traits under study regarding parents have been analyzed to validate whether the parental groups (tolerant and susceptible) performed significantly differently within and across normal and heat stress conditions. For this purpose, interaction profilers were plotted using the least square fit model with the help of the statistical software package SAS-JMP Pro 16 (SAS Institute Inc., Cary, NC, USA, 1989-2011). After confirming the existence of a reasonable amount of differences among parents under normal and heat stress conditions, all collected data were subjected to standard least-square effect fit to all the variables under study, considering the replications, genotypes and treatments. The data was then averaged for replications regarding treatments and genotypes and the mean data were then subjected to multivariate analyses, including correlation, PCA, and cluster analyses. Pairwise correlations have been calculated to make a Correlogram through Pearson correlation method. The means for genotype and treatments were used for correlation-based principal component analysis using default functions for standardization/scaling by JMP software: SAS-JMP Pro 16 (SAS Institute Inc., Cary, NC, USA, 1989-2011). Cluster analysis was performed using the Hierarchical clustering approach with two-way clustering to generate a tree diagram based on Euclidean distances by Ward's method. To decide on a number of clusters, a distance graph was followed representing the threshold for selection of the number of clusters by determining when the distances between clusters no longer appear to be of practical importance. Default standardization option in JMP Pro 16 SAS-JMP Pro 16 (SAS Institute Inc., Cary, NC, USA, 1989-2021) has been selected for analysis to standardize the range of the continuous initial variables so that each one contributes equally to the analysis.

Statistical Analysis
Initially, the data collected for all traits under study regarding parents have been analyzed to validate whether the parental groups (tolerant and susceptible) performed significantly differently within and across normal and heat stress conditions. For this purpose, interaction profilers were plotted using the least square fit model with the help of the statistical software package SAS-JMP Pro 16 (SAS Institute Inc., Cary, NC, USA, 1989-2011). After confirming the existence of a reasonable amount of differences among parents under normal and heat stress conditions, all collected data were subjected to standard least-square effect fit to all the variables under study, considering the replications, genotypes and treatments. The data was then averaged for replications regarding treatments and genotypes and the mean data were then subjected to multivariate analyses, including correlation, PCA, and cluster analyses. Pairwise correlations have been calculated to make a Correlogram through Pearson correlation method. The means for genotype and treatments were used for correlation-based principal component analysis using default functions for standardization/scaling by JMP software: SAS-JMP Pro 16 (SAS Institute Inc., Cary, NC, USA, 1989-2011). Cluster analysis was performed using the Hierarchical clustering approach with two-way clustering to generate a tree diagram based on Euclidean distances by Ward's method. To decide on a number of clusters, a distance graph was followed representing the threshold for selection of the number of clusters by determining when the distances between clusters no longer appear to be of practical importance. Default standardization option in JMP Pro 16 SAS-JMP Pro 16 (SAS Institute Inc., Cary, NC, USA, 1989-2021) has been selected for analysis to standardize the range of the continuous initial variables so that each one contributes equally to the analysis.

Results
Results presented as interaction profilers ( Figure S1) indicated a reasonable variation among the tolerant and susceptible genotypes across normal and heat stress conditions. All studied traits have shown a varying trend for parents across different conditions ( Figure S1). A sufficient amount of variation presented in results validated the hypothesis for the use of reference susceptible and tolerant parental accessions to be included for the current study.
After initial validation, the standard least-square effects fit applied on all test accessions across different conditions revealed that all genotypes performed sigSuppnificantly differently for studied characters under normal and heat stress conditions (Tables S2-S10). Summary of effect tests as presented in Table 1 and Figure 2 indicating the least square effects of heat stress on all parents and their hybrids were significant regarding all studied traits. The genotype × treatment interaction was also significant for all studied traits ( Table 1).
The correlation analysis was performed to assess the relationship among nine agronomical and biochemical traits. The Correlogram showed scatter plot, variable distribution and correlation coefficient. Overall, the combined correlation coefficient revealed that the PH, TNB and chlorophyll contents have a significant negative relationship with H 2 O 2 , POD and CAT. The SCY and BW positively correlated with PH, TNB, chlorophyll contents (a & b), while negatively associated with CAT. The SCY also exhibited a negative association with POD. The POD and CAT both revealed a significant negative association with H 2 O 2 ( Figure 3). The plants responded to increased H 2 O 2 by up-regulating defensive enzymes activities such as POD and CAT to metabolize H 2 O 2 . Under heat stress, the box plots showed that H 2 O 2 contents, CAT and POD activities were increased, but PH, TNB, BW, SCY and chlorophyll contents were decreased and revealed a negative association with H 2 O 2 , CAT and POD. The prediction profiler plot was also used to find out the most desirable genotypes for both environments. The desirability graph showed that the genotype FH-458 × FH-313 performed well as compared to all other genotypes for agronomical and biochemical traits under both conditions ( Figure 3).

Results
Results presented as interaction profilers ( Figure S1) indicated a reasonable variation among the tolerant and susceptible genotypes across normal and heat stress conditions. All studied traits have shown a varying trend for parents across different conditions (Figure S1). A sufficient amount of variation presented in results validated the hypothesis for the use of reference susceptible and tolerant parental accessions to be included for the current study.
After initial validation, the standard least-square effects fit applied on all test accessions across different conditions revealed that all genotypes performed sigSuppnificantly differently for studied characters under normal and heat stress conditions (Tables S2-S10). Summary of effect tests as presented in Table 1 and Figure 2 indicating the least square effects of heat stress on all parents and their hybrids were significant regarding all studied traits. The genotype × treatment interaction was also significant for all studied traits (Table  1).

Figure 2.
The least square mean plots for nine studied traits representing a comparison of variables across Normal and Heat stress conditions. The mean plots are depicting spectacular differences regarding traits performance under both conditions.

Cluster Analysis
Hierarchical clustering analysis was executed to clusters the 8 parents and their 15F1 hybrid based on agronomical and ionic traits under normal and heat stress conditions. Under normal conditions, 23 genotypes were clustered into 3 groups. Groups 1, 2 and 3 consists of 12, 3 and 8 genotypes, respectively (Figure 4). Under heat stress, 23 genotypes were also clustered in 3 groups. Group 1, 2 and 3 contained 7, 12 and 4 genotypes, respec-

Cluster Analysis
Hierarchical clustering analysis was executed to clusters the 8 parents and their 15F1 hybrid based on agronomical and ionic traits under normal and heat stress conditions. Under normal conditions, 23 genotypes were clustered into 3 groups. Groups 1, 2 and 3 consists of 12, 3 and 8 genotypes, respectively (Figure 4). Under heat stress, 23 genotypes were also clustered in 3 groups. Group 1, 2 and 3 contained 7, 12 and 4 genotypes, respectively (Figure 4). Under normal and heat stress conditions, different clusters were represented by different colors. Generally, the genotypes of cluster 3 have shown unchanged agronomic performance along with unchanged or higher values of biochemical traits across normal and stress conditions have been designated as heat tolerant. The second cluster have shown a nominal decrease in agronomic traits along with increased biochemical performance have been denoted as moderately tolerant. The third cluster with red color based on declined agronomic performance along with a decreased biochemical performance under stress have been designated as susceptible. Based on clustering, under normal and heat stress conditions, among parents and their hybrids, NIAB-1048 × FH-313, VH-329, NIAB-1048, IUB-65 × CRS-2, and FH-458 × FH-313 were superior genotypes for most of the studied traits. The cluster analysis also validated the prediction of desirability plots and showed that FH-458 × FH-313 performed well in yield and physiological traits under both conditions. The genotypes in obtained clusters under both conditions are given in (Table S11).  (Table S11).

Principal Component Analysis (PCA)
PCA is a multivariate statistical approach to examine complex and large datasets. Based on the association among the studied characters and extracted clusters, the variation pattern in cotton genotypes was also examined by PCA to determine the relationship between different variables that lead to variability in studied data. Under normal and heat stress conditions, 9 PCs were found, and from these 9 PCs, 3 PCs had shown commutative variation >80% (Table S12). Hence dividing explainable variation into 3 PCs and it validates 3 groups already suggested by AHMC clustering method. These 3 PCs shared

Principal Component Analysis (PCA)
PCA is a multivariate statistical approach to examine complex and large datasets. Based on the association among the studied characters and extracted clusters, the variation pattern in cotton genotypes was also examined by PCA to determine the relationship between different variables that lead to variability in studied data. Under normal and heat stress conditions, 9 PCs were found, and from these 9 PCs, 3 PCs had shown commutative variation >80% (Table S12). Hence dividing explainable variation into 3 PCs and it validates 3 groups already suggested by AHMC clustering method. These 3 PCs shared 86.93% of the total variability, while the rest of all PCs contribute 13.07% to the total variability. The PC1 had (55.83%), followed by PC2 (22.81%), PC3 (8.29%), PC4 (4.55%), and PC5 (3.41%) variability among the genotypes for studied traits, respectively ( Figure 5).  (Figure 5). In biplot, traits were depicted in the form of vectors. The relative distance of the traits from origin revealed their contribution to the total variation. A PCA biplot analysis could select characters classified into main sets and subsets based on similarity and dissimilarity. In our data set, two groups of traits were identified in the PCA biplot considering PC1 and PC2. The PC-I and PC-II explained 78.64% of the total variation. The Chla, Chlb, SCY, TNB and PH had long vectors and exhibited higher variation. The characters Chla, Chlb, SCY, and TNB lay close to each other in the direction of PC1; clustered in group I and exhibited higher association.
Interestingly the H2O2, POD and CAT lied close to each other towards the direction of PC2; they clustered in group II and revealed a high correlation among themselves. Interestingly, the PCA biplot revealed that group I traits, the major contributors in PC1, were strongly associated with cluster 3, while the traits of group II, the contributors in PC1, were associated with the genotypes of cluster 3. The biplot of PC1 and PC2 showed that the genotypes FH-458 × FH-313 and IUB-65 × CRS-2 are superior genotypes for studied traits.
The biplot of PC1 and PC3 explained 64.09% variability of the total variability. In the biplot of PC1 and PC3, the traits Chla, Chlb, SCY, and TNB are more discriminating traits and significantly correlated. This biplot showed that the FH-458 × FH-313, IUB-65 × CRS-2 and NIAB-1048 × FH-313 are the most desirable genotypes. The biplot of PC2 and PC3 exhibited the least variability (31.09%) than the biplot of PC1 with PC2 and PC3 independently. In the biplot of PC2 and PC3, CAT, POD and H2O2 were more discriminating traits and positively correlated among themselves ( Figure 5).  (Figure 5). In biplot, traits were depicted in the form of vectors. The relative distance of the traits from origin revealed their contribution to the total variation. A PCA biplot analysis could select characters classified into main sets and subsets based on similarity and dissimilarity. In our data set, two groups of traits were identified in the PCA biplot considering PC1 and PC2. The PC-I and PC-II explained 78.64% of the total variation. The Chla, Chlb, SCY, TNB and PH had long vectors and exhibited higher variation. The characters Chla, Chlb, SCY, and TNB lay close to each other in the direction of PC1; clustered in group I and exhibited higher association.
Interestingly the H 2 O 2 , POD and CAT lied close to each other towards the direction of PC2; they clustered in group II and revealed a high correlation among themselves. Interestingly, the PCA biplot revealed that group I traits, the major contributors in PC1, were strongly associated with cluster 3, while the traits of group II, the contributors in PC1, were associated with the genotypes of cluster 3. The biplot of PC1 and PC2 showed that the genotypes FH-458 × FH-313 and IUB-65 × CRS-2 are superior genotypes for studied traits.
The biplot of PC1 and PC3 explained 64.09% variability of the total variability. In the biplot of PC1 and PC3, the traits Chla, Chlb, SCY, and TNB are more discriminating traits and significantly correlated. This biplot showed that the FH-458 × FH-313, IUB-65 × CRS-2 and NIAB-1048 × FH-313 are the most desirable genotypes. The biplot of PC2 and PC3 exhibited the least variability (31.09%) than the biplot of PC1 with PC2 and PC3 independently. In the biplot of PC2 and PC3, CAT, POD and H 2 O 2 were more discriminating traits and positively correlated among themselves ( Figure 5).
The PCA results showed that the first three components explained 86.931% variability to the total variation. The yield and biochemical characters are illustrated by factor map squared cosines ( Figure 6). It is deduced that the higher values of squared cosines offer an adequate contribution to the particular variable. The PCA results showed that the first three components explained 86.931% variability to the total variation. The yield and biochemical characters are illustrated by factor map squared cosines ( Figure 6). It is deduced that the higher values of squared cosines offer an adequate contribution to the particular variable.

Discussion
Genetic diversity in germplasm allows the plant breeders to select favorable traits for developing new cotton cultivars conferring resistance to abiotic and biotic stresses. In recent years, cotton yield is decreasing alarmingly due to heat stress in hot climatic regions. Heat stress at reproductive phases alters cotton morphology and physiology, greatly reducing yield [10]. Significant effects of heat stress were observed in all genotypes for studied traits. This study observed a significant reduction in PH, TNB, BW, SCY and chlorophyll contents under heat stress [24,46,47] while H2O2, POD and CAT activities were increased. These results align with the previous findings [8,11]. Under heat stress, the increased production of H2O2 contents causes oxidative stress and results in electrolyte leakages [47]. All the studied traits revealed significant correlations in our study, which suggested that change in any one of those characters correspondingly changes the other characters. It indicated that cotton's morpho-physiological traits play an important role under heat stress conditions to determine high temperature response. It means that if one reliable trait is picked in heat stress and used as a selection criterion that will lead to affect other

Discussion
Genetic diversity in germplasm allows the plant breeders to select favorable traits for developing new cotton cultivars conferring resistance to abiotic and biotic stresses. In recent years, cotton yield is decreasing alarmingly due to heat stress in hot climatic regions. Heat stress at reproductive phases alters cotton morphology and physiology, greatly reducing yield [10]. Significant effects of heat stress were observed in all genotypes for studied traits. This study observed a significant reduction in PH, TNB, BW, SCY and chlorophyll contents under heat stress [24,46,47] while H 2 O 2 , POD and CAT activities were increased. These results align with the previous findings [8,11]. Under heat stress, the increased production of H 2 O 2 contents causes oxidative stress and results in electrolyte leakages [47]. All the studied traits revealed significant correlations in our study, which suggested that change in any one of those characters correspondingly changes the other characters. It indicated that cotton's morpho-physiological traits play an important role under heat stress conditions to determine high temperature response. It means that if one reliable trait is picked in heat stress and used as a selection criterion that will lead to affect other traits under heat stress conditions. Different studies reported that the chlorophyll contents measurement is an important marker for identifying the tolerant genotypes in a stress breeding program and suggested that the higher levels of chlorophyll content are associated with heat tolerance [48,49]. In our study, chlorophyll contents revealed a strong positive association with TNB, BW and SCY. Reduction in chlorophyll contents under heat stress can be attributed to increased production of H 2 O 2 . The significant negative correlation between chlorophyll contents and H 2 O 2 established the oxidative stress-induced decrease in chlorophyll synthesis under heat. So, the reduction in chlorophyll contents is associated with decreased SCY. The heat stress destabilizes the photosynthesis process due to the excessive production of reactive oxygen species that speeds up the oxygen-induced cellular damage [50]. Due to heat stress, plants lose their cellular homeostasis by producing reactive oxygen species (ROS) [51]. The plant responded to ROS by up-regulating antioxidants defensive enzymes such as peroxidase (POD) and catalase (CAT) to metabolize ROS [52]. Peroxidase (POD) detoxifies H 2 O 2 inside the cytosol and chloroplasts of the plant cells [53]. The catalase (CAT) converts the toxic H 2 O 2 into water and oxygen. In this study, with the rise of temperature, H 2 O 2 production increased; however, due to the scavenging activity of catalase, its negative effects were averted, as reported in the earlier findings [54]. The CAT and POD activities were higher under heat stress conditions as compared with the control condition. Previous reports have also observed a high level of CAT activity in cotton leaves under heat stress conditions [46]. The genotypes exhibiting a higher level of POD and CAT were found to be optimal for H 2 O 2 and declared as heat-tolerant [47]. The correlation analysis revealed that the SCY was positively correlated with TNB, BW, POD, CAT, and chlorophyll contents (a & b) under normal and heat stress conditions. Our results were in accordance with reported research previously by different studies [55][56][57].
Cluster analysis was performed to explore the genetic variation among 23 cotton genotypes for biochemical and yield-related traits, leading to breeding programs [58]. Hierarchical cluster analysis revealed that genotypes of cluster 3 have the highest ability to tolerate heat stress as judged by almost all studied traits. Our PCA results covering 86.93% of total variation in first three PCs for yield, yield contributing traits, fiber quality traits, and biochemical traits. It can affirm that discrimination of cotton genotypes can be made across the experimental conditions based on this information. Moreover, these selected varieties can also be further utilized in the cotton breeding program against heat stress. Many scientists have confirmed the presence of genetic diversity in plants species through this efficient statistical technique. Our PCA results agreed with other researchers who used this technique to find variability in cotton germplasm [59][60][61]. We observed that the first two PCs added maximum variation existent in experimental material. It was inconsistent with previous reports [62,63]. Under both conditions, among all studied traits in cotton genotypes, TNB, SCY, POD, CAT and chlorophyll contents (a & b) were the major contributing traits observed in PC1. The PC2 was characterized by BW and H 2 O 2 [64]. This technique could help exploit and conserve the germplasm by dividing genetic variability into different components and being used in future breeding programs to increase the efficiency of specific traits. Our results showed that clustering of cotton genotypes varying across heat stress tolerance based on the hierarchical clustering have been validated and confirmed through the results obtained from PCA. Altogether, it is evident that 23 cotton genotypes exhibited significant variation in the heat stress changes in studied traits and the multivariate analyses could be effectively used for the identification of the cotton genotypes of desirable traits for heat tolerance.
Depending upon clustering, under normal and heat stress conditions, among parents and their hybrids, NIAB-1048 × FH-313, VH-329, NIAB-1048, IUB-65 × CRS-2, and FH-458 × FH-313 were superior genotypes for most of the studied traits. Different studies revealed that heat-tolerant species exhibited higher antioxidant enzyme activities in response to heat stress, while heat-sensitive species could not do so. The facts accumulated from data revealed that plants' intrinsic antioxidant resistance system might lay out a strategy to increase heat tolerance.

Conclusions
The constantly warming surroundings of every organism make them change their phenotypes accordingly for their endurance. The rapidly changing environmental conditions across the globe are enabling every genotype to make certain changes in its genome to transcribe subsequent related phenotypes. To sustain these fluctuations, there is a need to hunt for such genotypes in the available genetic stock of cotton and every other crop. keeping in context, the study was conducted to explore the heat resilient genotypes from available assets. The current research results revealed the substantial impacts of heat stress on available cotton genetic resources regarding agronomic and biochemical traits. Based on different biometrical approaches, we found that the genotype FH-458 × FH-313 performed superior to the rest of the experimental genotypes under normal and heat stress conditions. This genotype could efficiently be utilized in future cotton improvement programs to enhance heat stress endurance in high-yielding cultivars to sustain abrupt climatic calamities.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10.3 390/agronomy11091885/s1, Figure S1: Interaction profiler for preliminary performance comparison of tolerant and susceptible parents across normal and stress conditions. Table S1: List of accessions used for evaluation across normal and heat stress conditions. Table S2: Detailed effect report for parameter estimates related to Plant Height under normal and heat stress conditions. Table S3: Detailed effect report for parameter estimates related to TNB under normal and heat stress conditions. Table S4: Detailed effect report for parameter estimates related to Boll weight under normal and heat stress conditions. Table S5: Detailed effect report for parameter estimates related to Seed Cotton Yield under normal and heat stress conditions. Table S6: Detailed effect report for parameter estimates related to H 2 O 2 under normal and heat stress conditions. Table S7: Detailed effect report for parameter estimates related to POD under normal and heat stress conditions. Table S8: Detailed effect report for parameter estimates related to CAT under normal and heat stress conditions. Table S9: Detailed effect report for parameter estimates related to Chlorophyll-A under normal and heat stress conditions. Table S10: Detailed effect report for parameter estimates related to Chlorophyll-B under normal and heat stress conditions. Table S11: Distribution of cotton genotypes in different clusters based on studied traits.   Data Availability Statement: Data supporting reported results is available and provided at reasonable request.