Sebastien Gagneux and colleagues analyze a global collection of Mycobacterium tuberculosis clinical isolates to classify sublineages by phylogeography. They find globally distributed ‘generalist’ and geographically restricted ‘specialist’ sublineages of lineage 4, indicating that different evolutionary strategies were adopted to succeed in various ecological niches. Generalist and specialist species differ in the breadth of their ecological niches. Little is known about the niche width of obligate human pathogens. Here we analyzed a global collection of Mycobacterium tuberculosis lineage 4 clinical isolates, the most geographically widespread cause of human tuberculosis. We show that lineage 4 comprises globally distributed and geographically restricted sublineages, suggesting a distinction between generalists and specialists. Population genomic analyses showed that, whereas the majority of human T cell epitopes were conserved in all sublineages, the proportion of variable epitopes was higher in generalists. Our data further support a European origin for the most common generalist sublineage. Hence, the global success of lineage 4 reflects distinct strategies adopted by different sublineages and the influence of human migration.

Multidrug-resistance is a substantial threat to global elimination of tuberculosis. Understanding transmission patterns is crucial for control of the disease. We used a genomic and epidemiological approach to assess recent transmission of multidrug-resistant (MDR) tuberculosis and identify potential risk factors for transmission. We did a population-based, retrospective study of patients who tested positive for tuberculosis between Jan 1, 2009, and Dec 31, 2012, in Shanghai, China. We did variable-number-of-tandem-repeat genotyping and whole-genome sequencing of isolates. We measured strain diversity within and between genomically clustered isolates. Genomic and epidemiological data were combined to construct transmission networks. 367 (5%) of 7982 patients with tuberculosis had MDR tuberculosis and 324 (88%) of these had isolates available for genomic analysis. 103 (32%) of the 324 MDR strains were in 38 genomic clusters that differed by 12 or fewer single nucleotide polymorphisms (SNPs), indicating recent transmission of MDR strains. Patients who had delayed diagnosis or were older than 45 years had high risk of recent transmission. 235 (73%) patients with MDR tuberculosis probably had transmission of MDR strains. Transmission network analysis showed that 33 (87%) of the 38 clusters accumulated additional drug-resistance mutations through emergence or fixation of mutations during transmission. 68 (66%) of 103 clustered MDR strains had compensatory mutations of rifampicin resistance. Recent transmission of MDR tuberculosis strains, with increasing drug-resistance, drives the MDR tuberculosis epidemic in Shanghai, China. Whole-genome sequencing can measure of the heterogeneity of drug-resistant mutations within and between hosts and help to determine the transmission patterns of MDR tuberculosis. National Science and Technology Major Project, National Natural Science Foundation of China, and US National Insitutes of Health.

Significance: Mycobacterium tuberculosis Beijing family is a group of globally emerging bacterial strains that are responsible for more than a quarter of the global tuberculosis epidemic. Here, we combine whole-genome sequencing and large-scale genotyping to map the temporal and spatial changes of the genetic diversity within this strain family. We reveal a southern East Asia origin and a parallel evolution of this bacterial genotype with modern humans in East Asia during the last 30,000 years. The recently globally emerged Beijing strains mainly belong to a hypervirulent subtype that most likely has initially been selected for adaption to increased population densities during the agricultural transition in northern China.

Combination therapy is one of the most effective tools for limiting the emergence of drug resistance in pathogens. Despite the widespread adoption of combination therapy across diseases, drug resistance rates continue to rise, leading to failing treatment regimens. The mechanisms underlying treatment failure are well studied, but the processes governing successful combination therapy are poorly understood. We address this question by studying the population dynamics of Mycobacterium tuberculosis within tuberculosis patients undergoing treatment with different combinations of antibiotics. By combining very deep whole genome sequencing (~1000-fold genome-wide coverage) with sequential sputum sampling, we were able to detect transient genetic diversity driven by the apparently continuous turnover of minor alleles, which could serve as the source of drug-resistant bacteria. However, we report that treatment efficacy has a clear impact on the population dynamics: sufficient drug pressure bears a clear signature of purifying selection leading to apparent genetic stability. In contrast, M. tuberculosis populations subject to less drug pressure show markedly different dynamics, including cases of acquisition of additional drug resistance. Our findings show that for a pathogen like M. tuberculosis, which is well adapted to the human host, purifying selection constrains the evolutionary trajectory to resistance in effectively treated individuals. Nonetheless, we also report a continuous turnover of minor variants, which could give rise to the emergence of drug resistance in cases of drug pressure weakening. Monitoring bacterial population dynamics could therefore provide an informative metric for assessing the efficacy of novel drug combinations.