Mycobacterium bovis bacillus Calmette-Guerin (BCG), the only vaccine currently used against tuberculosis, is an attenuated derivative of M. bovis that has been propagated in vitro for more than 40 years. We have previously reported that the experimentally-verified human T cell epitopes of the M. tuberculosis complex (MTBC) are the most conserved elements of the genome; whether immune recognition is the force driving the conservation of epitopes in the MTBC is unknown. Therefore, we sequenced the genomes of 12 BCG strains to determine whether T cell epitopes were under selection pressure during BCG in vitro evolution. We constructed a genome-wide phylogeny and refined the previously-determined BCG phylogeny. Notably, we identified a new cluster between BCG Japan and BCG Russia, and repositioned the relationships of several strains within the lineage. We also compared the sequence diversity of 1530 experimentally verified human T cell epitopes in the BCG vaccines with those in the MTBC. We found 23% of the known T cell epitopes are absent, and that the majority (82%) of the absent epitopes in BCG are contained in 6 proteins encoded in 2 regions of difference (RD) unique to BCG strains. We also found that T cell epitope sequences in BCG are more conserved than non-epitope sequences in the same gene. Finally, we find evidence that epitope sequence variation in BCG potentially affects human T cell recognition. These findings provide new insight into sequence variation in a slow-growing bacterium closely related to the MTBC that has been subjected to prolonged passage outside of a mammalian host, and indicate little difference in the extent of variation in vivo and in vitro.

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.

Unlabelled

The Mycobacterium tuberculosis genome includes the large family of pe_pgrs genes, whose functions are unknown. Because of precedents in other pathogens in which gene families showing high sequence variation are involved in antigenic variation, a similar role has been proposed for the pe_pgrs genes. However, the impact of immune selection on pe_pgrs genes has not been examined. Here, we sequenced 27 pe_pgrs genes in 94 clinical strains from five phylogenetic lineages of the M. tuberculosis complex (MTBC). We found that pe_pgrs genes were overall more diverse than the remainder of the MTBC genome, but individual members of the family varied widely in their nucleotide diversity and insertion/deletion (indel) content: some were more, and others were much less, diverse than the genome average. Individual pe_pgrs genes also differed in the ratio of nonsynonymous to synonymous mutations, suggesting that different selection pressures act on individual pe_pgrs genes. Using bioinformatic methods, we tested whether sequence diversity in pe_pgrs genes might be selected by human T cell recognition, the major mechanism of adaptive immunity to MTBC. We found that the large majority of predicted human T cell epitopes were confined to the conserved PE domain and experimentally confirmed the antigenicity of this domain in tuberculosis patients. In contrast, despite being genetically diverse, the PGRS domains harbored few predicted T cell epitopes. These results indicate that human T cell recognition is not a significant force driving sequence diversity in pe_pgrs genes, which is consistent with the previously reported conservation of human T cell epitopes in the MTBC.

Importance

Recognition of Mycobacterium tuberculosis antigens by T lymphocytes is known to be important for immune protection against tuberculosis, but it is unclear whether human T cell recognition drives antigenic variation in M. tuberculosis. We previously discovered that the known human T cell epitopes in the M. tuberculosis complex are highly conserved, but we hypothesized that undiscovered epitopes with naturally occurring sequence variants might exist. To test this hypothesis, we examined the pe_pgrs genes, a large family of genes that has been proposed to function in immune evasion by M. tuberculosis. We found that the pe_pgrs genes exhibit considerable sequence variation, but the regions containing T cell epitopes and the regions of variation are distinct. These findings confirm that the majority of human T cell epitopes of M. tuberculosis are highly conserved and indicate that selection forces other than T cell recognition drive sequence variation in the pe_pgrs genes.

Pathogens that evade adaptive immunity typically exhibit antigenic variation. By contrast, it appears that although the chronic human tuberculosis (TB)-causing pathogen Mycobacterium tuberculosis needs to counter host T cell responses, its T cell epitopes are hyperconserved. Here we present an extensive analysis of the T cell epitopes of M. tuberculosis. We combined population genomics with experimental immunology to determine the number and identity of T cell epitope sequence variants in 216 phylogenetically diverse strains of M. tuberculosis. Antigen conservation is indeed a hallmark of M. tuberculosis. However, our analysis revealed a set of seven variable antigens that were immunogenic in subjects with active TB. These findings suggest that M. tuberculosis uses mechanisms other than antigenic variation to evade T cells. T cell epitopes that exhibit sequence variation may not be subject to the same evasion mechanisms, and hence vaccines that include such variable epitopes may be more efficacious.

Molecular epidemiological assessments, drug treatment optimization, and development of immunological interventions all depend on understanding pathogen adaptation and genetic variation, which differ for specific pathogens. Mycobacterium tuberculosis is an exceptionally successful human pathogen, yet beyond knowledge that this bacterium has low overall genomic variation but acquires drug resistance mutations, little is known of the factors that drive its population genomic characteristics. Here, we compared the genetic diversity of the bacteria that established infection to the bacterial populations obtained from infected tissues during murine M. tuberculosis pulmonary infection and human disseminated M. bovis BCG infection. We found that new mutations accumulate during in vitro culture, but that in vivo, purifying selection against new mutations dominates, indicating that M. tuberculosis follows a dominant lineage model of evolution. Comparing bacterial populations passaged in T cell-deficient and immunocompetent mice, we found that the presence of T cells is associated with an increase in the diversity of the M. tuberculosis genome. Together, our findings put M. tuberculosis genetic evolution in a new perspective and clarify the impact of T cells on sequence diversity of M. tuberculosis.

Recent studies suggest that M. tuberculosis lineage and host genetics interact to impact how active tuberculosis presents clinically. We determined the phylogenetic lineages of M. tuberculosis isolates from participants enrolled in the Tuberculosis Trials Consortium Study 28, conducted in Brazil, Canada, South Africa, Spain, Uganda and the United States, and secondarily explored the relationship between lineage, clinical presentation and response to treatment. Large sequence polymorphisms and single nucleotide polymorphisms were analyzed to determine lineage and sublineage of isolates. Of 306 isolates genotyped, 246 (80.4%) belonged to the Euro-American lineage, with sublineage 724 predominating at African sites (99/192, 51.5%), and the Euro-American strains other than 724 predominating at non-African sites (89/114, 78.1%). Uneven distribution of lineages across regions limited our ability to discern significant associations, nonetheless, in univariate analyses, Euro-American sublineage 724 was associated with more severe disease at baseline, and along with the East Asian lineage was associated with lower bacteriologic conversion after 8 weeks of treatment. Disease presentation and response to drug treatment varied by lineage, but these associations were no longer statistically significant after adjustment for other variables associated with week-8 culture status.

Although small molecules shed from pathogens are widely used to diagnose infection, such tests have not been widely implemented for tuberculosis. Here we show that the recently identified compound, 1-tuberculosinyladenosine (1-TbAd), accumulates to comprise >1% of all Mycobacterium tuberculosis lipid. In vitro and in vivo, two isomers of TbAd were detected that might serve as infection markers. Using mass spectrometry and nuclear magnetic resonance, we established the structure of the previously unknown molecule, N(6)-tuberculosinyladenosine (N(6)-TbAd). Its biosynthesis involves enzymatic production of 1-TbAd by Rv3378c followed by conversion to N(6)-TbAd via the Dimroth rearrangement. Intact biosynthetic genes are observed only within M. tuberculosis complex bacteria, and TbAd was not detected among other medically important pathogens, environmental bacteria, and vaccine strains. With no substantially similar known molecules in nature, the discovery and in vivo detection of two abundant terpene nucleosides support their development as specific diagnostic markers of tuberculosis.

The transcontinental spread of multidrug-resistant (MDR) tuberculosis is poorly characterized in molecular epidemiologic studies. We used genomic sequencing to understand the establishment and dispersion of MDR Mycobacterium tuberculosis within a group of immigrants to the United States. We used a genomic epidemiology approach to study a genotypically matched (by spoligotype, IS6110 restriction fragment length polymorphism, and mycobacterial interspersed repetitive units-variable number of tandem repeat signature) lineage 2/Beijing MDR strain implicated in an outbreak of tuberculosis among refugees in Thailand and consecutive cases within California. All 46 MDR M. tuberculosis genomes from both Thailand and California were highly related, with a median difference of 10 single-nucleotide polymorphisms (SNPs). The Wat Tham Krabok (WTK) strain is a new sequence type distinguished from all known Beijing strains by 55 SNPs and a genomic deletion (Rv1267c) associated with increased fitness. Sequence data revealed a highly prevalent MDR strain that included several closely related but distinct allelic variants within Thailand, rather than the occurrence of a single outbreak. In California, sequencing data supported multiple independent introductions of WTK with subsequent transmission and reactivation within the state, as well as a potential super spreader with a prolonged infectious period. Twenty-seven drug resistance-conferring mutations and 4 putative compensatory mutations were found within WTK strains. Genomic sequencing has substantial epidemiologic value in both low- and high-burden settings in understanding transmission chains of highly prevalent MDR strains.

The Beijing family of Mycobacterium tuberculosis strains is part of lineage 2 (also known as the East Asian lineage). In clinical studies, we have observed that isolates from the sublineage RD207 of lineage 2 were more readily transmitted among humans. To investigate the basis for this difference, we tested representative strains with the characteristic Beijing spoligotype from four of the five sublineages of lineage 2 in the guinea pig model and subjected these strains to comparative whole-genome sequencing. The results of these studies showed that all of the clinical strains were capable of growing and causing lung pathology in guinea pigs after low-dose aerosol exposure. Differences between the abilities of the four sublineages to grow in the lungs of these animals were not overt, but members of RD207 were significantly more pathogenic, resulting in severe lung damage. The RD207 strains also induced much higher levels of markers associated with regulatory T cells and showed a significant loss of activated T cells in the lungs over the course of the infections. Whole-genome sequencing of the strains revealed mutations specific for RD207 which may explain this difference. Based on these data, we hypothesize that the sublineages of M. tuberculosis are associated with distinct pathological and clinical phenotypes and that these differences influence the transmissibility of particular M. tuberculosis strains in human populations.

Background

Bedaquiline is a life-saving tuberculosis drug undergoing global scale-up. People at risk of weak tuberculosis drug regimens are a priority for novel drug access despite the potential source of Mycobacterium tuberculosis-resistant strains. We aimed to characterise bedaquiline resistance in individuals who had sustained culture positivity during bedaquiline-based treatment.

Methods

We did a retrospective longitudinal cohort study of adults (aged ≥18 years) with culture-positive pulmonary tuberculosis who received at least 4 months of a bedaquiline-containing regimen from 12 drug-resistant tuberculosis treatment facilities in Cape Town, South Africa, between Jan 20, 2016, and Nov 20, 2017. Sputum was programmatically collected at baseline (ie, before bedaquiline initiation) and each month to monitor treatment response per the national algorithm. The last available isolate from the sputum collected at or after 4 months of bedaquiline was designated the follow-up isolate. Phenotypic drug susceptibility testing for bedaquiline was done on baseline and follow-up isolates in MGIT960 media (WHO-recommended critical concentration of 1 μg/mL). Targeted deep sequencing for Rv0678, atpE, and pepQ, as well as whole-genome sequencing were also done.

Findings

In total, 40 (31%) of 129 patients from an estimated pool were eligible for this study. Overall, three (8%) of 38 patients assessable by phenotypic drug susceptibility testing for bedaquiline had primary resistance, 18 (47%) gained resistance (acquired or reinfection), and 17 (45%) were susceptible at both baseline and follow-up. Several Rv0678 and pepQ single-nucleotide polymorphisms and indels were associated with resistance. Although variants occurred in Rv0676c and Rv1979c, these variants were not associated with resistance. Targeted deep sequencing detected low-level variants undetected by whole-genome sequencing; however, none were in genes without variants already detected by whole-genome sequencing. Patients with baseline fluoroquinolone resistance, clofazimine exposure, and four or less effective drugs were more likely to have bedaquiline-resistant gain. Resistance gain was primarily due to acquisition; however, some reinfection by resistant strains occurred.

Interpretation

Bedaquiline-resistance gain, for which we identified risk factors, was common in these programmatically treated patients with sustained culture positivity. Our study highlights risks associated with implementing life-saving new drugs and shows evidence of bedaquiline-resistance transmission. Routine drug susceptibility testing should urgently accompany scale-up of new drugs; however, rapid drug susceptibility testing for bedaquiline remains challenging given the diversity of variants observed.

Funding

Doris Duke Charitable Foundation, US National Institute of Allergy and Infectious Diseases, South African Medical Research Council, National Research Foundation, Research Foundation Flanders, Stellenbosch University Faculty of Medicine Health Sciences, South African National Research Foundation, Swiss National Science Foundation, and Wellcome Trust.