Genetic effects conferring heat tolerance in upland cotton (Gossypium hirsutum L.)

Abstract

Background

Climate change and particularly global warming has emerged as an alarming threat to the crop productivity of field crops and exerted drastic effects on the cropping patterns. Production of cotton has been dropped down to one million bales from 1.4 million bales since 2015 in Pakistan due to the increase in temperature at critical growth stages, i.e., flowering and boll formation. Keeping in view the importance of cotton in the country, this study was conducted to investigate the genetic effects conferring heat tolerance in six populations (P₁, P₂, F₁, F₂, BC₁ and BC₂) developed from cross-1 and cross-2, i.e., VH-282 × FH-142 and DNH-40 × VH-259.

Results

The results revealed that cross-1 performed better in heat stress as compared with cross-2 for majority of the traits recorded. Boll weight and ginning outturn (GOT) were highly effected under heat stress and had negative correlation with Relative cell injury (RCI). Boll weight, fiber length, fiber strength and fiber fineness were under the control of non-additive gene action, whereas RCI was controlled by additive gene effects. Lower values of genetic advance coupled with higher values of broad sense heritability for these traits except RCI confirmed the role of non-additive genetic effects. Duplicate types of epistasis were recorded for fiber strength in cross-1 in normal condition. However, complementary type of non-allelic interaction was recorded for fiber strength under normal condition, fiber fineness and RCI under heat stressed condition in cross-1. Likewise, boll weight, GOT and fiber length in populations derived from cross-2 in normal condition were also under the influence of complementary type of non-allelic interaction. Significant differences among values of mid parent and better parent heterosis for boll weight in both normal and heat stress condition provided the opportunity to cotton breeders for utilization of this germplasm for improvement of this trait through exploitation of heterosis breeding.

Conclusion

Cross-1 performed better in heat stress and could be utilized for development of heat tolerant cultivar. RCI was under the influence of additive gene action, so one can rely on this trait for screening of large number of accessions of cotton for heat stress. While other traits were predominantly controlled by non-additive gene action and selection based on these should be delayed in later generations.

Background

Cotton belongs to family Malvaceae containing more than 200 genera and about 2 300 species. There are more than 50 species of Gossypium reported till now, which are native to Africa, Australia, Central and South America and Asia, respectively (Fryxell 1992; Wendel and Grover 2015). Out of 50 species, only four are domesticated and widespread. Two diploid (2n = 26) species, namely G. arboreum and G. herbaceum belong to Old World cotton produce only 1% of the total cotton production in the world, whereas two tetraploid (2n = 52) species, namely G. barbadense and G. hirsutum belong to New World cotton produce 94% of the total world cotton production. G. barbadense produces 4%, while G. hirsutum also known as upland cotton produces about 90% of the total cotton production in the world (Lu et al. 1997; McCarty et al. 2004).

Upland cotton is a key source of spinnable fiber and cultivated in more than 61 countries in the world on an area of 29.3 million hectares (ICAC 2018). Cotton and cotton-based industry has a pivoting role in the economy of Pakistan. Pakistan ranks the fourth in terms of area and production in the world after India, China and USA, 3rd in consumption and 2nd in yarn production in the world. Cotton contributes 1% share in GDP, while 55% in total foreign exchange earnings of Pakistan. Cotton was planted on an area of 2.7 million hectares in 2017, showing an increase of 10% over the previous year. About 8% more cotton production, i.e., 11.54 million bales was recorded during 2017/2018 as compared with 2016/2017 where 10.72 million bales was recorded (PCCC 2017). However, in terms of per acre yield (679 kg·hm^− 2), Pakistan is lagging far behind from the major cotton producing countries like Australia (1 816 kg·hm^− 2), China (1 719 kg·hm^− 2), Turkey (1 826 kg·hm^− 2) and USA (985 kg·hm^− 2) (ICAC 2018).

A loss of about one-third of cotton produce was recorded in Pakistan during 2015/2016 due to adverse climatic conditions particularly heavy rains during reproductive phase of crop. But high temperature with dry weather conditions favored the spread of whitefly in 2016 and 2017 which affected the productivity of cotton crop on a wide range of area in Punjab province. In recent times besides drought, salinity, insect pests, diseases and seed quality: high temperature has emerged as a major threat to cotton productivity. It is estimated that the global temperature is increasing by 0.4~0.8 °C/year (PMD 2016). The consequences of high temperature in cotton could be low germination, higher fruit shedding (≥ 30 °C/22 °C), pollen sterility and abortion (Guilioni et al. 1997; Ismail and Hall 1999), unavailability of macro and micro nutrients due to increase in soil pH, higher levels of CO₂ in the air will increase photosynthetic activity resulting in enhanced nutrient requirement of cotton plants.

Keeping in view the importance of emerging threat of climate change, it is need of the day to develop climate smart varieties of cotton which could withstand harsh climatic conditions particularly heat stress due to significant adverse effects on yield of seed cotton. So, this experiment was conducted to explore and understand the genetic mechanisms controlling resistance to high temperature and to identify the potential germplasm having tolerance against heat stress which could be used in breeding programs for the introgression and development of new germplasm of upland cotton.

Materials and methods

Screening of germplasm for heat tolerance

The germplasm consisting of 80 accessions of cotton was collected from various Agricultural Research Institutes and Centers of Pakistan to determine heat tolerant and susceptible parental genotypes. Relative cell injury (RCI) percentage was calculated by using the following formula proposed by Sullivan (1972).

$$ \mathrm{RCI}\%=\left[1-\left[\left(1-\left(\mathrm{T}1-\mathrm{T}2\right)\right)\right]/\left(1-\left(\mathrm{C}1/\mathrm{C}2\right)\right)\right]\Big]\times 100 $$

Where, “T” is EC of heat treated and “C” is EC of controlled samples, subscripts 1 and 2 represent initial and final EC readings, respectively.

Based on means from RCI, two heat tolerant, namely VH-259 and FH-142, and two susceptible genotypes, namely VH-282 and DNH-40 were identified against high temperature (Table 1). First part of this study about the details of 80 accessions and screening procedure has already been published (Salman et al. 2016).

Table 1 List of identified heat tolerant and susceptible genotypes of G. hirsutum L.

Development of populations

The four genotypes were hybridized which are named as cross-1 (VH-282 × FH-142) and cross-2 (DNH-40 × VH-259) in the manuscript. A crossing scheme was designed for the development of various populations, i.e., F₁, F₂, BC₁ and BC₂ to fulfill the criteria of generation mean analysis. BC₁ populations were developed by keeping F₁ as female parent and parent 1 as male parent, whereas BC₂ was developed by using F₁ as female parent and parent 2 as a male parent. Some of flowers were self-pollinated for the development of seed for F₂ population. These populations were developed by using greenhouse and field facilities of the Department.

Assessment of populations for heat stress

Average daily temperature during summer season of last five years was collected from AgriMet to determine the duration of maximum heat stress during cotton crop. By having this information all of populations from cross-1 and cross-2 alongwith parents were planted during 2016–2017 on two different sowing times, i.e., early and late. These two sowings were planned based on temperature data of last ten years. The flowering stage in early sowing coincides with maximum annual temperature whereas late sowing coincides with optimal temperature (Ahamed et al. 2010; Abro et al. 2015). The plant material was sown in the experimental area of the department in triplicate by following randomized complete block design. During planting, plant to plant and row to row distance were maintained at 30 cm and 75 cm, respectively, for optimal supply of nutrition and plant protection practices to get good population except effects of heat stress. At the time of reproductive stage RCI %, boll weight and fiber traits were determined.

Statistical analysis

Analysis of variance among the generations was conducted according to Steel et al. (1997). The populations showing significant differences for certain traits were used to conduct generation mean analysis by following the method described by Mather and Jinks (2013).

Results

Assessment of populations

Mean values of F₁ were higher than F₂, BC₁ and BC₂ populations for all of the traits included in this study except fiber fineness and RCI for both crosses under normal and stress conditions (Table 2). The range of boll weight was found to be 2.08~4 g, GOT 37.11~39%, fiber length 23.31~28.05 mm, fiber strength 23.52~25.17 g·tex^-1 and fiber fineness 4.21~4.68 μg·inch^-1 for all of population derived from two crosses under normal and heat stress conditions.

Table 2 Generation means for boll weight, gin turn out, fibre length, fibre strength, fibre fineness and relative cell injury in two crosses VH-282 × FH-142 (1) and DNH-40 × VH-259 (2) under normal (N) and heat stress (H) conditions

Genetic effects

The additive [d) and dominant [h] effects were statistically significant for boll weight and fiber strength in normal conditions while d and h were involved in the inheritance of fiber fineness and relative cell injury in heat stressed condition in cross-1 (Tables 3 and 4). In cross-2, additive and dominant effects were significant for boll weight and GOT in normal and heat stressed conditions, for fiber length in normal condition, for fiber fineness and relative cell injury in heat stressed condition. It indicates that both additive and dominant genes played an important role in inheritance of these traits.

Table 3 Genetic effects for boll weight, ginning outturn, fibre length, fibre strength, fibre fineness and relative cell injury in cross VH-282 × FH-142 (1) under normal (N) and heat stress (H)

Table 4 Genetic effects for boll weight, ginning outturn, fibre length, fibre strength, fibre fineness and relative cell injury in cross DNH-40 × VH-259 (2) under normal (N) and heat stress (H)

Dominance [h], additive × dominance [j] and dominance × dominance [l] variances referred as non-additive gene action, were significant for RCI under heat stress in cross-1, and this pattern of inheritance was found for boll weight under heat stress and RCI under both normal and heat stress condition in case of cross-2. This indicated that these traits were affected by dominance as main affect and epistasis as interallelic interaction.

Correlation

Genotypic correlation was lower than phenotypic correlation that showed involvement of environmental × genotypic interaction. The correlation analysis revealed that boll weight was significantly but negatively correlated with cell membrane stability at phenotypic level grown in heat stress condition for cross-1 (Table 5 and Fig. 1). Likewise, in cross-2, GOT was also negatively and significantly correlated with relative cell injury under heat stress condition (Table 6 and Fig. 1).

Table 5 Phenotypic (lower diagonal) and genetic correlation (upper diagonal) matrix for boll weight, ginning outturn, fibre length, fibre strength, fibre fineness and relative cell injury in cross VH-282 × FH-142 (1) under normal (N) and heat stress (H) conditions

Fig. 1

Correlation coefficients under normal(N) and heat stressed(H) conditions for cross-1 and cross-2

Table 6 Phenotypic (lower diagonal) and genetic correlation (upper diagonal) matrix for boll weight, ginning outturn, fibre length, fibre strength, fibre fineness and relative cell injury in cross DNH-40 × VH-259 (2) under normal (N) and heat stress (H) conditions

Heritability and genetic advance

Narrow sense heritability was moderate (0.43–0.74) whereas broad sense heritability was found high (0.76–0.96) in both crosses (Table 7 & Fig. 2). Narrow sense heritability was lower than broad sense heritability for all the traits under study in both crosses. Genetic advance was low to moderate for both the crosses under both normal and heat stress conditions and ranged from 0.52 to 16.91 (Table 7 & Fig. 2).

Table 7 Narrow sense heritability (h²ns), broad sense heritability (h²bs), Genetic Advance (GA), Heterosis (Ht) and Better parent heterosis (Hbt) for boll weight, ginning outturn, fibre length, fibre strength, fibre fineness and relative cell injury in VH-282 × FH-142

Fig. 2

Heritability, genetic advance and heterosis for cross-1 and cross-2

Heterosis

Heterosis and Heterobeltiosis was statistically significant for boll weight under heat stress condition for cross-1, while it was significant for cross-2 in normal condition. The values of heterosis and heterobeltiosis were ranged from − 0.2 to 17.47 and − 0.24 to 16.73, respectively, for both crosses under normal and heat stress conditions (Table 7).

Discussion

Cotton production is facing several biotic and abiotic challenges including CLCV (Cotton leaf curl virus), wilting disease, sucking and chewing insect pests, drought and elevated temperature. In recent years, high temperature has been reported as a serious threat to crop productivity (Zafar et al. 2018). So, when cotton is exposed to high temperature for longer duration lead to wilting of leaf (Ahuja et al. 2010; Zahid et al. 2016), shedding of fruiting bodies, i.e., squares, buds, flowers and bolls (Cao et al. 2008; Iqbal et al. 2017), and decreased rate of photosynthesis is reported by Marchand et al. (2005). Therefore, it is need of the time to focus on development of new germplasm which can cope with high temperature.

The biometrical analysis indicated that values of dominance or epistasis were many times greater than those of additive effects which shows that these traits were governed by non-additive gene action (Ahmad et al. 2009; Batool et al. 2013; Iqbal et al. 2013). In contrary, these traits had higher values of broad sense heritability but low values of genetic advance which further elaborated the role of non-additive gene action (Singh and Verma 2018). According to Jagtap (1986) when dominant effects are larger than the additive ones then the intensive selection is required for improvement of these traits and selection may be delayed in later generations. Lower values of narrow sense heritability than broad sense heritability for these traits showed that the environmental component was contributing significantly. Low heritability under heat stress condition validate the role of environmental component as well as genotypic × environmental interaction (Murtaza 2006; Desalegn et al. 2009; Batool et al. 2010).

Correlation study revealed that the boll weight (Farooq et al. 2014) and GOT (Azhar et al. 1999; Farooq et al. 2014) had a significant and negative phenotypic correlation with RCI under heat stress conditions. Under heat stress RCI was increased which resulted in increased transpiration and less assimilation of photosynthates which had an adverse effect on boll weight and GOT. Consequently, injury of cell membrane lead to the disturbance of normal functioning of cell which exerts adverse effects on synthesis of fiber. The information from heterosis help plant breeders to identify the superior parental and certain combinations for the development of hybrids. The same data was also exploited for heterosis where it was known that boll weight had significant and positive values over mid parent and better parent under heat stress condition. This revealed that heterozygosity could increase the weight of boll in cotton. A significant gain in boll weight due to heterosis has been reported by several researchers (Abd-El-Haleem et al. 2010; Panni et al. 2012; El-Refaey and El-Razek 2013). Th progeny from cross-1 exhibited high heterosis under heat stressed condition due to the involvement of heat tolerant donor parent, i.e., FH-142. The germplasm derived from cross-1 had higher values for majority of the traits under heat stress which showed that this particular combination was comparatively more heat tolerant than cross-2. On the basis of these results one could conclude that FH-142 and VH-282 could be the desirable parents for their utilization in the breeding programs for the development of heat tolerant germplasm.

Conclusion

All of traits studied in this experiment were predominantly controlled by non-additive gene action except RCI where additive gene action was involved. Therefore, selection based on RCI could be reliable for development of heat tolerant varieties. It is determined that heat stress had role in reduction of boll weight and GOT, which are important yield contributing parameters. The parental lines VH-282 and FH-142 performed better under normal and heat stress conditions and could be utilized for the development of new germplasm for high temperature areas in addition fiber related parameters can be improved.