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

Analysis of fine-scale phylogeny of Burkholderia pseudomallei in relation to regional geography and drug susceptibility in Thailand

Melioidosis caused by Burkholderia pseudomallei (Bp) is a public health threat. Genomic-epidemiology research on this deadly disease is scarce. We investigated whole-genome sequences of Bp isolates in relation to environmental source and drug susceptibility. In total, 563 Bp isolates were collected from 11 Northeast Thai provinces during the period 2004-2021. Patients (n = 530 isolates), infected animals (n = 8), and environmental sources (n = 25) provided samples. Phylogenetic analysis revealed genetic diversity among the Bp isolates, including numerous well-supported clusters of varying sizes. Through in-depth analysis of 38 monophyletic clades (MCs), we found eleven associated with province of origin (p-value < 0.001). Closely related clusters (CRCs) within MCs resembled MLST-identified "sequence types" (STs). We found 102 known and 52 novel STs. ST-70 was the most prevalent in this area (n = 78; 13.85%). Sample type (human/environmental) and sampling time intervals were not correlated with genetic distance among clonal Bp isolates. Some members of 12 CRCs had acquired resistance to co-trimoxazole and one against amoxicillin-clavulanic acid. Within Northeast Thailand, there is an association between Bp genotype and geographical origin.

Transmission dynamics of tuberculosis in a high-burden area of China: An 8-year population-based study using whole genome sequencing

OBJECTIVES

This study investigated the transmission patterns of tuberculosis (TB) and its associated risk factors in Hunan province to inform the development of prevention and control strategies in the region.

METHODS

An 8-year retrospective population-based genomic epidemiological study was conducted. Genomic clusters were defined using distance thresholds of 12-single-nucletide-polymorphisms. Risk factors associated with TB transmission were analyzed using logistic regression model. Kernel Density analysis was used to locate hotspots where transmission occurred.

RESULTS

Among 2649 TB cases included in this study, 275 clusters were identified, with an overall clustering rate of 24.7% (654/2649). Nearly 95% (620/654) of clustered strains were isolated from the same county. Of the 275 clusters, 23 (8.4%, 23/275) had differences in drug-resistant profiles, with FQs resistance mutations occurring most frequently (52.2%, 12/23). Multivariate analysis identified male TB patients, those aged 30-60 years, ethnic minorities, nonfarmers, retreated TB patients, and individuals infected with MDR/RR-TB as independent risk factors for TB transmission (P < 0.05). Kernel density analysis showed that among the 5 drug-resistant surveillance sites, Leiyang had the highest clustering rate, followed by Yongshun, Qidong, Hecheng, and Taojiang.

CONCLUSION

Recent transmission in the region is predominantly occurring within counties. The risk factors related to TB transmission and the hotspots where transmission occurs can provide a scientific basis for the formulation of targeted TB prevention and control strategies.

Transmission of drug-resistant Mycobacterium tuberculosis isolates between Finnish- and foreign-born cases, 2014-2021: A molecular epidemiological study

BACKGROUND

Data on the molecular epidemiology and transmission of drug-resistant Mycobacterium tuberculosis (MTB) in low-incidence settings with immigration from high-incidence settings is limited.

METHOD

We included 115 drug-resistant (DR) MTB isolates with whole-genome sequencing data isolated in Finland between 2014 and 2021. Potential transmission clusters were identified using a threshold of 12 single-nucleotide polymorphisms (SNPs). Highly related clusters were identified using a threshold of 5 SNPs.

RESULT

Of the 115 DR MTB isolates, 31 (27.0%) isolates were from Finnish-born cases and 84 (73.0%) were from foreign-born cases. The proportion of multidrug-resistant (MDR) MTB isolates (30/84, 35.7%) from foreign-born cases was higher than that of MDR MTB isolates from Finnish-born cases (8/31, 25.8%). Lineage 2 (40/115, 34.8%) and lineage 4 (40/115, 34.8%) were the most prevalent lineages. A total of 25 (21.7%) isolates were classified into eight potential transmission clusters (≤12 SNPs). Furthermore, five highly related clusters (≤5 SNPs) were identified, including three DR MTB isolates from Finnish-born cases and 14 DR isolates from foreign-born cases.

CONCLUSION

The risk of DR MTB transmission between Finnish- and foreign-born persons is not negligible. Further research on clustering analysis in drug-susceptible MTB is worth to inform tuberculosis management and control in low-incidence settings with increasing immigration.

Multidrug-resistant Mycobacterium tuberculosis transmission in Shandong, China

Multidrug-resistant tuberculosis (MDR-TB) has imposed a significant economic and health burden worldwide, notably in China. Using whole genome sequence, we sought to understand the mutation and transmission of MDR-TB in Shandong. A retrospective study of patients diagnosed with pulmonary tuberculosis in Shandong from 2009 to 2018 was conducted. To explore transmission patterns, we performed whole genome sequencing on MDR-TB isolates, identified genomic clusters, and assessed the drug resistance of TB isolates. Our study analyzed 167 isolates of MDR-TB, finding that 100 were clustered. The predominant lineage among MDR-TB isolates was lineage 2, specifically with a notable 88.6% belonging to lineage 2.2.1. Lineage 4 constituted a smaller proportion, accounting for 4.2% of the isolates. We discovered that Shandong has a significant clustering percentage for MDR-TB, with Jining having the highest percentage among all Shandong cities. The clustering percentages of MDR-TB, pre-extensively drug-resistant tuberculosis, and extensively drug-resistant tuberculosis were 59.9%, 66.0%, and 71.4%, respectively, and the clustering percentages increased with the expansion of the anti-TB spectrum. Isolates from genomic clusters 1 and 3 belonged to lineage 2.2.1 and showed signs of cross-regional transmission. The distribution of rrs A1401G and katG S315T mutations in lineage 2.2.1 and 2.2.2 strains differed significantly (P < .05). MDR-TB isolates with rpoB I480V, embA-12C > T, and rrs A1401G mutations showed a higher likelihood of clustering (P < .05). Our findings indicate a significant problem of local transmission of MDR-TB in Shandong, China. Beijing lineage isolates and some drug-resistant mutations account for the MDR-TB transmission in Shandong.

Co-resistance to isoniazid and second-line anti-tuberculosis drugs in isoniazid-resistant tuberculosis at a tertiary care hospital in Thailand

Isoniazid-resistant tuberculosis (Hr-TB) is an important drug-resistant tuberculosis (TB). In addition to rifampicin, resistance to other medications for Hr-TB can impact the course of treatment; however, there are currently limited data in the literature. In this study, the drug susceptibility profiles of Hr-TB treatment and resistance-conferring mutations were investigated for Hr-TB clinical isolates from Thailand. Phenotypic drug susceptibility testing (pDST) and genotypic drug susceptibility testing (gDST) were retrospectively and prospectively investigated using the Mycobacterium Growth Indicator Tube (MGIT), the broth microdilution (BMD) method, and whole-genome sequencing (WGS)-based gDST. The prevalence of Hr-TB cases was 11.2% among patients with TB. Most Hr-TB cases (89.5%) were newly diagnosed patients with TB. In the pDST analysis, approximately 55.6% (60/108) of the tested Hr-TB clinical isolates exhibited high-level isoniazid resistance. In addition, the Hr-TB clinical isolates presented co-resistance to ethambutol (3/161, 1.9%), levofloxacin (2/96, 2.1%), and pyrazinamide (24/118, 20.3%). In 56 Hr-TB clinical isolates, WGS-based gDST predicted resistance to isoniazid [katG S315T (48.2%) and fabG1 c-15t (26.8%)], rifampicin [rpoB L430P and rpoB L452P (5.4%)], and fluoroquinolones [gyrA D94G (1.8%)], but no mutation for ethambutol was detected. The categorical agreement for the detection of resistance to isoniazid, rifampicin, ethambutol, and levofloxacin between WGS-based gDST and the MGIT or the BMD method ranged from 80.4% to 98.2% or 82.1% to 100%, respectively. pDST and gDST demonstrated a low co-resistance rate between isoniazid and second-line TB drugs in Hr-TB clinical isolates.

IMPORTANCE

The prevalence of isoniazid-resistant tuberculosis (Hr-TB) is the highest among other types of drug-resistant tuberculosis. Currently, the World Health Organization (WHO) guidelines recommend the treatment of Hr-TB with rifampicin, ethambutol, pyrazinamide, and levofloxacin for 6 months. The susceptibility profiles of Hr-TB clinical isolates, especially when they are co-resistant to second-line drugs, are critical in the selection of the appropriate treatment regimen to prevent treatment failure. This study highlights the susceptibility profiles of the WHO-recommended treatment regimen in Hr-TB clinical isolates from a tertiary care hospital in Thailand and the concordance and importance of using the phenotypic drug susceptibility testing or genotypic drug susceptibility testing for accurate and comprehensive interpretation of results.

Analysis of Drug-Resistance Characteristics and Genetic Diversity of Multidrug-Resistant Tuberculosis Based on Whole-Genome Sequencing on the Hainan Island, China

Purpose

Given the high burden of Tuberculosis (TB) in China, the prevalence of multidrug-resistant tuberculosis (MDR-TB) is significant. Whole-genome sequencing (WGS) of Mycobacterium tuberculosis (MTB) enables the identification of lineages, drug-resistant mutations, and transmission patterns, offering valuable insights for TB control, clinical diagnosis, and treatment.

Methods

We collected 202 MDR-MTB strains from 3519 suspected pulmonary TB patients treated at The Second Affiliated Hospital of Hainan Medical University between July 2019 and June 2021. Proportional drug-susceptibility testing was performed using 8 common anti-tuberculosis drugs. Subsequently, the genotypic drug resistance and genetic characteristics were analyzed by the WGS.

Results

Lineages are identified by TB-profiler revealed 202 MDR-MTB strains, showcasing three predominant lineages, with lineage 2 being the most prevalent. Close genomic relatedness analysis and evidence of MTB transmission led to the formation of 15 clusters comprising 42 isolates, resulting in a clustering rate of 20.8%. Novelty, lineage 2.1 (non-Beijing) accounted for 27.2% of the MDR-MTB strains, which is rare in China and Neighboring countries. Regarding first-line anti-TB drugs, genes associated with rifampicin resistance, primarily the rpoB gene, were detected in 200 strains (99.0%). Genes conferring resistance to isoniazid, ethambutol, and streptomycin were identified in 191 (94.5%), 125 (61.9%), and 100 (49.5%) strains, respectively. Among the second-line drugs, 97 (48.0%) strains exhibited genes encoding resistance to fluoroquinolones. Comparing the results to phenotypic drug susceptibility-based testing, the sensitivity of WGS for detecting resistance to each of the six drugs (rifampicin, isoniazid, ethambutol, ofloxacin, kanamycin, capreomycin) was 90% or higher. With the exception of ethambutol, the specificity of WGS prediction for the remaining drugs exceeded 88%.

Conclusion

Our study provides crucial insights into genetic mutation types, genetic diversity, and transmission of MDR-MTB on Hainan Island, serving as a significant reference for MDR-MTB surveillance and clinical decision-making.

Whole Genomic Analysis Revealed High Genetic Diversity and Drug-Resistant Characteristics of Mycobacterium tuberculosis in Guangxi, China

Background

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb), is a major public health issue in China. Nevertheless, the prevalence and drug resistance characteristics of isolates vary in different regions and provinces. In this study, we investigated the population structure, transmission dynamics and drug-resistant profiles of Mtb in Guangxi, located on the border of China.

Methods

From February 2016 to April 2017, 462 clinical M. tuberculosis isolates were selected from 5 locations in Guangxi. Drug-susceptibility testing was performed using 6 common anti-tuberculosis drugs. The genotypic drug resistance and transmission dynamics were analyzed by the whole genome sequence.

Results

Our data showed that the Mtb in Guangxi has high genetic diversity including Lineage 1 to Lineage 4, and mostly belong to Lineage 2 and Lineage 4. Novelty, 9.6% of Lineage 2 isolates were proto-Beijing genotype (L2.1), which is rare in China. About 12.6% of isolates were phylogenetically clustered and formed into 28 transmission clusters. We observed that the isolates with the high resistant rate of isoniazid (INH, 21.2%), followed by rifampicin (RIF, 13.2%), and 6.7%, 12.1%, 6.7% and 1.9% isolates were resistant to ethambutol (EMB), streptomycin (SM), ofloxacin (OFL) and kanamycin (KAN), respectively. Among these, 6.5% and 3.3% of isolates belong to MDR-TB and Pre-XDR, respectively, with a high drug-resistant burden. Genetic analysis identified the most frequently encountered mutations of INH, RIF, EMB, SM, OFL and KAN were katG_Ser315Thr (62.2%), rpoB_Ser450Leu (42.6%), embB_Met306Vol (45.2%), rpsL_Lys43Arg (53.6%), gyrA_Asp94Gly (29.0%) and rrs_A1401G (66.7%), respectively. Additionally, we discovered that isolates from border cities are more likely to be drug-resistant than isolates from non-border cities.

Conclusion

Our findings provide a deep analysis of the genomic population characteristics and drug-resistant of M. tuberculosis in Guangxi, which could contribute to developing effective TB prevention and control strategies.

Metabolomic analysis of Mycobacterium tuberculosis reveals metabolic profiles for identification of drug-resistant tuberculosis

The detection of pre-extensively (pre-XDR) and extensively drug-resistant tuberculosis (XDR-TB) is challenging. Drug-susceptibility tests for some anti-TB drugs, especially ethambutol (ETH) and ethionamide (ETO), are problematic due to overlapping thresholds to differentiate between susceptible and resistant phenotypes. We aimed to identify possible metabolomic markers to detect Mycobacterium tuberculosis (Mtb) strains causing pre-XDR and XDR-TB. The metabolic patterns of ETH- and ETO-resistant Mtb isolates were also investigated. Metabolomics of 150 Mtb isolates (54 pre-XDR, 63 XDR-TB and 33 pan-susceptible; pan-S) were investigated. Metabolomics of ETH and ETO phenotypically resistant subgroups were analyzed using UHPLC-ESI-QTOF-MS/MS. Orthogonal partial least-squares discriminant analysis revealed distinct separation in all pairwise comparisons among groups. Two metabolites (meso-hydroxyheme and itaconic anhydride) were able to differentiate the pre-XDR and XDR-TB groups from the pan-S group with 100% sensitivity and 100% specificity. In comparisons of the ETH and ETO phenotypically resistant subsets, sets of increased (ETH = 15, ETO = 7) and decreased (ETH = 1, ETO = 6) metabolites specific for the resistance phenotype of each drug were found. We demonstrated the potential for metabolomics of Mtb to differentiate among types of DR-TB as well as between isolates that were phenotypically resistant to ETO and ETH. Thus, metabolomics might be further applied for DR-TB diagnosis and patient management.

Discrepancy in the transmissibility of multidrug-resistant Mycobacterium tuberculosis in urban and rural areas in China

The fitness of multidrug-resistant tuberculosis (MDR-TB) is thought to be an important determinant of a strain's ability to be transmitted and cause outbreaks. Studies in the laboratory have demonstrated that MDR-TB strains have reduced fitness but the relative transmissibility of MDR-TB versus drug-susceptible (DS) TB strains in human populations remains unresolved. We used data on genomic clustering from our previous molecular epidemiological study in Songjiang (2011-2020) and Wusheng (2009-2020), China, to compare the relative transmissibility of MDR-TB versus DS-TB. Genomic clusters were defined with a threshold distance of 12-single-nucleotide-polymorphisms and the risk for MDR-TB clustering was analyzed by logistic regression. In total, 2212 culture-positive pulmonary TB patients were enrolled in Songjiang and 1289 in Wusheng. The clustering rates of MDR-TB and DS-TB strains were 19.4% (20/103) and 26.3% (509/1936), respectively in Songjiang, and 43.9% (29/66) and 26.0% (293/1128) in Wusheng. The risk of MDR-TB clustering was 2.34 (95% CI 1.38-3.94) times higher than DS-TB clustering in Wusheng and 0.64 (95% CI 0.38-1.06) times lower in Songjiang. Neither lineage 2, compensatory mutations nor rpoB S450L were significantly associated with MDR-TB transmission, and katG S315T increased MDR-TB transmission only in Wusheng (OR 5.28, 95% CI 1.42-19.21). MDR-TB was not more transmissible than DS-TB in either Songjiang or Wusheng. It appears that the different transmissibility of MDR-TB in Songjiang and Wusheng is likely due to differences in the quality of the local TB control programs. These results suggest that the most effective way to control MDR-TB is by improving local TB control programs.

Molecular Characterization of Mutations in Isoniazid- and Rifampicin-Resistant Mycobacterium tuberculosis Isolated in Thailand

The control of drug-resistant tuberculosis (TB) is a major challenge. The frequency and mutation characteristics indicate the efficiency of molecular tests for the rapid detection of TB drug resistance. This study examined the existence of katG and inhA mutations for isoniazid (INH) resistance and rpoB mutations for rifampicin (RFP) resistance. In total, 178 drug-resistant Mycobacterium tuberculosis (MTB) isolates were analyzed. Mutations in katG encoding and inhA regulatory regions were detected in 136/168 (81.0%) and 29/168 (17.3%), respectively, with the most prominent mutation of Ser315Thr substitution in katG in 126/168 (75.0%), and -15 C to T substitution in the regulatory region of the inhA (26/168; 15.5%). Two distinct katG mutations (Tyr337Cys, 1003InsG) were identified. Of 125 RFP-resistant isolates, 118 (94.4%) carried mutations affecting the 81-bp RFP resistance-determining region, with the most commonly affected codons 450, 445, and 435 identified in 74 (59.2%), 26 (20.8%), and 12 (9.6%) isolates, respectively. Genetic mutations were highly associated with phenotypic INH and RFP resistance, and the majority shared similarities with those reported in previous studies in Thailand and other Asian countries. These data are useful for guiding the use and improvement of molecular tests for TB drug resistance.

High clustering rate and genotypic drug-susceptibility screening for the newly recommended anti-tuberculosis drugs among global extensively drug-resistant Mycobacterium tuberculosis isolates

Multidrug-resistant tuberculosis (MDR-TB) and extensively drug-resistant tuberculosis (XDR-TB) make TB difficult to control. Global susceptibility data for six newly recommended anti-TB drugs against M/XDR-TB are still limited. Using publicly available whole-genome sequences, we determined the proportion of 513 phenotypically XDR-TB isolates that carried mutations associated with resistance against these drugs (bedaquiline, clofazimine, linezolid, delamanid, pretomanid and cycloserine). Mutations of Rv0678 and Rv1979c were detected in 69/513 isolates (13.5%) for bedaquiline resistance and 79/513 isolates (15.4%) for clofazimine resistance with additional mmpL5 mutations. Mutations conferring resistance to delamanid were detected in fbiB and ddn genes for 11/513 isolates (2.1%). For pretomanid, a mutation was detected in the ddn gene for 3/513 isolates (0.6%). Nineteen mutations of pykA, cycA, ald, and alr genes, conferring resistance to cycloserine, were found in 153/513 isolates (29.8%). No known mutations associated with linezolid resistance were detected. Cluster analysis showed that 408/513 isolates fell within 99 clusters and that 354 of these isolates were possible primary drug-resistant TB (292 XDR-TB, 57 pre-XDR-TB and 5 MDR-TB). Clonal transmission of primary XDR isolates might contribute significantly to the high prevalence of DR-TB globally.

Whole genome sequencing of multidrug-resistant Mycobacterium tuberculosis isolates collected in the Czech Republic, 2005-2020

The emergence and spread of resistant tuberculosis (TB) pose a threat to public health, so it is necessary to diagnose the drug-resistant forms in a clinically short time frame and closely monitor their transmission. In this study, we carried out a first whole genome sequencing (WGS)-based analysis of multidrug resistant (MDR) M. tuberculosis strains to explore the phylogenetic lineages diversity, drug resistance mechanisms, and ongoing transmission chains within the country. In total, 65 isolates phenotypically resistant to at least rifampicin and isoniazid collected in the Czech Republic in 2005-2020 were enrolled for further analysis. The agreement of the results obtained by WGS with phenotypic drug susceptibility testing (pDST) in the determination of resistance to isoniazid, rifampicin, pyrazinamide, streptomycin, second-line injectables and fluoroquinolones was more than 80%. Phylogenetic analysis of WGS data revealed that the majority of MDR M. tuberculosis isolates were the Beijing lineage 2.2.1 (n = 46/65; 70.8%), while the remaining strains belonged to Euro-American lineage. Cluster analysis with a predefined cut-off distance of less than 12 single nucleotide polymorphisms between isolates showed 19 isolates in 6 clusters (clustering rate 29.2%), located mainly in the region of the capital city of Prague. This study highlights the utility of WGS as a high-resolution approach in the diagnosis, characterization of resistance patterns, and molecular-epidemiological analysis of resistant TB in the country.

Drug-Resistant Characteristics, Genetic Diversity, and Transmission Dynamics of Rifampicin-Resistant Mycobacterium tuberculosis in Hunan, China, Revealed by Whole-Genome Sequencing

To gain a deep insight into the additional drug-resistant profiles, genetic diversity, and transmission dynamics of rifampicin-resistant tuberculosis (RR-TB) circulating in Hunan province, drug susceptibility testing and whole-genome-sequencing were performed among RR-TB strains collected from Jan. 2013 to Jun. 2018 in Hunan province. A total of 124 RR-TB strains were recovered successfully and included into the final analysis. Lineage 2.2.1 was the dominant sublineage, accounting for 72.6% (90/124), followed by lineage 4.5 (11.3%, 14/124), lineage 4.4 (8.1%, 10/124), lineage 4.2 (6.5%, 8/124) and lineage 2.2.2 (1.6%, 2/124). Overall, 83.1% (103/124) and 3.2% (4/124) of RR-TB were MDR-TB and XDR-TB, respectively. Nearly 30% of RR-TB isolates were resistant to fluoroquinolones, and 26.6% (33/124) were pre-XDR-TB. Moreover, 30.6% (38/124) of RR-TB strains were identified as phenotypically resistance to pyrazinamide. Totally, 17 clusters containing 48 (38.7%, 48/124) RR-TB strains were identified, ranging in size from 2 to 10 isolates. No significant difference was detected in clustering rate between lineage 2 and lineage 4 (χ² = 0.027, P = 0.870). Our study revealed the complexity of RR-TB strains circulating in Hunan province with complex additional drug-resistant profile and relatively higher clustering rates. Comprehensive information based on WGS should be used to guide the design of treatment regimens and tailor public interventions. IMPORTANCE Comprehensive information such as genetic background and drug-resistant profile of MTB strains could help to tailor public interventions. However, these data are limited in Hunan province, one of the provinces with high-TB burden in China. So, this study aimed to provide us with deep insight into the molecular epidemiology of RR-TB isolates circulating in Hunan province by combining phenotypic drug susceptibility testing and whole-genome sequencing. To our knowledge, this is the first study to use whole-genome sequencing data of RR-TB strains spanning more than 5 years for molecular epidemiology analysis in Hunan province, which allows us to identify genetic background information and clustered strains more accurately. Our study revealed the complexity of RR-TB strains circulating in Hunan province with complex additional drug-resistant profile and relatively higher clustering rates. Comprehensive information based on WGS should be used to guide the design of treatment regimens and tailor public interventions.

Molecular Epidemiology of Drug-Resistant Mycobacterium Tuberculosis in Japan

Clinical isolates of drug-resistant (isoniazid and/or rifampicin-resistant) Mycobacterium tuberculosis were obtained from 254 patients diagnosed with drug-resistant tuberculosis in Japan from April 2015 to March 2017 in National Hospital Organization hospitals. The 254 patients were approximately 32% of all 795 patients who were diagnosed with culture-confirmed drug-resistant tuberculosis from 2015 to 2016 nationwide in Japan. The whole-genome sequences of all the isolates from the 254 patients and the lineages of these isolates were determined, and phylogenetic trees were constructed based on single nucleotide polymorphism concatemers. Of these patients, 202 (79.5%) were born in Japan and 52 (20.5%) were born elsewhere. Of the 254 drug-resistant isolates, 54 (21.3%) were multidrug resistant, being resistant to both isoniazid and rifampicin. The percentages of multidrug-resistant isolates were significantly higher in foreign-born (38.5% [20/52]) than Japanese-born patients (16.8% [34/202]). Of the 54 multidrug-resistant isolates, nine were extensively drug resistant, which were all obtained from Japanese-born patients. Five extensively drug-resistant isolates were obtained from patients with incipient tuberculosis. A significant number of multidrug-resistant M. tuberculosis strains were isolated from foreign-born patients from Asian countries that have a high tuberculosis burden. Foreign-derived isolates affect the nationwide genetic diversity of drug-resistant M. tuberculosis in Japan. Extensively drug-resistant M. tuberculosis isolates were transmitted among the Japanese population.

IMPORTANCE The incidence rate of tuberculosis (TB) in Japan was 11.5 per 100,000 of the population in 2019. Of TB patients in Japan, 61.1% were aged >70 years, and 10.7% were born outside Japan, mostly in Asian countries with a high burden of tuberculosis. Of the tuberculosis patients in the present study, 5.4% and 1.0% showed resistance to isoniazid and rifampicin, respectively, and 0.7% were multidrug resistant. The objective of this study was to clarify the molecular epidemiological properties of drug-resistant tuberculosis in Japan. Molecular epidemiology provides several clues to inform potential measures to control drug-resistant tuberculosis in Japan.