Melhoramento Genético 4d6r

Introdução

A pecuária leiteira no Paraná é reconhecida por seus rebanhos de alto potencial genético para produção e composição de leite. Para que possam ser aplicadas práticas efetivas de seleção para essas características, é necessário que parâmetros genéticos sejam estimados na própria população sob seleção. A correta estimação dos parâmetros genéticos é fundamental em programas de melhoramento genético, pois permite predizer o valor genético dos animais e identificar animais geneticamente superiores.

No Paraná, estimativas para as características de produção (Ribas et al., 1983; Ribas et al., 1993; Almeida et al., 1997; Falcão et al., 2006) e composição do leite (Almeida et al., 1997) têm sido obtidas, mas, em virtude das diferenças genéticas da população, de ambiente, do tipo de análise e do método de estimação de componentes de co-variância, essas estimativas têm variado consideravelmente.

Métodos Bayesianos têm sido utilizados para estimação dos componentes de co-variância e dos parâmetros genéticos para avaliação genética de animais (Gianola & Fernando, 1986; Rekaya et al., 2003; Falcão et al., 2006). Esses métodos permitem a análise de grandes conjuntos de dados, não requerem soluções para as equações de modelos mistos e propiciam estimativas diretas e acuradas dos componentes de co-variância e dos parâmetros genéticos das características de interesse.

Os objetivos neste trabalho foram estimar as herdabilidades e as correlações genéticas para as produções de leite (PL), gordura (PG) e proteína (PP) e para os percentuais de gordura (%G) e proteína (%P) do leite em bovinos da raça Holandesa no estado do Paraná utilizando-se inferência Bayesiana.

Material e Métodos

Foram utilizados 195.886 registros de lactações encerradas de vacas da raça Holandesa, cujo controle foi realizado pelo Programa de Análise de Rebanhos Leiteiros do Paraná (PARLPR), desenvolvido em convênio com a Associação Paranaense de Criadores de Bovinos da Raça Holandesa (APCBRH), localizada em Curitiba, Paraná. Os animais foram controlados oficialmente pelo Serviço de Controle Leiteiro Mensal da APCBRH entre janeiro de 1992 e dezembro de 2003.

Para garantir a consistência das informações, o conjunto de dados original foi editado no programa computacional SAS (Statistical Analysis System) (2000), impondo-se restrições que resultaram na eliminação de observações inconsistentes. Foram eliminadas informações de vacas cujos pais tinham menos de oito filhas e de rebanhos com menos de dez lactações. Definiram-se duas estações de parto, uma de abril a setembro e outra de outubro a março, e grupos de contemporâneos (GC) formados por lactações de mesmo rebanho, ano e estação de parição. Apenas informações de grupos de contemporâneos com dez observações ou mais foram mantidas.

Após a aplicação de todas as restrições, o conjunto de dados final permaneceu com 117.082 registros de primeira à décima lactação para as produções de leite (PL), gordura (PG) e proteína (PP), ajustados para os 305 dias de lactação de 49.676 vacas provenientes de 308 rebanhos, e com 116.747 registros de primeira à décima lactação para os percentuais de gordura (%G) e de proteína (%P) ajustados para os 305 dias de lactação, de 49.612 vacas provenientes de 308 rebanhos. Os rebanhos estavam distribuídos em sete bacias leiteiras do estado do Paraná: Castro, Carambeí, Witmarsum, Arapoti, Sul, Norte e Oeste.

Tanto para as análises unicaracteres como para as tricaracteres, foram geradas cadeias de 120.000 iterações, com descarte inicial de 20.000 iterações e intervalo de amostragem de 100 iterações, resultando em 1.000 amostras dos componentes de co-variância para cada análise realizada. A monitoração da convergência das cadeias geradas pelo amostrador de Gibbs foi realizada por meio de análise gráfica e dos testes de diagnóstico de Geweke e de Heidelberger & Welch, disponíveis no CODA (Convergence Diagnosis and Output Analysis), implementado no programa R (2004).

Foram construídos os intervalos de credibilidade para todos os componentes de co-variância e parâmetros genéticos estimados, a 90% de probabilidade.

Resultados e Discussão

As médias e os desvios-padrão para as produções de leite (PL), de gordura (PG) e de proteína (PP), em kg, ajustados para os 305 dias de lactação foram 8.181,23 ± 1.903,51; 270,88 ± 67,58; e 249,01 ± 55,65, respectivamente, e para os percentuais de gordura (%G) e de proteína (%P) foram 3,33 ± 0,45 e 3,06 ± 0,19, respectivamente.

Houve indicação de convergência para todas as cadeias por meio de análise gráfica e da utilização dos testes de diagnóstico nas análises uni e tricaracteres.

As médias posteriores das estimativas de herdabilidade para PL, PG, PP, %G e %P obtidas nas análises unicaracteres foram iguais a 0,26; 0,28; 0,25; 0,60 e 0,58, respectivamente, semelhantes às obtidas nas análises tricaracteres (Tabela 1).

As médias foram obtidas com grande precisão, pois apresentaram baixo erro-padrão e se mantiveram no intervalo de credibilidade calculado. Estes resultados indicam que as características estudadas possuem variação genética aditiva suficiente para boa resposta à seleção. Estimativas semelhantes foram obtidas por Chauhan & Hayes (1991) utilizando o método da máxima verossimilhança restrita. Outros autores estimaram os parâmetros genéticos para características produtivas de bovinos da raça Holandesa utilizando REML e também encontraram valores semelhantes aos encontrados neste trabalho. As estimativas para PL, PG e %G encontradas por Ribas et al. (1993) no estado do Paraná foram de 0,23 ± 0,05; 0,28 ± 0,05; e 0,47 ± 0,08, respectivamente. Chauhan & Hayes (1991) encontraram, para vacas em primeira lactação, valores de herdabilidade para PL, PG, PP, %G e %P de 0,29; 0,31; 0,25; 0,65 e 0,61, respectivamente. Os valores de herdabilidade estimados por Urioste et al. (2003) para PL foram de 0,23; 0,24 e 0,25 para a primeira, segunda e últimas lactações.

Os valores encontrados neste estudo para as características PL e PG foram superiores ao descritos por Ferreira & Fernandes (2000), que estimaram os parâmetros genéticos para PL e PG utilizando o método da máxima verossimilhança restrita sem derivadas (DFREML) e encontraram valores de herdabilidade iguais a 0,18, tanto para PL como para a PG.

As correlações genéticas foram altas e positivas para PL e PG, PL e PP, PG e PP e para %G e %P e altas e negativas para PL e %G e PL e %P (Tabela 2), indicando a existência de associação genética entre essas características. A seleção para aumentar a PL promoverá aumento na PG e PP e redução nos %G e %P.

Correlações genéticas altas e positivas entre PL e PG também foram observadas por Chauhan & Hayes (1991), 0,45, Ribas et al. (1993), 0,78, Almeida et al. (1997), 0,52, Ferreira & Fernandes (2000), 0,98, Freitas et al. (2001), 0,86 a 1,00, e Weber et al. (2005), 0,91 a 1,00. Para PL e %G, Chauhan & Hayes (1991) e Ribas et al. (1993) também relataram correlações genéticas negativas, de -0,49 e -0,13, respectivamente.

Conclusões

As produções de leite, gordura e proteína e os percentuais de gordura e proteína possuem variabilidade genética aditiva suficiente para justificar sua inclusão em programas de melhoramento de bovinos da raça Holandesa no estado do Paraná. A seleção para aumentar a produção de leite nesta população promoverá aumento nas produções de gordura e proteína e redução nos percentuais de gordura e proteína.