Genetic filtering and optimal sampling of Schistosoma mansoni populations

Compatibility Polymorphism Based on Long-Term Host-Parasite Relationships: Cross Talking Between Biomphalaria glabrata and the Trematode Schistosoma mansoni From Endemic Areas in Brazil

Sympatric snail populations have been kept in the laboratory since the isolation of the parasite from the field. To evaluate the influence of the intermediate host in the infectivity of S. mansoni, this allopatric strain was compared to two sympatric strains, from different geographical origins, and with different time of maintenance in the laboratory. Snail–trematode compatibility was accessed for a total of nine possible combinations (three snail populations, three schistosome strains), using different charges of parasite: 1, 5, 10, and 15 miracidia/snail. Each S. mansoni strain was characterized according to its infectivity phenotype that reflects the efficiency of their infection mechanism and all B. glabrata populations were characterized according to its (in)compatible phenotype that reflects the level of (un)susceptibility they display. For all host-parasite combinations tested the dose-response relation indicated a trend for an increase in the infectivity of S. mansoni when higher miracidial doses were used. SmRES-2 presented the highest overall infectivity rate, especially in the SmRES-2/BgRES interaction with 15 miracidia/snail. However, SmRES was more infective to BgBAR than SmRES-2, indicating that SmRES strain was more infective at the first contact with this new host than after 2 years of interaction (SmRES-2). BgBAR presented the highest susceptibility to infection. SmRES and SmRES-2 are the same parasite strains. It seems that during these 2 years of interaction, BgBAR acted like a filter and shifted the compatibility polymorphism of the strain SmRES. SmRES-2 became more infective to BgRES (sympatric) than to BgBAR (allopatric), and conversely, SmRES was more infective to BgBAR (allopatric) than to BgRES (sympatric). This interplay suggests that epigenetic mechanisms are prompting these changes. This study concerns with infection of B. glabrata snails from different Brazilian localities with S. mansoni in allopatric and sympatric associations that will partially help in understanding the natural epidemiology of schistosomiasis within natural snail populations in watercourses. This work demonstrates that there is a shift on the compatibility polymorphism profile resulting from sympatric and allopatric interactions of B. glabrata and S. mansoni that constantly change during the time of interaction.

The effect of sample size on estimates of genetic differentiation and effective population size for Schistosoma mansoni populations

Eradication or local extinction of the human parasite Schistosoma mansoni is a goal for many control programs. Population genetic analyses are helping to evaluate and guide these efforts, yet what to sample, how to sample and how densely to sample is not well established. We determined the S. mansoni allele frequency profile of nearly all infected inhabitants in two small Brazilian communities and created sub-samples representing 5-50% of all detected human infections (infrapopulations). Samples were selected at random with replacement, and each size class was replicated 100 times. Mean pairwise differentiation for all infrapopulations (Di) and the variance effective population size (Ne) were calculated for each sample. Prior to community-wide treatment, the true mean Di was moderate (0.095-0.123) and Ne large (>30,000). Most samples of <50% of those infected produced estimates outside of 5% of the true value. For estimates within 10%, sample sizes of >15% of all infrapopulations were required. At the 3 year follow-up after treatment, the Di increased and Ne was reduced by >15 fold. At this time sampling of >30-45% was needed to achieve the same accuracy. Following a second treatment and 4 years from baseline, the Di further increased and Ne decreased with little change in the sampling effort required. Extensive sampling is required for accurate estimates of these important population parameters. Characteristics such as population census size, infection prevalence, the community's treatment history and the degree of infrapopulation differentiation should be taken into account. The intensity of infection was weakly correlated with the ability of a single infrapopulation to represent the component population (Dic), indicating a tendency toward random acquisition of parasite genotypes. This also suggests that targeted sampling from those most heavily infected will better represent the genetic diversity of the whole community than a random sample of infrapopulations.

High genetic variability of Schistosoma haematobium in Mali and Nigeria

Schistosoma haematobium is one of the most prevalent parasitic flatworms, infecting over 112 million people in Africa. However, little is known about the genetic diversity of natural S. haematobium populations from the human host because of the inaccessible location of adult worms in the host. We used 4 microsatellite loci to genotype individually pooled S. haematobium eggs directly from each patient sampled at 4 endemic locations in Africa. We found that the average allele number of individuals from Mali was significantly higher than that from Nigeria. In addition, no significant difference in allelic composition was detected among the populations within Nigeria; however, the allelic composition was significantly different between Mali and Nigeria populations. This study demonstrated a high level of genetic variability of S. haematobium in the populations from Mali and Nigeria, the 2 major African endemic countries, suggesting that geographical population differentiation may occur in the regions.

Field-derived Schistosoma mansoni and Biomphalaria pfeifferi in Kenya: a compatible association characterized by lack of strong local adaptation, and presence of some snails able to persistently produce cercariae for over a year

Background

Schistosoma mansoni is widely distributed in sub-Saharan Africa with Biomphalaria pfeifferi being its most widespread and important snail intermediate host. Few studies have examined the compatibility of field-derived B. pfeifferi snails with S. mansoni miracidia derived from human hosts. We investigated compatibility (as defined by shedding of cercariae following exposure to miracidia) of two isolates of S. mansoni from school children from Asao (western Kenya) and Mwea (central Kenya) with B. pfeifferi collected directly from Asao stream or the Mwea rice fields.

Methods

We exposed snails from both regions to four different doses of miracidia (1, 5, 10 and 25) from sympatric or allopatric S. mansoni, and maintained them in a shaded, screened out-of-doors rearing facility in Kisian, in western Kenya. Both snail survival and the number of snails that became infected were monitored weekly. This was done for 25 weeks post-exposure (PE). Those infected snails which survived beyond this period were monitored until they all died.

Results

Although overall survival of Mwea snails maintained in western Kenya was generally low, both sympatric and allopatric combinations of parasites and snails exhibited high compatibility (approximately 50% at a dose of one miracidium per snail), with an increase in infection rates as the miracidial dose was increased (P < 0.002). Schistosomes were no more compatible with sympatric than allopatric snails, nor were snails less compatible with sympatric than allopatric schistosomes. Snail mortality increased significantly with dose of miracidia (P < 0.05). Approximately 3% of Asao snails exposed to a low dose of sympatric miracidia (1 or 5) continued to shed cercariae for as long as 58 weeks post exposure.

Conclusions

There were no significant local adaptation effects for either schistosomes or snails. Also, the existence of “super-survivor” snails is noteworthy for its implications for current control initiatives that mostly rely on mass drug administration (MDA). Long-term shedders could provide an ongoing source of cercariae to initiate human infections for many months, suggesting care is required in considering how human MDA treatments are timed. Future control programs should incorporate means to eliminate infected snails to complement chemotherapy interventions in controlling schistosomiasis.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-014-0533-3) contains supplementary material, which is available to authorized users.

Isolation and characterization of the first polymorphic microsatellite markers for Schistosoma haematobium and their application in multiplex reactions of larval stages

The ability of microsatellite loci to reveal genetic diversity within the trematode Schistosoma haematobium is demonstrated for the first time. Nine novel polymorphic microsatellite markers were isolated and their viability assessed on 36 S. haematobium adult worm individuals from three geographical populations. Allelic diversity and gene diversity ranged from two to seven and from 0.29 to 0.76, respectively, suggesting high variability between individuals and between unrelated populations. Three primers also amplified Schistosoma mansoni and two Schistosoma japonicum. The results suggest these primers are useful for population genetic analyses of S. haematobium.

Population genetics of Schistosoma haematobium: development of novel microsatellite markers and their application to schistosomiasis control in Mali

The recent implementation of mass drug administration (MDA) for control of uro-genital schistosomiasis has identified an urgent need for molecular markers to both directly monitor the impact of MDA, for example to distinguish re-infections from uncleared infections, as well as understand aspects of parasite reproduction and gene flow which might predict evolutionary change, such as the development and spread of drug resistance. We report the development of a novel microsatellite tool-kit allowing, for the first time, robust genetic analysis of individual S. haematobium larvae collected directly from infected human hosts. We genotyped the parasite populations of 47 children from 2 schools in the Ségou region of Mali, the first microsatellite study of this highly neglected parasite. There was only limited evidence of population subdivision between individual children or between the two schools, suggesting that few barriers to gene flow exist in this population. Complex relationships between parasite reproductive success, infection intensity and host age and gender were identified. Older children and boys harboured more diverse infections, as measured by the number of unique adult genotypes present. Individual parasite genotypes had variable reproductive success both across hosts, a pre-requisite for evolutionary selection, and, phenotypically, in hosts of different ages and genders. These data serve as a baseline against which to measure the effect of treatment on parasite population genetics in this region of Mali, and the tools developed are suitable to further investigate this important pathogen, and its close relatives, throughout their range.

Optimal sample storage and extraction procotols for reliable multilocus genotyping of the human parasite Schistosoma mansoni

Genotyping individual larval stages and eggs of natural parasite populations is complicated by the difficulty of obtaining reliable genotypes from low quantity DNA template. A suitable storage and extraction protocol, together with a thorough quantification of genotyping errors are therefore crucial for molecular epidemiological studies. Here we test the robustness, handling time, ease of use, cost effectiveness and success rate of various fixation (Whatman FTA(®) Classic and Elute Cards, 70% EtOH and RNAlater(®)) and subsequent DNA extraction methods (commercial kits and proteinase K protocol). None of these methods require a cooling chain and are therefore suitable for field collection. Based on a multiplex microsatellite PCR with nine loci the success and reliability of each technique is evaluated by the proportion of samples with at least eight scored loci and the proportion of genotyping errors. If only the former is taken into account, FTA(®) Elute is recommended (83% success; 44% genotyping error; 0.2 €/sample; 1h 20 m handling time). However, when also considering the genotyping errors, handling time and ease of use, we opt for 70% EtOH with the 96-well plate technology followed by a simple proteinase K extraction (73% success; 0% genotyping error; 0.2 €/sample; 15m handling time). For eggs we suggest (1) to pool all eggs per person in 1.5 ml tubes filled with 70% EtOH for transport and (2) to identify each egg to species level prior to genotyping. To this end we extended the Rapid diagnostic PCR developed by Webster et al. (2010) with a S. mansoni-specific primer to discriminate between S. mansoni, S. haematobium and S. bovis in a single PCR reaction. The success rate of genotyping eggs was 75% (0% genotyping error). This is the first study to incorporate genotyping errors through re-amplification for the evaluation of schistosome sampling protocols and the identification of error-prone loci.

Transmission dynamics of two strains of Schistosoma mansoni utilizing novel intermediate and definitive hosts

The intimate host-parasite relationship mandates adaptation to the genetic and phenotypic variability of their counterparts. Here, inbred and outcrossed strains of Schistosoma mansoni were challenged with "local" and "novel" intermediate and definitive hosts to examine effects of genetic variability and novelty on infection success and dynamics. Genetically distinct lines of Biomphalaria glabrata intermediate hosts exposed to inbred and outcrossed S. mansoni larvae were assessed for differences in both snail and parasite life-history parameters. Cercariae from each parasite-snail treatment were used to infect "local" and "novel" Mus musculus definitive hosts to assess parasite infectivity and fitness. Outcrossed parasites significantly reduced snail growth, were more productive, and induced greater host mortality than inbred parasites. Mouse strain did not influence parasite infectivity or reproduction, but parasite and snail host genetic background did, affecting both sex-specific infectivity and parasite productivity. Overall, genetic background of S. mansoni and its intermediate snail host altered life history traits and transmission dynamics of the parasite throughout its life cycle.

Applying evolutionary genetics to schistosome epidemiology

We review how molecular markers and evolutionary analysis have been applied to the study of schistosome parasites, important pathogens that infect over 200 million people worldwide. Topics reviewed include phylogenetics and biogeography, hybridization, infection within snails, mating systems, and genetic structure. Some interesting generalizations include that schistosome species hybridize frequently and have switched definitive hosts repeatedly in evolutionary time. We show that molecular markers can be used to infer epidemiologically relevant processes such as spatial variation in transmission, or to reveal complex patterns of mate choice. Analysis of genetic structure data shows that transmission foci can be structured by watershed boundaries, habitat types, and host species. We also discuss sampling and analytical problems that arise when using larvae to estimate genetic parameters of adult schistosome populations. Finally, we review pitfalls in methodologies such as genotyping very small individuals, statistical methods for identifying clonemates or for identifying sibling groups, and estimating allele frequencies from pooled egg samples.

Analysis of Schistosoma mansoni population structure using total fecal egg sampling

Many parasite populations are difficult to sample because they are not uniformly distributed between several host species and are often not easily collected from the living host, thereby limiting sample size and possibly distorting the representation of the population. For the parasite Schistosoma mansoni, we investigated the use of eggs, in aggregate, from the stools of infected individuals as a simple and representative sample. Previously, we demonstrated that microsatellite allele frequencies can be accurately estimated from pooled DNA of cloned S. mansoni adults. Here, we show that genotyping of parasite populations from reproductively isolated laboratory strains can be used to identify these specific populations based on characteristic patterns of allele frequencies, as observed by polyacrylamide gel electrophoresis and automated sequencer analysis of fluorescently labeled PCR products. Microsatellites used to genotype aggregates of eggs collected from stools of infected individuals produced results consistent with the geographic distribution of the samples. Preferential amplification of smaller alleles, and stutter PCR products, had negligible effect on measurement of genetic differentiation. Direct analysis of total stool eggs can be an important approach to questions of population genetics for this parasite by increasing the sample size to thousands per infected individual and by reducing bias.

Efficient genotyping of Schistosoma mansoni miracidia following whole genome amplification

Small parasites and larval stages pose a problem for molecular analyses because limited amounts of DNA template are available. Isothermal methods for faithfully copying DNA have the potential to revolutionize studies of such organisms. We evaluated the fidelity of multiple displacement amplification (MDA) for amplifying DNA extracted from a single miracidium of Schistosoma mansoni. To do this we genotyped DNA extracted from 28 F1 miracidia following MDA using 56 microsatellite markers. Because these miracidia were obtained from a cross between a male and female worm of known genotypes, we were able to predict the alleles present in the progeny and quantify the genotyping error rate. We found just 8/1568 genotypes deviated from Mendelian expectations. Furthermore, because 1 of these resulted from a genuine mutation, the error rate due to MDA is 7/1568 (0.45%). We conclude that many hundreds of microsatellites or other genetic markers can be accurately genotyped from a single miracidium using this method, greatly expanding the scope of population genetic, epidemiological and evolutionary studies on this parasite.

A new approach to characterize populations of Schistosoma mansoni from humans: development and assessment of microsatellite analysis of pooled miracidia

OBJECTIVES

To develop and assess a microsatellite technique to characterize populations of Schistosoma mansoni from humans.

METHODS

For each of five patients, we calculated the allele count and frequency at 11 loci for several pools of miracidia (50 and 100), and compared these to population values, determined by amplifying microsatellites from 186 to 200 individual miracidia per patient.

RESULTS

We were able to detect up to 94.5% of alleles in pools. Allele count and frequency strongly and significantly correlated between singles and pools; marginally significant differences (P < 0.05) were detected for one patient (pools of 50) for allele frequencies and for two patients (pools of 100) for allele counts. Kato-Katz egg counts and number of alleles per pool did not co-vary, indicating that further direct comparisons of the results from these two techniques are needed.

CONCLUSIONS

Allele counts and frequency profiles from pooling provide important information about infection intensity and complexity, beyond that obtained from traditional methods. Although we are not advocating use of pooling to replace individual genotyping studies, it can potentially be useful in certain applications as a rapid and cost effective screening method for studies of S. mansoni population genetics, or as a more informative way to quantify and characterize human worm populations.

Schistosome monogamy: who, how, and why?

Schistosomes represent a unique animal model for comparative analyses of monogamy. Indeed, schistosomes are classified at the lowest taxonomical level of monogamous species and lack complex social interactions, which could alter our understanding of their unusual mating system. Elements discussed here include the fact that monogamy in schistosomes could be an ancestral state between hermaphroditism and polygyny or polygynandry and the occurrence of mate changes. In addition, hypotheses are proposed to explain monogamy in schistosomes (e.g. female dispersion, the need for paternal care, oviposition site limitation or aggressiveness, and mate guarding). We also propose future experimental and analytical approaches to improve our understanding of the schistosomes' mating system.

Effects of intermediate host genetic background on parasite transmission dynamics: a case study using Schistosoma mansoni

For parasites that require multiple hosts to complete their development, genetic interplay with one host may impact parasite transmission and establishment in subsequent hosts. In this study, we used microsatellite loci to address whether the genetic background of snail intermediate hosts influences life-history traits and transmission patterns of dioecious trematode parasites in their definitive hosts. We performed experimental Schistosoma mansoni infections utilizing two allopatric populations of Biomphalaria glabrata snails and assessed intensities and sex ratios of adult parasites in mouse definitive hosts. Our results suggest that the genetic background of hosts at one point in a parasite's life cycle can influence the intensities and sex ratios of worms in subsequent hosts.

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.

Making the most of mitochondrial genomes--markers for phylogeny, molecular ecology and barcodes in Schistosoma (Platyhelminthes: Digenea)

An increasing number of complete sequences of mitochondrial (mt) genomes provides the opportunity to optimise the choice of molecular markers for phylogenetic and ecological studies. This is particularly the case where mt genomes from closely related taxa have been sequenced; e.g., within Schistosoma. These blood flukes include species that are the causative agents of schistosomiasis, where there has been a need to optimise markers for species and strain recognition. For many phylogenetic and population genetic studies, the choice of nucleotide sequences depends primarily on suitable PCR primers. Complete mt genomes allow individual gene or other mt markers to be assessed relative to one another for potential information content, prior to broad-scale sampling. We assess the phylogenetic utility of individual genes and identify regions that contain the greatest interspecific variation for molecular ecological and diagnostic markers. We show that variable characters are not randomly distributed along the genome and there is a positive correlation between polymorphism and divergence. The mt genomes of African and Asian schistosomes were compared with the available intraspecific dataset of Schistosoma mansoni through sliding window analyses, in order to assess whether the observed polymorphism was at a level predicted from interspecific comparisons. We found a positive correlation except for the two genes (cox1 and nad1) adjoining the putative control region in S. mansoni. The genes nad1, nad4, nad5, cox1 and cox3 resolved phylogenies that were consistent with a benchmark phylogeny and in general, longer genes performed better in phylogenetic reconstruction. Considering the information content of entire mt genome sequences, partial cox1 would not be the ideal marker for either species identification (barcoding) or population studies with Schistosoma species. Instead, we suggest the use of cox3 and nad5 for both phylogenetic and population studies. Five primer pairs designed against Schistosoma mekongi and Schistosoma malayensis were tested successfully against Schistosoma japonicum. In combination, these fragments encompass 20-27% of the variation amongst the genomes (average total length approximately 14,000bp), thus providing an efficient means of encapsulating the greatest amount of variation within the shortest sequence. Comparative mitogenomics provides the basis of a rational approach to molecular marker selection and optimisation.