Genetic diversity of a population of Schistosoma haematobium derived from schoolchildren in east central Zimbabwe

Dcifer: an IBD-based method to calculate genetic distance between polyclonal infections

An essential step toward reconstructing pathogen transmission and answering epidemiologically relevant questions from genomic data is obtaining pairwise genetic distance between infections. For recombining organisms such as malaria parasites, relatedness measures quantifying recent shared ancestry would provide a meaningful distance, suggesting methods based on identity by descent (IBD). While the concept of relatedness and consequently an IBD approach is fairly straightforward for individual parasites, the distance between polyclonal infections, which are prevalent in malaria, presents specific challenges, and awaits a general solution that could be applied to infections of any clonality and accommodate multiallelic (e.g. microsatellite or microhaplotype) and biallelic [single nucleotide polymorphism (SNP)] data. Filling this methodological gap, we present Dcifer (Distance for complex infections: fast estimation of relatedness), a method for calculating genetic distance between polyclonal infections, which is designed for unphased data, explicitly accounts for population allele frequencies and complexity of infection, and provides reliable inference. Dcifer’s IBD-based framework allows us to define model parameters that represent interhost relatedness and to propose corresponding estimators with attractive statistical properties. By using combinatorics to account for unobserved phased haplotypes, Dcifer is able to quickly process large datasets and estimate pairwise relatedness along with measures of uncertainty. We show that Dcifer delivers accurate and interpretable results and detects related infections with statistical power that is 2–4 times greater than that of approaches based on identity by state. Applications to real data indicate that relatedness structure aligns with geographic locations. Dcifer is implemented in a comprehensive publicly available software package.

Genetic Diversity of Schistosoma haematobium in Qena Governorate, Upper Egypt

Introduction

Schistosomiasis is an important neglected tropical disease (NTD) in several developing countries. Praziquantel is the principle and efficacious chemotherapeutic agent that has been used to treat schistosomiasis for decades. Unfortunately, emerging resistance to praziquantel with accompanying reduced efficacy is reported in some localities. Hence, genetic diversity among parasite populations is of significant interest in assessing the effects of selective pressure generated by praziquantel therapy that might result in encouraging the emergence of new genotypes that are either non-susceptible or drug-resistant. The present study aimed to investigate the genetic diversity of Schistosoma haematobium among human populations using the RAPD technique to help clarify disease epidemiology and transmission.

Materials and Methods

S. haematobium eggs were isolated from 50 of 134 patients from four different localities in Qena Governorate, Upper Egypt. These patients complained of terminal hematuria and burning micturition. Samples were used for molecular analysis using RAPD-PCR primers (A02, A07, A09, A10).

Results

Twenty S. haematobium isolates (40%) were amplified using the selected RAPD primers. Amplification patterns of these isolates showed distinct variation in the size and number of amplified fragments, indicating high genetic variation among these isolates.

Conclusion

To the best of our knowledge, this study is the first to characterize the genetic diversity of S. haematobium in human populations in Upper Egypt. Future studies on a larger geographic scale involving many districts in Upper Egypt should be encouraged. Information from such a study would provide better insight into clonal lineages of S. haematobium in this endemic area. In turn, understanding transmission of the parasite may have a major role in establishing control strategies for urogenital schistosomiasis in Upper Egypt.

Genetic diversity among natural populations of Schistosoma haematobium might contribute to inconsistent virulence and diverse clinical outcomes

There is an evident difference in the intensity of morbidity caused by Schistosoma haematobium in North-African zones compared to Sub-Saharan ones. Clinical outcome dichotomy corresponds to two geographically distinct intermediate host snail species that are only infected by the related strain of the parasite. In concert, there is a manifest hybridization of the parasite with other Schistosoma species confined to certain regions of Africa. This raises a reasonable suggestion that S. haematobium has no less than two phylogenetic clusters that have different virulence. The aim of the study was to examine the possible diversity among S. haematobium using simultaneous amplification of genomic DNA of selected isolates. Random amplified polymorphic DNA-polymerase chain reaction markers were used to study the genetic diversity among S. haematobium natural isolates from selected regions of Africa (Egypt, Zimbabwe, and South Africa) that represent different ecological conditions, different species of intermediate host, and different possibilities of field hybridization with other schistosomes. A moderate to high level of genetic diversity was evident among the three isolates. More bands were shared by the isolates from Zimbabwe and South Africa (similarity index = 0.721) than those shared by each with the Egyptian isolate (similarity index = 0.551 and 0.566, respectively), suggesting that at least two phylogenetic groups of S. haematobium do exist in distinct geographic regions of Africa. The elucidation of the possible genetic diversity among S. haematobium parasites may explain many ambiguous aspects of the biology of the parasite-like virulence, immune evasion and drug resistance.

Genetic Diversity of Schistosoma haematobium Eggs Isolated from Human Urine in Sudan

The genetic diversity of Schistosoma haematobium remains largely unstudied in comparison to that of Schistosoma mansoni. To characterize the extent of genetic diversity in S. haematobium among its definitive host (humans), we collected S. haematobium eggs from the urine of 73 infected schoolchildren at 5 primary schools in White Nile State, Sudan, and then performed a randomly amplified polymorphic DNA marker ITS2 by PCR-RFLP analysis. Among 73 S. haematobium egg-positive cases, 13 were selected based on the presence of the S. haematobium satellite markers A4 and B2 in their genomic DNA, and used for RFLP analysis. The 13 samples were subjected to an RFLP analysis of the S. haematobium ITS2 region; however, there was no variation in size among the fragments. Compared to the ITS2 sequences obtained for S. haematobium from Kenya, the nucleotide sequences of the ITS2 regions of S. haematobium from 4 areas in Sudan were consistent with those from Kenya (> 99%). In this study, we demonstrate for the first time that most of the S. haematobium population in Sudan consists of a pan-African S. haematobium genotype; however, we also report the discovery of Kenyan strain inflow into White Nile, Sudan.

Genetic diversity of Schistosoma haematobium parasite IS NOT associated with severity of disease in an endemic area in Sudan

Background

Over 650 million people globally are at risk of schistosomiasis infection, while more than 200 million people are infected of which the higher disease rates occur in children. Eighty three students between 6-20 years (mean 12.45 ± 3.2) from Quran School for boys in Radwan village, Gezira state were recruited to investigate for the relationship between the genetic diversity of Schistosoma haematobium strains and the severity of the disease.

Method

Schistosoma haematobium infection was detected by filtration of urine. Ultrasonography was done on each study subject, while PCR technique was used for genotyping via random amplified polymorphic DNA (RAPD) with A01, A02, A12, Y20 and A13 primers. A01 primer gave three different genotypes (A01-1, A01-2 and A01-3).

Results

About 54.2% (45/83) were S. haematobium egg positive by urine filtration. On assessment of the upper and lower urinary tract by ultrasound technique, 61.4% (51/83) were positiveand73.3% (60/83) samples were PCR positive. No significant difference was found when comparing the three different genotypes with severity of the disease.

Conclusion

This study concludes that no association was found between the different genotypes of S.haemtobium and the severity of the disease. Examination of more samples from different areas to identify any possible differences between the parasites genes and disease severity was recommended.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2334-14-469) contains supplementary material, which is available to authorized users.

Elucidating the temporal and spatial dynamics of Biomphalaria glabrata genetic diversity in three Brazilian villages

OBJECTIVE

The freshwater snail Biomphalaria glabrata is the principal intermediate host for the parasite Schistosoma mansoni within Brazil. We assessed the potential effects of snail population dynamics on parasite transmission dynamics via population genetics.

METHODS

We sampled snail populations located within the confines of three schistosome-endemic villages in the state of Minas Gerais, Brazil. Snails were collected from individual microhabitats following seasonal periods of flood and drought over the span of 1 year. Snail spatio-temporal genetic diversity and population differentiation of 598 snails from 12 sites were assessed at seven microsatellite loci.

RESULTS

Average genetic diversity was relatively low, ranging from 4.29 to 9.43 alleles per locus, and overall, subpopulations tended to exhibit heterozygote deficits. Genetic diversity was highly spatially partitioned among subpopulations, while virtually, no partitioning was observed across temporal sampling. Comparison with previously published parasite genetic diversity data indicated that S. mansoni populations are significantly more variable and less subdivided than those of the B. glabrata intermediate hosts.

DISCUSSION

Within individual Brazilian villages, observed distributions of snail genetic diversity indicate temporal stability and very restricted gene flow. This is contrary to observations of schistosome genetic diversity over the same spatial scale, corroborating the expectation that parasite gene flow at the level of individual villages is likely driven by vertebrate host movement.

Genetic diversity and population structure of Schistosoma mansoni within human infrapopulations in Mwea, central Kenya assessed by microsatellite markers

A recently developed high-throughput technique that allows multi-locus microsatellite analysis of individual miracidia of Schistosoma mansoni was used to assess the levels of genetic diversity and population structure in 12 infrapopulations of the parasite, each infrapopulation derived from an infected school child from the Mwea area, central Kenya. The mean number of alleles per locus was in the range 8.22-10.22, expected heterozygosity in Hardy-Weinberg equilibrium was 0.68-0.70, and pairwise F(ST) values ranged from 0.16% to 3.98% for the 12 infrapopulations. Although the genetic diversity within each infrapopulation of S. mansoni in this area was generally high, low levels of genetic structure were observed, suggestive of high levels of gene flow among infrapopulations. Private alleles were found in 8 of the 12 infrapopulation, the highest number of private alleles recorded per infrapopulation was 3. Our data suggest that the level of gene flow among infrapopulations of S. mansoni in Mwea is extremely high, thus providing opportunity for spread of rare alleles, including those that may confer character traits such as drug resistance and virulence.

Efficacy and side effects of praziquantel treatment against Schistosoma haematobium infection among primary school children in Zimbabwe

We examined the efficacy of praziquantel against Schistosoma haematobium among primary school children during a school-based deworming programme in the Burma Valley commercial farming area and the Nyamaropa rural areas in Zimbabwe, where the disease is highly endemic. Among 767 individuals infected with S. haematobium, 675 (88.0%) received treatment. Two single oral doses of 40mg/kg praziquantel were given 6 weeks apart. Of the 675 participants, heavy infection intensity was more common in males than females (chi(2)=6.61, P=0.010). Six weeks later, 624 participants (92.4%) were successfully followed up. The overall cure rate was 88.5% and the egg reduction rate was 98.2%. The highest cure rate was among those individuals with light infection. Seventy-two individuals remained infected at 6 weeks post treatment, among which 3 and 69 individuals had heavy and light infection, respectively. Forty-six of these children resolved following a second round of treatment at 6 weeks follow-up. Of the remaining children successfully followed-up, 22 resolved after a third round of treatment 6 months later. A wide range of observed mild and transient side effects were not associated with egg intensity. The parasitological cure rate was not associated with gender or age. Our study demonstrates that praziquantel is efficacious against S. haematobium in Zimbabwe, although low levels of persistent infection warrant further investigation.

Multiple near-identical genotypes of Schistosoma japonicum can occur in snails and have implications for population-genetic analyses

We genotyped (using 16 or 17 microsatellite loci) numerous adult Schistosoma japonicum raised in rabbits exposed to pooled cercariae from small numbers of naturally infected snails from several localities in China. As expected, duplicate multi-locus genotypes (MLGs) were found among these worms. Additionally, many more MLGs, often near-identical, were found than snails used as sources of cercariae. Explanations for these results include (i) genotyping errors, (ii) development within each infected snail of multiple sibling miracidia and (iii) somatic mutation producing genetically varied cercariae from a single miracidium. To control for genotyping errors we re-analysed samples from many individual worms, including repeating the initial PCR. Explanations invoking the development of multiple sibling miracidia within a single snail are not likely to be correct because almost all duplicate MLGs fell within same-sex clusters in a principal coordinates analysis. We would expect both sexes to be represented in a multi-miracidium infection. In addition, we exposed several snails to infection by a single miracidium. One such snail, via an experimentally infected mouse, yielded 48 adult worms. The presence of at least nine near-identical MLGs among these worms was confirmed by re-genotyping. We regard somatic mutation as the most likely explanation for our results. The implications of multiple MLGs for population-genetic studies in S. japonicum are discussed.

An analysis of genetic diversity and inbreeding in Wuchereria bancrofti: implications for the spread and detection of drug resistance

Estimates of genetic diversity in helminth infections of humans often have to rely on genotyping (immature) parasite transmission stages instead of adult worms. Here we analyse the results of one such study investigating a single polymorphic locus (a change at position 200 of the beta-tubulin gene) in microfilariae of the lymphatic filarial parasite Wuchereria bancrofti. The presence of this genetic change has been implicated in benzimidazole resistance in parasitic nematodes of farmed ruminants. Microfilariae were obtained from patients of three West African villages, two of which were sampled prior to the introduction of mass drug administration. An individual-based stochastic model was developed showing that a wide range of allele frequencies in the adult worm populations could have generated the observed microfilarial genetic diversity. This suggests that appropriate theoretical null models are required in order to interpret studies that genotype transmission stages. Wright's hierarchical F-statistic was used to investigate the population structure in W. bancrofti microfilariae and showed significant deficiency of heterozygotes compared to the Hardy-Weinberg equilibrium; this may be partially caused by a high degree of parasite genetic differentiation between hosts. Studies seeking to quantify accurately the genetic diversity of helminth populations by analysing transmission stages should increase their sample size to account for the variability in allele frequency between different parasite life-stages. Helminth genetic differentiation between hosts and non-random mating will also increase the number of hosts (and the number of samples per host) that need to be genotyped, and could enhance the rate of spread of anthelmintic resistance.

Genetic diversity and population structuring of Schistosoma mansoni in a Brazilian village

The digenean trematode Schistosoma mansoni is responsible for chronic schistosomiasis worldwide, and in Brazil alone an estimated 35 million people are at risk. To evaluate epidemiological patterns among human definitive hosts, we assessed genetic diversity and population subdivision of S. mansoni infrapopulations in human hosts from the highly endemic village of Virgem das Graças in the state of Minas Gerais, Brazil. We believe this is the largest such survey to date. Genetic diversity of parasites, measured over eight polymorphic microsatellite loci, was relatively high and standard measures of inbreeding indicated that the population was panmictic. Furthermore, there was no significant isolation-by-distance of parasite infrapopulations, and measures of population subdivision indicated significant but low to moderate levels of population differentiation. We conclude that patients within this village sample from a broad range of schistosome genetic diversity and effectively act as "genetic mixing bowls" for the parasites. These results contrast with those previously observed in the Brazilian village of Melqui ades and thus provide the opportunity for comparisons of environmental and epidemiological differences that are likely to influence host-parasite coevolution and parasite transmission.

Clonal diversity of the marine trematode Maritrema novaezealandensis within intermediate hosts: the molecular ecology of parasite life cycles

We quantified the clonal diversity of the New Zealand marine trematode Maritrema novaezealandensis (n = 1250) within Zeacumantus subcarinatus snail (n = 25) and Macrophthalmus hirtipes crab (n = 25) intermediate hosts using four to six microsatellite loci, and investigated the potential biological and physical factors responsible for the observed genetic patterns. Individual snails harboured one to five trematode genotypes and 48% of snails were infected by multiple parasite genotypes. Overall, the number of parasite genotypes did not increase with snail size, but was highest in intermediate-sized snails. Significantly larger numbers of parasite genotypes were detected in crabs (relative to snails; P < 0.001), with 16-25 genotypes recovered from individual crabs. Although crabs are typically infected by small numbers of cercariae sourced from many snails, they are occasionally infected by large numbers of cercariae sourced from single snails. The latter cases explain the significant genetic differentiation of trematode populations detected among their crab hosts (F(ST) = 0.009, P < 0.001). Our results suggest that the timing of infection and/or intraspecific competition among parasite clones within snails determine(s) the diversity of parasite clones that snails harbour. The presence of a large number of infected snails and tidal mixing of cercariae prior to infection results in crabs potentially harbouring hundreds of parasite genotypes despite the crabs' territorial behaviour.

Development and application of an ethically and epidemiologically advantageous assay for the multi-locus microsatellite analysis of Schistosoma mansoni

Non-availability of adult worms from living hosts remains a key problem in population genetic studies of schistosomes. Indirect sampling involving passage through laboratory animals presents significant ethical and practical drawbacks, and may result in sampling biases such as bottlenecking processes and/or host-induced selection pressures. The novel techniques reported here for sampling, storage and multi-locus microsatellite analysis of larval Schistosoma mansoni, allowing genotyping of up to 7 microsatellite loci from a single larva, circumvent these problems. The utility of these assays and the potential problems of laboratory passage, were evaluated using 7 S. mansoni population isolates collected from school-children in the Hoima district of Uganda, by comparing the associated field-collected miracidia with adult worms and miracidia obtained from a single generation in laboratory mice. Analyses of laboratory-passaged material erroneously indicated the presence of geographical structuring in the population, emphasizing the dangers of indirect sampling for population genetic studies. Bottlenecking and/or other sampling effects were demonstrated by reduced variability of adult worms compared to their parent field-collected larval samples. Patterns of heterozygote deficiency were apparent in the field-collected samples, which were not evident in laboratory-derived samples, potentially indicative of heterozygote advantage in establishment within laboratory hosts. Genetic distance between life-cycle stages in the majority of isolates revealed that adult worms and laboratory-passaged miracidia clustered together whilst segregating from field miracidia, thereby further highlighting the utility of this assay.

Molecular ecology of parasites: elucidating ecological and microevolutionary processes

We review studies that have used molecular markers to address ecological and microevolutionary processes in parasites. Our goal is to highlight areas of research that may be of particular interest in relation to the parasitic lifestyle, and to draw attention to areas that require additional study. Topics include species identification, phylogeography, host specificity and speciation, population genetic structure, modes of reproduction and transmission patterns, and searching for loci under selection.

Molecular ecology of Schistosoma mansoni transmission inferred from the genetic composition of larval and adult infrapopulations within intermediate and definitive hosts

We investigated the genotypic composition of the digenetic parasite Schistosoma mansoni for its adult stages within the definitive host (the wild rat, Rattus rattus) and for the larval stages within the intermediate host (the snail, Biomphalaria glabrata) both collected at the same transmission site. Our analyses are based upon the recognition and distribution of 200 different multilocus genotypes generated by RAPD markers. While intramolluscan larval infrapopulations are characterized by a low infection rate (0.6 % on average) and low intra-host genetic diversity (1.1 genotype on average per infected snail), adult infrapopulations within rats showed a high infection rate (94%) and a substantial intra-host genetic diversity (34 genotypes on average) linked to high intensities (160 worms per host on average). A single definitive host bearing 105 different genotypes harboured 52 % of the total genetic diversity detected within the whole parasite population. Analysis of the genetic data allowed the identification of various ecological, behavioural and immunological factors which are likely to enhance transmission of multiple parasite genotypes towards the vertebrate hosts. From the distribution of repeated identical multilocus genotypes within the parasite population and among the hosts, we have inferred different parameters of the cercarial transmission efficiency as well as patterns and processes by which vertebrate hosts acquire infection in the field.

Detection of Genetic Variation in Taenia solium

Genetic variability among Taenia solium isolates was studied in 160 cysticerci from 6 pigs, 4 from Mexico, 1 from Honduras, and 1 from Tanzania. Random amplified polymorphic DNA (RAPD) analysis performed with 4 commercial primers showed 88% polymorphic loci and an average heterozygosity of 0.077; however, several alleles were fixed within each isolate. Linkage disequilibrium analysis indicated that 3 of the 6 isolates had a random association of alleles, whereas the other 3 had a clonal structure. These results suggest the existence of local lineages in T. solium, with events of genetic recombination within them.

Schistosome genetic diversity: the implications of population structure as detected with microsatellite markers

Blood flukes in the genus Schistosoma are important human parasites in tropical regions. A substantial amount of genetic diversity has been described in populations of these parasites using molecular markers. We first consider the extent of genetic variation found in Schistosoma mansoni and some factors that may be contributing to this variation. Recently, though, attempts have been made to analyze not only the genetic diversity but how that diversity is partitioned within natural populations of schistosomes. Studies with non-allelic molecular markers (e.g. RAPDs and mtVNTRs) have indicated that schistosome populations exhibit varying levels of gene flow among component subpopulations. The recent characterization of microsatellite markers for S. mansoni provided an opportunity to study schistosome population structure within a population of schistosomes from a single Brazilian village using allelic markers. Whereas the detection of population structure depends strongly on the type of analysis with a mitochondrial marker, analyses with a set of seven microsatellite loci consistently revealed moderate genetic differentiation when village boroughs were used to define parasite subpopulations and greater subdivision when human hosts defined subpopulations. Finally, we discuss the implications that such strong population structure might have on schistosome epidemiology.