|
Memórias do Instituto Oswaldo Cruz, Vol. 102, No. 7, November, 2007, pp. 879-882 Short Communication A preliminary analysis of the population genetics and molecular phylogenetics of Onchocerca volvulus (Nematoda: Filarioidea) using nuclear ribosomal second internal transcribed spacer sequences Ramiro Morales-Hojas/*/+, Robert A Cheke*, Rory J Post** Laboratório
de Evolução Molecular, Instituto de Biologia Molecular
e Celular, Universidade do Porto, Rua do Campo Alegre 823, 4150-180
Porto, Portugal *Natural Resources Institute, University of Greenwich
at Medway, Kent, UK **Department of Entomology, The Natural History
Museum, London, UK and Department of Infectious & Tropical Diseases,
London School of Hygiene & Tropical Medicine, London, UK Received
24 April 2007 Code Number: oc07158 Nuclear internal transcribed spacer 2 (ITS2) rDNA sequences were used for a molecular phylogenetics analysis of five Onchocerca species. The sister species of the human parasite O. volvulus was found to be the cattle parasite O. ochengi and not O. gibsoni, contrary to chromosomal evidence. The genetic differentiation of two African populations (representing the two African strains) and a Brazilian population of O. volvulus was also studied. Phylogenetic and network reconstruction did not show any clustering of ITS2 alleles on geographic or strain grounds. Furthermore, population genetics tests showed no indication of population differentiation but suggested gene flow among the three populations. Key words: Onchocerca volvulus - onchocerciasis - strains - molecular phylogenetics - population differentiation The parasite Onchocerca volvulus (Nematoda: Filarioidea) is the causative agent of onchocerciasis (river blindness), a major public health problem in tropical Africa and in isolated foci in Yemen, South and Central America. There is evidence for different pathogenicity between "strains" from the savanna and the forest of West Africa (Duke 1981) with reported levels of blindness among infected people in the savanna being higher than in the forest (up to a maxima of 15% compared to 2%, respectively) and in the forested areas severe ocular complications such as sclerosing keratitis are mostly absent, with dermatological symptoms predominant. There are also important differences in the competence of different vector species to support the different "strains" (Duke 1981), and geographic variation in morphology (Eichner & Renz 1990), isoenzymes (Cianchi et al. 1985), endosymbionts (Higazi et al. 2005), and O-150 repetitive DNA (Zimmerman et al. 1992) has been interpreted as correlating with forest/savanna "strain" designation. Furthermore, the O-150 repetitive DNA suggested that the New World populations of O. volvulus were derived from African savanna populations (Zimmerman et al. 1994). However, the morphological and isoenzymatic data were never followed-up, and subsequent data on the O-150 repetitive DNA has not given such a clear correlation in Africa (Fischer et al. 1996). The nuclear ITS rDNA sequence has been extensively used to define genetic markers for different species of nematodes (Powers et al. 1997), and the aims of this study were to use the ITS2 to investigate the levels of genetic variation within O. volvulus and determine whether it is related to strain differentiation. This is information of importance for understanding the epidemiology of onchocerciasis and the design of control programmes. Thus, if there are normal levels of gene flow between populations (or strains), these populations are likely to be similar in their overall genetic characteristics (and response to control), genes for important characteristics (such as filaricide resistance) are likely to spread easily, and disease control at single locations or targeted at a single strain will be countered by immigration. ITS2 rDNA sequences were already available from Brazilian and African forest samples of O. volvulus (Morales-Hojas et al. 2001, GenBank accession numbers AF228565-AF228570 and AF228572-AF228576). For samples of savanna O. volvulus (Kati district, Mali), O. ochengi (Bamako, Mali), O. lienalis (Shrewsbury, UK), O. gutturosa (Bamako, Mali) and O. gibsoni (Australia), amplification of the ITS2 was carried out using primers 5.8SOv-A 5-AGCGGTGGATCACTTGGCTC-3 (located in position 9-29 of the 5.8S gene) and NC2 5-TTAGT-TTCTTTTCCTCCGCT-3 (located 70 bp downstream of the 5 end of the 28s rRNA gene). The origin of this material is described in Morales-Hojas et al. (2006). Reactions were performed in 1x Buffer (Promega), 3 mM MgCl2 (Promega), 200 mM of each dNTP (Promega), 0.2 mM of each primer, 1 unit of Taq DNA Polymerase (Promega) and 1.5 ml of the DNA extraction. Amplifications consisted of a first denaturation step at 94oC for 2 min, followed by 30 cycles of 30 s at 94oC, 30 s at 50oC (annealing) and 30 s at 72oC, with a final extension step of 5 min at 72oC. Amplification products were extracted directly from the gel using the Geneclean II kit (Anachem) and cloned using the TOPO TA Cloning Kit for Sequencing (Invitrogen). Positive clones were extracted from transformed bacteria using the S.N.A.P. MiniPrep kit (Invitrogen) and sent to the Advanced Biotechnology Centre, Imperial College School of Medicine, for sequencing in both directions, to correct for possible reading errors, using Big Dye v3.1 chemistry (Perkin Elmer/Applied Biosystems) and the universal primers T3 and T7. It is possible that the sequences include PCR-induced errors. Nevertheless, these would be difficult to correct for because the DNA used is a pool of genomes of the different individuals and thus amplification products are not expected to be homogeneous. Furthermore, the ribosomal genes are a family of tandemly repeated units that despite concerted evolution do not necessarily show complete homoge-nisation. However, the expected number of errors introduced by non-proofreading polymerases is very small, in the range of 1 x 10-4 to 10-5 errors per bp. The new ITS2 sequences have been deposited in GenBank and have accession numbers DQ317639-DQ317666. Amplification products showed the expected size, ranging between 533 bp and 568 bp. The length of the ITS2 ranged between 318-319 bp in O. lienalis, 333-337 bp in O. gutturosa, 341-351 bp in O. gibsoni, 347-352 bp in O. ochengi, and 340-352 bp in O. volvulus. Phylogenetic results using both Maximum Likelihood and Maximum Parsimony resulted in similar trees (Fig. 1). They supported the mitochondrial analyses that placed O. ochengi and O. volvulus as sister species and not O. gibsoni and O. volvulus, as suggested by karyotype evolution data (Post 2005). It is interesting to note that O. ochengi ITS2 alleles are paraphyletic with respect to O. volvulus, which is consistent with a recent cladogenesis of O. volvulus after the spread of domestic cattle in Africa (5000 to 2500 years ago) as discussed by Morales-Hojas et al. (2006). However, this is inconsistent with mitochondrial DNA analyses which resulted in a tree where O. volvulus and O. ochengi showed reciprocal monophyly (Morales-Hojas et al. 2006). This difference might be explained by the faster rate of coalescence of mitochondrial alleles (four times faster than that of the nuclear alleles; Ballard & Whitlock 2004). The evolutionary relationships between O. volvulus alleles were further explored with statistical parsimony networks (Templeton et al. 1992). The reconstructed network showed allele A (from the forest strain) as the most recent common ancestral allele (Fig. 2), and branches in a star-like shape with no geographical or "strain" clustering of sequences. Consistent with these results, population differentiation estimates Fst, Kst* and Snn [estimated with DnaSP 4.00 (Rozas et al. 2004) and Proseq 2.9 (Filatov 2002)] were all non-significant indicating no genetic differentiation in the ITS2 sequences between the three populations (Table I). Furthermore, the estimated number of migrant individuals per generation based on Fst values was higher than 1, which would be sufficient to prevent population differentiation. The high haplotype diversity (H) observed, combined with the low sequence diversity (p) in the sample of O. volvulus (Table II), is often considered indicative of a population expansion after a reduced ancestral effective population size (Avise 2000). This would be consistent with a host switch from cattle to humans in which only a few individuals were involved, thus resulting in a bottleneck. Following the host switch, an increase in population size in the new (human) host must have occurred to give the current millions of infected people. Although a small number of errors may have been introduced in the sequences due to non-proofreading polymerase mistakes (around 1x10-4 to 10-5 per bp), the main consequence of these would be to increase the genetic and haplotype diversities. This could affect the interpretation that after a bottleneck the species population has increased its effective size. Nevertheless, despite this possibility of an artificial increase of the diversity the estimated levels of genetic diversity remain low suggesting that the PCR-introduced mutations are few, and not affecting the historical demographic interpretation. Furthermore, this is otherwise supported by the evolutionary history of the parasite, with a recent host-switch and an increase in population size alongside its human host. Other potential artifact of these PCR mistakes could be to increase the number of nucleotide differences between populations. Despite this, still non-significant levels of genetic differentiation between populations were obtained with this marker. Thus, the possible PCR-induced errors introduced in the sequences appear to have little consequences in the interpretation of the results. In conclusion, the rDNA ITS2 shows no intraspecific structure in O. volvulus and the suggested gene flow between the African populations is consistent with what would be expected as a result of the known annual migrations of the West African vector species (Boakye et al. 1998). However, further samples and additional (more sensitive) nuclear and mitochondrial genes should be investigated in order to corroborate the results presented in this study. This will be of great importance for control strategies. Poor response to the microfilaricidal drug Ivermectin has been reported in some populations of O. volvulus (e.g. Awadzi et al. 2004). Results presented here could suggest that resistance genes would spread freely among strains and populations in Africa if resistance to Ivermectin resulted in a re-emergence of parasite transmission at significant levels. Thus, strategies for onchocerciasis control should take into account the fact that populations may not be isolated and that "strain" delimitation is not necessarily indicative of genetic isolation and generalised genetic differentiation. REFERENCES
Copyright 2007 Instituto Oswaldo Cruz - Fiocruz The following images related to this document are available:Photo images[oc07158f1.jpg] [oc07158t1.jpg] [oc07158t2.jpg] [oc07158f2.jpg] |
|