Acessibilidade / Reportar erro

Métodos de agrupamento em estudo de divergência genética de pimentas

Clustering methods in a study of genetic diversity of peppers

Resumos

O objetivo deste trabalho foi comparar três métodos para determinação do número de grupos em estudos com aplicação de métodos hierárquicos de agrupamentos, baseando-se em dados obtidos a partir da caracterização de acessos de Capsicum, de modo a identificar aquele com maior poder de discriminação. Os métodos de Mojena, de Tocher e o método RMSSTD foram aplicados com a finalidade de determinar o número ideal de grupos formados na fase final do procedimento de agrupamento com o método UPGMA. Foram analisados 49 acessos da espécie Capsicum chinense do Banco de Germoplasma de Hortaliças da Universidade Federal de Viçosa, em relação a dez características morfológicas com o intuito de identificar e agrupar os acessos mais similares, tornando possível a seleção de genótipos superiores, ou seja, com as características comerciais de interesse. Os resultados mostraram que o método RMSSTD permitiu concluir sobre a existência de sete grupos, evidenciando um maior poder de discriminação para este método, em relação ao método de otimização de Tocher e ao método de Mojena, que formaram respectivamente, quatro e três grupos.

Capsicum chinense; método de Mojena; método de Tocher; RMSSTD


The objective of this study was to compare three methods to determine number of groups in studies with hierarquical cluster analysis, based at data from characterization of Capsicum accessions, in order to identify those with the greatest power of discrimination. Mojena method, Tocher method and RMSSTD method were applied to determine the optimal number of groups formed in final stage of the UPGMA method. Forty nine Capsicum chinense accessions from the Vegetable Germplasm Bank of the Universidade Federal de Viçosa were analyzed, in relation to ten morphological characters for identifying the most similar accessions group, making possible the selection of superior genotypes, i.e., of commercial interest. The RMSSTD method allowed to conclude on the existence of seven groups, with a greater power of discrimination for this method, compared to the Tocher optimization method and the Mojena method, which formed, respectively, four and three groups.

Capsicum chinense; Mojena method; Tocher method; RMSSTD


PESQUISA RESEARCH

Priscila N FariaI; Paulo R CeconII; Anderson R da SilvaII; Fernando L FingerII; Fabyano F e SilvaII; Cosme D CruzII; Filipe L SávioI

IUSP-ESALQ, 13418-900 Piracicaba-SP

IIUFV, 36570-000 Viçosa-MG; anderson.rodrigo@ufv.br

RESUMO

O objetivo deste trabalho foi comparar três métodos para determinação do número de grupos em estudos com aplicação de métodos hierárquicos de agrupamentos, baseando-se em dados obtidos a partir da caracterização de acessos de Capsicum, de modo a identificar aquele com maior poder de discriminação. Os métodos de Mojena, de Tocher e o método RMSSTD foram aplicados com a finalidade de determinar o número ideal de grupos formados na fase final do procedimento de agrupamento com o método UPGMA. Foram analisados 49 acessos da espécie Capsicum chinense do Banco de Germoplasma de Hortaliças da Universidade Federal de Viçosa, em relação a dez características morfológicas com o intuito de identificar e agrupar os acessos mais similares, tornando possível a seleção de genótipos superiores, ou seja, com as características comerciais de interesse. Os resultados mostraram que o método RMSSTD permitiu concluir sobre a existência de sete grupos, evidenciando um maior poder de discriminação para este método, em relação ao método de otimização de Tocher e ao método de Mojena, que formaram respectivamente, quatro e três grupos.

Palavras-chave:Capsicum chinense, método de Mojena, método de Tocher, RMSSTD.

ABSTRACT

The objective of this study was to compare three methods to determine number of groups in studies with hierarquical cluster analysis, based at data from characterization of Capsicum accessions, in order to identify those with the greatest power of discrimination. Mojena method, Tocher method and RMSSTD method were applied to determine the optimal number of groups formed in final stage of the UPGMA method. Forty nine Capsicum chinense accessions from the Vegetable Germplasm Bank of the Universidade Federal de Viçosa were analyzed, in relation to ten morphological characters for identifying the most similar accessions group, making possible the selection of superior genotypes, i.e., of commercial interest. The RMSSTD method allowed to conclude on the existence of seven groups, with a greater power of discrimination for this method, compared to the Tocher optimization method and the Mojena method, which formed, respectively, four and three groups.

Keywords:Capsicum chinense, Mojena method, Tocher method, RMSSTD.

A diversidade genética pode ser avaliada de forma simultânea em relação às várias características, e para isso recomenda-se a utilização de medidas de dissimilaridade (Cruz & Carneiro, 2003). Uma forma prática e eficiente de se obter essas medidas é por meio da análise de agrupamentos (ou análise de cluster), a qual tem por finalidade reunir os indivíduos em grupos, de forma que exista a máxima homogeneidade dentro do grupo e a máxima heterogeneidade entre os grupos (Johnson & Wichern, 1992; Cruz & Regazzi, 2001). Dos métodos de agrupamento, os mais utilizados são os de otimização e os hierárquicos (Cruz et al., 2011).

Uma dificuldade na etapa final dos estudos que se utilizam de algoritmos hierárquicos de agrupamento é a falta de critérios objetivos para identificar o número ideal de grupos formados, uma vez que na prática este número é dado simplesmente por uma inspeção gráfica visual ou estabelecido em pontos de alta mudança de nível dos dendrogramas (Milligan, 1981). Milligan & Cooper (1985) apresentaram um estudo comparativo de 30 critérios de corte para determinação do número de agrupamentos, e utilizando dados artificiais com número conhecido de agrupamentos, mostraram que critérios diferentes podem conduzir a resultados muito diferentes.

Um dos métodos considerados como "objetivos", dentre os poucos existentes, é o de Mojena (1977). Este método recorre a um critério estatístico e propõe um procedimento de cálculo baseado no tamanho relativo dos níveis de fusões ou distâncias no dendrograma. Dessa forma, o método não necessita do conhecimento prévio da conformação dos grupos, ao contrário do método subjetivo, com base na inspeção visual das ramificações do dendrograma.

Vários trabalhos têm utilizado a análise de agrupamentos no estudo da diversidade genética e relatado a subjetividade no estabelecimento do ponto de corte no dendrograma. Dentre eles, Sudré et al. (2005) analisaram a divergência genética entre acessos de pimenta e pimentão utilizando técnicas multivariadas, obtendo a separação dos acessos em sete grupos distintos, utilizando pontos de mudança brusca de nível como referência para estabelecer o corte. Karasawa et al. (2005) relataram que a utilização do procedimento subjetivo, baseado no exame visual do dendrograma, pode gerar alguma dificuldade na tomada de decisão quanto ao número de grupos gerados, uma vez que qualquer inferência rígida sobre este número pode não ser produtiva.

Na tentativa de solucionar esse problema, Faria (2009) propôs que o Método da Máxima Curvatura Modificado (Lessman & Atkins, 1963) fosse utilizado juntamente com o índice RMSSTD (Root Mean Square Standard Deviation). A metodologia proposta é fundamentada na obtenção do ponto de máxima curvatura para a trajetória do índice RMSSTD dada em função do número de grupos, onde o referido ponto é que determina o número ideal de grupos.

O objetivo deste trabalho foi comparar três métodos aplicados com a finalidade de determinar o número de grupos de agrupamentos hierarquicos, quais sejam: Método de Mojena, Método de otimização de Tocher e o método que utiliza o índice RMSSTD, a fim de se detectar aquele com maior poder em discriminar grupos homogêneos em relação a caracteres do fruto de pimentas Capsicum chinense.

MATERIAL E MÉTODOS

O Banco de Germoplasma de Hortaliças da Universidade Federal de Viçosa (BGH/UFV) possui cerca de 100 acessos de C. chinense com grande potencial para contribuir em programas de melhoramento, uma vez que grandes quantidades de informações a respeito de características agronômicas, químicas, morfológicas e moleculares já foram coletadas. Neste trabalho foram utilizados resultados de caracterização feita para 49 desses acessos por Faria (2009).

O experimento foi conduzido em área experimental da UFV no ano de 2009 sob delineamento em blocos casualizados com três repetições, utilizando espaçamento de 1x1 m, sendo cada linha constituída por três plantas de cada acesso. Os tratos culturais seguiram as recomendações feitas por Filgueira (2000). Foram utilizados dados de caracteres quantitativos de fruto, os quais foram colhidos no estádio maduro e imediatamente caracterizados de acordo com os descritores de Capsicum (IPGRI, 1995). Dez características foram avaliadas: comprimento do fruto (mm); maior diâmetro do fruto (mm); espessura do pericarpo (mm); porcentagem de matéria seca (g/100g); massa seca do fruto (g); massa fresca do fruto (g); capsaicina total (mg/gMS); teor de sólidos solúveis (ºBrix); Vitamina C (mg/100 g de fruto fresco); cor extraível (unidades ASTA de cor). A identificação dos referidos acessos pode ser verificada em Faria (2009).

Realizou-se análise de variância univariada para se verificar a existência de variabilidade genética significativa entre os 49 acessos, em relação às características avaliadas com base na média das parcelas.

Foram calculadas as médias de cada acesso para cada variável analisada. A matriz de correlações residuais entre as variáveis foi estimada por ocasião da realização de análise de variância univariada. Foi aplicado o teste de esfericidade de Bartlet (Rencher, 2002) para testar a hipótese de que esta é uma matriz diagonal, isto é, se as dez variáveis apresentadas são independentes. Tendo encontrado resultado não significativo (p>0,05) pela aproximação de qui-quadrado com 45 graus de liberdade, a distância generalizada de Mahalanobis coincidiu com a distância euclidiana como medida de dissimilaridade. Pode-se concluir então, que a utilização de quaisquer das distâncias supracitadas na análise de agrupamento realizada neste trabalho não altera significativamente os resultados obtidos. Desta forma, optou-se pelo uso da distância Euclidiana. A diversidade genética entre os acessos foi avaliada utilizando-se o método da ligação média entre grupos ou UPGMA (Unweighted Pair Group Method with Arithmetic Mean).

Três métodos foram utilizados para determinação do número ótimo de grupos no dendrograma, a saber: o Método de Mojena, o Método de Tocher e o método que utiliza o índice RMSSTD.

O método de Mojena (1977) é um procedimento baseado no tamanho relativo dos níveis de fusões (distâncias) no dendrograma. Este método consiste em selecionar o número de grupos no estágio j que, primeiramente, satisfizer à seguinte inequação: αj > θk em que αj é o valor de distâncias dos níveis de fusão correspondentes ao estágio j (j = 1, 2, ..., n) e θk é o valor referencial de corte, expresso por: , em que α e são, respectivamente, as estimativas não viesadas da média e do desvio padrão dos valores de α; k é uma constante. Adotou-se k= 1,25 como regra de parada na definição do número de grupos, como sugerem Milligan & Cooper (1985).

O método de otimização de Tocher é um método de agrupamento que se baseia na formação de grupos cujas distâncias dentro dos grupos sejam menores que as distâncias entre grupos. Ao final do processo obtém-se o número de grupos e os acessos contidos em cada grupo. A aplicação deste método foi feita conforme sugerem Cruz & Carneiro (2003).

A determinação do número de grupos pelo índice RMSSTD foi realizada por meio de análise gráfica (determinação geométrica) dos valores da raiz quadrada do desvio padrão médio (valores de RMSSTD). De acordo com Sharma (1996) o valor do índice deve ser tão baixo quanto possível para um agrupamento, isto é, quanto menor o RMSSTD, mais homogêneo ou compacto é o agrupamento formado a um determinado passo.

Juntamente com o índice RMSSTD, utilizou-se o Método da Máxima Curvatura Modificado (Lessman & Atkins, 1963), que se baseia na modificação do método de Smith (1938). Relacionaram-se então os valores de RMSSTD em função do número de grupos. Assim, o ponto que estabelece declínio acentuado do RMSSTD indica no eixo das abscissas o número ótimo de grupos. A metodologia pode ser expressa pela seguinte equação,

em que X representa o número de clusters (grupos), a e b são constantes estimadas. A partir desta função de curvatura, determinou-se o valor da abscissa onde ocorre o ponto de máxima curvatura ou ponto crítico da função de curvatura, conforme apresentado por Meier & Lessman (1971), por meio da seguinte equação,

em que XMC é o número de grupos que representa a máxima curvatura em RMSSTD; a e b são constantes estimadas.

As análises foram realizadas utilizando-se o módulo 'Diversidade Genética' do programa GENES (Cruz, 2006), para obtenção do número de grupos pelo Método de Mojena e pelo Método de Tocher; e SAS® (SAS, 2003), para o método de agrupamento UPGMA e para empregar a metodologia proposta por Faria (2009), que utiliza o método da máxima curvatura modificado e o índice RMSSTD.

RESULTADOS E DISCUSSÃO

Pela análise de variância (p<0,01), foi possível detectar diferenças significativas entre os acessos para todas as características avaliadas, indicando a presença de variabilidade genética. Os coeficientes de variação CV (%) do experimento variaram de 5,34% (capsaicina total) a 36,82% (cor extraível).

De acordo com a Tabela 1 é possível verificar que os acessos 5 e 28 foram os mais similares geneticamente, possuindo a menor distância (9,47); entre os acessos 1 e 27 houve a maior magnitude (463,59), sendo portanto os acessos mais dissimilares. Observa-se ainda que é viável um corte no dendrograma na altura de θ = 135,27 (critério geral: ), indicando que o número ideal de grupos deve ser igual a três. Sendo assim, o primeiro grupo foi composto somente pelo acesso 27, o segundo grupo pelos acessos 2, 12, 14, 45, 13, 29, 37 e o terceiro grupo contendo o restante dos acessos (Figura 1).


Pelo método de Tocher, foi detectada a formação de quatro grupos distintos, sendo os grupos II, III e IV formados por apenas um acesso cada (Tabela 2). O grupo I foi composto por 46 acessos, aproximadamente 94%. Verifica-se, portanto, a dificuldade em analisar a divergência entre os acessos, visto que a maioria deles encontra-se em apenas um grupo. Em estudo sobre dissimilaridade genética entre 17 acessos de pimentas Capsicum spp., Neitzke et al. (2010) encontraram resultados semelhantes, quando verificaram a formação de quatro grupos distintos pelo método de Tocher para dados quantitativos. Os mesmos autores relatam ainda que o primeiro grupo também reuniu a grande maioria dos acessos (13 acessos, cerca de 76%), e que, quando o método foi aplicado a dados qualitativos, foi maior o número de grupos formados (oito grupos distintos).

A informação do índice RMSSTD permitiu a identificação do número ótimo de grupos por meio da máxima curvatura (Figura 2). De posse desta informação, os 49 acessos foram agrupados, e a utilização do índice RMSSTD, juntamente com o método da máxima curvatura modificado possibilitou a formação de sete grupos (Tabela 3). Posteriormente foi possível identificar um ponto de corte no dendrograma, obtido pela utilização do método UPGMA com base na distância Euclidiana (Figura 2). Observou-se que o maior grupo formado foi composto por 24 acessos e os grupos IV, V e VII formados apenas por um acesso cada: o acesso 45 pertencente ao grupo IV, o acesso 13 pertencente ao grupo V e o acesso 27 pertencente ao grupo VII (Tabela 3). O acesso 27 mostrou-se bastante divergente dos demais, visto que formou um grupo exclusivo e permaneceu isolado dos demais acessos no dendrograma. Sua divergência genética pode ser explorada em programas de melhoramento com base em descritores quantitativos para obtenção de novas cultivares.


Utilizando os métodos de agrupamento de Tocher e UPGMA em estudo sobre diversidade genética entre 23 acessos de espécies cultivada de pimentas Capsicum spp., Monteiro et al. (2010) relataram que a espécie C. chinense foi a mais divergente em todos os agrupamentos, com as maiores distâncias intragrupo. Este fato está de acordo com os resultados obtidos no presente trabalho, evidenciando o elevado número de acessos num mesmo grupo. Nesse caso, tendo a espécie C. chinense apresentado grande divergência genética, é esperado que haja a formação de mais grupos e até mesmo com alguns dos grupos contendo muitos acessos, uma vez que apenas poucos destes acessos apresentam grande dissimilaridade. Assim, um método mais sensível para determinação do número de grupos consegue captar mais eficientemente a existência de acessos discrepantes.

A aplicação da metodologia do RMSSTD aos dados de pimenta permitiu concluir sobre a existência de sete grupos. Este número de grupos foi superior ao obtido pelo método de Tocher (quatro grupos) e ao método de Mojena (três grupos), sendo, portanto de maior poder de discriminação, possibilitando a identificação de mais grupos contendo acessos similares.

Recebido para publicação em 20 de abril de 2011

Aceito em 25 de junho de 2012

  • CRUZ CD. 2006. Programa Genes: Biometria. Viçosa: Editora UFV. 382p.
  • CRUZ CD; CARNEIRO PCS. 2003. Modelos biométricos aplicados ao melhoramento genético Volume 2. Viçosa: UFV. 585p.
  • CRUZ CD; FERREIRA FM; PESSONI LA. Biometria aplicada ao estudo da diversidade genética Visconde do Rio Branco: Suprema, 2011. 620p.
  • CRUZ CD; REGAZZI AJ. 2001. Modelos biométricos aplicados ao melhoramento genético. 2 ed. Viçosa: UFV. 390p.
  • FARIA PN. 2009. Avaliação de métodos para determinação do número ótimo de clusters em estudo de divergência genética entre acessos de pimenta Viçosa: UFV. 67p. (Dissertação mestrado).
  • FILGUEIRA FAR. 2000. Novo Manual de Olericultura Viçosa: UFV. 402p.
  • IPGRI AVRDC; CATIE. 1995. Descriptors for Capsicum (Capsicum spp.). International Plant Genetic Resources Institute, Rome, Italy; the Asian Vegetable Research and Development Center, Taipei, Taiwan, and the Centro Agronómico Tropical de Investigación y Enseñanza, Turrialba, Costa Rica. 110p.
  • JOHNSON RA; WICHERN DW. 1992. Applied Multivariate Statistical Analysis. New Jersey: Englewood Cliffs. 642p.
  • KARASAWA M; RODRIGUES R; SUDRÉ CP; SILVA MP; RIVA EM; AMARAL JÚNIOR AT. 2005. Aplicação de métodos de agrupamento na quantificação da divergência genética entre acessos de tomateiro. Horticultura Brasileira 23: 1000-1005.
  • LESSMAN KJ; ATKINS RE. 1963. Optimum plot size and relative efficiency of lattice designs for grain sorghum yield tests. Crop Science 3: 477-481.
  • MEIER VD; LESSMAN RJ. 1971. Estimation of optimum field plot shape and size for testing yield in Crambe abyssinica Hochst. Crop Science 11: 648-645.
  • MILLIGAN GW. 1981. A Monte Carlo study of thirty internal criterion measures for cluster analysis. Psychometrika 46: 187-199.
  • MILLIGAN GW; COOPER MC. 1985. An examination of procedures for determining the number of clusters in a data set. Psychometrika 50: 159-179.
  • MOJENA R. 1977. Hierarchical grouping methods and stopping rules: an evaluation. The Computer Journal 20: 359-363.
  • MONTEIRO ER; BASTOS EM; LOPES ACA; GOMES RLF; NUNES JAR. 2010. Diversidade genética entre acessos de espécies cultivadas de pimentas. Ciência Rural 40: 288-283.
  • NEITZKE RS; BARBIERI RL; RODRIGUES WF; CORRÊA IV; CARVALHO FIF. 2010. Dissimilaridade genética entre acessos de pimenta com potencial ornamental. Horticultura Brasileira 28: 47-53.
  • RENCHER AC. 2002. Methods of multivariate analysis New York: John Wiley. 708p.
  • SUDRÉ CP; RODRIGUES R; RIVA EM; KARASAWA M; AMARAL JÚNIOR AT. 2005. Divergência genética entre acessos de pimenta e pimentão utilizando técnicas multivariadas. Horticultura Brasileira 23: 22-27.
  • SAS INSTITUTE. 2003. SAS System: SAS/STAT version 9.1 (software). Cary.
  • SHARMA S. 1996. Applied multivariate techniquies. NewYork: John Wiley & Sons. 493p.
  • SMITH HF. 1938. An empirical law describing heterogeneity in the yields of agricultural crops. Journal of Agricultural Science 28: 1-23.
  • Métodos de agrupamento em estudo de divergência genética de pimentas

    Clustering methods in a study of genetic diversity of peppers
  • Datas de Publicação

    • Publicação nesta coleção
      04 Out 2012
    • Data do Fascículo
      Set 2012

    Histórico

    • Recebido
      20 Abr 2011
    • Aceito
      25 Jun 2012
    Associação Brasileira de Horticultura Embrapa Hortaliças, C. Postal 218, 70275-970 Brasília-DF, Tel. (61) 3385 9099, Tel. (81) 3320 6064, www.abhorticultura.com.br - Vitoria da Conquista - BA - Brazil
    E-mail: associacaohorticultura@gmail.com