ABSTRACT
Seed color, brightness and texture are important aspects for the consumption of immature beans. This study aimed to estimate the phenotypic variability among thirty cowpea genotypes for immature seed harvesting, in order to guide breeding programs. The experiment was carried out in a randomized blocks design, with three replications. The evaluated traits were immature seeds yield, pods and dry seeds production, number of immature pods, 100-immature seeds weight, immature pods length, number of seeds per immature pod, cooking time, harvesting period of immature pods, seed color, pod color, plant growth habit, opening of immature pods, seed threshing from immature pods, cultivation value and seed browning after threshing. Anova and three clustering methods were applied: UPGMA using the Jaccard coefficient of similarity; modified Tocher method based on the Gower algorithm; and k-means. All the analyses were performed based on the variable means of experiments evaluated in six environments. A statistically significant phenotypic variability was observed among all the quantitative traits. Discrepancies occurred in the accessions grouping and number of clusters. The accession BRS Acauã, in one group, and P508 and PCCR3F6L15, in other groups, showed a good quality for immature grain harvesting, were consistently grouped into different clusters and are recommended for use as parents in the development of new cultivars.
KEYWORDS: Vigna unguiculata ; cluster analysis; genetic variability
RESUMO
Cor, brilho e textura das sementes são importantes para o consumo de feijão imaturo. Objetivou-se estimar a variabilidade fenotípica entre trinta genótipos de feijão-caupi destinados à colheita de sementes imaturas, para orientar programas de melhoramento. O experimento foi conduzido em delineamento de blocos casualizados, com três repetições. As características avaliadas foram rendimento de sementes imaturas, produção de vagens e sementes secas, número de vagens imaturas, peso de 100 sementes imaturas, comprimento de vagens imaturas, número de sementes por vagem imatura, tempo de cozimento, período de colheita de vagens imaturas, cor da semente, cor da vagem, hábito de crescimento da planta, abertura de vagens imaturas, debulha de sementes de vagens imaturas, valor de cultivo e escurecimento das sementes após a debulha. A Anova e três análises de agrupamento foram adotadas: UPGMA, com o coeficiente de similaridade de Jaccard; método Tocher modificado, com o algoritmo de Gower; e k-means. Todas as análises foram realizadas com base em médias das variáveis de experimentos avaliados em seis ambientes. Variabilidade fenotípica estatisticamente significativa foi observada entre todas as variáveis quantitativas. Discrepâncias ocorreram no agrupamento de acessos e número de agrupamentos. O acesso BRS Acauã, em um grupo, e P508 e PCCR3F6L15, em outros grupos, apresentaram boa qualidade para colheita imatura, foram consistentemente agrupados em diferentes grupos e são indicados para utilização como genitores no desenvolvimento de novas cultivares.
PALAVRAS-CHAVE: Vigna unguiculata ; análise de agrupamento; variabilidade genética
INTRODUCTION
In northeastern Brazil, cowpea [Vigna unguiculata (L.) Walp] is mainly grown for dry seed production, although the commercialization of immature seeds is also expressive. Immature seeds are obtained from pods around maturity, i.e., immediately before or shortly after the plant stops accumulating photosynthates and begins the natural dehydration process, which occurs with seed moisture close to 60 or 70 % (Sousa et al. 2015).
Cowpea has a high genetic variability for almost all morphological and agronomic traits, though little has been investigated in this regard (Teixeira et al. 2010). However, information on the phenotypic variability of cowpea germplasm and on how much of it is due to genetic differences is fundamental in breeding programs (Cruz et al. 2004). According to Cruz et al. (2004), diversity denotes the genetic distance within species, populations or accessions measured through morphological, agronomic, biochemical, physiological and molecular traits.
In the Rio Grande do Norte state (Brazil), Silva Filho et al. (2013) evaluated the morphological traits of dry cowpea seeds and observed variability, highlighting the selection potential of several accessions that can be incorporated into breeding programs. Moreover, Andrade (2010) studied immature seed harvesting in 20 cowpea accessions and found genetic variability for number of days to flowering, number of days to immature pod maturity, immature pod length, cooking time and other traits related to immature seeds. However, an extensive characterization of cowpea germplasm for immature seed harvesting has not been conducted to date in the São Francisco Valley region, to assist in developing cowpea genotypes suitable for cultivation. Therefore, this study aimed to estimate the phenotypic variability of cowpea genotypes for immature seed harvesting and select accessions with traits of agronomic interest for breeding programs in the semi-arid region of Brazil.
MATERIAL AND METHODS
Thirty cowpea genotypes were selected to evaluate the crop phenotypic variability, including landraces, cultivars and some breeding lines with potential for immature seed harvesting. These genotypes consisted of 14 lines provided by the Embrapa Semiárido, six commercial cultivars and ten landraces collected in cowpea-producing areas of Juazeiro (Bahia state) and Petrolina (Pernambuco state), Brazil (Table 1). The experiments were conducted in experimental fields located in Bebedouro (Petrolina) and Mandacaru (Juazeiro), during the second half of 2013, 2014 and 2015.
A randomized blocks design, with three replications, was used. The plots measured 6 m2 and had two rows spaced 1.0 m and plants spaced 0.1 m, representing a population of 100,000 plants ha-1. As a common practice in the region, soil fertilization did not occur. Moreover, a micro-sprinkler system was used for irrigation twice a week, to avoid soil waterlogging. Weeding and pest control with insecticides were performed whenever necessary.
The evaluated agronomic traits were immature seed yield (kg ha-1), dry seed yield (kg ha-1), immature pod yield (kg ha-1), number of immature pods, 100-immature seed weight (g), immature pod length (cm), number of seeds per immature pod, cooking time (Matteson® cooker; min) and harvesting period of immature pods (days).
The evaluated morphological traits were seed and pod color (1 = cream; 2 = green; 3 = purple); plant growth habit (1 = erect; 2 = semi-erect; 3 = semi-prostrate; 4 = prostrate); opening of immature pods by manual threshing and seed release from immature pods by manual threshing (1 = very difficult; 2 = difficult; 3 = easy; 4 = very easy); cultivation value considering the appropriate traits for commercial cultivation (1 = none - inappropriate; 2 = little - poor; 3 = many - adequate; 4 = most - good; 5 = all - very good; 6 = all - excellent); and seed browning at threshing and 24 and 48 h after threshing (1 = light; 2 = little dark; 3 = medium dark; 4 = dark; 5 = very dark).
The analysis of variance (Anova) of the six experiments considered only the agronomic traits, and all means were grouped by the Scott-Knott test at 5 % of probability, using the Genes software (Cruz 2016). The means of immature seed yield, dry seed yield and immature pod yield for each environment were adjusted to the mean number of plants in the experimental units using the covariance method with the aid of a SAS code (Vencovsky & Barriga 1992).
Origin of the cowpea genotypes [Vigna unguiculata (L.) Walp] evaluated for immature seed harvesting [Pernambuco (PE) and Bahia (BA) states, Brazil].
The means of the cowpea genotypes were subjected to three cluster analyses: unweighted pair group method with arithmetic mean (UPGMA), using the NTSYs software (Rohlf 2000), considering the Jaccard similarity coefficient with nine quantitative variables categorized into five classes; presence (1) or absence (0) of traits evaluated in each class using the modified Tocher method based on the Gower algorithm with the Genes software (Cruz 2016); k-means clustering with qualitative data transformed to √x + 1, considering the formation of three groups (cultivars, lines and varieties from local producers). The k-means analysis was performed using the fastclus procedure, which employs the nearest centroid sorting method, and the Euclidean distance to evaluate the dissimilarity of the means (SAS Institute 1989).
RESULTS AND DISCUSSION
Seed and pod color were different, with a predominant cream color for seed color. The pod color of the genotypes BRS-Guariba, PC951016D01E, CPCR3F6L15, C1J, C2S, C3S and P303 was purple. In contrast, the genotypes PC951015D01E, P508, PJM22 and PL23 had a cream pod color, while the remaining genotypes were green (Table 2). The landraces had semi-prostrate and prostrate plant growth habits, while the cowpea lines and cultivars from research institutes had semi-erect and prostrate plant growth habits (Table 2).
The genotypes differed with regard to the opening of immature pods and seed release from immature pods by manual threshing. The genotypes BRS-Acauã, C3Q, PJJ21 and PD26 showed an easier opening of immature pods and seed release from immature pods (Table 2). The cultivation values of all genotypes were close to the maximum value required for commercial cultivation, except for the genotype C3F, which showed plants with many appropriate traits (Table 2). The seed browning after threshing showed differences among all the genotypes, with the BRS-Acauã and PJN29 showing the lowest seed browning at 48 h after threshing (Table 2).
The cowpea genotypes showed statistically significant differences by the F-test (p < 0.05), suggesting a high phenotypic variability for all the agronomic variables. The immature seed yield of the cowpea genotypes ranged from 931 (BRS-Marataoã) to 2,532 kg ha-1 (P508) (Table 3). On the other hand, the BRS-Guariba showed the highest immature seed yield (2,029 kg ha-1) among the cultivars, whereas, among the cowpea lines, P508, PCCR3F6L15 and P303 showed the highest immature seed yield (Table 3). The immature pod yield ranged from 1,783 (genotype PD26) to 4,715 kg ha-1 (genotype P508) (Table 3).
Means of nine morphological traits of thirty cowpea genotypes [Vigna unguiculata (L.) Walp] evaluated for immature seed harvesting.
The immature pod yield values of the genotypes PC951015D01E, PC950409D02E, PCCR3FL15 and P508 were higher than 3,700 kg ha-1. Furthermore, the dry seed yield ranged from 557 (genotype PLP30) to 1,350 kg ha-1 (genotype P508). The genotypes BRS-Patativa, PCCR3F6L15 and PJM22 showed dry seed yield values higher than 1,200 kg ha-1. On the other hand, the cooking time in a Matteson® cooker ranged from 10 to 15.3 min for the genotypes PLP30, BRS-Marataoã and C3Q (Table 3). In addition, the cowpea genotypes PD26, BRS-Pujante, BRS-Rouxinol, C2M, PJM22 and PL23 showed immature pod length values higher than 20 cm. On the other hand, the harvesting period of immature pods ranged from 15.6 to 22.6 days, and the genotypes BRS-Acauã, PCCR3F6F15, C2M and C3F showed the longest harvesting period of immature pods (Table 3).
The cowpea cultivars, lines and landraces formed a large group of 24 genotypes by the k-means clustering (Table 4). However, another large group consisting of 19 genotypes was formed when the Tocher method was used (Table 5). All 19 genotypes of group 1, formed using the Tocher method, were clustered in group 2 using the k-means method. However, the cowpea lines C3F and C3S and the genotype PD26 were clustered in the same group by the two methods, while the other cowpea genotypes showed inconsistencies in clusters developed by both methods (Tables 4 and 5).
Mean square of treatments (MST), coefficient of variation (CV) and means of nine agronomic variables of thirty cowpea genotypes [Vigna unguiculata (L.) Walp] evaluated for immature seed harvesting.
The cophenetic value by the UPGMA method was 0.69, indicating limitations in the developed dendrogram. The coefficient of similarity in this dendrogram ranged from 18 to 54 %, suggesting a high variability among the cowpea genotypes. The cut-off point at the similarity distance of 0.21 formed seven groups. The genotypes were clustered as group 1 (BRS-Acauã, BRS-Guariba, BRS-Pujante and CPCR3F6L17), group 2 (C3Q), group 3 (C3F, PJM22, C3S, PD26 and PLP30), group 4 (BRS-Marataoã, P508, PCCR3F6L15 and C2M), group 5 (the remaining 14 genotypes), group 6 (PJJ21) and group 7 (PC25). The cultivar BRS-Patativa and the line C2S showed the highest similarity (54 %), and the genotype PC25 showed the greatest divergence among the evaluated genotypes, justifying its location at the base of the dendrogram (Figure 1).
Due to its genetic variability, high adaptability, yield potential, excellent nutritional value and good behavior, in relation to water deficit, cowpea is considered a valuable crop for semi-arid regions. However, although highly adaptable to the most diverse environments, Santos (2011) reported that cowpea still reaches low dry seed yields, estimated at around 300 kg ha-1 in northeastern Brazil. According to Silva et al. (2010), climatic adversities, poor agricultural technology and the use of unimproved seeds are the leading causes of low yields in the semi-arid region of Brazil. Moreover, the growing conditions of different environments may result in unequal genotype performances. In order to estimate the phenotypic variability, the comprehensive evaluation of cowpea genotypes intended for immature seed harvesting is a crucial step toward developing new cowpea cultivars by crossing divergent genotypes with good agronomic performance, as proposed in this study.
Number and identification of thirty cowpea genotypes [Vigna unguiculata (L.) Walp] evaluated for eight qualitative and nine agronomic variables and distributed into three groups according to the k-means clustering method.
Number and identification of thirty cowpea genotypes [Vigna unguiculata (L.) Walp] evaluated for eight qualitative and nine agronomic variables and distributed into six groups according to the Tocher clustering method based on the Gower algorithm.
Cowpea production for immature seed or pod commercialization is often carried out by farmers who prefer genotypes with an indeterminate growth habit and uneven maturation, resulting in several harvests during the crop cycle and high yields. The consumer preference is predominantly for easily cooked, purple and immature pods and seeds with good taste and pleasant smell. Moreover, traders prefer genotypes that are easy to thresh and with a good post-harvest quality (Rocha et al. 2006).
The overall average of plant growth habit among the genotypes was 1.90, which was similar to the average found by Rocha et al. (2006) (1.93). The overall average of the cultivation value for the tested genotypes was 2.18, higher than the average found by Andrade et al. (2010) (1.89). The overall average of pod length for the tested genotypes was 20.1 cm, similar to the value reported by Andrade et al. (2010) (19.63 cm).
The overall average of 100-immature seed weight was 45.4, above the weight required by the immature seed market (30 g). This average was also higher than that obtained by Rocha et al. (2006) (36 g) and Andrade et al. (2010) (37 g). According to Rocha et al. (2007), genotypes for immature seed production include a longer flowering period, high pod yield, attractive and well-seeded pods, immature seed to pod weight ratio over 60 % and seeds capable of preserving good post-harvest and post-threshing aspects. In this scenario, fast darkening seeds are not well received by the market.
The number of groups differed according to the three applied methods, with the k-means clustering resulting in the smallest number of groups, since, in this method, there is a previously defined number of groups. The same discrepancies were observed for the accessions, except for C3F, C3S and PD26, which were consistently grouped by the k-means (group 1) and Tocher (group 2) methods. Another exception was observed for the accession PC25, consistently grouped by the Tocher method and the dendrogram based on the Jaccard coefficient.
Based on two clustering validation methods (Dunn and DB index), Adarsh et al. (2020) concluded that the non-hierarchical k-means clustering was superior to the Tocher method and other hierarchical methods adopted in the present study, e.g., UPGMA. The criticism against the k-means clustering is because the number of clusters is defined before the analysis (Cruz et al. 2004), which could occur based on some biological evidence, as adopted in the present study.
Unweighted pair group method with arithmetic mean dendrogram obtained from the evaluation of seventeen traits in thirty cowpea genotypes [Vigna unguiculata (L.) Walp]. Cophenetic value = 0.69.
Despite the preference for commercialization as immature seeds, there are no cowpea cultivars identified for this purpose, and studies on the evaluation of cowpea accessions for breeding programs are scarce. The thirty cowpea genotypes showed a high phenotypic variability for the agronomic traits evaluated, highlighting their potential for breeding programs involving immature seed harvesting in the near future.
The genotypes P508 and PCCR3F6L15 stood out with the highest immature seed yields, a harvesting period of about 20 days and cooking times of 12 and 13 min, respectively. They also had semi-erect to prostrate plant growth habits, easily opened pods and easy seed release by manual threshing. Nonetheless, these two genotypes showed dark seeds at 48 h after threshing, which is undesirable for the commercialization of immature seeds.
The genotype BRS-Acauã showed light-dark seeds at 48 h after threshing, with an immature seed yield of 2,000 kg ha-1 and a harvesting period of more than 22 days. Moreover, this genotype was clustered into different groups based on the cluster analysis conducted in this study. Therefore, BRS-Acauã is a good option for crosses with P508 and PCCR3F6L15 to develop new immature seed lines.
CONCLUSIONS
-
The thirty cowpea genotypes showed statistically significant differences (p < 0.05), suggesting a high phenotypic variability for all the agronomic variables evaluated for immature seed harvesting;
-
The number of groups differed according to the three methods applied, with the k-means clustering resulting in the smallest number of groups. The accession groups showed discrepancies, with only two grouping exceptions;
-
The accessions BRS Acauã, in one group, and P508 and PCCR3F6L15, in other groups, showed a good quality for immature grain harvest and were consistently grouped into different clusters by the three methods applied. These accessions are recommended as parents to develop new cultivars.
ACKNOWLEDGMENTS
The authors would like to thank the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (Capes), for the scholarship granted.
REFERENCES
- ADARSH, V. S.; JOSEPH, B.; GOPINATH, P. P. Comparison of clustering techniques used in divergence analysis. International Journal of Statistics and Applied Mathematics, v. 5, n. 4, p. 19-21, 2020.
- ANDRADE, F. N. Avaliação e seleção de linhagens de tegumento e cotilédone para o mercado de feijão-caupi verde 2010. Dissertação (Mestrado em Agronomia) - Universidade Federal do Piauí, Teresina, 2010.
- ANDRADE, F. N.; ROCHA, M. M.; GOMES, R. L. F.; FREIRE FILHO, F. R.; RAMOS, S. R. R. Estimativas de parâmetros genéticos em genótipos de feijão-caupi avaliados para feijão fresco. Agronomic Science Magazine, v. 41, n. 2, p. 253-258, 2010.
- CRUZ, C. D. Genes software: extended and integrated with the R, Matlab and Selegen. Version 2015.5.0. Acta Scientiarum, v. 38, n. 4, p. 547-552, 2016.
- CRUZ, C. D.; REGAZZI, A. J.; CARNEIRO, P. C. S. Modelos biométricos aplicados ao melhoramento genético 4. ed. Viçosa: UFV, 2004.
- ROCHA, M. M.; FREIRE FILHO, F. R.; RAMOS, S. R. R.; RIBEIRO, V. Q.; ANDRADE, F. N.; GOMES, R. L. F. Avaliação agronômica de genótipos de feijão-caupi para produção de grãos verdes Teresina: Embrapa Meio-Norte, 2006. (Boletim de pesquisa e desenvolvimento, 67).
- ROCHA, M. M.; FREIRE FILHO, F. R.; RIBEIRO, V. Q.; CARVALHO, H. W. L. C.; BELARMINO FILHO, J.; RAPOSO, J. A. A.; ALCÂNTARA, J. P.; RAMOS, S. R. R.; MACHADO, C. F. Adaptabilidade e estabilidade produtiva de genótipos de feijão-caupi de porte semiereto na região Nordeste do Brasil. Pesquisa Agropecuária Brasileira, v. 42, n. 9, p. 1283-1289, 2007.
- ROHLF, F. J. NTSYS-pc numerical taxonomy and multivariate analysis system Version 2.1. New York: Exeter, 2000.
- SANTOS, C. A. F. Cultivares de feijão-caupi para o Vale do São Francisco Petrolina: Embrapa Semiárido, 2011. (Circular técnica, 94).
- SAS INSTITUTE. SAS/STAT user’s guide Version 6. Cary: SAS Institute, 1989.
- SILVA FILHO, A. J. R.; ANTONIO, R. P.; SILVA, P. S. L.; SILVEIRA, L. M.; ALBURQUERQUE, L. B. Avaliação morfológica e agronômica de sementes de acessos de caupi coletados no Rio Grande do Norte. Agropecuária Científica no Semiárido, v. 9, n. 4, p. 102-106, 2013.
- SILVA, V. P. R.; CAMPOS, J. H. B. C.; SILVA, M. T.; AZEVEDO, P. V. Impact of global warming on cowpea bean cultivation in northeastern Brazil. Agricultural Water Management, v. 97, n. 11, p. 1760-1768, 2010.
- SOUSA, J. L. M.; ROCHA, M. de M.; SILVA, K. J. D. e; NEVES, A. C. das; SOUSA, R. R. de. Potencial de genótipos de feijão-caupi para o mercado de vagens e grãos verdes. Pesquisa Agropecuária Brasileira, v. 50, n. 5, p. 392-398, 2015.
- TEIXEIRA, I. R.; SILVA, G. C.; OLIVEIRA, J. P. R.; SILVA, A. G.; PELÁ, A. Desempenho agronômico e qualidade de semente de cultivares de feijão-caupi na região do Cerrado. Revista Ciência Agronômica, v. 41, n. 2, p. 300-307, 2010.
- VENCOVSKY, R.; BARRIGA, P. Genética biométrica no fitomelhoramento Ribeirao Preto: Sociedade Brasileira de Genética, 1992.
Publication Dates
-
Publication in this collection
14 Jan 2022 -
Date of issue
2021
History
-
Received
09 July 2021 -
Accepted
04 Nov 2021 -
Published
26 Nov 2021