Para o estudo da diversidade genética com base em marcadores RAPD, utilizaram-se 43 acessos, sendo 42 do gênero Abelmoschus e um do gênero Hibiscus. Dos 42 acessos de Abelmoschus spp, quinze são A. esculentus, quinze A. caillei, dois A. manihot, oito A. moschatus, um A. ficulneus e um A. tetraphyllus. Estes genótipos foram fornecidos pela Estação Regional da Geórgia de Introdução de Plantas do Departamento de Agricultura dos Estados Unidos (USDA), os quais foram solicitados através do Centro Nacional de Recursos Genéticos (CENARGEN/EMBRAPA) e atualmente estão sendo mantidos pelo Banco de Germoplasma da UENF. Os genótipos constam de espécies cultivadas e silvestres, de diversas procedências, conforme a Tabela 1.

As plantas foram cultivadas em vasos, em condição de casa de vegetação na Unidade de Apoio à Pesquisa da UENF. As amostras para análise dos marcadores de DNA foram constituídas de discos foliares obtidos de cinco plantas de cada acesso, que foram embaladas, identificadas e mergulhadas em nitrogênio líquido. A extração de moléculas a partir de folhas de quiabeiro foi realizada no Laboratório de Melhoramento Genético Vegetal do CCTA, com base no protocolo desenvolvido por Doyle & Doyle (1990).

Para o isolamento do DNA, foram utilizados cerca de 500 mg da mistura de tecidos macerados obtidos a partir de discos foliares das cinco plantas. Ao tecido foi adicionado 700 ml de CTAB com forte agitação em vortex por 30 segundos. A mistura foi incubada a 60ºC por 1 hora e centrifugada para a retirada de impurezas brutas. O sobrenadante passou por duas extrações com clorofórmio-álcool isoamílico (24:1), ambas seguidas por centrifugação (14.000 rpm, por 5 minutos em centrífuga Eppendorf modelo 5415 C). Para a precipitação do DNA, foi adicionado 2/3 do volume da amostra de isopropanol seguido por uma lavagem em etanol a 80%. Após a secagem do pellet, o DNA foi ressuspenso em tampão TE (0,01M Tris-HCl; 0,001M EDTA) pH 8,00. Em seguida foi adicionado RNAse (10 mg ml^-1) mantendo-se incubado por 30 min a 37ºC. Passou por uma nova lavagem com etanol (80%) e, após, nova ressuspensão em TE. A concentração do DNA foi quantificada pelo fluorômetro Hoefer modelo DyNA Quant 200.

Para verificar a integridade do DNA, géis de agarose a 1,4% foram preparados para eletroforese e, a cada poço do gel foi aplicada uma amostra de 3 ml de DNA. Após corrida a 100 V por 2 horas, os géis foram, então, banhados em brometo de etídio (0,5 mg/litro) por 30 a 45 min à temperatura ambiente. A descoloração foi realizada imergindo o gel em água durante 20 min também a temperatura ambiente e o DNA visualizado pelo fotodocumentador Stratagene modelo Eagle Eye II.

Um acesso de cada espécie foi testado para 60 iniciadores para avaliar o polimorfismo de DNA, com o propósito de selecionar aqueles que amplificassem simultaneamente para todas as espécies. O procedimento de amplificação foi similar ao descrito por Williams et al. (1990). O procedimento de PCR («Polymerase Chain Reaction») foi realizado com 30 ng de DNA, 0,1 mM de cada dNTP (Perkin-Elmer, USA), 2,4 mM de MgCl₂ (Perkin-Elmer, USA), 1U (unidade) de Taq DNA-Polimerase em tampão de PCR (Pharmacia Biotech, USA), 10mM/50mM de Tris-KCl e 0,4 mM de iniciador (decâmero – Operon Technologies). As concentrações de cada componente foram ajustadas para o volume final de 25 ml de reação com H₂O ultra pura estéril. As condições da reação foram as seguintes: 1 min a 90ºC, seguida por 55 ciclos, cada ciclo constituindo as etapas de desnaturação (1 min, 94ºC), anelamento do iniciador (1 min, 36ºC) e extensão do iniciador (2 min, 72ºC). Após os 55 ciclos, foi inserido um ciclo de 7 min a 72ºC para extensão final. As amplificações foram desenvolvidas em termociclador GeneAmp 9600 (Perkin-Elmer, USA).

Após a amplificação, o volume final da reação foi de 25 ml. Entretanto, adicionou-se 5 ml de corante «Blue Juice» (6X) e o volume final passou para 30 ml. O DNA amplificado foi aplicado em gel de agarose a 1,4% em cuba de eletroforese horizontal e submetido a 100 V por 2 horas. Após a eletroforese, o gel foi banhado em brometo de etídio (0,5 ml/litro) por 30 a 45 min a temperatura ambiente. A descoloração foi realizada imergindo o gel em água durante 20 min também a temperatura ambiente e o DNA visualizado sob luz ultravioleta, para, a seguir, ser fotodocumentado e selecionados os iniciadores polimórficos a serem utilizados nos ensaios de RAPD.

Após a seleção dos iniciadores, procedeu-se o ensaio de polimorfismo de DNA amplificado ao acaso (RAPD), utilizando-se a técnica de PCR descrita anteriormente (Williams et al., 1990). A leitura das bandas no gel foi feita pela tabulação dos resultados em 0 ou 1, respectivamente para ausência ou presença de banda e 2 para os dados perdidos ou dúbios. A tabela de dados binários resultante foi utilizada para análise da divergência genética dos 43 acessos de quiabeiro.

Na interpretação dos resultados, gerou-se uma matriz de dissimilaridade genética a partir do complemento aritmético do Índice de Jaccard (Sneath & Sokal, 1973), representado por:

a = número de casos em que a banda está presente nos dois genótipos, simultaneamente; b = número de casos em que a banda está presente somente no genótipo i; c = número de casos em que a banda está presente somente no genótipo j.

Posteriormente, adotou-se o método hierárquico do «vizinho mais próximo» para o agrupamento dos acessos. No procedimento analítico, partindo-se da matriz de distâncias genéticas, procedeu-se a sucessivas identificações dos genótipos mais próximos, a partir do par mais semelhante, até que se estabeleceu um diagrama de árvore.

A distância entre um acesso k e um grupo formado pelos acessos i e j, é dada por: D²_(ij)k = min (D²_ik; D²_jk),

em que:

D²_(ij)k é a distância entre o grupo ij e o acesso k, e min (D²_ik; D²_jk) é a menor distância entre os grupos de acessos ik e jk.

O agrupamento dos acessos foi feito, também, pela utilização do método de Tocher (Rao, 1952), da forma seguinte: a partir da matriz de distâncias genéticas identificou-se o par mais próximo que constituiu o primeiro grupo; concomitantemente, obteve-se q, que é o maior valor do conjunto de menores distâncias entre os acessos. A partir de então, avaliou-se a possibilidade de inclusão de outros acessos neste grupo, respeitando-se o critério da distância média intragrupo inferior a quaisquer distâncias intergrupos. Para tanto, procedeu-se à comparação entre o acréscimo no valor médio da distância dentro do grupo e o máximo permitido, comumente denominado q. A inclusão de um acesso em um grupo foi possível quando:

D²_(ij)k = D²_ik + D²_jk, em que:

D²_ik = distância entre os acessos i e k; D²_jk = distância entre os acessos j e k; e n = número de acessos que constitui o grupo original.

A partir da análise de agrupamento pelo método de Tocher, os genótipos pertencentes ao primeiro grupo foram reagrupados pelo método hierárquico do vizinho mais próximo, com o propósito de se verificar o novo perfil da divergência genética.

As análises genético-estatísticas foram realizadas com o auxílio do Programa Genes (Cruz, 1997).

Foram utilizados 31 iniciadores, previamente selecionados por exibirem pelo menos uma banda polimórfica no gel para os 43 genótipos de quiabeiro. Nenhum destes oligonucleotídios mostrou monomorfismo, provavelmente por estarem envolvidas diversas espécies do gênero Abelmoschus e também um acesso do gênero Hibiscus. Os 31 iniciadores amplificaram um total de 103 fragmentos polimórficos, na proporção de 3,3 RAPDs por iniciador. Até 6 fragmentos foram polimórficos para diversos iniciadores.

Os acessos de número 1293, 1294, 1295, 1300, 1305, 1322, 1323 e 1324, identificados na Tabela 1 como A. manihot, mostraram-se tanto fenotipica como genotipicamente muito similares a A. caillei, portanto serão considerados como pertencentes à última espécie, o que é fundamentado pelo fato de que estes acessos são provenientes do oeste africano, onde são, muitas vezes, tratados como A. manihot, embora a descrição correta seja A. caillei. Existe, ainda, a possibilidade de ocorrência de erro na descrição dos acessos no Banco de Germoplasma, pois A. caillei é considerada uma espécie distinta e não mais uma variedade de A. manihot.

O agrupamento dos 43 genótipos com base no método hierárquico do vizinho mais próximo, utilizando marcadores RAPD, separou os acessos, de modo geral, conforme as espécies botânicas, formando seis grupos distintos (Figura 1). O primeiro grupo foi composto pelos acessos de A. esculentus e A. caillei, por se mostrarem muito semelhantes. Os acessos com maior proximidade genética foram 4 e o 13 (A. esculentus var. Itaocara I e A. esculentus de Camarões), seguidos por 2 e 10 (A. esculentus, cv. Santa Cruz e A. esculentus da Algéria). Por outro lado, os acessos de A. esculentus que mostraram-se mais distantes foram os de número 6 e 15 (A. esculentus da Síria e A. esculentus cv. Parbhani Kranti). Com relação aos acessos de A. caillei, os genótipos 16 com 29 e 19 com 32 formaram, respectivamente, os dois pares mais semelhantes.

O segundo grupo é composto pelos acessos de A. moschatus, o terceiro pelos acessos de A. manihot, o quarto envolveu a espécie A. tetraphyllus e o quinto, A. ficulneus. A espécie H. sabdariffa formou, separadamenrte o grupo 6, mostrando ser a mais distante dentre os 43 acessos avaliados.

A projeção das distâncias genéticas no plano bidimensional (Figura 2), foi eficaz na separação dos acessos. Por esta figura, pode-se observar que os acessos de A. esculentus e de A. caillei, formaram um único grupo. Este resultado está de acordo com o agrupamento pelo método hierárquico (Figura 1), quando procede-se um corte vertical a uma distância de ligação de aproximadamente 0,35. Assim, as espécies A. esculentus e A. caillei revelaram semelhanças genéticas em elevada magnitude, constituindo genomas muito aproximados.

Martinello et al. (1996), numa análise por componentes principais aplicada a características morfológicas em quiabeiro, verificaram que a espécie A. caillei foi a única das silvestres que fez parte do grupo em que todos os outros acessos eram A. esculentus. Estes autores verificaram, ainda, que hibridações interespecíficas entre A. caillei (2n= 185-189) e A. esculentus (2n= 72-144) foram eficazes e que, os híbridos F₁ são férteis e produzem grandes quantidades de sementes viáveis, apesar da diferença do número de cromossomos entre os parentais. Isto indica uma grande afinidade entre os genomas, que podem fazer parte de um mesmo «pool» gênico. Estas constatações denotam que, no processo de evolução, estas espécies divergiram ou bifurcaram mais recentemente. Talvez A. caillei se diferenciou de A. esculentus, ou vice-versa, por um processo evolucionário reticulado, isto é, anfiploidização, poliploidização, hibridação ou outro sistema genético (Baum & Gupta, 1996). Convém salientar que este tipo de situação envolvendo evolução reticulada para as plantas superiores, ainda é muito pouco estudado. Siemonsma (1982) hipotizou que A. caillei é um anfiplóide natural de A. esculentus (2n= 130-140) e A. manihot (2n= 60-68).

Ainda com relação à análise gráfica das distâncias genéticas no plano bidimensional (Figura 2), as espécies representadas por apenas um acesso formaram grupos individuais, como é o caso de A. tetraphyllus , A. ficulneus e Hibiscus sabdariffa. Os acessos de A. moschatus e A. manihot também formaram grupos separados. Sendo assim, os marcadores RAPD foram capazes de classificar e também diferenciar corretamente as espécies, com exceção do genótipo 31 (A. caillei) que se assemelhou aos acessos de A. esculentus.

A aplicação do método de agrupamento por otimização de Tocher, utilizando a matriz de distâncias obtidas pelo complemento aritmético do Índice de Jaccard, permitiu a formação de quatro grupos distintos (Tabela 2). O grupo 1 é composto pelos acessos de A. caillei, A. esculentus, A. moschatus e A. manihot, sendo que os acessos das espécies A. caillei e A. esculentus mostraram-se altamente similares e estão representados pelos genótipos de número 4, 13, 5, 31, 9, 3, 1, 2, 12, 10, 11, 14, 7, 8, 15, 6, 29, 16, 18, 19, 26, 21, 17, 32, 25, 30, 28, 22, 24 e 20. Neste grupo, os acessos de número 4 e 13 são os mais próximos. Os genótipos de números 40, 39, 36, 33, 34, 35 e 38 pertencem a espécie A. moschatus e os de número 27 e 23, à espécie A. manihot que, por sua vez, mostraram-se mais distantes em relação às espécies A. caillei e A. esculentus. O grupo 1 compreendeu aproximadamente 90% do germoplasma avaliado, indicando baixa amplitude da variabilidade genética nos genótipos em estudo.

A análise de agrupamento por otimização de Tocher, mostrou-se coerente apenas em parte à análise de agrupamento hierárquico, vez que o grupo 1 conteve, além das espécies A. esculentus e A. caillei, A. moschatus e A. manihot. O grupo 2 compreendeu A. tetraphyllus e A. ficulneus. Por sua vez, A. moschatus compôs o grupo 3 e H. sabdariffa o grupo 4.

Quando aplicou-se o método hierárquico do vizinho mais próximo apenas sobre o grupo I de Tocher (Tabela 2), houve, consequentemente, o surgimento de seis fragmentos monomórficos, resultando 88 locos polimórficos. Esse polimorfismo foi reanalisado pelo complemento aritmético do Índice de Jaccard. A redistribuição hierárquica dos 39 genótipos que compunham o grupo I de Tocher (Figura 3), mostrou-se bastante semelhante ao apresentado pela Figura 1, sendo que na Figura 3 houve a formação de quatro grupos, quando passa-se uma linha vertical a uma distância de 0,35. O primeiro grupo reuniu mais uma vez os acessos de A. esculentus e A. caillei. Os acessos 2 e 10 mostraram-se, desta vez, os mais próximos, porém o acesso 13 mostrou-se mais distante com relação aos demais do mesmo grupo, inclusive afastando-se bastante do acesso 4, que na Figura 1 apresentou-se intimamente associado àquele. O grupo 2, da Figura 3, permaneceu idêntico ao grupo 2 da Figura 1. O grupo 3 reuniu os acessos de A. manihot que, por sua vez, mostraram-se os mais distantes dentre os 39 avaliados. Assim, o reagrupamento dos acessos revelou a existência de uma heterogeneidade genética no grupo I que, desta forma, tornou os resultados da metodologia de Tocher semelhantes aos obtidos com o método hierárquico do vizinho mais próximo.

O reagrupamento tem sido utilizado por alguns autores, como forma alternativa de tornar os resultados mais elucidativos. Amaral Júnior (1996), por exemplo, trabalhando com tomateiro, subdividiu o grupo I de Tocher, em dois subgrupos para uma melhor comparação entre a diversidade obtida com marcadores RAPD e características morfoagronômicas. Vasconcellos (1995), por sua vez, trabalhando com feijoeiro, também procedeu à redistribuição intragrupo para uma melhor visualização das distâncias entre os indivíduos que foram separados com base nos padrões de faseolina e nos tamanhos de sementes.

Com esta pesquisa verificou-se que os marcadores RAPD foram eficientes no estudo da diversidade genética do quiabeiro. Verificou-se ainda que as técnicas estatísticas empregadas foram capazes de estabecer padrões e classificações da diversidade genética, revelando a formação de seis grupos. O agrupamento pelo método hierárquico do vizinho mais próximo, com base nos marcadores RAPD, revelou que os acessos de A. esculentus, apesar das diferentes procedências, são bastante homogêneos e que estão muito próximos de A. caillei. O reagrupamento dos acessos, por sua vez, revelou a existência de uma heterogeneidade genética no grupo I que, desta forma, tornou os resultados da metodologia de Tocher semelhantes aos obtidos com o método hierárquico do vizinho mais próximo. Este trabalho apresentou informações valiosas que podem ser utilizadas na escolha de progenitores mais adequados para iniciar um programa de melhoramento da cultura em questão.

Agradeço à Fundação Estadual Norte Fluminense (FENORTE) pela concessão da bolsa de estudos.

AMARAL Jr, A.T. Análise dialélica de betacaroteno, vitamina C, sólidos solúveis e produção e variabilidade em cultivares de tomateiro (Lycopersicon esculentum Mill.) via marcadores RAPD. Viçosa: UFV. Imprensa Universitária, 1996, 198 p. (Tese doutorado).

BAUM, B.R.; GUPTA, P.K. Numerical taxonomy and genome analysis. In: JAUHAR, P.P. (ed.) Methods of genome analysis in plants. North Dakota: CRC Press, 1996. p. 345-356.

BRETTING, P.K.; WIDRLECHENER, M.P. Genetic markers and horticultural germplasm management. HortScience, v. 30, n. 7, p. 1349-1355, 1995.

CRUZ, C.D. Programa Genes: aplicativo computacional em genética e estatística. Viçosa, MG: Editora UFV, 1997. 442 p.

DOYLE, J.J.; DOYLE, J.L. Isolation of plant DNA from fresh tissue. Focus, v. 12, p. 13-15, 1987.

FERREIRA. M.E.; GRATTAPAGLIA, D. Introdução ao uso de marcadores RAPD e RFLP em análise genética. Brasília: MAARA/EMBRAPA/CENARGEN, 1995. 220 p.

HAMON, S.; CHARRIER, A.; KOEKLIN, J.; VAN SLOTEN, D.H. Potential contributions to okra breeding through the study of their genetic resources. In: REPORT OF AN INTERNATIONAL WORKSHOP ON OKRA GENETIC RESOURCES. 1990, New Delhi. Anais…; Roma: IBPGR, 1991. p. 77-88.

MARTINELLO, G.E.; ABBOUD, A.C.S.; LEAL, N.R. Análise em componentes principais e de agrupamentos aplicada a características morfológicas e agronômicas em quiabeiro. Horticultura Brasileira, v. 14, n. 2, p. 200-203, 1996.

PENNER, G.A. RAPD analysis of plant genomes. In: JAHUAR, P.P. (ed.) Methods of genome analysis in plants. Dakota: CRC Press, 1996. p. 251-264.

RAO, C.R. Advanced statistical methods in biometric research. New York: John-Wiley and Sons, 1952. 389 p.

SIEMONSMA, J.S. West-African okra morphological and cytogenetical indications for the existence of a natural amphidiploid of Abelmoschus esculentus (L.) Moench and A. manihot (L.) Medikus. Euphytica, v. 31, n. 1, p. 241-252, 1982.

SNEATH, P.H.; SOKAL, R.R. Numerical taxonomy: The principles and pratice of numerical taxonomy. San Francisco: W.H. Freeman, 1973. 573 p.

VASCONCELOS, M.J.V. Avaliação da variabilidade genética de cultivares de feijão (Phaseolus vulgaris L.) pelo uso de marcadores moleculares RAPD. Viçosa: UFV. Imprensa Universitária, 1995, 54 p. (Tese mestrado).

WILLIAMS, J.G.K.; KUBELIK, A.R.; LIVAK, K.J.; RAFALSKI, J.A.; TINGEY, S.V. DNA polymorphisms amplified by arbitrary primers are useful as genetic markers. Nucleic acid Research, v. 18, p. 6531-6535, 1990.

UENF, CCTA, LMGV, Av. Alberto Lamego, 2000 – Horto, 28.0150620 Campos dos Goytacazes-RJ.