Genomic diversity among drug sensitive and multidrug resistant isolates of Mycobacterium tuberculosis with identical DNA fingerprints

Applications and advances in molecular diagnostics: revolutionizing non-tuberculous mycobacteria species and subspecies identification

Non-tuberculous mycobacteria (NTM) infections pose a significant public health challenge worldwide, affecting individuals across a wide spectrum of immune statuses. Recent epidemiological studies indicate rising incidence rates in both immunocompromised and immunocompetent populations, underscoring the need for enhanced diagnostic and therapeutic approaches. NTM infections often present with symptoms similar to those of tuberculosis, yet with less specificity, increasing the risk of misdiagnosis and potentially adverse outcomes for patients. Consequently, rapid and accurate identification of the pathogen is crucial for precise diagnosis and treatment. Traditional detection methods, notably microbiological culture, are hampered by lengthy incubation periods and a limited capacity to differentiate closely related NTM subtypes, thereby delaying diagnosis and the initiation of targeted therapies. Emerging diagnostic technologies offer new possibilities for the swift detection and accurate identification of NTM infections, playing a critical role in early diagnosis and providing more accurate and comprehensive information. This review delineates the current molecular methodologies for NTM species and subspecies identification. We critically assess the limitations and challenges inherent in these technologies for diagnosing NTM and explore potential future directions for their advancement. It aims to provide valuable insights into advancing the application of molecular diagnostic techniques in NTM infection identification.

K-mer applied in Mycobacterium tuberculosis genome cluster analysis

Abstract According to studies carried out, approximately 10 million people developed tuberculosis in 2018. Of this total, 1.5 million people died from the disease. To study the behavior of the genome sequences of Mycobacterium tuberculosis (MTB), the bacterium responsible for the development of tuberculosis (TB), an analysis was performed using k-mers (DNA word frequency). The k values ranged from 1 to 10, because the analysis was performed on the full length of the sequences, where each sequence is composed of approximately 4 million base pairs, k values above 10, the analysis is interrupted, as consequence of the program's capacity. The aim of this work was to verify the formation of the phylogenetic tree in each k-mer analyzed. The results showed the formation of distinct groups in some k-mers analyzed, taking into account the threshold line. However, in all groups, the multidrug-resistant (MDR) and extensively drug-resistant (XDR) strains remained together and separated from the other strains.

Clustering minimal inhibitory concentration data through Bayesian mixture models: An application to detect Mycobacterium tuberculosis resistance mutations

Antimicrobial resistance is becoming a major threat to public health throughout the world. Researchers are attempting to contrast it by developing both new antibiotics and patient-specific treatments. In the second case, whole-genome sequencing has had a huge impact in two ways: first, it is becoming cheaper and faster to perform whole-genome sequencing, and this makes it competitive with respect to standard phenotypic tests; second, it is possible to statistically associate the phenotypic patterns of resistance to specific mutations in the genome. Therefore, it is now possible to develop catalogues of genomic variants associated with resistance to specific antibiotics, in order to improve prediction of resistance and suggest treatments. It is essential to have robust methods for identifying mutations associated to resistance and continuously updating the available catalogues. This work proposes a general method to study minimal inhibitory concentration distributions and to identify clusters of strains showing different levels of resistance to antimicrobials. Once the clusters are identified and strains allocated to each of them, it is possible to perform regression method to identify with high statistical power the mutations associated with resistance. The method is applied to a new 96-well microtiter plate used for testing Mycobacterium tuberculosis.

Early alveolar macrophage response and IL-1R-dependent T cell priming determine transmissibility of Mycobacterium tuberculosis strains

Mechanisms underlying variability in transmission of Mycobacterium tuberculosis strains remain undefined. By characterizing high and low transmission strains of M.tuberculosis in mice, we show here that high transmission M.tuberculosis strain induce rapid IL-1R-dependent alveolar macrophage migration from the alveolar space into the interstitium and that this action is key to subsequent temporal events of early dissemination of bacteria to the lymph nodes, Th1 priming, granulomatous response and bacterial control. In contrast, IL-1R-dependent alveolar macrophage migration and early dissemination of bacteria to lymph nodes is significantly impeded in infection with low transmission M.tuberculosis strain; these events promote the development of Th17 immunity, fostering neutrophilic inflammation and increased bacterial replication. Our results suggest that by inducing granulomas with the potential to develop into cavitary lesions that aids bacterial escape into the airways, high transmission M.tuberculosis strain is poised for greater transmissibility. These findings implicate bacterial heterogeneity as an important modifier of TB disease manifestations and transmission.

High prevalence of multiple antibiotic resistance in clinical E. coli isolates from Bangladesh and prediction of molecular resistance determinants using WGS of an XDR isolate

Multi-drug-resistance (MDR) is a severe public health concern worldwide, and its containment is more challenging in developing countries due to poor antimicrobial resistance (AMR) surveillance and irrational use of antibiotics. The current study investigated 100 clinical E. coli isolates and revealed that 98% of them were MDR. PCR analysis using 25 selected isolates showed the predominance of metallo-β-lactamase gene bla_NDM (80%) and ESBL genes bla_OXA (48%) and bla_CTX-M-15 (32%). The AmpC gene was detected in 68% of the isolates, while 32% was tetC positive. Notably, 34% of the isolates were resistant to carbapenem. Whole genome sequence (WGS) analysis of an extensively drug-resistant (XDR) isolate (L16) revealed the presence of the notorious sequence type 131 responsible for multi-drug-resistant infections, multiple antibiotic resistance genes (ARGs), virulence genes, and mobile genetic elements that pose risks to environmental transmission. Our results indicate that MDR is alarmingly increasing in Bangladesh that critically limits the treatment option against infections and contributes to further aggravation to the prevailing situation of MDR worldwide. The findings of this study will be valuable in designing sustainable strategies to contain MDR in the region.

Whole genome analysis of extensively drug resistant Mycobacterium tuberculosis strains in Peru

Peru has the highest burden of multidrug-resistant tuberculosis in the Americas region. Since 1999, the annual number of extensively drug-resistant tuberculosis (XDR-TB) Peruvian cases has been increasing, becoming a public health challenge. The objective of this study was to perform genomic characterization of Mycobacterium tuberculosis strains obtained from Peruvian patients with XDR-TB diagnosed from 2011 to 2015 in Peru. Whole genome sequencing (WGS) was performed on 68 XDR-TB strains from different regions of Peru. 58 (85.3%) strains came from the most populated districts of Lima and Callao. Concerning the lineages, 62 (91.2%) strains belonged to the Euro-American Lineage, while the remaining 6 (8.8%) strains belonged to the East-Asian Lineage. Most strains (90%) had high-confidence resistance mutations according to pre-established WHO-confident grading system. Discordant results between microbiological and molecular methodologies were caused by mutations outside the hotspot regions analysed by commercial molecular assays (rpoB I491F and inhA S94A). Cluster analysis using a cut-off ≤ 10 SNPs revealed that only 23 (34%) strains evidenced recent transmission links. This study highlights the relevance and utility of WGS as a high-resolution approach to predict drug resistance, analyse transmission of strains between groups, and determine evolutionary patterns of circulating XDR-TB strains in the country.

Molecular diversity of Mycobacterium tuberculosis complex in Sikkim, India and prediction of dominant spoligotypes using artificial intelligence

In India, tuberculosis is an enormous public health problem. This study provides the first description of molecular diversity of the Mycobacterium tuberculosis complex (MTBC) from Sikkim, India. A total of 399 Acid Fast Bacilli sputum positive samples were cultured on Lőwenstein-Jensen media and genetic characterisation was done by spoligotyping and 24-loci MIRU-VNTR typing. Spoligotyping revealed the occurrence of 58 different spoligotypes. Beijing spoligotype was the most dominant type constituting 62.41% of the total isolates and was associated with Multiple Drug Resistance. Minimum Spanning tree analysis of 249 Beijing strains based on 24-loci MIRU-VNTR analysis identified 12 clonal complexes (Single Locus Variants). The principal component analysis was used to visualise possible grouping of MTBC isolates from Sikkim belonging to major spoligotypes using 24-MIRU VNTR profiles. Artificial intelligence-based machine learning (ML) methods such as Random Forests (RF), Support Vector Machines (SVM) and Artificial Neural Networks (ANN) were used to predict dominant spoligotypes of MTBC using MIRU-VNTR data. K-fold cross-validation and validation using unseen testing data set revealed high accuracy of ANN, RF, and SVM for predicting Beijing, CAS1_Delhi, and T1 Spoligotypes (93-99%). However, prediction using the external new validation data set revealed that the RF model was more accurate than SVM and ANN.

Whole genome sequencing analysis of multidrug-resistant tuberculosis in Singapore, 2006-2018

There were 290 multidrug-resistant (MDR)-TB cases diagnosed in Singapore from 2006 to 2018. Eighty-one percent were foreign-born. Spoligotyping and MIRU-VNTR methods identified 108 patients in 24 clusters. The Beijing spoligotype accounted for 22 clusters. Whole genome sequencing (WGS) analysis reduced the number of clustered patients and clusters to 43 and nine respectively. One MIRU cluster was redefined into three WGS clusters. All the clusters had foreign-born source cases. Forty percent of local-born, versus 9% of foreign-born, MDR-TB cases belonged to WGS clusters. WGS more accurately elucidated potential MDR-TB transmission which was overestimated by conventional genotyping methods in Singapore.

Multi-Omics Technologies Applied to Tuberculosis Drug Discovery

Multi-omics strategies are indispensable tools in the search for new anti-tuberculosis drugs. Omics methodologies, where the ensemble of a class of biological molecules are measured and evaluated together, enable drug discovery programs to answer two fundamental questions. Firstly, in a discovery biology approach, to find new targets in druggable pathways for target-based investigation, advancing from target to lead compound. Secondly, in a discovery chemistry approach, to identify the mode of action of lead compounds derived from high-throughput screens, progressing from compound to target. The advantage of multi-omics methodologies in both of these settings is that omics approaches are unsupervised and unbiased to a priori hypotheses, making omics useful tools to confirm drug action, reveal new insights into compound activity, and discover new avenues for inquiry. This review summarizes the application of Mycobacterium tuberculosis omics technologies to the early stages of tuberculosis antimicrobial drug discovery.

A retrospective genomic analysis of drug-resistant strains of M. tuberculosis in a high-burden setting, with an emphasis on comparative diagnostics and reactivation and reinfection status

Background

Recurrence of drug-resistant tuberculosis (DR-TB) after treatment occurs through relapse of the initial infection or reinfection by a new drug-resistant strain. Outbreaks of DR-TB in high burden regions present unique challenges in determining recurrence status for effective disease management and treatment. In the Republic of Moldova the burden of DR-TB is exceptionally high, with many cases presenting as recurrent.

Methods

We performed a retrospective analysis of Mycobacterium tuberculosis from Moldova to better understand the genomic basis of drug resistance and its effect on the determination of recurrence status in a high DR-burden environment. To do this we analyzed genomes from 278 isolates collected from 189 patients, including 87 patients with longitudinal samples. These pathogen genomes were sequenced using Illumina technology, and SNP panels were generated for each sample for use in phylogenetic and network analysis. Discordance between genomic resistance profiles and clinical drug-resistance test results was examined in detail to assess the possibility of mixed infection.

Results

There were clusters of multiple patients with 10 or fewer differences among DR-TB samples, which is evidence of person-to-person transmission of DR-TB. Analysis of longitudinally collected isolates revealed that many infections exhibited little change over time, though 35 patients demonstrated reinfection by divergent (number of differences > 10) lineages. Additionally, several same-lineage sample pairs were found to be more divergent than expected for a relapsed infection. Network analysis of the H3/4.2.1 clade found very close relationships among 61 of these samples, making differentiation of reactivation and reinfection difficult. There was discordance between genomic profile and clinical drug sensitivity test results in twelve samples, and four of these had low level (but not statistically significant) variation at DR SNPs suggesting low-level mixed infections.

Conclusions

Whole-genome sequencing provided a detailed view of the genealogical structure of the DR-TB epidemic in Moldova, showing that reinfection may be more prevalent than currently recognized. We also found increased evidence of mixed infection, which could be more robustly characterized with deeper levels of genomic sequencing.

Computational pan-genome mapping and pairwise SNP-distance improve detection of Mycobacterium tuberculosis transmission clusters

Next-generation sequencing based base-by-base distance measures have become an integral complement to epidemiological investigation of infectious disease outbreaks. This study introduces PANPASCO, a computational pan-genome mapping based, pairwise distance method that is highly sensitive to differences between cases, even when located in regions of lineage specific reference genomes. We show that our approach is superior to previously published methods in several datasets and across different Mycobacterium tuberculosis lineages, as its characteristics allow the comparison of a high number of diverse samples in one analysis—a scenario that becomes more and more likely with the increased usage of whole-genome sequencing in transmission surveillance.

Tuberculosis still is a threat to global health. It is essential to detect and interrupt transmissions to stop the spread of this infectious disease. With the rising use of next-generation sequencing methods, its application in the surveillance of Mycobacterium tuberculosis has become increasingly important in the last years. The main goal of molecular surveillance is the identification of patient-patient transmission and cluster detection. The mutation rate of M. tuberculosis is very low and stable. Therefore, many existing methods for comparative analysis of isolates provide inadequate results since their resolution is too limited. There is a need for a method that takes every detectable difference into account. We developed PANPASCO, a novel approach for comparing pairs of isolates using all genomic information available for each pair. We combine improved SNP-distance calculation with the use of a pan-genome incorporating more than 100 M. tuberculosis reference genomes representing lineages 1-4 for read mapping prior to variant detection. We thereby enable the collective analysis and comparison of similar and diverse isolates associated with different M. tuberculosis strains.

Genotyping of Mycobacterium tuberculosis spreading in Hanoi, Vietnam using conventional and whole genome sequencing methods

Hanoi is the capital of Vietnam, one of the 30 countries with a high tuberculosis (TB) burden. Fundamental data on the molecular epidemiology of the disease is required for future TB management. To identify lineages and genotypes of Mycobacterium tuberculosis (Mtb), conventional genotyping data from clinical isolates of the Hanoi area was compared with whole genome sequencing (WGS) analysis from 332 of 470 samples. It was obtained from lineage-specific single nucleotide variants (SNVs), large sequence polymorphisms, spoligotyping, and variable number of tandem repeats (VNTR) analysis using mycobacterial interspersed repetitive unit (MIRU) and Japan anti-tuberculosis association (JATA) locus sets. This information was directly compared with results obtained from WGS. Mini-satellite repeat unit variants were identified using BLAST search against concatenated short read sequences, the RepUnitTyping tool. WGS analysis revealed that the Mtb strains tested are diverse and classified into lineage (L) 1, 2 and 4 (24.7, 57.2 and 18.1% respectively). The majority of the L2 strains were further divided into ancient and modern Beijing genotypes, and most of the L1 group were EAI4_VNM strains. Although conventional PCR-based genotyping results were mostly consistent with information obtained through WGS analysis, in-depth analysis identified aberrant deletions and spacers that may cause discordance. JATA-VNTR sets, including hypervariable loci, separated large Beijing genotypic clusters generated by MIRU15 into smaller groups. The distribution of repeat unit variants observed within 33 VNTR loci showed clear variation depending on the three lineages. WGS-based pairwise-SNV differences within VNTR-defined genotypic clusters were greater in L1 than in L2 and L4 (P = .001). Direct comparisons between results of PCR-based genotyping and in silico analysis of WGS data would bridge a gap between classical and modern technologies during this transition period, and provide further information on Mtb genotypes in specific geographical areas.

Whole genome sequencing provides additional insights into recurrent tuberculosis classified as endogenous reactivation by IS6110 DNA fingerprinting

Recurrent tuberculosis (TB) after successful TB treatment occurs due to endogenous reactivation (relapse) or exogenous reinfection. We revisited the conclusions of relapse in a high TB incidence setting that were drawn on the basis of IS6110 restriction fragment length polymorphism (RFLP) analysis in a large retrospective cohort study in suburban Cape Town, South Africa. Using whole genome sequencing (WGS), we undertook pair-wise genome comparison of Mycobacterium tuberculosis strains cultured from diagnostic sputum samples collected at the index and recurrent TB episode for 25 recurrent TB cases who had been classified as relapse based on identical DNA fingerprint patterns in the earlier study. We found that paired strain genome sequences were identical or showed minimal variant differences in 22 of 25 recurrent TB cases, consistent with relapse. One showed 20 variant differences, suggestive of exogenous reinfection. Two of the 25 had mixed infections, each with the index episode strain detected as the dominant strain at recurrence in one of these patients, the minority strain harboured drug-resistance conferring mutations (rpoB, katG). In conclusion, our study highlights the additional value of WGS for investigating recurrent TB in settings with high infection pressure and closely related circulating strains, where the extent of re- and mixed infection may be underestimated.

Genomic characterization of MDR/XDR-TB in Kazakhstan by a combination of high-throughput methods predominantly shows the ongoing transmission of L2/Beijing 94-32 central Asian/Russian clusters

Background

Kazakhstan remains a high-burden TB prevalence country with a concomitent high-burden of multi-drug resistant tuberculosis. For this reason, we performed an in depth genetic diversity and population structure characterization of Mycobacterium tuberculosis complex (MTC) genetic diversity in Kazakhstan with both patient and community benefit.

Methods

A convenience sample of 700 MTC DNA cultures extracts from 630 tuberculosis patients recruited from 12 out of 14 regions in Kazakhstan, between 2010 and 2015, was independently studied by high-throughput hybridization-based methods, TB-SPRINT (59-Plex, n = 700), TB-SNPID (50-Plex, n = 543). DNA from 391 clinical isolates was successfully typed by two methods. To resolve the population structure of drug-resistant clades in more detail two complementary assays were run on the L2 isolates: an IS6110-NTF insertion site typing assay and a SigE SNP polymorphism assay.

Results

Strains belonged to L2/Beijing and L4/Euro-American sublineages; L2/Beijing prevalence totaled almost 80%. 50% of all samples were resistant to RIF and to INH., Subtyping showed that: (1) all L2/Beijing were “modern” Beijing and (2) most of these belonged to the previously described 94–32 sublineage (Central Asian/Russian), (3) at least two populations of the Central Asian/Russian sublineages are circulating in Kazakhstan, with different evolutionary dynamics.

Conclusions

For the first time, the global genetic diversity and population structure of M. tuberculosis genotypes circulating in Kazakhstan was obtained and compared to previous local studies. Results suggest a region-specific spread of a very limited number of L2/Beijing clonal complexes in Kazakhstan many strongly associated with an MDR phenotype.

Electronic supplementary material

The online version of this article (10.1186/s12879-019-4201-2) contains supplementary material, which is available to authorized users.

Clonality and genetic profiles of drug-resistant Mycobacterium tuberculosis in the Eastern Cape Province, South Africa

In this study, we investigated the diversity of drug-resistant Mycobacterium tuberculosis isolates from families who own cattle in the Eastern Cape Province of South Africa using spoligotyping and mycobacterial interspersed repetitive-unit-variable number tandem repeat (MIRU-VNTR) typing. The Mycobacterium tuberculosis was investigated using MIRU-VNTR and the Mycobacterium tuberculosis families were evaluated using spoligotyping. Spoligotyping grouped 91% of the isolates into seven clusters, while 9% of the deoxyribonucleic acid (DNA) from TB isolates were unclustered from a total of 154 DNA used. Previously described shared types were observed in 89.6% of the isolates, with the Beijing family, SIT1, the principal genotype in the province, while the families T, SIT53 and X1, SIT1329 were the least detected genotypes. MIRU-VNTR grouped 81% of the isolates in 23 clusters while 19% were unclustered. A combination of the VNTR and spoligotyping grouped 79% of the isolates into 23 clusters with 21% unclustered. The low level of diversity and the clonal spread of drug-resistant Mycobacterium tuberculosis isolates advocate that the spread of TB in this study may be instigated by the clonal spread of Beijing genotype. The results from this study provide vital information about the lack of TB control and distribution of Mycobacterium tuberculosis complex strain types in the Eastern Cape Province of South Africa.

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.

Antibiotics Resistance Genes Screening and Comparative Genomics Analysis of Commensal Escherichia coli Isolated from Poultry Farms between China and Sudan

Escherichia coli (E. coli) strains, from the gut of animals and humans, harbor wide range of drug resistance genes. A comparative study is conducted on the intestinal E. coli from fecal samples of healthy chicken from China and Sudan in order to monitor the antimicrobial sensitivity pattern. A number of 250 E. coli isolates from chicken farms, including 120 from China and 130 from Sudan, were isolated and identified. All isolates were subjected to susceptibility tests against 10 antibiotics and the distribution of antibiotic resistant genes was confirmed by PCR amplification, involving genes such as ampC, tetA, pKD13, acrA, ermA, ermB, ermC, tetB, mphA, aadA14, aadA1, aac3-1, and aac3- III. Many isolates were found to exhibit resistance against more than one antibiotic. However, the Chinese isolates showed more antibiotics resistance and resistance genes compared to the Sudanese isolates. For better understanding of the multidrug resistance factors, we conducted whole genome analyses of E. coli D107 isolated from China, which revealed that the genome possesses multiple resistance genes including tetracycline, erythromycin, and kanamycin. Furthermore, E. coli D4 isolate from Sudan was more sensitive to antibiotics such as erythromycin, tetracycline, and gentamicin. After analysis by RAST and MAUVE, the two strains showed 89% average nucleotide identity. However, the genomes mostly differed at the number of antibiotics-related genes, as the genome of D107 revealed a considerable number of antibiotics resistance genes such as ermA and mphD which were found to be absent in D4 genome. These outcomes provided confirmation that the poultry farms environment in different countries (China and Sudan) may serve as a potential reservoir of antimicrobial resistance genes and also indicated the evolutionary differences of strains in terms of resistant genes expression.

Comparing Two Mycobacterium tuberculosis Genomes from Chinese Immigrants with Native Genomes Using Mauve Alignments

BACKGROUND

The number of immigrants with tuberculosis (TB) increases each year in South Korea. Determining the transmission dynamics based on whole genome sequencing (WGS) to cluster the strains has been challenging.

METHODS

WGS, annotation refinement, and orthology assignment for the GenBank accession number acquisition were performed on two clinical isolates from Chinese immigrants. In addition, the genomes of the two isolates were compared with the genomes of Mycobacterium tuberculosis isolates, from two native Korean and five native Chinese individuals using a phylogenetic topology tree based on the Multiple Alignment of Conserved Genomic Sequence with Rearrangements (Mauve) package.

RESULTS

The newly assigned accession numbers for two clinical isolates were CP020381.2 (a Korean-Chinese from Yanbian Province) and CP022014.1 (a Chinese from Shandong Province), respectively. Mauve alignment classified all nine TB isolates into a discriminative collinear set with matched regions. The phylogenetic analysis revealed a rooted phylogenetic tree grouping the nine strains into two lineages: strains from Chinese individuals and strains from Korean individuals.

CONCLUSION

Phylogenetic trees based on the Mauve alignments were supposed to be useful in revealing the dynamics of TB transmission from immigrants in South Korea, which can provide valuable information for scaling up the TB screening policy for immigrants.

Mycobacteriosis and Tuberculosis: Laboratory Diagnosis

Background

Tuberculosis is one of the most important infectious diseases that has claimed its victims throughout much of known human history. With Koch's discovery of the tubercle bacillus as the etiologic agent of the disease, his sanitary and hygienic measures, which were based on his discovery and the development of a vaccine against tuberculosis by Albert Calmette and Camille Guérin in 1921, an attenuated Mycobacterium bovis strain, bacilli Calmette-Guérin (BCG), and the discovery of the first antibiotic against tuberculosis, streptomycin by Selman Waksman in 1943, soon led to the opinion that appropriate control measures had become available for tuberculosis and it had been assumed that the disease could ultimately be eradicated.The emergence of resistant strains of this bacteria and widespread distribution of the disease in the world, and the emergence of the AIDS epidemic destroyed any possibility of global control of tuberculosis in the foreseeable future.

Objectives

The purpose of this review is to highlight the current scientific literature on mycobacterial infections and provide an overview on the laboratory diagnosis of tuberculosis and non-tuberculosis infections based on conventional phenotypic and modern molecular assays.

Method

In this study, a number of 65 papers comprising 20 reviews, 9 case reports, and 36 original research in association with mycobacteriosis and the laboratory diagnosis of mycobacterial infections, were reviewed.

Results

Based on our analysis on the published documents methods applied for the laboratory diagnosis of tuberculosis are continually assessed and developed in order to achieve more rapid, less expensive, and accurate results. Acid-fast staining and culture for mycobacteria remain at the core of any diagnostic algorithm with the sensitivity of 20-70% and specificity of 95-98% for AFB microscopy and the sensitivity of 95% and the specificity of 98% for culture based diagnosis. Following growth in culture, molecular tests such as nucleic acid hybridization probes and DNA sequencing may be used for definitive species identification. Nucleic acid amplification methods provide the means for direct detection of Mycobacterium tuberculosis in respiratory specimens without the prerequisite to isolate or culture the organism, leading to more rapid diagnosis and better patient care.

Conclusion

As the researchers in a developing country, we strongly believe that despite significant advances in laboratory capacity, in many countries reliable confirmation of suspected mycobacterial diseases is hindered by a lack of knowledge on proper standardized methods, sufficient funds, suitably trained staff and laboratory supplies.

Low Frequency of Acquired Isoniazid and Rifampicin Resistance in Rifampicin-Susceptible Pulmonary Tuberculosis in a Setting of High HIV-1 Infection and Tuberculosis Coprevalence

Background

We estimated the incidence of acquired isoniazid and rifampicin resistance in rifampicin-susceptible tuberculosis in a setting of high human immunodeficiency virus type 1 (HIV-1) infection and tuberculosis coprevalence.

Methods

GeneXpert MTB/RIF–confirmed patients with rifampicin-susceptible tuberculosis were recruited at antituberculosis treatment initiation in Khayelitsha, South Africa. Liquid culture and adherence assessment were performed at 2 and 5–6 months. MTBDRplus was performed on mycobacteria-positive cultures to ascertain acquired drug resistance (ADR). Spoligotyping and whole-genome sequencing were performed to ascertain homogeneity between baseline isolates and isolates with ADR. Baseline isolates were retrospectively tested for isoniazid monoresistance. An electronic database review was performed to ascertain tuberculosis recurrences.

Results

A total of 306 participants (62% with HIV-1 coinfection, of whom 71% received antiretroviral therapy) were recruited. Ascertainment of outcomes was complete for 284 participants. Five acquired a resistant Mycobacterium tuberculosis strain during or subsequent to treatment. One strain was confirmed to have ADR, 2 were confirmed as causing exogenous reinfection, and 2 were unrecoverable for genotyping. Incident ADR was estimated to have ranged from 0.3% (95% confidence interval [CI], .1%–1.9%; 1 of 284 participants) to 1% (95% CI, .2%–3%; 3 of 284 participants). Seventeen of 279 baseline isolates (6.1%; 95% CI, 3.6%–9.6%) had isoniazid monoresistance (13 of 17 had an inhA promoter mutation), but 0 of 17 had amplified resistance.

Conclusions

Treatment with standardized antituberculosis regimens dosed daily throughout, high uptake of antiretroviral therapy, and low prevalence of isoniazid monoresistance were associated with a low frequency of ADR.

Shaping the niche in macrophages: Genetic diversity of the M. tuberculosis complex and its consequences for the infected host

Pathogenic mycobacteria of the Mycobacterium tuberculosis complex (MTBC) have co-evolved with their individual hosts and are able to transform the hostile environment of the macrophage into a permissive cellular habitat. The impact of MTBC genetic variability has long been considered largely unimportant in TB pathogenesis. Members of the MTBC can now be distinguished into three major phylogenetic groups consisting of 7 phylogenetic lineages and more than 30 so called sub-lineages/subgroups. MTBC genetic diversity indeed influences the transmissibility and virulence of clinical MTBC isolates as well as the immune response and the clinical outcome. Here we review the genetic diversity and epidemiology of MTBC strains and describe the current knowledge about the host immune response to infection with MTBC clinical isolates using human and murine experimental model systems in vivo and in vitro. We discuss the role of innate cytokines in detail and portray two in our group recently developed approaches to characterize the intracellular niches of MTBC strains. Characterizing the niches and deciphering the strategies of MTBC strains to transform an antibacterial effector cell into a permissive cellular habitat offers the opportunity to identify strain- and lineage-specific key factors which may represent targets for novel antimicrobial or host directed therapies for tuberculosis.

Evolution of Mycobacterium tuberculosis: New Insights into Pathogenicity and Drug Resistance

The tuberculosis agent Mycobacterium tuberculosis has undergone a long and selective evolution toward human infection and represents one of the most widely spread pathogens due to its efficient aerosol-mediated human-to-human transmission. With the availability of more and more genome sequences, the evolutionary trajectory of this obligate pathogen becomes visible, which provides us with new insights into the molecular events governing evolution of the bacterium and its ability to accumulate drug-resistance mutations. In this review, we summarize recent developments in mycobacterial research related to this matter that are important for a better understanding of the current situation and future trends and developments in the global epidemiology of tuberculosis, as well as for possible public health intervention possibilities.

A systematic review of East African-Indian family of Mycobacterium tuberculosis in Brazil

Introduction

The Mycobacterium tuberculosis East African-Indian (EAI) spoligotyping family (belonging to lineage 1, Indo-Oceanic, defined by the region of deletion RD239) is distributed worldwide, but is more prevalent in Southeast Asia, India, and East Africa. Studies in Latin America have rarely identified EAI. In this study, we describe the occurrence of the EAI family in Brazil.

Methods

EAI was identified in a systematic literature review of genetic diversity studies pertaining to M. tuberculosis in Brazil, as well as in a survey conducted in Salvador, Bahia, located in the northeastern region of this country.

Results

The EAI6-BGD1 spoligotyping family and the EAI5 Spoligotype International Type (SIT) 1983 clade were the most frequently reported, with wide distribution of this particular clade described in Brazil. The distribution of other EAI spoligotyping patterns with broader worldwide distribution was restricted to the southeastern region of the country.

Conclusions

EAI may be endemic at a low frequency in Brazil, with some clades indicating increased fitness with respect to this population.

A joint cross-border investigation of a cluster of multidrug-resistant tuberculosis in Austria, Romania and Germany in 2014 using classic, genotyping and whole genome sequencing methods: lessons learnt

Molecular surveillance of multidrug-resistant tuberculosis (MDR-TB) using 24-loci MIRU-VNTR in the European Union suggests the occurrence of international transmission. In early 2014, Austria detected a molecular MDR-TB cluster of five isolates. Links to Romania and Germany prompted the three countries to investigate possible cross-border MDR-TB transmission jointly. We searched genotyping databases, genotyped additional isolates from Romania, used whole genome sequencing (WGS) to infer putative transmission links, and investigated pairwise epidemiological links and patient mobility. Ten isolates from 10 patients shared the same 24-loci MIRU-VNTR pattern. Within this cluster, WGS defined two subgroups of four patients each. The first comprised an MDR-TB patient from Romania who had sought medical care in Austria and two patients from Austria. The second comprised patients, two of them epidemiologically linked, who lived in three different countries but had the same city of provenance in Romania. Our findings strongly suggested that the two cases in Austrian citizens resulted from a newly introduced MDR-TB strain, followed by domestic transmission. For the other cases, transmission probably occurred in the same city of provenance. To prevent further MDR-TB transmission, we need to ensure universal access to early and adequate therapy and collaborate closely in tuberculosis care beyond administrative borders.

Some Synonymous and Nonsynonymous gyrA Mutations in Mycobacterium tuberculosis Lead to Systematic False-Positive Fluoroquinolone Resistance Results with the Hain GenoType MTBDRsl Assays

In this study, using the Hain GenoType MTBDRsl assays (versions 1 and 2), we found that some nonsynonymous and synonymous mutations in gyrA in Mycobacterium tuberculosis result in systematic false-resistance results to fluoroquinolones by preventing the binding of wild-type probes. Moreover, such mutations can prevent the binding of mutant probes designed for the identification of specific resistance mutations. Although these mutations are likely rare globally, they occur in approximately 7% of multidrug-resistant tuberculosis strains in some settings.

Whole-genome analysis of mycobacteria from birds at the San Diego Zoo

Methods

Mycobacteria isolated from more than 100 birds diagnosed with avian mycobacteriosis at the San Diego Zoo and its Safari Park were cultured postmortem and had their whole genomes sequenced. Computational workflows were developed and applied to identify the mycobacterial species in each DNA sample, to find single-nucleotide polymorphisms (SNPs) between samples of the same species, to further differentiate SNPs between as many as three different genotypes within a single sample, and to identify which samples are closely clustered genomically.

Results

Nine species of mycobacteria were found in 123 samples from 105 birds. The most common species were Mycobacterium avium and Mycobacterium genavense, which were in 49 and 48 birds, respectively. Most birds contained only a single mycobacterial species, but two birds contained a mixture of two species. The M. avium samples represent diverse strains of M. avium avium and M. avium hominissuis, with many pairs of samples differing by hundreds or thousands of SNPs across their common genome. By contrast, the M. genavense samples are much closer genomically; samples from 46 of 48 birds differ from each other by less than 110 SNPs. Some birds contained two, three, or even four genotypes of the same bacterial species. Such infections were found in 4 of 49 birds (8%) with M. avium and in 11 of 48 birds (23%) with M. genavense. Most were mixed infections, in which the bird was infected by multiple mycobacterial strains, but three infections with two genotypes differing by ≤ 10 SNPs were likely the result of within-host evolution. The samples from 31 birds with M. avium can be grouped into nine clusters within which any sample is ≤ 12 SNPs from at least one other sample in the cluster. Similarly, the samples from 40 birds with M. genavense can be grouped into ten such clusters. Information about these genomic clusters is being used in an ongoing, companion study of mycobacterial transmission to help inform management of bird collections.

Structure and variation of CRISPR and CRISPR-flanking regions in deleted-direct repeat region Mycobacterium tuberculosis complex strains

Background

CRISPR and CRISPR-flanking genomic regions are important for molecular epidemiology of Mycobacterium tuberculosis complex (MTBC) strains, and potentially for adaptive immunity to phage and plasmid DNA, and endogenous roles in the bacterium. Genotyping in the Israel National Mycobacterium Reference Center Tel-Aviv of over 1500 MTBC strains from 2008–2013 showed three strains with validated negative 43-spacer spoligotypes, that is, with putatively deleted direct repeat regions (deleted-DR/CRISPR regions). Two isolates of each of three negative spoligotype MTBC (a total of 6 isolates) were subjected to Next Generation Sequencing (NGS). As positive controls, NGS was performed for three intact-DR isolates belonging to T3_Eth, the largest multiple-drug-resistant (MDR)-containing African-origin cluster in Israel. Other controls consisted of NGS reads and complete whole genome sequences from GenBank for 20 intact-DR MTBC and for 1 deleted-DR MTBC strain recognized as CAS by its defining RD deletion.

Results

NGS reads from negative spoligotype MTBC mapped to reference H37Rv NC_000962.3 suggested that the DR/CRISPR regions were completely deleted except for retention of the middle IS6110 mobile element. Clonally specific deletion of CRISPR-flanking genes also was observed, including deletion of at least cas2 and cas1 genes. Genomic RD deletions defined lineages corresponding to the major spoligotype families Beijing, EAI, and Haarlem, consistent with 24 loci MIRU-VNTR profiles. Analysis of NGS reads, and analysis of contigs obtained by manual PCR confirmed that all 43 gold standard DR/CRISPR spacers were missing in the deleted-DR genomes.

Conclusions

Although many negative spoligotype strains are recorded as spoligotype-international-type (SIT) 2669 in the SITVIT international database, this is the first time to our knowledge that it has been shown that negative spoligotype strains are found in at least 4 different 24 loci MIRU-VNTR and RD deletion families. We report for the first time negative spoligotype-associated total loss of CRISPR region spacers and repeats, with accompanying clonally specific loss of flanking genes, including at least CRISPR-associated genes cas2 and cas1. Since cas1 deleted E.coli shows increased sensitivity to DNA damage and impaired chromosomal segregation, we discussed the possibility of a similar phenotype in the deleted-DR strains and Beijing family strains as both lack the cas1 gene.

Electronic supplementary material

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

The Evolution of Strain Typing in the Mycobacterium tuberculosis Complex

Tuberculosis (TB) is a contagious disease with a complex epidemiology. Therefore, molecular typing (genotyping) of Mycobacterium tuberculosis complex (MTBC) strains is of primary importance to effectively guide outbreak investigations, define transmission dynamics and assist global epidemiological surveillance of the disease. Large-scale genotyping is also needed to get better insights into the biological diversity and the evolution of the pathogen. Thanks to its shorter turnaround and simple numerical nomenclature system, mycobacterial interspersed repetitive unit-variable-number tandem repeat (MIRU-VNTR) typing, based on 24 standardized plus 4 hypervariable loci, optionally combined with spoligotyping, has replaced IS6110 DNA fingerprinting over the last decade as a gold standard among classical strain typing methods for many applications. With the continuous progress and decreasing costs of next-generation sequencing (NGS) technologies, typing based on whole genome sequencing (WGS) is now increasingly performed for near complete exploitation of the available genetic information. However, some important challenges remain such as the lack of standardization of WGS analysis pipelines, the need of databases for sharing WGS data at a global level, and a better understanding of the relevant genomic distances for defining clusters of recent TB transmission in different epidemiological contexts. This chapter provides an overview of the evolution of genotyping methods over the last three decades, which culminated with the development of WGS-based methods. It addresses the relative advantages and limitations of these techniques, indicates current challenges and potential directions for facilitating standardization of WGS-based typing, and provides suggestions on what method to use depending on the specific research question.

Molecular Strain Typing of Mycobacterium tuberculosis: a Review of Frequently Used Methods

Tuberculosis, caused by the bacterium Mycobacterium tuberculosis, remains one of the most serious global health problems. Molecular typing of M. tuberculosis has been used for various epidemiologic purposes as well as for clinical management. Currently, many techniques are available to type M. tuberculosis. Choosing the most appropriate technique in accordance with the existing laboratory conditions and the specific features of the geographic region is important. Insertion sequence IS6110-based restriction fragment length polymorphism (RFLP) analysis is considered the gold standard for the molecular epidemiologic investigations of tuberculosis. However, other polymerase chain reaction-based methods such as spacer oligonucleotide typing (spoligotyping), which detects 43 spacer sequence-interspersing direct repeats (DRs) in the genomic DR region; mycobacterial interspersed repetitive units-variable number tandem repeats, (MIRU-VNTR), which determines the number and size of tandem repetitive DNA sequences; repetitive-sequence-based PCR (rep-PCR), which provides high-throughput genotypic fingerprinting of multiple Mycobacterium species; and the recently developed genome-based whole genome sequencing methods demonstrate similar discriminatory power and greater convenience. This review focuses on techniques frequently used for the molecular typing of M. tuberculosis and discusses their general aspects and applications.

Methodological and Clinical Aspects of the Molecular Epidemiology of Mycobacterium tuberculosis and Other Mycobacteria

Molecular typing has revolutionized epidemiological studies of infectious diseases, including those of a mycobacterial etiology. With the advent of fingerprinting techniques, many traditional concepts regarding transmission, infectivity, or pathogenicity of mycobacterial bacilli have been revisited, and their conventional interpretations have been challenged. Since the mid-1990s, when the first typing methods were introduced, a plethora of other modalities have been proposed. So-called molecular epidemiology has become an essential subdiscipline of modern mycobacteriology. It serves as a resource for understanding the key issues in the epidemiology of tuberculosis and other mycobacterial diseases. Among these issues are disclosing sources of infection, quantifying recent transmission, identifying transmission links, discerning reinfection from relapse, tracking the geographic distribution and clonal expansion of specific strains, and exploring the genetic mechanisms underlying specific phenotypic traits, including virulence, organ tropism, transmissibility, or drug resistance. Since genotyping continues to unravel the biology of mycobacteria, it offers enormous promise in the fight against and prevention of the diseases caused by these pathogens. In this review, molecular typing methods for Mycobacterium tuberculosis and nontuberculous mycobacteria elaborated over the last 2 decades are summarized. The relevance of these methods to the epidemiological investigation, diagnosis, evolution, and control of mycobacterial diseases is discussed.

The genomic basis of adaptation to the fitness cost of rifampicin resistance in Pseudomonas aeruginosa

Antibiotic resistance carries a fitness cost that must be overcome in order for resistance to persist over the long term. Compensatory mutations that recover the functional defects associated with resistance mutations have been argued to play a key role in overcoming the cost of resistance, but compensatory mutations are expected to be rare relative to generally beneficial mutations that increase fitness, irrespective of antibiotic resistance. Given this asymmetry, population genetics theory predicts that populations should adapt by compensatory mutations when the cost of resistance is large, whereas generally beneficial mutations should drive adaptation when the cost of resistance is small. We tested this prediction by determining the genomic mechanisms underpinning adaptation to antibiotic-free conditions in populations of the pathogenic bacterium Pseudomonas aeruginosa that carry costly antibiotic resistance mutations. Whole-genome sequencing revealed that populations founded by high-cost rifampicin-resistant mutants adapted via compensatory mutations in three genes of the RNA polymerase core enzyme, whereas populations founded by low-cost mutants adapted by generally beneficial mutations, predominantly in the quorum-sensing transcriptional regulator gene lasR. Even though the importance of compensatory evolution in maintaining resistance has been widely recognized, our study shows that the roles of general adaptation in maintaining resistance should not be underestimated and highlights the need to understand how selection at other sites in the genome influences the dynamics of resistance alleles in clinical settings.

Rapid Drug Susceptibility Testing of Drug-Resistant Mycobacterium tuberculosis Isolates Directly from Clinical Samples by Use of Amplicon Sequencing: a Proof-of-Concept Study

Increasingly complex drug-resistant tuberculosis (DR-TB) is a major global health concern and one of the primary reasons why TB is now the leading infectious cause of death worldwide. Rapid characterization of a DR-TB patient's complete drug resistance profile would facilitate individualized treatment in place of empirical treatment, improve treatment outcomes, prevent amplification of resistance, and reduce the transmission of DR-TB. The use of targeted next-generation sequencing (NGS) to obtain drug resistance profiles directly from patient sputum samples has the potential to enable comprehensive evidence-based treatment plans to be implemented quickly, rather than in weeks to months, which is currently needed for phenotypic drug susceptibility testing (DST) results. In this pilot study, we evaluated the performance of amplicon sequencing of Mycobacterium tuberculosis DNA from patient sputum samples using a tabletop NGS technology and automated data analysis to provide a rapid DST solution (the Next Gen-RDST assay). One hundred sixty-six out of 176 (94.3%) sputum samples from the Republic of Moldova yielded complete Next Gen-RDST assay profiles for 7 drugs of interest. We found a high level of concordance of our Next Gen-RDST assay results with phenotypic DST (97.0%) and pyrosequencing (97.8%) results from the same clinical samples. Our Next Gen-RDST assay was also able to estimate the proportion of resistant-to-wild-type alleles down to mixtures of ≤1%, which demonstrates the ability to detect very low levels of resistant variants not detected by pyrosequencing and possibly below the threshold for phenotypic growth methods. The assay as described here could be used as a clinical or surveillance tool.

Does Choice Matter? Reference-Based Alignment for Molecular Epidemiology of Tuberculosis

When using genome sequencing for molecular epidemiology, short sequence reads are aligned to an arbitrary reference strain to detect single nucleotide polymorphisms. We investigated whether reference genome selection influences epidemiological inferences of Mycobacterium tuberculosis transmission by aligning sequence reads from 162 closely related lineage 4 (Euro-American) isolates to 7 different genomes. Phylogenetic trees were consistent with use of all but the most divergent genomes, suggesting that reference choice can be based on considerations other than M. tuberculosis lineage.

Interpreting whole genome sequencing for investigating tuberculosis transmission: a systematic review

BACKGROUND

Whole genome sequencing (WGS) is becoming an important part of epidemiological investigations of infectious diseases due to greater resolution and cost reductions compared to traditional typing approaches. Many public health and clinical teams will increasingly use WGS to investigate clusters of potential pathogen transmission, making it crucial to understand the benefits and assumptions of the analytical methods for investigating the data. We aimed to understand how different approaches affect inferences of transmission dynamics and outline limitations of the methods.

METHODS

We comprehensively searched electronic databases for studies that presented methods used to interpret WGS data for investigating tuberculosis (TB) transmission. Two authors independently selected studies for inclusion and extracted data. Due to considerable methodological heterogeneity between studies, we present summary data with accompanying narrative synthesis rather than pooled analyses.

RESULTS

Twenty-five studies met our inclusion criteria. Despite the range of interpretation tools, the usefulness of WGS data in understanding TB transmission often depends on the amount of genetic diversity in the setting. Where diversity is small, distinguishing re-infections from relapses may be impossible; interpretation may be aided by the use of epidemiological data, examining minor variants and deep sequencing. Conversely, when within-host diversity is large, due to genetic hitchhiking or co-infection of two dissimilar strains, it is critical to understand how it arose. Greater understanding of microevolution and mixed infection will enhance interpretation of WGS data.

CONCLUSIONS

As sequencing studies have sampled more intensely and integrated multiple sources of information, the understanding of TB transmission and diversity has grown, but there is still much to be learnt about the origins of diversity that will affect inferences from these data. Public health teams and researchers should combine epidemiological, clinical and WGS data to strengthen investigations of transmission.

TGS-TB: Total Genotyping Solution for Mycobacterium tuberculosis Using Short-Read Whole-Genome Sequencing

Whole-genome sequencing (WGS) with next-generation DNA sequencing (NGS) is an increasingly accessible and affordable method for genotyping hundreds of Mycobacterium tuberculosis (Mtb) isolates, leading to more effective epidemiological studies involving single nucleotide variations (SNVs) in core genomic sequences based on molecular evolution. We developed an all-in-one web-based tool for genotyping Mtb, referred to as the Total Genotyping Solution for TB (TGS-TB), to facilitate multiple genotyping platforms using NGS for spoligotyping and the detection of phylogenies with core genomic SNVs, IS6110 insertion sites, and 43 customized loci for variable number tandem repeat (VNTR) through a user-friendly, simple click interface. This methodology is implemented with a KvarQ script to predict MTBC lineages/sublineages and potential antimicrobial resistance. Seven Mtb isolates (JP01 to JP07) in this study showing the same VNTR profile were accurately discriminated through median-joining network analysis using SNVs unique to those isolates. An additional IS6110 insertion was detected in one of those isolates as supportive genetic information in addition to core genomic SNVs. The results of in silico analyses using TGS-TB are consistent with those obtained using conventional molecular genotyping methods, suggesting that NGS short reads could provide multiple genotypes to discriminate multiple strains of Mtb, although longer NGS reads (≥300-mer) will be required for full genotyping on the TGS-TB web site. Most available short reads (~100-mer) can be utilized to discriminate the isolates based on the core genome phylogeny. TGS-TB provides a more accurate and discriminative strain typing for clinical and epidemiological investigations; NGS strain typing offers a total genotyping solution for Mtb outbreak and surveillance. TGS-TB web site: https://gph.niid.go.jp/tgs-tb/.

Genomics and Machine Learning for Taxonomy Consensus: The Mycobacterium tuberculosis Complex Paradigm

Infra-species taxonomy is a prerequisite to compare features such as virulence in different pathogen lineages. Mycobacterium tuberculosis complex taxonomy has rapidly evolved in the last 20 years through intensive clinical isolation, advances in sequencing and in the description of fast-evolving loci (CRISPR and MIRU-VNTR). On-line tools to describe new isolates have been set up based on known diversity either on CRISPRs (also known as spoligotypes) or on MIRU-VNTR profiles. The underlying taxonomies are largely concordant but use different names and offer different depths. The objectives of this study were 1) to explicit the consensus that exists between the alternative taxonomies, and 2) to provide an on-line tool to ease classification of new isolates. Genotyping (24-VNTR, 43-spacers spoligotypes, IS6110-RFLP) was undertaken for 3,454 clinical isolates from the Netherlands (2004-2008). The resulting database was enlarged with African isolates to include most human tuberculosis diversity. Assignations were obtained using TB-Lineage, MIRU-VNTRPlus, SITVITWEB and an algorithm from Borile et al. By identifying the recurrent concordances between the alternative taxonomies, we proposed a consensus including 22 sublineages. Original and consensus assignations of the all isolates from the database were subsequently implemented into an ensemble learning approach based on Machine Learning tool Weka to derive a classification scheme. All assignations were reproduced with very good sensibilities and specificities. When applied to independent datasets, it was able to suggest new sublineages such as pseudo-Beijing. This Lineage Prediction tool, efficient on 15-MIRU, 24-VNTR and spoligotype data is available on the web interface “TBminer.” Another section of this website helps summarizing key molecular epidemiological data, easing tuberculosis surveillance. Altogether, we successfully used Machine Learning on a large dataset to set up and make available the first consensual taxonomy for human Mycobacterium tuberculosis complex. Additional developments using SNPs will help stabilizing it.

Genetic markers, genotyping methods & next generation sequencing in Mycobacterium tuberculosis

Molecular epidemiology (ME) is one of the main areas in tuberculosis research which is widely used to study the transmission epidemics and outbreaks of tubercle bacilli. It exploits the presence of various polymorphisms in the genome of the bacteria that can be widely used as genetic markers. Many DNA typing methods apply these genetic markers to differentiate various strains and to study the evolutionary relationships between them. The three widely used genotyping tools to differentiate Mycobacterium tuberculosis strains are IS6110 restriction fragment length polymorphism (RFLP), spacer oligotyping (Spoligotyping), and mycobacterial interspersed repeat units - variable number of tandem repeats (MIRU-VNTR). A new prospect towards ME was introduced with the development of whole genome sequencing (WGS) and the next generation sequencing (NGS) methods, where the entire genome is sequenced that not only helps in pointing out minute differences between the various sequences but also saves time and the cost. NGS is also found to be useful in identifying single nucleotide polymorphisms (SNPs), comparative genomics and also various aspects about transmission dynamics. These techniques enable the identification of mycobacterial strains and also facilitate the study of their phylogenetic and evolutionary traits.

A microbiological revolution meets an ancient disease: improving the management of tuberculosis with genomics

Tuberculosis (TB) is an ancient disease with an enormous global impact. Despite declining global incidence, the diagnosis, phenotyping, and epidemiological investigation of TB require significant clinical microbiology laboratory resources. Current methods for the detection and characterization of Mycobacterium tuberculosis consist of a series of laboratory tests varying in speed and performance, each of which yields incremental information about the disease. Since the sequencing of the first M. tuberculosis genome in 1998, genomic tools have aided in the diagnosis, treatment, and control of TB. Here we summarize genomics-based methods that are positioned to be introduced in the modern clinical TB laboratory, and we highlight how recent advances in genomics will improve the detection of antibiotic resistance-conferring mutations and the understanding of M. tuberculosis transmission dynamics and epidemiology. We imagine the future TB clinic as one that relies heavily on genomic interrogation of the M. tuberculosis isolate, allowing for more rapid diagnosis of TB and real-time monitoring of outbreak emergence.

Comparative whole-genome analysis of clinical isolates reveals characteristic architecture of Mycobacterium tuberculosis pangenome

The tubercle complex consists of closely related mycobacterium species which appear to be variants of a single species. Comparative genome analysis of different strains could provide useful clues and insights into the genetic diversity of the species. We integrated genome assemblies of 96 strains from Mycobacterium tuberculosis complex (MTBC), which included 8 Indian clinical isolates sequenced and assembled in this study, to understand its pangenome architecture. We predicted genes for all the 96 strains and clustered their respective CDSs into homologous gene clusters (HGCs) to reveal a hard-core, soft-core and accessory genome component of MTBC. The hard-core (HGCs shared amongst 100% of the strains) was comprised of 2,066 gene clusters whereas the soft-core (HGCs shared amongst at least 95% of the strains) comprised of 3,374 gene clusters. The change in the core and accessory genome components when observed as a function of their size revealed that MTBC has an open pangenome. We identified 74 HGCs that were absent from reference strains H37Rv and H37Ra but were present in most of clinical isolates. We report PCR validation on 9 candidate genes depicting 7 genes completely absent from H37Rv and H37Ra whereas 2 genes shared partial homology with them accounting to probable insertion and deletion events. The pangenome approach is a promising tool for studying strain specific genetic differences occurring within species. We also suggest that since selecting appropriate target genes for typing purposes requires the expected target gene be present in all isolates being typed, therefore estimating the core-component of the species becomes a subject of prime importance.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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