Publications

Pooled sequencing provides a rapid cost-effective approach to assess genetic variation segregating within populations of organisms. However, such studies are typically limited to single nucleotide variants and small indels (≤ 50bp), and have not been used for structural variants (SVs; >50bp) which impact large portions of most genomes and may significantly impact phenotype. Here, we examined SVs circulating in five laboratory populations of the human parasite Schistosoma mansoni by generating long-read sequences from pools of worms (92 -152 per population). We were able identify and genotype 17,446 SVs, representing 6.5% of the genome despite challenges in identifying low frequency variants. SVs included deletions (n=8,525), duplications (n=131), insertions (n=8,410), inversions (n=311), and translocations (n=69) and were enriched in repeat regions. More than half (59%) of the SVs were shared between ≥4 populations, but 12% were found in only one of the five populations. Within this subset, we identified 168 population-specific SVs that were at-or-near fixation (>95% alternate allele frequency) in one population but missing (<5%) in the other four populations. Five of these variants impact the coding sequence of 6 genes. We also identified 8 SVs with extreme allele frequency differences between populations within quantitative trait loci for biomedically important pathogen phenotypes (drug resistance, larval stage production) identified in prior genetic mapping studies. These results demonstrate that long-read sequence data from pooled individuals is a viable method to quickly catalogue SVs circulating within populations. Furthermore, some of these variants may be responsible for, or linked to, regions experiencing, population-specific directional selection.

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.