Population genomics reveals multi-scale mechanisms sustaining schistosomiasis re-emergence in a near-elimination setting

In China, sustained snail control, environmental management, and mass drug administration with praziquantel reduced schistosomiasis to near-elimination levels, yet re-emergence in Sichuan Province during the early 2000s exposed vulnerabilities in late-stage control. We conducted a novel investigation of Schistosoma japonicum re-emergence in Sichuan in order to identify the processes underlying re-emergence and to showcase how genomic data can be used in such investigations. We sequenced whole genomes from 270 miracidia collected from 53 human hosts across 17 villages in 2007 – the year after re-emergence was documented. Population genomic analyses revealed a broadly cohesive regional parasite population with weak geographic structuring. Genome-wide diversity was substantial, and demographic reconstructions revealed no recent decline in effective population size, demonstrating that parasite populations were not demographically fragmented at the onset of re-emergence, hinting at maintenance in reservoir populations. At finer spatial scales, several villages exhibited low diversity and elevated inbreeding, consistent with localized transmission maintained by small founding populations. Estimates of pairwise genetic relatedness revealed dense within-village sibling clusters alongside second- and third-degree relationships spanning villages, and rare first- and second-degree cross-village links. These findings are consistent with highly focal transmission with episodic parasite dispersal across villages, leading to a regional transmission network. Inferred minimum worm burden across hosts varied from 1 to 11 adult worm pairs, indicating heterogeneity in within-host parasite diversity, although uneven sampling limited inference. Together, these findings indicate that schistosomiasis re-emergence in this near-elimination setting was likely facilitated by a diverse parasite population maintained in reservoir populations, and that transmission, while predominantly local was occurring across a network of connected villages. This work illustrates how population genomics can reveal mechanisms driving re-emergence in late-stage elimination in complex, multi-host transmission systems.