African Journal of
Agricultural Research

  • Abbreviation: Afr. J. Agric. Res.
  • Language: English
  • ISSN: 1991-637X
  • DOI: 10.5897/AJAR
  • Start Year: 2006
  • Published Articles: 6918

Full Length Research Paper

Genetic parameters in Stylosanthes using different statistical methods

Ronaldo Simão de Oliveira
  • Ronaldo Simão de Oliveira
  • Universidade Estadual de Feira de Santana - UEFS, Brasil.
  • Google Scholar
Manoel Abílio de Queiróz
  • Manoel Abílio de Queiróz
  • Universidade do Estado da Bahia - UNEB, Brasil.
  • Google Scholar
Roberto Lisboa Romão
  • Roberto Lisboa Romão
  • Universidade Estadual de Feira de Santana - UEFS, Brasil.
  • Google Scholar
Bruno Augusto de Souza Almeida
  • Bruno Augusto de Souza Almeida
  • Universidade do Estado da Bahia - UNEB, Brasil.
  • Google Scholar
Cláudio Mistura
  • Cláudio Mistura
  • Universidade do Estado da Bahia - UNEB, Brasil.
  • Google Scholar
Luciano Paganucci de Queiróz
  • Luciano Paganucci de Queiróz
  • Universidade Estadual de Feira de Santana - UEFS, Brasil.
  • Google Scholar


  •  Received: 21 August 2015
  •  Accepted: 25 September 2015
  •  Published: 12 November 2015

 ABSTRACT

Genetic parameters in Stylosanthes accessions were estimated through (ANOVA) and REML/BLUP (Restricted Maximum Likelihood /Best Linear Unbiased Prediction), to compare them for the genetic values in order to select superior accessions. Twenty five genotypes were evaluated in two environments in a randomized blocks experimental design with four replications. The genetic parameters were estimated for 12 descriptors by the two methods. Both methods indicating that the accessions presented genetic variability for the descriptors, but, the ANOVA and the REML/BLUP presented divergent values for the unbalanced data and for descriptors with high environmental influence, as it happens with mass descriptors, what, in turn, indicates that the method REML/BLUP leads to a more accurate predictions and allows for the selection of Stylosanthes accessions. It allows the inclusion of characters, even if they present heterogeneity of the residual variance. The indirect selection may be used, because the primary stem length descriptor revealed itself as a good option due to the correlation with the total dry mass. The accessions BGF-016 and BGF-015 are the most promising ones to be taken up by a Stylosanthes plant breeding program for fodder in the Brazilian Semiarid region.

 

Key words: Semiarid, forages, breeding


 INTRODUCTION

The Stylosanthes genus calls our attention due to its number of species with extraordinary fodder potential with excellent nutritional quality and their easy adaptation to different environmental conditions. Even though Brazil is the centre of diversity of these species, Australia, African and Asian countries have released cultivars with higher impact on the production system than the releases made in Brazil (Resende et al., 2006). In Brazil the most used cultivar is the  Estilosantes  Campo  Grande,  which belongs to the species Stylosanthes capitata (Embrapa, 2007), however this variety was not developed for the Semiarid. A survey done in the herbariums from the universities of Bahia and Embrapa Genetic Resources and Biotechnology indicated the existence of different species in the State of Bahia, revealing the great diversity for the genus (Costa, 2006), and thus demonstrating that there is a great potential for the study of this germplasm, so far quite neglected, in order to develop cultivars for the Brazilian Semiarid.
 
The procedures normally used for the selection of fodder plants is done by the evaluation of accessions and the selection of the superior ones (Chakraborty, 2004). However, this selection is more efficient when it is based on the genetic values estimate and in association with variables of economic importance (Assis et al., 2010). Due to the lack of Stylosanthes cultivars developed for the region, it is desirable to establish a plant breeding program for this fodder plant, and thus the genetic parameters estimation is needed to identify the best strategies to be used in the selection of Stylosanthes individuals which are tolerant to droughts, has high fodder productivity and high nutritional value.
 
For a long time, the variance analysis (ANOVA) method has been the most used in plant breeding to estimate variance components, but, in situations where there are unbalanced data, environmental variance, variance heterogeneity among experiments, competition between genotypes due to the difference in aggressively and sensibility of different genetic materials (Resende, 2007) its use is limited.
 
This study had the objective to estimate the genetic parameters through the statistic models (ANOVA and REML/BLUP), to select the best method of identifying the most prominent accessions in the selection for superior materials for starting a breeding program of this fodder in the Semiarid of Bahia, which is also very important for in the Brazilian semiarid.


 MATERIALS AND METHODS

Two field experiments were carried out between July 2012 and January 2013. The first one was done at the experimental station Horto Florestal (12° 16’087’’S; 38° 56’346”W; 243 m) which belongs to the State University of  Feira de Santana - UEFS in Feira de Santana, BA (A1) and the other one at the experimental field of the Technology and Social Sciences Department (09° 24’ 50”S; 40°30’ 10”W; 368 m) of the University of Bahia - UNEB in Juazeiro, BA (A2).
 
For the experimental set up, 25 accessions of Stylosanthes spp., from collections acquired between the years of 2008-2011 in three semiarid regions in the state of Bahia (Sisaleira, the Middle Lower São Francisco River Basin and the micro region of Feira de Santana) were used along with a commercial variety, Estilosantes Campo Grande, as a control (Table 1).The fruits of each accession obtained in the collecting expeditions were sent to the Evolutive Ecology Laboratory - LEE from the Estate University of Feira de Santana, Feira de Santana, BA for manual processing, consisting of the removal of loment by friction using rubber devices of four mm of thickness. The seeds from each accession were put in labelled envelopes, and kept in air tight bottles with silica gel as a humidity indicator and were kept in Lab conditions (temperature of 25 °C). Between the months of April and May 2012, 54 seeds of each accession were mechanically scarified with wood sandpaper n° 150 and sown in polyethylene tubes with dimensions of 6 x 20 cm filled with the commercial substrate. They were kept under greenhouse conditions, and watered twice daily. The seedlings in both experiments were produced in green house of the State University of Feira de Santana (UEFS) with the purpose of standardizing the procedures for obtaining the seedlings and after three  months  they were transplanted to the field.
 
 
The soil, in both experiments, was prepared by disk harrow and the marking of the plots were marked and, then pitting was performed with the aid of a manual digger. A complete randomized block design with four replications and four plants per useful area of each plot was used, with a total of 16 plants per treatment and the spacing of 3.0 m between rows and 0.8 m between plants was adopted. There was a basal fertilization, in both trials, with the application of 30 kg ha-1 P, 30 kg ha-1 K and 20 kg N, and after around 35 to 40 days of the transplanting, topdressing fertilization with 30 kg ha-1 K and 20 kg ha-1N was used. Both experiments were irrigated using the drip irrigation system whenever necessary.
 
The evaluations were performed between four and five months after transplanting. For morphological descriptors the measurements were made on all the plants of useful area from each treatment. They were: PD - Plant Diameter (mm) measured at the base of the plant; PH – Plant Height (cm), measured from the ground level to the highest leaf on the stem; PBL - Primary branch length (cm), its measurement was taken from the insertion of the primary branch at the bottom of the stem up to the last leaf; PAL- Length of the central axis of the plants (cm), measured from the ground level to the highest leaf of the main stem; NS - Number of stems (units), counting of the number of stems from ground level to the last stem inserted in the main axis; CLL - central leaflet length (mm) and CLW - width of the central leaflet (mm), measured in the longitudinal part of the central leaflet of the third definitive leaf  of the plant, inserted in the first ten centimeters of the central axis; SLL - Side leaflet length (mm) and LLW- Lateral width of leaves (mm), measured in the longitudinal and latitudinal part, respectively, of the right of the third leaf stage of plant, inserted in the first ten centimeters from the central axis lateral leaflet. The TDM descriptors (total dry mass - g), DSM (dry stalk mass - g) and DLM (dry leaf mass - g) were evaluated in one plant per plot.
 
Initially the individual variance analysis was done for each place to check if the accessions differed significantly among themselves. Then, the homogeneity of variance was tested (F maximum - ratio between the largest and smallest mean squared residual for each descriptor) following the test of Hartley (1950) recommended by Cruz et al. (2004), excluding descriptors that presented F maximum superior to 7.
 
In this study, each accession was  formed by a group of individuals which represented a given population, and therefore, the genotypes effects were regarded as random and the environment effects as fixed using the following model (Cruz et al., 2004): Yijk = μ + (B / A)jk + Gi + Aj + GAij + Eijk where: Yijk = observation in the k-the block, reported in the i-the genotype and j-the environment; μ is the general mean of the test; (B/A) the block effect k within the environment j; Gi = effect of genotype i; Aj = effect of environment j; GAij = effect of the interaction between genotype i and environment j; and Eijk = random error associated with the ijk observation.
 
From the joint analysis, the genetic parameters were estimated (Cruz et al., 2004), thus the decomposition of the mean squares estimate of the interaction in simple and complex parts was measured (Cruz and Castoldi, 1991).
 
Genetic parameters were also estimated using the methodology of mixed models REML/BLUP (Restricted Maximum Likelihood /Best Linear Unbiased Prediction) to observe possible differences from ANOVA (Resende, 2007). In this analysis each individual has an individual genotypic value allowing for a more accurate estimate of the genetic value and more adequate ordering of superior individuals, leading to a better selection (Martinez et al., 2011).
 
In the analysis of mixed models with unbalanced data such as the present one, the effects of the model were tested via LRT (likelihood ratio test) over the F test which is used in the analysis of variance and in replacement of the framework of ANOVA a similar framework called analysis of deviance (ANADEV) was performed. All these effects along with the estimates of the variance components   and   genetic   parameters   were   obtained   by    the SELEGEN-REML/BLUP software (model 23 - selection based on various replications, taking into consideration one or more plants per plot, evaluation in more than one location and experiments in completely randomized blocks) .
According to Resende (2007) using the following statistical model y = Xr + Zg + Wi + e, where:y= data vector;ris the vector of the effects assumed to be fixed (repetition + overall mean);gis the vector of genotypic effects,iis the effects vector of the genotype x environment interaction, where: eis the vector of errors or residues (random), where: , and Ware the incidence matrices for the referred effects. Based on the genetic parameters it was possible to estimate the selection genetic gain calculated by the formulae where is the gain in environment i, with selection based on environment j; dsi (j) is the selection differential in the environment i in which the individuals selected have the best performance in the environment j; and is the heritability of the character in the environment i, adopting a selection intensity of 20% for each analysis. To estimate the gain selection percentage it was used the following formulae: .


 RESULTS AND DISCUSSION

The joint analysis, using ANOVA, showed that the mean squares of  the  effects  of  genotypes  and  environments presented high (P≤0.01) for almost all variables except lateral leaflet width (LLW), for environment (P≤0.05). The effects of G x E interaction were also great (P≤0.01), except plant diameter (PD) and lateral leaflet width (LLW) (Table 2). When considering the values between the two environments, most of the descriptors presented F maximum below seven (Table 2), except the total dry mass (TDM), dry stalk mass (DSM) and dry leaf mass (DLM) descriptors and the genetic variances (Vg) were higher in comparison to variances of the interaction GxE (Vgxe) and environmental variance (Ve) for almost all descriptors except for side leaflet length (SLL) (Table 3) indicating genetic diversity among the accessions evaluated and also demonstrated that the genetic variability among the accessions  is partially due to differences among individuals within each accession and, also, due to characteristics of different species, since the 25 treatments comprised four species (Stylosanthes scabra, Stylosanthes viscosa, Stylosanthes capitata and Stylosanthes humilis).
 
 
 
The relative coefficient of variation (CVr) was almost equal to 1.0 side leaflet length (SLL) and higher than 1.0 for   the    remaining    variables.    As    for    broad-sense heritability (h2g), the values ranged from 0.41 to 0.95. The results also showed that all characters expressed interaction of the simple type, except for central leaflet width (CLW) and plant diameter (PD) that were not significant and, therefore, the correlation between environments (r12) were relatively high (Table 3).
On the other hand, the values of ANADEV (Analysis of Deviance) show that seven descriptors plant diameter (PD), plant height (PH), primary branch length (PBL), length of the central axis of the plants (PAL), number of stems (NS), central leaflet length (CLL) and central leaflet width (CLW) presented effects of genetic variances (Vg) and broad-sense heritability (h2g) high (P≤0.01) and the heritability ranged from 0.59 to 0.83. Regarding the descriptors side leaflet length (SLL), lateral leaflet width (LLW), total dry mass (TDM), dry stalk mass (DSM) and dry leaf mass (DLM), the effects of genotypes were not significant via the LRT test. For the interaction effects of GxE and for the coefficient of determination (r2int) the values of the LRT test for all descriptors were high (P≤0.01), except for plant diameter (PD) and central leaflet width (CLW) (Tables 4 and 5). 
 
 
 
When comparing values for the variance components obtained by ANOVA with the ones acquired by REML/BLUP, there was divergence in value, once for the ANOVA method the central leaflet length (CLL) and lateral leaflet width (LLW) were high (P≤0.01), whereas for the REML/BLUP these descriptors showed no significant differences. The mean values of the descriptors that were possible to assess through ANOVA were higher when compared to those obtained by REML/BLUP, thus indicating elevation of means by ANOVA in the joint analysis. As such, this may induce a wrong selection of the best individuals.
 
It is noteworthy that the descriptors mass total dry mass (TDM), dry stalk mass (DSM) and dry leaf mass (DLM), the most important ones from the forage point of view, had to be excluded from the joint analysis by ANOVA due to this method`s limitations on heterogeneity of residual variances (Cruz et al., 2004), thus preventing the evaluation of their genetic parameters, the potential in relation to the genetic gains and consequently the indication of possible accessions for the selection of superior materials for the forage production.
 
Therefore, the methodology REML/BLUP, despite being a method little used for annual plants, can be used to detect variability among accessions in different species of Stylosanthes particularly for descriptors that exhibit heterogeneity of variance and unbalanced data, thus showing, that these factors are not limiting to estimate genetic predictions for descriptors that have a strong environmental influence. This reassures this methodology as essential for breeding programs of forages that need to be grown in different environments, as it is the case for Stylosanthes in the Semiarid of Bahia. The methodology of mixed models was also used efficiently to quantify the variation and estimate genetic predictions in some annual species (Borges et al., 2010; Neto et al., 2013; Resende et al., 2006).
 
It is also worth mentioning that apart from mass descriptors, other relevant traits in Stylosanthes need to be considered since in a breeding program the accessions selected by genetic merits of mass production should also be assessed in the association with grasses, produce enough seeds and present tolerance to diseases such as anthracnose – Colletotrichum gloeosporioides (Penz.) Penz. and Sacc. (Resende et al., 2006), although not yet reported in the state of Bahia. On the other hand, the mass descriptors were not significant among the accessions evaluated using REML/BLUP, but the broad-sense heritability and mean of dry matter found in this study was similar to the values found for Stylosanthes (Resende et al., 2006; Santana et al., 2013).
 
The CVg% ranged from 17.38 to 44.90%, with the lowest value observed for the lateral leaflet width (LLW) descriptor, while the highest variation was for dry stalk mass (DSM); the environment coefficient of variation (CVe) were superior for all descriptors, except for the mass ones; the relative  variation  coefficient  (CVr),  were above a unit for all descriptors, except for mass descriptors (Table 5). Regarding the mean heritability of the genotype (h2mg), the values obtained were high for all the descriptors, although lower for mass descriptors. A very high precision, the selective accuracy (rgg, Resende and Duarte, 2007), was observed for most descriptors of the plant and only two of them (side leaflet length -SLL and lateral leaflet width -LLW) presented only high accuracy, however, the mass descriptors showed moderate accuracy.
 
The genetic correlations (rg) between the accessions and the environments showed that the plant diameter (PD), plant height (PH), primary branch length (PBL), length of the central axis of the plants (PAL), number of stems (NS), central leaflet length (CLL) and central leaflet width (CLW) descriptors showed higher coefficients than 0.78 and plant diameter (PD), central leaflet width (CLW) and primary branch length (PBL) were highlighted. On the other hand, the characters side leaflet length (SLL), lateral leaflet width (LLW), total dry mass (TDM), dry stalk mass (DSM) and dry leaf mass (DLM) displayed genetic correlation between both environments(rg) values below 0.50 but the lowest levels were observed for total dry mass (TDM), dry stalk mass (DSM) and dry leaf mass (DLM) (Table 5).
 
From the values of genetic variances, heritability and relative coefficient of variation obtained in this work for mass descriptors, it can be inferred that the genotypes demonstrated potential for obtaining genetic progress with the selection in the environments of study. By comparing the results found here with other annual forages commonly used in the Brazilian Semiarid region, the potential of the genus Stylosanthes is confirmed (Assis et al., 2008, 2010; Cunha and Lima, 2010).
 
Moreover, the values of heritability in the broad sense were also overestimated by the method ANOVA, because nearly all descriptors were superior to the REML/BLUP values, except for the central leaflet length (CLL), confirming that the method of ANOVA for more complex situations leads to inaccurate estimates of the variance components and consequently inaccurate predictions of breeding values (Resende, 2007).
 
The behavior of the variance components when decomposed to G x E interaction by the two methods (ANOVA and REML/BLUP) were similar to some descriptors except side leaflet length (SLL) and lateral leaflet width (LLW) that showed simple interaction type by ANOVA and complex interaction by REML/BLUP. Then, the choice of the best accessions should happen through the REML/BLUP, because it was more viable for descriptors that exhibit different behavior in different environments (side leaflet length - SLL), lateral leaflet width - LLW, total dry mass - TDM, dry stalk mass - DSM and dry leaf mass - DLM). This method facilitates the evaluation and selection in this situation, similar to what happens with other  forages (Cargnin et al., 2006; Luz et al., 2010;  Pereira  et  al.,  2010)  and  studies  using  this methodology in selecting genotypes (Arantes et al., 2013; Bastos et al., 2007; Zeni-Neto et al., 2008; Verardi et al., 2009).
 
The estimate of the simple correlations revealed that the majority of descriptors showed significant correlations (P≤0.01), but of moderate to low magnitude from one character to the other, except for plant height (PH) and length of the central axis of the plants (PAL) that showed no significance with the majority of the characters. For the primary branch length (PBL) it was observed a positive association from moderate to low intensity for most descriptors, however, there was a significant and moderate correlation between PBL with the mass descriptors (Table 6).
 
 
Therefore, a possibility to improve the selection of genotypes with efficiency for higher dry matter would be through indirect selection. Moreover, indirect selection can lead to superior gains over direct  selection,   when   the   auxiliary   character displays heritability above the main one, a positive correlation and of high magnitude (Falconer and Mackay, 1996). Therefore, the indirect selection can be used, once the primary branch length (PBL) is a good option for selection due to the correlation with total dry mass. This is a feature of great significance to the breeding of Stylosanthes in the Brazilian Semiarid region with the purpose of increasing the total dry mass, since there is a great variation in the environmental conditions.
 
The selection of the best five accessions by ANOVA generated total gains ranging from the minimum and the maximum equal to 20.81%, and total dry mass (TDM) presented the smallest increase, while the largest was observed in central leaflet length (CLL). By the REML/BLUP method the amplitude of variation was approximately 38.02%, where the central leaflet length (CLL) had the highest and the lowest value was for dry  leaf  mass   (DLM).   As   for   primary branch length (PBL), it was observed that the two methods presented similar high values to gain selection and the estimation of genetic gains by the two methods showed that the BGF-016 and BGF-015 accessions were the top among the materials evaluated (Table 7).
 
Thus, the REML/BLUP method allowed the elimination of accessions that had negative genotypic effects (values below the overall mean of the experiment), and this increased the probability of selecting superior individuals; it provided the estimate for the number of individuals to be selected by accession for each variable and also enabled the identification of the plots of the best individuals and the number to be selected in each one (Martinez et al., 2011; Resende et al., 2006). 
 
 
Besides all these advantages presented here, it is known that the REML/BLUP method is little used  for  annual  species  in  the   area   of   plant

genetic resources, but the results presented here show that this methodology can be used routinely in the selection of accessions since it is easily handled, it estimates the genetic values with greater accuracy, increases the efficiency of selection, and consequently decreases the cost ofbreeding programs to increase mass descriptors for the species of Stylosanthes.
 

 


 CONCLUSIONS

The methods (ANOVA and REML/BLUP) presented divergent values, indicating that the REML/BLUP method estimated the genetic values with greater accuracy, increases the efficiency of selection and decreases the cost of breeding programs. The indirect selection may be used, because the primary stem length descriptor revealed itself as a good option due to the correlation with the total dry mass. The accessions BGF-016 and BGF-015 are the most promising ones to be taken up by a Stylosanthes plant breeding program for fodder in the Brazilian Semiarid region.


 CONFLICT OF INTERESTS

The authors have not declared any conflict of interests.


 ACKNOWLEDGEMENTS

The Coordination of Improvement of Higher Education Personnel (CAPES) for the scholarship grant; the State University of Feira de Santana (UEFS) and the Department of Technology and Social Sciences (DTCS/UNEB), for the support in carrying out the field experiments.



 REFERENCES

Arantes FC, Junior EJS, Gonçalves OS, Moraes MLT, Gonçalves ECP, Resende MDV (2013). Adaptability and stability in rubber tree progenies under different environmental conditions. Pesquisa Florestal Brasileira 33:37-44. doi: 10.4336/2013.pfb.33.73.436
Crossref

 

Assis GML, Valentim JF, Júnior JMC, Azevedo JMA, Ferreira AS (2008). Seleção de genótipos de amendoim forrageiro para cobertura do solo e produção de biomassa aérea no período de estabelecimento utilizando-se metodologia de modelos mistos. Revista Brasileira de Zootecnia 37:1905-1911. 
Crossref

 
 

Assis LCSL, Lira MA, Santos MVF, Júnior JCBD, Cunha MV (2010). Estimativa de parâmetros genéticos sob duas estratégias de avaliação em híbridos intra e interespecíficos de capim-elefante. Revista Brasileira de Zootecnia 39:2589-2597. 
Crossref

 
 

Bastos IT, Barbosa MHP, Resende MDV, Peternelli LA, Silveira LCI and Donda, LR, Fortunato AA, Costa PMA and Figueiredo ICR (2007). Avaliação da interação genótipo x ambiente em cana-de-açúcar via modelos mistos. Pesquisa Agropecuária Tropical 37:195-203.

 
 

Borges V, Soares AA, Reis MS, Resende MDV, Cornélio VMO, Leite NA, Vieira AR (2010). Desempenho genotípico de linhagens de arroz de terras altas utilizando metodologia de modelos mistos. Bragantia 69:833-841. 
Crossref

 
 

Cargnin A, Souza MA, Carneiro PCS, Sofiatti V (2006). Interação entre genótipos e ambientes e implicações em ganhos com seleção em trigo. Pesquisa Agropecuária Brasileira 41:987-993. 
Crossref

 
 

Chakraborty S (2004). High-yielding anthracnose resistant Stylosanthes for agriculturalsystems. (Australian Centre for International Agricultural Research, ACIAR: Canberra).

 
 

Costa NMS (2006). Revisão do gênero Stylosanthes Sw. Tese (Doutorado em Engenharia Agronômica) - Universidade Técnica de Lisboa, Lisboa, Portugal. 469p.

 
 

Cruz CD, Castoldi FL (1991). Decomposição da interação genótipos x ambientes em partes simples e complexa. Revista Ceres 38:422-430.

 
 

Cruz CD, Regazzi AJ, Carneiro PCS (2004). Modelos biométricos aplicados ao melhoramento genético. Editora UFV: Viçosa, MG.

 
 

Cunha EE, Lima JMP (2010). Caracterização de genótipos e estimativa de parâmetros genéticos de características produtivas de sorgo forrageiro. Revista Brasileira de Zootecnia 39:701-706.
Crossref

 
 

Embrapa (2007). "Cultivo e uso do Estilosantes-campo-grande." Embrapa Gado de Corte, Comunicado técnico No.105, Campo Grande, MS.

 
 

Falconer DS, Mackay TFC (1996). Introduction to quantitative genetics.Essex: Longman.

 
 

Luz LN, Santos RC, Filho JLS and Melo Filho PA (2010). Estimativas de parâmetros genéticos em linhagens de amendoim baseadas em descritores associados ao ginóforo. Revista Ciência Agronômica 41:132-138.

 
 

Martinez DT, Resende MDV, Higa AR, Costa RB (2011). Procedimentos de predição e efeitos de heterogeneidade de variâncias residuais dentro de tratamentos genéticos. Pesquisa Florestal Brasileira 31:193-202. 
Crossref

 
 

Neto AR, Junior EUR, Gallo PB, Freitas JG, Azzini LE (2013). Comportamento de genótipos de arroz de terras altas no estado de São Paulo. Revista Ciência Agronômica 44:512-519. 
Crossref

 
 

Pereira HS, Melo LC, Faria LC, Peloso MJDP, Wendland A (2010). Estratificação ambiental na avaliação de genótipos de feijoeiro-comum tipo Carioca em Goiás e no Distrito Federal. Pesquisa Agropecuária Brasileira 45:554-562. 
Crossref

 
 

Resende MDV (2007). "SELEGEN-REML/BLUP". Embrapa Florestas: Colombo, PR.

 
 

Resende MDV de, Duarte JB (2007). Precisão e controle de qualidade em experimentos de avaliação de cultivares. Pesquisa Agropecuária Tropical 37:182-194.

 
 

Resende RMS, Resende MDV, Laura VA (2006).Genotypic evaluation of accessions and individual selection in Stylosanthes spp. by simulated BLUP method. Crop Breed. Appl. Biotechnol. 6:253-260.
Crossref

 
 

Santana AS, Oliveira RS, Romão RL, Brasileiro BP, Generozo DB (2013). Divergência genética entre acessos de Stylosanthes Sw. (Fabaceae) coletados no Semiárido Baiano. Magistra 24:304-313.

 
 

Verardi CK, Resende MDV, Costa RB, Gonçalves PS (2009). Adaptabilidade e estabilidade da produção de borracha e seleção em progênies de seringueira. Pesquisa Agropecuária Brasileira 44:1277-1282. 
Crossref

 
 

Zeni-Neto H, Oliveira RA, Daros E, Bespalhok-Filho JC, Zambon JLC, Ido OT, Weber H (2008). Seleção para produtividade, estabilidade e adaptabilidade de clones de cana-de-açúcar em três ambientes no estado do Paraná via modelos mistos. Scientia Agraria 09:425-430.
Crossref

 

 




          */?>