Efeitos genéticos e ambientais que afetam a produção de leite e duração da lactação de cabras mestiças no estado da Paraíba

Genetic and environmental effects that affect milk production and lactation length of crossbred goats in the state of Paraiba

MELHORAMENTO, GENÉTICA E REPRODUÇÃO

Efeitos genéticos e ambientais que afetam a produção de leite e duração da lactação de cabras mestiças no estado da Paraíba

Genetic and environmental effects that affect milk production and lactation length of crossbred goats in the state of Paraiba

Edgard Cavalcanti Pimenta Filho^I; José Lindenberg Rocha Sarmento^II; Maria Norma Ribeiro^{III, IV}

^IProfessor do Departamento de Zootecnia/CCA/UFPB, 58397-000, Areia - PB. E.mail: [e-mail]

^IIDoutorando em Genética e Melhoramento. Bolsista da CAPES/UFV, Viçosa-MG. E.mail: [e-mail]

^IIIProfessora do Departamento de Zootecnia/UFRPE, Recife - PE. E.mail: [e-mail]

^IVBolsista do CNPq

RESUMO

Utilizaram-se dados de 689 lactações de 225 cabras Pardo Alpina x Gurguéia filhas de 18 bodes, controladas de 1988 a 1996, criadas na Fazenda Carnaúba, Taperoá-PB, com os objetivos de estudar os efeitos ambientais e estimar parâmetros genéticos para as características produção de leite (PL) e duração da lactação (DL). Os efeitos ambientais foram estimados por meio de um modelo estatístico que incluiu o efeito aleatório de bode, os efeitos fixos de grupo contemporâneo, tipo de parto e idade da cabra ao parto e duração da lactação, para PL, como covariáveis. Os parâmetros genéticos foram estimados por intermédio do programa MTDFREML, sob modelo animal, contendo os efeitos fixos e covariáveis utilizados no modelo anterior. As características foram influenciadas pelo efeito fixo de grupo contemporâneo. A idade da cabra influenciou apenas a característica PL. As médias obtidas para PL e DL foram 257,71 ± 54,73 kg e 185,56 ± 46,87 dias, respectivamente, equivalendo à produção média diária de 1,38 kg de leite/dia. As estimativas de herdabilidade e repetibilidade foram, respectivamente, de 0,23 e 0,35 para PL e 0,20 e 0,24 para DL. A correlação genética estimada entre as características foi de 0,96. Em razão de as herdabilidades das características serem de média magnitude e limitadas ao sexo, sugere-se o teste de progênie como método de seleção mais eficaz para os bodes.

Palavras-chave: caprinos leiteiros, nordeste, parâmetros genéticos, semi-árido

ABSTRACT

Records of 689 lactations of 225 Alpine x Gurguéia goats, progenies of 18 buck, controlled from 1988 to 1996, raised at Carnaúba Farm, Taperoá-PB, were analyzed with the objective of study environmental effects and estimate genetics parameters for Milk Production (MP) and Length of Lactation (LL). The environmental effects were analyzed by an statistical model that included buck random effects and the fixed effects of contemporary groups, type of kidding and age of dam at kidding and length lactation as covariate. The genetic parameters were estimated using MTDFREML program, under animal model, including fixed effects and covariates through a described model before. The traits were influenced by contemporary groups. The age of dam at kidding influenced only MP trait. The means for MP and LL were 257.72 ± 54.73 kg and 185.54± 46.87 days, respectively, giving a daily mean production of 1.38 kg of milk/day. The heritability and repeatability estimates were, respectively, 0.23 and 0.35 for MP and 0.20 and 0.24 for LL. The genetic correlation estimate between MP and LL was 0.96. Due to the traits heritability are of average magnitude and sex-limited, it is suggested to use the progeny test as method of more effective selection for bucks.

Key words: genetic parameters, goats milk, Northeast, semi-arid

Introdução

O nordeste detém o maior rebanho caprino do país, com, aproximadamente, 93,7% do rebanho nacional. Em termos de tamanho populacional, a Paraíba é o quinto estado nordestino, com cerca de 403.801 cabeças, o equivalente a 6,1% da população caprina desta região. Quanto à densidade populacional, a Paraíba apresenta uma das maiores concentrações de caprinos do nordeste, destacando-se a microrregião do Cariri Ocidental, com densidade de 20 cabras/km² (IBGE, 1996).

Apesar de numericamente expressivo, o rebanho caprino da região nordeste mantém índices produtivos ainda baixos, principalmente em razão do baixo padrão tecnológico empregado na região, que é fruto do sistema de produção adotado pela maioria dos criadores (cultura de subsistência) e da falta de apoio governamental.

A partir da iniciativa de alguns criadores dotados de maior visão empresarial, a caprinocultura leiteira vem se desenvolvendo na região e tem-se mostrado uma atividade promissora. Recentemente, o trabalho integrado da Associação de Criadores de Caprinos do Rio Grande do Norte, amparado por efetivo apoio governamental, por meio do Programa do Leite, fez com que, em apenas três anos, a produção de leite caprino naquele estado alcançasse 2.200.000 litros anuais (Cordeiro, 2001).

Faltam ainda, no entanto, estudos que definam melhor os caminhos para a consolidação do negócio do leite caprino. Na área de melhoramento animal, fundamental para a sustentação do sistema de produção, pouco foi feito. A porcentagem de criadores que adotam a prática de controle leiteiro no Brasil e, especialmente, no Nordeste, é aquém do mínimo necessário para iniciar qualquer programa de seleção (Pimenta Filho et al., 2000). São escassas as estimativas de herdabilidade das características de produção de leite referentes à população regional, essenciais o melhor procedimento em programas de seleção. Dessa forma, os poucos trabalhos aqui realizados são referenciados com informações de outras regiões de clima tropical e temperado, havendo, então, a necessidade de realização de pesquisas na área de melhoramento genético com a população regional.

Estudos com raças nativas em regiões tropicais registraram valores de herdabilidade para produção de leite entre 0,15 e 0,35 (Mavrogenis et al., 1984; Constantinou et al., 1985). No Brasil, a maioria das pesquisas tem sido realizada com cabras originadas de clima temperado, como os trabalhos realizados por Gonçalves et al. (2001), com cabras das raças Saanen, Alpina e Toggenburg, que obtiveram estimativa de herdabilidade para produção de leite de 0,29 ± 0,07, e Soares Filho et al. (2001), que, ao trabalharem com as mesmas raças na região centro-oeste do país, encontraram herdabilidade de 0,06; 0,12 e 0,11, respectivamente.

Embora a duração da lactação seja uma característica de produção, são poucos os autores que apresentam estimativas de sua herdabilidade: Ribeiro (1997), em animais da raça Saanen, verificou herdabilidade de 0,06 e Soares Filho et al. (2001), em cabras Saanen, Alpinas e Toggenburg, encontraram, respectivamente, 0,07; 0,07 e 0,03 de herdabilidade.

A correlação entre duas características de importância econômica deve ser estimada quando se deseja conhecer o grau de dependência entre as características. Essas estimativas, para caprinos, são escassas na literatura. As correlações genética e fenotípica estimadas por Kennedy et al. (1982) entre produção de leite e duração da lactação foram 0,78 e 0,66 respectivamente, para cabras das raças Alpina, Saanen e Toggenburg. Para essas raças, no Brasil, Gonçalves (1996) obteve 0,81 e 0,62 para as correlações genética e fenotípica, enquanto Ribeiro (1997) registrou 0,43 e -0,50, respectivamente, em rebanho Saanen.

Objetivou-se, com este estudo, analisar a influência de alguns efeitos de ambiente e estimar parâmetros genéticos para as características produção de leite e duração da lactação de caprinos mestiços criados no semi-árido paraibano, utilizando a metodologia da Máxima Verossimilhança Restrita sob modelo animal.

Material e Métodos

Foram utilizados dados referentes a lactações de cabras mestiças Pardo Alpina x Gurguéia, controladas no período de 1988 a 1998 na Fazenda Carnaúba, em Taperoá, Estado da Paraíba. A fazenda está localizada na microrregião do Cariri Ocidental, a 7º 12' 23'' de latitude sul e 36º 49' 25'' de longitude W.G., a 500 m de altitude.

O clima da região está incluído no subtipo B, sub-desértico quente, de tendência tropical, com precipitação média de 300 mm anuais, e apresenta curta estação chuvosa, com estação seca prolongada (maior que oito meses) e umidade relativa do ar em torno de 65%.

Os animais foram criados em sistema semi-extensivo e alimentados com capim-buffel (Cenchrus ciliaris L.) e pasto nativo, na estação chuvosa, e capim-elefante (Pennisetum purpureum schum.), palma forrageira (Opuntia sp.), raspa de mandioca (Manihot sativa), bagaço de cana hidrolisado (Sacharum oficcinarum) e concentrado protéico, na estação seca, além de sal mineral à vontade. Utilizou-se o sistema de monta natural controlada e controle sanitário, sistematicamente. O controle leiteiro foi realizado a cada 28 dias, em duas ordenhas diárias, com intervalo de 12 horas. A partir das fichas de produção e reprodução, foi editado um arquivo contendo o número da cabra, paternidade e maternidade, data de nascimento, datas dos partos, produção de leite, duração da lactação, tipo e ordem de parto e intervalo de partos.

Inicialmente, o arquivo continha 1.939 registros de produção de leite. Para gerar o arquivo utilizado nas análises, foram feitas algumas restrições, permanecendo no arquivo cabras paridas entre 1988 e 1996, animais com pai e mãe conhecidos, bodes com, no mínimo, cinco filhas no rebanho, grupos de contemporâneos com, no mínimo, dois animais (cabras que produziram mais de 90 e menos de 500 kg de leite por lactação, com durações de mais de 70 e menos de 400 dias e máximo de 12 parições). Após os descartes, ficaram disponíveis para a análise 689 registros de produção de leite e duração da lactação de 255 cabras filhas de 18 reprodutores e 142 matrizes.

O estudo dos fatores de meio que podem influenciar as características produção de leite (PL) e duração da lactação (DL) foi feito pelo método dos quadrados mínimos, por meio do PROC GLM do SAS^® (1999), utilizando o comando RANDOM, que considera o efeito de reprodutor como aleatório, a partir do seguinte modelo estatístico:

em que: Y_ijkl = PL do animal l, filha do bode i, pertencente ao grupo contemporâneo j, e com tipo de parição k; m = média geral da população; R_i = efeito aleatório do bode i; GC_k = efeito fixo de grupo contemporâneo k; T_k = efeito fixo de tipo de parto, sendo 1 correspondente ao parto simples e 2 ao múltiplo; I_ijkl = efeito da idade da cabra ao parto, como covariável; = média da covariável; b₁ = coeficiente de regressão linear para idade da cabra ao parto; b₂ = coeficiente de regressão quadrático para idade da cabra ao parto; DL_ijkl = efeito da duração da lactação, como covariável; b₃ = coeficiente de regressão linear para duração da lactação; b₄ = coeficiente de regressão quadrática para duração da lactação; = média da covariável; erro aleatório independentemente distribuído com média zero e variância s².

Os grupos de contemporâneos foram formados pelo ano de parição (1988, ..., 1996) mais a estação de parição, sendo a estação 1 correspondente ao período chuvoso (fevereiro a junho) e a estação 2, ao período seco (julho a janeiro do ano seguinte).

O modelo para análise da DL difere do anterior apenas, por não incluir a covariável duração de lactação.

As estimativas dos componentes de (co)variância necessárias para estimação dos parâmetros genéticos foram obtidas pela metodologia da Máxima Verossimilhança Restrita Livre de Derivadas, sob modelo animal em análise bicaracterística, utilizando o programa MTDFREML, descrito por Boldman et al. (1995). O programa utiliza o algoritmo simplex para localizar o mínimo de -2 log_e L (L = função de verossimilhança) e os componentes de (co)variância que minimizam -2 log_e L são estimativas de máxima verossimilhança, que maximizam a função de verossimilhança.

O critério de convergência utilizado foi a variância dos valores do simplex inferior a 10^-9. Após a convergência, o programa era reiniciado usando as estimativas obtidas anteriormente como valores iniciais. Este procedimento foi repetido até que a diferença entre as estimativas das duas últimas convergências não diferisse. Os valores iniciais foram os obtidos com base na análise de variância efetuada pelo método dos quadrados mínimos.

O modelo de análise pode ser representado, matricialmente, como segue:

Y = X_b + Z₁d + Z₂p + e

em que: Y = vetor das observações n x 1; X = matriz n x f de incidência, relacionando as observações aos efeitos fixos; b = vetor f x 1 de efeitos fixos; Z₁ = matriz n x N de incidência, relacionando as observações aos efeitos genéticos aditivos diretos; d = vetor N x 1 de efeitos genéticos aditivos diretos, ~ N(0, ); Z₂ = matriz n x N de incidência, relacionando as observações aos efeitos de ambientes permanentes; p = vetor N x 1 de efeitos de ambientes permanentes, ~ N(0, ); e = vetor n x 1 de efeitos residuais, ~ N(0, ). Sendo: A = matriz de parentesco; I = matriz identidade de ordem n; = variância genética aditiva direta; = variância de ambiente permanente; = variância residual; n = número de observações de produção de leite e duração da lactação; f = número de níveis de efeitos fixos; e N = número de animais com observação.

Os efeitos fixos e covariáveis considerados nesste modelo foram os mesmos utilizados no modelo das análises dos efeitos de ambiente.

Resultados e Discussão

As médias encontradas para produção de leite e duração de lactação foram, respectivamente, 257,7 ± 54,7 kg de leite e 185,5 ± 46,9 dias de duração da lactação, resultando em produção média diária de 1,38 kg de leite. Esse resultado é satisfatório, considerando-se a região onde a fazenda está localizada e o período estudado.

O resumo da análise de variância para as características estudadas é apresentado na Tabela 1. Todos os efeitos foram significativos (P<0,05), com exceção do tipo de parto para PL e DL e da idade da cabra ao parto para DL.

A influência de grupo contemporâneo já era esperada, pois as flutuações do clima, especificamente a variação de pluviosidade entre e dentro dos anos estudados causa variações no desempenho por meio da estacionalidade na disponibilidade de alimentos, em termos qualiquantitativos, mesmo com o manejo alimentar utilizado na Fazenda.

A não-significância do tipo de parto sobre as características estudadas pode ser atribuída ao manejo dos cabritos, que são aleitados artificialmente. Dessa forma, fica evidente a falta de estímulo diferenciado, em que as cabras, mesmo de tipos de partos diferentes, passam, após a parição, receber o mesmo estímulo (ordenha), corroborando os relatos de Soares Filho et al. (2001).

O resultado obtido para idade da cabra ao parto em relação à PL é biologicamente correto, uma vez que a cabra aumenta sua produção até certa idade e tende a decrescer, descrevendo uma função quadrática, o que está de acordo com a maioria dos trabalhos encontrados na literatura, como os de Mourad (1992) e Gonçalves et al. (2001). A equação e respectiva curva de regressão da produção de leite, em função da idade ao parto, encontram-se na Figura 1. A idade em que as cabras apresentaram máxima produção foi, aproximadamente, aos 56,7 meses (4,7 anos).

A duração da lactação influenciou a produção de leite de forma quadrática (Tabela 1). A regressão e a respectiva equação são apresentadas na Figura 2. Houve relação crescente entre produção de leite e duração da lactação. Entretanto, esta relação, acontece até certo limite, evidenciando que lactações mais longas não favorecem maior produção de leite, ou seja, a produção de leite não aumenta linearmente com o avanço do período de lactação.

As estimativas de herdabilidade para PL e DL foram de magnitude média (Tabela 2). A herdabilidade estimada para PL encontra-se dentro do intervalo das estimativas observadas na literatura para cabras de origem tropical, próxima à relatada por Gonçalves et al. (2001) e superior às encontradas por Soares Filho et al. (2001). A estimativa de herdabilidade para DL encontra-se acima dos valores obtidos para cabras especializadas, criadas em condições brasileiras. Essas estimativas de herdabilidade indicam correlação média entre o fenótipo e genótipo dos indivíduos. Portanto, a seleção de fêmeas pelo fenótipo promoverá lento incremento genético na população. Assim, a partir dos resultados obtidos e pelo fato de as características serem limitadas ao sexo, para os reprodutores, o teste de progênie é o método indicado como o mais eficaz para seleção dos animais destinados à reprodução.

A seleção para DL só seria interessante em programas de melhoramento se a meta almejada fosse maximizar a eficiência produtiva e reprodutiva do rebanho. Assim, a seleção seria praticada com base nas cabras que apresentassem maior produção em menor período de lactação possível. Como foi constatado neste estudo, a relação entre DL e PL não é linear, porém a relação entre DL e intervalo de parto, constatada por Sarmento et al. (2003), é linear, ou seja, cada dia a mais na DL correspondeu a um aumento de 0,59 dias no intervalo de parto. Dessa forma, a seleção de cabras com maiores períodos lactantes refletirá em maiores intervalos de partos, conseqüentemente menor pressão de seleção e ganho genético anual.

As estimativas de repetibilidade para as características estudadas foram de média magnitude (Tabela 2), de modo que uma única medida no animal pode não representar a real capacidade de produção, sendo necessário, portanto, mais de uma medida para aumentar a acurácia das estimativas. Estimativas inferiores foram registradas por Ribeiro (1997), com cabras de raças especializadas no Brasil, para ambas as características.

A correlação genética estimada entre PL e DL foi alta (Tabela 2), indicando que grande parte dos efeitos genéticos aditivos age sinergicamente em ambas as características, resultado que sugere que a seleção de cabras com maior produção de leite levaria a lactações mais longas. Entretanto, a resposta correlacionada não depende só da associação genética entre as características, como também das herdabilidades das características. Dessa forma, a resposta correlacionada pela seleção indireta, embora eficiente, em decorrência da alta correlação genética, resultará em baixos ganhos genéticos. A correlação estimada neste estudo encontra-se próxima às verificadas por Kennedy et al. (1982) e Gonçalves (1996), em cabras de raças exóticas.

Conclusões

Melhorias no manejo alimentar devem ser implementadas, principalmente no que tange ao armazenamento de forragem, visando diminuir os efeitos indiretos do clima.

Em razão de as herdabilidades das características PL e DL serem de média magnitude e por estas serem limitadas ao sexo, indica-se o teste de progênie como o método mais eficaz para seleção de bodes.

Apesar da alta correlação estimada entre produção de leite e duração da lactação, o ganho genético pela seleção indireta é pequeno, pois as herdabilidades das características são relativamente baixas.

Literatura Citada

Recebido em: 09/12/02

Aceito em: 17/03/03

• BOLDMAN, K.G.; KRIESE, L.A.; Van VLECK, D.L. et al. A manual for use of MTDFREML A set of programs to obtain estimates of variances and covariances [DRAFT] Lincon: USDA/ARS, 1995. 120p.
• CONSTANTINOU, A.; BEUING, R.; MAVROGENIS, A.P. Genetic and phenotypic parameters for some reproductive and milk prodution characters of the Damascus goat. Zeitschrift fur Tierzuchtung und Zuchtungsbiologie, v.102, p.301-307, 1985.
• CORDEIRO, P.R.C. Produção de leite d cabra no Brasil. In: REUNIÃO ANUAL DA SOCIEDADE BRASILEIRA DE ZOOTECNIA, 38., 2001, Piracicaba. Anais... Piracicaba: Sociedade Brasileira de Zootecnia, 2001. p.497.
• GONÇALVES, H.C.; SILVA, M.A.; WECHSLER, F.S. et al. Fatores genéticos e de meio na produção de leite de caprinos leiteiros. Revista Brasileira de Zootecnia, v.30, n.3, p.719-729, 2001.
• GONÇALVES, H.C. Fatores genéticos e de meio em algumas características produtivas e reprodutivas de caprinos Viçosa, MG: Universidade Federal de Viçosa, 1996. 141p. Tese (Doutorado em Zootecnia) - Universidade Federal de Viçosa, 1996.
• KENNEDY, B.W.; FINLEY, C.M.; BRADFORD, G.E. Phenotypic and genetic relationships between reproduction and milk production in dairy goats. Journal of Dairy Science, v.65, p.2373-83, 1982.
• MAVROGENIS, A.P.; CONSTATINOU, A.; LOUCA, A. Environmental and genetic causes of variation in production traits of Damascus goats. Animal Production Science, v.33, p.99-104, 1984.
• MOURAD, M. Effects of month of kidding, parity and litter size on milk yield of Alpine goats in Egypt. Small Ruminants Research, v.8, p.545-554, 1992.
• PIMENTA FILHO, E.C.; RIBEIRO, M.N.; SOUSA, W.H. Melhoramento genético de pequenos ruminantes para carne e leite. In: Congresso Nordestino de Produção Animal, 2., 2000, Teresina. Anais... Teresina: Sociedade Nordestina de Produção Animal, 2000. p.107-116.
• RIBEIRO, A.C. Estudo dos efeitos genéticos e de ambiente sobre características de importância econômica em caprinos da raça Saanen Jaboticabal: Universidade Estadual Paulista, 1997. 85p. Dissertação (Mestrado em Zootecnia) - Universidade Estadual Paulista, 1997.
• SARMENTO, J.L.R.; PIMENTA FILHO, E.C.; RIBEIRO, M.N. et al. Fatores genéticos e ambientais sobre o intervalo de partos de cabras leiteiras no semi-árido nordestino. Revista Brasileira de Zootecnia, v.32, n.4, p.875-879, 2003.
• STATISTICAL ANALYSIS SYSTEM - SAS. User's guide: Statistics. Version 8.0, Cary: 1999.
• SOARES FILHO, G.; MCMANUS, C.; MARIANTE, A.S. Fatores genéticos e ambientais que influenciam algumas características de reprodução e produção de leite em cabras no Distrito Federal. Revista Brasileira de Zootecnia, v.30, n.1, p.53-59, 2001.

Datas de Publicação

• Publicação nesta coleção
  22 Ago 2005
• Data do Fascículo
  Dez 2004

Histórico

• Aceito
  17 Mar 2003
• Recebido
  09 Dez 2002