MIRU-profiler: a rapid tool for determination of 24-loci MIRU-VNTR profiles from assembled genomes of Mycobacterium tuberculosis.
ABSTRACT: Background:Tuberculosis (TB) resulted in an estimated 1.7 million deaths in the year 2016. The disease is caused by the members of Mycobacterium tuberculosis complex, which includes Mycobacterium tuberculosis, Mycobacterium bovis and other closely related TB causing organisms. In order to understand the epidemiological dynamics of TB, national TB control programs often conduct standardized genotyping at 24 Mycobacterial-Interspersed-Repetitive-Units (MIRU)-Variable-Number-of-Tandem-Repeats (VNTR) loci. With the advent of next generation sequencing technology, whole-genome sequencing (WGS) has been widely used for studying TB transmission. However, an open-source software that can connect WGS and MIRU-VNTR typing is currently unavailable, which hinders interlaboratory communication. In this manuscript, we introduce the MIRU-profiler program which could be used for prediction of MIRU-VNTR profile from WGS of M. tuberculosis. Implementation:The MIRU-profiler is implemented in shell scripting language and depends on EMBOSS software. The in-silico workflow of MIRU-profiler is similar to those described in the laboratory manuals for genotyping M. tuberculosis. Given an input genome sequence, the MIRU-profiler computes alleles at the standard 24-loci based on in-silico PCR amplicon lengths. The final output is a tab-delimited text file detailing the 24-loci MIRU-VNTR pattern of the input sequence. Validation:The MIRU-profiler was validated on four datasets: complete genomes from NCBI-GenBank (n = 11), complete genomes for locally isolated strains sequenced using PacBio (n = 4), complete genomes for BCG vaccine strains (n = 2) and draft genomes based on 250 bp paired-end Illumina reads (n = 106). Results:The digital MIRU-VNTR results were identical to the experimental genotyping results for complete genomes of locally isolated strains, BCG vaccine strains and five out of 11 genomes from the NCBI-GenBank. For draft genomes based on short Illumina reads, 21 out of 24 loci were inferred with a high accuracy, while a number of inaccuracies were recorded for three specific loci (ETRA, QUB11b and QUB26). One of the unique features of the MIRU-profiler was its ability to process multiple genomes in a batch. This feature was tested on all complete M. tuberculosis genome (n = 157), for which results were successfully obtained in approximately 14 min. Conclusion:The MIRU-profiler is a rapid tool for inference of digital MIRU-VNTR profile from the assembled genome sequences. The tool can accurately infer repeat numbers at the standard 24 or 21/24 MIRU-VNTR loci from the complete or draft genomes respectively. Thus, the tool is expected to bridge the communication gap between the laboratories using WGS and those using the conventional MIRU-VNTR typing.
Project description:BACKGROUND:Tuberculosis (TB) caused an estimated 1.4 million deaths and 10.4 million new cases globally in 2015. TB rates in the United States continue to steadily decline, yet rates in the State of Hawaii are perennially among the highest in the nation due to a continuous influx of immigrants from the Western Pacific and Asia. TB in Hawaii is composed of a unique distribution of genetic lineages, with the Beijing and Manila families of Mycobacterium tuberculosis (Mtb) comprising over two-thirds of TB cases. Standard fingerprinting methods (spoligotyping plus 24-loci Mycobacterial Interspersed Repetitive Units-Variable Number Tandem Repeats [MIRU-VNTR] fingerprinting) perform poorly when used to identify actual transmission clusters composed of isolates from these two families. Those typing methods typically group isolates from these families into large clusters of non-linked isolates with identical fingerprints. Next-generation whole-genome sequencing (WGS) provides a new tool for molecular epidemiology that can resolve clusters of isolates with identical spoligotyping and MIRU-VNTR fingerprints. METHODS:We performed WGS and SNP analysis and evaluated epidemiological data to investigate 19 apparent TB transmission clusters in Hawaii from 2003 to 2017 in order to assess WGS' ability to resolve putative Mtb clusters from the Beijing and Manila families. This project additionally investigated MIRU-VNTR allele prevalence to determine why standard Mtb fingerprinting fails to usefully distinguish actual transmission clusters from these two Mtb families. RESULTS:WGS excluded transmission events in seven of these putative clusters, confirmed transmission in eight, and identified both transmission-linked and non-linked isolates in four. For epidemiologically identified clusters, while the sensitivity of MIRU-VNTR fingerprinting for identifying actual transmission clusters was found to be 100%, its specificity was only 28.6% relative to WGS. We identified that the Beijing and Manila families' significantly lower Shannon evenness of MIRU-VNTR allele distributions than lineage 4 was the cause of standard fingerprinting's poor performance when identifying transmission in Beijing and Manila family clusters. CONCLUSIONS:This study demonstrated that WGS is necessary for epidemiological investigation of TB in Hawaii and the Pacific.
Project description: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.
Project description:<h4>Background</h4>Tuberculosis (TB) is a serious health problem in Tibet where Tibetans are the major ethnic group. Although genotyping of Mycobacterium tuberculosis (M. tuberculosis) isolates is a valuable tool for TB control, our knowledge of population structure of M. tuberculosis circulating in Tibet is limited.<h4>Methodology/principal findings</h4>In our study, a total of 576 M. tuberculosis isolates from Tibetans in Tibet, China, were analyzed via spoligotyping and 24-locus MIRU-VNTR. The Beijing genotype was the most prevalent family (90.63%, n = 522). Shared-type (ST) 1 was the most dominant genotype (88.89%, n = 512). We found that there was no association between the Beijing genotype and sex, age and treatment status. In this sample collection, 7 of the 24 MIRU-VNTR loci were highly or moderately discriminative according to their Hunter-Gaston discriminatory index. An informative set of 12 loci had similar discriminatory power with 24 loci set.<h4>Conclusions/significance</h4>The population structure of M. tuberculosis isolates in Tibetans is homogeneous and dominated by Beijing genotype. The analysis of 24-locus MIRU-VNTR data might be useful to select appropriate VNTR loci for the genotyping of M. tuberculosis.
Project description:Epidemiological contact tracing complemented with genotyping of clinical Mycobacterium tuberculosis isolates is important for understanding disease transmission. In Sweden, tuberculosis (TB) is mostly reported in migrant and homeless where epidemiologic contact tracing could pose a problem. This study compared epidemiologic linking with genotyping in a low burden country. Mycobacterium tuberculosis isolates (n?=?93) collected at Scania University Hospital in Southern Sweden were analysed with the standard genotyping method mycobacterial interspersed repetitive units-variable number tandem repeats (MIRU-VNTR) and the results were compared with whole genome sequencing (WGS). Using a maximum of twelve single nucleotide polymorphisms (SNPs) as the upper threshold of genomic relatedness noted among hosts, we identified 18 clusters with WGS comprising 52 patients with overall pairwise genetic maximum distances ranging from zero to nine SNPs. MIRU-VNTR and WGS clustered the same isolates, although the distribution differed depending on MIRU-VNTR limitations. Both genotyping techniques identified clusters where epidemiologic linking was insufficient, although WGS had higher correlation with epidemiologic data. To summarize, WGS provided better resolution of transmission than MIRU-VNTR in a setting with low TB incidence. WGS predicted epidemiologic links better which could consolidate and correct the epidemiologically linked cases, avoiding thus false clustering.
Project description:Yukon Territory (YT) is a remote region in northern Canada with ongoing spread of tuberculosis (TB). To explore the utility of whole genome sequencing (WGS) for TB surveillance and monitoring in a setting with detailed contact tracing and interview data, we used a mixed-methods approach. Our analysis included all culture-confirmed cases in YT (2005-2014) and incorporated data from 24-locus Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats (MIRU-VNTR) genotyping, WGS and contact tracing. We compared field-based (contact investigation (CI) data + MIRU-VNTR) and genomic-based (WGS + MIRU-VNTR + basic case data) investigations to identify the most likely source of each person's TB and assessed the knowledge, attitudes and practices of programme personnel around genotyping and genomics using online, multiple-choice surveys (n = 4) and an in-person group interview (n = 5). Field- and genomics-based approaches agreed for 26 of 32 (81%) cases on likely location of TB acquisition. There was less agreement in the identification of specific source cases (13/22 or 59% of cases). Single-locus MIRU-VNTR variants and limited genetic diversity complicated the analysis. Qualitative data indicated that participants viewed genomic epidemiology as a useful tool to streamline investigations, particularly in differentiating latent TB reactivation from the recent transmission. Based on this, genomic data could be used to enhance CIs, focus resources, target interventions and aid in TB programme evaluation.
Project description:In South Africa and other high prevalence countries, transmission is a significant contributor to rising rates of multidrug resistant tuberculosis (MDR-TB). Thus, there is a need to develop an early detection system for transmission clusters suitable for high burden settings. We have evaluated the discriminatory power and clustering concordance of a novel and simple genotyping approach, combining spoligotyping with pncA sequencing (SpoNC), against two well-established methods: IS6110-RFLP and 24-loci MIRU-VNTR.A total of 216 MDR-TB isolates collected from January to June 2010 from the NHLS Central TB referral laboratory in Braamfontein, Johannesburg, representing a diversity of strains from South Africa, were included. The isolates were submitted for genotyping, pncA sequencing and analysis to the Centre for Tuberculosis in South Africa and the Public Health Research Institute Tuberculosis Center at Rutgers University in the United States. Clustering rates, Hunter-Gaston Discriminatory Indexes (HGI) and Wallace coefficients were compared between the methods.Overall clustering rates were high by both IS6110-RFLP (52.8%) and MIRU-VNTR (45.8%), indicative of on-going transmission. Both 24-loci MIRU-VNTR and IS6110-RFLP had similar HGI (0.972 and 0.973, respectively), with close numbers of unique profiles (87 vs. 70), clustered isolates (129 vs. 146), and cluster sizes (2 to 26 vs. 2 to 25 isolates). Spoligotyping alone was the least discriminatory (80.1% clustering, HGI 0.903), with 28 unique types. However, the discriminatory power of spoligotyping was improved when combined with pncA sequencing using the SpoNC approach (61.8% clustering, HGI 0.958). A high proportion of MDR-TB isolates had mutations in pncA (68%, n = 145), and pncA mutations were significantly associated with clustering (p = 0.007 and p = 0.0013 by 24-loci MIRU-VNTR and IS6110-RFLP, respectively), suggesting high rates of resistance to pyrazinamide among all MDR-TB cases and particularly among clustered cases.We conclude that SpoNC provides good discrimination for MDR-TB surveillance and early identification of outbreaks in South Africa, with 24-loci MIRU-VNTR applied for pncA wild-type strains as needed.
Project description:Tuberculosis (TB) remains an immense public health problem in the Republic of Korea despite a more than fivefold decrease in the prevalence of the disease over the last 3 decades. The rise in drug-resistant TB has compounded the situation. We analyzed 208 clinical isolates of M. tuberculosis from the National Masan Tuberculosis Hospital by spoligotyping, IS6110 restriction fragment length polymorphism (RFLP), and 24-locus-based mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) typing to assess the diversity and transmission dynamics of the tubercle bacilli in the Republic of Korea. The majority of the isolates (97.1%) belonged to the Beijing genotype. Cluster analysis by MIRU-VNTR yielded a low clustering rate of 22.3%, with most of the clusters comprising isolates with diverse drug resistance patterns. The discriminatory capacity of the typing methods was high for RFLP and MIRU-VNTR (allelic diversity [h] = 0.99) but low for spoligotyping (h = 0.31). Although analysis of 19 MIRU-VNTR loci was needed to achieve maximum discrimination, an informative set of 8 loci (960, 1955, 2163b, 2165, 2996, 3192, 4052, and 4348) (h = 0.98) that was able to differentiate most of the closely related strains was identified. These findings suggest that 24-locus-based MIRU-VNTR typing is a likely suitable alternative to RFLP to differentiate clinical isolates in this setting, which is dominated by M. tuberculosis Beijing strains. Within the study limits, our results also suggest that the problem of drug-resistant TB in the Republic of Korea may be largely due to acquired resistance as opposed to transmission.
Project description:INTRODUCTION:Recurrent tuberculosis (TB) is one of the main challenges in TB control. Genotyping based on Mycobacterial Interspersed Repetitive Units-Variable Tandem Repeats (MIRU-VNTR) has been widely used to differentiate between relapse and reinfection, which are the two main causes of recurrent TB. There is a lack of data regarding the causes of TB recurrence in Georgia, and while differentiating between relapse and reinfection plays a key role in defining appropriate interventions, the required genotyping methodologies have not been implemented. The objective of this study was to implement MIRU-VNTR genotyping at the National Center for Tuberculosis and Lung Diseases (NCTBLD) and differentiate between relapse and reinfection in multidrug resistant (MDR-) TB patients from Tbilisi, Georgia. METHODS:Recurrent MDR tuberculosis cases from 2014-2016 diagnosed at NCTLD were included in the study when bacterial samples from both episodes were available. Genotyping based on the MIRU-VNTR 24 loci was implemented and used for differentiating between relapse and reinfection. Paired samples showing the same MIRU-VNTR pattern or one locus difference were classified as relapse, while two and more loci differences were treated as reinfection. Exact logistic regression was used to identify predictors of recurrence. RESULTS:Thirty two MDR-TB patients (64 samples) were included and MIRU-VNTR 24 typing was performed on the corresponding paired samples. Of the 32 patients, 25 (83.3%) were identified as relapse while 5 (16.7%) were due to re-infection. Patients with a history of incarceration were significantly associated with TB reinfection (p< 0.05). CONCLUSION:Recurrent TB in MDR patients in Georgia are mainly caused by relapse, raising concerns on the efficacy of the TB control program. An association between incarceration and reinfection likely reflects high levels of ongoing TB transmission in prisons, indicating the need for better TB infection control measures in these settings. Our results add to the rationale for implementing genotypic surveillance of TB more broadly to support TB control in Georgia.
Project description:RIFAQUIN was a tuberculosis chemotherapy trial in southern Africa including regimens with high-dose rifapentine with moxifloxacin. Here, the application of whole-genome sequencing (WGS) is evaluated within RIFAQUIN for identifying new infections in treated patients as either relapses or reinfections. WGS is further compared with mycobacterial interspersed repetitive units-variable number tandem repeats (MIRU-VNTR) typing. This is the first report of WGS being used to evaluate new infections in a completed clinical trial for which all treatment and epidemiological data are available for analysis.DNA from 36 paired samples of Mycobacterium tuberculosis cultured from patients before and after treatment was typed using 24-loci MIRU-VNTR, in silico spoligotyping and WGS. Following WGS, the sequences were mapped against the reference strain H37Rv, the single-nucleotide polymorphism (SNP) differences between pairs were identified, and a phylogenetic reconstruction was performed.WGS indicated that 32 of the paired samples had a very low number of SNP differences (0-5; likely relapses). One pair had an intermediate number of SNP differences, and was likely the result of a mixed infection with a pre-treatment minor genotype that was highly related to the post-treatment genotype; this was reclassified as a relapse, in contrast to the MIRU-VNTR result. The remaining three pairs had very high SNP differences (>750; likely reinfections).WGS and MIRU-VNTR both similarly differentiated relapses and reinfections, but WGS provided significant extra information. The low proportion of reinfections seen suggests that in standard chemotherapy trials with up to 24 months of follow-up, typing the strains brings little benefit to an analysis of the trial outcome in terms of differentiating relapse and reinfection. However, there is a benefit to using WGS as compared to MIRU-VNTR in terms of the additional genotype information obtained, in particular for defining the presence of mixed infections and the potential to identify known and novel drug-resistance markers.
Project description:<h4>Background</h4>There is a lack of information on the clinical characteristics of multidrug-resistant (MDR) tuberculosis (TB) and extensively drug-resistant (XDR) TB in the Jiangxi Province of China; furthermore, data have not been reported on the utility of mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) analyses in genotyping Mycobacterium tuberculosis strains isolated from this region. The aim of this study was to analyse the clinical features of patients with MDR and XDR TB from Jiangxi Province and to evaluate the discriminatory power of the 15-loci MIRU-VNTR method.<h4>Methods</h4>A retrospective study was conducted on patients diagnosed with MDR and XDR TB at the Jiangxi Chest Hospital from July 2010 to June 2011. The RD105 deletion-targeted multiplex PCR (DTM-PCR) and the 15-loci MIRU-VNTR method were used to determine the genetic background of the identified MDR and XDR M. tuberculosis clinical isolates.<h4>Results</h4>Of 804?M. tuberculosis clinical isolates, 159 (159/804, 19.8%) of the isolates were identified as MDR with first-line drug susceptibility testing. Of the 123 available MDR isolates, 13 (13/123, 10.6%) were XDR. The RD105 deletion-targeted multiplex PCR method identified 85 (85/110, 77.3%) MDR and 12 (12/13, 92.3%) XDR isolates as the Beijing genotype. MIRU-VNTR cluster analysis demonstrated that 101 MDR and 13 XDR strains had unique genotype patterns; the remaining 9 MDR strains were in 4 clusters, namely 1 cluster with 3 strains and 3 clusters with 2 strains, resulting in a low clustering rate (4.06%). The Hunter-Gaston discriminatory index (HGDI) of the 15-loci MIRU-VNTR method was as high as 0.992. In addition, clinical surveys showed that 87 (87/110, 79.1%) MDR TB patients and 10 (10/13, 76.9%) XDR TB patients had been previously treated. Diabetes mellitus was the most common comorbidity in both MDR TB (16/110, 14.5%) and XDR TB (2/13, 15.4%) patients.<h4>Conclusions</h4>Based on our preliminary data, the MDR and XDR M. tuberculosis clinical isolates identified at the Jiangxi Chest Hospital were genetically diverse and clustered at a low frequency. The 15-loci MIRU-VNTR method showed high discriminatory power and may be used as a first-line genotyping tool in investigating the molecular epidemiology of M. tuberculosis in Jiangxi, China. Decisive measures are urgently needed to effectively prevent and manage MDR and XDR tuberculosis in this province.