a total of 22 mycobacterial strains were included in this work . to study the sequence diversity of t cell antigens in mtbc , we used illumina next - generation dna sequencing to generate nearly complete genome sequences from 20 strains representative of the six main human mtbc lineages , and one strain of mycobacterium canettii which is the closest known outgroup of mtbc3,22 ( table 1 ) . in addition , we used the published genome sequence of the h37rv laboratory strain of m. tuberculosis as a common reference23 . for each of the 21 strains newly sequenced , a mean of 6.8 million sequence reads with a mean length of 51 base pairs were generated and mapped to the h37rv reference genome . on average , the reads covered 98.9% of the 4.4 mb reference genome ( table 1 ) . the regions not covered primarily included members of the highly gc - rich and repetitive pe / ppe gene families24 . a total of 32,745 snps were identified , corresponding to an average of 1 snp call for every 3 kb of sequence generated . we used a total of 9,037 unique snps ( i.e. snps that occurred in one or several strains ) to derive a genome - wide phylogeny of 22 strains ( fig . these lineages were completely congruent to the strain groupings previously defined based on genomic deletion analysis and multilocus sequencing3,7,10 . the perfect congruence between these different phylogenetic markers further corroborates the highly clonal population structure of mtbc and lack of ongoing horizontal gene transfer in this organism25 . because of the comprehensive nature of genome - scale data , a higher degree of phylogenetic resolution the brown and green lineages ( also known as mycobacterium africanum ) are the most basal groups when compared to the m. canettii outgroup . however , the fact that the two m. africanum lineages represent the most ancestral forms of human mtbc reinforces the notion that human mtbc originated in africa3,7 . we used these genome sequence data and the phylogeny derived from them to compare the genetic diversity in antigens and other experimentally determined gene classes . for comparisons across different gene categories antigens were defined based on the presence of 491 experimentally confirmed human t cell epitopes ( supplementary table 3 ) , which were compiled through the immune epitope database ( iedb ) initiative26 . essential gene category was defined based on genome - wide analyses of transposon insertion mutants that were defective for the ability to grow on middlebrook 7h11 agar , or in the spleens of intravenously - infected mice , published previously27 - 28 . we excluded from this analysis genes belonging to the pe / ppe gene family24 and those related to mobile elements as they are difficult to study using current next - generation dna sequencing technologies ( total genes excluded : 273/3,990 ( 6.8% ) genes annotated in the h37rv reference genome ; supplementary table 4 ) . based on evolutionary theory and findings in other bacteria29 , one would expect that in contrast to non - essential genes , the essential genes in mtbc will be under stronger purifying selection and thus more evolutionary conserved . in support of this notion , we observed that on average essential genes harboured less nucleotide diversity than non - essential genes ( fig . we then compared the rates of synonymous and non - synonymous snps in the essential and non - essential gene categories . the synonymous and non - synonymous changes were derived by comparison to the most likely recent common ancestor of mtbc , which we inferred based on our new genome - wide phylogeny ( fig . 1 ) . because mtbc harbours little sequence diversity , it was necessary to analyze the distribution of synonymous and non - synonymous snps based on gene concatenates rather than individual genes . the two measures of distribution we used were based on the number of non - redundant snps across all 21 mtbc strains ( dn / ds based on measure a in table 2 and fig . 3 ) , and on the individual pairwise comparisons between each strain and the inferred most likely recent common ancestor ( dn / ds based on measures b in table 2 ) . from these analyses , we found that the dn / ds measures for essential genes were significantly lower than for non - essential genes ( measure a in fig . 3 ; measure b in table 2 , mann - whitney u test p<0.0001 ) . taken together , these data show that in mtbc essential genes are more evolutionary conserved than non - essential genes . because mtbc interacts with humans through antigen - specific cd4 or cd8 t - cells , we would expect t cell antigens to be among the most diverse genes in the genome . particularly when invoking a co - evolutionary arms - race and associated immune evasion , we would anticipate these antigens to be under diversifying selection and to be more variable than other genes in order to escape t cell recognition . however , when we analyzed the nucleotide diversity in 78 experimentally confirmed human t cell antigens ( supplementary table 2 ) , we found that they were on average not more diverse than essential genes ( fig . moreover , we found that the dn / ds measures in these antigens also resembled those of essential genes ( measure a in fig . 3 ; measure b in table 2 , mann - whitney u test p=0.77 ) . thus , human t cell antigens in mtbc do not appear to be under diversifying selection . t cell antigens consist of epitope regions that interact with human t cells , and non - epitope regions which are not targets of t cell recognition . , we generated a separate concatenate of the epitope regions and another concatenate of all corresponding non - epitope regions . because little data is currently available in the iedb with respect to whether these 491 epitopes are recognized by cd4 or cd8 t cells if immune escape was driving antigen evolution to evade t cell recognition in mtbc , we would expect non - synonymous changes to accumulate in epitope regions , leading to a high dn / ds . contrary to this expectation however , the overall dn / ds of the epitope regions was 0.53 , which was still similar to essential genes and lower than non - essential genes ( table 2, fig . 3 ) . moreover , when we analyzed the distribution of amino acid replacements in individual epitopes we found that the large majority ( 95% ) of the 491 epitopes showed no amino acid change ( fig . only five epitopes , contained in esxh , psts1 , and rv1986 , harboured more than one variable position ( supplementary table 5 ) . the higher number of amino acid substitutions in these five epitopes may reflect ongoing immune evasion , but further investigation is needed to determine whether the observed changes are due to immune pressure , other selection pressure(s ) , or mere random genetic drift3 . because these five epitopes were clear outliers compared to the large majority of t cell epitopes analyzed here , we repeated our dn / ds analysis after excluding the three antigens harbouring the five outlier epitopes . our analysis revealed that the epitope regions had the lowest dn / ds of all gene categories ( table 2 , fig . furthermore , when we compared the proportion of non - redundant non - synonymous changes in epitope and non - epitope regions , we found that epitopes were less likely than non - epitopes to harbour changes at non - synonymous sites ( measure a in table 2 ,  , p<0.05 ) , whereas no difference was observed at synonymous sites ( table 2 ,  , p=0.89 ) . to further corroborate our finding of hyperconservation of human t cell epitopes in mtbc , we repeated our analysis using a data set from a previous study in which 89 individual genes were sequenced in 99 human - adapted strains representative of the six major global lineages of mtbc3 . sixteen of these 89 genes belonged to the t cell antigens analyzed here , including two of the three outlier antigens esxh and psts1 . analysis of this additional dataset of 16 antigens in 99 mtbc strains revealed an overall dn / ds for the epitope regions of 0.74 . however , after excluding the two outlier antigens , the dn / ds dropped to 0.46 , which was again lower than the genome - based dn / ds values for essential and non - essential genes ( fig . , our findings strongly suggest that a large proportion of the mtbc genome known to interact with human t cells is highly conserved and under as strong , or perhaps even stronger , purifying selection than essential genes . a total of 22 mycobacterial strains were included in this work . to study the sequence diversity of t cell antigens in mtbc , we used illumina next - generation dna sequencing to generate nearly complete genome sequences from 20 strains representative of the six main human mtbc lineages , and one strain of mycobacterium canettii which is the closest known outgroup of mtbc3,22 ( table 1 ) . in addition , we used the published genome sequence of the h37rv laboratory strain of m. tuberculosis as a common reference23 . for each of the 21 strains newly sequenced , a mean of 6.8 million sequence reads with a mean length of 51 base pairs were generated and mapped to the h37rv reference genome . on average , the reads covered 98.9% of the 4.4 mb reference genome ( table 1 ) . the regions not covered primarily included members of the highly gc - rich and repetitive pe / ppe gene families24 . a total of 32,745 snps were identified , corresponding to an average of 1 snp call for every 3 kb of sequence generated . we used a total of 9,037 unique snps ( i.e. snps that occurred in one or several strains ) to derive a genome - wide phylogeny of 22 strains ( fig . these lineages were completely congruent to the strain groupings previously defined based on genomic deletion analysis and multilocus sequencing3,7,10 . the perfect congruence between these different phylogenetic markers further corroborates the highly clonal population structure of mtbc and lack of ongoing horizontal gene transfer in this organism25 . because of the comprehensive nature of genome - scale data , a higher degree of phylogenetic resolution the brown and green lineages ( also known as mycobacterium africanum ) are the most basal groups when compared to the m. canettii outgroup . however , the fact that the two m. africanum lineages represent the most ancestral forms of human mtbc reinforces the notion that human mtbc originated in africa3,7 . we used these genome sequence data and the phylogeny derived from them to compare the genetic diversity in antigens and other experimentally determined gene classes . for comparisons across different gene categories essential genes ,  non - essential genes , and  antigens ( supplementary fig . antigens were defined based on the presence of 491 experimentally confirmed human t cell epitopes ( supplementary table 3 ) , which were compiled through the immune epitope database ( iedb ) initiative26 . the  essential gene category was defined based on genome - wide analyses of transposon insertion mutants that were defective for the ability to grow on middlebrook 7h11 agar , or in the spleens of intravenously - infected mice , published previously27 - 28 . we excluded from this analysis genes belonging to the pe / ppe gene family24 and those related to mobile elements as they are difficult to study using current next - generation dna sequencing technologies ( total genes excluded : 273/3,990 ( 6.8% ) genes annotated in the h37rv reference genome ; supplementary table 4 ) . based on evolutionary theory and findings in other bacteria29 , one would expect that in contrast to non - essential genes , the essential genes in mtbc will be under stronger purifying selection and thus more evolutionary conserved . in support of this notion , we observed that on average essential genes harboured less nucleotide diversity than non - essential genes ( fig . we then compared the rates of synonymous and non - synonymous snps in the essential and non - essential gene categories . the synonymous and non - synonymous changes were derived by comparison to the most likely recent common ancestor of mtbc , which we inferred based on our new genome - wide phylogeny ( fig . 1 ) . because mtbc harbours little sequence diversity , it was necessary to analyze the distribution of synonymous and non - synonymous snps based on gene concatenates rather than individual genes . the two measures of distribution we used were based on the number of non - redundant snps across all 21 mtbc strains ( dn / ds based on measure a in table 2 and fig . 3 ) , and on the individual pairwise comparisons between each strain and the inferred most likely recent common ancestor ( dn / ds based on measures b in table 2 ) . from these analyses , we found that the dn / ds measures for essential genes were significantly lower than for non - essential genes ( measure a in fig . 3 ; measure b in table 2 , mann - whitney u test p<0.0001 ) . taken together , these data show that in mtbc essential genes are more evolutionary conserved than non - essential genes . because mtbc interacts with humans through antigen - specific cd4 or cd8 t - cells , we would expect t cell antigens to be among the most diverse genes in the genome . particularly when invoking a co - evolutionary arms - race and associated immune evasion , we would anticipate these antigens to be under diversifying selection and to be more variable than other genes in order to escape t cell recognition . however , when we analyzed the nucleotide diversity in 78 experimentally confirmed human t cell antigens ( supplementary table 2 ) , we found that they were on average not more diverse than essential genes ( fig . moreover , we found that the dn / ds measures in these antigens also resembled those of essential genes ( measure a in fig . 3 ; measure b in table 2 , mann - whitney u test p=0.77 ) . thus , human t cell antigens in mtbc do not appear to be under diversifying selection . instead , purifying selection appears to be the driving selection pressure in these genes . t cell antigens consist of epitope regions that interact with human t cells , and non - epitope regions which are not targets of t cell recognition . hence , we decided to study these regions separately . to this end , we generated a separate concatenate of the epitope regions and another concatenate of all corresponding non - epitope regions . because little data is currently available in the iedb with respect to whether these 491 epitopes are recognized by cd4 or cd8 t cells , we analyzed them as one class . if immune escape was driving antigen evolution to evade t cell recognition in mtbc , we would expect non - synonymous changes to accumulate in epitope regions , leading to a high dn / ds . contrary to this expectation however , the overall dn / ds of the epitope regions was 0.53 , which was still similar to essential genes and lower than non - essential genes ( table 2, fig . 3 ) . moreover , when we analyzed the distribution of amino acid replacements in individual epitopes we found that the large majority ( 95% ) of the 491 epitopes showed no amino acid change ( fig . only five epitopes , contained in esxh , psts1 , and rv1986 , harboured more than one variable position ( supplementary table 5 ) . the higher number of amino acid substitutions in these five epitopes may reflect ongoing immune evasion , but further investigation is needed to determine whether the observed changes are due to immune pressure , other selection pressure(s ) , or mere random genetic drift3 . because these five epitopes were clear outliers compared to the large majority of t cell epitopes analyzed here , we repeated our dn / ds analysis after excluding the three antigens harbouring the five outlier epitopes . our analysis revealed that the epitope regions had the lowest dn / ds of all gene categories ( table 2 , fig . furthermore , when we compared the proportion of non - redundant non - synonymous changes in epitope and non - epitope regions , we found that epitopes were less likely than non - epitopes to harbour changes at non - synonymous sites ( measure a in table 2 ,  , p<0.05 ) , whereas no difference was observed at synonymous sites ( table 2 ,  , p=0.89 ) . to further corroborate our finding of hyperconservation of human t cell epitopes in mtbc , we repeated our analysis using a data set from a previous study in which 89 individual genes were sequenced in 99 human - adapted strains representative of the six major global lineages of mtbc3 . sixteen of these 89 genes belonged to the t cell antigens analyzed here , including two of the three outlier antigens esxh and psts1 . analysis of this additional dataset of 16 antigens in 99 mtbc strains revealed an overall dn / ds for the epitope regions of 0.74 . however , after excluding the two outlier antigens , the dn / ds dropped to 0.46 , which was again lower than the genome - based dn / ds values for essential and non - essential genes ( fig . , our findings strongly suggest that a large proportion of the mtbc genome known to interact with human t cells is highly conserved and under as strong , or perhaps even stronger , purifying selection than essential genes . in this study of 22 mtbc genomes , we demonstrate that , as expected , essential genes are more conserved than non - essential genes . surprisingly , however , we found that the large majority of the currently known t cell antigens are as conserved as essential genes . furthermore , the epitope regions of these antigen genes are the most highly conserved regions we studied . this observation , that the regions of the genome that interact with the human adaptive immune system appear to be under even stronger purifying selection than essential genes , is inconsistent with a classical model of an evolutionary arms - race . it is possible that the known human t cell epitopes that we found to be hyperconserved represent a select subset of all of the human t cell epitopes encoded in the genome , and that certain approaches to epitope identification have favoured discovery of hyperconserved epitopes in mtbc . for example , since most , if not all of the epitope discovery efforts to date have utilized proteins and/or peptide sequences of strains from one lineage ( lineage 4 ) and t cells from humans that are likely to have been infected by strains of other lineages , the assays used may have been especially suited to identification of hyperconserved and/or cross - reactive epitopes . while further investigation using alternative approaches to epitope discovery may reveal that variable epitopes that exhibit evidence of positive selection exist in the mtbc , it is likely that the large number of epitopes that we examined will remain a significant subset of the total , and that future vaccine development efforts will need to account for the possibility that immune recognition of certain epitopes may actually provide a net benefit to the bacteria . lack of antigenic variation and immune evasion has been reported for a number of other human pathogens , including rna viruses such as measles , mumps , rubella , and influenza type c31 . theoretical studies have suggested that the absence of immune escape variants in these viruses might be due to structural constraints in viral proteins or negative mutational effects leading to reduced infectivity or transmission31 . while we can not exclude the possibility that structural and functional constraints that are independent of t cell recognition contribute to hyperconservation of the regions encoding mtbc peptides recognized by human t cells , one important characteristic of the aformentioned viral pathogens is that they spread among young and immunologically naive hosts , which might eliminate the need for immune evasion31 . moreover , infection by these viruses usually results in acute disease , followed by elimination of the infection through adaptive immunity , and acquisition of lifelong immunity against re - infection . by contrast , mtbc causes chronic and often lifelong infections , and adaptive immunity is usually unable to completely clear the infection18 . furthermore , tuberculosis patients are prone to re - infection32 , and mixed infections are also increasingly recognized33 . these observations suggest that the biological basis for the lack of antigenic variation in mtbc reported here differs from what has been proposed for antigenically homogeneous rna viruses31 . in addition , we determined that the fraction of hyperconserved t cell epitopes of the mtbc that are derived from essential genes is indistinguishable from the frequency of essential genes in the mtbc genome as a whole ( 18% versus 21% , respectively ;  = 0.28 , p = 0.59 ) , indicating that our results were not skewed by over - representation of t cell epitopes in essential genes . moreover , the t cell epitopes that we analyzed are present in genes from diverse gene ontologies , and the representation of five main gene categories ( defined based on the ncbi categories of orthologous groups ( cog ) ) was no different in the t cell antigens when compared to the genome overall (  with 4 degrees of freedom = 5.8 , p = 0.21 ; supplementary table 6 ) . hence the only identifiable common property of these regions is their recognition by human t lymphocytes . these findings suggest that t lymphocyte recognition is an important factor in hyperconservation of these sequences , and that other structural or functional constraints are unlikely to fully account for the lack of sequence variation in these domains . our data suggest that t cell epitopes in mtbc are under strong selection pressure to be maintained , perhaps because the immune response they elicit in humans , which are essential for survival of an individual host , might partially work towards the pathogen s benefit . one potential mechanism of benefit to mtbc from human t cell recognition is that human t cell responses are essential for mtbc to establish latent infection . this notion is supported by the fact that cd4 t cell - deficient hiv - positive individuals progress rapidly to active disease after infection , rather than to sustain prolonged periods of latent tuberculosis34 . latent infection mediated by host t cell responses , with subsequent reactivation to active disease often occurring decades after initial infection , is a key characteristic of human tuberculosis , and might have evolved as a way for mtbc to transmit to later generations of susceptible hosts35 . in addition , there is evidence that t cell responses may contribute directly to human - to - human transmission of mtbc . in particular , cavitary tuberculosis , which generates secondary cases more efficiently than other disease forms36 , rarely occurs in cd4 t cell - deficient hiv - positive individuals , and the frequency of cavitary lung lesions in hiv - infected patients with tuberculosis is directly correlated with the number of peripheral cd4 t cells37 . while the mechanisms of lung cavitation in tuberculosis are poorly understood , these observations suggest that cd4 t cells directly or indirectly mediate tissue damage in tuberculosis , and together with our finding of epitope hyperconservation indicate that certain t cell responses may be detrimental to the host and beneficial to the pathogen . hence our findings suggest that mtbc takes advantage of host adaptive immunity to increase its likelihood of spread , and that the benefits of enhanced transmission exceed the costs of within - host cellular immune responses to these epitopes . in this manner , mtbc may resemble hiv , for which there is evidence that virulence has evolved , not to maximize replication of the virus within individual hosts , but to maximize the likelihood of its transmission38 . whether t cell responses to other epitopes , or whether specific t cell subsets ( e.g. th17 versus th1 ) that benefit the host and not the bacteria can be identified will require additional studies in humans . one limitation of this study was the exclusion of pe / ppe genes because of technical reasons . some of these genes are known to vary and to be cell - surface exposed , which has lead to the hypothesis they might be involved in antigenic variation24 . however , no direct evidence for this has yet been presented . future work will need to clarify the function and evolution of pe / ppe genes . by contrast , all the t cell antigens included in this study have been experimentally confirmed26 . furthermore , some of them are being targeted by new tuberculosis diagnostics and vaccines39 . our findings thus have important implications for the development of these new tools . on the one hand , the fact that mtbc harbours little sequence diversity in t cell antigens will facilitate the development of diagnostics that are universally applicable across geographical regions where mtbc strains differ8 . on the other hand , the possibility that the immune responses induced by vaccine antigens might partially benefit the pathogen suggests current efforts in vaccine research should be broadened . most disturbing is the suggestion that vaccine induced immunity against these conserved epitopes may perversely increase transmission . in this respect , it is interesting to note that the currently available tuberculosis vaccine bacille - calmette - guerin ( bcg ) , which is a live vaccine based on an attenuated from of m. bovis , offers no protection against pulmonary tuberculosis in adults40 . more importantly , some clinical trials of bcg have even reported an increased risk of tuberculosis in vaccinees compared to unvaccinated individuals41 . thus , in contrast to standard reverse vaccinology , in which the least variable antigens of a genome are targeted42 , research into new tuberculosis vaccines should explore more variable regions of the mtbc genome . while most of the t cell epitopes anlyzed here were highly conserved , five epitopes in three antigens harboured a larger number of amino acid changes . the fact that the dn / ds measure dropped sharply after excluding these outlier antigens from the analysis further supports the notion that they are indeed outliers compared to the other antigens . one of these outlier antigens , esxh ( rv0288 , also known as tb10.4 ) is a member of a gene family known to encode a type vii secretion system43 . thus even though most of the other vaccine antigens analyzed here are conserved , our finding that this particular vaccine antigen harbours a comparatively high number of amino acid substitutions across a panel of global mtbc isolates , suggests that strain diversity should be considered during further development of the new vaccine candidates containing esxh8 . we detected significant differences in dn / ds between essential , non - essential , and antigenic genes . however , the individual dn / ds values remain high when compared to most other bacteria44 . such a high dn / ds was reported previously for mtbc , and has been linked to reduced selective constraint against slightly deleterious mutations3 . it was proposed that the serial transmission bottlenecks associated with patient - to - patient transmission in mtbc could lead to an increase in random genetic drift compared to the forces of natural selection . our new data show that even though the strength of purifying selection in mtbc might be reduced overall compared to other bacteria , it clearly is still acting on and capable of differentiating between gene categories . in summary , we show that t cell epitopes of mtbc are highly conserved , and do not reflect any ongoing evolutionary arms - race or immune - evasion . instead , the patterns observed might be indicative of a distinct evolutionary strategy of immune - subversion developed by this highly successful pathogen . other intracellular bacteria such as salmonella enterica serovar typhi exhibit a similar lack of antigenic variation45 , suggesting