Acessibilidade / Reportar erro

Genetic correlation between traits in the ESALQ-PB1 maize population divergently selected for tassel size and ear height

Correlação genética entre caracteres na população de milho ESALQ-PB1 após seleção divergente para tamanho do pendão e altura da espiga

Abstracts

Full-sib and selfed (S1) progenies were obtained from sub-populations of ESALQ-PB1, divergently selected for tassel size (T+ and T-) and ear height (E+ and E-), and used for estimating genetic and phenotypic correlation coefficients between traits. The analyzed traits were: EW- total ear weight (g/plant), PH- plant height (cm), EH- ear height (cm), TB- tassel branch number and TL- tassel length. The highest genetic (rG) and phenotypic (rF) correlation was observed for the combination PH x EH, as expected, with average of 0.800 and 0.778, respectively over sub-populations and locations. It is apparent that divergent selection for tassel size did not affect greatly the correlation between PH and EH in the full sib progenies, but in the inbred progenies the correlation was smaller in the sub-population selected for larger tassels. Genetic correlation between PH and EH with tassel traits was always positive but ranged from 0.020 to 0.668 in Piracicaba and from 0.06 to 0.309 in Rio Verde. Genetic correlation between PH and EH with yield (EW) also was positive in the range of 0.087 to 0.503. EH showed higher correlation with EW in relation to PH x EW and differences were larger in the sub-populations divergently selected for ear height. Correlation between tassel traits with other traits was positive in most of instances and a lack of consistency was observed among sub-populations. Generally the coefficients of genetic and phenotypic correlation differed substantially from the estimates in the base population ESALQ-PB1 before divergent selection for tassel size and ear placement. Divergent selection affected the correlation between traits under unpredicted and varying magnitudes.

Zea mays; genetic correlation; divergent selection; tassel size; ear height


Progênies de irmãos germanos e de autofecundação (S1) foram obtidas de sub-populações de ESALQ-PB1 após seleção divergente para tamanho do pendão (T+ e T-) e altura da espiga (E+ e E-) e utilizadas para estimar os coeficientes de correlação genética e fenotípica entre caracteres. Foram analisados os caracteres: EW- peso de espigas (g/planta), PH- altura da planta (cm), EH- altura da espiga (cm), TB- número de ramificações do pendão e TL- comprimento do pendão. Como esperado, as maiores correlações genéticas (rG) e fenotípicas (rF) foram observadas para a combinação PH x EH, com médias de 0,800 e 0,778 , respectivamente, em várias sub-populações e locais. Aparentemente, a seleção divergente para tamanho do pendão não afetou fortemente a correlação entre PH e EH nas progênies de irmãos germanos, mas nas progênies endogâmicas a correlação foi menor na sub-população selecionada para pendão grande. A correlação genética entre PH e EH com os caracteres do pendão foram sempre positivas mas variaram de 0,020 a 0,068 em Piracicaba e de 0,060 a 0,309 em Rio Verde. A correlação genética entre PH e EH com produção (EW) também foi positiva, na faixa de 0,087 a 0,503. EH mostrou maior correlação com EW em relação a PH x EW e as diferenças foram maiores nas sub-populações que sofreram seleção divergente para altura da espiga. A correlação entre os caracteres do pendão com outros caracteres foi positiva na maioria dos casos e uma falta de consistência foi observada entre sub-populações. Em geral, os coeficientes de correlação genética e fenotípica diferiram substancialmente das estimativas na população base ESALQ-PB1 antes da seleção divergente para tamanho do pendão e posição da espiga. A seleção divergente afetou a correlação entre caracteres sob magnitudes imprevisíveis e variáveis.

Zea mays; correlação genética; seleção divergente; tamanho do pendão; altura da planta


GENETIC CORRELATION BETWEEN TRAITS IN THE ESALQ-PB1 MAIZE POPULATION DIVERGENTLY SELECTED FOR TASSEL SIZE AND EAR HEIGHT

Austeclínio Lopes Farias Neto1; José Branco de Miranda Filho2*

1Embrapa Cerrados, Rodovia BR 020, km 17 - CEP: 73301-970 - Planaltina, DF.

2Depto. de Genética - USP/ESALQ, C.P. 83 - CEP: 13400-970 - Piracicaba, SP.

*Autor correspondente <jbmirand@carpa.ciagri.usp.br>

ABSTRACT: Full-sib and selfed (S1) progenies were obtained from sub-populations of ESALQ-PB1, divergently selected for tassel size (T+ and T–) and ear height (E+ and E–), and used for estimating genetic and phenotypic correlation coefficients between traits. The analyzed traits were: EW- total ear weight (g/plant), PH- plant height (cm), EH- ear height (cm), TB- tassel branch number and TL- tassel length. The highest genetic (rG) and phenotypic (rF) correlation was observed for the combination PH x EH, as expected, with average of 0.800 and 0.778, respectively over sub-populations and locations. It is apparent that divergent selection for tassel size did not affect greatly the correlation between PH and EH in the full sib progenies, but in the inbred progenies the correlation was smaller in the sub-population selected for larger tassels. Genetic correlation between PH and EH with tassel traits was always positive but ranged from 0.020 to 0.668 in Piracicaba and from 0.06 to 0.309 in Rio Verde. Genetic correlation between PH and EH with yield (EW) also was positive in the range of 0.087 to 0.503. EH showed higher correlation with EW in relation to PH x EW and differences were larger in the sub-populations divergently selected for ear height. Correlation between tassel traits with other traits was positive in most of instances and a lack of consistency was observed among sub-populations. Generally the coefficients of genetic and phenotypic correlation differed substantially from the estimates in the base population ESALQ-PB1 before divergent selection for tassel size and ear placement. Divergent selection affected the correlation between traits under unpredicted and varying magnitudes.

Key words:Zea mays, genetic correlation, divergent selection, tassel size, ear height

CORRELAÇÃO GENÉTICA ENTRE CARACTERES NA POPULAÇÃO DE MILHO ESALQ-PB1 APÓS SELEÇÃO DIVERGENTE PARA TAMANHO DO PENDÃO E ALTURA DA ESPIGA

RESUMO: Progênies de irmãos germanos e de autofecundação (S1) foram obtidas de sub-populações de ESALQ-PB1 após seleção divergente para tamanho do pendão (T+ e T–) e altura da espiga (E+ e E–) e utilizadas para estimar os coeficientes de correlação genética e fenotípica entre caracteres. Foram analisados os caracteres: EW- peso de espigas (g/planta), PH- altura da planta (cm), EH- altura da espiga (cm), TB- número de ramificações do pendão e TL- comprimento do pendão. Como esperado, as maiores correlações genéticas (rG) e fenotípicas (rF) foram observadas para a combinação PH x EH, com médias de 0,800 e 0,778 , respectivamente, em várias sub-populações e locais. Aparentemente, a seleção divergente para tamanho do pendão não afetou fortemente a correlação entre PH e EH nas progênies de irmãos germanos, mas nas progênies endogâmicas a correlação foi menor na sub-população selecionada para pendão grande. A correlação genética entre PH e EH com os caracteres do pendão foram sempre positivas mas variaram de 0,020 a 0,068 em Piracicaba e de 0,060 a 0,309 em Rio Verde. A correlação genética entre PH e EH com produção (EW) também foi positiva, na faixa de 0,087 a 0,503. EH mostrou maior correlação com EW em relação a PH x EW e as diferenças foram maiores nas sub-populações que sofreram seleção divergente para altura da espiga. A correlação entre os caracteres do pendão com outros caracteres foi positiva na maioria dos casos e uma falta de consistência foi observada entre sub-populações. Em geral, os coeficientes de correlação genética e fenotípica diferiram substancialmente das estimativas na população base ESALQ-PB1 antes da seleção divergente para tamanho do pendão e posição da espiga. A seleção divergente afetou a correlação entre caracteres sob magnitudes imprevisíveis e variáveis.

Palavras-chave:Zea mays,correlação genética, seleção divergente, tamanho do pendão, altura da planta

INTRODUCTION

Maize (Zea mays L.) is, among the cultivated species, one that has undergone a fantastic improvement in yield through artificial selection (Hallauer & Miranda Filho, 1995). Cultural practices also have greatly contributed to increase yield in the maize crops during this century. More recently, attention has been given to plant architecture toward a better efficiency in the utilization of solar energy and photosynthetic capacity (Donald, 1968; Wittner, 1974; Paterniani, 1974, 1981; Miranda Filho, 1974; Sampaio, 1986). Also, intentionally or not, changes have occurred in tassel size for a better utilization of photosynthates (Mock & Pearce, 1975).

The association between traits is an important aspect to deal with in breeding programs, because genetic change in a given trait may change positively or negatively other traits (Vencovsky & Barriga, 1992). In addition, in most breeding programs the strategy is based on selection for several traits simultaneously and, therefore, knowledge on the genetic association between traits is inevitably useful for the establishment of selection criteria. The basic causes of genetic correlation are pleiotropy, and linkage disequilibrium (Falconer, 1964; Vencovsky, 1978; Hallauer & Miranda Filho, 1995).

One important application of the genetic correlation in breeding programs refers to indirect selection for traits of low heritability and consequently low direct response to selection. Selection for another trait may result in indirect response in the low heritable trait, provided the following conditions are satisfied: i) traits under consideration must be highly correlated genetically; ii) heritability of the secondary trait must be higher than the trait of higher interest (Falconer, 1964; Vencovsky & Barriga, 1992). In maize, for example, two secondary traits were identified for indirect selection to increase yield: prolificacy (Paterniani, 1981) and tassel size (Geraldi et al., 1985).

Although many authors have referred to a negative association between tassel size and yield potential (Leonard & Kiesselbach, 1932; Hunter et al., 1973; Fakorede & Mock, 1978; Lordelo & Miranda Filho, 1981; Geraldi et al., 1985), positive association of those traits also have been reported (Lordelo & Miranda Filho, 1981; Sampaio, 1986; Soares Filho, 1987). There are evidences that the correlation between yield and tassel size tends to be higher and negative under stresses caused by unfavorable environments (Soares Filho, 1987). Brunini et al. (1983) emphasized that environmental factors such as photoperiod, solar radiation and rainfall affect decisively the yield potential of a corn crop and consequently the association between traits may change if there is a differentiated variety response to the environmental factors.

Souza Júnior et al. (1985) reported on a negative correlation between tassel size and prolificacy which were explained by a larger amount of indol-acetic-acid (IAA) produced by larger tassels and causing inhibition of prolificacy, or vice-versa (Anderson, 1967). Other studies referring to the correlation of tassel size with other characters are sometimes contradictories. The correlation between tassel branch number with plant height, ear height and ear placement was reported as positive by Obilana & Hallauer (1974), Ayala Osuna et al. (1986); Miranda Filho & Andrade (2000). However, the same associations were positive in some cases and negative in other cases, as reported by Lordelo (1982), Aguilar Morán (1984), Sampaio (1986), Soares Filho (1987), and Araújo (1992). Some estimates of additive and phenotypic correlation between tassel branch number with plant height, ear height, and grain yield are summarized in TABLE 1. The objectives of this work was to study changes in the genetic correlation between traits after divergent selection for tassel size and ear height in the maize population ESALQ-PB1.

MATERIAL AND METHODS

Four sub-populations were derived from ESALQ-PB1, after six cycles of divergent selection for tassel size and ear height. Sub-populations were designated by

ESALQ-PB1 (E+): selection for increasing ear height end/or ear placement

ESALQ-PB1 (E–): selection for decreasing ear height end/or ear placement

ESALQ-PB1 (T+): selection for increasing tassel size

ESALQ-PB1 (T–): selection for decreasing tassel size

Full sib progenies were obtained through biparental crosses from each sub-population. In ESALQ-PB1 (T+) and ESALQ-PB1 (T–), S1 progenies were also obtained by selfing individual plants. The six sets totaled 921 progenies that were evaluated in 18 experiments following the completely randomized block design with three replications, at Piracicaba, SP, and Rio Verde, GO. The numbers of progenies of each set in each experiment are shown in TABLE 2.

Experimental plots were represented by a single row 4.0 m long with 20 plants per plot (0.20 m between plants) after thinning. The following traits were analyzed: PH-plant height (cm), EH-ear height (cm), TB-tassel branch number, TL-tassel length (cm), and EW-total ear weight (g/plant). Except for EW, the experimental units were represented by means of three plants per plot. In Rio Verde (GO), only PH, EH and TB of full sib progenies were evaluated.

The statistical model for variance and covariance analyses following the randomized complete block design is Yij = m + pi + bj + eij , where Yij is the observed mean of the ith progeny in the jth replication, m is the general mean, pi is the random effect of the ith progeny, bj is the random effect of the jth replication, and eij is the error term. In the analysis of variance the mean squares for progenies (MP) and Error (ME) have the following expectations: E(MP) = s2 + 3 and E(ME) = s2 . In the same way, the corresponding expected mean products (PP and PE) in the analysis of covariance have the expectations E(PP) = cov + 3 covP and E(PE) = cov. In the formulations, s2 and cov refer to the variance and covariance of effects relative to the error term. In the same way, and covP refer to the genetic variance and covariance relative to progeny effects. The coefficients of correlation were then calculated by

Genetic correlation coefficient -------- rG =

Phenotypic correlation of progeny means - rF =

X and Y refer to pairs of characters

RESULTS AND DISCUSSION

Observed means for six traits in sub-populations of ESALQ-PB1, divergently selected for tassel size (T+ and T–) and ear height (E+ and E–) are shown in TABLE 3. The direct effect of divergent selection on means of selected traits and the correlated response on other non-selected traits can be visualized in TABLE 3. Also, Farias Neto & Miranda Filho (2000) already presented comparisons involving inbred and non-inbred progenies, leading to estimates of inbreeding depression. Therefore, our attention here focuses on results of genetic (rG ) and phenotypic (rF) correlation (TABLE 4).

Correlation between plant height (PH) and ear height (EH)

The correlation coefficients (rG and rF) between PH and EH were the highest among the combinations of traits in this study. Estimates of rG and rF averaged 0.800 and 0.778, ranging from 0.755 to 0.870 and 0.675 to 0.911, respectively, over sub-populations and locations. It is apparent that divergent selection for tassel size did not affect greatly the correlation between PH and EH in the full sib progenies, but in the inbred progenies the correlation was smaller in the sub-population selected for larger tassels. Also, rG and rF were smaller in the sub-population selected for decreasing ear height. In general, a high correlation between PH and EH has been reported by many authors. Estimates of the additive genetic correlation (rA) and phenotypic correlation summarized by Moraes (1989) averaged 0.82 and 0.78, respectively. Also, Hallauer & Miranda Filho (1995) presented a high average genetic correlation of the order of 0.81. Miranda Filho & Andrade (2000) reported on estimates of rA and rF of 0,842 and 0.803 between PH and EH in the base population ESALQ-PB1.

Correlation between plant height (PH) and ear height (EH) with tassel traits (TB, and TL)

The correlation coefficients between PH and EH with tassel branch number were always positive but not consistent among sub-populations. The highest correlations of EH with TB were observed with inbred progenies. Also, for both inbred and non-inbred progenies in Piracicaba the genetic correlation was higher in the sub-population selected for larger tassels. However, in Rio Verde, the highest positive correlation (0.301) was for EH x TB in the sub-population (T–). In general, correlation coefficients between EH and TB were higher than for PH and TB. Estimates of rA and rF in the base population ESALQ-PB1 were 0.379 and 0.339 for PH x TB and 0.436 and 0.382 for EH x TB, respectively (Miranda Filho & Andrade, 2000). The correlation of PH and EH with tassel length (TL) was estimated in only three and two sub-populations, respectively, and the estimates of rG were all positive and above 0.40, except for PH x TL in the sub-population ESALQ-PB1 (T+) with inbred progenies (rG = 0.127). Estimates of rA in the base population ESALQ-PB1 were 0.273 and 0.101 for PH x TL and EH x TB, respectively (Miranda Filho & Andrade, 2000).

Correlation between plant height (PH) and ear height (EH) with ear yield (EW)

The association of plant height and ear height with yield tends to be always positive. Hallauer & Miranda Filho (1995) reported genetic correlation of 0.26 and 0.31 for the combinations PH x EW and EH x EW, respectively, on the average of 23 experiments. The association of these traits in our experiments comprised only non-inbred progenies, where rG varied from 0.087 to 0.503 among sub-populations. EH showed higher correlation with EW than PH, and differences were larger in the sub-populations divergently selected for ear height. When looking to the sub-populations divergently selected for tassel size, a higher genetic correlation of PH and EH with EW was observed in ESALQ-PB1 (T+). Estimates of rA in the base population ESALQ-PB1 were 0.594 and 0.506 for PH x EW and EH x EW, respectively (Miranda Filho & Andrade, 2000).

Correlation between tassel traits (TB, TL) and EW

The genetic correlation coefficients between TB and TL were low but negative (rG = -0.015) in ESALQ-PB1 (T+) and positive (rG = 0.232) in ESALQ-PB1 (T–), the two sub-populations representing the divergent selection for tassel size. The later was closer to the rA estimate (0.217) in the base population ESALQ-PB1 (Miranda Filho & Andrade, 2000). The genetic correlation TBxEW was calculated only for full-sib progenies and was small negative (-0.022) in ESALQ-PB1 (T+) and positive (0.137 to 0.456) in the other three sub-populations. In a study with the same sub-populations in the preceding cycle of divergent selection, Araújo (1992) found the same pattern of correlation between TB and EW. The only rG estimate for TL x EW was in the sub-population ESALQ-PB1 (T+) and was small positive (0.027). Apparently, divergent selection for tassel size and ear height affected the original pattern of correlation between tassel traits and ear yield because all the estimates depart from the rA estimates of 0.207 and 0.216 for TB x EW and TL x EW, respectively, as reported by Miranda Filho & Andrade (2000) for the base population ESALQ-PB1. Falconer (1964) emphasized that selection may change significantly the pattern of correlation between traits. However, the pattern of correlation between tassel size and yield depends also on the effect of environment and tends to be positive under favorable environments (Soares Filho, 1987).

Generally the coefficients of genetic and phenotypic correlation differed substantially from the estimates in the base population ESALQ-PB1 before divergent selection for tassel size and ear placement (Miranda Filho & Andrade,2000). Divergent selection affected the correlation between traits under unpredicted and varying magnitudes.

MIRANDA FILHO, J.B.; ANDRADE, J.A.C. Quantitative variation in the maize population ESALQ-PB1. Piracicaba: USP/ESALQ. Depto de Genética, 2000. n.p. /não publicado/

Received October 13, 1999

  • AGUILAR MORÁN, J.F. Avaliaçăo do potencial genético de linhagens e respectivos testadores obtidos de duas populaçőes de milho (Zea mays L.). Piracicaba, 1984. 118p. Dissertaçăo (Mestrado) - Escola Superior de Agricultura Luiz de Queiroz, Universidade de Săo Paulo.
  • ANDERSON, I.C. Plant characteristics that affect yield. In: ANNUAL HYBRID CORN INDUSTRY-RESEARCH CONFERENCE, 22., Washington, 1967. Proceedings Washington, 1967. p.71-73.
  • ARAÚJO, P.M. Variabilidade genética em sub-populaçőes de milho (Zea mays L.) obtidas por seleçăo divergente. Piracicaba, 1992. 153p. Dissertaçăo (Mestrado) - Escola Superior de Agricultura Luiz de Queiroz, Universidade de Săo Paulo.
  • AYALA OSUNA, J.; ARAÚJO, S.M.C.; BANZATO, D.A. Avaliaçăo do tamanho do pendăo em relaçăo a alguns caracteres agronômicos em duas populaçőes de milho. In: CONGRESSO NACIONAL DE MILHO E SORGO, 15., Maceió, 1984. Anais Brasília: EMBRAPA, CNPMS, 1986.
  • BRUNINI, O.; ALFONSI, R.R.; PAES DE CAMARGO, M.B. Efeito dos elementos climáticos no desenvolvimento da cultura do milho. In: SIMPÓSIO SOBRE A PRODUTIVIDADE DO MILHO, Londrina, 1983. Anais Londrina: IAPAR, 1983.
  • DONALD, C.M. The breeding of crop ideotypes. Euphytica, v.17, p.385-403, 1968.
  • FAKOREDE, M.A.B.; MOCK, J.J. Changes in morphological and physiological traits associated with recurrent selection for grain yield in maize. Euphytica, v.27, p.397-409, 1978.
  • FALCONER, D.S. Introduction to quantitative genetics 2.ed. New York: Ronald Press, 1964. 365p.
  • FARIAS NETO, A.L.; MIRANDA FILHO, J.B. Inbreeding in sub-populations of maize divergently selected for tassel size and ear height. Scientia Agricola, v.57, p.487-490, 2000.
  • GERALDI, I.O.; MIRANDA FILHO, J.B.; VENCOVSKY, R. Estimates of genetic parameters for tassel characters in maize (Zea mays L.) and breeding perspectives. Maydica, v.30, p.1-14, 1985.
  • HALLAUER, A.R.; MIRANDA FILHO, J.B. Quantitative genetics in maize breeding 2.ed. Ames: Iowa State University Press, 1995. 468p.
  • HUNTER, R.B.; DAYNARD, T.B.; HUME, D.J.; TANNER, J.W.; CURTIS, J.O; KANNEMBERG, L.W. Effect of tassel removal on grain yield of corn (Zea mays L.). Crop Science, v.9, p.405-406, 1973.
  • LEONARD, W.A.; KIESSELBACH, T.A. The effect of the removal of tassels on the yield of corn. Agronomy Journal, v.24, p.514-516, 1932.
  • LORDELO, J.A.C. Parâmetros genéticos nas populaçőes de milho Piranăo-VD2 e Piranăo-VF1 Piracicaba: ESALQ/USP, 1982.
  • LORDELO, J.A.C.; MIRANDA FILHO, J.B. Correlaçőes genéticas e fenotípicas entre caracteres em duas populaçőes de milho braquítico. Relatório Científico do Departamento de Genética, ESALQ/USP, v.15, p.104-108, 1981.
  • MIRANDA FILHO, J.B. Cruzamentos dialélicos e síntese de compostos de milho (Zea mays L) com ęnfase na produtividade e no porte da planta. Piracicaba, 1974. p. 164. Tese (Doutorado) - Escola Superior de Agricultura Luiz de Queiroz, Universidade de Săo Paulo.
  • MOCK, J.J.; PEARCE, R.B. An ideotype of maize. Euphytica, v.24, p.613-623, 1975.
  • MORAES, M.L.T. Estimativas de parâmetros genéticos para os caracteres: produçăo de grăos, altura da planta e altura da espiga em populaçőes de milho. Piracicaba: ESALQ/USP, Depto. de Genética, 1989. 40p.
  • OBILANA, A.T.; HALLAUER, A.R. Estimation of variability of quantitative traits in BSSS by using unselected maize inbred lines. Crop Science, v.14, p. 99-103, 1974.
  • PATERNIANI, E. Origem e comportamento do milho Piranăo. Relatório Científico do Departamento de Genética, ESALQ/USP, v.8, p.104-108, 1974.
  • PATERNIANI, E. Influence of tassel size on ear placement in maize (Zea mays L.). Maydica, v.26, p.85-91, 1981.
  • SAMPAIO, N.F. Propriedades genéticas e potencial para o melhoramento dos compostos de milho (Zea mays L) ESALQ-PB4 e ESALQ-PB5. Piracicaba, 1986. 105p. Dissertaçăo (Mestrado) - Escola Superior de Agricultura Luiz de Queiroz, Universidade de Săo Paulo.
  • SOARES FILHO, W.S. Características fenotípicas e genéticas das populaçőes de milho (Zea mays L.) braquítico Piranăo-VD2B e Piranăo-VF1B. Piracicaba, 1987. 185p. Tese (Doutorado) - Escola Superior de Agricultura Luiz de Queiroz, Universidade de Săo Paulo.
  • SOUZA JR., C.L.; GERALDI, I.O; ZINSLY, J.R. Influence of tassel size on the expression of prolificacy in maize (Zea mays L.) Maydica, v.30, p.321-328, 1985.
  • VENCOVSKY, R. Herança quantitativa. In: PATERNIANI, E.; VIÉGAS, P. (Ed.) Melhoramento e produçăo do milho. Campinas: Fundaçăo Cargill, 1978.
  • VENCOVSKY, R.; BARRIGA, P. Genética biométrica no fitomelhoramento Ribeirăo Preto: SBG, 1992. 496p.
  • WITTNER, S.H. Maximum production capacity of food crops. Bioscience, v.24, p.216-224, 1974.

Publication Dates

  • Publication in this collection
    09 Feb 2001
  • Date of issue
    Mar 2001

History

  • Received
    13 Oct 1999
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