Interação genótipo x ambiente e estimativas de parâmetros genéticos dos pesos aos 205 e 365 dias de idade de bovinos Nelore

Genotype by environment interaction and genetic parameter estimates for 205 and 365 body weight of Nelore cattle

ZOOTECNIA E TECNOLOGIA E INSPEÇÃO DE PRODUTOS DE ORIGEM ANIMAL

Interação genótipo x ambiente e estimativas de parâmetros genéticos dos pesos aos 205 e 365 dias de idade de bovinos Nelore

Genotype by environment interaction and genetic parameter estimates for 205 and 365 body weight of Nelore cattle

A.B. Fridrich^I; M.A. SilvaI, ^*; B.D. Valente^I; J.E.R. Sousa^I; G.S.S. Corrêa^I; I.C. Ferreira^I; R.V. Ventura^I; L.O.C. Silva^II

^IEscola de Veterinária - UFMG - Caixa Postal 567- 30123-970 - Belo Horizonte, MG

^IIEmbrapa Gado de Corte - Campo Grande, MS

RESUMO

Dados de pesos aos 205 (P205) e 365 (P365) dias de idade, de 46.408 animais Nelore, nascidos entre 1976 e 2000, provenientes de 530 rebanhos dos diversos estados brasileiros, foram utilizados para avaliar as interações genótipo x ambiente (IGA) e estimar herdabilidades direta e materna dos pesos pelo método da máxima verossimilhança restrita. O modelo estatístico utilizado incluiu efeitos fixos de grupo contemporâneo, idade da vaca ao parto (covariável) e efeitos aleatórios genéticos aditivo direto e materno. As correlações genéticas, estimadas para as características P205 e P365, cada uma considerada como característica distinta em cada uma das regiões Sul (S), Sudeste (SE), Centro-Oeste (CO), Norte (N) e Nordeste (NE), foram, respectivamente, 0,86 e 0,84, 0,64 e 0,35, 0,75 e 0,42, 0,79 e 0,86, 0,92 e 0,81, 0,95 e 0,83, 0,83 e 1,00, 0,88 e 0,81, 0,33 e 0,85, 0,63 e 0,99 para S e SE, S e CO, S e N, S e NE, SE e CO, SE e N, SE e NE, CO e N, CO e NE e N e NE. As baixas correlações genéticas indicaram que há efeito da interação genótipo x ambiente nas combinações que envolvem as regiões S/CO, S/N, S/NE, CO/NE e N/NE para P205 dias de idade e S/CO e S/N para P365 dias de idade, havendo, portanto, necessidade de avaliação genética regional quando se consideram regiões bastante distintas.

Palavras-chave: bovino, Nelore, herdabilidade, interação genótipo x ambiente

ABSTRACT

Body weight records at 205 (205BW) and 365 (365BW) days of age of 46,408 Nelore animals, born from 1976 to 2000 period in 530 Nelore herds of several states of Brazil, were used to evaluate genotype by environment interactions and to estimate genetic and maternal heritability by restricted maximum likelihood methodology. The statistical model included the fixed effects of contemporary group and age of cow (covariate), and the random direct and maternal additive genetic effects. The genetic correlation for 205BW and 365BW, each one considered as different traits in each of the South (S), Southeast (SE), Central west (CO), North (N) and Northeast (NE) regions were 0.86 and 0.84, 0.64 and 0.35, 0.75 and 0.42, 0.79 and 0.86, 0.92 and 0.81, 0.95 and 0.83, 0.83 and 1.00, 0.88 and 0.81, 0.33 and 0.85, 0.63 and 0.99 for S/SE, S/CO, S/N, S/NE, SE/CO, SE/N, SE/NE, CO/N, CO/NE and N/NE, respectively. There is a significant genotype by environment interaction for 205BW in the combination involving the S/CO, S/N, S/NE, CO/NE, and N/NE regions. There are significant genotype by environment interaction effects in the combination involving the S/CO and S/N regions for 365BW. Based on genetic x environment interaction results, regional genetic evaluation is recommend for the very distinct regions.

Keywords: bovine, Nelore, heritability, genotype by environment interaction

INTRODUÇÃO

A investigação da existência de interação genótipo x ambiente (IGA) é essencial para os programas de melhoramento genético dos animais. Os estudos deste tipo de interação começaram a partir do século XX e estão basicamente relacionados ao ambiente no qual a seleção deve ser praticada, se no ambiente em que os animais serão criados ou no ambiente melhorado que permitiria que os animais expressassem todo o seu potencial genético.

Hammond (1947) propunha que os animais deveriam ser selecionados em ambiente que permitisse sua máxima expressão, mas Falconer (1952) sugeria que o conjunto de genes responsáveis pela expressão de determinada característica poderia variar em diferentes ambientes de criação, caso esta característica apresentasse IGA.

O Brasil, entre os países que apresentam grande extensão territorial, o clima, o manejo e a alimentação dos animais podem variar muito entre as regiões. O valor genético de diferentes animais pode apresentar diferentes sensibilidades a esta variação ambiental, o que caracteriza a existência de IGA. Dessa maneira, diferentes regiões podem apresentar diferentes ordens de classificação dos valores genéticos preditos dos animais. Isso faz com que avaliações genéticas com base em modelos que consideram características mensuradas em diferentes ambientes como uma só característica, como os modelos utilizados por Marques et al. (2000) e Pereira et al. (2001), possam apresentar menor eficiência na predição de valores genéticos para determinados ambientes.

Os resultados obtidos por Alencar et al. (2005), por meio de componentes principais, indicam interação touro x época de nascimento significativa para pesos e ganhos de pesos de animais Canchim até os 12 meses de idade, e sugerem que essa interação deva ser considerada nos estudos de avaliação genética.

Fridrich et al. (2005), ao estudarem pesos aos 205 e 365 dias de idade de animais Tabapuã nas diversas regiões brasileiras e considerarem a expressão da mesma característica em ambientes diferentes como sendo expressão de diferentes características, obtiveram resultados que indicam existência de IGA, principalmente para as combinações de regiões que envolvem a região nordeste.

Ribeiro (2006) avaliou o peso aos 205 e 550 dias de idade, além do ganho de peso entre 205 e 550 dias de idade de animais Nelore de diferentes rebanhos do estado de São Paulo e um rebanho do estado do Mato Grosso do Sul. Ao considerar, informações provenientes de um universo de ambientes de criação mais estreito, o autor estimou correlações genéticas entre observações da mesma característica coletadas em ambientes diferentes que variaram entre 0,82 a 1,00. Esses resultados sugerem pequeno ou nenhum efeito da IGA nos rebanhos estudados. Ferreira et al. (2001) estimaram altas correlações genéticas para peso aos 120 dias e para peso aos 205 dias de idade em animais Nelore criados em duas fazendas localizadas em regiões bem distintas do estado de Minas Gerais. As correlações obtidas foram, respectivamente, 0,96 e 0,95, estimativas que sugerem ausência de efeito importante de IGA. Por outro lado, para o peso aos 550 dias, a correlação estimada foi 0,53, indicando efeito de IGA.

Portanto, o assunto é ainda controverso e merece estudos mais detalhados que considerem maior amplitude de ambientes e envolvam maior número de rebanhos distribuídos nos vários estados brasileiros, para melhor avaliação genética dos animais, o que constitui a base de qualquer programa de melhoramento.

Neste trabalho, investigou-se a existência de interação genótipo x ambiente nos pesos aos 205 e 365 dias de idade de bovinos Nelore, considerando peso (205 ou 365 dias de idade) como característica distinta em cada região, bem como estimaram-se os parâmetros genéticos desses pesos.

MATERIAL E MÉTODOS

Os dados utilizados neste estudo são de 46.408 animais Nelore, nascidos entre 1976 e 2000, provenientes de 530 rebanhos dos diversos estados brasileiros incluídos no Controle de Desenvolvimento Ponderal da Associação Brasileira de Criadores de Zebu (ABCZ) e cedidos pela Embrapa - Gado de Corte, localizada em Campo Grande, MS.

As características estudadas foram os pesos aos 205 (P205) e aos 365 dias de idade (P365), sendo os pesos padronizados por interpolação. Foram formados três arquivos para as análises. O primeiro arquivo continha a identificação do animal, do pai e da mãe; o segundo e o terceiro arquivo, continham a identificação de animal e mãe, grupo contemporâneo, idade da mãe ao parto (covariável) e os pesos aos 205 e 365 dias de idade, respectivamente, considerados como características distintas em cada região.

Para avaliar o efeito da IGA sobre as características analisadas nas diversas regiões brasileiras, os arquivos foram constituídos apenas com informações de animais criados exclusivamente em pasto. Os estados foram agrupados em cinco regiões distintas: 1-Sul (S), 2-Sudeste (SE), 3-Centro-Oeste (CO), 4-Norte (N) e 5-Nordeste (NE). Foram consideradas duas estações de nascimento.

Foram também definidos os "touros de conexão", isto é, touros com número mínimo de três progênies em cada uma de duas ou mais regiões. Os grupos contemporâneos (GC) foram formados por ano de nascimento, propriedade, estação de nascimento e sexo, totalizando 4.490 grupos contemporâneos.

As análises estatísticas dos pesos foram realizadas com o auxílio do procedimento general linear models (GLM), contido no SAS (Statistical..., 1997). Na Tab. 1, encontram-se as estatísticas descritivas das características da amostra estudada nas regiões.

Para estimação dos componentes de (co)variância e predição dos valores genéticos, utilizou-se o método de máxima verossimilhança restrita aplicado a um modelo animal, usando o programa Multiple Trait Derivative-Free Restricted Maximum Likelihood (MTDFREML), desenvolvido por Boldman et al. (1995). Para estimação dos parâmetros genéticos, a saber, componentes de variância genéticos aditivos, maternos, covariância e correlação entre efeitos diretos e maternos, foram feitas análises unicaracterísticas com o emprego do modelo que incluiu o efeito fixo de GC, a idade da vaca ao parto como covariável (efeito linear e quadrático) e os efeitos genéticos aditivos direto e materno, bem como a covariância entre efeitos genético aditivo e o materno .

Os efeitos da IGA foram obtidos pelas correlações genéticas entre cada peso, considerando a mesma característica (P205 ou P365) como característica distinta nas várias regiões. Os dados foram analisados considerando-se duas características (análise bicaracterística) em cada análise, segundo o modelo:

y = Xb + Zu + Wm + ε;

em que:

y = vetor de pesos padronizados para 205 ou 365 dias de idade dos animais.

X, Z e W= matrizes de incidência dos efeitos fixos e aleatórios;

b = vetor de efeito fixo de grupo contemporâneo (GC) e da idade da vaca ao parto como covariável;

u = vetor dos efeitos aleatórios de valores genéticos aditivos diretos do animal;

m = vetor dos efeitos aleatórios de valores genéticos aditivos maternos e;

e = vetor dos erros aleatórios associados às observações.

Outra estatística utilizada para o estudo da interação genótipo x ambiente foi a correlação de posto (Spearman) dos valores genéticos preditos para cada par de combinações das regiões estudadas, para P205 e P365.

RESULTADOS E DISCUSSÃO

As variâncias, covariâncias, herdabilidades (h²) e correlação genética entre os efeitos genéticos aditivos direto e materno (r_am) dos pesos aos 205 (P205) e aos 365 (P365) dias de idade, estimadas para cada região, são apresentadas, respectivamente, nas Tab. 2 e 3.

As estimativas de parâmetros genéticos variaram muito em todas as regiões para ambas as características. Esta grande variação pode ter ocorrido em razão do número diferente de informações em cada região. Na região CO, a segunda com maior número de animais, foram observados valores mais elevados da σ²_a e σ²_m para ambas as características P205 (125,92 e 116,69) e P365 (189,76 e 111,75). Grande heterogeneidade de variância genética aditiva direta e materna foi observada para P205 e P365 entre as regiões SE e CO, as regiões com maior número de animais. Esta heterogeneidade não foi observada por Valente et al. (2008), que estimaram estruturas de covariância de peso em função da idade de bovinos Nelore das regiões SE e CO. Na região S foi observada menor σ²_a, tanto para P205 (2,69) quanto para P365 (26,94). As análises na região S forneceram valores atípicos em relação às outras regiões para os parâmetros genéticos dos animais, mas essa foi a região que apresentou o menor número de dados, automaticamente, reflexo dos menores números de rebanhos, touros e de matrizes.

Valores positivos de correlação entre efeitos genético aditivo direto e materno foram estimados somente para P365 nas regiões SE e N, o que indica que ao se utilizar, nessas regiões, o valor genético aditivo direto no estabelecimento do critério de seleção não haverá impacto negativo no mérito genético aditivo materno dos animais. Os valores estimados da correlação na região S, referente ao P205 e P365, devem ser encarados com reserva, em razão do baixo número de dados existentes para essa região. As r_am foram negativas para P205 em todas as regiões e para P365 nas regiões CO e NE, o que evidencia um antagonismo entre os efeitos genéticos direto e materno, ou seja, ao se utilizar, como critério de seleção, o valor genético aditivo direto, resultará em impacto negativo no mérito genético aditivo materno dos animais

Correlações entre efeito genético aditivo direto e materno negativas para P205 (-0,70) e P365 (-0,74) foram encontradas por Ribeiro et al. (2001), ao trabalharem com bovinos Nelore no estado da Paraíba. Na região S, a r_am para P365 foi igual a zero, o que indica dissociação entre os efeitos genéticos aditivos direto e materno. Fridrich et al. (2005) encontraram r_am iguais a zero nas regiões CO e NE para a característica P205, e na região S para P365, ao trabalharem com bovinos Tabapuã.

De modo geral, as estimativas de herdabilidade direta (h²_a) foram menores do que as herdabilidades maternas (h²_m) para P205, e o inverso ocorreu para P365, o que indica grande influência do efeito materno sobre P205, uma vez que existe maior dependência do bezerro em relação às suas mães até a desmama.

As estimativas de h²_a para P205 foram semelhantes às encontradas por Souza et al. (2003) que estudaram apenas o peso à desmama de bovinos Nelore nas diversas regiões brasileiras e por Fridrich et al. (2005) para P205 nas regiões S e SE em bovinos Tabapuã. O valor da h²_a foi menor nas regiões S para P205 e S e NE para P365, em razão do aumento da σ²_e, que foi proporcionalmente maior do que a σ²_a para essas regiões. Para P205, as regiões SE e CO apresentaram maiores estimativas de σ²_a e h²_a, sugerindo que a seleção baseada no peso à desmama seria mais efetiva nessas regiões.

Apesar da variação observada nas estimativas de parâmetros genéticos para P205 e P365, os valores encontram-se na faixa de variação retratada na literatura para bovinos Nelore.

As estimativas das (co)variâncias genéticas aditivas direta e materna, fenotípicas e residuais de P205 e P365, consideradas características distintas nas diversas regiões, são apresentadas nas Tab. 4, 5, 6 e 7, respectivamente.

Em geral, ocorreu ampla variação nas estimativas de variâncias e covariâncias genéticas direta, materna e fenotípica, tanto para P205 quanto para P365 em análises conjuntas, considerando-se cada uma das características como característica distinta nas combinações de pares de regiões estudadas, o que pode ser atribuído à expressiva variação no número de informações disponíveis em cada região, principalmente quando são cotejados os números de dados disponíveis nas regiões S, N e NE com os disponíveis nas regiões SE e CO que representam 87,17% de todos os dados.

As estimativas de herdabilidade e correlações genéticas dos pesos aos 205 dias de idade (P205) e aos trezentos e sessenta e cinco dias de idade (P365), cada um deles considerado característica distinta nas diversas regiões, obtidas em análises bicaracterísticas, são apresentadas nas Tab. 8 e 9, respectivamente.

As correlações genéticas encontradas para P205 e P365, ambos considerados características distintas nas diversas regiões, indicam que, na desmama (P205), o efeito da IGA não foi observado somente nas combinações que envolveram a região SE com as demais regiões e em CO/N. Para as demais regiões o efeito da IGA foi importante. Considerando a característica P365, o efeito da interação foi evidenciado somente nas combinações S/CO e S/N. Nas demais regiões não houve efeito da IGA. As estimativas de correlações genéticas para a maioria das combinações de regiões foram altas, sugerindo que grande parte dos genes de ação aditiva que influencia uma característica também influencia a outra. De acordo com Robertson (1959), quando o valor da correlação genética for menor que 0,8, a interação genótipo x ambiente passa a assumir papel importante na expressão da característica.

Fridrich et al. (2005), ao avaliarem pesos de bovinos Tabapuã nas diversas regiões brasileiras, também encontraram interação genótipo ambiente para P205 nas combinações que envolviam as regiões S e NE e SE e NE e, para P365, em todas as combinações que envolviam a região NE

As altas correlações observadas nas Tab. 8 (0,95) para SE/N e Tab. 9 para SE/NE (1,00) e N/NE (0,99), próximas ou iguais a um, mesmo considerando número diferente de informações nas regiões, indicam que os mesmos grupos de genes são responsáveis pela expressão das duas características estudadas.

De maneira geral, as estimativas das herdabilidades (moderadas) dos pesos P205 e P365 foram semelhantes e indicam que a raça Nelore apresenta variabilidade genética aditiva suficiente para que ocorra resposta à seleção baseada nestes pesos. Ferraz Filho et al. (2002), entretanto, observaram diminuição da magnitude das herdabilidades com o avanço da idade dos animais, em decorrência do menor número de observações com o avanço da idade, o que redunda na obtenção de estimativas de herdabilidades menos precisas e acuradas.

Os altos valores das herdabilidades maternas para P205 sugerem que o efeito materno de e ser incluído nos modelos de avaliação genética em razão da sua expressiva participação nesta fase.

As correlações de posto entre os valores genéticos dos 62 melhores reprodutores com mais de 150 filhos em pelo menos três regiões, obtidas em análise bicaracterística, do P205 e P365 dias de idade, são apresentadas na Tab.10.

Os resultados das correlações de posto, ao se considerar o P205 ou o P365, cada uma como característica distinta nas diversas regiões, indicam que os reprodutores, em análise envolvendo pares de região, foram ordenados (geneticamente) de forma semelhante quando se consideraram as combinações das regiões S/SE, S/NE, SE/N, SE/NE e CO/N. Porém, quando se consideraram as combinações S/CO e S/N, o mesmo não foi verificado. Os resultados indicaram baixas correlações de posto 0,71 (P205) e 0,73 (P365) para S/CO e 0,78 (P205) e 0,78 (P365) para S/N. Quando se considerou P205 para as combinações S/CO, S/N, SE/CO, SE/NE, CO/NE e N/NE, os valores das correlações variaram de 0,71 a 0,79, o que mostra que a classificação genética de um reprodutor em uma região seria diferente da realizada em outra região. Quando se considerou P365, somente as combinações que envolveram as regiões S/CO, S/N e SE/CO apresentaram correlações de posto com valores de 0,73, 0,78 e 0,76 que são menores que 0,80. Estes resultados reforçam a idéia de que os reprodutores deveriam ser testados por região, principalmente quando se consideram regiões muito distintas. As grandes diferenças climáticas e nutricionais em que os animais são submetidos em cada região podem ser responsáveis pelos diferentes desempenhos.

Fridrich et al. (2005) encontraram correlação de posto para P205 em bovinos Tabapuã iguais a -0,81 quando consideraram a combinação das regiões S e NE, para P365 dias de idade as correlações de posto foram iguais a -0,51 para as regiões S e NE, 0,76 para SE e NE e 0,53 para as regiões CO e NE. Esses resultados indicam que os reprodutores com melhores resultados nas regiões S, SE e CO não apresentam o mesmo desempenho da região NE.

Buchanan e Nielsen (1979) e Bertrand et al. (1985) observaram efeitos das interações touro x região, touro x rebanho/região e touro x GC/região sobre a ordem de classificação dos animais entre as regiões, rebanhos e GC dentro de região, quando foram analisados os pesos ao nascimento e ao desmame dos animais, respectivamente. Matos et al. (1996) observaram, em bovinos Hereford, que a interação reprodutor x região dentro de país (Uruguai) foi grande o suficiente para causar mudanças na ordem de classificação dos reprodutores entre regiões.

A existência da IGA encontrada nas combinações que envolveram as regiões S, SE, CO, N e NE indicam que a utilização de modelos de avaliação que consideram a IGA seria de grande valia para maior qualidade das predições de valores genéticos em diferentes regiões.

CONCLUSÕES

Não se observa interação genótipo x ambiente para pesos aos 205 e 365 dias de idade quando se consideram somente as regiões Sudeste e Centro-Oeste, que representam 87,17% dos dados. Os resultados obtidos sinalizam a necessidade de avaliação genética regional quando se consideram observações de regiões S, N ou NE.

Recebido em 5 de outubro de 2007

Aceito em 3 de junho de 2008

• ALENCAR, M.M.; MASCIOLI, A.S.; FREITAS, A.R. Evidências de interação genótipo x ambiente sobre características de crescimento de bovinos de corte. Rev. Bras. Zootec, v.34, p.489-495, 2005.
• BERTRAND, J.K.; BERGER, P.J.; WILHAM, R.L. Sire x environment interactions in beef cattle weaning weight field data. J. Anim. Sci, v.60, p.1397-1402, 1985.
• BOLDMAN, K.G.; KRIESE, L.A.; VAN VLECK, L.D. et al. A manual for use of MTDFREML : a set of programs to obtain estimates of variance and covariance. Lincoln : Agricultural Research Service, 1995. 120p.
• BUCHANAN, D.S.; NIELSEN, M.K. Sire by environment interactions in beef cattle field data. J. Anim. Sci, v.48, p.307-312, 1979.
• FALCONER, D.S. The problem of environment and selection. Am. Nat, v.86, p.293-298, 1952.
• FERRAZ FILHO, P.B.; RAMOS, A.A.; SILVA, L.O.C. et al. Tendência genética dos efeitos direto e materno sobre os pesos à desmama e pós-desmama de bovinos da raça Tabapuã. Rev. Bras. Zootec, v.31, p.635-640, 2002.
• FERREIRA, V.C.P.; PENNA, V.M.; BERGMANN, J.A.G. et al. Interação genótipo-ambiente em algumas características produtivas de gado de corte no brasil. Arq. Bras. Med. Vet. Zootec, v.53, p.385-392, 2001.
• FRIDRICH, A.B.; SILVA, M.A.; FRIDRICH, D. et al. Interação genótipo x ambiente e estimativas de parâmetros genéticos de características ponderais de bovinos Tabapuã. Arq. Bras. Med. Vet. Zootec, v.57, p.663-672, 2005.
• HAMMOND, J. Animal breeding in relation to nutrition and environment conditions. J. Anim. Sci, v.22, p.195-213, 1947.
• MARQUES, L.F.A.; PEREIRA, J.C.C.; OLIVEIRA, H.N. et al. Análise de características de crescimento da raça Simental. Arq. Bras. Med. Vet. Zootec, v.52, p.527-533, 2000.
• MATOS, S.D.; BERTRAND, J.K.; HERRING, W.O. et al. Sire and maternal grandsire by environment interaction for weaning weight in a Hereford beef cattle population in Uruguay. J. Anim. Sci, v.74, suppl.1, p.106, 1996. (Abstract).
• PEREIRA, E.; ELER, J.P.; FERRAZ, J.B.S. Análise genética de algumas características reprodutivas e suas relações com o desempenho ponderal na raça Nelore. Arq. Bras. Med. Vet. Zootec, v.53, p.720-727, 2001.
• RIBEIRO, M.N.; PIMENTA FILHO, E.C.; MARTINS, G.A. et al. Herdabilidade para efeitos diretos e maternos de características de crescimento de bovinos Nelore no estado da Paraíba. Rev. Bras. Zootec, v.30, p.1224- 1227, 2001.
• RIBEIRO, S. Estudo da interação genótipos x ambiente em algumas características produtivas na raça Nelore. 2006. 72f. Tese (Doutorado) - Faculdade de Zootecnia e Engenharia de Alimentos, Universidade de São Paulo, Pirassununga, SP.
• ROBERTSON, A. The sampling variance of the genetic correlation coefficient. Biometrics, v.15, p.469-485, 1959.
• SOUZA, J.C.; GADINI, C.H.; SILVA, L.O.C. et al. Estimates of genetic parameters and evaluation os genotype x environment interaction for weaning weight in Nellore cattle. Arch. Latinoam. Prod. Anim., v.11, p.94-100, 2003.
• STATISTICAL analysis system. Versão 6.12 para windows. Cary, NC: SAS Institute, 1997. 1290p.
• VALENTE, B.D.; SILVA, M.A.; SILVA, L.O.C. et al. Estruturas de Covariância de peso em função da idade de animais Nelore das regiões Sudeste e Centro-Oeste do Brasil. Arq. Bras. Med. Vet. Zootec, v.60, p.389-400, 2008.

Datas de Publicação

• Publicação nesta coleção
  15 Ago 2008
• Data do Fascículo
  Ago 2008

Histórico

• Recebido
  05 Out 2007
• Aceito
  03 Jun 2008