Saifullah Abro1
Muhammad Rizwan- 1Plant Breeding and Genetics Division, Nuclear Institute of Agriculture (NIA), Tando Jam, Pakistan
- 2Technical Services Division, Nuclear Institute of Agriculture (NIA), Tando Jam, Pakistan
Heat stress in cotton reduces its productivity. The development of heat-tolerant cotton varieties having resilience against changing climate is feasible. The purpose of this study was to probe the genetic variability in upland cotton for heat tolerance, the association of cell membrane thermostability (CMT), stomata, and trichome size with cotton adaptation to high temperature and effective breeding strategy to advance the valued traits. Relative cell injury percentage (RCI%) in studied genotypes ranged from 39 to 86%. Seventeen genotypes were found heat tolerant on the basis of low RCI%, heat susceptibility index (HSI<1), higher number of boll/plant, and seed cotton yield (SCY). Scanning electron microscopy (SEM) of heat-tolerant genotypes revealed the presence of different size of stomata (21.57 to 105.04 μm2) and trichomes (177 to 782.6 μm) on leaves of selected genotypes. The regression analysis showed a strong and negative association of RCI% and stomata size with SCY. However, no association was observed between the trichome size, yield, and fiber traits. On the overall location basis, a significant genotype × environment interaction was observed. All selected genotypes produced a higher SCY as compared with check varieties. But the stability analysis showed that the high yielding genotypes NIA-M-30, NIA-80, NIA-83, and CRIS-342 were also wide adaptive with unit regression (bi∼1) and non-significant deviation from the regression line (S2d∼0). The ability for the combination of some heat-tolerant genotypes was estimated by using the line × tester method among nine hybrids along with their 3 testers (i.e., male) and 3 lines (i.e., females). Genotypes, CRIS-342 and NIA-Perkh, were observed as best general combiners for SCY with a negative general combining ability effects for RCI%. Five hybrids showed a positive specific combining ability and heterotic effects for studied traits and also found lowest for HSI. RCI% and SCY/plant displayed higher estimates of heritability and genetic advance, indicating the heritability due to additive gene effects and chances of effective selection. The identified heat-tolerant and wide adaptive germplasm can be further advanced and utilized in cotton breeding programs for developing heat-tolerant cultivars. Selection criteria involving CMT and stomata size concluded to be an effective strategy for the screening of heat-tolerant cotton.
Introduction
Cotton (Gossypium hirsutum L.) is an international agricultural product of which quality and quantity are subject to various whims of nature. It occupies an important position in global status of commercial crops with annual impact of >US$50 billion in the economy of the world (Organisation for Economic Co-operation Development/Food and Agriculture Organization [OECD/FAO], 2019). The lint quality in general, while the quantity of produce, i.e., the seed cotton yield (SCY) particularly is highly sensitive to climatic conditions. It can be seen in case of Pakistan where it was grown across 2.3 million hectares during 2018–2019 with an average per hectare yield of approximately 707 kg/ha (Anonymous, 2018–2019)) compared with 2,320 and 1,765 kg/ha for Australia and China, respectively (International Cotton Advisory Committee [ICAC], 2018). The quality of lint produced is also inferior, having a short fiber length, coarse fiber fineness, and lower uniformity, resulting in a higher import of longer fiber and lower price of locally produced cotton lint (Pakistan Central Cotton Committee [PCCC], 2014). The cotton belt of Pakistan is mainly located in the high-temperature zone where mean maximum temperature often exceed 48°C during the cotton-growing season. The optimum temperature for a successful cotton production is about 35°C, which may vary among cultivars. However, temperatures above this threshold level badly affect the peak time of cotton flowers and boll setting, resulting in excessive evapo-transpiration and abscission/no boll set at lower half of the plant (Bibi et al., 2003). Singh et al. (2007) reported a strong negative correlation between high temperature and the lint yield, resulting in a yield decrease of about 110 kg/ha for each unit increase in maximum temperature.
Identification of genotypes that have a greater ability to withstand the peaks of heat stress coupled with their limited water use is important to enhance the cotton productivity (Cottee et al., 2010). However, suitable selection standards are required for measuring resilience of cotton germplasm against heat stress. Plant phenological traits especially flowering and boll retention capacity in high temperature environments are effective in repeatable heat stress screening environments (Ismail and Hall, 1999). But genotypes with significant G × Y (genotype × year) and G × E (genotype × environment) interaction in the field often compromise genetic advance. Therefore, genetic gain in the cotton yield can be exploited by targeting those traits that are closely associated with plant adaptation to high-temperature environment (Singh et al., 2007). Cell membrane thermostability (CMT) has been reported as most relevant and suitable selection criteria for measuring heat tolerance (Abro et al., 2015). The ratio of trichomes to stomata is also found associated with the efficiency of water use (Galdon-Armero et al., 2018; Lei et al., 2018). Slow progress in development of heat tolerance in cotton cultivars due to its polygenic inheritance has been reported in earlier findings (Jones et al., 2014). A low heritability and time reliance of adaptive traits in high-temperature environments have also been reported as major obstacles (Singh et al., 2007; Jones et al., 2014). However, combining morphological and plant physiological traits for the elucidation of thermo tolerance in cotton is important and expected to result in improved heat tolerance in cotton (Rana et al., 2011). Better understanding of genetic, morphological, and physiological traits conferring heat tolerance in cotton is necessary, if climate change resilience, high yielding, and fine fiber quality of cotton cultivars are to be developed. An ability for the combination of genotypes in general (general combining ability [GCA]) and specific combining ability (SCA), heterotic effects and heritability of the traits coupled with genetic advance are also important in devising an effective breeding strategy for improving heat stress tolerance of cotton genotypes (Azhar et al., 2005).
The present study investigates the genetic variability in upland cotton for heat tolerance and identifies the association of CMT, stomata, and trichome size with cotton adaptation to high temperature and effects on yield and fiber quality parameters. The line × tester analysis was also performed for deciding the effective breeding strategy to advance the valued traits.
Materials and Methods
Study Site and Experimental Material
The research work was conducted at Plant Breeding and Genetics Division, Nuclear Institute of Agriculture (NIA), Tando Jam, Pakistan, between the years 2015 and 2019. The experimental material consisted of fifty-eight diverse cotton genotypes including advance lines developed at NIA using hybridization and radiation-induced mutagenesis. History of breeding material used in the present study is given in Supplementary Table 1.
Evaluation of Cotton Genotypes for Heat Stress Tolerance
Heat tolerance of the experimental material in terms of high CMT, boll retention, SCY, ginning out turn percentage (GOT%), and staple length was tested in two sowing dates (i.e., early and normal) in a randomized complete block design (RCBD) with three replications. The RCBD was used to evaluate the large set of 58 genotypes because of the non-significant within-block variability and also because we were using the same piece of land every season for cotton breeding trials only. Further, the plot size used was smaller than the normal plot size being used for advance line yield trials. One set of experimental material was sown early (15 March, 2015) to expose the genotypes to high temperature (up to 45°C in the months June and July) at the peak flowering and boll setting stage. Another set of same genotypes was planted in normal sowing (15 May, 2015) to experience mild temperature at the flowering and boll setting stage during late August. Cultivars CRIS-134 and CIM-469 were included as heat tolerant and susceptible check, respectively. Each genotype was planted in a plot size of 6.1 × 3.0 m in four rows with 30 cm plant-to-plant and 75 cm row-to-row distance. Standard agronomic and plant protection measures were performed throughout the growing season. From each genotype, five random plants were selected for data recording on number of bolls/plant. At maturity, seed cotton was picked, and the data were recorded on SCY. Further, seed cotton from each repeat of each genotype was ginned using the ginning machine. The staple length (in mm) was measured using the fibrograph by taking forty grams of lint from each sample. GOT% was calculated by following formula.
Analysis of Cell Membrane Thermostability
Relative cell injury percentage (RCI%), a measurement for CMT under a high temperature stress, was determined according to Sullivan (1972). RCI% was measured from fully developed green leaves of 20–22 days during the peak flowering season. A steel punch with a 10 mm inner diameter was used for collection of two sets of leaf disc from both sides of the midrib of selected leaves. One set was used for heat treatment, and other was kept as control. Leaf discs were collected at about 1 pm to 3 pm, and samples were immediately kept in glass vials having 2 ml of deionized water. Leaf discs were thoroughly washed two times with deionized water to rinse any electrolyte adhering on the outside. After the final rinse, 2 ml of deionized water was added in sterilized glass vials and covered with a lid to avoid evaporation during heat treatment. One set of glass vials was kept at room temperature and other at 50°C in a water bath for 1 h. After heat treatment, 10 ml of deionized water was added in vials and kept at 10°C for 24 h to allow diffusion of electrolytes. On the next day, vials were kept at room temperature and shaken three times for mixing of electrolytes. Initial electrical conductivity (EC) was measured using the EC meter (Model HI-933300, Hanna Instrument, United States). After the measurement of EC, the samples were autoclaved for 10 min at 0.10 MPa pressure for releasing all electrolytes. The same procedure for measurement of EC was repeated. RCI% was calculated by using the following formula:
RCI% = 1-[{1-(T1/T2)}/{1-(C1/C2)}] × 100
Where as
T1 and T2 = EC of sap at 50°C before and after autoclaving, respectively.
C1 and C2 = EC of sap at 25°C before and after autoclaving, respectively.
Measurement of Heat Susceptibility Index
Heat susceptibility index (HSI), a measurement of a relative decrease in yield influenced by non-favorable vs. favorable environments, was determined according to Fischer and Maurer (1978).
Heat susceptibility index = [(1-yh/yp)/H]
Where:
Yh = Average yield of genotypes under heat stress
Yp = Average yield of genotypes under non-stress condition or potential yield
H = Heat stress intensity = 1-(average yh of all genotypes/average yp of all genotypes).
Scanning Electron Microscopy of Leaves for Stomata and Trichome Size
For scanning electron microscopy (SEM) study of selected heat-tolerant genotypes along with the susceptible check, prepared leaf specimens were mounted on the stubs of a double-sided cellophane tape. The sputter was covered with Jeol, JFC-1500 ion Sputter device with C-30–50 mm gold. The sample specimens were observed, and pictures were taken by SEM, Jeol, JSM T-200, and Jeol-T6380 at voltage of 5–15 KV with a distinct magnification (Hayat, 1989). The stomata and trichome measurements with 4–5 repetitions were carried from an individual leaf using an SEM measuring bar. The SEM study was conducted at Central Research Laboratory, University of Karachi, Pakistan.
Multienvironment Field Trial of Selected Heat-Tolerant Genotypes for the Stability Analysis
Seventeen selected heat-tolerant genotypes along with one susceptible check (CIM-469) were evaluated for their stability analysis across five different locations of Sindh (i.e., Tando Jam, Matiari, Shaheed Benazirabad, Dadu, and Khairpur) during 2016–2017. The multienvironment trials were conducted in the RCBD with four replications. The plot size was 18 m2 (6 × 3 m) with a row-to-row and plant-to-plant distance of 75 and 30 cm, respectively. Standard agronomic and plant protection measures were adopted throughout the growing season. At maturity, seed cotton was picked, and data were recorded on SCY (kg/plot).
Population Development for Genetic Studies of Heat Tolerance
To determine the pattern of inheritance and genetic variability, three heat-tolerant cotton genotypes (i.e., CRIS-342, NIA-Perkh, and CRIS-134) and three heat-susceptible genotypes (i.e., HariDost, Sohni, and NIA-148) were selected as parents for population development. The genetic population was developed by following the line × tester (3 × 3) method. Heat-tolerant genotypes were used as testers (i.e., male), whereas heat-susceptible genotypes were kept as lines (i.e., female). For crossing, unopened flowers commonly known as candles or buds were hand emasculated in the evening and covered with a butter paper bag. Male parent flowers were also covered with butter paper bags to avoid insect-intervened pollen contamination. Pollination of flowers was carried out in the morning (10:00 AM), and pollinated flowers were again covered with butter paper bags. At maturity, seed cotton was picked from all cross combinations. Ginning was performed, and F0 seeds were collected for evaluation in F1 generation along with parents.
Evaluation of F1 and F2 Generations
Two sets of F1 genetic material (i.e., 9 hybrids) along with six parents were planted in the field (i.e., normal) and glasshouse (i.e., elevated temperature) in the RCBD with 3 replications. The data were recorded on five random plants for RCI% and SCY (g/plant) in each cross and analyzed for combining ability and heterotic effects. Seeds harvested from F1 generation were again evaluated in F2 generation along with parents. Again two sets were planted, one in the field (i.e., normal) and other in the glasshouse (i.e., elevated temperature) conditions. In F2, 10 and 20 plants were maintained in each cross under the glasshouse and field experiment, respectively. The data on five random plants in each cross along with parents were recorded for RCI% and SCY (g/plant). Genetic variance, heritability, and genetic advance were estimated from F2 data.
Statistical Analysis
The data collected were analyzed separately for each parameter and subjected to analysis of variance (ANOVA) following Steel et al. (1997). The means comparison (honest significant difference [HSD] test at alpha 0.05) and correlation matrix were computed using the STATISTIX® VERSION 8.1 software. The values presented are mean of three replicates ± standard error (SE). The stability analysis of heat-tolerant genotypes was carried out for SCY data from replicated trials at multienvironments. The analysis comprised location-wise or environmental ANOVA and the combined ANOVA for any place/environment (Pooled ANOVA). Stability parameters were calculated following the model by Eberhart and Russell (1966). For the combining analysis, the mean squares of the line × tester design, the overall ability of the combination (i.e., GCA), and SCA variances and effects were calculated according to the method proposed by Kempthorne (1957). Mid-parent heterosis and better parent heterosis were worked out as percent mean deviation of the mean F1 performance over the mean performance of the mid-parent and better parent, respectively. Genetic variance was calculated by following Singh and Chaudhary (1999). Heritability was calculated by following Cochran and Cox (1957), and genetic advance was estimated according to Allard (1960).
Results
Response of Different Genotypes of Cotton to Heat Stress
The analysis of variance showed highly significant differences among studied genotypes for RCI%, SCY, and number of bolls/plant. Differences were also observed among sowing dates for number of bolls/plant and SCY. Significant differences were present among sowing dates for RCI%, staple length, and GOT%. The sowing dates × genotypes interaction was also found highly significant for number of bolls/plant and SCY (Supplementary Table 2). The RCI% values of 58 genotypes are presented in Table 1. RCI% in studied genotypes ranged from 39 to 86%. The RCI level was noted >70% in 22 genotypes, while 60–70% in 20 genotypes and 50–60% in 08 genotypes. Eight genotypes showed < 50% level of cell damage, i.e., 40–50%. The tolerant genotypes, 50–60% levels of cell damage, were Sadori, NIA-Perkh, NIAB-111, CRIS-134, and CRIS-342. Whereas, some advance lines such as NIA-80, NIA-82, NIA-83, NIA-84, NIA-Bt1, NIA-Bt2, NIA-M-30, NIA-M31, NIA-HM-48, NIA-H-32, NIA-HM-327, and NIA-HM-2-1 had the lowest RCI level < 50% (i.e., higher CMT). The genotypes NIA-86, NIA-M2, NIA-M32, NIA-HM-1, NIA-H67, CRIS-121, and CIM-496 showed the highest RCI level (i.e., lower CMT). Out of 58, 17 genotypes were identified as more heat-tolerant because of having a higher CMT.
Table 1. Mean performance of 58 cotton genotypes for relative cell injury percentage (RCI%), yield, and fiber traits.
The SCY of 58 genotypes are presented in Table 1. SCY in studied genotypes ranged from 2,215 to 4,115 kg/ha under heat stress and non-heat stress conditions. The selected seventeen heat-tolerant genotypes performed well in both regimes. NIA-M-30 had the highest SCY (4,115 kg/ha) followed by NIA-HM-2 (3,936 kg/ha) and NIA-80 (3,912 kg/ha). NIA-HM-320 produced lowest SCY (2,215 kg/ha) under heat stress conditions. The HSI based on SCY varied in different genotypes (Table 2). In the present study, HSI values ranged between 0.32 and 10.85. Superior heat-tolerant genotypes gave the least values of HSI < 1 and high yield under heat stress conditions. Whereas, other genotypes, which were found as heat susceptible, had a higher HSI.
Table 2. Heat susceptible index (HSI) of 58 cotton genotypes for yield and fiber traits.
For the determination of heat tolerance, the data for bolls/plant in cotton genotypes were also considered. The mean values of number of bolls/plant showed that both heat stress and non-stress conditions had differently affected the boll retention capacity of fifty-eight genotypes. The heat-tolerant genotypes produced maximum bolls/plant in both heat stress and non-stress regimes. While other genotypes produced less number of bolls as a result of high temperature (Table 1). Superior heat-tolerant genotypes produced maximum bolls/plant under both sowing condition indicating that these genotypes could be less affected by heat stress (Table 2).
Table 1 represents the mean values of GOT% of studied genotypes, which ranged from 30 to 40.56% as a result of early and normal sowing dates. Less effects of high temperature stress were recorded in terms of GOT% in NIA-M-30, NIA-HM-48, NIA-H-32, NIA-Bt1, NIA-Perkh, CRIS-342, NIAB-111, CRIS-134, and Sadori. Non-significant differences in mean values were observed for the staple length (in mm). However, the staple length for all genotypes ranged from 27 to 30 mm. A decrease in fiber length can be attributed to high temperature. But in the present study, eighteen genotypes showed less reduction in the staple length and gave the lowest value of HSI < 1 (Table 2).
Scanning Electron Microscopy Study of Stomata and Trichome Size of Heat-Tolerant Cotton Genotypes
Highly significant and significant differences among genotypes were observed for the stomata and trichome size, respectively (Supplementary Table 3). The stomata and trichome size of different genotypes examined with SEM are presented in Table 3, and the difference in stomata and trichome size of heat-tolerant and heat-susceptible genotypes is shown in Figures 1, 2. The minimum value of the stomata size on the surface of cotton leaf was recorded as 21.57 μm2, and the maximum stomata size was recorded as 105 μm2 with ± 23.3 standard deviation. The large size of stomata was observed in genotype CIM-469 (105.0 μm2) followed by genotype NIAB-111 (87.0 μm2). Whereas, two genotypes, namely, NIA-80 (21.6 μm2) and NIA-M30 (22.6 μm2), gave small size stomata. The results of the electron micrograph obtained from the leaf epidermis tissues showed that all genotypes were non-glandular single-cell trichome. All the selected genotypes showed a high density of trichomes and large trichome size as compared with susceptible check. Among them, the highest average size and density of trichome were recorded (782.7 μm) in genotype NIA-81 followed by Sadori (468 μm) and NIA-80 (465 μm). Whereas, the average lowest density and trichome size (177 μm) were observed in the susceptible check (i.e., CIM-469). The standard deviation value for the minimum trichome length of all genotypes was ± 156.40, and maximum trichomes lengths were ± 144.12. The average standard deviation was observed to be ± 136.31.
Table 3. Mean stomata and trichome size of selected 18 cotton genotypes.
Figure 1. Difference in stomata size of heat-tolerant (NIA-80) and heat-susceptible (CIM-469) cotton genotypes revealed through scanning electron microscopy (SEM).
Figure 2. Difference in trichome size of heat-susceptible (CIM-469) and heat-tolerant (NIA-81) cotton genotypes revealed through SEM.
Association of Relative Cell Injury Percentage, Stomata, and Trichome Size With Seed Cotton Yield and Ginning Out Turn Percentage
Coefficient of determination (R2) depicted that the variation in SCY was due to the effects of RCI% and stomata size (Figure 3). This shows that a higher SCY among genotypes was significantly associated with a lower RCI% and stomata size. Therefore, the regression analysis illustrated that the RCI% and stomata size had a strong negative effect on SCY, with coefficients of determination (R2) 0.61 and 0.51, respectively. The trichome size exhibited positive but non-significant relationships with SCY. Relationships between RCI%, stomata size, and trichome size with GOT% were observed to be weak or non-significant.
Figure 3. Relationship of physiological traits with seed cotton yield (kg ha– 1) and ginning out turn percentage (GOT%) of selected heat-tolerant cotton genotypes.
G × E Interaction Study of Selected Heat Tolerant Cotton Genotypes
Mean squares showed highly significant (P < 0.01) differences for mean SCY among genotypes, environments, and the G × E interaction (Supplementary Table 4). Highly significant differences among locations indicated the presence of the G × E interaction in this area. The stability analysis showed the presence of significant G × E interactions for SCY in this region. As per results, all selected genotypes showed higher yield and stability as compare to susceptible check. Among the genotypes, NIA-M30, NIA-80, and NIA-83 produced a high average yield (3.47, 3.47, and 3.31 tons/ha) and wide adaptability with unit regression (bi∼1) and non-significant deviation from the regression line (S2d∼0). These genotypes are high yielding and stable (Table 4) and may be tested over all environments to attain an additional stable performance. Moreover, three genotypes, namely, NIA-M31, NIA-HM48, and NIA-Bt2, were considered suitable for the specific environment by way of expected performance as they possess a high seed yield (ton/ha) with increased regression average (bi > 1). Five other genotypes viz., NIA-Perkh, NIA-Bt2, Sadori, CRIS-134, and CIM-469 were found suitable to poor environments as they showed good performance in the seed yield with regression (bi < 1) and above average non-significant deviation from the regression (Table 4).
Table 4. Seed cotton yield (tons/ha) and stability parameters for SCY of 18 selected genotypes tested over different locations during 2017–2018.
Genetic Improvement Through Hybridization for Heat Tolerance and Genetic Analysis of Relative Cell Injury and Seed Cotton Yield in F1 and F2 Generations
This research work was carried out to find out combining ability variance effects of parents and their F1 hybrids of cotton and their efficient use in the improvement of breeding program. The ANOVA (Supplementary Table 5) reveals that the mean squares for the interaction of temperature regimes were highly significant (P < 0.01) for parents and their crosses. A significant variation was identified for RCI% in parents as well as their respective crosses under field regimes. Parents vs. crosses mean square, an indicator of average heterosis, was also significant under field regimes (P < 0.01), indicating that the mean RCI% among crosses deviated substantially from the mean of their parental cultivars. ANOVA of combining ability for RCI% and SCY (g/plant) was computed separately using two temperature regimes (March and May) indicated a significant variation for RCI% due to lines and testers (P < 0.01) under both regimes in the field.
Mean Performance of Parents and Their Hybrids for Relative Cell Injury and Seed Cotton Yield in F1
In the presents study, lower RCI was interpreted as higher CMT due to heat sensitivity. The variation among parents as well as hybrids was significant (P < 0.01) for RCI under field and glasshouse regimes. The mean performance of parents and their hybrids for RCI% and SCY (g/plant) is presented in Table 5. Among parents, RCI ranged between 43 and 78% in the glasshouse and between 42 and 75% in the field, whereas tester NIA-Perkh, CRIS-342, and CRIS-134 showed minimum RCI in glasshouse (i.e., heat stress) and field (i.e., non-stress) regimes, respectively; however, line NIA-148 showed maximum RCI in both glasshouse and field temperature regimes. Among F1 hybrids, RCI ranged between 42 and 70% in the glasshouse and 42 and 67% in the field regime. NIA-Perkh, CRIS-342, and CRIS-134 among parents had lowest, while NIA-148, Haridost, and Sohni had the highest relative cell injury percent. The crosses NIA-148 × NIA-Perkh, Sohni × CRIS-134, and NIA-148 × CRIS-342 had lowest RCI%, which showed stability in the cell membrane among hybrids was below average in both glasshouse and field regimes. Hybrid NIA-148 × CRIS-134 had maximum RCI%.
Table 5. Mean performance of parents and their hybrids for RCI% and SCY (g/plant) in F1 and F2 generations.
The results regarding stocktickerSCY (g/plant) are presented in Table 5, which showed that the testers stocktickerCRIS-324, NIA-Perkh, and stocktickerCRIS-134 produced the highest SCY, whereas, three lines gave the lowest yield under both temperature regimes. Among the crosses, HariDost × stocktickerCRIS-342, NIA-148 × NIA-Perkh, NIA-148 × stocktickerCRIS-342, Sohni × stocktickerCRIS-342, and Sohni × stocktickerCRIS-134 showed better performance under non-heat stress and heat stress regimes. The HSI for parents and F1 hybrids was ranged from 0.32 to 2.0%. The parents and hybrids with a higher stocktickerSCY performance have also shown lowest values of HSI. The superior genotypes for heat tolerance gave the least values of HSI (HSI < 1) and high yield under heat stress conditions (Figure 4).
Figure 4. Heat susceptibility index (HSI) for SCY (g/plant) among parents and their progenies in F1 and F2 generations.
Combining Ability Variance and Effects
The effects of general combining ability associated with parental varieties were computed for RCI% under glasshouse and field conditions. The scatter diagram (Figure 5) revealed that CRIS-342 and NIA-Perkh could be categorized in high CMT and high GCA group under both temperature regimes, therefore, could be used as a suitable parent in breeding programs aimed to improve CMT in cotton. Cultivars CRIS-134, CRIS-342, NIA-Perkh, and Sohni were recognized in low CMT and low GCA group under both regimes. The variety CRIS-134, NIA-148, and Haridost showed positive GCA effects. The varieties NIA-148 and Haridost were recognized in low CMT and low GCA group under both temperature regimes, indicated that these cultivars will not be suitable as donor parents for CMT studies. For SCY, the parents NIA-Perkh, CRIS-342, and Sohni displayed the positive GCA (Figure 6). Three tolerant (based on CMT) genotypes expressed differential patterns of GCA effects and also show differential response (i.e., SCA effects) in different cross combinations. The combinations NIA-148 × CRIS-342, Haridost × CRIS-342, Sohni × CRIS-134, and NIA-148 × NIA-Perkh had a negative RCI% value, which proved to be the best indicator showing tolerance to heat stress in hybrids (based on CMT) in both regimes. Three hybrids Sohni × CRIS-342, Sohni × NIA-Perkh, and NIA-148 × CRIS-134 showed a positive RCI%, indicating susceptibility to heat stress. Sohni × NIA-Perkh and NIA-148 × CRIS-342 crosses proved as the best combiners for SCY (Figure 7).
Figure 5. Scatter diagram between RCI% and general combining ability (GCA) effects in upland cotton genotypes (parents) in the glasshouse (up) and field (down) data based on F1 generation.
Figure 6. Estimates of GCA effects for SCY (g/plant) in six parents data based on F1 generation.
Figure 7. Estimates of SCA effects for RCI% and SCY (g/plant) in F1 generation of cotton genotypes in heat-stress (glasshouse) and non-heat-stress (field) conditions.
Heterotic Effects in F1 Generation
The heterosis values for RCI% and SCY (g/plant) expressed as the percentage increase or decrease over mid-parent (i.e., relative heterosis or mid parent heterosis) and over better parent (i.e., hetrobeltiosis) are presented in Table 6. Results revealed that there was a general tendency of increase in the cell membrane stability of hybrids over mid- and better parents, respectively. Negative heterotic effects for RCI (%) indicated less damage to cellular membranes and vice versa. In the present study, maximum negative mid-parent and better parent heterosis were observed in HariDost × CRIS-342, Sohni × CRIS-342, NIA-148 × CRIS-342, and NIA-148 × NIA-Perkh, while the maximum positive heterotic effects were produced by hybrids NIA-148 × CRIS-134 and HariDost × NIA-Perkh under both temperature regimes. The results of heterosis for SCY (g/plant) have shown maximum significant positive heterosis and heterobeltiosis under both heat stress and non-heat stress regimes in the crosses HariDost × CRIS-342, Sohni × CRIS-342, Sohni × NIA-Perkh, NIA-148 × CRIS-342, and NIA-148 × NIA-Perkh. However, maximum negative heterotic effects were observed in cross NIA-148 × CRIS-134 in heat stress conditions (i.e., glasshouse).
Table 6. Heterotic effects for RCI% and SCY (g/plant) in F1 generation.
Performance of Parents and Their Hybrids for Relative Cell Injury and Seed Cotton Yield in F2 Generation
The results for RCI and SCY in parents and F2 progenies are presented in Table 5. Among parents, RCI ranged between 40 and 80% in glasshouse and between 38 and 82% in field conditions. Tester NIA-Perkh, CRIS-342, and CRIS-134 showed minimum RCI in both regimes, whereas line NIA-148 showed maximum RCI in both glasshouse and field conditions. Among F2progenies, RCI ranged between 38 and 72% in the glasshouse and 37 to 72% in field conditions. NIA-Perkh, CRIS-342 and CRIS-134 among parents had lowest, while NIA-148, Haridost, and Sohni had the highest relative cell injury percent. The crosses NIA-148 × NIA-Perkh, HariDost × CRIS-342, Sohni × CRIS-342, and NIA-148 × CRIS-342 had a lowest relative cell injury level percentage, which showed the stability in cell membrane among hybrids was below average RCI in both heat stress (i.e., glasshouse) and non-stress (i.e., field) regimes. For SCY, the results showed that the testers CRIS-324, NIA-Perkh, and CRIS-134 produced the highest SCY, whereas, three lines gave the lowest SCY under both temperature regimes. Among the crosses, HariDost × CRIS-342, HariDost × NIA-Perkh, Sohni × CRIS-342, Sohni × NIA-Perkh, and NIA-148 × NIA-Perkh showed better performance under both temperature regimes.
Values of HSI in parents and F2 progenies (Figure 4) were ranged between 0.2 and 2.3%. Genotypes NIA-Perkh, CRIS-342, and CRIS-134 gave the lowest value for HSI. Among the nine crosses, five crosses HariDost × CRIS-342, Sohni × NIA-CRIS-342, Sohni × NIA-Perkh, NIA-148 × NIA-Perkh, and NIA-148 × CRIS-342 gave the lowest value for HSI. The superior genotypes for heat tolerance gave the least values (HSI < 1) and high yield under heat stress conditions.
Genetic Parameters of Relative Cell Injury and Seed Cotton Yield in F2 Generation
As far as heritability and genetic advance for RCI% are concerned, crosses HariDost × CRIS-342, Sohni × CRI-342, Sohni × NIA-Perkh, and NIA-148 × NIA-Perkh displayed a high heritability and moderate genetic advance under both temperature regimes, while other crosses displayed a low heritability and low genetic advance. Lowest heritability and genetic advance were recorded in cross NIA-148 × CRIS-134 (Table 7).
Table 7. Genetic studies of RCI% and SCY (g/plant) in F2 progenies.
For SCY, out of nine crosses, five crosses showed a high heritability along with high genetic advance. Sohni × CRIS-342, Sohni × NIA-Perkh, NIA-148 × CRIS-342, and HariDost × CRIS-342 displayed the highest heritability coupled with genetic advance under both temperature regimes, which is indicative of the additive type of gene action. However, NIA-148 × CRIS-134 and Sohni × CRIS-134 displayed the low heritability estimates coupled with low genetic advance (Table 7).
Discussion
In the present study, the temperature at peak flowering in early sowing (March) was high (44.2/32.2°C) as compared with normal sowing (39.8/28.4°C) during the month of May. Both temperature regimes in the field were able to expose the studied genotypes to differential responses, resulting in significant sowing the date × genotypes interaction. Previously, Rahman et al. (2004) reported a highly significant genotype × temperature interaction in field. The present study reports seventeen heat-tolerant genotypes, identified on the basis of low RCI% (40–50) and HSI (<1), higher boll retention, and SCY. Differential responses of genotypes to heat in terms of RCI were consistent under heat stress and non-stress regimes. RCI is an indicator of cellular heat tolerance, i.e., high RCI (>70%) means low CMT and vice versa (Khan et al., 2008). Heat-tolerant genotypes produced a high SCY with high GOT% and staple length, while susceptible genotypes were low yielding with inferior fiber traits. Our results are in confirmation with the reports of Azhar et al. (2009) and Karademir et al. (2012). HSI was also used as an extent of relative decrease in yield influenced by non-favorable vs. favorable environments. Results indicated the values of HSI < 1 for heat-tolerant genotypes. Previously, Hafeez-ur-Rahman (2006) measures the heat tolerance in cotton at the whole-plant level using the heat-tolerant index. The SEM study of heat-tolerant genotypes along with susceptible check revealed a wide range of stomata and trichome size present on leaves of these genotypes. Galdon-Armero et al. (2018) reported that ratio of trichomes to stomata in cotton is associated with the efficiency of water use. Lei et al. (2018) observed the linkage of the efficiency of water use with variation between veins and stomata under drought stress. In the present study, a strong negative association of RCI% and stomata size with SCY suggested that low RCI% and small stomata are reliable indicators of genotypic heat tolerance response. This analysis indicated that simultaneous selection of high SCY and heat tolerance would be effective in the presence of heat stress. However, density and length of trichomes on the leaf surface, which protect against excessive radiation and reduce the effect of high temperature, were found non-significant, and these results are in contradiction with Celka et al. (2006) and Passos et al. (2009) who emphasized that trichomes play an important role in reducing biotic and abiotic effects. Another study was conducted for evaluating the heat tolerance and stability in multienvironments. The results showed that all selected genotypes were high yielding and stable as compared with the susceptible check and may be tested over all environments to attain an additional stable performance. These findings are in confirmation with some previous reports (Killi and Harem, 2006; Khan et al., 2007; Riaz et al., 2013), which observed the stability and G × E interaction in genotypes of cotton and the effect on SCY. The ability to the combination of some heat-tolerant genotypes was also estimated by a line × tester method among nine hybrids along with their 3 testers (i.e., male) and 3 lines (i.e., females). Genotypes CRIS-342 and NIA-Perkh were observed as best general combiners for SCY with negative GCA effects for RCI%. This finding is in confirmation with the reports of Jatoi et al. (2010) and Srinivas et al. (2014). Five hybrids showed positive SCA and negative heterotic effects for RCI% and also found lowest for HSI. This is in accordance with some previous studies of Abro et al. (2009) and Nidagundi et al. (2012). These findings suggested that the expression of RCI% under field regimes was predominantly controlled by combining ability because of additive genetic variability. Relative contribution of lines or testers individually as compared with lines × testers was also high for RCI and SCY under both temperature regimes, depicting a comparatively significant role of additive genetic variability in the expression of RCI and SCY. In cotton heat tolerance breeding programs, these results are important for selection point of view from segregating populations (Baloch and Bhutto, 2003). In the present study, RCI% and SCY/plant also displayed higher estimates of heritability and genetic advance, indicating the heritability due to additive gene effects and chances of effective selection. Alike degree of heritability and genetic advance were also got by Khan (2002) and Ahsan et al. (2015). The combination of high heritability and genetic advance provides clear picture of the trait in the selection process. The results for additive gene action are in confirmation with the results of Kumaresan et al. (2000) and Larik et al. (2000).
Conclusion
Selection criteria involving CMT and stomata size concluded to be an effective strategy for the screening of heat-tolerant cotton. Low RCI% and small stomata are good indicators of heat tolerance in optimum as well as heat stress environment. These traits have a good association with SCY. The line × tester analysis revealed two best general combiners for SCY with negative GCA effects for RCI% and five hybrids with positive SCA and heterotic effects for studied traits. Vigilant and rigorous selection of the developed segregating material recognized the heat-tolerant progenies with improved reproductive traits. The identified heat tolerant and wide adaptive germplasm can be further advanced for the release of commercial heat-tolerant cotton varieties.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author Contributions
SaA conceived the idea, executed experiments, recorded the data in the field and laboratory, and prepared the draft manuscript. MR helped in the execution of field and lab experiments, data recording, and analysis and was a major contributor in writing and finalization of the manuscript. ZD helped in execution of field experiments including multilocation trials. ShA performed the data analysis. MS validated the idea, supervised the work, and provided technical inputs in the execution of field experiments. All authors have read and approved the manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.804315/full#supplementary-material
References
Abro, S., Kandhro, M. M., Laghari, S., Arain, M. A., and Deho, Z. A. (2009). Combining ability and heterosis for yield contributing traits in upland cotton (Gossypium hirsutum L.). Pak. J. Bot. 41, 1769–1774.
Abro, S., Rajput, M. T., Khan, M. A., Sial, M. A., and Tahir, S. S. (2015). Screening of cotton (Gossypium hirsutum L.) genotypes for heat tolerance. Pak. J. Bot. 47, 2085–2091.
Ahsan, M. Z., Majidano, M. S., Bhutto, H., Soomro, A. W., Panhwar, F. H., Channa, A. R., et al. (2015). Genetic variability, coefficient of variance, heritability and genetic advance of some Gossypium hirsutum L. accessions. J. Agric. Sci. 7, 147–153. doi: 10.5539/jas.v7n2p147
Anonymous. (2018–2019)). Economic Survey of Pakistan. Islamabad: Ministry of Finance, Govt. of Pakistan.
Azhar, F. M., Ali, Z., Akhtar, M. M., Khan, A. A., and Trethowan, R. (2009). Genetic variability of heat tolerance and its effect on yield and fibre quality traits in upland cotton (Gossypium hirsutum L.). Plant Breed. 128, 356–362. doi: 10.1111/j.1439-0523.2008.01574.x
Azhar, M. T., Khan, A. A., and Khan, I. A. (2005). Combining ability analysis of heat tolerance in Gossypium hirsutum L. Czech J. Genet. Plant Breed. 41, 23–28. doi: 10.17221/3678-CJGPB
Baloch, M. J., and Bhutto, H. U. (2003). Design–II analysis of estimating general and specific combining ability effects of cotton leaf curl virus resistant inbred parents. Zagazig J. Agric. Res. 30, 635–639.
Bibi, A. C., Oosterhuis, D. M., Brown, R. S., Gonias, E. D., and Bourland, F. M. (2003). “The physiological response of cotton to high temperatures for germplasm screening,” in Summaries of Arkansas Cotton Research, Vol. 521, ed. D. M. Oosterhuis (Fayetteville, AR: Arkansas Agricultural Experiment Station), 87–93.
Celka, Z., Szkudlarz, P., and Biereżnoj, U. (2006). Morphological variation of hairs in Malvaalcea L. (Malvaceae). Biodiv. Res. Conserv. 3, 258–261.
Cochran, W. G., and Cox, G. M. (1957). Experimental Designs. New York, NY: John Wiley and Sons. Inc, 661.
Cottee, N. S., Tan, D. K. Y., Bange, M. P., Cothren, J. T., and Campbell, L. C. (2010). Multi-level determination of heat tolerance in cotton (Gossypium hirsutum L.) under field conditions. Crop Sci. 50, 2553–2564. doi: 10.2135/cropsci2010.03.0182
Eberhart, S. A., and Russell, W. A. (1966). Stability parameters for comparing varieties. Crop Sci. 6, 36–40. doi: 10.2135/cropsci1966.0011183X000600010011x
Fischer, R. A., and Maurer, R. (1978). Drought resistance in spring wheat cultivars. I. Grain yield responses. Aust. J. Agric. Res. 29, 897–912. doi: 10.1071/AR9780897
Galdon-Armero, J., Fullana-Pericas, M., Mulet, P. A., Conesa, M. A., Martin, C., and Galmes, J. (2018). The ratio of trichomes to stomata is associated with water use efficiency in Solanum lycopersicum (tomato). Plant J. 96, 607–619. doi: 10.1111/tpj.14055
Hafeez-ur-Rahman. (2006). Environmental interaction, additive and non-additive genetic variability is involved in the expression of tissue and whole-plant heat tolerance in upland cotton (Gossypium hirsutum L). Genet. Mol. Biol. 29, 525–532. doi: 10.1590/S1415-47572006000300022
Hayat, M. A. (1989). Principles and Techniques of Electron Microscopy. Biological Applications, 3rd Edn. Houndmills: Macmillan Press. doi: 10.1007/978-1-349-09857-6
International Cotton Advisory Committee [ICAC] (2018). ‘Cotton: Review of the World Situation’: The 2018/19 Season Featured Rollercoaster Prices, with Decreasing Production, Area and Yields. Available online at: https://icac.org/News/NewsDetails?NewsId=2314&YearId=2019 (accessed April 10, 2021).
Ismail, A. M., and Hall, A. E. (1999). Reproductive stage heat tolerance, leaf membrane thermostability and plant morphology in cowpea. Crop Sci. 39, 1762–1768. doi: 10.2135/cropsci1999.3961762x
Jatoi, W. A., Baloch, M. J., Khan, N. U., Veesar, N. F., and Batool, S. (2010). Identification of potential parents and hybrids in intraspecific crosses of upland cotton. Sarhad J. Agric. 26, 25–30.
Jones, D. C., Gore, M. A., Sanchez, P. A., Fang, D., Hunsaker, D., White, J., et al. (2014). “High- throughput phenotyping of drought-adaptive traits in cotton,” in Proceedings of the Abstracts Belt Wide Cotton Conference. New Orleans, LA.
Karademir, E., Karademir, Ç, Ekinci, R., and Basbag, S. (2012). Screening cotton varieties (Gossypium hirsutum L.) for heat tolerance under field conditions. Afr. J. Agric. Res. 7, 6335–6342.
Kempthorne, O. (1957). An Introduction to Genetic Statistics. New York, NY: John Wiley & sons, 191–200.
Khan, A. I., Khan, I. A., and Sadaqat, H. A. (2008). Heat tolerance is variable in cotton (Gossypium hirsutum L.) and can be exploited for breeding of better yielding cultivars under high temperature regimes. Pak. J. Bot. 40, 2053–2058.
Khan, N. G., Naveed, M., Islam, N. U., and Iqbal, M. S. (2007). Assessment of new upland cotton genotypes (Gossypium hirsutum L.) for yield stability and adaptability. Asian J. Plant Sci. 6, 1012–1017. doi: 10.3923/ajps.2007.1012.1015
Khan, U. Q. (2002). Study of heterosis in fiber quality traits of upland cotton. Asian J. Plant Sci. 1, 593–595. doi: 10.3923/ajps.2002.593.595
Killi, F., and Harem, E. (2006). Genotype × environment interaction and stability analysis of cotton yield in Aegean region of Turkey. J. Environ. Biol. 37, 427–430.
Kumaresan, D., Ganesan, J., and Ashok, S. (2000). Genetic analysis of qualitative characters in cotton (Gossypium hirsutum L.). Crop Res. 19, 481–484.
Larik, A. S., Malik, S. I., Kakar, A. A., and Naz, M. A. (2000). Assessment of heritability and genetic advance for yield and yield components in Gossypium hirsutum L. Sci. Khyb. 13, 39–44.
Lei, Z. Y., Han, J. M., Yi, X. P., Zhang, W. F., and Zhang, Y. L. (2018). Coordinated variation between veins and stomata in cotton and its relationship with water-use efficiency under drought stress. Photosynthetica 56, 1326–1335. doi: 10.1007/s11099-018-0847-z
Nidagundi, J. M., Patil, S. S., Salimath, P. M., Kajjidoni, S. T., Patil, B. C., and Hegde, M. G. (2012). Heterobeltiosis in multiple environments for seed cotton yield and yield attributes in cotton (Gossypium hirsutum L.). Karnataka J. Agric. Sci. 25, 301–304.
Organisation for Economic Co-operation Development/Food and Agriculture Organization [OECD/FAO] (2019). OECD-FAO Agricultural Outlook, OECD Agriculture statistics (database). Available online at: http://dx.doi.org/10.1787/agr-outl-data-en (accessed December 1, 2020).
Passos, J. L., Meira, R. M. S. A., and Barbosa, L. C. A. (2009). Foliar anatomy of the species Lantana camara and L. radula (Verbenaceae). Planta Danin. 27, 689–700. doi: 10.1590/S0100-83582009000400007
Pakistan Central Cotton Committee [PCCC] (2014). Summary Progress Report of National Coordinated Varietal Trials, 2013-14. Karachi: PCCC.
Rahman, H., Malik, S. A., and Saleem, M. (2004). Heat tolerance of upland cotton during the fruiting stage evaluated using cellular membrane thermostability. Field Crop Res. 85, 149–158. doi: 10.1016/S0378-4290(03)00159-X
Rana, R. M., Khan, S. H., Ali, Z., Khan, A. I., and Khan, I. A. (2011). Elucidation of thermo tolerance diversity in cotton (Gossypium hirsutum L.) using physio-molecular approaches. Genet. Mol. Res. 10, 1156–1167. doi: 10.4238/vol10-2gmr1180
Riaz, M., Naveed, M., Farooq, J., Farooq, A., Mahmood, A., Rafiq, C. M., et al. (2013). AMMI analysis for stability, adaptability and GE interaction studies in cotton (Gossypium hirsutum L.). J. Anim. Plant Sci. 23, 865–871.
Singh, R. K., and Chaudhary, B. D. (1999). Biometrical Methods in Quantitative Genetic Analysis. New Delhi: Kalyani Publisher, 1–303.
Singh, R. P., Prasad, P. V., Sunita, K., Giri, S. N., and Reddy, K. R. (2007). Influence of high temperature and breeding for heat tolerance in cotton: a review. Adv. Agron. 93, 313–385. doi: 10.1016/S0065-2113(06)93006-5
Srinivas, B., Bhadru, D., Rao, M. B., and Gopinath, M. (2014). Combining ability studies for yield and fibre quality traits in upland cotton (Gossypium hirsutum L.). SABRAO J. Breed. Genet. 46, 313–318.
Steel, R. G. D., Torrie, J. H., and Dickey, D. A. (1997). Principles and Procedures of Statistics: A Biometrical Approach. New York, NY: McGraw Hill Book Co. Inc, 400–428.
Keywords: breeding, cotton, high temperature, germplasm, physiology, stress tolerance
Citation: Abro S, Rizwan M, Deho ZA, Abro SA and Sial MA (2022) Identification of Heat Tolerant Cotton Lines Showing Genetic Variation in Cell Membrane Thermostability, Stomata, and Trichome Size and Its Effect on Yield and Fiber Quality Traits. Front. Plant Sci. 12:804315. doi: 10.3389/fpls.2021.804315
Received: 29 October 2021; Accepted: 06 December 2021;
Published: 05 January 2022.
Edited by:
Raul Antonio Sperotto, Universidade do Vale do Taquari, BrazilReviewed by:
Zhengwen Sun, Hebei Agricultural University, ChinaPangirayi Tongoona, University of KwaZulu-Natal, South Africa
Copyright © 2022 Abro, Rizwan, Deho, Abro and Sial. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Muhammad Rizwan, cnppX3JpendhbkB5YWhvby5jb20=