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.

Background

The HLA (human leukocyte antigen) molecules that present pathogen-derived epitopes to T cells are highly diverse. Correspondingly, many pathogens such as HIV evolve epitope variants in order to evade immune recognition. In contrast, another persistent human pathogen, Mycobacterium tuberculosis, has highly conserved epitope sequences. This raises the question whether there is also a difference in the ability of these pathogens' epitopes to bind diverse HLA alleles, referred to as an epitope's binding promiscuity. To address this question, we compared the in silico HLA binding promiscuity of T cell epitopes from pathogens with distinct infection strategies and outcomes of human exposure.

Methods

We used computer algorithms to predict the binding affinity of experimentally-verified microbial epitope peptides to diverse HLA-DR, HLA-A and HLA-B alleles. We then analyzed binding promiscuity of epitopes derived from HIV and M. tuberculosis. We also analyzed promiscuity of epitopes from Streptococcus pyogenes, which is known to exhibit epitope diversity, and epitopes of Bacillus anthracis and Clostridium tetani toxins, as these bacteria do not depend on human hosts for their survival or replication, and their toxin antigens are highly immunogenic human vaccines.

Results

We found that B. anthracis and C. tetani epitopes were the most promiscuous of the group that we analyzed. However, there was no consistent difference or trend in promiscuity in epitopes contained in HIV, M. tuberculosis, and S. pyogenes.

Conclusions

Our results show that human pathogens with distinct immune evasion strategies and epitope diversities exhibit equivalent levels of T cell epitope promiscuity. These results indicate that differences in epitope promiscuity do not account for the observed differences in epitope variation and conservation.

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.

One hundred healthy infants enrolled as controls in a tuberculosis vaccine study in Nyanza Province, Kenya provided anonymized samples for DNA sequence-based typing at the HLA-A, -B, -C, -DPB1, -DQA1, -DQB1, -DRB1, and -DRB3/4/5 loci. The purpose of the study was to characterize allele frequencies in the local population, to support studies of T cell immunity against pathogens, including Mycobacterium tuberculosis. There are no detectable deviations from Hardy Weinberg proportions for the HLA-B, -C, -DRB1, -DPB1, -DQA1 and -DQB1 loci. A minor deviation was detected at the HLA-A locus due to an excess of HLA-A*02:02, 29:02, 30:02, and 68:02 homozygotes. The genotype data are available in the Allele Frequencies Net Database under identifier 3393.

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.

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.