Drug-resistance profiling and transmission dynamics of multidrug-resistant Mycobacterium tuberculosis in Saudi Arabia revealed by whole genome sequencing

Genetic Characterization and Population Structure of Drug-Resistant Mycobacterium tuberculosis Isolated from Brazilian Patients Using Whole-Genome Sequencing

The present study aimed to determine the genetic diversity of isolates of Mycobacterium tuberculosis (Mtb) from presumed drug-resistant tuberculosis patients from several states of Brazil. The isolates had been submitted to conventional drug susceptibility testing for first- and second-line drugs. Multidrug-resistant (MDR-TB) (54.8%) was the most frequent phenotypic resistance profile, in addition to an important high frequency of pre-extensive resistance (p-XDR-TB) (9.2%). Using whole-genome sequencing (WGS), we characterized 298 Mtb isolates from Brazil. Besides the analysis of genotype distribution and possible correlations between molecular and clinical data, we determined the performance of an in-house WGS pipeline with other online pipelines for Mtb lineages and drug resistance profile definitions. Sub-lineage 4.3 (52%) was the most frequent genotype, and the genomic approach revealed a p-XDR-TB level of 22.5%. We detected twenty novel mutations in three resistance genes, and six of these were observed in eight phenotypically resistant isolates. A cluster analysis of 170 isolates showed that 43.5% of the TB patients belonged to 24 genomic clusters, suggesting considerable ongoing transmission of DR-TB, including two interstate transmissions. The in-house WGS pipeline showed the best overall performance in drug resistance prediction, presenting the best accuracy values for five of the nine drugs tested. Significant associations were observed between suffering from fatal disease and genotypic p-XDR-TB (p = 0.03) and either phenotypic (p = 0.006) or genotypic (p = 0.0007) ethambutol resistance. The use of WGS analysis improved our understanding of the population structure of MTBC in Brazil and the genetic and clinical data correlations and demonstrated its utility for surveillance efforts regarding the spread of DR-TB, hopefully helping to avoid the emergence of even more resistant strains and to reduce TB incidence and mortality rates.

Mycobacterium tuberculosis Sub-Lineage 4.2.2/SIT149 as Dominant Drug-Resistant Clade in Northwest Ethiopia 2020-2022: In-silico Whole-Genome Sequence Analysis

Introduction

Drug resistance (DR) in Mycobacterium tuberculosis complex (MTBC) is mainly associated with certain lineages and varies across regions and countries. The Beijing genotype is the leading resistant lineage in Asia and western countries. M. tuberculosis (Mtb) (sub) lineages responsible for most drug resistance in Ethiopia are not well described. Hence, this study aimed to identify the leading drug resistance sub-lineages and characterize first-line anti-tuberculosis drug resistance-associated single nucleotide polymorphisms (SNPs).

Methods

A facility-based cross-sectional study was conducted in 2020-2022 among new and presumptive multidrug resistant-TB (MDR-TB) cases in Northwest Ethiopia. Whole-genome sequencing (WGS) was performed on 161 isolates using Illumina NovaSeq 6000 technology. The SNP mutations associated with drug resistance were identified using MtbSeq and TB profiler Bioinformatics softwares.

Results

Of the 146 Mtb isolates that were successfully genotyped, 20 (13.7%) harbored one or more resistance-associated SNPs. L4.2.2.ETH was the leading drug-resistant sub-lineage, accounting for 10/20 (50%) of the resistant Mtb. MDR-TB isolates showed extensive mutations against first-line anti-TB drugs. Ser450Leu/(tcg/tTg) for Rifampicin (RIF), Ser315Thr/(agc/aCc) for Isoniazid (INH), Met306Ile/(atg/atA(C)) for Ethambutol (EMB), and Gly69Asp for Streptomycin (STR) were the leading resistance associated mutations which accounted for 56.5%, 89.5%, 47%, and 29.4%, respectively. The presence of both clustered and non-clustered drug resistance (DR) isolates indicated that the epidemics is driven by both new DR development and acquired resistance.

Conclusion

The high prevalence of drug-resistant TB due to geographically restricted sub-lineages (L4.2.2.ETH) indicates the ongoing local micro epidemics. The Mtb drug resistance surveillance system must be improved. Further evolutionary analysis of L4.2.2.ETH strain is highly desirable to understand evolutionary forces that leads L4.2.2.ETH in to high level DR and transmissible sub-lineage.

Borderline rpoB mutations transmit at the same rate as common rpoB mutations in a tuberculosis cohort in Bangladesh

The spread of multidrug-resistant tuberculosis (MDR-TB) is a growing problem in many countries worldwide. Resistance to one of the primary first-line drugs, rifampicin, is caused by mutations in the Mycobacterium tuberculosis rpoB gene. So-called borderline rpoB mutations confer low-level resistance, in contrast to more common rpoB mutations which confer high-level resistance. While some borderline mutations show lower fitness in vitro than common mutations, their in vivo fitness is currently unknown. We used a dataset of 394 whole genome sequenced MDR-TB isolates from Bangladesh, representing around 44 % of notified MDR-TB cases over 6 years, to look at differences in transmission clustering between isolates with borderline rpoB mutations and those with common rpoB mutations. We found a relatively low percentage of transmission clustering in the dataset (34.8 %) but no difference in clustering between different types of rpoB mutations. Compensatory mutations in rpoA, rpoB, and rpoC were associated with higher levels of transmission clustering as were lineages two, three, and four relative to lineage one. Young people as well as patients with high sputum smear positive TB were more likely to be in a transmission cluster. Our findings show that although borderline rpoB mutations have lower in vitro growth potential this does not translate into lower transmission potential or in vivo fitness. Proper detection of these mutations is crucial to ensure they do not go unnoticed and spread MDR-TB within communities.

Molecular epidemiology and transmission dynamics of multi-drug resistant tuberculosis strains using whole genome sequencing in the Amhara region, Ethiopia

BACKGROUND

Drug resistant Mycobacterium tuberculosis prevention and care is a major challenge in Ethiopia. The World health organization has designated Ethiopia as one of the 30 high burden multi-drug resistant tuberculosis (MDR-TB) countries. There is limited information regarding genetic diversity and transmission dynamics of MDR-TB in Ethiopia.

OBJECTIVE

To investigate the molecular epidemiology and transmission dynamics of MDR-TB strains using whole genome sequence (WGS) in the Amhara region.

METHODS

Forty-five MDR-TB clinical isolates from Amhara region were collected between 2016 and 2018, and characterized using WGS and 24-loci Mycobacterium Interspersed Repetitive Units Variable Number of Tandem Repeats (MIRU-VNTR) typing. Clusters were defined based on the maximum distance of 12 single nucleotide polymorphisms (SNPs) or alleles as the upper threshold of genomic relatedness. Five or less SNPs or alleles distance or identical 24-loci VNTR typing is denoted as surrogate marker for recent transmission.

RESULTS

Forty-one of the 45 isolates were analyzed by WGS and 44% (18/41) of the isolates were distributed into 4 clusters. Of the 41 MDR-TB isolates, 58.5% were classified as lineage 4, 36.5% lineage 3 and 5% lineage 1. Overall, TUR genotype (54%) was the predominant in MDR-TB strains. 41% (17/41) of the isolates were clustered into four WGS groups and the remaining isolates were unique strains. The predominant cluster (Cluster 1) was composed of nine isolates belonging to lineage 4 and of these, four isolates were in the recent transmission links.

CONCLUSIONS

Majority of MDR-TB strain cluster and predominance of TUR lineage in the Amhara region give rise to concerns for possible ongoing transmission. Efforts to strengthen TB laboratory to advance diagnosis, intensified active case finding, and expanded contact tracing activities are needed in order to improve rapid diagnosis and initiate early treatment. This would lead to the interruption of the transmission chain and stop the spread of MDR-TB in the Amhara region.

Origin and Global Expansion of Mycobacterium tuberculosis Complex Lineage 3

Simple Summary

Tuberculosis still causes 1.5 million deaths annually and is mainly caused by Mycobacterium tuberculosis complex strains belonging to three evolutionary modern lineages (Lineages 2–4). While Lineage 2 and Lineage 4 virtually conquered the world, Lineage 3 is particularly successful in Northern and Eastern Africa, as well as in Southern Asia, the suspected evolutionary origin of these strains. Here, we sought to understand how Lineage 3 strains came to the African continent. To this end, we performed routine genotyping to characterize over 2500 clinical isolates from 38 countries. We then selected a representative collection of 373 isolates for a whole-genome analysis and a modeling approach to infer the geographic origin of different sublineages. In fact, the origin of Lineage 3 could be located in India, and we found evidence for independent introductions of four distinct sublineages into North/East Africa, in line with known ancient exchanges and migrations between both world regions. Our study illustrates that the evolutionary history of humans and their pathogens are closely connected and further provides a systematic understanding of the genomic diversity of Lineage 3, which could be important for the development of new tuberculosis vaccines or new therapeutics.

Abstract

Mycobacterium tuberculosis complex (MTBC) Lineage 3 (L3) strains are abundant in world regions with the highest tuberculosis burden. To investigate the population structure and the global diversity of this major lineage, we analyzed a dataset comprising 2682 L3 strains from 38 countries over 5 continents, by employing 24-loci mycobacterial interspersed repetitive unit-variable number of tandem repeats genotyping (MIRU-VNTR) and drug susceptibility testing. We further combined whole-genome sequencing (WGS) and phylogeographic analysis for 373 strains representing the global L3 genetic diversity. Ancestral state reconstruction confirmed that the origin of L3 strains is located in Southern Asia and further revealed multiple independent introduction events into North-East and East Africa. This study provides a systematic understanding of the global diversity of L3 strains and reports phylogenetic variations that could inform clinical trials which evaluate the effectivity of new drugs/regimens or vaccine candidates.

Emerging Status of Multidrug-Resistant Bacteria and Fungi in the Arabian Peninsula

Simple Summary

The incidence and developing status of multidrug-resistant bacteria and fungi, as well as their related mortality, is reviewed by a systematic published literature search from nine countries in the Arabian Peninsula. In order to analyse the emerging status and mortality, a total of 382 research articles were selected from a comprehensive screening of 1705 papers. More than 850 deaths reported since 2010 in the Arabian Peninsula due to the infection of multidrug-resistant bacteria and fungi. Multidrug-resistant bacteria Acinetobacter baumannii, Mycobacterium tuberculosis, Staphylococcus aureus, and fungi Candida auris are the most prevalent and causing high deaths. To control these infections and associated deaths in the Arabian Peninsula, continuous preventive measures, accurate methods for early diagnosis of infection, active surveillance, constant monitoring, developing vaccines, eradicating multidrug resistance modulators, and data sharing among countries are required.

Abstract

We aimed to identify the prevalence and emerging status of multidrug-resistant bacteria and fungi and their associated mortality in nine countries in the Arabian Peninsula. Original research articles and case studies regarding multidrug-resistant bacteria and fungi in the Arabian Peninsula, published during the last 10 years, were retrieved from PubMed and Scopus. A total of 382 studies were included as per the inclusion and exclusion criteria, as well as the PRISMA guidelines, from a thorough screening of 1705 articles, in order to analyse the emerging status and mortality. The emerging nature of >120 multidrug-resistant (MDR) bacteria and fungi in the Arabian Peninsula is a serious concern that requires continuous monitoring and immediate preventive measures. More than 50% (n = 453) of multidrug-resistant, microbe-associated mortality (n = 871) in the Arabian Peninsula was due to MDR Acinetobacter baumannii, Mycobacterium tuberculosis and Staphylococcus aureus infection. Overall, a 16.51% mortality was reported among MDR-infected patients in the Arabian Peninsula from the 382 articles of this registered systematic review. MDR A. baumannii (5600 isolates) prevailed in all the nine countries of the Arabian Peninsula and was one of the fastest emerging MDR bacteria with the highest mortality (n = 210). A total of 13,087 Mycobacterium tuberculosis isolates were reported in the region. Candida auris (580 strains) is the most prevalent among the MDR fungal pathogen in the Arabian Peninsula, having caused 54 mortalities. Active surveillance, constant monitoring, the development of a candidate vaccine, an early diagnosis of MDR infection, the elimination of multidrug resistance modulators and uninterrupted preventive measures with enhanced data sharing are mandatory to control MDR infection and associated diseases of the Arabian Peninsula. Accurate and rapid detection methods are needed to differentiate MDR strain from other strains of the species. This review summarises the logical relation, prevalence, emerging status and associated mortality of MDR microbes in the Arabian Peninsula.

A Saudi Arabian Public Health Perspective of Tuberculosis

Tuberculosis is a global health challenge due to its spreading potential. The Kingdom of Saudi Arabia (KSA) faces a challenge in the spread of tuberculosis from migrant workers, but the foremost threat is the huge number of pilgrims who travel to visit sacred sites of the Islamic world located in the holy cities of Makkah and Al Madina. Pilgrims visit throughout the year but especially in the months of Ramadan and Zul-Hijah. The rise of resistance in Mycobacterium tuberculosis is an established global phenomenon that makes such large congregations likely hotspots in the dissemination and spread of disease at a global level. Although very stringent and effective measures exist, the threat remains due to the ever-changing dynamics of this highly pathogenic disease. This overview primarily highlights the current public health challenges posed by this disease to the Saudi health system, which needs to be highlighted not only to the concerned authorities of KSA, but also to the concerned global quarters since the pilgrims and migrants come from all parts of the world with a majority coming from high tuberculosis-burdened countries.

Clusters of Drug-Resistant Mycobacterium tuberculosis Detected by Whole-Genome Sequence Analysis of Nationwide Sample, Thailand, 2014-2017

Multidrug-resistant tuberculosis (MDR TB), pre-extensively drug-resistant tuberculosis (pre-XDR TB), and extensively drug-resistant tuberculosis (XDR TB) complicate disease control. We analyzed whole-genome sequence data for 579 phenotypically drug-resistant M. tuberculosis isolates (28% of available MDR/pre-XDR and all culturable XDR TB isolates collected in Thailand during 2014–2017). Most isolates were from lineage 2 (n = 482; 83.2%). Cluster analysis revealed that 281/579 isolates (48.5%) formed 89 clusters, including 205 MDR TB, 46 pre-XDR TB, 19 XDR TB, and 11 poly–drug-resistant TB isolates based on genotypic drug resistance. Members of most clusters had the same subset of drug resistance-associated mutations, supporting potential primary resistance in MDR TB (n = 176/205; 85.9%), pre-XDR TB (n = 29/46; 63.0%), and XDR TB (n = 14/19; 73.7%). Thirteen major clades were significantly associated with geography (p<0.001). Clusters of clonal origin contribute greatly to the high prevalence of drug-resistant TB in Thailand.

Highly transmitted M. tuberculosis strains are more likely to evolve MDR/XDR and cause outbreaks, but what makes them highly transmitted?

Multi-Drug-Resistant strains of Mycobacterium tuberculosis (MDR-TB) are a serious obstacle to global TB eradication. While most MDR-TB strains are infrequently transmitted, a few cause large transmission clusters that contribute substantially to local MDR-TB burdens. Here we examine whether the known mutations in these strains can explain their success. Drug resistance mutations differ in fitness costs and strains can also acquire compensatory mutations (CM) to restore fitness, but some highly transmitted MDR strains have no CM. The acquisition of resistance mutations that maintain high transmissibility seems to occur by chance and are more likely in strains that are intrinsically highly transmitted and cause many cases. Modern Beijing lineage strains have caused several large outbreaks, but MDR outbreaks are also caused by ancient Beijing and lineage 4 strains, suggesting the lineage is less important than the characteristics of the individual strain. The development of fluoroquinolone resistance appears to represent another level of selection, in which strains must surmount unknown fitness costs of gyrA mutations. The genetic determinants of high transmission are poorly defined but may involve genes encoding proteins involved in molybdenum acquisition and the Esx systems. In addition, strains eliciting lower cytokine responses and producing more caseating granulomas may have advantages for transmission. Successful MDR/XDR strains generally evolve from highly transmitted drug sensitive parent strains due to selection pressures from deficiencies in local TB control programs. Until TB incidence is considerably reduced, there will likely be highly transmitted strains that develop resistance to any new antibiotic.

Drug resistance, fitness and compensatory mutations in Mycobacterium tuberculosis

For tuberculosis to be eradicated, the transmission of Multi-Drug-Resistant and eXtensively Drug Resistant strains of Mycobacterium tuberculosis (MDR and XDR-TB) must be considerably reduced. Drug resistant strains were initially thought to have reduced fitness, and the majority of resistant strains may actually have compromised fitness because they are found in only one or a few patients. In contrast, some MDR/XDR-TB strains are highly transmitted and cause large outbreaks. Most antibiotics target essential bacterial functions and the mutations that confer resistance to anti-TB drugs can incur fitness costs manifested as slower growth and reduced viability. The fitness costs vary with different resistance mutations and the bacilli can also accumulate secondary mutations that compensate for the compromised functions and partially or fully restore lost fitness. The compensatory mutations (CM) are different for each antibiotic, as they mitigate the deleterious effects of the specific functions compromised by the resistance mutations. CM are generally more common in strains with resistance mutations incurring the greatest fitness costs, but for RIF resistance, CM are most frequent in strains with the mutation carrying the least fitness cost, Ser450Leu. Here, we review what is known about fitness costs, CM and mechanisms of resistance to the drugs that define a strain as MDR or XDR-TB. The relative fitness costs of the resistance mutations and the mitigating effects of CM largely explain why certain mutations are frequently found in highly transmitted clusters while others are less frequently, rarely or never found in clinical isolates. The CM illustrate how drug resistance affects bacteria and how bacteria evolve to overcome the effects of the antibiotics, and thus a paradigm for how mycobacteria can evolve in response to stress.

Controlling the drug-resistant tuberculosis epidemic in India: challenges and implications

India has a higher tuberculosis (TB) burden than any other country, accounting for an estimated one-fourth of the global burden. Drug-resistant tuberculosis (DR-TB) presents a major public health problem in India. Patients with DR-TB often require profound changes in their drug regimens, which are invariably linked to poor treatment adherence and sub-optimal treatment outcomes compared to drug-sensitive TB. The challenge of addressing DR-TB is critical for India, as India contributes over 27% of global DR-TB cases. In recent decades, India has been proactive in its battle against TB, even implementing a revised National Strategic Plan to eliminate TB by 2025. However, to achieve this ambitious goal, the country will need to take a multifaceted approach with respect to its management of DR-TB. Despite concerted efforts made by the National TB Elimination Program, India faces substantial challenges with regard to DR-TB care, especially in peripheral and resource-limited endemic zones. This article describes some of the major challenges associated with mitigating the growing DR-TB epidemic in India and their implications.

Prevalence of Multidrug-Resistant Tuberculosis in Dalian, China: A Retrospective Study

Purpose

Multidrug-resistant tuberculosis (MDR-TB) is the cause of serious health and economic burdens worldwide. The present study aimed to explore the initial and acquired drug-resistance rates among TB patients from 2012 to 2019 in Dalian, China. The effectiveness of MDR-TB prevention and control strategies were then evaluated.

Patients and Methods

Drug susceptibility testing (DST) was performed for 6429 diagnosed, culture-positive, Mycobacterium tuberculosis (MTB) strains, including 4661 new cases and 1768 previously treated cases. Descriptive statistics were employed to calculate the frequencies and percentages of TB strains, and the average annual growth rates (AAGRs) for each strain were calculated. The Chi-square test was applied to examine the significance of linear drug-resistance trends over time during the study period.

Results

Over the eight-year study period, the percentages of both initial (from 9.01% to 4.82%) and acquired (from 40.85% to 9.09%) MDR-TB cases decreased significantly, AAGRs of 8.55% and 19.32%, respectively. Among new and previously treated TB patients, significant downtrends were observed for the rates of both initial and acquired MDR-TB among young and middle-aged individuals (P < 0.05). Additionally, among both new and previously treated TB patients, the percentages of individuals with drug resistance against isoniazid (INH), rifampicin (RFP), ofloxacin (OFX), and amikacin (AMK) decreased significantly (P < 0.05) from 2012 to 2019 in Dalian, China.

Conclusion

The initial and acquired multidrug resistance rates exhibited significantly decreasing trends from 2012 to 2019, suggesting that MDR-TB prevalence has been controlled effectively in Dalian, China. The MDR-TB epidemic was reversed in the short term by establishing feasible strategies for detection, diagnosis, treatment, and infection control.

Reporting practices for genomic epidemiology of tuberculosis: a systematic review of the literature using STROME-ID guidelines as a benchmark

BACKGROUND

Pathogen genomics have become increasingly important in infectious disease epidemiology and public health. The Strengthening the Reporting of Molecular Epidemiology for Infectious Diseases (STROME-ID) guidelines were developed to outline a minimum set of criteria that should be reported in genomic epidemiology studies to facilitate assessment of study quality. We evaluate such reporting practices, using tuberculosis as an example.

METHODS

For this systematic review, we initially searched MEDLINE, Embase Classic, and Embase on May 3, 2017, using the search terms "tuberculosis" and "genom* sequencing". We updated this initial search on April 23, 2019, and also included a search of bioRxiv at this time. We included studies in English, French, or Spanish that recruited patients with microbiologically confirmed tuberculosis and used whole genome sequencing for typing of strains. Non-human studies, conference abstracts, and literature reviews were excluded. For each included study, the number and proportion of fulfilled STROME-ID criteria were recorded by two reviewers. A comparison of the mean proportion of fulfilled STROME-ID criteria before and after publication of the STROME-ID guidelines (in 2014) was done using a two-tailed t test. Quasi-Poisson regression and tobit regression were used to examine associations between study characteristics and the number and proportion of fulfilled STROME-ID criteria. This study was registered with PROSPERO, CRD42017064395.

FINDINGS

976 titles and abstracts were identified by our primary search, with an additional 16 studies identified in bioRxiv. 114 full texts (published between 2009 and 2019) were eligible for inclusion. The mean proportion of STROME-ID criteria fulfilled was 50% (SD 12; range 16-75). The proportion of criteria fulfilled was similar before and after STROME-ID publication (51% [SD 11] vs 46% [14], p=0·26). The number of criteria reported (among those applicable to all studies) was not associated with impact factor, h-index, country of affiliation of senior author, or sample size of isolates. Similarly, the proportion of criteria fulfilled was not associated with these characteristics, with the exception of a sample size of isolates of 277 or more (the highest quartile). In terms of reproducibility, 100 (88%) studies reported which bioinformatic tools were used, but only 33 (33%) reported corresponding version numbers. Sequencing data were available for 86 (75%) studies.

INTERPRETATION

The reporting of STROME-ID criteria in genomic epidemiology studies of tuberculosis between 2009 and 2019 was low, with implications for assessment of study quality. The considerable proportion of studies without bioinformatics version numbers or sequencing data available highlights a key concern for reproducibility.

Whole-Genome Sequencing of a Mycobacterium tuberculosis Strain Belonging to Lineage 1 (Indo-Oceanic) and the East African Indian Spoligotype, Isolated in Jazan, Saudi Arabia

We describe here the draft genome sequence of AY1MRC, a Mycobacterium tuberculosis strain belonging to lineage 1 (Indo-Oceanic) and the East African Indian spoligotype, isolated from a patient with tuberculosis in Jazan, Saudi Arabia.

Molecular characterization of multidrug-resistant Mycobacterium tuberculosis (MDR-TB) isolates identifies local transmission of infection in Kuwait, a country with a low incidence of TB and MDR-TB

Background

Increasing incidence of multidrug-resistant Mycobacterium tuberculosis infections is hampering global tuberculosis control efforts. Kuwait is a low-tuberculosis-incidence country, and ~ 1% of M. tuberculosis strains are resistant to rifampicin and isoniazid (MDR-TB). This study detected mutations in seven genes predicting resistance to rifampicin, isoniazid, pyrazinamide, ethambutol and streptomycin in MDR-TB strains. Sequence data were combined with spoligotypes for detecting local transmission of MDR-TB in Kuwait.

Methods

Ninety-three MDR-TB strains isolated from 12 Kuwaiti and 81 expatriate patients and 50 pansusceptible strains were used. Phenotypic drug susceptibility was determined by MGIT 460 TB/960 system. Mutations conferring resistance to rifampicin, isoniazid, pyrazinamide, ethambutol and streptomycin were detected by genotype MTBDRplus assay and/or PCR sequencing of three rpoB regions, katG codon 315 (katG315) + inhA regulatory region, pncA, three embB regions and rpsL + rrs-500–900 regions. Spoligotyping kit was used, spoligotypes were identified by SITVIT2, and phylogenetic tree was constructed by using MIRU-VNTRplus software. Phylogenetic tree was also constructed from concatenated sequences by MEGA7 software. Additional PCR sequencing of gidB and rpsA was performed for cluster isolates.

Results

Pansusceptible isolates contained wild-type sequences. Mutations in rpoB and katG and/or inhA were detected in 93/93 and 92/93 MDR-TB strains, respectively. Mutations were also detected for pyrazinamide resistance, ethambutol resistance and streptomycin resistance in MDR-TB isolates in pncA, embB and rpsL + rrs, respectively. Spoligotyping identified 35 patterns with 18 isolates exhibiting unique patterns while 75 isolates grouped in 17 patterns. Beijing genotype was most common (32/93), and 11 isolates showed nine orphan patterns. Phylogenetic analysis of concatenated sequences showed unique patterns for 51 isolates while 42 isolates grouped in 16 clusters. Interestingly, 22 isolates in eight clusters by both methods were isolated from TB patients typically within a span of 2 years. Five of eight clusters were confirmed by additional gidB and rpsA sequence data.

Conclusions

Our study provides the first insight into molecular epidemiology of MDR-TB in Kuwait and identified several potential clusters of local transmission of MDR-TB involving 2–6 subjects which had escaped detection by routine surveillance studies. Prospective detection of resistance-conferring mutations can identify possible cases of local transmission of MDR-TB in low MDR-TB settings.