A genomic epidemiology investigation of yaws re-emergence and bacterial drug resistance selection
Background: In a longitudinal study assessing the WHO strategy for yaws eradication using mass azithromycin treatment, we observed resurgence of yaws cases with dominance of a single JG8 sequence type and emergence of azithromycin-resistant Treponema pallidum. Here, we analyse genomic changes in the bacterial population using samples collected during the study. Methods: We performed whole bacterial genome sequencing directly on DNA extracted from 37 lesion swabs collected from patients on Lihir Island, Papua New Guinea, between 2013 and 2016. We produced phylogenies and correlated these with temporo-spatial information to investigate the source of new cases and the emergence of five macrolide-resistant cases. We used deep amplicon sequencing of surveillance samples to assess the presence of minority macrolide resistance populations. Findings: We recovered 20 whole Treponema genomes, and phylogenetic analysis showed that the re-emerging JG8 sequence type was composed of three bacterial sub-lineages characterised by distinct temporo-spatial patterns. Of five patients with resistant Treponema, all epidemiologically linked, we recovered genomes from three and found no variants. Deep sequencing showed that pre-treatment, the index patient harboured fixed macrolide-sensitive Treponema, while the post-treatment sample harboured a fixed resistant genotype, as did three of four contact cases. We also found no evidence of pre-existing minority Treponema drug-resistance variants in the general population. Interpretation: In this study, re-emergence of yaws cases was polyphyletic, indicating multiple epidemiological sources. However, given the genomic and epidemiological linkage of resistant cases and the rarity of resistance alleles in the general population, it is likely that azithromycin resistance evolved only once in this study, followed by onward dissemination.