Genetic diversity, multiplicity of infection and population structure of Schistosoma mansoni isolates from human hosts in Ethiopia

MOIRE: A software package for the estimation of allele frequencies and effective multiplicity of infection from polyallelic data

Malaria parasite genetic data can provide insight into parasite phenotypes, evolution, and transmission. However, estimating key parameters such as allele frequencies, multiplicity of infection (MOI), and within-host relatedness from genetic data has been challenging, particularly in the presence of multiple related coinfecting strains. Existing methods often rely on single nucleotide polymorphism (SNP) data and do not account for within-host relatedness. In this study, we introduce a Bayesian approach called MOIRE (Multiplicity Of Infection and allele frequency REcovery), designed to estimate allele frequencies, MOI, and within-host relatedness from genetic data subject to experimental error. Importantly, MOIRE is flexible in accommodating both polyallelic and SNP data, making it adaptable to diverse genotyping panels. We also introduce a novel metric, the effective MOI ( e M O I ) , which integrates MOI and within-host relatedness, providing a robust and interpretable measure of genetic diversity. Using extensive simulations and real-world data from a malaria study in Namibia, we demonstrate the superior performance of MOIRE over naive estimation methods, accurately estimating MOI up to 7 with moderate sized panels of diverse loci (e.g. microhaplotypes). MOIRE also revealed substantial heterogeneity in population mean MOI and mean relatedness across health districts in Namibia, suggesting detectable differences in transmission dynamics. Notably, e M O I emerges as a portable metric of within-host diversity, facilitating meaningful comparisons across settings, even when allele frequencies or genotyping panels are different. MOIRE represents an important addition to the analysis toolkit for malaria population dynamics. Compared to existing software, MOIRE enhances the accuracy of parameter estimation and enables more comprehensive insights into within-host diversity and population structure. Additionally, MOIRE's adaptability to diverse data sources and potential for future improvements make it a valuable asset for research on malaria and other organisms, such as other eukaryotic pathogens. MOIRE is available as an R package at https://eppicenter.github.io/moire/.

Interleukin-13 (IL-13)-A Pleiotropic Cytokine Involved in Wound Healing and Fibrosis

The liver, as a central metabolic organ, is systemically linked to metabolic-inflammatory diseases. In the pathogenesis of the metabolic syndrome, inflammatory and metabolic interactions between the intestine, liver, and adipose tissue lead to the progression of hepatic steatosis to metabolic-dysfunction-associated steatohepatitis (MASH) and consecutive MASH-induced fibrosis. Clinical and animal studies revealed that IL-13 might be protective in the development of MASH through both the preservation of metabolic functions and Th2-polarized inflammation in the liver and the adipose tissue. In contrast, IL-13-associated loss of mucosal gut barrier function and IL-13-associated enhanced hepatic fibrosis may contribute to the progression of MASH. However, there are only a few publications on the effect of IL-13 on metabolic diseases and possible therapies to influence them. In this review article, different aspects of IL-13-associated effects on the liver and metabolic liver diseases, which are partly contradictory, are summarized and discussed on the basis of the recent literature.

Geographical Influence on Morphometric Variability of Genetically “Pure” Schistosoma haematobium Eggs from Sub-Saharan Migrants in Spain

Schistosome eggs play a key role in schistosomiasis diagnosis and research. The aim of this work is to morphogenetically study the eggs of Schistosoma haematobium found in sub-Saharan migrants present in Spain, analyzing their morphometric variation in relation to the geographical origin of the parasite (Mali, Mauritania and Senegal). Only eggs considered “pure” S. haematobium by genetic characterization (rDNA ITS-2 and mtDNA cox1) have been used. A total of 162 eggs obtained from 20 migrants from Mali, Mauritania and Senegal were included in the study. Analyses were made by the Computer Image Analysis System (CIAS). Following a previously standardized methodology, seventeen measurements were carried out on each egg. The morphometric analysis of the three morphotypes detected (round, elongated and spindle) and the biometric variations in relation to the country of origin of the parasite on the egg phenotype were carried out by canonical variate analysis. Mahalanobis distances, when all egg measurements were analyzed, showed differences between: (i) Mali-Mauritania, Mali-Senegal and Mauritania-Senegal in the round morphotype; (ii) Mali-Mauritania and Mauritania-Senegal in the elongated morphotype; and (iii) Mauritania-Senegal in the spindle morphotype. Mahalanobis distances, when spine variables were analyzed, showed differences between Mali-Senegal in the round morphotype. In conclusion, this is the first phenotypic study performed on individually genotyped “pure” S. haematobium eggs, allowing the assessment of the intraspecific morphological variations associated with the geographical origin of the schistosome eggs.

Population genetic structure of Schistosoma haematobium and Schistosoma haematobium × Schistosoma bovis hybrids among school-aged children in Côte d'Ivoire

While population genetics of Schistosoma haematobium have been investigated in West Africa, only scant data are available from Côte d’Ivoire. The purpose of this study was to analyze both genetic variability and genetic structure among S. haematobium populations and to quantify the frequency of S. haematobium × S. bovis hybrids in school-aged children in different parts of Côte d’Ivoire. Urine samples were subjected to a filtration method and examined microscopically for Schistosoma eggs in four sites in the western and southern parts of Côte d’Ivoire. A total of 2692 miracidia were collected individually and stored on Whatman^® FTA cards. Of these, 2561 miracidia were successfully genotyped for species and hybrid identification using rapid diagnostic multiplex mitochondrial cox1 PCR and PCR Restriction Fragment Length Polymorphism (PCR-RFLP) analysis of the nuclear ITS2 region. From 2164 miracidia, 1966 (90.9%) were successfully genotyped using at least 10 nuclear microsatellite loci to investigate genetic diversity and population structure. Significant differences were found between sites in all genetic diversity indices and genotypic differentiation was observed between the site in the West and the three sites in the East. Analysis at the infrapopulation level revealed clustering of parasite genotypes within individual children, particularly in Duekoué (West) and Sikensi (East). Of the six possible cox1-ITS2 genetic profiles obtained from miracidia, S. bovis cox1 × S. haematobium ITS2 (42.0%) was the most commonly observed in the populations. We identified only 15 miracidia (0.7%) with an S. bovis cox1 × S. bovis ITS2 genotype. Our study provides new insights into the population genetics of S. haematobium and S. haematobium × S. bovis hybrids in humans in Côte d’Ivoire and we advocate for researching hybrid schistosomes in animals such as rodents and cattle in Côte d’Ivoire.

Population Genetic Structure and Hybridization of Schistosoma haematobium in Nigeria

Background: Schistosomiasis is a major poverty-related disease caused by dioecious parasitic flatworms of the genus Schistosoma with a health impact on both humans and animals. Hybrids of human urogenital schistosome and bovine intestinal schistosome have been reported in humans in several of Nigeria’s neighboring West African countries. No empirical studies have been carried out on the genomic diversity of Schistosoma haematobium in Nigeria. Here, we present novel data on the presence and prevalence of hybrids and the population genetic structure of S. haematobium. Methods: 165 Schistosoma-positive urine samples were obtained from 12 sampling sites in Nigeria. Schistosoma haematobium eggs from each sample were hatched and each individual miracidium was picked and preserved in Whatman® FTA cards for genomic analysis. Approximately 1364 parasites were molecularly characterized by rapid diagnostic multiplex polymerase chain reaction (RD-PCR) for mitochondrial DNA gene (Cox1 mtDNA) and a subset of 1136 miracidia were genotyped using a panel of 18 microsatellite markers. Results: No significant difference was observed in the population genetic diversity (p > 0.05), though a significant difference was observed in the allelic richness of the sites except sites 7, 8, and 9 (p < 0.05). Moreover, we observed two clusters of populations: west (populations 1−4) and east (populations 7−12). Of the 1364 miracidia genotyped, 1212 (89%) showed an S. bovis Cox1 profile and 152 (11%) showed an S. haematobium cox1 profile. All parasites showed an S. bovis Cox1 profile except for some at sites 3 and 4. Schistosoma miracidia full genotyping showed 59.3% of the S. bovis ITS2 allele. Conclusions: This study provides novel insight into hybridization and population genetic structure of S. haematobium in Nigeria. Our findings suggest that S. haematobium x S. bovis hybrids are common in Nigeria. More genomic studies on both human- and animal-infecting parasites are needed to ascertain the role of animals in schistosome transmission.

Cryptic population structure and transmission dynamics uncovered for Schistosoma mansoni populations by genetic analyses

Patterns of diversity in pathogen genomes provide a window into the spatiotemporal spread of disease. In this study, we tested the hypothesis that Schistosoma mansoni parasites form genetic clusters that coincide with the communities of their human hosts. We also looked for genetic clustering of parasites at the sub-community level. Our data consists of 14 microsatellite DNA markers, typed from pooled DNA samples from [Formula: see text] infected individuals living in three Brazilian communities. We found a one-to-one correspondence between genetic clusters found by K-means cluster analysis and communities when [Formula: see text]. These clusters are also easily identified in a neighbor-joining tree and principal coordinates plots. K-means analysis with [Formula: see text] also reveals genetic clusters of parasites at the sub-community level. These sub-clusters also appear on the neighbor-joining tree and principal coordinates plots. A surprising finding is a genetic relationship between subgroups in widely separated human communities. This connection suggests the existence of common transmission sites that have wide influence. In summary, the genetic structure of S. mansoni in Brazil juxtaposes local isolation that is occasionally broken by long-range migration. Permanent eradication of schistosomes will require both local efforts and the identification of regional infection reservoirs.

A Schistosoma mansoni tri- and tetramer microsatellite catalog for genetic population diversity and differentiation

All Schistosoma mansoni tri- and tetranucleotide repeat microsatellites published as of December 2018 were identified. All 52 were evaluated for autosomal location, strength of amplification, scorability and behavior as single-copy loci by polyacrylamide and capillary gel electrophoresis. Of these, 27 were unique, autosomal, polymorphic, easily scored and single copy as assessed on pooled adult worm DNA from two different continental origins and adult worm clones. These microsatellites were distributed across all seven autosomal chromosomes. On laboratory strains their heterozygosity ranged from 0.22 to 0.77. Individual markers had 5-13 alleles, allelic richness of 2-10 and an effective allele number of 1.3-8.14. Those infected by Schistosoma mansoni carry many genetically distinct, sexually reproducing parasites, therefore, for an individual infection the complete allele frequency profile of their progeny consists of a pool of DNA from multiple diploid eggs. Using a set of 25 microsatellites, we calculated allele frequency profiles of eggs in fecal samples from people in two Brazilian communities separated by 6 km: Jenipapo (n = 80) and Volta do Rio (n = 38). There were no a priori characteristics that could predict the performance of markers in natural infections based on their performance with laboratory strains. Increasing marker number did not change accuracy for differentiation and diversity but did improve precision. Our data suggest that using a random set of 10-20 microsatellites appears to result in values that exhibit low standard deviations for diversity and differentiation indices. All identified microsatellites as well as PCR conditions, allele size, primer sequences and references for all tri- and tetramer microsatellites markers presented in this work are available at: https://sites.google.com/case.edu/cwru-and-fiocruz-wdrc/home.

Whole-genome sequencing of Schistosoma mansoni reveals extensive diversity with limited selection despite mass drug administration

Control and elimination of the parasitic disease schistosomiasis relies on mass administration of praziquantel. Whilst these programmes reduce infection prevalence and intensity, their impact on parasite transmission and evolution is poorly understood. Here we examine the genomic impact of repeated mass drug administration on Schistosoma mansoni populations with documented reduced praziquantel efficacy. We sequenced whole-genomes of 198 S. mansoni larvae from 34 Ugandan children from regions with contrasting praziquantel exposure. Parasites infecting children from Lake Victoria, a transmission hotspot, form a diverse panmictic population. A single round of treatment did not reduce this diversity with no apparent population contraction caused by long-term praziquantel use. We find evidence of positive selection acting on members of gene families previously implicated in praziquantel action, but detect no high frequency functionally impactful variants. As efforts to eliminate schistosomiasis intensify, our study provides a foundation for genomic surveillance of this major human parasite.

Revisiting density-dependent fecundity in schistosomes using sibship reconstruction

The stability of parasite populations is regulated by density-dependent processes occurring at different stages of their life cycle. In dioecious helminth infections, density-dependent fecundity is one such regulatory process that describes the reduction in egg production by female worms in high worm burden within-host environments. In human schistosomiasis, the operation of density-dependent fecundity is equivocal and investigation is hampered by the inaccessibility of adult worms that are located intravascularly. Current understanding is almost exclusively limited to data collected from two human autopsy studies conducted over 40 years ago, with subsequent analyses having reached conflicting conclusions. Whether egg production is regulated in a density-dependent manner is key to predicting the effectiveness of interventions targeting the elimination of schistosomiasis and to the interpretation of parasitological data collected during monitoring and evaluation activities. Here, we revisit density-dependent fecundity in the two most globally important human Schistosoma spp. using a statistical modelling approach that combines molecular inference on the number of parents/adult worms in individual human hosts with parasitological egg count data from mainland Tanzania and Zanzibar. We find a non-proportional relationship between S. haematobium egg counts and inferred numbers of female worms, providing the first clear evidence of density-dependent fecundity in this schistosome species. We do not find robust evidence for density-dependent fecundity in S. mansoni because of high sensitivity to some modelling assumptions and the lower statistical power of the available data. We discuss the strengths and limitations of our model-based analytical approach and its potential for improving our understanding of density dependence in schistosomiasis and other human helminthiases earmarked for elimination.

Schistosomiasis is a devastating disease of poverty currently estimated to infect over 220 million people. It is caused by parasitic worms (blood flukes) that live for, on average, 5–7 years inside the blood vessels of infected hosts and produce hundreds of eggs daily. Whether egg production is regulated in a density-dependent manner, and if so under what conditions, has been controversial for schistosomiasis, and investigation is hampered due to the inaccessible location of adult worms. Resolving this fundamental question is important because density dependencies determine the resilience of helminthiases to interventions. Here, we have revisited this longstanding and unresolved question of density-dependent fecundity in human schistosomes using a novel statistical modelling approach that combines information from molecular and parasitological data. We report the first clear evidence of density-dependent fecundity in S. haematobium, the causative agent of millions of cases of urogenital schistosomiasis. Our findings are of critical importance both to mathematical modellers predicting the impact of interventions and to public health policy makers striving to meet the 2030 elimination targets for schistosomiasis. This study also serves to illustrate a new biostatistical approach that could be applied to investigate density dependencies in other helminthiases where adult parasites are inaccessible.

Schistosomiasis Morbidity Hotspots: Roles of the Human Host, the Parasite and Their Interface in the Development of Severe Morbidity

Schistosomiasis is the second most important human parasitic disease in terms of socioeconomic impact, causing great morbidity and mortality, predominantly across the African continent. For intestinal schistosomiasis, severe morbidity manifests as periportal fibrosis (PPF) in which large tracts of macro-fibrosis of the liver, visible by ultrasound, can occlude the main portal vein leading to portal hypertension (PHT), sequelae such as ascites and collateral vasculature, and ultimately fatalities. For urogenital schistosomiasis, severe morbidity manifests as pathology throughout the urinary system and genitals, and is a definitive cause of squamous cell bladder carcinoma. Preventative chemotherapy (PC) programmes, delivered through mass drug administration (MDA) of praziquantel (PZQ), have been at the forefront of schistosomiasis control programmes in sub-Saharan Africa since their commencement in Uganda in 2003. However, despite many successes, 'biological hotspots' (as distinct from 'operational hotspots') of both persistent high transmission and morbidity remain. In some areas, this failure to gain control of schistosomiasis has devastating consequences, with not only persistently high infection intensities, but both "subtle" and severe morbidity remaining prevalent. These hotspots highlight the requirement to revisit research into severe morbidity and its mechanisms, a topic that has been out of favor during times of PC implementation. Indeed, the focality and spatially-structured epidemiology of schistosomiasis, its transmission persistence and the morbidity induced, has long suggested that gene-environmental-interactions playing out at the host-parasite interface are crucial. Here we review evidence of potential unique parasite factors, host factors, and their gene-environmental interactions in terms of explaining differential morbidity profiles in the human host. We then take the situation of schistosomiasis mansoni within the Albertine region of Uganda as a case study in terms of elucidating the factors behind the severe morbidity observed and the avenues and directions for future research currently underway within a new research and clinical trial programme (FibroScHot).

Population genetics of African Schistosoma species

Blood flukes within the genus Schistosoma (schistosomes) are responsible for the major disease, schistosomiasis, in tropical and sub-tropical areas. This disease is predominantly present on the African continent with more than 85% of the human cases. Schistosomes are also parasites of veterinary importance infecting livestock and wildlife. Schistosoma population genetic structure and diversity are important characteristics that may reflect variations in selection pressures such as those induced by host (mammalian and snail) environments, habitat change, migration and also treatment/control interventions, all of which also shape speciation and evolution of the whole Schistosoma genus. Investigations into schistosome population genetic structure, diversity and evolution has been an area of important debate and research. Supported by advances in molecular techniques with capabilities for multi-locus genetic analyses for single larvae schistosome genetic investigations have greatly progressed in the last decade. This paper aims to review the genetic studies of both animal and human infecting schistosome. Population genetic structures are reviewed at different spatial scales: local, regional or continental (i.e. phylogeography). Within species genetic diversities are discussed compared and the compounding factors discussed, including the effect of mass drug administration. Finally, the ability for intra-species hybridisation questions species integrities and poses many questions in relation to the natural epidemiology of co-endemic species. Here we review molecularly confirmed hybridisation events (in relation to human disease) and discuss the possible impact for ongoing and future control and elimination.

Genetic evidence for the role of non-human primates as reservoir hosts for human schistosomiasis

Background

Schistosomiasis is a chronic parasitic disease, that affects over 207 million people and causes over 200,000 deaths annually, mainly in sub-Saharan Africa. Although many health measures have been carried out to limit parasite transmission, significant numbers of non-human primates such as Chlorocebus aethiops (Ch. aethiops) (vervet) and Papio anubis (baboon) are infected with S. mansoni, notably in Ethiopia, where they are expected to have potentially significant implications for transmission and control efforts.

Objective

The objective of this study was to assess and compare the genetic diversity and population structure of S. mansoni isolates from human and non-human primates free-ranging in close proximity to villages in selected endemic areas of Ethiopia.

Methods

A cross-sectional study was conducted in three transmission sites: Bochesa, Kime and Fincha. A total of 2,356 S. mansoni miracidia were directly isolated from fecal specimens of 104 hosts (i.e. 60 human hosts and 44 non-human primates). We performed DNA extraction and PCR amplification using fourteen microsatellite loci.

Results

At population scale we showed strong genetic structure between the three sample sites. At the definitive host scale, we observed that host factors can shape the genetic composition of parasite infra-populations. First, in male patients, we observed a positive link between parasite genetic diversity and the age of the patients. Second, we observed a difference in genetic diversity which was high in human males, medium in human females and low in non-human primates (NHPs). Finally, whatever the transmission site no genetic structure was observed between human and non-human primates, however, there appears to be little barriers, if any, host specificity of the S. mansoni populations with cross-host infections.

Conclusion

Occurrence of infection of a single host with multiple S. mansoni strains and inter- and intra-host genetic variations was observed. Substantial genetic diversity and gene flow across the S. mansoni population occurred at each site and non-human primates likely play a role in local transmission and maintenance of infection. Therefore, public health and wildlife professionals should work together to improve disease control and elimination strategies.

Schistosomiasis is a chronic disease caused by flukes (trematodes). The definitive host spectrum of schistosomes, whether human, non-human primates (NHPs) or other mammals, is highly dependent on the schistosome species concerned. Genetic diversity and population structure studies of S. mansoni have provided insights into the variation of natural populations. Understanding S. mansoni genetic diversity and population structure of isolates from human and non-human primate hosts living in close proximity showed the occurrence of infection of a single host with multiple S. mansoni strains and inter- and intra-host genetic variations. In this article, the researchers assert the fact that genetic approach reveals that parasites from the three different sites are independent. Thus, we could consider the three sites as geographical replicates showing the influence of NHPs in parasitic transmission in Ethiopia. This study provides insights into the epidemiology, genetic diversity and population structure of S. mansoni in human and non-human primates in Ethiopia, all of which are crucial for the control of schistosomiasis.

Population Structure and Dynamics of Helminthic Infection: Schistosomiasis

While disease and outbreaks are mainly clonal for bacteria and other asexually reproducing organisms, sexual reproduction in schistosomes and other helminths usually results in unique individuals. For sexually reproducing organisms, the traits conserved in clones will instead be conserved in the group of organisms that tends to breed together, the population. While the same tools are applied to characterize DNA, how results are interpreted can be quite different at times (see another article in this collection, http://www.asmscience.org/content/journal/microbiolspec/10.1128/microbiolspec.AME-0002-2018). It is difficult to know what the real effect any control program has on the parasite population without assessing the health of this population, how they respond to the control measure, and how they recover, if they do. This review, part of the Microbiology Spectrum Curated Collection: Advances in Molecular Epidemiology of Infectious Diseases, concentrates on one approach using pooled samples to study schistosome populations and shows how this and other approaches have contributed to our understanding of this parasite family's biology and epidemiology. *This article is part of a curated collection.

Two-year longitudinal survey reveals high genetic diversity of Schistosoma mansoni with adult worms surviving praziquantel treatment at the start of mass drug administration in Uganda

BACKGROUND

A key component of schistosomiasis control is mass drug administration with praziquantel. While control interventions have been successful in several endemic regions, mass drug administration has been less effective in others. Here we focus on the impact of repeated praziquantel treatment on the population structure and genetic diversity of Schistosoma mansoni.

METHODS

We examined S. mansoni epidemiology, population genetics, and variation in praziquantel susceptibility in parasites isolated from children across three primary schools in a high endemicity region at the onset of the Ugandan National Control Programme. Children were sampled at 11 timepoints over two years, including one week and four weeks post-praziquantel treatment to evaluate short-term impacts on clearance and evidence of natural variation in susceptibility to praziquantel.

RESULTS

Prevalence of S. mansoni was 85% at baseline. A total of 3576 miracidia larval parasites, isolated from 203 individual children, were genotyped at seven loci. Overall, genetic diversity was high and there was low genetic differentiation, indicating high rates of parasite gene flow. Schistosome siblings were found both pre-treatment and four weeks post-treatment, demonstrating adult worms surviving treatment and natural praziquantel susceptibility variation in these populations at the beginning of mass drug administration. However, we did not find evidence for selection on these parasites. While genetic diversity decreased in the short-term (four weeks post-treatment), diversity did not decrease over the entire period despite four rounds of mass treatment. Furthermore, within-host genetic diversity was affected by host age, host sex, infection intensity and recent praziquantel treatment.

CONCLUSIONS

Our findings suggest that praziquantel treatments have short-term impacts on these parasite populations but impacts were transient and no long-term reduction in genetic diversity was observed. High gene flow reduces the likelihood of local adaptation, so even though parasites surviving treatment were observed, these were likely to be diluted at the beginning of the Ugandan National Control Programme. Together, these results suggest that MDA in isolation may be insufficient to reduce schistosome populations in regions with high genetic diversity and gene flow.

Estimating helminth burdens using sibship reconstruction

BACKGROUND

Sibship reconstruction is a form of parentage analysis that can be used to identify the number of helminth parental genotypes infecting individual hosts using genetic data on only their offspring. This has the potential to be used for estimating individual worm burdens when adult parasites are otherwise inaccessible, the case for many of the most globally important human helminthiases and neglected tropical diseases. Yet methods of inferring worm burdens from sibship reconstruction data on numbers of unique parental genotypes are lacking, limiting the method's scope of application.

RESULTS

We developed a novel statistical method for estimating female worm burdens from data on the number of unique female parental genotypes derived from sibship reconstruction. We illustrate the approach using genotypic data on Schistosoma mansoni (miracidial) offspring collected from schoolchildren in Tanzania. We show how the bias and precision of worm burden estimates critically depends on the number of sampled offspring and we discuss strategies for obtaining sufficient sample sizes and for incorporating judiciously formulated prior information to improve the accuracy of estimates.

CONCLUSIONS

This work provides a novel approach for estimating individual-level worm burdens using genetic data on helminth offspring. This represents a step towards a wider scope of application of parentage analysis techniques. We discuss how the method could be used to assist in the interpretation of monitoring and evaluation data collected during mass drug administration programmes targeting human helminthiases and to help resolve outstanding questions on key population biological processes that govern the transmission dynamics of these neglected tropical diseases.

Population genetic structure of Schistosoma bovis in Cameroon

Background

Schistosomiasis is neglected tropical parasitic disease affecting both humans and animals. Due to the human health impact, population genetic studies have focused on the three main human-infecting schistosome species: Schistosoma mansoni, S. haematobium and S. japonicum. Here we present novel data on the population genetic structure of Schistosoma bovis, a highly widespread and prevalent schistosome infecting ruminants, and therefore of veterinary importance.

Methods

Adult S. bovis were sampled in the two main abattoirs of Cameroon (Yaoundé and Douala). Twenty-two cows originating from four distinct localities were sampled and a total of 218 parasites were recovered. All parasites were genotyped using a panel of 14 microsatellite markers and a sub-sample of 91 parasites were sequenced and characterized with the mitochondrial (cox1) and nuclear (ITS) genetic markers.

Results

No significant difference in allelic richness, heterozygosity, nucleotide diversity and haplotype diversity was observed between the populations. Additionally, no strong genetic structure was observed at the country scale. Our data also show that S. bovis is more polymorphic than its sister species, S. haematobium, and that the haplotype diversity is similar to that of S. mansoni while the nucleotide diversity does not significantly differ from that of S. haematobium. The resulting negative Tajima’s D* and Fu and Li’s D* indices could be a signature of population demographic expansion. No S. haematobium/S. bovis hybrids were observed in our populations, thus all samples were considered as pure S. bovis.

Conclusions

This study provides novel insights into genetic diversity and population genetic structure of S. bovis. No strong genetic structure was observed at the country scale but some genetic indices could be associated as a signature of population demographic expansion.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3307-0) contains supplementary material, which is available to authorized users.

Whole genome amplification and exome sequencing of archived schistosome miracidia

Adult schistosomes live in the blood vessels and cannot easily be sampled from humans, so archived miracidia larvae hatched from eggs expelled in feces or urine are commonly used for population genetic studies. Large collections of archived miracidia on FTA cards are now available through the Schistosomiasis Collection at the Natural History Museum (SCAN). Here we describe protocols for whole genome amplification of Schistosoma mansoni and Schistosome haematobium miracidia from these cards, as well as real time PCR quantification of amplified schistosome DNA. We used microgram quantities of DNA obtained for exome capture and sequencing of single miracidia, generating dense polymorphism data across the exome. These methods will facilitate the transition from population genetics, using limited numbers of markers to population genomics using genome-wide marker information, maximising the value of collections such as SCAN.

Phenotypic and genotypic monitoring of Schistosoma mansoni in Tanzanian schoolchildren five years into a preventative chemotherapy national control programme

BACKGROUND

Schistosoma mansoni is a parasite of profound medical importance. Current control focusses on mass praziquantel (PZQ) treatment of populations in endemic areas, termed Preventative Chemotherapy (PC). Large-scale PC programmes exert prolonged selection pressures on parasites with the potential for, direct and/or indirect, emergence of drug resistance. Molecular methods can help monitor genetic changes of schistosome populations over time and in response to drug treatment, as well as estimate adult worm burdens through parentage analysis. Furthermore, methods such as in vitro drug sensitivity assays help phenotype in vivo parasite genotypic drug efficacy.

METHODS

We conducted combined in vitro PZQ efficacy testing with population genetic analyses of S. mansoni collected from children from two schools in 2010, five years after the introduction of a National Control Programme. Children at one school had received four annual PZQ treatments and the other school had received two mass treatments in total. We compared genetic differentiation, indices of genetic diversity, and estimated adult worm burden from parasites collected in 2010 with samples collected in 2005 (before the control programme began) and in 2006 (six months after the first PZQ treatment). Using 2010 larval samples, we also compared the genetic similarity of those with high and low in vitro sensitivity to PZQ.

RESULTS

We demonstrated that there were individual parasites with reduced PZQ susceptibility in the 2010 collections, as evidenced by our in vitro larval behavioural phenotypic assay. There was no evidence, however, that miracidia showing phenotypically reduced susceptibility clustered together genetically. Molecular analysis also demonstrated a significant reduction of adult worm load over time, despite little evidence of reduction in parasite infection intensity, as measured by egg output. Genetic diversity of infections did not reduce over time, despite changes in the genetic composition of the parasite populations.

CONCLUSIONS

Genotypic and phenotypic monitoring did not indicate a selective sweep, as may be expected if PZQ treatment was selecting a small number of related "resistant" parasites, but there was evidence of genetic changes at the population level over time. Genetic data were used to estimate adult worm burdens, which unlike parasite infection intensity, showed reductions over time, suggesting the relaxation of negative density-dependent constraints on parasite fecundity with PZQ treatment. We thereby demonstrated that density-dependence in schistosome populations may complicate evaluation and monitoring of control programmes.

Geographic strain differentiation of Schistosoma japonicum in the Philippines using microsatellite markers

BACKGROUND

Microsatellites have been found to be useful in determining genetic diversities of various medically-important parasites which can be used as basis for an effective disease management and control program. In Asia and Africa, the identification of different geographical strains of Schistosoma japonicum, S. haematobium and S. mansoni as determined through microsatellites could pave the way for a better understanding of the transmission epidemiology of the parasite. Thus, the present study aims to apply microsatellite markers in analyzing the populations of S. japonicum from different endemic areas in the Philippines for possible strain differentiation.

METHODOLOGY/ PRINCIPAL FINDINGS

Experimental mice were infected using the cercariae of S. japonicum collected from infected Oncomelania hupensis quadrasi snails in seven endemic municipalities. Adult worms were harvested from infected mice after 45 days of infection and their DNA analyzed against ten previously characterized microsatellite loci. High genetic diversity was observed in areas with high endemicity. The degree of genetic differentiation of the parasite population between endemic areas varies. Geographical separation was considered as one of the factors accounting for the observed difference between populations. Two subgroups have been observed in one of the study sites, suggesting that co-infection with several genotypes of the parasite might be present in the population. Clustering analysis showed no particular spatial structuring between parasite populations from different endemic areas. This result could possibly suggest varying degrees of effects of the ongoing control programs and the existing gene flow in the populations, which might be attributed to migration and active movement of infected hosts from one endemic area to another.

CONCLUSIONS/ SIGNIFICANCE

Based on the results of the study, it is reasonable to conclude that genetic diversity could be one possible criterion to assess the infection status in highly endemic areas. Genetic surveillance using microsatellites is therefore important to predict the ongoing gene flow and degree of genetic diversity, which indirectly reflects the success of the control program in schistosomiasis-endemic areas.

Genetic differences among Haplorchis taichui populations in Indochina revealed by mitochondrial COX1 sequences

Haplorchis taichui is an intestinal heterophyid fluke that is pathogenic to humans. It is widely distributed in Asia, with a particularly high prevalence in Indochina. Previous work revealed that the lack of gene flow between three distinct populations of Vietnamese H. taichui can be attributed to their geographic isolation with no interconnected river basins. To test the hypothesis that interconnected river basins allow gene flow between otherwise isolated populations of H. taichui, as previously demonstrated for another trematode, Opisthorchis viverrini, we compared the genetic structures of seven populations of H. taichui from various localities in the lower Mekong Basin, in Thailand and Laos, with those in Vietnam, using the mitochondrial cytochrome c oxidase subunit 1 (COX1) gene. To determine the gene flow between these H. taichui populations, we calculated their phylogenetic relationships, genetic distances and haplotype diversity. Each population showed very low nucleotide diversity at this locus. However, high levels of genetic differentiation between the populations indicated very little gene flow. A phylogenetic analysis divided the populations into four clusters that correlated with the country of origin. The negligible gene flow between the Thai and Laos populations, despite sharing the Mekong Basin, caused us to reject our hypothesis. Our data suggest that the distribution of H. taichui populations was incidentally associated with national borders.