Download PDF
ads:
UNIVERSIDADE FEDERAL DE MINAS GERAIS
INSTITUTO DE CIÊNCIAS BIOLÓGICAS
DEPARTAMENTO DE BIOLOGIA GERAL
PROGRAMA DE PÓS-GRADUAÇÃO EM GENÉTICA
DISSERTAÇÃO DE MESTRADO
Estrutura e diversidade genética de populações de minhocuçu
Rhinodrilus alatus, Righi 1971 (Glossoscholecidae, Clitellata) do
estado de Minas Gerais, Brasil
ORIENTADA: Flávia de Faria Siqueira
ORIENTADORA: Maria Raquel Santos Carvalho
BELO HORIZONTE
Junho - 2009
ads:
Livros Grátis
http://www.livrosgratis.com.br
Milhares de livros grátis para download.
UNIVERSIDADE FEDERAL DE MINAS GERAIS
INSTITUTO DE CIÊNCIAS BIOLÓGICAS
DEPARTAMENTO DE BIOLOGIA GERAL
PROGRAMA DE PÓS-GRADUAÇÃO EM GENÉTICA
Estrutura e diversidade genética de populações de minhocuçu
Rhinodrilus alatus, Righi 1971 (Glossoscholecidae, Clitellata) do
estado de Minas Gerais, Brasil
Flávia de Faria Siqueira
Dissertação submetida ao programa de Pós-
graduação em Genética rea de
Concentração em Genética Evolutiva e de
Populações) da Universidade Federal de Minas
Gerais como requisito parcial para a obtenção
do grau de Mestre em Genética.
ads:
Aos meus pais e ao meu companheiro Eduardo
Charge publicada na “Punch Magazineno dia 6 de dezembro de 1881, logo após o
lançamento do livro de Charles Darwin, “The formation of vegetable mould through the
action of worms, with observations on their habits. Sua obra sobre minhocas foi a
publicação de maior sucesso durante sua vida, vendendo cerca de 3500 cópias em poucos
dias.
v
Agradecimentos
À minha orientadora Pro Dra. Maria Raquel Santos Carvalho, pela orientação e
ensinamentos que contribuíram para a realização deste trabalho.
À banca examinadora, por ter aceito o convite.
Aos professores do programa de Pós-graduação em Genética pela formação ao longo do
curso, em especial ao Prof. Evanguedes Kalapothakis, pelo auxílio nos sequenciamentos; à
Profª Cleusa Graça da Fonseca, ao Profº Eduardo Tarazona Santos, à Pro Maria
Bernadete Lovato, pelas sugestões dadas nas análises dos dados e interpretação.
Aos coordenadores do programa de Pós-graduação em Genética, Prof. Vasco Azevedo e
Prof. Evanguedes Kalapothakis, pelo contínuo incentivo e aos funcionários da secretaria, por
cuidar da parte burocrática.
Ao aluno de iniciação científica Sávio Henrique de Cicco Sandes, estagiário presente e
atuante, que trabalhou muito e também me divertiu nos momentos de tensão. À técnica
Marlene de Miranda, pelo preparo de soluções, organização do laboratório permitindo que
os experimentos fluíssem bem e, acima de tudo, pela preocupação e amizade. Aos colegas
do LGHM, de tempos atrás - Lu Lara, Joana, Daiane, Viviane - e os atuais - Latife,
Gutemberg, Raphael, Gustavo Lacorte e Bastos (saudades!!!), Leo, Izinara, e tantos outros -
pela convivência, carinho e amizade.
Ao ProDr. Rogério Parentoni Martins, à Maria Auxiliadora Drumond e a toda equipe do
Projeto Minhocuçu, em especial à Silvia Campos, pela colaboração e coleta das amostras.
Aos colegas de outros laboratórios, em especial: à Tatiana Moura Barroca e demais
responsáveis pelo ABI3130 que passaram pelo LBMM pela paciência e atenção com meus
sequenciamentos, quase que diários; ao Wagner Carlos Magalhães e Moara Machado do
LDGH pela boa vontade e suporte com análises de bioinformática e por ensinar passo a
passo o Pipeline desenvolvido por eles;
Aos colegas e amigos encontrados na Pós-Graduação, sempre presentes, especialmente ao
“Bonde do E3” e à Flavinha e ao Gui, amizade herdada após organização do simpósio.
vi
Às amigas das antigas: Carolzinha e Rê; e aos amigos da PUC, que compreenderam minha
ausência em alguns momentos e mesmo assim ficaram sempre do meu lado.
Acima de tudo, à minha amada família: Pai e Mãe pelo amor incondicional e pelo incentivo
constante, meus irmãos pela torcida e força, e meu grande amor Eduardo, pelo
companheirismo, amor, estímulo e ainda por ajudar com as imagens.
A Deus, força que me guia em todos os momentos da minha vida.
Enfim, a todos que de alguma forma contribuíram para a conclusão de mais essa etapa da
minha vida. “Se muito vale o já feito, mais vale o que será!
vii
Sumário
Lista de Figuras viii
Lista de Tabelas viii
Lista de Abreviaturas, Siglas e Unidades x
RESUMO 12
INTRODUÇÃO GERAL 13
Sobre a espécie Rhinodrilus alatus 13
Sobre a família Glossoscolecidae e o Gênero Rhinodrilus 14
Genética de Populações e da Conservação 15
Estudos genéticos em Annelida 16
Marcadores moleculares escolhidos 17
Manuscript for Molecular Ecology 19
Genetic diversity and population structure in the brazilian giant earthworm Rhinodrilus
alatus (Annelida, Clitellata, Lumbricina, Glossoscolecidae) 19
ABSTRACT 20
INTRODUCTION 21
MATERIALS AND METHODS 22
Ecological data and distribution of species 22
DNA Extraction and Polymerase Chain Reaction 22
Sequencing 24
Data analysis 24
RESULTS 26
Genetic diversity and population genetic structure with COI 26
Genetic diversity and structure in the ITS1-5.8S-ITS2 segment 28
DISCUSSION 30
ACKNOWLEDGMENTS 35
REFERENCES 35
FIGURE LEGENDS 43
CONCLUSÕES 54
REFERÊNCIAS BIBLIOGRÁFICAS 55
viii
Lista de Figuras
Figure 1 50
Figure 2 51
Figure 3 52
Figura 4 53
ix
Lista de Tabelas
Table 1 45
Table 2 47
Table 3 48
Table 4 49
x
Lista de Abreviaturas, Siglas e Unidades
- Diversidade nucleotídica
µL - Microlitros
µM - Micromolar
AMOVA - Analysis of Molecular Variance, Análise de variância molecular
ºC - Graus Celsius
cm - Centímetros
COI - Subunidade I da Citocromo c Oxidase
DMSO - Dimetilsulfóxido
DNA - Deoxyribonucleic acid, Ácido desoxiribonucléico
dNTP - desoxinucleotídeo tri-fosfato
GPS - Global Position System
GTR - General Time Reversible
h - Diversidade haplotípica
HCl - Ácido Clorídrico
IBAMA - Instituto Brasileiro do Meio Ambiente
Ile - Isoleucina
INDEL - Inserção/Deleção
ITS - Internal Transcribed Spacer, Espaço Interno Transcrito
IUCN - União Internacional para a Conservação da Natureza
k - Número médio de diferenças nucleotídicas
KCl - Cloreto de Sódio
Leu - Leucina
Met - Metionina
mg/mL - Miligramas por microlitros
MgCl
2
- Cloreto de Magnésio
Mm - Milimolar
DNAmt – DNA mitocondrial
mv - Median vector
NCBI - National Center for Biotechonology Information
nDNA - DNA nuclear
ng - Nanogramas
NT - Near threaned, Quase ameaçado
Pb - Pares de Base
PCR - Polymerase Chain Reaction, reação em cadeia da polimerase
xi
PEG - Polietilenoglicol
pH - Potencial Hidrogeniônico
Phe - Fenilalanina
RAPD - Randomly Amplified Polymorphic DNA
rRNA - Ácido ribonuclécico Ribossômico
SAMOVA – Spatial Analysis of Molecular Variance
SNP - Single nucleotide polymorphism, Polimorfismo de nucleotídeo único
SPR - Subtree Pruning and Regrafting
U - unidades
UFMG - Universidade Federal de Minas Gerais
12
RESUMO
O minhocuçu Rhinodrilus alatus (Righi 1971) é um oligoqueta gigante endêmico do cerrado
central da região sudeste do Brasil, muito usado como isca no mercado da pesca. A
atividade extrativista dessa espécie inclui muitas pessoas, e acarreta impactos sociais e
ambientais. O objetivo desse trabalho foi caracterizar a diversidade genética e estrutura
populacional genética e geográfica de R. alatus através do sequenciamento do gene de
DNA mitocondrial (DNAmt) Subunidade I da Citocromo c Oxidase (COI) e do gene nuclear
de RNAr 5.8S e os espaços internos transcritos (ITS) 1 e 2. Foram coletados 75 indivíduos
R. alatus e 6 R. motucu, representando 21 pontos de coleta. A diversidade genética foi
avaliada baseada na diversidade haplotípica, nucleotídica e número médio de diferenças
nucleotídicas. Estrutura populacional foi caracterizada por AMOVA, Teste de Mantel,
SAMOVA e análises filogenéticas (máxima verossimilhança, inferência Bayesiana e median
joining). Uma extraordinária diversidade genética foi encontrada para R. alatus, com 53
haplótipos para o COI, sugerindo a existência de seis linhagens, e 22 haplótipos para o
segmento ITS1-5.8S-ITS2, sugerindo quatro linhagens. A linhagem mais comum para o
segmento ITS1-5.8S-ITS2 coincide com três linhagens diferentes de COI e três linhagens
menos freqüentes foram concordantes para ambos marcadores. Moderada estruturação
genética das populações foi sugerida pela distribuição espacial das linhagens, valores de
Φ
ST
(COI = 0,65, ITS1-5.8S-ITS2 = 0,63, p < 0,00001) e os resultados do Teste de Mantel
(COI = 0,524, ITS1-5.8S-ITS2 = 0,416, p < 0.001). Estes resultados podem sugerir uma
baixa taxa de dispersão para essa espécie, como descrito para outros oligoquetas. Este
estudo contribuirá no delineamento de ações de manejo e uso sustentável da espécie.
13
INTRODUÇÃO GERAL
Sobre a espécie Rhinodrilus alatus
O minhocuçu Rhinodrilus alatus (Righi, 1971) (Annelida, Clitellata,
Glossoscolecidae), é um oligoqueta gigante endêmico do cerrado central do estado de
Minas Gerais, muito apreciado como isca para pesca. Embora não incluído em listas
internacionais de espécies ameaçadas, como a da IUCN, entre 1995 e 2003 o minhocuçu foi
considerado ameaçado de extinção no Estado de Minas Gerais (na categoria “em perigo”,
por meio da publicação da Deliberação Normativa do Conselho de Política Ambiental
41/1995) e no Brasil (Instrução Normativa do Ministério do Meio Ambiente 03/2003).
Atualmente, após uma revisão, a espécie será incluída sob o status de quase ameaçada
(NT) (Machado et al. 2008).
Os animais atingem cerca de 60 centímetros de comprimento e seu ciclo anual
parece ser diretamente influenciado pelas alterações climáticas, vinculado à sazonalidade
dos períodos chuvosos e secos nos locais de sua ocorrência. De acordo com Drumond
(2008), a estação chuvosa (setembro a março) corresponde à fase de reprodução e
forrageamento, peodo no qual ocorre dispersão, ao passo que na estação seca (abril a
agosto), os indivíduos permanecem enrolados em câmaras de quiescência, sem muito
deslocamento. O minhocuçu é um animal hermafrodita que copula através de transferência
mútua de espermatozóides, sendo a fecundação recíproca e cruzada. Segundo os
extratores, os indivíduos reproduzem uma vez ao ano e cada ovo possui dois filhotes.
A área de ocorrência original da espécie era restrita aos municípios de Paraopeba e
Sete Lagoas (Righi 1971), entretanto uma distribuição mais ampla foi confirmada por
Drumond (2008), abrangendo 17 municípios da região central do estado. Esses animais
ocorrem em diferentes fisionomias vegetais como cerrado, cerradão, veredas, e estão
presentes também em locais modificados pelo homem, como plantações de braquiárias,
canaviais e em eucaliptais.
A atividade extrativista e a comercialização ocorrem desde a década de 30 (Miranda
1987) e traz como característica conflitos legais e sociais constantes entre extratores e
proprietários rurais. Além disso, impactos ambientais - como o uso de fogo, revolvimento do
solo, remoção da vegetação e invasão de reservas ambientais - também contribuem para a
problemática dessa atividade.
Por outro lado, a cadeia produtiva dessa atividade envolve centenas de moradores
da região, que participam na extração, comercialização, e até mesmo na produção de
panelas de barro onde são armazenados os animais até que sejam vendidos. Dessa forma,
14
a extração e a comercialização do minhocuçu representam uma alternativa importante com
fonte de renda para populações humanas da região.
Nesse contexto, surgiu o Projeto Minhocuçu, que buscou elucidar aspectos da
ecologia e biodiversidade desse oligoqueto e ecologia humana das comunidades da região,
assim como propor estratégias de manejo e uso sustentável de R. alatus. O projeto
envolveu diversas instituições como a Universidade Federal de Minas Gerais (UFMG), o
IBAMA, o Ministério Público, Prefeituras locais, dentre outras. Na UFMG, uma parceria foi
firmada entre o Laboratório de Ecologia e Comportamento de Insetos - sob a coordenação
do Prof. Dr. Rogério Parentoni Martins e da Dra. Maria Auxiliadora Drumond - e o
Laboratório de Genética Humana e Médica sob coordenação da Prof. Dra. Maria Raquel
Carvalho -, ambos do Departamento de Biologia Geral do Instituto de Ciências Biológicas.
Além da candidata, participaram também desse projeto a Profª Dra. Cleusa Graça da
Fonseca, a bióloga Sílvia Helena Campos, o biólogo Arthur Queiroz, o biólogo Javan Tarsis
e o aluno de iniciação cienfica Sávio Henrique Cicco de Sandes. O presente trabalho foi
uma das vertentes desse projeto maior e buscou caracterizar a estrutura genético-
populacional e contribuir com informações importantes para futuras ações de proteção da
espécie, uma vez que as populações remanescentes são reservatórios de material genético
a ser iminentemente preservado.
Sobre a família Glossoscolecidae e o Gênero Rhinodrilus
Os oligoquetos da família Glossoscolecidae o endêmicos na América do Sul,
podendo ser encontrados em quase todos os habitats terrestres da região compreendida
entre o curso do Rio Juramento Salado na Argentina até o Paralelo de 15º Norte na
Guatemala. Uma das características do solo onde se encontram as espécies dessa família é
a presença de certa quantidade de umidade e de humus e que a acidez não seja excessiva
(Righi 1971). A família foi inicialmente descrita por Michaelsen em 1900, e após algumas
discussões, Stephenson (1930) estabeleceu a existência de cinco subfamílias:
Glossoscolecinae, Sparganophilinae, Microchaetinae, Homorgastrinae e Criodrilinae.
O gênero Rhinodrilus (Perrier, 1872) possui 44 espécies e subespécies, sendo que a
maioria desses táxons foi descrito por Gilberto Righi (Christoffersen 2007). Entretanto, são
poucos os estudos relacionados à identificação e distribuição da biodiversidade desses
oligoquetos.
15
Genética de Populações e da Conservação
A biodiversidade do planeta vem sofrendo uma rápida diminuição como
consequência de ações antrópicas, seja de maneira direta ou indireta. Dentre essas ações,
destacam-se a destruição de habitats e fragmentação, a super exploração, a poluição e a
introdução de espécies exóticas. Em meio a esse declínio, os tamanhos reduzidos de
populações contribuem com a perda da diversidade genética, podendo, assim, limitar a
capacidade de adaptação às mudanças futuras no ambiente.
A diversidade genética, produzida ao longo dos 3,5 bilhões de anos de evolução do
planeta, é um dos três níveis abordados pela Biologia da Conservação que devem ser
protegidos, sendo os demais a diversidade de espécies e a de ecossistemas (Allendorf et al
2007). A partir daí, atua a Genética da Conservação, que faz uso de teorias e técnicas da
genética para tentar minimizar o risco de extinção de espécies ameaçadas. As principais
atividades propostas por essa disciplina são: o manejo genético de populações reduzidas,
afim de se preservar a diversidade genética e minimizar endocruzamentos; a resolução de
incertezas taxonômicas e identificação de unidades de manejo; e o uso da genética
molecular na análise forense e na elucidação da biologia das espécies (Frankham et al
2002; O’Brien 1994). Nesse contexto, os estudos de estrutura genética revelam como a
variação genética está distribuída dentro e entre populações e espécies, afim de se
compreender melhor a dinâmica de populações da espécie na natureza ou monitorar
populações em cativeiro (Avise 2004).
A maioria das populações é agrupada em subpopulações menores dentro das quais
o cruzamento entre os indivíduos geralmente ocorre. Quando há essa estruturação, é
comum encontrar diferenciação genética entre as subpopulações, muitas vezes devido aos
efeitos da endogamia, deriva genética e/ou fluxo gênico (Hartl & Clark 2007; Templeton
2006). Assim, estudos com marcadores moleculares são capazes de inferir sobre tais
fatores evolutivos e os processos históricos que modularam a estrutura genética da
população. Desse modo, o conhecimento da variação intra-específica pode contribuir para a
criação de estratégias de conservação de espécies, que é possível a identificação de
unidades de manejo (O’Brien 1994).
Outra abordagem que utiliza a variação existente dentro de espécies é a
Filogeografia, que estuda os princípios e processos que geraram as distribuições
geográficas de linhagens genealógicas. Essa é uma disciplina relativamente recente, que
busca integrar diferentes áreas como a genética molecular, genética de populações,
etologia, demografia, filogenia, paleontologia, geologia e geografia histórica (Avise 2000). O
uso da informação existente em árvores filogenéticas, juntamente com dados de distribuição
e ocorrência de determinadas espécies, pode fornecer importantes insight sobre a história
16
da Terra. James (2004) mostra algumas razões e exemplos para a aplicação da sistemática
e da biogeografia de minhocas na compreensão da história geológica da Terra,
particularmente em relação ao movimento da crosta terrestre. Por serem considerados
organismos com baixa dispersão, a distribuição de minhocas pode ser justificada, muitas
vezes, somente por conexões passadas entre regiões e por eventos de vicariância.
Estudos genéticos em Annelida
O papel fundamental das minhocas para o solo foi inicialmente descrito por Darwin,
em 1881, através do livro The formation of vegetable mould through the action of worms,
with observations on their habits”. As minhocas tem grande importância para a formação do
solo através da capacidade de alterar seu habitat e criar novos ambientes para outros
organismos. Dentre as modificações pode-se citar o aumento da porosidade, formação de
matéria orgânica, alteração da microbiota e fornecimento de nutrientes por meio de suas
fezes (Römbke et al. 2005). Do ponto de vista ecológico, são organismos fundamentais em
algumas redes alimentares, sendo presas para vários vertebrados e invertebrados
(Symondson 2000). Além disso, sua importância é ainda mais valorizada, uma vez que são
bioindicadores ambientais, reforçando a necessidade de estudos de estrutura e diversidade
genética, pois estes são capazes de fornecer informações valiosas sobre o táxon estudado.
Embora muito relevante, estudos moleculares com esse xon são poucos. De modo
geral, em Annelida, eles se baseiam na filogenia molecular (McHugh 2005; Erseus 2005;
Struck et al. 2007; Huang et al. 2007) e em filogeografia e genética de populações (Jolly et
al. 2005; Bastrop & Blank 2006; Jolly et al. 2006; King et al. 2008; Chang et al. 2008).
Entretanto, a maioria desses estudos são realizados com espécies de poliquetas. Estudos
moleculares em nível intraespecífico com oligoquetos já foram desenvolvidos utilizando
isoenzimas (Haimi et al. 2007), RAPD (Kautenburger 2006; Lentzsch & Golldack 2006),
microssatélites (Harper et al. 2006; Velavan et al. 2007) e haplótipos de mtDNA (Heethoff et
al. 2004; Field et al. 2007 e Cameron et al. 2008).
O estudo da diversidade genética, estrutura de populações e filogeografia de
populações de minhocuçus é inédito, não havendo seqüências disponíveis dessa espécie
em banco de dados. Recentemente, foram depositadas no Barcode of Life Data System
sequências do gene da Subunidade I da Citocromo c Oxidase para três indivíduos do
gênero Rhinodrilus coletados na Amazônia (www.barcodinglife.org). Entretanto, tais
sequências não estão disponíveis para a comunidade. A espécie mais próxima com
sequências disponíveis é a Lumbricus rubellus (Hoffmeister, 1843) com mais de 20.000
sequências (NCBI 2009). Desse modo, esse é o primeiro esforço em descrever a estrutura e
diversidade genética de populações de minhocuçu da região central de Minas Gerais
17
Marcadores moleculares escolhidos
O surgimento da reação em cadeia da polimerase (PCR) (Saiki et al. 1985) marcou o
início da revolução na biologia molecular. Aliado a isso, o uso de iniciadores universais
permitiu o sequenciamento do DNA de espécies nunca antes estudadas. O sequenciamento
de algumas regiões genômicas é cada vez mais comum para estudos de diversidade
genética e para inferir modelos de seleção e demográficos (Schlötterer 2004)
Um único locus ou região de DNA não é capaz de capturar totalmente a estrutura
das populações de espécies e sua história evolutiva, sendo ideal o exame de mais de um
loci e regiões gênicas (Templeton 2006). Estudos genético-populacionais e a reconstituição
dos processos demográficos quando baseados em diferentes sistemas moleculares, trazem
uma interpretação mais real dos fatos. Essa abordagem multiloci permite acessar a
evolução molecular de acordo com diferentes taxas evolutivas.
O DNA mitocondrial animal (DNAmt), com poucas exceções, é um molécula circular
com cerca de 15 a 20 kilobases de tamanho composto por 37 genes, sendo que 22 deles
codificam RNAt, 2 RNAr e 13 codificam proteínas relacionadas no transporte de elétrons e
na fosforilação oxidativa (Avise 2004). Algumas características típicas do DNAmt o tornam
ferramentas valiosas para análises genéticas: ocorrência de milhares de cópias dentro das
células, fácil manipulação e amplificação, herança uniparental materna, frequência mínima
de recombinação, e evolução rápida a nível de sequências, devido ao ineficiente mecanismo
de reparo do DNA (Simon et al. 1994; Avise 2004).
Em relão a outras regiões do DNAmt, postula-se que aquelas que codificam
proteínas são regiões úteis para diversos estudos pois em geral não apresentam
deleções/inserções. O gene da subunidade I da Citocromo c Oxidase apresenta vantagens,
como o uso consolidado de primers universais testado em alguns filos de invertebrados
(Folmer et al. 1994), e um grande alcance de sinal filogenético, sendo utilizado tanto na
discriminação de espécies relacionadas, quanto na diferenciação de grupos filogeográficos
dentro da mesma espécie (Hebert 2003). Assim como nos demais genes DNAmt
codificantes a maior incidência de substituições nucleotídicas ocorrem na terceira posição
do códon, sendo que muitas vezes as substituições são ditas silenciosas, pois não
provocam mudança de aminoácido (Simon et al. 1994).
Os genes de RNA ribossômico no núcleo de células eucarióticas geralmente ocorrem
em elementos repetidos em tandem. Cada unidade repetitiva é composta por sequências
altamente conservadas que codificam os genes ribossomais (subunidades 18S, 5.8S e 28S)
e por regiões espaçadoras menores e mais variáveis. O espaço interno transcrito 1 está
localizado entre os genes de RNAr 18S e 5.8S, enquanto o espaço interno transcrito 2 está
18
presente entre o RNAr 5.8S e 28S. Ambos variam de tamanho entre e, às vezes, dentro de
espécies (Lewin 2004). Por conseguinte, esse segmento tem como característica taxas
diferentes de divergência entre as regiões codificadoras e os espaçadores, e uma evolução
em conjunto (Concerted Evolution). Esse processo evolutivo sugere que unidades repetidas
tendem ser homogeneizadas em sequências dentro das espécies ao invés de evoluir
independente entre as regiões (Page e Holmes 2007).
Os espaçadores tem sido aplicados em diferentes abordagens, tanto na distinção
filogenética de espécies proximamente relacionadas (Chen 2002; Lee & Foighill 2003; Meyer
et al. 2008) quanto para comparação entre populações e estudos filogeográficos (Vogler &
DeSalle 1994; Bargues et al. 2006) .
O presente trabalho gerou um artigo científico que será submetido à revista
Molecular Ecology, conforme descrito a seguir.
19
Manuscript for Molecular Ecology
Genetic diversity and population structure in the brazilian
giant earthworm Rhinodrilus alatus (Annelida, Clitellata,
Lumbricina, Glossoscolecidae)
Flávia F. Siqueira,vio H.C. Sandes, Maria A. Drumond, Sílvia H. Campos, Rogério
P. Martins, Cleusa G. Fonseca, Maria Raquel S. Carvalho
Departamento de Biologia Geral, Instituto de Ciências Biológicas, Universidade
Federal de Minas Gerais Av. Antônio Carlos, 6627, Pampulha, Belo Horizonte,
Minas Gerais, Brazil
KEYWORDS
Earthworm, Rhinodrilus alatus, genetic diversity, population structure
RUNNING TITLE
Genetic diversity and population structure of Rhinodrilus alatus
Correspondence: Dra Maria Raquel S. Carvalho, Fax: +55 31 3409-2570, E-mail:
ma.raquel.carvalho@gmail.com
20
ABSTRACT 1
Rhinodrilus alatus (Righi 1971) is a giant worm (“minhocuçu”) endemic of the central 2
savannah of Southeast Brazil, much appreciated as bait for fishing. R. alatus 3
extractive activity includes many people and results in social and environmental 4
impacts. The aims of this work were to characterize R. alatus genetic diversity and 5
population genetic and geographic structure by sequencing of the mitochondrial gene 6
cytochrome oxidase I (COI) and a nuclear segment including 5.8S rRNA gene and 7
the internal transcriber spacers 1 (ITS1) and 2 (ITS2). Sample was composed by 75 8
R. alatus individuals and 6 R. motucu ones collected at 21 sampling sites. Genetic 9
diversity was evaluated based on haplotype diversity, nucleotide diversity and 10
average number of nucleotides differences. Population structure was characterized 11
with AMOVA, Mantel tests, SAMOVA, and phylogenetic analysis (maximum 12
likelihood, Bayesian inference and median joining). Extraordinary genetic diversity 13
was found for this species, with 53 haplotypes for the COI, suggesting the existence 14
of six lineages, and 22 haplotypes for the ITS1-5.8S-ITS2 segment, suggesting four 15
lineages. Most common ITS1-5.8S-ITS2 lineage segregates with three different COI 16
lineages and the three less frequent lineages were coincident for both markers. 17
Moderate population structure was suggested by the spatial distribution of lineages, 18
Φ
ST
values (COI = 0.65, ITS1-5.8S-ITS2 = 0.63, P < 0.00001), and Mantel tests 19
results (COI = 0.524, ITS1-5.8S-ITS2 = 0.416, P < 0.001). These results can be 20
taken as suggestive of low dispersal rate for the species as described for other 21
earthworms. This study will contribute in the delineation of actions for management 22
and sustainable use for the species. 23
24
21
INTRODUCTION 1
Earthworms are important for soil formation because of its ability to change 2
characteristics such as porosity, amount of organic matter, composition of the 3
microbiota, providing nutrients through their faeces (Römbke et al. 2005). This way, 4
they modify their environments and create new habitats for other organisms. Also, 5
they are involved in some fundamental ecological interactions and are prey for many 6
vertebrates and invertebrates (Symondson 2000). They can also be directly applied 7
in bioremediation strategies to promote biodegradation of organic contaminants 8
(Ceccanti et al. 2006; Hickman & Reid 2008). Earthworms had also great prominence 9
in Charles Darwin studies (1881), who regarded them as one of the most important 10
groups of animals in the planet. 11
Brazil is among the countries with greatest earthworm biodiversity in the world, 12
with about 305 species and subspecies described (James & Brown 2007). 13
Earthworm populations have been found in various ecosystems of the country, both 14
in native environments as well as in those modified by man. Giant earthworms longer 15
than 30 cm and with a diameter larger than 1 cm are collectively called “minhocuçus”. 16
It was estimated that more than 50 species of “minhocuçus” inhabit the country 17
(Brown & James 2007; Christoffersen 2007a,b). 18
The minhocuçu” Rhinodrilus alatus (Righi 1971; Clitellata: Glossoscolecidae) 19
is a terrestrial giant oligochaete endemic to the savannah (cerrado) of Southeast 20
Brazil. The genus Rhinodrilus (Perrier, 1872) has at least 44 species and subspecies, 21
most of them described by Gilberto Righi (Christoffersen 2007b). The family 22
Glossoscolecidae has a geographic distribution restricted to South America. 23
R. alatus has been highly exploited as bait for fishing since the 1930s and 24
therefore, was included in local lists of threatened species (Machado et al. 2008). In 25
addition to the presumed impact of the high extraction on the genetic diversity of the 26
22
species, there are constant legal and social conflicts between extractors and 1
landowners. Moreover, environmental impacts such as the use of fire, soil 2
disturbance, removal of vegetation, and invasion of environmental reserves are 3
additional problems related to minhocuçu collection activity. On the other hand, 4
collecting and marketing the animals is an important alternative source of income for 5
human populations in the region. 6
Worldwide, few studies on earthworm diversity or population genetic structure 7
have been carried out (Kautenburger 2004; King et al. 2008; Cameron et al. 2008), 8
none of them in Brazil. The aim of this study was to characterize the genetic structure 9
and phylogeography of the populations of the giant worm R. alatus. 10
MATERIALS AND METHODS 11
Ecological data and distribution of species 12
The collection of individuals was made with the help of the locals. For the 13
population genetic studies, 75 R. alatus individual were sampled in Southeast Brazil. 14
Additionally, 6 R. motucu individuals have been collected from swamps in Midwest 15
Brazil. Sampling points were geocoding (Table 1). In the field, animals were 16
anesthetized with 10% ethanol and thereafter set on 70% ethanol. Species 17
identification was based on the criteria proposed by Righi (1971). 18
DNA Extraction and Polymerase Chain Reaction 19
In the laboratory, muscle fragments from the posterior region of the body wall 20
were dissected for use in molecular analysis and then the specimens were fixed in 21
10% formaldehyde. Muscle fragments were kept in 70% ethanol at 4 °C until DNA 22
extraction, which was performed following a salting out protocol (Miller et al. 1988), 23
adapted according to the following description: 0.4 X 0.4 cm fragments of muscle 24
tissue were transferred to 1.5 mL tubes containing 500 L of SE (75 mM NaCl, 25 25
mM EDTA, pH 8.0) and 25 L 20 % sodium dodecyl sulfate, and 20 L proteinase K 26
23
(20 mg/mL). After an overnight incubation at 56 °C, 50 L 5M NaCl was added. 1
Tubes were vortexed and centrifuged for 10 minutes at 13000 rpm. Then, the 2
supernatant was transferred to another 1.5 mL tube and precipitated with 2 volumes 3
100% ethanol at -20 °C. Tubes were incubated for 20 minutes at -20 °C and 4
centrifuged for 15 minutes at 13000 rpm. The supernatant was discarded and the 5
pellet was washed twice with 70% ethanol at -20 °C. After removing the 70% ethanol 6
of the second wash, tubes remained open at room temperature for 15 minutes in 7
order to let evaporate ethanol rests. DNA was resuspended in 300 L of TE (10 mM 8
TRIS, 1 mM EDTA, pH 7.4). 9
A 1050 bp fragment containing the internal transcribed spacer region 1 (ITS1), 10
the 5.8S rDNA gene, and the internal transcribed spacer region 2 (ITS2) was 11
amplified with primers: ETTS1: 5' TGCTTAAGTTCAGCGGGT 3', ETTS2: 5' 12
TAACAAGGTTTCCGTAGGTGAA 3' (Kane & Rollinson 1994). ETTS1 primer was 13
anchored at the 3' end of 18S rRNA gene and the ETTS2 primer at the 5' region of 14
28S rRNA. PCR reactions were performed in a 50 L final volume with 10 mM Tris-15
HCl pH 8.4, 50 mM KCl, 1.5 mM MgCl2, 0.1% Triton X-100, 0.2 mM of each dNTP, 16
5% DMSO, 1 U Taq DNA polymerase, 0.3 µM of each primer, and approximately 100 17
ng of genomic DNA. Amplification process consisted of an initial denaturation step at 18
95 ºC for 5 minutes followed by 32 cycles consisting of denaturation at 94 ºC for 45 19
seconds, annealing at 60 ºC for 1 minute and elongation at 72 ºC for 2 minutes, 20
followed by a last extension step at 72 ºC for 2 minutes. 21
A 650 bp fragment of COI gene was amplified with universal primers 22
HCO2198: 5 'TAAACTTCAGGGTGACCAAAAAATCA 3' and LCO1490: 5 23
'GGTCAACAAATCATAAAGATATTGG 3' (Folmer et al. 1994). PCR reactions were 24
set in the same concentrations described above except for primers (1 µM each). The 25
24
PCR program consisted of an initial denaturation step at 94 ºC for 4 minutes, 25 1
cycles with denaturation at 94 ºC for 1 minute, annealing at 50 ºC for 1 minute and 2
extension at 72 ºC for 90 seconds, followed by a 5 minutes last extension step at 72 3
°C. 4
The same individuals were used for both ITS1-5.8S-ITS2 and COI PCR 5
amplification. However, some individuals amplified for only one marker, even after 6
multiple attempts, so the sample size differs between the molecular markers used. 7
Sequencing 8
Amplicons were purified with PEG 8000 (Sambrook & Russel 2001). For each 9
individual, two to four sequences were produced, from independent PCR reactions, 10
so that at least one complete sequence was obtained from each direction for all the 11
individuals. DYEnamic ET Dye Terminator Sequencing (GE HealthCare, Little 12
Chalfont, UK) and ABI PRISM ® BigDye ® Terminator v.3.1 Cycle Sequencing 13
(Applied Biosystems, Foster City, USA) kits were used to produce sequences in the 14
MegaBACE ® 1000 (GE Healthcare, USA) or ABI 3130 (Applied Biosystems, USA) 15
automatic sequencer. In the sequencing, the same primers were used as in the PCR 16
reactions. 17
Data analysis 18
Sequences were analyzed with Polyphred software of the Phred-Phrap-19
Consed package (Ewin et al. 1998, Gordon et al. 1998, Nickerson et al. 1997). In 20
order to confirm the polymorphisms, all chromatogram peaks were conferred by 21
visual inspection and compared to a reference sequence, which was the first one 22
obtained for the species by us. 23
Identification of 5.8S rRNA gene was performed by alignment with the Eisenia 24
fetida sequence (accession number EF534709) available in the database of NCBI, 25
using the program ClustalW implemented in MEGA software version 4.0 (Tamura et 26
25
al. 2007). After confirmation of the SNPs, data were subjected to a pipeline to 1
generate input files (W. Magalhães & E. Tarazona-Santos, personal communication) 2
for PHASE v.2.1 software (Stephens et al. 2001) and DnaSP v.4.50.3 package 3
(Rozas et al. 2003). With the PHASE v.2.1 software, phase and haplotypes were 4
inferred with 10,000 numbers of iterations, 100 thinning interval and 1,000 burn-in. 5
Haplotype diversity (h), nucleotide diversity () and average number of 6
nucleotide differences (k) were calculated for the total sample and for each lineage 7
inferred by phylogenetic analysis (see below) using DnaSP 4.50.3 (Rozas et al. 8
2003). Analysis of molecular variance (AMOVA) was based on the index Φ
ST
and 9
was implemented by the program Arlequin v.3.0, using 1,000 permutations (Excoffier 10
et al. 2005). 11
Pairwise comparison with the evolutionary model of Tamura and Nei (1993) 12
was used to establish genetic distances matrixes between all haplotypes (MEGA 4.0 13
software) and between the lineages identified in phylogenetic studies (Arlequin 3.0 14
software). Other evolutionary models were also tested with similar results. 15
Phylogenetic relationships between haplotypes were reconstructed with three 16
different methods: maximum likelihood (Felsenstein 1981) using PhyML v.3.0 17
software (Guindon & Gascuel 2003), Bayesian inference using MrBayes 3.0 18
(Ronquist & Huelsenbeck 2003), and median joining (Bandelt et al. 1999) with 19
Network v.4.5.1.0 software (Available at: http://www.fluxus-technology.com). Best 20
evolutionary model was selected with MrModelTest V.2.3 software (Nylander 2004) 21
based on Akaike Information Criterion. Maximum likelihood trees were constructed by 22
subtree pruning and regrafting (SPR) methods and the branches were supported by 23
approximate likelihood ratios. Parameters for Bayesian inference were based on two 24
independent runs and four Markov chains with 1x10
6
generations for COI and 1.5 25
26
x10
6
generations for the segment ITS1-5.8S-ITS2, with sampling every 300 1
generations. The number of generations was set corresponding to P < 0.01. The first 2
170 trees of COI and the first 250 trees of the segment ITS1-5.8S-ITS2 were 3
discarded as burn-in, retaining for analysis only those trees present in the stationary 4
phase of the runs. Trees were rooted with R. motucu as out group. 5
Spatial analysis of molecular variance (SAMOVA) was used in order to 6
ascertain population genetic structure and identification of potential barriers to 7
dispersal with SAMOVA v.1.0 software (Dupanloup et al. 2002). This software 8
defines groups of geographically homogeneous populations and maximize the 9
proportion of genetic variance due to differences between populations (F
CT
), based 10
on the number of groups (K) defined a priori through simulated annealing 11
procedures. SAMOVA was performed using 100 simulated annealing procedures for 12
K values from 2 to 9. Correlation between genetic and geographic distances was 13
accessed with Mantel test (Mantel 1967), using 1,000 permutations and logarithmic 14
transformations for both genetic and geographical distances with Alleles in Space 15
software (Miller 2005). 16
RESULTS 17
Genetic diversity and population genetic structure with COI 18
Good quality, partial sequences of COI gene (593 bp) were obtained for 69 19
individuals R. alatus and 6 R. motucu. There were 203 nucleotide substitutions and 20
no insertions or deletions. 21
Considering the 203 (34.2%) polymorphic positions found in both species, 183 22
were parsimoniously informative and 59 haplotypes were identified. Most haplotypes 23
had low frequencies and occurred only in one collection point (Tab.1). Through 24
phylogenetic analysis we identified six distinct lineages of R. alatus. Genetic diversity 25
levels for each lineage and for the whole sample are shown in table 2. Estimates for 26
27
h varied from 0.933 (lineage 4) to 1.00 (lineage 2 and 6), for from 0.003 (lineage 6) 1
to 0.040 (lineage 2) and for k from 2 (lineage 6) to 24 (lineage 2). The most frequent 2
amino acid substitution observed was a methionine to isoleucine, which was present 3
in some individuals in all lineages. 4
Three lineages were the most divergent: lineages 4, 5 and 6. The greatest 5
number of amino acid substitutions per lineage was observed in lineage 5 (six 6
substitutions), and 3 of them were unique to this lineage (2 Ile → Met; 1 Leu → Phe). 7
Lineage 6 has 3 non-synonymous substitutions, the same number as the lineage 4. 8
However, an asparagine to lysine substitution was found only in the last lineage. 9
Moreover, substitutions were also observed in other lineages, for example, in the 10
lineage 1 there was an individual with a replacement of a tryptophan by a serine and 11
in the lineage 2 an individual had a substitution of an isoleucine for a tyrosine. 12
AMOVA test showed a fixation index Φ
ST
equal to 0.65 (P < 0.00001). 13
Difference between the 53 haplotypes according to the model of Tamura and Nei 14
(1993) ranged from 0.17% to 23.08%, with an average of 8.82%. Pairwise distance 15
matrix for the lineages (Table 3) showed that most of the observed distances are 16
significant (P < 0.01), while the smallest difference was 59.3% between lineages 1 17
and 2 and the highest was 87.6% between lineages 2 and 4, considering only the 18
significant values. 19
Phylogenetic trees generated with both the maximum likelihood analysis and 20
Bayesian inference showed similar topologies, with few branches not being 21
recovered by both methods (Fig. 1). Approximate likelihood ratios and a posteriori 22
Bayesian probabilities were higher than 0.5 in all branches and showed close 23
similarity between the two methods of phylogenetic reconstruction. The lineages 24
28
identified in the phylogenetic trees were also supported by median joining haplotype 1
network (Fig. 2). 2
It was possible to identify some patterns of geographic distribution of lineages 3
(Fig. 3). Lineage 1 is widely distributed, occurring throughout the center of the 4
sampled area, from north to south and limited by Três Marias Reservoir/Pará River 5
and Das Velhas River. Lineage 2 was present at three points close to each other on 6
both riversides of Das Velhas River. More evolutionarily divergent lineages, such as 7
3, 4, 5, and 6, had restricted distributions, respectively, to points 20, 1, 13, and 11/13. 8
The highest F
CT
value identified with SAMOVA was observed under K = 2 9
model, separating individuals from point 1 (UTM1 470914, UTM2 7831039) from the 10
other (F
CT
= 0.715, P < 0.05). F
CT
values diminished as K values increased (Table 4) 11
and the proposed groupings were not congruent with the lineages identified. Thus, a 12
possible genetic or geographical barrier should be considered separating point 1 13
from the other sampling sites. A significant and positive correlation was detected 14
between genetic and geographic distances by Mantel test with r equal to 0.524 (P < 15
0.001). 16
Genetic diversity and structure in the ITS1-5.8S-ITS2 segment 17
In order to investigate a nuclear DNA segment, a 935 bp fragment was PCR 18
amplified, which included 492 bp of the ITS1, 157 bp of the 5.8S rRNA gene, and 286 19
bp of ITS2. In this segment, nucleotide substitution and insertion/deletion (INDEL) 20
polymorphisms were identified. For example, ITS1 size ranged from 447 to 492 bp 21
(INDEL blocks varied from 1 to 8 nucleotides) and ITS2 size went from 282 to 286 bp 22
(with only one, 4 nucleotides INDEL block). There were no INDELS in 5.8 rRNA 23
gene. 24
29
Nucleotide substitutions had the following distribution: in ITS1 there were 48 1
transitions and 31 transversions, in ITS2 there were 19 transitions and 18 2
transversions and in 5.8S rRNA there was only one transition. 3
The most common ITS1-5.8S-ITS2 haplotype (H11) is distributed over 7 4
sampling points, while the other ones had restricted distributions or were unique to a 5
sampling point (Table 1). Based on the genetic diversity present in this gene segment 6
and phylogenetic analysis, R. alatus was subdivided into four lineages. Lineage 1 7
showed a greater number of sequences (126), but has only 18 haplotypes. h 8
remained between 0.0005 (Lineage 3) and 0.818 (Lineage 1), ranged from 0.0005 9
(Lineage 3) to 0.003 (Lineage 1) and k was 0.5 (Lineage 3) to 3.470 (Lineage 1) 10
(Table 2). 11
With AMOVA, most of the variance was located between populations with Φ
ST
12
equal to 0.637 (P < 0.05). Genetic divergence between the 22 ITS1-5.8S-ITS2 13
haplotypes observed in R. alatus ranged from 1x10
-6
% and 6.94%, with an average 14
of 2.4%. Most of the distances observed in the pairwise distance matrix between the 15
lineages were significant (P <0.01) (Table 3). The smallest genetic divergence 16
occurred between lineages 1 and 2 (91.7%) and highest between 2 and 4 (100%), 17
considering only the significant values. 18
Phylogenetic reconstructions with both maximum likelihood and Bayesian 19
inference, produced highly similar results for ITS1-5.8S-ITS2 considering topologies, 20
with few branches not being recovered by one of the methodologies used (Fig. 1). 21
Approximate likelihood ratios and a posteriori Bayesian probabilities rendered similar 22
values (greater than 0.5) by phylogenetic reconstruction with the ITS1-5.8S-ITS2. 23
Five lineages were obtained in the phylogenetic reconstructions as well as in the 24
median joining networking of the ITS1-5.8S-ITS2 haplotypes (Fig. 4). 25
30
Lineage 1 showed a very wide distribution, being absent only on the left 1
riverside of the Pará River (Fig. 3). Other lineages had more restricted spatial 2
distribution, e. g., lineage 2 present only at sampling site 1. 3
Similarly to COI results, SAMOVA indicated a higher value of F
CT
(= 0.821, P < 4
0.05) to K equal to 2, separating site 1 from the others. Lower but still significant F
CT
5
values were obtained for K values 3 to 9 (Table 4). Thus, the same putative genetic 6
and/or geographical barrier was observed with both COI and ITS1-5.8S-ITS2. A 7
significant and positive correlation between genetic and geographic distances was 8
identified by Mantel test (r = 0.416, P <0.001). 9
Besides, the 6 COI and the 4 ITS1-5.8S-ITS2 lineages were compared in 10
terms of intra-individual and spatial distribution (for a better understanding the 11
nomenclature correlating lineages of both markers will be: lineage-marker/lineage-12
marker, for example 6-COI/3-ITS). At intra-individual level and spatial distribution 13
there was a perfect coincidence between 4-COI/2-ITS (sampling site 1), 5-COI/4-ITS 14
(sampling site 13), and 6-COI/3-ITS (sampling sites 11 and 13). ITS lineage 1 comes 15
together in individuals with COI lineages 1, 2, and 3 (Fig. 3). 16
DISCUSSION 17
This work represents the first effort to describe relevant aspects about the 18
biology of a giant oligochaete with great local ecological, economic and social 19
importance. There are reports of the use of “minhocuçu” R. alatus as bait for fishing 20
since the 1930s. Consequently, it was included in regional and national lists of 21
endangered species. A better understanding of the role of these invertebrates in 22
ecosystem processes will help delineating conservation strategies for the species 23
(Lewinsohn et al 2005; Dupont 2009). 24
Currently, just few studies have investigated genetic/geographic population 25
structure based on sequencing results in oligochaete (Heethoff et al. 2004, Field et 26
31
al. 2007, Cameron et al. 2008, King et al. 2008, review by Dupont 2009). In the 1
present study, molecular analysis revealed a high intra-specific genetic diversity in R. 2
alatus with 31.87% polymorphic sites in COI and 15.4% in the segment ITS1-5.8S-3
ITS2. Similar values were described for COI gene in Allolobophora chlorotica with 4
30.2% of polymorphic sites (King et al. 2008). Lower diversity was observed for the 5
partenogenetic species Dendrobaena octahedron, with 10.7% of polymorphic 6
positions in COI (Cameron et al. 2008). In Octolasion tyrtaeum, 32.53% variable sites 7
were identified in cytochrome oxidase II gene (Heethoff et al. 2004) and in Lumbricus 8
terrestris, 12.44% polymorphic sites were found in subunits 2 and 4 of the NADH 9
gene (Field et al. 2007). This is the first work describing genetic diversity in the whole 10
ITS1-5.8S-ITS2 segment for an oligochaete species. This region has already been 11
investigated in polychaetes of the genus Perinereis, including four species, in which 12
27.6% of variation was found (Chen et al. 2002). 13
The identification of cryptic lineages within populations of R. alatus is 14
consistent with results found in phylogenetic studies of other earthworms. Six 15
different lineages were determined in Tubifex tubifex, by the sequencing of 16S 16
mitochondrial ribosomal gene (Beauchamp 2001) and 5 lineages for A. chlorotica 17
with COI (King 2008). Heethoff (2004) found evidence for two genetic and 18
morphologically distinct lineages of O. tyrtaeum. One cryptic species was identified in 19
a species complex Metaphire formosae in Taiwan (Chang et al. 2008). In other 20
groups of Annelida, two and three lineages cryptic were identified in populations of 21
Pectinaria koreni and Owenia fusiformis, respectively, in the Northeast Atlantic (Jolly 22
et al. 2005; Jolly et al. 2006). Accordingly, the existence of lineages in R. alatus was 23
supported by the following evidences: a) the high genetic divergence observed in the 24
matrices of pairwise distance (Table 3), b) pattern of haplotypes distribution in the 25
32
phylogenetic reconstructions and their branch sizes (Fig. 1) and, c) the large number 1
of mutational steps between the lineages in the haplotype networks (Figs. 2 and 4). 2
We observed six lineages with COI and four with ITS1-5.8S-ITS2. The most 3
frequent lineage was ITS1-5.8S-ITS2 lineage 1, which occurred in individuals bearing 4
COI lineages 1-COI, 2-COI, or 3-COI. Moreover, there were three less frequent 5
lineages for both systems that always coincided at intra-individual level (5-COI/4-ITS, 6
4-COI/2-ITS, and 6-COI/3-ITS). Two of them (5-COI/4-ITS and 6-COI/3-ITS) were 7
identified at collection points where the most common ones were also detected. 8
These two lineages should be treated with caution because they contain only, 9
respectively, one and two individuals each. 10
At AMOVA, significant and high values (63% of variance between populations) 11
were obtained for both markers, suggesting a genetic structure. Higher 12
interpopulational variation was also described for D. octahedron with COI (F
ST
= 13
0.717; Cameron et al. 2008), for L. terrestris based on RAPD (Φ
ST
= 0.424; 14
Kautenburger et al. 2006); and Sabella spallanzanii an invader polichaete with ITS2 15
(F
ST
values from 0.533 to 0.838; Patti & Gambi 2001). It has been suggested that 16
FST values higher than 0.15 evidence significant differentiation between fragments 17
(Frankham et al. 2002). High values of divergence among fragments or populations 18
suggest low dispersal ability. 19
Generally, earthworms are considered low dispersion and highly localized 20
reproduction animals (James, 2004). However, dispersal rates may differ among 21
earthworm species due to many reasons, such as reproduction systems and 22
ecological niches. The results reported here suggest that R. alatus most frequent 23
lineages extend throughout the species geographical distribution area, while those 24
with higher genetic divergence are restricted and peripheral (Fig.3). Highly positive 25
33
and significant correlation values between genetic and geographic distances 1
obtained for both COI and ITS1-5.8S-ITS2 markers at Mantel tests suggest distance 2
as an important factor for genetic differentiation of R. alatus populations. The 3
restricted distribution of genetic diversity observed in some lineages reinforces the 4
idea of low dispersion for this species. These findings contradict the opinions the 5
locals, who report that the “minhocuçu” R. alatus is able to travel long distances over 6
the soil surface. According to Brown and James (2007), except for some species with 7
wide distribution (Pontoscolex corethrurus and Urobenus brasiliensis), about 80% of 8
Brazilian species of earthworms have restricted dispersals, occurring in only a few, 9
nearby locations. Similarly, King et al. (2008) showed that the relatively sedentary 10
behavior of A. chlorotica contributed to genetic differentiation of lineages and 11
subsequent sympatric speciation. 12
In addition to the dispersion data, the inference of possible barriers to gene 13
flow may also be useful in understanding the distribution patterns of the species. The 14
results of SAMOVA (Table 4) for both markers showed the existence of a potential 15
barrier between the Point 1 and other sampling points. Other simulated values for K 16
were also significant, however presented lower F
CT
values than those found for K=2. 17
Plotting the sampling points over different maps of the region (vegetation, 18
topography, hydrology, and geological formation), points to the Pará River as a 19
possible geographical barrier, that contributed to the genetic differentiation of this 20
population (Fig. 3). Therefore, SAMOVA with K=2 reinforces the identification of the 21
lineage 4-COI/2-ITS. 22
In this study, a multiloci approach was adopted, using nuclear and 23
mitochondrial markers, for addressing different evolutionary rates. Representing the 24
nuclear genome, it was used a segment including the coding region of the 5.8S rRNA 25
34
gene and parts of the internal transcribed spacers 1 and 2 (ITS1 and ITS2). RNA 1
ribosomal gene complex (rDNA) is a nuclear unit of tandem repeats, with one to 2
many hundreds of copies (Hillis & Dixon 1991). Due to its moderate rate of molecular 3
evolution, ITS sequences are important tools for population genetics and 4
phylogenetic studies (Hillis & Davis 1986, Vogler & DeSalle 1994, Reed et al. 2000, 5
Chen et al 2002). In order to ascertain the mitochondrial genome mutation rate, the 6
region chosen was the gene for subunit I of cytochrome C oxidase (COI). 7
Mitochondrial DNA has proven useful in the reconstruction of historical models of 8
population demography, migration, biogeography and speciation (Brown et al. 1979, 9
Moore 1995, Hebert et al 2003, Hurst & Jiggins 2005). 10
Higher genetic diversity was detected with COI (6 lineages), when compared 11
to ITS1-5.8S-ITS2 segment (4 lineages). COI location in the mtDNA explains its high 12
mutation rate, in spite of being a protein coding gene (Simon et al. 1994). ITS1-5.8S-13
ITS2 segment, due to its large number of copies in tandem, is subject to gene 14
conversion mechanisms typical of such sequences, which reduce their genetic 15
diversity. Furthermore, these sequences are subject to greater selective pressure 16
because they contain signals required for the processing of RNA transcribed (Gerbi 17
1985, Page & Holmes 1998). The results reported here reinforce the usefulness of 18
associating genetic markers with different evolutionary dynamics. 19
An amazing aspect of this study is the high amount of genetic diversity 20
observed in individuals sampled over a relatively small geographic area, supporting 21
the current idea of large genetic diversity among earthworms. Our study intended to 22
elucidate important features of life history and conservation of a species with high 23
local importance. As recommended by Machado et al. (2008), basic research is 24
essential for conservation biology of R. alatus, and to encourage management 25
35
projects for its population. High genetic diversity values observed for this species in 1
such a small area reiterate the necessity of protection of the species. The present 2
results will help delineating conservation. The initial characterization of genetic 3
diversity and geographical distribution of species that will be useful in the follow up of 4
projects for the conservation and sustainable use R. alatus. 5
ACKNOWLEDGMENTS 6
Thank Samuel Davis Wooster, of University of Kansas / USA and George Brown, of 7
EMBRAPA Florestas (Empresa Brasileira de Pesquisa Agropecuária), Curitiba, 8
Brazil, for training in the taxonomic classification of terrestrial oligochaetes. This work 9
received grants from FAPEMIG (Fundação de Amparo à Pesquisa de Minas Gerais), 10
CNPq (National Research Council), Conservation Internacional, and IEF-MG 11
(Instituto Estadual de Florestas de Minas Gerais) from Brazil. 12
REFERENCES 13
Bandelt HJ, Forster P, Rohl A (1999) Median-joining networks for inferring 14
intraspecific phylogenies. Molecular Biology and Evolution, 16, 37–48. 15
Beauchamp KA, Kathman RD, McDowell TS, Hedrick RP (2001) Molecular 16
phylogeny of Tubicid Oligochaetes with special emphasis on Tubifex tubifex 17
(Tubicidae). Molecular Phylogenetics and Evolution, 19, 216–224. 18
Brown GG, James SW (2007) Ecologia, Biodiversidade e biogeografia das 19
minhocas no Brasil. In: Minhocas na América Latina: Biodiversidade e Ecologia (ed. 20
Brown GG, Fragoso C.), pp. 297-381. Embrapa Soja, Londrina, Paraná. 21
Brown WM, George MJ, Wilson AC (1979) Rapid evolution of animal 22
mitochondrial DNA. Proceedings of the National Academy of Sciences of the United 23
States of America, 76, 1967–1971 24
36
Cameron EK, Bayne EM, Coltman DW (2008) Genetic structure of invasive 1
earthworms Dendrobaena octaedra in the boreal forest of Alberta: insights into 2
introduction mechanisms. Molecular Ecology, 17, 1189-1197. 3
Ceccanti B, Masciandaro G, Garcia C, Macci C, Doni S (2006) Soil 4
bioremediation: combination of earthworms and compost for the ecological 5
remediation of a hydrocarbon polluted soil. Water, Air, and Soil Pollution, 177, 383–6
397 7
Chang CH, Lin SM, Chen JH. (2008) Molecular systematics and 8
phylogeography of the gigantic earthworms of the Metaphire formosae species group 9
(Clitellata, Megascolecidae). Molecular Phylogenetics and Evolution, 49, 958-968. 10
Chen CA, Chen CP, Fan TY, Yu JK, Hsieh HL (2002) Nucleotide sequences of 11
ribosomal internal transcribed spacer and their utility in distinguishing closely related 12
Perinereis Polychaets (Annelida; Polychaeta; Nereididae). Marine Biotechnology, 4, 13
17-29. 14
Christoffersen ML (2007)a Biodiversity and distribution of the Glossoscolex 15
genus-group (Annelida, Clitellata, Lumbricina, Glossoscolecidae) in South America. 16
Neotropical Biology and Conservation, 2, 155-160. 17
Christoffersen ML (2007)b Distribution and species diversity of Rhinodrilus 18
Perrier, 1872 (Annelida, Clitellata, Lumbricina, Glossoscolecidae) in South America. 19
Neodiversity, 2, 1-6. 20
Darwin C (1881) The Formation of Vegetable Mould Through the Actions of 21
Worms. John Murray, London. 22
37
Dupanloup I, Schneider S, Excofer L (2002) A simulated annealing approach 1
to dene the genetic structure of populations. Molecular Ecology, 11, 2571–2582. 2
Dupont L (2009) Perspectives on the application of molecular genetics to 3
earthworm ecology. Pedobiologia, 52, 191-205. 4
Excoffier L, Laval G, Schneider S (2005) Arlequin (version 3.0): An integrated 5
software package for population genetics data analysis. Evolutionary Bioinformatics 6
Online, 1, 47–50. 7
Ewing B, Hillier L, Wendl MC, Green P (1998) Base-calling of automated 8
sequencer traces using phred. I. Accuracy assessment. Genome Research, 8,175-9
185. 10
Felsenstein J (1981) Evolutionary trees from DNA sequences: A maximum 11
likelihood approach. Journal of Molecular Evolution, 17, 368-376. 12
Field SG, Lange M, Schulenburg H, Velavan TP, Michiels NK. (2007) Genetic 13
diversity and parasite defense in a fragmented urban metapopulation of earthworms. 14
Animal Conservation, 10, 162-175. 15
Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R (1994) DNA primers for 16
amplification of mitochondrial cytochrome c oxidase subunit I from diverse 17
metazoan invertebrates. Molecular Marine Biology and Biotechnology, 3, 294299. 18
Frankham R, Ballou JD, Briscoe DA (2002) Introduction to conservation 19
genetics. Cambridge University Press, Cambridge. 20
Gerbi S (1985) Evolution of ribosomal RNA. In: Molecular Evolutionary 21
Genetics, MacIntyre, R.J. (ed). New York, N.Y.: Plenum, 419–517. 22
38
Gordon D, Abajian C, Green P (1998) Consed: A graphical tool for sequence 1
finishing. Genome Research, 8, 195-202. 2
Guindon S, Gascuel O (2003) A Simple, fast, and accurate algorithm to 3
estimate large phylogenies by maximum likelihood. Systematic Biology, 52, 696-704. 4
Hebert PDN, Cywinska A, Ball SL, DeWaard JR (2003) Biological identifications 5
through DNA barcodes. Proceedings of the Royal Society B: Biological Sciences, 6
270, 313–321. 7
Heethoff M, Etzold K, Scheu S (2004) Mitochondrial COII sequences indicate 8
that the parthenogenetic earthworm Octolasion tyrtaeum (Savigny 1826) constitutes 9
of two lineages differing in body size and genotype. Pedobiologia, 48, 9–13. 10
Hickman ZA, Reid BJ (2008) Earthworm assisted bioremediation of organic 11
contaminants. Environment International, 34, 1072-1081. 12
Hillis DM, Davis SK (1986) Evolution of ribosomal DNA: fifty million years of 13
recorded history in the frog genus Rana. Evolution, 40, 1275-1288. 14
Hillis DM, Dixon MT (1991) Ribosomal DNA: molecular evolution and 15
phylogenetic inference. Quarterly Review of Biology, 66, 411–453. 16
Hurst GDD, Jiggins FM (2005) Problems with mitochondrial DNA as a marker in 17
population, phylogeographic and phylogenetic studies: the effects of inherited 18
symbionts. Proceedings of the Royal Society B: Biological Sciences, 272, 1525-1534. 19
James SW (2004) Planetary process and their interactions with earthworm 20
distribution and ecology. In: Earthworm Ecology (ed. Edwards CA), pp. 53–62. CRC 21
Press, Boca Raton, Florida. 22
39
Jolly MT, Jollivet D, Gentil F, Thiébaut E, Viard F (2005) Sharp genetic break 1
between Atlantic and English Channel populations of the polychaete Pectnaria 2
koreni, along the north coast of France. Heredity, 94, 23-32. 3
Jolly MT, Viard F, Gentil F, Thiébaut E, Jollivet D. (2006) Comparative 4
phylogeography of two coastal polychaete tubeworms in the Northeast Atlantic 5
supports shared history and vicariant events. Molecular Ecology, 15, 1841-1855. 6
Kane RA, Rollinson D (1994) Repetitive sequences in the ribosomal DNA 7
internal transcribed spacer of Schistosoma haematobium, Schistosoma intercalatum 8
and Schistosoma mattheei. Molecular and Biochemical Parasitology, 63, 153-156. 9
Kautenburger R (2006) Genetic structure among earthworms (Lumbricus 10
terrestris L.) from different sampling sites in western Germany based on random 11
amplied polymorphic DNA. Pedobiologia, 50, 257-266. 12
King RA, Tibble AL, Symondson WOC (2008) Opening a can of worms: 13
unprecedented sympatric cryptic diversity within British lumbricid earthworms. 14
Molecular Ecology, 17, 4684-4698. 15
Lewinsohn TM, Freitas AVL, Prado PI (2005) Conservação de invertebrados 16
terrestres e seus habitats no Brasil. Megadiversidade, 1, 62-69. 17
Machado ABM, Drummond GM, Paglia AP (2008) Livro Vermelho da Fauna 18
Brasileira Ameaçada de Extinção. Fundação Biodiversitas, Belo Horizonte, Minas 19
Gerais. 20
Mantel N (1967) The detection of disease clustering and a generalized 21
regression approach. Cancer Research, 27, 209-220. 22
40
Miller SA, Dykes DD, Polesky HF (1988) A simple salting out procedure for 1
extracting DNA from human nucleated cells. Nucleic Acids Research, 16, p.1215. 2
Miller MP (2005) Alleles In Space (AIS): Computer software for the joint 3
analysis of interindividual spatial and genetic information. Journal of Heredity, 96, 4
722-724. 5
Moore WS (1995) Inferring phylogenies from mtDNA variation: mitochondrial 6
gene trees versus nuclear gene trees. Evolution, 49, 718–726. 7
Nickerson DA, Tobe VO, Taylor SL (1997) PolyPhred: Automating the detection 8
and genotyping of single nucleotide substitutions using fluorescence-based 9
resequencing. Nucleic Acids Research, 25, 2745-2751. 10
Nylander JAA (2004) MrModeltest v2. Program distributed by the author. 11
Evolutionary Biology Centre, Uppsala University,Sweden. 12
Page RDM, Holmes EC (1998) Molecular Evolution: A Phylogenetic Approach. 13
WileyBlackwell Science Inc., Oxford. 14
Patti FP, Gambi MC (2001) Phylogeography of the invasive polychaete Sabella 15
spallanzanii (Sabellidae) based on the nucleotide sequence of internal transcribed 16
spacer 2 (ITS2) of nuclear rDNA. Marine Ecology Progress Series, 215, 169-177. 17
Righi G (1971) Sobre a família Glossoscolecidae (Oligochaeta) no Brasil. 18
Arquivos de Zoologia, 65, 1-95. 19
Reed KM, Hackett JD, Phillips RB (2000) Comparative analysis of intra-20
individual and inter-species DNA sequence variation in salmonid ribosomal DNA 21
cistrons. Gene, 249, 115-125. 22
41
Römbke J, Jänsch S, Didden W (2005) The use of earthworms in ecological soil 1
classification and assessment concepts. Ecotoxicology and Environmental Safety, 2
62, 249-265. 3
Ronquist F, Huelsenbeck JP (2003) MrBayes 3: Bayesian phylogenetic 4
inference under mixed models. Bioinformatics, 19, 1572–1574. 5
Rozas J, Sánchez-DelBarrio JC, Messeguer X, Rozas R (2003) DnaSP, DNA 6
polymorphism analyses by the coalescent and other methods. Bioinformatics, 19, 7
2496-2497. 8
Sambrook J, Russel DW (2001) Molecular cloning: a laboratory manual. Cold 9
Spring Harbor Laboratory, New York. 10
Simon C, Frati F, Beckenbach A, Crespi B, Liu H, Flook P (1994) Evolution, 11
weighting, and phylogenetic utility of mitochondrial gene sequences and a 12
compilation of conserved polymerase chain reaction primers. Annals of the 13
Entomological Society of America, 87, 651-701. 14
Stephens M, Smith NJ, Donnelly P (2001) A new statistical method for 15
haplotype reconstruction from population data. American Journal of Human Genetics, 16
68, 978-989. 17
Symondson WOC , Glen DM, Erickson ML, Liddell JE, Langdon CJ (2000) Do 18
earthworms help to sustain the slug predator Pterostichus melanarius (Coleoptera: 19
Carabidae) within crops? Investigations using monoclonal antibodies. Molecular 20
Ecology, 9, 1279-1292 21
42
Tamura K, Nei M (1993) Estimation of the number of nucleotide substitutions in 1
the control region of mitochondrial DNA in humans and chimpanzees. Molecular 2
Biology and Evolution, 10, 512-526. 3
Tamura K, Dudley J, Nei M, Kumar S (2007) MEGA4: Molecular evolutionary 4
genetics analysis (MEGA) software version 4.0. Molecular Biology and Evolution, 24, 5
1596-1599. 6
Vogler AP, DeSalle R (1994) Evolution and Phylogenetic Information Content of 7
the ITS-l Region in the Tiger Beetle Cicindela dorsalis. Molecular Biology and 8
Evolution, 11, 393-405. 9
10
43
FIGURE LEGENDS 1
2
FIGURE 1 Phylogenetic tree of 53 haplotypes of COI and the 22 haplotypes of ITS1-3
5.8S-ITS2. Topologies shown were based on the maximum likelihood method and 4
the branches not reconstituted by Bayesian inference have an asterisk. Respectively, 5
the values of Bayesian posterior probabilities and the approximate likelihood ratios 6
are indicated next to branch. The scale represents 0.1 substitutions per site for COI 7
tree and 0.02 for segment ITS1-5.8S-ITS2 tree. 8
9
FIGURE 2 Haplotype network of COI. Colored ellipses emphasize the lineages 10
identified and their numbers are next to them. The white circles represent 53 11
haplotypes found and their size is proportional to the frequency of them. Small gray 12
circles (mv = mediam vector) represent hypothetical haplotypes. The size of the 13
branches is not proportional to genetic divergence and the numbers above them 14
indicate the number of mutational steps between the different lineages. 15
16
FIGURE 3 Map of Brazil showing the location of Southeast region. R. alatus sampling 17
points are shown in the rectangles. COI lineages are shown in the top and ITS1-18
5.8S-ITS2 lineages in botton rectangle. Numbers refer to sampling sites (Table 1). 19
The colors used to represent the lineages in both cases are the same used in figure 20
1. Observe that ITS1-5.8S-ITS2 lineage 1-ITS (in violet) occurs over the central 21
geographic distribution area and coincided at intra-individual level with COI lineages 22
1-COI (orange), 2-COI (grey), and 3-COI (pink). Less frequent lineages are 23
represented by colors blue (5-COI/4-ITS), green (4-COI/2-ITS), and yellow (6-COI/3-24
ITS) for both markers. Note the most periphery distribution of these lineages. 25
26
44
FIGURE 4 Haplotype network of ITS1-5.8S-ITS2. Colored ellipses illustrate the 1
lineages identified and their numbers next to them. Lineages 2, 3 and 4 have the 2
same color of the lineages 4, 6 and 5 of the COI, respectively, because they 3
represent the same individuals. The white circles represent the 22 haplotypes found 4
and their size is proportional to the frequency of them. Small gray circles (mv = 5
mediam vector) represent hypothetical haplotypes. The size of the branches is not 6
proportional to genetic divergence and the numbers above them represent the 7
number of mutational steps between the different lineages. 8
45
TABLE 1
Characterization of the sampling sites, haplotypes and lineages present. Part a represents COI and b ITS1-5.8S-ITS2
segment. Ncr = number of chromosomes, NH = number of haplotypes, Lin. Pres. = Lineages present.
a) COI
Sampling
Site UTM1 UTM2 Ncr
NH
Lin.
Pres.
Lin. 1 Lin. 2
Lin.
3
Lin. 4
Lin.
5
Lin.
6
1 470914
7831039
6
5
4
H49,50,51,52,53
2 585530
7848179
4
4
1,2
H40,21 H35,36
3
590033
7847313
1
1
2
H34
4 535291
7849567
5
3
1
H42,43,44
5 539341
7857588
5
4
1
H9,10,41,46
6 554312
7856446
1
1
1
H39
7 551789
7877140
4
4
1
H2,3,11,16
8 559229
7870722
6
5
1
H7,8,13,14,15
9 562251
7877894
1
1
1
H12
10 522661
7878691
4
2
1
H5,38
11 569464
7899894
11
8
1,6
H22,23,24,25,26,28,30 H48
12 535526
7887775
4
4
1
H17,31,32,33
13 563591
7905440
4
3
1,5,6
H28 H1 H47
14 546328
7910163
1
1
1
H29
15 540239
7940493
1
1
1
H4
16 525726
7986906
3
2
1
H20,28
17 515194
7977096
3
2
1
H18,19
18 489648
7962718
3
2
1
H27,45
19 617279
7858220
1
1
2
H6
20
750849
7810197
1
1
3
H37
21*
700813
8514271
6
6
-
-
-
-
-
-
-
46
-
5.8S
-
ITS2
Sampling
site UTM1 UTM2 Ncr
NH
Lin.
Pres.
Lin. 1 Lin. 2
Lin.
3
Lin. 4
1 470914
7831039
12
1
2
H20
2 585530
7848179
8
2
1
H10, 12
3
590033
7847313
4
1
1
H12
4 535291
7849567
10
2
1
H6,8
5 539341
7857588
10
5
1
H6,7,9,10,11
6 554312
7856446
2
1
1
H10
7 551789
7877140
8
1
1
H11
8 559229
7870722
14
5
1
H10,11,14,17,18
9 562251
7877894
-
-
-
10 522661
7878691
8
1
1
H19
11 569464
7899894
26
4
1,4
H11,15 H21,22
12 535526
7887775
8
2
1
H11,19
13 563591
7905440
8
4
1,3,4
H11,13 H22
H1
14 546328
7910163
-
-
-
15 540239
7940493
4
1
1
H11
16 525726
7986906
6
3
1
H3,4,5
17 515194
7977096
6
1
1
H3
18 489648
7962718
6
1
1
H2
19 617279
7858220
2
1
1
H12
20
750849
7810197
2
1
1
H16
21*
700813
8514271
10
1
-
-
-
-
-
47
TABLE 2
Summary details of the number of sequences (Nseq), haplotypes (Nhap), segregating sites (S), and values of haplotype (h) and
nucleotide (π) diversity, and average number of nucleotide differences (k). At the top (a), data of 6 COI lineages and bottom (b)
data of the 4 ITS1-5.8S-ITS2 lineages. Standard deviation is in parentheses.
a) COI
Total R. alatus Total R. motucu
Lin. 1 Lin. 2 Lin. 3 Lin. 4 Lin. 5 Lin. 6
Nseq 69 6 55 4 1 6 1 2
Nhap 53 6 40 4 1 5 1 2
S 189 42 67 47 - 12 - 2
h 0.988(±0.005) 1.000(±0.092) 0.983(±0.008) 1.000(±0.176) - 0.933(±0.121) - 1.000(±0.500)

0.064(±0.031) 0.026(±0.016) 0.025(±0.012) 0.040(±0.027) - 0.007(±0.004) - 0.003(±0.004)
k 38.403(±16.894)
15.933(±8.309) 14.923(±6.777)
24.00(±13.468)
- 4.200(±2.424) - 2.0001.732)
b) ITS1-5.8S-ITS2
Total R. alatus Total R. motucu
Lin. 1 Lin. 2 Lin. 3 Lin. 4
Nseq 144 10 126 12 4 2
Nhap 19 1 18 1 2 1
S 84 - 20 - 1 -
h 0.854(±0.023) - 0.818(±0.029) - 0.5(±0.265) -

0.019(±0.009) - 0.003(±0.002) - 0.0005(±0.0006)
-
k 17.779(±7.939) - 3.470(±1.782) - 0.5(±0.519) -
48
TABLE 3
Matrixes of pairwise differences between lineages, according to model of Tamura
and Nei. At the top (a), distance between the 6 COI lineages and bottom (b) distance
between of 4 ITS1-5.8S-ITS2. The significance level adopted was P < 0.01. NS: not
significant.
a) COI
Lin. 1 Lin. 2 Lin. 3 Lin. 4 Lin. 5 Lin. 6
Lin. 1
-
Lin. 2 0.693
-
Lin. 3 0.738(NS)
0.586(NS)
-
Lin. 4 0.856 0.876 0.950(NS)
-
Lin. 5 0.846(NS)
0.765(NS)
1.000(NS)
0.953(NS)
-
Lin. 6 0.847 0.803(NS)
0.978(NS)
0.953(NS)
0.976(NS)
-
b) ITS1
-
5.8S
-
ITS2
Lin. 1 Lin. 2 Lin. 3 Lin. 4
Lin. 1
-
Lin. 2 0.917
-
Lin. 3 0.952 0.997
-
Lin. 4 0.956 1.000
0.994
(NS)
-
49
TABLE 4
Simulations of SAMOVA. Respectively, values of K, FSC (variation among
populations within groups), FST (variation within populations), and FCT (variation
between groups), are showed. The significance level adopted was P < 0.05 and are
in parentheses. At the top (a), values of COI and bottom (b) values of ITS1-5.8S-
ITS2.
a) COI
K FSC (P) FST (P) FCT (P)
2 0.497(0.000) 0.856(0.000) 0.715(0.041)
3 0.477(0.000) 0.848(0.000) 0.709(0.001)
4 0.458(0.000) 0.840(0.000) 0.705(0.000)
5 0.435(0.000) 0.833(0.000) 0.705(0.000)
6 0.436(0.000) 0.821(0.000) 0.683(0.000)
7 0.440(0.000) 0.809(0.000) 0.659(0.000)
8 0.096(0.000) 0.674(0.000) 0.640(0.000)
9 0.040(0.000) 0.675(0.000) 0.661(0.000)
b) ITS1-5.8S-ITS2
K FSC (P) FST (P) FCT (P)
2 0.372(0.000) 0.887(0.000) 0.821(0.030)
3 0.370(0.000) 0.873(0.000) 0.798(0.001)
4 0.269(0.000) 0.84(0.000) 0.781(0.016)
5 0.267(0.000) 0.829(0.000) 0.767(0.000)
6 0.271(0.000) 0.819(0.000) 0.752(0.000)
7 0.027(0.003) 0.742(0.000) 0.735(0.000)
8 0.025(0.000) 0.74(0.000) 0.733(0.000)
9 0.028(0.000) 0.739(0.000) 0.731(0.000)
50
Figure 1
51
Figure 2
52
Figure 3
53
Figura 4
54
CONCLUSÕES
As populações do minhocuçu R. alatus revelaram uma alta diversidade genética,
sendo possível a identificação de linhagens distintas para ambos marcadores.
Houve concordância de indivíduos para ambos marcadores em três linhagens, as
quais foram menos freqüentes. Por outro lado, a linhagem mais comum para o segmento
ITS1-5.8S-ITS2 coincidiu com três linhagens diferentes de COI.
A distribuição espacial da diversidade genética da espécie no estado de Minas
Gerais se mostrou moderada, o que pode ser justificada pela baixa capacidade de
dispersão, comum entre minhocas.
O presente estudo elucidou características relevantes para o conhecimento da
história de vida e para a conservação de uma espécie com grande importância local.
Pesquisas básicas em qualquer âmbito da biologia da conservação de R. alatus são úteis
para o delineamento de medidas de manejo adequado. Desse modo, entende-se que esse
seja apenas o primeiro passo como contribuição para a preservão e uso sustentável do
minhocuçu.
55
REFERÊNCIAS BIBLIOGRÁFICAS
Allendorf FW, Luikart G (2007) Conservation and the Genetics of Populations. Blackwell
Publishing, Malden.
Avise JC (2000) Phylogeography: The History and Formation of Species. Harvard University
Press, Cambridge, Massachusetts.
Avise JC (2004) Molecular Markers, Natural History, and Evolution. 2
nd
ed. Sinauer
Associates, Inc. Publishers, Sunderland, Massachusetts.
BARCODE OF LIFE DATA SYSTEM. Disponível em: www.barcodinglife.org. Acesso em: 18
de maio de 2009.
Bargues MD, Klisiowicz DR, Panzera F, Noireau F, Marcilla A, Perez R, Rojas MG, O’connor
JE, Gonzalez CF, Galvão C, Jurberg J, Carcavallo RU, Dujardin JP, Mascoma S (2006)
Origin and phylogeography of the Chagas disease main vector Triatoma infestans based on
nuclear rDNA sequences and genome size. Infections, Genetics and Evolution, 6, 46-62.
Bastrop R, Blank M (2006) Multiple invasions a polychaete genus enter the Baltic sea.
Biological invasions, 8, 1195-1200.
Cameron E K, Bayne EM, Coltman DW (2008) Genetic structure of invasive earthworms
Dendrobaena octaedra in the boreal forest of Alberta: insights into introduction mechanisms.
Molecular Ecology, 17, 1189-1197.
Chen CA, Chen CP, Fan TY, Yu JK, Hsieh HL (2002) Nucleotide sequences of ribosomal
internal transcribed spacer and their utility in distinguishing closely related Perinereis
polychaets (Annelida; Polychaeta; Nereididae). Mar. Biotechnol., 4, 17-29.
Chang C, Lin S, Chen J. (2008) Molecular systematics and phylogeography of the gigantic
earthworms of the Metaphire formosae species group (Clitellata, Megascolecidae). Molecular
Phylogenetics and Evolution, 49, 958-968.
Christoffersen ML (2007) Distribution and species diversity of Rhinodrilus Perrier, 1872
(Annelida, Clitellata, Lumbricina, Glossoscolecidae) in South America. Neodiversity, 2, 1-6.
Darwin C (1881) The Formation of Vegetable Mould Through the Actions of Worms. John
Murray, London.
Drumond MA (2008) Manejo adaptativo do minhocuçu Rhinodrilus alatus. Tese de
Doutorado. Universidade Federal de Minas Gerais.
Erseus C (2005) Phylogeny of oligochaetous Clitellata. Hydrobiologia, 535/536, 357-372
Field SG, Lange M, Schulenburg H, Velavan T P, Michiels NK. (2007) Genetic diversity and
parasite defense in a fragmented urban metapopulation of earthworms. Animal
Conservation, 10, 162-175.
56
Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R (1994) DNA primers for the amplification
of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates.
Molecular Marine Biology and Biotechnology, 3, 294–299.
Frankham R, Ballou JD, Briscoe DA (2002) Introduction to conservation genetics. Cambridge
University Press, Cambridge.
Haimi J, Knott KE, Yliniemi I (2007) Does metal contamination affect clonal diversity of the
parthenogenetic earthworm Dendrobaena octaedra? European Journal of Soil Biology, 43,
268-272.
Harper GL, Cesarini S, Casey SP, Morgan AJ, Kille P, Bruford MW (2006) Microsatellite
markers for the earthworm Lumbricus rubellus. Molecular Ecology Notes, 6, 325-327.
Hartl DL, Clark AG (2007) Principles of population genetics. Sinauer Associations,
Massachusetts.
Hebert PDN, Cywinska A, Ball SL, DeWaard JR (2003) Biological identifications through
DNA barcodes. Proceedings of the Royal Society B: Biological Sciences, 270, 313–321.
Heethoff M, Etzold K, Scheu S (2004) Mitochondrial COII sequences indicate that the
parthenogenetic earthworm Octolasion tyrtaeum (Savigny 1826) constitutes of two lineages
differing in body size and genotype. Pedobiologia, 48, 9–13.
Huang J, Xu Q, Sun ZJ, Tang GL, Su ZY (2007) Identifying earthworms through DNA
barcodes. Pedobiologia, 51, 301-309.
James SW (2004) Planetary process and their interactions with earthworm distribution and
ecology. In: Earthworm Ecology (ed. Edwards CA), pp. 53–62. CRC Press, Boca Raton,
Florida.
Jolly MT, Jollivet D, Gentil F, Thiébaut E, Viard F (2005) Sharp genetic break between
Atlantic and English Channel populations of the polychaete Pectnaria koreni, along the North
coast of France. Heredity, 94, 23-32.
Jolly MT, Viard F, Gentil F, Thiébaut E, Jollivet D. (2006) Comparative phylogeography of
two coastal polychaete tubeworms in the Northeast Atlantic supports shared history and
vicariant events. Molecular Ecology, 15, 1841-1855.
Kautenburger R (2006) Genetic structure among earthworms (Lumbricus terrestris L.) from
different sampling sites in western Germany based on random amplified polymorphic DNA.
Pedobiologia, 50, 257-266.
King RA, Tibble AL, Symondson WOC (2008) Open a can or worms: unprecedented
sympatric cryptic diversity within British lumbricid earthworms. Molecular Ecology, 17, 4684-
4698.
Lee T, Foighil DO (2003) Phylogenetic structure of the Sphaeriinae, a global clade of
freshwater bivalves mollusks, inferred from nuclear (ITS-1) and mitochondrial (16S)
ribosomal gene sequences. Zoological Journal of the Linnean Society, 137, 245-260.
57
Lentzsch P, Golldack J (2006) Genetic diversity of Aporrectodea caliginosa from agricultural
sites in Northeast Brandenburg, Germany. Pedobiologia, 50, 369–376.
Lewin B (2004) Genes VIII. Pearson Prentice Hall, Upper Saddle River.
Machado ABM, Drummond GM, Paglia AP (2008) Livro Vermelho da Fauna Brasileira
Ameaçada de Extinção. Fundação Biodiversitas, Belo Horizonte, Brasil.
McHugh D (2005) Molecular Systematics of polychaetes (Annelida). Hydrobiologia, 535/536,
309-318.
Meyer A, Bleidorn C, Rouse GW, Hausen H (2008) Morphological and molecular data
suggest a cosmopolitan distribution of the polychaete Proscoloplos cygnochaetus, Day, 1954
(Annelida, Orbiniidae). Mar. Biol., 153, 789-889.
Miranda AL (1987) Era uma vez os rios. Belo Horizonte, Minas Gerais, Brasil.
NATIONAL CENTER FOR BIOTECHNOLOGY INFORMATION. Disponível em:
<www.ncbi.nlm.nih.gov >. Acesso em: 08 de junho 2009.
O’brien SJ (1994) A role for molecular genetics in biological conservation. Proc. Natl. Acad.
Sci. USA, 91, 5748-5755.
Page RDM, Holmes EC (1998) Molecular Evolution: A Phylogenetic Approach. Blackwell
Science Inc., Oxford.
Righi G (1971) Sobre a família Glossoscolecidae (Oligochaeta) no Brasil. Arquivos de
Zoologia, 65, 1-95.
Römbke J, Jänsch S, Didden W (2005) The use of earthworms in ecological soil
classification and assessment concepts. Ecotoxicology and Environmental Safety, 62, 249-
265.
Saiki RK, Scharff S, Faloona F, Mullis KB, Horn GT, Erlich HA, Arnheim N (1985) Enzimatic
amplification of beta-globin genomic sequences and restriction site analysis for diagnosis of
sickle cell anemia. Science, 230, 1350-1354.
Schlötterer C (2004) The evolution of molecular markers just a matter of fashion? Nature
Review Genetics, 5, 63-69.
Simon C, Frati F, Beckenbach A, Crespi B, Liu H, Flook P (1994) Evolution, Weighting, and
Phylogentic utility of mitochondrial gene sequences and a compilation of conserved
Polymerase Chain Reaction primers. Annals of the Entomological Society of America, 87,
651-700.
Struck TH, Schult N, Kusen T, Hickman E, Bleidorn C, McHugh D, Halanyck KM (2007)
Annelid Phylogeny and the status of Sipuncula and Echiura. BMC Evolutionary Biology, 7,
57.
Symondson WOC, Glend M, Erickson ML, Liddell JE, Langdon CJ (2000) Do earthworms
help to sustain the slug predator Pterostichus melanarius (Coleoptera: Carabidae) within
crops? Investigations using monoclonal antibodies. Molecular Ecology, 9, 1279-1292.
58
Templeton AR (2006) Population genetics and microevolutionary theory. John Wiley & Sons,
New Jersey.
Velavan TP, Schulenburg H, Michiels NK (2007) Development and characterization of novel
microsatellite markers for the common earthworm (Lumbricus terrestris L.). Molecular
Ecology Notes, 7, 160-162.
Vogler AP, DeSalle R (1994) Evolution and phylogenetic information content of the ITS-1
region in the Tiger Beetle Cicindela dorsalis. Molecular Biology Evolution, 11, 393-405.
Livros Grátis
( http://www.livrosgratis.com.br )
Milhares de Livros para Download:
Baixar livros de Administração
Baixar livros de Agronomia
Baixar livros de Arquitetura
Baixar livros de Artes
Baixar livros de Astronomia
Baixar livros de Biologia Geral
Baixar livros de Ciência da Computação
Baixar livros de Ciência da Informação
Baixar livros de Ciência Política
Baixar livros de Ciências da Saúde
Baixar livros de Comunicação
Baixar livros do Conselho Nacional de Educação - CNE
Baixar livros de Defesa civil
Baixar livros de Direito
Baixar livros de Direitos humanos
Baixar livros de Economia
Baixar livros de Economia Doméstica
Baixar livros de Educação
Baixar livros de Educação - Trânsito
Baixar livros de Educação Física
Baixar livros de Engenharia Aeroespacial
Baixar livros de Farmácia
Baixar livros de Filosofia
Baixar livros de Física
Baixar livros de Geociências
Baixar livros de Geografia
Baixar livros de História
Baixar livros de Línguas
Baixar livros de Literatura
Baixar livros de Literatura de Cordel
Baixar livros de Literatura Infantil
Baixar livros de Matemática
Baixar livros de Medicina
Baixar livros de Medicina Veterinária
Baixar livros de Meio Ambiente
Baixar livros de Meteorologia
Baixar Monografias e TCC
Baixar livros Multidisciplinar
Baixar livros de Música
Baixar livros de Psicologia
Baixar livros de Química
Baixar livros de Saúde Coletiva
Baixar livros de Serviço Social
Baixar livros de Sociologia
Baixar livros de Teologia
Baixar livros de Trabalho
Baixar livros de Turismo