Genome-wide analysis of Schistosoma mansoni reveals population structure and praziquantel drug selection pressure within Ugandan hot-spot communities

Abstract Populations within schistosomiasis control areas, especially those in Africa, are recommended to receive regular mass drug administration (MDA) with praziquantel (PZQ) as the main strategy for controlling the disease. The impact of PZQ treatment on schistosome genetics remains poorly understood, and is limited by a lack of high-resolution genetic data on the population structure of parasites within these control areas. We generated whole-genome sequence data from 174 individual miracidia collected from both children and adults from fishing communities on islands in Lake Victoria in Uganda that had received either annual or quarterly MDA with PZQ over four years, including samples collected immediately before and four weeks after treatment. Genome variation within and between samples was characterised and we investigated genomic signatures of natural selection acting on these populations that could be due to PZQ treatment. The parasite population on these islands was more diverse than found in nearby villages on the lake shore. We saw little or no genetic differentiation between villages, or between the groups of villages with different treatment intensity, but slightly higher genetic diversity within the pre-treatment compared to post-treatment parasite populations. We identified classes of genes significantly enriched within regions of the genome with evidence of recent positive selection among post-treatment and intensively treated parasite populations. The differential selection observed in post-treatment and pre-treatment parasite populations could be linked to any reduced susceptibility of parasites to praziquantel treatment. Author summary Schistosomiasis is caused by parasitic helminths of the genus Schistosoma . Schistosoma mansoni is the primary cause of intestinal schistosomiasis, a devastating and widespread parasitic infection that causes morbidity, death and socio-economic impact on endemic communities across the world and especially sub-Saharan Africa. Using whole-genome sequencing, we were able to elucidate the parasite population within Lake Victoria island fishing communities in Uganda which are among the major hotspots for schistosomiasis. We further assessed genetic markers that might be linked to recent observations concerning reduced susceptibility to praziquantel, the major drug used in the treatment of this disease. Whole-genome data on the population genetics of S. mansoni in an African setting will provide a strong basis for future functional genomics or transcriptomic studies that will be key to identifying drug targets, improving existing drugs or developing new therapeutic interventions.