Acessibilidade / Reportar erro

Desempenho e Predição de Híbridos e Análise de Agrupamento de Características de Matrizes de Frangos de Corte

Performance and Prediction and Clustering Analyses of Characteristics of Parental Broiler Lines of Production Traits

Resumos

Este trabalho foi conduzido com cruzamentos (dialelo parcial) entre as linhas do grupo I (LL, TL, ZL, TT, LT, ZT, ZZ, LZ e TZ) e as linhas do grupo II (PP, KP, VP, KK, PK, VK, VV, PV e KV), além de duas linhas (LLc e PPc) controle. O delineamento experimental foi em blocos casualizados consistindo de três blocos, 82 cruzamentos e nove fêmeas por cruzamento. As características de produção, no total de onze, foram: fertilidade na 32ª, 38ª, 42ª e 49ª semanas, eclodibilidade na 32ª semana, peso inicial dos pintinhos, peso no 42º dia, conversão no 42º dia, rendimento de carcaça, peito e gordura abdominal no 42º dia de idade dos frangos. As médias e predição dos híbridos mostraram ser fonte importante de informações para o melhorista. Existe similaridade entre os cruzamentos oriundos da utilização da linha 1, como progenitor.

análise de agrupamento; desempenho; híbridos; matrizes de frango de corte; predição


This study was conducted using data from single crosses (partial diallel) of male lines, group I (LL, TL, ZL, TT, LT, ZT, ZZ, LZ, TL), female lines, group II (PP, KP, VP, KK, PK, VK, VV, PV, KV), besides two control lines (LLc and PPc). A completely randomized block design, with three replicales, 82 single crosses treatments with nine females for experimental unit was used in the analyses. The production traits analyzed, eleven total, were: fertility at 32nd, 38th, 42nd and 49th weeks of age, hatchability at 32nd, body weight at day of age, body weight and feed:gain ratio at 42nd days of age, and carcass yield, breast yield and fast pad yield at 42nd day of age. The averages and prediction of the hybrid ones showed to be important source of information. Similarity exists among the crossings originating from of the use of the line 1, as progenitor.

clustering analyses; broiler parental lines; hybrids; performance; prediction


Desempenho e Predição de Híbridos e Análise de Agrupamento de Características de Matrizes de Frangos de Corte

Valéria Maria Nascimento Abreu2 1 Parte do trabalho de tese de D.S. do primeiro autor, financiada pela CAPES/UFV/EMBRAPA. , Martinho de Almeida e Silva3 1 Parte do trabalho de tese de D.S. do primeiro autor, financiada pela CAPES/UFV/EMBRAPA. , Cosme Damião Cruz4 1 Parte do trabalho de tese de D.S. do primeiro autor, financiada pela CAPES/UFV/EMBRAPA. , Élsio Antônio Pereira de Figueiredo5 1 Parte do trabalho de tese de D.S. do primeiro autor, financiada pela CAPES/UFV/EMBRAPA. , Paulo Giovanni de Abreu6 1 Parte do trabalho de tese de D.S. do primeiro autor, financiada pela CAPES/UFV/EMBRAPA.

RESUMO - Este trabalho foi conduzido com cruzamentos (dialelo parcial) entre as linhas do grupo I (LL, TL, ZL, TT, LT, ZT, ZZ, LZ e TZ) e as linhas do grupo II (PP, KP, VP, KK, PK, VK, VV, PV e KV), além de duas linhas (LLc e PPc) controle. O delineamento experimental foi em blocos casualizados consistindo de três blocos, 82 cruzamentos e nove fêmeas por cruzamento. As características de produção, no total de onze, foram: fertilidade na 32a, 38a, 42a e 49a semanas, eclodibilidade na 32a semana, peso inicial dos pintinhos, peso no 42o dia, conversão no 42o dia, rendimento de carcaça, peito e gordura abdominal no 42o dia de idade dos frangos. As médias e predição dos híbridos mostraram ser fonte importante de informações para o melhorista. Existe similaridade entre os cruzamentos oriundos da utilização da linha 1, como progenitor.

Palavras-chave: análise de agrupamento, desempenho, híbridos, matrizes de frango de corte, predição

Performance and Prediction and Clustering Analyses of Characteristics of Parental Broiler Lines of Production Traits

ABSTRACT - This study was conducted using data from single crosses (partial diallel) of male lines, group I (LL, TL, ZL, TT, LT, ZT, ZZ, LZ, TL), female lines, group II (PP, KP, VP, KK, PK, VK, VV, PV, KV), besides two control lines (LLc and PPc). A completely randomized block design, with three replicales, 82 single crosses treatments with nine females for experimental unit was used in the analyses. The production traits analyzed, eleven total, were: fertility at 32nd, 38th, 42nd and 49th weeks of age, hatchability at 32nd, body weight at day of age, body weight and feed:gain ratio at 42nd days of age, and carcass yield, breast yield and fast pad yield at 42nd day of age. The averages and prediction of the hybrid ones showed to be important source of information. Similarity exists among the crossings originating from of the use of the line 1, as progenitor.

Key Words: clustering analyses, broiler parental lines, hybrids, performance, prediction

Introdução

O estudo do desempenho de híbridos por meio da média de características permite avaliar a possibilidade de uso per si dos híbridos avaliados, considerando para isto, fundamentalmente, a média fenotípica e o grau de complementariedade genética de seus progenitores, expressos pelos efeitos da capacidade específica de combinação. Ao se selecionarem progenitores para compor programa de cruzamentos, aproveitando as informações providas pela capacidade específica de combinação, inicialmente devem-se observar dentro dos progenitores os que apresentam maior média e, depois, os de maior capacidade geral de combinação.

De acordo com Cruz & Regazzi (1994), a predição de híbridos, sejam duplos, triplos ou múltiplos, tem sido de grande utilidade, sendo oriunda das informações originadas pelos cruzamentos em dialelos parciais.

O híbrido triplo é o resultado do cruzamento entre o híbrido simples e uma terceira linhagem. Na situação mais freqüente, o híbrido simples é obtido com as linhagens de um mesmo grupo e o triplo do cruzamento desse híbrido simples com uma linhagem de outro grupo.

O híbrido duplo é o resultado do cruzamento entre dois híbridos simples. Geralmente, os híbridos simples são obtidos do cruzamento entre linhagens de determinado grupo, e o híbrido duplo, do cruzamento entre híbridos simples de grupos distintos.

A análise de agrupamento tem por finalidade reunir, por algum critério de classificação, os progenitores (ou qualquer outro tipo de unidade amostral) em vários grupos, de tal forma que existam homogeneidade dentro do grupo e heterogeneidade entre grupos. Alternativamente, as técnicas de análise de agrupamento têm por objetivo, ainda, dividir um grupo original de observações em vários grupos, segundo algum critério de similaridade ou dissimilaridade. O processo de agrupamento envolve, basicamente, duas etapas, a primeira relaciona-se com a estimação de uma medida de similaridade (ou dissimilaridade) entre os progenitores e a segunda, com a adoção de uma técnica de agrupamento para a formação dos grupos (Cruz & Regazzi, 1994).

Segundo Freitas (1996), a estimação das medidas de dissimilaridade consiste na conversão da matriz de observações n ´ p (sendo n o número de indivíduos ou progenitores e p o de variáveis), em uma matriz quadrada e simétrica, n ´ n, de similaridade individuais, que são medidas da distância entre pares de indivíduos, dado o valor de um conjunto de p variáveis. Na interseção da i-ésima linha e k-ésima coluna dessa matriz, coloca-se a distância entre o i-ésimo e o k-ésimo indivíduo. A alta similaridade (pequena distância) indica que os indivíduos são próximos em relação ao conjunto de variáveis, enquanto que a dissimilaridade (grande distância) indica o contrário.

Nos estudos da divergência genética destinados à identificação de progenitores para hibridação, têm sido de uso mais rotineiro a distância euclidiana média ou a distância generalizada de Mahalanobis, sendo essa última a preferida, entretanto, possível de ser estimada apenas quando se dispõe da matriz de covariâncias residuais estimada a partir de ensaios experimentais com repetições (Cruz & Regazzi, 1994).

Segundo Cruz (1990), a distância euclidiana, ou o seu quadrado, é bastante utilizada quando as unidades para cálculo são escores de componentes principais ou de variáveis canônicas. Quando é estimada a partir de variáveis originais, apresenta a inconveniência de ser influenciada, além do número e da escala, pelas correlações entre aquelas variáveis. Para contornar a influência do número de variáveis sobre as estimativas da distância euclidiana, recomenda-se a utilização da distância euclidiana média; problemas inerentes à escala são contornados pela padronização prévia das variáveis em estudo. Assim, o maior inconveniente é mesmo a pressuposição da independência entre variáveis. Entretanto, Arunachalam (1981) argumenta que, nos estudos sobre divergência genética, vários caracteres com graus significativos de correlação são avaliados, desaconselhando a quantificação da divergência pela distância euclidiana e citando, como maneira alternativa para avaliar a divergência, a estatística D2 de Mahalanobis, que leva em consideração as correlações entre características analisadas por meio da matriz de covariâncias residuais entre variáveis.

A distância generalizada de Mahalanobis considera a variabilidade de cada unidade experimental e não somente a medida de tendência amostral, sendo recomendada para dados provenientes de delineamentos e quando as características são correlacionadas. Considera, também, a correlação entre as características analisadas por meio da matriz de dispersão dessas. Quando as correlações entre as características são nulas e as características padronizadas, a estatística D2 é equivalente à distância euclidiana (Cruz, 1990).

Segundo Cruz & Regazzi (1994), quando se dispõe de várias características, o valor de D2 pode ser, alternativamente, estimado a partir das médias dos dados originais e da matriz de covariâncias residuais, ou a partir de dados transformados. Na transformação dos dados, equivalente ao processo de rotação de eixos, são obtidas combinações lineares das características originais, cujos coeficientes de ponderação são derivados de um processo denominado condensação pivotal. A condensação pivotal é aplicada na matriz de dispersão e evita a inversão de matrizes, que, quando de ordem elevada, proporciona grandes erros numéricos. Com isso, torna os cálculos mais simples (Rao, 1952). Com o avanço dos computadores, já é possível a inversão dessas matrizes sem a necessidade da condensação pivotal. E quando se aplica a condensação pivotal, nas variáveis originais, são estimadas as médias das variáveis transformadas que têm variâncias residuais unitárias e covariâncias residuais nulas, que são variáveis não-correlacionadas, permitindo-se, assim, a estimação de D2 de maneira análoga ao quadrado da distância euclidiana.

A segunda etapa do processo de agrupamento consiste na escolha de um método de agrupamento. Dentre os métodos de agrupamento, os mais utilizados, em estudos de divergência genética, são os hierárquicos e os de otimização, sendo que o método de otimização apresentado por Tocher vem sendo extensivamente usado nesses estudos (Piassi, 1994).

Segundo Cruz & Regazzi (1994), no método de Tocher, adota-se o critério de que a média das medidas de dissimilaridade, dentro de cada grupo, deve ser menor que as distâncias médias entre quaisquer grupos. O método requer a obtenção da matriz de dissimilaridade, sobre a qual é identificado o par de progenitores mais similar. Esses progenitores formarão o grupo inicial e, a partir daí, será avaliada a possibilidade de inclusão de novos progenitores, adotando-se o critério anteriormente citado. A entrada de um progenitor em um grupo sempre aumenta o valor médio da distância dentro do grupo. Assim, pode-se tomar a decisão de incluir o progenitor em um grupo por meio da comparação entre o acréscimo no valor médio da distância dentro do grupo e um nível máximo permitido, que pode ser estabelecido arbitrariamente, ou adotar, como tem sido geralmente feito, o valor máximo da medida de dissimilaridade encontrado das menores distâncias envolvendo cada progenitor.

Este trabalho teve como objetivo obter resultados de desempenho e predição de híbridos resultantes do cruzamento entre linhagens de matrizes de frango de corte, bem como estudar a diversidade genética, por meio da análise de agrupamento, destes cruzamentos.

Material e Métodos

Os dados foram provenientes de experimento realizado no Setor de Melhoramento Genético de Aves (SMGA) do Centro Nacional de Pesquisa de Suínos e Aves (CNPSA), da Empresa Brasileira de Pesquisa Agropecuária (EMBRAPA), em Concórdia, Santa Catarina, de janeiro de 1994 a dezembro de 1995.

Os pintos de um dia de idade deste experimento foram obtidos de cruzamentos (dialelo parcial) entre as linhas do grupo I (LL, TL, ZL, TT, LT, ZT, ZZ, LZ e TZ) e as linhas do grupo II (PP, KP, VP, KK, PK, VK, VV, PV e KV), além de duas linhas (LLc e PPc) controle.

A codificação das linhagens foi feita da seguinte maneira:

As fêmeas foram alojadas em gaiolas individuais perfazendo um total de 2.214 fêmeas. Os machos foram alojados em 60 gaiolas individuais, com piso ripado, para coleta de sêmen, destinado à fertilização posterior, via inseminação artificial, em delineamento experimental em blocos casualizados, consistindo de três blocos, 82 cruzamentos e nove fêmeas por cruzamento.

A inseminação artificial foi efetuada com 0,05 mL de um "pool" de sêmen fresco, não diluído, uma vez por semana, diretamente na vagina da galinha (2 cm de profundidade), com seringa graduada, conforme procedimento padrão utilizado no CNPSA.

Os ovos foram coletados diariamente e armazenados, a 12oC, por sete dias antes da incubação semanal normal para a venda de pintos, com a identificação da galinha, do tratamento e da data de postura. Todos os ovos foram submetidos à ovoscopia aos 19 dias de incubação

As 11 características de produção estudadas, na análise de agrupamento e estudo do desempenho dos híbridos foram: fertilidade na 42a (FERT42); 49a (FERT49); 32a (FERT32) e 38a (FERT38); eclodibilidade na 32a (ECLO32); peso inicial do pintinho (PI); peso dos frangos no 42o dia de idade (P42); conversão alimentar dos frangos no 42o dia de idade (CA); rendimento de carcaça dos frangos no 42o dia de idade (RCAR); rendimento de peito dos frangos no 42o dia de idade (RPEIT); rendimento de gordura abdominal dos frangos no 42o dia de idade (RGA).

Para a predição de híbridos foram utilizadas apenas as características: fertilidade na 38a (FERT38); peso dos frangos no 42o dia de idade (P42) e rendimento de carcaça dos frangos no 42o dia de idade (RCAR).

As análises foram realizadas utilizando-se o programa Genes (Cruz, 1997).

Este trabalho foi desenvolvido seguindo as etapas descritas:

a) Avaliação das médias

Este estudo foi feito levando-se em consideração as médias dos cruzamentos;

b) Predição de híbridos triplos e duplos

Um dos processos de predição baseia-se na expressão geral, de acordo com Cruz & Regazzi (1994):

em que: M = média a ser predita; [X] = contribuição dos progenitores para o conjunto gênico que é repassado pelos gametas masculinos; [Y] = contribuição dos progenitores para o conjunto gênico que é repassado pelos gametas femininos.

c) Análise multivariada

A análise de agrupamento foi feita, adotando-se a distância euclidiana média como medida de dissimilaridade, sobre a qual empregou-se o método de agrupamento de otimização de Tocher.

c1) Medida de dissimilaridade

A distância euclidiana média foi utilizada para contornar a influência do número de caracteres e, para contornar o problema da escala dos dados, foi feita padronização dos dados da seguinte maneira:

em que: Xij é a observação no progenitor i (i = 1, 2,..., p), em referência ao caráter j (j = 1, 2,..., n); S(Xj) é o desvio-padrão dos dados do caráter j.

Então, a distância euclidiana média entre dois progenitores i e i' foi dada por:

em que: di i' é distância euclidiana média baseada em dados padronizados; n é o número de caracteres analisados.

c2) Método de otimização de Tocher

Um dos métodos de agrupamento mais indicados para estudar o grau de semelhança existente entre genótipos é o método de otimização de Tocher, que é utilizado na formação dos conjuntos de genótipos com maior semelhança entre si (Piassi, 1994).

Para o agrupamento segundo o método de Tocher, foi adotado o critério do maior entre os mínimos, ou seja, o acréscimo médio na distância intracluster (intragrupo), pela inclusão de um novo cruzamento a um dado cluster (grupo), deveria ser menor ou igual ao limite adotado. Estabeleceu-se que esse limite seria o valor máximo da medida de dissimilaridade encontrada no conjunto composto pelas menores distâncias de cada cruzamento a qualquer outro, na matriz de distâncias original. O acréscimo médio no valor da distância intragrupo foi calculado, dividindo-se a distância entre o cruzamento a ser inserido e o grupo que poderia recebê-lo pelo número de cruzamentos do grupo.

A matriz de dissimilaridade adotada foi a de distância euclidiana média, na qual se identificou, inicialmente, o par de cruzamentos que apresentava o menor valor de distância, para formar o primeiro grupo, quando essa distância não superava o limite estabelecido. Na seqüência, de acordo com o critério adotado, avaliou-se a possibilidade de inclusão de outros cruzamentos no grupo ou se haveria necessidade de formação de outros grupos, seguindo o mesmo critério.

Esquematicamente, mostra-se o critério para formação do agrupamento:

então, ou em que: a = limite de acréscimo adotado; j, k, l = cruzamentos do grupo; i = cruzamento a ser incluído, ou não, no grupo; n = número de cruzamentos do grupo que poderá receber i.

Resultados e Discussão

Na Tabela 1, estão apresentadas as médias das características em estudo, de acordo com os cruzamentos, incluindo-se, também, os valores do cruzamento testemunha.

A média é uma variável importante, pois, ao se selecionarem progenitores para compor programa de cruzamentos, aproveitando as informações providas pela capacidade especifica de combinação, inicialmente devem-se observar dentro dos progenitores os que apresentam maior média e depois os de maior capacidade geral de combinação.

Segundo Costa (1994), são padrões de produtividade de uma boa matriz:

· peso médio, às 66 semanas, da fêmea: 3.900g

· ovos produzidos/ave alojada: 184,5

· eclosão média: 85,55%

· fertilidade: 90%

De acordo com Torres (1996), a carcaça ideal nas aves deve ser bem musculosa, livre de defeitos e com alta produção de carne. Carcaças com maior produção de carne estão associadas com menores custos de processamento e de comercialização. Frangos com menos deposição de gordura têm mostrado pequena superioridade na produção de carcaça.

Ainda segundo o mesmo autor, boa produção de ovos com alta fertilidade e eclodibilidade é importante em todos os níveis de produção de frangos de corte. Nas linhas puras, estas características permitem intensidade de seleção ótima; em avós e matrizes, permitem menor custo de ovos e de pintos.

Pelos resultados apresentados, nota-se que as características de fertilidade e eclodibilidade estão dentro do esperado, excetuando-se um ou outro cruzamento que apresentou baixas fertilidade e eclodibilidade, assim como alguns valores de 100% para fertilidade, que não condizem com a realidade, mas que podem ter surgido por problemas na amostragem. Para as características do frango de corte, em geral, os resultados dos cruzamentos são melhores quando comparados com o cruzamento controle, mostrando a superioridade dos híbridos.

O híbrido de maior valor da característica P42 é o que envolve as linhas 5, do grupo I, e 4' do grupo II. Estas linhas também haviam sido selecionadas pela maior capacidade geral de combinação em trabalho anterior. Tanto para a característica P42, como para as características RCAR e FERT38, os 20 melhores híbridos envolvem cruzamentos que contêm as linhas selecionadas, em trabalho anterior pela capacidade geral de combinação, no grupo I e II.

Dos 648 híbridos triplos possíveis, ou seja pq (p + q - 2)/2 híbridos apresentam-se, neste trabalho, os 20 melhores (Tabelas 2 e 3), para as características FERT38, P42 e RCAR.

O híbrido de maior valor da característica P42 é o que envolve as linhas 5, do grupo I, e 4' do grupo II. Estas linhas também haviam sido selecionadas pela maior capacidade geral de combinação em trabalho anterior. Tanto para a característica P42, como para as características RCAR e FERT38, os 20 melhores híbridos envolvem cruzamentos que contêm as linhas selecionadas, em trabalho anterior pela capacidade geral de combinação, nos grupos I e II.

Dos 1.296 híbridos duplos possíveis neste estudo, ou seja, dos pq (p - 1)(q - 1)/4 híbridos, são apresentados somente os 40 melhores (Tabela 4).

Também para estes híbridos, os 40 melhores híbridos duplos são compostos por cruzamentos que envolvem as linhagens anteriormente selecionadas pela capacidade geral de combinação. Estes resultados mostram que a predição de híbridos oriunda de informações dos dialelos parciais é mais uma fonte importante de informação para o melhorista.

Na análise de agrupamento, estes agrupamentos foram feitos usando-se como medida de dissimilaridade a distância euclidiana média. O maior entre os mínimos estabelecido pelo limite de distância - Método de Tocher - foi de 0,912 (a = 0,912). Dessa forma, foram formados vinte e nove grupos de cruzamentos (Tabela 5), sendo que os cruzamentos foram distribuídos uniformemente dentro dos grupos, não existindo grupos com grandes concentrações de cruzamentos. Os grupos 2 e 1 foram os que apresentaram maior número de cruzamentos (7 e 6, respectivamente). O grupo 1 mostra a similaridade entre os cruzamentos obtidos pelo uso da linha 1 (linha macho) como progenitor, o que não é verificado para mais nenhum grupo.

Conclusões

As médias das características são fundamentais em programa de seleção para compor cruzamentos entre linhas, para melhoria de características tanto na matrizes, como no frango de corte.

A predição de híbridos é fonte importante de informação complementar para o melhorista.

Neste estudo, existe similaridade entre os cruzamentos onde se utiliza a linha 1 (linha macho) como progenitor, o que não acontece com as demais linhas.

Literatura Citada

Recebido em: 17/01/01

Aceito em: 08/11/01

2 Pesquisadora EMBRAPA/CNPSA, bolsista da CAPES. Embrapa Suínos e Aves, CP. 21. Concórdia ¾ SC. CEP. 89700-000. E.mail: valeria@cnpsa.embrapa.br

3 Professor da Universidade Federal de Minas Gerais.

4 Professor da Universidade Federal de Viçosa.

5 Pesquisador EMBRAPA/CNPSA, bolsista do CNPq. E.mail: elsio@cnpsa.embrapa.br

6 Pesquisador EMBRAPA/CNPSA. E.mail: pabreu@cnpsa.embrapa.br

  • ARUNACHALAM, V. Genetic distance in plant breeding. The Indian Journal of Genetic and Plants Breeding, v.41, n.2, p.226-236. 1981.
  • COSTA, C.A. Pontos críticos do manejo de matrizes. In: Manejo de matrizes Campinas: Fundação APINCO de Ciência e Tecnologia Avícolas, 1994. 10p.
  • CRUZ, C.D.; REGAZZI, A.J. Modelos biométricos aplicados ao melhoramento genético Viçosa, MG: Universidade Federal de Viçosa, 1994. 390p.
  • CRUZ, C.D. Aplicação de algumas técnicas multivariadas no melhoramento de plantas. Piracicaba: Escola Superior de Agricultura "Luiz de Queiroz", 1990. 188p. Tese (Doutorado em Genética e Melhoramento) - Escola Superior de Agricultura "Luiz de Queiroz", 1990.
  • CRUZ, C.D. Programa GENES - Aplicativo computacional em genética e estatística. Viçosa, MG: Universidade Federal de Viçosa, 1997. 442p.
  • FREITAS, R.T.F. Estudo da divergência genética de suínos em cruzamentos, utilizando-se técnica de análise multivariada. Viçosa, MG: Universidade Federal de Viçosa, 1996. 152p. Tese (Doutorado em Genética e Melhoramento) - Universidade Federal de Viçosa, 1996.
  • PIASSI, M. Avaliação do desempenho de linhagens de postura mantidas na Universidade Federal de Viçosa, em competição com marcas comerciais Viçosa, MG: Universidade Federal de Viçosa, 1994. 86p. Dissertação (Mestrado em Genética e Melhoramento) - Universidade Federal de Viçosa, 1994.
  • RAO, C.R. Advanced statistical methods in biometric research New York: John Wiley & Sons, 1952. 389p.
  • TORRES, R.A. Conseqüências do melhoramento genético. Cadernos Técnicos da Escola de Veterinária da Universidade Federal de Minas Gerais, n.17, p.73-77, 1996.
  • 1
    Parte do trabalho de tese de D.S. do primeiro autor, financiada pela CAPES/UFV/EMBRAPA.
  • Datas de Publicação

    • Publicação nesta coleção
      12 Ago 2002
    • Data do Fascículo
      2002

    Histórico

    • Aceito
      08 Nov 2001
    • Recebido
      17 Jan 2001
    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