Genetic variation of Mycobacterium tuberculosis circulating in Kharkiv Oblast, Ukraine

Phylogenetic relationships of Mycobacterium tuberculosis isolates in Poland: The emergence of Beijing genotype among multidrug-resistant cases

Introduction

The epidemiological situation of tuberculosis (TB) in Poland urges for its continuous and scrupulous monitoring. The objective of this study was to explore the genetic diversity of multidrug-resistant (MDR) and drug-susceptible (DS) Mycobacterium tuberculosis isolates from Poland with a combination of spoligotyping and high-resolution mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) analysis. The results were placed in the Northern and Eastern Europe context.

Methods

The study included 89 (39 MDR and 50 DS) M. tuberculosis isolates collected from as many patients between 2018 and 2021 in Poland. The analysis was done using spoligotyping, and MIRU-VNTR typing at 24 standard loci. The data were compared to those available on Poland and neighbors and global M. tuberculosis datasets.

Results

The main identified families were Beijing (28.1%) and Haarlem (16.8%) while 34.8% of isolates were in the heterogeneous L4-unclassified group. Although the Beijing family was the most prevalent (61.5%) among MDR-TB cases, it accounted for only 2% of DS isolates. Among foreign-born patients, a higher ratio of MDR isolates were observed when compared with those who Poland-born (64.3% vs. 40%). Furthermore, all patients from the Former Soviet Union (FSU) countries were infected with MDR-TB.

Discussion

Whereas DS M. tuberculosis population in Poland is dominated by L4 isolates, MDR isolates are mostly of the Beijing genotype. The rise in the prevalence of the Beijing isolates in Poland, coupled with high proportion of the Beijing genotype among foreign-born TB patients may reflect an ongoing transmission of this family, imported to Poland mainly from FSU countries.

Whole genomic sequencing as a tool for diagnosis of drug and multidrug-resistance tuberculosis in an endemic region in Mexico

Background

Whole genome sequencing (WGS) has been proposed as a tool for diagnosing drug resistance in tuberculosis. However, reports of its effectiveness in endemic countries with important numbers of drug resistance are scarce. The goal of this study was to evaluate the effectiveness of this procedure in isolates from a tuberculosis endemic region in Mexico.

Methods

WGS analysis was performed in 81 tuberculosis positive clinical isolates with a known phenotypic profile of resistance against first-line drugs (isoniazid, rifampin, ethambutol, pyrazinamide and streptomycin). Mutations related to drug resistance were identified for each isolate; drug resistant genotypes were predicted and compared with the phenotypic profile. Genotypes and transmission clusters based on genetic distances were also characterized.

Findings

Prediction by WGS analysis of resistance against isoniazid, rifampicin, ethambutol, pyrazinamide and streptomycin showed sensitivity values of 84%, 96%, 71%, 75% and 29%, while specificity values were 100%, 94%, 90%, 90% and 98%, respectively. Prediction of multidrug resistance showed a sensitivity of 89% and specificity of 97%. Moreover, WGS analysis revealed polymorphisms related to second-line drug resistance, enabling classification of eight and two clinical isolates as pre- and extreme drug-resistant cases, respectively. Lastly, four lineages were identified in the population (L1, L2, L3 and L4). The most frequent of these was L4, which included 90% (77) of the isolates. Six transmission clusters were identified; the most frequent was TC6, which included 13 isolates with a L4.1.1 and a predominantly multidrug-resistant condition.

Conclusions

The results illustrate the utility of WGS for establishing the potential for prediction of resistance against first and second line drugs in isolates of tuberculosis from the region. They also demonstrate the feasibility of this procedure for use as a tool to support the epidemiological surveillance of drug- and multidrug-resistant tuberculosis.

Sequencing of the entire rpob gene and characterization of mutations in isolates of Mycobacterium tuberculosis circulating in an endemic tuberculosis setting

OBJECTIVE

To evaluate the use of a sequencing procedure for the entire rpoB gene of Mycobacterium tuberculosis to identify mutations pre-rifampicin resistance determining region (RRDR), within RRDR, and post-RRDR in isolates circulating in a region affected by tuberculosis (TB).

METHODS

Five primers were designed, with which five DNA fragments of rpoB were obtained, sequenced by Sanger, and analysed in silico in order to identify mutations over the entire rpoB gene in rifampicin-sensitive and rifampicin-resistant TB.

RESULTS

It was possible to analyse the entire rpoB gene in five rifampicin-sensitive and 15 rifampicin-resistant isolates. Thirty-six mutations were identified. Two mutations were found pre-RRDR, nine within-RRDR and 25 post-RRDR. The most frequent mutations within RRDR were S531L (53%), followed by S512T (20%), all of which were found in rifampicin-resistant isolates. Of the 25 mutations found post-RRDR, 14 were only in resistant isolates, and the most frequent was D853N, which was present in 85% of isolates. Mutations E818K, D836N and T882P were observed in 80% of the rifampicin-resistant and rifampicin-sensitive isolates.

CONCLUSIONS

The proposed sequencing method allowed identification of mutations in the entire rpoB gene. This procedure represents a useful tool for diagnosing rifampicin resistance. The number of mutations that were found raises new questions about the diversity of mutations in the rpoB gene and their role in rifampicin resistance in regions where TB is endemic.

Genetic diversity of drug and multidrug-resistant Mycobacterium tuberculosis circulating in Veracruz, Mexico

Background

Mexico is one of the most important contributors of drug and multidrug-resistant tuberculosis in Latin America; however, knowledge of the genetic diversity of drug-resistant tuberculosis isolates is limited.

Methods

In this study, the genetic structure of 112 Mycobacterium tuberculosis strains from the southeastern Mexico was determined by spoligotyping and 24-loci MIRU-VNTRs.

Findings

The results show eight major lineages, the most of which was T1 (24%), followed by LAM (16%) and H (15%). A total of 29 (25%) isolates were identified as orphan. The most abundant SITs were SIT53/T1 and SIT42/LAM9 with 10 isolates each and SIT50/H3 with eight isolates. Fifty-two spoligotype patterns, twenty-seven clusters and ten clonal complexes were observed, demonstrating an important genetic diversity of drug and multidrug-resistant tuberculosis isolates in circulation and transmission level of these aggravated forms of tuberculosis. Being defined as orphan or as part of an orphan cluster, was a risk factor for multidrug resistant-tuberculosis (OR 2.5, IC 1.05–5.86 and OR 3.3, IC 1–11.03, respectively). Multiple correspondence analyses showed association of some clusters and SITs with specific geographical locations.

Conclusions

Our study provides one of the most detailed description of the genetic structure of drug and multidrug-resistant tuberculosis strains in southeast Mexico, establishing for the first time a baseline of the genotypes observed in resistant isolates circulating, however further studies are required to better elucidate the genetic structure of tuberculosis in region and the factors that could be participating in their dispersion.

Genetic Diversity of Mycobacterium tuberculosis Isolates from Assam, India: Dominance of Beijing Family and Discovery of Two New Clades Related to CAS1_Delhi and EAI Family Based on Spoligotyping and MIRU-VNTR Typing

Tuberculosis (TB) is one of the major public health concerns in Assam, a remote state located in the northeastern (NE) region of India. The present study was undertaken to explore the circulating genotypes of Mycobacterium tuberculosis complex (MTBC) in this region. A total of 189 MTBC strains were collected from smear positive pulmonary tuberculosis cases from different designated microscopy centres (DMC) from various localities of Assam. All MTBC isolates were cultured on Lowenstein-Jensen (LJ) media and subsequently genotyped using spoligotyping and 24-loci mycobacterial interspersed repetitive units-variable number of tandem repeats (MIRU-VNTR) typing. Spoligotyping of MTBC isolates revealed 89 distinct spoligo patterns. The most dominant MTBC strain belonged to Beijing lineage and was represented by 35.45% (n = 67) of total isolates, followed by MTBC strains belonging to Central Asian-Delhi (CAS/Delhi) lineage and East African Indian (EAI5) lineage. In addition, in the present study 43 unknown spoligo patterns were detected. The discriminatory power of spoligotyping was found to be 0.8637 based on Hunter Gaston Discriminatory Index (HGDI). On the other hand, 24-loci MIRU-VNTR typing revealed that out of total 189 MTBC isolates from Assam 185 (97.9%) isolates had unique MIRU-VNTR profiles and 4 isolates grouped into 2 clusters. Phylogenetic analysis of 67 Beijing isolates based on 24-loci MIRU-VNTR typing revealed that Beijing isolates from Assam represent two major groups, each comprising of several subgroups. Neighbour-Joining (NJ) phylogenetic tree analysis based on combined spoligotyping and 24-loci MIRU-VNTR data of 78 Non-Beijing isolates was carried out for strain lineage identification as implemented by MIRU-VNTRplus database. The important lineages of MTBC identified were CAS/CAS1_Delhi (41.02%, n = 78) and East-African-Indian (EAI, 33.33%). Interestingly, phylogenetic analysis of orphan (23.28%) MTBC spoligotypes revealed that majority of these orphan isolates from Assam represent two new sub-clades Assam/EAI and Assam/CAS. The prevalence of multidrug resistance (MDR) in Beijing and Non-Beijing strains was found to be 10.44% and 9.01% respectively. In conclusion, the present study has shown the predominance of Beijing isolates in Assam which is a matter of great concern because Beijing strains are considered to be ecologically more fit enabling wider dissemination of M. tuberculosis. Other interesting finding of the present study is the discovery of two new clades of MTBC isolates circulating in Assam. More elaborate longitudinal studies are required to be undertaken in this region to understand the transmission dynamics of MTBC.

The Association between Mycobacterium Tuberculosis Genotype and Drug Resistance in Peru

Background

The comparison of Mycobacterium tuberculosis bacterial genotypes with phenotypic, demographic, geospatial and clinical data improves our understanding of how strain lineage influences the development of drug-resistance and the spread of tuberculosis.

Methods

To investigate the association of Mycobacterium tuberculosis bacterial genotype with drug-resistance. Drug susceptibility testing together with genotyping using both 15-loci MIRU-typing and spoligotyping, was performed on 2,139 culture positive isolates, each from a different patient in Lima, Peru. Demographic, geospatial and socio-economic data were collected using questionnaires, global positioning equipment and the latest national census.

Results

The Latin American Mediterranean (LAM) clade (OR 2.4, p<0.001) was significantly associated with drug-resistance and alone accounted for more than half of all drug resistance in the region. Previously treated patients, prisoners and genetically clustered cases were also significantly associated with drug-resistance (OR's 2.5, 2.4 and 1.8, p<0.001, p<0.05, p<0.001 respectively).

Conclusions

Tuberculosis disease caused by the LAM clade was more likely to be drug resistant independent of important clinical, genetic and socio-economic confounding factors. Explanations for this include; the preferential co-evolution of LAM strains in a Latin American population, a LAM strain bacterial genetic background that favors drug-resistance or the "founder effect" from pre-existing LAM strains disproportionately exposed to drugs.

Map the gap: missing children with drug-resistant tuberculosis

BACKGROUND

The lack of published information about children with multidrug-resistant tuberculosis (MDR-TB) is an obstacle to efforts to advocate for better diagnostics and treatment.

OBJECTIVE

To describe the lack of recognition in the published literature of MDR-TB and extensively drug-resistant TB (XDR-TB) in children.

DESIGN

We conducted a systematic search of the literature published in countries that reported any MDR- or XDR-TB case by 2012 to identify MDR- or XDR-TB cases in adults and in children.

RESULTS

Of 184 countries and territories that reported any case of MDR-TB during 2005-2012, we identified adult MDR-TB cases in the published literature in 143 (78%) countries and pediatric MDR-TB cases in 78 (42%) countries. Of the 92 countries that reported any case of XDR-TB, we identified adult XDR-TB cases in the published literature in 55 (60%) countries and pediatric XDR-TB cases for 9 (10%) countries.

CONCLUSION

The absence of publications documenting child MDR- and XDR-TB cases in settings where MDR- and XDR-TB in adults have been reported indicates both exclusion of childhood disease from the public discourse on drug-resistant TB and likely underdetection of sick children. Our results highlight a large-scale lack of awareness about children with MDR- and XDR-TB.

Mycobacterium tuberculosis phylogeography in the context of human migration and pathogen's pathobiology: Insights from Beijing and Ural families

Here, I review the population structure and phylogeography of the two contrasting families of Mycobacterium tuberculosis, Beijing and Ural, in the context of strain pathobiology and human history and migration. Proprietary database (12-loci MIRU-VNTR profiles of 3067 Beijing genotype isolates) was subjected to phylogenetic and statistical analysis. The highest rate (90%) and diversity (HGI 0.80-0.95) of the Beijing genotype in North China suggest it to be its area of origin. Under VNTR-based MDS analysis the interpopulation genetic distances correlated with geography over uninterrupted landmasses. In contrast, large water distances together with long time generated remarkable outliers. Weak and less expected affinities of the distant M. tuberculosis populations may reflect hidden epidemiological links due to unknown migration. Association with drug-resistance or increased virulence/transmissibility along with particular human migration flows shape global dissemination of some Beijing clones. The paucity of data on the Ural genotype prevents from high-resolution analysis that was mainly based on the available spoligotyping data. The North/East Pontic area marked with the highest prevalence of the Ural family may have been the area of its origin and primary dispersal in Eurasia. Ural strains are not marked by increased pathogenic capacities, increased transmissibility and association with drug resistance (but most recent reports describe an alarming increase of MDR Ural strains in some parts of eastern Europe and northwestern Russia). Large-scale SNP or WGS population-based studies targeting strains from indigenous populations and, eventually, analysis of ancient DNA will better test these hypotheses. Host genetics factors likely play the most prominent role in differential dissemination of particular M. tuberculosis genotypes.

Effect of study design and setting on tuberculosis clustering estimates using Mycobacterial Interspersed Repetitive Units-Variable Number Tandem Repeats (MIRU-VNTR): a systematic review

OBJECTIVES

To systematically review the evidence for the impact of study design and setting on the interpretation of tuberculosis (TB) transmission using clustering derived from Mycobacterial Interspersed Repetitive Units-Variable Number Tandem Repeats (MIRU-VNTR) strain typing.

DATA SOURCES

MEDLINE, EMBASE, CINHAL, Web of Science and Scopus were searched for articles published before 21st October 2014.

REVIEW METHODS

Studies in humans that reported the proportion of clustering of TB isolates by MIRU-VNTR were included in the analysis. Univariable meta-regression analyses were conducted to assess the influence of study design and setting on the proportion of clustering.

RESULTS

The search identified 27 eligible articles reporting clustering between 0% and 63%. The number of MIRU-VNTR loci typed, requiring consent to type patient isolates (as a proxy for sampling fraction), the TB incidence and the maximum cluster size explained 14%, 14%, 27% and 48% of between-study variation, respectively, and had a significant association with the proportion of clustering.

CONCLUSIONS

Although MIRU-VNTR typing is being adopted worldwide there is a paucity of data on how study design and setting may influence estimates of clustering. We have highlighted study design variables for consideration in the design and interpretation of future studies.

Molecular characterisation of extensively drug-resistant Mycobacterium tuberculosis isolates in China

The emergence of extensively drug-resistant tuberculosis (XDR-TB) in China is a great threat to TB control. To determine the molecular characterisation of XDR-TB isolates from China and the correlations between specific drug resistance-associated mutations and different genotype strains, 58 XDR-TB isolates were sequenced in eight drug loci, including katG, inhA, oxyR-ahpC intergenic region, rpoB, eis, rrs, gyrA and gyrB, and were genotyped using spoligotyping and analysis of the noise transfer function region. Compared with the phenotypic data, the sensitivities and specificities for DNA sequencing were 87.9% and 100.0% for isoniazid (INH), 91.4% and 98.3% for rifampicin (RIF), 60.4% and 100.0% for kanamycin (KAN) and 81.0% and 100.0% for ofloxacin (OFX), respectively. A combination of eight drug loci predicted XDR-TB phenotypes with 53.4% sensitivity (31/58 isolates) and 100.0% specificity. The most frequent mutations among these XDR-TB isolates were katG315 and inhA-15 (for INH), 531, 526 and 516 in rpoB (for RIF), rrs1401 and eis-10 (for KAN) and 94, 90 and 91 in gyrA (for OFX). Also, among these XDR-TB isolates, 44 (75.9%) were identified as Beijing genotype strain, of which 31 (70.5%) belonged to the modern Beijing sublineage. inhA-8, rpoB526 and rpoB531 mutations demonstrated significant statistical associations with ancient and modern Beijing family sublineage (P<0.05). However, Beijing and non-Beijing genotypes showed no association with specific resistance-conferring mutations. These results will be helpful in designing new molecular biology-based techniques to diagnose XDR-TB in China.

Molecular characterization of multidrug-resistant Mycobacterium tuberculosis isolates from China

To investigate the molecular characterization of multidrug-resistant tuberculosis (MDR-TB) isolates from China and the association of specific mutations conferring drug resistance with strains of different genotypes, we performed spoligotyping and sequenced nine loci (katG, inhA, the oxyR-ahpC intergenic region, rpoB, tlyA, eis, rrs, gyrA, and gyrB) for 128 MDR-TB isolates. Our results showed that 108 isolates (84.4%) were Beijing family strains, 64 (59.3%) of which were identified as modern Beijing strains. Compared with the phenotypic data, the sensitivity and specificity of DNA sequencing were 89.1% and 100.0%, respectively, for isoniazid (INH) resistance, 93.8% and 100.0% for rifampin (RIF) resistance, 60.0% and 99.4% for capreomycin (CAP) resistance, 84.6% and 99.4% for kanamycin (KAN) resistance, and 90.0% and 100.0% for ofloxacin (OFX) resistance. The most prevalent mutations among the MDR-TB isolates were katG315, inhA15, rpoB531, -526, and -516, rrs1401, eis-10, and gyrA94, -90, and -91. Furthermore, there was no association between specific resistance-conferring mutations and the strain genotype. These findings will be helpful for the establishment of rapid molecular diagnostic methods to be implemented in China.

Insights into the origin, emergence, and current spread of a successful Russian clone of Mycobacterium tuberculosis

Mycobacterium tuberculosis variant Beijing B0/W148 is regarded as a successful clone of M. tuberculosis that is widespread in the former Soviet Union and respective immigrant communities. Understanding the pathobiology and phylogeography of this notorious strain may help to clarify its origin and evolutionary history and the driving forces behind its emergence and current dissemination. I present the first review and analysis of all available data on the subject. In spite of the common perception of the omnipresence of B0/W148 across post-Soviet countries, its geographic distribution shows a peculiar clinal gradient. Its frequency peaks in Siberian Russia and, to a lesser extent, in the European part of the former Soviet Union. In contrast, the frequency of B0/W148 is sharply decreased in the Asian part of the former Soviet Union, and it is absent in autochthonous populations elsewhere in the world. Placing the molecular, clinical, and epidemiological features in a broad historical, demographic, and ecological context, I put forward two interdependent hypotheses. First, B0/W148 likely originated in Siberia, and its primary dispersal was driven by a massive population outflow from Siberia to European Russia in the 1960s to 1980s. Second, a historically recent, phylogenetically demonstrated successful dissemination of the Beijing B0/W148 strain was triggered by the advent and wide use of modern antituberculosis (anti-TB) drugs and was due to the remarkable capacity of this strain to acquire drug resistance. In contrast, there is some indication, but not yet systematic proof, of an enhanced virulence of this strain.

Multidrug-resistant Myobacterium tuberculosis caused by the Beijing genotype and a specific T1 genotype clone (SIT No. 266) is widely transmitted in Minsk

SETTING

This study was performed in the city of Minsk in Belarus, where a very severe problem with MDR-TB was demonstrated in a recent drug resistant survey.

OBJECTIVE

The aim of this study was to use molecular typing of MDR and pan-susceptible clinical isolates of Mycobacterium tuberculosis to increase the understanding of the transmission patterns and possible differences between the strains causing susceptible and drug-resistant tuberculosis.

STUDY POPULATION AND METHODS

Consecutive isolates from pulmonary TB patients in Minsk were collected at the Belarusian National Reference Laboratory. Isolates found to be either pan-susceptible or MDR were included in the study, which totally comprised 81 MDR and 82 pan-susceptible clinical isolates. All isolates were characterized by spoligotyping. The major clusters were characterized using sequencing of the pncA gene.

RESULTS

Three out of four MDR cases were caused by one out of two drug-resistant clones of M. tuberculosis belonging to the Beijing and T1 genotypes, respectively. A single T1 clone, SIT No. 266, found exclusively in the MDR cohort, was shown to cause no less than 30% of all MDR-TB cases.

DISCUSSION

The findings indicate that the major cause of MDR-TB in Minsk is an ongoing transmission of certain already resistant M. tuberculosis strains.

CONCLUSION

The significant transmission of MDR-TB in Minsk underlines the urgent need for strengthened infection control measures to limit the transmission in order to better control MDR-TB.

Molecular snapshot of Mycobacterium tuberculosis population structure and drug-resistance in Kyrgyzstan

Kyrgyzstan is a post-Soviet country in Central Asia marked with high incidence and mortality rates of tuberculosis (TB). The present study provided first assessment of Mycobacterium tuberculosis population structure and drug-resistance in civilian population here. The collection included 103 M. tuberculosis DNA samples subjected to the analysis of rifampin and isoniazid resistance mutations and spoligotyping. The major spoligotype-defined families were Beijing (n = 62), T (n = 14), LAM (n = 9), Ural-2 (n = 6) and Ural-1 (n = 3). Genotypically, 20 isolates were RIF-resistant, 28 were INH-resistant, 17 were multidrug-resistant. Drug resistant isolates were more prevalent among Beijing than non-Beijing groups (P = 0.03). The predominance of the mainly "Russian" spoligotypes among the non-Beijing strains (LAM-RUS and Ural-1) in this study along with previously demonstrated prevalence of the Russia-specific subtype of the Beijing family in Kyrgyz prison (Mokrousov et al., 2009) suggest that the current population structure of M. tuberculosis in Kyrgyzstan has been mainly formed within the course of the 20th century when the country was a part of the Russian Empire and Soviet Union. On the other hand, a prevalence of the Asia-specific Ural-2 type in the oldest age group (68-85 years old; P < 0.0001) may present a heritage of the more distant historical events. In summary, we suggest: (i) a clear shift of the local M. tuberculosis population structure during the last 100 years and (ii) a critical impact of the Beijing genotype on the current situation with drug resistant TB in Kyrgyzstan.

Mycobacterium tuberculosis population in northwestern Russia: an update from Russian-EU/Latvian border region

This study aimed to characterize the population structure of Mycobacterium tuberculosis in Pskov oblast in northwestern Russia, to view it in the geographical context, to compare drug resistance properties across major genetic families. Ninety M. tuberculosis strains from tuberculosis (TB) patients, permanent residents in Pskov oblast were subjected to LAM-specific IS6110-PCR and spoligotyping, followed by comparison with SITVITWEB and MIRU-VNTRplus databases. The Beijing genotype (n = 40) was found the most prevalent followed by LAM (n = 18), T (n = 13), Haarlem (n = 10), Ural (n = 5), and Manu2 (n = 1); the family status remained unknown for 3 isolates. The high rate of Beijing genotype and prevalence of LAM family are similar to those in the other Russian settings. A feature specific for M. tuberculosis population in Pskov is a relatively higher rate of Haarlem and T types. Beijing strains were further typed with 12-MIRU (followed by comparison with proprietary global database) and 3 hypervariable loci QUB-3232, VNTR-3820, VNTR-4120. The 12-MIRU typing differentiated 40 Beijing strains into 14 types (HGI = 0.82) while two largest types were M2 (223325153533) prevalent throughout former USSR and M11 (223325173533) prevalent in Russia and East Asia. The use of 3 hypervariable loci increased a discrimination of the Beijing strains (18 profiles, HGI = 0.89). Both major families Beijing and LAM had similar rate of MDR strains (62.5 and 55.6%, respectively) that was significantly higher than in other strains (21.9%; P = 0.001 and 0.03, respectively). The rpoB531 mutations were more frequently found in Beijing strains while LAM drug resistant strains mainly harbored rpoB516 and inhA -15 mutations. Taken together with a high rate of multidrug resistance among Beijing strains from new TB cases (79.3% versus 44.4% in LAM), these findings suggest the critical impact of the Beijing genotype on the current situation with MDR-TB in the Pskov region in northwestern Russia.

New insight into the molecular characterization of isoniazid and rifampicin resistant Mycobacterium tuberculosis strains from Saudi Arabia

Data on the genetic variation of isolates of Mycobacterium tuberculosis and spectrum of mutations determining resistance to principal anti-tuberculosis drugs isoniazid (INH) and rifampicin (RIF) have not yet been studied in Saudi Arabia. One hundred and fifty-one clinical isolates of M. tuberculosis from different regions in the country showing resistance to RIF and INH were subjected to drug susceptibility testing, characterization of mutations conferring drug resistance and genotyping. Phenotypically 17 (11.3%) isolates were resistance to RIF, 75 (49.6%) were resistant to INH and 59 (39.1%) were resistant to both RIF and INH, respectively. Sixteen (10.6%), 74 (49%) and 56 (37.1%) were determined as resistant to RIF, INH and to both by line probe assay. High frequency of rpoB 531 mutations (67.1%) in RIF resistant strains and katG 315 mutations (65.2%) in INH resistant strains were found. Mutations responsible for INH resistance, katG 315 (P value<0.001, odds ratio: 1.81, 95% CI [1.51, 2.18]) and inhA-15 (P value - 0.004, odds ratio: 1.48, 95% CI [1.22, 1.8]) were predominant among the newly diagnosed cases. Beijing strains were significantly associated with multi drug resistance and mutations in combination of rpoB531 and katG315 (P value - <0.001, odds ratio: 6.83, 95% CI [2.65, 17.58]). In addition multi drug resistance was significantly associated with treatment history (P value<0.001, odds ratio: 3.16, 95% CI [2.14, 4.67]). Furthermore, a higher rate (39.3%) of clustering among the multidrug resistant strains particularly with Beijing family (52.9%) was observed. Saudi Arabia harbors highly diverse drug resistant M. tuberculosis population with an ongoing transmission which needs to be immediately managed.