Mycobacterium tuberculosis lineages and anti-tuberculosis drug resistance in reference hospitals across Viet Nam

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.

Distribution and origins of Mycobacterium tuberculosis L4 in Southeast Asia

Molecular and genomic studies have revealed that Mycobacterium tuberculosis Lineage 4 (L4, Euro-American lineage) emerged in Europe before becoming distributed around the globe by trade routes, colonial migration and other historical connections. Although L4 accounts for tens or hundreds of thousands of tuberculosis (TB) cases in multiple Southeast Asian countries, phylogeographical studies have either focused on a single country or just included Southeast Asia as part of a global analysis. Therefore, we interrogated public genomic data to investigate the historical patterns underlying the distribution of L4 in Southeast Asia and surrounding countries. We downloaded 6037 genomes associated with 29 published studies, focusing on global analyses of L4 and Asian studies of M. tuberculosis. We identified 2256 L4 genomes including 968 from Asia. We show that 81 % of L4 in Thailand, 51 % of L4 in Vietnam and 9 % of L4 in Indonesia belong to sub-lineages of L4 that are rarely seen outside East and Southeast Asia (L4.2.2, L4.4.2 and L4.5). These sub-lineages have spread between East and Southeast Asian countries, with no recent European ancestor. Although there is considerable uncertainty about the exact direction and order of intra-Asian M. tuberculosis dispersal, due to differing sampling frames between countries, our analysis suggests that China may be the intermediate location between Europe and Southeast Asia for two of the three predominantly East and Southeast Asian L4 sub-lineages (L4.2.2 and L4.5). This new perspective on L4 in Southeast Asia raises the possibility of investigating host population-specific evolution and highlights the need for more structured sampling from Southeast Asian countries to provide more certainty of the historical and current routes of dispersal.

Tuberculosis Treatment Monitoring and Outcome Measures: New Interest and New Strategies

Despite the advent of new diagnostics, drugs and regimens, tuberculosis (TB) remains a global public health threat. A significant challenge for TB control efforts has been the monitoring of TB therapy and determination of TB treatment success. Current recommendations for TB treatment monitoring rely on sputum and culture conversion, which have low sensitivity and long turnaround times, present biohazard risk, and are prone to contamination, undermining their usefulness as clinical treatment monitoring tools and for drug development. We review the pipeline of molecular technologies and assays that serve as suitable substitutes for current culture-based readouts for treatment response and outcome with the potential to change TB therapy monitoring and accelerate drug development.

Mycobacterium tuberculosis typing using Allele-specific oligonucleotide multiplex PCR (ASO-PCR) method

Mycobacterium tuberculosis (M. tuberculosis) genotyping provides valuable information related to the origin and the evolution of the isolates. This study aimed to evaluate the applicability of single-nucleotide polymorphisms (SNPs) technique for lineages identification of M. tuberculosis and compare it with mycobacterial interspersed repetitive units-variable number of tandem repeats (MIRU-VNTR) method. The lineages of 162 clinically isolates were evaluated using six pair primers by Multiplex-PCR based on SNPs. Among 162 isolates, 70 (43.2%) isolates were lineage 4, following that 62 (38.3%) and 22 (13.6%) isolates were lineage 3 and 2, respectively. The method could not type 8 (4.9%) isolates. Moreover, we could identify 71 out of 79 unknown isolates resulted from the MIRU-VNTR method. The results showed that the SNP typing method has the potential to determine the lineages of M. tuberculosis as a rapid laboratory screening test.

Genotypic and phenotypic diversity of Mycobacterium tuberculosis complex genotypes prevalent in West Africa

Findings from previous comparative genomics studies of the Mycobacterium tuberculosis complex (MTBC) suggest genomic variation among the genotypes may have phenotypic implications. We investigated the diversity in the phenotypic profiles of the main prevalent MTBC genotypes in West Africa. Thirty-six whole genome sequenced drug susceptible MTBC isolates belonging to lineages 4, 5 and 6 were included in this study. The isolates were phenotypically characterized for urease activity, tween hydrolysis, Thiophen-2-Carboxylic Acid Hydrazide (TCH) susceptibility, nitric oxide production, and growth rate in both liquid (7H9) and solid media (7H11 and Löwenstein-Jensen (L-J)). Lineage 4 isolates showed the highest growth rate in both liquid (p = 0.0003) and on solid (L-J) media supplemented with glycerol (p<0.001) or pyruvate (p = 0.005). L6 isolates optimally utilized pyruvate compared to glycerol (p<0.001), whereas L5 isolates grew similarly on both media (p = 0.05). Lineage 4 isolates showed the lowest average time to positivity (TTP) (p = 0.01; Average TTP: L4 = 15days, L5 = 16.7days, L6 = 29.7days) and the highest logCFU/mL (p = 0.04; average logCFU/mL L4 = 5.9, L5 = 5.0, L6 = 4.4) on 7H11 supplemented with glycerol, but there was no significant difference in growth on 7H11 supplemented with pyruvate (p = 0.23). The highest release of nitrite was recorded for L5 isolates, followed by L4 and L6 isolates. However, the reverse was observed in the urease activity for the lineages. All isolates tested were resistant to TCH except for one L6 isolate. Comparative genomic analyses revealed several mutations that might explain the diverse phenotypic profiles of these isolates. Our findings showed significant phenotypic diversity among the MTBC lineages used for this study.

Phenotypic and genotypic features of the Mycobacterium tuberculosis lineage 1 subgroup in central Vietnam

Mycobacterium tuberculosis (Mtb) has different features depending on different geographic areas. We collected Mtb strains from patients with smear-positive pulmonary tuberculosis in Da Nang, central Vietnam. Using a whole genome sequencing platform, including genome assembly complemented by long-read-sequencing data, genomic characteristics were studied. Of 181 Mtb isolates, predominant Vietnamese EAI4_VNM and EAI4-like spoligotypes (31.5%), ZERO strains (5.0%), and part of EAI5 (11.1%) were included in a lineage-1 (L1) sublineage, i.e., L1.1.1.1. These strains were found less often in younger people, and they genetically clustered less frequently than other modern strains. Patients infected with ZERO strains demonstrated less lung infiltration. A region in RD2bcg spanning six loci, i.e., PE_PGRS35, cfp21, Rv1985c, Rv1986, Rv1987, and erm(37), was deleted in EAI4_VNM, EAI4-like, and ZERO strains, whereas another 118 bp deletion in furA was specific only to ZERO strains. L1.1.1.1-sublineage-specific deletions in PE_PGRS4 and PE_PGRS22 were also identified. RD900, seen in ancestral lineages, was present in majority of the L1 members. All strains without IS6110 (5.0%) had the ZERO spoligo-pattern. Distinctive features of the ancestral L1 strains provide a basis for investigation of the modern versus ancestral Mtb lineages and allow consideration of countermeasures against this heterogeneous pathogen.

Lineage-Specific Proteomic Signatures in the Mycobacterium tuberculosis Complex Reveal Differential Abundance of Proteins Involved in Virulence, DNA Repair, CRISPR-Cas, Bioenergetics and Lipid Metabolism

Despite the discovery of the tubercle bacillus more than 130 years ago, its physiology and the mechanisms of virulence are still not fully understood. A comprehensive analysis of the proteomes of members of the human-adapted Mycobacterium tuberculosis complex (MTBC) lineages 3, 4, 5, and 7 was conducted to better understand the evolution of virulence and other physiological characteristics. Unique and shared proteomic signatures in these modern, pre-modern and ancient MTBC lineages, as deduced from quantitative bioinformatics analyses of high-resolution mass spectrometry data, were delineated. The main proteomic findings were verified by using immunoblotting. In addition, analysis of multiple genome alignment of members of the same lineages was performed. Label-free peptide quantification of whole cells from MTBC lineages 3, 4, 5, and 7 yielded a total of 38,346 unique peptides derived from 3092 proteins, representing 77% coverage of the predicted proteome. MTBC lineage-specific differential expression was observed for 539 proteins. Lineage 7 exhibited a markedly reduced abundance of proteins involved in DNA repair, type VII ESX-3 and ESX-1 secretion systems, lipid metabolism and inorganic phosphate uptake, and an increased abundance of proteins involved in alternative pathways of the TCA cycle and the CRISPR-Cas system as compared to the other lineages. Lineages 3 and 4 exhibited a higher abundance of proteins involved in virulence, DNA repair, drug resistance and other metabolic pathways. The high throughput analysis of the MTBC proteome by super-resolution mass spectrometry provided an insight into the differential expression of proteins between MTBC lineages 3, 4, 5, and 7 that may explain the slow growth and reduced virulence, metabolic flexibility, and the ability to survive under adverse growth conditions of lineage 7.

The recombination-cold region as an epidemiological marker of recombinogenic opportunistic pathogen Mycobacterium avium

Background

The rapid identification of lineage remains a challenge in the genotyping of clinical isolates of recombinogenic pathogens. The chromosome of Mycobacterium avium subsp. hominissuis (MAH), an agent of Mycobacterium avium complex (MAC) lung disease, is often mosaic and is composed of chromosomal segments originating from different lineages. This makes it difficult to infer the MAH lineage in a simple experimental set-up. To overcome this difficulty, we sought to identify chromosomal marker genes containing lineage-specific alleles by genome data mining.

Results

We conducted genetic population structure analysis, phylogenetic analysis, and a survey of historical recombination using data from 125 global MAH isolates. Six MAH lineages (EA1, EA2, SC1, SC2, SC3, and SC4) were identified in the current dataset. One P-450 gene (locus_tag MAH_0788/MAV_0940) in the recombination-cold region was found to have multiple alleles that could discriminate five lineages. By combining the information about allele type from one additional gene, the six MAH lineages as well as other M. avium subspecies were distinguishable. A recombination-cold region of 116 kb contains an insertion hotspot and is flanked by a mammalian cell-entry protein operon where allelic variants have previously been reported to occur. Hence, we speculate that the acquisition of lineage- or strain-specific insertions has introduced homology breaks in the chromosome, thereby reducing the chance of interlineage recombination.

Conclusions

The allele types of the newly identified marker genes can be used to predict major lineages of M. avium. The single nucleotide polymorphism typing approach targeting multiallelic loci in recombination-cold regions will facilitate the epidemiological study of MAC, and may also be useful for equivalent studies of other nontuberculous mycobacteria potentially carrying mosaic genomes.

First insights into the genetic characteristics and drug resistance of Mycobacterium tuberculosis population collected during the first national tuberculosis prevalence survey of Lao PDR (2010-2011)

Background

In Lao People’s Democratic Republic (PDR), tuberculosis (TB) prevalence was estimated at 540/100,000 in 2011. Nevertheless, little is known about the genetic characteristics and anti-TB drug resistance of the Mycobacterium tuberculosis population. The main objective of this work was to study the genetic characteristics and drug resistance of M. tuberculosis population collected during the first National TB Prevalence Survey (TBPS) of Lao PDR (2010–2011).

Methods

Two hundred and twenty two isolates collected during TBPS (2010–2011) were analyzed with the GenoType MTBDRplus test for M. tuberculosis identification and drug resistance detection. Then, 206 of the 222 isolates were characterized by spoligotyping and MIRU-VNTR typing.

Results

Among the 222 M. tuberculosis isolates, 11 were mono-resistant to isoniazid and 2 were resistant to isoniazid and rifampicin (MDR-TB), using the GenoType MTBDRplus test. Among the 202 genetically characterized isolates, the East African-Indian (EAI) family was predominant (76.7%) followed by the Beijing (14.4%) and T (5.5%) families. EAI isolates came from all the country provinces, whereas Beijing isolates were found mainly in the northern and central provinces. A higher proportion of Beijing isolates was observed in people younger than 35 years compared to EAI. Moreover, the percentage of drug resistance was higher among Beijing (17.2%) than EAI (5.2%) isolates, and the two MDR-TB isolates belonged to the Beijing family. Combined analysis of the MIRU-VNTR and spoligotyping results (n = 202 isolates) revealed an estimated clustering rate of 11% and the occurrence of mini-outbreaks of drug-resistant TB caused by Beijing genotypes.

Conclusions

The EAI family, the ancient and endemic family in Asia, is predominant in Lao PDR whereas the prevalence of Beijing, the most harmful M. tuberculosis family for humans, is still low, differently from neighboring countries. However, its association with drug resistance, its presence in young patients and its potential association with recent transmission suggest that the Beijing family could change TB epidemiological pattern in Lao PDR. Therefore, efficient TB control and surveillance systems must be maintained and reinforced to prevent the emergence of highly transmissible and drug-resistant strains in Lao PDR, as observed in neighboring countries.

Electronic supplementary material

The online version of this article (10.1186/s12879-019-4435-z) contains supplementary material, which is available to authorized users.

Association of gyrA and rrs gene mutations detected by MTBDRsl V1 on Mycobacterium tuberculosis strains of diverse genetic background from India

There is limited data on the use of Genotype MTBDRslVersion 1 (MTBDRsl V1) as an initial rapid screening test to rule out XDR-TB and most importantly its performance in various genotypes of Mycobacterium tuberculosis is scarcely studied. A total of 359 MDR-TB isolates were tested for gene mutations representing second line drug resistance, using the MTBDRsl_V.1 and the results were compared with phenotypic method (Bactec MGIT-960 system) for second-line drug (SLD) susceptibility testing. Genetic lineages of all these isolates were also determined using spoligotyping and SITVIT2 WEB database. The MTBDRsl V1 detected mutations in the gyrA, rrs, and emb genes in 108 (30%), 2 (0.5%) and 129 (35.9%) isolates, respectively. Remaining 120 (33.4%) had no second line drug (SLD) resistance. In 17 (4.7%) isolates mutations were detected in both gyrA and rrs genes. Its concordance with MGIT-960 culture drug susceptibility testing (DST) was 97% and 94.1%, 93.5%, 60.5% and 50% for the detection of XDR-TB, pre-XDR, Ethambutol, and Aminoglycosides/Cyclopeptides resistance. The Beijing lineage was predominant (46%) between both the pre-XDR/XDR-TB isolates. We conclude that MTBDRsl is useful for rapid detection of SLD resistance. Also in pre-XDR and XDR-TB isolates the frequency of relevant genetic mutations was significantly higher in the Beijing strains.

Frequent transmission of the Mycobacterium tuberculosis Beijing lineage and positive selection for the EsxW Beijing variant in Vietnam

To examine the transmission dynamics of Mycobacterium tuberculosis (Mtb) isolated from tuberculosis patients in Ho Chi Minh City, Vietnam, we sequenced the whole genomes of 1,635 isolates and compared these with 3,144 isolates from elsewhere. The data identify an underlying burden of disease caused by the endemic Mtb lineage 1 associated with the activation of long-term latent infection, and a threefold higher burden associated with the more recently introduced Beijing lineage and lineage 4 Mtb strains. We find that Beijing lineage Mtb is frequently transferred between Vietnam and other countries, and detect higher levels of transmission of Beijing lineage strains within this host population than the endemic lineage 1 Mtb. Screening for parallel evolution of Beijing lineage-associated SNPs in other Mtb lineages as a signal of positive selection, we identify an alteration in the ESX-5 type VII-secreted protein EsxW, which could potentially contribute to the enhanced transmission of Beijing lineage Mtb in Vietnamese and other host populations.

Frequent transmission of the Mycobacterium tuberculosis Beijing lineage and positive selection for the EsxW Beijing variant in Vietnam

To examine the transmission dynamics of Mycobacterium tuberculosis (Mtb) isolated from tuberculosis patients in Ho Chi Minh City, Vietnam, we sequenced the whole genomes of 1,635 isolates and compared these with 3,144 isolates from elsewhere. The data identify an underlying burden of disease caused by the endemic Mtb lineage 1 associated with the activation of long-term latent infection, and a threefold higher burden associated with the more recently introduced Beijing lineage and lineage 4 Mtb strains. We find that Beijing lineage Mtb is frequently transferred between Vietnam and other countries, and detect higher levels of transmission of Beijing lineage strains within this host population than the endemic lineage 1 Mtb. Screening for parallel evolution of Beijing lineage-associated SNPs in other Mtb lineages as a signal of positive selection, we identify an alteration in the ESX-5 type VII-secreted protein EsxW, which could potentially contribute to the enhanced transmission of Beijing lineage Mtb in Vietnamese and other host populations.

Importance of differential identification of Mycobacterium tuberculosis strains for understanding differences in their prevalence, treatment efficacy, and vaccine development

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb), remains a serious global health problem in the 21^st century because of its high mortality. Mtb is an extremely successful human-adapted pathogen that displays a multifactorial ability to control the host immune response and to evade killing by drugs, resulting in the breakdown of BCG vaccine-conferred anti-TB immunity and development of multidrug-resistant (MDR) and extensively drug-resistant (XDR) Mtb. Although genetic components of the genomes of the Mtb complex strains are highly conserved, showing over 99% similarity to other bacterial genera, recently accumulated evidence suggests that the genetic diversity of the Mtb complex strains has implications for treatment outcomes, development of MDR/XDR Mtb, BCG vaccine efficacy, transmissibility, and epidemiological outbreaks. Thus, new insights into the pathophysiological features of the Mtb complex strains are required for development of novel vaccines and for control of MDR/XDR Mtb infection, eventually leading to refinement of treatment regimens and the health care system. Many studies have focused on the differential identification of Mtb complex strains belonging to different lineages because of differences in their virulence and geographical dominance. In this review, we discuss the impact of differing genetic characteristics among Mtb complex strains on vaccine efficacy, treatment outcome, development of MDR/XDR Mtb strains, and epidemiological outbreaks by focusing on the best-adapted human Mtb lineages. We further explore the rationale for differential identification of Mtb strains for more effective control of TB in clinical and laboratory settings by scrutinizing current diagnostic methods.

Bacterial risk factors for treatment failure and relapse among patients with isoniazid resistant tuberculosis

Background

Drug resistant tuberculosis (TB) is increasing in prevalence worldwide. Treatment failure and relapse is known to be high for patients with isoniazid resistant TB treated with standard first line regimens. However, risk factors for unfavourable outcomes and the optimal treatment regimen for isoniazid resistant TB are unknown. This cohort study was conducted when Vietnam used the eight month first line treatment regimen and examined risk factors for failure/relapse among patients with isoniazid resistant TB.

Methods

Between December 2008 and June 2011 2090 consecutive HIV-negative adults (≥18 years of age) with new smear positive pulmonary TB presenting at participating district TB units in Ho Chi Minh City were recruited. Participants with isoniazid resistant TB identified by Microscopic Observation Drug Susceptibility (MODS) had extended follow-up for 2 years with mycobacterial culture to test for relapse. MGIT drug susceptibility testing confirmed 239 participants with isoniazid resistant, rifampicin susceptible TB. Bacterial and demographic factors were analysed for association with treatment failure and relapse.

Results

Using only routine programmatic sputum smear microscopy for assessment, (months 2, 5 and 8) 30/239 (12.6%) had an unfavourable outcome by WHO criteria. Thirty-nine patients were additionally detected with unfavourable outcomes during 2 year follow up, giving a total of 69/239 (28.9%) of isoniazid (INH) resistant cases with unfavourable outcome by 2 years of follow-up. Beijing lineage was the only factor significantly associated with unfavourable outcome among INH-resistant TB cases during 2 years of follow-up. (adjusted OR = 3.16 [1.54–6.47], P = 0.002).

Conclusion

One third of isoniazid resistant TB cases suffered failure/relapse within 2 years under the old eight month regimen. Over half of these cases were not identified by standard WHO recommended treatment monitoring. Intensified research on early identification and optimal regimens for isoniazid resistant TB is needed. Infection with Beijing genotype of TB is a significant risk factor for bacterial persistence on treatment resulting in failure/relapse within 2 years. The underlying mechanism of increased tolerance for standard drug regimens in Beijing genotype strains remains unknown.

Spoligotyping and whole-genome sequencing analysis of lineage 1 strains of Mycobacterium tuberculosis in Da Nang, Vietnam

Background

Spacer oligonucleotide typing (spoligotyping), a widely used, classical genotyping method for Mycobacterium tuberculosis complex (MTBC), is a PCR-based dot-blot hybridization technique to detect the genetic diversity of the direct repeat (DR) region. Of the seven major MTBC lineages in the world, lineage 1 (Indo-Oceanic) mostly corresponds to the East African–Indian (EAI) spoligotype family in East Africa and Southeast Asia.

Objectives

We investigated the genomic features of Vietnamese lineage 1 strains, comparing spoligotype patterns using whole-genome sequencing (WGS) data.

Methods

M. tuberculosis strains isolated in Da Nang, Vietnam were subjected to conventional spoligotyping, followed by WGS analysis using a high-throughput sequencer. Vietnamese lineage 1 strains were further analyzed with other lineage 1 strains obtained from a public database.

Results

Indicating a major spoligotype in Da Nang, 86 (46.2%) of the 186 isolates belonged to the EAI family or lineage 1. Although typical EAI4-VNM strains are characterized by the deletion of spacers 26 and 27, 65 (75.6%) showed ambiguous signals on spacer 26. De novo assembly of the entire DR region and in silico spoligotyping analysis suggested the absence of spacer 26, and direct sequencing revealed that the 17^th spacer sequence not used for conventional typing, was cross-hybridized to the spacer 26 probe. Vietnamese EAI4-VNM, other EAI-like strains, and those showing a non-EAI pattern lacking many spacers formed a monophyletic group separate from other EAI families in the world.

Conclusion

Information about the alignment of spacers in the entire DR region obtained from WGS data provides a clue for the determination of experimentally ambiguous spoligo patterns. WGS data also helped to analyze the hidden relationships between apparently distinct spoligo patterns.

Molecular analysis of pyrazinamide resistance in Mycobacterium tuberculosis in Vietnam highlights the high rate of pyrazinamide resistance-associated mutations in clinical isolates

Pyrazinamide (PZA) is a key antibiotic in current anti-tuberculosis regimens. Although the WHO has stressed the urgent need to obtain data on PZA resistance, in high tuberculosis burden countries, little is known about the level of PZA resistance, the genetic basis of such resistance or its link with Mycobacterium tuberculosis families. In this context, this study assessed PZA resistance through the molecular analysis of 260 Vietnamese M. tuberculosis isolates. First-line drug susceptibility testing, pncA gene sequencing, spoligotyping and mycobacterial interspersed repetitive units-variable number of tandem repeats (MIRU-VNTR) typing were performed. Overall, the pncA mutation frequency was 38.1% (99 out of 260 isolates) but was higher than 72% (89 out of 123 isolates) in multidrug and quadruple-drug resistant isolates. Many different pncA mutations (71 types) were detected, of which 55 have been previously described and 50 were linked to PZA resistance. Among the 16 novel mutations, 14 are likely to be linked to PZA resistance because of their mutation types or codon positions. Genotype analysis revealed that PZA resistance can emerge in any M. tuberculosis cluster or family, although the mutation frequency was the highest in Beijing family isolates (47.7%, 62 out of 130 isolates). These data highlight the high rate of PZA resistance-associated mutations in M. tuberculosis clinical isolates in Vietnam and bring into question the use of PZA for current and future treatment regimens of multidrug-resistant tuberculosis without PZA resistance testing.

Combining molecular typing and spatial pattern analysis to identify areas of high tuberculosis transmission in a moderate-incidence county in Taiwan

In total, 303 randomly selected clinical Mycobacterium tuberculosis (MTB) isolates from 303 patients (collected January to December 2012) in central Taiwan were examined. The major lineages found were Beijing (N = 114, 37.62%), Haarlem (N = 76, 25.08%) and East African-Indian (EAI) (N = 42, 13.86%). Notably, younger persons (≤30 years old) were 6.58 times more likely to be infected with a Beijing genotype compared to older persons (>70 years) (p < 0.05). Combining molecular typing methods and geographical information system (GIS) analysis, we uncovered a twofold higher incidence of Beijing strains in a hotspot area (33%) compared to non-hotspot areas (17%). By 24 MIRU-VNTR typing, persons in clustered groups were 1.96 times more likely to be infected with a Beijing strain compared with non-clustered persons, suggesting recent spread and emergence of MTB. Finally, we observed a trend in which TB incidence increased as the density/concentration of analyzed environmental factors increased, suggesting that environmental factors are associated with TB transmission; however, only population density was found to be significantly associated with increased risk of TB (p < 0.05). Molecular typing methods combined with spatial analysis suggest possible TB transmission. Early intervention to interrupt transmission may be most effective if targeted to hot zones of TB.