Acessibilidade / Reportar erro

Adaptability and stability of fruit set and production of peach trees in a subtropical climate

Abstract

Brazilian peach breeding programs have been established to improve peach [Prunus persica (L.) Batsch] production, yield consistency, quality, and disease resistance. Every year several genotypes are selected and their traits must be assessed. This study aimed to evaluate adaptability and stability of fruit set and production of peach genotypes in a subtropical climate, using the GGE biplot methodology. The experimental design was completely randomized with three replicates (trees) in a factorial arrangement of 29 × 3 for genotype and growing season, respectively. The genotypes 'Conserva 1129', 'Rubimel', 'Kampai', 'Tropic Beauty', and 'Cascata 967' had the greatest adaptability and stability for fruit set. The genotypes 'Conserva 681', 'Santa Áurea', 'Atenas', 'Kampai', 'Cascata 962', 'Tropic Beauty' and 'Cascata 967' had the greatest production adaptability and stability. The GGE-biplot methodology classified the peach tree genotypes with regard to adaptability and stability of fruit set and production.


GENETICS AND PLANT BREEDING

Adaptability and stability of fruit set and production of peach trees in a subtropical climate

Idemir CitadinI,* * Corresponding author < idemir@utfpr.edu.br> ; Silvia ScariottoI; Marcos Robson SachetI; Fábio José RosaI; Maria do Carmo Bassols RaseiraII; Américo Wagner JuniorI

IFederal University of Technology-Paraná – Dept. of Agronomy, Via do Conhecimento, km 01 – 85503-390 – Pato Branco, Paraná – Brazil

IIEmbrapa Temperate Agriculture, BR 392, km 78 – 96010- 971 – Pelotas, RS – Brazil

ABSTRACT

Brazilian peach breeding programs have been established to improve peach [Prunus persica (L.) Batsch] production, yield consistency, quality, and disease resistance. Every year several genotypes are selected and their traits must be assessed. This study aimed to evaluate adaptability and stability of fruit set and production of peach genotypes in a subtropical climate, using the GGE biplot methodology. The experimental design was completely randomized with three replicates (trees) in a factorial arrangement of 29 × 3 for genotype and growing season, respectively. The genotypes 'Conserva 1129', 'Rubimel', 'Kampai', 'Tropic Beauty', and 'Cascata 967' had the greatest adaptability and stability for fruit set. The genotypes 'Conserva 681', 'Santa Áurea', 'Atenas', 'Kampai', 'Cascata 962', 'Tropic Beauty' and 'Cascata 967' had the greatest production adaptability and stability. The GGE-biplot methodology classified the peach tree genotypes with regard to adaptability and stability of fruit set and production.

Introduction

The spread of commercial peaches [Prunus persica (L.) Batsch] to subtropical and tropical regions has increased rapidly in recent years. This increase is especially noticeable in Brazil, where climatic conditions are highly variable. Subtropical humid zones, located in Southern Brazil, have hot, humid, rainy summers that favor disease development. Furthermore, inconsistent winter dormancy conditions, caused by conflicting air masses of tropical and polar origins, result in both insufficient chill accumulation in some years or sites and a late frost risk during bloom. Temperatures above 20 ºC during the endodormancy period, considered undesirable (Erez et al., 1979), are also frequent. The cultivars better adapted to Brazil have low chilling requirements (0 to 400 chilling hours below 7.2 ºC) (Byrne et al., 2000; Topp et al., 2008). Therefore, Brazilian peach breeding programs have been introduced designed to improve peach production, yield consistency, quality, and disease resistance (Raseira et al., 2003; Medeiros et al., 2011; Sachet et al., 2013).

Mild winters with a wider range of temperature conditions during dormancy have been related to several physiological and biological disturbances leading to the appearance of a high level of floral anomalies and flower bud drops, which compromise blooming (Lamp et al., 2001; Bassi et al., 2006; Oukabli and Mahhou, 2007; Nava et al., 2009; Couto et al., 2010). Chilling, when insufficient, poses problems for trees that have to overcome endodormancy and can consequently result in low fruit set and yield for genotypes not suitably adapted.

The study of genotype × environment interaction is important for success in temperate fruit cultivation (Gil, 1997). In this sense, knowledge of the adaptability and stability of fruit set and production in different locations or years is crucial when the objective is to increase crop production potential. Thus, the aim of this work was to evaluate adaptability and stability for fruit set and production of peach genotypes in a subtropical climate using the GGE biplot methodology.

Materials and Methods

The evaluations were made from 2008 to 2010 in Pato Branco, in the state of Paraná, Brazil (26º10' S; 52º41' W, 764 m a.s.l.). The local climate is classified as subtropical humid (Cfa – by Köeppen Classification). The number of hours with temperatures below 12 ºC or 7.2 ºC or above 20 ºC was recorded from May to Aug (Figure 1). Temperatures below 7.2 ºC are traditionally used to determine chilling accumulation (Weinberger, 1950). However, nowadays temperatures below 12 ºC are also effective for chilling accumulation (Erez and Couvillion, 1987; Citadin et al., 2002; Chavarria et al., 2009), especially for low chilling genotypes. Temperatures above 20 ºC are undesirable during the dormant period, as they work against the accumulation of chilling (Erez et al., 1979).


Twenty-nine peach genotypes (Figure 2), most of which originated in Pelotas, in the state of Rio Grande do Sul, Brazil (31º52' S; 52º21' W; 13 m a.s.l.), were evaluated during the 2008/2009, 2009/2010 and 2010/2011 growing seasons. Each genotype was represented by three plants. Trees were grown in a 4 m × 5 m spacing and pruned twice a year (winter and summer). The experiment used a completely randomized design. Five one-year-old shoots, 25 to 30 cm long, per plant were randomly selected from around the canopy for fruit set. The trees were five years old in 2008. The orchard received standard fungicide and insecticide sprays, pruning and fertilization, similar to the treatments used in commercial orchards. No chemical means were used to break dormancy.


The evaluation of fruit set was carried out between 30 and 35 days after full bloom. The evaluation of tree production was made by multiplying the total number of fruits by the average fruit weight of each genotype and the results were expressed in kg per tree. For fruit numbering, all fruits in three trees of each genotype were counted close to harvest. The weight average was measured collecting 30 fruits at the beginning, at full and at the end of harvest totaling 90 fruits per genotype per year.

To analyze the adaptability and stability of fruit set and tree production, the GGE (genotype main effect plus genotype by environment interaction effect) biplot methodology (Yan and Rajcan, 2002) was used, based on the following model:

Yij - j = λ1ξi1ηj1 + λ2ξi2ηj2 + εij

where Yij represents the average of genotype i in the year j; j is the mean of all genotypes in the year j; λ1ξi1ηj1 the first principal component (PC1); λ2ξi2ηj2 the second principal component (PC2); λ1 and λ2 are the self values associated with PC1 and PC2, respectively; ξi1 and ξi2 the scores of PC1 and PC2, respectively, for genotype i; ηj1 and ηj2 the self values associated with PC1 and PC2, respectively, for the year j; and εij is the error ij associated with the model (Yan and Kang, 2003). In the GGE biplot method, only the genetic effect and the genotype × environment interaction are considered to be relevant, and both must be considered simultaneously in the evaluation of the cultivars. The two main axes represent most of the variation in the data, considering the environment effect as fixed, i.e., the variation in fruit set or tree production would be due only to genotype and the genotype × environment interaction (Yan and Rajcan, 2002).

In each graph, a polygon was constructed to join the points that represent the most distant genotypes in relation to the origin of the axes. Later, perpendicular lines were designed for each polygon edge passing through the origin, separating it into sections. The genotypes in each sector showed the best performance in environments/ years included in that sector (Yan and Kang, 2003).

A PC1 value near the origin indicates that the genotypes have means close to the general mean (represented by the origin of the lines). As the value becomes more distant and to the right of the origin, the greater the value of the variable can be considered and the more adapted the genotypes (in this case, for fruit set or tree production). A PC2 value near zero indicates more stable genotypes. The graphical biplot may also identify the ideal environment (year), indicated by a year that has a high value for PC1 and a value near zero for PC2 (Yan and Kang, 2003).

The normality and homogeneity of the data were checked by Shapiro-Wilk and Bartlett's tests, respectively. Then, an analysis of variance (ANOVA) was conducted, with genotypes grouped by Scott-Knott's grouping test (p < 0.05) using GENES Program Software for PC (Cruz, 2006). The graphs were generated using R Program Software for PC (R Development Core Team, 2010 package "GGEBiplotGUI").

Results and Discussion

In 2008, chilling accumulation was 424 h below 12 ºC and 96 h below 7.2 ºC, while the temperatures above 20 ºC amounted to 659 h. This year was the one with the highest accumulation of temperatures above 20 ºC among the years studied. May was mild, with night temperatures below 12 ºC and short periods with diurnal temperatures above 20 ºC. Chilling accumulation was observed from the end of May until mid-June, when frost occurred. After this, a tendency towards heat accumulation was observed. Of the years studied, this was the year with the greatest sum of temperatures above 20 ºC compared with temperatures below 12 ºC (Figure 1).

In 2009, chilling accumulation was 744 h below 12 ºC and 157 h below 7.2 ºC, while the temperatures above 20 ºC amounted to 484 h. No long period of heat, with temperatures above 20 ºC, was observed, particularly in June. Although not intense, the chilling was continuous, without abrupt temperature fluctuations (no periods of heat). Dangerous frosts were observed at the end of July (Figure 1). In 2010, chilling accumulation was 732 h below 12 ºC and 158 h below 7.2 ºC, while temperatures above 20 ºC amounted to 531 h. In this year, chilling accumulation was observed only from mid- July, with the occurrence of no dangerous frost (Figure 1).

Genotypes, growing season and the interaction of genotype × growing season were significant for fruit set and production (Table 1). During the three growing seasons, the response of fruit set and production differed among genotypes, which indicates the existence of a genotype-environment interaction. This interaction is expressed by the GGE biplot model.

The SREG (sites regression) obtained by genotype and genotype-environment interaction (GGE biplot) was performed to provide an easy visualization of the complete data set in a reduced dimension plot. Principal components analyses performed in this work showed that 81 % and 84 % of the total variability observed for fruit set and tree production, respectively, could be explained by the first two components: PC1 and PC2. Thus, PC1 and PC2 are expressive and can explain the genotype and genotype-environment interaction (Gauch Jr., 1988). In the GGE biplot model the effects of genotype and genotype-environment interaction are multiplicative.

The PC1 was directly connected to the genotypic effect allowing for almost perfect correlation between the variable analyzed and the PC1 score. Therefore, the PC1 represented 49 % and 60 % of the original data variability in the experiment for fruit set (Figure 2) and tree production (Figure 3), respectively. The difference (51 %, for fruit set, and 40 %, for tree production) is related to interaction genotypes × growing seasons (PC2: 31 %, for fruit set, and 24 %, for tree production) and error (20 %, for fruit set, and 16 %, for tree production).


For fruit set, the year 2009 was the environment that could be considered the closest to the ideal because it was the most stable (near the PC2 axis origin) and promoted good fruit set (Table 2). The years 2008 and 2010 were the most unstable (far from the PC2 axis origin) with poor fruit set (Table 2), which allowed for better discrimination among the genotypes (Figure 2).

The SREG-GGE biplot graph for fruit set was divided into seven sectors. The principal sectors grouped the years 2008 and 2009, and the following genotypes: 'Conserva 1129' (22), 'Cascata 967'(18), 'Rubimel'(2), 'Tropic Beauty' (1), 'Cascata 962' (15); 'Conserva 655' (27), 'Conserva 681' (25) and especially 'Conserva 1129' (22) with the greatest fruit set. The other sector grouped the year 2010 and the genotypes 'Kampai' (13), 'Santa Áurea' (28), 'Atenas' (10), 'Conserva 1223' (16), 'Cascata 1055' (17), 'Conserva 977' (26), with 'Kampai' being the genotype with the greatest fruit set for this year. Genotypes 'Conserva 681' (25), 'Cascata 1065' (4) and 'Âmbar' (9) had fruit set near the general average (50 %) and remained stable during the three years of evaluation. The genotypes 'Conserva 1153' (8) 'Bonão' (4) and 'Cascata 587' (23) had the lowest fruit set (Figure 2 and Table 2).

For tree production, the year 2008 was the most stable (near the PC2 axis origin) and promoted reasonable production for most of the genotypes, with an average of 45.6 kg per tree. The years 2009 and 2010 were the most unstable (far from the PC2 axis origin), which allowed for better discrimination among the genotypes, 2010 being the year with the greatest average production in the years studied (Figure 3 and Table 2).

The SREG-GGE biplot graph for tree production (Figure 3) was divided into six sectors. The principal sectors grouped the years 2008 and 2009, and the following genotypes: 'Santa Áurea' (28), 'Cascata 962' (15), 'Âmbar' (9), 'Conserva 1186' (21), 'Conserva 1129' (22) and 'Kampai' (13). The other sector grouped the year 2010 and the genotypes 'Conserva 681' (25), 'Atenas' (10), 'Tropic Beauty' (1), 'Cascata 967' (18), 'Conserva 655' (27) and 'Conserva 977' (26).

The genotypes 'Rubimel' (2), 'Bonão' (4), 'Conserva 871' (20), 'Cascata 1055' (17), 'Conserva 688' (29), 'Cascata 1070' (19), 'Libra' (12), 'Conserva 1187' (7) and 'Tropic Snow' (5) had a similar level of production when compared with each other, but were below the overall average (48.5 kg per tree). 'Cascata 587' (23) and 'Conserva 1153' were the most unstable (Figure 3). The first had the greatest chilling requirement among the genotypes studied, and the second had serious problems with floral bud abortion (data not shown).

Temperatures above 20 ºC during the endodormancy period and temperatures above 24 ºC near the bloom period, as observed in 2008 (Figure 1), may explain the low fruit set and production in that year (Table 2). Dangerous frosts after fruit set, as observed in 2009 (Figure 1), were the cause of instability and low production observed in this year (Figure 3 and Table 2). In Southern Brazil, the occurrence of frost usually coincides with the period of flowering and fruit development (Raseira et al., 1992) and is one of the serious problems found in peach cultivation. On the other hand, chilling negation by high temperatures is a cause of low budbreak in peach (Richardson et al., 1974; Erez et al. 1979). Temperatures above 24 ºC during pre-blooming and blooming are related to abnormal development of pollen and ovules (Nava et al., 2009; Couto et al., 2010), leading to low flower fertilization and, consequently, low fruit set.

In 2010, despite having the same chill accumulation as in 2009 (Figure 1), it occurred later, and was not as effective in overcoming dormancy of flower buds, with heat accumulation during pre-blooming period, from the beginning of to mid July (Figure 1), causing lower fruit set than in 2009, but more than in 2008 (Table 2). Frost in late July was not dangerous, but caused fruit drop in a few genotypes that had already reached full bloom. The heat, during pre-blooming, and the frost, after blooming, were responsible for instability observed in both fruit set and production in 2010 (Figure 2 and Figure 3).

Rodrigo and Herrero (2002) observed that an increase in temperature accelerated flower bud development causing a hastening in flowering time and following hand-pollination, reduced fruit set in apricot. At anthesis, flowers that had developed in warmer conditions weighed less and showed less development of the pistil than control flowers. 'Tropic Beauty' (1), 'Kampai' (13), 'Cascata 967' (18), 'Santa Áurea' (28), 'Atenas' and 'Conserva 681' (25) appear to tolerate higher temperatures during the dormancy period, thus maintaining good fruit set (Figure 2) and production (Figure 3) even in years with low chill and high heat accumulation, as in 2008.

The development of varieties with low chilling requirement with good tolerance to high temperatures during flowering, make production more consistent (Sherman and Lyrene, 1984), especially in subtropical regions in low latitudes, in which temperatures above 20 ºC often occur. 'Conserva 1153' (8) and 'Cascata 587' (23) had low fruit set (Figure 2) and tree production (Figure 3) in the years studied. Ensuring good fruit set is important, but this alone is not sufficient for good production. Production is a function of orchard density, flower density, index of fruit set and fruit size for all species, varieties and climatic conditions (Dennis Jr., 1981).

Excellent yields were obtained when there is good fruit set and fruit size. However, after fruit set, the number of fruits per plant is generally standardized by thinning practices; a fact that occurred in this experiment, therefore, from the thinning, production will be defined by fruit size of each genotype. Yamaguchi et al. (2002) found both cell number and cell size differences were significantly correlated with peach fruit size, although the correlation coefficients between cell number and fruit weight were consistently higher.

In sweet cherry, differences in mesocarp thickness were correlated with both cell number and cell size (Yamaguchi et al., 2004). Mesocarp cell number was remarkably stable and virtually unaffected by the environment as neither growing location nor physiological factors that reduced final fruit size altered the cell numbers. Cell length was also different among the cultivars (p < 0.05), but failed to contribute to the overall difference in fruit size. Cell length was influenced by the environment, indicating that cultural practices that maximize mesocarp cell size should be used to achieve a fruit size potential of the cultivar (Olmstead et al., 2007).

Other related genotype-dependent floral biology factors affect fruiting in temperate fruit species, such as flower bud production, flower production, flower bud drop, flowering time, percentage pollen germination, height difference between the stigma and the superior plane of the anthers, percentage of aborted pistils, autogamy level, fruit set and fruit production (Ruiz and Egea, 2008). Some of these characteristics may have contributed to the reduction of the average fruit set and production in the genotypes studied, even those genotypes that had sufficient chilling accumulation. These causes should be elucidated in future work.

In conclusion, genotypes 'Conserva 1129', 'Rubimel', 'Kampai', 'Tropic Beauty' and 'Cascata 967' demonstrated the greatest adaptability and stability for fruit set. The genotypes 'Conserva 681', 'Santa Áurea', 'Atenas', 'Kampai', 'Cascata 962', 'Tropic Beauty' and 'Cascata 967' demonstrated the greatest adaptability and stability of production (kg per tree). The GGE-biplot methodology was successful in classifying peach tree genotypes with regard to adaptability and stability of fruit set and production.

Acknowledgments

This research was supported in part by, CNPq (Brazilian National Council for Scientific and Technological Development), the Araucária Foundation, Embrapa Temperate Agriculture and CAPES (Coordination for the Improvement of Higher Level Personnel).

Received June 12, 2013

Accepted October 31, 2013

Edited by: Antonio Costa de Oliveira

  • Bassi, D.; Bartolini, S.; Viti, R. 2006. Recent advances on environmental and physiological challenges in apricot growing. Acta Horticulturae 717: 23-31.
  • Byrne, D.H.; Wayne, B.; Bacon, T.A. 2000. Stone fruit genetic pool and its exploitation for growing under warm winter conditions. p. 157-230. In: Erez, A., ed. Temperate fruit crops in warm climate. Kluwer Academic, Dordrecht, The Netherlands.
  • Chavarria, G.; Herter, F.G.; Raseira, M.C.B.; Rodrigues, A.C.; Reisser, C.; Silva, J.B. 2009. Mild temperatures on bud breaking dormancy in peaches. Ciência Rural 39: 2016-2021.
  • Citadin, I.; Raseira, M.C.B.; Herter, F.G.; Silveira, C.A.P. 2002. Avaliation of chilling requirement in peach. Revista Brasileira de Fruticultura 24: 703-706 (in Portuguese, with abstract in English).
  • Couto, M.; Raseira, M.C.B.; Herter, F.G.; Silva, J.B. 2010. influence of high temperatures at blooming time on pollen production and fruit set of peach cvs. Maciel and Granada. Acta Horticulturae 872: 225-230.
  • Cruz, C.D. 2006. Genes Program: Biometry = Programa Genes: Biometria. Editora UFV, Viçosa, MG, Brazil (in Portuguese).
  • Dennis Jr., F.G. 1981. Limiting factors in fruit set of 'Delicious' apple. Acta Horticulturae 120: 119-124.
  • Erez, A.; Couvillon, G.A. 1987. Characterization of the influence of moderate temperatures on rest completion in peach. Journal of the American Society for Horticultural Science 112: 677-680.
  • Erez, A.; Couvillon, G.A.; Hendershott, C.H. 1979. Quantitative chilling enhancement and negation in peach buds by high temperatures in a daily cycle. Journal of the American Society for Horticultural Science 104: 536-540.
  • Gauch Jr., H.G. 1988. Model selection and validation for yield trials with interaction. Biometrics 44: 705-715.
  • Gil, G. 1997. Fruticulture: The Productive Potential = Fruticultura: El Potencial Productivo. Ediciones Universidad Católica de Chile, Santiago, Chile (in Spanish).
  • Lamp, B.M.; Connell, J.H.; Duncan, R.A.; Viveros, M.; Polito, V.S. 2001. Almond flower development: floral initiation and organogenesis. Journal of the American Society for Horticultural Science 126: 689-696.
  • Medeiros, J.G.S.; Citadin, I.; Santos, I.; Assmann, A.P. 2011. Reaction of peach tree genotypes to bacterial leaf spot caused by Xanthomonas arboricola pv. pruni Scientia Agricola 68: 57- 61.
  • Nava, G.A.; Dalmago, G.A.; Bergamaschi, H.; Paniz, R.; Santos, R.P.; Marodin, G.A.B. 2009. Effect of high temperatures in the pre-blooming and blooming periods on ovule formation, pollen grains and yield of 'Granada' peach. Scientia Horticulturae 122: 37-44.
  • Olmstead, J.W.; Iezzoni, M.F.; Whiting, M.D. 2007. Genotypic differences in sweet cherry fruit size are primarily a function of cell number. Journal of the American Society for Horticultural Science 132: 697-703.
  • Oukabli, A.; Mahhou, A. 2007. Dormancy in sweet cherry (Prunus avium L.) under Mediterranean climatic conditions. Biotechnology, Agronomy, Society and Environment 11: 133- 139.
  • Raseira, M.C.B.; Herter, F.; Silveira, C.A.P. 2003. The Embrapa Clima Temperado peach breeding program and adaptation to subtropical regions. Acta Horticulturae 606: 45-50.
  • Raseira, M.C.B.; Baptista, J.S.; Herter, F.G.; Peters, J.A. 1992. Peach tree, Prunus persica L. Batsch, flower buds cold sensibility. Revista Brasileira de Fruticultura 14: 167-172 (in Portuguese, with abstract in English).
  • Richardson, E.A.; Seeley, Sd.; Walker, D.R. 1974. A model for estimating the completion of rest for 'Redhaven' and 'Elberta' peach trees. HortScience 1: 331-332.
  • Rodrigo, J.E.; Herrero, M. 2002. Effects of pre-blossom temperatures on flower development and fruit set in apricot. Scientia Horticulturae 92: 125-135.
  • Ruiz, D.; Egea, J. 2008. Analysis of the variability and correlations of floral biology factors affecting fruit set in apricot in a Mediterranean climate. Scientia Horticulturae 115: 154-163.
  • Sachet, M.R.; Citadin, I.; Scariotto, S.; Santos, I.; Zydek, P.H.; Raseira, M.C.B. 2013. Reaction of peach genotypes to bacterial leaf spot: correlations with environmental conditions, leaf phenology and morphology. HortScience 48: 28-33.
  • Sherman, W.B.; Lyrene, P.M. 1984. Biannual peaches in the tropics. Fruit Varieties Journal 38: 37-39.
  • Topp, B.L.; Sherman, W.B.; Raseira, M.C.B. 2008. Low-chill cultivar development. p. 106-138. In: Layne, D.R.; Bassi, D., eds. The peach: botany, production and uses. CAB international, Wallingford, UK.
  • Weinberger, J.H. 1950. Chilling requirements of peach varieties. Proceedings of the American Society for Horticultural Science 56: 122-128.
  • Yamaguchi, M.; Haji, T.; Miyake, M.; Yaegaki, H. 2002. Studies on the varietal differences and yearly deviation of mesocarp cell numbers and lengths and fruit weight among commercial peach [Prunus persica (L.) Batsch] cultivars and selections, wild types, and their hybrids. Journal of the Japanese Society for Horticultural Science 71: 459-466.
  • Yamaguchi, M.; Sato, I.; Takase, K.; Watanabe, A.; Ishiguro, M. 2004. Differences and yearly variation in number and size of mesocarp cells in sweet cherry (Prunus avium L.) cultivars and related species. Journal of the Japanese Society for Horticultural Science 73: 12-18.
  • Yan, W.; Kang, M.S. 2003. GGE Biplot analysis: a graphical tool for breeders, geneticists, and agronomists. CRC Press LLC, Boca Raton. FL, USA.
  • Yan, W.; Rajcan, I. 2002. Biplot analysis of test sites and trait relations of soybean in Ontario. Crop Science 42: 11-20.
  • *
    Corresponding author <
  • Publication Dates

    • Publication in this collection
      14 Apr 2014
    • Date of issue
      Apr 2014

    History

    • Accepted
      31 Oct 2013
    • Received
      12 June 2013
    Escola Superior de Agricultura "Luiz de Queiroz" USP/ESALQ - Scientia Agricola, Av. Pádua Dias, 11, 13418-900 Piracicaba SP Brazil, Phone: +55 19 3429-4401 / 3429-4486 - Piracicaba - SP - Brazil
    E-mail: scientia@usp.br