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

Mycobacterium caprae Infection of Red Deer in Western Austria-Optimized Use of Pathology Data to Infer Infection Dynamics

Austria is officially bovine tuberculosis (TB) free, but during the last decade the west of the country experienced sporadic TB cases in cattle. Free-ranging red deer are known to be the maintenance host of Mycobacterium (M.) caprae in certain areas in Austria, where cattle can become infected on alpine pastures shared with deer. The epidemiology of TB in deer in alpine regions is still poorly understood. To inform decisions on efficient interventions against TB in deer, a method is needed to better capture the infection dynamics on population level. A total of 4,521 free-ranging red deer from Austria's most western Federal state Vorarlberg were TB-tested between 2009 and 2018. M. caprae was confirmed in samples from 257 animals. Based on descriptions of TB-like lesions, TB positive animals were categorized with a newly developed lesion score called “Patho Score.” Analyses using this Patho Score allowed us to distinguish between endemic, epidemic and sporadic TB situations and revealed different roles of subgroups of infected deer in infection dynamics. Overall, deer in poor condition, deer of older age and stags were the subgroups that were significantly more often TB positive (p = 0.02 or smaller for all subgroups). Deer in poor condition (p < 0.001) and stags (p = 0.04) also showed more often advanced lesions, indicating their role in mycobacterial spread. TB was never detected in fawns, while hinds were the subgroup that showed the fewest advanced lesions. Analysis of outbreaks of TB and lesion development in yearlings provided some evidence for the role of winter feeding as a source for increased infection transmission. Sporadic cases in TB-free areas appear to precede outbreaks in these areas. These currently TB-free areas should receive particular attention in sampling schemes to be able to detect early spreading of the infection. The Patho Score is a quick, easy-to-apply and reproducible tool that provides new insights on the epidemiology of TB in deer at population level and is flexible enough to relate heterogeneous wildlife monitoring data collected following different sampling plans. This lesion score was used for systematic assessment of infection dynamics of mycobacterial infections.

Acquired Resistance to Antituberculosis Drugs in England, Wales, and Northern Ireland, 2000-2015

Among tuberculosis (TB) patients, acquired resistance to anti-TB drugs represents a failure in the treatment pathway. To improve diagnosis and care for patients with drug-resistant TB, we examined the epidemiology and risk factors associated with acquired drug resistance during 2000–2015 among TB patients in England, Wales, and Northern Ireland. We found acquired resistance in 0.2% (158/67,710) of patients with culture-confirmed TB. Using multivariate logistic regression, we identified the following factors associated with acquired drug resistance: having pulmonary disease; initial resistance to isoniazid, rifampin, or both; a previous TB episode; and being born in China or South Africa. Treatment outcomes were worse for patients with than without acquired resistance. Although acquired resistance is rare in the study area, certain patient groups are at higher risk. Identifying these patients and ensuring that adequate resources are available for treatment may prevent acquisition of resistance, thereby limiting transmission of drug-resistant strains of mycobacteria.

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.

Cross-border outbreak of extensively drug-resistant tuberculosis linked to a university in Romania

Extensively drug-resistant (XDR) tuberculosis (TB) poses a threat to public health due to its complicated, expensive and often unsuccessful treatment. A cluster of three XDR TB cases was detected among foreign medical students of a Romanian university. The contact investigations included tuberculin skin testing or interferon gamma release assay, chest X-ray, sputum smear microscopy, culture, drug susceptibility testing, genotyping and whole-genome sequencing (WGS), and were addressed to students, personnel of the university, family members or other close contacts of the cases. These investigations increased the total number of cases to seven. All confirmed cases shared a very similar WGS profile. Two more cases were epidemiologically linked, but no laboratory confirmation exists. Despite all the efforts done, the source of the outbreak was not identified, but the transmission was controlled. The investigation was conducted by a team including epidemiologists and microbiologists from five countries (Finland, Israel, Romania, Sweden and the UK) and from the European Centre for Disease Prevention and Control. Our report shows how countries can collaborate to control the spread of XDR TB by exchanging information about cases and their contacts to enable identification of additional cases and transmission and to perform the source investigation.

Macro-geographical specificities of the prevailing tuberculosis epidemic as seen through SITVIT2, an updated version of the Mycobacterium tuberculosis genotyping database

In order to provide a global overview of genotypic, epidemiologic, demographic, phylogeographical, and drug resistance characteristics related to the prevailing tuberculosis (TB) epidemic, we hereby report an update of the 6th version of the international genotyping database SITVIT2. We also make all the available information accessible through a dedicated website (available at http://www.pasteur-guadeloupe.fr:8081/SITVIT2). Thanks to the public release of SITVIT2 which is currently the largest international multimarker genotyping database with a compilation of 111,635 clinical isolates from 169 countries of patient origin (131 countries of isolation, representing 1032 cities), our major aim is to highlight macro- and micro-geographical cleavages and phylogeographical specificities of circulating Mycobacterium tuberculosis complex (MTBC) clones worldwide. For this purpose, we retained strains typed by the most commonly used PCR-based methodology for TB genotyping, i.e., spoligotyping based on the polymorphism of the direct repeat (DR) locus, 5-loci Exact Tandem Repeats (ETRs), and MIRU-VNTR minisatellites used in 12-, 15-, or 24-loci formats. We describe the SITVIT2 database and integrated online applications that permit to interrogate the database using easy drop-down menus to draw maps, graphics and tables versus a long list of parameters and variables available for individual clinical isolates (year and place of isolation, origin, sex, and age of patient, drug-resistance, etc.). Available tools further allow to generate phylogenetical snapshot of circulating strains as Lineage-specific WebLogos, as well as minimum spanning trees of their genotypes in conjunction with their geographical distribution, drug-resistance, demographic, and epidemiologic characteristics instantaneously; whereas online statistical analyses let a user to pinpoint phylogeographical specificities of circulating MTBC lineages and conclude on actual demographic trends. Available associated information on gender (n = 18,944), age (n = 16,968), drug resistance (n = 19,606), and HIV serology (n = 2673), allowed to draw some important conclusions on TB geo-epidemiology; e.g. a positive correlation exists between certain Mycobacterium tuberculosis lineages (such as CAS and Beijing) and drug resistance (p-value<.001), while other lineages (such as LAM, X, and BOV) are more frequently associated with HIV-positive serology (p-value<.001). Besides, availability of information on the year of isolation of strains (range 1759-2012), also allowed to make tentative correlations between drug resistance information and lineages - portraying probable evolution trends over time and space. To conclude, the present approach of geographical mapping of predominant clinical isolates of tubercle bacilli causing the bulk of the disease both at country and regional level in conjunction with epidemiologic and demographic characteristics allows to shed new light on TB geo-epidemiology in relation with the continued waves of peopling and human migration.

Mycobacterial biomaterials and resources for researchers

There are many resources available to mycobacterial researchers, including culture collections around the world that distribute biomaterials to the general scientific community, genomic and clinical databases, and powerful bioinformatics tools. However, many of these resources may be unknown to the research community. This review article aims to summarize and publicize many of these resources, thus strengthening the quality and reproducibility of mycobacterial research by providing the scientific community access to authenticated and quality-controlled biomaterials and a wealth of information, analytical tools and research opportunities.

Genetic diversity of Mycobacterium tuberculosis in south coastal Karnataka, India, using spoligotyping

Background & objectives:

Despite high occurrence of tuberculosis in India very little information is available about the genetic diversity of Mycobacterium tuberculosis isolates prevailing in coastal Karnataka, India. Thus, the present study was undertaken to explore the genetic biodiversity of M. tuberculosis isolates prevailing in south coastal region of Karnataka (Udupi District), India.

Methods:

A total of 111 Mycobacterial isolates were cultured in Lowenstein Jensen (LJ) medium and after obtaining growth, DNA was extracted and spoligotyping was performed. SITVIT WEB database was used to locate families of spoligotypes.

Results:

On analyzing the hybridization results of all 111 isolates on SITVIT WEB database 57 (51.35%) isolates were clustered into 11 Spoligotype International Types (SIT). The largest cluster of 14 (12.61%) isolates was SIT-48 (EAI1-SOM), followed by SIT-1942 (CAS1-Delhi) with 11 isolates (9.9%) and SIT-11 with seven (6.30%). Moreover, 23 isolates (20.72%) had unique spoligotypes and 31 (27.92%) were orphans. Spotclust analysis revealed that majority (67%) of orphan isolates were variants of CAS (37%) and EAI-5 (34%).

Interpretation & conclusions:

The present study revealed high biodiversity among the circulating isolates of M. tuberculosis in this region with the presence of mixed genotypes earlier reported from north and south India along with certain new genotypes with unique SITs. The study highlights the need for further longitudinal studies to explore the genetic diversity and to understand the transmission dynamics of prevailing isolates.

Early ancient sublineages of Mycobacterium tuberculosis Beijing genotype: unexpected clues from phylogenomics of the pathogen and human history

OBJECTIVE

The Mycobacterium tuberculosis Beijing family is an epidemiologically important lineage subdivided into large-scale phylogenetic sublineages: ancient, endemic in East Asia, and global modern. Here, we analysed ancient sublineages of the Beijing genotype in the Omsk region of southwestern Siberia, an intriguing area at the intersection of European Russia, Siberia, and Central Asia.

METHODS

The study included 423 M. tuberculosis strains isolated in 2013-2017 and subjected to drug susceptibility testing, genotyping, and whole genome sequencing.

RESULTS

The Beijing genotype constituted 280 out of 423 strains. Forty Beijing strains belonged to the early ancient sublineage (wild type mutT4-48). Of these, 11 belonged to the 14717-15 MIRU-VNTR cluster and had intact RD181, 29 belonged to the 1071-32 cluster and had the RD181 deletion. Thirty-nine ancient strains were multidrug-resistant (MDR) and 20 pre-extensively drug resistant (XDR)/XDR. Comparison with global data demonstrated that these clones circulate mainly in Asian Russia with certain phylogenetic affinity to strains from Japan, Korea, and northeastern China. The genome-wide analysis revealed 29-37 single nucleotide polymorphism distances between isolates from different Russian regions within these two clusters.

CONCLUSIONS

Based on phylogenetic, phylogeographic, genomic, and historical data, we hypothesize that these two clones or their direct ancestors were probably brought to Russia ∼70 years ago after the Second World War with Japanese prisoners of war and, until recently, were mainly circulating in Siberia and the Far East. Their elevated prevalence in Omsk along with the extremely strong association with not only MDR but also pre-XDR/XDR also observed in other locations highlight their epidemic potential and the need for monitoring and attention from health authorities.

Genomic sequencing is required for identification of tuberculosis transmission in Hawaii

Background

Tuberculosis (TB) caused an estimated 1.4 million deaths and 10.4 million new cases globally in 2015. TB rates in the United States continue to steadily decline, yet rates in the State of Hawaii are perennially among the highest in the nation due to a continuous influx of immigrants from the Western Pacific and Asia. TB in Hawaii is composed of a unique distribution of genetic lineages, with the Beijing and Manila families of Mycobacterium tuberculosis (Mtb) comprising over two-thirds of TB cases. Standard fingerprinting methods (spoligotyping plus 24-loci Mycobacterial Interspersed Repetitive Units-Variable Number Tandem Repeats [MIRU-VNTR] fingerprinting) perform poorly when used to identify actual transmission clusters composed of isolates from these two families. Those typing methods typically group isolates from these families into large clusters of non-linked isolates with identical fingerprints. Next-generation whole-genome sequencing (WGS) provides a new tool for molecular epidemiology that can resolve clusters of isolates with identical spoligotyping and MIRU-VNTR fingerprints.

Methods

We performed WGS and SNP analysis and evaluated epidemiological data to investigate 19 apparent TB transmission clusters in Hawaii from 2003 to 2017 in order to assess WGS’ ability to resolve putative Mtb clusters from the Beijing and Manila families. This project additionally investigated MIRU-VNTR allele prevalence to determine why standard Mtb fingerprinting fails to usefully distinguish actual transmission clusters from these two Mtb families.

Results

WGS excluded transmission events in seven of these putative clusters, confirmed transmission in eight, and identified both transmission-linked and non-linked isolates in four. For epidemiologically identified clusters, while the sensitivity of MIRU-VNTR fingerprinting for identifying actual transmission clusters was found to be 100%, its specificity was only 28.6% relative to WGS. We identified that the Beijing and Manila families’ significantly lower Shannon evenness of MIRU-VNTR allele distributions than lineage 4 was the cause of standard fingerprinting’s poor performance when identifying transmission in Beijing and Manila family clusters.

Conclusions

This study demonstrated that WGS is necessary for epidemiological investigation of TB in Hawaii and the Pacific.

Electronic supplementary material

The online version of this article (10.1186/s12879-018-3502-1) contains supplementary material, which is available to authorized users.

Genetic diversity of Mycobacterium tuberculosis isolates causing pulmonary and extrapulmonary tuberculosis in the capital of Iran

OBJECTIVES

Evaluation of the genetic diversity of Mycobacterium tuberculosis (M.tb) and determining if the association between a specific genotype and the site of infection is crucial. Accordingly, the current study aimed at comparing predominant M.tb genotypes in pulmonary (PTB) and extrapulmonary tuberculosis (EPTB) isolates circulating in the capital of Iran.

METHODS

The genetic diversity of culture-confirmed PTB and EPTB isolates were evaluated by Spoligotyping and MIRU-VNTR (mycobacterial interspersed repetitive-unit-variable-number tandem-repeat) typing methods. Genotyping data were analyzed with SITVIT, MIRU-VNTRplus, and TBminer databases. To assess adjusted associations, chi-square/the Fisher exact test and multiple logistic regression model were applied.

RESULTS

URAL2 (NEW-1) (28/88; 31.8%) and CAS1-DELHI (25/84; 29.8%) genotypes were predominant in EPTB and PTB strains, respectively. Based on MIRU-VNTR typing, 158 different MIRU-VNTR patterns were identified. Clustering rate and minimum estimate of the proportion of TB caused by recent transmission was 4.1% and 8.1%, respectively.

CONCLUSIONS

The current study provided new insight into circulating genotypes of M.tb in PTB and EPTB patients in Tehran, Iran. This low percentage of TB transmission rate, demonstrated that mode of TB transmission was mainly associated with reactivation of latent TB rather than recently transmitted infection in this region. There was no significant difference in the association between the genotypes of M.tb strains and the site of the disease.

Genetic diversity of drug resistant Mycobacterium Tuberculosis in local area of Southwest China: a retrospective study

BACKGROUND

By 2014 although tuberculosis (TB) incidence had fallen by an average of 1.5% per year since 2000 and was 18% lower than the level of in 2000, 1.5 million people died for TB in that year. One of reason was that drug resistant Mycobacterium tuberculosis (DRTB) spread. This study aims to determine drug resistant characteristics and genotype of DRTB that isolated from patients in a tuberculosis referral hospital of southwest China.

METHODS

Five hundred thirty-eight drug resistant tuberculosis samples were collected from July 2013 to March 2015. All the isolates were identified by genomic deletions in region of difference 105 (RD105) and genotyped by mycobacterial interspersed repetitive unit-variable number tandem repeat typing (MIRU-VNTR). Polymorphism and cluster analysis of each locus was carried out using Bionumerics Version 3.0 and phyloviz software.

RESULTS

Five hundred thirty-eight TB strains included 503 Mycobacterium tuberculosis (MTB) isolates and 35 non Mycobacterium tuberculosis (NMTB) isolates. Of 503 isolates Beijing family type was 447 (88.9%, 447/503) and non-Beijing family type was 56 (11.1%, 56/503). Five hundred three DRTB isolates were divided into 345 genotypes, of which 265 isolates were single genotype and the remaining 238 strains were classified into 80 clusters with cluster rate of 47.3% and cluster ratio of 31.4%. Sixty-nine clusters belonged to Beijing family with cluster rate was 48.3% and clustering ratio was 32.9%. The non - Beijing family had 11 clusters with a cluster rate of 39.3% and the clustering ratio of 19.6%. Beijing genotype had a significant correlation with the age (P < 0.05), the retreatment patients (P < 0.05) and the city of Chongqing (P < 0.05), not with gender (P > 0.05). In the 9 Beijing genotype clusters each cluster contained some patients who lived in the same region.

CONCLUSIONS

Beijing genotype was the predominant in the patients with DRTB in our hospital. In Chongqing retreatment patients with Beijing genotype MTB may be patient with DRTB. Drug resistance test (DST), regular medication and strict follow-up are very important for patients with Beijing genotype MTB. In Chongqing control and treatment of DRTB should be paid attention. Transmission and relations of patients with DRTB need to be further research.

Mixed Mycobacterium tuberculosis-Strain Infections Are Associated With Poor Treatment Outcomes Among Patients With Newly Diagnosed Tuberculosis, Independent of Pretreatment Heteroresistance

Background

Heteroresistant Mycobacterium tuberculosis infections (defined as concomitant infection with drug-resistant and drug-susceptible strains) may explain the higher risk of poor tuberculosis treatment outcomes observed among patients with mixed-strain M. tuberculosis infections. We investigated the clinical effect of mixed-strain infections while controlling for pretreatment heteroresistance in a population-based sample of patients with tuberculosis starting first-line tuberculosis therapy in Botswana.

Methods

We performed 24-locus mycobacterial interspersed repetitive unit-variable number tandem-repeat analysis and targeted deep sequencing on baseline primary cultured isolates to detect mixed infections and heteroresistance, respectively. Drug-sensitive, micro-heteroresistant, macro-heteroresistant, and fixed-resistant infections were defined as infections in which the frequency of resistance was <0.1%, 0.1%-4%, 5%-94%, and ≥95%, respectively, in resistance-conferring domains of the inhA promoter, the katG gene, and the rpoB gene.

Results

Of the 260 patients with tuberculosis included in the study, 25 (9.6%) had mixed infections and 30 (11.5%) had poor treatment outcomes. Micro-heteroresistance, macro-heteroresistance, and fixed resistance were found among 11 (4.2%), 2 (0.8%), and 11 (4.2%), respectively, for isoniazid and 21 (8.1%), 0 (0%), and 10 (3.8%), respectively, for rifampicin. In multivariable analysis, mixed infections but not heteroresistant infections independently predicted poor treatment outcomes.

Conclusions

Among patients starting first-line tuberculosis therapy in Botswana, mixed infections were associated with poor tuberculosis treatment outcomes, independent of heteroresistance.

The relationship between transmission time and clustering methods in Mycobacterium tuberculosis epidemiology

BACKGROUND

Tracking recent transmission is a vital part of controlling widespread pathogens such as Mycobacterium tuberculosis. Multiple methods with specific performance characteristics exist for detecting recent transmission chains, usually by clustering strains based on genotype similarities. With such a large variety of methods available, informed selection of an appropriate approach for determining transmissions within a given setting/time period is difficult.

METHODS

This study combines whole genome sequence (WGS) data derived from 324 isolates collected 2005-2010 in Kinshasa, Democratic Republic of Congo (DRC), a high endemic setting, with phylodynamics to unveil the timing of transmission events posited by a variety of standard genotyping methods. Clustering data based on Spoligotyping, 24-loci MIRU-VNTR typing, WGS based SNP (Single Nucleotide Polymorphism) and core genome multi locus sequence typing (cgMLST) typing were evaluated.

FINDINGS

Our results suggest that clusters based on Spoligotyping could encompass transmission events that occurred almost 200 years prior to sampling while 24-loci-MIRU-VNTR often represented three decades of transmission. Instead, WGS based genotyping applying low SNP or cgMLST allele thresholds allows for determination of recent transmission events, e.g. in timespans of up to 10 years for a 5 SNP/allele cut-off.

INTERPRETATION

With the rapid uptake of WGS methods in surveillance and outbreak tracking, the findings obtained in this study can guide the selection of appropriate clustering methods for uncovering relevant transmission chains within a given time-period. For high resolution cluster analyses, WGS-SNP and cgMLST based analyses have similar clustering/timing characteristics even for data obtained from a high incidence setting.

Mycobacterium simiae pulmonary disease in Iran: systematic review and meta-analysis

Mycobacterium simiae is one of the most common nontuberculous mycobacteria (NTM) microorganisms causing lung disease in many countries in the world. A reliable estimate of the extent of M. simiae pulmonary disease has not been well investigated in Iran. We systematically searched multiple databases to identify relative studies. Studies were excluded if they did not use the American Thoracic Society (ATS) and Infectious Diseases Society of America (IDSA) diagnostic criteria for NTM diseases. Data were extracted independently and in duplicate. We assessed pooled estimate by using a random model effect, and sources of heterogeneity were assessed by using Cochran's Q and the I 2 statistic. The potential for publication bias was explored by using Begg's and Egger's tests. All analyses were conducted with Stata 14.0 (StataCorp, College Station, TX, USA). Of 172 articles identified, seven met the inclusion criteria. Of 355 patients who were culture positive for NTM, 82 had M. simiae pulmonary disease according to the ATS/IDSA diagnostic criteria. The pooled frequency of M. simiae pulmonary disease among patients with NTM was 25.0% (95% confidence interval, 16.8-33.2). No evidence of publication bias was observed among the included studies (p >0.05 for Begg's and Egger's tests). Clinical isolates of M. simiae are increasingly being recognized as a cause of pulmonary disease in Iran and need further attention by health authorities.

Prevalence and Characterization of Heterogeneous Variable-Number Tandem-Repeat Clusters Comprising Drug-Susceptible and/or Variable Resistant Mycobacterium tuberculosis Complex Isolates in the Netherlands from 2004 to 2016

The variable-number tandem-repeat (VNTR) typing method is used to study tuberculosis (TB) transmission. Clustering of Mycobacterium tuberculosis isolates with identical VNTR patterns is assumed to reflect recent transmission. Hence, clusters are thought to be homogeneous regarding antibiotic resistance. In practice, however, heterogeneous clusters are also identified. This study investigates the prevalence and characteristics of heterogeneous VNTR clusters and assesses whether isolates in these clusters remain clustered when subjected to whole-genome sequencing (WGS). In the period from 2004 to 2016, 9,072 isolates were included. Demographic and epidemiological linkage data were obtained from the Netherlands Tuberculosis Register. VNTR clusters were defined as homogeneous when isolates shared identical resistance profiles or as heterogeneous if both susceptible and (variable) resistant isolates were found. Multivariate logistic regression analysis was performed to identify factors associated with heterogeneous clustering. Isolates from 2016 were subjected to WGS, and a genetic distance of 12 single nucleotide polymorphisms (SNPs) was used as the cutoff for WGS clustering. In total, 4,661/9,072 (51%) isolates were clustered into 985 different VNTR clusters, of which 217 (22%) were heterogeneous. Patient characteristics associated with heterogeneous clustering were non-Dutch ethnicity (odds ratio [OR], 1.46 [95% confidence interval {CI}, 1.22 to 1.75]), asylum seeker (OR, 1.51 [95% CI, 1.24 to 1.85]), extrapulmonary TB (OR, 1.26 [95% CI, 1.09 to 1.46]), previous TB diagnosis (OR, 1.38 [95% CI, 1.04 to 1.82]), and not being a contact of a TB patient (OR, 1.35 [95% CI, 1.08 to 1.69]). With WGS, 34% of heterogeneous and 78% of homogeneous isolates from 2016 remained clustered. Heterogeneous VNTR clusters are common but seem to be explained by a substantial degree of false clustering by VNTR typing compared to WGS.

Set-up and validation of mycobacterial interspersed repetitive unit-variable number of tandem repeat (MIRU-VNTR) analysis of Mycobacterium tuberculosis using BioNumerics software

The objective was to describe and validate a new and alternative software procedure for 24-locus mycobacterial interspersed repetitive unit-variable number-tandem repeat (MIRU-VNTR) typing of Mycobacterium tuberculosis (Mtb) based on the multipurpose BioNumerics software. DNA from randomly selected isolates of Mtb from two European laboratories, including external control samples for MIRU-VNTR typing, were analysed. Samples were genotyped using the commercial 24-locus VNTR typing kit from GenoScreen. The PCR amplified fragments were separated by capillary electrophoresis. For the subsequent analyses, the currently used software GeneMapper was compared with BioNumerics. The endpoint was the level of concordance when comparing genotyping results obtained from BioNumerics with results obtained from GeneMapper and the ECDC proficiency study reference results. Also, the number of necessary manual standard size corrections and allele assignments in the two different software methods were compared. In total, 272 DNA samples, including the ECDC proficiency panel, were analysed. For all samples, there were 100% concordance of results. For a randomly selected set of 96 samples the numbers of manual corrections needed for size standards were 199 with GeneMapper versus zero for BioNumerics. The numbers of manual corrections for allele assignments were 122 with GeneMapper versus 16 with BioNumerics. In conclusion, we have validated the multipurpose software BioNumerics for standard 24-locus MIRU-VNTR typing and the software shows promising benefits in terms of simplification and minimization of hand-on time.

Mutations of Mycobacterium tuberculosis induced by anti-tuberculosis treatment result in metabolism changes and elevation of ethambutol resistance

Selective pressure from antibiotic use is one of the most important risk factors associated with the development of drug resistance in Mycobacterium tuberculosis (MTB). However, the mechanisms underlying drug resistance at the molecular level remain partly unclear. Therefore, the purpose of this study was to investigate the potential functional effect of novel mutations arising from anti-tuberculosis treatment. We analyzed two multidrug-resistant TB (MDR-TB) isolates from the same patient; one collected before and one almost a year after commencing MDR-TB treatment. The post-treatment isolate exhibited elevated ethambutol resistance. We sequenced the whole genomes of the two clinical isolates and detected six novel polymorphisms affecting the genes Rv1026, nc0021, Rv2155c, Rv2437, and Rv3696c, and the intergenic region between Rv2764c and Rv2765. Metabolomics approach was used to reveal the effect of the found variation on the metabolic pathways of MTB. Partial least squares-discriminant analysis showed a clear differentiation between the two isolates, involving a total of 175 metabolites. Pathway analysis showed that these metabolites are mainly involved in amino sugar and nucleotide sugar metabolism, β-alanine metabolism, sulfur metabolism, and galactose metabolism. The increased ethambutol resistance exhibited by the post-treatment MDR-TB strain could speculatively be linked to the identified genetic variations, which affected the synthesis of a number of metabolites associated with sources of carbon and energy. This may have been the main factor underlying the increased ethambutol resistance of this isolate.

Role of Molecular Epidemiology on Tuberculosis Control in the Middle East Countries: a Systematic Review and Meta-Analysis

BACKGROUND

Tuberculosis (TB) is a public health problem in developing countries and yet the numbers of people with the disease are abundant. Early detection of transmission sources and effective treatment of the cases is essential to control the disease which will be possible by application of molecular epidemiology approaches. Studies conducted based on Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats (MIRU-VNTR) method in Muslim Middle East countries were evaluated to determine their role in TB control.

MATERIALS AND METHODS

All studies from January 2005 to April 2016 were systematically reviewed in four electronic databases and finally 16 articles were found eligible to be included in this study. The mean clustering rate was determined as 44% and the recent transmission rate was 12.3 to 78.8% with average of 33%.

RESULTS

The results showed that both reactivation and recent transmission were important in developing new cases of TB in Middle East countries; but, reactivation plays a more critical role.

CONCLUSION

Regarding to ongoing war and immigration in the region along with the increasing of drug-resistant TB, in the case of improper supervision in the future, the disease, especially drug- resistant TB, will be problematic.

Longitudinal genotyping surveillance of Mycobacterium tuberculosis in an area with high tuberculosis incidence shows high transmission rate of the modern Beijing subfamily in Japan

Tuberculosis (TB) is a severe and wide-spread infectious disease worldwide. The modern Beijing subfamily, one lineage of M. tuberculosis, reportedly has high pathogenicity and transmissibility. This study used a molecular epidemiological approach to investigate the transmissibility of the modern Beijing subfamily in the Airin area of Osaka City, Japan. During 2006-2016, we collected 596 M. tuberculosis clinical isolates in the Airin area, Osaka city, Japan. We analyzed the 24-locus variable number of tandem repeats typing optimized for the Beijing family of isolates, M. tuberculosis lineage, and patient epidemiological data. The proportion of the modern Beijing subfamily was significantly higher not only than previously obtained data for the Airin area: it was also higher than the nationwide in Japan. The rate of recent clusters, defined as a variable number of tandem repeats profile identified within two years, of the modern Beijing subfamily was significantly higher than that the rate of recent clusters of the ancient Beijing subfamily. Results suggest that TB control measures formulated with attention to the modern Beijing subfamily might be an important benchmark to understanding recent TB transmission in the area.

Tuberculosis caused by Mycobacterium bovis infection in a captive-bred American bullfrog (Lithobates catesbeiana)

BACKGROUND

Tuberculosis is widely known as a progressive disease that affects endothermic animals, leading to death and/or economical losses, while mycobacterial infections in amphibians are commonly due to nontuberculous mycobacteria. To the authors' knowledge, this report describes the first case of bovine tuberculosis in a poikilothermic animal.

CASE PRESENTATION

An adult female captive American bullfrog (Lithobates catesbeianus Shaw, 1802) died in a Brazilian aquarium. Multiple granulomas with acid-fast bacilli were observed in several organs. Identification of Mycobacterium bovis was accomplished by culture and PCR methods. The other animals from the same enclosure were euthanized, but no evidence of mycobacterial infection was observed.

CONCLUSIONS

The American bullfrog was introduced in several countries around the world as an alternative husbandry, and its production is purposed for zoological and aquarium collections, biomedical research, education, human consumption and pet market. The present report warns about an episode of bovine tuberculosis in an amphibian, therefore further studies are necessary to define this frog species' role in the epidemiology of M. bovis.

Molecular epidemiology and drug resistance patterns of Mycobacterium tuberculosis complex isolates from university students and the local community in Eastern Ethiopia

BACKGROUND

Previous studies suggest the burden of pulmonary tuberculosis (PTB) in Ethiopia may be greater in university students relative to the overall population. However, little is known about the transmission dynamics of PTB among students and members of the communities surrounding university campuses in Eastern Ethiopia.

METHODS

A cross sectional study was conducted in Eastern Ethiopia among prevalent culture-confirmed PTB cases from university students (n = 36) and community members diagnosed at one of four hospitals (n = 152) serving the surrounding area. Drug susceptibility testing (DST) was performed on Mycobacterium tuberculosis complex (MTBC) isolates using BD Bactec MGIT 960 and molecular genotyping was performed using spoligotyping and 24-loci MIRU-VNTR. MTBC strains with Identical genotyping patterns were assigned to molecular clusters as surrogate marker for recent transmission and further contact tracing was initiated among clustered patients.

RESULTS

Among all study participants, four MTBC lineages and 11 sub-lineages were identified, with Ethiopia_3 (Euro-American lineage) being most common sub-lineage (29.4%) in both cohorts and associated with strain clustering (P = 0.016). We further identified 13 (8.1%) strains phylogenetically closely related to Ethiopia_3 but with a distinct Spoligotyping pattern and designated as Ethiopia_4. The clustering rate of MTBC strains was 52.9% for university students and 66.7% for community members with a Recent Transmission Index (RTI) of 17.6% and 48.4%, respectively. Female gender, urban residence, and new TB cases were significantly associated with strain clustering (P<0.05). Forty-eight (30%) of the study participants were resistant to one or more first line anti TB drugs, three patients were classified as multidrug resistant (MDR).

CONCLUSION

We found evidence for recent transmission of PTB among Ethiopian university students and the local community in Eastern Ethiopia, mainly linked to strains classified as Ethiopia_3 sub lineage. Drug resistance didn't have a major impact on recent transmission but comprehensive molecular surveillance in combination with drug resistance profiling of MTBC strains is desirable to better characterize TB transmission dynamics in high risk congregate living environments such as university campuses and guide regional TB control programs.

Role of MIRU-VNTR and spoligotyping in assessing the genetic diversity of Mycobacterium tuberculosis in Henan Province, China

Background

Tuberculosis remains a serious threat to human health as an infectious disease in China. Henan, a most populated province in China, has a high incidence of tuberculosis (TB). Though the genetic diversity of Mycobacterium tuberculosis (MTB) has been investigated in many regions, there have been only a few studies on the molecular characteristics and drug resistance phenotypes in Henan. This is the first study on the genetic profile of MTB from Henan.

Methods

A total of 668 MTB isolates from various areas were genotyped with spoligotyping and 26-locus MIRU-VNTR (classical 24-locus MIRU-VNTR and 2 other loci). The association between TB spoligotype signatures and drug-resistant profiles was analysed.

Results

Our data revealed that MTB isolates circulating in Henan had a high degree of genetic variation. The Beijing family was the most predominant genotype (83.53%,n = 558), and the typical Beijing type(ST1) was the major sublineage (81.73%,n = 546). In total,668 isolates were divided into 567 different types, forming 38 clusters (2–15 isolates per cluster), and 529 unique types by 26-locus MIRU-VNTR analysis. There was no correlation between the Beijing family and gender, age at diagnosis or treatment history, whereas the Beijing family was significantly associated with all four first-line drug resistance and multidrug-resistant phenotypes. For these samples, 15 of 26 MIRU-VNTR loci had high or moderate discriminatory power according to the Hunter-Gaston discriminatory index. A combination of the 10 most polymorphic loci had similar discriminatory power as the 26-locus set.

Conclusion

The Beijing genotype is the most prevalent family. Ten-locus MIRU-VNTR in combination with spoligotyping can efficiently classify the molecular type of MTB in Henan Province.

Electronic supplementary material

The online version of this article (10.1186/s12879-018-3351-y) contains supplementary material, which is available to authorized users.

Mycobacterium bovis Persistence in Two Different Captive Wild Animal Populations in Germany: a Longitudinal Molecular Epidemiological Study Revealing Pathogen Transmission by Whole-Genome Sequencing

Bovine tuberculosis (bTB) caused by Mycobacterium bovis is a transmissible disease notifiable to the World Organization for Animal Health and to the European Union, with ongoing efforts of surveillance and eradication in every EU member state. In Germany, a country which has been declared officially free from bovine tuberculosis since 1997 by the EU, M. bovis infections still occur sporadically in cattle and other mammals, including humans. Here, the transmission routes of a bTB outbreak in a wildlife park in Germany affecting different cervid species, bison, lynx, and pot-bellied pigs were followed by employing whole-genome sequencing (WGS) combined with spoligotyping and mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) typing. One single M. bovis strain persisted from 2002 to 2015, and transmission between the park and a distantly located captive cervid farm was verified. The spoligotyping patterns remained identical, while MIRU-VNTR typing of 24 loci of the standardized panel and locus 2163a as an additional locus revealed one change at locus 2165 in a strain from a fallow deer and one at locus 2461 in isolates from red deer over the whole time period. WGS analysis confirmed the close relatedness of the isolates, with a maximum of 12 single nucleotide polymorphisms (SNPs) detected between any two sequenced isolates. In conclusion, our data confirm a longitudinal outbreak of M. bovis in a German wildlife park and provide the first insights into the dynamics of different genotyping markers in M. bovis.

Relapse Versus Reinfection of Recurrent Tuberculosis Patients in a National Tuberculosis Specialized Hospital in Beijing, China

Tuberculosis (TB) recurrence can result from either relapse of an original infection or exogenous reinfection with a new strain of Mycobacterium tuberculosis (MTB). The aim of this study was to assess the roles of relapse and reinfection among recurrent TB cases characterized by a high prevalence rate of drug-resistant TB within a hospital setting. After 58 paired recurrent TB cases were genotyped to distinguish relapse from reinfection, 37 (63.8%) were demonstrated to be relapse cases, while the remaining 21 were classified as reinfection cases. Statistical analysis revealed that male gender was a risk factor for TB reinfection, odds ratios and 95% confidence interval (OR [95% CI]: 4.188[1.012–17.392], P = 0.049). Of MTB isolates obtained from the 37 relapse cases, 11 exhibited conversion from susceptible to resistance to at least one antibiotic, with the most frequent emergence of drug resistance observed to be levofloxacin. For reinfection cases, reemergence of rifampicin-resistant isolates harboring double gene mutations, of codon 531 of rpoB and codon 306 of embB, were observed. In conclusion, our data demonstrate that relapse is a major mechanism leading to TB recurrence in Beijing Chest Hospital, a national hospital specialized in TB treatment. Moreover, male patients are at higher risk for reinfection. The extremely high rate of multidrug-resistant tuberculosis (MDR-TB) among reinfection cases reflects more successful transmission of MDR-TB strains versus non-resistant strains overall.

Harmonized Genome Wide Typing of Tubercle Bacilli Using a Web-Based Gene-By-Gene Nomenclature System

Background

Global tuberculosis (TB) control is challenged by uncontrolled transmission of Mycobacterium tuberculosis complex (Mtbc) strains, esp. of multidrug (MDR) or extensively resistant (XDR) variants. Precise analysis of transmission networks is the basis to trace outbreak M/XDR clones and improve TB control. However, classical genotyping tools lack discriminatory power due to the high similarity of strains of particular successful lineages, e.g. Beijing or outbreak strains. This can be overcome by whole genome sequencing (WGS) approaches, but these are not yet standardized to facilitate larger investigations encompassing different laboratories or outbreak tracing across borders.

Methods

We established and improved a whole genome gene-by-gene multi locus sequence typing approach encompassing a stable set of core genome genes (cgMLST) and linked it to a web-based nomenclature server (cgMLST.org) facilitating assignment and storage of allele numbers.

Findings

We evaluated and refined a previously suggested cgMLST schema by using a reference strain set (n = 251) reflecting the global diversity of the Mtbc. A set of 2891 genes showed excellent performance with at least 97% of the genes reliably identified in strains of all Mtbc lineages and in discriminating outbreak strains. cgMLST allele numbers were automatically retrieved from and stored at cgMLST.org.

Interpretation

The refined cgMLST schema provides high resolution genome-based typing of clinical strains of all Mtbc lineages. Combined with a web-based nomenclature server, it facilitates rapid, high-resolution, and harmonized tracing of clinical Mtbc strains needed for prospective local and global surveillance.

A Major Mycobacterium tuberculosis outbreak caused by one specific genotype in a low-incidence country: Exploring gene profile virulence explanations

Denmark, a tuberculosis low burden country, still experiences significant active Mycobacterium tuberculosis (Mtb) transmission, especially with one specific genotype named Cluster 2/1112-15 (C2), the most prevalent lineage in Scandinavia. In addition to environmental factors, antibiotic resistance, and human genetics, there is increasing evidence that Mtb strain variation plays a role for the outcome of infection and disease. In this study, we explore the reasons for the success of the C2 genotype by analysing strain specific polymorphisms identified through whole genome sequencing of all C2 isolates identified in Denmark between 1992 and 2014 (n = 952), and the demographic distribution of C2. Of 234 non-synonymous (NS) monomorphic SNPs found in C2 in comparison with Mtb reference strain H37Rv, 23 were in genes previously reported to be involved in Mtb virulence. Of these 23 SNPs, three were specific for C2 including a NS mutation in a gene associated with hyper-virulence. We show that the genotype is readily transmitted to different ethnicities and is also found outside Denmark. Our data suggest that strain specific virulence factor variations are important for the success of the C2 genotype. These factors, likely in combination with poor TB control, seem to be the main drivers of C2 success.

Genetic diversity of Mycobacterium tuberculosis complex isolates circulating in an area with high tuberculosis incidence: Using 24-locus MIRU-VNTR method

We aimed to determine the genetic diversity, phylogenetic relationship and transmission dynamics of Mycobacterium tuberculosis complex (MTBC) genotypes in an area with high tuberculosis (TB) incidence. A set of 164 MTBC isolates from new TB patients of Golestan province, Iran, were subjected to genotyping using the standard 24-locus MIRU-VNTR method. Recent TB transmission was evaluated and phylogenetic relationships were analysed by minimum spanning tree and cluster-graph methods. Among the 164 isolates, 132 distinct patterns were detected. The 48 clustered isolates (29.3%) were distributed into 16 clusters ranging in size from 2 to 12 isolates. The most frequent genotype was Central Asian Strain/Delhi (CAS/Delhi) (n = 67, 40.8%), followed by NEW-1 (n = 53, 32.3%) and Beijing (n = 19, 11.6%) genotypes. Thirty five (72.9%) of NEW-1 isolates were recovered from immigrant patients and 84.2% (n = 16) of Beijing genotypes recovered from native cases. Statistically significant association was found between clustering and smoking (p = 0.047), drug addiction (p = 0.01) and prison history (p = 0.003). The estimated proportion of recent transmission was 19.5%. Presence of highly diverse MTBC isolates circulating in this province without a dominant genotype might be a consequence of importation of various genotypes in this area.

Investigating spillover of multidrug-resistant tuberculosis from a prison: a spatial and molecular epidemiological analysis

BACKGROUND

Congregate settings may serve as institutional amplifiers of tuberculosis (TB) and multidrug-resistant tuberculosis (MDR-TB). We analyze spatial, epidemiological, and pathogen genetic data prospectively collected from neighborhoods surrounding a prison in Lima, Peru, where inmates experience a high risk of MDR-TB, to investigate the risk of spillover into the surrounding community.

METHODS

Using hierarchical Bayesian statistical modeling, we address three questions regarding the MDR-TB risk: (i) Does the excess risk observed among prisoners also extend outside the prison? (ii) If so, what is the magnitude, shape, and spatial range of this spillover effect? (iii) Is there evidence of additional transmission across the region?

RESULTS

The region of spillover risk extends for 5.47 km outside of the prison (95% credible interval: 1.38, 9.63 km). Within this spillover region, we find that nine of the 467 non-inmate patients (35 with MDR-TB) have MDR-TB strains that are genetic matches to strains collected from current inmates with MDR-TB, compared to seven out of 1080 patients (89 with MDR-TB) outside the spillover region (p values: 0.022 and 0.008). We also identify eight spatially aggregated genetic clusters of MDR-TB, four within the spillover region, consistent with local transmission among individuals living close to the prison.

CONCLUSIONS

We demonstrate a clear prison spillover effect in this population, which suggests that interventions in the prison may have benefits that extend to the surrounding community.

Unexpected Genomic and Phenotypic Diversity of Mycobacterium africanum Lineage 5 Affects Drug Resistance, Protein Secretion, and Immunogenicity

Mycobacterium africanum consists of Lineages L5 and L6 of the Mycobacterium tuberculosis complex (MTBC) and causes human tuberculosis in specific regions of Western Africa, but is generally not transmitted in other parts of the world. Since M. africanum is evolutionarily closely placed between the globally dispersed Mycobacterium tuberculosis and animal-adapted MTBC-members, these lineages provide valuable insight into M. tuberculosis evolution. Here, we have collected 15 M. africanum L5 strains isolated in France over 4 decades. Illumina sequencing and phylogenomic analysis revealed a previously underappreciated diversity within L5, which consists of distinct sublineages. L5 strains caused relatively high levels of extrapulmonary tuberculosis and included multi- and extensively drug-resistant isolates, especially in the newly defined sublineage L5.2. The specific L5 sublineages also exhibit distinct phenotypic characteristics related to in vitro growth, protein secretion and in vivo immunogenicity. In particular, we identified a PE_PGRS and PPE-MPTR secretion defect specific for sublineage L5.2, which was independent of PPE38. Furthermore, L5 isolates were able to efficiently secrete and induce immune responses against ESX-1 substrates contrary to previous predictions. These phenotypes of Type VII protein secretion and immunogenicity provide valuable information to better link genome sequences to phenotypic traits and thereby understand the evolution of the MTBC.

A Quantitative Evaluation of MIRU-VNTR Typing Against Whole-Genome Sequencing for Identifying Mycobacterium tuberculosis Transmission: A Prospective Observational Cohort Study

Background

Mycobacterial Interspersed Repetitive Unit-Variable Number Tandem Repeat (MIRU-VNTR) typing is widely used in high-income countries to determine Mycobacterium tuberculosis relatedness. Whole-genome sequencing (WGS) is known to deliver greater specificity, but no quantitative prospective comparison has yet been undertaken.

Methods

We studied isolates from the English Midlands, sampled consecutively between 1 January 2012 and 31 December 2015. In addition to routinely performed MIRU-VNTR typing, DNA was extracted from liquid cultures and sequenced using Illumina technology. Demographic and epidemiological data for the relevant patients were extracted from the Enhanced Tuberculosis Surveillance system run by Public Health England. Closely related samples, defined using a threshold of five single nucleotide variants (SNVs), were compared to samples with identical MIRU-VNTR profiles, to samples from individuals with shared epidemiological risk factors, and to those with both characteristics.

Findings

1999 patients were identified for whom at least one M. tuberculosis isolate had been MIRU-VNTR typed and sequenced. Comparing epidemiological risk factors with close genetic relatedness, only co-residence had a positive predictive value of over 5%. Excluding co-resident individuals, 18.6% of patients with identical MIRU-VNTR profiles were within 5 SNVs. Where patients also shared social risk factors and ethnic group, this rose to 48%. Only 8% of MIRU-VNTR linked pairs in lineage 1 were within 5 SNV, compared to 31% in lineage 4.

Interpretation

In the setting studied, this molecular epidemiological study shows MIRU-VNTR typing and epidemiological risk factors are poorly predictive of close genomic relatedness, assessed by SNV. MIRU-VNTR performance varies markedly by lineage.

Funding

Public Health England, Health Innovation Challenge Fund, NIHR Health Protection Research Unit Oxford, NIHR Oxford Biomedical Research Centre.

•

M. tuberculosis relatedness was measured by MIRU-VNTR typing and single nucleotide variants (SNV).

•

In consecutive isolates in England, only 20% of samples with identical MIRU-VNTR profiles were within 5 SNVs.

•

MIRU-VNTR is more predictive of close relatedness in lineage 4 (Euro-American) lineages than in other lineages.

Geospatial distribution of Mycobacterium tuberculosis genotypes in Africa

Objective

To investigate the distribution of Mycobacterium tuberculosis genotypes across Africa.

Methods

The SITVIT2 global repository and PUBMED were searched for spoligotype and published genotype data respectively, of M. tuberculosis from Africa. M. tuberculosis lineages in Africa were described and compared across regions and with those from 7 European and 6 South-Asian countries. Further analysis of the major lineages and sub-lineages using Principal Component analysis (PCA) and hierarchical cluster analysis were done to describe clustering by geographical regions. Evolutionary relationships were assessed using phylogenetic tree analysis.

Results

A total of 14727 isolates from 35 African countries were included in the analysis and of these 13607 were assigned to one of 10 major lineages, whilst 1120 were unknown. There were differences in geographical distribution of major lineages and their sub-lineages with regional clustering. Southern African countries were grouped based on high prevalence of LAM11-ZWE strains; strains which have an origin in Portugal. The grouping of North African countries was due to the high percentage of LAM9 strains, which have an origin in the Eastern Mediterranean region. East African countries were grouped based on Central Asian (CAS) and East-African Indian (EAI) strain lineage possibly reflecting historic sea trade with Asia, while West African Countries were grouped based on Cameroon lineage of unknown origin. A high percentage of the Haarlem lineage isolates were observed in the Central African Republic, Guinea, Gambia and Tunisia, however, a mixed distribution prevented close clustering.

Conclusions

This study highlighted that the TB epidemic in Africa is driven by regional epidemics characterized by genetically distinct lineages of M. tuberculosis. M. tuberculosis in these regions may have been introduced from either Europe or Asia and has spread through pastoralism, mining and war. The vast array of genotypes and their associated phenotypes should be considered when designing future vaccines, diagnostics and anti-TB drugs.

Reduced transmission of Mycobacterium africanum compared to Mycobacterium tuberculosis in urban West Africa

OBJECTIVE

Understanding transmission dynamics is useful for tuberculosis (TB) control. A population-based molecular epidemiological study was conducted to determine TB transmission in Ghana.

METHODS

Mycobacterium tuberculosis complex (MTBC) isolates obtained from prospectively sampled pulmonary TB patients between July 2012 and December 2015 were characterized using spoligotyping and standard 15-locus mycobacterial interspersed repetitive unit variable number tandem repeat (MIRU-VNTR) typing for transmission studies.

RESULTS

Out of 2309 MTBC isolates, 1082 (46.9%) unique cases were identified, with 1227 (53.1%) isolates belonging to one of 276 clusters. The recent TB transmission rate was estimated to be 41.2%. Whereas TB strains of lineage 4 belonging to M. tuberculosis showed a high recent transmission rate (44.9%), reduced recent transmission rates were found for lineages of Mycobacterium africanum (lineage 5, 31.8%; lineage 6, 24.7%).

CONCLUSIONS

The study findings indicate high recent TB transmission, suggesting the occurrence of unsuspected outbreaks in Ghana. The observed reduced transmission rate of M. africanum suggests other factor(s) (host/environmental) may be responsible for its continuous presence in West Africa.

Molecular characterization and drug susceptibility profile of Mycobacterium tuberculosis isolates from Northeast Bangladesh

Tuberculosis (TB) remains a major public health problem worldwide including in Bangladesh. Molecular epidemiological tools provide genotyping profiles of Mycobacterium tuberculosis (M. tuberculosis) strains that can give insight into the transmission of TB in a specific region. The objective of the study was to identify the genetic diversity and drug susceptibility profile of M. tuberculosis strains circulating in the northeast Bangladesh. A total of 244 smear-positive sputum specimens were collected from two referral hospitals in Mymensingh and Netrakona districts. The isolated strains were genotyped by deletion analysis, spoligotyping, and MIRU-VNTR typing. We also analyzed the distributions of drug susceptibility pattern and demographic data among different genotypes. All isolates were identified as M. tuberculosis and among them 167 strains (68.44%) were 'ancestral' and the remaining 77 (31.56%) were 'modern' type. Spoligotyping analysis yielded 119 distinct patterns, among them, 86 isolates had unique patterns and the remaining 158 were grouped into 33 distinct clusters containing 2 to 18 isolates. The predominant spoligotypes belong to the EAI lineage strains, comprising 66 (27.04%) isolates followed by Beijing (7.38%), T1 (6.15%), CAS1-Delhi (5.33), LAM9 (3.28%), MANU-2 and X2. MIRU-VNTR analysis revealed 167 isolates (68%) had unique patterns, whereas 77 (32%) were grouped into 26 clusters and the rate of recent transmission was 20.9%, suggesting that the majority of TB cases in this region are caused by the reactivation of previous TB infections rather than recent transmission. About 136 (55.7%) isolates were sensitive to four anti-TB drugs, 69 (28.3%) were resistant to one or more (except rifampicin and isoniazid combination) drugs and 39 (15.9%) were MDR. In conclusion, our study provides a first insight into molecular characterization and drug resistance profile of M. tuberculosis strains in northeast Bangladesh which will ultimately contribute to the national TB control program.

Characterization of pncA mutations in multi-drug and pyrazinamide resistant Mycobacterium tuberculosis isolates cultured from Queensland migrants and Papua New Guinea residents

Outbreak of drug resistant tuberculosis in the Western province, Papua New Guinea is a concern to Queensland, Australia due to migration. We performed pncA mutation analysis and genotyping of multi-drug/pyrazinamide (MDR/PZA) resistant isolates from 18 Queensland (Qld) migrants and 81 Papua New Guinea (PNG) residents, to compare with phenotypic evidence of PZA resistance and to evaluate the genotypes obtained from the two countries. Seven different mutations were seen from Qld isolates of which 2 have not been described previously. A cluster of mutations were found between amino acids L35 and S65. Amongst the PNG isolates, 10 mutations were identified, of which 6 were unique and have not been described previously. Majority of the mutations formed 2 clusters, between amino acids Q10 to A20 and W68 to W119. Mutations identified at nucleotide (nt) position 202 and 307 were found to be the most common types, occurring in 25% and 51% of the PNG isolates respectively. The majority of the mutations were seen in MDR/PZA resistant isolates. These mutations could be utilized for direct screening of PZA resistance from PNG patient samples. Genotypic analysis of the isolates showed strong clustering amongst the PNG isolates as opposed to Qld isolates. A diversity of mutations and genotypes were seen amongst the Qld migrant isolates. Majority of PNG isolates had one genotype with two distinct pncA mutation patterns (T202C and T307G) which highlight on-going transmission. pncA mutation analysis provided a satisfactory alternative to PZA culture DST with high positive predictive value and an improved result turnaround time.

MIRU-profiler: a rapid tool for determination of 24-loci MIRU-VNTR profiles from assembled genomes of Mycobacterium tuberculosis

Background

Tuberculosis (TB) resulted in an estimated 1.7 million deaths in the year 2016. The disease is caused by the members of Mycobacterium tuberculosis complex, which includes Mycobacterium tuberculosis, Mycobacterium bovis and other closely related TB causing organisms. In order to understand the epidemiological dynamics of TB, national TB control programs often conduct standardized genotyping at 24 Mycobacterial-Interspersed-Repetitive-Units (MIRU)-Variable-Number-of-Tandem-Repeats (VNTR) loci. With the advent of next generation sequencing technology, whole-genome sequencing (WGS) has been widely used for studying TB transmission. However, an open-source software that can connect WGS and MIRU-VNTR typing is currently unavailable, which hinders interlaboratory communication. In this manuscript, we introduce the MIRU-profiler program which could be used for prediction of MIRU-VNTR profile from WGS of M. tuberculosis.

Implementation

The MIRU-profiler is implemented in shell scripting language and depends on EMBOSS software. The in-silico workflow of MIRU-profiler is similar to those described in the laboratory manuals for genotyping M. tuberculosis. Given an input genome sequence, the MIRU-profiler computes alleles at the standard 24-loci based on in-silico PCR amplicon lengths. The final output is a tab-delimited text file detailing the 24-loci MIRU-VNTR pattern of the input sequence.

Validation

The MIRU-profiler was validated on four datasets: complete genomes from NCBI-GenBank (n = 11), complete genomes for locally isolated strains sequenced using PacBio (n = 4), complete genomes for BCG vaccine strains (n = 2) and draft genomes based on 250 bp paired-end Illumina reads (n = 106).

Results

The digital MIRU-VNTR results were identical to the experimental genotyping results for complete genomes of locally isolated strains, BCG vaccine strains and five out of 11 genomes from the NCBI-GenBank. For draft genomes based on short Illumina reads, 21 out of 24 loci were inferred with a high accuracy, while a number of inaccuracies were recorded for three specific loci (ETRA, QUB11b and QUB26). One of the unique features of the MIRU-profiler was its ability to process multiple genomes in a batch. This feature was tested on all complete M. tuberculosis genome (n = 157), for which results were successfully obtained in approximately 14 min.

Conclusion

The MIRU-profiler is a rapid tool for inference of digital MIRU-VNTR profile from the assembled genome sequences. The tool can accurately infer repeat numbers at the standard 24 or 21/24 MIRU-VNTR loci from the complete or draft genomes respectively. Thus, the tool is expected to bridge the communication gap between the laboratories using WGS and those using the conventional MIRU-VNTR typing.

Mycobacterium tuberculosis cording in alveolar macrophages of patients with pulmonary tuberculosis is likely associated with increased mycobacterial virulence

Mycobacterium tuberculosis (Mtb) is an infectious agent that causes tuberculosis (TB) in humans. A study of the volume of Mtb population and the detection of Mtb virulence in the lungs of patients with pulmonary TB are of great importance for understanding the infectious process and the outcome of the disease. We analyzed the functional state of Mtb and their number in alveolar macrophages obtained from the resected lungs of patients with TB in ex vivo culture and determined that the number of Mtb, referred mainly to the Beijing genotype family (A0 and B0/W148 clusters), were significantly different in cells between different patients. Only single Mtb were found in alveolar macrophages of some patients, while Mtb were actively replicated in colonies in alveolar macrophages of other patients, including cord morphology of Mtb growth (the indicator of Mtb virulence). Our data demonstrated association between the formation of Mtb cording in alveolar macrophages of patients and increased virulence of Mtb from the lungs of these patients in guinea pig TB model. The find of cording formation by replicating Mtb in human alveolar macrophages may be used for preliminary quick estimation of increased Mtb virulence in individual patients with pulmonary TB.

Insight into multidrug-resistant Beijing genotype Mycobacterium tuberculosis isolates in Myanmar

OBJECTIVES

Myanmar is a World Health Organization high tuberculosis (TB) burden country with a high multidrug-resistant (MDR)-TB burden. Of significance, a high prevalence of the Beijing genotype of Mycobacterium tuberculosis (MTB) among MDR-MTB has been reported previously. A detailed genetic characterization of TB clinical isolates was performed in order to explore whether there is an association between the prevalence of the Beijing MTB genotype and MDR-TB in Myanmar.

METHODS

A total of 265 MDR-MTB clinical isolates collected in 2010 and 2012 were subjected to spoligotyping, mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) analysis, single nucleotide polymorphism (SNP) typing, and drug resistance-associated gene sequencing, including rpoC to detect potential compensatory evolution.

RESULTS

Of the total MDR-MTB isolates, 79.2% (210/265) were of the Beijing genotype, the majority of which were the 'modern' subtype. Beijing genotype isolates were differentiated by 15-locus MIRU-VNTR and a high clustering rate (53.0%) was observed in the modern subtype. These MIRU-VNTR patterns were similar to Beijing genotype clones spreading across Russia and Central Asia. A high prevalence of katG Ser315Thr, and genetic evidence of extensive drug resistance (XDR) and pre-XDR and compensatory mutations in rpoC were observed among clustered isolates.

CONCLUSIONS

MDR-MTB strains of the Beijing genotype might be spreading in Myanmar and present a major challenge to TB control in this country.

Polyclonal Pulmonary Tuberculosis Infections and Risk for Multidrug Resistance, Lima, Peru

Because within-host Mycobacterium tuberculosis diversity complicates diagnosis and treatment of tuberculosis (TB), we measured diversity prevalence and associated factors among 3,098 pulmonary TB patients in Lima, Peru. The 161 patients with polyclonal infection were more likely than the 115 with clonal or the 2,822 with simple infections to have multidrug-resistant TB.

Recurrence of tuberculosis in a low-incidence setting: a retrospective cross-sectional study augmented by whole genome sequencing

Background

The recurrence of tuberculosis (TB) disease in treated patients can serve as a marker of the efficacy of TB control programs. Recurrent disease represents either endogenous reactivation with the same strain of Mycobacterium tuberculosis due to non-compliance or inadequate therapy or exogenous reinfection with a new strain. Genotyping or whole genome sequencing (WGS) of M. tuberculosis isolates from initial and recurrent cases can differentiate between reinfection and reactivation. This study examined cases of recurrent TB in New South Wales, Australia, using genotyping and WGS.

Methods

Culture-confirmed TB cases diagnosed at least 12 months apart between January 2011 and December 2016 were included. Isolates of M. tuberculosis from patients were compared using 24-locus Mycobacterial Interspersed Repetitive Unit Variable Number Tandem Repeat (MIRU-24) typing and WGS.

Results

Eighteen cases of recurrent disease were identified but isolates from only 15 (83%) were available for study. MIRU-24 findings classified 13 (13/15; 87%) as reactivation and two (13%), as reinfection. Sequencing 13 cultivable paired isolates demonstrated 11 reactivations and two reinfections. There was genomic similarity in 10 out of 13 pairs while one case (1/13; 8%) had 12 SNPS differences. Two other cases (2/13;15%) had > 200 SNPs differences and were classified as reinfection. No phenotypic or genomic evidence of drug resistance was observed.

Conclusion

TB control programs can achieve consistently low rates of recurrent disease in low incidence settings. WGS of implicated isolates augments the differentiation between reactivation and reinfection and indicates that the majority of recurrences are due to reactivation rather than reinfection. Predominance of reactivation over reinfection indicates high-quality public health practices and a low risk of local transmission.

Trial registration

This study was approved by the Western Sydney Local Health District (WSLHD) Human Research Ethics Committee (HREC Ref: AU RED LNR/17/WMEAD/190; SSA Ref: LNR SSA/17/WMEAD/191).

A molecular epidemiological analysis of tuberculosis trends in South Korea

Molecular epidemiological data are needed to assess tuberculosis (TB)-management policy outcomes in South Korea. IS6110 restriction fragment-length polymorphism (IS6110-RFLP) and mycobacterial interspersed repetitive unit-variable-number tandem repeat (MIRU-VNTR) analyses are major molecular epidemiological tools for investigating the transmission or reactivation of active TB. Here, we determined trends in the clustering rate (i.e., the prevalence of Mycobacterium tuberculosis isolates with identical genotype patterns) of active TB and related differences between the 1990s and 2000s in Korea. M. tuberculosis isolates (1,007) of nationwide origins were analyzed by IS6110-RFLP and 24-locus standardized MIRU-VNTR genotyping. The clustering rate was measured by IS6110-RFLP, 24-locus MIRU-VNTR, and both analytical methods in combination. IS6110-RFLP, 24-locus MIRU-VNTR typing, and the combined method revealed 882, 754, and 983 distinct profiles; 809, 651, and 961 unique isolates; and 198, 356, and 46 clustered isolates grouped into 73, 103, and 22 clusters, respectively. In addition, we confirmed that the clustering rates in the 2000s decreased by 11.2%, 2.1%, and 3.1% relative to that in the 1990s using the three methods, respectively. Furthermore, in multivariate analysis, the younger-age group (<30) clustered more frequently than the older-age group (>50), based on all the three methods. Our study is the first report to provide nationwide molecular epidemiological information on TB in Korea.

Population Structure and Local Adaptation of MAC Lung Disease Agent Mycobacterium avium subsp. hominissuis

Mycobacterium avium subsp. hominissuis (MAH) is one of the most common nontuberculous mycobacterial species responsible for chronic lung disease in humans. Despite increasing worldwide incidence, little is known about the genetic mechanisms behind the population evolution of MAH. To elucidate the local adaptation mechanisms of MAH, we assessed genetic population structure, the mutual homologous recombination, and gene content for 36 global MAH isolates, including 12 Japanese isolates sequenced in the present study. We identified five major MAH lineages and found that extensive mutual homologous recombination occurs among them. Two lineages (MahEastAsia1 and MahEastAsia2) were predominant in the Japanese isolates. We identified alleles unique to these two East Asian lineages in the loci responsible for trehalose biosynthesis (treS and mak) and in one mammalian cell entry operon, which presumably originated from as yet undiscovered mycobacterial lineages. Several genes and alleles unique to East Asian strains were located in the fragments introduced via recombination between East Asian lineages, suggesting implication of recombination in local adaptation. These patterns of MAH genomes are consistent with the signature of distribution conjugative transfer, a mode of sexual reproduction reported for other mycobacterial species.

Direct genotyping of Mycobacterium tuberculosis from Xpert® MTB/RIF remnants

Genotyping of Mycobacterium tuberculosis (MTB) isolates has markedly improved our knowledge of its transmission dynamics. MIRU-VNTR is considered the reference molecular tool for MTB fingerprinting. However, the dependence of this technique on cultured isolates means that we lack molecular epidemiology data from many settings where culture facilities have not been implemented. Efforts have been made to adapt the MIRU-VNTR procedure to direct analysis of clinical specimens, although implementation of these efforts has not proven successful. The large-scale roll-out of Xpert MTB/RIF (Xpert) technology, which is now in almost every TB-endemic country, including many where MTB is not cultured, provides us with a new opportunity to explore whether MTB genotyping could be performed from the remnants of the Xpert cartridge. We ran a pilot study in Mozambique in which the remnants of 24 positive Xpert assays for detection of MTB were used as template material for the 15-locus or the more discriminatory 24-locus MIRU-VNTR technique. MTB fingerprinting was possible in specimens with a high bacterial burden, according to the Xpert load categories, and within the first week after Xpert was performed. Given the wide availability, simple processing, and rapid reporting of results with Xpert, our findings suggest that MIRU-VNTR-based fingerprinting from remnants of Xpert could play a major role in extending MTB molecular epidemiology studies to settings where information on the transmission dynamics of this pathogen is lacking.

High genetic diversity among Mycobacterium tuberculosis strains in Tehran, Iran

INTRODUCTION

Tuberculosis (TB) still remains an important public health problem in Iran. The genotyping of Mycobacterium tuberculosis isolates is expected to lead to a better understanding of M. tuberculosis transmission in Tehran, the most populated city of Iran.

MATERIALS AND METHODS

A total of 2300 clinical specimens were obtained from TB suspected patients who were referred to a TB center in Tehran from Jan 2014 to Dec 2016. Identification was performed using both conventional and molecular methods. The presence of resistance to rifampicin was examined by the GeneXpert MTB/RIF. The standard 15-locus mycobacterial interspersed repetitive units-variable number of tandem repeats (MIRU-VNTR) typing method was applied to genotype of clinical isolates.

RESULTS

Of 2300 specimens, 80 isolates were identified as M. tuberculosis by using biochemical and molecular tests. Of 80 M. tuberculosis isolates, 76 (95%) had unique genotypic profiles and 4 (5%) shared a profile with one or more other strains. Based on single loci variation (SLV) 4 clonal complexes were observed. NEW-1 was found to be the most predominant lineage (22.5%) followed by West African (1.25%), Central Asian (CAS)/Delhi (1.25%), Bovis (1.25%), H37Rv (1.25%) and multiple matches (1.25%). Loci MIRU10, MIRU26, MTUB21 and QUB26 were found as highly discriminative. No mutation was detected in the hotspot region of rifampicin by using GeneXpert MTB/RIF.

CONCLUSIONS

Our study findings show that there was considerable genotypic diversity among M. tuberculosis isolates in Tehran. The 15-locus MIRU-VNTR showed high HGDI and could be used as a first-line genotyping method for epidemiological studies.

Prevalence of pyrazinamide resistance and Wayne assay performance analysis in a tuberculosis cohort in Lima, Peru

BACKGROUND

Multidrug-resistant tuberculosis (MDR-TB) regimens often contain pyrazinamide (PZA) even if susceptibility to the drug has not been confirmed. This gap is due to the limited availability and reliability of PZA susceptibility testing.

OBJECTIVES

To estimate the prevalence of PZA resistance using the Wayne assay among TB patients in Lima, Peru, to describe characteristics associated with PZA resistance and to compare the performance of Wayne with that of BACTEC™ MGIT™ 960.

METHODS

PZA susceptibility using the Wayne assay was tested in patients diagnosed with culture-positive pulmonary TB from September 2009 to August 2012. Factors associated with PZA resistance were evaluated. We compared the performance of the Wayne assay to that of MGIT 960 in a convenience sample.

RESULTS

The prevalence of PZA resistance was 6.6% (95%CI 5.8-7.5) among 3277 patients, and 47.7% (95%CI 42.7-52.6) among a subset of 405 MDR-TB patients. In multivariable analysis, MDR-TB (OR 86.0, 95%CI 54.0-136.9) and Latin American-Mediterranean lineage (OR 3.40, 95%CI 2.33-4.96) were associated with PZA resistance. The Wayne assay was in agreement with MGIT 960 in 83.9% of samples (κ 0.66, 95%CI 0.56-0.76).

CONCLUSION

PZA resistance was detected using the Wayne assay in nearly half of MDR-TB patients in Lima. This test can inform the selection and composition of regimens, especially those dependent on additional resistance.

Universal Genotyping for Tuberculosis Prevention Programs: a 5-Year Comparison with On-Request Genotyping

Prospective universal genotyping of tuberculosis (TB) isolates is used by many laboratories to detect clusters of cases and inform contact investigations. Prior to universal genotyping, most TB prevention programs genotyped isolates on request only, relying on requests from public health professionals whose knowledge of a patient's clinical, demographic, and epidemiological characteristics suggested potential transmission. To justify the switch from on-request to universal genotyping—particularly in the public health domain, with its limited resources and competing priorities—it is important to demonstrate the additional benefit provided by a universal genotyping program. We compared the clustering patterns revealed by retrospective 24-locus mycobacterial interspersed repetitive unit–variable-number tandem repeat genotyping of all culture-positive isolates over a 5-year period to the patterns previously established by our genotyping-on-request program in the low-incidence setting of British Columbia, Canada. We found that 23.8% of isolates were requested during the study period, and while requested isolates had increased odds of belonging to a genotype cluster (adjusted odds ratio, 2.3; 95% confidence interval, 1.5 to 3.3), only 54.6% clustered with the requested comparator strain. Universal genotyping revealed 94 clusters ranging in size from 2 to 53 isolates (mean = 5) and involving 432 individuals. On-request genotyping missed 54 (57.4%) of these clusters and 130 (30.1%) clustered individuals. Our results underscore that TB patient networks are complex, with unrecognized linkages between patients, and a prospective province-wide universal genotyping program provides an informative, bias-free tool to explore transmission to a degree not possible with on-request genotyping.

The Epidemiological Significance and Temporal Stability of Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats-Based Method Applied to Mycobacterium tuberculosis in China

This study aimed to validate the epidemiological significance and temporal stability of Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats (MIRU-VNTR) typing in a genetically and geographically diverse set of clinical isolates from patients diagnosed with pulmonary tuberculosis in China. Between 2010 and 2013, a total of 982 Mycobacterium tuberculosis isolates were collected from four population-based investigations in China. Apart from the currently applied 24-locus MIRU-VNTR, six additional hypervariable loci were analyzed in order to validate the MIRU-VNTR combinations in terms of their epidemiological links, clustering time span, and paired geographic distance. In vitro temporal stability was analyzed for both individual MIRU-VNTR loci, and for several combinations of loci. In the present study, four MIRU-VNTR combinations, including the hypervariable loci 3820, 3232, 2163a, and 4120, were evaluated. All of these combinations obtained a Hunter-Gaston discriminatory index (HGDI) value over 0.9900 with a reduced clustering proportion (from 32.0% to 25.6%). By comparing epidemiological links, clustering time span, and paired geographic distance, we found that the performances of the four MIRU-VNTR combinations were comparable to the insertion sequence 6110 restriction fragment length polymorphism (IS6110-RFLP), and significantly better than that of 24-locus MIRU-VNTR genotyping alone. The proportion of temporally stable loci ranged from 90.5% to 92.5% within the combined MIRU-VNTR genotyping, which is higher than IS6110-RFLP (85.4%). By adding four hypervariable loci to the standard 24-locus MIRU-VNTR genotyping, we obtained a high discriminatory power, stability and epidemiological significance. This algorithm could therefore be used to improve tuberculosis transmission surveillance and outbreak investigation in China.

Phylogenetic uniqueness of Mycobacterium avium subspecies hominissuis isolated from an abnormal pulmonary bovine case

Mycobacterium avium subspecies hominissuis (MAH) is an important cause of infection in human pulmonary and swine intestinal cases. Although MAH is isolated from environmental sources frequently, infections of other animals have rarely been analysed. Recently, we detected granulomatous inflammation in bovine lung as an abnormal postmortem inspection case. To ascertain its genetic profile, we conducted a variable numbers of tandem repeats (VNTR) analysis and genomic characterization using deep sequencing. The VNTR type was a unique profile that differed from reported genotypes, but it was assigned within a broad genotypic complex of isolates from human patients and bathrooms. Genomic comparison with 116 registered genome sequences of the subspecies revealed that the strain was separate from five major genetic population groups proposed previously. Although the infection source remains unclear, its isolation from various resources such as animal infection cases should be elucidated more extensively to reveal its genetic diversity and ecological context.

New insights into the transposition mechanisms of IS6110 and its dynamic distribution between Mycobacterium tuberculosis Complex lineages

The insertion Sequence IS6110, only present in the pathogens of the Mycobacterium tuberculosis Complex (MTBC), has been the gold-standard epidemiological marker for TB for more than 25 years, but biological implications of IS6110 transposition during MTBC adaptation to humans remain elusive. By studying 2,236 clinical isolates typed by IS6110-RFLP and covering the MTBC, we remarked a lineage-specific content of IS6110 being higher in modern globally distributed strains. Once observed the IS6110 distribution in the MTBC, we selected representative isolates and found a correlation between the normalized expression of IS6110 and its abundance in MTBC chromosomes. We also studied the molecular regulation of IS6110 transposition and we found a synergistic action of two post-transcriptional mechanisms: a -1 ribosomal frameshift and a RNA pseudoknot which interferes translation. The construction of a transcriptionally active transposase resulted in 20-fold increase of the transposition frequency. Finally, we examined transposition in M. bovis and M. tuberculosis during laboratory starvation and in a mouse infection model of TB. Our results shown a higher transposition in M. tuberculosis, that preferably happens during TB infection in mice and after one year of laboratory culture, suggesting that IS6110 transposition is dynamically adapted to the host and to adverse growth conditions.

Since the pioneering discovery of transposition by Barbara McClintock in eukaryotes and later in prokaryotes by Robert W. Hedges and Alan E. Jacob, it has become clear the key role of mobile genetics elements in chromosome remodelling, microbial evolution and host adaptation. The insertion sequence IS6110 is widely recognized for its utility in TB diagnosis and epidemiology because it is only present in the M. tuberculosis Complex (MTBC) and its transposition provides an excellent chromosomal polymorphic variability allowing the study of recent TB transmission. This inherent feature of IS6110 leads us to hypothesize that IS6110 plays a crucial role during the TB infectious cycle. However, the biological significance of IS6110 has been hindered by its almost exclusive use as an epidemiological marker. Here, we study the regulatory mechanisms and the distribution of IS6110 in the different MTBC lineages. We discuss the potential biological implications of IS6110, that is much more than an excellent TB epidemiological tool. Since IS6110 could play an important role in the adaptation of MTBC to the host, this study opens new avenues to decipher the biological roles of IS6110 in TB pathogenesis.