Resumos

GENÉTICA

Diversidade genética estimada com marcadores ISSR em populações brasileiras de Zabrotes subfasciatus

Genetic diversity estimated through ISSR markers in populations of Zabrotes subfasciatus

Giselle Anselmo de Souza^I; Márcia Rodrigues de Oliveira Carvalho^I; Ernane Ronie Martins^II; Raul Narciso Carvalho Guedes^III; Luiz Orlando de Oliveira^I

^IUniversidade Federal de Viçosa (UFV), Instituto de Biotecnologia Aplicada à Agropecuária, Avenida P.H. Rolfs, s/nº, Campus Universitário, CEP 365700-000 Viçosa, MG. E-mail: gansouza@gmail.com, mrocarvalho@terra.com.br, luiz.ufv@hotmail.com

^IIUniversidade Federal de Minas Gerais, Instituto de Ciências Agrárias, Setor de Fitotecnia, Avenida Osmani Barbosa, s/nº, Bairro JK, CEP 39404-006 Montes Claros, MG. E-mail: ernane-martins@ufmg.br

^IIIUFV, Departamento de Biologia Animal. E-mail: guedes@ufv.br

RESUMO

O objetivo deste trabalho foi estimar a diversidade genética de populações de Zabrotes subfasciatus, por meio de marcadores moleculares ISSR. Foram avaliadas 12 populações, provenientes de oito estados brasileiros, no total de 269 indivíduos. Cinco iniciadores ISSR permitiram a obtenção do total de 51 fragmentos polimórficos. A percentagem média de locos polimórficos, dentro de cada população, foi de 83,8%. A heterozigosidade corrigida de Nei esperada variou de 0,23 a 0,33, com média de 0,29, e o índice de diversidade genética de Shannon e Weaver variou de 0,29 a 0,48, com média de 0,42. No nível de espécie, estes dois índices apresentaram valores de 0,36 e 0,54, respectivamente. A análise de variância molecular mostrou que 66% da variância molecular total pode ser atribuída a diferenças intrapopulacionais e que os 34% restantes podem ser atribuídos a diferenças interpopulacionais. O teste de Mantel mostrou baixa correlação entre: distância geográfica e diferenciação genética, identidade genética e diferenciação genética e, distância genética de Nei e diferenciação genética. As populações brasileiras de Z. subfasciatus possuem baixa diferenciação genética e fraca estruturação geográfica.

Termos para indexação: caruncho-do-feijão, diferenciação genética, índice de diversidade genética, marcadores ISSR.

ABSTRACT

The objective of this work was to evaluate the genetic diversity of populations of Zabrotes subfasciatus through ISSR molecular markers. Twelve populations from eight Brazilian states were sampled in the total of 269 individuals. Five ISSR primers were used and a total of 51 polymorphic bands were obtained. The percentage of polymorphism within population averaged 83.8%. Nei's unbiased expected heterozygosity ranged from 0.23 to 0.33, with an average of 0.29, and the Shannon & Weaver's index of gene diversity varied from 0.29 to 0.48, with an average of 0.42. At species level, these two indexes were 0.36 and 0.54, respectively. The analysis of molecular variance showed that 66% of the total molecular variance can be attributed to differences within population, while 34% was apportioned among populations. Mantel's test showed a low correlation between: geographic distance and genetic differentiation, genetic identity and genetic differentiation, and Nei's genetic distance and genetic differentiation. The genetic differentiation in Brazilian populations of Z. subfasciatus is low and geographic structuring is weak.

Index terms: bean bruchid, genetic differentiation, genetic diversity index, ISSR markers.

Introdução

O Brasil se destaca mundialmente como o maior produtor de feijão (Phaseolus vulgaris), entretanto, a cultura ainda apresenta baixa produtividade, quando comparada à de outros países (Abreu, 2005; Borém & Carneiro, 2006), o que resulta de uma combinação de fatores, entre os quais estão o baixo nível tecnológico empregado no cultivo e a falta de cuidados fitossanitários adequados, tanto durante o ciclo da cultura, quanto no armazenamento dos grãos (Rosolem & Marubayashi, 1994; Borém & Carneiro, 2006). Estima-se que entre 20 e 30% das perdas na produção brasileira possam ser atribuídas às pragas que atacam os grãos armazenados (Magalhães & Carvalho, 1988).

Zabrotes subfasciatus é uma praga que se destaca quanto às perdas durante o armazenamento de grãos de feijão. A espécie é provavelmente originária da América Central e se tornou praga agrícola, quando se estabeleceu e passou a se reproduzir continuamente em grãos armazenados. Difundiu-se pelas regiões tropicais e subtropicais, por meio do comércio de sementes e, atualmente, é cosmopolita e presente nas Américas Central e do Sul, além de estar amplamente distribuída na África, Índia e Região Mediterrânea (Southgate, 1978; Hill, 1990). Os insetos adultos reproduzem-se no campo durante a formação das sementes e, em grãos, durante o período de armazenagem. Com a perfuração causada pelos insetos, os grãos adquirem sabor desagradável e têm seu valor comercial reduzido. Além disso, os insetos consomem as reservas dos cotilédones e comprometem a germinação e o vigor das sementes (Pimbert & Pierre, 1983; Hill, 1990; Abate & Ampofo, 1996; Yokoyama, 2006).

Para ser efetivo, o controle de pragas requer entendimento sobre a variação genética dentro e entre populações e sobre a distribuição geográfica desta variação. As interações bioquímicas e moleculares de Z. subfasciatus com seus hospedeiros têm sido intensamente investigadas (Pereira et al., 1995; Mazzoneto, 2002; Sales et al., 2005), no entanto, pouco se sabe sobre a extensão da subdivisão genética dentro de suas populações e entre suas populações. Estudos sobre a variabilidade genética da espécie são relevantes, pois restrições no fluxo gênico entre populações têm papel importante na evolução da adaptação local, o que inclui desenvolvimento de resistência a inseticidas e estabelecimento de raças hospedeiras (González-Rodríguez et al., 2002).

Os marcadores moleculares podem ser utilizados como ferramenta para estudos de diversidade genética entre indivíduos, dentro e entre populações ou espécies relacionadas. Aebi et al. (2004) conduziram um estudo com Z. subfasciatus e Z. sylvestris, em que utilizaram uma das classes mais polimórficas de marcadores atualmente disponíveis, os microssatélites, ou "simple sequence repeat" - SSR. Entretanto, os autores detectaram uma deficiência significativa de heterozigotos com os iniciadores avaliados, o que sugere a presença de alelos nulos. Alelos nulos de microssatélites são quaisquer alelos cuja amplificação via PCR falha consistentemente e cuja presença é problemática em análises de dados desta classe de marcadores, pois podem levar à falsa interpretação dos resultados nas análises genéticas (Dakin & Avise, 2004). O estudo aqui apresentado foi conduzido por meio dos marcadores dominantes ISSR ("inter simple sequence repeat") que, embora sejam reconhecidos como menos polimórficos, foram escolhidos para evitar os problemas que os alelos nulos de microssatélites poderiam trazer para as análises desejadas. Marcadores ISSR não requerem informações prévias de seqüências de DNA da espécie-alvo, produzem fragmentos com grande reprodutibilidade, quando comparados a outros marcadores com base em PCR não-específico como RAPD (Wolfe & Liston, 1998), por exemplo, e requerem pouca infra-estrutura em termos de equipamento de laboratório para execução dos experimentos.

O objetivo deste trabalho foi estimar a diversidade genética de populações de Z. subfasciatus, por meio de marcadores moleculares ISSR. Os objetivos específicos foram avaliar o nível de variação de ISSRs em Z. subfasciatus, estimar os níveis de partição da variação genética dentro de populações e entre populações, e avaliar se haveria associação significativa entre a distribuição geográfica atual da praga no Brasil e a distância genética no nível de populações.

Material e Métodos

Amostras de feijão (Phaseolus vulgaris) não-expurgado e com sinais de infestação por Zabrotes subfasciatus - presença de ovos ou orifício de emergência de adultos presentes no tegumento dos grãos - foram obtidas de produtores rurais, com histórico de produção própria de sementes e baixa escala de produção, por intermédio das empresas de extensão rural dos estados brasileiros amostrados. Cerca de 250 a 500 g de grãos foram trazidos do campo para o Laboratório de Biologia Molecular II, no Instituto de Biotecnologia Aplicada à Agropecuária (Bioagro), da Universidade Federal de Viçosa, onde foram mantidas individualizadas e em condições de multiplicação do inseto. Após a primeira eclosão dos ovos em laboratório, os insetos ainda vivos foram fixados em etanol 70%, para etapa posterior de extração de DNA. A identificação do inseto foi confirmada por observação visual de características de morfologia externa, principalmente coxa posterior (Pereira & Almeida, 2001), sob microscópio estereoscópico.

Foram avaliados 269 indivíduos de Z. subfasciatus, provenientes de 12 populações, distribuídas em oito estados brasileiros. As coordenadas geográficas dos locais de coleta das amostras (Tabela 1) foram estimadas com base na localização da cidade de origem, por meio do aplicativo Google Earth versão 3.0. Este aplicativo também foi utilizado para o cálculo das distâncias geográficas entre as populações. Os dados de distâncias geográficas foram utilizados para construir uma matriz e correlacioná-la à matriz de diferenciação genética (F_ST) das populações.

Com base em modificações do protocolo adotado por Clark et al. (2001), foi extraído o DNA genômico total de cada inseto individualmente e, após a etapa de extração, o DNA foi ressuspenso em tampão TE. A concentração final da solução de DNA extraído de cada amostra foi obtida pela comparação de alíquotas de cada solução com amostras de DNA do fago λ, cuja concentração era previamente conhecida. O DNA de cada amostra foi então mantido a -20ºC para uso subseqüente.

A escolha dos iniciadores de ISSR começou por meio de um teste preliminar de amplificação realizado em 26 deles, provenientes do conjunto 9 da University of British Columbia. A amplificação foi conduzida com DNA amostrado de quatro indivíduos, escolhidos aleatoriamente e provenientes de duas populações distintas. Somente iniciadores que permitissem a obtenção de fragmentos de alta intensidade, quando avaliados em gel de agarose, foram selecionados para a próxima fase de avaliação. Foram selecionados 15 iniciadores entre os 26 testados; os 11 iniciadores restantes foram descartados do estudo, pois não apresentaram produtos de amplificação ou a amplificação foi de intensidade muito fraca.

Em seguida, foram testadas condições de reação, de forma a otimizar as amplificações para cada um dos 15 iniciadores anteriormente escolhidos. Foram testadas diferentes concentrações dos vários componentes de PCR, inclusive: temperatura de hibridação daqueles pré-selecionados (43 até 56°C), concentração de MgCl₂ (1,5 mM ou 2 mM), concentração de formamida (1 ou 2%) e quantidade de DNA. Os padrões de fragmentos que cada iniciador produziu, em cada combinação de fatores, foram inspecionados visualmente, tendo-se observado a reprodutibilidade e a intensidade dos fragmentos e a presença de polimorfismo. Finalmente, foram escolhidos cinco iniciadores entre os 15 testados (Tabela 2).

As reações de amplificação foram realizadas em volume total de 15 µL e consistiram de: 1% de formamida (iniciadores 808, 811 e 891) ou 2% de formamida (iniciadores 807, 809); 2 mM de MgCl₂ (iniciadores 807, 808, 811 e 891) ou 2,5 mM de MgCl₂ (iniciador 809); 0,5 U de Taq polimerase (Phoneutria); 10 a 20 ng de DNA molde; 0,2 mM de cada dNTP; 0,2 M de iniciador; 10 mM Tris-HCl pH 8; 50 mM KCl; 0,1% de Tween-20; e água ultrapura.

As amplificações foram realizadas em termociclador Modelo 9700. As condições de amplificação foram adaptadas de Huang & Sun (2000): etapa inicial de desnaturação de 94ºC por 5 min, seguida por 35 ciclos de 94ºC por 45 s, 45 s para hibridação do iniciador (com temperatura de hibridação específica para cada um), 72ºC por 1 min e 50 s, e uma etapa adicional de extensão de 72ºC por 7 min. Os produtos da amplificação foram separados por eletroforese em gel de agarose 1,5% em tampão TBE (90 mM Tris-base, 90 mM ácido bórico e 2 mM EDTA), com tensão constante de 110 V. Posteriormente, os géis foram corados em solução de brometo de etídio 0,6 ng mL^-1. Os fragmentos foram visualizados por meio de luz ultravioleta, e as imagens foram digitalizadas em sistema de captura de imagens. O tamanho dos fragmentos amplificados foi estimado por comparação com marcador molecular de 100 pb DNA Ladder.

Os fragmentos de ISSR foram tratados como marcadores dominantes e julgados como caracteres binários: presente (1) ou ausente (0). Em análises subseqüentes, somente foram considerados os fragmentos polimórficos que puderam ser julgados de forma reproduzível em todos os indivíduos analisados. Fragmentos cuja presença era inconsistente ou que foram observados como fragmentos fracos e de pequena reprodutibilidade, em todos os indivíduos, foram excluídos das análises finais. Uma matriz binária (presença ou ausência do fragmento) foi obtida e com ela estimaram-se os índices de diversidade genética dentro de e entre populações. Marcadores dominantes não permitem que os indivíduos heterozigotos sejam distinguidos na amostra estudada. Assim, a freqüência dos alelos em populações com acasalamentos ao acaso, como é o caso das populações de Z. subfasciatus, pode ser inferida pelo equilíbrio de Hardy-Weinberg. O programa POPGENE versão 1.3 (Yeh et al., 1999) foi utilizado para se calcular a estatística da variação genética de cada população pelo pressuposto de que as populações estão em equilíbrio de Hardy-Weinberg. As estimativas incluíram a percentagem de locos polimórficos (P), heterozigosidade esperada (H_E) e índice de diversidade genética de Shannon (I).

Estimou-se o valor de diferenciação genética, F_ST, que representa um índice de diferenciação genética da população e determina as proporções da variabilidade genética total que se referem às diferenças genéticas entre populações e às diferenças entre indivíduos (Mork, 2000).

A análise de variância molecular (AMOVA) foi realizada por meio do programa Arlequin versão 3.1 (Excoffier et al., 2005). A árvore filogenética foi obtida pelo método UPGMA ("unweighted pair group method with arithmetic mean") e pelo programa MEGA versão 3.1 (Kumar et al., 2004), a partir dos valores de F_ST. O teste de Mantel (Mantel, 1967), implementado pela rotina MXCOMP no programa NTSYS-pc, versão 2.2K, com 9.000 permutações ao acaso, foi usado para investigar a correlação entre F_ST, a identidade genética de Nei e a diversidade genética de Nei, com base nas matrizes cofenéticas.

Resultados e discussão

Os cinco iniciadores ISSR selecionados permitiram a obtenção de 52 fragmentos, dos quais apenas um foi monomórfico e, portanto, excluído da análise. O número de fragmentos por iniciador variou entre 7 e 13, mas, em média, foram obtidos 10 fragmentos. O polimorfismo dentro de cada população variou de 74,51% (populações PAR e RIB) a 92,16% (população ALF), com média de 83,82% (Tabela 3). Estes resultados indicam que os cinco iniciadores utilizados foram satisfatórios para a detecção de polimorfismo, nas populações de Z. subfasciatus analisadas.

Pelo acasalamento ao acaso e o equilíbrio de Hardy-Weinberg, os valores da heterozigosidade esperada corrigida de Nei (1978) (H_E) variaram de 0,1281 a 0,3281, com média de 0,270. Os índices de diversidade genética de Shannon & Weaver (1949) (I) variaram de 0,1893 a 0,4805, com média de 0,3870 (Tabela 3); no nível de espécie, estes dois índices apresentaram valores médios de 0,3626 e 0,5381, respectivamente.

Por meio da AMOVA, foi realizada a partição da variação genética em dois níveis: dentro de populações e entre populações. Essa partição da AMOVA mostrou maior variação genética dentro da população, com 66% da variação total, e os 34% da variação genética restante foram observados entre populações. Os valores de variação genética foram significativos para as duas fontes de variação, com probabilidades calculadas por 1.000 permutações ao acaso. O valor de diferenciação genética das populações (F_ST) foi 0,3358, obtido pela AMOVA, e foi significativo (p<0,001) com probabilidades calculadas por 1.000 permutações ao acaso. Os valores de F_ST entre pares de populações variaram de 0,0680 a 0,6165. Todos os valores encontrados entre pares de populações foram significativos a 5% de probabilidade.

A identidade genética de Nei, calculada entre pares de populações, variou de 0,7226 (MOC e PPO) a 0,9452 (PRA e MON), valores estes que, por serem próximos a 1, mostram o alto grau de similaridade genética entre as populações avaliadas. A distância genética de Nei, calculada para inferir a divergência genética entre populações, variou de 0,0640 a 0,3250, e seu menor valor foi observado entre os pares de populações MOC e PRA, enquanto o maior valor foi observado entre os pares MOC e PPO. MOC é uma população originária de Montes Claros, MG, e PRA é originária de Prado, BA, distante quase 500 km de MOC. PPO está localizada em Ponta Porá, MS, a aproximadamente 1.700 km de MOC.

Os valores de distância genética de Nei foram empregados para a construção do dendrograma com o método UPGMA (Figura 1). A população ALF, amostrada em Alta Floresta, MT, mostrou-se completamente diferenciada das demais populações analisadas e formou um grupo exclusivo. Com exceção de ALF, a análise do dendrograma identificou dois subgrupos principais, cujas composições mostram que não há relação clara entre distância genética e distância geográfica. O primeiro subgrupo agregou a população amostrada no Acre (RIB), juntamente com populações encontradas nos estados do Paraná (COL), Mato Grosso do Sul (PPO), Goiás (ORD e ITA), Bahia (CRA) e Espírito Santo (LAT). O segundo subgrupo apresentou discreta relação entre distribuição geográfica das populações e composição do subgrupo, uma vez que ele contém populações amostradas em localidades geograficamente próximas, das quais três são do Norte e Noroeste de Minas Gerais (MOC, MON, PAR) e uma do Sul da Bahia (PRA). O teste de Mantel mostrou baixa correlação entre: os valores da distância geográfica e F_ST, (r = -0,000708, p[Z_rdmd<Z_obs] = 0,5332); identidade genética de Nei e F_ST (r = -0,24190, p[Z_rdmd<Z_obs] = 0,0888); e distância genética de Nei e F_ST (r = 0,25470, p[Z_rdmd<Z_obs] = 0,9223).

Os insetos amostrados para o presente estudo foram obtidos em propriedades rurais com produção própria de sementes e baixa escala de produção. Deliberadamente, não foram amostrados insetos provenientes de produtores com grandes áreas de plantio. A escolha desta estratégia de coleta teve como objetivo maximizar as chances de amostrar insetos provenientes de populações localmente adaptadas e que, potencialmente, poderiam estar sob menor influência de fluxo gênico associado a dispersões recentes, mediadas por atividades humanas. Foram utilizados somente insetos recebidos vivos ou oriundos da primeira eclosão dos ovos, que vieram juntamente com os grãos de feijão das áreas de cultivo. Insetos que surgiram da multiplicação em condições de laboratório foram descartados do estudo, para não alterar as freqüências gênicas originalmente presentes nas populações multiplicadas em condições de campo.

Estudos que tratam da capacidade de dispersão de populações naturais de bruquídeos, especialmente Z. subfasciatus, não estão disponíveis na literatura. Dispersão a longa distância, apenas com auxílio de correntes aéreas, parece ser bastante improvável para bruquídeos, pelo menos assim foi considerado para Acanthoscelides pallidipennis (Tuda et al., 2001). Entretanto, dispersão a longa distância mediada por atividades humanas, como por exemplo, a livre movimentação de grãos não-expurgados pelo território nacional, pode ser um mecanismo facilitador, que permite a Z. subfasciatus colonizar de maneira eficaz novas áreas geográficas, que não poderiam ser atingidas facilmente por meios naturais de dispersão. Os resultados obtidos no presente estudo apóiam a conclusão de que as populações de Z. subfasciatus possuem estruturação genética fraca. Eles indicam que populações próximas, em termos de geografia, não são necessariamente aquelas mais similares quanto à genética, e que a diversidade genética intrapopulacional é pequena e, portanto, há um extenso fluxo gênico contemporâneo entre populações geograficamente distantes.

Vários fatores podem ter contribuído para a baixa diversidade genética e fraca estruturação geográfica detectada neste estudo. Z. subfasciatus não é uma espécie nativa do Brasil, e sua introdução no país tem provavelmente origem recente. É comum encontrar espécies introduzidas com baixos níveis de diversidade genética, pois a diversidade encontrada em áreas recém-colonizadas é composta apenas de uma pequena parcela da diversidade genética existente nos centros de origem, em razão do efeito do fundador (Kim & Sappington, 2004). Z. subfasciatus se reproduz continuadamente durante o ano, inclusive em hospedeiros alternativos a P. vulgaris e P. lunatus, seus hospedeiros primários (Hill, 1990). Embora os mecanismos de dispersão a longa distância não sejam conhecidos, uma alta mobilidade local permitiria que indivíduos de Z. subfasciatus estabelecessem migrações entre lavouras próximas, o que impediria, assim, o aparecimento de populações geneticamente diferenciadas. Embora este estudo tenha usado somente insetos provenientes de pequenas propriedades rurais, com produção própria de sementes, é possível que outros aspectos não avaliados, como por exemplo, a fonte original das sementes e a proximidade da lavoura das vias de transporte de grãos ou de outras lavouras com histórico de produção distinto tenham influenciado a amostragem de populações geneticamente isoladas.

Conclusões

1. A variabilidade genética entre as populações brasileiras de Zabrotes subfasciatus é baixa.

2. A variação genética intrapopulacional (66%) é maior do que a variação genética interpopulacional (34%).

3. As populações de Zabrotes subfasciatus avaliadas não se apresentam geograficamente estruturadas no Brasil.

Agradecimentos

À Fundação de Amparo à Pesquisa do Estado de Minas Gerais, pelo apoio financeiro concedido a Luiz Orlando de Oliveira, por meio do PPM CAG APQ-0673-3.08/07; à Coordenação de Aperfeiçoamento de Pessoal de Nível Superior e ao Conselho Nacional de Desenvolvimento Científico e Tecnológico, pela concessão de bolsas de estudos de pós graduação a Giselle Anselmo de Souza e Márcia Rodrigues de Oliveira Carvalho, respectivamente.

Recebido em 14 de fevereiro de 2008 e aprovado em 17 de junho de 2008

Datas de Publicação

Histórico