The genome and transcriptome of the snail Biomphalaria sudanica s.l.: immune gene diversification and highly polymorphic genomic regions in an important African vector of Schistosoma mansoni

BACKGROUND

Control and elimination of schistosomiasis is an arduous task, with current strategies proving inadequate to break transmission. Exploration of genetic approaches to interrupt Schistosoma mansoni transmission, the causative agent for human intestinal schistosomiasis in sub-Saharan Africa and South America, has led to genomic research of the snail vector hosts of the genus Biomphalaria. Few complete genomic resources exist, with African Biomphalaria species being particularly underrepresented despite this being where the majority of S. mansoni infections occur. Here we generate and annotate the first genome assembly of Biomphalaria sudanica sensu lato, a species responsible for S. mansoni transmission in lake and marsh habitats of the African Rift Valley. Supported by whole-genome diversity data among five inbred lines, we describe orthologs of immune-relevant gene regions in the South American vector B. glabrata and present a bioinformatic pipeline to identify candidate novel pathogen recognition receptors (PRRs).

RESULTS

De novo genome and transcriptome assembly of inbred B. sudanica originating from the shoreline of Lake Victoria (Kisumu, Kenya) resulted in a haploid genome size of ~ 944.2 Mb (6,728 fragments, N50 = 1.067 Mb), comprising 23,598 genes (BUSCO = 93.6% complete). The B. sudanica genome contains orthologues to all described immune genes/regions tied to protection against S. mansoni in B. glabrata, including the polymorphic transmembrane clusters (PTC1 and PTC2), RADres, and other loci. The B. sudanica PTC2 candidate immune genomic region contained many PRR-like genes across a much wider genomic region than has been shown in B. glabrata, as well as a large inversion between species. High levels of intra-species nucleotide diversity were seen in PTC2, as well as in regions linked to PTC1 and RADres orthologues. Immune related and putative PRR gene families were significantly over-represented in the sub-set of B. sudanica genes determined as hyperdiverse, including high extracellular diversity in transmembrane genes, which could be under pathogen-mediated balancing selection. However, no overall expansion in immunity related genes was seen in African compared to South American lineages.

CONCLUSIONS

The B. sudanica genome and analyses presented here will facilitate future research in vector immune defense mechanisms against pathogens. This genomic/transcriptomic resource provides necessary data for the future development of molecular snail vector control/surveillance tools, facilitating schistosome transmission interruption mechanisms in Africa.

The genome and transcriptome of the snail Biomphalaria sudanica s.l.: Immune gene diversification and highly polymorphic genomic regions in an important African vector of Schistosoma mansoni

Background

Control and elimination of schistosomiasis is an arduous task, with current strategies proving inadequate to break transmission. Exploration of genetic approaches to interrupt Schistosoma mansoni transmission, the causative agent for human intestinal schistosomiasis in sub-Saharan Africa and South America, has led to genomic research of the snail vector hosts of the genus Biomphalaria. Few complete genomic resources exist, with African Biomphalaria species being particularly underrepresented despite this being where the majority of S. mansoni infections occur. Here we generate and annotate the first genome assembly of Biomphalaria sudanica sensu lato, a species responsible for S. mansoni transmission in lake and marsh habitats of the African Rift Valley. Supported by whole-genome diversity data among five inbred lines, we describe orthologs of immune-relevant gene regions in the South American vector B. glabrata and present a bioinformatic pipeline to identify candidate novel pathogen recognition receptors (PRRs).

Results

De novo genome and transcriptome assembly of inbred B. sudanica originating from the shoreline of Lake Victoria (Kisumu, Kenya) resulted in a haploid genome size of ~944.2 Mb (6732 fragments, N50=1.067 Mb), comprising 23,598 genes (BUSCO=93.6% complete). The B. sudanica genome contains orthologues to all described immune genes/regions tied to protection against S. mansoni in B. glabrata. The B. sudanica PTC2 candidate immune genomic region contained many PRR-like genes across a much wider genomic region than has been shown in B. glabrata, as well as a large inversion between species. High levels of intra-species nucleotide diversity were seen in PTC2, as well as in regions linked to PTC1 and RADres orthologues. Immune related and putative PRR gene families were significantly over-represented in the sub-set of B. sudanica genes determined as hyperdiverse, including high extracellular diversity in transmembrane genes, which could be under pathogen-mediated balancing selection. However, no overall expansion in immunity related genes were seen in African compared to South American lineages.

Conclusions

The B. sudanica genome and analyses presented here will facilitate future research in vector immune defense mechanisms against pathogens. This genomic/transcriptomic resource provides necessary data for the future development of molecular snail vector control/surveillance tools, facilitating schistosome transmission interruption mechanisms in Africa.

Multi-strain compatibility polymorphism between a parasite and its snail host, a neglected vector of schistosomiasis in Africa

Interactions between Schistosoma mansoni and its snail host are understood primarily through experimental work with one South American vector species, Biomphalaria glabrata. However, 90% of schistosomiasis transmission occurs in Africa, where a diversity of Biomphalaria species may serve as vectors. With the long-term goal of determining the genetic and ecological determinants of infection in African snail hosts, we developed genetic models of Biomphalaria sudanica, a principal vector in the African Great Lakes. We determined laboratory infection dynamics of two S. mansoni lines in four B. sudanica lines. We measured the effects of the following variables on infection success and the number of cercariae produced (infection intensity): (i) the combination of parasite and snail line; (ii) the dose of parasites; and (iii) the size of snail at time of exposure. We found one snail line to be almost completely incompatible with both parasite lines, while other snail lines showed a polymorphism in compatibility: compatible with one parasite line while incompatible with another. Interestingly, these patterns were opposite in some of the snail lines. The parasite-snail combination had no significant effect on the number of cercariae produced in a successful infection. Miracidia dose had a strong effect on infection status, in that higher doses led to a greater proportion of infected snails, but had no effect on infection intensity. In one of the snail-schistosome combinations, snail size at the time of exposure affected both infection status and cercarial production in that the smallest size class of snails (1.5-2.9 mm) had the highest infection rates, and produced the greatest number of cercariae, suggesting that immunity increases with age and development. The strongest predictor of the infection intensity was the size of snail at the time of shedding: 1 ​mm of snail growth equated to a 19% increase in cercarial production. These results strongly suggest that infection status is determined in part by the interaction between snail and schistosome genetic lines, consistent with a gene-for-gene or matching allele model. This foundational work provides rationale for determining the genetic interactions between African snails and schistosomes, which may be applied to control strategies.

Genotypic-specific heat shock response of vector susceptibility to Schistosoma mansoni

Living organisms are vulnerable to thermal stress which causes a diversity of physiological outcomes. Previous work has shown that the snail vectors (Biomphalaria glabrata) of an important human pathogen, Schistosoma mansoni, revert from resistant to susceptible after short exposure to a heat stress as low as 31oC; however, due to lack of replicability among labs and genetic lines of snails, it has been hypothesized that this effect is genotype dependent. We examined the effects of heat shock on resistance of two species of snail vectors including B. glabrata and B. sudanica. We used 3 different inbred laboratory snail lines in addition to the F1 generation of field collected snails from Lake Victoria, Kenya, an area with high levels of schistosomiasis transmission. Our results showed marginal effects of heat shock on prevalence of infection in B. glabrata, and that this response was genotype specific. We found no evidence of a heat shock effect on prevalence of infection in B. sudanica or on intensity of infection (number of infectious stages shed) in either snail species. Such environmentally influenced defense responses stress the importance of considering this unique interaction between snail and parasite genotypes in determining infection dynamics under climate changes.

Host preference of field‐derived Schistosoma mansoni is influenced by snail host compatibility and infection status

Schistosome parasites cause a chronic inflammatory disease in humans, and recent studies have emphasized the importance of control programs for understanding the aquatic phases of schistosomiasis transmission. The host-seeking behavior of larval schistosomes (miracidia) for their snail intermediate hosts plays a critical role in parasite transmission. Using field-derived strains of Kenyan snails and parasites, we tested two main hypotheses: (1) Parasites prefer the most compatible host, and (2) parasites avoid hosts that are already infected. We tested preference to three Biomphalaria host snail taxa (B. pfeifferi, B. sudanica, and B. choanomphala), using allopatric and sympatric Schistosoma mansoni isolates and two different nonhost snail species that co-occur with Biomphalaria, Bulinus globosus, and Physa acuta. We also tested whether schistosomes avoid snail hosts that are already infected by another trematode species and whether competitive dominance played a role in their behavior. Preference was assessed using two-way choice chambers and by visually counting parasites that moved toward competing stimuli. In pairwise comparisons, we found that S. mansoni did not always prefer the more compatible snail taxon, but never favored an incompatible host over a compatible host. While parasites preferred B. pfeifferi to the nonhost species B. globosus, they did not significantly prefer B. pfeifferi versus P. acuta, an introduced species in Kenya. Finally, we demonstrated that parasites avoid infected snails if the resident parasite was competitively dominant (Patagifer sp.), and preferred snails infected with subordinates (xiphidiocercariae) to uninfected snails. These results provide evidence of “fine tuning” in the ability of schistosome miracidia to detect hosts; however, they did not always select hosts that would maximize fitness. Appreciating such discriminatory abilities could lead to a better understanding of how ecosystem host and parasite diversity influences disease transmission and could provide novel control mechanisms to improve human health.

Genomic analysis of a parasite invasion: Colonization of the Americas by the blood fluke Schistosoma mansoni

Schistosoma mansoni, a snail-borne, blood fluke that infects humans, was introduced into the Americas from Africa during the Trans-Atlantic slave trade. As this parasite shows strong specificity to the snail intermediate host, we expected that adaptation to South American Biomphalaria spp. snails would result in population bottlenecks and strong signatures of selection. We scored 475,081 single nucleotide variants in 143 S. mansoni from the Americas (Brazil, Guadeloupe and Puerto Rico) and Africa (Cameroon, Niger, Senegal, Tanzania, and Uganda), and used these data to ask: (i) Was there a population bottleneck during colonization? (ii) Can we identify signatures of selection associated with colonization? (iii) What were the source populations for colonizing parasites? We found a 2.4- to 2.9-fold reduction in diversity and much slower decay in linkage disequilibrium (LD) in parasites from East to West Africa. However, we observed similar nuclear diversity and LD in West Africa and Brazil, suggesting no strong bottlenecks and limited barriers to colonization. We identified five genome regions showing selection in the Americas, compared with three in West Africa and none in East Africa, which we speculate may reflect adaptation during colonization. Finally, we infer that unsampled populations from central African regions between Benin and Angola, with contributions from Niger, are probably the major source(s) for Brazilian S. mansoni. The absence of a bottleneck suggests that this is a rare case of a serendipitous invasion, where S. mansoni parasites were pre-adapted to the Americas and able to establish with relative ease.

Comparative Vectorial Competence of Biomphalaria sudanica and Biomphalaria choanomphala, Snail Hosts of Schistosoma mansoni, From Transmission Hotspots In Lake Victoria, Western Kenya

Schistosoma mansoni, which causes human intestinal schistosomiasis, continues to be a major public health concern in the Lake Victoria basin in western Kenya, with Biomphalaria sudanica (a shoreline inhabiting snail) and Biomphalaria choanomphala (a deep-water snail) playing roles in transmission. A recent study showed that B. sudanica was abundantly present near all study villages on the lakeshore, but B. choanomphala was significantly more abundant near villages known to be persistent transmission hotspots. The present study investigated the relative compatibility of B. sudanica and B. choanomphala with S. mansoni. A reciprocal cross-infection experiment used young adult F1 generation B. sudanica and B. choanomphala that were exposed to either 1, 5, or 10 sympatric or allopatric human-derived S. mansoni miracidia. Three weeks post-exposure (PE) and weekly thereafter, the snails were counted and screened for schistosome cercariae, and at 7 wk PE, total cercariae shed during a 2 hr period by each infected snail was determined. Pre-patent periods for S. mansoni in both B. sudanica and B. choanomphala were similar, and most snails in all exposure combinations started shedding cercariae 5 wk PE. Prevalences were significantly higher in B. choanomphala (12.2-80.9%) than in B. sudanica (5.2-18.6%) at each dose, regardless of whether miracidia were of an allopatric or a sympatric source (P < 0.0001). Overall, the odds of a snail becoming infected with 5 or 10 miracidia were significantly higher than the odds of being infected with 1 miracidium, (P < 0.0001), and fewer cercariae were produced by snails exposed to single as compared to 5 or 10 miracidia. On average, B. choanomphala produced more cercariae ( = 458, SD = 414) than B. sudanica ( = 238, SD = 208) (P < 0.0001). These results suggest that B. choanomphala is more compatible with S. mansoni than B. sudanica. Though B. choanomphala can be found in shallow shoreline waters, it is, for the most part, a deeper-water taxon. Because dredging is a relatively inefficient means of sampling, B. choanomphala is likely underestimated with respect to its population size, the number of S. mansoni-positive snails, and its role in maintaining transmission.

A Search for Snail-Related Answers to Explain Differences in Response of Schistosoma mansoni to Praziquantel Treatment among Responding and Persistent Hotspot Villages along the Kenyan Shore of Lake Victoria

Following a 4-year annual praziquantel (PZQ) treatment campaign, the resulting prevalence of Schistosoma mansoni was seen to differ among individual villages along the Kenyan shore of Lake Victoria. We have investigated possible inherent differences in snail-related aspects of transmission among such 10 villages, including six persistent hotspot (PHS) villages (≤ 30% reduction in prevalence following repeated treatments) located along the west-facing shore of the lake and four PZQ-responding (RESP) villages (> 30% prevalence reduction following repeated treatment) along the Winam Gulf. When taking into account all sampling sites, times, and water hyacinth presence/absence, shoreline-associated Biomphalaria sudanica from PHS and RESP villages did not differ in relative abundance or prevalence of S. mansoni infection. Water hyacinth intrusions were associated with increased B. sudanica abundance. The deeper water snail Biomphalaria choanomphala was significantly more abundant in the PHS villages, and prevalence of S. mansoni among villages both before and after control was positively correlated with B. choanomphala abundance. Worm recoveries from sentinel mice did not differ between PHS and RESP villages, and abundance of non-schistosome trematode species was not associated with S. mansoni abundance. Biomphalaria choanomphala provides an alternative, deepwater mode of transmission that may favor greater persistence of S. mansoni in PHS villages. As we found evidence for ongoing S. mansoni transmission in all 10 villages, we conclude that conditions conducive for transmission and reinfection occur ubiquitously. This argues for an integrated, basin-wide plan for schistosomiasis control to counteract rapid reinfections facilitated by large snail populations and movements of infected people around the lake.

Acanthocephalan Parasites of the Oarfish, Regalecus russelii (Regalecidae), With A Description of A New Species of Gymnorhadinorhynchus (Acanthocephala: Gymnorhadinorhynchidae)

Oarfish are rarely seen and seldom studied, which makes their parasite fauna even more enigmatic. Necropsy of 12 oarfish, Regalecus russelii (Regalecidae) (Cuvier, 1816), from Japan yielded 2 species of acanthocephalans. One species was found in 2 oarfish and a total of 76 specimens was collected, but only a single, immature specimen of the second species was found. The former represents an undescribed species from the order Echinorhynchida and is described here. Morphological and phylogenetic analyses of the small subunit ( SSU) rDNA place this species in the family Gymnorhadinorhynchidae, and genus Gymnorhadinorhynchus which is characterized by a cylindrical proboscis with longitudinal rows of hooks, basal circle of enlarged hooks, asymmetry of hook shape, 4 cement glands, and a spineless trunk. Diagnostic characters of this species within the genus are the number of longitudinal rows of hooks (14), smaller body size (males: 4.8-6.6 mm and females: 5.3-6.3 mm) and a number of molecular autapomorphies including a number of long insertions in both the SSU and large subunit rDNA (LSU). A single immature female of Bolbosoma sp. (Palaeacanthocephala: Plagiorhynchidae) was also found with its anterior end embedded in the mucosa of the pyloric ceca. The characters of this specimen are not consistent with any other known species of Bolbosoma; however, because only 1 immature specimen with a partially invaginated proboscis was recovered, it was not designated as a new species.

Host immunity, nutrition and coinfection alter longitudinal infection patterns of schistosomes in a free ranging African buffalo population

Schistosomes are trematode parasites of global importance, causing infections in millions of people, livestock, and wildlife. Most studies on schistosomiasis, involve human subjects; as such, there is a paucity of longitudinal studies investigating parasite dynamics in the absence of intervention. As a consequence, despite decades of research on schistosomiasis, our understanding of its ecology in natural host populations is centered around how environmental exposure and acquired immunity influence acquisition of parasites, while very little is known about the influence of host physiology, coinfection and clearance in the absence of drug treatment. We used a 4-year study in free-ranging African buffalo to investigate natural schistosome dynamics. We asked (i) what are the spatial and temporal patterns of schistosome infections; (ii) how do parasite burdens vary over time within individual hosts; and (iii) what host factors (immunological, physiological, co-infection) and environmental factors (season, location) explain patterns of schistosome acquisition and loss in buffalo? Schistosome infections were common among buffalo. Microgeographic structure explained some variation in parasite burdens among hosts, indicating transmission hotspots. Overall, parasite burdens ratcheted up over time; however, gains in schistosome abundance in the dry season were partially offset by losses in the wet season, with some hosts demonstrating complete clearance of infection. Variation among buffalo in schistosome loss was associated with immunologic and nutritional factors, as well as co-infection by the gastrointestinal helminth Cooperia fuelleborni. Our results demonstrate that schistosome infections are surprisingly dynamic in a free-living mammalian host population, and point to a role for host factors in driving variation in parasite clearance, but not parasite acquisition which is driven by seasonal changes and spatial habitat utilization. Our study illustrates the power of longitudinal studies for discovering mechanisms underlying parasite dynamics in individual animals and populations.

Genotypic variation in host response to infection affects parasite reproductive rate

Parasite fitness is largely influenced by a variation in host response due to the host's genetic background. Here we investigated the impact of host genotype on pathogen success in the snail vector of its castrating parasite, Schistosoma mansoni. We infected five inbred lines of Biomphalaria glabrata with two infection doses and followed their growth, reproductive output and parasite production throughout the course of infection. There was no difference in resistance to infection among inbred lines, but lines varied in their responses to infection and the numbers of parasites produced. Snails did not compensate for castration by increasing their fecundity during the early phase of infection (fecundity compensation). However, some lines were able to delay parasite shedding for up to 30 weeks, thus prolonging reproduction before the onset of castration. Here we propose this strategy as a novel defense against castrating pathogens in snails. Gigantism, a predicted outcome of castration due to energy reallocation, occurred early in infection (<15 weeks) and was not universal among the snail lines. Lines that did not show gigantism were also characterised by a high parasite production rate and low survivorship, perhaps indicating energy reallocation into parasite production and costly immune defense. We observed no differences in total parasite production among lines throughout the entire course of infection, although lines differed in their parasite reproductive rate. The average rate of parasite production varied among lines from 1300 to 2450 cercariae within a single 2h shedding period, resulting in a total production of 6981-29,509 cercariae over the lifetime of a single snail. Regardless of genetic background, snail size was a strong predictor of parasite reproduction: each millimetre increase in snail size at the time of the first shed resulted in up to 3500 more cercariae over the lifetime of the snail. The results of this study provide a detailed picture of variation in hosts' responses to infection and the resulting impacts on parasite fitness, further defining the intricacies of snail-schistosome compatibility.

Revision of Leptorhynchoides thecatus (Acanthocephala: Illiosentidae), with morphometric analysis and description of six new species

Six new species of Leptorhynchoides from the southeastern United States are described. These new species were once part of the Leptorhynchoides thecatus complex of species that was previously recognized on the basis of DNA sequence data. Multivariate morphometric analysis including discriminant function analysis and decision tree analysis indicated that each of the species is morphologically distinct. Both analyses classified more than 90% of specimens correctly and most misclassifications occurred between members of 2 pairs of species that are morphologically similar. The most discriminating continuous characters were: trunk length, number of longitudinal rows of hooks, length of the longest hook, and testes width. Hook asymmetry and missing hooks on the proboscis were also important taxonomic characters. The discriminant function and the decision tree generated from the data were used to classify new specimens, yielding a 96% and 84% correct classification rate, respectively. The new taxonomic designations account for much of the previously recognized variability in host use, habitat use, and development as determined by survey data. With the addition of these 6 new taxa, 10 species currently are recognized within the genus.

Monsters of the sea serpent: parasites of an oarfish, Regalecus russellii

Examination of a small portion of the viscera of an oarfish ( Regalecus russellii ) recovered from Santa Catalina Island, southern California, revealed numerous tetraphyllidean tapeworm plerocercoids, Clistobothrium cf. montaukensis; 2 juvenile nematodes, Contracaecum sp.; and a fragment of an adult acanthocephalan, family Arhythmacanthidae. This suggests that the fish was relatively heavily parasitized. The presence of larval and juvenile worms suggests that oarfish are preyed upon by deep-swimming predators such as the shortfin mako shark, Isurus oxyrinchus , known to be a definitive host for the adult tapeworm, and also by diving mammals such as sperm whales, Physeter catodon L., hosts of Contracaecum spp. nematodes.

Non-invasive sampling of schistosomes from humans requires correcting for family structure

For ethical and logistical reasons, population-genetic studies of parasites often rely on the non-invasive sampling of offspring shed from their definitive hosts. However, if the sampled offspring are naturally derived from a small number of parents, then the strong family structure can result in biased population-level estimates of genetic parameters, particularly if reproductive output is skewed. Here, we document and correct for the strong family structure present within schistosome offspring (miracidia) that were collected non-invasively from humans in western Kenya. By genotyping 2,424 miracidia from 12 patients at 12 microsatellite loci and using a sibship clustering program, we found that the samples contained large numbers of siblings. Furthermore, reproductive success of the breeding schistosomes was skewed, creating differential representation of each family in the offspring pool. After removing the family structure with an iterative jacknifing procedure, we demonstrated that the presence of relatives led to inflated estimates of genetic differentiation and linkage disequilibrium, and downwardly-biased estimates of inbreeding coefficients (F_IS). For example, correcting for family structure yielded estimates of F_ST among patients that were 27 times lower than estimates from the uncorrected samples. These biased estimates would cause one to draw false conclusions regarding these parameters in the adult population. We also found from our analyses that estimates of the number of full sibling families and other genetic parameters of samples of miracidia were highly intercorrelated but are not correlated with estimates of worm burden obtained via egg counting (Kato-Katz). Whether genetic methods or the traditional Kato-Katz estimator provide a better estimate of actual number of adult worms remains to be seen. This study illustrates that family structure must be explicitly accounted for when using offspring samples to estimate the genetic parameters of adult parasite populations.

Genetic epidemiology uses genetic data to uncover patterns of disease processes. To acquire data for these analyses, individual pathogens are collected and scored at genetic markers, and the resultant data are analyzed to infer biological patterns about the pathogen populations. In lieu of invasive sampling of adult pathogens in humans, researchers have relied on non-invasive sampling of parasite offspring (often shed in fecal samples). One potential problem with this approach is that analyses using the offspring data will be biased because many of the offspring are related and family sizes are likely to be unequal. We show that this sampling issue is relevant in a natural transmission zone in western Kenya and that it yields biases in three important parameters: genetic differentiation, inbreeding coefficients, and estimates of the amount of non-random association between loci (linkage disequilibrium). We also develop a method to remove these biases by removing the sibling structure present in the dataset. Finally, we suggest that our measure of family number, as well as other genetic measures, may be useful measures of the worm burdens in patients.

Host susceptibility is altered by light intensity after exposure to parasites

Translating research advances to natural systems using experimental laboratory studies is often difficult because of the variability between the natural environment and experimental conditions. Because environmental conditions have a large effect on an organism's physiology, responses to stressors like nutrient limitation, temperature, oxygen deprivation, predation, and parasite/pathogen infection are likely to be context dependent. Therefore, it is essential to examine the impact the study environment has on the experimental outcome. Here, we explored the effect of light exposure on susceptibility to parasite infection. The Biomphalaria glabrata / Schistosoma mansoni study system is a well-established model for studying schistosomiasis. It has been general practice to maintain the vector, B. glabrata, in dark conditions after exposure to miracidia of the human pathogen S. mansoni. We evaluated susceptibility of B. glabrata to S. mansoni under 3 different light conditions during the prepatent period, light (125 lx) on a 12-12 cycle, dim light (3 lx) on a 12-12 cycle, and no light (24 hr at 0 lx). We hypothesized that stress due to photoperiod disruption (24 hr of darkness) would result in compromised immune function and lead to higher susceptibility to infection. Prevalence of infected snails differed significantly between the light conditions, and higher susceptibility was observed in the full light and complete dark conditions compared with the low light conditions. The dim conditions are representative of current methods for evaluating susceptibility in this system. Our results indicate that light exposure during the prepatent period can affect infection outcomes, and environmental conditions must therefore be considered when assessing fitness and immune response due to interactions between host genotype and environment.

Polymorphism associated with the Schistosoma mansoni tetraspanin-2 gene

A vaccine against schistosomiasis would contribute significantly to reducing the 3-70 million disability-adjusted life years lost annually to the disease. Towards this end, inoculation with the large extracellular loop (EC-2) of Schistosoma mansoni tetraspanin-2 protein (Sm-TSP-2) has proved effective in reducing worm and egg burdens in S. mansoni-infected mice. The EC-2 loop of Schistosoma japonicum TSP-2, however, has been found to be highly polymorphic, perhaps diminishing the likelihood that this antigen can be used for vaccination against this species. Here, we examine polymorphism of the EC-2 of Sm-TSP-2 in genetically unique worms derived from six individuals from Kisumu, Kenya.

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.

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.

Reduced susceptibility to praziquantel among naturally occurring Kenyan isolates of Schistosoma mansoni

Background

The near exclusive use of praziquantel (PZQ) for treatment of human schistosomiasis has raised concerns about the possible emergence of drug-resistant schistosomes.

Methodology/Principal Findings

We measured susceptibility to PZQ of isolates of Schistosoma mansoni obtained from patients from Kisumu, Kenya continuously exposed to infection as a consequence of their occupations as car washers or sand harvesters. We used a) an in vitro assay with miracidia, b) an in vivo assay targeting adult worms in mice and c) an in vitro assay targeting adult schistosomes perfused from mice. In the miracidia assay, in which miracidia from human patients were exposed to PZQ in vitro, reduced susceptibility was associated with previous treatment of the patient with PZQ. One isolate (“KCW”) that was less susceptible to PZQ and had been derived from a patient who had never fully cured despite multiple treatments was studied further. In an in vivo assay of adult worms, the KCW isolate was significantly less susceptible to PZQ than two other isolates from natural infections in Kenya and two lab-reared strains of S. mansoni. The in vitro adult assay, based on measuring length changes of adults following exposure to and recovery from PZQ, confirmed that the KCW isolate was less susceptible to PZQ than the other isolates tested. A sub-isolate of KCW maintained separately and tested after three years was susceptible to PZQ, indicative that the trait of reduced sensitivity could be lost if selection was not maintained.

Conclusions/Significance

Isolates of S. mansoni from some patients in Kisumu have lower susceptibility to PZQ, including one from a patient who was never fully cured after repeated rounds of treatment administered over several years. As use of PZQ continues, continued selection for worms with diminished susceptibility is possible, and the probability of emergence of resistance will increase as large reservoirs of untreated worms diminish. The potential for rapid emergence of resistance should be an important consideration of treatment programs.

The emergence of drug resistant pathogens is a great challenge to the control of infectious diseases. Schistosomiasis is one of the world's greatest neglected tropical diseases, and it is primarily controlled with the drug praziquantel. This drug is often used by repeatedly treating patients to maintain reduced worm burdens, an ideal situation to encourage the evolution of resistant worms. Although drug based control programs are increasing, monitoring efforts for drug resistance remain rare. We measured drug susceptibility of schistosomes from a cohort of patients in Kenya who are enrolled in a longitudinal study in which they are repeatedly treated with praziquantel. We found that schistosomes from previously treated patients were significantly less susceptible than those that were not. Also, schistosomes derived from a single patient who had been treated with praziquantel 18 times showed marked resistance. Although the findings of this study indicated that reduced drug susceptibility occurs in this population of schistosomes, this trait does not seem to be spreading widely or creating clinical levels of resistance. We hypothesize that the trait remains at low frequency because of the large population of schistosomes that are not exposed to the drug and/or potential fitness costs associated with reduced susceptibility.

The sex lives of parasites: investigating the mating system and mechanisms of sexual selection of the human pathogen Schistosoma mansoni

The mating systems of internal parasites are inherently difficult to investigate although they have important implications for the evolutionary biology of the species, disease epidemiology, and are important considerations for control measures. Using parentage analyses, three topics concerning the mating biology of Schistosoma mansoni were investigated: the number of mates per adult male and female, variance in reproductive success among individuals, and the potential role for sexual selection on male body size and also mate choice for genetically dissimilar individuals. Results indicated that schistosomes were mostly monogamous, and evidence of only one mate change occurred over a period of 5-6 weeks. One male was polygynous and contained two females in its gynecophoral canal although offspring were only detected for one of the females. Even though they were primarily monogamous and the sex ratio near even, reproductive success was highly variable, indicating a potential role for sexual selection. Male body size was positively related to reproductive success, consistent with sexual selection via male-male competition and female choice for large males. However, relatedness of pairs was not associated with their reproductive success. Finally, genetically identical individuals differed significantly in their reproductive output and identical males in their body size, indicating important partner and environmental effects on these traits.

Genetic structure of Schistosoma mansoni in western Kenya: The effects of geography and host sharing

We examined the spatial structure of Schistosoma mansoni, a parasite of humans, from natural infections at two levels: across the Lake Victoria basin of Kenya and among snail hosts. Using 20 microsatellite markers we examined geographic patterns of relatedness and population structure of cercariae and found weak, but significant structure detected by some, but not all analyses. We hypothesise structure created by aggregations of clonal individuals or adherence of hosts to local transmission sites is eroded by high amounts of gene flow in the region. This finding also supports previous hypotheses concerning the evolution of drug resistance in the region. Intrasnail dynamics were investigated in the context of aggregation and kin selection theory to determine how relatedness and also sex influence host sharing and host exploitation. Cercarial production did not differ significantly between snails with one or two genotypes suggesting that mixed infections resulted in decreased individual fitness and provides a framework for reproductive competition. Coinfection patterns in snails were independent of parasite relatedness indicating that schistosomes were not aggregated according to their relatedness and that kin selection was not influencing host sharing. Additionally, host exploitation in coinfections (measured by cercarial production) was not negatively correlated with relatedness, as predicted by classical models due to increased competition and thus exploitation when parasites are unrelated. Because of the low levels of relatedness within the population, schistosomes may rarely encounter close relatives and kin selection mechanisms that influence the distribution of individuals within snails or the virulence mode of the parasites may simply have not evolved.

Impact of intense, longitudinal retreatment with praziquantel on cure rates of schistosomiasis mansoni in a cohort of occupationally exposed adults in western Kenya

OBJECTIVE

To investigate trends in the efficacy of praziquantel (PZQ) suggestive of the emergence of drug resistance against Schistosoma mansoni infection after 12.5 years of intense, repeated use in a small geographic area along the shores of Lake Victoria.

METHODS

As part of a longitudinal study, 178 men occupationally exposed to schistosomes were repeatedly tested for S. mansoni infection at 4- to 6-week intervals and treated with PZQ at each reinfection. We compared cure rates by year of study and examined factors associated with cure failure in a multivariate logistic regression model.

RESULTS

Overall, the cure rate after a single dose of PZQ was 66%, ranging annually from 36% to 82%. In multivariate analysis, failure to cure after 1 PZQ dose was significantly associated with high intensity of infection and having fewer previous exposures to dying worms. Even after adjustment for these factors, treatments administered in 2006 were significantly more likely to result in cure failures than treatments administered in 2004, the year in which PZQ efficacy was highest. While cure rates varied over the course of 12 years, there was no consistent downward trend towards decreased efficacy over time. In years for which malacological data were available, periods of low PZQ efficacy coincide with high rates of S. mansoni infection in nearby snail populations.

CONCLUSION

We did not find a pattern of cure failures consistent with development of clinical resistance to PZQ in our intensely treated cohort.

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.

Introgressive hybridization of human and rodent schistosome parasites in western Kenya

Hybridization and introgression can have important consequences for the evolution, ecology and epidemiology of pathogenic organisms. We examined the dynamics of hybridization between a trematode parasite of humans, Schistosoma mansoni, and its sister species, S. rodhaini, a rodent parasite, in a natural hybrid zone in western Kenya. Using microsatellite markers, rDNA and mtDNA, we showed that hybrids between the two species occur in nature, are fertile and produce viable offspring through backcrosses with S. mansoni. Averaged across collection sites, individuals of hybrid ancestry comprised 7.2% of all schistosomes collected, which is a large proportion given that one of the parental species, S. rodhaini, comprised only 9.1% of the specimens. No F1 individuals were collected and all hybrids represented backcrosses with S. mansoni that were of the first or successive generations. The direction of introgression appears highly asymmetric, causing unidirectional gene flow from the rodent parasite, S. rodhaini, to the human parasite, S. mansoni. Hybrid occurrence was seasonal and most hybrids were collected during the month of September over a 2-year period, a time when S. rodhaini was also abundant. We also examined the sex ratios and phenotypic differences between the hybrids and parental species, including the number of infective stages produced in the snail host and the time of day the infective stages emerge. No statistical differences were found in any of these characteristics, and most of the hybrids showed an emergence pattern similar to that of S. mansoni. One individual, however, showed a bimodal emergence pattern that was characteristic of both parental species. In conclusion, these species maintain their identity despite hybridization, although introgression may cause important alterations of the biology and epidemiology of schistosomiasis in this region.

Interactions between natural populations of human and rodent schistosomes in the Lake Victoria region of Kenya: a molecular epidemiological approach

BACKGROUND

Schistosoma mansoni exists in a complex environmental milieu that may select for significant evolutionary changes in this species. In Kenya, the sympatric distribution of S. mansoni with S. rodhaini potentially influences the epidemiology, ecology, and evolutionary biology of both species, because they infect the same species of snail and mammalian hosts and are capable of hybridization.

METHODOLOGY/PRINCIPAL FINDINGS

Over a 2-year period, using a molecular epidemiological approach, we examined spatial and temporal distributions, and the overlap of these schistosomes within snails, in natural settings in Kenya. Both species had spatially and temporally patchy distributions, although S. mansoni was eight times more common than S. rodhaini. Both species were overdispersed within snails, and most snails (85.2% for S. mansoni and 91.7% for S. rodhaini) only harbored one schistosome genotype. Over time, half of snails infected with multiple genotypes showed a replacement pattern in which an initially dominant genotype was less represented in later replicates. The other half showed a consistent pattern over time; however, the ratio of each genotype was skewed. Profiles of circadian emergence of cercariae revealed that S. rodhaini emerges throughout the 24-hour cycle, with peak emergence before sunrise and sometimes immediately after sunset, which differs from previous reports of a single nocturnal peak immediately after sunset. Peak emergence for S. mansoni cercariae occurred as light became most intense and overlapped temporally with S. rodhaini. Comparison of schistosome communities within snails against a null model indicated that the community was structured and that coinfections were more common than expected by chance. In mixed infections, cercarial emergence over 24 hours remained similar to single species infections, again with S. rodhaini and S. mansoni cercarial emergence profiles overlapping substantially.

CONCLUSIONS/SIGNIFICANCE

The data from this study indicate a lack of obvious spatial or temporal isolating mechanisms to prevent hybridization, raising the intriguing question of how the two species retain their separate identities.

Cryptic speciation and patterns of phenotypic variation of a highly variable acanthocephalan parasite

An investigation of a parasite species that is broadly host- and habitat-specific and exhibits alternative transmission strategies was undertaken to examine intraspecific variability and if it can be attributed to cryptic speciation or environmentally induced plasticity. Specimens of an acanthocephalan parasite, Leptorhynchoides thecatus, collected throughout North America were analysed phylogenetically using sequences of the cytochrome oxidase I gene and the internal transcribed spacer region. Variation in host use, habitat use, and transmission were examined in a phylogenetic context to determine if they were more likely phylogenetically based or due to environmental influences. Results indicated that most of the variation detected can be explained by the presence of cryptic species. The majority of these species have narrow host and microhabitat specificities although one species, which also may comprise a complex of species, exhibits broad host and habitat specificity. Alternate transmission pathways only occurred in two of the cryptic species and correlate with host use patterns. Taxa that mature in piscivorous piscine hosts use a paratenic fish host to bridge the trophic gap between their amphipod intermediate host and piscivorous definitive host. One potential example of environmentally induced variation was identified in three populations of these parasites, which differ on their abilities to infect different host species.

GEOGRAPHIC ANALYSIS OF HOST USE, DEVELOPMENT, AND HABITAT USE OF AN ACANTHOCEPHALAN SPECIES, LEPTORHYNCHOIDES THECATUS

Leptorhynchoides thecatus (Linton, 1891), an acanthocephalan parasite of freshwater fishes, varies in host use, development, and habitat use throughout North America. Spatial structure of these characteristics was examined from data extracted from the literature. Geographic patterns were inferred from point comparisons using correllograms and then tested with Moran's I statistic for global and local significance, and visually from regional means within major river drainages. Species of Micropterus Lacepède, 1802 (black basses) were common hosts in most regions, except the Lower Mississippi and South Atlantic regions where species of Lepomis Rafinesque, 1819 (sunfishes) were common hosts. Development, described as the proportions of adults relative to cystacanths (extraintestinal juveniles), decreased with latitude. Habitat use of L. thecatus showed marked geographic patterns. Leptorhynchoides thecatus occurred in the intestine of sunfishes in the South Atlantic and Lower Mississippi regions, in the ceca in fish of all species included in the study in the Missouri and Texas-Gulf regions, and both in ceca and intestines in fish of all species in northern regions. Leptorhynchoides thecatus showed geographic patterning within the variable traits across the range of the species. These patterns may be the result of ecological factors or of genetic differences that might indicate L. thecatus comprises multiple cryptic species.

First Sequenced Mitochondrial Genome from the Phylum Acanthocephala(Leptorhynchoides thecatus) and Its Phylogenetic Position Within Metazoa

The complete sequence of the mitochondrial genome of Leptorhynchoides thecatus (Acanthocephala) was determined, and a phylogenetic analysis was carried out to determine its placement within Metazoa. The genome is circular, 13,888 bp, and contains at least 36 of the 37 genes typically found in animal mitochondrial genomes. The genes for the large and small ribosomal RNA subunits are shorter than those of most metazoans, and the structures of most of the tRNA genes are atypical. There are two significant noncoding regions (377 and 294 bp), which are the best candidates for a control region; however, these regions do not appear similar to any of the control regions of other animals studied to date. The amino acid and nucleotide sequences of the protein coding genes of L. thecatus and 25 other metazoan taxa were used in both maximum likelihood and maximum parsimony phylogenetic analyses. Results indicate that among taxa with available mitochondrial genome sequences, Platyhelminthes is the closest relative to L. thecatus, which together are the sister taxon of Nematoda; however, long branches and/or base composition bias could be responsible for this result. The monophyly of Ecdysozoa, molting organisms, was not supported by any of the analyses. This study represents the first mitochondrial genome of an acanthocephalan to be sequenced and will allow further studies of systematics, population genetics, and genome evolution.

Effect of cystacanth body size on adult success

Laboratory-reared cystacanths of Leptorhynchoides thecatus (Acanthocephala: Rhadinorhynchidae) were used to study the effect of cystacanth size on adult success and the factors that influence cystacanth size within the intermediate host. To assess how host size and intensity of infection influence cystacanth size, infected amphipods (Hyalella azteca) were measured, and sex, length, and width of cystacanths were determined. After a subset of cystacanths was measured, small- and large-size classes of cystacanths were designated. To determine how cystacanth size relates to adult size, green sunfish (Lepomis cyanellus) were fed 10 large or small cystacanths. Fish were dissected 6 wk after infection, and worms were removed. After adult worms were permanently mounted on slides, their length and width were measured. Intensity of infection and amphipod size significantly influenced cystacanth size in that large amphipods harbored larger cystacanths than did small amphipods and heavy infections produced smaller cystacanths than did light infections. Adult worms from the small and large cystacanth-size classes showed no significant difference in size; however, large cystacanths had a significantly higher establishment and survival than did small cystacanths: 40% of large worms and 14% of small worms were recovered. The results of this study indicate that host size and host sharing influence cystacanth size and that cystacanth size is an important factor in determining adult success.

Seasonal dynamics of the helminths of bluegill (Lepomis macrochirus) in a subtropical region

Seasonal prevalence and abundance of the helminths of bluegill sunfish in a Louisiana oligohaline bayou were measured by a survey and a field "live-box" experiment. The survey took place from the spring of 1997 to the summer of 1998 and examined fish that were <7 cm. Three trematode species, Phagicola nana (Heterophyidae), Ascocotyle tenuicollis (Heterophyidae), and Posthodiplostomum minimum (Diplostomatidae), and 2 nematode species, Camallanus oxycephalus (Camallanidae) and Spinitectus carolini (Cystidicolidae), were examined. Camallanus oxycephalus was the only helminth that showed a distinct seasonal pattern. Abundance and prevalence peaked in summer, which was likely driven by concerted reproductive cycles of females. The survey data indicated that the other helminths fluctuated over time but did not seem to follow a distinct seasonal pattern. The fluctuations could be attributed to the variable nature of the habitat or to the stochastic events that influence transmission dynamics, particularly isolated events such as tropical storms, which cause dramatic changes in salinity. The live-box experiment successfully measured recruitment of A. tenauicollis and suggested a seasonal component of the distribution of this species, which was possibly caused by temperature-dependent emergence of cercariae from the snail intermediate host.