Abstract
The objective of this work was to compare different methods to estimate heritability in 30 pigeon pea families from the F3, F4, F5, and F6 generations, for nine variables. The experimental design was a randomized complete block with three replicates and 20 plants per plot. Broad-sense heritability was estimated by the analysis of variance ( A N O V A) [ h 2b-E(MS)], restricted maximum likelihood/best linear unbiased prediction (REML/BLUP) (h2b-REML), parent-offspring regression (h2PO), and standard deviation unit (h2UP). The h2b-E(MS) and h2b-REML estimates were similar for seven of the analyzed variables. For a higher genetic control and easier selection, values of h2b-E(MS) and h2b-REML >0.70 were estimated for two variables in four generations, two variables in three generations, three variables in two generations, and one variable in one generation. Values of h2UP and h2PO >0.70 were obtained for four and five variables, respectively. The estimates via regression or parent-offspring correlation showed some values outside the expected range of 0 to 1. The ANOVA [h2b-E(MS)] and REML/BLUP [h2b-REML] methods are the best to estimate pigeon pea heritability.
Index terms:
Cajanus cajan
; expected mean square; genetic parameters; mixed models; parent-offspring regression
Resumo
O objetivo deste trabalho foi comparar diferentes métodos para estimar a herdabilidade em 30 linhagens de guandu das gerações F3, F4, F5 e F6, para nove variáveis. O delineamento experimental foi em blocos ao acaso, com três repetições e 20 plantas por parcela. Estimou-se a herdabilidade no sentido amplo por meio de análise de variância (ANOVA) [h2b-E(MS)], máxima verossimilhança restrita/melhor predição linear não viciada (REML/BLUP) (h2b-REML), regressão pai-filho (h2PO) e unidade do desvio-padrão (h2UP). As estimativas de h2b-E(MS) e h2b-REML foram de magnitude próxima para sete das variáveis analisadas. Para maior controle genético e facilidade na seleção, valores de h2b-E(MS) e h2b-REML >0,70 foram estimados para duas variáveis em quatro gerações, duas variáveis em três gerações, três variáveis em duas gerações e uma variável em uma geração. Valores de h2UP e h2PO >0,70 foram obtidos para quatro e cinco variáveis, respectivamente. As estimativas via regressão ou correlação pai-filho mostraram alguns valores fora da variação esperada de 0 a 1. Os métodos ANOVA [h2b-E(MS)] e REML/BLUP [h2b-REML] são os melhores para estimar a herdabilidade em guandu.
Termos para indexação:
Cajanus cajan
; quadrado médio esperado; parâmetros genéticos; modelos mistos; regressão pai-filho
Introduction
Pigeon pea [Cajanus cajan (L.) Millsp.] is an important pulse in some countries, especially for the populations of tropical and subtropical regions, due to its protein content and ability to withstand drought (Cullis & Kunert, 2017CULLIS, C.; KUNERT, K.J. Unlocking the potential of orphan legumes. Journal of Experimental Botany, v.68, p.1895-1903, 2017. DOI: https://doi.org/10.1093/jxb/erw437.
https://doi.org/10.1093/jxb/erw437...
). The species is mainly grown in India, which concentrates 75% of its worldwide production, estimated at 5.6 million per hectare (FAO, 2019FAO. Food and Agriculture Organization of the United Nations. Faostat. 2019. Available at: <https://www.fao.org/faostat/en/#home>. Accessed on: Nov. 7 2021.
https://www.fao.org/faostat/en/#home...
).
Despite the advances in global pigeon pea production due to breeding, so far yield gains have been small, and, therefore, increasing them is still one of the main objectives and challenges of the breeders of the species (Saxena et al., 2021aSAXENA, K.; BOHRA, A.; CHOUDHARY, A.K.; SULTANA, R.; SHARMA, M.; PAZHAMALA, L.T.; SAXENA, R.K. The alternative breeding approaches for improving yield gains and stress response in pigeonpea (Cajanus cajan). Plant Breeding, v.140, p.74-86, 2021a. DOI: https://doi.org/10.1111/pbr.12863.
https://doi.org/10.1111/pbr.12863...
). According to Silva et al. (2018)SILVA, F.M. da; PEREIRA, E. de M.; VAL, B.H.P.; PERECIN, D.; DI MAURO, A.O.; UNÊDA-TREVISOLI, S.H. Strategies to select soybean segregating populations with the goal of improving agronomic traits. Acta Scientiarum. Agronomy, v.40, e39324, 2018. DOI: https://doi.org/10.4025/actasciagron.v40i1.39324.
https://doi.org/10.4025/actasciagron.v40...
, the estimates and analyses of important genetic parameters to predict yield gains in the selection of superior genotypes have allowed significant advances for several quantitative agronomic traits, with heritability acting as an essential parameter in classical breeding.
The main methods used to estimate heritability (h2) fall into three categories (Brown & Caligari, 2008BROWN, J.; CALIGARI, P.D.S. An introduction to plant breeding. Oxford: Wiley-Blackwell, 2008. 226p.). The first consists of methods based on variance components, such as the expected mean squares of the analysis of variance (ANOVA) and the restricted maximum likelihood (REML) (Hill, 2014HILL, W.G. Applications of population genetics to animal breeding, from Wright, Fisher and Lush to genomic prediction. Genetics, v.196, p.1-16, 2014. DOI: https://doi.org/10.1534/genetics.112.147850.
https://doi.org/10.1534/genetics.112.147...
). The second comprises methods based on the regression between the offspring and one or both parents, used to estimate narrow-sense heritability (Smith & Kinman, 1965SMITH, J.D.; KINMAN, M.L. The use of parent-ofspring regression as an estimator of heritability. Crop Science, v.5, p.595-596, 1965. DOI: https://doi.org/10.2135/cropsci1965.0011183X000500060035x.
https://doi.org/10.2135/cropsci1965.0011...
) and the standard unit heritability in terms of standard deviation (Frey & Horner, 1957FREY, K.J.; HORNER, T. Heritability in standard units. Agronomy Journal, v.49, p.59-62, 1957. DOI: https://doi.org/10.2134/agronj1957.00021962004900020001x.
https://doi.org/10.2134/agronj1957.00021...
). The third includes some methods that effectively express heritability achieved as a gain from selection, determining how much of the selection differential applied in previous generations was observed as a response in the progeny (Brown & Caligari, 2008BROWN, J.; CALIGARI, P.D.S. An introduction to plant breeding. Oxford: Wiley-Blackwell, 2008. 226p.).
Obala et al. (2018)OBALA, J.; SAXENA, R.K.; SINGH, V.K.; VECHALAPU, S.; DAS, R.; RATHORE, A.; SAMEER-KUMAR, C.V.; SAXENA, K.; TONGOONA, P.; SIBIYA, J.; VARSHNEY, R.K. Genetic variation and relationships of total seed protein content with some agronomic traits in pigeonpea (Cajanus cajan (L.) Millsp.). Australian Journal of Crop Science, v.12, p.1859-1865, 2018. DOI: https://doi.org/10.21475/ajcs.18.12.12.p1138.
https://doi.org/10.21475/ajcs.18.12.12.p...
used ANOVA to estimate broad-sense heritability (h2b) in pigeon pea, finding values from 0.519 for number of pods per plant to 0.999 for days to f lowering, with the h2b of 0.519 having par ticular importance for grain protein content. Studying five pigeon pea generations, Ajay et al. (2012)AJAY, B.C.; GNANESH, B.N.; GANAPATHY, K.N.; BYRE GOWDA, M.; PRASAD, P.S.; VEERAKUMAR, G.N.; VENKATESHA, S.C.; ABDUL FIYAZ, R.; RAMYA, K.T. Genetic analysis of yield and quantitative traits in pigeonpea (Cajanus cajan L. Millsp.). Euphytica, v.186, p.705-714, 2012. DOI: https://doi.org/10.1007/s10681-011-0556-1.
https://doi.org/10.1007/s10681-011-0556-...
reported a high heritability and genetic advance in yield-related traits, which are indicative of the high effectiveness of the performed selection. Furthermore, Sharma et al. (2019)SHARMA, S.; PAUL, P.J.; SAMEER KUMAR, C.V.; JAGANMOHAN RAO, P.; PRASHANTI, L.; MUNISWAMY, S.; SHARMA, M. Evaluation and identification of promising introgression lines derived from wild Cajanus species for broadening the genetic base of cultivated pigeonpea [Cajanus cajan (L.) Millsp.]. Frontiers in Plant Science, v.10, art.1269, 2019. DOI: https://doi.org/10.3389/fpls.2019.01269.
https://doi.org/10.3389/fpls.2019.01269...
used mixed models to evaluate pigeon pea lineages, observing a high broad-sense heritability (>70%) for most traits. However, to date, there are no known studies on heritability estimates for improved pigeon pea populations in Brazil or using different methods to estimate this parameter in the species.
The objective of this work was to compare different methods to estimate heritability in 30 pigeon pea families from the F3, F4, F5, and F6 generations, for nine variables.
Materials and Methods
The experiment was conducted at the Caatinga experimental field of Embrapa Semiárido, located in the municipality of Petrolina, in the state of Pernambuco, Brazil (09º09'S, 40º22'W, at 365.5 m above sea level). According to Köppen’s classification, the climate of the region is of the BSwh’ type, i.e., a hot semiarid. The soil of the experimental area was classified as an Ultisol and Oxisol (Soil Survey Staff, 2014SOIL SURVEY STAFF. Keys to soil taxonomy. 12th ed. Washington: USDA, 2014. 360p.), with the following properties according to Claessen (1997)CLAESSEN, M.E.C. (Org.). Manual de métodos de análise de solo. 2.ed. rev. e atual. Rio de Janeiro: Embrapa-CNPS, 1997. (EMBRAPA-CNPS. Documentos, 1).: pH 5.0, 3.28 mg dm-3 P, 0.22 cmolc dm-3 K, 0.03 cmolc dm-3 Na, 0.60 cmolc dm-3 Mg, 0.00 cmolc dm-3 Al, 1.5 cmolc dm-3 H + Al, sum of bases of 2.5 cmolc dm-3, cation exchange capacity of 4.3 cmolc dm-3, porosity of 38.27%, 48.8 g kg-1 clay, 225.5 g kg-1 silt, and 725.7 g kg-1 sand.
For the study, 30 families from different pigeon pea generations from the ICPL 90045 × UW10 cross were randomly selected and evaluated: six families from F3 and eight families from each of the F4, F5, and F6 generations. Accession ICPL 90045 was developed by International Crops Research Institute for the Semi-Arid Tropics (ICRISAT, Patancheru, India) as an extra-early plant, with maturation under 90 days, plant height lower than 100 cm, a 100-grain weight of ~9.0 g, three grains per pod, and pods measuring ~6.0 cm in length. Accession UW10 was developed at The University of the West Indies (St. Augustine, Trinidad and Tobago), with a cycle longer than 100 days, a 100-grain weight of ~11 g, five grains per pod, and pods measuring ~6.2 cm (Santos et al., 2000SANTOS, C.A.F.; ARAÚJO, F.P. de; MENEZES, E.A. Avaliação de genótipos de guandu de diferentes ciclos e portes no Sertão pernambucano. Magistra, v.12, p.31-40, 2000.).
To obtain the seeds of the families evaluated in the F3, F4, F5, and F6 generations, the plants from the previous generation were protected before anthesis with non-woven fabric bags (TNT) to avoid possible crosses. Two seeds were sown per cell in 200-cell polystyrene trays containing commercial substrate (Maxfertil, Pouso Redondo, SC, Brazil), and the seedlings were transplanted to the field after developing the first true leaves 15 days after sowing.
The experiment was conducted in randomized complete blocks, with three replicates and 20 plants per plot, distributed along two rows with 2.5 m of length, containing five holes each. During the experiment, the plants were irrigated and crop management practices were carried out whenever necessary according to Saxena et al. (2019)SAXENA, R.K.; SAXENA, K.B.; VARSHNEY, R.K. Pigeonpea (Cajanus cajan L. Millsp.): an ideal crop for sustainable agriculture. In: AL-KHAYRI, J.M.; JAIN, S.M.; JOHNSON, D.V. (Ed.). Advances in plant breeding strategies: legumes. Cham: Springer, 2019. v.7, p.409-429. DOI: https://doi.org/10.1007/978-3-030-23400-3_11.
https://doi.org/10.1007/978-3-030-23400-...
. No fertilizers were applied.
The following nine traits were evaluated: number of days to start flowering (NIF); number of days to 50% flowering (NDF); number of days to pod maturation (NDM), with 95% of the pods in the plot showing a typical mature pod color; mature plant height (PH), measured from the base of the plant to the tip of the main stem; mean pod length (PL) of five random pod samples, in millimeters; number of pods per plant (NPP); mean number of grains per pod (NGP) of five randomly sampled pods; 100-seed weight (P100), in grams; and grain yield, defined as the total weight of grains harvested per plant, in grams or kilograms. Selection was performed for all traits evaluated in each generation and was considered: positive, when an increasing mean was obtained for NPP, NGP, P100, and grain yield; and negative, when a decreasing mean was found for NIF, NDF, NDM, PH, and PL.
The data were tested for normality of residuals using the Kolmogorov-Smirnov test and the PROC GLM procedure of the SAS software (SAS Institute Inc., Cary, NC, USA). The variables NPP and grain yield did not meet the normality assumptions and were transformed to square root of x.
For PH, NPP, PL, NGP, P100, and grain yield, the variance components were calculated with the PROC VARCOMP procedure (SAS Institute Inc., Cary, NC, USA) using the expected mean square according to the following model: where Yijk is the value observed in the k-th plant, in the j-th replicate of the i-th treatment; μ is the overall mean; ti is the random effect of the i-th treatment (i=1, 2, ..., n;); bj is the random effect of the j-th block (j=1, 2, ..., r); eij is the random effect of the variation between plots; and πijk is the random effect of the variation between plants within the plot.
For the remaining traits – NIF, NDF, and NDM –, the parameters were estimated by the expected mean square using the following model: yij =m+bj+ti+eij, , where yij is the value of the studied trait in treatment i (i =1, 2, ..., I) and in block (or replicate) j (j =1, 2, ..., J); m is the overall mean (of all observations) of the experiment; ti is the effect of treatment i; and eij is the error associated with observation yij or the effect of uncontrolled factors. Individual ANOVAs by generations and a pooled ANOVA with all generations were performed to estimate genetic parameters and to test for phenotypic variance, respectively.
Broad-sense heritability by ANOVA (h2a-E(MS)), the genetic coefficient of variation between families (CVge), and the environmental coefficient of variation (CVe) were estimated as described by Cruz et al. (2014)CRUZ, C.D.; CARNEIRO, P.C.S.; REGAZZI, A.J. Modelos biométricos aplicados ao melhoramento genético. 3.ed. rev. e ampl. Viçosa: UFV, 2014. v.2, 668p.. According to Vencovsky & Barriga (1992)VENCOVSKY, R.; BARRIGA, P. Genética biométrica no fitomelhoramento. Ribeirão Preto: Sociedade Brasileira de Genética, 1992. 496p., when the CVge/CVe ratio is equal to or greater than 1.0, there is a favorable condition for selection.
Genetic gain as a percentage of the mean (GG%) was obtained and categorized as: low, from 0–10%; moderate, 10–20%; and high, >20% (Obala et al., 2018OBALA, J.; SAXENA, R.K.; SINGH, V.K.; VECHALAPU, S.; DAS, R.; RATHORE, A.; SAMEER-KUMAR, C.V.; SAXENA, K.; TONGOONA, P.; SIBIYA, J.; VARSHNEY, R.K. Genetic variation and relationships of total seed protein content with some agronomic traits in pigeonpea (Cajanus cajan (L.) Millsp.). Australian Journal of Crop Science, v.12, p.1859-1865, 2018. DOI: https://doi.org/10.21475/ajcs.18.12.12.p1138.
https://doi.org/10.21475/ajcs.18.12.12.p...
), using the following equation:
where is the mean of the evaluated phenotypes and k is the constant for selection intensity (1.27 at a selection intensity of 25%).
The estimates were obtained by the restricted maximum likelihood/best unbiased linear prediction (REML/BLUP) procedure. The genetic parameters were estimated for: PH, NPP, PL, NGP, P100, and grain yield, which are variables with plant-level data, using model 18 of the SELEGEN-REML/BLUP software (Resende, 2016RESENDE, M.D.V. de. Software Selegen-REML/BLUP: a useful tool for plant breeding. Crop Breeding and Applied Biotechnology, v.16, p.330-339, 2016. DOI: https://doi.org/10.1590/1984-70332016v16n4a49.
https://doi.org/10.1590/1984-70332016v16...
); and NIF, NDF, and NDM, with model 21 of the same software.
The several selections performed during generation advance made it impossible to estimate the additive and dominance variances and, therefore, narrow-sense heritability, as discussed by Vencovsky & Barriga (1992)VENCOVSKY, R.; BARRIGA, P. Genética biométrica no fitomelhoramento. Ribeirão Preto: Sociedade Brasileira de Genética, 1992. 496p.. Continuous self-fertilization resulted in an increased homozygosity and equal broad- and narrow-sense heritability values (Clo et al., 2019CLO, J.; GAY, L.; RONFORT, J. How does selfing affect the genetic variance of quantitative traits? An updated meta-analysis on empirical results in angiosperm species. Evolution, v.73, p.1578-1590, 2019. DOI: https://doi.org/10.1111/evo.13789.
https://doi.org/10.1111/evo.13789...
).
The parent-offspring regression method was used to estimate the heritability values (h2PO) between F6 and F3, offspring and parent, respectively. The estimates were obtained according to Smith & Kinman (1965)SMITH, J.D.; KINMAN, M.L. The use of parent-ofspring regression as an estimator of heritability. Crop Science, v.5, p.595-596, 1965. DOI: https://doi.org/10.2135/cropsci1965.0011183X000500060035x.
https://doi.org/10.2135/cropsci1965.0011...
: , where b is the regression coefficient and 2rop is the parent-offspring correlation coefficient (0.5 for F3). The analyses were performed using the PROC REG procedure (SAS Institute Inc., Cary, NC, USA). Narrow-sense heritability estimates were obtained by the standard deviation unit method (h2UP) of Frey & Horner (1957)FREY, K.J.; HORNER, T. Heritability in standard units. Agronomy Journal, v.49, p.59-62, 1957. DOI: https://doi.org/10.2134/agronj1957.00021962004900020001x.
https://doi.org/10.2134/agronj1957.00021...
, which shows the correlation between the means of the families of generations F6 and F3, offspring and parent, respectively. Heritability was estimated using the PROC CORR procedure (SAS Institute Inc., Cary, NC, USA).
Results and Discussion
All variables showed statistical differences, except PL and NGP, for pooled generations (Table 1). The CVs ranged from 8.93% for PL to 32.32% for grain yield, indicating a good to reasonable experimental control for most variables (Zavinon et al., 2019ZAVINON, F.; ADOUKONOU-SAGBADJA, H.; BOSSIKPONNON, A.; DOSSA, H.; AHANHANZO, C. Phenotypic diversity for agro-morphological traits in pigeon pea landraces [(Cajanus cajan L.) Millsp.] cultivated in southern Benin. Open Agriculture, v.4, p.487-499, 2019. DOI: https://doi.org/10.1515/opag-2019-0046.
https://doi.org/10.1515/opag-2019-0046...
).
Mean squares of number of days to start flowering in the plot (NIF), number of days to 50% flowering in the plot (NDF), number of days to pod maturation (NDM), mature plant height (PH), number of pods per plant (NPP), pod length (PL), number of grains per pod (NGP), 100-seed weight (P100), and grain yield (Yield) for 30 families of the F3, F4, F5 and F6 generations from the ICPL 90045 × UW10 pigeon pea (Cajanus cajan) cross.
The means of NIF, NDF, PH, NPP, and grain yield were lower in generation F3, with no significant differences for PL and NGP (Figure 1). With the exception of P100, none of the means of the analyzed variables differed significantly between the F4, F5, and F6 generations. With a high selection intensity, it is expected that the genetic gain of selection will be reduced after a few generations; therefore, as observed for the NIF and NDF variables, there is a greater chance that selection will be more effective when practiced in the initial generations for traits controlled by a greater number of genes, such as grain yield (Hallauer et al., 2010HALLAUER, A.R.; CARENA, M.J.; MIRANDA FILHO, J.B. Quantitative genetics in maize breeding. 3rd ed. New York: Springer, 2010. 663p. (Handbook of Plant Breeding, 6). DOI: https://doi.org/10.1007/978-1-4419-0766-0.
https://doi.org/10.1007/978-1-4419-0766-...
; Saxena et al., 2021bSAXENA, K.B.; WALLIS, E.S.; CHAUHAN, Y.S.; BYTH, D.E. Inheritance of photo-sensitivity in pigeonpea. Indian Journal of Genetics, v.81, p.1-7, 2021b. DOI: https://doi.org/10.31742/IJGPB.81.1.6.
https://doi.org/10.31742/IJGPB.81.1.6...
).
Heritability means obtained for the following variables evaluated in generations F3, F4, F5, and F6 from the ICPL 90045 × UW10 pigeon pea (Cajanus cajan) cross: number of days to start fowering in the plot (NIF), number of days to 50% fowering in the plot (NDF), number of days to 95% pod maturation (NDM), mature plant height (PH), number of pods per plant (NPP), pod length (PL), number of grains per pod (NGP), 100-seed weight (P100), and grain yield (Yield). Means followed by equal letters do not difer by the Bonferroni test, at 5% probability.
The present study highlights the high convergence of heritability in the broad-sense (h2b) estimated by ANOVA and REML-BLUP in generations F3, F4, F5, and F6 from the ICPL 90045 × UW10 pigeon pea cross, with similar values for 23 of the 36 estimates obtained by both methods (Figure 2). Among the various methodologies available to estimate the parameter, the one based on ANOVA and mixed models, such as REML-BLUP (You et al., 2016YOU, F.M.; JIA, G.; CLOUTIER, S.; BOOKER, H.M.; DUGUID, S.D.; RASHID, K.Y. A method of estimating broad-sense heritability for quantitative traits in the type 2 modifed augmented design. Journal of Plant Breeding and Crop Science, v.8, p.257-272, 2016. DOI: https://doi.org/10.5897/JPBCS2016.0614.
https://doi.org/10.5897/JPBCS2016.0614...
), helps the breeder to choose appropriate strategies to select superior genotypes.
Estimates of broad-sense heritability based on the expected mean square analysis of variance [E(MS)] and restricted maximum likelihood (REML) for the following variables evaluated in generations F3, F4, F5, and F6 from the ICPL 90045 × UW10 pigeon pea (Cajanus cajan) cross: number of d ays t o st ar t fowe r i ng i n t he plot ( NIF), nu mber of days to 50% fowering in the plot (NDF), number of days to 95% pod maturation (NDM), mature plant height (PH), number of pods per plant (NPP), pod length (PL), number of grains per pod (NGP), 100-seed weight (P100), and grain yield (Yield).
The heritability values obtained for the studied variables differed among generations. For NIF, heritability decreased from generation F3 to F4, but showed an increasing trend in F6, reaching a value close to 1.0. For NDF, heritability values increased from generation F3 to F4, decreasing in F5 and F6. For NDM, the estimated values of heritability oscillated, decreasing from F3 to F4 and increasing from F4 to F5, followed by another decrease from F5 to F6. The h2b estimates, except those obtained for NIF and NDM in F4, are in agreement with those reported in the literature on the crop (Ajay et al., 2014AJAY, B.C.; BYREGOWDA, M.; PRASHANTH BABU, H.; VEERA KUMAR, G.N.; REENA, M. Variability and transgressive segregation for yield and yield contributing traits in pigeonpea crosses. Electronic Journal of Plant Breeding, v.5, p.786-791, 2014.; Vohra et al., 2015VOHRA, A.; BAJPAI, G.C.; VERMA, S.K. Yield factor analysis in F4 and F5 progenies derived from interspecific hybridization between cultivated and wild pigeonpea [Cajanus cajan (L.) Millsp]. Legume Research: An International Journal, v.38, 2015. DOI: https://doi.org/10.5958/0976-0571.2015.00114.9.
https://doi.org/10.5958/0976-0571.2015.0...
). In pigeon pea, NDF is considered a trait with a high heritability value and has been used in crop improvement programs for germplasm classification (Zavinon et al., 2019ZAVINON, F.; ADOUKONOU-SAGBADJA, H.; BOSSIKPONNON, A.; DOSSA, H.; AHANHANZO, C. Phenotypic diversity for agro-morphological traits in pigeon pea landraces [(Cajanus cajan L.) Millsp.] cultivated in southern Benin. Open Agriculture, v.4, p.487-499, 2019. DOI: https://doi.org/10.1515/opag-2019-0046.
https://doi.org/10.1515/opag-2019-0046...
).
PH estimates increased from generation F3 to F4, with the highest values of 0.83 and 0.84 obtained by ANOVA and REML-BLUP, respectively, followed by a decrease in F5, with the lowest h2b values for this variable, and then a slight increase in F6. For NPP, there was an increase in the h2b values from generation F4 to F5, progressing since F3, which showed a negative value considered equal to zero via the ANOVA and to 0.03 using REML-BLUP, reaching 0.57 in generation F5 and later decreasing in F6. When evaluating the parents and progenies of the F1, F2, and F3 generations of four pigeon pea crosses, Ajay et al. (2012)AJAY, B.C.; GNANESH, B.N.; GANAPATHY, K.N.; BYRE GOWDA, M.; PRASAD, P.S.; VEERAKUMAR, G.N.; VENKATESHA, S.C.; ABDUL FIYAZ, R.; RAMYA, K.T. Genetic analysis of yield and quantitative traits in pigeonpea (Cajanus cajan L. Millsp.). Euphytica, v.186, p.705-714, 2012. DOI: https://doi.org/10.1007/s10681-011-0556-1.
https://doi.org/10.1007/s10681-011-0556-...
found high heritability values. According to these authors, the association of high h2b and genetic advance values is indicative of the possibility of a successful selection in later generations. In the present work, high heritability estimates in the broad sense were obtained for PH and NPP, highlighting the presence of the additive gene action.
Among the analyzed variables, only PL showed a continuous increase in h2b throughout the advance of all evaluated generations, with the lowest value in F3 and peak of 0.83 in F6; the greatest difference between the obtained estimates was due to the used model, i.e., the ANOVA and REML-BLUP, resulting in values of 0.24 and 0.05, respectively. In addition, NGP presented the smallest variation in h2b estimates between generations, with the highest value in F4.
For P100, the highest value among all h2b estimates was observed in generation F4, decreasing in F5 and F6, which showed the lowest value for this variable. For grain yield, the value for h2b was the lowest in generation F3 and the highest in F4, followed by successive decreases in F5 and F6. These values are higher than those found by Zavinon et al. (2019)ZAVINON, F.; ADOUKONOU-SAGBADJA, H.; BOSSIKPONNON, A.; DOSSA, H.; AHANHANZO, C. Phenotypic diversity for agro-morphological traits in pigeon pea landraces [(Cajanus cajan L.) Millsp.] cultivated in southern Benin. Open Agriculture, v.4, p.487-499, 2019. DOI: https://doi.org/10.1515/opag-2019-0046.
https://doi.org/10.1515/opag-2019-0046...
for pigeon pea grain yield.
The CVge estimated by the ANOVA and REML were similar (Figure 3 A and B). CVge values higher than those of CVe were only found for: NDF in generations F4, F5, and F6; PH in F3, F4, and F6; and P100 in F3, F4, and F5. A CVge/CVe ratio equal to or greater than 1.0 was obtained for NDM, grain yield, NIF, and PL in generations F5, F4, F5, and F6, respectively. For NPP and NGP, the ratio was <1.0 and >1.0, respectively. A CVge/CVe ratio greater than 1.0 favors selection for trait improvement in general, allowing the breeder to make inferences about genetic variability and expected genetic gain with selection (Lopes et al., 2017LOPES, K.V.; TEODORO, P.E.; SILVA, F.A.; SILVA, M.T.; FERNANDES, R.L.; RODRIGUES, T.C.; FARIA T.C.; CORRÊA, A.M. Genetic parameters and path analysis in cowpea genotypes grown in the Cerrado/Pantanal ecotone. Genetics and Molecular Research, v.16, gmr16029559, 2017. DOI: https://doi.org/10.4238/gmr16029559.
https://doi.org/10.4238/gmr16029559...
).
Estimates of the genetic coefficient of variation (CVge) and environmental coefficient of variation (CVe) for the following variables evaluated in generations F3, F4, F5, and F6 from the ICPL 90045 × UW10 pigeon pea (Cajanus cajan) cross: number of days to start fowering in the plot (NIF), number of days to 50% fowering in the plot (NDF), number of days to 95% pod maturation (NDM), mature plant height (PH), number of pods per plant (NPP), pod length (PL), number of grains per pod (NGP), 100-seed weight (P100), and grain yield (Yield), using: A, the analysis of variance (ANOVA); B, restricted maximum likelihood (REML); and C, genetic gain as a percentage of the mean (GG%).
GG% was considered low for most variables (Obala et al., 2018OBALA, J.; SAXENA, R.K.; SINGH, V.K.; VECHALAPU, S.; DAS, R.; RATHORE, A.; SAMEER-KUMAR, C.V.; SAXENA, K.; TONGOONA, P.; SIBIYA, J.; VARSHNEY, R.K. Genetic variation and relationships of total seed protein content with some agronomic traits in pigeonpea (Cajanus cajan (L.) Millsp.). Australian Journal of Crop Science, v.12, p.1859-1865, 2018. DOI: https://doi.org/10.21475/ajcs.18.12.12.p1138.
https://doi.org/10.21475/ajcs.18.12.12.p...
), being <10% for all of them in generation F3 (Figure 3 C) and also for NIF, NDF, NDM, and PL in all generations; high values of >20% were only obtained for grain yield in F5. A moderate GG% of 10–20% was found for: NPP in generation F5, grain yield in F4, NGP in F4 and F6, P100 in F4 and F5, and PH in F4 to F6. According to Obala et al. (2018)OBALA, J.; SAXENA, R.K.; SINGH, V.K.; VECHALAPU, S.; DAS, R.; RATHORE, A.; SAMEER-KUMAR, C.V.; SAXENA, K.; TONGOONA, P.; SIBIYA, J.; VARSHNEY, R.K. Genetic variation and relationships of total seed protein content with some agronomic traits in pigeonpea (Cajanus cajan (L.) Millsp.). Australian Journal of Crop Science, v.12, p.1859-1865, 2018. DOI: https://doi.org/10.21475/ajcs.18.12.12.p1138.
https://doi.org/10.21475/ajcs.18.12.12.p...
, the association of h2b with other parameters, such as CVge and GG%, results in a greater reliability for estimating the expected genetic gain in phenotypic selection. In the present study, grain yield showed a high genetic gain with selection and low heritability estimates, stressing the importance of parameter association in the breeding process.
The heritability values estimated by the standard deviation unit method (h2UP) were higher than those obtained by the parent-offspring regression (h2PO) for six of the nine variables evaluated in generations F3 and F6 (Table 2), as observed for NPP and PL. The values estimated by h2UP and h2PO for NDF and PH were similar via the ANOVA and REML. However, unlike for h2b, the proportion of h2UP and h2PO estimates higher than 0.70 was low. The greatest difference between the values estimated by the two used methods was found for NDM, whose h2UP was approximately nine times higher than that of h2PO. Among the studied variables, PH showed the highest value for h2UP and h2PO.
Heritability estimates obtained by the standard deviation unit method (h2UP) and parent-offspring regression (h2PO) using data from generations F3 and F6 from the ICPL 90045 × UW10 pigeon pea (Cajanus cajan) cross(1).
For NGP, the h2UP estimate was high, above 1.0 (Table 2). Grain yield also showed a h2PO higher than 1.0, with a h2UP estimate of 0.80. For Cavalli-Sforza & Feldman (1976)CAVALLI-SFORZA, L.L.; FELDMAN, M.W. Evolution of continuous variation: direct approach through joint distribution of genotypes and phenotypes. Proceedings of the National Academy of Sciences, v.73, p.1689-1692, 1976. DOI: https://doi.org/10.1073/pnas.73.5.1689.
https://doi.org/10.1073/pnas.73.5.1689...
, heritability values higher than 1.0 occur when the intensity of selection is so high that selection variance becomes small and phenotypic variance becomes smaller than genetic variance. In this context, Frey & Horner (1957)FREY, K.J.; HORNER, T. Heritability in standard units. Agronomy Journal, v.49, p.59-62, 1957. DOI: https://doi.org/10.2134/agronj1957.00021962004900020001x.
https://doi.org/10.2134/agronj1957.00021...
proposed a modification in the procedure to estimate heritability by regression, according to which, regardless of the environmental effects, the estimated value will be close to 1.0, as found in the present study.
Similarly to NPP in generation F3, two negative h2PO values were estimated for P100 using the ANOVA; REML was the only method that did not result in negative estimates (Tables 1 and 2). According to Maia et al. (2016)MAIA, M.C.C.; MACEDO, L.M.; VASCONCELOS, L.F.L.; AQUINO, J.P.A.; OLIVEIRA, L.C.; RESENDE, M.D.V. Estimates of genetic parameters using RELM/BLUP for intra-populational genetic breeding of Platonia insignis Mart. Revista Árvore, v.40, p.561-573, 2016. DOI: https://doi.org/10.1590/0100-67622016000300020.
https://doi.org/10.1590/0100-67622016000...
, this is a characteristic of the REML methodology, which shows a non-negativity restriction for the estimates of the variance components. Moreover, PH and NPP presented h2b estimates higher than 80 and 90%, respectively. These results suggest a greater probability of success in early selection for traits with high h2b values, e.g., PH, NGP, P100, and NDF in all evaluated generations (Figure 2), estimated both by the ANOVA and REML/BLUP. This observation was confirmed by the high h2UP and h2PO estimates obtained for NDF, NGP, and PH (Table 2).
Conclusions
-
The broad-sense heritability estimates obtained by the analysis of variance (ANOVA) and restricted maximum likelihood (REML) methods for the evaluated families of different pigeon pea (Cajanus cajan) generations are similar in terms of magnitude.
-
The estimates obtained by the regression or parent-offspring correlation are precise for estimating heritability in pigeon pea.
-
The broad-sense heritability values estimated for the four studied pigeon pea generations are similar or close via the ANOVA and REML/best linear unbiased prediction, showing that both methods are appropriate for studies on this species.
Acknowledgements
To Embrapa Semiárido, for financial support; to Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (Capes), for financing, in part, this study (Finance Code 001); and to Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq), for granted fellowship.
References
- AJAY, B.C.; BYREGOWDA, M.; PRASHANTH BABU, H.; VEERA KUMAR, G.N.; REENA, M. Variability and transgressive segregation for yield and yield contributing traits in pigeonpea crosses. Electronic Journal of Plant Breeding, v.5, p.786-791, 2014.
- AJAY, B.C.; GNANESH, B.N.; GANAPATHY, K.N.; BYRE GOWDA, M.; PRASAD, P.S.; VEERAKUMAR, G.N.; VENKATESHA, S.C.; ABDUL FIYAZ, R.; RAMYA, K.T. Genetic analysis of yield and quantitative traits in pigeonpea (Cajanus cajan L. Millsp.). Euphytica, v.186, p.705-714, 2012. DOI: https://doi.org/10.1007/s10681-011-0556-1
» https://doi.org/10.1007/s10681-011-0556-1 - BROWN, J.; CALIGARI, P.D.S. An introduction to plant breeding Oxford: Wiley-Blackwell, 2008. 226p.
- CAVALLI-SFORZA, L.L.; FELDMAN, M.W. Evolution of continuous variation: direct approach through joint distribution of genotypes and phenotypes. Proceedings of the National Academy of Sciences, v.73, p.1689-1692, 1976. DOI: https://doi.org/10.1073/pnas.73.5.1689
» https://doi.org/10.1073/pnas.73.5.1689 - CLAESSEN, M.E.C. (Org.). Manual de métodos de análise de solo 2.ed. rev. e atual. Rio de Janeiro: Embrapa-CNPS, 1997. (EMBRAPA-CNPS. Documentos, 1).
- CLO, J.; GAY, L.; RONFORT, J. How does selfing affect the genetic variance of quantitative traits? An updated meta-analysis on empirical results in angiosperm species. Evolution, v.73, p.1578-1590, 2019. DOI: https://doi.org/10.1111/evo.13789
» https://doi.org/10.1111/evo.13789 - CRUZ, C.D.; CARNEIRO, P.C.S.; REGAZZI, A.J. Modelos biométricos aplicados ao melhoramento genético 3.ed. rev. e ampl. Viçosa: UFV, 2014. v.2, 668p.
- CULLIS, C.; KUNERT, K.J. Unlocking the potential of orphan legumes. Journal of Experimental Botany, v.68, p.1895-1903, 2017. DOI: https://doi.org/10.1093/jxb/erw437
» https://doi.org/10.1093/jxb/erw437 - FAO. Food and Agriculture Organization of the United Nations. Faostat 2019. Available at: <https://www.fao.org/faostat/en/#home>. Accessed on: Nov. 7 2021.
» https://www.fao.org/faostat/en/#home - FREY, K.J.; HORNER, T. Heritability in standard units. Agronomy Journal, v.49, p.59-62, 1957. DOI: https://doi.org/10.2134/agronj1957.00021962004900020001x
» https://doi.org/10.2134/agronj1957.00021962004900020001x - HALLAUER, A.R.; CARENA, M.J.; MIRANDA FILHO, J.B. Quantitative genetics in maize breeding 3rd ed. New York: Springer, 2010. 663p. (Handbook of Plant Breeding, 6). DOI: https://doi.org/10.1007/978-1-4419-0766-0
» https://doi.org/10.1007/978-1-4419-0766-0 - HILL, W.G. Applications of population genetics to animal breeding, from Wright, Fisher and Lush to genomic prediction. Genetics, v.196, p.1-16, 2014. DOI: https://doi.org/10.1534/genetics.112.147850
» https://doi.org/10.1534/genetics.112.147850 - LOPES, K.V.; TEODORO, P.E.; SILVA, F.A.; SILVA, M.T.; FERNANDES, R.L.; RODRIGUES, T.C.; FARIA T.C.; CORRÊA, A.M. Genetic parameters and path analysis in cowpea genotypes grown in the Cerrado/Pantanal ecotone. Genetics and Molecular Research, v.16, gmr16029559, 2017. DOI: https://doi.org/10.4238/gmr16029559
» https://doi.org/10.4238/gmr16029559 - MAIA, M.C.C.; MACEDO, L.M.; VASCONCELOS, L.F.L.; AQUINO, J.P.A.; OLIVEIRA, L.C.; RESENDE, M.D.V. Estimates of genetic parameters using RELM/BLUP for intra-populational genetic breeding of Platonia insignis Mart. Revista Árvore, v.40, p.561-573, 2016. DOI: https://doi.org/10.1590/0100-67622016000300020
» https://doi.org/10.1590/0100-67622016000300020 - OBALA, J.; SAXENA, R.K.; SINGH, V.K.; VECHALAPU, S.; DAS, R.; RATHORE, A.; SAMEER-KUMAR, C.V.; SAXENA, K.; TONGOONA, P.; SIBIYA, J.; VARSHNEY, R.K. Genetic variation and relationships of total seed protein content with some agronomic traits in pigeonpea (Cajanus cajan (L.) Millsp.). Australian Journal of Crop Science, v.12, p.1859-1865, 2018. DOI: https://doi.org/10.21475/ajcs.18.12.12.p1138
» https://doi.org/10.21475/ajcs.18.12.12.p1138 - RESENDE, M.D.V. de. Software Selegen-REML/BLUP: a useful tool for plant breeding. Crop Breeding and Applied Biotechnology, v.16, p.330-339, 2016. DOI: https://doi.org/10.1590/1984-70332016v16n4a49
» https://doi.org/10.1590/1984-70332016v16n4a49 - SANTOS, C.A.F.; ARAÚJO, F.P. de; MENEZES, E.A. Avaliação de genótipos de guandu de diferentes ciclos e portes no Sertão pernambucano. Magistra, v.12, p.31-40, 2000.
- SAXENA, K.; BOHRA, A.; CHOUDHARY, A.K.; SULTANA, R.; SHARMA, M.; PAZHAMALA, L.T.; SAXENA, R.K. The alternative breeding approaches for improving yield gains and stress response in pigeonpea (Cajanus cajan). Plant Breeding, v.140, p.74-86, 2021a. DOI: https://doi.org/10.1111/pbr.12863
» https://doi.org/10.1111/pbr.12863 - SAXENA, K.B.; WALLIS, E.S.; CHAUHAN, Y.S.; BYTH, D.E. Inheritance of photo-sensitivity in pigeonpea. Indian Journal of Genetics, v.81, p.1-7, 2021b. DOI: https://doi.org/10.31742/IJGPB.81.1.6
» https://doi.org/10.31742/IJGPB.81.1.6 - SAXENA, R.K.; SAXENA, K.B.; VARSHNEY, R.K. Pigeonpea (Cajanus cajan L. Millsp.): an ideal crop for sustainable agriculture. In: AL-KHAYRI, J.M.; JAIN, S.M.; JOHNSON, D.V. (Ed.). Advances in plant breeding strategies: legumes. Cham: Springer, 2019. v.7, p.409-429. DOI: https://doi.org/10.1007/978-3-030-23400-3_11
» https://doi.org/10.1007/978-3-030-23400-3_11 - SHARMA, S.; PAUL, P.J.; SAMEER KUMAR, C.V.; JAGANMOHAN RAO, P.; PRASHANTI, L.; MUNISWAMY, S.; SHARMA, M. Evaluation and identification of promising introgression lines derived from wild Cajanus species for broadening the genetic base of cultivated pigeonpea [Cajanus cajan (L.) Millsp.]. Frontiers in Plant Science, v.10, art.1269, 2019. DOI: https://doi.org/10.3389/fpls.2019.01269
» https://doi.org/10.3389/fpls.2019.01269 - SILVA, F.M. da; PEREIRA, E. de M.; VAL, B.H.P.; PERECIN, D.; DI MAURO, A.O.; UNÊDA-TREVISOLI, S.H. Strategies to select soybean segregating populations with the goal of improving agronomic traits. Acta Scientiarum. Agronomy, v.40, e39324, 2018. DOI: https://doi.org/10.4025/actasciagron.v40i1.39324
» https://doi.org/10.4025/actasciagron.v40i1.39324 - SMITH, J.D.; KINMAN, M.L. The use of parent-ofspring regression as an estimator of heritability. Crop Science, v.5, p.595-596, 1965. DOI: https://doi.org/10.2135/cropsci1965.0011183X000500060035x
» https://doi.org/10.2135/cropsci1965.0011183X000500060035x - SOIL SURVEY STAFF. Keys to soil taxonomy 12th ed. Washington: USDA, 2014. 360p.
- VENCOVSKY, R.; BARRIGA, P. Genética biométrica no fitomelhoramento Ribeirão Preto: Sociedade Brasileira de Genética, 1992. 496p.
- VOHRA, A.; BAJPAI, G.C.; VERMA, S.K. Yield factor analysis in F4 and F5 progenies derived from interspecific hybridization between cultivated and wild pigeonpea [Cajanus cajan (L.) Millsp]. Legume Research: An International Journal, v.38, 2015. DOI: https://doi.org/10.5958/0976-0571.2015.00114.9
» https://doi.org/10.5958/0976-0571.2015.00114.9 - YOU, F.M.; JIA, G.; CLOUTIER, S.; BOOKER, H.M.; DUGUID, S.D.; RASHID, K.Y. A method of estimating broad-sense heritability for quantitative traits in the type 2 modifed augmented design. Journal of Plant Breeding and Crop Science, v.8, p.257-272, 2016. DOI: https://doi.org/10.5897/JPBCS2016.0614
» https://doi.org/10.5897/JPBCS2016.0614 - ZAVINON, F.; ADOUKONOU-SAGBADJA, H.; BOSSIKPONNON, A.; DOSSA, H.; AHANHANZO, C. Phenotypic diversity for agro-morphological traits in pigeon pea landraces [(Cajanus cajan L.) Millsp.] cultivated in southern Benin. Open Agriculture, v.4, p.487-499, 2019. DOI: https://doi.org/10.1515/opag-2019-0046
» https://doi.org/10.1515/opag-2019-0046
Publication Dates
-
Publication in this collection
21 Oct 2022 -
Date of issue
2022
History
-
Received
01 Mar 2022 -
Accepted
20 June 2022