Indian Journal of Agricultural Research

  • Chief EditorV. Geethalakshmi

  • Print ISSN 0367-8245

  • Online ISSN 0976-058X

  • NAAS Rating 5.60

  • SJR 0.217, CiteScore: 0.595

Frequency :
Monthly (January, February, March, April, May, June, July, August, September, October, November, December)
Indexing Services :
BIOSIS Preview, ISI Citation Index, Biological Abstracts, Elsevier (Scopus and Embase), AGRICOLA, Google Scholar, CrossRef, CAB Abstracting Journals, Chemical Abstracts, Indian Science Abstracts, EBSCO Indexing Services, Index Copernicus

Heat Stress Tolerance Indices and Their Contribution in Identifying Heat Tolerant Bottle Gourd [Lagenaria siceraria (Molina) Standl.] Genotypes

Md. Shafiqul Islam1,2, M. Hasanuzzaman1,*, Md. Arifuzzaman1
  • https://orcid.org/0009-0006-4340-9916, https://orcid.org/0000-0003-3037-8509, https://orcid.org/0000-0003-2786-8339
1Department of Genetics and Plant Breeding, Hajee Mohammad Danesh Science and Technology University, Dinajpur, Bangladesh.
2Vegetable Research Center, BRAC Seed and Agro Enterprise, Dinajpur, Bangladesh.

Background: Bottle gourd, a popular winter vegetable in Bangladesh, is experiencing a year-round increase in demand. High temperature stress disrupts physiological structure, fruit quality and crop yields, limiting the year-round cultivation of bottle gourd. The research aimed to identify heat stress-tolerant genotypes to extend growing seasons and increase production among varieties and breeding materials.

Methods: The study conducted at Vegetable Research Center, Dinajpur of BRAC Seed and Agro Enterprise, Bangladesh consisting 58 bottle gourd genotypes under normal and heat stress conditions using augmented block design. Meteorological data showed heat stress resulted ~5-7oC higher air temperatures than normal sowing environments.

Result: Heat stress tolerance indices indicated that the genotypes with a combination of greater values of MP, GMP, HM, YI and STI and lower values of SSI and TOL are high performed and productive breeding material. Yield based stress indices SSI and TOL had a negative association with fruit yield under heat stress (Ys), while strong positive correlations with STI, MP, GMP, HM and YI. The cluster analysis classified bottle gourd genotypes as tolerant, medium tolerant and susceptible to heat stress. Heat stress tolerance indices could be effective tools for identifying heat tolerant bottle gourd genotypes in field. The genotypes BG-45, BG-02, BG-57, BG-39 and BG-43 determined as high performing and stable breeding material under both normal and heat stress environments and might potentially be employed in future breeding programs.

Bottle gourd [Lagenaria siceraria (Molina.) Standl.] is a popular winter vegetable crop grown in Bangladesh, belong to the Cucurbitaceae family (Bose and Som, 1986). Its demand for fresh young fruits, shoots and succulent leaves is increasing year-round for their high water content, fibers, carbs, vitamins, minerals and amino acids (Milind and Satbir, 2011). It have a wide range of genetic diversity, heritability and genetic advance (Sharma and Sengupta, 2012).
       
The global climate has steadily changed during the last century (Zhang et al., 2019). Climate change, global warming and weather uncertainty have raised serious concerns regarding agriculture and crop productivity (Lesk et al., 2016). Climate change is elevating global surface temperatures, inducing heat stress in plants and having a negative impact on their development, growth and yield, altered rainfall patterns and shorter winter seasons (Sultan et al., 2019). Rising global temperatures are threatening human survival by decreasing crop productivity and increasing the global population, which is predicted to exceed nine billion by 2050, resulting in increased food demand (Cole et al., 2018).
               
Cucurbits is a perishable food item, serve various functions in meals but are sensitive to unstable climatic changes. It is predicted to reduce cucurbit production by 15 % due to climate change (Olarewaju et al., 2023). The ideal temperature for bottle gourd growing is 25-30oC (Kumar and Reddy, 2021). High temperatures (>30oC) disrupt cucurbits physiological structure and fruit shapes, crucial for market quality (Saroj and choudhary, 2020). Yield levels are vital role for assessing heat tolerance and growing under normal and heat-stress conditions is a standard method for selecting heat-tolerant genotypes. A number of stress tolerance indices are effective for identifying heat-tolerant bottle gourd genotypes, making them the most credible indicators (Lamba et al., 2023; Poudel et al., 2024). Therefore, present research was carried out to assess the stress indices in identifying high yielding with heat-tolerance bottle gourd genotypes.
The research was carried out at Vegetable research center, Dinajpur of BRAC Seed and Agro Enterprise in the old Himalayan piedmont plain agro-ecological zone of Bangladesh. The study assessed 58 diversified bottle gourd genotypes including registered varieties, advanced breeding lines and local land races collected from Agro-ecological zones of Bangladesh, India, Thailand, China and seed companies (Fig 1). Among the examined genotypes, registered varieties Khetlau (BG-50), BARI Lau-4 (BG-51) and BARI Lau-5 (BG-52), BU Lau-1 (BG-53) and BU Lau-2 (BG-54) were also employed. The experimental genotypes were cultivated for non-stress (winter season) during October 2022-February 2023 and heat stress (summer) in April 2023-August 2023. The genotypes were planted using an augmented block design. Ten plants of each genotype were planted with 2-meter space between rows and plants in each setups. With the exception of sowing time, every other aspect of the two trials (agronomic package and practices) were similar.

Fig 1: Morphological variations of studied 58 bottle gourd genotypes.


 
Data collection
 
Weather data were collected from the Bangladesh Meteorological Department, Dinajpur. The high temperature stress was recorded nearly 5-7oC higher than normal, with mean maximum temperatures reaching 29.66oC in Rabi 2022-23 while 34.62oC in summer, 2023 (Fig 2). Various morphological and yield data were collected from each genotype under both heat stressed and normal environments.

Fig 2: Monthly temperature and relative humidity (%) during non-stress and heat stress period.


 
Statistical analysis
 
MS Excel 2016 were utilized for data entry, processing and stress index calculation. The stress indices, PCA analysis and biplot were analysed using Plant Abiotic Stress Index Calculator (Pour-Aboughadareh et al., 2019). Correlation and cluster analysis were performed using Python using Jupyter 7.0.8.
       
The yield data were utilized in equations for following indices to calculate genotypic mean values for fruit yield under heat stress (Ys) and non-stress (Yp). Xs and Xp represent the average yield of all genotypes under heat stress and non-stress environments, respectively (Porch, 2006).
 
Yield based indices of heat tolerance
 
Yield-based stress indices TOL, MP, GMP, HM, SSI, STI, YI, YSI and RSI were determined for studied bottle gourd genotypes using yield under both normal (Yp) and heat stress (Ys) conditions (Table 1). The genetic variation among genotypes defined by different heat stress indices are depicted in figures (Fig 3-5). Nineteen bottle gourd genotypes possessed higher TOL, SSI, YSI and RSI and were marked as heat-susceptible genotypes because they recorded high yield under normal environments but zero yield when the plant died under heat stress and hence, these genotypes are only appropriate for normal environments.

Table 1: Fruit yield (Yp and Ys t ha-1) and stress tolerance indices, their mean rank and standard deviation of 58 bottle gourd genotypes.



Fig 3: Variability of YI under non-stress and heat stress environments.



Fig 4: Variability of STI under non-stress and heat stress environments.



Fig 5: Variability of YSI under non-stress and heat stress environments.


       
The lowest value of TOL, SSI, YSI and RSI was exhibited in genotypes BG-42 followed by BG-18, BG-35 and BG-39 denoted as tolerant to heat stress. These genotypes performed less under both conditions due to less yield differential, implying that low values do not always indicate high performance and emphasising the relevance of yield as a crucial factor. Lamba et al., (2023) suggested to determined heat tolerant genotypes considering lower TOL, SSPI, YSI and RSI which supported the current findings. Superior genotype selection has also been reported based on lower SSI and TOL in chickpeas (Jha et al., 2018), common bean (Porch, 2006) and bottle gourd (Mashilo et al., 2017). Based on greater values of MP, GMP, HM, STI and YI  and lower TOL and SSI combinations genotype BG-02 were considered to be most stable and productive genotypes followed by BG-43, BG-57, BG-01 and BG-45 identified among all the genotypes under both conditions as heat tolerant. Sofi et al., (2018) suggested classifying genotypes into four groups based on their resilience and productivity in both stress and non-stress environments. These indices were effective in identifying high-yielding cum heat-tolerant, which identical to the findings reported by Poudel et al., (2024) and Jha et al., (2018).
 
Correlation between fruit yield and stress tolerance indices
 
The indices that showed a strong relationship with fruit yield in both environments were chosen as best in selecting genotypes with high yield in both conditions (Table 2). Fruit yield under non-stress (Yp) was significant positive association (0.369) with fruit yield under heat stress (Ys) suggested that these might be utilized to select high yielding bottle gourd genotypes. Fruit yield (Ys) had positive and significant relationship with MP (0.904), GMP (0.983), HM (0.991), STI (0.968), YI (1) and YSI (0.952) but highly negative significant correlation with TOL (-0.796) and SSI (-0.952). Whereas, Yp demonstrated a significant and positive relationship with TOL (0.268), MP (0.732), GMP (0.405), HM (0.396), STI (0.514) and YI (0.369). Fruit yield (Ys) was negatively correlated with TOL (-0.796) while positively correlated (0.268) under normal conditions; thus, selection based on these indices will yield more fruit under normal conditions but less under heat stress (Lamba et al., 2023). Yield (Yp and Ys) were determined positive association with MP, GMP, HM, STI and negatively correlated with TOL and SSI for heat tolerance in wheat (Poudel et al., 2024), chickpeas (Jha et al., 2018) and water stress in cow pea (Gull et al., 2019) which matched with current findings. Kumar et al., (2020) reported that yield (both normal and stress environments) in mungbean was positively associated with all stress tolerance indices except SSI and Superiority measure. Both TOL and SSI had substantial negative associations with all stress tolerance indices, whereas SSI and TOL demonstrated a strong positive connection (0.886). High positive associations with SSI and TOL but negative correlations with other heat stress indices were also reported by Jha et al., (2018). The findings indicated that these parameters were more effective in selecting genotypes with high yields under different environments.

Table 2: Correlation coefficient between fruit yield (Yp and Ys) and stress tolerance indices.


 
Principal component and biplot analysis
 
The first two principal components (PCs) with an Eigen value >1.0 generated variations by 83.89% and 14.50% of PC1 and PC2, respectively, for heat tolerance indices (Table 3) which contributed highest variation (98.39%) of the six components. Seepal et al., (2025), Mashilo et al., (2017) and Farshadfar et al., (2013) reported 97.46%, 99.72% and 99.60% of total variations in field pea, bottle gourd and wheat respectively. PC1 explained significant positive association with Ys, MP, GMP, HM, STI, YI and YSI but had a negative relationship with TOL and SSI. Thus, PC1 regarded as stress tolerance component. PC2 were identified as heat susceptibility component since it explained a greater correlation with TOL, SSI and Yp. This method also utilized by Poudel et al., (2024) and Devi et al., (2021) to categorize the components of wheat under heat stress. Ullah et al., (2022) suggested that most appropriate criteria are to select stable genotypes with low PC2 and greater PC1 values and vice versa. Thus, genotypes BG-45, BG-57, BG-02 and BG-43 were found higher PC1 values with lower PC2 values regarded as superior with stable and tolerant genotypes under both heat stress and non-stress environments and might recommended for cultivation, Lamba et al., (2023) employed this approach to identify stable heat tolerant genotypes.

Table 3: Principal component analysis showing eigenvectors, eigenvalues and variability.


       
The biplot depicted that Yp and Ys were positively connected with MP, GMP, HM, STI, YI and YSI, while Ys was negatively associated with TOL and SSI, as evidenced by the obtuse and acute angles between their vectors respectively (Fig 6). Lamba et al., (2023) and Jha et al., (2018) displayed these types of relationship among indices using biplots for heat tolerance in wheat and chickpeas, respectively. Devi et al., (2021) similarly illustrated positive relationship with MP, GMP and STI as acute angles between their vectors while negative with TOL and SSI as an obtuse angle in Ys.

Fig 6: Biplot displayed relationship among heat stress indices of 58 genotypes.


 
Cluster analysis
 
The cluster analysis categorized genotypes into three clusters based on stress tolerance indices trends. Among the clusters, cluster 2 possessed the twenty-one genotypes with highest MP, GMP, STI, YI and YSI, indicating strong performance and desirable traits categorized as tolerant. Similarly, cluster 3 exhibited medium values of heat tolerance indices and consisted 18 genotypes classified as semi-tolerant. Consequently, cluster 1 had the highest TOL and SSI values but least MP, GMP and YI values, with 19 genotypes suited only for normal condition, indicating that they were susceptible genotypes (Fig 7). Naghavi et al., (2013) categorized genotypes into three clusters for drought tolerance which validated selection through clustering of studied genotypes.

Fig 7: Dendrogram showing classification of genotypes based on heat tolerance.


 
Ranking of genotypes
 
The ranking system was employed to determine the overall performance of genotypes. In order to identify the most heat-tolerant genotype based on all indicators, mean rank (R) and SDR of all indices were estimated. The genotypes BG-02, BG-45 and BG-57 were identified as most heat-tolerant with low rank values based on all heat stress tolerance indices (Table 1). Mashilo et al., (2017) and Naghavi et al., (2013) utilized ranking method for screening corn and bottle gourd genotypes, respectively. 
The correlation coefficient analysis of heat stress indices revealed that MP, GMP, HM and STI had positive association with Yp and Ys, but TOL and SSI had negative correlation with Ys. Thus, these indices might be utilized for identification heat tolerant bottle gourd genotypes. PCA analysis revealed 98.39% of total genetic variation employed in PC1 and PC2, which classified genotypes as tolerant and susceptible, respectively. Cluster analysis illustrated the genotypes more clearly for tolerant, mid-tolerant and susceptible grouped for heat stress conditions. The genotypes BG-45, BG-02, BG-57, BG-39 and BG-43 were determined to be heat tolerant by evaluating their contributions of stress indices by PCA, biplot and ranking. Further, the heat stress tolerance indices TOL, SSI, MP, GMP, STI and YI are the most appropriate indicators to determine heat tolerance genotypes.
The authors are grateful to BRAC Seed and Agro Enterprise, BRAC at 66, Dhaka- 1212, Bangladesh and its authorities for financial support of this study. The work is based on Ph.D research of the first author.
The authors declare no conflict of interest exists.

  1. Bidinger, F.R., Mahalakshmi, V. and Rao, G.D. (1987). Assessment of drought resistance in pearl millet [Pennisetum americanum (L.) Leeke]. I. Factors Affecting Yields under Stress. Australian Journal of Agricultural Research. 38: 49-59.

  2. Bose, T.K. and Som, M.G. (1986). Vegetable Crops in India. Naya Prokash, Calcutta, India. pp: 91-164.

  3. Bouslama, M. and Schapaugh, W.T. (1984). Stress tolerance in soybean. Part 1: Evaluation of three screening techniques for heat and drought tolerance. Crop Science. 24: 933-937. 

  4. Cole, M.B., Augustin, M.A., Robertson, M.J., Manners, J.M. (2018). The science of food security. NPJ Science of Food. 2: 14. 

  5. Devi, K., Chahal, S., Singh, S., Venkatesh, K., Mamrutha, H.M., Raghav, N., Singh, G., Singh, G.P., Tiwari, R. (2021). Assessment of wheat genotypes based on various indices under different heat stress conditions. Indian Journal of Genetics and Plant Breeding. 81(3): 376-382. 

  6. Farshadfar, E., Poursiahbidi, M.M. and Safavi, S.M. (2013). Assessment of drought tolerance in land races of bread wheat based on resistance/ tolerance indices. International Journal of Advanced Biological and Biomedical Research. 1(2): 143-158.

  7. Fernandez, G.C.J. (1992). Effective Selection Criteria for Assessing Plant Stress Tolerance. In: Asian Vegetable Research and Development Center, [Kuo, C.G. (ed.).] Taiwan. pp: 257-270. 

  8. Fischer, R.A. and Maurer, R. (1978). Drought resistance in spring wheat cultivars. I. Grain yield responses. Australian Journal of Agricultural Research. 29: 897-912.

  9. Fischer, R.A. and Wood, T. (1979). Drought resistance in spring wheat cultivars III. Yield association with morphological traits. Australian Journal of Agricultural Research. 30: 1001-1020. 

  10. Gavuzzi, P., Rizza, F., Palumbo, M., Campaline, R.G., Ricciardi, G.L., Borghi, B. (1997). Evaluation of field and laboratory predictors of drought and heat tolerance in winter cereals. Canandian Journal of Plant Science. 77: 523-531.

  11. Gull, M., Sofi, P.A., Mir, R.R., Ars, A., Zargar, S.M. (2019). Productivity and resilience based indices for identification of water stress resilient genotypes in cowpea (Vigna unguiculata L.). Indian Journal of Agricultural Research. 53(4): 391- 397. doi: 10.18805/IJARe.A-5139.

  12. Jha, U.C., Jha, R., Singh, N.R., Shil, S., Kole, P.C. (2018). Heat tolerance indices and their role in selection of heat stress tolerant chickpea (Cicer arietinum) genotypes. Indian Journal of Agricultural Sciences. 88(2): 260-267.

  13. Kumar, R., Singh, C.M., Arya, M., Kumar, R., Mishra, S.B., Singh, U.K., Paswan, S. (2020). Investigating stress indices to discriminate the physiologically efficient heat tolerant genotypes of mungbean [Vigna radiata (L.) Wilczek]. Legume Research. 43(1): 43-49. doi: 10.18805/LR-3950.

  14. Kumar, R. and Reddy, K.M. (2021). Impact of Climate Change on Cucurbitaceous Vegetables in Relation to Increasing Temperature and Drought. Springer, Cham. pp: 175-195.

  15. Lamba, K., Kumar, M., Singh, V., Chaudhary, L., Sharma, R., Yashveer, R., Dalal, M.S. (2023). Heat stress tolerance indices for identification of the heat tolerant wheat genotypes. Scientific Reports. 13: 10842. 

  16. Lesk, C., Rowhani, P. and Ramankutty, N. (2016). Influence of extreme weather disasters on global crop production. Nature: 529(7584): 84-87.

  17. Mashilo, J., Shimelis, H. and Odindo, A. (2017). Yield-based selection indices for drought tolerance evaluation in selected bottle gourd [Lagenaria siceraria (Molina) Standl.] landraces. Acta Agriculturae Scandinavica, Soil and Plant Science. 67(1): 43-50.

  18. Milind, P. and Satbir, K. (2011). Is bottle gourd a natural guard? International Research Journal of Pharmacy. 2: 13-17.

  19. Naghavi, M.R., Pour Aboughadareh, A. and Khalili, M. (2013). Evaluation of drought tolerance indices for screening some of corn (Zea mays L.) cultivars under environmental conditions. Notulae Scientia Biologicae. 5(3): 388-393.

  20. Olarewaju, O.O., Fajinmi, O.O., Arthur, G.D., Coopoosamy, R.M., Naidoo, K. (2023). Effect of climate change on the production of Cucurbitaceae species in North African countries. Journal of Agriculture and Food Research.14: 100-742 

  21. Porch, T.G. (2006). Application of stress indices for heat tolerance screening of common bean. Journal of Agronomy and Crop Science.192: 390-394.

  22. Poudel, R.M., Dhakal, A., Bhandari, R.K., Nyaupane, S., Poudel, S., Panthi, B., Neupane, P. (2024). Evaluation of wheat genotypes based on heat stress tolerance. Annals Agri Bio Research. 29(1): 50-56. 

  23. Pour-Aboughadareh, A., Yousefian, M., Moradkhani, H., Vahed, M.M., Poczai, P., Siddique, K.H.M. (2019). iPASTIC:  An online toolkit to estimate plant abiotic stress indices. Applications in Plant Science. 7(7): e11278.

  24. Ramirez, P.V. and Kelly, J.D. (1998). Traits related to drought resistance in common bean. Euphytica. 99: 127-136.

  25. Rosielle, A.A. and Hamblin, J. (1981). Theoretical aspects of selection for yield in stress and non stress environments. Crop Science. 21: 943-946.

  26. Saroj, P.L. and Choudhary, B.R. (2020). Improvement in cucurbits for drought and heat stress tolerance- A review. Current Horticulture. 8(2): 3-13. 

  27. Seepal, Y.S., Sharma, V., Singh, C.M., Shukla, G., Gangwar, V., Kamaluddin and Singh, S.K. (2025). Application of stress indices to identify terminal heat tolerance genotype in field pea (Pisum sativum var. arvense). Legume Research.  48(1): 20-25. doi: 10.18805/LR-4888.

  28. Sharma, A. and Sengupta, S.K. (2012). Evaluation of genetic variability in bottle gourd genotypes. Vegetable Science. 39(1): 83-85. 

  29. Sofi, P.A., Rehman, K., Ara, A. Gull, M. (2018). Stress tolerance indices based on yield, phenology and biomass partitioning: A review. Agricultural Reviews. 39(4): 292-299. doi: 10.18805/ag.R-1822

  30. Sultan, B., Defrance, D. and Iizumi, T. (2019). Evidence of crop production losses in West Africa due to historical global warming in two crop models. Scientific Reports. 9: 12834.

  31. Ullah, A., Shakeel, A., Ahmed, H.G.M.D., Naeem, M., Ali, M., Shah, A.N., Wang, L., Jaremko, M., Abdelsalam, N.R., Ghareeb, R.Y., Hasan, M.E. (2022). Genetic basis and principal component analysis in cotton (Gossypium hirsutum L.) grown under water deficit condition. Frontiers in Plant Science.13: 981-369. 

  32. Zhang, X., Li, X., Chen, D., Cui, H., Ge, Q. (2019). Overestimated climate warming and climate variability due to spatially homogeneous CO2 in climate modeling over the Northern Hemisphere since the mid-19th century. Scientific Reports. 9: 17426.

Editorial Board

View all (0)