Search

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.

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., Competing Interests: Declarations Competing interests The authors declare that they have no competing interests.

The trematodes that cause schistosomiasis in humans require aquatic snails as intermediate hosts. Identifying the genes in snails at which allelic variation controls resistance to infection by schistosomes could lead to novel ways to break the cycle of transmission. We therefore mapped genetic variation within the BS90 population of Biomphalaria glabrata snails that controls their resistance to infection by the SmLE population of Schistosoma mansoni. A marker in the PTC2 genomic region strongly associates with variation in resistance. The S-haplotype, which confers increased susceptibility, appears to be almost completely dominant to the R-haplotype, which confers increased resistance. This result suggests a model in which the parasite must match a molecule on the host side to successfully infect. The genomic region surrounding our marker shows high structural and sequence variability between haplotypes. It is also highly enriched for genes that code for single-pass transmembrane (TM1) genes. Several of the TM1 genes present on the S-haplotype lack orthologs on the R-haplotype, which makes them intriguing candidate genes in a model of dominant susceptibility. These results add to a growing body of work that suggests TM1 genes, especially those in this exceptionally diverse genomic region, may play an important role in snail-schistosome compatibility polymorphisms., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Blouin et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The human parasitic fluke, Schistosoma haematobium hybridizes with the livestock parasite S. bovis in the laboratory, but the extent 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) demonstrate 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 interspecific hybridization between schistosomes of medical and veterinary importance., Competing Interests: Competing interests The authors declare no competing interests.

Piperaquine (PPQ) is widely used in combination with dihydroartemisinin as a first-line treatment against malaria. Multiple genetic drivers of PPQ resistance have been reported, including mutations in the Plasmodium falciparum chloroquine resistance transporter ( pfcrt ) and increased copies of plasmepsin II/III ( pm2/3 ). We generated a cross between a Cambodia-derived multidrug-resistant KEL1/PLA1 lineage isolate (KH004) and a drug-susceptible Malawian parasite (Mal31). Mal31 harbors a wild-type (3D7-like) pfcrt allele and a single copy of pm2/3 , while KH004 has a chloroquine-resistant (Dd2-like) pfcrt allele with an additional G367C substitution and multiple copies of pm2/3 . We recovered 104 unique recombinant parasites and examined a targeted set of progeny representing all possible combinations of variants at pfcrt and pm2/3 . We performed a detailed analysis of competitive fitness and a range of PPQ susceptibility phenotypes with these progenies, including PPQ survival assay, area under the dose response curve, and a limited point IC 50 . We find that inheritance of the KH004 pfcrt allele is required for reduced PPQ sensitivity, whereas copy number variation in pm2/3 further decreases susceptibility but does not confer resistance in the absence of additional mutations in pfcrt . A deep investigation of genotype-phenotype relationships demonstrates that progeny clones from experimental crosses can be used to understand the relative contributions of pfcrt , pm2/3 , and parasite genetic background to a range of PPQ-related traits. Additionally, we find that the resistance phenotype associated with parasites inheriting the G367C substitution in pfcrt is consistent with previously validated PPQ resistance mutations in this transporter.IMPORTANCEResistance to piperaquine, used in combination with dihydroartemisinin, has emerged in Cambodia and threatens to spread to other malaria-endemic regions. Understanding the causal mutations of drug resistance and their impact on parasite fitness is critical for surveillance and intervention and can also reveal new avenues to limiting the evolution and spread of drug resistance. An experimental genetic cross is a powerful tool for pinpointing the genetic determinants of key drug resistance and fitness phenotypes and has the distinct advantage of quantifying the effects of naturally evolved genetic variation. Our study was strengthened since the full range of copies of KH004 pm2/3 was inherited among the progeny clones, allowing us to directly test the role of the pm2/3 copy number on resistance-related phenotypes in the context of a unique pfcrt allele. Our multigene model suggests an important role for both loci in the evolution of this multidrug-resistant parasite lineage., Competing Interests: The authors declare no conflict of interest.

Pathogen genomics is a powerful tool for tracking infectious disease transmission. In malaria, identity-by-descent (IBD) is used to assess the genetic relatedness between parasites and has been used to study transmission and importation. In theory, IBD can be used to distinguish genealogical relationships to reconstruct transmission history or identify parasites for genotype-to-phenotype quantitative-trait-locus experiments. MalKinID (Malaria Kinship Identifier) is a new likelihood-based classification model designed to identify genealogical relationships among malaria parasites based on genome-wide IBD proportions and IBD segment distributions. MalKinID was calibrated to the genomic data from three 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 accurately 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 of 0.76 (0.56, 0.92) and 0.23 (0.0, 0.4) for PC and FS and <0.05 for second-degree and third-degree relatives. Genealogical inference is most powered 1) when outcrossing is the norm or 2) when multi-sample comparisons based on a predefined pedigree are used. MalKinID lays the foundations for using IBD to track parasite transmission history and for separating progeny for quantitative-trait-locus experiments., Competing Interests: Conflicts of Interests The authors declare no relevant conflicts of interest.

Background: Artemisinin resistance in Plasmodium falciparum threatens global malaria elimination efforts. To contain and then eliminate artemisinin resistance in Eastern Myanmar a network of community-based malaria posts was instituted and targeted mass drug administration (MDA) with dihydroartemisinin-piperaquine (three rounds at monthly intervals) was conducted. The prevalence of artemisinin resistance during the elimination campaign (2013-2019) was characterized., Methods: Throughout the six-year campaign Plasmodium falciparum positive blood samples from symptomatic patients and from cross-sectional surveys were genotyped for mutations in kelch-13-a molecular marker of artemisinin resistance., Result: The program resulted in near elimination of falciparum malaria. Of 5162 P. falciparum positive blood samples genotyped, 3281 (63.6%) had K13 mutations. The prevalence of K13 mutations was 73.9% in 2013 and 64.4% in 2019. Overall, there was a small but significant decline in the proportion of K13 mutants (p < 0.001). In the MDA villages there was no significant change in the K13 proportions before and after MDA. The distribution of different K13 mutations changed substantially; F446I and P441L mutations increased in both MDA and non-MDA villages, while most other K13 mutations decreased. The proportion of C580Y mutations fell from 9.2% (43/467) before MDA to 2.3% (19/813) after MDA (p < 0.001). Similar changes occurred in the 487 villages where MDA was not conducted., Conclusion: The malaria elimination program in Kayin state, eastern Myanmar, led to a substantial reduction in falciparum malaria. Despite the intense use of artemisinin-based combination therapies, both in treatment and MDA, this did not select for artemisinin resistance., (© 2024. The Author(s).)

Background: The role of pathogen genotype in determining disease severity and immunopathology has been studied intensively in microbial pathogens including bacteria, fungi, protozoa and viruses but is poorly understood in parasitic helminths. The medically important blood fluke Schistosoma mansoni is an excellent model system to study the impact of helminth genetic variation on immunopathology. Our laboratory has demonstrated that laboratory schistosome populations differ in sporocyst growth and cercarial production in the intermediate snail host and worm establishment and fecundity in the vertebrate host. Here, we (i) investigate the hypothesis that schistosome genotype plays a significant role in immunopathology and related parasite life history traits in the vertebrate mouse host and (ii) quantify the relative impact of parasite and host genetics on infection outcomes., Methods: We infected BALB/c and C57BL/6 mice with four different laboratory schistosome populations from Africa and the Americas. We quantified disease progression in the vertebrate host by measuring body weight and complete blood count (CBC) with differential over a 12-week infection period. On sacrifice, we assessed parasitological (egg and worm counts, fecundity), immunopathological (organ measurements and histopathology) and immunological (CBC with differential and cytokine profiles) characteristics to determine the impact of parasite and host genetics., Results: We found significant variation between parasite populations in worm numbers, fecundity, liver and intestine egg counts, liver and spleen weight, and fibrotic area but not in granuloma size. Variation in organ weight was explained by egg burden and intrinsic parasite factors independent of egg burden. We found significant variation between infected mouse lines in cytokine levels (IFN-γ, TNF-α), eosinophils, lymphocytes and monocyte counts., Conclusions: This study showed that both parasite and host genotype impact the outcome of infection. While host genotype explains most of the variation in immunological traits, parasite genotype explains most of the variation in parasitological traits, and both host and parasite genotypes impact immunopathology outcomes., (© 2024. The Author(s).)

More than 60 zoonoses are linked to small mammals, including some of the most devastating pathogens in human history. Millions of museum-archived tissues are available to understand natural history of those pathogens. Our goal was to maximize the value of museum collections for pathogen-based research by using targeted sequence capture. We generated a probe panel that includes 39,916 80-bp RNA probes targeting 32 pathogen groups, including bacteria, helminths, fungi, and protozoans. Laboratory-generated, mock-control samples showed that we are capable of enriching targeted loci from pathogen DNA 2,882‒6,746-fold. We identified bacterial species in museum-archived samples, including Bartonella, a known human zoonosis. These results showed that probe-based enrichment of pathogens is a highly customizable and efficient method for identifying pathogens from museum-archived tissues.

Malaria parasites break down host haemoglobin into peptides and amino acids in the digestive vacuole for export to the parasite cytoplasm for growth: interrupting this process is central to the mode of action of several antimalarial drugs. Mutations in the chloroquine (CQ) resistance transporter, pfcrt, located in the digestive vacuole membrane, confer CQ resistance in Plasmodium falciparum, and typically also affect parasite fitness. However, the role of other parasite loci in the evolution of CQ resistance is unclear. Here we use a combination of population genomics, genetic crosses and gene editing to demonstrate that a second vacuolar transporter plays a key role in both resistance and compensatory evolution. Longitudinal genomic analyses of the Gambian parasites revealed temporal signatures of selection on a putative amino acid transporter (pfaat1) variant S258L, which increased from 0% to 97% in frequency between 1984 and 2014 in parallel with the pfcrt1 K76T variant. Parasite genetic crosses then identified a chromosome 6 quantitative trait locus containing pfaat1 that is selected by CQ treatment. Gene editing demonstrated that pfaat1 S258L potentiates CQ resistance but at a cost of reduced fitness, while pfaat1 F313S, a common southeast Asian polymorphism, reduces CQ resistance while restoring fitness. Our analyses reveal hidden complexity in CQ resistance evolution, suggesting that pfaat1 may underlie regional differences in the dynamics of resistance evolution, and modulate parasite resistance or fitness by manipulating the balance between both amino acid and drug transport., (© 2023. The Author(s).)

Background: The trematode parasite Schistosoma mansoni uses an aquatic snail intermediate and a vertebrate definitive host to complete its life cycle. We previously showed that a key transmission trait-the number of cercariae larvae shed from infected Biomphalaria spp. snails-varies significantly within and between different parasite populations and is genetically controlled by five loci. We investigated the hypothesis that the success of parasite genotypes showing high propagative fitness in the intermediate snail host may be offset by lower reproductive fitness in the definitive vertebrate host., Methods: We investigated this trade-off hypothesis by selecting parasite progeny producing high or low number of larvae in the snail and then comparing fitness parameters and virulence in the rodent host. We infected inbred BALB/c mice using two Schistosoma mansoni parasite lines [high shedder (HS) and low shedder (LS) lines] isolated from F2 progeny generated by genetic crosses between SmLE (HS parent) and SmBRE (LS parent) parasites. We used the F3 progeny to infect two populations of inbred Biomphalaria glabrata snails. We then compared life history traits and virulence of these two selected parasite lines in the rodent host to understand pleiotropic effects of genes determining cercarial shedding in parasites infecting the definitive host., Results: HS parasites shed high numbers of cercariae, which had a detrimental impact on snail physiology (measured by laccase-like activity and hemoglobin rate), regardless of the snail genetic background. In contrast, selected LS parasites shed fewer cercariae and had a lower impact on snail physiology. Similarly, HS worms have a higher reproductive fitness and produced more viable F3 miracidia larvae than LS parasites. This increase in transmission is correlated with an increase in virulence toward the rodent host, characterized by stronger hepato-splenomegaly and hepatic fibrosis., Conclusions: These experiments revealed that schistosome parasite propagative and reproductive fitness was positively correlated in intermediate and definitive host (positive pleiotropy). Therefore, we rejected our trade-off hypothesis. We also showed that our selected schistosome lines exhibited low and high shedding phenotype regardless of the intermediate snail host genetic background. ​., (© 2023. The Author(s).)

The human malaria parasite Plasmodium falciparum is globally widespread, but its prevalence varies significantly between and even within countries. Most population genetic studies in P. falciparum focus on regions of high transmission where parasite populations are large and genetically diverse, such as sub-Saharan Africa. Understanding population dynamics in low transmission settings, however, is of particular importance as these are often where drug resistance first evolves. Here, we use the Pacific Coast of Colombia and Ecuador as a model for understanding the population structure and evolution of Plasmodium parasites in small populations harboring less genetic diversity. The combination of low transmission and a high proportion of monoclonal infections means there are few outcrossing events and clonal lineages persist for long periods of time. Yet despite this, the population is evolutionarily labile and has successfully adapted to changes in drug regime. Using newly sequenced whole genomes, we measure relatedness between 166 parasites, calculated as identity by descent (IBD), and find 17 distinct but highly related clonal lineages, six of which have persisted in the region for at least a decade. This inbred population structure is captured in more detail with IBD than with other common population structure analyses like PCA, ADMIXTURE, and distance-based trees. We additionally use patterns of intra-chromosomal IBD and an analysis of haplotypic variation to explore past selection events in the region. Two genes associated with chloroquine resistance, crt and aat1, show evidence of hard selective sweeps, while selection appears soft and/or incomplete at three other key resistance loci (dhps, mdr1, and dhfr). Overall, this work highlights the strength of IBD analyses for studying parasite population structure and resistance evolution in regions of low transmission, and emphasizes that drug resistance can evolve and spread in small populations, as will occur in any region nearing malaria elimination., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2022 Carrasquilla et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

In 2005 Richard Carter's group surprised the malaria genetics community with an elegant approach to rapidly mapping the genetic basis of phenotypic traits in rodent malaria parasites. This approach, which he termed "linkage group selection", utilized bulk pools of progeny, rather than individual clones, and exploited simple selection schemes to identify genome regions underlying resistance to drug treatment (or other phenotypes). This work was the first application of "bulk segregant" methodologies for genetic mapping in microbes: this approach is now widely used in yeast, and across multiple recombining pathogens ranging from Aspergillus fungi to Schistosome parasites. Genetic crosses of human malaria parasites (for which Richard Carter was also a pioneer) can now be conducted in humanized mice, providing new opportunities for exploiting bulk segregant approaches for a wide variety of malaria parasite traits. We review the application of bulk segregant approaches to mapping malaria parasite traits and suggest additional developments that may further expand the utility of this powerful approach., Competing Interests: Declaration of Competing Interest None., (Copyright © 2022 The Authors. Published by Elsevier B.V. All rights reserved.)

Oxamniquine (OXA) is a prodrug activated by a sulfotransferase ( SULT ) that was only active against Schistosoma mansoni . We have reengineered OXA to be effective against S. haematobium and S. japonicum . Three derivatives stand out, CIDD-0066790, CIDD-0072229, and CIDD-0149830 as they kill all three major human schistosome species. However, questions remain. Is the OXA mode of action conserved in derivatives? RNA-interference experiments demonstrate that knockdown of the SmSULT , ShSULT , and SjSULT results in resistance to CIDD-0066790. Confirming that the OXA-derivative mode of action is conserved. Next is the level of expression of the schistosome SULTs in each species, as well as changes in SULT expression throughout development in S. mansoni . Using multiple tools, our data show that SmSULT has higher expression compared to ShSULT and SjSULT . Third, is the localization of SULT in the adult, multicellular eucaryotic schistosome species. We utilized fluorescence in situ hybridization and uptake of radiolabeled OXA to determine that multiple cell types throughout the adult schistosome worm express SULT . Thus, we hypothesize the ability of many cells to express the sulfotransferase accounts for the ability of the OXA derivatives to kill adult worms. Our studies demonstrate that the OXA derivatives are able to kill all three human schistosome species and thus will be a useful complement to PZQ.

What genes determine in vitro growth and nutrient utilization in asexual blood-stage malaria parasites? Competition experiments between NF54, clone 3D7, a lab-adapted African parasite, and a recently isolated Asian parasite (NHP4026) reveal contrasting outcomes in different media: 3D7 outcompetes NHP4026 in media containing human serum, while NHP4026 outcompetes 3D7 in media containing AlbuMAX, a commercial lipid-rich bovine serum formulation. To determine the basis for this polymorphism, we conducted parasite genetic crosses using humanized mice and compared genome-wide allele frequency changes in three independent progeny populations cultured in media containing human serum or AlbuMAX. This bulk segregant analysis detected three quantitative trait loci (QTL) regions [on chromosome (chr) 2 containing aspartate transaminase AST ; chr 13 containing EBA-140; and chr 14 containing cysteine protease ATG4 ] linked with differential growth in serum or AlbuMAX in each of the three independent progeny pools. Selection driving differential growth was strong ( s = 0.10 - 0.23 per 48-hour lifecycle). We conducted validation experiments for the strongest QTL on chr 13: competition experiments between ΔEBA-140 and 3D7 wildtype parasites showed fitness reversals in the two medium types as seen in the parental parasites, validating this locus as the causative gene. These results (i) demonstrate the effectiveness of bulk segregant analysis for dissecting fitness traits in P. falciparum genetic crosses, and (ii) reveal intimate links between red blood cell invasion and nutrient composition of growth media. Use of parasite crosses combined with bulk segregant analysis will allow systematic dissection of key nutrient acquisition/metabolism and red blood cell invasion pathways in P. falciparum., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Kumar, Li, McDew-White, Reyes, Delgado, Sayeed, Haile, Abatiyow, Kennedy, Camargo, Checkley, Brenneman, Button-Simons, Duraisingh, Cheeseman, Kappe, Nosten, Ferdig, Vaughan and Anderson.)

Aquatic snails, the intermediate hosts of schistosomes, harbor a diverse unexplored microbiome. We speculate that this may play a critical role in host-parasite interactions. We summarize our current knowledge of snail microbiomes and highlight future research priorities., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 Elsevier Ltd. All rights reserved.)

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., (© 2022 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Classical malaria parasite genetic crosses involve isolation, genotyping, and phenotyping of progeny parasites, which is time consuming and laborious. We tested a rapid alternative approach-bulk segregant analysis (BSA)-that utilizes sequencing of bulk progeny populations with and without drug selection for rapid identification of drug resistance loci. We used dihydroartemisinin (DHA) selection in two genetic crosses and investigated how synchronization, cryopreservation, and the drug selection regimen impacted BSA success. We detected a robust quantitative trait locus (QTL) at kelch13 in both crosses but did not detect QTLs at four other candidate loci. QTLs were detected using synchronized, but not unsynchronized progeny pools, consistent with the stage-specific action of DHA. We also successfully applied BSA to cryopreserved progeny pools, expanding the utility of this approach. We conclude that BSA provides a powerful approach for investigating the genetic architecture of drug resistance in Plasmodium falciparum ., Competing Interests: The authors declare no competing interests., (© 2022 The Authors.)

Mass drug administration with praziquantel (PZQ) monotherapy is considered the mainstay for control and elimination of the parasites causing schistosomiasis in humans. This drug shows imperfect cure rates in the field, and parasites showing reduced PZQ response can be selected in the laboratory, but the extent of resistance in Schistosoma mansoni populations is unknown. We examined the genetic basis of the variation in response in a PZQ-selected S. mansoni population (SmLE-PZQ-R) in which 35% of the parasitic worms survive high-dose PZQ (73 micrograms per milliliter) treatment. We used genome-wide association to map loci underlying PZQ response and identified a transient receptor potential ( Sm. TRPM PZQ ) channel (Smp&#95;246790) within the major chromosome 3 peak that is activated by nanomolar concentrations of PZQ. The PZQ response showed recessive inheritance and marker-assisted selection of parasites at a single Sm. TRPM PZQ SNP that produced populations of PZQ-enriched resistant (PZQ-ER) and PZQ-enriched sensitive (PZQ-ES) parasites, exhibiting >377-fold difference in PZQ response. The PZQ-ER parasites survived treatment in rodents at higher frequencies compared with PZQ-ES, and resistant parasites exhibited 2.25-fold lower expression of Sm. TRPM PZQ relative to sensitive parasites. Specific chemical blockers of Sm. TRPM PZQ enhanced PZQ resistance, whereas Sm. TRPM PZQ activators increased sensitivity. We surveyed Sm. TRPM PZQ sequence variations in 259 parasites from different global sites and identified one nonsense mutation that resulted in a truncated protein with no PZQ binding site. Our results demonstrate that Sm. TRPM PZQ underlies variation in PZQ responses in S. mansoni and provides an approach for monitoring emerging PZQ-resistant alleles in schistosome elimination programs.

Both theory and experimental data from pathogens suggest that the production of transmission stages should be strongly associated with virulence, but the genetic bases of parasite transmission/virulence traits are poorly understood. The blood fluke Schistosoma mansoni shows extensive variation in numbers of cercariae larvae shed and in their virulence to infected snail hosts, consistent with expected trade-offs between parasite transmission and virulence. We crossed schistosomes from two populations that differ 8-fold in cercarial shedding and in their virulence to Biomphalaria glabrata snail hosts, and determined four-week cercarial shedding profiles in F0 parents, F1 parents and 376 F2 progeny from two independent crosses in inbred snails. Sequencing and linkage analysis revealed that cercarial production is polygenic and controlled by five QTLs (i.e. Quantitative Trait Loci). These QTLs act additively, explaining 28.56% of the phenotypic variation. These results demonstrate that the genetic architecture of key traits relevant to schistosome ecology can be dissected using classical linkage mapping approaches.

Population genomics of bulk malaria infections is unable to examine intrahost evolution; therefore, most work has focused on the role of recombination in generating genetic variation. We used single-cell sequencing protocol for low-parasitaemia infections to generate 406 near-complete single Plasmodium vivax genomes from 11 patients sampled during sequential febrile episodes. Parasite genomes contain hundreds of de novo mutations, showing strong signatures of selection, which are enriched in the ApiAP2 family of transcription factors, known targets of adaptation. Comparing 315 P. falciparum single-cell genomes from 15 patients with our P. vivax data, we find broad complementary patterns of de novo mutation at the gene and pathway level, revealing the importance of within-host evolution during malaria infections., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 Elsevier Inc. All rights reserved.)

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed.  Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination., Competing Interests: No competing interests were disclosed., (Copyright: © 2021 MalariaGEN et al.)

Schistosomiasis control efforts in Nigeria received a boost in 2016 when Merck Group made the largest single donation of praziquantel to an African country. We examined urine samples from 2,023 school age children from 15 locations in 10 states and an Internally Displaced Person's (IDP) camp in Nigeria. We recorded an overall Schistosoma haematobium prevalence of 10.4% in the 10 states that ranged between 6 - 37%, while prevalence in the IDP camp was 2.9%. The highest infection prevalence (37%) recorded was from the population in Wasai Dam area in Minjibir (Kano State), while five locations had no positive urine samples. We observed heavy intensity of infection (≥ 50 eggs/10 ml urine) in 87.9% of infected samples and co-occurrence of the eggs of S. haematobium and S. mansoni in urine for two participants. The overall prevalence we recorded is slightly above the national average (9.5%) reported in 2015. Our findings indicate that despite the ongoing administration of praziquantel in Nigeria, urogenital schistosomiasis is still prevalent with heavy intensity of infection. Large-scale epidemiological monitoring is required to monitor the efficacy of schistosomiasis control in Nigeria., (Copyright © 2021 Elsevier B.V. All rights reserved.)

Genetic crosses are most powerful for linkage analysis when progeny numbers are high, parental alleles segregate evenly and numbers of inbred progeny are minimized. We previously developed a novel genetic crossing platform for the human malaria parasite Plasmodium falciparum, an obligately sexual, hermaphroditic protozoan, using mice carrying human hepatocytes (the human liver-chimeric FRG NOD huHep mouse) as the vertebrate host. We report on two genetic crosses-(1) an allopatric cross between a laboratory-adapted parasite (NF54) of African origin and a recently patient-derived Asian parasite, and (2) a sympatric cross between two recently patient-derived Asian parasites. We generated 144 unique recombinant clones from the two crosses, doubling the number of unique recombinant progeny generated in the previous 30 years. The allopatric African/Asian cross has minimal levels of inbreeding and extreme segregation distortion, while in the sympatric Asian cross, inbred progeny predominate and parental alleles segregate evenly. Using simulations, we demonstrate that these progeny provide the power to map small-effect mutations and epistatic interactions. The segregation distortion in the allopatric cross slightly erodes power to detect linkage in several genome regions. We greatly increase the power and the precision to map biomedically important traits with these new large progeny panels.

Non-human primate models will expedite therapeutics and vaccines for coronavirus disease 2019 (COVID-19) to clinical trials. Here, we compare acute severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in young and old rhesus macaques, baboons and old marmosets. Macaques had clinical signs of viral infection, mild to moderate pneumonitis and extra-pulmonary pathologies, and both age groups recovered in two weeks. Baboons had prolonged viral RNA shedding and substantially more lung inflammation compared with macaques. Inflammation in bronchoalveolar lavage was increased in old versus young baboons. Using techniques including computed tomography imaging, immunophenotyping, and alveolar/peripheral cytokine response and immunohistochemical analyses, we delineated cellular immune responses to SARS-CoV-2 infection in macaque and baboon lungs, including innate and adaptive immune cells and a prominent type-I interferon response. Macaques developed T-cell memory phenotypes/responses and bystander cytokine production. Old macaques had lower titres of SARS-CoV-2-specific IgG antibody levels compared with young macaques. Acute respiratory distress in macaques and baboons recapitulates the progression of COVID-19 in humans, making them suitable as models to test vaccines and therapies.

The microbiome - the microorganism community that is found on or within an organism's body - is increasingly recognized to shape many aspects of its host biology and is a key determinant of health and disease. Microbiomes modulate the capacity of insect disease vectors (mosquitoes, tsetse flies, sandflies) to transmit parasites and disease. We investigate the diversity and abundance of microorganisms within the hemolymph (i.e. blood) of Biomphalaria snails, the intermediate host for Schistosoma mansoni, using Illumina MiSeq sequencing of the bacterial 16S V4 rDNA. We sampled hemolymph from five snails from six different laboratory populations of B. glabrata and one population of B. alexandrina. We observed 279.84 ± 0.79 amplicon sequence variants per snail. There were significant differences in microbiome composition at the level of individual snails, snail populations and species. Snail microbiomes were dominated by Proteobacteria and Bacteroidetes while water microbiomes from snail tank were dominated by Actinobacteria. We investigated the absolute bacterial load using qPCR: hemolymph samples contained 2784 ± 339 bacteria/μl. We speculate that the microbiome may represent a critical, but unexplored intermediary in the snail-schistosome interaction as hemolymph is in very close contact with the parasite at each step of its development., (© 2020 Society for Applied Microbiology and John Wiley & Sons Ltd.)

Vaccine and antiviral development against SARS-CoV-2 infection or COVID-19 disease would benefit from validated small animal models. Here, we show that transgenic mice expressing human angiotensin-converting enzyme 2 (hACE2) by the human cytokeratin 18 promoter (K18 hACE2) represent a susceptible rodent model. K18 hACE2 transgenic mice succumbed to SARS-CoV-2 infection by day 6, with virus detected in lung airway epithelium and brain. K18 ACE2 transgenic mice produced a modest TH1/2/17 cytokine storm in the lung and spleen that peaked by day 2, and an extended chemokine storm that was detected in both lungs and brain. This chemokine storm was also detected in the brain at day 6. K18 hACE2 transgenic mice are, therefore, highly susceptible to SARS-CoV-2 infection and represent a suitable animal model for the study of viral pathogenesis, and for identification and characterization of vaccines (prophylactic) and antivirals (therapeutics) for SARS-CoV-2 infection and associated severe COVID-19 disease.

Characterising connectivity between geographically separated biological populations is a common goal in many fields. Recent approaches to understanding connectivity between malaria parasite populations, with implications for disease control efforts, have used estimates of relatedness based on identity-by-descent (IBD). However, uncertainty around estimated relatedness has not been accounted for. IBD-based relatedness estimates with uncertainty were computed for pairs of monoclonal Plasmodium falciparum samples collected from five cities on the Colombian-Pacific coast where long-term clonal propagation of P. falciparum is frequent. The cities include two official ports, Buenaventura and Tumaco, that are separated geographically but connected by frequent marine traffic. Fractions of highly-related sample pairs (whose classification using a threshold accounts for uncertainty) were greater within cities versus between. However, based on both highly-related fractions and on a threshold-free approach (Wasserstein distances between parasite populations) connectivity between Buenaventura and Tumaco was disproportionally high. Buenaventura-Tumaco connectivity was consistent with transmission events involving parasites from five clonal components (groups of statistically indistinguishable parasites identified under a graph theoretic framework). To conclude, P. falciparum population connectivity on the Colombian-Pacific coast abides by accessibility not isolation-by-distance, potentially implicating marine traffic in malaria transmission with opportunities for targeted intervention. Further investigations are required to test this hypothesis. For the first time in malaria epidemiology (and to our knowledge in ecological and epidemiological studies more generally), we account for uncertainty around estimated relatedness (an important consideration for studies that plan to use genotype versus whole genome sequence data to estimate IBD-based relatedness); we also use threshold-free methods to compare parasite populations and identify clonal components. Threshold-free methods are especially important in analyses of malaria parasites and other recombining organisms with mixed mating systems where thresholds do not have clear interpretation (e.g. due to clonal propagation) and thus undermine the cross-comparison of studies., Competing Interests: The authors have declared that no competing interests exist.

Currently there is only one method of treatment for human schistosomiasis, the drug praziquantel. Strong selective pressure has caused a serious concern for a rise in resistance to praziquantel leading to the necessity for additional pharmaceuticals, with a distinctly different mechanism of action, to be used in combination therapy with praziquantel. Previous treatment of Schistosoma mansoni included the use of oxamniquine (OXA), a prodrug that is enzymatically activated in S. mansoni but is ineffective against S. haematobium and S. japonicum. The oxamniquine activating enzyme was identified as a S. mansoni sulfotransferase (SmSULT-OR). Structural data have allowed for directed drug development in reengineering oxamniquine to be effective against S. haematobium and S. japonicum. Guided by data from X-ray crystallographic studies and Schistosoma worm killing assays on oxamniquine, our structure-based drug design approach produced a robust SAR program that tested over 300 derivatives and identified several new lead compounds with effective worm killing in vitro. Previous studies resulted in the discovery of compound CIDD-0066790, which demonstrated broad-species activity in killing of schistosome species. As these compounds are racemic mixtures, we tested and demonstrate that the R enantiomer CIDD-007229 kills S. mansoni, S. haematobium and S. japonicum better than the parent drug (CIDD-0066790). The search for derivatives that kill better than CIDD-0066790 has resulted in a derivative (CIDD- 149830) that kills 100% of S. mansoni, S. haematobium and S. japonicum adult worms within 7 days. We hypothesize that the difference in activation and thus killing by the derivatives is due to the ability of the derivative to fit in the binding pocket of each sulfotransferase (SmSULT-OR, ShSULT-OR, SjSULT-OR) and to be efficiently sulfated. The purpose of this research is to develop a second drug to be used in conjunction with praziquantel to treat the major human species of Schistosoma. Collectively, our findings show that CIDD-00149830 and CIDD-0072229 are promising novel drugs for the treatment of human schistosomiasis and strongly support further development and in vivo testing., Competing Interests: The authors have declared that no competing interests exist.

Human schistosomiasis is a disease which globally affects over 229 million people. Three major species affecting humans are Schistosoma mansoni, S. haematobium and S. japonicum. Previous treatment of S. mansoni includes the use of oxamniquine (OXA), a prodrug that is enzymatically activated in S. mansoni but is ineffective against S. haematobium and S. japonicum. The OXA activating enzyme was identified and crystallized, as being a S. mansoni sulfotransferase (SmSULT). S. haematobium and S. japonicum possess homologs of SmSULT (ShSULT and SjSULT) begging the question; why does oxamniquine fail to kill S. haematobium and S. japonicum adult worms? Investigation of the molecular structures of the sulfotransferases indicates that structural differences, specifically in OXA contact residues, do not abrogate OXA binding in the active sites as previously hypothesized. Data presented argue that the ability of SULTs to sulfate and thus activate OXA and its derivatives is linked to the ability of OXA to fit in the binding pocket to allow the transfer of a sulfur group., Competing Interests: Declaration of competing interest The authors declare that they have no conflicts of interest., (Copyright © 2020 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Background: Tracking and understanding artemisinin resistance is key for preventing global setbacks in malaria eradication efforts. The ring-stage survival assay (RSA) is the current gold standard for in vitro artemisinin resistance phenotyping. However, the RSA has several drawbacks: it is relatively low throughput, has high variance due to microscopy readout, and correlates poorly with the current benchmark for in vivo resistance, patient clearance half-life post-artemisinin treatment. Here a modified RSA is presented, the extended Recovery Ring-stage Survival Assay (eRRSA), using 15 cloned patient isolates from Southeast Asia with a range of patient clearance half-lives, including parasite isolates with and without kelch13 mutations., Methods: Plasmodium falciparum cultures were synchronized with single layer Percoll during the schizont stage of the intraerythrocytic development cycle. Cultures were left to reinvade to early ring-stage and parasitaemia was quantified using flow cytometry. Cultures were diluted to 2% haematocrit and 0.5% parasitaemia in a 96-well plate to start the assay, allowing for increased throughput and decreased variability between biological replicates. Parasites were treated with 700 nM of dihydroartemisinin or 0.02% dimethyl sulfoxide (DMSO) for 6 h, washed three times in drug-free media, and incubated for 66 or 114 h, when samples were collected and frozen for PCR amplification. A SYBR Green-based quantitative PCR method was used to quantify the fold-change between treated and untreated samples., Results: 15 cloned patient isolates from Southeast Asia with a range of patient clearance half-lives were assayed using the eRRSA. Due to the large number of pyknotic and dying parasites at 66 h post-exposure (72 h sample), parasites were grown for an additional cell cycle (114 h post-exposure, 120 h sample), which drastically improved correlation with patient clearance half-life compared to the 66 h post-exposure sample. A Spearman correlation of - 0.8393 between fold change and patient clearance half-life was identified in these 15 isolates from Southeast Asia, which is the strongest correlation reported to date., Conclusions: eRRSA drastically increases the efficiency and accuracy of in vitro artemisinin resistance phenotyping compared to the traditional RSA, which paves the way for extensive in vitro phenotyping of hundreds of artemisinin resistant parasites.

In high-transmission regions, we expect parasite lineages within complex malaria infections to be unrelated due to parasite inoculations from different mosquitoes. This project was designed to test this prediction. We generated 485 single-cell genome sequences from fifteen P. falciparum malaria patients from Chikhwawa, Malawi-an area of intense transmission. Patients harbored up to seventeen unique parasite lineages. Surprisingly, parasite lineages within infections tend to be closely related, suggesting that superinfection by repeated mosquito bites is rarer than co-transmission of parasites from a single mosquito. Both closely and distantly related parasites comprise an infection, suggesting sequential transmission of complex infections between multiple hosts. We identified tetrads and reconstructed parental haplotypes, which revealed the inbred ancestry of infections and non-Mendelian inheritance. Our analysis suggests strong barriers to secondary infection and outbreeding amongst malaria parasites from a high transmission setting, providing unexpected insights into the biology and transmission of malaria., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Do mutations required for adaptation occur de novo, or are they segregating within populations as standing genetic variation? This question is key to understanding adaptive change in nature, and has important practical consequences for the evolution of drug resistance. We provide evidence that alleles conferring resistance to oxamniquine (OXA), an antischistosomal drug, are widespread in natural parasite populations under minimal drug pressure and predate OXA deployment. OXA has been used since the 1970s to treat Schistosoma mansoni infections in the New World where S. mansoni established during the slave trade. Recessive loss-of-function mutations within a parasite sulfotransferase (SmSULT-OR) underlie resistance, and several verified resistance mutations, including a deletion (p.E142del), have been identified in the New World. Here we investigate sequence variation in SmSULT-OR in S. mansoni from the Old World, where OXA has seen minimal usage. We sequenced exomes of 204 S. mansoni parasites from West Africa, East Africa and the Middle East, and scored variants in SmSULT-OR and flanking regions. We identified 39 non-synonymous SNPs, 4 deletions, 1 duplication and 1 premature stop codon in the SmSULT-OR coding sequence, including one confirmed resistance deletion (p.E142del). We expressed recombinant proteins and used an in vitro OXA activation assay to functionally validate the OXA-resistance phenotype for four predicted OXA-resistance mutations. Three aspects of the data are of particular interest: (i) segregating OXA-resistance alleles are widespread in Old World populations (4.29-14.91% frequency), despite minimal OXA usage, (ii) two OXA-resistance mutations (p.W120R, p.N171IfsX28) are particularly common (>5%) in East African and Middle-Eastern populations, (iii) the p.E142del allele has identical flanking SNPs in both West Africa and Puerto Rico, suggesting that parasites bearing this allele colonized the New World during the slave trade and therefore predate OXA deployment. We conclude that standing variation for OXA resistance is widespread in S. mansoni., Competing Interests: The authors have declared that no competing interests exist.

Background: Parasite traits associated with transmission success, such as the number of infective stages released from the host, are expected to be optimized by natural selection. However, in the trematode parasite Schistosoma mansoni, a key transmission trait, i.e. the number of cercariae larvae shed from infected Biomphalaria spp. snails, varies significantly within and between different parasite populations and selection experiments demonstrate that this variation has a strong genetic basis. In this study, we compared the transmission strategies of two laboratory schistosome population and their consequences for their snail host., Methods: We infected inbred Biomphalaria glabrata snails using two S. mansoni parasite populations (SmBRE and SmLE), both isolated from Brazil and maintained in the laboratory for decades. We compared life history traits of these two parasite populations by quantifying sporocyst growth within infected snails (assayed using qPCR), output of cercaria larvae and impact on snail host physiological response (i.e. hemoglobin rate, laccase-like activity) and survival., Results: We identified striking differences in virulence and transmission between the two studied parasite populations. SmBRE (low shedder (LS) parasite population) sheds very low numbers of cercariae and causes minimal impact on the snail physiological response (i.e. laccase-like activity, hemoglobin rate and snail survival). In contrast, SmLE (high shedder (HS) parasite population) sheds 8-fold more cercariae (mean ± SE cercariae per shedding: 284 ± 19 vs 2352 ± 113), causes high snail mortality and has strong impact on snail physiology. We found that HS sporocysts grow more rapidly inside the snail host, comprising up to 60% of cells within infected snails, compared to LS sporocysts, which comprised up to 31%. Cercarial production is strongly correlated to the number of S. mansoni sporocyst cells present within the snail host tissue, although the proportion of sporocyst cells alone does not explain the low cercarial shedding of SmBRE., Conclusions: We demonstrated the existence of alternative transmission strategies in the S. mansoni parasite consistent with trade-offs between parasite transmission and host survival: a "boom-bust" strategy characterized by high virulence, high transmission and short duration infections and a "slow and steady" strategy with low virulence, low transmission but long duration of snail host infections.

Determining the genetic basis of fitness is central to understanding evolution and transmission of microbial pathogens. In human malaria parasites (Plasmodium falciparum), most experimental work on fitness has focused on asexual blood stage parasites, because this stage can be easily cultured, although the transmission of malaria requires both female Anopheles mosquitoes and vertebrate hosts. We explore a powerful approach to identify the genetic determinants of parasite fitness across both invertebrate and vertebrate life-cycle stages of P. falciparum. This combines experimental genetic crosses using humanized mice, with selective whole genome amplification and pooled sequencing to determine genome-wide allele frequencies and identify genomic regions under selection across multiple lifecycle stages. We applied this approach to genetic crosses between artemisinin resistant (ART-R, kelch13-C580Y) and ART-sensitive (ART-S, kelch13-WT) parasites, recently isolated from Southeast Asian patients. Two striking results emerge: we observed (i) a strong genome-wide skew (>80%) towards alleles from the ART-R parent in the mosquito stage, that dropped to ~50% in the blood stage as selfed ART-R parasites were selected against; and (ii) repeatable allele specific skews in blood stage parasites with particularly strong selection (selection coefficient (s) ≤ 0.18/asexual cycle) against alleles from the ART-R parent at loci on chromosome 12 containing MRP2 and chromosome 14 containing ARPS10. This approach robustly identifies selected loci and has strong potential for identifying parasite genes that interact with the mosquito vector or compensatory loci involved in drug resistance., Competing Interests: The authors have declared that no competing interests exist.