Skip to main content

Drought mitigation in cocoa (Theobroma cacao L.) through developing tolerant hybrids



Cocoa, being a shade loving crop cannot withstand long periods of water stress. Breeding for drought tolerance is the need of the hour due to change in climatic condition and extension of crop to non-traditional areas. Hybrids were produced by crossing four tolerant genotypes in all possible combination. The cross GV1 55 x M 13.12 didn’t yield any fruit due to cross incompatibility between these genotypes. Various biochemical parameters act as the true indicators to select tolerant and susceptible types. The major biochemical parameters considered after imposing stress included proline, nitrate reductase activity, superoxide dismutase content and glycine betaine.


The drought tolerant hybrids were having high amount of proline, superoxide dismutase enzyme and glycine betaine content. Normally, plants having drought stress show low amount of nitrate reductase activity. However, in case of hybrids, the drought tolerant hybrids were having higher NR activity than the susceptible hybrids. The highest amount of NR was found in the control plants kept at fully irrigated conditions.


This experiment showed the role of different biochemical enzymes and osmolytes in giving tolerance to plants during drought stress. Logistic regression analysis selected proline and nitrate reductase as the two biochemical markers for identifying efficient drought tolerant genotypes in the future breeding programmes.

Peer Review reports


Cocoa is a high industrial valued crop that originated in humid tropics of rain forest and produce cocoa beans which is the only source for chocolate. Cocoa plays a significant role in economics of both producing and consuming countries. In the whole world cocoa is cultivated by around five million farmers and contributes to the subsistence of 40–50 million people [1]. The continuous climate change and limitations in availability of water in several regions of the world, creates a serious threat to plantation industry, especially cocoa as the crop is drought sensitive. Even though cocoa is cultivated under rainfed condition in many cocoa growing countries, inconsistency in rainfall pattern forms the major constrain. It is estimated that a minimum of 24 l of water/ 4–5 days interval is required for the plant to express its maximum potential [2]. The main reason behind this is that it has a very shallow tap root system which enables it to absorb water from surface layers only. The growth and yield of cocoa is influenced by a number of environmental factors, particularly rainfall, temperature and water stress [3].

Drought is defined as a decrease in water inputs or precipitation in an agro/ecosystem over time that is sufficient to result in soil water deficit [4]. Drought stress is the most common environmental stress causing threat to successful production of crop plants. There is a growing concern about the global increase in temperature and simultaneous increase in potential evapo- transpiration [5]. Plant water demand may result in increased drought stress during the day and a further deterioration of climatic condition for cocoa. Drought stress normally decreases plant growth by stressing various physiological and biochemical processes in plant, so these parameters can be used as an indicator of stress responses in plants. The plant response to drought stress are usually studied based on physiological parameters, but in recent studies biochemical and antioxidant responses also been proposed as appropriate indicators of drought stress in plants [6, 7].

Moreover, due to the increased rate of consumption of cocoa in chocolate industry, the cultivation is extending to non-traditional area. In such area acute water shortage is the most important problem to be addressed. Hence to cope up with the changing climatic scenario and cropping pattern it is essential to breed genotypes that can tolerate water stress. The screening of genotypes based on biochemical parameters helps in developing drought tolerant hybrids which can be used in further breeding programmes. The objective of the present study is to screen and identify cocoa genotypes which can tolerate drought and the role of biochemical parameters in drought stress so as to develop cocoa varieties suitable for drought prone areas.

Materials and methods

Experimental site

The hybridization work was conducted at Cocoa Research Centre (CRC), Kerala Agricultural University, India, using four genotypes identified as tolerant to drought in a preliminary study [8] as parents (Table 1). They were crossed in diallel out of which, the cross GV1 55 x M 13.12 yielded no fruits due to cross incompatibility. The pods matured approximately within 5–6 months; mature pods from each cross were harvested separately and raised in the nursery. Since parents are in heterozygous condition F1 population itself was segregating. Hence each individual was considered as separate hybrids.

Table 1 List of parents used for hybridisation

Drought stress imposition and management

This study was carried out during 2019–20 in a research greenhouse at Cocoa Research Centre, Vellanikkara located in Thrissur (10° 32′ N, 76° 17′ E). Based on initial vigour, hybrids were selected and subjected to moisture stress following gravimetric method [9]. Seedlings of 5 month old with an average of 15 leaves were selected for drought screening in green house. Each progeny is replicated three times and used in the experiment. Initial trials in CRC indicated that cocoa cannot tolerate water stress less than 40% field capacity [10]. Hence, 40% field capacity was maintained for 2 weeks.

A control was also kept at fully irrigated condition representing each of the crosses. Based on the percentage of leaves retained after the drought imposition for 2 weeks, morphological classification of these hybrids was done based on the score chart (Table 2). The humidity and the temperature of the mist chamber were recorded using Berlin’s psychrometer on daily basis. Percent of leaves retained was calculated using the formula.

$$\mathrm{Percentage}\ \mathrm{of}\ \mathrm{leaves}\ \mathrm{retained}=\frac{\mathrm{Number}\ \mathrm{of}\ \mathrm{leaves}\ \mathrm{retained}}{\mathrm{Total}\ \mathrm{number}\ \mathrm{of}\ \mathrm{leaves}}\times 100$$
Table 2 Score chart depicting the leaves retained in the hybrids

Analysis of biochemical parameters

Biochemical analysis was carried out using standard procedures after 2 weeks of stress imposition. The parameters considered were proline (μg/g) [11], Nitrate reductase activity (mmol nitrate/g/hr) [12], superoxide dismutase (units/mg protein/g) [13] and glycine betaine (μmol/g) [14]. Analysis of variance was done for biochemical analysis for all selected hybrids following completely randomised design (CRD).

Corelation and path coefficient analysis

The correlation coefficients were calculated to determine the degree of association of characters with percentage of leaves retained. The genotypic coefficients of correlation between character pairs were determined by using the variance and covariance components [15]. The direct and indirect effects of various biochemical parameters was estimated by path coefficient analysis by using the simple correlation coefficient [16, 17] and the direct and indirect effects were grouped into very high (> 1.00), high (0.30–0.99), medium (0.20–0.29), low (0.10–0.19) and negligible (0.09–0.00) [18].

Binary logistic regression analysis

Binary logistic regression was carried out to find the relationship between the dependent biochemical parameters over the independent variable, the number of leaves retained. Improvement in selection over the base population was found by the equation,

$$\mathrm{Per}\ \mathrm{cent}\ \mathrm{improvement}\ \mathrm{over}\ \mathrm{base}\ \mathrm{population}=\frac{\mathrm{Exp}\left(\mathrm{B}\right)}{1+\mathrm{Exp}\ \left(\mathrm{B}\right)}\times 100$$


Selection of hybrids and classification

A total of 1505 hybrid seedlings were raised in the nursery. Based on initial vigour 120 hybrids were selected representing 11 crosses. After imposition of water stress for 2 weeks the percentage of leaves withered was calculated and hybrids were classified as highly tolerant (HT), tolerant (T), susceptible (S) and highly susceptible (HS) as per the score chart given. Table 2. Different biochemical characters recorded from 120 hybrids are statistically analysed and depicted in Table 3.

Table 3 Biochemical parameters of cocoa hybrid seedlings

Effect of water stress on proline content

In all the 11 crosses, all the tolerant and highly tolerant hybrids indicated high proline content as compared to the susceptible hybrids. The control recorded minimum amount of proline under full irrigated conditions (Table 3).

In the cross M 13.12 x G I 5.9, the highest content was found in hybrid H27 (1105.64 μg/g) and the lowest value was found in H18 (167.18 μg/g). The cross M 13.12 x G II 19.5 had ten hybrids out of which H38 showed the maximum content of proline of about 2710.82 μg/g. Lowest value was found in H40 (224.46 μg/g). In cross M 13.12 x G VI 55, the highest content was found in H43 (2817.39 μg/g) which was a tolerant hybrid and lowest in H44 of only about 85.52 μg/g of proline indicating its vulnerability to drought stress. The cross G I 5.9 x M 13.12 had six hybrids where H52 had highest proline content of 536.70 μg/g and H54 (181.17 μg/g) was having lower proline content. In the cross G I 5.9 x G II 19.5, out of 13 hybrids, H63 had the highest proline content of 706.01 μg/g and the lowest was found in H66 (177.26 μg/g). The cross G I 5.9 x G VI 55 had eight hybrids, and the highest content was found in H71 (1749.05 μg/g) while the susceptible hybrid H77 (169.54 μg/g) had the lowest amount. In cross G II 19.5 x M 13.12, the highly tolerant hybrid H85 (2293.88 μg/g) was having the highest amount while the lowest values were found in H91 (139.87 μg/g) which was a susceptible hybrid. In the cross G II 19.5 x G I 5.9, high values for proline was observed in H101 (2726.81 μg/g) and the lowest in H98 (169.18 μg/g) which was classified as susceptible. In the cross G II 19.5 x G VI 55), only two hybrids were there, one of which was highly tolerant and the other one was susceptible. The H107 was having the proline content of 1984.83 μg/g and the susceptible one, H108 was having 339.69 μg/g of proline. The cross G VI 55 x G I 5.9 was having six hybrids and H113 was having the highest value of proline, 520.85 μg/g and lowest value was found in H114 (128.55 μg/g). In the cross G VI 55 x G II 19.5, highest value was observed in H118 (1354.74 μg/g) while the lowest value was observed in H116 which was susceptible to drought having only 182.50 μg/g of proline.

When the progenies of all crosses were compared, it was seen that the content ranged from 85.52 μg/g in H44 (M 13.12 x G VI 55) to 2817.39 μg/g in H43 (M 13.12 x G VI 55).

Effect of water stress on nitrate reductase activity (NRA)

Tolerant hybrids had high resistance to drought stress and were able to regulate the nitrate reduction activity even with less water. In this experiment, the cross G VI 55 x G II 19.5 recorded the maximum amount of enzyme activity of about 17.06 mmol nitrate/g/hr. in H115 which is a tolerant hybrid. However, the susceptible hybrids had very low content of NRA, the lowest being in H17 (0.51 mmol nitrate/g/hr) from the cross M13.12 x G I 5.9. The control plant on the other hand, had the highest value of NRA which was kept under 100% field capacity.

In the cross between M 13.12 x G I 5.9, 32 hybrids were obtained out of which the H11 showed the maximum nitrate reductase activity of about 12.59 mmol nitrate/g/hr. while the lowest value of nitrate reductase activity was observed in the susceptible hybrid H17.

The cross M 13.12 x G II 19.5 had hybrids having values as high as 11.74 mmol nitrate/g/hr. in H37 to as low as 2.07 mmol nitrate/g/hr. in H36.In cross M 13.12 x G VI 55, H48 was having the highest NRA content of about 10.89 mmol nitrate/g/hr. The lowest value was found in H50 with only 1.88 mmol nitrate/g/hr. In the cross between G I 5.9 x M 13.12, highest NRA value was found in H52 (8.45 mmol nitrate/g/hr) which was a tolerant hybrid while the lowest value was found in H54 (1.81 mmol nitrate/g/hr). The cross G I 5.9 x G II 19.5 had H64 with the highest NRA value of 9.48 mmol nitrate/g/hr. while the lowest value for NRA was 1.41 mmol nitrate/g/hr. found in H59.In the cross between G I 5.9 x G VI 55, the highest NR activity was found in H76 with 10.41 mmol nitrate/g/hr. which is a highly tolerant hybrid. The lowest NR activity was found in H73 having 3.65 mmol nitrate/g/hr. which is susceptible hybrid. In cross G II 19.5 x M 13.12, H81 showed the maximum NR activity of 14.27 mmol nitrate/g/hr. and the lowest value (3.98 mmol nitrate/g/hr) was observed inthe susceptible hybrid H89.The cross G II 19.5 x G I 5.9 had hybrid having the highest value in NR activity (H97) with a value of 15.45 mmol nitrate/g/hr. and the lowest (4.76 mmol nitrate/g/hr) was reported in H98.The cross G II 19.5 x G VI 55 had two hybrids, the high NR value was observed in the tolerant hybrid H107 (12.98 mmol nitrate/g/hr) and the low value of NR activity was observed in H108 (6.65 mmol nitrate/g/hr). In the cross G VI 55 x G I 5.9, the tolerant hybrids were having the high NR activity, the highest being H111 (6.68 mmol nitrate/g/hr), and lowest in H110 (3.57 mmol nitrate/g/hr). The cross G VI 55 x G II 19.5 had values as high as 17.06 mmol nitrate/g/hr. in H115 while the lowest value was observed in H116 having 4.10 mmol nitrate/g/hr. nitrate activity.

Effect of water stress on superoxide dismutase

In cross M13.12 x G I 5.9, the highest value was observed in H10 (0.381 units/mg protein/g) while the lowest value was observed in H16 (0.128 units/mg protein/g). In the cross M 13.12 x G II 19.5, highest SOD value was found in H38 (0.332 units/mg protein/g) which is a highly tolerant hybrid while the lowest was H41 (0.184 units/mg protein/g). In the cross M 13.12 x G VI 55, the highest SOD value was found in H46 (0.362 units/mg protein/g). The lowest SOD value was in hybrid H45 (0.166 units/mg protein/g) which was highly susceptible. The cross G I 5.9 x M 13.12 had the highest SOD values of about 0.259 units/mg protein/g expressed in hybrid H52 while the lowest values was found in H55 (0.156 units/mg protein/g). In cross G I 5.9 x G II 19.5, the highest value was present in hybrid H68 (0.276 units/mg protein/g). The lowest value obtained was in hybrid H60 (0.162 units/mg protein/g). The cross G I 5.9 x G VI 55 had high SOD activity in hybrid H75 (0.272 units/mg protein/g) whereas lowest SOD value was found in H77 (0.191 units/mg protein/g). In cross G II 19.5 x M 13.12, SOD values was highest in H85 (0.227 units/mg protein/g). The lowest SOD value was found in H93 (0.133 units/mg protein/g). In the cross G II 19.5 x G I 5.9, hybrid H101 (0.270 units/mg protein/g) had the highest value and the lowest was found in H106 (0.154 units/mg protein/g). In the cross G II 19.5 x G VI 55, two hybrids were obtained in which the hybrid H107 was highly tolerant hybrid having high SOD content of 0.163 units/mg protein/g whereas the susceptible hybrid was having 0.106 units/mg protein/g of SOD activity. In the cross G VI 55 x G I 5.9, high SOD value was found in hybridH112 (0.221 units/mg protein/g). The low SOD value was found in hybrid H110 (0.136 units/mg protein/g). In cross G VI 55 x G II 19.5, the highest value was found in hybrid H118 (0.227 units/mg protein/g) while the lowest value was found in hybrid H116 (0.140 units/mg protein/g) which is a susceptible hybrid.

H10 (M13.12 x G I 5.9) recorded the high SOD content of about 0.381 units/mg protein/g which was a tolerant hybrid and the lowest content was observed in a susceptible hybrid, H108 (G II 19.5 x G VI 55) of about 0.106 units/mg protein/g.

Effect of water stress on glycine betaine (GB)

When the hybrids of present study were analysed, the tolerant hybrid, H103 (G II 19.5 x G I 5.9) was having high amount of glycine betaine of about 13.79 μmol/g and the susceptible hybrid, H3 (4.42 μmol/g) of the cross M13.12 x G I 5.9 recorded the least. All the tolerant hybrids had high glycine betaine content as compared to the susceptible hybrids. However, the control which was kept under fully irrigated condition had least amount of glycine betaine among the crosses indicating the accumulation of glycine betaine under drought stress conditions. In cross M 13.12 x G I 5.9, the highest glycine betaine value was found in H26 (11.64 μmol/g) and the lowest value was observed in H3 (4.42 μmol/g). In cross M 13.12 x G II 19.5, the highest value was found in hybrid H34 (9.68 μmol/g) whereas the lowest value was observed in H35 (6.73 μmol/g). In cross M13.12 x G VI 55 the highest value was observed in H48 (9.87 μmol/g) and the H45 (5.29 μmol/g) recorded the lowest glycine betaine value.

In cross G I 5.9 x M 13.12, high GB values were observed in H52 (8.10 μmol/g) while H55 (5.36 μmol/g) had the lowest value. In cross G I 5.9 x G II 19.5, the highest value was observed in H63 (9.52 μmol/g) to as low as 6.45 μmol/g in H70.In cross G I 5.9 x G VI 55, the highest value was observed in hybrid H71 (9.56 μmol/g). Lowest values were observed in H73 (7.07 μmol/g). In cross G II 19.5 x M 13.12, highest GB value was observed in H88 (10.08 μmol/g) while the lowest value was observed in H80 (6.78 μmol/g). In the cross between G II 19.5 x G I 5.9, H103 (13.79 μmol/g) recorded the highest value while the lowest value was observed in H99 (7.20 μmol/g). The cross G II 19.5 x G VI 55 had two hybrids, the hybrid H107 had high GB value of 9.34 μmol/g and the other hybrid, H108 which is the susceptible one, had 7.18 μmol/g of GB. In cross G VI 55 x G I 5.9, H112 was having the highest GB content of about 10.61 μmol/g while the lowest value was observed in H109 (6.93 μmol/g). In the cross G VI 55 x G II 19.5, the highest value was observed in hybrid H119 (9.23 μmol/g) while the lowest value was observed in H116 (4.42 μmol/g).

Corelation and path coefficient analysis

All the biochemical parameters estimated showed positive corelation with dependant variable (percentage of leaves retained). Among the biochemical characters, proline showed maximum correlation with the dependent variable, percentage of leaves retained (0.777) (Table 4).

Table 4 Correlation among drought tolerant contributing characters of hybrids

High positive effect on the percentage of leaves retained in path analysis result is a direct measure of plant’s tolerance to drought. It was expressed by proline (0.386) and low direct effect was shown by nitrate reductase activity (0.166) (Fig. 1). Negligible direct effects were found in superoxide dismutase (0.063), and negative but negligible effects were shown by glycine betaine (−0.016).

Fig. 1
figure 1

Path diagram for the biochemical characters having direct effects on percent of leaves retained

Binary logistic regression analysis

The positive and comparable value of odds ratio Exp (B) and positive correlation indicated that proline content and NRA content had a positive correlation with drought tolerance and also these characters expressed a significant value of less than 0.05 which is the constant indicating the 95% accuracy with the results.

Based on the Exp (B) value from regression model, expressed percentage for drought tolerance over the base population was calculated and it was found that if selection is based on proline, a new population from the base population will express 51.55% improvement regarding the tolerance. In case of NRA, a new population will show 87.48% improvementover the base population (Table 5).

Table 5 Logistic estimate of characters influencing drought tolerance in cocoa


Drought stress is one of the important problems which adversely affects the performance of different crop plants worldwide. Various strategies are being adopted to cop up with adverse conditions. Development of drought tolerant cultivars is prerequisite to encounter the existing drought stress on sustainable basis. In the present study selection of hybrids were done based on initial vigor and stress imposition was carried out. In cocoa initial vigour is directly correlated to yield [19].

In general, plants employ various methods to promote tolerance to drought stress. To assess the degree of tolerance in plants, various drought-related indicators are used, such as morphological, physiological and biochemical parameters [20]. In the present study evaluation and screening of genotypes under drought stress were done based on biochemical parameters.

Proline is an important amino acid found in proteins. It has a significant contribution in drought tolerance. Under drought stress, a higher accumulation of proline was recorded in the study, which might be due to osmotic adjustment. Even within a single cross, highly tolerant as well as susceptible hybrids were found. This is because parents are heterozygous in nature [21]. The analysis of proline clearly indicated that resistant plants were high in proline content compared to susceptible ones. A clear-cut difference was observed between highly tolerant and susceptible genotypes indicating that proline plays an important role in drought tolerance of cocoa. High levels of proline enabled the plant to maintain low water potentials. Apart from acting as an osmolyte for osmotic adjustment, proline contributes to stabilizing sub-cellular structures (eg., membranes and proteins), scavenging free radicals and buffering cellular redox potential under stress conditions [22]. Similar studies were carried out in other crops. Under drought stress, there was a progressive increase in free proline in cotton plants [23] and a notable proline accumulation was observed in sugar beet leaves when drought stress became severe [24]. It has been already proved that proline is having direct correlation with drought stress [25, 26] as indicated in the present study. Proline in general was known to correlate with stress tolerance and has a direct effect on tolerance capability of plant [22] supporting our present findings.

Another important biochemical parameter contributes to drought stress is nitrate reductase activity. Water deficit induces an abrupt reduction in the uptake and nitrate flux rates from roots to leaves, preventing the mechanisms of NR protein synthesis induction and NR activity. The NR activity decline during water stress is mainly attributed to low NO3 - absorption and availability resulting from water uptake deprivation [27]. The present results indicated that the tolerant hybrids had high resistance to drought stress and were able to regulate the nitrate reduction activity even with less water. Usually, drought stress reduces the enzyme activity and that is the reason the amount of reductase enzyme was low in hybrids whereas in the control, it was more. NRase is closely associated with plant growth and development [28]. It is generally accepted that drought stress has a negative impact on plant’s photosynthetic activity, N concentrations, free amino acids or soluble protein contents accompanied with a decline of nitrate reductase activity in many plant species, such as maize [29], potato [30] and winter wheat [31]. The plants subjected to water stress produces less amount of total protein which results in a decrease in the synthesis of nitrate reductase activity caused by low nitrate flux [32]. The amount of NR enzyme generally decreases during drought stress and hence, the hybrids having more NR enzyme were more tolerant to drought stress and were able to regulate the nitrogen assimilation in plants [33]. Therefore, in the present study NR was found to be directly related to the dependent variable.

Reactive Oxygen Species (ROS) accumulation during stress greatly depends on the balance between ROS production and ROS scavenging mechanism [34]. When plants are subjected to any kind of stress, the cells have an increased production of reactive oxygen species (ROS) which in normal cases, is removed from time to time. Under stress, these become high in number and results in oxidative damage. These are removed by anti-oxidant systems which form the first line of defence which is superoxide dismutase. ROS-scavenging mechanisms were shown to have an important role in protecting plants against osmotic stresses [35].

In the experiment when plants were subjected to analysis, the tolerant and highly tolerant hybrids showed more amount of superoxide dismutase enzyme as compared to the susceptible hybrids and the control which was kept under fully irrigated condition had the least amount of SOD in it. This indicated that SOD will get accumulated under drought stress conditions and forms a defence system against the stress [36]. Superoxide dismutase (SOD) are the enzymes that forms the first line defense and catalyses the dismutation of O2- radicals to H2 O and O2. Hence, the amount ofSOD increases with increase in stress conditions indicating it’sdirect relation to drought stress [37, 38]. Similar studies on development of water stress in Curtilobumsolanum and Solanumtuberosum due to over production of SOD in chloroplasts were also reported [39] which supported our findings.

Many plants accumulate compounds, termed compatible solutes, to cope with stress conditions. One of the most extensively studied compatible solutes is glycine betaine [40] Not only GB acts as an osmoregulator, but also stabilizes the structures and activities of enzymes and protein complexes, and maintains the integrity of membranes against the damaging effects of stress [41]. The role of glycine betaine to drought tolerance has been reported in many cases [42, 43]. Genes associated with glycine betaine synthesis in higher plants and microbes have been transferred into plants which do not accumulate glycine betaine, such as Arabidopsis thaliana [44], Brassica napus [45] Persimmon [46] and rice [47]. The metabolic engineering of glycine betaine biosynthesis in these plants improved the tolerance of transgenic plants to salt, drought and extreme temperature stresses [48]. When the levels of glycine betaine was correlated with the extent of increased tolerance, the accumulation of glycine betaine was found to be induced under stress conditions [49, 50].

Hence, these biochemical parameters were reliable source to detect drought stress in cocoa. However, the correlation and path analysis revealed that proline and the nitrate reductase are the two important parameter that contribute to drought tolerance. To find out the actual relationship between the dependent variables (the biochemical parameters) and the independent variable (the number of leaves retained, a binary regression analysis was carried out. The results of present study indicated that the proline and nitrate reductase can be used as biochemical indicators in cocoa drought breeding programme (Fig. 2).

Fig. 2
figure 2

Biochemical phenes and their association with drought tolerance in cocoa


The parameters proline and glycine betaine represented the osmolyte group whereas nitrate reductase and superoxide dismutase represented the enzyme group. In all the crosses, the content of proline was high in highly tolerant and tolerant hybrids as compared to susceptible hybrids. The control which was fully irrigated condition was having the least amount of proline. The glycine betaine also followed the same trend, as tolerant and highly tolerant hybrids were having more amount of glycine betaine as compared to susceptible hybrids and the control indicating that these two osmolytes accumulated only during water stress. In case of superoxide dismutase, the highly tolerant and tolerant hybrids were having high amount of superoxide dismutase as compared to the susceptible hybrids whereas the control was having the least amount of superoxide dismutase enzyme. In case of nitrate reductase activity, the highly tolerant and tolerant hybrids were having high amount of enzyme as compared to the susceptible hybrids. The control kept at fully irrigated condition was having the highest amount of nitrate reductase enzyme. The hybrids having high amount of nitrate reductase were more tolerant because generally, this enzyme reduces under drought stress. Based on corelation, path analysis and binary logistic regression analysis, proline and nitrate reductase can be used as reliable parameters in drought breeding programmes in the future.

Availability of data and materials

All data generated or analysed during this study are included in this article and its supplementary files.



Nitrate reductase


Nitrate reductase activity


Glycine betaine


Super oxide dismutase


Highly tolerant






Highly susceptible


Reactive oxygen species


Cocoa Research Centre


Completely randomised design


  1. 1.

    WCF. Cocoa market update. World Cocoa Foundation; 2012. Available at:

  2. 2.

    Prasannakumari AS, Vikraman Nair R, Lalithabai EK, Mallika VK, Minimol JS, Abraham K, et al. Cocoa in India. Kerala: Kerala Agricultural University; 2009. p. 72.

    Google Scholar 

  3. 3.

    Minimol JS, Suma B, Shija TK, Shilpa KS. Genotypic and seasonal variations affecting yield of cocoa. J Agrometeorol. 2020;22:528–31.

    Google Scholar 

  4. 4.

    Gilbert ME, Medina V. Drought adaptation mechanisms should guide experimental design. Trends Plant Sci. 2016;21:639–47.

    CAS  PubMed  Google Scholar 

  5. 5.

    Laderach P, Martinez-Valle A, Schroth G, Castro N. Predicting the future climatic suitability for cocoa farming of the world’s leading producer countries, Ghana and Cote d’Ivoire. Clim Chang. 2013;119:841–54.

    Google Scholar 

  6. 6.

    Huang S, Van Aken O, Schwarzlander M, Belt K, Millar AH. The roles of mitochondrial reactive oxygen species in cellular signaling and stress response in plants. Plant Physiol. 2016;171:1551–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Martinez V, Nieves-Cordones M, Lopez-Delacalle M, Rodenas R, Mestre TC, Garcia Sanchez F, et al. Tolerance to stress combination in tomato plants: new insights in the protective role of melatonin. Molecules. 2018;23:535–40.

    PubMed Central  Google Scholar 

  8. 8.

    Binimol B. Identification of drought tolerant cocoa types. Thrissur: M.Sc. (Ag.) Thesis, Kerala Agricultural University; 2005. p. 95.

    Google Scholar 

  9. 9.

    Souza CC, Oliveira FA, Silva IF, Amorim Neto MS. Evaluation of methods of available water determination and irrigation management in “Terra Roxa” under cotton crop. Rev Bras Eng Agric. 2000;4:338–42.

    Google Scholar 

  10. 10.

    CCRP [Cadbury-KAU Co-operative Cocoa Research Project]. Annual report 2014–2015. Vellanikkara: Cocoa Research Centre; 2015. p. 153.

    Google Scholar 

  11. 11.

    Bates L, Waldren RP, Teare ID. Rapid determination of free proline for water-stress studies. Plant Soil. 1973;39:205–7.

    CAS  Google Scholar 

  12. 12.

    Evans HJ, Nason A. Pyridine nucleotide-nitrate reductase from extracts of higher plants. Plant Physiol. 1953;28:233.

    CAS  PubMed  PubMed Central  Google Scholar 

  13. 13.

    Dhindsa RS, Plumb DP, Thorpe TA. Leaf senescence: correlated with increased levels of membrane permeability and lipid peroxidation, and decreased levels of superoxide dismutase and catalase. J Exp Bot. 1981;32:93–101.

    CAS  Google Scholar 

  14. 14.

    Grieve CM, Grattan SR. Rapid assay for determination of water-soluble quaternary amino compounds. Plant Soil. 1983;70:303–7.

    CAS  Google Scholar 

  15. 15.

    Al-Jabouri RA, Miller PA, Robinson HF. Genotypic and environmental variance in upland cotton cross of interspecific origin. Agron J. 1958;50:633–7.

    Google Scholar 

  16. 16.

    Wright S. Correlation and causation. J Agric Res. 1921;20:557–85.

    Google Scholar 

  17. 17.

    Dewey DI, Lu KH. A correlation and path-coefficient analysis of components of crested wheatgrass seed production. Agron J. 1959;51:515–8.

    Google Scholar 

  18. 18.

    Lenka D, Mishra B. Path coefficient analysis of yield in rice varieties. Indian J Agric Sci. 1973;43:376–9.

    Google Scholar 

  19. 19.

    Enriquez GA. Early selection for vigour of hybrid seedlings. In: Proceedingsof the 8th international symposium on cocoa research; 1981. p. 535–9.

    Google Scholar 

  20. 20.

    Rampino P, Pataleo S, Gerardi C, Mita G, Perrotta C. Drought stress response in wheat: physiological and molecular analysis of resistant and sensitive genotypes. Plant Cell Environ. 2006;29:2143–52.

    CAS  PubMed  Google Scholar 

  21. 21.

    Minimol JS, Amma SP. Self incompatibility assessment in vascular streak dieback (VSD) disease resistant hybrids of cocoa (Theobrama cocoa L.). Asian J Hort. 2013;8:114–7.

    Google Scholar 

  22. 22.

    Ashraf M, Foolad MR. Roles of glycine betaine and proline in improving plant abiotic stress resistance. Environ Exp Bot. 2007;59:206–16.

    CAS  Google Scholar 

  23. 23.

    Yi XP, Zhanga YL, Yao HS, Luo HH, Gou L, Chow WS, et al. Rapid recovery of photosynthetic rate following soil water defcit and re-watering in cotton plants (Gossypium herbaceum L.) is related to the stability of the photosystems. J Plant Physiol. 2016;193:23–34.

    Google Scholar 

  24. 24.

    Monreal JA, Jimenez ET, Remesal E. Proline content of sugar beet storage roots: response to water defcit and nitrogen fertilization at feld conditions. Environ Exp Bot. 2007;60:257–67.

    CAS  Google Scholar 

  25. 25.

    Mali PC, Mehta SL. Effect of drought on enzymes and free proline in rice varieties. Phytochem. 1977;16:1355–7.

    CAS  Google Scholar 

  26. 26.

    Karamanos AJ, Drossopoulos JB, Niavis CA. Free proline accumulation during development of two wheat cultivars with water stress. J Agric Sci. 1983;100:429–39.

    CAS  Google Scholar 

  27. 27.

    Foyer CH, Valadier MH, Migge A, Becker TH. Drought-induced effects on nitrate reductase activity and mRNA and on the coordination of nitrogen and carbon metabolism in maize leaves. Plant Physiol. 1998;117:283292.

    Google Scholar 

  28. 28.

    Ferrario-Mery S, Valadier MH, Foyer CH. Over-expression of nitrate reductase in tobacco delays drought-induced decreases in nitrate reductase activity and mRNA. Plant Physiol. 1998;117:293–302.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Sinha SK, Nicholas JD. Nitrate reductase in relation to drought. In: Palag LG, Aspinal LD, editors. Physiology and biochemistry of drought resistance. Sydney: Academic Press; 1981. p. 145.

    Google Scholar 

  30. 30.

    Ghosh SC, Asanuma K, Kusutani K, Toyota M. Effects of moisture stress at different growth stages on the amount of total non-structural carbohydrate, nitrate reductase activity and yield of potato. Jpn J Trop Agr. 2000;44:158–66.

    Google Scholar 

  31. 31.

    Xu ZZ, Yu ZW. Nitrogen metabolism in flag leaf and grain of wheat in response to irrigation regimes. J Plant Nutr Soil Sci. 2006;169:118–26.

    CAS  Google Scholar 

  32. 32.

    Costa RCL, Lobato AKS, Oliveira Neto CF, Maia PSP, Alves GAR, Laughinghouse HD. Biochemical and physiological responses in two Vignauguiculata(L.) walp cultivars under water stress. J Agron. 2008;7:98–101.

    Google Scholar 

  33. 33.

    Xu ZZ, Zhou GS. Effects of water stress on photosynthesis and nitrogen metabolism in vegetative and reproductive shoots of Leymus chinensis. Photosynthetica. 2005;43:29–35.

    Google Scholar 

  34. 34.

    Mittler R, Vanderauwera S, Gollery M, Breusegem VF. Reactive oxygen gene network of plants. Trends Plant Sci. 2004;9:490–8.

    CAS  PubMed  Google Scholar 

  35. 35.

    Abbasi AR, Hajirezaei M, Hofius D, Sonnewald U, Voll L, M. Specific roles of alpha- and gamma-tocopherol in abiotic stress responses of transgenic tobacco. Plant Physiol. 2007;143:1720–38.

    CAS  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Rizhsky L, Liang H, Shuman J, Shulaev V, Davletova S, Mittler R. When defensepathways collide, the response of Arabidopsis to a combination of drought and heat stress. Plant Physiol. 2004;134:1683–96.

    CAS  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Bowler C, Montagu MV, Inze D. Superoxide dismutase and stress tolerance. Annu Rev Plant Physiol Plant Mol Biol. 1992;43:83–116.

    CAS  Google Scholar 

  38. 38.

    Luna CM, Pastori GM, Driscoll S, Groten K, Bernard S, Foyer CH. Drought controls on H2O2 accumulation, catalase (CAT) activity and CAT gene expression in wheat. J Exp Bot. 2004;56:417–23.

    PubMed  Google Scholar 

  39. 39.

    Martinez CA, Loureiro ME, Oliva MA, Maestri M. Differential responses of superoxide dismutase in freezing resistant Solanumcurtilobum and freezing Solanumsensitive tuberosum subjected to oxidative and water stress. Plant Sci. 2001;160:505–15.

    CAS  PubMed  Google Scholar 

  40. 40.

    Rathinasabapathi B. Metabolic engineering for stress tolerance: installing osmoprotectants synthesis pathways. Ann Bot. 2000;86:709–16.

    CAS  Google Scholar 

  41. 41.

    Sakamoto A, Murata N. The role of glycine betaine in the protection of plants from stress: clues from transgenic plants. Plant Cell Environ. 2002;25:163–71.

    CAS  PubMed  Google Scholar 

  42. 42.

    Robinson SP, Jones GP. Accumulation of glycine betaine in chloroplasts provides osmotic adjustment during salt stress. Funct Plant Biol. 1986;13:659–68.

    CAS  Google Scholar 

  43. 43.

    Genard H, Saos J, Billard JP, Tremolieres A, Boucaud J. Effect of salinity on lipid composition, glycine betaine content and photosynthetic activity in chloroplasts of Suaeda maritima. Plant Physiol Biochem. 1991;29:421–7.

    CAS  Google Scholar 

  44. 44.

    Alia H, Sakamoto A, Murata N. Enhancement of the tolerance of Arabidopsis to high temperatures by genetic engineering of the synthesis of glycine betaine. Plant J. 1998;16:155–61.

    CAS  PubMed  Google Scholar 

  45. 45.

    Huang J, Hirji R, Adam L, Rozwadowski KL, Hammerlind JK, Keller WA, et al. Genetic engineering of glycine betaine production towards enhancing stress tolerance in plants: metabolic limitations. Plant Physiol. 2000;122:747–56.

    CAS  PubMed  PubMed Central  Google Scholar 

  46. 46.

    Gao M, Sakamoto A, Miura K, Murata N, Sugiura A, Tao R. Transformation of Japanese persimmon (Diospyros kaki Thunb.) with a bacterial gene for choline oxidase. Mol Breed. 2000;6:501–10.

    CAS  Google Scholar 

  47. 47.

    Mohanty A, Kathuria H, Ferjani A, Sakamoto A, Mohanty P, Murata N, et al. Transgenics of an elite indica rice variety Pusa Basmati 1 harbouring the coda gene are highly tolerant to salt stress. Theor Appl Gen. 2002;106:51–7.

    CAS  Google Scholar 

  48. 48.

    Sulpice R, Tsukaya H, Nonaka H, Mustardy L, Chen TH, Murata N. Enhanced formation of flowers in salt-stressed Arabidopsis after genetic engineering of the synthesis of glycine betaine. Plant J. 2003;36:165–76.

    CAS  PubMed  Google Scholar 

  49. 49.

    Yang WJ, Rich PJ, Axtell JD, Wood KV, Bonham CC, Ejeta G, et al. Genotypic variation for glycine betaine in sorghum. Crop Sci. 2003;43:162–9.

    CAS  Google Scholar 

  50. 50.

    Park EJ, Jeknic Z, Sakamoto A, DeNoma J, Yuwansiri R, Murata N, et al. Genetic engineering of glycine betaine synthesis in tomato protects seeds, plants, and flowers from chilling damage. Plant J. 2004;40:474–87.

    CAS  PubMed  Google Scholar 

Download references


The author acknowledge Kerala Agricultural University and Mondelez International (Earlier Cadbury India Ltd).


This study was financially supported by Mondelez International (Earlier Cadbury India Ltd) and Kerala Agricultural University, Vellanikkara.

Author information




Juby Baby carried out the work as a part of her post graduate programme under the guidance of Dr. Janaki Seifudeen Minimol as major advisor and the programme was undertaken under a project funded by Mondelez Internationlal. Dr. Suma Basura, Dr. Jiji Joseph and Dr. Santhoshkumar Adiyodi Venugopal as advisory committee members whose whole hearted co-operation helped to complete the project sucessfully. Dr. Panchami Pottekkat Sidharthan is the research associate associated with the project and assisted Juby Baby to complete the project. All authors read and approved the final manuscript. This manuscript is not submitted elsewhere for publication.

Corresponding author

Correspondence to Janaki Seifudeen Minimol.

Ethics declarations

Ethics approval and consent to participate

The article does not contain any studies with human participants or animals performed by any of the authors. Plant materials used in the experiments were from Cocoa Research Centre, Vellanikkara, Thrissur. All methods in the experiment were performed in accordance with the relevant guidelines/regulations/legislation of Kerala Agricultural University, Vellanikkara, Thrissur.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Juby, B., Minimol, J.S., Suma, B. et al. Drought mitigation in cocoa (Theobroma cacao L.) through developing tolerant hybrids. BMC Plant Biol 21, 594 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Nitrate reductase
  • Proline
  • Superoxide dismutase
  • Logistic regression