Next Article in Journal
Effects of Particle Size Distribution on the Physicochemical, Functional, and Structural Properties of Alfalfa Leaf Powder
Previous Article in Journal
Current Flaxseed Dehulling Technology in China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Variability and Interrelationships of Grain, Cooking, and Nutritional Quality Traits in Cowpea: Implications for Cowpea Improvement

by
Michael M. Chipeta
1,*,
Esnart Yohane
2,
John Kafwambira
1 and
Jessica Kampanje-Phiri
3
1
Department of Crops and Soil Sciences, Lilongwe University of Agriculture and Natural Resources, P.O. Box 219, Lilongwe, Malawi
2
Department of Agricultural Research Services, Chitedze Research Station, P.O. Box 128, Lilongwe, Malawi
3
Department of Human Ecology, Lilongwe University of Agriculture and Natural Resources, P.O. Box 219, Lilongwe, Malawi
*
Author to whom correspondence should be addressed.
Agriculture 2024, 14(4), 633; https://doi.org/10.3390/agriculture14040633
Submission received: 14 February 2024 / Revised: 5 March 2024 / Accepted: 8 March 2024 / Published: 19 April 2024

Abstract

:
Grain quality, cooking quality, and nutritional quality traits are some of the major attributes that enhance the uptake and utilization of improved cowpea varieties. Therefore, there is a need for a better understanding of the genetic variation and inter-relationships among these quality traits in cowpeas to integrate them into cowpea breeding programs. This study was conducted to determine genetic variability among 306 cowpea genotypes for grain quality, cooking quality, and nutritional quality traits and to understand the interrelationships among these traits for exploitation in breeding programs. The results showed highly significant differences (p < 0.001) among genotypes for grain quality, cooking quality, and nutritional quality traits. The mean performance for these quality traits was also very variable. These results suggest that genetic variability exists in the cowpea genotypes studied, which can be exploited in breeding programs aimed at developing high-performing varieties for the said traits. Significant (p < 0.001) positive correlations were detected for protein content with iron and zinc. On the other hand, nutritional quality traits did not exhibit any association with grain quality or cooking quality traits. Cooking quality traits were also shown to be significantly and positively correlated with grain quality traits. This study has identified several genotypes with desirable quality-related traits that could be used in crossing programs to generate improved varieties with consumer-preferred traits to improve the food, income, and nutritional status of many smallholder farmers that largely depend on cowpeas.

1. Introduction

Cowpea (Vigna unguiculata (L.) Walp.) is an important legume crop produced and consumed in Sub-Saharan Africa by the majority of smallholder farmers, especially women. Cowpea is diploid, with 2n = 22 and a genome size of about 620 Mb [1,2,3]. The crop is autogamous, but outcrossing of up to 5% has been reported in the cultivated varieties, mainly due to insect activities [4]. It is a crucial food crop for meeting dietary protein needs and preventing micronutrient deficiencies [5], especially among women and children [6].
Cowpea whole grains contain 23–32% protein, 50–60% carbohydrates, and 1% fat [7], as well as considerable amounts of phytochemicals, dietary fiber, minerals (calcium, iron, and zinc), and vitamins. Iron and zinc are essential for human well-being. Their adequate supply helps prevent anemia and boost the immune system, two common problems in developing countries [8]. The composition of nutrients can vary due to several factors, including varietal differences, climatic conditions, and agronomic practices [9]. Some studies have reported mature grains to contain lower concentrations of most minerals [10] than immature green pods [11] and leaves [12].
Despite its enormous economic potential, the productivity of the crop is considerably low in Malawi, which is estimated at an average of less than 500 kg/ha among smallholder farmers [13]. This low productivity is attributed to low adoption levels of improved varieties (currently at 10%) and continued use of landraces, which tend to be low-yielding and highly susceptible to biotic and abiotic stress. Improved varieties, however, have been blamed for lacking important quality traits such as grain size, cooking time, and broth thickness. Low adoption of improved varieties and a lack of farmer and consumer-preferred varieties have been attributed to the failure of national breeding programs to involve farmers in the process of designing and developing varieties to meet their priorities and preferences [14].
Although cowpeas are nutritionally rich, the longer cooking time makes the legume unacceptable to some consumers. In a trait preference study conducted in Malawi in 2021 [15], short cooking time and broth thickness were among the priority traits for consumers to adopt and utilize a particular cowpea variety. It has also been reported that broth thickness is an important trait in the determination of cooking quality [16]. Traore et al. [17] stated that cooking qualities must be combined with other consumer preferences to encourage cowpea utilization. Other traits that influence cowpea marketing and utilization include seed coat color, seed shape, and grain size [18,19,20,21,22].
Considerable genetic variability has been reported for cowpea’s physical, cooking, and nutritional quality traits [19,20,23,24,25,26]. In Malawi, reports on the genetic variability for grain quality (grain color, seed size, seed shape), cooking quality (cooking time, broth thickness), and nutritional quality (protein, zinc, and iron) traits on a diverse panel of cowpea genotypes are rare and non-existent. The degree of association among these quality traits has neither been fully dissected nor exploited to inform cowpea breeding programs. The present study was therefore conducted to evaluate the genetic variability of cowpea genotypes for grain, cooking, and nutritional quality traits. Further, this study was conducted to understand the degree of associations among the cowpea quality traits so that they can be exploited in a breeding program to develop cowpea varieties with consumer-preferred, market-driven, and resilient inclusive traits to enhance the food, income, and nutrition security of smallholder farmers in East and Southern Africa.

2. Materials and Methods

2.1. Study Area Description

A field experiment was conducted at Lilongwe University of Agriculture and Natural Resources (LUANAR), Bunda College of Agriculture, in the 2021/2022 growing season. The site is located at about 14°11′ S latitude and 33°46′ E longitude. The site receives a mean annual rainfall of about 930 mm (Figure 1) with an average temperature range of 17.2 °C to 19.1 °C and predominantly clay-loam soils.

2.2. Plant Materials

A diversity panel constituting 306 genotypes (advanced breeding lines, released varieties, and landraces) was used for this study (Supplementary Table S1). These genotypes were sourced from the International Institute of Tropical Agriculture-Nigeria genebank (IITA-44), India (61), Mozambique (24), the United States of America (22), Tanzania (16), Nigeria (13), South Africa (9), Uganda (4), Zambia (4), Hungary (2), Ghana (2), Cameroon (2), Botswana (2), Senegal (2), Russia (2), Argentina (1), Benin (1), Italy (1), Niger (1), Zimbabwe (1), Malawi (77), and 14 genotypes had unknown origins. These genotypes represent a subsample of the Centre of Innovation for Crop Improvement for East and Southern Africa (CICI-ESA) cowpea germplasm at LUANAR.

2.3. Experimental Design and Field Layout

An alpha lattice design with three replications generated using a Breeding Management System (BMS) was used for field establishment on 18 January 2022. The plot size was 2 ridges, 2 m long, each spaced 75 cm apart. Two seeds were planted per station at a depth of 3 cm, spaced 20 cm apart. Weeds were controlled by manual hoeing and uprooting to ensure weed-free conditions throughout the season. The infestation of insect pests such as aphids was controlled by the weekly application of Snowcron 500EC, a broad-spectrum emulsifiable concentrate insecticide, three weeks after germination. Data on grain quality traits were collected using the Field Book App-v5.6 [27], installed on mobile tablets. For laboratory analysis, graded and clean grains of each genotype were sampled from the three replications in the field.

2.4. Determination of Grain Quality Traits

Data for qualitative (seed shape, seed coat color) and quantitative traits (100 seed weight, seed length, seed width) of grain quality were collected following the cowpea descriptor [28]. Seed length and width were measured in millimeters on five randomly selected seeds using a ruler, and the weight of 100 seeds was determined as the weight in grams for 100 randomly selected seeds using an electronic balance.

2.5. Determination of Cooking Quality Traits

2.5.1. Cooking Time

The cooking time of each dry cowpea genotype was determined according to the method of Akinyele et al. [29]. Briefly, 20 g of dry cowpea grains were added to 450 mL of boiling water (tap water) in an aluminum cooking pot (without a top cover) on an electric hotplate. The water level was maintained by adding more boiled water intermittently. An electric kettle was used to heat the extra water, which was used for re-filling the pots during cooking. Cooking heat intensity was also maintained by using the maximum heating level of the hotplate for all cooking trials. Cooking time logging began when the water returned to the boiling point after immersing the grains in the boiling water. The cowpeas were left to cook for 30 min; thereafter, a doneness test was conducted every 5 min using the tactile method [30,31] by pressing the cooked seeds between two fingers until they were mushy, i.e., no hard material was found, just as also conducted by Munthali et al. [32] for common beans.

2.5.2. Broth Viscosity

The cooked cowpeas were removed from the pot, and then all broth and residual solids were transferred into a jar and mixed with boiling water to a total volume of 365 mL. The viscosity (cP or mPa) of 365 mL samples at 25 °C was measured using a programmable Brookfield digital viscometer (Model LV DV-I, Brookfield Engineering Laboratories, Inc., Stoughton, MA, USA) fitted with an LV spindle No. 1 as recommended in the Brookfield Labs Inc. guide (Brookfield Manual No. M14-023). A 400 mL open jar was used as a sample holder. Viscosity readings were taken after one minute at a spindle speed of 60 rpm [33]. All samples were thoroughly mixed just before taking viscosity measurements.

2.6. Determination of Nutritional Quality Traits

2.6.1. Digestion Procedure

Sample digestion was conducted according to the method proposed by Sahrawat et al. [34], with modifications. Approximately 0.2 g of finely ground cowpea seed samples were transferred to 50 mL digestion tubes, and six blanks were prepared to be used for standards. Thereafter, 2.5 mL of a sulfuric acid/selenium mixture was added to each digestion tube and blank. The sulfuric acid/selenium mixture was prepared by dissolving 3.5 g of selenium powder in 1 L of sulfuric acid and heating the mixture on a hot plate at high temperature with occasional stirring using a glass rod until the mixture turned clear and was cooled. The digestion tubes containing samples and the sulfuric acid/selenium mixture were placed in an aluminum block on a hotplate and heated to 200 °C until sample fumes were seen. The tubes were then removed from the hotplate and allowed to cool for 10 min. Then, 1 mL and 2 mL of 30% hydrogen peroxide were added to the samples and standards. The samples were replaced on the hotplate with a heavy glass vial on top of each tube and heated to 330 °C until the mixture turned clear and colorless, indicating the completion of digestion. The mixture was then allowed to cool.
N stock solution (0, 0.2, 0.4, 0.6, 0.8, and 1 mL, respectively) was added to the standards. The N stock solution was prepared by diluting 4.714 g of oven-dried ammonium sulfate in a 100 mL volumetric flask to make 10,000 ppm N. The standards are equivalent to 0, 1, 2, 3, 4, and 5% N in the plant digests.

2.6.2. Determination of Protein

The colorimetry method was used for the determination of nitrogen [35]. The N1 reagent was prepared by dissolving 68 g sodium salicylate, 50 g sodium citrate (Tri-Sodium Citrate), and 50 g sodium tartrate in 500 mL deionized water. This solution dissolved 0.24 g of sodium nitroprusside and was diluted to 2000 mL. The N2 reagent was prepared by dissolving 60 g sodium hydroxide in about 1500 mL of deionized water. The solution was then mixed with 28.5 mL of a 3.5% sodium hypochlorite solution and diluted to 2000 mL.
To determine the total nitrogen in samples, 0.750 mL extracts from the digestion above were put into 25 mL glass vials, and to each of the vials containing the extracts and standards, 5 mL of N1 solution was added. After 5 min, 5 mL of the N2 solution was added. The vials were allowed to stand until the color developed (about 1 h), and then the absorbance of the samples and standards was read on a spectrophotometer at 655 nm. To obtain the protein percentage, the percentage of nitrogen was multiplied by 6.25.

2.6.3. Determination of Iron and Zinc

The nitrogen digestion method was used to determine iron and zinc contents. To prepare stock solutions, 0.498 g dried ferrous sulfate septahydrate (FeSO4·7H2O) and 0.440 g dried zinc sulfate septahydrate (ZnSO4·7H2O) were weighed for iron and zinc determination, respectively. These were then dissolved in about 200 mL dilute HCl (0.1 N) and made up to the mark with distilled water in a 1 L volumetric flask to contain 100 ppm Fe and 100 ppm Zn stock solutions, respectively. The original samples of the digestion extract were passed on the Atomic Absorption Spectrophotometer (AAS).

2.7. Data Analysis

An exploratory analysis using density plots was performed on quantitative data to check for ANOVA assumptions. Thereafter, all the analysis was performed in the R software environment [36]. All laboratory analyses were performed in duplicate. ANOVA was performed by fitting the models with the aov function of the agricolae package [37] in R, and Duncan’s multiple range test was used to separate means at a 0.05 significance level. Principal component analysis (PCA) was performed using the FactoMineR package [38], and Pearson correlation was performed using the corr. function of the stats package in R. The Euclidean distance metric was used in hierarchical cluster analysis, and the associations between groups were performed using the Ward D2 method in a cluster R package version 2.1.6 [39].

3. Results

3.1. Frequency Distribution of Quantitative Traits

The exploratory analysis showed that all the traits except cooking time, iron content, and zinc content exhibited a nearly normal distribution pattern (Figure 2).

3.2. Genetic Variability and Mean Performance for Grain Quality Traits

Highly significant differences (p < 0.001) were observed among the genotypes for all grain quality traits (Table 1 and Table S1). The seed length ranged from 3.73 mm (TVu-322) to 10.33 mm (MWcp03). The highest seed width was recorded on genotype MWcp44 (7.13 mm), while genotype TVu-322 (2.93 mm) registered the lowest seed width. The 100 seed weights ranged from 6.0 g to 26.13 g. Genotype MWcp36 had the largest seed size, while genotype TVu-10169 had the smallest seed size. In terms of seed coat color, genotypes were assorted into brown (43.14%), white (19.28%), reddish-brown (13.40%), purple (11.76%), purplish-brown (4.25%), black (1.63%), and red (1.31%). Some of the genotypes (5.23%) were either mottled or had other mixed colors. For seed shape, the ovoid shape was dominant with about 60.78%, followed by the rhomboid shape (32.03%) and the kidney shape (6.54%). Globose and crowder shapes were observed in 0.66% of the genotypes (Table 2 and Table S1).

3.3. Genetic Variability and Mean Performance for Cooking Quality Traits

Highly significant differences (p < 0.001) were observed among genotypes for cooking quality traits, cooking time, and broth viscosity (Table 3). Genotype TVu-15631 recorded the shortest cooking time (39 min), while genotype MWcp07 reported the longest cooking time (191 min) (Table 4 and Table S1). The minimum broth viscosity at 25 °C was recorded at 1.6 cP for IT10K-834-3, whereas the maximum was observed for TVu-3252 (5.5 cP).

3.4. Genetic Variability and Mean Performance for Nutritional Quality Traits

The genotypes exhibited highly significant differences (p < 0.001) for both protein (%) and mineral (iron and zinc) contents (mg/kg) (Table 5). The highest average protein content was recorded on genotype TVu-3243 (22.92%), while the lowest protein content was reported on genotype TVu-3263 (3.6%) (Table 6 and Table S1). The highest iron content (21.44 mg/kg) was recorded on genotype TVu-3533, whereas the lowest level (0.12 mg/kg) was reported on genotype MWcp26. Genotype MWcp37 had the highest zinc level (2.63 mg/kg), while Raha 1 and TVu-1472 had the lowest zinc level (0.08 mg/kg).

3.5. Correlation Coefficients among Grain, Cooking, and Nutritional Quality Traits of Cowpea Genotypes

The relationships among the grain quality, cooking quality, and nutritional quality traits are presented in Table 7. Protein content was significantly positively correlated with iron content (r = 0.28, p < 0.001) and zinc content (r = 0.21, p < 0.001). Iron content was positively correlated with zinc content (r = 0.24, p < 0.001). Cooking time was significantly and positively correlated with seed length (r = 0.21, p < 0.001), seed width (r = 0.26, p < 0.001), and 100 seed weight (r = 0.19, p < 0.001). Highly significant positive correlations were also observed for seed length and seed width with 100 seed weights (r = 0.71 and r = 0.75 at p < 0.001, respectively).

3.6. Genetic Diversity among the Cowpea Genotypes

3.6.1. Principal Component Analysis for Grain Quality, Cooking Quality, and Nutritional Quality Traits

Table 8 shows principal component analyses (PCA) for the quantitative quality traits among the genotypes. The genotypes were plotted on two dimensions based on the PCA results (Figure 3). All the traits were grouped into three principal components (PC) (eigenvalues ≥ 1), which accounted for 65.17% of the variability. The first PC explained about 33.36% of the total variability present and was mainly associated with seed width, seed length, and 100 seed weight. The second PC accounted for 18.81% of the total variability, which was correlated with iron, zinc, and protein contents. The third PC explained nearly 12.98% of the total variation and was mainly associated with broth viscosity.

3.6.2. Hierarchical Clustering of Genotypes Based on Grain Quality, Cooking Quality and Nutritional Quality Traits

As shown in the dendrogram (Supplementary Figure S1), cluster analysis grouped the genotypes into two major clusters based on the quality traits. Cluster I comprised 234 genotypes (gold color), while Cluster II comprised 72 genotypes (blue color). The genotypes in Cluster I were relatively small-seeded with shorter cooking times, while those in Cluster II were generally large-seeded with longer cooking times and higher protein and iron content (Table 9).

4. Discussion

Most resource-constrained smallholder farmers in sub-Saharan Africa, especially women, rely on cowpeas as their major source of protein and other mineral elements such as iron and zinc, thereby contributing to their improved nutritional status and health well-being. In dryland areas, cowpea is a major source of income to support their daily livelihoods. At the same time, cowpea is one of the crops used to adapt to climate-related changes due to its resilience in these areas. Current cowpea breeding strategies in the region must, therefore, take into consideration its grain, nutritional value, and cooking quality, which influence farmers/consumers’ choices and consumption.
This study observed significant variations in the genotypes for grain quality, cooking quality, and nutritional quality traits (Table 1, Table 3 and Table 5). This suggests that genetic variability exists in the current set of cowpea germplasm in Malawi, which could be exploited in breeding programs aimed at developing high-performing varieties for the said traits. Assessment of genetic variability for different traits is an important pre-breeding step as it provides an opportunity for plant breeders to develop new and improved varieties with desirable characteristics that are both farmer- and consumer-preferred. Earlier studies [6,40,41] reported wide genetic variation in grain quality, crude protein, and micronutrient contents in cowpeas, suggesting that the nutritional value of varieties could be improved.
In the study, seed coat color among the genotypes was quite variable, predominated by brown color (43.14%), followed by white (16.39%), and reddish-brown (13.40%). In cowpeas, seed coat color is an important aspect of consumer preference, which directly influences the marketability of the grain. This implies that variety development programs need to respond to specific and unique market preferences for improved uptake of the varieties. Trait preference studies in different regions and countries reported these unique preferences [19,21,22,42,43]. For instance, in Malawi, Mozambique, and Tanzania, through a trait preference survey [15] among farmers, consumers, traders, and processors, seed coat color was among the most preferred characteristics.
The genotypes in the study also varied for seed shape, with over 51.67% being ovoid and 32.03% rhomboid. This finding is similar to Kim et al. [44], who reported that out of the 245 cowpea genotypes in the study, 66.9% were egg-shaped (ovoid) and 24.9% were rectangular (rhomboid). This suggests that the ovoid shape is a dominant seed shape in cowpeas.
The other important trait in cowpea breeding is grain size. Seed size is directly correlated with grain yield, and as such, it is considered an essential market trait [18]. Seed size is measured as grams per 100 seeds, and Moses and Zibokere [24] grouped seed size into small (10–15 g), medium (15.1–20 g), large (20.1–25), and very large (greater than 25 g). In the present study, genotypes exhibited wide genetic variability for seed size, which ranged from 6 to 26.13 g. This range is within most of the reported seed sizes [23,26]. It was noted that most of the landraces were larger-seeded than the improved varieties. As noted before, most farmers prefer large-seeded varieties, and this might explain the reasons behind the low adoption levels of improved varieties in Malawi since all the improved varieties have small grain sizes.
With cowpeas being an important source of protein and other elements for many smallholder farmers, breeding programs must continue to develop and release more nutritious varieties. The genetic variability for nutritional quality properties such as protein, zinc, and iron exhibited in the study is a sure way to set the trajectory towards sustainable, improved variety development that meets the nutritional needs of farmers, especially women and children. Genetic variability for nutritional quality in cowpeas is a common phenomenon [6,12,20,40,41], which means that selection for these traits in cowpeas can increase genetic gain.
Cooking time is usually associated with energy sources. In sub-Saharan Africa, where the predominant energy source is fuel wood, most households struggle to prepare meals that require more energy. Shorter cooking varieties are usually a strategy for energy savings. The findings of the study established some considerable genetic variation among genotypes, with some taking a minimum of 39 min to cook. Some studies have reported a minimum of 35 min [45] to cook. Contrastingly, some genotypes in the study took 191 min, which suggests that they are less economical in terms of energy consumption. Related to cooking quality is broth viscosity, which very few studies have looked at in cowpeas as a genetically controlled trait. This study suggests that the variability recorded is genetically controlled, as considerable variation was observed in the genotypes, which ranged from 1.6 to 5.5 cP. In Malawi, it has been reported that broth thickness is one of the key determinants of the cooking quality of the legume [16], and if a variety lacks this trait, it can rarely be taken up by farmers/consumers.
Significant positive correlations were detected among the 306 genotypes for protein content with iron and zinc (Table 7). This is in agreement with the positive correlation between crude protein and Fe contents in 11 genotypes reported by Moura et al. [25] and Boukar et al. [6,40]. Since the presence of genetic variability has already been elucidated for these traits, it is, therefore, possible to improve protein, iron, and zinc contents without adverse interactions, as also espoused by Nielsen et al. [46] and Jean Baptiste et al. [47], who indicated some heritability for crude protein content in cowpea.
On the other hand, nutritional quality traits did not exhibit any association with grain quality or cooking quality traits. The non-significant associations between grain quality and nutritional quality traits suggest that desirable traits can be incorporated through trait introgression techniques. This is in agreement with several reported [25,41,48] results where cooking time was not associated with protein content in their study. However, cooking quality traits were shown to be significantly and positively correlated with grain quality traits. For example, cooking time is correlated with seed length, seed width, and 100 seed weight (cooking time increases with seed size). The implication is that as breeding programs strive to develop varieties that are large-seeded with short cooking times to meet market demands, there will be a trade-off between these two traits. There is a need to devise breeding strategies that maximize both of these traits, as well as a careful selection of parental materials. In addition to this, there is a need to dissect the genetic mechanisms of these traits and determine whether any of them can be amenable to trait integration techniques.
A genotype by trait biplot was produced using the first two PCs, which explained 52.19% of the variability (Figure 3). This was used to identify the best-performing and genetically dissimilar genotypes that could be further explored for breeding. The random scattering of the cowpea genotypes across the quadrant, as seen in the biplot, suggests that the genotypes were genetically different. The genotypes clustered based on their performance for the traits; thus, those genotypes located close together had relatively similar values for specific traits. The biplot shows that many genotypes clustered near the origin, so they are likely derived from the same parents. Genotypes located further from the biplot origin had more extreme values for a specific trait than genotypes closer to the origin. For example, genotypes MWcp03 and MWcp36 showed extreme values for seed length, width, and 100 seed weight, whereas TVu-3533 showed higher values for protein, iron, and zinc concentrations. On the other hand, TVu-322 had the highest broth viscosity, and MWcp46 had the longest cooking time. Furthermore, MWcp37 was distant from the rest of the genotypes, indicating that the genotype is genetically dissimilar from most of the genotypes. Concentrations of protein, iron, and zinc were independent of broth viscosity, cooking time, seed length, width, and 100 seed weight. This suggests that protein, iron, and zinc can be improved without altering the rest of the traits. Noticeably, broth viscosity was not associated with cooking time or the quantitative grain quality traits. However, cooking time was positively related to the quantitative grain quality traits. This result further suggests that increasing the seed size in cowpeas may increase the cooking time but not the broth viscosity. As already elucidated, the landraces from Malawi were generally large-seeded with longer cooking times, but these tend to be common among smallholder farmers. This was further confirmed through cluster analysis, which clustered the majority of the landraces from Malawi into Cluster II, which had large-seeded genotypes with relatively higher protein and iron content. Nevertheless, genotypes in Cluster I were small-seeded, and within the cluster were most of the improved varieties and advanced lines from the IITA, which is an indication that they were selected for similar traits, mostly small seed size, which is contrary to farmer preferences.

5. Conclusions

This study has identified genetic variability in grain quality traits, cooking quality traits, and nutritional quality traits among the cowpea genotypes that constitute breeding lines and landraces. This study has further identified relationships among the quality traits. Significant and positive correlations among protein, iron, and zinc contents suggest the possibility of improving the concentrations of these nutrients simultaneously. Non-significant associations between grain quality and nutritional quality traits suggest that desirable traits can be incorporated through trait introgression techniques. This study has also identified several potential genotypes with desirable quality-related traits that could be used in crossing programs to generate improved varieties with consumer-preferred traits to improve the food, income, and nutritional status of many smallholder farmers that largely depend on cowpeas.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/agriculture14040633/s1, Table S1: List of genotypes used in the study, Analysis of variance, summary statistics and Genotype mean values for grain, cooking, and nutritional quality traits; Figure S1. Hierarchical clustering dendrogram analysis.

Author Contributions

Conceptualization, M.M.C. and J.K.-P.; data curation, J.K.; formal analysis, J.K.; funding acquisition, M.M.C.; investigation, M.M.C. and E.Y.; methodology, M.M.C. and J.K.; project administration, J.K.; resources, M.M.C.; supervision, M.M.C., E.Y. and J.K.-P.; writing—original draft preparation, J.K.; writing—review and editing, M.M.C., E.Y. and J.K.-P. All authors have read and agreed to the published version of the manuscript.

Funding

This study was made possible by the generous support of the American people through the United States Agency of International Development (USAID) under the Feed the Future Initiative. The contents are the responsibility of the authors and do not necessarily reflect the views of USAID or the United States Government. Program activities are funded by the United States Agency for International Development (USAID) under Cooperative Agreement No. 7200AA-19LE-00005.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Data is contained within the article and supplementary material.

Acknowledgments

We thank the Centre of Innovation for Crop Improvement for East and Southern Africa (CICI-ESA) research team and Lilongwe University of Agriculture and Natural Resources for the smooth implementation of project activities. Department of Agriculture Research Services (DARS) through Chitedze Research Station, Instituto de Investigação Agrária de Moçambique (IIAM), and Tanzania Agricultural Research Institute (TARI) for sharing with us their cowpea germplasm.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study, in the collection, analysis, or interpretation of data, in the writing of the manuscript, or in the decision to publish the results.

References

  1. Boukar, O.; Belko, N.; Chamarthi, S.; Togola, A.; Batieno, J.; Owusu, E.; Haruna, M.; Diallo, S.; Umar, M.L.; Olufajo, O.; et al. Cowpea (Vigna unguiculata): Genetics, genomics and breeding. Plant Breed. 2018, 138, 415–424. [Google Scholar] [CrossRef]
  2. Lonardi, S.; Muñoz-Amatriain, M.; Liang, Q.; Shu, S.; Wanamaker, S.; Lo, S.L.; Tanskanen, J.; Schulman, A.H.; Zhu, T.; Luo, M.C.; et al. The genome of cowpea (Vigna unguiculata [L.] Walp.). Plant J. Cell Mol. Biol. 2019, 98, 767–782. [Google Scholar] [CrossRef]
  3. Timko, M.P.; Rushton, P.J.; Laudeman, T.W.; Bokowiec, M.T.; Chipumuro, E.; Cheung, F.; Town, C.D.; Chen, X. Sequencing and analysis of the gene-rich space of cowpea. BMC Genom. 2008, 9, 103. [Google Scholar] [CrossRef]
  4. Badiane, F.A.; Diouf, M.; Diouf, D. Cowpea. In Broadening the Genetic Base of Grain Legumes; Singh, M., Bisht, I., Dutta, M., Eds.; Springer: New Delhi, India, 2014. [Google Scholar]
  5. Dakora, F.D.; Belane, A.K. Evaluation of protein and micronutrient levels in edible cowpea (Vigna unguiculata L. Walp.). Front. Sustain. Food Syst. 2019, 3, 70. [Google Scholar] [CrossRef]
  6. Boukar, O.; Massawe, F.; Muranaka, S.; Franco, J.; Maziya-Dixon, B.; Singh, B.; Fatokun, C. Evaluation of cowpea germplasm lines for protein and mineral concentrations in grains. Plant Genet. Resour. 2011, 9, 515–522. [Google Scholar] [CrossRef]
  7. Kirse, A.; Karklina, D. Integrated evaluation of cowpea (Vigna unguiculata (L.) Walp.) and maple pea (Pisum sativum var. Arvense L.) spreads. Agron. Res. 2015, 13, 956–968. [Google Scholar]
  8. Blair, M.; Astudillo, C.; Grusak, M.; Graham, R.; Beebe, S. Inheritance of seed iron and zinc content in common bean (Phaseolus vulgaris L.). Mol. Breed. 2009, 23, 197–207. [Google Scholar] [CrossRef]
  9. Affrifah, N.S.; Phillips, R.D.; Saalia, F.K. Cowpeas: Nutritional profile, processing methods and products—A review. Legume Sci. 2022, 4, e131. [Google Scholar] [CrossRef]
  10. Belane, A.; Dakora, F. Elevated concentrations of dietarily-important trace elements and macronutrients in edible leaves and grain of 27 cowpea (Vigna unguiculata L. Walp.) genotypes: Implications for human nutrition and health. Food Nutr. Sci. 2012, 3, 377–386. [Google Scholar] [CrossRef]
  11. Gerrano, A.S.; van Rensburg, W.S.J.; Adebola, P.O. Nutritional composition of immature pods in selected cowpea (Vigna unguiculata’(L.) Walp.) genotypes in South Africa. Aust. J. Crop Sci. 2017, 11, 134–141. [Google Scholar] [CrossRef]
  12. Gerrano, A.; Adebola, P.; van Rensburg, W.J.; Venter, S. Genetic variability and heritability estimates of nutritional composition in the leaves of selected cowpea genotypes [Vigna unguiculata (L.) Walp.]. HortScience A Publ. Am. Soc. Hortic. Sci. 2015, 50, 1435–1440. [Google Scholar] [CrossRef]
  13. FAOSTAT. Statistical Database. Food and Agriculture Organization of the United Nations. 2023. Available online: https://www.fao.org/faostat/en/#data/QCL (accessed on 1 March 2023).
  14. Coulibaly, O.; Alene, A.D.; Abdoulaye, T.; Chianu, C.; Manyong, V.; Aitchedji, C.; Fatokun, D.; Kamara, A.; Ousmane, B.; Tefera, H. Baseline Assessment of Cowpea Breeding and Seed Delivery Efforts to Enhance Poverty Impacts in Sub-Saharan Africa; ICRISAT: Telengana, India, 2010. [Google Scholar]
  15. Chipeta, M.M.; Kampanje-Phiri, J.; Moyo, D.; Colial, H.; Tamba, M.; Belarmino, D.; Hella, J.; Yohane, E.; Mvula, N.; Kafwambira, J. Understanding specific gender dynamics in the cowpea value chain for key traits to inform cowpea breeding programs in Malawi, Mozambique and Tanzania. Front. Sociol. 2024, 9, 1254292. [Google Scholar] [CrossRef] [PubMed]
  16. Nkongolo, K.; Bokosi, J.; Malusi, M.; Vokhiwa, Z.; Mphepo, M. Agronomic, culinary, and genetic characterization of selected cowpea elite lines using farmers’ and breeder’s knowledge: A case study from Malawi. Afr. J. Plant Sci. 2009, 3, 147–156. [Google Scholar]
  17. Traore, K.; Sawadogo, P.; Batiéno, T.B.J.; Sawadogo, N.; Zongo, H.; Poda, S.L.; Tignegré, J.B.D.S.; Ouédraogo, T.J.; Sawadogo, M. Physical and cooking characteristics of six cowpea varieties seeds cultivated in Burkina Faso. Am. J. Plant Sci. 2022, 13, 7. [Google Scholar] [CrossRef]
  18. Egbadzor, K.F.; Dadoza, M.; Danquah, E.; Yeboah, M.; Offei, S.; Ofori, K. Genetic control of seed size in cowpea (Vigna unguiculata (L.) Walp). Int. J. Agric. Sci. 2013, 5, 367–371. [Google Scholar] [CrossRef]
  19. Egbadzor, K.F.; Danquah, E.Y.; Ofori, K.; Yeboah, M.; Offei, S.K. Diversity in 118 cowpea [Vigna unguiculate (L.) Walp] accessions assessed with 16 morphological traits. Int. J. Plant Breed. Genet. 2014, 8, 13–24. [Google Scholar] [CrossRef]
  20. Gondwe, T.; Alamu, E.O.; Mdziniso, P.; Maziya-Dixon, B.B. Cowpea (Vigna unguiculata (L.) Walp) for food security: An evaluation of end-user traits of improved varieties in Swaziland. Sci. Rep. 2019, 9, 15991. [Google Scholar] [CrossRef]
  21. Hella, J.P.; Chilongo, T.; Mbwag, A.M.; Bokosi, J.; Kabambe, V.; Riches, C.; Massawe, C.L. Participatory market-led cowpea breeding in Sub- Saharan Africa: Evidence pathway from Malawi and Tanzania. Merit Res. J. Agric. Sci. Soil Sci. 2013, 1, 11–18. [Google Scholar]
  22. Herniter, I.A.; Lo, R.; Muñoz-Amatriaín, M.; Lo, S.; Guo, Y.-N.; Huynh, B.-L.; Lucas, M.; Jia, Z.; Roberts, P.A.; Lonardi, S.; et al. coat pattern QTL and development in cowpea (Vigna unguiculata [L.] Walp.). Front. Plant Sci. 2019, 10, 1346. [Google Scholar] [CrossRef]
  23. Kuldeep Tripathi, P.G.G.; Ahlawat, S.P.; Tyagi, V.; Semwal, D.P.; Gautam, N.K.; Rana, J.C.; Kumar, A. Cowpea genetic resources and its utilization: Indian perspective-A review. Legume Res. 2019, 42, 437–446. [Google Scholar] [CrossRef]
  24. Moses, D.R.; Zibokere, D.S. Effect of moisture content on some physical and mechanical properties of three varieties of cowpea (Vigna Unguiculata (L.) Walp). Agric. Eng. Int. CIGR J. 2011, 13, 1. [Google Scholar]
  25. de Oliveira Moura, J.; Rocha, M.D.M.; Gomes, R.L.F.; Freire, F.R.F.; Silva, K.J.D.E.; Ribeiro, V.Q. Path analysis of iron and zinc contents and others traits in cowpea. Crop Breed. Appl. Biotechnol. 2012, 12, 245–252. [Google Scholar] [CrossRef]
  26. Nalawade, A.D.; Patil, S.M.; Rajwade, P.R.; Kauthale, V.K. Evaluation of cowpea germplasm by using agro-morphological characters. Indian J. Agric. Res. 2020, 55, 364–368. [Google Scholar] [CrossRef]
  27. Rife, T.; Poland, J. Field Book: An open-source application for field data collection on android. Crop Sci. 2014, 54, 1624–1627. [Google Scholar] [CrossRef]
  28. IPGRI. Cowpea Descriptors; International Board for Plant Genetic Resources, FAO: Rome, Italy, 1983. [Google Scholar]
  29. Akinyele, I.O.; Onigbinde, A.O.; Hussain, M.A.; Omololu, A. Physicochemical characteristics of 18 cultivars of Nigerian cowpeas (V. uniguiculafa) and their cooking properties. J. Food Sci. 1986, 51, 1483–1485. [Google Scholar] [CrossRef]
  30. Kinyanjui, P.; Njoroge, D.; Makokha, A.; Christiaens, S.; Sila, D.; Hendrickxm, M. Quantifying the effects of postharvest storage and soaking pretreatments on the cooking quality of common beans (Phaseolus vulgaris). J. Food Process. Preserv. 2017, 41, e13036. [Google Scholar] [CrossRef]
  31. Vindiola, O.L.; Seib, O.A.; Hoseney, R.C. Accelerated development of hard-to-cook state in beans. Cereal Foods World 1986, 31, 538–552. [Google Scholar]
  32. Munthali, J.; Nkhata, S.G.; Masamba, K.; Mguntha, T.; Fungo, R.; Chirwa, R. Soaking beans for 12 h reduces split percent and cooking time regardless of type of water used for cooking. Heliyon 2022, 8, e10561. [Google Scholar] [CrossRef]
  33. Gama, A.P.; Hung, Y.-C.; Adhikari, K. Optimization of Emulsifier and Stabilizer Concentrations in a Model Peanut-Based Beverage System: A Mixture Design Approach. Foods 2019, 8, 116. [Google Scholar] [CrossRef]
  34. Sahrawat, K.; Kumar, G.; Murthy, K. Sulfuric acid-selenium digestion for multi-element analysis in a single digest. Commun. Soil Sci. Plant Anal. 2002, 33, 3757–3765. [Google Scholar] [CrossRef]
  35. Okalebo, J.R.; Gathua, K.W. Laboratory methods of soil and plant analysis. A Work. Man. 2002, 2, 29–68. [Google Scholar]
  36. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021; Available online: https://www.R-project.org/ (accessed on 2 January 2023).
  37. Mendiburu, F. Agricolae: Statistical Procedures for Agricultural Research. In R Package Version; R Foundation for Statistical Computing: Vienna, Austria, 2010; Volume 1, pp. 1–8. [Google Scholar]
  38. Lê, S.; Josse, J.; Husson, F. FactoMineR: An R package for multivariate analysis. J. Stat. Softw. 2008, 25, 1–18. [Google Scholar] [CrossRef]
  39. Maechler, M.; Rousseeuw, P.; Struyf, A.; Hubert, M.; Hornik, K. Cluster: Cluster Analysis Basics and Extensions; R Package Version 2.1.6; R Foundation for Statistical Computing: Vienna, Austria, 2023. [Google Scholar]
  40. Boukar, O.; Muranaka, S.; Franco, J.; Fatokun, C.C. Protein and mineral composition in grains of elite cowpea lines. In Proceedings of the Fifth World Cowpea Conference—Improving Livelihoods in the Cowpea Value Chain through Advancement in Science, 26 September–1 October 2010, Saly, Senegal; Boukar, O., Coulibaly, O., Fatokun, C.A., Lopez, K., Tamo, M., Eds.; IITA: Ibadan, Nigeria, 2012; pp. 88–100. [Google Scholar]
  41. Muranaka, S.; Shono, M.; Myoda, T.; Takeuchi, J.; Franco, J.; Nakazawa, Y.; Boukar, O.; Takagi, H. Genetic diversity of physical, nutritional and functional properties of cowpea grain and relationships among the traits. Plant Genet. Resour. 2015, 1, 1–10. [Google Scholar] [CrossRef]
  42. Ishikawa, H.; Drabo, I.; Joseph, B.B.; Muranaka, S.; Fatokun, C.; Boukar, O. Characteristics of farmers’ selection criteria for cowpea (Vigna unguiculata) varieties differ between north and south regions of Burkina Faso. Exp. Agric. 2020, 56, 94–103. [Google Scholar] [CrossRef]
  43. Quaye, W.; Adofo, K.; Buckman, S.E.; Frempong, G.; Jongerden, J.; Ruivenkamp, G. A socio-economic assessment of cowpea diversity on the Ghanaian market: Implications for breeding. Int. J. Consum. Stud. 2011, 35, 679–687. [Google Scholar] [CrossRef]
  44. Kim, D.K.; Son, D.M.; Choi, J.K.; Shin, H.R.; Choi, K.J.; Lee, J.R.; Lee, K.D.; Rim, Y.S. Agronomic characteristics and seed quality of cowpea (Vigna unguiculata L.) Germplasm. Korean J. Crop Sci. 2013, 58, 1–7. [Google Scholar] [CrossRef]
  45. Olapade, A.A.; Okafor, G.I.; Ozumba, A.U.; Olatunji, O. Characterization of common Nigerian cowpea (Vigna unguiculata L. Walp) varieties. J. Food Eng. 2002, 55, 101–105. [Google Scholar] [CrossRef]
  46. Nielsen, S.S.; Brandt, W.E.; Singh, B.B. Genetic variability for nutritional composition and cooking time of improved cowpea lines. Crop Sci. 1993, 33, 469–472. [Google Scholar] [CrossRef]
  47. Jean Baptiste, N.T.; Bell, J.M.; Nassourou, A.M.; Njintang, N.Y.; Youmbi, E. Genetic analysis of seed proteins contents in cowpea (Vigna unguiculata L. Walp.). Afr. J. Biotechnol. 2011, 10, 3077–3086. [Google Scholar] [CrossRef]
  48. Sanni, S.A.; Adebowale, A.-R.; Yusuf, A. Physical, proximate, and cooking properties of selected Nigerian cowpea varieties. Niger. Food J. 2006, 24, 34–41. [Google Scholar] [CrossRef]
Figure 1. Monthly rainfall data at Bunda from November 2021 to April 2022.
Figure 1. Monthly rainfall data at Bunda from November 2021 to April 2022.
Agriculture 14 00633 g001
Figure 2. Frequency distribution of traits: (a) Protein content (%), (b) iron content (mg/kg), (c) zinc content (mg/kg), (d) cooking time (min), (e) broth viscosity (cP), (f) seed length (mm), (g) seed width (mm), (h) 100 seed weight (g).
Figure 2. Frequency distribution of traits: (a) Protein content (%), (b) iron content (mg/kg), (c) zinc content (mg/kg), (d) cooking time (min), (e) broth viscosity (cP), (f) seed length (mm), (g) seed width (mm), (h) 100 seed weight (g).
Agriculture 14 00633 g002
Figure 3. Genotype by trait biplot. A total of 306 cowpea genotypes measured for eight quality traits: Pr = Protein content, Fe = Iron content, Zn = Zinc content, CT = Cooking time, BV = Broth Viscosity, SDLEN = Seed length, SDWDT = Seed width, HSDWT = 100 seed weight, black dots represent the different genotypes evaluated.
Figure 3. Genotype by trait biplot. A total of 306 cowpea genotypes measured for eight quality traits: Pr = Protein content, Fe = Iron content, Zn = Zinc content, CT = Cooking time, BV = Broth Viscosity, SDLEN = Seed length, SDWDT = Seed width, HSDWT = 100 seed weight, black dots represent the different genotypes evaluated.
Agriculture 14 00633 g003
Table 1. Mean square values for seed length, seed width, and one hundred seed weights.
Table 1. Mean square values for seed length, seed width, and one hundred seed weights.
Source of VariationDFSeed LengthSeed WidthHSDWT
Rep2124.68 ***9.49 ***543.50 ***
Block:Rep994.67 ***2.14 ***10.90 **
Genotype30518.55 ***6.31 ***38.30 ***
Residuals41831.180.547.00
DF = Degrees of freedom, HSDWT = 100 seed weight, *** highly significant at p < 0.001, ** highly significant at p < 0.01.
Table 2. Mean performance of the top ten and bottom ten genotypes was ranked based on 100 seed weights.
Table 2. Mean performance of the top ten and bottom ten genotypes was ranked based on 100 seed weights.
GenotypeSeed Length
(mm)
Seed Width
(mm)
HSDWT
(g)
Seed Coat ColorSeed Shape
Top 10 genotypes
MWcp369.33 ± 1.296.53 ± 1.0626.13 ± 0.67BrownRhomboid
MWcp0310.33 ± 0.986.67 ± 1.0523.70 ± 1.56 PurpleOvoid
MWcp677.80 ± 0.945.53 ± 0.5222.57 ± 3.32 Reddish-brownRhomboid
MWcp608.13 ± 1.136.13 ± 0.7422.23 ± 1.53 Reddish-brownRhomboid
UAM14-126-19-28.60 ± 1.245.80 ± 0.6821.70 ± 7.01 WhiteRhomboid
MWcp458.67 ± 1.356.00 ± 1.0721.30 ± 0.70 Reddish-brownRhomboid
MWcp248.87 ± 0.746.20 ± 0.8621.03 ± 1.24 Reddish-brownRhomboid
MWcp538.73 ± 1.285.93 ± 0.5921.00 ± 4.75 Reddish-brownOvoid
MWcp048.00 ± 1.775.27 ± 0.7020.97 ± 2.78 Reddish-brownRhomboid
TVu-140049.80 ± 0.775.73 ± 0.7020.93 ± 4.16 WhiteKidney
Bottom 10 genotypes
TVu-35265.40 ± 0.744.07 ± 0.707.87 ± 0.50 PurpleOvoid
TVu-11775.67 ± 0.623.73 ± 0.467.77 ± 0.40 BrownOvoid
TVu-30636.87 ± 1.064.33 ± 0.907.30 ± 0.79 WhiteOvoid
TVu-815.4 ± 0.633.73 ± 0.467.23 ± 0.95 WhiteRhomboid
TVu-32286.27 ± 1.104.73 ± 0.707.00 ± 0.40 Purplish-brownOvoid
TVu-32295.20 ± 0.684.07 ± 0.467.00 ± 1.11 PurpleOvoid
TVu-9725.80 ± 0.943.87 ± 0.646.97 ± 0.32 BrownOvoid
TVu-170605.07 ± 0.803.87 ± 0.646.53 ± 0.65 PurpleOvoid
TVu-32176.13 ± 0.834.27 ± 0.596.13 ± 0.70 BrownRhomboid
TVu-101696.27 ± 1.283.8 ± 0.566.00 ± 1.42 BlackRhomboid
Grand mean7.174.9312.86
Pr (>F)0.000.000.00
CV (%)15.1414.9720.51
CV = Coefficient of variation, Values with the same within the column.
Table 3. Mean square values for cooking quality traits in the study.
Table 3. Mean square values for cooking quality traits in the study.
Source of VariationDFBroth ViscosityCooking Time
Genotype3031.17 ***1274.80 ***
Residuals3040.4510.70
DF = Degrees of freedom, *** highly significant at p < 0.001.
Table 4. Mean values for cooking quality traits (top ten and bottom ten genotypes were ranked based on cooking time).
Table 4. Mean values for cooking quality traits (top ten and bottom ten genotypes were ranked based on cooking time).
EntryCooking Time (min)Broth Viscosity (cP)SDCCOLSDSHP
Top 10 genotypes
TVu-1563139.00 ± 1.414.60 ± 0.28RedKidney
TVu-32845.50 ± 2.12 4.20 ± 0.14MottledOvoid
IT98K-131-246.50 ± 4.95 3.75 ± 1.34BrownRhomboid
TVu-2247.50 ± 3.54 3.90 ± 0.14MottledKidney
TVu-355048.50 ± 0.71 5.00 ± 0.00Reddish-brownKidney
IT00k-126-358.50 ± 2.12 4.60 ± 0.28BrownOvoid
TVu-1400449.50 ± 2.12 3.60 ± 1.41WhiteKidney
MWcp4049.50 ± 2.12 3.85 ± 1.34PurpleOvoid
TVu-1326549.50 ± 2.12 3.85 ± 1.34WhiteKidney
TVu-266149.50 ± 3.54 2.75 ± 0.64BrownOvoid
Bottom 10 genotypes
TVu-11674137.50 ± 0.71 2.95 ± 0.64BrownRhomboid
MWcp31145.00 ± 2.83 3.20 ± 0.28BrownOvoid
MWcp17159.50 ± 0.71 4.20 ± 0.14WhiteRhomboid
MWcp61160.50 ± 3.54 4.60 ± 0.28MottledOvoid
MWcp64167.50 ± 0.71 3.30 ± 0.28BrownOvoid
MWcp54171.00 ± 5.66 4.10 ± 0.14MottledRhomboid
MWcp29171.50 ± 6.36 4.00 ± 0.14PurpleOvoid
MWcp46180.50 ± 0.71 2.40 ± 0.0MottledRhomboid
MWcp43181.00 ± 10.41 3.65 ± 0.21WhiteRhomboid
MWcp07191.00 ± 7.07 2.65 ± 0.64MottledOvoid
Grand mean81.523.62
Pr (>F)0.000.00
CV (%)4.0118.48
CV = Coefficient of variation, SDCCOL = Seed coatcolor, SDSHP = Seed shape.
Table 5. Mean square values for nutritional quality traits in the study.
Table 5. Mean square values for nutritional quality traits in the study.
Source of VariationDFProteinIronZinc
Genotype30534.42 ***13.58 ***0.10 ***
Residuals3061.050.050.00
DF = Degrees of freedom, *** highly significant at p < 0.001.
Table 6. Mean values for nutritional quality traits (top ten and bottom ten genotypes were ranked based on protein content).
Table 6. Mean values for nutritional quality traits (top ten and bottom ten genotypes were ranked based on protein content).
EntryProtein (%)Iron (mg/kg)Zinc (mg/kg)SDCCOLSDSHP
Top 10 genotypes
TVu-324322.92 ± 0.25 8.39 ± 0.000.45 ± 0.00BrownOvoid
MWcp30522.21 ± 0.85 6.56 ± 0.000.34 ± 0.00BrownOvoid
TVu-270621.32 ± 0.10 3.20 ± 0.000.39 ± 0.00Reddish-brownKidney
TVu-309421.29 ± 0.15 1.29 ± 0.000.36 ± 0.00WhiteKidney
MWcp6921.21 ± 0.75 8.67 ± 0.000.39 ± 0.00PurpleOvoid
MZcp 02420.97 ± 0.60 5.00 ± 0.000.34 ± 0.00BrownOvoid
TVu-352420.40 ± 0.10 1.19 ± 0.000.29 ± 0.00BrownRhomboid
IT 99K-529-120.29 ± 0.15 7.15 ± 0.000.45 ± 0.00BrownOvoid
TVu-1016920.18 ± 0.70 1.93 ± 0.000.33 ± 0.00BlackRhomboid
TVu-148320.04 ± 0.20 1.64 ± 0.000.43 ± 0.00Purplish- brownOvoid
Bottom 10 genotypes
TVu-10154.29 ± 0.10 1.58 ± 0.000.27 ± 0.00BrownOvoid
MWcp6014.19 ± 0.03 1.89 ± 0.000.22 ± 0.00BrownOvoid
MWcp484.17 ± 0.15 1.90 ± 0.001.77 ± 0.02PurpleOvoid
MWcp464.09 ± 0.03 2.00 ± 0.000.23 ± 0.00MottledRhomboid
IT90K-764.06 ± 0.10 1.57 ± 0.010.31 ± 0.00BrownOvoid
TVx-32364.00 ± 0.03 2.67 ± 0.000.27 ± 0.00WhiteRhomboid
TZcp673.78 ± 0.20 1.94 ± 0.000.18 ± 0.00BrownOvoid
MWcp503.64 ± 0.10 1.67 ± 0.000.24 ± 0.00PurpleOvoid
MZcp 0043.63 ± 0.07 1.98 ± 0.000.17 ± 0.00BrownOvoid
TVu-32633.60 ± 0.15 1.24 ± 0.000.28 ± 0.00PurpleRhomboid
Grand mean11.053.210.36
Pr (>F)0.000.000.00
CV%9.277.101.23
CV = Coefficient of variation, SDCCOL = Seed coat color, SDSHP = Seed shape.
Table 7. The correlation matrix among grain, cooking, and nutritional quality traits of cowpea genotypes.
Table 7. The correlation matrix among grain, cooking, and nutritional quality traits of cowpea genotypes.
Trait 1Trait 2r95% CItdfp
ProteinIron0.28[0.17, 0.38]5.05304<0.001 ***
ProteinZinc0.21[0.10, 0.31]3.723040.005 **
ProteinBroth−0.07[−0.18, 0.04]−1.23302>0.999
ProteinCT−0.03[−0.14, 0.08]−0.53302>0.999
ProteinSeed.L−0.06[−0.17, 0.05]−1.09304>0.999
ProteinSeed.W−0.03[−0.15, 0.08]−0.59304>0.999
ProteinHSDWT−0.06[−0.17, 0.05]−1.08304>0.999
IronZinc0.24[0.13, 0.34]4.34304<0.001 ***
IronBroth−0.04[−0.15, 0.08]−0.61302>0.999
IronCT−0.08[−0.19, 0.03]−1.45302>0.999
IronSeed.L−7.24 × 10−03[−0.12, 0.10]−0.13304>0.999
IronSeed.W0.03[−0.09, 0.14]0.47304>0.999
IronHSDWT−0.06[−0.17, 0.05]−1.08304>0.999
ZincBroth0.03[−0.08, 0.14]0.5302>0.999
ZincCT−0.04[−0.15, 0.07]−0.75302>0.999
ZincSeed.L0.07[−0.04, 0.18]1.29304>0.999
ZincSeed.W0.06[−0.05, 0.17]1.06304>0.999
ZincHSDWT0.03[−0.09, 0.14]0.48304>0.999
BrothCT−0.04[−0.16, 0.07]−0.76302>0.999
BrothSeed.L6.74 × 10−03[−0.11, 0.12]0.12302>0.999
BrothSeed.W−0.03[−0.14, 0.08]−0.54302>0.999
BrothHSDWT−0.09[−0.20, 0.02]−1.61302>0.999
CTSeed.L0.21[0.10, 0.31]3.723020.005 **
CTSeed.W0.26[0.15, 0.36]4.67302<0.001 ***
CTHSDWT0.19[0.08, 0.30]3.443020.013 *
Seed.LSeed.W0.89[0.86, 0.91]33.49304<0.001 ***
Seed.LHSDWT0.71[0.65, 0.76]17.57304<0.001 ***
Seed.WHSDWT0.75[0.70, 0.80]19.88304<0.001 ***
r = Correlation coefficient, CI = Confidence interval, t = t-value, df = Degrees of freedom, p = p-value, CT = Cooking time, Seed.L = Seed length, Seed.W = Seed width, HSDWT = 100 seed weight, Broth = Broth viscosity. *, **, ***=Significant at 5%, 1% and 0.1%, respectively. (p-value adjustment method: Holm (1979)).
Table 8. Trait contributions, eigenvalues, and cumulative percentage of the components.
Table 8. Trait contributions, eigenvalues, and cumulative percentage of the components.
Traits Principal ComponentsTrait Contributions (%)
PC1PC2PC3PC1PC2PC3
Seed length0.570.040.1132.130.161.30
Seed width0.580.060.0533.590.340.26
100 seed weight0.53−0.00−0.0228.510.000.05
Cooking time0.22−0.13−0.295.081.628.65
Broth viscosity−0.04−0.080.890.160.6579.84
Protein−0.060.56−0.200.3131.923.84
Iron−0.030.600.010.0836.040.01
Zinc0.040.540.250.1436.046.05
Eigenvalue2.671.511.04
% of variance33.3618.8112.98
Cumulative % of variance33.3752.1965.17
Table 9. Descriptive statistics of the clusters.
Table 9. Descriptive statistics of the clusters.
ClusterNStatisticPr (%)Fe (mg/kg)Zn (mg/kg)CT (min)BV (cP)Seed.L (mm)Seed.W (mm)HSDWT (g)
I234M10.953.260.3576.983.686.744.6611.43
Min3.600.320.0839.001.603.732.936.00
Max22.9221.441.23145.005.508.935.7318.80
II72M11.393.060.4096.033.438.575.7917.58
Min3.640.120.0949.502.205.674.5312.80
Max21.2114.642.63191.004.9010.337.1326.13
Pr = Protein, Fe = Iron, Zn = Zinc, CT = Cooking time, BV = Broth Viscosity, Seed.L = Seed length, Seed.W = Seed width, HSDWT = 100 seed weight.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chipeta, M.M.; Yohane, E.; Kafwambira, J.; Kampanje-Phiri, J. Genetic Variability and Interrelationships of Grain, Cooking, and Nutritional Quality Traits in Cowpea: Implications for Cowpea Improvement. Agriculture 2024, 14, 633. https://doi.org/10.3390/agriculture14040633

AMA Style

Chipeta MM, Yohane E, Kafwambira J, Kampanje-Phiri J. Genetic Variability and Interrelationships of Grain, Cooking, and Nutritional Quality Traits in Cowpea: Implications for Cowpea Improvement. Agriculture. 2024; 14(4):633. https://doi.org/10.3390/agriculture14040633

Chicago/Turabian Style

Chipeta, Michael M., Esnart Yohane, John Kafwambira, and Jessica Kampanje-Phiri. 2024. "Genetic Variability and Interrelationships of Grain, Cooking, and Nutritional Quality Traits in Cowpea: Implications for Cowpea Improvement" Agriculture 14, no. 4: 633. https://doi.org/10.3390/agriculture14040633

APA Style

Chipeta, M. M., Yohane, E., Kafwambira, J., & Kampanje-Phiri, J. (2024). Genetic Variability and Interrelationships of Grain, Cooking, and Nutritional Quality Traits in Cowpea: Implications for Cowpea Improvement. Agriculture, 14(4), 633. https://doi.org/10.3390/agriculture14040633

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop