High resolution discrimination of clinical Mycobacterium tuberculosis complex strains based on single nucleotide polymorphisms

Recently, the diversity of the Mycobacterium tuberculosis complex (MTBC) population structure has been described in detail. Based on geographical separation and specific host pathogen co-evolution shaping MTBC virulence traits, at least 20 major lineages/genotypes have evolved finally leading to a clear influence of strain genetic background on transmissibility, clinical presentation/outcome, and resistance development. Therefore, high resolution genotyping for characterization of strains in larger studies is mandatory for understanding mechanisms of host-pathogen-interaction and to improve tuberculosis (TB) control. Single nucleotide polymorphisms (SNPs) represent the most reliable markers for lineage classification of clinical isolates due to the low levels of homoplasy, however their use is hampered either by low discriminatory power or by the need to analyze a large number of genes to achieve higher resolution. Therefore, we carried out de novo sequencing of 26 genes (approx. 20000 bp per strain) in a reference collection of MTBC strains including all major genotypes to define a highly discriminatory gene set. Overall, 161 polymorphisms were detected of which 59 are genotype-specific, while 13 define deeper branches such as the Euro-American lineage. Unbiased investigation of the most variable set of 11 genes in a population based strain collection (one year, city of Hamburg, Germany) confirmed the validity of SNP analysis as all strains were classified with high accuracy. Taken together, we defined a diagnostic algorithm which allows the identification of 17 MTBC phylogenetic lineages with high confidence for the first time by sequencing analysis of just five genes. In conclusion, the diagnostic algorithm developed in our study is likely to open the door for a low cost high resolution sequence/SNP based differentiation of the MTBC with a very high specificity. High throughput assays can be established which will be needed for large association studies that are mandatory for detailed investigation of host-pathogen-interaction during TB infection.     

Human macrophage responses to clinical isolates from the Mycobacterium tuberculosis complex discriminate between ancient and modern lineages

The aim of the present study was to determine whether there is a correlation between phylogenetic relationship and inflammatory response amongst a panel of clinical isolates representative of the global diversity of the human Mycobacterium tuberculosis Complex (MTBC). Measurement of cytokines from infected human peripheral blood monocyte-derived macrophages revealed a wide variation in the response to different strains. The same pattern of high or low response to individual strains was observed for different pro-inflammatory cytokines and chemokines, and was conserved across multiple human donors. Although each major phylogenetic lineage of MTBC included strains inducing a range of cytokine responses, we found that overall inflammatory phenotypes differed significantly across lineages. In particular, comparison of evolutionarily modern lineages demonstrated a significant skewing towards lower early inflammatory response. The differential response to ancient and modern lineages observed using GM-CSF derived macrophages was also observed in autologous monocyte-derived dendritic cells and murine bone marrow-derived macrophages, but not in human unfractionated peripheral blood mononuclear cells. We hypothesize that the reduced immune responses to modern lineages contribute to more rapid disease progression and transmission, which might be a selective advantage in the context of expanding human populations. In addition to the lineage effects, the large strain-to-strain variation in innate immune responses elicited by MTBC will need to be considered in tuberculosis vaccine development.     

Mycobacterium africanum Is Associated with Patient Ethnicity in Ghana

Mycobacterium africanum is a member of the Mycobacterium tuberculosis complex (MTBC) and an important cause of human tuberculosis in West Africa that is rarely observed elsewhere. Here we genotyped 613 MTBC clinical isolates from Ghana, and searched for associations between the different phylogenetic lineages of MTBC and patient variables. We found that 17.1% (105/613) of the MTBC isolates belonged to M. africanum, with the remaining belonging to M. tuberculosis sensu stricto. No M. bovis was identified in this sample. M. africanum was significantly more common in tuberculosis patients belonging to the Ewe ethnic group (adjusted odds ratio: 3.02; 95% confidence interval: 1.67–5.47, p<0.001). Stratifying our analysis by the two phylogenetic lineages of M. africanum (i.e. MTBC Lineages 5 and 6) revealed that this association was mainly driven by Lineage 5 (also known as M. africanum West Africa 1). Our findings suggest interactions between the genetic diversity of MTBC and human diversity, and offer a possible explanation for the geographical restriction of M. africanum to parts of West Africa.     

Genetic Diversity of Mycobacterium tuberculosis Isolates from Assam,India: Dominance of Beijing Family and Discovery of Two New Clades Related to CAS1_Delhi and EAI Family Based on Spoligotyping and MIRU-VNTR Typing

Tuberculosis (TB) is one of the major public health concerns in Assam, a remote state located in the northeastern (NE) region of India. The present study was undertaken to explore the circulating genotypes of Mycobacterium tuberculosis complex (MTBC) in this region. A total of 189 MTBC strains were collected from smear positive pulmonary tuberculosis cases from different designated microscopy centres (DMC) from various localities of Assam. All MTBC isolates were cultured on Lowenstein-Jensen (LJ) media and subsequently genotyped using spoligotyping and 24-loci mycobacterial interspersed repetitive units-variable number of tandem repeats (MIRU-VNTR) typing. Spoligotyping of MTBC isolates revealed 89 distinct spoligo patterns. The most dominant MTBC strain belonged to Beijing lineage and was represented by 35.45% (n = 67) of total isolates, followed by MTBC strains belonging to Central Asian-Delhi (CAS/Delhi) lineage and East African Indian (EAI5) lineage. In addition, in the present study 43 unknown spoligo patterns were detected. The discriminatory power of spoligotyping was found to be 0.8637 based on Hunter Gaston Discriminatory Index (HGDI). On the other hand, 24-loci MIRU-VNTR typing revealed that out of total 189 MTBC isolates from Assam 185 (97.9%) isolates had unique MIRU-VNTR profiles and 4 isolates grouped into 2 clusters. Phylogenetic analysis of 67 Beijing isolates based on 24-loci MIRU-VNTR typing revealed that Beijing isolates from Assam represent two major groups, each comprising of several subgroups. Neighbour-Joining (NJ) phylogenetic tree analysis based on combined spoligotyping and 24-loci MIRU-VNTR data of 78 Non-Beijing isolates was carried out for strain lineage identification as implemented by MIRU-VNTRplus database. The important lineages of MTBC identified were CAS/CAS1_Delhi (41.02%, n = 78) and East-African-Indian (EAI, 33.33%). Interestingly, phylogenetic analysis of orphan (23.28%) MTBC spoligotypes revealed that majority of these orphan isolates from Assam represent two new sub-clades Assam/EAI and Assam/CAS. The prevalence of multidrug resistance (MDR) in Beijing and Non-Beijing strains was found to be 10.44% and 9.01% respectively. In conclusion, the present study has shown the predominance of Beijing isolates in Assam which is a matter of great concern because Beijing strains are considered to be ecologically more fit enabling wider dissemination of M. tuberculosis. Other interesting finding of the present study is the discovery of two new clades of MTBC isolates circulating in Assam. More elaborate longitudinal studies are required to be undertaken in this region to understand the transmission dynamics of MTBC.     

Genotyping of Genetically Monomorphic Bacteria: DNA Sequencing in Mycobacterium tuberculosis Highlights the Limitations of Current Methodologies

Because genetically monomorphic bacterial pathogens harbour little DNA sequence diversity, most current genotyping techniques used to study the epidemiology of these organisms are based on mobile or repetitive genetic elements. Molecular markers commonly used in these bacteria include Clustered Regulatory Short Palindromic Repeats (CRISPR) and Variable Number Tandem Repeats (VNTR). These methods are also increasingly being applied to phylogenetic and population genetic studies. Using the Mycobacterium tuberculosis complex (MTBC) as a model, we evaluated the phylogenetic accuracy of CRISPR- and VNTR-based genotyping, which in MTBC are known as spoligotyping and Mycobacterial Interspersed Repetitive Units (MIRU)-VNTR-typing, respectively. We used as a gold standard the complete DNA sequences of 89 coding genes from a global strain collection. Our results showed that phylogenetic trees derived from these multilocus sequence data were highly congruent and statistically robust, irrespective of the phylogenetic methods used. By contrast, corresponding phylogenies inferred from spoligotyping or 15-loci-MIRU-VNTR were incongruent with respect to the sequence-based trees. Although 24-loci-MIRU-VNTR performed better, it was still unable to detect all strain lineages. The DNA sequence data showed virtually no homoplasy, but the opposite was true for spoligotyping and MIRU-VNTR, which was consistent with high rates of convergent evolution and the low statistical support obtained for phylogenetic groupings defined by these markers. Our results also revealed that the discriminatory power of the standard 24 MIRU-VNTR loci varied by strain lineage. Taken together, our findings suggest strain lineages in MTBC should be defined based on phylogenetically robust markers such as single nucleotide polymorphisms or large sequence polymorphisms, and that for epidemiological purposes, MIRU-VNTR loci should be used in a lineage-dependent manner. Our findings have implications for strain typing in other genetically monomorphic bacteria.     

First Molecular Epidemiology Study of Mycobacterium tuberculosis in Kiribati

Tuberculosis incidence rates in Kiribati are among the highest in the Western Pacific Region, however the genetic diversity of circulating Mycobacterium tuberculosis complex strains (MTBC) and transmission dynamics are unknown. Here, we analysed MTBC strains isolated from culture positive pulmonary tuberculosis (TB) cases from the main TB referral centre between November 2007 and October 2009. Strain genotyping (IS6110 typing, spoligotyping, 24-loci MIRU-VNTR and SNP typing) was performed and demographic information collected. Among 73 MTBC strains analysed, we identified seven phylogenetic lineages, dominated by Beijing strains (49%). Beijing strains were further differentiated in two main branches, Beijing-A (n = 8) and -B (n = 28), that show distinct genotyping patterns and are characterized by specific deletion profiles (Beijing A: only RD105, RD207 deleted; Beijing B: RD150 and RD181 additionally deleted). Many Kiribati strains (59% based on IS6110 typing of all strains) occurred in clusters, suggesting ongoing local transmission. Beijing-B strains and over-crowded living conditions were associated with strain clustering (likely recent transmission), however little evidence of anti-tuberculous drug resistance was observed. We suggest enhanced case finding amongst close contacts and continued supervised treatment of all identified cases using standard first-line drugs to reduce TB burden in Kiribati. Beijing strains can be subdivided in different principle branches that might be associated with differential spreading patterns in the population.     

SNP/RD typing of Mycobacterium tuberculosis Beijing strains reveals local and worldwide disseminated clonal complexes

The Beijing strain is one of the most successful genotypes of Mycobacterium tuberculosis worldwide and appears to be highly homogenous according to existing genotyping methods. To type Beijing strains reliably we developed a robust typing scheme using single nucleotide polymorphisms (SNPs) and regions of difference (RDs) derived from whole-genome sequencing data of eight Beijing strains. SNP/RD typing of 259 M. tuberculosis isolates originating from 45 countries worldwide discriminated 27 clonal complexes within the Beijing genotype family. A total of 16 Beijing clonal complexes contained more than one isolate of known origin, of which two clonal complexes were strongly associated with South African origin. The remaining 14 clonal complexes encompassed isolates from different countries. Even highly resolved clonal complexes comprised isolates from distinct geographical sites. Our results suggest that Beijing strains spread globally on multiple occasions and that the tuberculosis epidemic caused by the Beijing genotype is at least partially driven by modern migration patterns. The SNPs and RDs presented in this study will facilitate future molecular epidemiological and phylogenetic studies on Beijing strains.     

Tuberculous Lymphadenitis in Ethiopia Predominantly Caused by Strains Belonging to the Delhi/CAS Lineage and Newly Identified Ethiopian Clades of the Mycobacterium tuberculosis Complex

Background

Recently, newly defined clades of Mycobacterium tuberculosis complex (MTBC) strains, namely Ethiopia 1–3 and Ethiopia H37Rv-like strains, and other clades associated with pulmonary TB (PTB) were identified in Ethiopia. In this study, we investigated whether these new strain types exhibit an increased ability to cause TB lymphadenitis (TBLN) and raised the question, if particular MTBC strains derived from TBLN patients in northern Ethiopia are genetically adapted to their local hosts and/or to the TBLN.

Methods

Genotyping of 196 MTBC strains isolated from TBLN patients was performed by spoligotyping and 24-loci mycobacterial interspersed repetitive unit-variable number of tandem repeats (MIRU-VNTR) typing. A statistical analysis was carried out to see possible associations between patient characteristics and phylogenetic MTBC strain classification.

Results

Among 196 isolates, the majority of strains belonged to the Delhi/CAS (38.8%) lineage, followed by Ethiopia 1 (9.7%), Ethiopia 3 (8.7%), Ethiopia H37RV-like (8.2%), Ethiopia 2 and Haarlem (7.7% each), URAL (3.6%), Uganda l and LAM (2% each), S-type (1.5%), X-type (1%), and 0.5% isolates of TUR, EAI, and Beijing genotype, respectively. Overall, 15 strains (7.7%) could not be allocated to a previously described phylogenetic lineage. The distribution of MTBC lineages is similar to that found in studies of PTB samples. The cluster rate (35%) in this study is significantly lower (P = 0.035) compared to 45% in the study of PTB in northwestern Ethiopia.

Conclusion

In the studied area, lymph node samples are dominated by Dehli/CAS genotype strains and strains of largely not yet defined clades based on MIRU-VNTR 24-loci nomenclature. We found no indication that strains of particular genotypes are specifically associated with TBLN. However, a detailed analysis of specific genetic variants of the locally contained Ethiopian clades by whole genome sequencing may reveal new insights into the host-pathogen co-evolution and specific features that are related to the local host immune system.     

A multiplex bead-based suspension array assay for interrogation of phylogenetically informative single nucleotide polymorphisms for Bacillus anthracis

Single nucleotide polymorphisms (SNPs) are abundant in genomes of all species and represent informative DNA markers extensively used to analyze phylogenetic relationships between strains. Medium to high throughput, open methodologies able to test many SNPs in a minimum time are therefore in great need. By using the versatile Luminex® xTAG technology, we developed an efficient multiplexed SNP genotyping assay to score 13 phylogenetically informative SNPs within the genome of Bacillus anthracis. The Multiplex Oligonucleotide Ligation-PCR procedure (MOL-PCR) described by Deshpande et al., 2010 has been modified and adapted for simultaneous interrogation of 13 biallelic canonical SNPs in a 13-plex assay. Changes made to the originally published method include the design of allele-specific dual-priming-oligonucleotides (DPOs) as competing detection probes (MOLigo probes) and use of asymmetric PCR reaction for signal amplification and labeling of ligation products carrying SNP targets. These innovations significantly reduce cross-reactivity observed when initial MOLigo probes were used and enhance hybridization efficiency onto the microsphere array, respectively. When evaluated on 73 representative samples, the 13-plex assay yielded unambiguous SNP calls and lineage affiliation. Assay limit of detection was determined to be 2 ng of genomic DNA. The reproducibility, robustness and easy-of-use of the present method were validated by a small-scale proficiency testing performed between four European laboratories. While cost-effective compared to other singleplex methods, the present MOL-PCR method offers a high degree of flexibility and scalability. It can easily accommodate newly identified SNPs to increase resolving power to the canSNP typing of B. anthracis.     

Origin, spread and demography of the Mycobacterium tuberculosis complex

The evolutionary timing and spread of the Mycobacterium tuberculosis complex (MTBC), one of the most successful groups of bacterial pathogens, remains largely unknown. Here, using mycobacterial tandem repeat sequences as genetic markers, we show that the MTBC consists of two independent clades, one composed exclusively of M. tuberculosis lineages from humans and the other composed of both animal and human isolates. The latter also likely derived from a human pathogenic lineage, supporting the hypothesis of an original human host. Using Bayesian statistics and experimental data on the variability of the mycobacterial markers in infected patients, we estimated the age of the MTBC at 40,000 years, coinciding with the expansion of "modern" human populations out of Africa. Furthermore, coalescence analysis revealed a strong and recent demographic expansion in almost all M. tuberculosis lineages, which coincides with the human population explosion over the last two centuries. These findings thus unveil the dynamic dimension of the association between human host and pathogen populations.     

Genomic Diversity among Drug Sensitive and Multidrug Resistant Isolates of Mycobacterium tuberculosis with Identical DNA Fingerprints

Background

Mycobacterium tuberculosis complex (MTBC), the causative agent of tuberculosis (TB), is characterized by low sequence diversity making this bacterium one of the classical examples of a genetically monomorphic pathogen. Because of this limited DNA sequence variation, routine genotyping of clinical MTBC isolates for epidemiological purposes relies on highly discriminatory DNA fingerprinting methods based on mobile and repetitive genetic elements. According to the standard view, isolates exhibiting the same fingerprinting pattern are considered direct progeny of the same bacterial clone, and most likely reflect ongoing transmission or disease relapse within individual patients.

Methodology/Principal Findings

Here we further investigated this assumption and used massively parallel whole-genome sequencing to compare one drug-susceptible (K-1) and one multidrug resistant (MDR) isolate (K-2) of a rapidly spreading M. tuberculosis Beijing genotype clone from a high incidence region (Karakalpakstan, Uzbekistan). Both isolates shared the same IS6110 RFLP pattern and the same allele at 23 out of 24 MIRU-VNTR loci.We generated 23.9 million (K-1) and 33.0 million (K-2) paired 50 bp purity filtered reads corresponding to a mean coverage of 483.5 fold and 656.1 fold respectively. Compared with the laboratory strain H37Rv both Beijing isolates shared 1,209 SNPs. The two Beijing isolates differed by 130 SNPs and one large deletion. The susceptible isolate had 55 specific SNPs, while the MDR variant had 75 specific SNPs, including the five known resistance-conferring mutations.

Conclusions

Our results suggest that M. tuberculosis isolates exhibiting identical DNA fingerprinting patterns can harbour substantial genomic diversity. Because this heterogeneity is not captured by traditional genotyping of MTBC, some aspects of the transmission dynamics of tuberculosis could be missed or misinterpreted. Furthermore, a valid differentiation between disease relapse and exogenous reinfection might be impossible using standard genotyping tools if the overall diversity of circulating clones is limited. These findings have important implications for clinical trials of new anti-tuberculosis drugs.     

Construction of a virtual Mycobacterium tuberculosis consensus genome and its application to data from a next generation sequencer

Background

Although Mycobacterium tuberculosis isolates are consisted of several different lineages and the epidemiology analyses are usually assessed relative to a particular reference genome, M. tuberculosis H37Rv, which might introduce some biased results. Those analyses are essentially based genome sequence information of M. tuberculosis and could be performed in sillico in theory, with whole genome sequence (WGS) data available in the databases and obtained by next generation sequencers (NGSs). As an approach to establish higher resolution methods for such analyses, whole genome sequences of the M. tuberculosis complexes (MTBCs) strains available on databases were aligned to construct virtual reference genome sequences called the consensus sequence (CS), and evaluated its feasibility in in sillico epidemiological analyses.

Results

The consensus sequence (CS) was successfully constructed and utilized to perform phylogenetic analysis, evaluation of read mapping efficacy, which is crucial for detecting single nucleotide polymorphisms (SNPs), and various MTBC typing methods virtually including spoligotyping, VNTR, Long sequence polymorphism and Beijing typing. SNPs detected based on CS, in comparison with H37Rv, were utilized in concatemer-based phylogenetic analysis to determine their reliability relative to a phylogenetic tree based on whole genome alignment as the gold standard. Statistical comparison of phylogenic trees based on CS with that of H37Rv indicated the former showed always better results that that of later. SNP detection and concatenation with CS was advantageous because the frequency of crucial SNPs distinguishing among strain lineages was higher than those of H37Rv. The number of SNPs detected was lower with the consensus than with the H37Rv sequence, resulting in a significant reduction in computational time. Performance of each virtual typing was satisfactory and accorded with those published when those are available.

Conclusions

These results indicated that virtual CS constructed from genome sequence data is an ideal approach as a reference for MTBC studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1368-9) contains supplementary material, which is available to authorized users.     

A First Insight on the Population Structure of Mycobacterium tuberculosis Complex as Studied by Spoligotyping and MIRU-VNTRs in Santiago,Chile

Tuberculosis (TB) remains a significant public health problem worldwide, but the ecology of the prevalent mycobacterial strains, and their transmission, can vary depending on country and region. Chile is a country with low incidence of TB, that has a geographically isolated location in relation to the rest of South American countries due to the Andes Mountains, but recent migration from neighboring countries has changed this situation. We aimed to assess the genotypic diversity of Mycobacterium tuberculosis complex (MTBC) strains in Santiago, Chile, and compare with reports from other Latin-American countries. We analyzed MTBC isolates from pulmonary tuberculosis cases collected between years 2008 and 2013 in Central Santiago, using two genotyping methods: spoligotyping and 12-loci mycobacterial interspersed repetitive unit-variable number of tandem repeats (MIRU-VNTRs). Data obtained were analyzed and compared to the SITVIT2 database. Mean age of the patients was 47.5 years and 61% were male; 11.6% were migrants. Of 103 strains (1 isolate/patient) included, there were 56 distinct spoligotype patterns. Of these, 16 strains (15.5%) corresponded to orphan strains in the SITVIT2 database, not previously reported. Latin American and Mediterranean (LAM) (34%) and T (33%) lineages were the most prevalent strains, followed by Haarlem lineage (16.5%). Beijing family was scarcely represented with only two cases (1.9%), one of them isolated from a Peruvian migrant. The most frequent clustered spoligotypes were SIT33/LAM3 (10.7%), SIT53/T1 (8.7%), SIT50/H3 (7.8%), and SIT37/T3 (6.8%). We conclude that LAM and T genotypes are the most prevalent genotypes of MTBC in Santiago, Chile, and together correspond to almost two thirds of analyzed strains, which is similar to strain distribution reported from other countries of Latin America. Nevertheless, the high proportion of SIT37/T3, which was rarely found in other Latin American countries, may underline a specific history or demographics of Chile related to probable human migrations and evolutions.     

High Sequence Variability of the ppE18 Gene of Clinical Mycobacterium tuberculosis Complex Strains Potentially Impacts Effectivity of Vaccine Candidate M72/AS01E

The development of an effective vaccine is urgently needed to fight tuberculosis (TB) which is still the leading cause of death from a single infectious agent worldwide. One of the promising vaccine candidates M72/AS01E consists of two proteins subunits PepA and PPE18 coded by Rv0125 and Rv1196. However, preliminary data indicate a high level of sequence variability among clinical Mycobacterium tuberculosis complex (MTBC) strains that might have an impact on the vaccine efficacy. To further investigate this finding, we determined ppE18 sequence variability in a well-characterized reference collection of 71 MTBC strains from 23 phylogenetic lineages representing the global MTBC diversity. In total, 100 sequence variations consisting of 96 single nucleotide polymorphisms (SNPs), three insertions and one deletion were detected resulting in 141 variable positions distributed over the entire gene. The majority of SNPs detected were non-synonymous (n = 68 vs. n = 28 synonymous). Strains from animal adapted lineages, e.g., M. bovis, showed a significant higher diversity than the human pathogens such as M. tuberculosis Haarlem. SNP patterns specific for different lineages as well as for deeper branches in the phylogeny could be identified. The results of our study demonstrate a high variability of the ppE18 gene even in the N-terminal domains that is normally highly conserved in ppe genes. As the N-terminal region interacts with TLR2 receptor inducing a protective anti-inflammatory immune response, genetic heterogeneity has a potential impact on the vaccine efficiency, however, this has to be investigated in future studies.     

Can 15-locus mycobacterial interspersed repetitive unit-variable-number tandem repeat analysis provide insight into the evolution of Mycobacterium tuberculosis?

The phylogeny and evolution of the bacterium Mycobacterium tuberculosis is still poorly understood despite the application of a variety of molecular techniques. We analyzed 469 M. tuberculosis and 49 Mycobacterium bovis isolates to evaluate if the mycobacterial interspersed repetitive units-variable-number tandem repeats (MIRU-VNTR) commonly used for epidemiological studies can define the phylogeny of the M. tuberculosis complex. This population was characterized by previously identified silent single-nucleotide polymorphisms (sSNPs) or by a macroarray based on these sSNPs that was developed in this study. MIRU-VNTR phylogenetic codes capable of differentiating between phylogenetic lineages were identified. Overall, there was 90.9% concordance between the lineages of isolates as defined by the MIRU-VNTR and sSNP analyses. The MIRU-VNTR phylogenetic code was unique to M. bovis and was not observed in any M. tuberculosis isolates. The codes were able to differentiate between different M. tuberculosis strain families such as Beijing, Delhi, and East African-Indian. Discrepant isolates with similar but not identical MIRU-VNTR codes often displayed a stepwise trend suggestive of bidirectional evolution. A lineage-specific panel of MIRU-VNTR can be used to subdivide each lineage for epidemiological purposes. MIRU-VNTR is a valuable tool for phylogenetic studies and could define an evolutionarily uncharacterized population of M. tuberculosis complex organisms.     

Host-pathogen coevolution in human tuberculosis

Tuberculosis (TB) is a disease of antiquity. Yet TB today still causes more adult deaths than any other single infectious disease. Recent studies show that contrary to the common view postulating an animal origin for TB, Mycobacterium tuberculosis complex (MTBC), the causative agent of TB, emerged as a human pathogen in Africa and colonized the world accompanying the Out-of-Africa migrations of modern humans. More recently, evolutionarily 'modern' lineages of MTBC expanded as a consequence of the global human population increase, and spread throughout the world following waves of exploration, trade and conquest. While epidemiological data suggest that the different phylogenetic lineages of MTBC might have adapted to different human populations, overall, the phylogenetically 'modern' MTBC lineages are more successful in terms of their geographical spread compared with the 'ancient' lineages. Interestingly, the global success of 'modern' MTBC correlates with a hypo-inflammatory phenotype in macrophages, possibly reflecting higher virulence, and a shorter latency in humans. Finally, various human genetic variants have been associated with different MTBC lineages, suggesting an interaction between human genetic diversity and MTBC variation. In summary, the biology and the epidemiology of human TB have been shaped by the long-standing association between MTBC and its human host.     

Significance of the Identification in the Horn of Africa of an Exceptionally Deep Branching Mycobacterium tuberculosis Clade

Molecular and phylogeographic studies have led to the definition within the Mycobacterium tuberculosis complex (MTBC) of a number of geotypes and ecotypes showing a preferential geographic location or host preference. The MTBC is thought to have emerged in Africa, most likely the Horn of Africa, and to have spread worldwide with human migrations. Under this assumption, there is a possibility that unknown deep branching lineages are present in this region. We genotyped by spoligotyping and multiple locus variable number of tandem repeats (VNTR) analysis (MLVA) 435 MTBC isolates recovered from patients. Four hundred and eleven isolates were collected in the Republic of Djibouti over a 12 year period, with the other 24 isolates originating from neighbouring countries. All major M. tuberculosis lineages were identified, with only two M. africanum and one M. bovis isolates. Upon comparison with typing data of worldwide origin we observed that several isolates showed clustering characteristics compatible with new deep branching. Whole genome sequencing (WGS) of seven isolates and comparison with available WGS data from 38 genomes distributed in the different lineages confirms the identification of ancestral nodes for several clades and most importantly of one new lineage, here referred to as lineage 7. Investigation of specific deletions confirms the novelty of this lineage, and analysis of its precise phylogenetic position indicates that the other three superlineages constituting the MTBC emerged independently but within a relatively short timeframe from the Horn of Africa. The availability of such strains compared to the predominant lineages and sharing very ancient ancestry will open new avenues for identifying some of the genetic factors responsible for the success of the modern lineages. Additional deep branching lineages may be readily and efficiently identified by large-scale MLVA screening of isolates from sub-Saharan African countries followed by WGS analysis of a few selected isolates.     

Microarrays for public health: genomic epidemiology of tuberculosis

In response to a large local school-based outbreak of tuberculosis, we have been evaluating the utility of microarray bacterial genomic analysis in outbreak management. After initial comparison of the isolate from the index case with Mycobacterium tuberculosis H37Rv, it was possible to design robust PCRs directed towards strain-specific deletions. Rapid PCR analysis of isolates proved valuable in determining whether or not other isolates were compatible with the outbreak strain and further microarray studies revealed genetic markers that could be used to discriminate between locally circulating strains.We suggest that this approach forms the basis for developing rapid local genotyping schemes applicable to M. tuberculosis and that application to other pathogens warrants consideration.     

Combined Genotypic,Phylogenetic, and Epidemiologic Analyses of Mycobacterium tuberculosis Genetic Diversity in the Rh?ne Alpes Region,France

BackgroundThe present work relates to identification and a deep molecular characterization of circulating Mycobacterium tuberculosis complex (MTBC) strains in the Rhône-Alpes region, France from 2000 to 2010. It aimed to provide with a first snapshot of MTBC genetic diversity in conjunction with bacterial drug resistance, type of disease and available demographic and epidemiologic characteristics over an eleven-year period, in the south-east of France.MethodsMycobacterium tuberculosis complex (MTBC) strains isolated in the Rhône-Alpes region, France (n = 2257, 1 isolate per patient) between 2000 and 2010 were analyzed by spoligotyping. MIRU-VNTR typing was applied on n = 1698 strains (with full results available for 974 strains). The data obtained were compared with the SITVIT2 database, followed by detailed genotyping, phylogenetic, and epidemiologic analyses in correlation with anonymized data on available demographic, and epidemiologic characteristics, and location of disease (pulmonary or extrapulmonary TB).ResultsThe most predominant spoligotyping clusters were SIT53/T1 (n = 346, 15.3%) > SIT50/H3 (n = 166, 7.35%) > SIT42/LAM9 (n = 125, 5.5%) > SIT1/Beijing (n = 72, 3.2%) > SIT47/H1 (n = 71, 3.1%). Evolutionary-recent strains belonging to the Principal Genetic Group (PGG) 2/3, or Euro-American lineages (T, LAM, Haarlem, X, S) were predominant and represented 1768 or 78.33% of all isolates. For strains having drug resistance information (n = 1119), any drug resistance accounted for 14.83% cases vs. 1.52% for multidrug resistance (MDR); and was significantly more associated with age group 21–40 years (p-value<0.001). Extra-pulmonary TB was more common among female patients while pulmonary TB predominated among men (p-value<0.001; OR = 2.16 95%CI [1.69; 2.77]). Also, BOV and CAS lineages were significantly well represented in patients affected by extra-pulmonary TB (p-value<0.001). The origin was known for 927/2257 patients: 376 (40.6%) being French-born vs. 551 (59.4%) Foreign-born. French patients were significantly older (mean age: 58.42 yrs 95%CI [56.04; 60.80]) than Foreign-born patients (mean age: 42.38 yrs. 95%CI [40.75; 44.0]).ConclusionThe study underlined the importance of imported TB cases on the genetic diversity and epidemiologic characteristics of circulating MTBC strains in Rhône-Alpes region, France over a large time-period. It helps better understand intricate relationships between certain lineages and geographic origin of the patients, and pinpoints genotypic and phylogenetic specificities of prevailing MTBC strains. Lastly, it also demonstrated a slow decline in isolation of M. africanum lineage in this region between 2000 and 2010.     

Combined Species Identification, Genotyping, and Drug Resistance Detection of Mycobacterium tuberculosis Cultures by MLPA on a Bead-Based Array

The population structure of Mycobacterium tuberculosis is typically clonal therefore genotypic lineages can be unequivocally identified by characteristic markers such as mutations or genomic deletions. In addition, drug resistance is mainly mediated by mutations. These issues make multiplexed detection of selected mutations potentially a very powerful tool to characterise Mycobacterium tuberculosis. We used Multiplex Ligation-dependent Probe Amplification (MLPA) to screen for dispersed mutations, which can be successfully applied to Mycobacterium tuberculosis as was previously shown. Here we selected 47 discriminative and informative markers and designed MLPA probes accordingly to allow analysis with a liquid bead array and robust reader (Luminex MAGPIX technology). To validate the bead-based MLPA, we screened a panel of 88 selected strains, previously characterised by other methods with the developed multiplex assay using automated positive and negative calling. In total 3059 characteristics were screened and 3034 (99.2%) were consistent with previous molecular characterizations, of which 2056 (67.2%) were directly supported by other molecular methods, and 978 (32.0%) were consistent with but not directly supported by previous molecular characterizations. Results directly conflicting or inconsistent with previous methods, were obtained for 25 (0.8%) of the characteristics tested. Here we report the validation of the bead-based MLPA and demonstrate its potential to simultaneously identify a range of drug resistance markers, discriminate the species within the Mycobacterium tuberculosis complex, determine the genetic lineage and detect and identify the clinically most relevant non-tuberculous mycobacterial species. The detection of multiple genetic markers in clinically derived Mycobacterium tuberculosis strains with a multiplex assay could reduce the number of TB-dedicated screening methods needed for full characterization. Additionally, as a proportion of the markers screened are specific to certain Mycobacterium tuberculosis lineages each profile can be checked for internal consistency. Strain characterization can allow selection of appropriate treatment and thereby improve treatment outcome and patient management.