The results of the ANOVA suggested significant variations among parents and their F1 hybrids under salt stress conditions. These differences indicated the existence of genetic diversity in the studied germplasm for salinity tolerance, as shown in (Table 4). Moreover, the treatment and Genotypes × Treatment interaction for all characteristics was highly significant, indicating that all genotypes responded differently to salt stress (Table 4). All-cotton genotypes experienced a negative impact on their agronomic and yield-related traits under salt stress. The reduction was observed in agronomic traits such as PH, NBP, BW, SCY, lint%, FL, and FS as shown in (Table 5). Interestingly, fiber fineness (FF) was increased under saline environments. Furthermore, the mean values for H2O2, Na+, CAT, POD, and SOD were increased under salt stress, while the mean values for K+/Na+, K+, TSP, Chla and b were decreased under saline conditions.
The GAM was classified into three categories: low (0–10%), moderate (10–20%), and high (20% or higher). In terms of heritability, values exceeding 80% were considered very high, values ranging from 60–79% were moderately high, values between 40–59% were medium, and values lower than 40% were regarded as low (Johnson, Robinson, & Comstock, 1955). Under saline conditions, high genetic advance as a percentage of the mean (GAM) was observed for BW, CAT, Chla, Chlb, FF, NBP, POD, SCY, SOD and TSP whereas moderate GAM was observed for PH, LP, Na+, K+, FS, FL, and Car. The heritability estimates for BW, CAT, Car, Chlb, FL, FS, K+, K+/Na+, LP, NBP, PH, POD, SOD and TSP, were very high whilst moderate heritability estimates were observed for SCY, Na+, H2O2, FF, and Chla under saline conditions (Table 5).
Table 4
ANOVA split plot design for different agronomic and physiological traits
SOV | Replication | Treatments | Error A | Genotypes | Treatments × Genotypes | Error B |
DF | 2 | 1 | 2 | 23 | 23 | 92 |
BW | 0.14 | 66.59** | 0.05 | 2.04** | 9.94** | 0.14 |
CAT | 23.68 | 2218.17* | 34.99 | 577.3** | 9.37** | 10.30 |
Car | 0.04 | 4.62** | 0.00 | 0.16** | 5.74** | 0.02 |
Chla | 0.13 | 17.32** | 0.11 | 0.41** | 2.86** | 0.07 |
Chlb | 0.003 | 1.20** | 0.00 | 0.04** | 4.17** | 0.00 |
FF | 0.06 | 81.48** | 0.23 | 5.24** | 7.18** | 0.14 |
FL | 0.09 | 5335.82** | 0.28 | 151.9** | 14.64** | 4.89 |
FS | 2.43 | 3235.33** | 1.63 | 127.1** | 8.93** | 5.98 |
H2O2 | 0.00 | 2.84** | 0.02 | 0.2** | 1.32 | 0.02 |
K+ | 104.10 | 12511.4* | 23.50 | 1166.1** | 6.29** | 42.70 |
K+/Na+ | 1.00 | 76.07** | 0.01 | 1.17 ** | 10.94** | 0.14 |
LP | 9.70 | 11507.4** | 19.10 | 356.4** | 14.31** | 14.40 |
Na+ | 120.04 | 7863.99** | 5.16 | 89.38** | 7.23** | 20.32 |
NBP | 7.34 | 4841.84** | 3.55 | 108.73** | 8.11** | 9.96 |
PH | 62.80 | 11612.8** | 38.20 | 429.8** | 5.32** | 30.20 |
POD | 11.82 | 762.86** | 1.48 | 165.7** | 6.32** | 4.11 |
SCY | 23.92 | 5451.73** | 40.26 | 130.77** | 14.36** | 6.27 |
SOD | 0.83 | 304.39** | 0.87 | 22.7** | 17.96** | 0.77 |
TSP | 772.40 | 71679.4** | 53.40 | 73718.5** | 5.36** | 308.20 |
* Significant at the 0.05 level |
** Significant at the 0.01 level |
Table 5
Genetic components of variability, genetic advance percentage means and heritability (broad sense) estimates for studied traits across control and salt stress conditions
SOV | Treatment | Min | Max | Mean | H2b | GAM |
BW | Control | 2.0973 | 4.838 | 3.4328 | 30.61 | 31.39 |
| Stress | 0.121 | 4.129 | 2.0733 | 89.86 | 22.19 |
CAT | Control | 14.05 | 45.6091 | 29.5031 | 84.88 | 20.77 |
| Stress | 16.4609 | 59.7862 | 37.3531 | 93.98 | 44.39 |
Car | Control | 0.43 | 1.49 | 0.7058 | 47 | 39.57 |
| Stress | 0.0863 | 1.2997 | 0.3475 | 83.98 | 16.43 |
Chla | Control | 0.36 | 1.985 | 1.2565 | 14.37 | 40.35 |
| Stress | 0.0607 | 1.6396 | 0.5637 | 75.45 | 33.07 |
Chlb | Control | 0.2059 | 0.6768 | 0.3822 | 12.90 | 44.47 |
| Stress | 0.0145 | 0.533 | 0.1993 | 90.11 | 38.05 |
FF | Control | 0.3462 | 4.4903 | 2.1767 | 93.14 | 18.13 |
| Stress | 2.05 | 5.15 | 3.6812 | 73.79 | 20.48 |
FL | Control | 16.29 | 31.785 | 26.3596 | 78.50 | 14.58 |
| Stress | 2.93 | 28.42 | 14.1857 | 90.92 | 11.19 |
FS | Control | 17.672 | 29.814 | 23.7402 | 25.41 | 21.09 |
| Stress | 4.148 | 27.932 | 14.2603 | 93.25 | 12.27 |
H2O2 | Control | 0.303 | 1.166 | 0.6221 | 68.85 | 36.76 |
| Stress | 0.4579 | 1.3713 | 0.9033 | 62.50 | 38.25 |
K+ | Control | 135.34 | 179.78 | 156.89 | 73 | 13.87 |
| Stress | 100.29 | 178.67 | 138.25 | 88 | 26.30 |
K+/Na+ | Control | 2.77 | 7.22 | 4.12 | 60 | 21.41 |
| Stress | 1.68 | 5.09 | 2.67 | 86 | 30.74 |
LP | Control | 32.1954 | 57.5275 | 44.502 | 16.01 | 13.55 |
| Stress | 7.9245 | 49 | 26.6242 | 94.80 | 11.48 |
Na+ | Control | 20.4 | 56.06 | 38.9845 | 65.03 | 21.92 |
| Stress | 33.18 | 69.98 | 53.7644 | 59.49 | 17.39 |
NBP | Control | 13 | 32 | 24.6389 | 18.82 | 36.23 |
| Stress | 3 | 27 | 13.0417 | 88.20 | 22.56 |
PH | Control | 50.29 | 93.22 | 66.941 | 19.64 | 12.05 |
| Stress | 28.36 | 70.49 | 48.9806 | 90.56 | 14.89 |
POD | Control | 5.25 | 19.7774 | 11.9577 | 81.13 | 54.77 |
| Stress | 7.1406 | 31.7854 | 16.5612 | 90.14 | 82.40 |
SCY | Control | 22.32 | 32.658 | 27.8429 | 81.07 | 29.61 |
| Stress | 5.4634 | 30.3084 | 15.5381 | 73.27 | 61.97 |
SOD | Control | 3.04 | 7.904 | 5.5599 | 55.94 | 19.11 |
| Stress | 4.1312 | 16.845 | 8.4673 | 90.63 | 77.37 |
TSP | Control | 336.8 | 739.6 | 553.3943 | 97.07 | 37.43 |
| Stress | 288 | 736.8 | 508.7726 | 97.96 | 48.25 |
Mean performance of agronomic, and fiber quality traits
Under the influence of salt stress conditions (15dsm− 1), the growth and agronomic characteristics of the examined genotypes were negatively impacted. Notably, the effects of salt stress were particularly evident in terms of plant height, the number of bolls per plant, and boll weight, lint percentage, fiber fineness, and fiber strength for all the genotypes studied along with crosses, as depicted in (Figs. 1 & 2). Under salt stress conditions (15dsm− 1), certain genotypes and crosses demonstrated remarkable growth and yield performance compared to other studied genotypes. Genotypes FBS-Falcon, Barani-333, JSQ-White Hold, Ghauri, along with crosses FBS- FALCON × JSQ- White Hold, FBG-222 × FBG-333, FBG-222 × Barani-222, and Barani-333 × FBG-333 displayed notable resilience and exhibited superior performance for following traits such as PH, NBP, BW, SCY, lint%, FL, FF, and FS. Conversely, the genotypes (FBS-Smart, FBG-222, FBG-333, and Barani-222), and crosses (FBS- FALCON × Ghauri-3, FBS- SMART × JSQ- White Hold, FBS- SMART × Barani-222, FBS- SMART × Ghauri-3, FBG-222 × JSQ- White Hold, FBG-222 × Ghauri-3, Barani-333 × JSQ- White Hold, and Barani-333 × Ghauri-3) demonstrated lower performance for agronomic and fiber quality traits under salt stress conditions. Notably, among the studied crosses (FBS- FALCON × FBG-333, FBS- FALCON × Barani-222, FBS- SMART × FBG-333, and Barani-333 × Barani-222), revealed moderate performance for PH, NBP, BW, SCY, lint%, FL, FF, and FS under salt stress conditions. These specific crosses demonstrated moderate resilience to the adverse effects of salt stress, indicating their potential for further exploration and breeding programs (Figs. 1 & 2).
Biochemical traits
The biochemical traits such chlorophyll a (Chla), chlorophyll b (Chlb), total soluble proteins (TSP), and carotenoid (CAR) exhibited a statistically significant decrease in all genotypes under salt stress (15dsm− 1) (Fig. 3). The application of salt stress treatment led to a more significant reduction in Chla, Chlb, TSP, and CAR levels in specific genotypes and crosses, such as FBS-Smart, FBG-222, FBG-333, Barani-222, FBS- FALCON × Ghauri-3, FBS- SMART × JSQ- White Hold, FBS- SMART × Barani-222, FBS- SMART × Ghauri-3, FBG-222 × JSQ- White Hold, FBG-222 × Ghauri-3, Barani-333 × JSQ- White Hold, and Barani-333 × Ghauri-3. However, the genotypes FBS-Falcon, Barani-333, JSQ-White Hold, Ghauri, FBS- FALCON × JSQ- White Hold, FBG-222 × FBG-333, FBG-222 × Barani-222, and Barani-333 × FBG-333 exhibited the ability to maintain moderate level of Chla, Chlb, TSP, and CAR (Fig. 3). Furthermore, certain genotypes and crosses, such as FBS- FALCON × FBG-333, FBS- FALCON × Barani-222, FBS- SMART × FBG-333, and Barani-333 × Barani-222 exhibited less reduction in Chla, Chlb, TSP, and CAR level compared to susceptible genotypes (Fig. 3).
The application of NaCl treatment had a significant impact on the ionic homeostasis of plants, particularly in relation to Na+, K+, and the K+/Na+ ratio. Notably, the Na+ content in NaCl-treated plants was found to be significantly higher compared to the control group, indicating an increase in sodium accumulation as a result of the treatment. The mean graph of (Fig. 4) depicted that certain genotype and crosses, including FBS-Falcon, Barani-333, JSQ-White Hold, Ghauri, FBS- FALCON × JSQ- White Hold, FBG-222 × FBG-333, FBG-222 × Barani-222, and Barani-333 × FBG-333 lower increase in sodium contents compared other genotypes. Compared to other genotypes, these genotypes also revealed less reduction in potassium contents. Consequently, these genotypes and crosses maintained a higher ratio of potassium to sodium content, similar to that of the control group and performed well for biochemical, agronomic and fiber quality characters under saline environments (Fig. 4). However, specific crosses such as FBS- FALCON × FBG-333, FBS- FALCON × Barani-222, FBS- SMART × FBG-333, and Barani-333 × Barani-222 demonstrated a moderate accumulation of Na+, and exhibited a moderate potassium to sodium ratio. The mean graph reveals that the genotypes and crosses, including FBS-Smart, FBG-222, FBG-333, Barani-222, FBS- FALCON × Ghauri-3, FBS- SMART × JSQ- White Hold, FBS- SMART × Barani-222, FBS- SMART × Ghauri-3, FBG-222 × JSQ- White Hold, FBG-222 × Ghauri-3, Barani-333 × JSQ- White Hold, and Barani-333 × Ghauri-3 exhibited higher levels of Na+, along with a lower potassium to sodium ratio under salt stress conditions and declared as salt susceptible genotypes (Fig. 4).
Under salt treatment (15 dsm− 1), the cotton genotypes displayed a significant rise in the levels of reactive oxygen species (ROS). In contrast to the control group, all genotypes exhibited noticeable increases in the levels of H2O2, indicating the presence of oxidative stress induced by salt stress. However, it is noteworthy that susceptible genotypes and specific crosses, including FBS-Smart, FBG-222, FBG-333, Barani-222, FBS- FALCON × Ghauri-3, FBS- SMART × JSQ- White Hold, FBS- SMART × Barani-222, FBS- SMART × Ghauri-3, FBG-222 × JSQ- White Hold, FBG-222 × Ghauri-3, Barani-333 × JSQ- White Hold, and Barani-333 × Ghauri-3, exhibited greater levels of oxidative damage and showed higher value for H2O2 under saline conditions (Fig. 5A). This observation suggests that these particular cultivars possess lower tolerance to salt stress when compared to tolerant genotypes such as FBS-Falcon, Barani-333, JSQ-White Hold, Ghauri, FBS- FALCON × JSQ- White Hold, FBG-222 × FBG-333, FBG-222 × Barani-222, and Barani-333 × FBG-333. In contrast, moderate levels of H2O2 were observed in the following crosses such as FBS- FALCON × FBG-333, FBS- FALCON × Barani-222, FBS- SMART × FBG-333, and Barani-333 × Barani-222 (Fig. 5A). Theses genotypes also showed a moderate level of SOD, POD, and CAT activities under saline environmental conditions. Our findings indicate that under the salt stress condition of (15 dsm− 1), there was a higher level of SOD, CAT and POD activity observed compared to the control group. The tolerant genotypes exhibited optimal level of SOD, CAT and POD activity under saline conditions. The tolerant genotypes FBS-Falcon, Barani-333, JSQ-White Hold, Ghauri, FBS- FALCON × JSQ- White Hold, FBG-222 × FBG-333, FBG-222 × Barani-222, and Barani-333 × FBG-333 revealed high performance for SOD, POD and CAT as compared to moderately tolerant and susceptible genotypes. The susceptible genotypes FBS-Smart, FBG-222, FBG-333, Barani-222, FBS- FALCON × Ghauri-3, FBS- SMART × JSQ- White Hold, FBS- SMART × Barani-222, FBS- SMART × Ghauri-3, FBG-222 × JSQ- White Hold, FBG-222 × Ghauri-3, Barani-333 × JSQ- White Hold, and Barani-333 × Ghauri-3 showed lesser increase in SOD, POD and CAT activities under saline conditions (Fig. 5B, 5C, 5D).
Correlation analysis
The relationship between various morpho-physiological and biochemical traits is crucial in determining the most suitable genotype under both normal and saline environment. The Pearson correlation coefficients were computed separately for normal and saline conditions. Under normal conditions, FL, Na+, TSP, CAT, FS, PH, FF, POD, and SOD showed positive association with each other. Lint % revealed significant positive relationship with FL, TSP, CAT, FS, FF, POD, Chla, Chlb, BW, and CAR contents under control conditions (Fig. 6). Under salt stress conditions, H2O2 and Na+ revealed significant negative relationship with all morphological, agronomical (SCY, BW, NBP, and LP), fiber quality traits (FL, FS, & FF) and physiological characters (CAT, TSP, SOD, POD, Car, chlorophyll contents, K+/Na+, and K+ respectively (Fig. 2). Under salt stress environments all agronomic characters (SCY, BW, NBP, and LP) revealed positive relationship with physiological traits CAT, TSP, POD, Car, chlorophyll contents, K+/Na+, and K+ respectively. The fiber quality parameters showed positive association with all antioxidant and biochemical traits under salt stress conditions (Fig. 7).
Principal component analysis (PCA) for agro-physiological, biochemical and fiber quality traits
PCA is a statistical approach for analyzing and simplifying complex datasets. We used PCA to evaluate the correlation among variables and the genetic diversity in studied accessions and their association with studied characters under normal and saline conditions. The PC-1 and PC-2 explained 37.6% and 13% of the whole variation under control conditions, respectively (Fig. 8). The relative distance of variables from the origin in PC1 and PC2 showed how much each variable contributes to the total variation among the studied accessions. Under normal conditions, the biplot demonstrated that genotypes of quadrant-IV had high potential for SCY, LP, BW, chlorophyll contents, POD, FS, SOD, K+, and carotenoids contents. The genotypes of quadrant-III were performed better for NBP, TSP, CAT, PH, Na+, and FL (Fig. 8). The genotypes of quadrant-I revealed good performance for potassium to sodium ratio whereas the genotypes of quadrant-II showed optimal performance for H2O2 under control conditions (Fig. 8).
Under saline conditions, the biplot of PC-1 and PC-2, which together account for 90.20% of the total variation. In biplot, traits lied close to each other were positively associated with each other (Fig. 9). Under salt stress conditions, the biplot of PC1 and PC2, revealed positive association between Na+ and H2O2 and the genotypes FBG-222, FBS- SMART, FBG-333, Barani-333 × JSQ-White Hold, FBS- SMART × Ghauri-3, Barani-333 × Ghauri-3, FBG-222 × JSQ- White Hold, FBS- SMART × JSQ-White Hold, FBS- SMART × Ghauri-3, FBS- FALCON × Ghauri-3, Barani-222, and FBS- SMART × Barani-222 showed positive relationship with these traits and declared as salt susceptible. Whereas the genotypes FBS- FALCON × JSQ- White Hold, FBS- FALCON, JSQ- White Hold, FBG-222 × FBG-333, FBG-222 × Barani-222, Barani-333, Ghauri-3, and Barani-333 × FBG-333 revealed positive correlation with K+/Na+, K+, POD, SOD, CAT, TSP, Chla, Chlb, CAR, LP, FS, FL, FS, NBP, BW, and SCY and declared as salt tolerant (Fig. 9). The genotypes FBS- FALCON × FBG-333, FBS- FALCON × Barani-222, and FBS- SMART × FBG-333 showed significant association with FF under saline conditions. The insights obtained from PCA can be utilized to obtain valuable outcomes (Fig. 9).
Examination of the expression of ERF genes (Ghi-ERF-2D.6 and Ghi-ERF-7A.6) by qPCR
The quantitative measurement of mRNA expression in the leaf of genotypes at salt stress levels of control and 15 dsm− 1 was performed. As an internal control of gene normalization, the GAPDH was employed. In Fig. 10, the genes Ghi-ERF-2D.6 and Ghi-ERF-7A.6 showed relatively higher expression at 15 dsm− 1. In genotype JSQ-white Hold, the expression of Ghi-ERF-2D.6 is 1.5 folds higher at 15 dsm− 1 than control whereas the expression of Ghi-ERF-7A.6 is 2.2 folds higher at 15 dsm− 1 than control. The genotype Barani-333 × FBG-333 showed 2 folds higher expression of Ghi-ERF-2D.6 at 15 dsm− 1 than control whereas the expression of Ghi-ERF-7A.6 is 2.7 folds higher at 15 dsm− 1 than control. The expression of both Ghi-ERF-2D.6 and Ghi-ERF-7A.6 in genotype FBG-222 × Barani-222 was measured 2 folds higher at 15 dsm− 1 than control. In genotype Ghauri, the expression of Ghi-ERF-2D.6 and Ghi-ERF-7A.6 was 2.5 folds and 2.2 folds higher at 15 dsm− 1 than control respectively. The genotype FBG-222 × FBG-333 showed 1.5 folds and 3 folds higher expression of Ghi-ERF-2D.6 and Ghi-ERF-7A.6 at 15 dsm− 1 than control respectively (Fig. 10). In genotype FBS- FALCON × JSQ-white Gold, the expression of Ghi-ERF-2D.6 is 2 folds higher whereas the expression of Ghi-ERF-7A.6 is 2.7 folds higher at 15 dsm− 1 than control respectively. The genotype Barani-333 showed 2.5 folds and 3 folds higher expression of Ghi-ERF-2D.6 and Ghi-ERF-7A.6 at 15 dsm− 1 than control respectively. The expression of Ghi-ERF-2D.6 and Ghi-ERF-7A.6 in genotype FBS- FALCON was measured 2.1 folds and 2.3 folds higher at 15 dsm− 1 than control respectively. The expression level of both the genes was increased at 15 dsm− 1 and showed significant tolerance against salt stress (Fig. 10).