Population-level genome sequencing reveals distinct Mycobacterium tuberculosis intrahost mutational trajectories in simian immunodeficiency virus co-infected and antiretroviral treated non-human primates

Whole genome sequencing of Mycobacterium tuberculosis (Mtb) populations from clinical samples has increasingly identified genes undergoing selection within and between hosts that drive differential infection and treatment outcomes. However, the intrahost Mtb mutational landscape—especially in the context of human immunodeficiency virus coinfection and antiretroviral therapy (ART)—remains less clear, as do the potential impacts of such mutations on Mtb infection dynamics. Here, we performed whole genome sequencing of Mtb populations isolated from approximately 480 infected tissues across 20 non-human primates, including animals co-infected with simian immunodeficiency virus (SIV) with or without virological suppression by ART. We identified 116 mutations that emerged during infection, including those that are overrepresented within individual tissues and a subset that are shared across tissues during Mtb dissemination. We further find differential mutation trajectories across treatment groups, with higher mutation rate and bacterial outgrowth in SIV-infected hosts and increased prevalence of oxidative damage-associated mutations in coinfected animals on ART. Finally, we demonstrate a common pattern of mutation in Mtb lipid metabolism and polyketide synthase genes and identify a subset of NHP-derived mutations that have also independently arisen in human clinical isolates. Together, our population-based sequencing uncovers Mtb diversification during early infection, captures discrete bacterial dissemination events and infers differential immune pressures faced by Mtb in the setting of SIV-Mtb coinfection and ART suppression.