Acessibilidade / Reportar erro

Phenotypic, additive genetic and environment correlations of maize landraces populations in family farm systems

Correlações fenotípicas, genéticas aditivas e ambientais de populações de milho crioulo em sistemas de agricultura familiar

Abstracts

Knowledge of the association between characters is very important in the genetic breeding programs, but there is limited information about correlations between traits in maize landraces in the current literature. The objectives of this study were to estimate phenotypic, additive genetic and environment correlations among traits in maize landraces, which were cultivated in family farm systems, in order to guide the participatory breeding on maize landraces. Between 294 and 400 half-sib progenies from the six populations were evaluated, subdivided in triple 7X7 and 10X10 lattice design, respectively, with single row plots of 4 m long and spaced 1.00 × 0.20 m. Estimates of correlation coefficients are reported for nine traits: grain yield (GY), relation grain weight per ear weight (GE), ears per plant (EP), days to flower (DF), plant height (PH), ear height (EH), percentage of root loding (%L), percentage of stalk loding (%BS) and percentage of damaged ears (%DE). A high number of additive genetic correlation (rA) was obtained between grain yield and the other traits, in maize landraces populations, cycles and locations. For grain yield, the highest averages of the rA estimate were with EP (0.67), GE (0.47) and %DE (-0.63). In contrast to the current published researches, negative correlations between grain yield and days to flower were observed.

Zea mays L.; population; recurrent selection; participatory breeding


O conhecimento da associação entre caracteres é de grande importância nos programas de melhoramento genético, contudo a literatura apresenta limitadas informações sobre correlações em populações de milho crioulo. Os objetivos deste estudo foram estimar correlações fenotípicas, genéticas aditivas e ambientais entre caracteres de populações de milho crioulo, cultivados em sistemas de agricultura familiar, para orientar o programa de melhoramento genético participativo. Entre 294 e 400 de progênies de meios-irmãos de seis populações, foram avaliadas em látices triplos 7X7 e 10X10, respectivamente, com parcelas de 4 m de comprimento e espaçadas 1,00 × 0,20 m. Os coeficientes de correlações foram estimados para nove características: produtividade de grãos total (GY), relação peso de grãos por peso de espiga (GE), nº de espigas por planta (EP), nº de dias do florescimento (DF), altura da planta (PH), altura espiga (EH), porcentagem de acamamento (%L), porcentagem de colmos quebrados (%BS) e porcentagem de espigas danificadas (%DE). Houve um alto número de correlação genética aditiva (rA) entre produtividade de grãos e as demais características, nas populações de milho crioulo, ciclos e locais. As maiores médias das estimativas rA para o produtividade de grãos foram com EP (0,67), GE (0,47) e DE% (-0,63). Diferentemente dos trabalhos publicados, foram observadas correlações negativas entre produtividade de grãos e dias para florescimento.

Zea mays; população; seleção recorrente; melhoramento genético participativo


GENETICS AND PLANT BREEDING

Phenotypic, additive genetic and environment correlations of maize landraces populations in family farm systems

Correlações fenotípicas, genéticas aditivas e ambientais de populações de milho crioulo em sistemas de agricultura familiar

Luciane Bertoletti BarrosI; Rosângela Maria Pinto MoreiraII; Josué Maldonado FerreiraII,* * Corresponding author < josuemf@uel.br>

IUEL Programa de Pós-Graduação em Genética e Biologia Molecular, C.P. 6001 86051-990 Londrina, PR Brasil

IIUEL/Centro de Ciências Biológicas, Depto. de Biologia Geral, Lab. de Melhoramento Genético de Plantas, C.P 6001 86051-990 Londrina, PR - Brasil

ABSTRACT

Knowledge of the association between characters is very important in the genetic breeding programs, but there is limited information about correlations between traits in maize landraces in the current literature. The objectives of this study were to estimate phenotypic, additive genetic and environment correlations among traits in maize landraces, which were cultivated in family farm systems, in order to guide the participatory breeding on maize landraces. Between 294 and 400 half-sib progenies from the six populations were evaluated, subdivided in triple 7X7 and 10X10 lattice design, respectively, with single row plots of 4 m long and spaced 1.00 × 0.20 m. Estimates of correlation coefficients are reported for nine traits: grain yield (GY), relation grain weight per ear weight (GE), ears per plant (EP), days to flower (DF), plant height (PH), ear height (EH), percentage of root loding (%L), percentage of stalk loding (%BS) and percentage of damaged ears (%DE). A high number of additive genetic correlation (rA) was obtained between grain yield and the other traits, in maize landraces populations, cycles and locations. For grain yield, the highest averages of the rA estimate were with EP (0.67), GE (0.47) and %DE (-0.63). In contrast to the current published researches, negative correlations between grain yield and days to flower were observed.

Key words: Zea mays L., population, recurrent selection, participatory breeding

RESUMO

O conhecimento da associação entre caracteres é de grande importância nos programas de melhoramento genético, contudo a literatura apresenta limitadas informações sobre correlações em populações de milho crioulo. Os objetivos deste estudo foram estimar correlações fenotípicas, genéticas aditivas e ambientais entre caracteres de populações de milho crioulo, cultivados em sistemas de agricultura familiar, para orientar o programa de melhoramento genético participativo. Entre 294 e 400 de progênies de meios-irmãos de seis populações, foram avaliadas em látices triplos 7X7 e 10X10, respectivamente, com parcelas de 4 m de comprimento e espaçadas 1,00 × 0,20 m. Os coeficientes de correlações foram estimados para nove características: produtividade de grãos total (GY), relação peso de grãos por peso de espiga (GE), nº de espigas por planta (EP), nº de dias do florescimento (DF), altura da planta (PH), altura espiga (EH), porcentagem de acamamento (%L), porcentagem de colmos quebrados (%BS) e porcentagem de espigas danificadas (%DE). Houve um alto número de correlação genética aditiva (rA) entre produtividade de grãos e as demais características, nas populações de milho crioulo, ciclos e locais. As maiores médias das estimativas rA para o produtividade de grãos foram com EP (0,67), GE (0,47) e DE% (-0,63). Diferentemente dos trabalhos publicados, foram observadas correlações negativas entre produtividade de grãos e dias para florescimento.

Palavras-chave: Zea mays, população, seleção recorrente, melhoramento genético participativo

Introduction

Brazilian maize landraces populations were originated through Amerindian races, introduced from central and of North America and several commercial materials (Paterniani, 1998; Paterniani et al., 2000; Sánchez et al., 2007). These varieties have been broadly and independently cultivated throughout Brazilian regions and they have relevant socio-economic importance for the family farm systems. As a result, accessions were developed and selected for different environments and morphological characteristics (Paterniani et al., 2000; Warburton et al., 2008; Wietholter et al., 2008). However, more consistent agronomic and genetic knowledge about these collections is still lacking and it is a serious limitation to utilizing, managing, and conserving of maize landrace gene pools (Nass et al., 1993; Andrade et al., 2002; Teixeira et al., 2002).

The use of genotypic correlation helps evaluating the magnitude and direction of associations between characters facilitating the application of indirect selection, because genetic change in a given trait may change other traits, leading to faster and larger genetic gains in maize breeding programs (Vencovsky and Barriga, 1992). Therefore, the selection for another trait may result in indirect response in the low heritable trait, provided the following conditions are satisfied: the genetic correlation between them is substantial, and the heritability of the secondary trait is greater than that of the primary trait (Hallauer and Miranda Filho, 1995). Many authors have estimated genetic correlations among different maize traits in different maize types (Lemos et al., 1992; González et al., 1994; Pinto et al., 2000; Silva et al., 2001; Farias Neto and Miranda Filho, 2001; Saleh et al., 2002; Granate et al., 2002; Daros et al., 2004; Santos et al., 2005; Andrade and Miranda Filho, 2008). However, there is limited published research with genetic correlation estimates for maize landraces. Thus, the objectives of this study were to estimate phenotypic, genotypic and environment correlations among traits in maize landraces, that were cultivated in family farm systems.

Material and Methods

Six maize landraces populations, four with yellow endosperm (Caiano, Carioca, Macaco and Palha Roxa) and two with white endosperm (Cinquentinha and Maizena), that have been improved by the Participatory Breeding Program of Maize Landrace (Ferreira et al., 2006), were used in this study. These landraces were selected individually, through two or three cycles of half-sib family selection, using selection intensity at about 12.0% to 16.3%, according to the methodology showed by Hallauer and Miranda Filho (1995), in family farm systems.

The populations were evaluated with 300 and 400 half-sib progenies in three and four 10X10 lattice experiments, respectively, with three replications each, during the initial cycle of recurrent selection (C0). In the following cycles of the recurrent selection (C1 and C2) 294 half-sib progenies in six 7X7 lattice experiments with three replications (Table 1) were evaluated. Experimental plots were represented by a single row 4 m long, spaced 1 m apart, with 20 plants after thinning. The following traits were taken for analysis: GY - total grain yield corrected for 13.5% moisture (g per plant and t ha1); GE - relation grain weight per ear weight; EP - ears per plant; DF - days to flower, (average between days to silk extrusion and days to anthesis); PH - plant height (cm), EH - ear height (cm); %L - percentage of root loding; %BS - percentage of stalk loding; and %DE - percentage of damaged ears.

The half-sib progenies were obtained and evaluated in farms of Paraná and Santa Catarina states, Brazil, using organic systems of production. The climate of these local is subtropical humid (Table 1).

The basic statistical model for variance and covariance analyses following the lattice design is Yik/j = m + pi + rj + (b|r)k(j) + eik/j, where Yik/j is the observed mean of ith progeny in the kth block of the jth replication; m is the general mean; pi is the random effect of the ith progeny; rj is the random effect of the jth replication; (b|r)k(j) is the random effect of kth block of the jth replication; and eik/j is the error term. The individual variance and covariance analyses were grouped to compute the following estimates for each population, local and cycle of selection: environmental variance ; phenotypic variance of progenies means within experiments

;

additive genetic variance within experiments

;

environmental covariance

; phenotypic covariance of progenies means within experiments

;

genetic covariance among progenies within experiments

;

environmental correlation coefficient

;

phenotypic (means) correlation coefficient

;

additive correlation coefficient

;

where: ME is the mean squares of error; MP/Exp is the mean squares of progeny within experiments; PE is the mean products of error; PP/Exp is the mean products of progeny within experiments. X and Y refer to pairs of traits (Viana and Regazzi, 1999; Silva et al, 2000).

The Student's t test was used to determinate the significance of the phenotypic, additive and environmental correlation coefficients at the 5% and 1% probability levels.

Results and Discussion

A progeny variation (p < 0.01) was detected for total grain yield (GY), days to flower (DF), plant height (PH), ear height (EH) in all experiments. The others traits had significant progeny effects in the majority of experiments. So, there is genetic variability for the populations in the different cycles and local for the vast majority of evaluated traits (Table 2).

The maximum and minimum population average by local and cycles of selection, were, respectively, 8.18 and 1.86 t ha1 for grain yield; 0.85 and 0.62 for relation grain weight per ear weight; 1.12 and 0.62 ears per plant; 50.3% and 5.7% for percentage of damaged ears; 111 and 80 for days to flower; 309 and 171 cm for plant height; 197 and 85 cm for ear height; 24.3% and 0% for percentage of root lodging; and, 53.7% and 0.8% percentage of stalk lodging (Table 2). These results show a good performance of populations in family farm systems for the majority of locals and cycles of selection, when compared to the average of 5.35 t ha1 (density of 50,000 plants ha1) mean estimates for Brazilian populations (Miranda Filho, 1985), particularly when the rustic environmental conditions in which those experiments were conducted are considered. Poor performance for grain yield and other traits were observed in few locations.

The majority of the traits, evaluated in the cycles and local, exhibited additive correlation coefficients (p < 0.01 or p < 0.05) with grain yield. High additive correlations were frequently observed between grain yield and the traits: grain weight per ear weight (GE), ears per plant (EP) and percentage of damaged ears (%DE). High correlations between yield and ears per plant have also been reportedby Lemos et al. (1992), González et al. (1994) and Silva et al. (2001). For the traits days to flower (DF), plant height (PH), ear height (EH), percentage of root loding (%L), percentage of stalk loding (%BS), smaller correlation with grain yield were predominant, but significant for the majority. The high degrees of freedom included in the t test must have contributed to show significant correlation in low estimates. The additive correlation indicated that an increase of grain yield of these maize landraces in family farm systems, is associated with increase of grain weight per ear weight (GE), number of ears per plant (EP) and decrease for percentage of damaged ears (%DE) (Table 3).

For plant height and ear height, frequently, the significant additive correlations with grain yield were positive, ranging from 0.58 to -0.21 and 0.63 to -0.21, and only three cases had low negative estimates (Table 3). Hallauer and Miranda Filho (1995) presented genetic correlation estimates of 0.26 and 0.31 between plant height and ear height with ear yield, respectively, on the average of 23 experiments. Estimates of rA between plant height and ear height with ear yield were, respectively, 0.59 and 0.51 in the base population ESALQ-PB1 (Andrade and Miranda Filho, 2008).

Nineteen estimates of rA between grain yield and days to flower (Table 3) were negatives and significant (-0.52 to -0.10), as long as four were positive and significant (0.43 to 0.13), suggesting that the selection of earlier maturity plants, in these maize landraces populations, could result in greater yield. Average estimate of genetic correlations across 13 populations showed that the correlation between these traits presented low value (0.14) (Hallauer and Miranda Filho, 1995).

Although the estimates of rA between grain yield with percentage of root loding (%L) and with percentage of stalk loding (%BS) were significant, there was no clear pattern, with the estimated values oscillating between positive and negative (Table 3).

The average of environmental (rE), additive genetic (rA) e phenotypic (rF) correlation estimates, for six landraces in the cycles and local (Table 4). For the greater part of traits a concordance in the direction and magnitude of rE, rA e rF estimates were observed. In a small amount of cases that discordance existed, which can be attributed to the environmental modifier effects and physiologic mechanisms, which controlled the expression of the traits (Falconer and Mackay, 1996). In almost every association, the additive genetic correlations presented superior magnitudes when compared to the other types of correlation, except for the correlations of grain yield with days to flower, plant height, ear height and percentage of root loding, ears per plant with days to flower and percentage of root loding.

The average of additive genetic correlations between the trait grain yield and ears per plant was high and positive (0.67), suggesting that the selection of prolific plants resulted in a gain of selection for yield (Table 4). Similar results were shown by Lemos et al. (1992), Santos et al. (2005) and González et al. (1994). For grain yield ´ percentage of damaged ears a high and negative average of additive genetic correlations (-0.63) was observed, permitting the simultaneous selection of plants that have less damaged ears and bigger yields. Daros et al. (2004) reported an estimate of rA between grain yield ´ percentage of damaged equal to -0.11, in a research with pop corn.

Days to flower presented a high positive average of additive genetic correlation with plant height (0.56) and ear height (0.60). These results indicate that selection for earliness would product a reduction in plant and ear heights, in these maize landraces evaluated in family farm systems. High genetic correlations for plant height with days to anthesis (0.67) and days to silk extrusion (0.53) was evidenced by Saleh et al. (2002), in studies using ten corn synthetic populations. However, Santos et al. (2002) reported that phenotypic correlation between plant height and days to anthesis and days to silk extrusion, ear height and days to anthesis and days to silk extrusion, separately, were not significant.

As expected, a high average additive positive genetic correlation between plants and ears heights (0.90) and a low deviation for this trait (0.04) were observed, showing a high association between these traits, over landraces and cycles of selection and local. Similar high estimates were also reported by several authors: 0.89 and 0.95 by Pinto et al. (2000); 0.84 by Silva et al. (2001); 0.78 and 0.80 by Farias Neto and Miranda Filho (2001); 0.97 by Granate et al. (2002); 0.81 and 0.82 by Santos et al. (2005).

The averages of environmental correlations were frequently small, except for the association between plant height and ear height, with coefficient average equal to 0.75. The positive value of environmental correlation demonstrates that both characters were damaged or beneficiated by the same environmental variations and negatives values, which indicate that the environment benefited one of the characters in detriment of the other (Falconer and Mackay, 1996). The highest negative averages of environmental correlations were observed between grain yield and days to flower (-0.31) and grain yield with percentage of damaged ears (-0.35), indicating that the environmental conditions, that lead to increase the days to flower and damaged ears, promote a reduction of yield. This tendency can also be verified through the results showed in Table 2.

The standard errors of the estimates of environmental, additive genetic and phenotypic correlations are shown in Table 4. The standard errors of the estimates of additive genetic correlations between two traits were bigger than the environmental and phenotypic correlation standard deviations, except between plant height and ear height deviation, when standard deviation of the environmental correlation was higher.

Through the correlations and their standard errors of the estimates it was possible to observe that some associations between traits always presented the same positive (grain yield with relation grain weight per ear weight and ears per plant; percentage of damaged ears with days to flower; days to flower with plant height and ear height; and, plant height with ear height) and negative sign (grain yield with percentage of damaged ears; relation grain weight per ear weight with percentage of damaged ears and days to flower).

The selection in these maize landraces to increase yield is related to the selection of plants more prolific, earlier plants, smaller percentage of damaged ears and higher relation grain weight per ear weight. There is a positive association between taller plants and increased number of days to flower, but with magnitude values that allow the selection of plants with higher yield, small height and number of days to flower. With the exception of the trait days to flower, the correlations observed in landraces within family farm systems are the same presented in several other reports (Lemos et al., 1992, González et al., 1994; Hallauer and Miranda Filho, 1995; Farias Neto and Miranda Filho, 2001; Silva et al., 2001; Saleh et al., 2002; Daros et al., 2004; Santos et al., 2005; Andrade and Miranda Filho, 2008), but with different magnitudes in some cases.

Acknowledgements

To Fundação Araucária for financial support and to AS-PTA and the family farmers that work with us.

Received December 19, 2009

Accepted January 05, 2010

  • Andrade, R.V.; Santos, M.X.; Ferreira, A.S.; Oliveira, A.C. 2002. Evaluation of maize land races accession collected in the Central region of Brazil. Revista Brasileira de Milho e Sorgo 1: 67-74. (in Portuguese, with abstract in English).
  • Andrade, J.A.C.; Miranda Filho, J.B. 2008. Quantitative variation in the tropical maize population, ESALQ-PB1. Scientia Agricola 65: 174-182.
  • Daros, M.; Amaral Júnior, A.T.; Pereira, M.; Santos, F.S.; Scapim, C.A.; Freitas Júnior, S.P.; Daher, R.F.; Ávila, M.R. 2004. Correlations among agronomic traits in two recurrent selection cycles in popcorn. Ciência Rural 34:1389-1394. (in Portuguese, with abstract in English).
  • Falconer, D.S; Mackay, T.F.C. 1996. Introduction to Quantitative Genetics. 4ed. Longmans Green, Harlow, Essex, UK.
  • Farias Neto, A.L.; Miranda Filho, J.B. 2001.Genetic correlation between traits in the ESALQ-PB1 maize population divergently selected for tassel size and ear height. Scientia Agricola 58: 119-123.
  • Ferreira, J.M.; Moreira, R.M.P.; Tardin, J.M.; Jantara, A.E. 2006. Participatory breeding landraces of maize: an experience in the Southern Central Region of Paraná State, Brazil p. 197-219. In: Labrada, H.R.; Espinosa, M.M.H.; Sotomayor, J.C.R., eds. Participatory breeding: farmers improve crop. Ediciones INCA, La Habana, Cuba. (in Spanish).
  • González, P.A.; Lemos, M.A.; Ramalho Neto, C.E.; Reis, O.V.; Tabosa J.N.; Tavares, Filho J.J. 1994. Genetic, phenotypic and environmental in two cycles of selection in the dent corn composite. Pesquisa Agropecuária Brasileira 29: 419-425. (in Portuguese, with abstract in English).
  • Granate, M.J.; Cruz, C.D.; Pacheco, C.A.P. 2002. Prediction of genetic gains from selection on popcorn half-sib families of CMS 43 population. Ciência Agrotecnologica 26: 1228-1235. (in Portuguese, with abstract in English).
  • Hallauer, A.R.; Miranda Filho, J.B. 1995. Quantitative Genetics in Maize Breeding. Iowa State University Press, Ames, IA, USA.
  • Lemos, M.A; Gama, E.E.G.; Oliveira, A.C.; Araújo, M.R.A. 1992. Correlations genotypic, phenotypic and environmental progenies of maize. Pesquisa Agropecuária Brasileira 27: 1563-1569. (in Portuguese, with abstract in English).
  • Miranda Filho, J.B. 1985. Breeding methodologies for tropical maize. p. 177-206. In: Brandolini A.; Salamini F., eds. Breeding strategies for maize production improvement in the tropics. FAO, Rome, Italy.
  • Nass, L.L.; Pellicano, I.J.; Valois, A.C.C. 1993. Utilization of genetic resources for maize and soybean breeding in Brazil. Brazilian Journal of Genetics 16: 983-988.
  • Paterniani, E. 1998. Genetic diversity and races of maize in Brazil. p. 28-31. In: Soares, A.C.; Machado A.C.; Silva B.M.; Weid J.M., eds. Maize landrace: conservation and use of biodiversity. AS-PTA, Rio de Janeiro, RJ, Brazil. (in Portuguese).
  • Paterniani, E.; Nass, L.L.; Santos, M.X. 2000. The value of genetic resources of maize in Brazil. p. 11-42. In: Udry C.V.; Duarte W., eds. A history of Brazilian maize: the value of genetic resources. Paralelo15, Brasília, DF, Brazil. (in Portuguese).
  • Pinto, R.M.C; Lima, Neto, F.P.; Souza Júnior, C.L. 2000. Estimates of appropriate number of S1 progenies for recurrent selection in maize. Pesquisa Agropecuária Brasileira 35: 63-73. (in Portuguese, with abstract in English).
  • Saleh, G.B; Alawi, S.A.S.; Panjaitan, K. 2002. Performance, correlation and herdability studies on selected sweet corn synthetic populations. Pakistan Journal of Biological Sciences 5: 251-254.
  • Sánchez, G.J.J.; Goodman, M.M.; Stuber, C.W. 2007. Racial diversity of maize in Brazil and adjacent areas. Maydica 52: 13-30.
  • Santos, P.G.; Juliatti, F.C.; Buiatti, A.L.; Hamawaki, O.T. 2002. Evaluation of the agronomic performance of corn hybrids in Uberlândia, MG. Pesquisa Agropecuária Brasileira 37: 597-602. (in Portuguese, with abstract in English).
  • Santos, M.F.; Moro, G.V.; Aguiar, A.M.; Souza Júnior, C.L. 2005. Responses to reciprocal recurrent selection and changes in genetic variability in IG-1 and IG-2 maize populations. Genetics and Molecular Biology 28: 781-788.
  • Silva, E.S.; Silva, P.S.L.; Nunes, G.H.S; Silva, K.M.B. 2001. Estimation of genetic parameters in the compound of maize ESAM-1. Caatinga 14: 43-52. (in Portuguese, with abstract in English).
  • Silva, H.D.; Ferreira, D.F.; Pacheco, C.A.P. 2000. Evaluation of four alternative analysis of experiments in square lattice, with the estimation of variance components. Bragantia 59: 117-123. (in Portuguese, with abstract in English).
  • Teixeira, F.F.; Andrade, R.V.; Oliveira, A.C.; Ferreira, A.S.; Santos, M.X. 2002. Diversity in maize germplasm collected in northeastern Brazil. Revista Brasileira de. Milho Sorgo 1: 59-67 (in Portuguese, with abstract in English).
  • Vencovsky, R.; Barriga, P. 1992. Genetics biometrics in plant breeding. Sociedade Brasileira de Genética, Ribeirão Preto, SP, Brazil. (in Portuguese)
  • Viana, J.M.S.; Regazzi, A.J. 1999. Estimation of genetic parameters in the square lattice analysis. Bragantia 58: 185-193.
  • Warburton, M.L.; Reif, J.C.; Frisch, M.; Bohn, M.; Bedoya, C.; Xia, X.C.; Crossa, J.; Franco, J.; Hoisington, D.; Pixley, K.; Taba, S.; Melchinger,. A.E. 2008. Genetic diversity in CIMMYT nontemperate maize germoplasma: landraces, open pollinated varieties, and inbred lines. Crop Science 48: 617-624.
  • Wietholter, P.; Sereno, M.J.M.S.; Terra, T.F.; Silva, S.D.A.; Barbosa Neto, J.F. 2008. Genetic variability in corn landraces from Southern Brazil. Maydica 53: 151-159.
  • *
    Corresponding author <
  • Publication Dates

    • Publication in this collection
      23 May 2011
    • Date of issue
      Dec 2010

    History

    • Received
      19 Dec 2009
    • Accepted
      05 Jan 2010
    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