Análise da variabilidade genética aditiva de características de crescimento na raça Nelore

Additive genetic variability analysis in the growth characteristics of Nellore breed

Análise da variabilidade genética aditiva de características de crescimento na raça Nelore¹

Additive genetic variability analysis in the growth characteristics of Nellore breed

Roberta Lisboa Pontes Gestal de Siqueira^I; João Ademir de Oliveira^II; Raysildo Barbosa Lôbo^III; Luiz Antônio Framartino Bezerra^IV; Humberto Tonhati^V

^IDoutoranda, FMRP-USP/Professora e Coordenadora do Núcleo de Melhoramento Genético-UNIRP, São José do Rio Preto. E.mail: [e-mail]

^IIDepartamento de Ciências Exatas-FCAV-UNESP, Jaboticabal. E.mail: [e-mail]

^IIIProfessor Associado, Departamento de Genética e Matemática Aplicada à Biologia-FMRP-USP, Ribeirão Preto. E.mail: [e-mail]

^IVAnalista de Sistemas, Departamento de Genética e Matemática Aplicada à Biologia-FMRP-USP, Ribeirão Preto. E.mail: [e-mail]

^VDepartamento de Zootecnia - FCAV-UNESP, Jaboticabal. E.mail: [e-mail]

RESUMO

Foram utilizados dados de cinqüenta e um rebanhos participantes do Programa de Melhoramento Genético da Raça Nelore (PMGRN), distribuídos nos estados de Goiás (GO), Mato Grosso do Sul (MS), Mato Grosso (MT), Minas Gerais (MG), São Paulo (SP), Maranhão (MA) e Bahia (BA). Foram obtidas estimativas de parâmetros genéticos para os pesos padronizados aos 120 (P120), 455 (P455) e 550 (P550) dias de idade. Análises unicaráter e bicaráter foram realizadas por modelo animal usando o aplicativo MTDFREML. Para P120 foi utilizado um modelo que incluiu como efeitos fixos, grupo de contemporâneos e classe de idade da vaca ao parto, e como aleatórios, os efeitos genéticos direto, materno e de ambiente permanente da vaca. Para P455 e P550, o modelo utilizado incluiu os mesmos efeitos fixos e o efeito genético direto do animal. ANas análises unicaráter, as estimativas de herdabilidade direta foram 0,29, 0,51 e 0,47 para P120, P455 e P550, respectivamente. Nas análises bicaráter, observaram–se coeficientes de herdabilidade direta de 0,50 e 0,58 para P120, 0,50 e 0,53 para P455 e 0,44 e 0,49 para P550. As correlações genéticas estimadas entre P120 e P455, P120 e P550 e P455 e P550, foram 0,92, 0,93 e 0,96, respectivamente. As estimativas de herdabilidade obtidas para P455 e as correlações genéticas deste peso com P120 e P550 sugerem que a avaliação genética pode ser feita aos 15 meses de idade em substituição aos 18 meses.

Palavras-chave: avaliação genética, bovinos, características de crescimento, Nelore, parâmetros genéticos

ABSTRACT

The data were obtained from 51 herds to participate in the Nelore Catttle Breeding Program (NCBP) from the states of Goiás (GO), Mato Grosso do Sul (MS), Mato Grosso (MT), Minas Gerais (MG), São Paulo (SP), Maranhão (MA) and Bahia (BA). Were used to estimative genetic parameters for standardized weights at 120 (P120), 455 (P455) and 550 (550) days of age. Univariate and bivariate analysis were performed by animal model using MTDFREML program. For P120 was used a model that included contemporary groups and cow age at calving as fixed effects, and direct genetic, maternal genetic and permanent environment effects as random effects. For P455 and P550 were utilized the same model but without maternal direct and permanent environment effects. The estimates of heritability direct from univariate analysis were: 0.29, 0.51 and 0.47 for P120, P455 and P550, respectively. In the bivariate analyses the direct heritability values were of high magnitude. The genetic correlation between P120 and P455, P120 and P550 and P455 and P550 were 0.92, 0.93 and 0.96, respectively. The values of the heritability coefficients estimated for the trait P455 and genetic correlation that characteristic with others indicate that the genetic evaluation could be made at the 15 months of age.

Keywords: bovine, genetic evaluation, genetic parameters, growth traits, Nellore

Introdução

A pecuária de corte no Brasil está experimentando importantes transformações, ressaltando-se a crescente busca dos produtores por tecnologias que possam resultar no incremento da rentabilidade do sistema de exploração. As fazendas estão se tornando empresas, contando com equipes de trabalho e assessorias nas áreas de melhoramento genético, nutrição, manejo e sistemas de criação, com introdução da informática no meio rural.

Até pouco tempo, o produtor de gado de corte além de abater os animais em idades tardias e com pesos elevados, ofertava carne de baixa qualidade ao mercado. Ademais, não percebia o custo que isso representava. Atualmente, a busca pela precocidade tornou-se uma prioridade tanto no meio rural como no científico. Este conceito tem sido muito discutido e vem gerando projetos de pesquisa com o objetivo de identificar critérios de seleção para obtenção de animais precoces, favorecendo a rentabilidade do sistema, já que reduz o tempo de permanência dos animais nas fazendas, além de colocar um produto final de melhor qualidade para o mercado consumidor.

A redução na idade da pesagem dos animais participantes de programas de melhoramento genético e, conseqüentemente, no número de pesagens é uma proposta de grande importância, pois permite aos criadores comercialização antecipada dos machos e seleção das novilhas para entrarem em reprodução a menor idade.

A reprodução em gado de corte é de fundamental importância econômica, portanto, a precocidade sexual irá interferir no desempenho lucrativo do sistema de produção.

As estações de acasalamento estabelecidas pela maioria dos criadores do PMGRN vão de outubro a março, no entanto, estão sendo atrasadas para o início de dezembro, com o objetivo de aproveitar a época de melhor qualidade das pastagens. Dessa forma, os animais atingem os 15 meses de idade no meio do período das águas quando estão ganhando peso, podendo, então, serem comercializados antes da época seca, com bom estado corporal e sem custos de suplementação.

Em melhoramento animal, o conhecimento de estimativas de parâmetros genéticos de características de importância econômica é o primeiro pré - requisito para o estabelecimento de um programa de seleção. Na literatura, Eler et al. (1996), Lôbo et al. (1997), Reyes et al. (1997) e Garnero et al. (1998), em rebanhos da raça Nelore, obtiveram para o peso aos 120 dias de idade, estimativas do coeficiente de herdabilidade direta, variando de 0,20 a 0,27, e para a materna, variando de 0,16 a 0,31. Para o peso aos 550 dias de idade, as estimativas da herdabilidade direta variaram de 0,19 a 0,35. Com relação às correlações genéticas entre P120 e P550, as estimativas variaram de 0,68 a 0,78. Quanto ao peso aos 455 dias de idade, não foram encontradas referências na literatura sobre este parâmetro.

Assim, o objetivo do presente estudo foi estimar parâmetros genéticos para as características de crescimento, pesos padronizados aos 120, 455 e 550 dias de idade (P120, P455 e P550), incluídas na avaliação genética do Programa de Melhoramento Genético da Raça Nelore, visando uma proposta para redução de pesagens e na idade das mesmas.

Material e Métodos

Foram utilizados dados de 51 rebanhos Nelore oriundos de fazendas participantes do Programa de Melhoramento Genético da Raça Nelore, localizadas nos Estados de Goiás, Mato Grosso do Sul, Mato Grosso, Minas Gerais, São Paulo, Maranhão e Bahia

Todas as fazendas adotavam estação de acasalamento que variava entre 60 e 90 dias de duração, sendo que algumas possuíam duas estações conforme a região e o ano. Geralmente, utilizava-se a inseminação artificial como técnica para reprodução dos animais, podendo ter repasse das fêmeas com touros, ou adotar-se exclusivamente monta natural. A desmama dos animais ocorria entre 7 e 8 meses de idade; a vacinação e os demais cuidados sanitários eram feitos conforme normas estabelecidas pelos órgãos de controle sanitário de cada região. Os rebanhos passavam por controle periódico e criterioso das características produtivas e reprodutivas, pois realizavam pesagens dos machos e das fêmeas de três em três meses (janeiro-abril-julho-outubro), do nascimento até 21 meses de idade e na ocasião da desmama, vaca e bezerro. A cada pesagem, a situação dos animais (normal, doente..) era informada. No cadastro, além da identificação do animal e dos pais, datas de cobertura e ou inseminação artificial, data e ocorrência no parto (normal, ajuda, cesárea) e a situação ao nascer (normal, natimorto..) também eram informadas.

As características estudadas foram os pesos padronizados para 120 (P120), 455 (P455) e 550 (P550) dias de idade. O ajuste dos pesos para as idades-padrão foi feito por interpolação linear.

Os grupos de contemporâneos para as características estudadas foram constituídos por animais nascidos na mesma fazenda, do mesmo sexo, na mesma época e mesmo ano e no mesmo regime alimentar.

Para P120, foram consideradas quatro épocas na formação dos grupos de contemporâneos, ou seja, 1 = janeiro a março, 2 = abril a junho, 3 = julho a setembro e 4 = outubro a dezembro; e para P455 e P550, duas épocas: 1 = outubro a março e 2 = abril a setembro.

A razão pela qual se consideraram apenas duas épocas para P455 e P550 foi a diminuição no número de animais nos grupos de contemporâneos quando se consideraram quatro épocas, o que acarretaria diminuição no valor estimado da acurácia.

A condição do regime alimentar foi estabelecida de modo a combinar todas as formas distintas de alimentação recebidas pelo animal durante toda sua vida, em cada idade-padrão analisada. Assim, um animal que, por exemplo, permaneceu em pasto do nascimento até o desmame e depois passou para regime de confinamento e a característica de interesse é o P455, seu regime alimentar (ra) foi codificado como 1133 (1 = pasto e 3 = confinamento), sendo que para P120, ra=1; para P240 (desmame), ra = 11 e para P365, ra = 113. Ou seja, todos os manejos adotados foram considerados em cada idade padrão analisada.

A idade da vaca ao parto (CIVP) foi agrupada em seis classes: menor que 36, de 36 a 47, de 48 a 59, de 60 a 71, de 72 a 119, e maior que 120 meses de idade

Consta na Tabela 1 o número de observações utilizadas nas análises para cada característica, bem como o número de animais na matriz de parentesco, o número de grupos de contemporâneos e o número de classes de idade da vaca ao parto.

As análises para estimação dos parâmetros genéticos foram efetuadas pelo método de Máxima Verossimilhança Restrita (REML) utilizando-se o software MTDFREML (Multiple Trait Derivative Free Restricted Maximum Likelihood), apresentado por Boldman et al. (1993). Para tal, foram fornecidos valores iniciais, previamente calculados, dos componentes de variância e, a cada interação, novos valores foram calculados até que a função fosse maximizada, quando então a convergência era atingida (com valores da função -2 log L inferiores a 10^-9). A cada convergência, o programa foi reiniciado, usando como valores de (co)variâncias iniciais aqueles obtidos na análise anterior, checando, assim, possíveis máximos locais e assegurando que o valor alcançado fosse o máximo global. Este procedimento foi repetido entre 2 e 5 vezes, segundo as características estudadas, até que os componentes estimados não diferissem entre uma análise e outra.

A princípio, cada característica foi analisada individualmente. Para o peso aos 120 dias de idade, foi utilizado um modelo estatístico que incluiu como efeitos fixos o grupo contemporâneo e a classe de idade da vaca ao parto e, como aleatórios, os efeitos genéticos direto, materno e de ambiente permanente da vaca.

Para as características P455 e P550, foi utilizado um modelo que incluiu os mesmos efeitos do anteriormente explicitado, exceto os efeitos aleatórios genéticos materno e de ambiente permanente da vaca, devido ter sido verificado em análises prévias que a contribuição destes efeitos na variância fenotípica total foram inferiores a 5%.

Posteriormente, análises bicaráter foram executadas com o objetivo principal de estimar as covariâncias e correlações genéticas entre as características estudadas. Desse modo, foram combinadas duas a duas, assumindo as covariâncias entre os efeitos maternal e direto como zero.

Em termos matriciais, o modelo geral utilizado pode ser descrito como:

em que: y = vetor das observações da característica; X = matriz de incidência dos efeitos fixos; b = vetor dos efeitos fixos; Z₁ = matriz de incidência do efeito genético direto; a = vetor de efeitos genéticos diretos aleatórios; Z₂ = matriz de incidência do efeito genético materno; m = vetor de efeitos genéticos maternos aleatórios; Z₃ = matriz de incidência do efeito de ambiente permanente; p = vetor de efeitos aleatórios de ambiente permanente; e = vetor de efeitos residuais aleatórios.

As pressuposições acerca das distribuições de y, a, m, p e e podem ser descritas como:

em que

V = Z₁G₁'_`+ Z₂MZ₂'+ 2Z₁CZ₂'+ Z₃PZ₃'+R

G = A s_a² em análise unicaráter e A Å G₀ em análise bicaráter; C =A s_am em análise unicaráter e A Å C₀ em análise bicaráter; M = A s_m² em análise unicaráter e A Å M₀ em análise bicaráter; P = Im s_c² em análise unicaráter e Im Å P₀ em análise bicaráter; R = In s_e² em análise unicaráter e In Å R₀ em análise bicaráter; m = número total de mães; n = número total de observações; NMV = multivariada normal.

Nestas designações, s_a² é a variância genética aditiva direta; s_am, a covariância genética entre os efeitos direto e materno; s_m², a variância genética materna; s_c², a variância de ambiente materno permanente; s_e², a variância do erro; G₀, M₀, P₀ e R₀, as matrizes de (co)variâncias genética direta, genética materna, de ambiente permanente e residual, respectivamente; C₀, uma matriz de covariâncias dos efeitos genéticos direto e materno. Assume-se que P ou P₀ e R ou R₀ não são correlacionados com os demais efeitos; A é a matriz de relações aditivas entre os animais; I, uma matriz identidade e Å, o produto direto entre as matrizes.

Resultados e Discussão

O produto de sistemas de bovinos de corte é o crescimento, sendo o objetivo da pecuária de corte a otimização do crescimento pré e pós-natal (Lanna, 1997). Dessa forma, o objetivo da pesquisa nesta área é estudar o crescimento dos animais através dos pesos obtidos a diferentes idades.

Na Tabela 2 são apresentadas as médias, os desvios-padrão, os valores máximos e mínimos e o coeficiente de variação de cada característica estudada.

A característica P120 está presente neste estudo por ser estabelecer ligação com as outras, já que neste período o bezerro sofre pouca influência de manejo e grande influência materna.

As médias de P120 e P550 estão de acordo com as encontradas na literatura consultada para a raça Nelore (Salles, 1995; Lôbo et al., 1997). No entanto, o valor médio encontrado para a característica P550 é superior ao valor citado, também na raça Nelore no Estado do Ceará, por Martins Filho et al. (1997) que provavelmente foi influenciado pelas condições ambientais.

O valor médio encontrado na literatura (Sumário de Touros Nelore 2001) é de 280,8 kg para peso aos 450 dias, de 59.353 observações das fazendas da Agro-Pecuária CFM, sendo superior ao valor encontrado no presente estudo P455 (250kg). Estes dados não foram previamente ajustados para idade, pois estas são incluídas no modelo matemático que se encarrega de efetuar o ajuste.

Em análises mais recentes realizadas por Lôbo et al. (2001), incluindo 147 rebanhos participantes do PMGRN, com 74.215 registros para peso aos 450 dias, o valor médio encontrado foi de 290 kg.

Há tendência desses valores fenótipicos aumentarem anualmente, já que os criadores que participam de um programa de melhoramento são orientados a utilizarem nos acasalamentos das matrizes, touros que possuem valores genéticos superiores a média do seu rebanho.

Os efeitos ambientais, como ano e época de nascimento, fazenda, onde estão incluídos o clima, solo, índice pluviométrico, manejo etc, são importantes na determinação de uma característica, e são denominados de efeitos fixos em uma avaliação genética.

A correta definição dos grupos contemporâneos é um dos aspectos mais importantes de qualquer programa de melhoramento genético e avaliação genética, pois parte dos efeitos fixos é considerada nos grupos.

As estimativas dos componentes de variância e herdabilidades obtidas nas análises unicarater são apresentadas na Tabela 3.

Para a característica P120, observou-se estimativa do coeficiente de herdabilidade direta de 0,29, que é superior as estimativas encontradas na literatura consultada para a raça Nelore (Lôbo et al., 1997; Reyes et al., 1997; Garnero et al., 2001). Os últimos autores utilizaram o mesmo modelo usado neste estudo em análise unicarater, empregando-se também o MTDFREML. No entanto, o valor do coeficiente de herdabilidade materna (0,08) é inferior aos valores relatados pelos referidos autores.

O coeficiente de herdabilidade direta estimado para a característica P550, ou seja, 0,47, é um valor alto, sendo superior às estimativas citadas na literatura consultada (Eler et al., 1995; Reyes et al., 1997; Garnero et al., 2001), que variaram de 0,16 a 0,31.

A estimativa de 0,51, obtida para o coeficiente de herdabilidade direta de P455, foi maior que as estimadas para os outros pesos. Todavia, verificaram-se na literatura (Mercadante,1995) valores de estimativas de herdabilidade que variam de 0,03 a 0,80 para P450, um artigo de revisão sobre parâmetros genéticos para características de crescimento em zebuínos de carne.

Nas Tabelas 4 e 5 são apresentadas as estimativas dos parâmetros genéticos obtidos para as características estudadas, mediante análises bicaráter.

Os coeficientes de herdabilidade direta obtidos para P120 (0,58 e 0,53), são superiores aos estimados por Lôbo et al. (1997) e Reyes et al. (1997). No entanto, as estimativas da herdabilidade materna (0,02 e 0,03) foram inferiores. Para a característica P550, as estimativas de herdabilidade direta (0,44 e 0,49) foram de magnitude maior que as estimadas por Eler et al. (1995) e Reyes et al. (1997).

Verifica-se também, na Tabela 5, que as estimativas dos coeficientes de herdabilidade direta para P455 (0,53 e 0,50), tal como as estimadas para P120 e P550, foram altas, indicando que grande parte da variância fenotípica das características estudadas pode ser atribuída aos efeitos aditivos dos genes. Por outro lado, as estimativas das correlações genéticas de P120 x P455, P120 x P550 e P455 x P550 (0,92, 0,93 e 0,96, respectivamente), sugerem que existe uma associação linear positiva entre as características, sendo que grande parte dos genes que determinam uma característica, determina, também, as demais.

Os valores encontrados no presente estudo para as variâncias genética direta e herdabilidade direta para P455 e P550, tanto em análise unicarater como na bicarater, são substancialmente altos, isto provavelmente se deva ao fato do modelo matemático adotado para análise destas duas características não incluirem o efeito de ambiente permanente da vaca.

Mercadante (1997), fez um estudo com o objetivo de apresentar estimativas dos componentes de (co)variância e dos parâmetros genéticos para algumas características de crescimento de fêmeas de um rebanho da raça Nelore, no intuito de quantificar a importância dos efeitos genéticos direto e materno e efeito de ambiente permanente. Nota-se, no estudo em que se avaliaram quatro modelos diferentes, que a não inclusão dos efeitos genético materno e de ambiente permanente, pode ocasionar mudanças significativas nas estimativas dos componentes e das herdabilidades, podendo até levar à superestimação destes valores e, conseqüentemente, o valor genético dos animais para estas características. A autora trabalhou com apenas duas características, peso aos 240 e aos 365 dias.

Portanto, o objetivo do nosso estudo foi estimar componentes e parâmetros genéticos não com o intuito de se estimar os valores genéticos dos animais, e sim analisar a viabilidade de utilizar o peso padronizado aos 455 dias como critério de seleção em substituição ao peso padronizado aos 550 dias no PMGRN.

Sabemos que o aumento da eficiência do setor cria, que é responsável pela maior parte do sistema de bovinos de corte, pode ser atingido pela redução na idade ao primeiro parto, sendo alcançada quando elevada proporção de novilhas apresenta cio ao início da estação de monta, aos 15 meses de idade.

A aceleração da maturidade sexual pode ser obtida através da utilização de uma série de técnicas de manejo, mas é ultimamente limitada pelo potencial genético (Lanna, 1997). Assim, estimular e orientar os criadores a selecionar os animais à idades mais jovens, e desperta-los para outras características de interesse econômico, além dos pesos, é de grande responsabilidade dos Programas de Melhoramento Genético, já que os reprodutores produzidos nestes rebanhos servirão a base da pirâmide, que são os rebanhos comerciais.

Conclusões

A característica P455 possui variabilidade genética aditiva de magnitude tal, que permite ser modificada por seleção. Isto, aliado ao fato de serem altas as estimativas das correlações genéticas observadas entre P455 e as demais características, sugere que o peso aos 455 dias de idade pode ser adotado como critério de seleção. No entanto, novos estudos são necessários para investigar possíveis associações do P455 com outras características de importância econômica.

A característica P120, por ser pouca influenciada pelo manejo, é a utilizada para se avaliar a habilidade maternal da mãe e de crescimento pré-desmama do bezerro, sendo de fundamental importância para o PMGRN, pois é usada como âncora nas análises bicaracteres para minimizar os efeitos de descartes posteriores à desmama.

Para a estimativa do P550, os animais devem ser pesados pelo menos mais uma vez, fazendo com que permaneçam no rebanho por pelo menos três meses a mais. Dessa forma, a alta correlação genética com o P455 permite a redução do número de pesagens e a utilização do P455 como critério de seleção.

Agradecimento

Ao Programa de Melhoramento Genético da Raça Nelore (PMGRN), dos criadores participantes e do GEMAC, pelo apoio.

Literatura Citada

Recebido em: 19/06/01

Aceito em: 22/08/02

• BOLDMAN, K.G.; KRIESE, L.A.; Van VLECK, L.D. et al. A manual for use for MTDFREML. A set of programs to obtain of variance and covariances [DRAF]. Lincoln: Agricultural Research Service, 1995. 120p.
• CFM. Sumário de Touros Nelore 2001. Publicação da Agro-Pecuária CFM. São José do Rio Preto: Agro-Pecuária CFM, 2001, 60p.
• ELER, J.P.; Van VLECK, L.D.; FERRAZ, J.B.S. et al. Estimation of variances due to direct and maternal effects for growth traits of Nelore cattle. Journal of Animal Science, v.77, p.3253-3258, 1995.
• GARNERO, A. del V.; LÔBO, R.B.; BEZERRA, L.A.F. et al. Comparação entre alguns critérios de seleção para Crescimento na Raça Nelore. Revista Brasileira de Zootecnia, v.30, n.3, p.714-718, 2001.
• LANNA, D.P. Fatores condicionantes e predisponentes da puberdade e da idade de abate. In: SIMPÓSIO SOBRE PECUÁRIA DE CORTE, 4., 1996, Piracicaba. Anais...Piracicaba, 1997. p.41-78.
• LÔBO, R.B.; REYES, A. de los; BEZERRA, L.A.F. et al. Avaliação genética de animais jovens, touros e matrizes. Ribeirão Preto: GEMAC Departamento de Genética/Universidade de São Paulo, 1997. 65p.
• LÔBO, R.B.; BEZERRA, L.A.F.; OLIVEIRA, H.N. et al. Avaliação genética de animais jovens, touros e matrizes. Ribeirão Preto: GEMAC Departamento de Genética/Universidade de São Paulo, 2001. 60p.
• MARTINS FILHO, R.; LOBO, R.N.B.; LIMA, F.A.M. et al. Parâmetros genéticos e fenotípicos de pesos e ganhos em pesos de bovinos Zebus no estado do Ceará. In: REUNIÃO ANUAL DA SOCIEDADE BRASILEIRA DE ZOOTECNIA, 34, 1997, Juiz de Fora. Anais... Juiz de Fora,1997. 600p. p.248-250.
• MERCADANTE, M.E.Z., LÔBO, R.B Estimativas de (co)variâncias e parâmetros genéticos dos efeitos direto e materno de características de crescimento de fêmeas de um rebanho Nelore. Revista Brasileira de Zootecnia, v.26, n.6, p.1124-1133, 1997.
• REYES, A. de los; LÔBO, R.B.; BEZERRA, L.A.F. et al. Variabilidade genética de características de crescimento alternativas para a seleção em gado de corte. In: REUNIÃO ANUAL DA SOCIEDADE BRASILEIRA DE ZOOTECNIA, 34, 1997, Juiz de Fora. Anais... Juiz de Fora, 1997. 600p. p.245-247.
• SALLES, P.A. Critérios de seleção para características de crescimento em machos da raça Nelore. Ribeirão Preto: Universidade de São Paulo, 1995. 65p. Dissertação (Mestrado em Genética) - Universidade de São Paulo, 1995.
• SAS Institute. SAS user's guide for windows enviroment. 6.08. Cary: 1995.

1

Parte da dissertação de mestrado da 1

^a autora apresentada à FCAV-UNESP, Jaboticabal

Datas de Publicação

• Publicação nesta coleção
  03 Jul 2003
• Data do Fascículo
  Fev 2003

Histórico

• Aceito
  22 Ago 2002
• Recebido
  19 Jun 2001