Tracking Mycobacterium tuberculosis sequencing samples using unique spikes of random DNA.

In the Netherlands, local laboratories are involved in the primary diagnosis of tuberculosis. Positive Mycobacterium tuberculosis complex cultures are sent to the National Institute for Public Health and the Environment (RIVM) for species identification, epidemiological typing, and screening for resistance by Whole Genome Sequencing (WGS). Occasional sample-swaps and cross-contaminations are known to occur in the diagnostic procedures. Such errors may lead to incorrect diagnoses resulting in the unnecessary or sub-optimal treatment of patients. Internal controls throughout the process ideally allow the early detection of such mistakes.

WGS more accurately predicts susceptibility of Mycobacterium tuberculosis to first-line drugs than phenotypic testing.

BACKGROUND: Drug-susceptibility testing (DST) of Mycobacterium tuberculosis complex (MTBC) isolates by the Mycobacteria Growth Indicator Tube (MGIT) approach is the most widely applied reference standard. However, the use of WGS is increasing in many developed countries to detect resistance and predict susceptibility. We investigated the reliability of WGS in predicting drug susceptibility, and analysed the discrepancies between WGS and MGIT against the first-line drugs rifampicin, isoniazid, ethambutol and pyrazinamide. METHODS: DST by MGIT and WGS was performed on MTBC isolates received in 2016/2017. Nine genes and/or their promotor regions were investigated for resistance-associated mutations: rpoB, katG, fabG1, ahpC, inhA, embA, embB, pncA and rpsA. Isolates that were discrepant in their MGIT/WGS results and a control group with concordant results were retested in the MGIT, at the critical concentration and a lower concentration, and incubated for up to 45 days after the control tube became positive in the MGIT. RESULTS: In total, 1136 isolates were included, of which 1121 were routine MTBC isolates from the Netherlands. The negative predictive value of WGS was ≥99.3% for all four first-line antibiotics. The majority of discrepancies for isoniazid and ethambutol were explained by growth at the lower concentrations, and for rifampicin by prolonged incubation in the MGIT, both indicating low-level resistance. CONCLUSIONS: Applying WGS in a country like the Netherlands, with a low TB incidence and low prevalence of resistance, can reduce the need for phenotypic DST for ∼90% of isolates and accurately detect mutations associated with low-level resistance, often missed in conventional DST.

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.

Correction: Epidemiological links between tuberculosis cases identified twice as efficiently by whole genome sequencing than conventional molecular typing: A population-based study.

[This corrects the article DOI: 10.1371/journal.pone.0195413.].

Epidemiological links between tuberculosis cases identified twice as efficiently by whole genome sequencing than conventional molecular typing: A population-based study.

BACKGROUND: Patients with Mycobacterium tuberculosis isolates sharing identical DNA fingerprint patterns can be epidemiologically linked. However, municipal health services in the Netherlands are able to confirm an epidemiological link in only around 23% of the patients with isolates clustered by the conventional variable number of tandem repeat (VNTR) genotyping. This research aims to investigate whether whole genome sequencing (WGS) is a more reliable predictor of epidemiological links between tuberculosis patients than VNTR genotyping. METHODS: VNTR genotyping and WGS were performed in parallel on all Mycobacterium tuberculosis complex isolates received at the Netherlands National Institute for Public Health and the Environment in 2016. Isolates were clustered by VNTR when they shared identical 24-loci VNTR patterns; isolates were assigned to a WGS cluster when the pair-wise genetic distance was ≤ 12 single nucleotide polymorphisms (SNPs). Cluster investigation was performed by municipal health services on all isolates clustered by VNTR in 2016. The proportion of epidemiological links identified among patients clustered by either method was calculated. RESULTS: In total, 535 isolates were genotyped, of which 25% (134/535) were clustered by VNTR and 14% (76/535) by WGS; the concordance between both typing methods was 86%. The proportion of epidemiological links among WGS clustered cases (57%) was twice as common than among VNTR clustered cases (31%). CONCLUSION: When WGS was applied, the number of clustered isolates was halved, while all epidemiologically linked cases remained clustered. WGS is therefore a more reliable tool to predict epidemiological links between tuberculosis cases than VNTR genotyping and will allow more efficient transmission tracing, as epidemiological investigations based on false clustering can be avoided.

A Predominant Variable-Number Tandem-Repeat Cluster of Mycobacterium tuberculosis Isolates among Asylum Seekers in the Netherlands and Denmark, Deciphered by Whole-Genome Sequencing.

In many countries, Mycobacterium tuberculosis isolates are routinely subjected to variable-number tandem-repeat (VNTR) typing to investigate M. tuberculosis transmission. Unexpectedly, cross-border clusters were identified among African refugees in the Netherlands and Denmark, although transmission in those countries was unlikely. Whole-genome sequencing (WGS) was applied to analyze transmission in depth and to assess the precision of VNTR typing. WGS was applied to 40 M. tuberculosis isolates from refugees in the Netherlands and Denmark (most of whom were from the Horn of Africa) that shared the exact same VNTR profile. Cluster investigations were undertaken to identify in-country epidemiological links. Combining WGS results for the isolates (all members of the central Asian strain [CAS]/Delhi genotype), from both European countries, an average genetic distance of 80 single-nucleotide polymorphisms (SNPs) (maximum, 153 SNPs) was observed. The few pairs of isolates with confirmed epidemiological links, except for one pair, had a maximum distance of 12 SNPs. WGS divided this refugee cluster into several subclusters of patients from the same country of origin. Although the M. tuberculosis cases, mainly originating from African countries, shared the exact same VNTR profile, most were clearly distinguished by WGS. The average genetic distance in this specific VNTR cluster was 2 times greater than that in other VNTR clusters. Thus, identical VNTR profiles did not represent recent direct M. tuberculosis transmission for this group of patients. It appears that either these strains from Africa are extremely conserved genetically or there is ongoing transmission of this genotype among refugees on their long migration routes from Africa to Europe.

Occurrence and Nature of Double Alleles in Variable-Number Tandem-Repeat Patterns of More than 8,000 Mycobacterium tuberculosis Complex Isolates in The Netherlands.

Since 2004, variable-number tandem-repeat (VNTR) typing of Mycobacterium tuberculosis complex isolates has been applied on a structural basis in The Netherlands to study the epidemiology of tuberculosis (TB). Although this technique is faster and technically less demanding than the previously used restriction fragment length polymorphism (RFLP) typing, reproducibility remains a concern. In the period from 2004 to 2015, 8,532 isolates were subjected to VNTR typing in The Netherlands, with 186 (2.2%) of these exhibiting double alleles at one locus. Double alleles were most common in loci 4052 and 2163b. The variables significantly associated with double alleles were urban living (odds ratio [OR], 1.503; 95% confidence interval [CI], 1.084 to 2.084; P = 0.014) and pulmonary TB (OR, 1.703; 95% CI, 1.216 to 2.386; P = 0.002). Single-colony cultures of double-allele strains were produced and revealed single-allele profiles; a maximum of five single nucleotide polymorphisms (SNPs) was observed between the single- and double-allele isolates from the same patient when whole-genome sequencing (WGS) was applied. This indicates the presence of two bacterial populations with slightly different VNTR profiles in the parental population, related to genetic drift. This observation is confirmed by the fact that secondary cases from TB source cases with double-allele isolates sometimes display only one of the two alleles present in the source case. Double alleles occur at a frequency of 2.2% in VNTR patterns in The Netherlands. They are caused by biological variation rather than by technical aberrations and can be transmitted either as single- or double-allele variants.

Parallel reaction monitoring of clinical Mycobacterium tuberculosis lineages reveals pre-existent markers of rifampicin tolerance in the emerging Beijing lineage.

The spread of multidrug resistant Mycobacterium tuberculosis is one of the major challenges in tuberculosis control. In Eurasia, the spread of multidrug resistant tuberculosis is driven by the M. tuberculosis Beijing genotype. In this study, we examined whether selective advantages are present in the proteome of Beijing isolates that contribute to the emergence of this genotype. To this end, we compared the proteome of M. tuberculosis Beijing to that of M. tuberculosis H37Rv, both in the presence and absence of the first-line antibiotic rifampicin. During rifampicin exposure, both M. tuberculosis genotypes express proteins belonging to the DosR dormancy regulon, which induces a metabolically hypoactive-, drug tolerant phenotype. However, these markers of rifampicin tolerance were already more abundant in the M. tuberculosis Beijing isolate prior to drug exposure. To determine whether the a priori high abundance of specific proteins contribute to the formation of antibiotic resistance in M. tuberculosis Beijing, we quantified the abundance of 33 selected proteins in 27 clinical isolates from the five most common M. tuberculosis lineages using parallel reaction monitoring. The observed pre-existing high abundance of dormancy proteins in Beijing strains provides an evolutionary advantage that allows these strains to persist for prolonged periods during rifampicin treatment. SIGNIFICANCE: M. tuberculosis is the leading cause of death by a bacterial infection worldwide. Treatment-regimen to eradicate this pathogen make use of the first-line antibiotic rifampicin, which is considered to be the cornerstone of modern day anti-tuberculosis treatment. Despite the potency of rifampicin, there is an increasing occurrence of rifampicin resistant mutants in a specific cluster of M. tuberculosis, the Beijing genotype. Using both a data dependent acquisition and a targeted proteomic approach we identified markers of rifampicin tolerance to be high abundant in members of the M. tuberculosis Beijing genotype, already prior drug exposure. The identification of this M. tuberculosis Beijing specific trait will contribute to improved diagnostics and treatment of M. tuberculosis.

Thioridazine Alters the Cell-Envelope Permeability of Mycobacterium tuberculosis.

The increasing occurrence of multidrug resistant tuberculosis exerts a major burden on treatment of this infectious disease. Thioridazine, previously used as a neuroleptic, is active against extensively drug resistant tuberculosis when added to other second- and third-line antibiotics. By quantitatively studying the proteome of thioridazine-treated Mycobacterium tuberculosis, we discovered the differential abundance of several proteins that are involved in the maintenance of the cell-envelope permeability barrier. By assessing the accumulation of fluorescent dyes in mycobacterial cells over time, we demonstrate that long-term drug exposure of M. tuberculosis indeed increased the cell-envelope permeability. The results of the current study demonstrate that thioridazine induced an increase in cell-envelope permeability and thereby the enhanced uptake of compounds. These results serve as a novel explanation to the previously reported synergistic effects between thioridazine and other antituberculosis drugs. This new insight in the working mechanism of this antituberculosis compound could open novel perspectives of future drug-administration regimens in combinational therapy.

Mechanisms of Phenotypic Rifampicin Tolerance in Mycobacterium tuberculosis Beijing Genotype Strain B0/W148 Revealed by Proteomics.

The "successful" Russian clone B0/W148 of Mycobacterium tuberculosis Beijing is well-known for its capacity to develop antibiotic resistance. During treatment, resistant mutants can occur that have inheritable resistance to specific antibiotics. Next to mutations, M. tuberculosis has several mechanisms that increase their tolerance to a variety of antibiotics. Insights in the phenotypic mechanisms that contribute to drug tolerance will increase our understanding of how antibiotic resistance develops in M. tuberculosis. In this study, we examined the (phospho)proteome dynamics in M. tuberculosis Beijing strain B0/W148 when exposed to a high dose of rifampicin; one of the most potent first-line antibiotics. A total of 2,534 proteins and 191 phosphorylation sites were identified, and revealed the differential regulation of DosR regulon proteins, which are necessary for the development of a dormant phenotype that is less susceptible to antibiotics. By examining independent phenotypic markers of dormancy, we show that persisters of in vitro rifampicin exposure entered a metabolically hypoactive state, which yields rifampicin and other antibiotics largely ineffective. These new insights in the role of protein regulation and post-translational modifications during the initial phase of rifampicin treatment reveal a shortcoming in the antituberculosis regimen that is administered to 8-9 million individuals annually.

Whole genome sequencing as the ultimate tool to diagnose tuberculosis.

In the past two decades, DNA techniques have been increasingly used in the laboratory diagnosis of tuberculosis (TB). The (sub) species of the Mycobacterium tuberculosis complex are usually identified using reverse line blot techniques. The resistance is predicted by the detection of mutations in genes associated with resistance. Nevertheless, all cases are still subjected to cumbersome phenotypic resistance testing. The production of a strain-characteristic DNA fingerprint, to investigate the epidemiology of TB, is done by the 24-locus variable number tandem repeat (VNTR) typing. However, most of the molecular techniques in the diagnosis of TB can eventually be replaced by whole genome sequencing (WGS). Many international TB reference laboratories are currently working on the introduction of WGS; however, standardization in the international context is lacking. The European Centre for Infectious Disease Prevention and Control in Stockholm, Sweden organizes a yearly round of quality control on VNTR typing and in 2015 for the first time also WGS. In this first proficiency study, only three out of eight international TB laboratories produced WGS results in line with those of the reference laboratory. The whole process of DNA isolation, purification, quantification, sequencing, and analysis/interpretation of data is still under development. In this presentation, many aspects will be covered that influence the quality and interpretation of WGS results. The turn-around-time, analysis, and utility of WGS will be discussed. Moreover, the experiences in the use of WGS in the molecular epidemiology of TB in The Netherlands are detailed. It can be concluded that many difficulties still have to be conquered. The state of the art is that bacteria still have to be cultured to have sufficient quality and quantity of DNA for succesful WGS. The quality of sequencing has improved significantly over the past 7years, and the detection of mutations has, therefore, become more reliable. The resistance mutations detected in WGS are in line with the ones visualized in reverse line blot techniques. The turnover in the genome of M. tuberculosis is very low, ∼0.3-0.5 mutations per genome per year. However, there is a wide variation in the occurrence of mutations per strain and genotype. Still, the resolution of WGS in epidemiological typing is higher than that in VNTR typing; previously suggested epidemiological links by VNTR typing are sometimes refuted on the basis of WGS. Although WGS offers the highest resolution in typing, in a country like The Netherlands, there are many strains with a limited genetic distance up to 100 mutations, without an apparent epidemiological link between the respective cases. These lookalikes are presumably even more prevalent in settings where predominant genotypes of M. tuberculosis are circulating. In summary, WGS seems to yield a more reliable prediction of resistance by the (lack of) detection of mutations in all 25 genes ever associated with resistance. This may within a short while prevent the need for many phenotypic resistance tests. Although more robust algorithms need to be developed, the recognition of the (sub) species in the M. tuberculosis complex seems possible. The first detailed studies on the population structure of M. tuberculosis strains in The Netherlands provide more resolution in typing but also an interesting observation that a part of the strains are genetically so conserved that they are separated by less than 100 mutations. This demands a more extended and accurate validation and understanding of the utility of WGS in the epidemiology of TB.

Optimization of standard in-house 24-locus variable-number tandem-repeat typing for Mycobacterium tuberculosis and its direct application to clinical material.

Variable-number tandem-repeat (VNTR) typing with a panel of 24 loci is the current gold standard in the molecular typing of Mycobacterium tuberculosis complex isolates. However, because of technical problems, a part of the loci often cannot be amplified by multiplex PCRs. Therefore, a considerable number of single-locus PCRs have to be performed for the loci with missing results, which impairs the laboratory work flow. Therefore, the original in-house method described by Supply et al. in 2006 was reevaluated. We modified seven primers and the PCR master mixture and obtained a strongly optimized in-house 24-locus VNTR typing method. The percentage of instantly complete 24-locus VNTR patterns detected in the routine flow of typing activities increased to 84.7% from the 72.3% obtained with the typing conducted with the commercially available Genoscreen MIRU-VNTR typing kit. The analytical sensitivity of the optimized in-house method was assessed by serial dilutions of M. tuberculosis in bronchoalveolar lavage fluid. A 1:10 dilution of the different strains tested was the lowest dilution for the detection of a complete 24-locus VNTR pattern. The optimized in-house 24-locus VNTR typing method will reduce the turnaround time of typing significantly and also the financial burden of these activities.

Current trends of Mycobacterium tuberculosis molecular epidemiology in Saudi Arabia and associated demographical factors.

Data are scarce on demographical factors related to the population structure of Mycobacterium tuberculosis in Saudi Arabia. A study was conducted on 902 clinical isolates to explore current trends in the phylogeography and associated demographical factors of tuberculosis by using spoligotyping and 24 loci based MIRU-VNTR typing. Young male patients (aged 16-29 and 30-44) were predominant in this cohort. The phylogenetic diversity among M. tuberculosis isolates was found high, as almost all known genetic lineages were identified. Delhi/CAS (26.4%), EAI (13.7%) and Haarlem (11.3%) were the most common lineages observed, particularly among the low age groups (16-29 and 30-44 years), whereas elderly patients (>60 years) showed a predominance in the lineages S, Ghana, TUR and Uganda-I. A statistically significant association was observed between gender of the patients and lineages of EAI (p value 0.026) and LAM (p value 0.005). Overall, molecular strain cluster rate was 34.4% with an elevated rate among patients aged below 15 years (43.1%), while cases among the elderly (>60 years) showed the lowest degree of clustering (12.5%). The largest level of clustering was noticed among cases caused by strains of the lineages Haarlem (59.8%), Beijing (55.8%) and LAM (42.8%). The current population structure of M. tuberculosis in Saudi Arabia is highly diverse with significant associations to demography, transmission dynamics and origin of the patients. The difference in genotype distributions among low and high aged patients reflects the ongoing change in the strain population structure in the country.

Comparative study of IS6110 restriction fragment length polymorphism and variable-number tandem-repeat typing of Mycobacterium tuberculosis isolates in the Netherlands, based on a 5-year nationwide survey.

In order to switch from IS6110 and polymorphic GC-rich repetitive sequence (PGRS) restriction fragment length polymorphism (RFLP) to 24-locus variable-number tandem-repeat (VNTR) typing of Mycobacterium tuberculosis complex isolates in the national tuberculosis control program in The Netherlands, a detailed evaluation on discriminatory power and agreement with findings in a cluster investigation was performed on 3,975 tuberculosis cases during the period of 2004 to 2008. The level of discrimination of the two typing methods did not differ substantially: RFLP typing yielded 2,733 distinct patterns compared to 2,607 in VNTR typing. The global concordance, defined as isolates labeled unique or identically distributed in clusters by both methods, amounted to 78.5% (n = 3,123). Of the remaining 855 cases, 12% (n = 479) of the cases were clustered only by VNTR, 7.7% (n = 305) only by RFLP typing, and 1.8% (n = 71) revealed different cluster compositions in the two approaches. A cluster investigation was performed for 87% (n = 1,462) of the cases clustered by RFLP. For the 740 cases with confirmed or presumed epidemiological links, 92% were concordant with VNTR typing. In contrast, only 64% of the 722 cases without an epidemiological link but clustered by RFLP typing were also clustered by VNTR typing. We conclude that VNTR typing has a discriminatory power equal to IS6110 RFLP typing but is in better agreement with findings in a cluster investigation performed on an RFLP-clustering-based cluster investigation. Both aspects make VNTR typing a suitable method for tuberculosis surveillance systems.

Characterization of Mycobacterium orygis as M. tuberculosis complex subspecies.

The oryx bacilli are Mycobacterium tuberculosis complex organisms for which phylogenetic position and host range are unsettled. We characterized 22 isolates by molecular methods and propose elevation to subspecies status as M. orygis. M. orygis is a causative agent of tuberculosis in animals and humans from Africa and South Asia.

Validation of pncA gene sequencing in combination with the mycobacterial growth indicator tube method to test susceptibility of Mycobacterium tuberculosis to pyrazinamide.

Pyrazinamide is important in the treatment of tuberculosis. Unfortunately, the diagnosis of pyrazinamide resistance is hampered by technical difficulties. We hypothesized that mutation analysis combined with the mycobacterial growth indicator tube (MGIT) phenotypic method would be a good predictor of pyrazinamide resistance. We prospectively analyzed 1,650 M. tuberculosis isolates referred to our tuberculosis reference laboratory in 2008 and 2009. In our laboratory, the MGIT 960 system was used for pyrazinamide resistance screening. If a pyrazinamide-resistant strain was detected, we performed a pncA gene mutation analysis. A second MGIT 960 susceptibility assay was performed afterwards to evaluate the accuracy of the pncA mutation analysis to detect true- or false-positive MGIT results. We observed pyrazinamide resistance in 69 samples using the first MGIT 960 analysis. In a second MGIT 960 analysis, 47 of the 69 samples proved susceptible (68% false positivity). Sensitivity of nonsynonymous pncA mutations for detecting resistant isolates was 73% (95% confidence interval [CI], 61% to 73%), and specificity was 100% (95% CI, 95% to 100%). A diagnostic algorithm incorporating phenotypic and molecular methods would have a 100% positive predictive value for detecting pyrazinamide-resistant isolates, indicating that such an algorithm, based on both methods, is a good predictor for pyrazinamide resistance in routine diagnostics.

[Establishment and application of a standard IS6110-RFLP method in the study of molecular genotyping analysis on Mycobacterium tuberculosis].

OBJECTIVE: To develop a standardized IS6110-restriction fragment length polymorphism (RFLP) method, used for evaluating the capacity of genotyping. METHODS: IS6110-RFLP of 78 Mycobacterium (M.) tuberculosis strains were studied by bio-molecular techniques including DNA isolation, PCR, restriction endonuclease enzyme analysis, southern blotting, agarose gel electrophoresis, together with data analysis by software Gel-Pro analyzer 3.1 and BioNumerics (Version 5.0). RESULTS: IS6110-RFLP method was established and standardized successfully, including DNA isolation, PCR, restriction endonuclease enzyme analysis, southern blotting, agarose gel electrophoresis and usage of the analysis software with standard parameters. By this method, 78 M. tuberculosis isolates were classified into 75 genotypes which belonged to 11 different clusters. Of all the isolates, 66.7% (52/78) belonged to a main cluster. CONCLUSION: Standard IS6110-RFLP method was established successfully. This method had powerful capacity for genotyping and strain level identification and could be used for the surveillance on pathogens of M. tuberculosis in China.