Identification and Genetic Mapping of Potential QTLs Conferring Heat Tolerance in Cotton (Gossypium hirsutum L.) by Using Micro Satellite Marker’s Approach

High-temperature stress can cause serious abiotic damage that limits the yield and quality of cotton plants. Heat Tolerance (HT) during the different developmental stages of cotton can guarantee a high yield under heat stress. HT is a complex trait that is regulated by multiple quantitative trait loci (QTLs). In this study, the F2 population derived from a cross between MNH-886, a heattolerant cultivar, and MNH-814, a heat-sensitive variety, was used to map HT QTLs during different morphological stages in cotton. A genetic map covering 4402.7 cm, with 175 marker loci and 26 linkage groups, was constructed by using this F2 population (94 individuals). This population was evaluated for different 23 morpho-physiological HT contributing traits QTL analysis via composite interval mapping detected 17 QTLs: three QTLs each for Total Number of Sympodes (TNS), Length of Bract (LOB), and Length of Staminal-column (LOS); two QTLs for First Sympodial Node Height (FSH), and one QTL each for Sympodial Node Height (SNH), Percent Boll set on second position along Sympodia (PBS), Total Number of Nodes (TNN), Number of Bolls (NOB), Total Number of Buds (TNB), and Length of Petal (LOP). Individually, the QTLs accounted for 7.76%–36.62% of phenotypic variation. QTLs identified linked with heat tolerance traits can facilitate marker-assisted breeding for heat tolerance in cotton.


Introduction
Cotton is a miracle of the plant realm as it fulfills most of the vital needs and provides more than 90% of the world's total production of fiber for the textile industry and edible oil for almost half of the world's population [1]. It has been observed that more than 50% of cotton around the globe is affected by abiotic stress such as salinity, drought, and heat stress that lead to deficient production of this field crop, especially when affected at the seedling stage [2]. Cotton growth requires sufficient fresh water for better fiber quality, but if it faces drought or heat stress the fiber production is reduced [3]. Many new drought tolerant cultivars of cotton have been introduced with improved plant growth, and even other genetically engineered genotypes of cotton by breeding techniques are being cultivated that can tolerate many abiotic stresses [4]. However, the genetic basics and amendments behind these stresses need to be evaluated more to combat these problems from the genetic roots. Cotton is divided into eight genomes (groups) from A to G and K including 45 diploids and the basic seven tetraploid [5,6]. Evolutionary data based on DNA sequencing suggested that about six to seven million years, ago due to trans-oceanic dispersal, D genome divergence gave rise to the A genome and in America (primarily Mexico), it became a separate lineage [7,8]. An incredible diversification occurred over this time that resulted in the worldwide spread of the Gossypium species. Domestication of wild varieties of cotton by human beings resulted in lot of change in all phenotypic and genotypic characteristics.
In terms of production, Pakistan is at the fourth position among the cotton growers of the world; raw cotton exported from Pakistan holds third position in the world as per records of 2012-2013 [9]. Pakistan is more prone to climate changes due to its geographical location [10]. Heat stress is a combination of different intricate functions of intensity duration of temperature. Because of its geographical position, in Pakistan during the summer in some locations, the temperature reaches up to 50 • C and the scorching heat adversely affects cotton plants. Cotton is cultivated in hot areas in Pakistan [11]. High temperature affects growth and development of the plant as well as fiber quality traits [12,13]. Episodes of periodic heat stress and increase in average temperature for the full season enhances the detrimental effects on almost all the factors of plant growth, and that is the reason there is great reduction in the seed number, fiber quality, and content [14]. Cotton yield is suppressed when the plant faces heat and drought stress due to decreased plant transpiration and reduced biomass accumulation, resulting in an inadequate yield [15]; these stresses adversely affect cell elongation, differentiation, and division and also suppress stomatal conductance [16].
The cotton plant has a wide range of adaptability [17], but high temperature is one of the major constraints in cotton productivity and greatly reduces seed cotton yield and quality, which can be addressed by breeding methods. Marker-assisted selection fastens the breeding technology with an accurate approach towards the desired phenotypic traits among the breeding population [18], and it requires detection and analysis of genetic variations using advanced genetic approaches, leading to phenotypic traits of quantitative and agro-economic importance [19]. Genomic selection (GS) and MAS developed by molecular markers techniques has made it possible to map quantitative trait loci (QTL) and identifying QTLs for high-temperature stress and breeding heat-tolerant varieties is an effective way to address this issue. MAS methodology has been used globally to acquire ordered and swift ways for cotton improvement on large scales internationally, with both highly demanded attributes like high seed production and excellent quality of fiber [20]. For dissection of QTLs related to traits with complex genetic patterns of inheritance, molecular marker use has been an efficient tool and these markers have also facilitated MAS breeding [21].
Both agronomic and economically important traits are approached by researchers for obtaining the aim of better yield of cotton [22]. The main challenging goal for current cotton breeders is to further enhance cotton production. However, this aim is hindered by the use of locally available germplasm and extreme environmental fluctuations that influence yield attributing traits [23,24]. Certain different genes cause different expressions of characters regarding tolerance of heat stress at vegetative and reproductive growth stages [25]. Genes attributing to relative water content, stomatal conductance, especially along with Percent Boll set on the First Position along Sympodia (PBF), Percent Boll set on the second Position along Sympodia (PBS), Cell Injury (CIY), Boll Number (BON), Total number of Buds (TNB), Size of Petiole (SOP), Total number of Flowers (TNF), Length of Bract (LOB), Length of Petal (cm) (LOP), Length of Staminal Column (LOS), Length of Pistil (LPI), and Proline Con. (µg mL −1 ) (PCO) have been reported as crucial for heat stress determination [26,27]. Therefore, during the selection of heat tolerant varieties, both vegetative and reproductive traits should be considered equally.
Molecular genetic methods, especially molecular markers, have been applied widely in cotton in recent couple of decades. Recently, the development of molecular markers was accelerated with the release of assembled genome sequences of G. hirsutum [28,29]. Numerous genetic linkage maps including the intraspecific map of G. hirsutum have been constructed using restriction fragment length polymorphisms (RFLPs), simple sequence repeats (SSRs), and single nucleotide polymorphisms (SNPs). Thousands of quantitative trait loci (QTLs) for yield and fiber quality in cotton have been documented in Cotton QTLdb, Release 2.3 [30,31]. However, there are few studies about the simultaneous dissection of the genetic basis underlying complex traits and their genetic correlations in multiple upland cotton populations by QTL mapping. In the situation of changing weather and elevating temperature around the globe, it is of the utmost importance to recognize QTLs for morphological, architectural, and physiological traits that are directly or indirectly affected by high heat stress at some stages of cotton plant development. This study was conducted to identify and map quantitative trait loci (QTLs) conferring heat tolerance in an Intraspecific cross and used microsatellite markers to identify polymorphism between two upland cotton cultivars in the scorching heat of Multan (Pakistan) during summer. QTL identified in this project could be helpful for future cotton growers of high-temperature regions in the world.
In this study, F 2 populations were used, which were derived from hybridization of two G. hirsutum normal lines (MNH-886 and MNH-814). The corresponding genetic linkage map was constructed using 175 polymorphic SSR markers. QTL mapping was implemented with the integration of the genotypic and phenotypic data of twenty-three agronomic and economic traits contributing towards heat tolerance; the aim of this study was to (a) screen cotton cultivars for heat tolerance, (b) select diverse cultivars as parental lines and then their assessment by SSRs for parental survey, (c) develop the segregating/mapping population (F 2 ) of selected parents and collect phenotypic trait data at different time intervals, (d) survey the F 2 population by polymorphic markers obtained from the parental survey, (e) evaluate phenotypic traits with the association of genotypic markers (SSR) data, (f) identify QTLs directing heat tolerance by QTL cartographer software, and (g) construct a genetic linkage map of Gossypium from the obtained information. The outcomes of this study will help plant breeders to produce heat-resistant varieties that will help farmers and countries with agriculture-dependent economies, especially in high-temperature areas around the globe.

Materials and Methods
This study was conducted to identify and map QTLs conferring heat tolerance in an Intraspecific cross and used microsatellite markers to identify polymorphism between two upland cotton cultivars in the scorching heat of Multan (Pakistan) during summer. QTL identified in this project could be helpful for future cotton growers of high temperature regions in the world. The research was arranged at Cotton Research Station (CRS) Multan to coincide the reproductive phase with higher temperature. The field work encompassed 14 cultivars sown in Randomized Complete Block Design (RCBD) replicated three times during the year 2012. All fourteen cultivars were tagged randomly altogether to evaluate 23 morphological and physiological parameters contributing to heat tolerance for identifying the genomic regions under plant breeding techniques; F2 generation was observed for screening purposes.

Heat Stress Estimation
Heat stress was measured in plants that were sown late in the month of April and traits were compared to plants sown earlier in May because temperature in the latter was higher than 46-48 • C during the research time period and the heat stress-related 23 morphophysiological traits were observed to be affected by temperature in late-sown irrigated conditions. The heat was estimated by a weather forecast taken from the automated metrological station of cotton research station, Multan, as given in the Table 1: Heat-susceptible and -resistant varieties (MNH-814 and MNH-886 respectively) were selected on the basis of data for relative water content, osmotic potential, cell injury, and proline concentration. Relative water content was measured by the following [32] formula: Cell Injury (CIY) was measured when the crop was 55-60 days old, and a sufficient number of leaves was taken from the upper portion and stored in a paper bag. By the use of a punching machine, 15-25 discs of 1 cm diameter were cut. With distilled water, leaf discs were washed three times, were put in test tubes, and then the test tubes were filled up to 40 mL with distilled water. Three sets were made each of 14 test tubes containing leaf discs of 14 cultivars. The first set of test tubes was kept at room temperature as control and the electrical conductivity of the water was noted. The second set was heated at 48 • C for 45 min in a water bath. When water was cooled after 6 h, its electrical conductivity was recorded, while the third set of test tubes was autoclaved at 15 lbs (pressure) for 15 min and electrical conductivity was noted when water was cooled.
The greater the EC, greater the damage caused to plant cells due to heat stress as the maximum number of electrolytes came out of the cell due to cell injury. Consequently, cell injury was also greater. Cell injury was expressed in percentage. Proline is an organic compound synthesized from glutamine. It is located in cytoplasm under stressed conditions as nontoxic compatible organic solute to compensate for the dehydrating effects of high osmotic pressure in the vacuole and in the external media. The proline concentration at 700 mol m 3 was not inhibitory to enzymes and develops in consequences of poor plant growth under toxic effects. Therefore, its exogenous application should promote tolerance [33]. Different workers stated that upon heat stress, when starch and protein synthesis are inhibited, proline might be used by the plant for growth [34,35]. Proline from different tissues was measured by Spectrophotometry based on the method of ref. [36]. , and Proline Con. (µg mL −1 ) (PCO), contributing to heat tolerance to identify the genomic regions. Arithmetic means of three replicates were calculated for fourteen cultivars for each characteristic. The data were compared. The variance and standard deviation were also calculated. The computation of trait correlation was carried out using Minitab Inc., University Park, PA, USA and the following shortlisted traits had considerably varying phenotypes among two genotypes, i.e., MNH 886 and MNH 814 ( Figure 1). growth under toxic effects. Therefore, its exogenous application should promote tolerance [33]. Different workers stated that upon heat stress, when starch and protein synthesis are inhibited, proline might be used by the plant for growth [34,35]. Proline from different tissues was measured by Spectrophotometry based on the method of ref. [36]. , and Proline Con. (μg mL −1 ) (PCO), contributing to heat tolerance to identify the genomic regions. Arithmetic means of three replicates were calculated for fourteen cultivars for each characteristic. The data were compared. The variance and standard deviation were also calculated. The computation of trait correlation was carried out using Minitab Inc., University Park, PA, USA and the following shortlisted traits had considerably varying phenotypes among two genotypes, i.e., MNH 886 and MNH 814 ( Figure 1).

Mapping Population
Based on highly significant differences between two parental lines, the F2 population was developed by self-pollinating F1 plants from a cross between upland cotton line MNH-886 (a heat-tolerant cultivar), and MNH-814 (a heat-sensitive cultivar), and was used to map HT-QTLs during different morphological stages in cotton. Five plants were tagged at random in each line for recording physiological traits data. Ninety-four plants from the F2 population were selected to derive phenotypic and molecular data along with two parents. The experimental field area of Cotton Research Station Multan under natural conditions was selected for experiment to coincide the reproductive phase with higher temperature.

Phenotypic Data Collection Statistical Analysis
Selected parental lines and 94 F2 individuals' phenotypic data were collected from fields at different time intervals. Arithmetic means of 3 replicates were calculated for each parent for each characteristic. The data for heat characteristics were compared. The computation of trait correlation was carried out using Minitab Inc., University Park, PA, USA.

Microsatellite Analysis
Laboratory techniques for DNA extraction were performed as described by Peterson. Amplification reactions were carried out in 15 uL reaction volumes containing 30 mg genomic DNA, 1.0 μM each of SSR primers sequences, which were drawn from the following sources: BNL primers from the Research Genetics Co. (Huntsville, AL, USA, http://www.resgen.com, accessed on 7 April 2022); JESPR primers [37]; CIR primers [37]; and NAU primers [38,39], 100 uM each of dATP, dCTP, dGTP, and dTTP, 1 unit of Taq

Mapping Population
Based on highly significant differences between two parental lines, the F 2 population was developed by self-pollinating F 1 plants from a cross between upland cotton line MNH-886 (a heat-tolerant cultivar), and MNH-814 (a heat-sensitive cultivar), and was used to map HT-QTLs during different morphological stages in cotton. Five plants were tagged at random in each line for recording physiological traits data. Ninety-four plants from the F 2 population were selected to derive phenotypic and molecular data along with two parents. The experimental field area of Cotton Research Station Multan under natural conditions was selected for experiment to coincide the reproductive phase with higher temperature.

Phenotypic Data Collection Statistical Analysis
Selected parental lines and 94 F 2 individuals' phenotypic data were collected from fields at different time intervals. Arithmetic means of 3 replicates were calculated for each parent for each characteristic. The data for heat characteristics were compared. The computation of trait correlation was carried out using Minitab Inc., University Park, PA, USA.

QTL Mapping
Genetic mapping and QTL analysis were performed on each population separately and combined across populations. Linkage maps were constructed using MAPMAKER/Exp Version 3.0b software [41]. QTLs were identified by composite interval mapping [42] using Windows QTLs Cartographer 2.5 [43]. A LOD threshold of 3.0 was used [44]. Marker's order was confirmed with the "ripple" command. Recombination frequencies were converted into map distances (cm) using the Kosambi mapping function [45].
Tests for independence of QTLs were also conducted when 2 or more QTLs of a trait were located on the same chromosome [46]. QTLs were declared significant if the corresponding LR score were greater than 11.5 (equal to a LOD score of 2.5). The proportion of the phenotypic variation explained by each QTL was calculated as R 2 (%) = Phenotypic variability explained by QTL/all of the variation in the population × 100. The total phenotypic variance explained together by all the putative QTLs for each trait was estimated by fitting a multiple-QTL model in the Mapmaker/QTL program.  (Table 2).

Average Performance of Cotton Varieties Based on Morpho-Physiological Traits
Phenotypic distribution of F 2 population for morpho-physiological traits is shown in Figure 2a,b. The phenotypic values of morpho-physiological traits are presented in Table 3. Twenty-one morpho-physiological traits displayed a normal distribution (skewness did not exceed 1.0), while two traits, TNF and NOB, showed a non-normal distribution. These results indicated the trend of having major QTL involvement in this population and it was thus suitable for QTL analysis.

Stress Determining Physiological Traits
Physiological traits measure the response of plants to different phenomena taking place internally, such as cell injury and production of certain proteins, such as proline, in response to heat stress. MNH-886 showed less CIY while exposed to high temperatures, with an average value of 65, while MNH-814 was susceptible to extreme temperatures and CIY was higher with a value of 80. CIM-557 showed 64.2, while MNH-886 showed a significant value of 76.6 for proline content in heat stress. MNH-814 was found as the most susceptible among fourteen experimental cultivars and showed a proline content value of 5.3 under stress (Table 4).

Correlation
Correlation (Figure 3) was observed by OriginPro 8.5 software and it was observed that plant height showed a positive correlation with the number of fruiting branches per plant, total number of nodes, size of petiole and balls, length of bracts, length of petals, and length of pistil but it had no correlation with total number of flowers, whereas plant height was negatively correlated with total number of nodes, first sympodial node height, sympodial nose number bearing first effective boll, sympodia node height, bearing last effective boll, cell injury, total number of sympods, length of staminal column, and proline content. Fully dehiscent anther had a positive correlation with sympodial node number, percent boll set on first position, percent boll set on second position along sympodia, total number of sympods, total number of nodes, size of petiole, boll number, total number of bolls, and length of bract. Hence, the length of petiole, proline contents, sympodial node number, percent boll set on first and second position along sympodia, total number of nodes, size of petiole, branch number, total number of bolls, length of bract, and length of petiol were all positively correlated with each other and a significant effect was observed among the traits.

Correlation
Correlation (Figure 3) was observed by OriginPro 8.5 software and it was observed that plant height showed a positive correlation with the number of fruiting branches per plant, total number of nodes, size of petiole and balls, length of bracts, length of petals, and length of pistil but it had no correlation with total number of flowers, whereas plant height was negatively correlated with total number of nodes, first sympodial node height, sympodial nose number bearing first effective boll, sympodia node height, bearing last effective boll, cell injury, total number of sympods, length of staminal column, and proline content. Fully dehiscent anther had a positive correlation with sympodial node number, percent boll set on first position, percent boll set on second position along sympodia, total number of sympods, total number of nodes, size of petiole, boll number, total number of bolls, and length of bract. Hence, the length of petiole, proline contents, sympodial node number, percent boll set on first and second position along sympodia, total number of nodes, size of petiole, branch number, total number of bolls, length of bract, and length of petiol were all positively correlated with each other and a significant effect was observed among the traits.

Construction and Characterization of Intra Specific Linkage Map
Among the 1450 SSR primer pairs tested on parental lines, 175 markers were found to be polymorphic. These markers were applied on population. Using a LOD score > 3.0, these markers were assigned to 26 chromosomes for population based on the information on the cotton SSR map [47]. The linkage map was constructed for the F 2 population. Each linkage group was assigned to specific chromosome (Figure 4). The linkage maps covered approximately 4402.7 cm (Table 5) with an average distance of 20 cm within the markers which, according to the position of SSR markers, is common with the cotton map [48]. We estimate that we surveyed close to 70% of the cotton map, comparing the length of our map with that of the cotton map. The genetic map for the population was generated by MAPMAKER/version 3.0. Genotypic frequencies deviation from the expected segregation ratio of 1:2:1 for the co-dominant locus or 3:1 for the dominant locus was detected with the legitimacy of the additive-dominance model by means of the Chi square (χ 2 ) method [49].

QTLs Mapping for Traits Associated with Heat Tolerance in Cotton
A summary of statistically important QTLs is shown in Table 6. All QTLs for First Sympodial node height, Sympodial node height, percent boll set along sympodia at 2nd position along sympodia, Total number of Sympodes, Total no. of nodes, number of bolls, Total no of buds, Length of bract, Length of Staminal column, and and Length of petal are shown in Figure 3. A total of 17 regions were recognized that contain the QTLs with LOD value 3.0 and above. The most noteworthy QTLs are described in Table 6.  Two QTLs, qFSHa1 and qFSHa2, for first sympodial node height were detected on chromosome 15 with LOD ≥ 6.0, which explain 35% and 36% of phenotypic variance in F 2 respectively. These Loci detected 35%-36% of the PV value. When two QTLs were assembled together, they explained 71% of the PV value. The additive values for both QTLs were 0.59 and 0.58 respectively ( Table 6). Position of the QTLs on the linkage map is shown in Figure 4.

QTLs for Sympodial Node Height Bearing First Effective Boll Set (SNH)
One QTL, qSNH1, influencing Last Effective Boll Set with a LOD score of 3.42 was detected in the F 2 population and it was located on chromosome 6. Putative QTL in this region accounted for 17% of phenotypic variance. So, this QTL explained 17% of the total phenotypic variance (Table 6 and Figure 4).

QTLs for Percent Boll Set on Second Position along Sympodia (PBS)
In the F 2 population of one QTL, qPBS1, the total influencing number of nodes was identified with a LOD score of 18.21 and it was located on chromosome 26. Phenotypic variance in this region was 14.56% ( Table 6). The additive value for this QTL was 0.69.

QTLs for Total No of Sympodes (TNS)
Two QTLs, qTNSa1 and qTNSa2, on chromosome 03 were detected for a total number of sympodes with LOD values 3.59 and 3.71 respectively. Phenotypic variance was observed between 10.05% and 10.12%, and the additive effect was 6.00 and 6.27 respectively. A total of 22% of phenotypic variance was seen when two QTL were fitted simultaneously. The third QTLqTNSa3 was detected on chromosome 05 with LOD value 3.98. The additive effect was 2.89. Phenotypic variance observed was 16.93%. (Table 6).

QTL for Total No of Nodes (TNN)
On chromosome 23 single QTL qTNN1was detected for total no of nodes with LOD value 4.05. Positive additive effect was seen with value 0.18. Phenotypic variance seen was 12.91% (Table 6).

QTLs for Number of Bolls (NOB)
In experiment, one QTL, q NOB1, for Length of bract was identified on chromosome 26 with accumulative phenotypic variance of 21.52%. The LOD value was 3.80. So, this QTL showed phenotypic variance of 22%. Additive positive effect of q NOB1 was 4.25 (Table 6).

QTLs for Total No of Buds (TNB)
In the F 2 population, one QTL, qTNB1, influencing the Total No. of buds was identified with a LOD score of 3.79 and it was located on chromosome 18. Phenotypic variance in this region was 17.67%. Additive effect was positive with value 1.05 (Table 6).

QTLs for Length of Bract (LOB)
Three QTLs, qLOBa1, qLOBa2, and qLOBa3, for length of bract were detected during analysis. The first QTL was on chromosome 2 with LOD ≥ 3.24 and with a positive additive effect of 0.18. Phenotypic variance observed was 8.59%. The second QTL was detected on chromosome 16 with LOD ≥ 3.01 and a negative additive effect of 0.13 and phenotypic variance 7.76%. The third QTL was detected on chromosome 19 with a positive additive effect of 0.18 and phenotypic variance of 12.91%. When three QTLs were fitted together simultaneously the phenotypic variance was 28% (  (Table 6).

Discussion
The current study was carried out to identify the genetic basis responses of cotton plants under heat stress. The data collected were at the parental line and then after by F 2 generation from which heat-susceptible and heat-tolerant genotypes were selected for the screening process. Initially, the emergence of the first sympodial branch at lower nodes determined the early maturity of cotton plants. Theoretically, it is implicated for the 1st sympodial branch to appear on lower nodes as it is highly correlated with earliness and heat tolerance [50,51]. The strong relationship between early maturity and lower sympodial branch node number was reported in previous studies [52]. It was reported that there was a strong association of the 1st sympodial branch node number and heat tolerance. Highly significant differences were found in analysis of variance for the 1st sympodial node number (Table 4). The data of correlation ( Figure 3) showed a positive correlation of the 1st sympodial node number expressed with all the traits except sympodial node number, present boll set on first position along sympodia, cell injury, and length of pistil. Node number to set the initial fruiting sympodia is a reliable and realistic morphological trait of heat tolerance [53]. Minimum and maximum temperature significantly affected the first sympodial branch with 1st boll [54]. All genotypes under study differed significantly for this trait (Table 4). Hussain et al. (2000) revealed similar results for plant height under heat stress, presenting a familiar correlation among traits that plant height has a positive correlation with the morphological traits under study [55]. Boll development was affected by the high temperature stress as compared with vegetative phase and a similar reduction in boll weight was observed when the temperature fluctuated [56]. Morris (1964) also reported a reduction in cotton boll maturity time at high temperature stress [57]. After screening the genotypes on morphological parameters, one genotypes was selected as tolerant against heat stress and another one was selected as heat susceptible, among others, on the basis of physiological characteristics, i.e., relative water contents, water potential, osmotic potential, cell injury, and proline contents. Highly significant differences were perceived by analysis of variance for all the physiological traits among the genotypes, except photosynthesis rate, which is significant ( Table 3). The membrane structure of plant cells was distorted under severe temperature stress, which caused the increased permeability of membrane. As a result, electrolyte leakage increased and eventually led to cell death [58]. Azhar et al. (2009) measured the heat tolerance in term of relative cell injury percentage in cotton and found that thermal stress-tolerant genotypes were more stable for seed cotton yield and also maintained fiber quality as compared with heat-susceptible genotypes. A significant decrease was observed in leaf relative water content % (RWC) for heat-susceptible genotypes when exposed to heat stress, and similar findings were also obtained by Rahman et al. (2000), Siddique et al. (2000), and Parida et al. (2007) under stress conditions [59][60][61]. Higher leaf relative water content (RWC) could be a criterion for selection of a parent for hybridization to develop stress-tolerant genotypes [62,63]. On the basis of grand mean attained from normal and heat-stress situations, the protein contents was variable among genotypes and Raison et al. (1982) revealed that for temperature conditions above the optimum, significant reticence of photosynthesis takes place, resulting in substantial reduction in protein formation [64].
Finally, it was observed that high heat tolerance is a multigenic trait and its expression is controlled by many QTLs. Almost all the vegetative and floral characteristics of cotton plants were affected adversely because of this stress. The identification of QTLs activated to combat heat stress allowed the estimation of genetic architecture and improvement of heattolerance traits by molecular marker-assisted selection (MAS). A total of 1450 markers were applied, among which 175 SSR markers were observed to be polymorphic and were found to be significant; the observations were also in accordance with some other researchers [65]. In order to dissect the genetic basis of heat tolerance, two upland cotton cultivars (MNH-886 and MNH-884) were selected as parents and an F 2 population was developed. A high LOD (logarithm of odds) value provided strong evidence that the reported QTLs are actually associated with the respective traits. We only reported QTLs whose LOD score values were greater than three and which showed a significant additive or dominance genetic effect. A total of 17 QTLs with different effects on ten morphological and physiological traits such as First sympodial node height (FSH), sympodial node height (SNH), Percent boll set along sympodia on 2nd position (PBS), total no. of sympods (TNS), total no. of nodes (TNN), number of bolls (NOB), total no. of buds (TNB), length of bracts (LOB), length of staminal column (LOS), and length of petal (LOP) were detected in the present study. These QTLs were mapped on chromosome numbers 2, 3, 5, 6, 15, 16, 18, 19, 23, and 26. QTLs for length of petal and length of bracts were located on Chr. 2 while QTLs for total no. of buds and length of staminal column were located on Chr. 18 [66,67]. Likewise, QTLs for Boll no. and Percent boll set along sympodia on 2nd position were located on Chr. 26. Our findings are in accordance with work carried out by [68,69].

Conclusions
The purpose of cotton breeding is to boost and stabilize its yield in abiotic and biotic stress environments and to make cultivars with such physiological and architectural characteristics that can tolerate heat stress conditions. A low level of polymorphism is one of the major constraints for plant breeders and geneticists that can be attributed to the different processes like selection and domestication. It resulted in narrowing genetic shuffling in cotton. The use of an enormous number of SSRs can overcome the constraint of low polymorphism. In this study project, more than 1450 SSRs were assessed and the polymorphism rate was 12%, meaning the genetic diversity level was low owing to Intraspecific cross and segregation distortions. In spite of Intraspecific cross, 17 QTLs were detected by evaluating earlier-used and some novel traits. QTL detection can be attributed to a high rate of diversity in both parents. SSR markers were found best to deal with and easiest to assess polymorphism. The main goal of cotton breeding is to help increase and stabilize its productivity in stress environments and to develop cultivars with morphological traits which can withstand heat conditions. Our data suggest that favorable alleles for morphological traits can be combined to improve heat stress tolerance in cotton. Comparisons could be made to evaluate the consistency of QTL detection for the same trait in various backgrounds, which will help to determine the value of targeting these loci for selection in breeding programs.

Future Recommendation
Such coverage in the localization of QTLs controlling different quantitative traits suggested a close genotypic correlation among these traits or a pleiotropic effect of a single gene. It remains to be tested whether these common genomic regions have pleiotropic effects or there are clusters of tightly linked genes for some related traits in these regions. A more numerous mapping population and more closely spaced markers in the map are needed to determine whether the QTLs correspond to a gene with pleiotropic effects or to several separate but closely linked genes, each controlling a single character.