Acessibilidade / Reportar erro

PHENOTYPIC MODELS OF COMPETITION FOR Pinus taeda L GENETIC PARAMETERS ESTIMATION

MODELOS FENOTÍPICOS DE CONCORRÊNCIA PARA Pinus taeda L ESTIMAÇÃO DE PARÂMETROS GENÉTICOS

ABSTRACT

The objective of this study was to evaluate the efficiency of the phenotypic models of competition, through spatial analysis in the genetic evaluation of Pinus taeda L progenies. For this, four competition covariates were used to adjust the phenotypic values ​​in a P. taeda progeny test installed in four different locations in the state of Santa Catarina. The test was implemented in randomized block design, with seven repetitions, linear plots containing six plants per plot in 2.5 m x 2.0 m spacing. The test installed in sites A, B, and D present 63 families and site C 53 families. At nine years old, the diameter at the breast height was measured for all individuals. The presence or absence of competition was based on the residual autocorrelation coefficients, which had its significance tested by the Durbin-Watson test. In general, the use of covariates corrected the competition effect. The variances among and within plots, as well as the residual variation coefficient, were reduced. The classification by the genetic effect of the individuals in the progeny test was extremely altered for this data set with and without the use of covariates for sites A and D, as well as the genotype x environment interaction. The use of these two tools is of great importance in the analysis of data in P. taeda progeny tests, since the effects of competition can lead to mistakes in the selection of individuals and in the definition of improvement zones.

Keywords:
Spatial analysis; Covariate; REML/BLUP; Components of variance

RESUMO

O objetivo foi avaliar a eficiência dos modelos fenotípicos de competição, via análise espacial, na avaliação genética de testes de progênies de Pinus taeda L. Para isso foram utilizadas quatro covariáveis de competição para ajuste dos valores fenotípicos em um teste de progênies de P. taeda plantado em quatro diferentes locais no estado de Santa Catarina. O teste foi implantado em delineamento de blocos casualizados, com sete repetições, parcelas lineares de seis plantas e espaçamento de 2,5 m x 2,0 m. O teste instalado nos locais A, B e D possuem 63 famílias e no local C 53 famílias. Aos nove anos de idade foi efetuada a mensuração da variável diâmetro à altura do peito de todos os indivíduos. A presença ou não da competição foi baseada nos coeficientes de autocorrelação residual, que teve sua significância testada pelo teste de Durbin-Watson. No geral, o uso de covariáveis corrigiu o efeito da competição e as variâncias entre e dentro de parcela e o coeficiente de variação residual foram reduzidos. A análise espacial foi apropriada para validar a eficiência das covariáveis. A classificação pelo efeito genético dos indivíduos no teste de progênies foi extremamente alterada para o conjunto de dados com e sem o uso de covariáveis para os locais A e D, assim como a interação genótipo x ambiente. O uso dessas duas ferramentas é de suma importância na análise de dados em testes de progênies de P. taeda, pois os efeitos de competição podem conduzir a erros na seleção de indivíduos e na definição de zonas de melhoramento.

Palavras chave:
Análise espacial; Covariável; REML/BLUP; Componentes de variância

INTRODUCTION

Biological variation is the basis for the breeder’s work genetic improvement programs of forest species, the most widely used way to distinguish between environmental and genetic variability is progeny testing. This allows to estimate genetic and phenotypic parameters and to predict the genetic value of a specific set of individuals. These tests are commonly established in small linear plots, ranging from four to nine plants (STOEHR et al., 2010). This way, there is an increase in statistical efficiency, since they reduce the variance within the block, but increase the interaction between neighboring individuals due to the competition (LEONARDECZ-NETO et al., 2003LEONARDECZ- NETO, E.; VENCOVSKY, R.; SEBBENN, A. M. Ajuste para a competição entre plantas em teste de progênies e procedências de essências florestais. Scientia Forestalis, n. 63, p. 136-149, 2003.).

Competition is caused by genetic and environmental factors and is due to genotypes and the spatial arrangement of neighboring trees (CAPPA et al., 2016CAPPA, E. P., STOEHR, M. U., XIE, C. Y., & YANCHUK, A. D. Identification and joint modeling of competition effects and environmental heterogeneity in three Douglas-fir (Pseudotsuga menziesii var. menziesii) trials. Tree Genetics & Genomes, v. 12, n. 6, p. 102, 2016.). In forest field trials, it is mainly caused by the lack of space induced by high intraspecific density (LEONARDECZ-NETO et al., 2003LEONARDECZ- NETO, E.; VENCOVSKY, R.; SEBBENN, A. M. Ajuste para a competição entre plantas em teste de progênies e procedências de essências florestais. Scientia Forestalis, n. 63, p. 136-149, 2003.). In which the growth of the individuals of the plot is correlated to its neighbors, due to the common micro environment (PAVAN et al., 2012PAVAN, B. E.; DE PAULA, R. C.; PERECIN, D.; CANDIDO, L. S.; SCARPINATI, E. A. Efeito competicional em testes de progênies de eucalipto. Ciência Florestal, v. 22, n. 3, p. 433-443, 2012.). In linear plots, this competition can occur within the plots (inter-genotypic competition) or between the experimental plots (intra-genotypic competition).

However, in the models used for analysis of progeny tests, it is not common to consider the effects of competition, as there are no established routines; and the definition of specific methods to correct this type of interference (PAVAN et al., 2012PAVAN, B. E.; DE PAULA, R. C.; PERECIN, D.; CANDIDO, L. S.; SCARPINATI, E. A. Efeito competicional em testes de progênies de eucalipto. Ciência Florestal, v. 22, n. 3, p. 433-443, 2012.), this implies in the reduction of the ability to identify the genetic variation and the precision of the genetic value prediction (YE and JAYAWICKRAMA, 2008YE, T. Z.; JAYAWICKRAMA, K. J. S. Efficiency of using spatial analysis in first-generation coastal Douglas-fir progeny tests in the US Pacific Northwest. Tree Genetics & Genomes , v. 4, n. 4, p. 677-692, 2008.).

This interference may occur in the favoring of certain plants due to the lack of competition as there may be faults in the experiment or by more aggressive genotypes that tend to have their performances overestimated once they compete with more sensitive genotypes (RESENDE et al., 2005RESENDE, M.D. V.; STRINGER, J.; CULLIS, B.; THOMPSON, R. joint modelling of competition and spatial variability in forest field trials. Revista Matemática e Estatística, v.23, n.2, p.7-22, 2005.). Progeny selection is also influenced, as the genetic effects may be underestimated or overestimated (RESENDE et al., 2016).

Overall, two competition models can be used: the phenotypic model, when competition or interference does not have genetic control, thus, treating the neighbors phenotypic value as a covariate; and the genotypic model, when there is genetic control in the competition, relating the genetic effects of the neighbors to the residual value of the plant (RESENDE et al., 2005RESENDE, M.D. V.; STRINGER, J.; CULLIS, B.; THOMPSON, R. joint modelling of competition and spatial variability in forest field trials. Revista Matemática e Estatística, v.23, n.2, p.7-22, 2005.).

Genetic models of competition can be seen in Resende et al. (2016RESENDE, R. T.; MARCATTI, G. E.; PINTO, D. S.; TAKAHASHI, E. K.; CRUZ, C.D.; RESENDE, M. D. V. Intra-genotypic competition of Eucalyptus clones generated by environmental heterogeneity can optimize productivity in forest stands. Forest Ecology and Management, v. 380, p. 50-58, 2016.) testing intra-genotypic competition in Eucalyptus clones, Costa et al. (2013). Using simulated modeling data and Cappa and Cantet (2008CAPPA, E. P.; CANTET, R. J. Direct and competition additive effects in tree breeding: Bayesian estimation from an individual tree mixed model. Silvae Genetica, v. 57, n. 2, p. 45-55, 2008.) in P. taeda tests at 13 years of age using models that incorporate competition for both, genetic and non-genetic levels. Phenotypic models were studied by Leonardecz-Neto et al. (2003), who used Hegyi competition index for DBH data in several forest species and Pavan et al. (2012PAVAN, B. E.; DE PAULA, R. C.; PERECIN, D.; CANDIDO, L. S.; SCARPINATI, E. A. Efeito competicional em testes de progênies de eucalipto. Ciência Florestal, v. 22, n. 3, p. 433-443, 2012.), using besides this index other six covariates in Eucalyptus progenies tests for volume. The DBH and volume variables are more influenced by spacing and more susceptible to competition effects than the height (YE and JAYAWICKRAMA, 2008YE, T. Z.; JAYAWICKRAMA, K. J. S. Efficiency of using spatial analysis in first-generation coastal Douglas-fir progeny tests in the US Pacific Northwest. Tree Genetics & Genomes , v. 4, n. 4, p. 677-692, 2008.; STOEHR et al., 2010). Since DBH is a readily available variable through direct measurement and because volume estimates are dependent on this variable, it is preferable to adjust phenotypic competition models for DBH.

Thus, with the use of competition models, it is expected to isolate in a more efficient way the variation that occurs due to the environment, causing the analysis to express more consistently the experimental reality (LEONARDECZ-NETO et al., 2003LEONARDECZ- NETO, E.; VENCOVSKY, R.; SEBBENN, A. M. Ajuste para a competição entre plantas em teste de progênies e procedências de essências florestais. Scientia Forestalis, n. 63, p. 136-149, 2003.). However, knowing which covariant and how it fits to faithfully represent, the biological effect, is often unknown (PAVAN et al., 2012PAVAN, B. E.; DE PAULA, R. C.; PERECIN, D.; CANDIDO, L. S.; SCARPINATI, E. A. Efeito competicional em testes de progênies de eucalipto. Ciência Florestal, v. 22, n. 3, p. 433-443, 2012.). One way to diagnose the presence of competition effects in field trials is from the residual correlation coefficient in which competition effects between neighboring plants cause negative autocorrelation and spatial heterogeneity causes positive correlations (RESENDE et al., 2005RESENDE, M.D. V.; STRINGER, J.; CULLIS, B.; THOMPSON, R. joint modelling of competition and spatial variability in forest field trials. Revista Matemática e Estatística, v.23, n.2, p.7-22, 2005.; CAPPA et al., 2016CAPPA, E. P., STOEHR, M. U., XIE, C. Y., & YANCHUK, A. D. Identification and joint modeling of competition effects and environmental heterogeneity in three Douglas-fir (Pseudotsuga menziesii var. menziesii) trials. Tree Genetics & Genomes, v. 12, n. 6, p. 102, 2016.). Thus, the objective of the present work is to evaluate the efficiency of the phenotypic models of competition using four distinct covariates, through spatial analysis, in the genetic evaluation of a P. taeda progeny test.

MATERIAL AND METHODS

Experimental area and genetic material

The experimental area is located in two cities: Rio Negrinho and Lages in a plateau region of Santa Catarina state. The climate is classified as Cfb according to Köppen, subtropical humid with mild summer, without dry season (ALVARES et al., 2013ALVARES, C. A.; STAPE, J. L.; SENTELHAS, P. C.; de MORAES, G.; LEONARDO, J.; SPAROVEK, G. Köppen’s climate classification map for Brazil. Meteorologische Zeitschrift , v. 22, n. 6, p. 711-728, 2013.). The regions present frost incidence in the colder months. Table 1 shows the other characteristics and location of the experiment.

TABLE 1
Location and characterization of the experimental area of the P. taeda progeny test.

The P. taeda L. genetic material comes from a 1st generation Clonal Orchard established in 1994, with material introduced from Georgia and Zimbabwe provenance, and from a seed production area with material from Zimbabwe provenance. The open pollinated progeny test was conducted in four different environments (Sites A, B, C and D) and the installation occurred in 2006 in randomized block design with seven repetitions, six plant linear plots, in 2.5 m x 2.0 m spacing and double border rows. The test installed in sites A, B, and D has 63 families and site C 53 families. The 53 families present at site C are common to all sites. At the age of nine, a direct measurement of the variable Diameter at Breast height (DBH, cm) of all individuals at all progeny test sites was performed.

Estimation of the spatial autocorrelation function

The autocorrelation coefficients between residues following the direction of the rows or between columns (ρL ), following the direction of the columns or between rows (ρc ) and the value of the d statistic for the Durbin-Watson residue independence test were obtained using SELEGEN-REML / BLUP® software statistical model number 113 (RESENDE, 2007aRESENDE, M. D. V. Selegen-Reml/Blup: Sistema Estatístico e Seleção Genética Computadorizada via Modelos Lineares Mistos. Colombo: Embrapa Florestas, 2007a. 361p.). The value of the statistic d is defined as described in equation 1, where L is the order of the plot positioning in the experiment and n refers to the residues as ordered in the space. The residues (r i ) were obtained via model 1 SELEGEN - REML/BLUP® software (equation 3) considering zero in the analysis. The relation between d and ρ is given by equation 2.

d = i = 1 n ( r i r i 1 ) 2 / L = 1 n r i 2 (1)

d = 2 ( 1 ρ ) (2)

y = X r + Z a + W p + ε (3)

The calculated value of d was compared with theoretical values with (n-1) degrees of freedom, where n is the number of data. The calculated value of dl (lower d) and du (upper d) are given in tables associated to certain significance levels. For hypothesis Ha: ρ>0 the decision rule about d is given by: if ddL, it rejects H0 and autocorrelation is positive; If ddu, it accepts H0 and the autocorrelation is null and if dl<d<du the test is inconclusive. For hypothesis Ha: ρ < 0 the decision rule is equivalent, but using (4-d) instead of d (RESENDE, 2007bRESENDE, M.D. V. Matemática e estatística na análise de experimentos e no melhoramento genético. Colombo: Embrapa Florestas, 2007b. 362p.).

Estimation of genetic parameters

Genetic analysis of progeny tests were performed using SELEGEN - REML / BLUP® software. For the selection of individuals and progenies from in open-pollinated parents, in the design used and performing the selection one site at a time, the statistical model 1 (equation 3) was used; and for the selection in the four sites (joint analysis) model 4 (equation 4) was used. In this model, y is the vector of data, r is the vector of the effects of repetition (assumed as fixed) added to the general mean; a is the vector of the individual additive genetic effects (assumed as random); p is the vector of plot effects (assumed as random); i is vector of the effects of genotype x environment interaction (random) and ε is the vector of (random) errors or residues. Capital letters represent the incidence matrices for these effects (RESENDE, 2007aRESENDE, M. D. V. Selegen-Reml/Blup: Sistema Estatístico e Seleção Genética Computadorizada via Modelos Lineares Mistos. Colombo: Embrapa Florestas, 2007a. 361p.). For the estimates of the genetic parameters the zeros were not considered in the analysis.

y = X r + Z a + W p + T i + ε (4)

Analysis of covariance

The analysis of covariance of the progeny tests were performed with SELEGEN - REML / BLUP® software using statistical model 131 (equation 5). In this, y is the data vector; r is the vector of the repetition effects (assumed as fixed values) added to the general mean; a is the vector of the individual additive genetic effects (assumed as random), p is the vector of the plot effects (assumed as random), ε is the vector of (random) errors or residues. The coefficient β refers to the regression associated with the covariate Cov Capital letters represent the incidence matrices for these effects (RESENDE, 2007aRESENDE, M. D. V. Selegen-Reml/Blup: Sistema Estatístico e Seleção Genética Computadorizada via Modelos Lineares Mistos. Colombo: Embrapa Florestas, 2007a. 361p.). For the analysis of covariance, zeros were not considered.

y = X r + β C o v + Z a + W p + ε (5)

The assumptions for analysis of covariance are: (I) the covariate values ​​are fixed, measured without error and independent of the treatments. (II) The regression between the covariate and the variable of interest, after removal of block and treatment differences, is linear and independent of blocks and treatments. (III) The residue is normally and independently distributed, with zero mean and common variance (RESENDE, 2002RESENDE, M. D. V. Genética biométrica e estatística no melhoramento de plantas perenes. Brasília: Embrapa. Informação Tecnológica, 2002. 975p.; PALLANT, 2010PALLANT, J. SPSS survival manual: A step by step guide to data analysis using SPSS. Maidenhead, 2010. 354 p. ; RAMALHO et al., 2012RAMALHO, M. A. P.; FERREIRA, D. F.; OLIVEIRA, A. C. Experimentação em genética e melhoramento de plantas 3ª Ed. Lavras: Editora UFLA, 2012. 305p.).

Three covariates were tested based on the performance of neighboring trees and one covariate based on the number of neighboring plants.

Hegyi Competition Index (IC)

Competition Mean (IC) through the distances and performance of the eight potential neighbors of the target plant i (LEONARDECZ-NETO et al., 2003LEONARDECZ- NETO, E.; VENCOVSKY, R.; SEBBENN, A. M. Ajuste para a competição entre plantas em teste de progênies e procedências de essências florestais. Scientia Forestalis, n. 63, p. 136-149, 2003.). The expression for the calculation is given by equation 6, where IC i is the competition index of the target plant i, Y i is the observed value of the target plant i; Y j is the observed value of competing plant j and D ij is the distance between plants i and j.

I C i = ( j = 1 8 Y j / Y i ) / D i j (6)

Self-competition mean (MAT)

Arithmetic mean of the silvicultural characteristic of the auto competitive trees (PAVAN et al., 2012PAVAN, B. E.; DE PAULA, R. C.; PERECIN, D.; CANDIDO, L. S.; SCARPINATI, E. A. Efeito competicional em testes de progênies de eucalipto. Ciência Florestal, v. 22, n. 3, p. 433-443, 2012.), belonging to the same genetic material (equation 7), where Y j is the observed value of the self-competing plant j and n the number of auto competitive plants.

M A T = ( j = 1 2 Y j ) / n (7)

Alo competition Mean (MAL)

Arithmetic mean of the silvicultural characteristic of the alo competitive trees (PAVAN et al., 2012PAVAN, B. E.; DE PAULA, R. C.; PERECIN, D.; CANDIDO, L. S.; SCARPINATI, E. A. Efeito competicional em testes de progênies de eucalipto. Ciência Florestal, v. 22, n. 3, p. 433-443, 2012.), belonging to different genetic material (equation 8), where Y j is the observed value of the alo competitive plant j and n the number of alo competitive plants.

M A L = ( j = 1 7 Y j ) / n (8)

Number of neighboring plants (NV)

Algebraic sum of the number of surviving neighboring plants (equation 9), where n is the surviving neighboring plant (RESENDE, 2002RESENDE, M. D. V. Genética biométrica e estatística no melhoramento de plantas perenes. Brasília: Embrapa. Informação Tecnológica, 2002. 975p.).

N V = j = 1 8 ( n ) (9)

Adjustment of phenotypic values

This is based on the covariance value of each individual (x ij ), the mean value of the covariate () and the residual regression coefficient (β) between the covariate and the DBH data. The phenotypic values for each surviving individual were adjusted according to equation 10 for the covariates IC, MAT and MAL and equation 11 (RESENDE, 2002RESENDE, M. D. V. Genética biométrica e estatística no melhoramento de plantas perenes. Brasília: Embrapa. Informação Tecnológica, 2002. 975p.) for covariate NV, where Y ij is the phenotypic value measured for each individual and Y ijc the adjusted phenotypic value for each individual.

Y i j c = Y i j β ( x i j x ¯ ) (10)

Y i j c = Y i j β ( x i j 8 ) (11)

RESULTS AND DISCUSSION

The major cause of the interference found in the four progeny tests was due to experiment failures, which caused absence of competition and the favoring of isolated individuals. At test site A the favoring was intensified by the presence of three families that had much inferior performance to the others throughout the test. Mortality had genetic control for the tests at site A and C (TABLE 2). However, the competition effect was random in all tests, showing no tendency for either family.

TABLE 2
Deviance analysis for survival of P. taeda progenies at 9 years of age in Santa Catarina state.

Assumptions for analysis of covariates

The residual regression coefficients (β) for the covariates alo competition mean (MAL) and number of neighboring plants (NV) showed no interaction between covariate and DBH in any of the sites. The covariate self-competition mean (MAT) showed significant interaction between covariate x treatment for sites A, B and D. The same occurred for Hegyi competition index (IC) for site B. There was significance of the β coefficients for all covariates, and these were negative for the covariates tested (TABLE 3).

TABLE 3
Residual regression coefficient (β) between covariate and DBH data and mean value of covariate () in P. taeda progeny test in four test sites in the state of Santa Catarina.

As in Leonardecz-Neto et al. (2003) and Pavan et al. (2012PAVAN, B. E.; DE PAULA, R. C.; PERECIN, D.; CANDIDO, L. S.; SCARPINATI, E. A. Efeito competicional em testes de progênies de eucalipto. Ciência Florestal, v. 22, n. 3, p. 433-443, 2012.), for the covariates calculated based on the performance of neighboring trees, the assumptions for analysis of covariance are not satisfied in full, since they do not fully meet the first condition of the analysis of covariance. To minimize this implication, for covariates where residual regression coefficients that were not significant or showed no significant interaction with the treatments they were not used (TABLE 3). The same authors say that the use of indexes containing estimates and data-dependent errors should not be the only or the main concern for a reliable analysis. The proper choice of the site for the allocation of the experiments, the reduction of the variance of the experimental and blocks errors are as or more important than the error associated with the covariate.

Phenotypic models of competition for sites A and D

For DBH data without adjustment, the coefficient of autocorrelation between residues following the direction of the columns (ρc ) was negative, significant and the statistic value d was higher than two for sites A and D. The same trend happened for the direction of the lines (ρL ) (TABLE 4). Negative and significant values of ρc and ρL are indicative of the presence of competition in the progeny tests studied, since competition effects between neighboring plants cause negative autocorrelation and d statistic values greater than two. Spatial heterogeneity causes positive correlations (RESENDE et al., 2005RESENDE, M.D. V.; STRINGER, J.; CULLIS, B.; THOMPSON, R. joint modelling of competition and spatial variability in forest field trials. Revista Matemática e Estatística, v.23, n.2, p.7-22, 2005.; CAPPA et al., 2016CAPPA, E. P., STOEHR, M. U., XIE, C. Y., & YANCHUK, A. D. Identification and joint modeling of competition effects and environmental heterogeneity in three Douglas-fir (Pseudotsuga menziesii var. menziesii) trials. Tree Genetics & Genomes, v. 12, n. 6, p. 102, 2016.).

TABLE 4
Variance components and genetic parameters, autocorrelation coefficients between residues in the directions of the lines (ρl), columns (ρc) and value of the statistic d for rows and columns; with and without the use of covariates for DBH in a progeny test of P. taeda planted in two different sites in the state of Santa Catarina.

The sites showed the same trend in progeny behavior after the covariate analysis. The coefficients ρc and ρL were negative and significant for sites A and D when the data were adjusted via NV covariate. This covariant practically did not alter the components of variance and genetic parameters when compared with the original data (TABLE 4). The adjustment using the covariate MAL presented negative and significant ρc , increase in σ²e , reduction in σ²f and in parc . The value of this covariate considers only the competition among plots, thus not correcting intra-genotypic competition. The ρL indicated absence of inter-genotypic competition, indicating that this covariate corrected the data between plots for the tests of sites A and D, which can also be verified by the reduction in σ²parc (TABLE4).

The covariate IC resulted in ρc not significant and ρL significant and positive, indicating absence of competition between lines and environmental heterogeneity between columns, respectively. Ye and Jayawickrama (2008YE, T. Z.; JAYAWICKRAMA, K. J. S. Efficiency of using spatial analysis in first-generation coastal Douglas-fir progeny tests in the US Pacific Northwest. Tree Genetics & Genomes , v. 4, n. 4, p. 677-692, 2008.) obtained average autocorrelation coefficients of 0.75, indicating strong environmental variation, with environmental heterogeneity being the major source of interference in the progeny tests evaluated by the authors.

This environmental heterogeneity reflects changes in the variance components after the adjustment when compared to the original data, where there was a small increase in environmental variance between plots (σ²parc ) and a large reduction in individual phenotypic variance (σ²f ), resulting in an increase in the coefficient of determination of plot effects ( parc ) at both sites. The absence of correlation between the lines corresponds to an expressive reduction in the residual variance within the plot (σ²e ) (TABLE 4). The reduction of σ²parc and σ²e leads to an environmental improvement providing more reliable data for the selection process (PAVAN et al., 2012PAVAN, B. E.; DE PAULA, R. C.; PERECIN, D.; CANDIDO, L. S.; SCARPINATI, E. A. Efeito competicional em testes de progênies de eucalipto. Ciência Florestal, v. 22, n. 3, p. 433-443, 2012.). There was also a reduction in the additive genetic variance (σ2 a ), that is, the adjusted data presented lower amplitude between the best and worst genetic material (TABLE 4). As in the work of Pavan et al. (2012PAVAN, B. E.; DE PAULA, R. C.; PERECIN, D.; CANDIDO, L. S.; SCARPINATI, E. A. Efeito competicional em testes de progênies de eucalipto. Ciência Florestal, v. 22, n. 3, p. 433-443, 2012.), the IC was the covariate that mostly altered the components of variance at sites A, C and D, indicating great environmental influence on progeny tests due to competition between and within families.

Martinez et al. (2012MARTINEZ, D.T.; RESENDE, M. D. V.; COSTA, R. B. da; HIGA, A. R.; SANTOS, G. A. dos; FIER, I. S. N. Estudo da interação genótipo x ambiente em progênies de Pinus taeda por meio da análise de parâmetros genéticos. Floresta, v. 42, n. 3, p. 539 - 552, 2012.) obtained an individual heritability average (h2a ) of close to 0.3 and parc considered as low magnitude. As mentioned in the work of Sturion and Resende (2005STURION, J. A. e RESENDE, M. D. V. Eficiência do delineamento experimental e capacidade de teste no melhoramento genético da erva-mate (Ilex paraguariensis St. Hil). Pesquisa Floresta l Brasileira, n. 50, p. 3-10, 2005.) it is considered as low parc values lower than 0.10 for heritability values close to 0.30. The test at site A, the adjusted IC data, presented a high value of h 2 a (0.42) and c 2 parc of 0.15 and a coefficient of residual variation (CV e% ) of 6%, and site B showed moderate h 2 a (0.24) and parc of 0.23 and a residual coefficient of variation (CV e% ) of 5.5% (TABLE 4). Although parc is slightly above recommended, the improvement in CV e% and absence of competition effects justify the use of IC covariate in the selection for the two sites.

Phenotypic models of competition for sites B and C

For the DBH data without adjustment, ρc was negative, significant and statistic value d higher than two for sites B and C. The local B presented ρL null, that is, absence of correlation and the negative and significant C site (TABLE 5).

TABLE 5
Parameters estimates of the omnidirectional and cross semivariograms Vol = volume, Exp = Exponential.

For the tests installed at these sites, the NV covariate corrected the data for the competition effects, presenting non-significant ρc and ρL and the genetic parameters were similar to the data set parameters without adjustment. The covariate MAL exhibited ρL not significant and ρc negative and significant, and that according to what occurred in sites A and D, increase in σ²e and reduction in σ²parc . The same occurred with the use of this covariate at site C (TABLE 5).

The MAT covariate also corrected the data for the competition effects for site C, provided increase in σ²parc , reduction in σ²e , σ²f and a slight increase in parc . For this location, two covariates corrected the data for competition effects with similar variance components, so the choice of the covariate to be used in the selection can be made based on the one that offers greater individual heritability ( a ), in the data set adjusted by MAT (TABLE 5).

The CV e% remained the same or reduced in the tests in which the non-competition condition was met, showing that the data adjustment provided good experimental accuracy. As in the works of Leonardecz-Neto et al. (2003), Ye and Jayawickrama (2008YE, T. Z.; JAYAWICKRAMA, K. J. S. Efficiency of using spatial analysis in first-generation coastal Douglas-fir progeny tests in the US Pacific Northwest. Tree Genetics & Genomes , v. 4, n. 4, p. 677-692, 2008.) and Pavan et al. (2012PAVAN, B. E.; DE PAULA, R. C.; PERECIN, D.; CANDIDO, L. S.; SCARPINATI, E. A. Efeito competicional em testes de progênies de eucalipto. Ciência Florestal, v. 22, n. 3, p. 433-443, 2012.), in general, coefficients of experimental variation were reduced with the use of covariates.

The values of individual heritability ( a ), mean heritability of progenies ( mp ), and additive heritability within plot (h 2 ad ), were not significantly changed by analysis without adjustment or with adjustment by covariates (TABLE 5). Similar results were seen by Resende et al. (2005RESENDE, M.D. V.; STRINGER, J.; CULLIS, B.; THOMPSON, R. joint modelling of competition and spatial variability in forest field trials. Revista Matemática e Estatística, v.23, n.2, p.7-22, 2005.) testing several models of competition in experiment with Eucalyptus maculata. The author did not obtain reductions in the estimates of heritability between the traditional analysis and several models of competition tested, because in the dataset studied the competition was purely environmental.

Individual selection results

The genetic effects predicted for DBH using adjusted data suffered an intense reduction in site A and a slightly smaller reduction for site D, that is, gains on selection would be overestimated for these sites. For these sites, only eight individuals in the ranking of the twenty with the highest genetic values are present in both sets of data (TABLE 6). Some individuals which would be selected without data adjustment, mainly the closest to the gaps, were no longer selected after data adjustment. However, some individuals which were selected previously the adjustment of data and that were not close to the gaps, remained in the selection (FIGURE 1).

TABLE 6
Classification of the twenty best individuals by genetic effects for DBH in a progeny test of P. taeda planted at two sites in the state of Santa Catarina.

FIGURE 1
Detail of the progenies test installed on the site A which demonstrates the individuals that would be selected before the adjustment, after the adjustment and which remains selected after the data adjustment.

For site B there was a small reduction in the genetic effect and in site C the adjustment for the effects of competition increased the value of the predicted genetic effect. At site B only 5% of individuals are not the same for the set and unadjusted data sets and at site C 30% (TABLE 7). In this way, genetic materials have different responses to bioclimatic conditions, and may have different aptitudes, with some presenting or not increased growth, due to the lack of competition caused by the failure of neighboring plants (PAVAN et al., 2012PAVAN, B. E.; DE PAULA, R. C.; PERECIN, D.; CANDIDO, L. S.; SCARPINATI, E. A. Efeito competicional em testes de progênies de eucalipto. Ciência Florestal, v. 22, n. 3, p. 433-443, 2012.).

TABLE 7
Classification of the twenty best individuals by genetic effects for DBH in a progeny test of P. taeda planted in two sites in the state of Santa Catarina.

Correlation among sites

The variance components and genetic parameters for joint analysis, considering all four sites, were estimated for DBH data with and without adjustment (TABLE 8). The value of a for the unadjusted data set was higher than for the adjusted data (TABLE 8), since the estimation of the average individual heritability in the restricted sense in the environment is influenced by the genotype x environment interaction.

TABLE 8
Variance components and genetic parameters estimated by joint analysis for DBH in a progeny test of P. taeda planted in four sites in the state of Santa Catarina.

Genetic correlation between progeny performance in the four environments was high (0.80) for unadjusted DBH data, the same value found in P. taeda tests at four sites (distributed in the cities Otacílio Costa / SC And Telêmaco Borba / PR) in the research presented by Martinez et al. (2012MARTINEZ, D.T.; RESENDE, M. D. V.; COSTA, R. B. da; HIGA, A. R.; SANTOS, G. A. dos; FIER, I. S. N. Estudo da interação genótipo x ambiente em progênies de Pinus taeda por meio da análise de parâmetros genéticos. Floresta, v. 42, n. 3, p. 539 - 552, 2012.). As in the above, the variance of the genotype x environment interaction was of low magnitude (0.21), which resulted in a low coefficient of determination for the effects of genotype x environment interaction (0.03). For the adjusted DBH data, the correlation between the four sites was low (0.52), indicating loss in gain from performing selection at one site for use in distinct locations, as the best individuals in the environment may not be in another one (RESENDE, 2007bRESENDE, M.D. V. Matemática e estatística na análise de experimentos e no melhoramento genético. Colombo: Embrapa Florestas, 2007b. 362p.) (TABLE 8).

It is emphasized that if the data were not corrected for the competition effect, only one breeding zone for the tested sites would erroneously be established. When correcting the data for competition effects, for better efficiency of indirect gains, two breeding areas should be established, ne for test sites B and C and another for sites A and D (TABLE 9).

TABLE 9
Genetic correlation between sites and heritabilities estimated by means of joint analysis for DBH in a progeny test of P. taeda in the state of Santa Catarina.

CONCLUSION

The phenotype models of competition using covariates are efficient for correcting the data for competition effects. The covariate IC was more efficient for the test at sites A and D, the covariate NV for the test at site B and the covariate MAT for the test site C. That is, there is not a single model that fits all locations, making it necessary to examine individually which covariate best fits the data set based on autocorrelation coefficients and variance components.

The residual autocorrelation coefficient and Durbin-Watson test can be used to identify the presence or absence of competition in P. taeda progeny tests and the effectiveness of the covariate used in the model.

The use of these two tools is of great importance in the analysis of data in P. taeda progeny tests, since there is a great discrepancy in the individual classification of the progenies, in the genotypic effects and in the definition of breeding areas between the unadjusted data and the adjusted for the effects of competition.

ACKNOWLEDGMENT

We thank CAPES / CNPQ for the financial support, Mobasa Reforestation S / A and the Laboratory of Genetics and Forest Improvement of UFPR.

REFERENCES

  • ALVARES, C. A.; STAPE, J. L.; SENTELHAS, P. C.; de MORAES, G.; LEONARDO, J.; SPAROVEK, G. Köppen’s climate classification map for Brazil. Meteorologische Zeitschrift , v. 22, n. 6, p. 711-728, 2013.
  • CAPPA, E. P.; CANTET, R. J. Direct and competition additive effects in tree breeding: Bayesian estimation from an individual tree mixed model. Silvae Genetica, v. 57, n. 2, p. 45-55, 2008.
  • CAPPA, E. P., STOEHR, M. U., XIE, C. Y., & YANCHUK, A. D. Identification and joint modeling of competition effects and environmental heterogeneity in three Douglas-fir (Pseudotsuga menziesii var. menziesii) trials. Tree Genetics & Genomes, v. 12, n. 6, p. 102, 2016.
  • COSTA, J. S.; KERR, R. Accounting for competition in genetic analysis, with particular emphasis on forest genetic trials. Tree Genetics & Genomes , v. 9, n. , p. 1-17, 2013.
  • LEONARDECZ- NETO, E.; VENCOVSKY, R.; SEBBENN, A. M. Ajuste para a competição entre plantas em teste de progênies e procedências de essências florestais. Scientia Forestalis, n. 63, p. 136-149, 2003.
  • MARTINEZ, D.T.; RESENDE, M. D. V.; COSTA, R. B. da; HIGA, A. R.; SANTOS, G. A. dos; FIER, I. S. N. Estudo da interação genótipo x ambiente em progênies de Pinus taeda por meio da análise de parâmetros genéticos. Floresta, v. 42, n. 3, p. 539 - 552, 2012.
  • PALLANT, J. SPSS survival manual: A step by step guide to data analysis using SPSS Maidenhead, 2010. 354 p.
  • PAVAN, B. E.; DE PAULA, R. C.; PERECIN, D.; CANDIDO, L. S.; SCARPINATI, E. A. Efeito competicional em testes de progênies de eucalipto. Ciência Florestal, v. 22, n. 3, p. 433-443, 2012.
  • RAMALHO, M. A. P.; FERREIRA, D. F.; OLIVEIRA, A. C. Experimentação em genética e melhoramento de plantas 3ª Ed. Lavras: Editora UFLA, 2012. 305p.
  • RESENDE, M. D. V. Genética biométrica e estatística no melhoramento de plantas perenes. Brasília: Embrapa. Informação Tecnológica, 2002. 975p.
  • RESENDE, M. D. V. Selegen-Reml/Blup: Sistema Estatístico e Seleção Genética Computadorizada via Modelos Lineares Mistos Colombo: Embrapa Florestas, 2007a. 361p.
  • RESENDE, M.D. V. Matemática e estatística na análise de experimentos e no melhoramento genético. Colombo: Embrapa Florestas, 2007b. 362p.
  • RESENDE, M.D. V.; STRINGER, J.; CULLIS, B.; THOMPSON, R. joint modelling of competition and spatial variability in forest field trials. Revista Matemática e Estatística, v.23, n.2, p.7-22, 2005.
  • RESENDE, R. T.; MARCATTI, G. E.; PINTO, D. S.; TAKAHASHI, E. K.; CRUZ, C.D.; RESENDE, M. D. V. Intra-genotypic competition of Eucalyptus clones generated by environmental heterogeneity can optimize productivity in forest stands. Forest Ecology and Management, v. 380, p. 50-58, 2016.
  • STOEHR, M.; BIRD, K.; NIGH, G.; WOODS, J.; YANCHUK, A. Realized genetic gains in coastal Douglas-fir in British Columbia: implications for growth and yield projections. Silvae Genetica , v. 59, n. 5, p. 223-233, 2010.
  • STURION, J. A. e RESENDE, M. D. V. Eficiência do delineamento experimental e capacidade de teste no melhoramento genético da erva-mate (Ilex paraguariensis St. Hil). Pesquisa Floresta l Brasileira, n. 50, p. 3-10, 2005.
  • YE, T. Z.; JAYAWICKRAMA, K. J. S. Efficiency of using spatial analysis in first-generation coastal Douglas-fir progeny tests in the US Pacific Northwest. Tree Genetics & Genomes , v. 4, n. 4, p. 677-692, 2008.

Publication Dates

  • Publication in this collection
    Jul-Sep 2017

History

  • Received
    23 May 2017
  • Accepted
    11 Sept 2017
UFLA - Universidade Federal de Lavras Universidade Federal de Lavras - Departamento de Ciências Florestais - Cx. P. 3037, 37200-000 Lavras - MG Brasil, Tel.: (55 35) 3829-1706, Fax: (55 35) 3829-1411 - Lavras - MG - Brazil
E-mail: cerne@dcf.ufla.br