Acessibilidade / Reportar erro

Estrutura genética da raça Sindi no Brasil

Genetic structure of the red Sindhi breed in Brazil

Resumos

O objetivo deste trabalho foi descrever a estrutura genética da raça Sindi no Brasil, utilizando dados do registro genealógico de animais nascidos entre 1955 e 1998. O banco de dados foi separado nos seguintes períodos: 1979-1983, 1984-1988, 1989-1993 e 1994-1998. A endogamia total aumentou de 0,38 para 10,13%; a esperada sob acasalamento ao acaso, de 0,07 para 5,65%; e a endogamia atribuída à subdivisão populacional, de 0,30 para 4,74%, considerando os períodos inicial e final, indicando a existência de subdivisão genética na raça Sindi. O tamanho efetivo populacional estimado pelo coeficiente total de endogamia decresceu vertiginosamente de 161 para 9. Tendo-se como base a probabilidade de origem do gene, foram calculados os números efetivos de fundadores, de ancestrais e de genomas remanescentes, os quais decresceram, ao longo do período total, atingindo valores de 16, 9 e 7, respectivamente.

estrutura genética; probabilidade de origem do gene; Sindi; tamanho efetivo populacional


The aim of this research was to describe the genetic structure of the Brazilian Red Sindhi breed using pedigree records from registered animals born from 1955 to 1998. The data file was separated in the following periods: 1979-1983, 1984-1988, 1989-1993 and 1994-1998. The total inbreeding increased from 0.38 to 10.13%, the expected inbreeding under random mating increased from 0.07 to 5.65%, and the inbreeding due to population subdivision increased from 0.30 to 4.74%, from the first to the final period, indicating that genetic subdivision in the Red Sindhi breed is relevant. The effective population size was estimated from the total inbreeding and decreased vertiginously from 161 to 9. Based on probabilities of gene origin the effective numbers of founders, ancestors and remaining genomes were calculated and decreased over the total period, reaching values of 16, 9 and 7, respectively.

genetic structure; probability of gene origin; Red Sindhi; effective population size


MELHORAMENTO, GENÉTICA E REPRODUÇÃO

Estrutura genética da raça Sindi no Brasil

Genetic structure of the red Sindhi breed in Brazil

Fábio José Carvalho FariaI; Anibal Eugênio Vercesi FilhoII; Fernando Enrique MadalenaIII; Luiz Antônio JosahkianIV

IProfessor Visitante do Depto. de Produção Animal da UFMS, Cidade Universitária s/n, Cx. P. 549, CEP 79070-900, Campo Grande - MS. E-mail: cifal_ms@yahoo.com.br

IIAluno de Doutorado em Ciência Animal, Departamento de Zootecnia, Escola de Veterinária da UFMG, Av. Antônio Carlos, 6627, Cx. P. 567, CEP 30123-970, Belo Horizonte - MG

IIIProfessor do Departamento de Zootecnia, Escola de Veterinária da UFMG, Av. Antônio Carlos 6627, Cx. P. 567, CEP 30123-970, Belo Horizonte - MG. E-mail: fermadal@dedalus.lcc.ufmg.br

IVSuperintendente técnico da Associação Brasileira de Criadores de Zebu, Praça Vicente Rodrigues da Cunha, 110 Bl 1, CEP 38022-330, Uberaba - MG. Professor de Melhoramento Animal da Faculdade de Agronomia e Zootecnia de Uberaba - FAZU

RESUMO

O objetivo deste trabalho foi descrever a estrutura genética da raça Sindi no Brasil, utilizando dados do registro genealógico de animais nascidos entre 1955 e 1998. O banco de dados foi separado nos seguintes períodos: 1979-1983, 1984-1988, 1989-1993 e 1994-1998. A endogamia total aumentou de 0,38 para 10,13%; a esperada sob acasalamento ao acaso, de 0,07 para 5,65%; e a endogamia atribuída à subdivisão populacional, de 0,30 para 4,74%, considerando os períodos inicial e final, indicando a existência de subdivisão genética na raça Sindi. O tamanho efetivo populacional estimado pelo coeficiente total de endogamia decresceu vertiginosamente de 161 para 9. Tendo-se como base a probabilidade de origem do gene, foram calculados os números efetivos de fundadores, de ancestrais e de genomas remanescentes, os quais decresceram, ao longo do período total, atingindo valores de 16, 9 e 7, respectivamente.

Palavras-chave: estrutura genética, probabilidade de origem do gene, Sindi, tamanho efetivo populacional

ABSTRACT

The aim of this research was to describe the genetic structure of the Brazilian Red Sindhi breed using pedigree records from registered animals born from 1955 to 1998. The data file was separated in the following periods: 1979-1983, 1984-1988, 1989-1993 and 1994-1998. The total inbreeding increased from 0.38 to 10.13%, the expected inbreeding under random mating increased from 0.07 to 5.65%, and the inbreeding due to population subdivision increased from 0.30 to 4.74%, from the first to the final period, indicating that genetic subdivision in the Red Sindhi breed is relevant. The effective population size was estimated from the total inbreeding and decreased vertiginously from 161 to 9. Based on probabilities of gene origin the effective numbers of founders, ancestors and remaining genomes were calculated and decreased over the total period, reaching values of 16, 9 and 7, respectively.

Key Words: genetic structure, probability of gene origin, Red Sindhi, effective population size

Introdução

A raça Sindi é originária do Estado de Sind no Paquistão e, graças à sua reputação em rusticidade e tolerância ao calor, se espalhou pela Ásia, Oceania, África e Américas. Os animais Sindi caracterizam-se por pequeno porte e pelagem vermelha, que pode variar em tonalidades, de acordo com o sexo do animal. As fêmeas produzem, em média, 1.700 kg de leite por lactação. Entretanto, sob condições ótimas, algumas chegam a ultrapassar 4.000 kg de leite por lactação (Leite et al., 2001).

Os primeiros exemplares dessa raça chegaram ao Brasil em 1952 e, atualmente, a maior parte do rebanho está presente no nordeste brasileiro (Leite et al., 2001). Desde sua introdução, o gado Sindi brasileiro manteve-se concentrado em poucos rebanhos, não apresentando a evolução numérica verificada em outras raças, mantendo o rebanho nacional reduzido. De acordo com Faria et al. (2001), o número, estimado em menos de 1.000 fêmeas puras em reprodução, classifica a raça no patamar crítico em relação ao risco de extinção, conforme critérios adotados pela FAO - Food and Agriculture Organization of the United Nations.

Este trabalho foi conduzido com o objetivo de quantificar a variabilidade genética da raça Sindi, por meio de informações de seu registro genealógico, com o intuito de fornecer subsídios para programas de preservação e expansão da raça.

Material e Métodos

Foram utilizados os dados do registro genealógico da raça Sindi, disponíveis sob processamento eletrônico, da Associação Brasileira de Criadores de Zebu (ABCZ). O arquivo de dados consistiu de 4.884 observações entre os anos de 1955 a 1998, referentes a amostragem de 8.080 registros (ABCZ, 2001). As informações referentes a cada animal consistiam de: pai, mãe, sexo, categoria de registro (PO = puro de origem / LA = livro aberto), dia, mês e ano de nascimento, criador e proprietário. Inicialmente, foram geradas estatísticas descritivas de cada amostra populacional e, em seguida, procedeu-se à separação do banco de dados em períodos, compreendendo os anos de 1979-1983, 1984-1988, 1989-1993 e 1994-1998. O tamanho efetivo populacional (Ne) foi estimado por meio da taxa de variação dos coeficientes médios de endogamia entre gerações sucessivas, (DF), conforme Falconer & Mackay (1996):

em que: Ft= coeficiente médio de endogamia estimado na geração atual; Ft-1= coeficiente médio de endogamia estimado na geração anterior.

A estimação dos parâmetros baseados na probabilidade de origem do gene e o coeficiente de endogamia foram obtidos por meio da utilização do pacote PEDIGÒ, que apresenta vários programas escritos em linguagem FORTRAN, gentilmente cedido pelo Dr. Didier Boichard (INRA - Station de Génétique Quantitative et Appliquée). O número efetivo de fundadores representa o número de animais com igual contribuição que produziria a mesma variabilidade genética encontrada na população estudada (Lacy, 1989). A aproximação clássica na determinação do número efetivo de fundadores foi obtida por:

em que Nfun = número efetivo de fundadores;

= somatório da contribuição esperada do número de progênies (q) do fundador k na população.

O número efetivo de ancestrais representa o número mínimo de animais (fundadores ou não) necessário para explicar a total diversidade genética da população estudada (Boichard et al., 1997). A determinação do número efetivo de ancestrais foi obtida computando-se a contribuição marginal de cada ancestral por:

em que Nanc = número efetivo de ancestrais;

= contribuição marginal do ancestral k na população.

O número efetivo de genomas remanescentes representa o número de fundadores com igual contribuição que não teria perda de alelos fundadores por deriva genética na progênie e que produziria a mesma diversidade genética que a encontrada na população estudada (Lacy, 1989). O programa simulou segregações para cada progênie, gerando genótipos em função da amostragem aleatória de um alelo em cada progenitor, e as freqüências fk foram determinadas por contagem na população referência. O Ngen foi determinado em função da média de 1.000 simulações de segregações gênicas.

em que Ngen = número efetivo de genomas remanescentes; Na = número efetivo de alelos (Crow & Kimura, 1970).

Foram calculados o número médio de ancestrais, a proporção de ancestrais identificados por geração e o número equivalente de gerações completas, definido pela soma entre todas as gerações da proporção de ancestrais identificados em cada geração. Todas essas estatísticas foram computadas para indivíduos com pelo menos um ancestral identificado. O número equivalente de gerações completas foi dado por:

em que nj = total de ancestrais do animal j;gij = número de gerações entre o indivíduo j e seu ancestral i; N = total de animais na população-referência.

As estatísticas de F (FIT, FST e FIS) de Wright (1965) também foram calculadas. FIT é o coeficiente médio de endogamia da população, FST designa o coeficiente médio de endogamia esperado se os reprodutores de cada período fossem acasalados aleatoriamente e foi estimado em populações hipotéticas produzidas pelo acasalamento ao acaso dos pais dentro de cada período por meio do programa escrito em linguagem FORTRAN, gentilmente cedido pelo Dr. Tetsuro Nomura (Kyoto Sangyo University), e FIS expressa o desvio da casualidade nos acasalamentos, obtido como:

Se FIS > 0, então a endogamia real ultrapassa o nível esperado sob acasalamentos ao acaso, implicando que os acasalamentos estão acontecendo entre animais mais aparentados que a média, conseqüentemente a população estará dividida em subpopulações. Quando FIS < 0, evita-se a endogamia ou predomina o acasalamento entre subpopulações.

Resultados e Discussão

Determinaram-se o coeficiente de endogamia, a taxa de incremento (DF) e o Ne, que podem ser visualizados na Figura 1. Os valores calculados do primeiro ao último período foram de 0,38; 3,16; 6,20 e 10,13% para o FIT, que expressa o coeficiente de endogamia; de 0,3; 2,3; 3,5 e 5,1% para o DF, que expressa a taxa de incremento da endogamia; e de 161, 20, 14 e 9 para o Ne. A endogamia e a taxa de incremento foram crescentes do primeiro ao último período, ao passo que o Ne decresceu vertiginosamente, o que traz preocupações quanto a manutenção da variabilidade genética dessa raça.


Dos 4.981 animais do arquivo de genealogia, 1.611 apresentaram coeficiente de endogamia positivo, com média de 13,23% e máximo de 42,96%. Na Tabela 1, são apresentados os números de animais distribuídos por classe de endogamia.

O coeficiente total de endogamia apresentou-se em nível crítico, atingindo no último período 10,13%, com crescimento da ordem de 33 vezes em relação ao valor do primeiro período. O coeficiente de regressão linear do F médio anual, ponderado pelo número de animais, foi estimado em 0,007 (P<0,05). O gráfico dos valores observados e preditos do coeficiente de endogamia podem ser visualizados na Figura 2. Os valores dos parâmetros que expressam a variabilidade genética descritos para essa raça são críticos, e puderam ser confirmados em seguida por meio das análises baseadas na probabilidade de origem do gene.


Goddard & Smith (1990) sugerem Ne mínimo de 40 por geração para maximizar o ganho genético em população selecionada para produção de leite, e Meuwissen & Woolliams (1994) recomendam valores entre 30 e 250 para prevenir decréscimo no valor adaptativo. Nicholas (1980) propõe valor mínimo de 50 por geração para o Ne. Nesse contexto, esse grupamento apresenta valores assaz reduzidos em comparação às recomendações da literatura.

Na Figura 3, pode-se observar a evolução anual do número de indivíduos endogâmicos em relação ao ano de nascimento na raça Sindi. Nesse período, pode-se verificar que, em 1979, não foram registrados animais endogâmicos. Entretanto, nos últimos três anos a quase totalidade dos animais é endogâmica.


A evolução das estatísticas de F pode ser vista na Figura 4. O FST passou de 0,07%, no primeiro período, para 5,65% no último, aumentando em 7.300% seu valor, sendo este incremento resultado do aumento da relação de parentesco entre os progenitores, o que reflete acasalamentos preferenciais ocorridos no passado. O FIS passou de 0,30% para 4,74%, valor que indica a existência de subdivisão na raça. Vercesi Filho et al. (2002) também relataram a existência de subdivisão populacional na raça Tabapuã, embora os valores encontrados estejam aquém dos aqui obtidos. Resultados opostos foram encontrados por Nomura et al. (2001), na raça Japanese Black, e por Faria (2002) nas raças Nelore e Nelore Mocho. Esses autores relataram a ausência de subdivisão nessas raças, ocasionada por alguns fatores, dentre esses o uso mais intenso de IA.


O arquivo de genealogia permitiu traçar o número médio de ancestrais por indivíduo e o número de gerações. Embora o arquivo trabalhado contivesse dados a partir de 1955, as informações citadas foram obtidas de indivíduos com pelo menos um ancestral conhecido. Em 1963 foram registradas duas fêmeas com média de 0,75 ancestral identificado, enquanto em 1998 esse número foi de 85 fêmeas registradas com cerca de 28 ancestrais e 3,7 gerações equivalentes, em média, identificadas por animal. No caso dos machos, essas estimativas não foram muito diferentes; em 1965, foi registrado um animal com apenas um ancestral identificado, ao passo que, no ano de 1998, 83 machos com número médio de gerações equivalentes e ancestrais similares aos encontrados nas fêmeas. O número máximo de gerações identificadas foi de seis para ambos os sexos. Faria (2002) encontrou raças zebuínas em que se pode traçar até 13 gerações ascendentes.

Finalizando a análise estrutural, seguem descritos na Tabela 2 os números efetivos de fundadores, de ancestrais e de genomas remanescentes calculados para essa população. Segundo Boichard et al. (1997), as análises realizadas dos parâmetros baseados na probabilidade de origem do gene são mais eficientes na predição da variabilidade genética porque são menos sensíveis a falhas ou perdas de informações dos pedigrees. Constata-se, na Tabela 2, que houve queda na magnitude dos parâmetros estimados entre os períodos, exprimindo redução na variabilidade, com a seguinte ordem de magnitude Nfun > Nanc > Ngen. O Nfun avalia como o balanço esperado das contribuições dos fundadores é mantido entre as gerações, levando em consideração a taxa de seleção e a variação de tamanho de família. Uma das causas de perda de genes nas populações é a diminuição do número de reprodutores em determinado período, o que é levado em consideração quando da estimação do Nanc, ao passo que o Ngen avalia o número de genes fundadores que é mantido na população para determinado loco. Desse modo, o Nfun apresentou os maiores valores, por expressar somente a contribuição de fundadores, e o Ngen, as menores estimativas, por considerar todas as possibilidades de perdas de alelos fundadores.

O Nfun estimado no presente trabalho é o menor descrito para raças zebuínas criadas no Brasil, descritas por Faria (2002), o que pode ser atribuído ao pequeno número de fundadores e machos em uso. O número efetivo de ancestrais Sindi revela a predominância na utilização de alguns reprodutores, explicada por sua contribuição marginal junto à população-referência (Tabela 3).

A razão Nfun/Nanc expressa o efeito "gargalo" e resulta da diminuição do número de reprodutores entre períodos. No rebanho Sindi, esse valor passa de 1,03 para 1,77, com valor médio de 1,52. De acordo com Vercesi Filho et al. (2002), na raça Tabapuã esse valor esteve entre 1,02 e 1,44 no primeiro e último períodos. Boichard et al. (1997) relataram razão igual a 3,0 na raça Normanda, relacionada à intensa utilização de poucos animais por meio de inseminação artificial. Embora não possam ser consideradas comparações diretas, o Nanc do último período na raça Sindi está entre os menores valores descritos na literatura (Boichard et al., 1997; Sölkner et al., 1998; Vercesi Filho et al., 2002). A proporção entre Ngen/Nfun indica a extensão do processo de deriva genética e os valores reduziram de 0,81 a 0,43, indicando aumento nesse processo ao longo dos períodos. Vercesi Filho et al. (2002) relataram para a raça Tabapuã o valor de 0,54 entre 1994 e 1998.

Na Tabela 3, evidencia-se maior contribuição de certos animais. No último qüinqüênio, 50 animais respondem pela quase totalidade dos alelos presentes na população. Tornou-se crecente, com o passar dos qüinqüênios, a contribuição alélica de cada ancestral, exemplificada entre o primeiro e último períodos nos quais o ancestral que mais contribuiu incrementou sua participação em 500%, detendo, no último qüinqüênio, quase 30% dos alelos na população-referência sendo que, nesse último, os cinco ancestrais mais importantes responderam por cerca de 60% dos genes.

Os resultados encontrados nesta pesquisa corroboram as descrições de redução na variabilidade genética já descritos por Faria et al. (2001) no rebanho Sindi. No referido trabalho, estes autores apresentaram os efeitos da variância de tamanho de progênie sobre o Ne, determinando valor crítico quanto ao número de fêmeas puras em reprodução. No presente estudo, os parâmetros determinados por meio das estatísticas de F e da probabilidade de origem do gene indicaram redução bastante acentuada da variabilidade genética da raça Sindi.

Conclusões

Os valores do coeficiente de consangüinidade (FIT) e da taxa de incremento da endogamia (DF) aumentaram entre os períodos, com concomitante redução do tamanho efetivo populacional (Ne). A relação de parentesco entre os animais elevou-se conforme os valores determinados do coeficiente médio de endogamia esperado sob acasalamentos ao acaso (FST). A magnitude do índice de subdivisão genética (FIS) indicou a existência de subpopulações. Os parâmetros estimados do número efetivo de fundadores (Nfun), de ancestrais (Nanc) e de genomas remanescentes (Ngen) decresceram entre os períodos e as relações entre eles (Nfun/Nanc e Ngen/Nfun) indicaram o aumento do efeito gargalo e do processo de deriva genética na população avaliada.

Agradecimento

À ABCZ, pela concessão dos dados para este trabalho, realizado com apoio do CNPq, CAPES e FAPEMIG.

Literatura Citada

Recebido em: 22/04/03

Aceito em: 09/10/03

  • ABCZ. "Estatísticas", http://www.abcz.org.br, acessado em 10/04/2001.
  • BOICHARD, D.; MAIGNEL, L.; VERRIER, É. The value of using probabilities of gene origin to measure genetic variability in a population. Genetics, Selection and Evolution, v.29, p.5-23, 1997.
  • CROW, J.F.; KIMURA, M. An introduction to population genetics theory New York: Harper and Row, 1970. 591p.
  • FALCONER, D.S.; MACKAY, T. Introduction to quantitative genetics 4.ed. London: Longman Group Ltda, 1996. 464p.
  • FARIA, F.J.C. Estrutura genética das populações zebuínas brasileiras registradas Belo Horizonte: Universidade Federal de Minas Gerais, 2002. 177p. Tese (Doutorado em Ciência Animal) - Universidade Federal de Minas Gerais, 2002.
  • FARIA, F.J.C.; VERCESI FILHO, A.E.; MADALENA, F.E. et al. Parâmetros populacionais do rebanho Sindi registrado no Brasil. Revista Brasileira de Zootecnia, v.30(6S), p.1989-1994, 2001.
  • GODDARD, M.G.; SMITH, C. Optimum number of bull sires in dairy cattle breeding. Journal of Dairy Science, v.73, p.1113-1122, 1990.
  • LACY, R.C. Analysis of founder representation in pedigrees: founder equivalents and founder genome equivalents. Zoo Biology, v.8, p.111-123, 1989.
  • LEITE, P.R.M.; SANTIAGO, A.A.; NAVARRO FILHO, H.R. et al. Sindi: Gado vermelho para o semi-árido. João Pessoa: EMEPA - PB / Banco do Nordeste, 2001. 147p.
  • MEUWISSEN, T.H.E.; WOOLLIAMS, J. Effective sizes of livestock populations to prevent a decline in fitness. Theoretical and Applied Genetics, v.89, p.1019-1026, 1994.
  • NICHOLAS, F.W. Size of population required for artificial selection. Genetical Research, v.35, p.85-105, 1980.
  • NOMURA, T.; HONDA, T.; MUKAI, F. Inbreeding and effective population size of Japanese Black cattle. Journal of Animal Science, v.79, p.366-370, 2001.
  • SÖLKNER, J.; FILIPCIC, L.; HAMPSHIRE, N. Genetic variability of populations and similarity of subpopulations in Austrian cattle breeds determined by analysis of pedigrees. Animal Science, v.67, p.249-256, 1998.
  • VERCESI FILHO, A.E.; FARIA, F.J.C.; MADALENA, F.E. et al. Estrutura populacional do rebanho Tabapuã registrado no Brasil. Arquivo Brasileiro de Medicina Veterinária e Zootecnia, v.54, n.6, p.609-617, 2002.
  • WRIGHT, S. The interpretation of population structure by F-statistics with special regards to systems of mating. Evolution, v.19, p.395-420, 1965.

Datas de Publicação

  • Publicação nesta coleção
    17 Nov 2004
  • Data do Fascículo
    Ago 2004

Histórico

  • Recebido
    22 Abr 2003
  • Aceito
    09 Out 2003
Sociedade Brasileira de Zootecnia Universidade Federal de Viçosa / Departamento de Zootecnia, 36570-900 Viçosa MG Brazil, Tel.: +55 31 3612-4602, +55 31 3612-4612 - Viçosa - MG - Brazil
E-mail: rbz@sbz.org.br