FITNESS OF INDIRECTLY TRANSMITTED PATHOGENS: RESTRAINT AND CONSTRAINT

Abundant genetic variation is retained in many laboratory schistosome populations

Schistosomes are obligately sexual blood flukes that can be maintained in the laboratory using freshwater snails as intermediate and rodents as definitive hosts. The genetic composition of laboratory schistosome populations is poorly understood: whether genetic variation has been purged due to serial inbreeding or retained is unclear. We sequenced 19 - 24 parasites from each of five laboratory Schistosoma mansoni populations and compared their genomes with published exome data from four S. mansoni field populations. We found abundant genomic variation (0.897 - 1.22 million variants) within laboratory populations: these retained on average 49% (π = 3.27e-04 - 8.94e-04) of the nucleotide diversity observed in the four field parasite populations (π = 1.08e-03 - 2.2e-03). However, the pattern of variation was very different in laboratory and field populations. Tajima's D was positive in all laboratory populations except SmBRE, indicative of recent population bottlenecks, but negative in all field populations. Current effective population size estimates of laboratory populations were lower (2 - 258) compared to field populations (3,174 - infinity). The distance between markers at which linkage disequilibrium (LD) decayed to 0.5 was longer in laboratory populations (59 bp - 180 kb) compared to field populations (9 bp - 9.5 kb). SmBRE was the least variable; this parasite also shows low fitness across the lifecycle, consistent with inbreeding depression. The abundant genetic variation present in most laboratory schistosome populations has several important implications: (i) measurement of parasite phenotypes, such as drug resistance, using laboratory parasite populations will determine average values and underestimate trait variation; (ii) genome-wide association studies (GWAS) can be conducted in laboratory schistosome populations by measuring phenotypes and genotypes of individual worms; (iii) genetic drift may lead to divergence in schistosome populations maintained in different laboratories. We conclude that the abundant genetic variation retained within many laboratory schistosome populations can provide valuable, untapped opportunities for schistosome research.

Rapid phenotypic and genotypic change in a laboratory schistosome population

Background

Genomic analysis has revealed extensive contamination among laboratory-maintained microbes including malaria parasites, Mycobacterium tuberculosis and Salmonella spp. Here, we provide direct evidence for recent contamination of a laboratory schistosome parasite population, and we investigate its genomic consequences. The Brazilian Schistosoma mansoni population SmBRE has several distinctive phenotypes, showing poor infectivity, reduced sporocysts number, low levels of cercarial shedding and low virulence in the intermediate snail host, and low worm burden and low fecundity in the vertebrate rodent host. In 2021 we observed a rapid change in SmBRE parasite phenotypes, with a ~10x increase in cercarial production and ~4x increase in worm burden.

Methods

To determine the underlying genomic cause of these changes, we sequenced pools of SmBRE adults collected during parasite maintenance between 2015 and 2023. We also sequenced another parasite population (SmLE) maintained alongside SmBRE without phenotypic changes.

Results

While SmLE allele frequencies remained stable over the eight-year period, we observed sudden changes in allele frequency across the genome in SmBRE between July 2021 and February 2023, consistent with expectations of laboratory contamination. (i) SmLE-specific alleles rose in the SmBRE population from 0 to 41-46% across the genome between September and October 2021, documenting the timing and magnitude of the contamination event. (ii) After contamination, strong selection (s = ~0.23) drove replacement of low fitness SmBRE with high fitness SmLE alleles. (iii) Allele frequency changed rapidly across the whole genome, except for a region on chromosome 4 where SmBRE alleles remained at high frequency.

Conclusions

We were able to detect contamination in this case because SmBRE shows distinctive phenotypes. However, this would likely have been missed with phenotypically similar parasites. These results provide a cautionary tale about the importance of tracking the identity of parasite populations, but also showcase a simple approach to monitor changes within populations using molecular profiling of pooled population samples to characterize fixed single nucleotide polymorphisms. We also show that genetic drift results in continuous change even in the absence of contamination, causing parasites maintained in different labs (or sampled from the same lab at different times) to diverge.

Genetic and molecular basis of drug resistance and species-specific drug action in schistosome parasites

Oxamniquine resistance evolved in the human blood fluke (Schistosoma mansoni) in Brazil in the 1970s. We crossed parental parasites differing ~500-fold in drug response, determined drug sensitivity and marker segregation in clonally derived second-generation progeny, and identified a single quantitative trait locus (logarithm of odds = 31) on chromosome 6. A sulfotransferase was identified as the causative gene by using RNA interference knockdown and biochemical complementation assays, and we subsequently demonstrated independent origins of loss-of-function mutations in field-derived and laboratory-selected resistant parasites. These results demonstrate the utility of linkage mapping in a human helminth parasite, while crystallographic analyses of protein-drug interactions illuminate the mode of drug action and provide a framework for rational design of oxamniquine derivatives that kill both S. mansoni and S. haematobium, the two species responsible for >99% of schistosomiasis cases worldwide.

Stage-specific differences in fecundity over the life-cycle of two characterized isolates of the liver fluke, Fasciola hepatica

The variability inherent in different isolates of Fasciola hepatica has been evident from reports in the literature but to date there has been no systematic examination of the relationship between these differences and the fecundity of the parasite. In this study we have attempted to remedy this situation by comparing the relative efficiencies with which 2 well-characterized isolates of the liver fluke (Oberon and Fairhurst) progress through both their definitive and intermediate hosts. We did not observe a reduction in fitness in the Oberon isolate which has been reported to be triclabendazole-resistant, compared to the triclabendazole-susceptible Fairhurst isolate, but considerable inter- and intra-isolate variability at different life-cycle stages was recorded. Thus the Oberon isolate gave 4-fold the number of cercariae when 100 snails were each challenged with a single miracidium and was more successful in establishing productive infections in rats. Fairhurst metacercariae excysted at a higher rate than those from the Oberon isolate and Fairhurst flukes produced 4-fold more eggs. The extent of the intra- and inter-isolate variability revealed in this work will provide a basis for the development of models of population dynamics aimed at predicting the response of the liver fluke to changing environmental conditions such as the use of anthelmintics or climatic change.