Acessibilidade / Reportar erro

GREEN BEAN YIELD AND PATH ANALYSIS IN COWPEA LANDRACES1 1 Paper extracted from the masters thesis of the first author.

RENDIMENTO DE GRÃOS VERDES E ANÁLISE DE TRILHA EM VARIEDADES TRADICIONAIS DE FEIJÃO-CAUPI

ABSTRACT

Cowpea (Vigna unguiculata (L.) Walp.) is widely cultivated in northeast Brazil for production of both green grains and dry grains. Green grains are cowpea grains with a 60 to 70% moisture content that are very appreciated by the local people. The cultivation of cowpea in northeast Brazil is performed mostly using landraces and the low yields achieved in some areas of that region are generally attributed to the cultivation of less productive landraces. The objectives of this research were to identify the best landraces in terms of green bean yield and the characters that hold the most positive direct effects on this yield via path analysis. Twelve landraces of undetermined growth were evaluated in two experiments conducted in Mossoró, Brazil. The first experiment (E1) was sprinkler-irrigated and the second (E2) was conducted under rainfed conditions, but was irrigated whenever necessary. Pod yield and dimensions and green grains were evaluated, in addition to the main grain yield components. A randomized block design with five replications was used. The data presented here identified eight landraces with the highest yields. In addition, cowpea yield was higher in the dry season (E1) than in the rainy season (E2). The number of pods per plant had the largest direct positive effect on green grain yield.

Keywords:
Vigna unguiculata; Green grain; Landraces.

RESUMO

O feijão-caupi (Vigna unguiculata (L.) Walp.) é amplamente cultivado na região nordeste do Brasil para a produção de grãos secos e de feijão verde. O feijão verde, isto é, grãos de feijão-caupi com 60% a 70% de umidade, é muito apreciado pelos nordestinos. O cultivo do feijão-caupi no nordeste brasileiro é feito principalmente com variedades tradicionais. Presume-se que os baixos rendimentos da cultura em algumas áreas da referida região são devidos ao uso de variedades tradicionais pouco produtivas. Os objetivos do trabalho foram identificar as melhores variedades tradicionais quanto ao rendimento de grãos verdes e os caracteres que possuem os maiores efeitos positivos diretos sobre esse rendimento, via análise de trilha. Doze variedades tradicionais de crescimento indeterminado foram avaliadas em dois experimentos realizados em Mossoró-RN. O primeiro experimento (E1) foi realizado com irrigação por aspersão e o segundo, sob condições de sequeiro, mas irrigado quando necessário. A avaliação contemplou os rendimentos e dimensões de vagens e de grãos verdes e os componentes principais do rendimento de grãos. Utilizou-se o delineamento de blocos casualizados com cinco repetições. Oito variedades foram as mais produtivas. O rendimento das cultivares no experimento realizado fora da estação chuvosa (E1) foi superior ao do outro experimento. O número de vagens por planta foi a característica que apresentou o maior efeito positivo direto sobre o rendimento de grãos.

Palavras-chave:
Vigna unguiculata; Feijão-verde; Variedades crioulas.

INTRODUCTION

Cowpea (Vigna unguiculata (L.) Walp.) is a multipurpose crop used mainly for the production of green grains and dry grains. The cultivation of cowpea in northeast Brazil is performed mostly using landraces and the low cowpea yields achieved in some areas of the region are generally attributed to the cultivation of less productive landraces. In addition, there are few research institutions in northeast Brazil working on improving cowpea cultivars, and even fewer companies focused on producing improved seeds.

According to some authors (BRESEGHELLO; COELHO, 2013BRESEGHELLO, F.; COELHO, A. S. G. Traditional and modern plant breeding methods with examples in rice (Oriza sativa L.). Journal of Agricultural and Food Chemistry, Washington, v. 61, n. 35, p. 8277-8286, 2013.), landraces have lower yield, but others have shown that some landraces may be more productive than improved cultivars (MKOI; CHIMPHANGO; DAKORA, 2009; PEKSEN, 2004PEKSEN, A. Fresh pod yield and some pod characteristics of cowpea (Vigna unguiculata L. Walp.) genotypes from Turkey. Asian Journal of Plant Sciences, Deira, v. 3, n. 3, p. 269-273, 2004.). Thus, studies of traditional cultivars enable the identification of cultivars that can be used directly by farmers or used in breeding programs (BOUKAR; FATOKUN, 2009BOUKAR, O.; FATOKUN, C. Strategies in cowpea breeding. In: ZERIHUN, T. (Ed.). New approaches to plant breeding of orphan crops in Africa. Rome: FAO, 2009. cap. 6, p. 69-92.). A landrace of a seed-propagated crop is a variable population, which is identifiable and usually has a local name; it lacks ''formal'' crop improvement, is characterized by specific adaptations to the environmental conditions of the cultivation area, and is closely associated with the uses of the people who have developed and continue to grow it (POLEGRI; NEGRI, 2010POLEGRI, L.; NEGRI, V. Molecular markers for promoting agro-biodiversity conservation: a case study from Italy. How cowpea landraces were saved from extinction. Genetic Resources and Crop Evolution , Dordrecht, v. 57, n. 6, p. 867-880, 2010.).

Studies of cowpea landraces have been carried out in various regions of the world (GBAGUIDI et al., 2013GBAGUIDI, A. A. et al. Diversity and agronomic performances of the cowpea (Vigna unguiculata Walp.) landraces in Southern Benin. International Research Journal of Agricultural Science and Soil Science, Aliero, v. 3, n. 4, p. 121-133, 2013.; GHALMI et al., 2010GHALMI, N. et al. Morphological and molecular diversity within Algerian cowpea (Vigna unguiculata (L.) Walp.) landraces. Genetic Resources and Crop Evolution, Dordrecht, v. 57, n. 3, p. 371-386, 2010.) and have examined several characteristics other than grain yield. Knowledge of the association between characteristics of distinct landraces is important for breeding, especially when selection of a trait is hampered by low heritability and/or difficulties with measurement and identification (CRUZ et al., 2012CRUZ, C. D.; REGAZZI, A. J.; CARNEIRO, P. C. S. Modelos biométricos aplicados ao melhoramento genético. 4. ed. Viçosa, MG: Editora UFV, 2012. v. 1, 514p.). Moreover, path analysis has been shown to be the best method for comparisons between grain yield and landrace characteristics (OGUNBODEDE, 1989OGUNBODEDE, B. A. Comparison between three methods of determining the relationships between yield and eight of its components in cowpea, Vigna unguiculata L. Walp. Scientia Horticulturae, Amsterdam, v. 38, n. 3-4, p. 201-205, 1989.).

Path analysis was developed in 1918 by Sewall Wright and detailed later (WRIGHT, 1960WRIGHT, S. Path coefficients and path analysis: alternative or complementary. Biometrics, Washington, v. 16, n. 2, p. 189-202, 1960.). This analysis is based on the study of direct and indirect effects of characters on a relevant variable (overall yield), whose estimates are obtained by regression equations. Although correlation is an intrinsic characteristic of two characters under a given experimental condition, its decomposition is dependent on the set of studied characters (CRUZ et al., 2012CRUZ, C. D.; REGAZZI, A. J.; CARNEIRO, P. C. S. Modelos biométricos aplicados ao melhoramento genético. 4. ed. Viçosa, MG: Editora UFV, 2012. v. 1, 514p.).

Path analysis studies of green grain yield in cowpea showed that the number of pods per plant had the greatest direct effect on yield (HITIKSHA et al., 2014HITIKSHA, K. P.; ACHARYA, R. R.; SHEETAL, R. P. Interrelationship among green pod yield and its yield contributing characters and biochemical characters vegetable cowpea [Vigna unguiculata (L.) Walp.]. Trends in Biosciences, New Delhi, v. 7, n. 23, p. 3972-3976, 2014.; SAPARA; JAVIA, 2014SAPARA, G. K.; JAVIA, R. M. Correlation and path analysis in vegetable cowpea (Vigna unguiculata L.). International Journal of Plant Sciences, Amsterdam, v. 9, n. 1, p. 138-141, 2014.). In addition, these studies also confirmed that other characteristics such as pod length had important direct effects on yield (SANTOS et al. 2014SANTOS, A. et al. Correlations and path analysis of yield components in cowpea. Crop Breeding and Applied Biotechnology, Viçosa, v. 14, n. 2, p. 82-87, 2014., SAPARA; JAVIA, 2014).

In 2009, the 73 municipalities in the state of Rio Grande do Norte, Brazil with the largest cultivated area of cowpea were visited. In each municipality, a farmer was contacted and provided seeds. Landraces were collected and six preliminary assessments were conducted in small borderless plots. The 12 most promising landraces in terms of grain yield were evaluated in this study. The objectives were to identify the best cowpea landraces in terms of green grain yield and the characters that had the greatest direct positive effects over yield, via path analysis. The best landraces could thus be used directly by farmers or selected for breeding programs.

MATERIAL AND METHODS

Two experiments were carried out at the Rafael Fernandes Experimental Farm (5°03'S and 37°24' W, 87 m elevation) of the Federal Rural University of the Semi-Arid (UFERSA), located 20 km away from the city of Mossoró, state of Rio Grande do Norte (RN), Brazil. The first experiment was carried out between August and December 2013, and the second experiment was carried out from January to April 2014. According to Gaussen's bioclimatic classification, the climate in the Mossoró region is type 4aTH or distinctly xerothermic, i.e., tropical hot with a pronounced dry season lasting from seven to eight months and with a xerothermic index between 150 and 200. The average maximum air temperature in the region ranges between 32.1 and 34.5 °C, with June and July being the coldest months, and an average annual precipitation of approximately 825 mm. Insolation increases from March to October with a monthly average of 241.7 h. The maximum relative humidity is 78% in April and the minimum is 60% in September (CARMO FILHO; OLIVEIRA, 1989CARMO FILHO, F.; OLIVEIRA, O. F. Mossoró: um município do semi-árido nordestino. Mossoró: ESAM, 1989. 62p. (Coleção mossoroense, 672).). Table 1 shows climate data for each experimental period.

The soil in the experimental area is classified as red-yellow argisol according to the Brazilian Soil Classification System (EMBRAPA, 2006) and as ferric lixisol according to the Soil Map of the World (FAO, 1988). The results of the chemical analysis of the soil for the two experiments are shown in Table 2. Phosphorus (P) and potassium (K) were extracted with Mehlich 1 extracting solution (0.05M HC1 in 0.0125M H2SO4). P was colorimetrically determined from the extracted solution and K was determined by photometry. Calcium (Ca), magnesium (Mg), and aluminum (Al) were extracted with 1M of KC1. Aluminum was determined by titration with NaOH and Ca and Mg by titration with EDTA.

Before sowing, 30 kg ha-1 of ammonium sulfate (1/3 of the nitrogen dose), 100 kg ha-1 of P2O5 (simple superphosphate), and 50 kg ha-1 of K2O (potassium chloride) were applied as fertilizers in furrows located below and alongside the sowing furrows. The remaining N was blanket applied, in equal parts, after each of two weedings. A spacing of 1.0 m was used between rows, with the pits of each row spaced at 1.0 m. Sowing was done manually with four seeds per pit. Thinning was done 20 days after sowing, leaving the two largest plants in each pit. Hence, after thinning, each experiment had a planting density of 20000 plants ha-1.

The first experiment was sprinkler-irrigated, with the experimental plots arranged in parallel with the sprinklers. The required water depth for beans (5.3 mm) was calculated considering the effective depth of the root system as 0.40 m. The plots were irrigated every two days based on the water retained in the soil at a tension of 0.04 Mpa. Irrigation began after sowing and was suspended five days before the last harvest. The second experiment was carried out under rainfed conditions, but was irrigated whenever necessary. Irrigation was performed in the second experiment as described above. Pest control was conducted with two Deltamethrin sprays (250 ml ha-1), at seven and 14 days after sowing. Weed control was performed by hoe weeding at 20 and 40 days after sowing.

A randomized block experimental design with five replicates was used in the two experiments. The landraces evaluated were: Baraúna, Campo Grande, Carnaubais, Ceará Mirim, Itaú, José da Penha, Lagoa de Pedra, Lagoa Salgada, Macaíba, São Tomé, Umarizal, and Upanema. Landraces were named after the municipalities from where they were collected. Each plot consisted of three rows of plants, each row measuring 6.0 m in length. The usable area was occupied by the central row and plants in pits on either extremity of the row were not considered at harvesting.

Six harvests were performed at every two to four days for the evaluation of green grain yield. Green grain yield was estimated based on the weight of pods and green grains. In addition, the number of pods per plant (based on the pods collected in the usable area); number of grains per pod (in 20-pod samples); 100-green grain weight (in five samples); and the length, width, and thickness of pods and grains (in 10 sample units) were also determined.

Individual and joint analysis of variance were performed using the landrace data. The effects of landraces and experiments were considered as fixed effects in these analyses. Means were compared by the Scott-Knott test (1974) when the value of the F-test in the analysis of variance was significant (P < 0.05). The data were tested for homogeneity of variances before the analyses of variance, which were performed using Statistical Analysis System (SAS) software version 9.1 (SAS Institute 2003). Multicollinearity diagnostics and path analysis were performed using the Genes Program - Computational Genetics and Statistics Application (CRUZ, 2013CRUZ, C. D. GENES - A software package for analysis in experimental statistics and quantitative genetics. Acta Scientiarum Agronomy, Maringá, v. 35, n. 3, p. 271-276, 2013.).

In path analysis, the degree of multicollinearity of the matrix X'X was established based on its condition number (CN, the ratio between the largest and smallest eigenvectors of the matrix) and the determinant value of the correlation matrix between the characters analyzed. Multicollinearity does not cause serious problems in path analysis when the CN is smaller than 100 (TOEBE; CARGNELUTTI FILHO, 2013TOEBE, M.; CARGNELUTTI FILHO, A. Multicollinearity in path analysis of maize (Zea mays L.). Journal of Cereal Science, London, v. 57, n. 3, p. 453-462, May 2013.) and determinant values close to zero indicate strong associations between the studied characters, which may lead to bias in the estimates. Previous analyses were conducted and multicollinearity was detected in the first (CN = 204.76 and determinant value = 5.67 x 10-6) and second (CN = 289.34 and determinant value = 7.34 x 10-7) experiments, which characterizes moderate to strong multicollinearity according to Montgomery and Peck (1981MONTGOMERY, D. C.; PECK, E. A. Introduction to linear regression analysis. New York: John Wiley & Sons, 1981. 504p.). For this reason, path analysis was conducted under multicollinearity. In this method, a procedure similar to ridge regression analysis was used (CARVALHO; CRUZ, 1996CARVALHO, S. P.; CRUZ, C. D. Diagnosis of multicollinearity: assesssment of the condition of correlation matrices used in genetic studies. Revista Brasileira de Genética, Ribeirão Preto, v. 19, n. 3, p. 479-484, 1996.). Contrary to traditional path analysis, path analysis under multicollinearity includes a constant k in the X'X correlation matrix to reduce the variance associated with the least squares estimator of the path analysis (CARVALHO; CRUZ, 1996). Thus, the system of normal X'Xβ = X'Y equations becomes (X'X + k) β + X'Y. Twenty-one k values (k = 0.00; 0.05 ...1.00) were tested. The smallest values that determine the path coefficients (5.83 x 10-2 for the first experiment and 6.4 x 10-2 for the second experiment) were chosen as recommended by Carvalho and Cruz (1996CARVALHO, S. P.; CRUZ, C. D. Diagnosis of multicollinearity: assesssment of the condition of correlation matrices used in genetic studies. Revista Brasileira de Genética, Ribeirão Preto, v. 19, n. 3, p. 479-484, 1996.).

RESULTS AND DISCUSSION

Green pod yield was affected by landraces (L) and experiments (E), but not by the interaction between L and E (Table 3). Green pod yield was not different in the landraces Macaíba, São Tomé, Umarizal, and Upanema, but was lower than in the other eight landraces (Table 4). The average pod yield in experiment 1 (second semester 2013) was lower than the average yield in experiment 2 (first semester 2014). In addition, pod length was not affected by E and by the interaction between L and E (Table 3). Four of the most productive landraces (Baraúna, Campo Grande, Carnaubais, and Ceará Mirim) had the longest pods (Table 4). Santos et al. (2014SANTOS, A. et al. Correlations and path analysis of yield components in cowpea. Crop Breeding and Applied Biotechnology, Viçosa, v. 14, n. 2, p. 82-87, 2014.) found a positive correlation between green grain yield and pod length and large pods may have a larger number of grains and/or heavier grains, resulting in higher grain yield.

Pod width and thickness were affected by landraces (L) and experiments (E) and pod thickness was also affected by the L x E interaction (Table 3). Pod width and thickness were highest in the landrace Baraúna, but thicker pods were found in the landraces Campo Grande and Lagoa de Pedra in experiment 2 (Table 4). On average, pod width was larger in experiment 2 than in experiment 1. In addition, pod thickness in experiment 2 was larger than or equal to in experiment 1. The fact that pod width and thickness were larger in experiment 2 had no effect on pod yield, which was higher in experiment 1 (Table 4). Large grains generally are heavier than small grains. In experiment 1, grain dimensions and the number of pods per plant, on average, were greater than in experiment 2 (Table 6). Pericarp weight apparently did not contribute to the higher average pod yield in experiment-1 because pod yield (Table 4) and grain yield (Table 6) were similar across landraces.

Green grain yield and grain length were affected by landraces (L) and experiments (E), but not by the L x E interaction (Table 5). Grain yield was significantly lower in the landraces Macaíba, São Tomé, Umarizal, and Upanema than in the other eight landraces (Table 6). In addition, landraces with the highest pod yield (Table 4) also had the highest grain yield (Table 6). Grain length was largest in the landrace Baraúna (Table 6). On average, grain yield and grain length were higher in experiment 1 than in experiment 2. Grain width and thickness were affected by landraces (L) and experiments (E), but not by the interaction between L and E (Table 5). Grain width and thickness were largest in the landraces Baraúna and Campo Grande, but the landraces Carnaubais, Ceará Mirim, and Itaú were among those with the thickest grains (Table 6). On average, grain width and thickness were larger in experiment 1 than in experiment 2 (Table 6).

Information on green cowpea grain dimensions is not available. In Brazil and in various international markets, there is a preference for large mature grains (see review in FREIRE FILHO et al. 2012FREIRE FILHO, F. R. et al. Production, breeding and potential of cowpea crop in Brazil. Teresina: Embrapa Mid-North, 2012. 54 p. (Documentos 216).). However, cowpea grain classification scales are based on weight and not on grain dimensions (FREIRE FILHO et al. 2012). The following classification for a 100-grain weight of mature grains (in g) was proposed by Freire Filho et al. (2012): extra-small: ≤ 10; small: 10.2‒15; medium: 15.1‒25; large: 25.1‒30; extra-large: ˃ 30 g. On average, grain yield was higher and bean dimensions were larger in experiment 1 than in experiment 2 (Table 6), indicating that there is a positive correlation between green and mature bean dimensions. Consumers prefer green beans with large grains. Some landraces also had wide and thick grains, but only the Baraúna landrace was superior when considering the three grain dimensions (Table 6). This fact, coupled with its high yield, makes the Baraúna landrace of special interest to be used by farmers or in breeding programs.

The number of pods per plant, number of grains per pod, and 100-grain weight were affected by landraces (L) and experiments (E), and the latter two traits were also affected by the L x E interaction (Table 5). Landraces that produced the largest numbers of pods (Table 3) and grains (Table 6) also had the largest number of pods per plant (Table 6). Of the eight landraces with the highest green bean yield, only two (Baraúna and Campo Grande) also had the largest number of grains per pod and the largest 100-grain weight (Table 7). For number of grains per pod, the experimental conditions in the first semester (experiment 2) were unfavorable for some varieties, resulting in a significant interaction between landraces and experiments. For 100-grain weight, the interaction between E and L was caused by a variation in performance across landraces between the first and second experiments (Table 7).

The Baraúna landrace was one of the most productive in terms of pod yield (Table 4) and green grain yield (Table 6), and also had the largest grains (Table 6). These characteristics, in addition to ease of threshing, are among the most desirable traits of green pods for the production of fresh grains (ANDRADE et al., 2010). Thus, the Baraúna landrace is the most promising of all landraces evaluated in this study for direct use by farmers or for use in breeding programs. However, other landraces had similar yield and similarly wide (Campo Grande and Carnaubais) and thick (Campo Grande, Carnaubais, Ceará Mirim, and Itaú) grains as the Baraúna landrace. These landraces could also be potentially used by farmers or for breeding. Several characteristics other than those examined in our study must be considered before a landrace is selected by farmers. However, because the landraces evaluated have already been grown by farmers, they are more likely to be accepted due to their other desirable characteristics.

Landraces can be used for crop production and in breeding programs. Different traits are selected in breeding programs, including grain yield, drought resistance, and photosynthetic rate (PADI; EHLERS, 2007; HALL, GRANTZ, 1980HALL, A. E.; GRANTZ, D. A. Drought resistance of cowpea improved by selecting for early appearance of mature pods. Crop Science , Madison, v. 21, n. 3, p. 461-461, 1980.; BIRADAR; SALIMATH; PATIL, 2010BIRADAR, S.; SALIMATH, P. M.; PATIL, B. C. Land race as a source for improving photosynthetic rate and productivity in cowpea (Vigna unguiculata W.). Eletronic Journal of Plant Breeding, Coimbatore, v. 1, n. 4, p. 872-877, 2010.). In these programs, landraces can be used for hybridization with improved cultivars (BIRADAR; SALIMATH; PATIL, 2010; PADI; EHLERS, 2007) or subjected to selection (HALL; GRANTZ, 1980; HEDGE; MISHRA, 2009MISHRA, H. N.; DASH, S. R. A study on variability, character association and path analysis for green pod yield in yard long bean [Vigna unguiculata subsp sesquipedalis (L.) Verdcourt]. Vegetable Science, New Delhi, v. 36, n. 3, p. 319-322, 2009.) because they tend to exhibit large variability resulting from mutations, natural crossing, and the mechanical mixing of seeds. Mass selection and individual plant selection using the progeny test are potentially interesting for selecting commercial landraces (FREIRE FILHO et al., 2012FREIRE FILHO, F. R. et al. Production, breeding and potential of cowpea crop in Brazil. Teresina: Embrapa Mid-North, 2012. 54 p. (Documentos 216).). Padi and Ehlers (2007) concluded that the bulk-population and the single-seed descent methods would be more efficient than the pedigree method in developing cowpea landraces using traditional landraces.

The two experiments were conducted in nearby areas, but some chemical characteristics were more favorable in the soil of the second experimental area (Table 2). The fact that bean yield in the first experiment was higher than (10 landraces) or equal to (two landraces) in the second experiment suggests that climatic conditions were more favorable to cowpea growth in the first experiment (Table 1). Belane and Dakora (2010BELANE, A. K.; DAKORA, F. D. Symbiotic N2 fixation in 30 field -grown cowpea (Vigna unguiculata L. Walp.) genotypes in the Upper West Region of Ghana measured using 15N natural abundance. Biology and Fertility of Soils, Heildeberg, v. 46, n. 2, p. 191-198, 2010.) found that genotypes with the highest N2 fixation rates also had the largest biomass and the highest grain yield, indicating that several abiotic factors can affect symbiotic nitrogen fixation (MOHAMMADI et al., 2012MOHAMMADI, K. et al. Effective factors on biological nitrogen fixation. African Journal of Agricultural Research, Nairobi, v. 7, n. 12, p. 1782-1788, 2012.). In addition, N2 fixation depends on genotype and year of evaluation (BELANE; DAKORA, 2010BELANE, A. K.; DAKORA, F. D. Symbiotic N2 fixation in 30 field -grown cowpea (Vigna unguiculata L. Walp.) genotypes in the Upper West Region of Ghana measured using 15N natural abundance. Biology and Fertility of Soils, Heildeberg, v. 46, n. 2, p. 191-198, 2010.). However, the climatic conditions during the first experiment may have been favorable to the crop, the nitrogen-fixing bacteria, or both.

In path analysis, k values of 5.83 x 10-2 and 6.4 x104 x 10-2 were used in the first and second experiments, respectively. In the first and second experiments, the following variables were estimated, respectively: coefficient of determination (0.91 and 0.96), condition number (96.9 and 67.0), residual effects (0.31 and 0.2), and X'X matrix determinants (1.59 x 10-3 and 3.49 x 10-3). However, the coefficients of determination of the path analysis were measured, and they were even higher than those reported in some studies (ALMEIDA et al., 2014ALMEIDA, W. S. et al. Correlation and path analysis in components of grain yield of cowpea genotypes. Revista Ciência Agronômica, Fortaleza, v. 45, n. 4, p. 726-736, 2014.; SANTOS et al., 2014SANTOS, A. et al. Correlations and path analysis of yield components in cowpea. Crop Breeding and Applied Biotechnology, Viçosa, v. 14, n. 2, p. 82-87, 2014.). The high coefficients of determination indicate that the variation of the main character (green grain yield) was determined in large part by the explanatory variables.

In path analysis, no variance inflator factors (VIFs) > 10 were observed in the first and second experiments (Tables 8-10). According to Neter et al. (1983NETER, J.; WASSERMAN, W.; KUTNER, M. H. Applied linear regression models. Homewood: Richard D. Irwin Inc., 1983. 547p.), VIFs exceeding 10 are often regarded as indicating multicollinearity; thus, the regression coefficients may be associated with these values and some estimates may be very influenced by multicollinearity. The results of this study indicate that multicollinearity caused a minor bias in regression analysis and that reliable estimates were obtained (CARVALHO; CRUZ, 1996CARVALHO, S. P.; CRUZ, C. D. Diagnosis of multicollinearity: assesssment of the condition of correlation matrices used in genetic studies. Revista Brasileira de Genética, Ribeirão Preto, v. 19, n. 3, p. 479-484, 1996.).

In general, the explanatory variables had similar effects on the response variable in the two experiments, which were conducted in different times of the year (Table 8, table 9, table 10). The number of pods per plant was the major determinant of green grain yield, regardless of time of year (Table 8): green grain yield was positively affected by the number of pods per plant. Similar results were reported in other studies (MISHRA; DASH 2009MISHRA, H. N.; DASH, S. R. A study on variability, character association and path analysis for green pod yield in yard long bean [Vigna unguiculata subsp sesquipedalis (L.) Verdcourt]. Vegetable Science, New Delhi, v. 36, n. 3, p. 319-322, 2009., NWOFIA et al. 2013NWOFIA, G. E.; OGBONNA, N. D.; AGBO, C. U. Path analysis and heritability estimates of yield and yield components in vegetable cowpea as influenced by planting season. American-Eurasian Journal of Agriculture & Environmental Sciences, Dubai, v. 13, n. 9, p. 1283-1289, 2013., HITIKSHA et al. 2014HITIKSHA, K. P.; ACHARYA, R. R.; SHEETAL, R. P. Interrelationship among green pod yield and its yield contributing characters and biochemical characters vegetable cowpea [Vigna unguiculata (L.) Walp.]. Trends in Biosciences, New Delhi, v. 7, n. 23, p. 3972-3976, 2014., SAPARA; JAVIA, 2014SAPARA, G. K.; JAVIA, R. M. Correlation and path analysis in vegetable cowpea (Vigna unguiculata L.). International Journal of Plant Sciences, Amsterdam, v. 9, n. 1, p. 138-141, 2014.). This result is particularly interesting for selection purposes because the number of pods is relatively easier to determine compared to grain yield, besides the fact that the latter requires that green pods be shelled, which can be a labor-intensive process.

The high correlation between yield and number of pods per plant suggests that indirect selection for grain yield via number of pods per plant may be an interesting alternative for cowpea cultivation (UMAHARAN et al., 1997UMAHARAN, P.; ARIANAYAGAM, R. P.; HAQUE, S. Q. Genetic analysis of yield and its componentes in vegetable cowpea (Vigna unguiculata L. Walp.). Euphytica, Dordrecht, v. 96, n. 2, p. 207-213, 1997.). However, some studies reported low heritability estimates for number of pods per plant (ALIDU et al., 2013; SHIMELIS; SHIRINGANI, 2010SHIMELIS, H.; SHIRINGANI, R. Variance components and heritabilities of yield and agronomic traits among cowpea genotypes. Euphytica, Dordrecht, v. 176, n. 3, p. 383-389, 2010.; MOHAMMED et al., 2009MOHAMMED, M. S.; RUSSOM, Z.; ABDUL, S. D. Inheritance of hairiness and pod shattering, heritability and correlation studies in crosses between cultivated cowpea (Vigna unguiculata (L.) Walp.) and its wild (var. pubescens) relative. Euphytica, Dordrecht, v. 171, n. 3, p. 397-407, 2009.). Thus, direct selection for cowpea yield has been proposed as the best method for selecting cowpea landraces (SHIMELI; SHIRINGANI, 2010).

However, heritability for grain yield may be relatively low and close to zero in some cases, which may be a problem for direct selection (PADI; EHLERS, 2007; ALIDU et al., 2013). There are no studies comparing the effectiveness of direct and indirect selection for cowpea. However, in other crops, direct selection can be more, less, or as effective as indirect selection (ATLIN; FREY, 1989ATLIN, G. N.; FREY, K. J. Predicting the relative effectiveness of direct versus indirect selection for oat yield in three types of stress environments. Euphytica, Dordrecht, v. 44, n. 1, p. 137-142, 1989.; BRANCOURT-HULMEL et al., 2005BRANCOURT-HULMEL, M. et al. Indirect versus direct selection of winter wheat for low-input or high-input levels. Crop Science, Madison, v. 45, n. 4, p. 1427-1431, 2005.; KUMAR; BAHL, 1992KUMAR, J.; BAHL, P. N. Direct and indirect selection for yield in chickpea. Euphytica, Dordrecht, v. 60, n. 3, p. 197-199, 1992.; GRAVOIS; McNEW, 1992GRAVOIS, K. A.; McNEW, R. W. Genetic relationships among and selection for rice yield and yield components. Crop Science , Madison, v. 33, n. 2, p. 249-252, 1992.). An alternative approach is to combine selection for grain yield with selection for number of pods per plant using selection indexes (ROSIELLE; FREY, 1975ROSIELLE, A. A.; FREY, K. J. Estimates of selection parameters associated with harvest index in oat lines derived from a bulk population. Euphytica, Dordrecht, v. 24, n. 1, p. 121-131, 1975.; GRAVOIS, 1998).

Estimates of indirect effects for the number of pods per plant were, in some cases, better than estimates of direct effects observed for other traits (Table 8, Table 9 and Table 10). Even though 100-grain weight had a greater direct effect than the number of grains per pod, it was the only trait that was negatively correlated with most yield components. Other traits with moderate effects (> 0.5) such as pod and grain length and thickness negatively affected green grain yield. This effect was probably due to the well-known compensation that occurs between the main components of yield, i.e., varieties with larger pods tend to have fewer seeds per pod and/or lower grain weight.

CONCLUSIONS

Landraces collected in the municipalities of Baraúna, Campo Grande, Carnaubais, Ceará Mirim, Itaú, José da Penha, Lagoa de Pedra, and Lagoa Salgada had the highest green pod and green grain yields. The number of pods per plant had the greatest direct positive effect on yield. Average crop yield was higher in the dry season (experiment 1) than in the rainy season (experiment 2).

REFERENCES

  • ALÇIDU, M. S.; ATOKPLE, I. D. K.; AKROMAH, R. Genetic analysis of vegetative-stage drought tolerance in cowpea. Greener Journal of Agricultural Sciences, Geneva, v. 3, n. 6, p. 481-496, 2013.
  • ALMEIDA, W. S. et al. Correlation and path analysis in components of grain yield of cowpea genotypes. Revista Ciência Agronômica, Fortaleza, v. 45, n. 4, p. 726-736, 2014.
  • ATLIN, G. N.; FREY, K. J. Predicting the relative effectiveness of direct versus indirect selection for oat yield in three types of stress environments. Euphytica, Dordrecht, v. 44, n. 1, p. 137-142, 1989.
  • BELANE, A. K.; DAKORA, F. D. Symbiotic N2 fixation in 30 field -grown cowpea (Vigna unguiculata L. Walp.) genotypes in the Upper West Region of Ghana measured using 15N natural abundance. Biology and Fertility of Soils, Heildeberg, v. 46, n. 2, p. 191-198, 2010.
  • BERTINI, C. H. C. M. et al. Análise multivariada e índice de seleção na identificação de genótipos superiores de feijão-caupi. Acta Scientiarum: Agronomy, Maringá, v. 32, n. 4, p. 613-619, 2010.
  • BIRADAR, S.; SALIMATH, P. M.; PATIL, B. C. Land race as a source for improving photosynthetic rate and productivity in cowpea (Vigna unguiculata W.). Eletronic Journal of Plant Breeding, Coimbatore, v. 1, n. 4, p. 872-877, 2010.
  • BOUKAR, O.; FATOKUN, C. Strategies in cowpea breeding. In: ZERIHUN, T. (Ed.). New approaches to plant breeding of orphan crops in Africa. Rome: FAO, 2009. cap. 6, p. 69-92.
  • BRANCOURT-HULMEL, M. et al. Indirect versus direct selection of winter wheat for low-input or high-input levels. Crop Science, Madison, v. 45, n. 4, p. 1427-1431, 2005.
  • BRESEGHELLO, F.; COELHO, A. S. G. Traditional and modern plant breeding methods with examples in rice (Oriza sativa L.). Journal of Agricultural and Food Chemistry, Washington, v. 61, n. 35, p. 8277-8286, 2013.
  • CARMO FILHO, F.; OLIVEIRA, O. F. Mossoró: um município do semi-árido nordestino. Mossoró: ESAM, 1989. 62p. (Coleção mossoroense, 672).
  • CARVALHO, S. P.; CRUZ, C. D. Diagnosis of multicollinearity: assesssment of the condition of correlation matrices used in genetic studies. Revista Brasileira de Genética, Ribeirão Preto, v. 19, n. 3, p. 479-484, 1996.
  • CRUZ, C. D. GENES - A software package for analysis in experimental statistics and quantitative genetics. Acta Scientiarum Agronomy, Maringá, v. 35, n. 3, p. 271-276, 2013.
  • CRUZ, C. D.; REGAZZI, A. J.; CARNEIRO, P. C. S. Modelos biométricos aplicados ao melhoramento genético. 4. ed. Viçosa, MG: Editora UFV, 2012. v. 1, 514p.
  • EMBRAPA. Centro Nacional de Pesquisa do Solo. Sistema brasileiro de classificação de solos. 2. ed. Brasília, DF: Serviço de Produção de Informação, 2006. 306p.
  • FAO. Soil map of the world; revised legend. Rome: UNESCO, 1988. 119p.
  • FREIRE FILHO, F. R. et al. Production, breeding and potential of cowpea crop in Brazil. Teresina: Embrapa Mid-North, 2012. 54 p. (Documentos 216).
  • GBAGUIDI, A. A. et al. Diversity and agronomic performances of the cowpea (Vigna unguiculata Walp.) landraces in Southern Benin. International Research Journal of Agricultural Science and Soil Science, Aliero, v. 3, n. 4, p. 121-133, 2013.
  • GHALMI, N. et al. Morphological and molecular diversity within Algerian cowpea (Vigna unguiculata (L.) Walp.) landraces. Genetic Resources and Crop Evolution, Dordrecht, v. 57, n. 3, p. 371-386, 2010.
  • GRAVOIS, K. A. Optimizing selection for rough rice yield, head rice, and total milled rice. Euphytica, Dordrecht, v. 101, n. 2, p. 151-156, 1998.
  • GRAVOIS, K. A.; McNEW, R. W. Genetic relationships among and selection for rice yield and yield components. Crop Science , Madison, v. 33, n. 2, p. 249-252, 1992.
  • HALL, A. E.; GRANTZ, D. A. Drought resistance of cowpea improved by selecting for early appearance of mature pods. Crop Science , Madison, v. 21, n. 3, p. 461-461, 1980.
  • HEDGE, V. S.; MISHRA, S. K. Landraces of cowpea, Vigna unguiculata (L.) Walp., as potential sources of genes for unique characters in breeding. Genetics resources and Crop Evolution, Dordrecht, v. 56, n. 5, p. 615-627, 2009.
  • HITIKSHA, K. P.; ACHARYA, R. R.; SHEETAL, R. P. Interrelationship among green pod yield and its yield contributing characters and biochemical characters vegetable cowpea [Vigna unguiculata (L.) Walp.]. Trends in Biosciences, New Delhi, v. 7, n. 23, p. 3972-3976, 2014.
  • KUMAR, J.; BAHL, P. N. Direct and indirect selection for yield in chickpea. Euphytica, Dordrecht, v. 60, n. 3, p. 197-199, 1992.
  • MAKOI, J. H. J. R.; CHIMPHANGO, S. B. M.; DAKORA, F. D. Effect of legume plant density and mixed culture on symbiotic N2 fixation in five cowpea (Vigna unguiculata L. Walp.) genotypes in South Africa. Symbiosis, Normal, v. 48, n. 1, p. 57-67, 2009.
  • MISHRA, H. N.; DASH, S. R. A study on variability, character association and path analysis for green pod yield in yard long bean [Vigna unguiculata subsp sesquipedalis (L.) Verdcourt]. Vegetable Science, New Delhi, v. 36, n. 3, p. 319-322, 2009.
  • MOHAMMADI, K. et al. Effective factors on biological nitrogen fixation. African Journal of Agricultural Research, Nairobi, v. 7, n. 12, p. 1782-1788, 2012.
  • MOHAMMED, M. S.; RUSSOM, Z.; ABDUL, S. D. Inheritance of hairiness and pod shattering, heritability and correlation studies in crosses between cultivated cowpea (Vigna unguiculata (L.) Walp.) and its wild (var. pubescens) relative. Euphytica, Dordrecht, v. 171, n. 3, p. 397-407, 2009.
  • MONTGOMERY, D. C.; PECK, E. A. Introduction to linear regression analysis. New York: John Wiley & Sons, 1981. 504p.
  • NETER, J.; WASSERMAN, W.; KUTNER, M. H. Applied linear regression models. Homewood: Richard D. Irwin Inc., 1983. 547p.
  • NWOFIA, G. E.; OGBONNA, N. D.; AGBO, C. U. Path analysis and heritability estimates of yield and yield components in vegetable cowpea as influenced by planting season. American-Eurasian Journal of Agriculture & Environmental Sciences, Dubai, v. 13, n. 9, p. 1283-1289, 2013.
  • OGUNBODEDE, B. A. Comparison between three methods of determining the relationships between yield and eight of its components in cowpea, Vigna unguiculata L. Walp. Scientia Horticulturae, Amsterdam, v. 38, n. 3-4, p. 201-205, 1989.
  • PADI, F. K.; EHLERS, J. D. Effectiveness of early generation selection in cowpea for grain yield and agronomic characteristics in semiarid West Africa. Crop Science , Madison, v. 48, n. 2, p. 533-540, 2008.
  • PEKSEN, A. Fresh pod yield and some pod characteristics of cowpea (Vigna unguiculata L. Walp.) genotypes from Turkey. Asian Journal of Plant Sciences, Deira, v. 3, n. 3, p. 269-273, 2004.
  • POLEGRI, L.; NEGRI, V. Molecular markers for promoting agro-biodiversity conservation: a case study from Italy. How cowpea landraces were saved from extinction. Genetic Resources and Crop Evolution , Dordrecht, v. 57, n. 6, p. 867-880, 2010.
  • ROSIELLE, A. A.; FREY, K. J. Estimates of selection parameters associated with harvest index in oat lines derived from a bulk population. Euphytica, Dordrecht, v. 24, n. 1, p. 121-131, 1975.
  • SANTOS, A. et al. Correlations and path analysis of yield components in cowpea. Crop Breeding and Applied Biotechnology, Viçosa, v. 14, n. 2, p. 82-87, 2014.
  • SAPARA, G. K.; JAVIA, R. M. Correlation and path analysis in vegetable cowpea (Vigna unguiculata L.). International Journal of Plant Sciences, Amsterdam, v. 9, n. 1, p. 138-141, 2014.
  • SHIMELIS, H.; SHIRINGANI, R. Variance components and heritabilities of yield and agronomic traits among cowpea genotypes. Euphytica, Dordrecht, v. 176, n. 3, p. 383-389, 2010.
  • TOEBE, M.; CARGNELUTTI FILHO, A. Multicollinearity in path analysis of maize (Zea mays L.). Journal of Cereal Science, London, v. 57, n. 3, p. 453-462, May 2013.
  • UMAHARAN, P.; ARIANAYAGAM, R. P.; HAQUE, S. Q. Genetic analysis of yield and its componentes in vegetable cowpea (Vigna unguiculata L. Walp.). Euphytica, Dordrecht, v. 96, n. 2, p. 207-213, 1997.
  • WRIGHT, S. Path coefficients and path analysis: alternative or complementary. Biometrics, Washington, v. 16, n. 2, p. 189-202, 1960.
  • 1
    Paper extracted from the masters thesis of the first author.

Publication Dates

  • Publication in this collection
    Oct-Dec 2016

History

  • Received
    21 Apr 2015
  • Accepted
    17 June 2016
Universidade Federal Rural do Semi-Árido Avenida Francisco Mota, número 572, Bairro Presidente Costa e Silva, Cep: 5962-5900, Telefone: 55 (84) 3317-8297 - Mossoró - RN - Brazil
E-mail: caatinga@ufersa.edu.br