Adaptability and stability of strawberry cultivars using a mixed model

Adaptabilidade e estabilidade de cultivares de morango via modelos mistos

ABSTRACT

Although strawberry crops have a strong socio-economic impact on the agricultural sector of Espírito Santo State, there are few studies on the performance of strawberry cultivars in different locations and years under a low tunnel management system (LT). Therefore, this study aimed to estimate the parameters of adaptability and stability of strawberry cultivars under LT-protected cultivation using the harmonic mean of the relative performance of genetically predicted values. Seven strawberry cultivars were assessed ('Dover', 'Camino Real', 'Ventana', 'Camarosa', 'Seascape', 'Diamante' and 'Aromas') in the agricultural years 2006/7, 2007/8 and 2008/9 in three locations in the mountainous region of Espírito Santo State. The experiment was arranged in a randomized block design, with 3 replications and 15 plants per plot. The analysis of deviance for the yield (ton. ha^-1) demonstrated that only the effects of genotype and the triple interaction genotype x location x year were significant, which indicates the presence of genetic variability among the cultivars and the inconsistency of the position among the genotypes for combinations of year and location. Considering the selection for yield, adaptability and stability, in LT-protected cultivation, cultivars Camarosa and Aromas are highlighted for expressing average values that were 22% higher than the overall mean of the cultivars (24.55 ton. ha^-1) in the environments assessed.

Keywords
Fragaria x ananassa; REML/BLUP; MHPRVG method; genotype x environment interaction

RESUMO

A cultura do morangueiro vem exercendo um importante papel sócio-econômico no setor agrícola do Estado do Espírito Santo, devido principalmente ao aumento de renda nas pequenas propriedades e na fixação de trabalhadores na área rural. Contudo, são escassos os estudos sobre o desempenho de cultivares de morango em diferentes locais e anos sob o sistema de manejo túnel baixo (TB). Portanto, o presente trabalho teve como objetivo avaliar a produtividade, adaptabilidade e estabilidade de cultivares de morango sob cultivo protegido TB, utilizando o método da média harmônica da performance relativa dos valores genéticos preditos. Foram avaliados sete cultivares de morango ('Dover', 'Camino Real', 'Ventana', 'Camarosa', 'Seascape', 'Diamante' e 'Aromas') nos anos agrícolas 2006/7, 2007/8 e 2008/9, em três locais da região Serrana do Estado do Espírito Santo. O delineamento experimental utilizado foi de blocos ao acaso com três repetições e 15 plantas por parcela. Pela análise de deviance em relação à produtividade (t ha^-1), apenas os efeitos de genótipos e da interação tripla genótipos x locais x anos foram significativas, evidenciando a presença de variabilidade genética entre as cultivares e a inconsistência da posição entre os genótipos para as combinações de anos e locais. Considerando s seleção para produtividade, adaptabilidade e estabilidade as duas cultivares que se destacaram foram 'Camarosa' e 'Aromas', tendo obtido valores superiores em média de 22% sobre a média geral das cultivares (24.55 t ha^-1) nos ambientes avaliados sob cultivo protegido TB.

Palavras-chave:
Fragaria x ananassa Duch.; REML/BLUP; método MHPRVG; interação genótipos x ambientes

Introduction

Strawberry (Fragaria x ananassa Duch.) is an accidental hybrid of two native species, F. chiloensis and F. virginiana, and is one of the most appreciated fruits in the world due to its organoleptic properties and health benefits (BOMBARELY et al., 2010;KURAS; KORBIN, 2010). In 2010, the world strawberry production reached 4.17 million tons, with the United States of America as the largest producer, accounting for 30% of the global production. In Brazil, the annual production is approximately of 100 thousand tons and is distributed in temperate and subtropical regions where it is used both for fresh and industrial consumption (RADIN et al., 2011).

In Espírito Santo State, strawberry cultivation has increased over the years and has played an important socio-economic role, mainly due to the increased income of small farms and the establishment of workers in rural areas because much manual labor is required during the harvest season. Approximately 7200 tons are produced per year, involving more than 240 hectares and generating approximately 2500 direct jobs (INCAPER, 2012).

Espírito Santo producers have gradually been using protected cultivation to increase production. This management system has a number of advantages, particularly the protection of the crop against abiotic (winds, hail, rain, frost and low temperatures) and biotic (pests and diseases) factors (ANTUNES et al., 2007;CALVETE et al., 2008;WITTER et al., 2012). However, there are few studies on the performance of strawberry cultivars in different locations and years under this management system.

Therefore, a study on the genotype x environment interaction (GE) plays an important role in determining the adaptation and fitness of genotypes to the physical environment. GE interaction can be classified as simple, such as when it results only from genetic variability in a genotype, and complex, such as when there is no correlation between the measures of the same genotype in different environments. The latter indicates inconsistency of genotype superiority under environmental variations, which impedes the recommendation of cultivars with wide adaptability (CRUZ et al., 2004).

To minimize the effects of GE interactions and to achieve a higher performance predictability, the identification of the most stable genotypes and cultivars adapted to the specific conditions of each environment is necessary. Several methods have been described for examining stability and adaptability, as based on the analysis of variance, linear regression, nonlinear regression, multivariate analysis and nonparametric statistics (BASTOS et al., 2007;PENA et al., 2012).

Within the context of mixed models, the method of the harmonic mean of the relative performance of genetically predicted values proposed byResende (2004) has shown good results for studies on stability and adaptability (BASTOS et al., 2007;CARBONELL et al., 2007;SILVA et al., 2011;SILVA et al., 2012). The advantages of this procedure include the following: i) allowing the simultaneous selection of yield, adaptability and stability; ii) considering the genotypic effects as random and, therefore, providing estimates of stability and adaptability for predicted and non-predicted genotypic values; iii) allowing the use of unbalanced data, non-orthogonal designs and heterogeneity of variance; iv) allowing the consideration of correlated errors within the locations in addition to stability and adaptability in the selection of individuals within progeny; v) providing genetic values already discounted (penalized) with regard to instability and allowing the application of any number of environments; and vi) generating results at the unit or scale of the evaluated trait, which can be directly interpreted as genetic values (RESENDE, 2002a).

Therefore, the present study aimed to estimate the adaptability and stability of strawberry cultivars under protected cultivation using the method of the harmonic mean of the relative performance of genetically predicted values.

Material and methods

The experiments were conducted during the agricultural years of 2006/7, 2007/8 and 2008/9 in three locations in the mountainous region of Espírito Santo State (municipalities Domingos Martins, Iúna and Muniz Freire). The elevations of the experimental sites in the municipalities were 700 m in Iúna, 950 m in Domingos Martins and 1100 m in Muniz Freire. The planting for the experiments was performed in late May, and the experiment was conducted until early January, with the exception of the experiment in 2007/8, which was conducted until March 2008. In all years, the harvests were performed twice a week, from the early harvest (usually in August) to the end of the cycle.

The experiment was arranged in a randomized block design, with 3 replications and 15 plants per plot. The study involved seven strawberry cultivars ('Dover', 'Camino Real', 'Ventana', 'Camarosa', 'Seascape', 'Diamante' and 'Aromas'). Plant spacing was 0.4 x 0.4 m, arranged in three rows on 30-cm high and 14-m long beds covered with black mulching. The environmental conditions for cultivation consisted of covering the plants with white plastic film, which was suspended over arcs of galvanized iron, at a height of approximately 1.0 m.

The values of productivity (ton. ha^-1) in the different environments were subjected to statistical analyses using the following linear model:

where:

Y_ijkn is the observation value relating to treatment i in repetition j in year k within location n; μ is the general average; g_i is the effect of genotype i; b_jkn is the effect of block j within year k within location n; a_k is the effect of the planting year k; l_n is the effect of the location n; ga_ik is the effect of the genotype x planting year interaction; gl_kn is the effect of the genotype x location interaction; gal_ikn is the effect of the genotype x year x location interaction; and є is the random error or residue.

In matrix form, the corresponding model is as follows:

where:

Y is the vector of the observations; b is the vector of the effects of the block-local-year combinations (fixed effects) plus the overall average; g is the vector of the genotypic effects, g (assumed as random); ga is the vector of random effects of the genotype x year interaction; gl is the vector of random effects of the genotype x location interaction; gla is the vector of random effects of the triple interaction genotype x location x year; and є is the vector of errors (random). X, Z, Q, T and U refer to the incidence matrices for these effects, respectively. The effects of the year, local and block-local-years are grouped into effect b, as they are environmental effects (RESENDE, 2007).

The values of the harmonic mean of the genotypic values (MHVG) to evaluate stability, the relative performance of the genotypic values (PRVG) for adaptability and the harmonic mean of the relative performance of the genotypic values (MHPRVG) for stability, adaptability and productivity, simultaneously, for all genotypes, were obtained according to the expressions:

where:

l is the number of locations, VG is the genotypic value, and j is the genotype.

The significance of the effects of the model was estimated by an analysis of deviance, as recommended byResende (2002a). The deviances were obtained through analyses with and without the effects of g, ga, gl and gla. The deviance without the mentioned effect was then subtracted from each full-model deviance and was compared to the value of chi-square, with one degree of freedom, at 1 and 5% probability. The Selegen software system (RESENDE, 2002b) was used for the analyses.

Results and discussion

The analysis of deviance in relation to productivity (ton. ha^-1) showed that the effects of genotype and the triple interaction genotype x location x year were statistically significant (p ≤ 0.05 and p ≤ 0.01, respectively). In contrast, no significance was observed for the genotype x location and genotype x year interactions (Table 1). These results demonstrate the presence of genetic variability among the strawberry cultivars tested and that the genotype x location and genotype x year interactions did not cause changes in the performance of the genotypes in the different locations and years. However, for the triple interaction genotype x location x year, an inconsistency in the position among the genotypes was verified for all the combinations of year and location. According toAllard and Bradshaw (1964), this type of interaction is caused by unpredictable environmental variations, such as precipitation, temperature, relative humidity and the occurrence of pests and diseases.

Table 1:
Deviance, variance components and coefficients of determination related to the global joint analysis involving seven strawberry cultivars evaluated in nine environments.

• ⁺Deviance of the adjusted model without the corresponding effects. ^* and ^**Significant by the chi-square test at 5% (3.84) and 1% (6.63), respectively

The coefficients of determination (for the genotype x location interaction,for the genotype x year interaction, andfor the genotype x location x year interaction) indicate how much each component contributes to the total phenotypic variance. Thus, the genotype x location, genotype x year, and genotype x location x year interactions contributed 0.24, 8.60 and 33.70%, respectively, and demonstrated the importance of the triple interaction in the total phenotypic variance. The heritability estimate of the broad-sense individual plotsfor productivity was 30.60%. This result is explained by the fact that productivity is a polygenic character and is therefore greatly influenced by the environment.

The genotypic correlation of the cultivars across the locations and years was 0.4183, showing a moderate level of complex interaction, which indicates that the ranking of the cultivars across environments (locations and years) will not necessarily be the same, i.e., a cultivar rated with excellent productivity in environment 1 may not be excellent in environment 2 or in any other environment.

Table 2shows the cultivars with the respective genotypic values (µ+g+ge) for each cultivation location in the different agricultural harvests. The Dover cultivar presented the lowest genotypic values for most environments evaluated, whereas cultivar Camarosa presented the best genotypic values for most environments and occupied the first or second position in eight of the nine environments assessed.Resende et al. (2010) evaluated the effect of three cultivation systems (high tunnel, low tunnel and field) on strawberry productivity in Guarapuava-PR and reported that the cultivars Camorosa and Oso Grande achieved the highest productivity rates (56.66 and 53.76 ton. ha^-1, respectively) and differed statistically from cultivars Sweet Charlie (45.70 ton. ha^-1) and Dover (36.98 ton. ha^-1) under low tunnel-protected cultivation.

According to the expected results (in ton. ha^-1) for stability (MHVG), adaptability (PRVG), and stability and adaptability (MHPRVG) for the set of locations, either harming or favoring cultivars, according to their performance, it was observed that cultivar Camorosa was superior, followed by cultivars Aromas, Diamante, Camino Real, Seascape, Ventana and Dover (Table 3). Therefore, the two best cultivars (Camarosa and Aromas) by the MHPRVG criterion were, on average, 22% superior to the overall average of the nine environments (24.55 kg ha^-1). The MHPRVG method has the advantage of providing results for the scale of trait measurement, which can be directly interpreted as predicted genetic values for productivity, stability and adaptability (BASTOS et al., 2007;CARBONELL et al., 2007;VERARDI et al., 2009).

Silva et al. (2011) compared the methods of adaptability and stability and verified that the REML/BLUP-MHPRVG analysis stood out from the others (AMMI and GGE biplot). Although the method does not provide information about the grouping of similar sites, the results are according to the genotypic values when considering the parameters of adaptability and stability.

Table 2:
Estimates of the genotypic values (µ + g + ge) related to the yield trait of seven strawberry cultivars evaluated in nine environments.

It is important to mention that the three cultivars with the highest values for MHVG, PRVG and MHPRVG have excellent performance for firmness and a good balance between the total soluble solids and acidity, such that they will be widely accepted by both consumers and industry. Within this context, such information on the stability and adaptability of strawberry cultivars under protected cultivation (low tunnel) are important because they allow recommendations for cultivation and the verification of the performance of these cultivars in different strawberry-producing regions in the state of Espírito Santo.

Table 3:
Stability of genotypic values (MHVG), adaptability of genotypic values (PRVG), stability and adaptability of genotypic values (MHPRVG) for the yield trait in seven strawberry cultivars evaluated in nine environments.

Conclusion

Based on selection for yield, adaptability and stability, the Camarosa and Aromas cultivars may be indicated for a low tunnel management system.

Acknowledgements

The author thanks FAPES (Fundação de Amparo à Pesquisa do Espírito Santo) for funding the project and Agricola Llahuen and Bioagro for donating a portion of the seedlings.

References

• ALLARD, R. W.; BRADSHAW, A. D. Implications of genotype environment interaction applied plant breeding. Crop Science, v. 4, n. 5, p. 503-508, 1964.
• ANTUNES, O. T.; CALVETE, E. O.; ROCHA, H. C.; NIENOW, A. A.; CECCHETTI, D.; RIVA, E.; MARAN, R. E. Produção de cultivares de morangueiro polinizadas pela abelha jataí em ambiente protegido. Horticultura Brasileira, v. 25, n. 1, p. 94-99, 2007.
• BASTOS, I. T.; BARBOSA, M. H. P.; RESENDE, M. D. V.; PETERNELLI, L. A.; SILVEIRA, L. C. I.; DONDA, L. R.; FORTUNATO, A. A.; COSTA, P. M. A.; FIGUEIREDO, I. C. R. Avaliação da interação genótipo x ambiente em cana-de-açucar via modelos mistos. Pesquisa Agropecuária Tropical, v. 37, n. 4, p. 195-203, 2007.
• BOMBARELY, A.; MERCHANTE, C.; CSUKASI, F.; CRUZ-RUS, E.; CABALLERO, J. L.; ESCOBAR N. M.; PORTALES, R. B.; BOTELLA, M. A.; BLANCO, J. M.; SEVILLA, J. F. S.; VALPUESTA, V. Generation and analysis of ESTs from strawberry (Fragaria x ananassa) fruits and evaluation of their utility in genetic and molecular studies. BMC Genomics, v. 11, n. 1, p. 1-17, 2010.
• CALVETE, E. O.; MARIANI, F.; WESP, C. L.; NIENOW, A. A.; CASTILHOS, T.; CECCHETTI, D. Fenologia, produção e teor de antocianinas de cultivares de morangueiro em ambiente protegido. Revista Brasileira de Fruticultura, v. 30, n. 2, p. 396-401, 2008.
• CARBONELL, S. A. M.; CHIORATO, A. F.; RESENDE, M. D. V.; DIAS, L. A. S.; BERALDO, A. L. A.; PERINA, E. F. Estabilidade de cultivares e linhagens de feijoeiro em diferentes ambientes no estado de São Paulo. Bragantia, v. 66, n. 2, p. 193-201, 2007.
• CRUZ, C. D.; REGAZZI, A. J.; CARNEIRO, P. C. S. Modelos biométricos aplicados ao melhoramento genético. 3. ed. Viçosa: UFV, 2004, v. 1, 2004.
• INCAPER-Instituto Capixaba de Pesquisa, Assistência Técnica e Extensão Rural. Available from: <Available from:http://www. incaper.es.gov.br >. Access on: Aug. 2, 2012.
  » http://www. incaper.es.gov.br
• KURAS, A.; KORBIN, M. PCo-A analysis of strawberry germplasm used in European breeding programs, based on evaluation of DNA polymorphism of investigated plants. Journal of Fruit and Ornamental Plant Research, v. 18, n. 2, p. 7-16, 2010.
• PENA, G. F.; AMARAL JÚNIOR, A. T.; GONÇALVES, L. S. A.; CANDIDO, L. S.; VITTORAZZI, C.; RIBEIRO, R. M.; FREITAS JÚNIOR, S. P. Stability and adaptability of popcorn genotypes in the state of Rio de Janeiro, Brazil. Genetics and Molecular Research, v. 11, n. 4, p. 1320-127, 2012.
• RADIN, B.; LISBOA, B. B.; WITTER, S.; BARNI, V.; REISSER JUNIOR, C.; MATZENAUER, R.; FERMINO, M. H. Desempenho de quatro cultivares de morangueiro em duas regiões ecoclimáticas do Rio Grande do Sul. Horticultura Brasileira v. 29, n. 3, p. 287-291, 2011.
• RESENDE, M. D. V. Genética biométrica e estatística no melhoramento de plantas perenes. Brasília: Embrapa Informação Tecnológica, 2002a.
• RESENDE, M. D. V. Software Selegen-REML/BLUP. Colombo: Embrapa Florestas, 2002b.
• RESENDE, M. D. V. Métodos estatísticos ótimos na análise de experimentos de campo., Colombo: Embrapa Florestas 2004.
• RESENDE, M. D. V. Matemática e estatística na análise de experimentos e no melhoramento genético., Colombo: Embrapa Florestas 2007.
• RESENDE, J. T. V.; MORALES, R. G. F.; FARIA, M. V.; RISSINI, A. L. L.; CAMARGO, C. K. Produtividade e teor de sólidos solúveis de frutos de cultivares de morangueiro em ambiente protegido. Horticultura Brasileira, v. 28, n. 2, p. 185-189, 2010.
• SILVA, G. O.; CARVALHO, A. D. F.; VIEIRA, J. V.; BENIN, G. Verificação da adaptabilidade e estabilidade de populações de cenoura pelos métodos AMMI, GGE biplot e REML/BLUP. Bragantia, v. 70, n. 3, p. 494-501, 2011.
• SILVA, G. O.; CARVALHO, A. D. F.; VIEIRA, J. V.; FRITSCHE-NETO, R. Adaptabilidade e estabilidade de populações de cenoura. Horticultura Brasileira, v. 30, n. 1, p. 80-83, 2012.
• VERARDI, C. K.; RESENDE, M.D.V.; COSTA, R. B.; GONÇALVES, P. S. Adaptabilidade e estabilidade da produção de borracha e seleção em progênies de seringueira. Pesquisa Agropecuária Brasileira, v. 44, n. 1, p. 1277-1282, 2009.
• WITTER, S.; RADIN, B.; LISBOA, B. B.; TEIXEIRA, J. S. G.; BLOCHTEIN, B.; IMPERATRIZ-FONSECA, V. C. Desempenho de cultivares de morango submetidas a diferentes tipos de polinização em cultivo protegido. Pesquisa Agropecuária Brasileira , v. 47, n. 1, p. 58-65, 2012.

Publication Dates

• Publication in this collection
  Dec 2015

History

• Received
  15 Aug 2012
• Accepted
  30 Oct 2012