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

A Guide to Phylogenomic Inference

Phylogenomics aims at reconstructing the evolutionary histories of organisms taking into account whole genomes or large fractions of genomes. Phylogenomics has significant applications in fields such as evolutionary biology, systematics, comparative genomics, and conservation genetics, providing valuable insights into the origins and relationships of species and contributing to our understanding of biological diversity and evolution. This chapter surveys phylogenetic concepts and methods aimed at both gene tree and species tree reconstruction while also addressing common pitfalls, providing references to relevant computer programs. A practical phylogenomic analysis example including bacterial genomes is presented at the end of the chapter.

Genomic epidemiological analysis identifies high relapse among individuals with recurring tuberculosis and provides evidence of recent household-related transmission of tuberculosis in Ghana

•

Unresolved previous infection as major cause of recurring tuberculosis (TB) in Ghana.

•

Genomic epidemiology identifies high relapse among recurrent TB cases in Ghana.

•

15-locus MIRU-VNTR typing is sufficient to predict the cause of TB recurrence.

•

Evidence of recent household-related TB transmission in Ghana.

•

Need for increased education by national TB control program.

Objective

To retrospectively investigate the cause of recurring tuberculosis (rcTB) among participants with pulmonary TB recruited from a prospective population-based study conducted between July 2012 and December 2015.

Methods

Mycobacterium tuberculosis complex isolates obtained from rcTB cases were characterized by standard mycobacterial genotyping tools, whole-genome sequencing, and phylogenetic analysis carried out to assess strain relatedness.

Results

The majority (58.3%, 21/36) of study participants with rcTB episodes had TB recurrence within 12 months post treatment. TB strains with isoniazid (INH) resistance were found in 19.4% (7/36) of participants at the primary episode, of which 29% (2/7) were also rifampicin-resistant. On TB recurrence, an INH-resistant strain was found in a larger proportion of participants, 27.8% (10/36), of which 40% (4/10) were MDR-TB strains. rcTB was attributed to relapse (same strain) in 75.0% (27/36) of participants and 25.0% (9/36) to re-infection.

Conclusion

Our findings indicate that previous unresolved infectiondue to inadequate treatment, may be the major cause of rcTB.

Mycobacterium tuberculosis transmission in an ethnically-diverse high incidence region in England, 2007-11

BACKGROUND

Transmission patterns in high tuberculosis incidence areas in England are poorly understood but need elucidating to focus contact tracing. We study transmission within and between age, ethnic and immigrant groups using molecular data from the high incidence West Midlands region.

METHODS

Isolates from culture-confirmed tuberculosis cases during 2007-2011 were typed using 24-locus Mycobacterial Interspersed Repetitive Unit-Variable Number Tandem Repeats (MIRU-VNTR). We estimated the proportion of disease attributable to recent transmission, calculated the proportion of isolates matching those from the two preceding years ("retrospectively clustered"), and identified risk factors for retrospective clustering using multivariate analyses. We calculated the ratio (RCR) between the observed and expected proportion clustered retrospectively within or between age, ethnic and immigrant groups.

RESULTS

Of the 2159 available genotypes (79% of culture-confirmed cases), 34% were attributed to recent transmission. The percentage retrospectively clustered decreased from 50 to 24% for 0-14 and ≥ 65 year olds respectively (p = 0.01) and was significantly lower for immigrants than the UK-born. Higher than expected clustering occurred within 15-24 year olds (RCR: 1.4 (95% CI: 1.1-1.8)), several ethnic groups, and between UK-born or long-term immigrants with the UK-born (RCR: 1.8 (95% CI: 1.1-2.4) and 1.6 (95% CI: 1.2-1.9) respectively).

CONCLUSIONS

This study is the first to consider "who clusters with whom" in a high incidence area in England, laying the foundation for future whole-genome sequencing work. The higher than expected clustering seen here suggests that preferential mixing between some age, ethnic and immigrant groups occurs; prioritising contact tracing to groups with which cases are most likely to cluster retrospectively could improve TB control.

Epidemiology of Mycobacterium bovis Disease in Humans in England, Wales, and Northern Ireland, 2002-2014

Despite slightly increased cases in these areas, human infection with this cattle pathogen remains rare.

Despite control efforts, Mycobacterium bovis incidence among cattle remains high in parts of England, Wales, and Northern Ireland, attracting political and public health interest in potential spread from animals to humans. To determine incidence among humans and to identify associated factors, we conducted a retrospective cohort analysis of human M. bovis cases in England, Wales, and Northern Ireland during 2002–2014. We identified 357 cases and observed increased annual case numbers (from 17 to 35) and rates. Most patients were >65 years of age and born in the United Kingdom. The median age of UK-born patients decreased over time. For 74% of patients, exposure to risk factors accounting for M. bovis acquisition, most frequently consumption of unpasteurized milk, was known. Despite the small increase in case numbers and reduction in patient age, M. bovis infection of humans in England, Wales, and Northern Ireland remains rare.

Phylogenomics

Phylogenomics aims at reconstructing the evolutionary histories of organisms taking into account whole genomes or large fractions of genomes. The abundance of genomic data for an enormous variety of organisms has enabled phylogenomic inference of many groups, and this has motivated the development of many computer programs implementing the associated methods. This chapter surveys phylogenetic concepts and methods aimed at both gene tree and species tree reconstruction while also addressing common pitfalls, providing references to relevant computer programs. A practical phylogenomic analysis example including bacterial genomes is presented at the end of the chapter.

Tuberculosis in migrants moving from high-incidence to low-incidence countries: a population-based cohort study of 519 955 migrants screened before entry to England, Wales, and Northern Ireland

BACKGROUND

Tuberculosis elimination in countries with a low incidence of the disease necessitates multiple interventions, including innovations in migrant screening. We examined a cohort of migrants screened for tuberculosis before entry to England, Wales, and Northern Ireland and tracked the development of disease in this group after arrival.

METHODS

As part of a pilot pre-entry screening programme for tuberculosis in 15 countries with a high incidence of the disease, the International Organization for Migration screened all applicants for UK visas aged 11 years or older who intended to stay for more than 6 months. Applicants underwent a chest radiograph, and any with results suggestive of tuberculosis underwent sputum testing and culture testing (when available). We tracked the development of tuberculosis in those who tested negative for the disease and subsequently migrated to England, Wales, and Northern Ireland with the Enhanced Tuberculosis Surveillance system. Primary outcomes were cases of all forms of tuberculosis (including clinically diagnosed cases), and bacteriologically confirmed pulmonary tuberculosis.

FINDINGS

Our study cohort was 519 955 migrants who were screened for tuberculosis before entry to the UK between Jan 1, 2006, and Dec 31, 2012. Cases notified on the Enhanced Tuberculosis Surveillance system between Jan 1, 2006, and Dec 31, 2013, were included. 1873 incident cases of all forms of tuberculosis were identified, and, on the basis of data for England, Wales, and Northern Ireland, the estimated incidence of all forms of tuberculosis in migrants screened before entry was 147 per 100 000 person-years (95% CI 140-154). The estimated incidence of bacteriologically confirmed pulmonary tuberculosis in migrants screened before entry was 49 per 100 000 person-years (95% CI 45-53). Migrants whose chest radiographs were compatible with active tuberculosis but with negative pre-entry microbiological results were at increased risk of tuberculosis compared with those with no radiographic abnormalities (incidence rate ratio 3·2, 95% CI 2·8-3·7; p<0·0001). Incidence of tuberculosis after migration increased significantly with increasing WHO-estimated prevalence of tuberculosis in migrants' countries of origin. 35 of 318 983 pre-entry screened migrants included in a secondary analysis with typing data were assumed index cases. Estimates of the rate of assumed reactivation tuberculosis ranged from 46 (95% CI 42-52) to 91 (82-102) per 100 000 population.

INTERPRETATION

Migrants from countries with a high incidence of tuberculosis screened before being granted entry to low-incidence countries pose a negligible risk of onward transmission but are at increased risk of tuberculosis, which could potentially be prevented through identification and treatment of latent infection in close collaboration with a pre-entry screening programme.

FUNDING

Wellcome Trust, UK National Institute for Health Research, UK Medical Research Council, Public Health England, and Department of Health Policy Research Programme.

Transmission of multidrug-resistant tuberculosis in the UK: a cross-sectional molecular and epidemiological study of clustering and contact tracing

BACKGROUND

Between 2000 and 2012 the number of multidrug-resistant (MDR) tuberculosis cases in the UK increased from 28 per year to 81 per year. We investigated the proportion of MDR tuberculosis cases arising from transmission in the UK and associated risk factors.

METHOD

We identified patients with MDR tuberculosis notified in England, Wales, and Northern Ireland between Jan 1, 2004, and Dec 31, 2007, by linking national laboratory and surveillance data. Data for laboratory isolates, including drug sensitivities and 24-mycobacterial interspersed repetitive-unit-variable-number tandem repeat (MIRU-VNTR) typing were obtained routinely from the National Tuberculosis Reference laboratories as part of national tuberculosis surveillance. We investigated clusters of cases with indistinguishable MIRU-VNTR profiles to identify epidemiological links. We calculated transmission using the n-1 method and established associated risk factors by logistic regression. We also assessed the likelihood of transmission to additional secondary active tuberculosis cases, identified through conventional contact tracing.

FINDINGS

204 patients were diagnosed with MDR tuberculosis in the study period; 189 (92·6%) had an MIRU-VNTR profile. We identified 12 clusters containing 40 individuals and 149 unique strains. The proportion of cases attributable to recent transmission, on the basis of molecular data, was 15% (40 cases clustered-12 clusters/189 with a strain type). The proportion of cases attributable to recent transmission (ie, transmission within the UK) after adjustment for epidemiological links was 8·5% (22 cases with epidemiological links-six clusters/189 cases with a strain type). Being UK born (odds ratio 4·81; 95% CI 2·03-11·36, p=0·0005) and illicit drug use (4·75; 1·19-18·96, p=0·026) were significantly associated with clustering. The most common transmission setting was the household but 21 of 22 of epidemiological links were missed by conventional contact tracing. 13 secondary active tuberculosis cases identified by conventional contact tracing were mostly contacts of patients with MDR tuberculosis from countries of high tuberculosis burden. 11 (85%) of 13 shared the same country of birth as the index case, of whom ten did not share a strain type or drug resistance pattern.

INTERPRETATION

Transmission of MDR tuberculosis in the UK is low and associated with being UK born or illicit drug use. MIRU-VNTR typing with cluster investigation was more successful at identifying transmission events than conventional contact tracing. Individuals with tuberculosis who have had contact with a known MDR tuberculosis source case from a country of high tuberculosis burden should have drug-sensitivity testing on isolates to ensure appropriate treatment is given.

FUNDING

Public Health England.

Molecular epidemiology, drug susceptibility and economic aspects of tuberculosis in Mubende district, Uganda

Background

Tuberculosis (TB) remains a global public health problem whose effects have major impact in developing countries like Uganda. This study aimed at investigating genotypic characteristics and drug resistance profiles of Mycobacterium tuberculosis isolated from suspected TB patients. Furthermore, risk factors and economic burdens that could affect the current control strategies were studied.

Methods

TB suspected patients were examined in a cross-sectional study at the Mubende regional referral hospital between February and July 2011. A questionnaire was administered to each patient to obtain information associated with TB prevalence. Isolates of M. tuberculosis recovered during sampling were examined for drug resistance to first line anti-TB drugs using the BACTEC-MGIT960^TMsystem. All isolates were further characterized using deletion analysis, spoligotyping and MIRU-VNTR analysis. Data were analyzed using different software; MIRU-VNTR plus, SITVITWEB, BioNumerics and multivariable regression models.

Results

M. tuberculosis was isolated from 74 out of 344 patients, 48 of these were co-infected with HIV. Results from the questionnaire showed that previously treated TB, co-infection with HIV, cigarette smoking, and overcrowding were risk factors associated with TB, while high medical related transport bills were identified as an economic burden. Out of the 67 isolates that gave interpretable results, 23 different spoligopatterns were detected, nine of which were novel patterns. T2 with the sub types Uganda-I and Uganda-II was the most predominant lineage detected. Antibiotic resistance was detected in 19% and multidrug resistance was detected in 3% of the isolates.

Conclusion

The study detected M. tuberculosis from 21% of examined TB patients, 62% of whom were also HIV positive. There is a heterogeneous pool of genotypes that circulate in this area, with the T2 lineage being the most predominant. High medical related transport bills and drug resistance could undermine the usefulness of the current TB strategic interventions.

Whole-genome sequencing to delineate Mycobacterium tuberculosis outbreaks: a retrospective observational study

Background

Tuberculosis incidence in the UK has risen in the past decade. Disease control depends on epidemiological data, which can be difficult to obtain. Whole-genome sequencing can detect microevolution within Mycobacterium tuberculosis strains. We aimed to estimate the genetic diversity of related M tuberculosis strains in the UK Midlands and to investigate how this measurement might be used to investigate community outbreaks.

Methods

In a retrospective observational study, we used Illumina technology to sequence M tuberculosis genomes from an archive of frozen cultures. We characterised isolates into four groups: cross-sectional, longitudinal, household, and community. We measured pairwise nucleotide differences within hosts and between hosts in household outbreaks and estimated the rate of change in DNA sequences. We used the findings to interpret network diagrams constructed from 11 community clusters derived from mycobacterial interspersed repetitive-unit–variable-number tandem-repeat data.

Findings

We sequenced 390 separate isolates from 254 patients, including representatives from all five major lineages of M tuberculosis. The estimated rate of change in DNA sequences was 0·5 single nucleotide polymorphisms (SNPs) per genome per year (95% CI 0·3–0·7) in longitudinal isolates from 30 individuals and 25 families. Divergence is rarely higher than five SNPs in 3 years. 109 (96%) of 114 paired isolates from individuals and households differed by five or fewer SNPs. More than five SNPs separated isolates from none of 69 epidemiologically linked patients, two (15%) of 13 possibly linked patients, and 13 (17%) of 75 epidemiologically unlinked patients (three-way comparison exact p<0·0001). Genetic trees and clinical and epidemiological data suggest that super-spreaders were present in two community clusters.

Interpretation

Whole-genome sequencing can delineate outbreaks of tuberculosis and allows inference about direction of transmission between cases. The technique could identify super-spreaders and predict the existence of undiagnosed cases, potentially leading to early treatment of infectious patients and their contacts.

Funding

Medical Research Council, Wellcome Trust, National Institute for Health Research, and the Health Protection Agency.

Molecular epidemiology of Mycobacterium tuberculosis, Buenos Aires, Argentina

To analyze the molecular epidemiology of Mycobacterium tuberculosis strains at a hospital in Buenos Aires, Argentina, and mutations related to multidrug-resistant and extensively drug-resistant tuberculosis, we conducted a prospective case-control study. Our findings reinforce the value of incorporating already standardized molecular methods for rapidly detecting resistance.

Molecular epidemiology of Mycobacterium tuberculosis clinical isolates in Southwest Ireland

Tuberculosis has had significant effects on Ireland over the past two centuries, causing persistently higher morbidity and mortality than in neighbouring countries until the last decade. This study describes the results of genotyping and drug susceptibility testing of 171 strains of Mycobacterium tuberculosis complex isolated between January 2004 and December 2006 in a region of Ireland centred on the city of Cork. Spoligotype comparisons were made with the SpolDB4 database and clustered 130 strains in 23 groups, forty-one strains showed unique Spoligotyping patterns. The commonest spoligotypes detected were ST0137 (X2) (16.9%), and ST0351 (15.8%) ('U' clade). The major spoligotype clades were X (26.2%), U (19.3%), T (15.2%), Beijing (5.9%), Haarlem (4.7%), LAM (4.1%), BOVIS (1.75%), with 12.9% unassigned strains. A 24-locus VNTR genotyping produced 15 clusters containing 49 isolates, with high discrimination index (HGDI>0.99). A combination of Spoligotyping and VNTR reduced the number of clustered isolates to 47 in 15 clusters (27.5%). This study identified ST351 as common among Irish nationals, and found a low rate of drug resistance with little evidence of transmission of drug resistant strains. Strain clustering was significantly associated with age under 55 years and Irish nationality. Only strains of Euro-American lineage formed clusters. Molecular typing did not completely coincide with the results of contact investigations.

Associations between Mycobacterium tuberculosis strains and phenotypes

This population-based study was used to investigate strong associations between phenotypes and genotypes.

To inform development of tuberculosis (TB) control strategies, we characterized a total of 2,261 Mycobacterium tuberculosis complex isolates by using multiple phenotypic and molecular markers, including polymorphisms in repetitive sequences (spoligotyping and variable-number tandem repeats [VNTRs]) and large sequence and single-nucleotide polymorphisms. The Beijing family was strongly associated with multidrug resistance (p = 0.0001), and VNTR allelic variants showed strong associations with spoligotyping families: >5 copies at exact tandem repeat (ETR) A, >2 at mycobacterial interspersed repetitive unit 24, and >3 at ETR-B associated with the East African–Indian and M. bovis strains. All M. tuberculosis isolates were differentiated into 4 major lineages, and a maximum parsimony tree was constructed suggesting a more complex phylogeny for M. africanum. These findings can be used as a model of pathogen global diversity.

Typing Mycobacterium tuberculosis using variable number tandem repeat analysis

DNA-based typing has contributed to the understanding of M. tuberculosis epidemiology and evolution. IS6110 RFLP was the first method described and has been used in many epidemiologic investigations. Technological difficulties have hampered the widespread establishment of this method, and it has been found to be of little use in evolutionary studies. PCR-based methods such as spoligotyping and variable number tandem repeat (VNTR) analysis largely overcome these difficulties. Spoligotyping alone is of limited value in epidemiologic investigations due to low discrimination but can be useful in evolutionary studies. Panels of VNTR loci selected from the 59 polymorphic VNTRs described to date have been shown to be useful in both epidemiologic and evolutionary studies. A VNTR type is identified by, first, amplifying a series of PCR fragments each encompassing a different VNTR locus and, second, determining the PCR fragment sizes from which the number of repeats present is calculated. The repeat number present at a series of loci is used as numerical code to describe a type. This chapter describes a high-throughput automated method for VNTR analysis at 15 loci using a capillary fragment analyzer and a manual method using agarose gel analysis.

Discordance between mycobacterial interspersed repetitive-unit-variable-number tandem-repeat typing and IS6110 restriction fragment length polymorphism genotyping for analysis of Mycobacterium tuberculosis Beijing strains in a setting of high incidence of tuberculosis

IS6110 restriction fragment length polymorphism (RFLP) genotyping is the most widely used genotyping method to study the epidemiology of Mycobacterium tuberculosis. However, due to the complexity of the IS6110 RFLP genotyping technique, and the interpretation of RFLP data, mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) genotyping has been proposed as the new genotyping standard. This study aimed to determine the discriminatory power of different MIRU-VNTR locus combinations relative to IS6110 RFLP genotyping, using a collection of Beijing genotype M. tuberculosis strains with a well-established phylogenetic history. Clustering, diversity index, clustering concordance, concordance among unique genotypes, and divergent and convergent evolution were calculated for seven combinations of 27 different MIRU-VNTR loci and compared to IS6110 RFLP results. Our results confirmed previous findings that MIRU-VNTR genotyping can be used to estimate the extent of recent or ongoing transmission. However, molecular epidemiological linking of cases varied significantly depending on the genotyping method used. We conclude that IS6110 RFLP and MIRU-VNTR loci evolve independently and at different rates, which leads to discordance between transmission chains predicted by the respective genotyping methods. Concordance between the two genotyping methods could be improved by the inclusion of genetic distance (GD) into the clustering formulae for some of the MIRU-VNTR loci combinations. In summary, our findings differ from previous reports, which may be explained by the fact that in settings of low tuberculosis incidence, the genetic distance between epidemiologically unrelated isolates was sufficient to define a strain using either marker, whereas in settings of high incidence, continuous evolution and persistence of strains revealed the weaknesses inherent to these markers.

Molecular epidemiology and prevalence of mutations conferring rifampicin and isoniazid resistance in Mycobacterium tuberculosis strains from the southern Ukraine

Understanding the molecular epidemiology of tuberculosis (TB) and mutations in genes associated with drug resistance may contribute to the development of appropriate interventions to improve tuberculosis control. A structured questionnaire was used to collect basic epidemiological data from 589 patients with radiologically confirmed TB in the Odessa and Nikolaev regions of the Ukraine in 2003-2004. A non-commercial reverse hybridisation assay and DNA sequencing were used to detect mutations associated with rifampicin and isoniazid resistance. Genotyping was performed using multilocus variable number tandem repeat (VNTR) typing and spoligotyping. Mutations conferring rifampicin and isoniazid resistance were detected in 32.9% and 44.0%, respectively, of 225 Mycobacterium tuberculosis isolates from individual consecutive patients. Mutations in codon 531 and codon 315 of the rpoB and katG genes, respectively, were predominant among drug-resistant isolates. Multidrug (MDR) resistance rates were significantly higher among former prison inmates compared with non-prisoners (54.8% vs. 27.3%; RR 2.01; 95% CI 1.35-2.97) and the prevalence of mutations was higher in Beijing strains sharing the VNTR signature 223325173533424 than in other Beijing strains (71.4% vs. 45.7%; RR 1.74; 95% CI 1.17-2.57), suggesting that this group may be responsible for rapid transmission of MDR TB in the southern Ukraine.

Proposal for standardization of optimized mycobacterial interspersed repetitive unit-variable-number tandem repeat typing of Mycobacterium tuberculosis

Molecular typing based on 12 loci containing variable numbers of tandem repeats of mycobacterial interspersed repetitive units (MIRU-VNTRs) has been adopted in combination with spoligotyping as the basis for large-scale, high-throughput genotyping of Mycobacterium tuberculosis. However, even the combination of these two methods is still less discriminatory than IS6110 fingerprinting. Here, we define an optimized set of MIRU-VNTR loci with a significantly higher discriminatory power. The resolution and the stability/robustness of 29 loci were analyzed, using a total of 824 tubercle bacillus isolates, including representatives of the main lineages identified worldwide so far. Five loci were excluded for lack of robustness and/or stability in serial isolates or isolates from epidemiologically linked patients. The use of the 24 remaining loci increased the number of types by 40%--and by 23% in combination with spoligotyping--among isolates from cosmopolitan origins, compared to those obtained with the original set of 12 loci. Consequently, the clustering rate was decreased by fourfold--by threefold in combination with spoligotyping--under the same conditions. A discriminatory subset of 15 loci with the highest evolutionary rates was then defined that concentrated 96% of the total resolution obtained with the full 24-locus set. Its predictive value for evaluating M. tuberculosis transmission was found to be equal to that of IS6110 restriction fragment length polymorphism typing, as shown in a companion population-based study. This 15-locus system is therefore proposed as the new standard for routine epidemiological discrimination of M. tuberculosis isolates and the 24-locus system as a high-resolution tool for phylogenetic studies.

Bottlenecks and broomsticks: the molecular evolution of Mycobacterium bovis

Mycobacterium bovis is the cause of tuberculosis in cattle and is a member of the Mycobacterium tuberculosis complex. In contrast to many other pathogenic bacterial species, there is little evidence for the transfer and recombination of genes between cells. The clonality of this group of organisms indicates that the population structure is dominated by reductions in diversity, caused either by population bottlenecks or selective sweeps as entire chromosomes become fixed in the population. We describe how these forces have shaped not only the phylogeny of this group but also, at a very local level, the population structure of Mycobacterium bovis in the British Isles. We also discuss the practical implications of applying this knowledge to understanding the spread of infection and the development of improved vaccines and diagnostic tests.

Mycobacterium tuberculosis complex genetic diversity: mining the fourth international spoligotyping database (SpolDB4) for classification, population genetics and epidemiology

Background

The Direct Repeat locus of the Mycobacterium tuberculosis complex (MTC) is a member of the CRISPR (Clustered regularly interspaced short palindromic repeats) sequences family. Spoligotyping is the widely used PCR-based reverse-hybridization blotting technique that assays the genetic diversity of this locus and is useful both for clinical laboratory, molecular epidemiology, evolutionary and population genetics. It is easy, robust, cheap, and produces highly diverse portable numerical results, as the result of the combination of (1) Unique Events Polymorphism (UEP) (2) Insertion-Sequence-mediated genetic recombination. Genetic convergence, although rare, was also previously demonstrated. Three previous international spoligotype databases had partly revealed the global and local geographical structures of MTC bacilli populations, however, there was a need for the release of a new, more representative and extended, international spoligotyping database.

Results

The fourth international spoligotyping database, SpolDB4, describes 1939 shared-types (STs) representative of a total of 39,295 strains from 122 countries, which are tentatively classified into 62 clades/lineages using a mixed expert-based and bioinformatical approach. The SpolDB4 update adds 26 new potentially phylogeographically-specific MTC genotype families. It provides a clearer picture of the current MTC genomes diversity as well as on the relationships between the genetic attributes investigated (spoligotypes) and the infra-species classification and evolutionary history of the species. Indeed, an independent Naïve-Bayes mixture-model analysis has validated main of the previous supervised SpolDB3 classification results, confirming the usefulness of both supervised and unsupervised models as an approach to understand MTC population structure. Updated results on the epidemiological status of spoligotypes, as well as genetic prevalence maps on six main lineages are also shown. Our results suggests the existence of fine geographical genetic clines within MTC populations, that could mirror the passed and present Homo sapiens sapiens demographical and mycobacterial co-evolutionary history whose structure could be further reconstructed and modelled, thereby providing a large-scale conceptual framework of the global TB Epidemiologic Network.

Conclusion

Our results broaden the knowledge of the global phylogeography of the MTC complex. SpolDB4 should be a very useful tool to better define the identity of a given MTC clinical isolate, and to better analyze the links between its current spreading and previous evolutionary history. The building and mining of extended MTC polymorphic genetic databases is in progress.