Publications

Benzimidazole resistance has evolved in a variety of organisms and typically results from mutations in the beta-tubulin locus at specific amino acid sites. Despite widespread treatment of human intestinal nematodes with benzimidazole drugs, there have been no unambiguous reports of resistance. However, since beta-tubulin mutations conferring resistance are generally recessive, frequencies of resistance alleles less than 30% would be difficult to detect on the basis of drug treatment failures. Here we investigate sequence variation in a 1079 bp segment of the beta-tubulin locus in the human whipworm Trichuris trichiura from 72 individual nematodes from seven countries. We did not observe any alleles with amino acid mutations indicative of resistance, and of 40 point mutations there were only four non-synonymous mutations all of which were singletons. Estimated effective population sizes are an order of magnitude lower than those from another nematode species in which benzimidazole resistance has developed (Haemonchus contortus). Both the lower diversity and reduced population sizes suggest that benzimidazole resistance is likely to evolve less rapidly in Trichuris than in trichostrongyle parasites of livestock. We observed moderate levels of population subdivision (Phi(ST)=0.26) comparable with that previously observed in Ascaris lumbricoides, and identical alleles were frequently found in parasites from different continents, suggestive of recent admixture. A particularly interesting feature of the data is the high nucleotide diversities observed in nematodes from the Caribbean. This genetic complexity may be a direct result of extensive admixture and complex history of human populations in this region of the world. These data should encourage (but not make complacent) those involved in large-scale benzimidazole treatment of human intestinal nematodes.

Molecular markers are used widely to discriminate between closely related species of parasites, and in many cases a single locus is used for this purpose. This article aims to show how molecular data derived from a single genetic marker or linkage group - in this case mitochondrial DNA - can lead to ambiguous conclusions and to illustrate how a multilocus approach has enhanced our understanding of the epidemiology of two closely related parasites, the nematodes Ascaris suum, which infects pigs, and Ascaris lumbicoides, which infects humans.

Infection with the endosymbiotic bacteria Wolbachia is widespread in filarial nematodes. Previous studies have suggested concordance between the phylogeny of Wolbachia with that of their nematode hosts. However, there is only one published molecular phylogenetic study of filarial species, based on the 5S rRNA gene spacer. The phylogeny proposed by this study is partially incongruent with previous classifications of filarial nematodes, based on morphological characters. Furthermore, both traditional classifications and molecular phylogenies are, in part, inconsistent with the phylogeny of Wolbachia. Here we report mitochondrial cytochrome oxidase I (COI) gene sequences for 11 species of filaria and for another spirurid nematode which was included as an outgroup. In addition, 16S rRNA, wsp and ftsZ gene sequences were generated for the Wolbachia of several filarial species, in order to complete the available data sets and further resolve the phylogeny of Wolbachia in nematodes. We used these data to evaluate whether nematode and Wolbachia phylogenies are concordant. Some of the possible phylogenetic reconstructions based on COI gene were congruent with the phylogeny of Wolbachia and supported the grouping of the rodent filaria Litomosoides sigmodontis with the lymphatic filariae (i.e. Brugia spp. and Wuchereria spp.) and the sister group relationship of Dirofilaria spp. and Onchocerca spp. However, the placement of the Wolbachia-free filaria Acanthocheilonema viteae is ambiguous and dependent on the phylogenetic methods used.

Microsatellite loci are generally assumed to evolve via a stepwise mutational process and a battery of statistical techniques has been developed in recent years based on this or related mutation models. It is therefore important to investigate the appropriateness of these models in a wide variety of taxa. We used two approaches to examine mutation patterns in the malaria parasite Plasmodium falciparum: (i) we examined sequence variation at 12 tri-nucleotide repeat loci; and (ii) we analysed patterns of repeat structure and heterozygosity at 114 loci using data from 12 laboratory parasite lines. The sequencing study revealed complex patterns of mutation in five of the 12 loci studied. Alleles at two loci contain indels of 24 bp and 57 bp in flanking regions, while in the other three loci, blocks of imperfect microsatellites appear to be duplicated or inserted; these loci essentially consist of minisatellite repeats, with each repeat unit containing four to eight microsatellites. The survey of heterozygosity revealed a positive relationship between repeat number and microsatellite variability for both di- and trinucleotides, indicating a higher mutation rate in loci with longer repeat arrays. Comparisons of levels of variation in different repeat types indicate that the mutation rate of dinucleotide-bearing loci is 1.6-2.1 times faster than trinucleotides, consistent with the lower mean number of repeats in trinucleotide-bearing loci. However, despite the evidence that microsatellite arrays themselves are evolving in a manner consistent with stepwise mutation model in P. falciparum, the high frequency of complex mutations precludes the use of analytical tools based on this mutation model for many microsatellite-bearing loci in this protozoan. The results call into question the generality of models based on stepwise mutation for analysing microsatellite data, but also demonstrate the ease with which loci that violate model assumptions can be detected using minimal sequencing effort.

Wolbachia endosymbiotic bacteria are widespread in arthropods and are also present in filarial nematodes. Almost all filarial species so far examined have been found to harbor these endosymbionts. The sequences of only three genes have been published for nematode Wolbachia (i.e., the genes coding for the proteins FtsZ and catalase and for 16S rRNA). Here we present the sequences of the genes coding for the Wolbachia surface protein (WSP) from the endosymbionts of eight species of filaria. Complete gene sequences were obtained from the endosymbionts of two different species, Dirofilaria immitis and Brugia malayi. These sequences allowed us to design general primers for amplification of the wsp gene from the Wolbachia of all filarial species examined. For these species, partial WSP sequences (about 600 base pairs) were obtained with these primers. Phylogenetic analysis groups these nematode wsp sequences into a coherent cluster. Within the nematode cluster, wsp-based Wolbachia phylogeny matches a previous phylogeny obtained with ftsZ gene sequences, with a good consistency of the phylogeny of hosts (nematodes) and symbionts (Wolbachia). In addition, different individuals of the same host species (Dirofilaria immitis and Wuchereria bancrofti) show identical wsp gene sequences.

Polymerase chain reaction (PCR)-based genotyping of oocysts dissected from mosquito midguts has previously been used to investigate overall levels of inbreeding within malaria parasite populations. We present a re-analysis of the population structure of Plasmodium falciparum malaria using diploid genotypes at three antigen-encoding loci in 118 oocysts dissected from 34 mosquitoes. We use these data to ask whether mating is occurring at random within the mosquito midgut, as is generally assumed. We observe a highly significant deficit of heterozygous oocysts within mosquitoes at all three loci, suggesting that fusion of gametes occurs non-randomly in the mosquito gut. A variety of biological explanations, such as interrupted feeding of mosquitoes, positive assortative mating and outcrossing depression, could account for this observation. However, an alternative artefactual explanation–the presence of non-amplifying or null alleles–can account for the observed data equally well, without the need to invoke non-random mating. To evaluate this explanation further, we estimate the frequencies of null alleles within the oocyst population using maximum likelihood, by making the assumption that non-amplifying oocysts at any of the three loci are homozygous for null alleles. Observed levels of visible heterozygotes fit closely with those expected under random mating when non-amplifying oocysts are accounted for. Other lines of evidence also support the artefactual explanation. Overall inbreeding coefficients have been recalculated in the light of this analysis, and may be considerably lower than those estimated previously. In conclusion, we suggest that the deficit of heterozygotes observed is unlikely to indicate non-random mating within the mosquito gut and is better explained by misscoring of heterozygotes as homozygotes.

Parasite genotyping by the polymerase chain reaction was used to distinguish recrudescent from newly acquired Plasmodium falciparum infections in a Karen population resident on the northwestern border of Thailand where malaria transmission is low (one infection/person/year). Plasmodium falciparum infections were genotyped for allelic variation in three polymorphic antigen loci, merozoite surface proteins-1 and -2 (MSP-1 and -2) and glutamaterich protein (GLURP), before and after antimalarial drug treatment. Population genotype frequencies were measured to provide the baseline information to calculate the probability of a new infection with a different or the same genotype to the initial pretreatment isolate. Overall, 38% of the infections detected following treatment had an identical genotype before and up to 121 days after treatment. These post-treatment genotypes were considered recrudescent because of the low (< 5%) probability of repeated occurrence by chance in the same patient. This approach allows studies of antimalarial drug treatment to be conducted in areas of low transmission since recrudescences can be distinguished confidently from newly acquired infections.

Multiple, selectively neutral genetic markers are the most appropriate tools for analysis of parasite population structure and epidemiology, but yet existing methods for characterization of malaria field samples utilize a limited number of antigen encoding genes, which appear to be under strong selection. We describe protocols for characterization of 12 microsatellite markers from finger-prick blood samples infected with Plasmodium falciparum. A two-step, heminested strategy was used to amplify all loci, and products were visualized by fluorescent end-labelling of internal primers. This procedure allows amplification from low levels of template, while eliminating the problem of spurious products due to primer carry over from the primary round of PCR. The loci can be conveniently multiplexed, while accurate sizing and quantification of PCR products can be automated using the GENOTYPER software. The primers do not amplify co-infecting malaria species such as P. vivax and P. malariae. To demonstrate the utility of these markers, we characterized 57 infected finger-prick blood samples from the village of Mebat in Papua New Guinea for all 12 loci, and all samples were genotyped a second time to measure reproducibility. Numbers of alleles per locus range from 4 to 10 in this population, while heterozygosities range from 0.21 to 0.87. Reproducibility (measured as concordance between predominant alleles detected in replicate samples) ranged from 92 to 98% for the 12 loci. The composition of PCR products from infections containing multiple malaria clones could also be defined using strict criteria and scored in a highly repeatable manner.

It has recently been shown using genetic markers that Ascaris in humans and pigs in Central America comprise reproductively isolated populations. We present a similar analysis for a region of China in which close association between pigs and humans has been the norm for thousands of years, and agricultural practices will result in frequent exposure to eggs from both sources. DNA fragments from selected regions of mitochondrial and ribosomal DNA were amplified by PCR and allelic forms identified following digestion with a panel of restriction enzymes, using DNA from a total of 115 individual worms from both people and pigs from 2 neighbouring villages. Significant frequency differences in both mtDNA haplotypes and the rDNA spacer were found between the 2 host-associated populations, indicating that they represented reproductively isolated populations. Mitochondrial haplotype frequencies were different from those observed in Guatemala and also from other Asian Ascaris populations, suggesting low levels of gene flow between populations. However, we found no evidence for significant heterogeneity in the genetic composition of Ascaris infrapopulations in either humans or pigs, possibly indicative of agricultural practices in China which have resulted in a random distribution of alleles within the parasite populations.