Skip to main content

Deciphering the molecular basis for photosynthetic parameters in Bambara groundnut (Vigna subterranea L. Verdc) under drought stress

Abstract

Background

Assessment of segregating populations for their ability to withstand drought stress conditions is one of the best approaches to develop breeding lines and drought tolerant varieties. Bambara groundnut (Vigna subterranea L. Verdc.) is a leguminous crop, capable of growing in low-input agricultural systems in semi-arid areas. An F4 bi-parental segregating population obtained from S19-3 × DodR was developed to evaluate the effect of drought stress on photosynthetic parameters and identify QTLs associated with these traits under drought-stressed and well-watered conditions in a rainout shelter.

Results

Stomatal conductance (gs), photosynthesis rate (A), transpiration rate (E) and intracellular CO2 (Ci) were significantly reduced (p < 0.05) while water use efficiency (WUE) was significantly increased (p < 0.05) under drought-stressed conditions. A strong linear correlation was observed between gs, WUE, A, E and Ci under both water regimes. The variability between different water treatment, among individual lines and the interaction between lines and environment for photosynthetic parameters provides resources for superior lines selection and drought resistant variety improvement. Significant QTL for gs and FV/FM under well-watered conditions were mapped on LG5 and LG3, respectively, with more than 20% of the PVE, which could be considered as the major QTL to control these traits. Five clustered QTLs for photosynthetic traits under drought-stressed and well-watered conditions were mapped on LG5, LG6A, LG10 and LG11, respectively.

Conclusions

Significant and putative QTLs associated with photosynthetic parameters and the effect of drought stress on these traits have been revealed by QTL linkage mapping and field experiment in the F4 segregating population derived from S19-3 × DodR in bambara groundnut. The study provides fundamental knowledge of how photosynthetic traits response to drought stress and how genetic features control these traits under drought-stressed and well-watered conditions in bambara groundnut.

Peer Review reports

Background

Drought is one of the major abiotic stresses, negatively impacting plant growth and reduce crop production worldwide [1, 2]. Drought stress caused significant changes in photosynthesis, relative water content, root and shoot dry weight, which are good indicators of drought monitoring in chickpea (Cicer arietinum L.) [3,4,5]. Three common bean (Phaseolus vulgaris L.) elite lines (NCB 226, SER 78, SER 125) showed superior levels of adaptation to drought stress conditions by remobilizing photosynthate to increase grain yield [6].

Bambara groundnut is an underutilised and drought-resistant leguminous crop with high protein content (16%-25%), which are mainly grown by subsistence farmers and served as an edible protein source in Africa [7,8,9,10,11]. It was shown that S19-3 landrace from Namibia experienced reduced respiration and stomata closure at a comparatively lower water threshold coupled with fast phenological development, short life cycle and early maturing proved to be among the mechanisms to ameliorate drought conditions [12, 13]. Three landraces of bambara groundnut collected from South Africa i.e., Brown, Red and Light Brown were reported to have reduced stomatal conductance of 1% – 8%, reduced chlorophyll content index (CCI), plant height, leaf number, reduced leaf area index and biomass accumulation of 5% – 8% and yield loss of 50% under water defict conditions [14]. Landrace Brown and Red showed higher emergence rate, gs, CCI and yielded more than Light Brown in response to water deficit conditions [14].

Similar to most of the underutilised and neglected crop species which have limited established breeding programmes due to lacking of commercial interest in breeding this crop and genetic improvement activities, landraces (mixture of genotypes) have remained as the main source of planting in bambara groundnut [12, 15,16,17,18]. Single plant descent (SPD) and single seed descent (SSD) has been highlighted to develop pure lines/genotypes of bambara groundnut [12]. Variation among genotypes with different drought response ability provides resources for breeders to select drought resistance varieties with high yield in bambara groundnut [13, 14, 19]. Strong genotypic variation was observed for many traits, i.e., 100-seed weight, harvest index, stomatal density and leaf area in the F5 segregating population derived from Tiga Nicuru × DipC, facilitating the identification of superior and drought tolerant lines for advancement [20]. Kendabie et al. [21] reported that five segregating population, i.e., Ankpa4 × IITA-686 (reciprocal), Ankpa4 × DodR, Ankpa4 × DipC, S19-3 × Ankpa4 and IITA-686 × LunT were developed to create genetic linkage map and enhance trait dissection in bambara groundnut to accelerate crop breeding process. A genetic linkage map covering 1,040.92 cM across 11 linkage groups was constructed using 234 DArTseq-based SNP markers in the F2 segregating population from S19-3 × DodR [22]. Significant QTLs associated with number of seeds per plant, number of double-seeded pod per plant, seed weight per plant and pod weight per plant were mapped on LG4 with overlapping confidence intervals under well-watered conditions in the F4 population, which could be considered as major QTL involved in the control of these traits [22]. QTLs associated with stomatal density, length and conductance were co-located on chromosome II under greenhouse conditions in faba bean (Vicia faba L.) [23]. Lopez et al. [24] reported different E values but similar A values at the same QTL in Sorghum (Sorghum bicolor L. Moench). However, few studies have been reported for genetic analysis and variety development in the structured populations of bambara groundnut.

The first genome sequence of bambara groundnut has been assembled with 513 Mb in size and predicted 31,707 protein-coding genes [25]. However, the current bambara groundnut genome information and QTL mapping does not afford adequate resolution to identify genes. High density genetic linkage maps and QTL detection are very useful tools to identify genomic regions that may be responsible for target traits for MAS breeding of bambara groundnut [25]. In the present study, we evaluated the effect of drought stress on photosynthetic parameters and mapped QTLs for these traits under drought-stressed and well-watered conditions in the F4 segregating population derived from S19-3 × DodR in bambara groundnut. The study provided critical insights into how genetic features control photosynthetic traits in bambara groundnut under drought-stressed and well-watered conditions, which is also essential for crop improvement of bambara groundnut in response to drought stress.

Results

Photosynthesis response to drought stress during plant growth

The average of total reduction of soil moisture content under drought-stressed conditionswas 36.15% from 47 to 74 DAS. On average, soil moisture content declined by 0.41% per day at depth 300 mm and 0.31% per day at depth 400 mm over 28 days of drought (Supplementary Fig. S1). Significant reduction (p < 0.01) in soil moisture by 5.66%, 7.04% and 11.51% was observed under drought-stressed conditions compared to well-watered conditions at depth 100 mm, 200 mm and 300 mm, respectively. However, there was no significant difference (p > 0.01) for soil moisture content at depth 400 mm, 600 mm and 1000 mm between drought-stressed and well-watered conditions (Supplementary Fig. S2).

Parental lines showed significant differences (p < 0.05) for gs and Ci at 64 DAS, E, Ci and WUE at 71 DAS, Ci and WUE at 78 DAS between drought-stressed and well-watered conditions (Fig. 1). S19-3 had significantly higher (p < 0.05) E at 57 DAS, A and E at 78 DAS, gs, A, E and Ci at 86 DAS but significantly lower (p < 0.05) WUE6 compared to DodR under drought-stressed conditions (Fig. 1). The interaction between genotype and environment was significant (p < 0.05) for E at 57 DAS, A and E at 78 DAS, gs, A, E, Ci and WUE at 86 DAS among parental lines (Supplementary Table S1).

Fig. 1
figure 1

Comparison of (a) photosynthesis rate, A (b) stomatal conductance, gs (c) transpiration rate, E (d) intracellular CO2, Ci and (E) water use efficiency, WUE between parental lines, S19-3 and DodR under drought-stressed (DS) and well-watered (WW) conditions. Mean and standard error are indicated at the time of measurement. n = 9. Arrow Irrigation was resumed at 74 DAS

On average, A declined from 36.24 μmol m−2 s−1 to 18.61 μmol m−2s−1 (by 48.6%) was observed under drought-stressed conditions from 47 to 71 DAS followed by recovery to 23.3 μmol m−2 s−1 (by 25.20%) at 78 DAS after irrigation was resumed, with significant difference (p < 0.01) observed between drought-stressed and well-watered treatments at 71 DAS (p < 0.01) (Fig. 2a). A significant difference (p < 0.05) was observed for A among the individual lines during drought period at 64 DAS and 71 DAS and at 78 DAS and 86 DAS after irrigation was resumed. The interaction between the individual lines and water treatment was significant (p < 0.05) during drought period at 61 DAS and 71 DAS.

Fig. 2
figure 2

Comparison of (a) photosynthesis rate, A (b) stomatal conductance, gs (c) transpiration rate, E (d) intracellular CO2, Ci and (e) water use efficiency, WUE between individual lines under drought-stressed (DS) and well-watered (WW) conditions. Mean and standard error are indicated at the time of measurement. n = 36. * = Significant at p = 0.05, ** = Significant at p = 0.01. Arrow Irrigation was resumed at 74 DAS

On average, gs declined significantly (p < 0.01) from 0.497 mol m−2 s−1 to 0.203 mol m−2 s−1 (by 59.2%) while WUE increased significantly (p < 0.01) from 79.97 μmol mol−1 to 142.7 μmol mol−1 (by 55.9%) under drought-stressed conditions from 47 to 71 DAS. On average, gs then was observed to recover to 0.276 mol m−2 s−1 (by 35.96%) and WUE to 116.5 μmol mol−1 (by 18.36%) at 78 DAS after irrigation was resumed (Fig. 2b and e). A significant difference was observed for gs among the individual lines before drought was imposed at 46 DAS (p < 0.01), during drought period from 57 to 71 DAS (p < 0.01), and at 78 DAS and 86 DAS after irrigation was resumed (p < 0.05). WUE exhibited significant difference (p < 0.05) among individual lines similar to gs, but they were not significantly different during drought period at 64 DAS (p = 0.079) and after irrigation was resumed at 86 DAS (p = 0.141). Stomatal conductance, gs exhibited significant interaction (p < 0.05) between individual lines and water conditions before drought conditions was imposed at 46 DAS, during drought period at 57 DAS, 64 DAS and at 78 DAS and 86 DAS after irrigation was resumed. WUE exhibited significant interaction (p < 0.05) between individual lines and water conditions similar to gs, but not significantly different during drought period at 64 DAS (p = 0.051).

Similar to gs, on average, Ci significantly (p < 0.01) declined by 50.37% from 273.6 μmol m−1 to 135.8 μmol m−1, after drought stress was imposed at 47 DAS and recovered to 186.6 μmol m−1 (by 37.41%) at 78 DAS after irrigation was resumed, with significant difference (p < 0.01) observed between drought-stressed and well-watered treatments during drought period from 57 to 71 DAS and at 78 DAS and 86 DAS after irrigation was resumed (Fig. 2d). A significant difference was observed for Ci among the individual lines during drought period at 71 DAS (p < 0.05) and at 78 DAS and 86 DAS after irrigation was resumed (p < 0.01). The interaction between individual lines and water conditions for Ci was significant (p < 0.01) before drought stress was imposed, at 46 DAS and during drought period, at 57 DAS.

On average, E significantly (p < 0.01) declined from 12.47 mol m−2 s−1 to 3.54 mol m−2 s−1, a reduction of 71.61%, under drought-stressed conditions from 46 to 71 DAS and recovered to 4.98 mol m−2 s−1 (by 40.68%) at 78 DAS after irrigation was resumed, with significant difference observed between drought-stressed and well-watered treatments during drought period from 57 to 71 DAS (p < 0.01) and at 78 DAS after irrigation was resumed (p < 0.05) (Fig. 2c). A significant difference (p < 0.05) for E was observed among individual lines and interaction between individual lines and water conditions at 71 DAS.

Parental lines showed significant differences (p < 0.05) for FV/FM at 46 DAS and 64 DAS, CCI at 57 DAS and RWC at 78 DAS between drought-stressed and well-watered conditions (Fig. 3). In the F4 segregating population, RWC was reduced by 7.55% under drought-stressed conditions from 81.98% at 57 DAS to 75.79% at 71 DAS, with significant difference observed between drought-stressed and well-watered conditions at 71 DAS (p < 0.05) and after irrigation was resumed at 86 DAS (p < 0.05) (Fig. 3b). A significant difference was observed among the individual lines during drought period at 57 DAS (p < 0.01) and after irrigation was resumed at 78 DAS (p < 0.05). The interaction between individual lines and water conditions was significant (p < 0.05) after drought stress was imposed at 57 DAS.

Fig. 3
figure 3

The effect of drought stress on (a) and (b) relative water content, RWC (c) and (d) chlorophyll content index, CCI (e) and (f) quantum yield of PSII photochemistry, FV/FM in parental lines and the F4 segregating population. Data represent mean values ± standard error. DS, drought-stressed conditions; WW, well-watered conditions. * = Significant at p = 0.05, ** = Significant at p = 0.01. Arrow Irrigation was resumed at 74 DAS

CCI declined by 17.9% over 28 days (p = 0.193) after drought stressed was imposed from 47 to 74 DAS (Fig. 3). CCI under drought-stressed conditions showed 2.5% significant reduction (p < 0.05) at 71 DAS compared to well-watered conditions (Fig. 3d). A significant difference was observed among the individual lines during drought period at 71 DAS (p < 0.01), after irrigation was resumed at 78 DAS (p < 0.05) and 86 DAS (p < 0.05). The interaction between individual lines and water conditions was significant (p < 0.01) at 71 DAS. FV/FM declined gradually (p = 0.208) by 8.48% during drought period from 0.66 at 57 DAS to 0.60 at 71 DAS (Fig. 3f). FV/FM had 7.69% significant reduction (p < 0.01) under drought-stressed conditions at 71 DAS compared to well-watered conditions (Fig. 3f). An increase in FV/FM value up to 5.62% was also observed after irrigation was resumed starting from 74 DAS, with significantly higher (p < 0.05) FV/FM observed at 78 DAS and 86 DAS in well-watered treatment. A significant difference was observed among the individual lines before drought conditions was imposed at 46 DAS (p < 0.01), during drought period at 64 DAS (p < 0.01) and after irrigation was resumed at 78 DAS (p < 0.05) and 86 DAS (p < 0.05). The interaction between individual lines and water conditions was significant (p < 0.05) after drought stress was imposed at 64 and 78 DAS.

The effect of drought stress on photosynthetic parameters

The average of A, E, gs, Ci and WUE, RWC, CCI and FV/FM under drought-stressed and well-watered conditions during drought period from 47 to 74 DAS were presented in Table 1. The reduction of 1.76% and 9.03% in A, 3.26% and 15.23% in Ci was observed under drought-stressed conditions in S19-3 and DodR, respectively, compared to well-watered conditions (Table 1).

Table 1 Effects of drought stress on photosynthetic parameters under drought stressed (DS) and well-watered (WW) conditions in the F4 segregating population derived from S19-3 × DodR and their parental lines

A significant reduction (p < 0.05) of 13.65% in A, 20.59% in E, 47.92% in gs, 26.36% in Ci and 1.58% in RWC was observed in the individual lines under drought-stressed conditions compared to well-watered conditions (Table 1). Compared to well-watered conditions, WUE significantly (p < 0.01) increased by 45.61% under drought-stressed conditions (Table 1). A, gs, WUE and CCI showed significant difference (p < 0.05) among individual lines. The interaction between conditions and individual lines was significant (p < 0.05) for A, WUE and CCI.

A positively correlated with gs (rWW = 0.55, p < 0.01; rDS = 0.53, p < 0.01), Ci (rWW = 0.66, p < 0.01; rDS = 0.46, p < 0.05) and E (rWW = 0.60, p < 0.05; rDS = 0.35, p = 0.07), and negatively correlated with RWC (rWW = –0.45, p < 0.05; rDS = –0.07, p = 0.73) (Supplementary Table S2). gs positively correlated with Ci (rWW = 0.49, p < 0.05; rDS = 0.32, p = 0.09) and E (rWW = 0.41, p < 0.05; rDS = 0.17, p = 0.53), and negatively correlated with WUE (rWW = –0.71, p < 0.01; rDS = –0.65, p < 0.01) (Supplementary Table S2).

Detection of QTLs associated with photosynthetic parameters under drought-stressed and well-watered conditions

Significant and putative QTLs for photosynthetic traits were detected under both water regimes in the F4 segregating population (Fig. 4). Most QTLs were distributed in LG5, LG6A and LG11. Significant QTL for gs5 (LOD: 3.09, 37.7% of the PVE) and significant QTL for FV/FM4 (LOD: 3.06, 36.5% of the PVE) under well-watered conditions were mapped on LG5 and LG3, respectively (Table 2). Putative QTL for gs3 (LOD: 2.12, 26.6% of the PVE) and putative QTL for WUE3 (LOD: 2.26, 28.7% of the PVE) under well-watered conditions were co-located on LG11 (35.83 cM, nearest marker: 4,177,456) with overlapping confidence intervals (Table 2). Putative QTL for Ci6 (LOD: 2.41, 29.8% of the PVE) and putative QTL for WUE5 (LOD: 2.39, 33.5% of the PVE) under well-watered conditions were co-located on LG11 (45.18 cM, nearest marker: 4,183,896) with overlapping confidence intervals (Table 2). Five clustered QTLs were found to have overlapping confidence intervals for photosynthetic traits, which included gs5 and A6 under well-watered conditions, A2 under drought-stressed conditions and FV/FM2 under well-watered conditions on LG5, A6 under drought-stressed conditions and RWC3 under well-watered conditions on LG6A, A5 and WUE2 under drought-stressed conditions on LG10, gs3, WUE3, Ci6 and WUE5 under well-watered conditions on LG11 (Fig. 4).

Fig. 4
figure 4

Map position of the quantitative trait loci (QTL) under drought-stressed (DS) and well-watered (WW) conditions in the F4 segregating population developed from S19-3 × DodR. Rectangular bars represent the 1- and 2-LOD QTL interval (inner and outer interval). Solid rectangular bars represent significant QTLs, while blank bars represent putative QTLs. LG1, LG6 and LG7 were divided into subgroups ‘1A’ and ‘1B’, respectively, based on the association observed in the maximum likelihood mapping (MLM) due to insufficient linkage to complete the map using regression mapping (RM). A Photosynthesis rate, gs Stomatal conductance, E Transpiration rate, Ci Intracellular CO2, WUE Water use efficiency, RWC Relative water content, FV/FM Quantum yield of PSII photochemistry

Table 2 Significant and putative QTLs for photosynthetic traits under drought-stressed (DS) and well-watered (WW) conditions in the F4 segregating population derived from S19-3 × DodR

Discussion

Stomatal closure usually happened during the initial stages of drought stress, which results in the reduction of transpiration in plant leaves, a decrease in CO2 flow into leaves, a decline in net photosynthesis, and ultimately reduced plant growth [14, 26, 27]. In the present study, WUE, calculated as A/gs, increased after drought stress was imposed, then declined gradually after irrigation was resumed under drought-stressed conditions in the F4 segregating population. Singh and Reddy [28] reported that WUE increased under drought stress in 15 cowpea genotypes, suggesting that stomatal regulation was a major limitation to photosynthesis and plant growth. WUE is regulated by gs and multiple factors including the available energy impinging on the leaf, vapour pressure deficit, and aerodynamic exchange [29]. The negative correlation between WUE and gs under well-watered conditions (r = –0.79, p < 0.05) and under drought-stressed conditions (r = –0.63, p < 0.01) suggests that gs decreases faster than A, leading to increased WUE under drought stress [30]. The QTL associated with gs3 and WUE3 were mapped on LG11 with overlapping confidence intervals in the F4 population, which may suggest that these traits are controlled by the same loci. Similar findings have been reported in Sorghum that the QTL for gs was associated with reduced E and increased WUE in Sorghum [24].

Genotypes with high stomatal conductance and WUE in response to drought stress were suggested to have good drought tolerance and adaptation ability [28]. For example, a drought-tolerant cowpea cultivar (PO) maintained higher photochemical activity and leaf gas exchange under water deficit and showed faster recovery of photosynthesis after irrigation was resumed than the drought-sensitive cultivar (SI), revealing possible mechanisms enabling plants to overcome stressful conditions [31, 32]. Plants maintain high water status by reducing stomatal conductance during periods of drought stress, which involves either drought avoidance or tolerance or both mechanisms [30, 33, 34]. The positive correlation between A and Ci under well-watered conditions (r = 0.67, p < 0.05) and under drought-stressed conditions (r = 0.42, p < 0.05) suggests that lower internal CO2 accumulation concentration during drought is responsible for the reduction in photosynthesis [35]. In the present study, the QTL for Ci6 and WUE5 were co-located on LG11 with overlapping confidence intervals, which may suggest that these traits are controlled by the same loci.

Chlorophyll content and FV/FM are non-stomatal limiting factors and capture light energy for plant photosynthesis [30, 36]. In the present study, drought stress significantly reduced (p < 0.05) CCI and FV/FM in the F4 segregating population, which suggests that the ability of bambara groundnut plants to capture light energy for plant photosynthesis is significantly curtailed by drought conditions. Similar to these findings, Mafakheri et al. [26] reported that drought significantly reduced total chlorophyll content (p < 0.05) under drought stress during vegetative growth in three chickpea cultivars. Rahbarian et al. [37] also reported that drought stress reduced FV/FM in two drought-tolerant genotypes and two drought-sensitive genotypes of Chickpea. Additionally, in a study involving several bambara groundnut landraces, CCI was lower under water-deficit compared to irrigated conditions [14, 38]. The FV/FM value was also reported to have declined by 25% at the end of drought stress trial involving three bambara groundnut landraces [39]. Significant QTL associated with FV/FM4 under well-watered conditions was mapped on LG3 with 36.5% of the PVE, while putative QTL for FV/FM4 was mapped on LG7A with reduced PVE (27.2 of the PVE) under drought-stressed conditions in the F4 population. Similar to gs5, a reduced PVE was detected under drought-stressed conditions compared to well-watered conditions, suggesting the traits identified under well-watered conditions were unable to fully express their potential trait values under drought conditions [22].

RWC is an indicator of plant water status revealing the stress intensity [40]. In the present study, RWC increased in the initial stage of drought, declined gradually until the end of drought period, suggesting some individual lines have the ability to adapt to drought stress. RWC was reported to have decreased by 21% – 24% with time after water deficit and increased by 13% – 17% after irrigation was resumed in bambara groundnut [39]. RWC was higher in well-watered plants than drought-stressed plants, although bambara groundnut accessions were still able to maintain high RWC despite the water stress [41]. Keyvan [35] reported wheat cultivars with high RWC under drought stress conditions to be resistant. The putative QTL for RWC3 and seeds weight per plant under well-watered conditions were mapped on LG6A with overlapping confidence intervals, suggesting RWC and seed yield may be controlled by the same gene [22]. Further validation of consensus markers, significant QTLs associated with various traits and candidate genes is required in different populations, across locations and seasons in bambara groudnnut. Individual lines with overall superior performance such as high gs, E, RWC, FV/FM and CCI than S19-3 under drought-stressed conditions are recommended for further field investigation to develop drought-tolerant varieties (Supplementary Table S3). The major QTLs identified in this study are essential to support the development of improved varieties of bambara groundnut in molecular-enabled breeding programmes.

Conclusions

The development of drought resistant materials is essential to cope with the effects of climate change, especially in the tropical arid and semi-arid areas where rainfall is scarce and erratic. Drought stress significantly reduced (p < 0.05) gs, A, E, Ci and RWC, while WUE significantly increased (p < 0.01) under drought-stressed conditions in the F4 segregating population. The linear correlation between photosynthetic parameters suggests the synergy of photosynthesis mechanisms when plant response to drought stress. Significant QTL for gs and FV/FM under well-watered conditions were mapped on LG5 and LG3, respectively. QTLs identified under well-watered conditions would reflect the intrinsic genetic mechanisms underlying photosynthetic parameters. Five clustered QTLs were found to have overlapping confidence intervals for photosynthetic traits, which included gs, A, Ci, WUE, FV/FM and RWC under well-watered conditions and A and WUE under drought-stressed conditions, suggesting these traits are controlled by the same major QTLs. The QTLs identified in this study are essential to identify candidate genes related to photosynthetic traits in response to drought stress in bambara groundnut.

Methods

Plant material and experimental design

A total of 36 individual lines of the F4 segregating population derived from a cross between S19-3 and DodR were evaluated in a rainout shelter at the University of Nottingham Malaysia (2°56′46.74"N; 101°52′24.35"E) with mean air temperature of 36 °C/25 °C day/night and relative humidity of 58%/91% day/night from April to July 2019. The experiment was carried out in a completely randomized design (CRD) with three replicates and two treatments, drought-stressed and well-watered treatments [22]. Each of the replicates was represented by one plant from each of the individual lines. Irrigation for the well-watered conditions was continued throughout the experiment while the drought-stressed conditions was imposed after 100% flowering was observed at 47 days after sowing (DAS) and no further irrigation was applied until early pod-filling stage at 74 DAS, at which irrigation of plants for the drought-stressed conditions was resumed [22].

Field management

A trickle irrigation system was set to irrigate the plants at 07:00 and 19:00 h for 10 min with a flow rate of 2 L/h, with each tube 6 m in length [22]. A distance of 40 cm × 30 cm was kept between the plants. NPK (nitrogen, phosphorus and potassium) fertiliser was applied at a rate of 20:40:60 kg/ha (133 kg/ha NPK (15:15:15), 44 kg/ha TSP (triple-super-phosphate) and 67 kg/ha MOP (muriate of potash) at sowing and after emergence [22]. All other agronomic procedures, such as weeding and spraying of pesticides, were carried out when necessary [22].

Soil moisture content

Two evenly spaced PR2 profile tubes (Delta-T Devices Ltd., Cambridge, UK) were inserted into the centre of each of each plot with a distance of 3 m between two profile tubes in each plot [22]. There were 12 access tubes in total [22]. Three PR2 readings %Vol (volumetric water content as a percentage) were taken twice a week between 0900 and 1100 h at soil depth of 100, 200, 300, 400, 600 and 1000 mm from seeds sowing until maturity [22].

Photosynthetic parameters

A, gs, Ci and E were measured by LI-6400XT Portable Photosynthesis System (Li-Cor, Lincoln, USA). Water use efficiency (WUE) was estimated as the ratio of A/gs [28].

Relative water content (RWC) was calculated as:

$$\mathrm{RWC }= [(\mathrm{Fw}-\mathrm{Dw}) / (\mathrm{Tw}-\mathrm{Dw})] \times 100$$

where FW = fresh weight of leaves, TW = turgid weight of leaves after incubating leaves in distilled water for 24 h, and DW = dry weight of leaves after oven drying at 80 °C for 48 h.

CCI was measured by chlorophyll meter SPAD-502 (Spectrum Technologies, Inc., Aurora, Illinois, USA). Three readings were taken per leaf, three leaves per plant, and averaged to give a final reading. Quantum yield of PSII photochemistry (FV/FM) was estimated from dark-adapted leaves for 30 min using FlourPen FP 100 (PSI, CZ, Czech Republic). Photosystem II quantum yield is equivalent to ratio of variable fluorescence/maximal fluorescence (FV/FM) in dark-adapted samples.

Readings were taken for A, gs, E, Ci, RWC, CCI and FV/FM on the middle leaflet of one most fully expanded leaf between 08:00 and 12:00 h for all individual lines, starting from 50% flowering observed at 46 DAS before drought conditions was imposed, during drought period at 57 DAS, 64 DAS and 71 DAS, and after irrigation was resumed, at 78 DAS and 86 DAS.

Genetic linkage map construction and QTL analysis

The genetic linkage map was constructed using an F2 individual population derived from the same parents [22]. Genetic linkage map and phenotypic data from drought-stressed and well-watered conditions were subjected to QTL analysis using MapQTL 6.0 software [42]. The significant threshold of the Genome-Wide (GM) LOD threshold was obtained from permutation test using 10,000 repetitions at p < 0.05 (5%). Interval mapping (IM) was carried out following the permutation test and the LOD values from IM was compared with GW LOD threshold at p < 0.05 from the permutation test. Significant QTLs were detected if the LOD score was equivalent or higher than GM LOD threshold. Putative QTLs were detected if the LOD score was lower than GM LOD threshold by up to a 1-LOD interval [20, 22]. MapChart 2.3.2 [43] was used to depict the linkage groups and QTLs. QTLs explaining more than 20.0% of the phenotypic variation (occurring at least once) or more than 10% of the phenotypic variation (occurring at least twice) were defined as major QTL, whereas QTL was defined as minor QTL [44, 45].

Data analysis

Normality of trait data was examined using Shapiro–Wilk normality test and data transformation was performed for non-normally distributed trait data. Two-way analysis of variance (ANOVA) and Pearson’s correlation coefficient analysis were conducted for all photosynthetic parameters using 18th edition of Genstat Statistical package (18th edition, VSN International, UK).

Availability of data and materials

The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.

Abbreviations

A:

Photosynthesis rate

A2:

Photosynthesis rate at 57 days after sowing

A3:

Photosynthesis rate at 64 days after sowing

A4:

Photosynthesis rate at 71 days after sowing

A5:

Photosynthesis rate at 78 days after sowing

A6:

Photosynthesis rate at 86 days after sowing

ANOVA:

Analysis of variance

CCI:

Chlorophyll content index

Ci:

Intracellular CO2

Ci1:

Intracellular CO2 at 46 days after sowing

Ci2:

Intracellular CO2 at 57 days after sowing

Ci6:

Intracellular CO2 at 86 days after sowing

CRD:

Completely randomized design

DArTseq-based SNP markers:

Diversity array technology sequencing based single nucleotide polymorphism markers

DAS:

Days after sowing

DS:

Drought-stressed

E:

Transpiration rate

E5:

Transpiration rate at 78 days after sowing

F V/F M :

Quantum yield of PSII photochemistry

FV/FM1:

Quantum yield of PSII photochemistry at 46 days after sowing

FV/FM2:

Quantum yield of PSII photochemistry at 57 days after sowing

FV/FM3:

Quantum yield of PSII photochemistry at 64 days after sowing

FV/FM4:

Quantum yield of PSII photochemistry at 71 days after sowing

FV/FM6:

Quantum yield of PSII photochemistry at 86 days after sowing

G*E:

Interaction between conditions and genotypes

gs:

Stomatal conductance

gs3:

Stomatal conductance at 64 days after sowing

gs5:

Stomatal conductance at 78 days after sowing

GW LOD:

Genome-Wide logarithm of odds

IM LOD:

Interval mapping logarithm of odds

MLM:

Maximum likelihood mapping

MOP:

Muriate of potash

NPK:

Nitrogen, phosphorus and potassium

PVE:

Phenotypic variation explanation

QTL:

Quantitative trait loci

RM:

Regression mapping

RWC:

Relative water content

RWC2:

Relative water content at 57 days after sowing

RWC3:

Relative water content at 64 days after sowing

RWC4:

Relative water content at 71 days after sowing

SD:

Standard deviation

TSP:

Triple-super-phosphate

WUE:

Water use efficiency

WUE2:

Water use efficiency at 57 days after sowing

WUE3:

Water use efficiency at 64 days after sowing

WUE5:

Water use efficiency at 78 days after sowing

WW:

Well-watered

References

  1. Zhu J. Salt and drought stress signal transduction in plants. Annu Rev Plant Biol. 2011;53:247–73. https://doi.org/10.1146/annurev.arplant.53.091401.143329.

    Article  CAS  Google Scholar 

  2. Farooq M, Gogoi N, Barthakur S, Baroowa B, Bharadwaj N, Alghamdi SS, Siddique KHM. Drought stress in grain legumes during reproduction and grain filling. J Agron Crop Sci. 2017;203(2):81–102. https://doi.org/10.1111/jac.12169.

    Article  Google Scholar 

  3. Farooq M, Wahid A, Kobayashi N, Fujita D, Basra SMA. Plant drought stress: effects, mechanisms and management. Agron Sustain Dev. 2009;29:185–212. https://doi.org/10.1051/agro:2008021.

    Article  Google Scholar 

  4. Pang J, Turner NC, Khan T, Du YL, Xiong JL, Colmer TD, Devilla R, Stefanova K, Siddique KHM. Response of Chickpea (Cicer Arietinum L.) to terminal drought: leaf stomatal conductance, pod abscisic acid concentration, and seed set. J Exp Bot. 2017;68(8):1973–85. https://doi.org/10.1093/jxb/erw153.

    Article  CAS  PubMed  Google Scholar 

  5. Khan H, Gul R, Khan NU, Naz R, Shah S, Asim N, Latif A. Role of selection indices in ascertaining high yielding drought stress tolerant chickpea (Cicerarietinum L.). J Anim Plant Sci. 2018;28(1):146–54.

    Google Scholar 

  6. Rao IM, Beebe SE, Polania J, Grajales M, Cajiao C, Ricaurte J, García R, Rivera M. Evidence for genotypic differences among elite lines of Common bean in the ability to remobilize photosynthate to increase yield under drought. J Agric Sci. 2017;155(6):857–75. https://doi.org/10.1017/S0021859616000915.

    Article  CAS  Google Scholar 

  7. Basu S, Roberts JA, Azam-Ali SN, Mayes S. Development of microsatellite markers for Bambara groundnut (Vigna Subterranea L. Verdc.) - an underutilized African legume crop species. Mol Ecol Notes. 2007;7(6):1326–8. https://doi.org/10.1111/j.1471-8286.2007.01870.x.

    Article  CAS  Google Scholar 

  8. Aliyu S, Massawe F, Mayes S. Beyond landraces: developing improved germplasm resources for underutilized species - a case for Bambara groundnu. Biotechnol Genet Eng Rev. 2015;30(2):127–41. https://doi.org/10.1080/02648725.2014.992625.

    Article  Google Scholar 

  9. Massawe F, Mayes S, Cheng A. Crop diversity: an unexploited treasure trove for food security. Trends Plant Sci. 2016;21(5):365–8. https://doi.org/10.1016/j.tplants.2016.02.006.

    Article  CAS  PubMed  Google Scholar 

  10. Atoyebi JO, Oyatomi O, Osilesi O, Adebawo O, Abberton M. Morphological characterisation of selected African accessions of Bambara groundnut (Vigna subterranea (L.) Verdc.). Int J Plant Res. 2017;7(2):29–35. https://doi.org/10.5923/j.plant.20170702.01.

    Article  Google Scholar 

  11. Halimi RA, Barkla BJ, Mayes S, King GJ. Characteristics of the underutilised pulse Bambara groundnut (Vigna Subterranea (L) Verdc.) relevant to food & nutritional security. Proceedings. 2020;36(1):199. https://doi.org/10.3390/proceedings2019036199.

    Article  Google Scholar 

  12. Massawe FJ, Mwale SS, Roberts JA. Breeding in Bambara groundnut (Vigna subterranea (L.) Verdc.): Strategic Considerations. Afr J Biotechnol. 2005;4(6):463–71. https://doi.org/10.5897/AJB2005.000-3084.

    Article  Google Scholar 

  13. Jørgensen ST, Liu F, Ouedraogo M, Ntundu WH, Sarrazin J, Christiansen JL. Drought responses of two Bambara groundnut (Vigna ubterranea L. Verdc.) landraces collected from a dry and a humid area of Africa. J Agron Crop Sci. 2010;196. https://doi.org/10.1111/j.1439-037X.2010.00435.x

  14. Mabhaudhi T, Modi AT. Growth, phenological and yield responses of a Bambara groundnut (Vigna subterranea (L.) Verdc.) landrace to imposed water stress under field conditions. S Afr J Plant Soil. 2013;30(2):69–79. https://doi.org/10.1080/02571862.2013.790492.

    Article  Google Scholar 

  15. Olukolu BA, Mayes S, Stadler F, Ng NQ, Fawole I, Dominique D, Azam-Ali SN, Abbott AG, Kole C. Genetic diversity in Bambara groundnut (Vigna subterranea (L.) Verdc.) as revealed by phenotypic descriptors and DArT marker analysis. Genet Resour Crop Evol. 2012;59(3):347–58. https://doi.org/10.1007/s10722-011-9686-5.

    Article  Google Scholar 

  16. Mayes S, Kendabie P, Ho WK, Massawe F, Azam-Ali SN. Increasing the contribution that underutilised crops could make to food security – Bambara groundnut as an example. Asp Appl Biol. 2015;124:1–8.

    Google Scholar 

  17. Mayes S, Ho WK, Chai HH, et al. Bambara groundnut: An exemplar underutilised legume for resilience under climate change. Planta. 2019;250(3):803–20. https://doi.org/10.1007/s00425-019-03191-6.

    Article  CAS  PubMed  Google Scholar 

  18. Muhammad I, Rafii MY, Ramlee SI. Exploration of Bambara groundnut (Vigna Subterranea (L.) Verdc, an underutilized crop, to aid global food security: varietal improvement, genetic diversity and processing. Agronomy. 2020;10(6):1–20. https://doi.org/10.3390/agronomy10060766.

    Article  CAS  Google Scholar 

  19. Nautiyal PC, Kulkarni G, Singh AL, Basu MS. Evaluation of water-deficit stress tolerance in Bambara groundnut land races for cultivation in sub-tropical environments in india evaluation of water-deficit stress tolerance in Bambara groundnut landraces for cultivation in sub-tropical environments. Indian J Plant Physiol. 2017;22:190–6. https://doi.org/10.1007/s40502-017-0296-x.

    Article  Google Scholar 

  20. Chai HH, Massawe F, Mayes S. Effects of mild drought stress on the morpho-physiological characteristics of a Bambara groundnut segregating population. Euphytica. 2016;208(2):225–36. https://doi.org/10.1007/s10681-015-1581-2.

    Article  CAS  Google Scholar 

  21. Kendabie P, Massawe F, Mayes S. Developing genetic mapping resources from landrace-derived genotypes that differ for photoperiod sensitivity in Bambara groundnut (Vigna Subterranea L.). Asp Appl Biol. 2015;124:1–8.

    Google Scholar 

  22. Gao X, Chai HH, Ho WK, Kundy AC, Mateva KI, Mayes S, Massawe F. Genetic linkage map construction and identification of QTLs associated with agronomic traits in Bambara groundnut (Vigna Subterranea (L.) Verdc.) using DArTseq-Based SNP markers. Food Energy Secur. 2022;00:e353. https://doi.org/10.1002/fes3.353.

    Article  CAS  Google Scholar 

  23. Khazaei H, O’Sullivan DM, Sillanpää MJ, Stoddard FL,. Use of synteny to identify candidate genes underlying QTL controlling stomatal traits in faba bean (Vicia faba L.). Theor Appl Genet. 2014;127(11):2371–85. https://doi.org/10.1007/s00122-014-2383-y.

    Article  PubMed  Google Scholar 

  24. Lopez JR, Erickson JE, Munoz P, Saballos A, Felderhoff TJ, Vermerris W. QTLs associated with crown root angle, stomatal conductance, and maturity in sorghum. Plant Genome. 2017;10(2). https://doi.org/10.3835/plantgenome2016.04.0038

  25. Chang Y, Liu H, Liu M, et al. The draft genomes of five agriculturally important African orphan crops. GigaScience. 2018;8(3):1–16. https://doi.org/10.1093/gigascience/giy152.

    Article  CAS  Google Scholar 

  26. Mafakheri A, Siosemardeh A, Bahramnejad B, Struik PC, Sohrabi E. Effect of drought Stress on yield, proline and chlorophyll contents in three chickpea cultivars. Aust J Crop Sci. 2010;4(8):580–5.

    CAS  Google Scholar 

  27. Ashraf M, Harris PJC. Photosynthesis under stressful environments : an overview. 2013;51(2):163–90. https://doi.org/10.1007/s11099-013-0021-6.

  28. Singh SK, Reddy KR. Biology regulation of photosynthesis, fluorescence, stomatal conductance and water-use efficiency of Cowpea (Vigna nguiculata [L.] Walp.) under drought. J Photochem Photobiol B Biol. 2011;105(1):40–50. https://doi.org/10.1016/j.jphotobiol.2011.07.001.

    Article  CAS  Google Scholar 

  29. Hatfield JL, Dold C. Water-use efficiency: advances and challenges in a changing climate. Front Plant Sci. 2019;10:1–14. https://doi.org/10.3389/fpls.2019.00103.

    Article  Google Scholar 

  30. Chaves MM, Oliveira MM. Mechanisms underlying plant resilience to water deficits: prospects for water-saving agriculture. Int J Curr Microbiol Appl Sci. 2004;55(407):2365–84. https://doi.org/10.1093/jxb/erh269.

    Article  CAS  Google Scholar 

  31. Bastos EA, Nascimento S, Silva EMD, Filho FFR, Gomide RL. Identification of cowpea genotypes for drought tolerance. Rev Ciência Agronômica. 2011;42:100–7. https://doi.org/10.1590/s1806-66902011000100013.

    Article  Google Scholar 

  32. Rivas R, Falcão HM, Ribeiro RV, Machado EC, Pimentel C, Santos MG. Drought tolerance in cowpea species is driven by less sensitivity of leaf gas exchange to water deficit and rapid recovery of photosynthesis after rehydration. South African J Bot. 2016;103:101–7. https://doi.org/10.1016/j.sajb.2015.08.008.

    Article  CAS  Google Scholar 

  33. Turner NC, Wright GC, Siddique KHM. Adaptation of grain legumes (Pulses) to water-limited environments. In: Advances in Agronomy. Amsterdam: Elsevier; 2001. p. 193–231.

    Google Scholar 

  34. Kavar T, Maras M, Kidrič M, Šuštar-Vozlič J, Meglič V. Identification of genes involved in the response of leaves of phaseolus vulgaris to drought stress. Mol Breed. 2008;21(2):159–72. https://doi.org/10.1007/s11032-007-9116-8.

    Article  CAS  Google Scholar 

  35. Lidon ZZ, Cebola. An overview on drought induced changes in plant growth, water relations and photosynthesis. Emi J Food Agric. 2012;24(1):57–72. https://doi.org/10.9755/ejfa.v24i1.10599.

    Article  Google Scholar 

  36. Keyvan S. The effects of drought stress on yield, relative water content, proline, soluble carbohydrates and chlorophyll of bread wheat cultivars. J Anim Plant Sci. 2010;8(3):1051–60.

    Google Scholar 

  37. Rahbarian R, Khavari-Nejad R, Ganjeali A, Bagheri A, Najafi F. Drought stress effects on photosynthesis, chlorophyll fluorescence and water relations in tolerant and susceptible Chickpea (Cicer Arietinum L.) genotypes. Acta Biol Crac Ser Bot. 2011;53(1):47–56. https://doi.org/10.2478/v10182-011-0007-2.

    Article  Google Scholar 

  38. Muhammad YY, Massawe F. Photosynthetic gas exchange and chlorophyll in Bambara groundnut (Vigna subterranea L. Verdc.) subjected to water deficit. Bayero J Pure Appl Sci. 2015;8(1):50. https://doi.org/10.4314/bajopas.v8i1.10.

    Article  Google Scholar 

  39. Muhammad YY, Mayes S, Massawe F. Effects of short-term water deficit stress on physiological characteristics of Bambara groundnut (Vigna subterranea (L.) Verdc.). S Afr J Plant Soil. 2016;33(1):51–8. https://doi.org/10.1080/02571862.2015.1056847.

    Article  Google Scholar 

  40. Lawlor DW, Cornic G. Photosynthetic carbon assimilation and associated metabolism in relation to water deficits in higher plants. Plant Cell Environ. 2002;25(2):275–94. https://doi.org/10.1046/j.0016-8025.2001.00814.x.

    Article  CAS  PubMed  Google Scholar 

  41. Abejide DR, Falusi OA, Adebola MO, Gana AS, Abubakar A, Dangana M. Effect of water stress on physiological parameters on Bambara groundnut (Vigna subterranea (L.) Verdc.) accessions. J Plant Dev. 2020;27:111–20. https://doi.org/10.33628/jpd.2020.27.1.111.

    Article  Google Scholar 

  42. Van Ooijen JW. MapQTL 6. Software for the Mapping of Quantitative Trait Loci in Experimental Populations of Diploid Species. Kyazma BV, Wageningen, Netherlands. 2009. https://www.kyazma.nl/index.php/MapQTL/.

  43. Voorrips RE. Computer note MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered. 2002;93(1):77–8. https://doi.org/10.1093/jhered/93.1.77.

    Article  CAS  PubMed  Google Scholar 

  44. Maccaferri M, Sanguineti MC, Corneti S, Ortega JLA, Salem MB, Bort J, DeAmbrogio E, Moral LFG, Demontis A, El-Ahmed A, Maalouf F, Machlab H, Martos V, Moragues M, Motawaj J, Nachit M, Nserallah N, Ouabbou H, Royo C, Slama A, Tuberosa R. Quantitative Trait Loci for Grain Yield and Adaptation of Durum Wheat (Triticum durum Desf.) Across a Wide Range of Water Availability. Genetics. 2008;178:489–511. https://doi.org/10.1534/genetics.107.077297.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Price AH. Believe it or not, QTLs are accurate! Trends Plant Sci. 2006;11:213–6.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank Dr. Aloyce Callist Kundy from The University of Nottingham Malaysia for helping with trait measurements.

Funding

This study was supported by School of Chemistry and Chemical Engineering, North University of China (304-98001816), China and Future Food Beacon, School of Biosciences, University of Nottingham Malaysia.

Author information

Authors and Affiliations

Authors

Contributions

XG, HHC, SM and FM conceived and designed research. XG conducted experiments. XG, HHC, WKH analyzed data. XG wrote the manuscript. All authors read and approved the manuscript.

Corresponding authors

Correspondence to Xiuqing Gao or Festo Massawe.

Ethics declarations

Ethics approval and consent to participate

The related plant was carried out with the permission of the University of Nottingham Malaysia, and all the methods were carried out in accordance with relevant guidelines and regulations.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Supplementary Table S1.

The effects of drought stress on photosynthetic parameters under drought-stressed (DS) and well-watered (WW) conditions in parental lines, S19-3 and DodR.

Additional file 2: Supplementary Table S2.

Correlation coefficient analysis of photosynthetic parameters under drought-stressed and well-watered conditions in the F4 segregating populations of bambara groundnut derived from S19-3 × DodR.

Additional file 3:

Supplementary Table S3. Potential with superior performance than S19-3 for advancement based on photosynthesis rate, stomatal conductance, transpiration rate, intracellular CO2, water use efficient, relative water content, chlorophyll content index and quantum yield of PSII photochemistry in the F4 segregating population derived from S19-3 × DodR.

Additional file 4: Supplementary Fig S1.

Soil moisture content measurements at depth 100 mm, 200 mm and 300 mm based on PR2 reading (% vol) under drought-stressed (DS) well-watered (WW) conditions. Data represent mean values of soil moisture content during plant growth season in 2019; n = 6. Data represent mean values ± standard error. (Gao et al. 2022 [22]).

Additional file 5: Supplementary Fig S2.

Soil moisture content measurements at depth of 400 mm, 600 mm and 1000 mm based on PR2 reading (% vol) under drought stress (DS) conditions plots and well-watered (WW) conditions plots. Data represent mean values of soil moisture content during plant growth season in 2019; n = 6. Data represent mean values ± standard error (Gao et al. 2022 [22]).

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 http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) 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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gao, X., Chai, H.H., Ho, W.K. et al. Deciphering the molecular basis for photosynthetic parameters in Bambara groundnut (Vigna subterranea L. Verdc) under drought stress. BMC Plant Biol 23, 287 (2023). https://doi.org/10.1186/s12870-023-04293-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-023-04293-w

Keywords