A single-step sequencing method for the identification of Mycobacterium tuberculosis complex species

Background

The Mycobacterium tuberculosis complex (MTC) comprises closely related species responsible for strictly human and zoonotic tuberculosis. Accurate species determination is useful for the identification of outbreaks and epidemiological links. Mycobacterium africanum and Mycobacterium canettii are typically restricted to Africa and M. bovis is a re-emerging pathogen. Identification of these species is difficult and expensive.

Methodology/Principal Findings

The Exact Tandem Repeat D (ETR-D; alias Mycobacterial Interspersed Repetitive Unit 4) was sequenced in MTC species type strains and 110 clinical isolates, in parallel to reference polyphasic identification based on phenotype profiling and sequencing of pncA, oxyR, hsp65, gyrB genes and the major polymorphism tandem repeat. Inclusion of M. tuberculosis isolates in the expanding, antibiotic-resistant Beijing clone was determined by Rv0927c gene sequencing. The ETR-D (780-bp) sequence unambiguously identified MTC species type strain except M. pinnipedii and M. microti thanks to six single nucleotide polymorphisms, variable numbers (1–7 copies) of the tandem repeat and two deletions/insertions. The ETR-D sequencing agreed with phenotypic identification in 107/110 clinical isolates and with reference polyphasic molecular identification in all isolates, comprising 98 M. tuberculosis, 5 M. bovis BCG type, 5 M. canettii, and 2 M. africanum. For M. tuberculosis isolates, the ETR-D sequence was not significantly associated with the Beijing clone.

Conclusions/Significance

ETR-D sequencing allowed accurate, single-step identification of the MTC at the species level. It circumvented the current expensive, time-consuming polyphasic approach. It could be used to depict epidemiology of zoonotic and human tuberculosis, especially in African countries where several MTC species are emerging.     

Distribution of spoligotyping defined genotypic lineages among drug-resistant Mycobacterium tuberculosis complex clinical isolates in Ankara, Turkey

Background

Investigation of genetic heterogeneity and spoligotype-defined lineages of drug-resistant Mycobacterium tuberculosis clinical isolates collected during a three-year period in two university hospitals and National Tuberculosis Reference and Research Laboratory in Ankara, Turkey.

Methods and Findings

A total of 95 drug-resistant M. tuberculosis isolates collected from three different centers were included in this study. Susceptibility testing of the isolates to four major antituberculous drugs was performed using proportion method on Löwenstein–Jensen medium and BACTEC 460-TB system. All clinical isolates were typed by using spoligotyping and IS6110-restriction fragment length polymorphism (RFLP) methods. Seventy-three of the 95 (76.8%) drug resistant M. tuberculosis isolates were isoniazid-resistant, 45 (47.4%) were rifampicin-resistant, 32 (33.7%) were streptomycin-resistant and 31 (32.6%) were ethambutol-resistant. The proportion of multidrug-resistant isolates (MDR) was 42.1%. By using spoligotyping, 35 distinct patterns were observed; 75 clinical isolates were grouped in 15 clusters (clustering rate of 79%) and 20 isolates displayed unique patterns. Five of these 20 unique patterns corresponded to orphan patterns in the SITVIT2 database, while 4 shared types containing 8 isolates were newly created. The most prevalent M. tuberculosis lineages were: Haarlem (23/95, 24.2%), ill-defined T superfamily (22/95, 23.2%), the Turkey family (19/95, 20%; previously designated as LAM7-TUR), Beijing (6/95, 6.3%), and Latin-America & Mediterranean (LAM, 5/95 or 5.3%), followed by Manu (3/95, 3.2%) and S (1/95, 1%) lineages. Four of the six Beijing family isolates (66.7%) were MDR. A combination of IS6110-RFLP and spoligotyping reduced the clustering rate from 79% to 11.5% among the drug resistant isolates.

Conclusions

The results obtained showed that ill-defined T, Haarlem, the Turkey family (previously designated as LAM7-TUR family with high phylogeographical specifity for Turkey), Beijing and LAM were predominant lineages observed in almost 80% of the drug-Resistant M. tuberculosis complex clinical isolates in Ankara, Turkey.     

Molecular-genetic method of identification of drug-resistant Mycobacterium tuberculosis

Above two million people annually die worldwide of tuberculosis. There has been a serious uptrend in the tuberculosis morbidity observed during the recent decade in Russia. The growing prevalence of drug-resistant strains of pathogens is highly dangerous. A timely diagnosis reduces significantly the risk of epidemics. The use of molecular methods cuts the time of tuberculosis diagnosis, including its resistant forms, from several weeks to several days. The presented survey contains a discussion of modem methods used for identification of Mycobacterium tuberculosis resistance to the main anti-TB drugs.     

人型结核杆菌全染色体DNA探针鉴别结核杆菌和其安分枝杆菌

The dot-blots containing DNA isolated from nonmycobacterial and mycobacterial microorganisms were hybridized with 32P-labeled M. tuberculosis whole chromosomal DNA at the various temperatures. The probe did not cross-hybridize to DNA of nonmycobacterial microorganisms (E. coli, Plasmid pUC19, Nocardia asteriodes), nor with DNA from all mycobacteria tested except M. bovis BCG under the higher temperature conditions. Microorganisms could also be directly spotted and lysed on nitrocellulose filters and used for hybridization thus making this technique suitable for clinical diagnosis.     

Genetic typing of Mycobacterium tuberculosis strains by spoligotyping and genome fingerprinting techniques

A total of 122 M. tuberculosis clinical drug-resistant strains isolated in Central Russia were studied by spoligotyping and genome fingerprinting techniques. According to spoligotyping results 77% of M. tuberculosis strains were distributed to 13 oligotypes, while 23% of these strains were found to form unique patterns. Most of them belonged to the families Beijing and Haarlem (43.4% and 13.9% respectively). The patterns of the strains of oligotype 12 (7F-7F-7E-0E-78-3E) were identical to those of the strains isolated in Brazil, France and the Netherlands. The strains of the spoligotype 22 (7F-1E-7F-7F-07-3E) had the patterns identical to those of the strains of group S13, also isolated in Brazil. According to genome fingerprinting 31.4% of the strains were found to belong to clusters with the similarity coefficient equal to 1. The strains belonging to genotypes W and A1 were found to prevail in the analyzed group.     

Nature of DNA polymorphism in the direct repeat cluster of Mycobacterium tuberculosis; application for strain differentiation by a novel typing method

Mycobacterium tuberculosis complex strains contain a unique chromosomal region, which consists of multiple 36bp direct repeats (DRs), which are interspersed by unique spacers 35 to 41 bp in length. In this study we investigated the nature of the DNA polymorphism of this DR cluster by sequencing part of this region in a large number of M. tuberculosis complex strains. Two types of genetic rearrangements were observed. One type consists of the variation in one or a few discrete, contiguous DRs plus spacer sequences. This variation is probably driven by homologous recombination between adjacent or distant DRs. The other type of polymorphism is probably driven by transpositional events of the insertion sequence, IS6110, which is almost invariably present in the DR cluster of M. tuberculosis complex strains. Based on the nature of the DNA polymorphism in the DR cluster, we developed a novel method of strain differentiation, direct variable repeat polymer chain reaction (DVR-PCR), which enables typing of individual M. tuberculosis strains in a single PCR. The method allows an excellent differentiation of epidemiologically unrelated isolates and, in principle, the DVR-PCR allows the detection of M. tuberculosis and strain differentiation at the same time.     

Biological properties of clinical isolates of Mycobacterium tuberculosis

The invasiveness of 2 grud sensitive and 8 holy-resistant M. tuberculosis clinical isolated was evaluated in experiments on BALB/c mice. Mycobacterial suspension was injected into the caudal vein of experimental animals. The results were evaluated by the degree of contamination of lungs and spleen of infected animals euthanized at different periods on time. The study revealed high variability in the degree of contamination of the organs of the animals infected with M. tuberculosis drug-resistant clinical isolates.     

Comparative genomic analysis of Mycobacterium tuberculosis clinical isolates

Background

Due to excessive antibiotic use, drug-resistant Mycobacterium tuberculosis has become a serious public health threat and a major obstacle to disease control in many countries. To better understand the evolution of drug-resistant M. tuberculosis strains, we performed whole genome sequencing for 7 M. tuberculosis clinical isolates with different antibiotic resistance profiles and conducted comparative genomic analysis of gene variations among them.

Results

We observed that all 7 M. tuberculosis clinical isolates with different levels of drug resistance harbored similar numbers of SNPs, ranging from 1409–1464. The numbers of insertion/deletions (Indels) identified in the 7 isolates were also similar, ranging from 56 to 101. A total of 39 types of mutations were identified in drug resistance-associated loci, including 14 previously reported ones and 25 newly identified ones. Sixteen of the identified large Indels spanned PE-PPE-PGRS genes, which represents a major source of antigenic variability. Aside from SNPs and Indels, a CRISPR locus with varied spacers was observed in all 7 clinical isolates, suggesting that they might play an important role in plasticity of the M. tuberculosis genome. The nucleotide diversity (Л value) and selection intensity (dN/dS value) of the whole genome sequences of the 7 isolates were similar. The dN/dS values were less than 1 for all 7 isolates (range from 0.608885 to 0.637365), supporting the notion that M. tuberculosis genomes undergo purifying selection. The Л values and dN/dS values were comparable between drug-susceptible and drug-resistant strains.

Conclusions

In this study, we show that clinical M. tuberculosis isolates exhibit distinct variations in terms of the distribution of SNP, Indels, CRISPR-cas locus, as well as the nucleotide diversity and selection intensity, but there are no generalizable differences between drug-susceptible and drug-resistant isolates on the genomic scale. Our study provides evidence strengthening the notion that the evolution of drug resistance among clinical M. tuberculosis isolates is clearly a complex and diversified process.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-469) contains supplementary material, which is available to authorized users.     

Molecular identification and characterization of rifampicin-resistant Mycobacterium tuberculosis isolates by line probe assay: an approach for rapid diagnosis of multidrug-resistant tuberculosis

Aim:  Early identification and characterization of rifampicin-resistant (R^r) Mycobacterium tuberculosis isolates recovered from the samples of tuberculosis (TB) patients in the Aegean (West Anatolian) Region was intended.
Methods and Results:  Sixty isolates [47 (78·3%) multidrug-resistant (MDR)], which were identified as M. tuberculosis complex and phenotypically resistant to rifampicin by both BACTEC mycobacteria growth indicator tube (MGIT) 960 and 460 systems were analysed by a commercial line probe assay (INNO-LiPA Rif TB). The concordance of LiPA with the in vitro susceptibility test was found as 98·3%. Among the isolates, S531L (R5 pattern; 46·7%) and L511P/R, S512T, Q513L/K (ΔS1 pattern; 11·7%) were the most frequent mutation patterns. As compared with the BACTEC systems and conventional techniques for cultivation, identification and in vitro susceptibility testing, INNO-LiPA Rif TB after cultivation in BACTEC MGIT 960 system provided an average of 20 days early diagnosis of R^r M. tuberculosis isolates.
Conclusions:  Rapid molecular identification and characterization of R^r M. tuberculosis isolates after BACTEC MGIT 960 cultivation would be useful for faster diagnosis, infection control and planning of accurate treatment in MDR-TB patients.
Significance and Impact of the Study:  Patients with MDR-TB need a specified treatment and efficient follow-up strategies. Rapid and practical methodologies to diagnose and follow these patients should be applied in routine use.     

A two-step strategy for molecular typing of multidrug-resistant Mycobacterium tuberculosis clinical isolates from Poland

Tuberculosis (TB) continues to be one of the most challenging public health problems in the world. An important contributor to the global burden of the disease is the emergence and spread of drug-resistant and particularly multidrug-resistant Mycobacterium tuberculosis strains (MDR), defined as being resistant to at least isoniazid and rifampicin. In recent years, the introduction of different DNA-based molecular typing methods has substantially improved the knowledge of the epidemiology of TB. The purpose of this study was to employ a combination of two PCR-based genotyping methods, namely spoligotyping and IS6110-Mtb1/Mtb2 PCR to investigate the clonal relatedness of MDR M. tuberculosis clinical isolates recovered from pulmonary TB patients from Poland. Among the 50 isolates examined, 28 (56%) were clustered by spoligotyping, whereas IS6110-Mtb1/Mtb2 PCR resulted in 16 (32%) clustered isolates. The isolates that clustered in both typing methods were assumed to be clonally related. A two-step strategy consisting of spoligotyping as a first-line test, performed on the entire pool of isolates, and IS6110-Mtb1/Mtb2 PCR typing as a confirmatory subtyping method, performed only within spoligotype-defined clusters, is an efficient approach for determining clonal relatedness among M. tuberculosis clinical isolates.     

A real-time PCR assay for detection of isoniazid resistance in Mycobacterium tuberculosis clinical isolates

A real-time PCR genotypic assay was developed for the detection of isoniazid (INH) resistance in Mycobacterium tuberculosis. The assay detects mutations C(-15)T and, possibly, G(-24)T in the regulatory region of the inhA gene and proved as sensitive and specific as nucleotide sequencing in all the clinical isolates tested. Our assays mapped the mutations efficiently in 10 out of 35 resistant isolates, thereby covering 29% of all resistant strains.     

Optimizing multiplex SNP-based data analysis for genotyping of Mycobacterium tuberculosis isolates

Background

Multiplex ligation-dependent probe amplification (MLPA) is a powerful tool to identify genomic polymorphisms. We have previously developed a single nucleotide polymorphism (SNP) and large sequence polymorphisms (LSP)-based MLPA assay using a read out on a liquid bead array to screen for 47 genetic markers in the Mycobacterium tuberculosis genome. In our assay we obtain information regarding the Mycobacterium tuberculosis lineage and drug resistance simultaneously. Previously we called the presence or absence of a genotypic marker based on a threshold signal level. Here we present a more elaborate data analysis method to standardize and streamline the interpretation of data generated by MLPA. The new data analysis method also identifies intermediate signals in addition to classification of signals as positive and negative. Intermediate calls can be informative with respect to identifying the simultaneous presence of sensitive and resistant alleles or infection with multiple different Mycobacterium tuberculosis strains.

Results

To validate our analysis method 100 DNA isolates of Mycobacterium tuberculosis extracted from cultured patient material collected at the National TB Reference Laboratory of the National Center for Tuberculosis and Lung Diseases in Tbilisi, Republic of Georgia were tested by MLPA. The data generated were interpreted blindly and then compared to results obtained by reference methods. MLPA profiles containing intermediate calls are flagged for expert review whereas the majority of profiles, not containing intermediate calls, were called automatically. No intermediate signals were identified in 74/100 isolates and in the remaining 26 isolates at least one genetic marker produced an intermediate signal.

Conclusion

Based on excellent agreement with the reference methods we conclude that the new data analysis method performed well. The streamlined data processing and standardized data interpretation allows the comparison of the Mycobacterium tuberculosis MLPA results between different experiments. All together this will facilitate the implementation of the MLPA assay in different settings.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-572) contains supplementary material, which is available to authorized users.     

The use of macroarrays for the identification of MDR Mycobacterium tuberculosis

The emergence of Mycobacterium tuberculosis (Mtb), resistant to both isoniazid (INH) and rifampicin (RIF) (MDR-TB), is an increasing threat to tuberculosis control programs. Susceptibility testing of Mtb complex isolates by phenotypic methods requires a minimum of 14 days from a primary specimen. This can be reduced significantly if molecular analysis is used. Low density oligonucleotide arrays (macroarrays) have been used successfully for the detection of RIF resistance in Mtb. We describe the use of macroarray technology to identify Mtb complex isolates resistant to INH and/or RIF. The macroarray MDR-Mtb screen has been designed to detect mutations in the RIF resistance determining region (RRDR) of Mtb rpoB and loci in katG and mabA-inhA associated with INH resistance. A panel of Mtb isolates containing 38 different RRDR genotypes, 4 different genotypes within codon 315 of katG and 2 genotypes at mabA-inhA was used to validate the macroarray. The wild type (WT) genotype was correctly identified at all three loci. Of the 37 mutant rpoB genotypes, 36 were correctly detected; the single mutant not detected contained a 9 base insertion. All mutations within katG and mabA-inhA were correctly identified. We conclude that this low cost, rapid system can usefully detect the mutations associated with the vast majority of MDR-Mtb.     

Molecular typing of the Mycobacterium tuberculosis strains circulating in central Russia: efficiency of spoligotyping and VNTR-MIRU

Current problems of molecular epidemiology of the Mycobacterium tuberculosis strains circulating in Samara Region, Russia are discussed. A total of 190 isolates of Mycobacterium tuberculosis were typed using two PCR-based molecular methods. The cultures were isolated from civil and prison patients with pulmonary tuberculosis recruited from different tuberculosis institutions across the Samara region. The usefulness of spoligotyping and 15-locii VNTR-MIRU was assessed for genotyping of Mycobacterium in population with high prevalence of Beijing strains (67.9%) using statistical analyses that included calculation of Hunter-Gaston index. The VNTR-MIRU method was demonstrated to be more efficient and was characterized by higher discrimination (index 0.747) compare to spoligotyping (index 0.572). VNTR-MIRU loci 10, 26, 31, 39, 40 and ETR-A were mostly polymorphic and therefore recommended for use in screening. It could be performed by manual electrophoresid, provided that automated sequencing is not available.     

A modified acid-fast staining method for rapid detection of Mycobacterium tuberculosis

A modified acid-fast staining method was developed for rapid detection of Mycobacterium tuberculosis and its L forms, wherein carbol fuchsin and dioxogen were mixed into the sputum smear. With this method, the dyeing time is shortened and heating is not required. The sensitivity, specificity, positive predictive value, negative predictive value, positive rate, and diagnostic efficiency of the new method were compared to those obtained by PCR using 50 clinical samples. Further, 468 clinical samples were analyzed using the new method, the modified intensified Kinyoun (IK) acid-fast staining method, and the traditional Ziehl-Neelsen acid-fast staining method. Differences among the positive detection rates of the three methods were analyzed using Student's t-test, and no significant differences were found between the new method and the modified IK acid-fast staining method, while the rates of both these methods were higher than that of the traditional acid-fast staining method. Additionally, the dyeing time in the new method was markedly less than that in the modified IK acid-fast staining method (5min and 24h, respectively).     

Genetic diversity of Iranian clinical isolates of Mycobacterium tuberculosis

In the current study we aimed to execute a rather less complicated molecular tying method, i.e. the random amplification of polymorphic DNA (RAPD) analysis to find the heterogeneity of Iranian strains of Mycobacterium tuberculosis. The isolates comprised a total of 96 strains of M. tuberculosis collected from clinical specimens of patients in Isfahan and Tehran. The isolates were assigned to the species M. tuberculosis by the key conventional and molecular methods. They were then subjected to RAPD analysis by four arbitrary primers, namely, the primers 27F, 1525R, MS- GF and INS-2. They were then evaluated for the number and intensity of the band patterns. The RAPD profiles of the Iranian isolates showed a degree of heterogeneity which varied based on the primer used. However, analysis of the isolates by primer INS-2 revealed the highest degree of diversity yielding 31 distinguishable RAPD types. RAPD analysis provides a rapid and easy means of identifying heterogeneity among the M. tuberculosis isolates. This typing system might be considered a valuable alternative molecular typing for countries with limited resources provided that the reproducibility and reliability of the method is carefully assured.     

Use of the insertion element IS6110 for DNA fingerprinting of Mycobacterium tuberculosis isolates presenting various profiles of drug susceptibility

Abstract IS 6100 is an insertion sequence of the IS3 family and it is present in multiple copies in the chromosome of Mycobacterium tuberculosis . Four to 15 copies are present in various strains of M. tuberculosis . In this study, the value of IS 6110 as an epidemiological marker of tuberculosis was examined. Unrelated clinical strains from Greek patients presented, in restriction fragment length polymorphism analysis, a high degree of polymorphism, whereas patterns of related clinical strains from familial outbreaks were identical. Since RFLP analysis with acetylaminofluorene labeled IS 6110 as the probe gave satisfactory results, it is suggested that this non-radioactive probe can be used in hospitals and health centres for the epidemiological survey of M. tuberculosis infections.     

Standardization of broth microdilution method for Mycobacterium tuberculosis

Indirect drug susceptibility tests of Mycobacterium tuberculosis was done to investigate the accuracy and feasibility of a broth microdilution method (BMM) for determining minimal inhibitory concentrations of conventional drugs against M. tuberculosis. Test drugs included isoniazid (H), rifampicin (R), ethambutol (E), streptomycin (S) and pyrazinamide (Z). Fifty isolates of M. tuberculosis from patients who had never received drug therapy, and H37Rv strain for control, were evaluated in the system. When comparing this method with the gold standard proportional method in Lowenstein-Jensen medium, sensitivity of 100% for all drugs and specifities of 91, 100, 96, 98 and 85% were observed respectively for H, R, E, S and Z. The BMM was read faster (14-20 days) than the proportional method (20-28 days). The microdilution method evaluated allows the testing of multiple drugs in multiple concentrations. It is easy to perform and does not require special equipment or expensive supplies. In contrast to radiometric method it does not use radioactive material.