Whole genome sequencing analysis of multidrug-resistant tuberculosis in Singapore, 2006-2018

Phylodynamic assessment of SNP distances from whole genome sequencing for determining Mycobacterium tuberculosis transmission

The global tuberculosis (TB) epidemic is driven by primary transmission. Pathogen genome sequencing is increasingly used in molecular epidemiology and outbreak investigations. Based on contact tracing and epidemiological links, Single Nucleotide Polymorphism (SNP) cut-offs, ranging from 3 to 12 SNPs, identify probable transmission clusters or exclude direct transmission. However, contact tracing can be limited by recall bias and inconsistent methodologies across TB settings. We propose phylodynamic models, i.e. methods to infer transmission processes from pathogen genomes and associated epidemiological data, as an alternative reference to infer transmission events. We analyzed 2,008 whole-genome sequences from Dutch TB patients collected from 2015 to 2019. Genetic clusters were defined within a 20-SNP range, and the phylodynamic model phybreak was employed to infer transmission. Probable transmission SNP cut-offs were assessed by the proportion of inferred transmission events with a SNP distance below these cut-offs. A total of 79 clusters were identified, with a median size of 4 isolates (IQR = 3-8). A SNP cut-off of 4 captured 98% of inferred transmission events while reducing pairs without transmission links. A cut-off beyond 12 SNPs effectively excluded transmission. Phylodynamic approaches provide a valuable alternative to contact tracing for defining SNP cut-offs, allowing for a more precise assessment of transmission events.

Implementation of national whole-genome sequencing of Mycobacterium tuberculosis, National Public Health Laboratory, Singapore, 2019-2022

The National Tuberculosis Programme (NTBP) monitors the occurrence and spread of tuberculosis (TB) and multidrug-resistant TB (MDR-TB) in Singapore. Since 2020, whole-genome sequencing (WGS) of Mycobacterium tuberculosis isolates has been performed at the National Public Health Laboratory (NPHL) for genomic surveillance, replacing spoligotyping and mycobacterial interspersed repetitive unit-variable number tandem repeats analysis (MIRU-VNTR). Four thousand three hundred and seven samples were sequenced from 2014 to January 2023, initially as research projects and later developed into a comprehensive public health surveillance programme. Currently, all newly diagnosed culture-positive cases of TB in Singapore are prospectively sent for WGS, which is used to perform lineage classification, predict drug resistance profiles and infer genetic relationships between TB isolates. This paper describes NPHL's operational and technical experiences with implementing the WGS service in an urban TB-endemic setting, focusing on cluster detection and genomic drug susceptibility testing (DST). Cluster detection: WGS has been used to guide contact tracing by detecting clusters and discovering unknown transmission networks. Examples have been clusters in a daycare centre, housing apartment blocks and a horse-racing betting centre. Genomic DST: genomic DST prediction (gDST) identifies mutations in core genes known to be associated with TB drug resistance catalogued in the TBProfiler drug resistance mutation database. Mutations are reported with confidence scores according to a standardized approach referencing NPHL's internal gDST confidence database, which is adapted from the World Health Organization (WHO) TB drug mutation catalogue. Phenotypic-genomic concordance was observed for the first-line drugs ranging from 2959/2998 (98.7 %) (ethambutol) to 2983/2996 (99.6 %) (rifampicin). Aspects of internal database management, reporting standards and caveats in results interpretation are discussed.

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.

The mutation rate of rpoB gene showed an upward trend with the increase of MIRU10, MIRU39 and QUB4156 repetitive number

BACKGROUND

Mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) is a frequently used typing method for identifying the Beijing genotype of Mycobacterium tuberculosis (Mtb), which is easily transformed into rifampicin (RIF) resistance. The RIF resistance of Mtb is considered to be highly related with the mutation of rpoB gene. Therefore, this study aimed to analyze the relationship between the repetitive number of MIRU loci and the mutation of rpoB gene.

METHODS

An open-source whole-genome sequencing data of Mtb was used to detect the mutation of rpoB gene and the repetitive number of MIRU loci by bioinformatics methods. Cochran-Armitage analysis was performed to analyze the trend of the rpoB gene mutation rate and the repetitive number of MIRU loci.

RESULTS

Among 357 rifampicin-resistant tuberculosis (RR-TB), 304 strains with mutated rpoB genes were detected, and 6 of 67 rifampicin susceptible strains were detected mutations. The rpoB gene mutational rate showed an upward trend with the increase of MIRU10, MIRU39, QUB4156 and MIRU16 repetitive number, but only the repetitive number of MIRU10, MRIU39 and QUB4156 were risk factors for rpoB gene mutation. The Hunter-Gaston discriminatory index (HGDI) of MIRU10 (0.65) and QUB4156 (0.62) was high in the overall sample, while MIRU39 (0.39) and MIRU16 (0.43) showed a moderate discriminatory Power.

CONCLUSION

The mutation rate of rpoB gene increases with the addition of repetitive numbers of MIRU10, QUB4156 and MIRU39 loci.

Quantifying Mycobacterium tuberculosis Transmission Dynamics Across Global Settings: A Systematic Analysis

The degree to which individual heterogeneity in the production of secondary cases ("superspreading") affects tuberculosis (TB) transmission has not been systematically studied. We searched for population-based or surveillance studies in which whole genome sequencing was used to estimate TB transmission and in which the size distributions of putative TB transmission clusters were enumerated. We fitted cluster-size-distribution data to a negative binomial branching process model to jointly infer the transmission parameters $R$ (the reproduction number) and the dispersion parameter, $k$, which quantifies the propensity of superspreading in a population (generally, lower values of $k$ ($<1.0$) suggest increased heterogeneity). Of 4,796 citations identified in our initial search, 9 studies from 8 global settings met the inclusion criteria (n = 5 studies of all TB; n = 4 studies of drug-resistant TB). Estimated $R$ values (range, 0.10-0.73) were below 1.0, consistent with declining epidemics in the included settings; estimated $k$ values were well below 1.0 (range, 0.02-0.48), indicating the presence of substantial individual-level heterogeneity in transmission across all settings. We estimated that a minority of cases (range, 2%-31%) drive the majority (80%) of ongoing TB transmission at the population level. Identifying sources of heterogeneity and accounting for them in TB control may have a considerable impact on mitigating TB transmission.

A university-clustered tuberculosis outbreak during the COVID-19 pandemic in eastern China

During the COVID-19 pandemic in 2020, a tuberculosis outbreak occurred in a university in eastern China, with 4,488 students and 421 staff on the campus. A 19-year-old student was diagnosed in August 2019. Later, the first round of screening was initiated among close contacts, but no active cases were found. Till September 2020, four rounds of screening were performed. Four rounds of screening were conducted on September 9, November 8, November 22-25 in 2019 and September 2020, with 0, 5, 0 and 43 cases identified, respectively. A total of 66 active tuberculosis were found in the same university, including 4 sputum culture-positive and 7 sputum smear-positive. The total attack rate of active tuberculosis was 1.34% (66/4909). The whole-genome sequencing showed that the isolates belonged to the same L2 sub-specie and were sensitive to all tested antituberculosis drugs. Delay detection, diagnosis and report of cases were the major cause of this university tuberculosis epidemic. More attention should be paid to the asymptomatic students in the index class. After the occurrence of tuberculosis cases in schools, multiple rounds of screening should be carried out, and preventive therapy should be applied in a timely manner.

Transmission of multidrug-resistant Mycobacterium tuberculosis in Wuhan, China: A retrospective molecular epidemiological study

How multidrug-resistant tuberculosis (MDR-TB) spreads and expands in Wuhan population is not clear. The study aimed to determine the transmission patterns of MDR-TB in Wuhan city, China, including 149 patients with MDR-TB.Tuberculosis isolates were genotyped by deletion-targeted multiplex polymerase chain reaction, mycobacterial interspersed repetitive unit-variable number tandem repeat typing, and sequencing of drug resistance-associated genes. The risk factors of genomic-clustering were analyzed with logistic regression. The genomic-clustering patients were deeply investigated.The analysis identified 111 unique and 11 clustered genotypes (38 isolates). The clustering rate was 25.50% and the minimum estimate proportion of recent transmission was 18.12%. Two clusters (5 isolates) shared the same mutation, the remain 9 clusters (33 isolates) had different mutation. Logistic regression showed that older than 60 years (adjusted OR 2.360, 95% CI:1.052-5.292) was an independent factor associated with the genomic-clustering of MDR-TB. Among the 38 genomic-clustering cases, 14 cases had epidemiological transmission links. The most common type of transmission link was social contact.The local transmission of MDR-TB in Wuhan was really an issue. The elderly population might be the high-risk groups for transmission of MDR-TB, and the community or public transportation might be the main transmission places.

Genomic-based surveillance reveals high ongoing transmission of multi-drug-resistant Mycobacterium tuberculosis in Southern Brazil

Genomic-based surveillance on the occurrence of drug resistance and its transmission dynamics has emerged as a powerful tool for the control of tuberculosis (TB). A whole-genome sequencing approach, phenotypic testing and clinical-epidemiological investigation were used to undertake a retrospective population-based study on drug-resistant (DR)-TB in Rio Grande do Sul, the largest state in Southern Brazil. The analysis included 305 resistant Mycobacterium tuberculosis strains sampled statewide from 2011 to 2014, and covered 75.7% of all DR-TB cases identified in this period. Lineage 4 was found to be predominant (99.3%), with high sublineage-level diversity composed mainly of 4.3.4.2 [Latin American and Mediterranean (LAM)/RD174], 4.3.3 (LAM/RD115) and 4.1.2.1 (Haarlem/RD182) sublineages. Genomic diversity was also reflected in resistance of the variants to first-line drugs. A large number of distinct resistance-conferring mutations, including variants that have not been reported previously in any other setting worldwide, and 22 isoniazid-monoresistant strains with mutations described as disputed in the rpoB gene but causing rifampicin resistance generally missed by automated phenotypic tests as BACTEC MGIT. Using a cut-off of five single nucleotide polymorphisms, the estimated recent transmission rate was 55.1%, with 168 strains grouped into 28 genomic clusters. The most worrying fact concerns multi-drug-resistant (MDR) strains, of which 73.4% were clustered. Different resistance profiles and acquisition of novel mutations intraclusters revealed important amplification of resistance in the region. This study described the diversity of M. tuberculosis strains, the basis of drug resistance, and ongoing transmission dynamics across the largest state in Southern Brazil, stressing the urgent need for MDR-TB transmission control state-wide.

Recurrent Tuberculosis Disease in Singapore

Background

Previously treated (ie, recurrent) tuberculosis (TB) cases account for approximately 7%-8% of incident TB globally and in Singapore. Molecular fingerprinting has enabled the differentiation of these patients into relapsed or reinfection cases.

Methods

Patient demographics, disease characteristics, and treatment information were obtained from the national TB notification registry and TB Control Unit. We performed a retrospective, case-control study to evaluate factors associated with recurrent TB disease in Singapore citizens and permanent residents with culture-positive TB from 2006 to 2013 and who developed a second episode of culture-positive TB up to 2016 using multivariable logistic regression analyses.

Results

Ninety-one cases with culture-positive first and recurrent TB disease episodes were identified. Recurrent TB was associated with age ≥60 years (adjusted odds ratio [aOR], 1.98 [95% confidence interval {CI}, 1.09-3.61), male sex (aOR, 2.29 [95% CI, 1.22-4.51]), having concomitant pulmonary and extrapulmonary TB (aOR, 3.10 [95% CI, 1.59-6.10]) and extrapulmonary TB alone (aOR, 3.82 [95% CI, 1.12-13.31]), and was less likely in non-Malays (aOR, 0.52 [95% CI, .27-.99]). DNA fingerprinting results for both episodes in 49 cases differentiated these into 28 relapsed and 21 reinfection cases. Relapse was associated with having concomitant pulmonary and extrapulmonary TB (aOR, 9.24 [95% CI, 2.50-42.42]) and positive sputum acid-fast bacilli smear (aOR, 3.95 [95% CI, 1.36-13.10]).

Conclusions

Relapse and reinfection contributed to 57% and 43%, respectively, of recurrent TB in Singapore. Our study highlights the underappreciated association of concomitant pulmonary and extrapulmonary TB as a significant risk factor for disease relapse.