Publications

The human parasitic fluke, Schistosoma haematobium hybridizes with the livestock parasite S. bovis in the laboratory, but the frequency of hybridization in nature is unclear. We analyzed 34.6 million single nucleotide variants in 162 samples from 18 African countries, revealing a sharp genetic discontinuity between northern and southern S. haematobium. We found no evidence for recent hybridization. Instead the data reveal admixture events that occurred 257-879 generations ago in northern S. haematobium populations. Fifteen introgressed S. bovis genes are approaching fixation in northern S. haematobium with four genes potentially driving adaptation. We identified 19 regions that were resistant to introgression; these were enriched on the sex chromosomes. These results (i) suggest strong barriers to gene flow between these species, (ii) indicate that hybridization may be less common than currently envisaged, but (iii) reveal profound genomic consequences of rare interspecific hybridization between schistosomes of medical and veterinary importance.

There are limited control measures for the disease schistosomiasis, despite the fact that infection with parasitic blood flukes affects hundreds of millions of people worldwide. The current treatment, praziquantel, has been in use since the 1980’s and there is a concern that drug resistance may emerge with continued monotherapy. Given the need for additional antischistosomal drugs, we have re-visited an old lead, meclonazepam. In comparison to praziquantel, there has been relatively little work on its antiparasitic mechanism. Recent findings indicate that praziquantel and meclonazepam act through distinct receptors, making benzodiazepines a promising chemical series for further exploration. Previous work has profiled the transcriptional changes evoked by praziquantel treatment. Here, we examine in detail schistosome phenotypes evoked by in vitro and in vivo meclonazepam treatment. These data confirm that meclonazepam causes extensive tegument damage and directly kills parasites, as measured by pro-apoptotic caspase activation. In vivo meclonazepam exposure results in differential expression of many genes that are divergent in parasitic flatworms, as well as several gene products implicated in blood feeding and regulation of hemostasis in other parasites. Many of these transcripts are also differentially expressed with praziquantel exposure, which may reflect a common schistosome response to the two drugs. However, despite these similarities in drug response, praziquantel-resistant parasites retain susceptibility to meclonazepam’s schistocidal effects. These data provide new insight into the mechanism of antischistosomal benzodiazepines, resolving similarities and differences with the current frontline therapy, praziquantel.

Pathogen genomics is a powerful tool for tracking infectious disease transmission. In malaria, identity-by-descent is used to assess the genetic relatedness between parasites and has been used to study transmission and importation. In theory, identity-by-descent can be used to distinguish genealogical relationships to reconstruct transmission history or identify parasites for QTL experiments. MalKinID (Malaria Kinship Identifier) is a new classification model designed to identify genealogical relationships among malaria parasites based on genome-wide identity-by-descent proportions and identity-by-descent segment distributions. MalKinID was calibrated to the genomic data from 3 laboratory-based genetic crosses (yielding 440 parent-child and 9060 full-sibling comparisons). MalKinID identified lab-generated F1 progeny with >80% sensitivity and showed that 0.39 (95% CI 0.28, 0.49) of the second-generation progeny of a NF54 and NHP4026 cross were F1s and 0.56 (0.45, 0.67) were backcrosses of an F1 with the parental NF54 strain. In simulated outcrossed importations, MalKinID reconstructs genealogy history with high precision and sensitivity, with F1-scores exceeding 0.84. However, when importation involves inbreeding, such as during serial co-transmission, the precision and sensitivity of MalKinID declined, with F1-scores (the harmonic mean of precision and sensitivity) of 0.76 (0.56, 0.92) and 0.23 (0.0, 0.4) for parent-child and full-sibling and <0.05 for second-degree and third-degree relatives. Disentangling inbred relationships required adapting MalKinID to perform multisample comparisons. Genealogical inference is most powered when (1) outcrossing is the norm or (2) multisample comparisons based on a predefined pedigree are used. MalKinID lays the foundations for using identity-by-descent to track parasite transmission history and for separating progeny for quantitative-trait-locus experiments.

Background: The microbiome is increasingly recognized to shape many aspects of its host biology and is a key determinant of health and disease. The microbiome may influence transmission of pathogens by their vectors, such as mosquitoes or aquatic snails. We previously sequenced the V4 region of the bacterial 16S rRNA gene from the hemolymph (blood) of Biomphalaria spp. snails, vectors of the human blood fluke schistosome. We showed that snail hemolymph harbored an abundant and diverse microbiome. This microbiome is distinct from the water environment and can discriminate snail species and populations. As hemolymph bathes snail organs, we then investigated the heterogeneity of the microbiome in these organs.

Results: We dissected ten snails for each of two different species (B. alexandrina and B. glabrata) and collected their hemolymph and organs (ovotestis, hepatopancreas, gut, and stomach). We also ground in liquid nitrogen four whole snails of each species. We sampled the water in which the snails were living (environmental controls). Sequencing the 16S rRNA gene revealed organ-specific microbiomes. These microbiomes harbored a lower diversity than the hemolymph microbiome, and the whole-snail microbiome. The organ microbiomes tend to cluster by physiological function. In addition, we showed that the whole-snail microbiome is more similar to hemolymph microbiome.

Conclusions: These results are critical for future work on snail microbiomes and show the necessity of sampling individual organ microbiomes to provide a complete description of snail microbiomes.

The microbiome is increasingly recognized to shape many aspects of its host biology and is a key determinant of health and disease. The microbiome may influence transmission of pathogens by their vectors, such as mosquitoes or aquatic snails. We previously sequenced the V4 region of the bacterial 16S rRNA gene from the hemolymph (blood) of Biomphalaria spp. snails, vectors of the human blood fluke schistosome. We showed that snail hemolymph harbored an abundant and diverse microbiome. This microbiome is distinct from the water environment and can discriminate snail species and populations. As hemolymph bathes snail organs, we then investigated the heterogeneity of the microbiome in these organs.

We dissected ten snails for each of two different species (B. alexandrina and B. glabrata) and collected their hemolymph and organs (ovotestis, hepatopancreas, gut, and stomach). We also ground in liquid nitrogen four whole snails of each species. We sampled the water in which the snails were living (environmental controls). Sequencing the 16S rRNA gene revealed organ-specific microbiomes. These microbiomes harbored a lower diversity than the hemolymph microbiome, and the whole-snail microbiome. The organ microbiomes tend to cluster by physiological function. In addition, we showed that the whole-snail microbiome is more similar to hemolymph microbiome.

These results are critical for future work on snail microbiomes and show the necessity of sampling individual organ microbiomes to provide a complete description of snail microbiomes.

The anthelmintic praziquantel (PZQ) has been used for decades as the clinical therapy for schistosomiasis, and remains the only available drug. As a cheap and effective drug therapy for all human disease-causing Schistosoma species, usage of PZQ underpins mass drug administration strategies aimed at eliminating schistosomiasis as a public health problem by 2030. Concern over the potential emergence of resistance to PZQ is therefore warranted, as it would constitute a major threat to this approach. In terms of molecular adaptations conferring PZQ resistance, variation in the sequence and/or expression of the drug target is an obvious mechanism and should be a priority for surveillance efforts. The target of PZQ is a transient receptor potential ion channel, TRPM PZQ , which is established as a locus that regulates schistosome sensitivity to PZQ. Here, we describe the establishment of a community resource, 'TRPtracker', which coalesces data on TRPM PZQ natural variants together with measurements of individual variant sensitivity to PZQ. A compendium of laboratory-generated mutants in TRPM PZQ is also compiled in TRPtracker to map regions within TRPM PZQ critical for PZQ sensitivity. Aggregation of data from multiple research groups into TRPtracker permits rapid community-wide exchange of data, cataloguing which TRPM PZQ variants have been functionally profiled, where geographically these variants have been found, their frequency within populations and their potential impact on PZQ sensitivity.

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 carried on average 62% (π = 1.52e-04 - 7.15e-04) less nucleotide diversity than the four field parasite populations (π = 9.06e-03 - 2.24e-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 - 271 kb) compared to field populations (9 bp - 17.1 kb). SmBRE was the least variable laboratory population; 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.

BACKGROUND: Artemisinin partial resistance (ART-R) has spread throughout Southeast Asia and mutations in Pfkelch13, the molecular marker of resistance, are widely reported in East Africa. Effective in vitro assays and robust phenotypes are crucial for monitoring populations for the emergence and spread of resistance. The recently developed extended Recovery Ring-stage Survival Assay used a qPCR-based readout to reduce the labour intensiveness for in vitro phenotyping of ART-R and improved correlation with the clinical phenotype of ART-R. Here, the assay is extended and refined to include measurements of parasite growth and recovery after drug exposure. Clinical isolates and progeny from two genetic crosses were used to optimize and validate the reliability of a straight-from-blood, SYBR Green-based qPCR protocol in a 96-well plate format to accurately measure phenotypes with this new Growth, Resistance, and Recovery assay (GRRA).

RESULTS: The assay determined growth between 6 and 96 h, resistance at 120 h, and recovery from 120 to 192 h. Growth can be accurately captured by qPCR and is shown by reproduction of previous growth phenotypes from HB3 × Dd2. Resistance measured at 120 h continually shows the most consistent phenotype for ring stage susceptibility. Recovery identifies an additional response to drug in parasites that are determined sensitive by replicative viability at 120 h. Comparison of progeny phenotypes for Growth versus Resistance showed a minor but significant correlation, whereas Growth versus Recovery and Resistance versus Recovery showed no significant correlation. Additionally, dried blood spot (DBS) samples matched replicative viability measured from liquid samples demonstrating Resistance can be easily quantified using either storage method.

CONCLUSIONS: The direct-from-blood qPCR-based methodology provides the throughput needed to quickly measure large numbers of parasites for multiple relevant phenotypes. Growth can reveal fitness defects and illuminate relationships between proliferation rates and drug response. Recovery serves as a complementary phenotype to resistance that quantifies the ability of sensitive parasites to tolerate drug exposure. All three phenotypes offer a comprehensive assessment of parasite-drug interaction each with potential independent genetic determinants of main effect and overlapping secondary effects. By adapting the method to include DBS, readouts can be easily extended to ex vivo surveillance applications.

Avian Influenza viruses (AIVs) present a public health risk, especially with seasonal vaccines offering limited protection. AIV H5N1 clade 2.3.4.4b has caused a multi-state outbreaks in the United States (US) poultry and cattle since March 2024, raising pandemic concerns. We developed a nonstructural protein 1 (NS1)-deficient mutant of a low pathogenic version of the cattle-origin human influenza A/Texas/37/2024 H5N1, namely LPhTXdNS1, and assessed its safety, immunogenicity, and protection efficacy. LPhTXdNS1 is attenuated in vitro, showing reduced replication efficiency in Vero cells and inability to control IFNβ promoter activation. The LPhTXdNS1-immunized C57BL/6 J mice exhibit significantly reduced viral replication and pathogenicity compared to those infected with the low pathogenic version expressing NS1, namely LPhTX. Notably, a single intranasal dose of LPhTXdNS1 elicited protective immune responses, providing robust protection against lethal wild-type H5N1 challenge. These results demonstrate that LPhTXdNS1 is safe and able to induce protective immune responses against H5N1.

Genomic analysis of parasites can deepen our understanding of their transmission, population structure, and important biological characteristics. Onchocerciasis (river blindness), caused by the parasitic nematode Onchocerca volvulus, involves adult worms residing in subcutaneous nodules that produce larval-stage microfilariae (mf), which are routinely detected in the skin for diagnosis. Whole-genome studies of mf are limited; most analyses have focused on the mitochondrial genome. We conducted a genome-wide analysis with 94% median nuclear genome coverage, analyzing 171, 37, and 98 mf from 16, 3, and 5 individuals from Ghana, Liberia, and the Democratic Republic of Congo, respectively. These data were used to investigate population differentiation, estimate the number of reproductive adult worms, and analyze genetic variation across chromosomes. Population genetic analyses across hosts and countries showed that nuclear genome diversity can reveal fine-scale genetic structure, even between geographically close countries, providing more resolution than mitochondrial haplotype data. By reconstructing maternal and paternal sibships, we estimated the number of reproductively active adult filariae. Comparisons between adult worm estimates from genetic data and nodule observations showed that genetics-based estimates were higher or equal to observed worm counts in 8 out of 9 hosts for female worms and 7 out of 9 hosts for male worms. Our analysis also revealed lower-than-expected X chromosome diversity, consistent with neo-X chromosome fusions in filarial species. This study represents an important step in using nuclear genome data from mf to support onchocerciasis elimination efforts and in developing genetic tools that could inform mass drug administration programs.