On-line version ISSN 1678-3921
Pesq. agropec. bras. vol.34 no.7 Brasília July 1999
Phenotypic stability and genetic gains in six-year girth growth of Hevea clones1
ABSTRACT - Rubber tree [Hevea brasiliensis (Willd. ex Adr. de Juss.) Müell. Arg.] budgrafts of seven clones were evaluated on five contrasting sites in the plateau region of the São Paulo State, Brazil. The objective of this work was to study the phenotypic stability for girth growth. The experimental design was a randomized block design with three replications and seven treatments. Analysis of variance of girth at six-year plant growth indicated a highly significant clone x site interaction. Only linear sites and clone x site components of clone x year interaction were significant, indicating that the performance of clones over sites for this trait could be predicted. The clones GT 1 and PB 235 showed the greatest stability in relation to girth growth, with foreseen responses to change, introduced in the sites. The clones PB 235 and IAN 873 showed significative difference in relation to regression coefficient, representing clones with specific adaptability on favorable and unfavorable sites respectively. The clone GT 1 became the most promissory one in the study of stability and adaptability even showing low girth growth. Expected genetic gains from planting sites, along with estimates of clonal variance and repeatability of clonal means are generally greatest or close to the greatest when selection is done at the same site.
Index terms: rubber tree, Hevea brasiliensis, heritability, repeatability.
Estabilidade fenotípica e ganhos genéticos de seis anos de avaliação do perímetro do caule em clones de seringueira
RESUMO - Enxertos de sete clones de seringueira [Hevea brasiliensis (Willd. ex Adr. de Juss.) Muell. Arg.] foram plantados em cinco diferentes locais da região do Planalto do Estado de São Paulo, objetivando o estudo da estabilidade do perímetro do caule. O delineamento experimental foi o de blocos casualizados, com três repetições e sete tratamentos. A análise de variância mostrou que a interação genótipos x locais foi altamente significativa. A maior parte da variação que se refletiu na estabilidade deveu-se à regressão linear das médias dos ambientes. Os clones GT 1 e PB 235 se apresentaram como os mais estáveis em relação ao perímetro do caule, com respostas previsíveis às mudanças introduzidas nos locais. Os clones PB 235 e IAN 873 mostraram diferenças significativas em relação ao coeficiente de regressão, correspondendo a clones de adaptabilidade específica para locais favoráveis e desfavoráveis, respectivamente. O clone GT 1 foi o mais promissor em termos de adaptabilidade e estabilidade, devido ao baixo desvio da regressão. Ganhos genéticos esperados nos locais, junto com as estimativas de variância clonal e repetibilidade, foram geralmente maiores quando a seleção foi feita no mesmo local.
Termos para indexação: Hevea brasiliensis, herdabilidade, repetibilidade.
Rubber tree [ Hevea brasiliensis ( Willd. ex Adr. de Juss.) Müell. Arg.] is one of the indigenous species from the tropical rain forest in the Amazon Basin. Their high commercial importance is well known for their latex production from which rubber is coagulated and processed and as a consequence the dry rubber yield determines the level of output.
Traditionally, rubber in Brazil has been cultivated in the equatorial region, in a zone lying between 10ºN and 10ºS of the Equator. Cultivation was concentrated in this location as these areas provided the ideal environmental conditions, i.e. mean annual temperature of 26ºC ± 2ºC and high rainfall (1,800 to 3,000 mm) distributed throughout the year with 120 - 240 rainy days per year. Therefore, regions outside conventional zones are beginning to establish rubber due to an increasing world demand for the commodity. In the São Paulo State, Brazil, the rubber growing region, considered as an area well recognized as hospitable for the crop, lies between 20ºS and 21ºS, far beyond the normal limits.
Girth is an economically important character and is considered to be a measure of vigour. It determines the age at which a clone can be exploited and defines the shortening of the uneconomic, immature period. Trees with large trunks have high timber values.
Therefore, one of the main tasks in Hevea breeding is the comparative evaluation of clones, either of those bred locally or imported from other countries, over a wide range of environments in the plateau region of São Paulo rubber planting areas. The testing of these planting materials in different environments is necessary due to the environmental influence on girth performance of rubber trees. In addition, it has also been observed that a particular clone which performs well in one environment may not perform as well in a different environment concerning a given set of clones. This phenomenon of clone genotype - environment (GE) interaction can therefore make conclusive assessment on the value of individual clones for final clonal recommendation to the planting industry.
No extensive studies have been performed on the phenotypic stability of vigour. In general, a large amount of variation in girth is observed. According to Pushparajah (1980), many agroclimatic elements, such as dryness, temperature and rainfall, that affect several components of the growth of Hevea, contribute to a large variability.
Clones with various degrees of adaptability/stability were identified in Malaysia (Tan, 1995), as a supplementary guide for Hevea clonal selection. In Sri Lanka, Jayasekera et al. (1977) estimated regression coefficients for first height, second height and survival rate of ten clones on eight sites. In Indonesia, Suhendi (1989) obtained regression estimates of sixteen clones tested in three locations. These authors illustrated the use of information as a guide in clonal planting recommendation .
Several methods have been proposed to analyze genotype x environment interactions (Freeman, 1973), but joining regression technique has been widely applied to agronomic stability (Finlay & Wilkinson, 1963; Eberhart & Russell, 1966; Perkins & Jinks, 1968). Although girth stability is one of the most important breeding objectives, a statistical measure of this trait is very seldom considered in a breeding programme of perennial crops like Hevea. The regression coefficient (bi) can be used to describe the general response to environmental conditions (Stringfield & Satter, 1934) and the mean squares deviations () from linear regression actually measure the vigour stability (Yates & Cochran, 1938; Eberhart & Russell, 1966).
In Brazil, Kalil Filho (1982) and Gonçalves et al. (1992) employed the method of Eberhart & Russell (1966) to study temporal yielding, where the environments were represented by different months and years through an environmental index obtained by the medium performance of all the genotype in each month and year.
The present study was conducted to examine the magnitude of the genotype - environment interaction on girth growth of rubber tree at the age of six-years, searching for clone stability and genetic gains.
MATERIAL AND METHODS
The genetic material originally consisted of seven Hevea clones i.e. GT 1, PR 261, PR 255, IAN 873, RRIM 701, PB 235, and RRIM 600 chosen to represent good clones available at the establishing of the experiments. They were budgrafted onto one and a half-year-old rootstock seedlings (Tjir 1 x Tjir 16) raised in nurseries. The successful budgrafts were uprooted and planted in polyethylene bags. After the emission of the first flush of leaves, the experiments were established in the field. The budgraftings were planted at five contrasting sites in the plateau region of the State of São Paulo (Table 1). These sites represent the most important non traditional area of continental climate of rubber production.
The experiment at each site was set in a randomized complete block design with three replications. In all locations, except Matão, thirty plants were used in each plot. At Matão, a larger plot size of forty plants was used, due to the local facilities. Each plot consisted of five rows of six or eight plants. In all sites the girth of each tree was measured at 1.20 m above the highest point of the bud union.
Data from the sixth year girth for all sites combined were submitted to analysis of variance (ANOVA) using the following statistical model:
Yijk is the observed value of the i th clone of jth site in the kth replication;
m is the general mean;
Ci is the fixed effect of the ith clone (i= 1, 2 ... c, c = 7);
Sj is the effects of the jth site (j = 1,2...s, s = 5);
(CS)ij is the interaction between ith clone with the jth site;
R k(j) is the effect of kth randomized block with the jth site;
Eijk is the experimental error associated with ijkth observation.
All effects were calculated through "Expected Mean Square" components given on the right side of the ANOVA (Table 2). Tukey's test was applied to find out entries different from each other.
Estimated variance components were used to calculate phenotypic variance, broad-sense heritability, repeatability of clonal means; and genetic gain was estimated from algebric manipulation of the mean square (Vencovsky & Barriga, 1992), worked out as follows:
, and are estimated variance components (Table 1);
is the phenotypic variance;
is the phenotypic variance of clonal means;
is the broad-sense heritability;
is the repeatability of clonal means;
s is the number of sites;
r is the number of replications per site.
Genetic gain was calculated as (Falconer, 1981), where i is the intensity of selection and is the phenotypic standard deviation. In this study, i was assumed to be 0.82, which corresponds to selection of three clones out of seven (Becker, 1984). As the clone x site variance component increases relative clonal variance, the repeatability of clonal means decreases, and thus, genetic gains are reduced. One solution is to select clones with high stability. The technique of Eberhart & Russell (1966) was used in this study to characterize genotypic stability. The following linear regression model was used:
Yij is the mean of the clone ith at the site j;
m is the general mean of clone i;
bi is the regression coefficient of the ith clone at the site index which measures the response of this clone to varying sites;
Ij is the environmental index which is defined as the mean deviation of all clones at a given site from the overall mean;
dij is the deviation from regression of the ith clone in the jth site;
is the mean of experimental error.
With this methodology two parameters of stability were estimated: a) the regression coefficient, which is the regression of the performance of each clone under different site on the annual mean over all the clones, and b) the variance of the deviation () from the regression line.
The contribution of sites in generating interactions was studied using a method developed by Burdon (1977). It considers growth in different sites as separate traits and calculates the genetic correlations between pairs of sites. In this evaluation, girth was used instead of height. Gain for selection at one site with planting at another is determined from the formula for correlated response to selection (Falconer, 1981). Estimated correlations can then be used to indicate which sites are the most similar to one another ( in terms of showing least interaction) and a matrix of expected gains can indicate which sites are best for testing.
The formula for expected gain at site y based on clonal selection at site x can be written as:
i is the intensity of selection;
is the repeatibility of clonal means at site x;
is the clonal component of variance at site y;
rxy is the correlation between clonal means at sites x and y.
Variance components and heritability of clonal means were estimated from analyses of variance performed separately for each site. The statistical model considered effects due to clone, repetition and error.
RESULTS AND DISCUSSION
Overall means and analyses of variance
All effects, including clone x site interaction, were highly statistically significant (p < 0.01) in the overall analyses of variance (Table 3). There were significant differences among test sites for six-year girths. The significant clone x site interaction indicates that clones perform differently between sites. Of greater interest other than statistical significance is the importance of the interaction in reducing gains. The magnitude of the interaction component was almost the same to the clonal variance component (Table 3). Site means are given in Table 4. Matão was the top-ranked site with a mean of nearly 43.00 centimeters, 38% better than the bottom-ranked site of Ribeirão Preto. The facilities of a more hard budgrafts selection conducted in Matão before setting the experiment in the field was one of the reasons of higher growth. Overall clonal means ranged from 34.82 cm for clone RRIM 701 to 39.34 cm for clone PB 235. Clone PB 235 was particularly outstanding, ranking first at two sites and second at the other three sites (Table 4).
Broad sense heritability () overall for six-year girth growth was 0.23, which is relatively low. However, repeatability of clonal means () overall was 0.79. The genetic gain based on a selection intensity of i= 0.82 was 5.55 cm, a gain of 15.26 percent above the overall mean. If the three top clones are selected and genetic gain is calculated based on the observed selection differential, the genetic gain is 14.40 cm, 40% above de overall mean. The difference between the two gain estimates results from the outstanding girth of clones PB 235 and GT 1. These gain estimates represent gains to be achieved from further clonal selection, and do not include gains to be achieved from past selection.
The combined analysis of variance conducted for six-year girth growth at five sites of the study are presented in Table 5, and revealed that the clones differ significantly. The clone x site interaction components were further partitioned into linear (sites and clone x sites) and non-linear (pooled deviations) components. Mean square for both these components were tested against pooled error mean square. Only the linear component was significant, indicating that the predictable components shared clone x site interactions. Preponderance of linear clone x site interaction is of great practical importance, implying that there are differences among the linear regression coefficient for each clone.
The stability parameters for all clones are given in Table 6. Eberhart & Russell (1966) emphasized the need of considering both linear (bi) and non-linear () components of genotype-environment interactions in judging the stability of a genotype. A wide adaptability genotype was defined as one with bi= 1.0 and high stability as one with = 0. Values for the regression coefficient (bi) range from 0.6868 (IAN 873) to 1.4121 (PB 235). The plot clonal means against site means with the corresponding regression lines for seven clones illustrates the differential reactions of clones to changing environments (Fig. 1). As defined by Eberhart & Russell (1966), clone GT 1, PR 255 and RRIM 701 represent clones of wide adaptability (bi= 1.0). Its performance is relatively equal on poor and good sites. Clone PB 235 represents a clone of specific adaptability (bi=1.4). It performs relatively better on favorable sites while clone IAN 873 represents one of specific adaptability to unfavorable sites (bi= 0.69).
The point is that regression coefficient measures "relative" performance of the clone (Clair & Kleinschmit, 1986). In rubber trees breeding, this information is useful to distinguish clones for specific environments, but if all environments tested are in one planting zone, and each one represents the same proportion of area to be planted, then this information is not important. Selection on the overall mean is all that is necessary to assure the largest overall gains. Therefore clone GT 1 could be considered superior in future breeding programmes in order to incorporate stability. According to Singh & Gupta (1988), it is possible that stable genotypes carry genes for stability, useful in breeding programmes because it facilitates economic production.
More important to rubber tree breeding is the predictability of clone vigor in various environments. This concept of stability may be measured by the mean deviations from regression line, (Breeze, 1969; Becker, 1981). Clone RRIM 600 performs well overall, but is very unstable as measured by = 2.15 (Table 6). On the other hand, clone PB 235 also performs well but is very stable as measured by = 0.08. This is apparent from the deviations of the individual points from the regression line in Fig.1.
A plot of the stability parameter against the clonal means is useful as an aid to selection (Figs. 2 and 3). Clones falling in the lower right-hand side are preferred. Clone GT 1 is the best clone no matter which stability parameter is used. Selection of clones based on girth growth and stability can then proceed using a selection index for multiple trait selection (Stonecypher & Arbes, 1976), but assigning weights based on relative economic value may prove difficult. This would require knowledge of the reduction in gain associated with using unstable versus stable genotypes. The use of independent culling levels, that is, setting an acceptable value for a stability parameter followed by selection based on girth, may prove more pratical.
Correlated gains among sites
The expected gains from various combinations of testing and planting sites, along with estimates of clonal variance and repeatability of clonal means are presented in Table 7. Expected gains at planting sites are generally greater or close to the greatest when selection is done at the same site. Where this is not true (Indiana and Ribeirão Preto), the clonal variance is small when compared to those of Mococa, Matão, and Votuporanga. Selection based on overall means results in gains second only to those when selection is done at the same site as planting.
Matão was the best single site for large testing because of the high correlations with most of the other sites (Table 8). Votuporanga shows fair correlation with Ribeirão Preto (r= 0.53). The estimated gains for selection at Matão and planting at Indiana, Mococa and Ribeirão Preto were, in all cases greater than for selection at either Mococa or Votuporanga.
The least similar sites are Ribeirão Preto and Mococa as indicated by a low correlation of clonal means and low gains from selection and planting between sites.
1. The clone GT 1 recommended to the plateau region of São Paulo State shows the second highest girth average and a value of bi = 1.0, representing a clone of wide adaptability, performing well on poor and good sites.
2. The clones PR 255 and PR 261 show high stability and wide adaptability, but the means of girth are low.
3. The clones PB 235 and IAN 873 are the most stable to girth growth and superior to the average, for favorable and unfavorable sites, respectively.
4. Expected gains from planting sites, along with estimates of clonal variance and repeatability of clonal means are generally greatest or close to the greatest when selection is done at the same site.
To Cambuhy Empreendimentos Agropecuários Ltda and Indiana Rubber Plantation for providing land support for the experiments in Matão and Indiana, respectively; to Mrs. Lígia Regina L. Gouvêa for compiling the data; to Dr. Maria Regina G. Ungaro for providing constructive criticism on the manuscript.
BECKER, H. Correlations among some statistical measures of phenotypic stability. Euphytica, Wageningen, v.30, p.835 -840, 1981. [ Links ]
BECKER, W.A. Manual of quantitative genetics. Pullman: Academic Enterprises, 1984. 290p. [ Links ]
BREEZE, E.L. The measurement and significance of genotype-environment interactions in grasses. Heredity, London, v.24, p.27-44, 1969. [ Links ]
BURDON, R.D. Genetic correlation as a concept for studying genotype-environment interaction in forest tree breeding. Silvae Genetica, Frankfurt, v.26, n.5/6, p.168-175, 1977. [ Links ]
CLAIR, J.B.; KLEINSCHMIT, J. Genotype-environment interaction and stability in ten-year height growth of Norway Spruce clones (Picea abies Karts). Silvae Genetica, Frankfurt, v.35, n.2, p.177-186, 1986. [ Links ]
EBERHART, A.A.; RUSSELL, W.A. Stability parameters for comparing varieties. Crop Science, Madison, v.6, p.36-40, 1966. [ Links ]
FALCONER, D.S. Introduction to quantitative genetics. London: Longman, 1981. 370p. [ Links ]
FINLAY, K.W.; WILKINSON, G.N. The analysis of adaptation in a breeding programme. Australian Journal of Agricultural Research, Melbourne, v.14, p.742-754, 1963. [ Links ]
FREEMAN, J.H. Statistical methods for the analysis of genotype-environment interactions. Heredity, London, v.31, p.339-354, 1973. [ Links ]
GONÇALVES, P. de S.; CARDOSO, M.; SANTOS, I.C.I. dos; MARTINS, A.L.M.; ORTOLANI, A.A.; COLOMBO, C.A. Selection of Hevea mother tree adapted to unpredictable annual climatic variability. Revista Brasileira de Genética, Ribeirão Preto, v.15, n.1, p.137-147, 1992. [ Links ]
JAYASEKERA, N.E.M.; SAMARANAYAKE, P.; KARUNASEKERA, K.B. Initial studies on the nature of genotype-environment interactions in some Hevea cultivars. Journal of Rubber Research Institute of Sri Lanka, Agalawata, v.54, p.33-42, 1977. [ Links ]
KALIL FILHO, A.N. Potencial de produtividade e estabilidade fenotípica na caracterização de clones de seringueira (Hevea spp.). Piracicaba: USP-ESALQ, 1982. 95p. Dissertação de Mestrado. [ Links ]
PERKINS, J.M.; JINKS, J.L. Environmental and genotype x environmental components of variability III: Multiple lines and cross. Heredity, London, v.23, p.339-356, 1968. [ Links ]
PUSHPARAJAH, E. Problems and potentials for establishing Hevea under difficult environmental conditions. Planter, Kuala Lumpur, v.50, p.242-251, 1980. [ Links ]
SINGH, S.V.; GUPTA, P.K. Phenotypic stability for tuber yield potato (Solanum tuberosum L.). Indian Journal of Genetics and Plant Breeding, New Delhi, v.48, p.9-12, 1988. [ Links ]
STONECYPHER, R.; ARBES, M. Methods of selection. In: IUFRO JOINT MEETING ON ADVANCED GENERATION BREEDING, 1976, Bordeaux. Proceedings... Bordeaux: INRA, 1976. p.92-99. [ Links ]
STRINGFIELD, G.H.; SATTER, R.M. Differential response of corn varieties to fertility levels and to seasons. Journal of Agricultural Research, Punjab, v.49, p.991-1000, 1934. [ Links ]
SUHENDI, D.A. A yield stability analysis of rubber clones. Menara Perkebunan, Bogor, v.57, n.1, p.10-13, 1989. [ Links ]
TAN, H. Genotype x environment interaction studies in rubber (Hevea) clones. Journal of Natural Rubber Research, Kuala Lumpur, v.10, n.1, p.63-76, 1995. [ Links ]
VENCOVSKY, R.; BARRIGA, P. Genética biométrica no fitomelhoramento. Ribeirão Preto: Sociedade Brasileira de Genética, 1992. 496p. [ Links ]
YATES, F.; COCHRAN, W.G. The analysis of groups of experiments. Journal of Agricultural Science, London, v.28, p.556-580, 1938. [ Links ]
1 Accepted for publication on July 7, 1998. Partially supported by FAPESP and Secretaria de Agricultura e Abastecimento do Estado de São Paulo.
2 Agronomist, Dr., Embrapa/Programa Integrado de São Paulo, Programa Seringueira, Centro de Café e Plantas Tropicais, Instituto Agronômico de Campinas (IAC), Caixa Postal 28, CEP 13020-970 Campinas, SP. E-mail: email@example.com
3 Forestry Engineer. Fellowship of the CNPq.
4 Agronomist, Dr., Centro de Ecofisiologia e Biofisica, IAC.
5 Agronomist, Dr., Centro de Solos e Recursos Agroam-bientais, IAC.
6 Agronomist, Núcleo de Agronomia do Noroeste, IAC, Votu-poranga, SP.
7 Agronomist, Gerente de Divisão Seringueira/Cambuhy Empreendimentos Agropecuários Ltda.