Applied phyloepidemiology: Detecting drivers of pathogen transmission from genomic signatures using density measures

Understanding the driving forces of an epidemic is key to inform intervention strategies against it. Correlating measures of the epidemic success of a pathogen with ancillary parameters such as its drug resistance profile provides a flexible tool to identify such driving forces. The recently described time-scaled haplotypic density (THD) method facilitates the inference of a pathogen's epidemic success from genetic data. Contrary to demogenetic approaches that define success in an aggregated fashion, the THD computes an independent index of success for each isolate in a collection. Modeling this index using multivariate regression, thus, allows us to control for various sources of bias and to identify independent predictors of success. We illustrate the use of THD to address key questions regarding three exemplary epidemics of multidrug-resistant (MDR) bacterial lineages, namely Mycobacterium tuberculosis Beijing, Salmonella Typhi H58, and Staphylococcus aureus ST8 (including ST8-USA300 MRSA), based on previously published, international genetic datasets. In each case, THD analysis allowed to identify the impact, or lack thereof, of various factors on the epidemic success, independent of confounding by population structure and geographic distribution. Our results suggest that rifampicin resistance drives the MDR Beijing epidemic and that fluoroquinolone resistance drives the S. aureus ST8/USA300 epidemic, in line with previous evidence of a lack of resistance-associated fitness cost in these pathogens. Conversely, fluoroquinolone resistance measurably hampered the success of S. Typhi H58 and non-H58. These findings illustrate how THD can help leverage the massive genomic datasets generated by molecular epidemiology studies to address new questions. THD implementation for the R platform is available at https://github.com/rasigadelab/thd.