Abstract

Cowpea (Vigna unguiculata (L.) Walp.) is an important commodity in West Africa. Its seeds are a valuable source of protein, vitamins, and income for humans. However, cowpea cultivation in Benin faces climatic constraints such as water stress caused by a prolonged absence of rain during the rainy season. Thus, this work aims at selecting cowpea varieties that can be cultivated in times of drought without compromising their yields and yield components. Twenty cowpea varieties were used, including 17 improved cultivars and 3 landraces. The experiment was conducted at the International Institute of Tropical Agriculture in Benin and laid at a split-plot design with four replicates. Each genotype was exposed to three water treatments: fully irrigated control, vegetative stress (when plants were 23 days old, drought stress was imposed for 30 days), and reproductive stress (once the first flowers were observed, water stress was imposed for 30 days). The results showed that photochemical yield, chlorophyll content, and relative water content were reduced under water deficit at the vegetative and reproductive stages. But there were no significant differences in proline content among cowpea varieties. Agronomic traits such as number of days to flowering, number of pods, yield per plant, the weight of 100 seeds, and harvest time showed significant differences under water stress. Overall, the landraces and cultivars including Kpodjiguegue, KVX 61-1, and IT 06-K-242-3 were the most tolerant to drought stress at the vegetative and reproductive stages and could potentially be used in breeding programs to improve drought tolerance of cowpeas.

1. Introduction

Cowpea (Vigna unguiculata (L.) Walp.) is one of the oldest crops known to man. It is native to Africa, and it was domesticated from a wild plant to a cultivated plant [1]. It is widely adapted and cultivated around the world on about 14.5 million hectares, and approximately 6.5 million metric tons of cowpea is yearly harvested [2]. Of the total area of around 14.5 million hectares planted with cowpeas worldwide, West Africa alone accounted for around 10.6 million hectares in 2017 [3]. Many cowpea producers in Africa are smallholder farmers. Many Africans consume cowpea as fresh pods, tender leaves, and grains [4, 5]. It is the most economically important indigenous African pulse crop [6] and is of vital importance to the livelihoods of millions of people in West and Central Africa. In Benin, cowpea is cultivated nationwide mainly for its edible seeds [7], and it occupies 7% of the areas sown for annual crops with a production of 11,224 tonnes [8, 9].

One of the key characteristics of cowpea is that it can withstand dry environmental conditions, making it the crop of choice in semiarid/arid areas of West and Central Africa. Additionally, when grown together, cowpea is the first crop harvested before cereals get ready and is therefore referred to as a “hungry season crop” [10]. Owing to its capacity for symbiotic fixation of atmospheric nitrogen, the inclusion of cowpea in crop rotations makes it possible to meet the nitrogen fertilizer needs of subsequent crops which help increase their yields and subsequently contributes to the sustainability of cropping systems [1113].

Due to its high protein (19–25%), carbohydrate, and mineral content, cowpea plays an important role in human nutrition and the fight against malnutrition [14]. With over 25% protein in the seeds as well as in the young leaves (in dry weight), cowpea is a major source of protein, minerals, and vitamins in the daily human diet and is just as important as nutritious fodder for livestock [15, 16]. Cowpea is Africa’s most versatile crop; it feeds people and their livestock.

Water deficit is a key limiting factor to cowpea production due to poor and irregular rainfall [17]. However, cowpea exhibits, to some extent, resistance to water stress when compared to other crops [1719]. Cowpea demonstrates tolerance to water stress conditions that could have been very disastrous to several other annual crops. It is also able to produce a higher yield under dry conditions than any other cultivated crop. When resuming watering after a drought spell, the yield of cowpea is similar to that of a normally irrigated control [20]. However, cowpea is heavily damaged by repeated water deficits in the savannah and Sahel [21]. Early maturing cowpea genotypes are susceptible to water deficit during flowering and fruit setting phases [22]. Cowpeas yield more than 1000 kg of grains/ha, but drought reduces this potential to around 360 kg/ha especially when stress sets in during preflowering [23, 24]. Similarly, in some drier locations in Niger, the average yield is about 20 times lower than the average yield in the United States where water is plentiful [25]. IPCC [26] has shown that predicted climate changes will further exacerbate extreme events. Cowpea resistance to water stress lies in its ability to maintain leaf turgor [27]. Despite its capability to resist water deficit more than any other legume grown in the tropics [28, 29], a significant difference exists between cowpea varieties adapted to water deficit conditions [30, 31]. “Morphophysiological characteristics can be used to select for drought-tolerant genotypes.” Therefore, the objectives of this study were to (1) study the effect of water stress at the vegetative and reproductive stages of twenty cowpea varieties, (2) identify the physiological and agronomic traits that explain drought tolerance in cowpea, and (3) identify high-yielding performing genotypes that can withstand drought conditions.

2. Materials and Methods

2.1. Plant Material and Growth Conditions

The experiment was conducted at the International Institute of Tropical Agriculture (IITA) in Benin, Africa. Twenty cowpea genotypes, including three local landraces and seventeen hybrids, were screened for their tolerance to drought. The landraces collected originally in Benin were Gboto (V7), Kpodjiguegue (V19), and Tawa (V20), and the improved varieties tested are commonly grown in Benin, as well as Burkina-Faso, and included IT 97K-556-6 (variety 1 = V1), KVX 61-1 (V2), KVX 396-18 (V3), IT 87 S-1390 (V4), IT 99K-573-2-1 (V5), IT 86 D-888 (V6), IT 84 S-2246-4 (V8), IT 07K-206-1-1 (V9), IT 99K-573-1-1 (V10), IT 84 D-449 (V11), IT 84 E-124 (V12), IT 93K-452-1 (V13), IT 86 D-1038 (V14), IT 07K-188-49 (V15), IT 06K-242-3 (V16), IT 07 K-211-1-8 (V17), and IT 86 D-1033 (V18). Each experimental pot was filled with 20 kg of soil with 500 g of poultry manure. Five percent diluted carbofuran was used to treat the soil against microorganisms. The seeds were sown on November 18, 2019. The seeds were watered daily. Plants were grown under natural conditions. During the experiment, weather data were collected and plotted (Figure 1). The pots were hand-weeded regularly and were sprayed weekly with abamectin (18 g/L) and cyhalothrin (25 g/L) to protect the plants from the attack of insect pests and plant pathogens. The drought stress at the vegetative stage was imposed on December 10, 2019, and ended on January 08, 2020, while the drought stress at the reproductive stage started on December 23, 2019, and it was terminated on January 21, 2020. There was no rainfall during the drought imposition at both vegetative and reproductive stages.

2.2. Drought Treatments and Experimental Design

A split plot with 4 replications and 2 factors was used. The main factor was water regime with 3 levels: RO = no drought, R1 = vegetative stress, and R2 = reproductive stress. That is, two drought treatments were initiated, one at the vegetative stage and one at the reproductive stage. For the vegetative stage, drought treatments were initiated 23 days after sowing during. Plants were either irrigated daily or exposed to drought (no water) for 30 days. For the reproductive stage, the drought stress was initiated 40 days after sowing, when the first flowers appeared. During this reproductive phase, plants were either watered daily or exposed to drought for 30 days. The subplot was the variety factor with 20 levels (V1 to V20). This generated 60 interactions repeated 4 times given 240 experimental units (Pots). The pots were spaced 0.6 m apart within the rows and 0.4 m apart between the rows. Replications were separated from each other by 2 m.

2.3. Measurements

A portable chlorophyll meter (Minolta SPAD-502, Soil Plant Analysis Development, Minolta Co., Osaka, Japan) was used to determine the chlorophyll content of the cowpea plants. Chlorophyll fluorescence was also measured with a portable chlorophyll fluorometer (Model OS-30P; Opti-Sciences, New Hampshire, USA). Data on fluorescence were collected from fully expanded leaflets clipped from the upper part of the plant canopy, after a 1-hour dark adaptation period. Dark fluorescence (Fo), maximal fluorescence (Fm), and photochemical yield (/Fm, where  = FmFo) were recorded. After stress, the relative water content (RWC) of different plants was determined. Here, 3 leaflets on the third leaf from the apex were taken per plant in each treatment and immediately wrapped in aluminum foil to avoid loss of water. The samples were brought to the laboratory to determine their fresh weight (FW). The leaflets were placed in bottles containing 150 ml of distilled water left in the dark for 24 hours at 16°C [32]. The leaf discs were then removed and weighed again to determine the weight of the full turgor (WFC). The dry weight (DW) was obtained after oven-drying at 80°C for 24 hours. The weighing of the leaf discs was carried out using a 0.01 g precision balance. The relative water content (RWC) of the leaves was calculated by the formula of Turner and Begg [33]: RWC = (FW − DW)/(WFC − DW) × 100. Proline was determined according to Bates and Hall [34]. Under drought stress at the reproductive stage, fresh tissue (0.25 g) was collected and homogenized in 10 ml of water with 3% of aqueous sulfosalicylic acid. After 3 hours, the mixture was centrifuged at 1500 g for 10 min. After removing 2 ml of the supernatant, 2 ml of glacial acetic acid and 2 ml of hydrofluoric acid were added to the supernatant and boiled at 100°C in a water bath for one hour. The reaction was stopped by placing it on ice. 4 ml of toluene was added and mixed vigorously using the vortex for 15–20 seconds, and the toluene containing the chromophore was separated using a separatory funnel and the absorbance was measured at 520 nm in a spectrophotometer against an appropriate toluene blank. The proline content was determined from a standard curve prepared with L-proline and expressed in mg/g MS. The unknown proline content is calculated from the samples using the standard graph. The proline concentration is then calculated withwhere e is the μg/ml proline, V1 is the toluene volume, V2 is the volume of sulfosalicylic acid, and g Fr.MF is the of the dry sample.

Agronomic traits measured included the number of days from sowing to flowering, number of pods per plant, yield per plant, the weight of 100 grains, and harvest time, i.e., from the first harvest to the last.

2.4. Statistical Analysis

Data were analyzed statistically using analysis of variance (ANOVA). Principal component analysis (PCA) was also carried out based on all measured parameters (physiological and agronomic traits measured) to identify the main characteristics contributing to the variability observed between the varieties. Based on these groups, an analysis was carried out to identify water stress yet high-yielding genotypes. All analyses were performed with R version 4.0 software. The separation of the means for the various measured parameters was carried out by the Student Newman Keuls test.

3. Results

3.1. Physiological Traits

Water stress significantly () reduced the /Fm, RWC, and chlorophyll content of genotypes. However, KVX 61–1 (V2), Kpodjiguegue (V19), and GBOTO (V7) recorded the highest values for /Fm (0.90825; 0.88225, and 0.80675, respectively) and RWC (85.3700, 90.8750, and 90.1350, respectively) among stressed plants at the end of vegetative stress (Table 1). At the end of reproductive stress, there was a significant () difference between genotypes for physiological parameters such as photosynthetic yield, chlorophyll content, and relative water content. Chlorophyll fluorescence decreased in stressed plants during the vegetative phase. Nevertheless, KVX 61-1 and Kpodjiguegue were not affected by the stress because their photosynthetic yield values obtained were optimal (0.800). At the end of the reproductive phase, we observed the same result for KVX 61-1 and Kpodjiguegue. But KVX 61-1 (V2), Gboto (V7), and Kpodjiguegue (V19) recorded the highest average mean values of /Fm, thus indicating the genotypes with the best performance under drought stress during the reproductive stage. Chlorophyll content increased in almost all the stressed plants at the vegetative phase; especially in KVX 61-1 (V2), IT 84E-124 (V12), IT 93K-452-1 (V13), IT 86 D-1038 (V14), Kpodjiguegue (V19), and Gboto (V7) whereas, at the reproductive stage, their chlorophyll contents were negatively affected by drought conditions when compared with the controls, indicating that reproductive stage was the most affected. Relative water content was significantly () decreased in stressed plants during the vegetative and reproductive phases. KVX 61-1, Gboto, and Kpodjiguegue recorded the highest relative water content during the vegetative phase, while KVX 61-1 (V2), IT 84 S-2246-4 (V8), IT 84 D-449 (V11), IT 07K-188-49 (V15), IT 06K-242-3 (V16), and Kpodjiguegue (V19) performed better than others in the reproductive stage. The proline content in plants under drought conditions was higher than that of the control plants. There were significant differences () in proline content among the varieties under drought stress during the reproductive stage (Figure 2). The varieties Kpodjiguegue, IT 06K-242-3 (V16), and Gboto (V7) had the highest proline contents under water deficit during the reproductive stage (Table 2).

3.2. Agronomic Traits

There was a significant () difference between the water regimes for all agronomic parameters (Table 3). For the stress during the vegetative phase, the number of pods increased as well as the number of days to flowering while for the stress at the reproductive phase, the weight of seeds significantly decreased. However, IT 97K-556-6, KVX 61-1, IT 07K211-1-8, and Kpodjiguegue gave the best performance among genotypes subjected to drought stress at the vegetative phase, which had the best yields with large seeds. IT 06K242-3, Kpodjiguegue, and KVX 61-1 demonstrated high yield with small seeds among the 20 genotypes under water deficit stress at the reproductive phase. IT 06K-242-3 had the biggest seeds among the genotypes subjected to drought conditions at the vegetative phase, while KVX 61-1 and Kpodjiguegue were more productive than others under drought at the reproductive phase (Table 3).

The varieties KVX 61-1, IT 87 S-1390, IT 06K-242-3, IT 07K-211-1-8, and Kpodjiguegue recorded the highest number of pods under drought stress at the vegetative stage, obtaining 24, 18, 28, 25, and 39 pods, respectively, while KVX 61-1, IT 87 S-1390, IT 99K-573-2-1, IT 84 S-2246-4, IT 07K-188-49, IT 84 S-2246-4, IT 07K-188-49, IT 06K-242-3, and Kpodjiguegue produced 18, 23, 19, 24, and 30 pods under drought stress at reproductive stage, respectively. On the other hand, IT 97K-556-6, IT 07K-206-1-1, and IT 86 D-1033 had fewer pods under water stress at the vegetative stage, producing 13, 6, and 13 pods, respectively, and KVX 396-18, IT 07K-206-1-1, IT 99K-573-1-1, IT 84E-124, and IT 07K-211-1-8 had 5, 6, 7, 8, and 7 pods, respectively, under water stress at reproductive stress (Figures 2 and 3).

3.3. Principal Component Analysis and Ascending Hierarchical Classification

Individual varieties are presented in Supplementary Figure 1, and physiological and agronomic variables are represented in Figure 4. Under water stress at the vegetative stage, the principal component analysis indicated that the PC1 explained most of the variation observed in the traits and accounted for 42.14% of the variation. The PC 1 comprises individual varieties KVX 61-1 (V2), IT 07K-188-49 (V15), and IT 99K-573- 2-1 (V5) on the lower right portion of the graph which is characterized by a strongly positive coordinate on the axis while varieties such as IT 84 S-2246-4 (V8), IT 99K-573-1-1 (V10), IT 07K-206-1-1 (V9), and IT 86 D-1033 (V18) on the lower left of the graph are characterized by a strongly negative coordinate on the axis (Figure 4). KVX 61–1 (V2), Gboto (V7), IT 07K-188-49 (V15), and IT 99K-573-2-1 (V5) are characterized by (i) high values for the variables _Fm, , RWC, Chl_C, and Fm and (ii) low values for the variable Fo while the individual varieties such as IT 84 S-2246-4 (V8), IT 99K-573-1-1 (V10), IT 07K-206-1-1 (V9), and IT 86 D-1033 (V18) are characterized by low values for the variables RWC, Chl_C, Har_T, and DFL.

The PC 2 particularly distinguishes individuals such as IT 86 D-888 (V6) V14, V4, and Kpodjiguegue (V19) on the upper side of the graph. The individual varieties such as IT 86 D-888 and IT 87 S-1390 are characterized by low values for the variables Fv_Fm and Fv_Fo, whereas Kpodjiguegue (V19) is characterized by high values for the variables number of pods and Yield.

Cluster 1 includes IT 86 D-888 (V6), IT 84 S-2246-4 (V8), IT 07K-206-1-1 (V9), IT 99 K—573-1-1 (V10), and IT 86 D-1033 (V18), which are characterized by (i) high values for the variable Fo and (ii) low values for the variables RWC, Fv_Fm, number of pods (No. P), Yield, Har_T, Chl_C, Fv_Fo, and days to flowering (DFL) (Supplementary Figure 2 and Figure 5). Cluster 2 comprises IT 87 S-1390 (V4), which is characterized by low values for the variables Fo and Fm (Supplementary Figure 2 and Figure 5). Cluster 3 is made of individuals such as KVX 61-1 (V2), IT 99K-573-2-1 (V5), Gboto (V7), IT 07K-188-49 (V15), and Kpodjiguegue (V19). This cluster is characterized by high values for the variables Fv, Fm, Yield, No. P, Fv_Fm, DFL, and Chl_C (Supplementary Figures 2 and Figure 5).

Under water stress at the reproductive stage, the principal component analysis indicated that the PC1 explained most of the variation observed in the traits and accounted for 54.28% (Supplementary Figure 3). Dimension 1 includes individuals such as IT 07K-188-49 (V15), Kpodjiguegue (V19), IT 06K-242-3 (V16), and KVX 61-1 (V2) on the right of the graph and individuals such as KVX 396-18 (V3), IT 07K-211-1-8 (V17), and IT 93K-452-1 (V13) on the left of the graph. IT 07K-188-49 (V15), Kpodjiguegue (V19), IT 06 K-24-3 (V16), and KVX 61-1 (V2) are characterized by high values for the variables DFL, Yield, Fv_Fo, Fv_Fm, Proli, Fv, RWC, and Fm while the individuals KVX 396-18 (V3), IT 07K-211-1-8 (V17), and IT 93K-452-1 (V13) are characterized by low values for the variables Fm, Fv, Fv_Fm, Har_T, Fv_Fo, RWC, and No.P (Figure 6).

Hierarchical tree and ascending hierarchical classification of the individuals under water stress at the reproductive stage are presented in Figure 7 and Supplementary Figure 4. Cluster 1 is made of individuals such as KVX 396-18 (V3), IT 93K-452-1 (V13), and IT 07K-211-1-8 (V17), which are characterized by low values for the variables Fm, Fv_Fm, Fv, Fv_Fo, Fo, and RWC. Cluster 2 comprises IT 87 S-1390 (V4), IT 86 D-888 (V6), IT 84 D-449 (V11), and Tawa (V20). This cluster is characterized by (i) high values for the variable Fo and (ii) low values for the variable Yield. Cluster 3 includes KVX 61-1 (V2), IT 07K-188-49 (V15), IT 06K-242-3 (V16), and Kpodjiguegue (V19). This group is characterized by high values for the variables Yield, DFL, Fv_Fo, RWC, Chl_C, Proli, No.P, Fv_Fm, and Fv.

4. Discussion

Water stress affects physiological parameters in vegetative and reproductive stages. Chlorophyll fluorescence represents the quantitative measure of the photosynthetic performance of plants under stress conditions. It is a good tool for the identification of genotypes tolerant to abiotic stress especially drought stress. Chlorophyll fluorescence measures the activity of photosystem II, which has a large measurement variability, i.e., water photolysis into gaseous oxygen and hydrogen ion, electron release, and production of energy compared to photosystem I [35, 36]. This study made it possible to compare the results of chlorophyll fluorescence obtained under water stress conditions and nonstressed conditions. The values of the chlorophyll fluorescence parameters recorded are lower under water deficit conditions than those under unstressed conditions as shown in Tables 1 and 2. There was a significant difference in the photosynthetic yield between the two regimes at vegetative and reproductive stages. This can be explained by the negative effect of water deficit on the photosynthetic apparatus of the genotypes studied. These results agree with those of Schreiber and Berry [37] and Al-Khatib and Paulsen [38] who indicated that the presence of stress causes a decrease in maximum fluorescence (Fm) and therefore a decrease in Fv/Fm, thus reducing photosynthesis. Bukhov and Mohanty [39], Guissé et al. [40], Lazár et al. [41], and Yamori et al. [42] reported that water stress causes a decline in the rate of O2 release, electron transport, and the Fv/Fm ratio which is explained by the increase of Fo due to the increase in leaf temperature and decrease in Fm. This corroborates our results. The reduction in the chlorophyll content in stressed plants can be explained by the fact that chlorophylls are more degraded than synthesized in the latter while in unstressed plants these pigments are synthesized. Moreover, according to Bousba et al. [43], the decline in chlorophyll content under water deficit is the consequence of the reduction in the opening of the stomata to avoid water losses through transpiration, and this reduction of stomata opening also leads to the decrease in the atmospheric CO2 necessary for photosynthesis.

The relative water content is a parameter indicating the water status of a plant. In the present study, water deficit had a significant effect on the RWC of stressed plants. The unavailability of water caused the stressed plants to lose their water and their turgidity which was illustrated through the lower values of the recorded RWC. Our results are in agreement with those on cowpea [4447]. Lobato et al. [45] studied the effect of the water deficit on cowpea and noticed that the RWC decreases by 25.7% while the control value remains between 88% and 91%. Anyia and Herzog [48] noticed an average reduction of 12.2% in leaf RWC in stressed plants compared to the control and that the magnitude of reduction was cultivar-dependent. Lobato et al. [46] observed a significant decrease of more than 31% in the RWC in two cowpea varieties under water stress. Meftah [47] noted a reduction of 14.41% under moderate stress and 23.03% under severe stress.

Under water stress conditions, the agronomic traits were significantly affected especially when there is water stress at the beginning of the inflorescence. These results confirm those of previous work by [48, 49] who showed that water stress negatively and significantly affected the flowering and filling of cowpea pods, which considerably reduced yield. The severity of this stress is because it affects parameters related to seed formation including photosynthesis and translocation of assimilates [4951]. Ahmed and Suliman [52] observed a significant reduction of over 40% under water deficit. The same result was obtained by Dadson et al. [29] in ten cultivars of cowpea under water constraints. Water stress affects the growth and development of plants and reduces seed yield. Lisar et al. [53] reported that, in plants under water stress, during the setting of flowers and pod formation stages, the yield reduction was 50% lower than the control plants. According to Chiulele et al. [54], drought during the formation of flowers and pods influences yields negatively. Water stress affects seed development, production, or distribution of photoassimilates [49]. Water deficit causes a delay in flowering depending on the species [55]. On the other hand, the stressed plants in the reproductive stage showed a similar trend in days to flowering with controlled plants [56]. Zombré et al. [57] observed that water deficit conditions cause a large number of flower abortions which reduces the yield of plants. Lalsaga et al. [56] explained the reduction in the weight of the seeds by a reduction in the synthesis of photoassimilates necessary for filling the seeds following the unavailability of water for photochemical reactions. This is similar to our observations during the two spells of water stress imposition. Thus, in our study, the more a plant had a long vegetative/development cycle especially the resistant ones [KVX 61–1 (V2), Kpodjiguegue (V19)], the more it showed a large number of pods and thus a high yield. This could also be due to the fact that with a long development cycle there was sufficient time for stressed plants recovery and the regaining of physiological functions indispensable for high productivity.

Proline is a very important osmoprotectant in the life of stressed plants. Proline can protect the integrity of proteins by preventing the aggregation of proteins and thus maintain the stability of the configuration of proteins and enzymes; its biosynthesis is activated during dehydration [58]. Our study showed an increase in the proline content in stressed plants, especially in stress-resistant varieties. Our results corroborate those obtained on beans [59] and cowpea [45, 60]. The increase in proline content and chlorophyll content results in keeping plant water status and chlorophyll activity.

The cowpea varieties fell into different clusters that equated their yield potentials and stress tolerance. Under drought imposition at the vegetative stage, it was observed that KVX 61-1 (V2), Gboto (V7), IT 07K-188-49 (V15), and IT 99 K-573-2-1 (V5) shared high values for Fv_Fm, Fv, RWC, Chl_C, and Fm as the most contributive variables for genotypes dispersion along the first axis. Similarly, in the second axis, IT 86 D-888 (V6) and Kpodjiguegue (V19) are characterized by a strongly positive coordinate on the axis. V6 was characterized by low values of Fv_Fm and Fv_Fo while Kpodjiguegue (V19) was characterized by high values for the variables number of pods and yield as the main contributive factor to drought tolerance in cowpea. The varieties KVX 61-1 (V2), Gboto (V7), IT 07K-188-49 (V15), IT 99K-573-2-1 (V5), IT 06K-242-3 (V16), and Kpodjiguegue (V19) were associated with Fv_Fm, Fv, RWC, Chl_C, Fm, and yield, which represent the cluster of higher-yielding and stress-tolerant varieties under water deficit at vegetative stage. Santos et al. [61] reported a strong contribution of the pod and grain-related components, namely, pod dry weight, along the PC1 axis, for the screening of drought-tolerant cowpea genotypes at the seedling stage. Under water deficit at reproductive stage, IT 07 K-188-49 (V15), Kpodjiguegue (V19), IT 06K-242-3 (V16), and KVX 61–1 (V2) were characterized by high values for DFL, Yield, Fv_Fo, Fv_Fm, Proli, Fv, RWC, and Fm. Thus, high yield, high proline content, high relative water content, Fv/Fm, and days to flowering significantly contributed to the drought tolerance of these varieties used in this study. The drought stress attributes such as yield, proline content, relative water content, and Fv/Fm were correlated with the varieties KVX 61-1 (V2), IT 07K-188-49 (V15), IT 06K-242-3 (V16), IT 99K-573-2-1 (V5), and Kpodjiguegue (V19), which represent the cluster of higher-yielding and stress-tolerant varieties under water deficit at reproductive stage. Overall, the dispersion of cowpea varieties in the PCA reflected the groupings based on their physiological and agronomic responses to drought stress effects.

5. Conclusion

The varieties Kpodjiguegue, Gboto, and KVX 61-1 gave the highest grain yields and were found to be the most drought-tolerant under water stress at the vegetative stage and the varieties Kpodjiguegue, IT 06K-242-3, IT 07K-188–49, and KVX 61-1 gave the highest grain yields and were found to be the most drought-tolerant under water deficit at the reproductive stage. To have more efficient genotypes in the condition of water stress, we must initiate crosses between the different tolerant genotypes and between tolerant genotypes and those sensitive to water stress to combine the different mechanisms of tolerance.

Abbreviations

CCI:Chlorophyll concentration index
ANOVA:Analysis of variance
PCA:Principal component analysis
RWC:Relative water content
FW:Fresh weight
DW:Dry weight
WFC:Weight of the full turgor
Fv/Fm:Photochemical yield.

Data Availability

Climatic, physiological, morphological, and agronomic data collected as the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

VE and AA designed the experiment. CT carried out the fieldwork. CT and IBC conducted laboratory work. CT, VE, and AA performed the statistical analysis and CT drafted the manuscript. VE, AA, and IBC coordinated the study and contributed to the write-up of the manuscript. VE and AA edited the manuscript. All authors read and approved the final manuscript.

Supplementary Materials

Supplementary Figure 1: individual varieties under drought stress at the vegetative stage. The labeled individuals are those with the higher contribution to the plane construction. Supplementary Figure 2: ascending hierarchical classification of the individuals under drought stress at the vegetative stage. Supplementary Figure 3: individual varieties under drought stress at the reproductive stage. The labeled individuals are those with the higher contribution to the plane construction. Supplementary Figure 4: ascending hierarchical classification of the individuals under water stress at the reproductive stage. (Supplementary Materials)