Sebastien Gagneux

Global phylogeography of Mycobacterium tuberculosis and implications for tuberculosis product development

New tools for controlling tuberculosis are urgently needed. Despite our emerging understanding of the biogeography of Mycobacterium tuberculosis, the implications for development of new diagnostics, drugs, and vaccines is unknown. M tuberculosis has a clonal genetic population structure that is geographically constrained. Evidence suggests strain-specific differences in virulence and immunogenicity in light of this global phylogeography. We propose a strain selection framework, based on robust phylogenetic markers, which will allow for systematic and comprehensive evaluation of new tools for tuberculosis control.

Human T cell epitopes of Mycobacterium tuberculosis are evolutionarily hyperconserved

Mycobacterium tuberculosis is an obligate human pathogen capable of persisting in individual hosts for decades. We sequenced the genomes of 21 strains representative of the global diversity and six major lineages of the M. tuberculosis complex (MTBC) at 40- to 90-fold coverage using Illumina next-generation DNA sequencing. We constructed a genome-wide phylogeny based on these genome sequences. Comparative analyses of the sequences showed, as expected, that essential genes in MTBC were more evolutionarily conserved than nonessential genes. Notably, however, most of the 491 experimentally confirmed human T cell epitopes showed little sequence variation and had a lower ratio of nonsynonymous to synonymous changes than seen in essential and nonessential genes. We confirmed these findings in an additional data set consisting of 16 antigens in 99 MTBC strains. These findings are consistent with strong purifying selection acting on these epitopes, implying that MTBC might benefit from recognition by human T cells.

Out-of-Africa migration and Neolithic coexpansion of Mycobacterium tuberculosis with modern humans

Tuberculosis caused 20% of all human deaths in the Western world between the seventeenth and nineteenth centuries and remains a cause of high mortality in developing countries. In analogy to other crowd diseases, the origin of human tuberculosis has been associated with the Neolithic Demographic Transition, but recent studies point to a much earlier origin. We analyzed the whole genomes of 259 M. tuberculosis complex (MTBC) strains and used this data set to characterize global diversity and to reconstruct the evolutionary history of this pathogen. Coalescent analyses indicate that MTBC emerged about 70,000 years ago, accompanied migrations of anatomically modern humans out of Africa and expanded as a consequence of increases in human population density during the Neolithic period. This long coevolutionary history is consistent with MTBC displaying characteristics indicative of adaptation to both low and high host densities.

The Influence of Host and Bacterial Genotype on the Development of Disseminated Disease with Mycobacterium tuberculosis

The factors that govern the development of tuberculosis disease are incompletely understood. We hypothesized that some strains of Mycobacterium tuberculosis (M. tuberculosis) are more capable of causing disseminated disease than others and may be associated with polymorphisms in host genes responsible for the innate immune response to infection. We compared the host and bacterial genotype in 187 Vietnamese adults with tuberculous meningitis (TBM) and 237 Vietnamese adults with uncomplicated pulmonary tuberculosis. The host genotype of tuberculosis cases was also compared with the genotype of 392 cord blood controls from the same population. Isolates of M. tuberculosis were genotyped by large sequence polymorphisms. The hosts were defined by polymorphisms in genes encoding Toll-interleukin 1 receptor domain containing adaptor protein (TIRAP) and Toll-like receptor-2 (TLR-2). We found a significant protective association between the Euro-American lineage of M. tuberculosis and pulmonary rather than meningeal tuberculosis (Odds ratio (OR) for causing TBM 0.395, 95% confidence intervals (C.I.) 0.193–0.806, P = 0.009), suggesting these strains are less capable of extra-pulmonary dissemination than others in the study population. We also found that individuals with the C allele of TLR-2 T597C allele were more likely to have tuberculosis caused by the East-Asian/Beijing genotype (OR = 1.57 [95% C.I. 1.15–2.15]) than other individuals. The study provides evidence that M. tuberculosis genotype influences clinical disease phenotype and demonstrates, for the first time, a significant interaction between host and bacterial genotypes and the development of tuberculosis.

Whole-genome sequencing of rifampicin-resistant Mycobacterium tuberculosis strains identifies compensatory mutations in RNA polymerase genes

Epidemics of drug-resistant bacteria emerge worldwide, even as resistant strains frequently have reduced fitness compared to their drug-susceptible counterparts. Data from model systems suggest that the fitness cost of antimicrobial resistance can be reduced by compensatory mutations; however, there is limited evidence that compensatory evolution has any significant role in the success of drug-resistant bacteria in human populations. Here we describe a set of compensatory mutations in the RNA polymerase genes of rifampicin-resistant M. tuberculosis, the etiologic agent of human tuberculosis (TB). M. tuberculosis strains harboring these compensatory mutations showed a high competitive fitness in vitro. Moreover, these mutations were associated with high fitness in vivo, as determined by examining their relative clinical frequency across patient populations. Of note, in countries with the worlds highest incidence of multidrug-resistant (MDR) TB, more than 30% of MDR clinical isolates had this form of mutation. Our findings support a role for compensatory evolution in the global epidemics of MDR TB.

Genomic Deletions Classify the Beijing/W Strains as a Distinct Genetic Lineage of Mycobacterium tuberculosis

ABSTRACT Beijing/W strains of Mycobacterium tuberculosis are geographically widespread and hypervirulent. To enhance our understanding of their origin and evolution, we sought phylogenetically informative large sequence polymorphisms (LSPs) within the Beijing/W family. Comparative whole-genome hybridization of Beijing/W strains revealed 21 LSPs, 7 of which were previously unreported. We show that some of these LSPs are unique event polymorphisms that can be used to define and subdivide the Beijing/W family. One LSP (RD105) was seen in all Beijing/W strains and thus serves as a useful marker for the identification of this family of strains. Additional LSPs (RD142, RD150, and RD181) further divided this family into four monophyletic subgroups, demonstrating a deeper population structure than previously appreciated. All Beijing/W strains were also observed to have an intact pks15/1 gene that is involved in the biosynthesis of a phenolic glycolipid, a putative virulence factor. A simple PCR assay using these Beijing/W strain-defining deletions will facilitate molecular epidemiological studies and may assist in the identification of the molecular basis of phenotypes associated with this important lineage of M. tuberculosis.

Progression to Active Tuberculosis, but Not Transmission, Varies by Mycobacterium tuberculosis Lineage in The Gambia

BACKGROUND There is considerable variability in the outcome of Mycobacterium tuberculosis infection. We hypothesized that Mycobacterium africanum was less likely than M. tuberculosis to transmit and progress to tuberculosis disease. METHODS In a cohort study of patients with tuberculosis and their household contacts in The Gambia, we categorized 1808 HIV-negative tuberculosis contacts according to exposure to M. tuberculosis or M. africanum. Positive skin test results indicated transmission, and development of tuberculosis during 2 years of follow-up indicated progression to disease. RESULTS Transmission rates were similar, but rates of progression to disease were significantly lower in contacts exposed to M. africanum than in those exposed to M. tuberculosis (1.0% vs. 2.9%; hazard ratio [HR], 3.1 [95% confidence interval {CI}, 1.1-8.7]). Within M. tuberculosis sensu stricto, contacts exposed to a Beijing family strain were most likely to progress to disease (5.6%; HR relative to M. africanum, 6.7 [95% CI, 2.0-22]). CONCLUSIONS M. africanum and M. tuberculosis transmit equally well to household contacts, but contacts exposed to M. africanum are less likely to progress to tuberculosis disease than those exposed to M. tuberculosis. The variable rate of progression by lineage suggests that tuberculosis variability matters in clinical settings and should be accounted for in studies evaluating tuberculosis vaccines and treatment regimens for latent tuberculosis infection.

Mycobacterium tuberculosis mutation rate estimates from different lineages predict substantial differences in the emergence of drug-resistant tuberculosis

A key question in tuberculosis control is why some strains of M. tuberculosis are preferentially associated with resistance to multiple drugs. We demonstrate that M. tuberculosis strains from lineage 2 (East Asian lineage and Beijing sublineage) acquire drug resistances in vitro more rapidly than M. tuberculosis strains from lineage 4 (Euro-American lineage) and that this higher rate can be attributed to a higher mutation rate. Moreover, the in vitro mutation rate correlates well with the bacterial mutation rate in humans as determined by whole-genome sequencing of clinical isolates. Finally, using a stochastic mathematical model, we demonstrate that the observed differences in mutation rate predict a substantially higher probability that patients infected with a drug-susceptible lineage 2 strain will harbor multidrug-resistant bacteria at the time of diagnosis. These data suggest that interventions to prevent the emergence of drug-resistant tuberculosis should target bacterial as well as treatment-related risk factors.

Genomic analysis identifies targets of convergent positive selection in drug-resistant Mycobacterium tuberculosis.

M. tuberculosis is evolving antibiotic resistance, threatening attempts at tuberculosis epidemic control. Mechanisms of resistance, including genetic changes favored by selection in resistant isolates, are incompletely understood. Using 116 newly sequenced and 7 previously sequenced M. tuberculosis whole genomes, we identified genome-wide signatures of positive selection specific to the 47 drug-resistant strains. By searching for convergent evolution—the independent fixation of mutations in the same nucleotide position or gene—we recovered 100% of a set of known resistance markers. We also found evidence of positive selection in an additional 39 genomic regions in resistant isolates. These regions encode components in cell wall biosynthesis, transcriptional regulation and DNA repair pathways. Mutations in these regions could directly confer resistance or compensate for fitness costs associated with resistance. Functional genetic analysis of mutations in one gene, ponA1, demonstrated an in vitro growth advantage in the presence of the drug rifampicin.

The W-beijing lineage of Mycobacterium tuberculosis overproduces triglycerides and has the DosR dormancy regulon constitutively upregulated

The Beijing family of Mycobacterium tuberculosis strains has been associated with epidemic spread and an increased likelihood of developing drug resistance. The characteristics that predispose this family to such clinical outcomes have not been identified, although one potential candidate, the phenolic glycolipid PGL-tb, has been shown to mediate a fulminant lethal disease in mice and rabbits due to lipid-mediated immunosuppression. However, PGL-tb is not uniformly expressed throughout the Beijing lineage and may not be the only unique virulence trait associated with this family. In an attempt to define phenotypes common to all Beijing strains, we interrogated a carefully selected set of isolates representing the five extant lineages of the Beijing family. Comparison of lipid production in this set revealed that all Beijing strains accumulated large quantities of triacylglycerides in in vitro aerobic culture. This accumulation was found to be coincident with upregulation of Rv3130c, whose product was previously characterized as a triacylglyceride synthase. Rv3130c is a member of the DosR-controlled regulon of M. tuberculosis, and further examination revealed that several members of this regulon were upregulated throughout this strain family. The upregulation of the DosR regulon may confer an adaptive advantage for growth in microaerophilic or anaerobic environments encountered by the bacillus during infection and thus may be related to the epidemiological phenomena associated with this important strain lineage.