Abstract:Background: Tuberculosis remains a global health threat, and the World Health Organization reports a limited reduction in disease incidence rates, including both new and relapse cases. Therefore, studies targeting tuberculosis transmission chains and recurrent episodes are crucial for developing the most effective control measures. Herein, multiple tuberculosis clusters were retrospectively investigated by integrating patients' epidemiological and clinical information with median-joining networks recreated based on whole genome sequencing (WGS) data of Mycobacterium tuberculosis isolates. Methods: Epidemiologically linked tuberculosis patient clusters were identified during the source case investigation for pediatric tuberculosis patients. Only M. tuberculosis isolate DNA samples with previously determined spoligotypes identical within clusters were subjected to WGS and further median-joining network recreation. Relevant clinical and epidemiological data were obtained from patient medical records. Results: We investigated 18 clusters comprising 100 active tuberculosis patients 29 of whom were children at the time of diagnosis; nine patients experienced recurrent episodes. M. tuberculosis isolates of studied clusters belonged to Lineages 2 (sub-lineage 2.2.1) and 4 (sub-lineages 4.3.3, 4.1.2.1, 4.8, and 4.2.1), while sub-lineage 4.3.3 (LAM) was the most abundant. Isolates of six clusters were drug-resistant. Within clusters, the maximum genetic distance between closely related isolates was only 5–11 single nucleotide variants (SNVs). Recreated median-joining networks, integrated with patients' diagnoses, specimen collection dates, sputum smear microscopy, and epidemiological investigation results indicated transmission directions within clusters and long periods of latent infection. It also facilitated the identification of potential infection sources for pediatric patients and recurrent active tuberculosis episodes refuting the reactivation possibility despite the small genetic distance of ≤5 SNVs between isolates. However, unidentified active tuberculosis cases within the cluster, the variable mycobacterial mutation rate in dormant and active states, and low M. tuberculosis genetic variability inferred precise transmission chain delineation. In some cases, heterozygous SNVs with an allelic frequency of 10–73% proved valuable in identifying direct transmission events. Conclusion: The complex approach of integrating tuberculosis cluster WGS-data-based median-joining networks with relevant epidemiological and clinical data proved valuable in delineating epidemiologically linked patient transmission chains and deciphering causes of recurrent tuberculosis episodes within clusters.

Characteristic SNPs defining the major multidrug-resistant Mycobacterium tuberculosis clusters identified by EuSeqMyTB to support routine surveillance, EU/EEA, 2017 to 2019

Emerging Challenges of Whole-Genome-sequencing–powered Epidemiological Surveillance of Globally Distributed Clonal Groups of Bacterial Infections, Giving Acinetobacter Baumannii ST195 As an Example

Towards standardisation: comparison of five whole genome sequencing (WGS) analysis pipelines for detection of epidemiologically linked tuberculosis cases

Define SNP thresholds for delineation of tuberculosis transmissions using whole-genome sequencing

Whole genome sequencing identifies bacterial factors affecting transmission of multidrug-resistant tuberculosis in a high-prevalence setting

Whole-genome Sequencing to Detect Recent Transmission of Mycobacterium Tuberculosis in Settings with a High Burden of Tuberculosis

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

Multidrug-resistant tuberculosis clusters and transmission in Taiwan: a population-based cohort study

Genomic epidemiology of multidrug-resistant Mycobacterium tuberculosis during transcontinental spread

Integrating standardized whole genome sequence analysis with a global Mycobacterium tuberculosis antibiotic resistance knowledgebase

Detection of multidrug-resistance in Mycobacterium tuberculosis by phenotype- and molecular-based assays

A Bioinformatics Whole-Genome Sequencing Workflow for Clinical Mycobacterium tuberculosis Complex Isolate Analysis, Validated Using a Reference Collection Extensively Characterized with Conventional Methods and In Silico Approaches

Comprehensive and accurate genetic variant identification from contaminated and low-coverage Mycobacterium tuberculosis whole genome sequencing data

[Application of Whole Genome Sequencing Technology in the Epidemiology of Tuberculosis].

Whole genome sequencing for the prediction of resistant tuberculosis strains from northern India

Whole-genome Sequencing of Mycobacterium Tuberculosis for Prediction of Drug Resistance.

Evaluation of whole-genome sequence data analysis approaches for short- and long-read sequencing of Mycobacterium tuberculosis

Detection of low-frequency resistance-mediating SNPs in next-generation sequencing data of Mycobacterium tuberculosis complex strains with binoSNP

Clinical use of whole genome sequencing for Mycobacterium tuberculosis

TBtypeR: Sensitive detection and sublineage classification of low-frequency Mycobacterium tuberculosis complex mixed infections

Unraveling tuberculosis patient cluster transmission chains: integrating WGS-based network with clinical and epidemiological insights