Extensive endemic transmission of multidrug-resistant Mycobacterium tuberculosis in Bhutan: A retrospective genomic-epidemiological study

Despite decreasing overall tuberculosis notifications, the proportion of multidrug-resistant tuberculosis (MDR-TB) cases are increasing in Bhutan. While most MDR-TB cases are diagnosed among patients in the bordering districts and the capital, current diagnostic tests are limited in their ability to differentiate between the recurrent introductions and local transmission. For the first time, we conducted a retrospective genomic-epidemiological study to provide insights into the population structure, genotypic resistance patterns, and explore recent transmission of drug-resistant TB in Bhutan. Whole genome sequencing was performed on randomly selected drug-resistant and drug-sensitive TB isolates from Bhutan, collected between 2018-2022 at Microbiological Diagnostic Unit Public Health Laboratory in Melbourne, Australia. We investigated drug resistance mutations, and genomic clustering of cases using different single nucleotide polymorphisms (SNPs) thresholds. Of the 203 sequences that passed the quality control, 126 (62.1%) were MDR-TB and 15 (7.4%) were isoniazid-resistant TB. There were four different circulating lineages, with most sequences belonging to lineage 2 (86.2%). Using a SNP-threshold of ≤12 SNPs, 71% of sequences formed 12 genomic clusters. Surprisingly, the largest cluster included of 88% of all MDR-TB sequences and spanned the entire study period. These cases were highly clonal (mean pairwise SNP-distance of 10, range 0-25). Phylogenetic analysis with publicly available international sequence data showed that this MDR-TB cluster formed a distinct clade. The major burden of MDR-TB in Bhutan appears to be due to recent local transmission of cases resulting in a large, single endemic cluster. This information will be critical for TB control program in Bhutan to tackle the major burden of MDR-TB through enhanced contact tracing of this MDR-TB clade. Additionally, this genomic data will be valuable to regional neighbours to monitor for dissemination of the strain. This study highlights the significant value of investing in TB genomics in resource limited settings globally to gain actionable insights into transmission dynamics.