Transmission dynamics study of tuberculosis isolates with whole genome sequencing in southern Sweden

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.

Impact of Whole-Genome Sequencing of Mycobacterium tuberculosis on Treatment Outcomes for MDR-TB/XDR-TB: A Systematic Review

The emergence and persistence of drug-resistant tuberculosis is a major threat to global public health. Our objective was to assess the applicability of whole-genome sequencing (WGS) to detect genomic markers of drug resistance and explore their association with treatment outcomes for multidrug-resistant/extensively drug-resistant tuberculosis (MDR/XDR-TB).

METHODS

Five electronic databases were searched for studies published in English from the year 2000 onward. Two reviewers independently conducted the article screening, relevant data extraction, and quality assessment. The data of the included studies were synthesized with a narrative method and are presented in a tabular format.

RESULTS

The database search identified 949 published articles and 8 studies were included. An unfavorable treatment outcome was reported for 26.6% (488/1834) of TB cases, which ranged from 9.7 to 51.3%. Death was reported in 10.5% (194/1834) of total cases. High-level fluoroquinolone resistance (due to gyrA 94AAC and 94GGC mutations) was correlated as the cause of unfavorable treatment outcomes and reported in three studies. Other drug resistance mutations, like kanamycin high-level resistance mutations (rrs 1401G), rpoB Ile491Phe, and ethA mutations, conferring prothionamide resistance were also reported. The secondary findings from this systematic review involved laboratory aspects of WGS, including correlations with phenotypic DST, cost, and turnaround time, or the impact of WGS results on public health actions, such as determining transmission events within outbreaks.

CONCLUSIONS

WGS has a significant capacity to provide accurate and comprehensive drug resistance data for MDR/XDR-TB, which can inform personalized drug therapy to optimize treatment outcomes.

Genomic transmission clusters and circulating lineages of Mycobacterium tuberculosis among refugees residing in refugee camps in Ethiopia

BACKGROUND

Understanding the transmission dynamics of Mycobacterium tuberculosis (Mtb) could benefit the design of tuberculosis (TB) prevention and control strategies for refugee populations. Whole Genome Sequencing (WGS) has not yet been used to document the Mtb transmission dynamics among refugees in Ethiopia. We applied WGS to accurately identify transmission clusters and Mtb lineages among TB cases in refugee camps in Ethiopia.

METHOD AND DESIGN

We conducted a cross-sectional study of 610 refugees in refugee camps in Ethiopia presenting with symptoms of TB. WGS data of 67 isolates was analyzed using the Maximum Accessible Genome for Mtb Analysis (MAGMA) pipeline; iTol and FigTree were used to visualize phylogenetic trees, lineages, and the presence of transmission clusters.

RESULTS

Mtb culture-positive refugees originated from South Sudan (52/67, 77.6%), Somalia (9/67, 13.4%). Eritrea (4/67, 6%), and Sudan (2/67, 3%). The majority (52, 77.6%) of the isolates belonged to Mtb lineage (L) 3, and one L9 was identified from a Somalian refugee. The vast majority (82%) of the isolates were pan-susceptible Mtb, and none were multi-drug-resistant (MDR)-TB. Based on the 5-single nucleotide polymorphisms cutoff, we identified eight potential transmission clusters containing 23.9% of the isolates. Contact investigation confirmed epidemiological links with either family or social interaction within the refugee camps or with neighboring refugee camps.

CONCLUSION

Four lineages (L1, L3, L4, and L9) were identified, with the majority of strains being L3, reflecting the Mtb L3 dominance in South Sudan, where the majority of refugees originated from. Recent transmission among refugees was relatively low (24%), likely due to the short study period. The improved understanding of the Mtb transmission dynamics using WGS in refugee camps could assist in designing effective TB control programs for refugees.

Whole-genome sequencing of carbapenem-resistant Enterobacterales isolates in southeast Louisiana reveals persistent genetic clusters spanning multiple locations

BACKGROUND

We investigated 51 g-negative carbapenem-resistant Enterobacterales (CRE) isolates collected from 22 patients over a five-year period from six health care institutions in the Ochsner Health network in southeast Louisiana.

METHODS

Short genomic reads were generated using Illumina sequencing and assembled for each isolate. Isolates were classified as Enterobacter spp. (n = 20), Klebsiella spp. (n = 30), and Escherichia coli (n = 1) and grouped into 19 different multi-locus sequence types (MLST). Species and patient-specific core genomes were constructed representing ∼50% of the chromosomal genome.

RESULTS

We identified two sets of patients with genetically related infections; in both cases, the related isolates were collected > 6 months apart, and in one case, the isolates were collected in different locations. On the other hand, we identified four sets of patients with isolates of the same species collected within 21 days from the same location; however, none had genetically related infections. Genes associated with resistance to carbapenem drugs (blaKPC and/or blaCTX-M-15) were found in 76% of the isolates. We found three blaKPC variants (blaKPC-2, blaKPC-3, and blaKPC-4) associated with four different Enterobacter MLST variants, and two blaKPC variants (blaKPC-2, blaKPC-3) associated with seven different Klebsiella MLST variants.

CONCLUSIONS

Molecular surveillance is increasingly becoming a powerful tool to understand bacterial spread in both community and clinical settings. This study provides evidence that genetically related infections in clinical settings do not necessarily reflect temporal associations, and vice versa. Our results also highlight the regional genomic and resistance diversity within related bacterial lineages.

Advancing tuberculosis management: the role of predictive, preventive, and personalized medicine

Tuberculosis is a major global health issue, with approximately 10 million people falling ill and 1.4 million dying yearly. One of the most significant challenges to public health is the emergence of drug-resistant tuberculosis. For the last half-century, treating tuberculosis has adhered to a uniform management strategy in most patients. However, treatment ineffectiveness in some individuals with pulmonary tuberculosis presents a major challenge to the global tuberculosis control initiative. Unfavorable outcomes of tuberculosis treatment (including mortality, treatment failure, loss of follow-up, and unevaluated cases) may result in increased transmission of tuberculosis and the emergence of drug-resistant strains. Treatment failure may occur due to drug-resistant strains, non-adherence to medication, inadequate absorption of drugs, or low-quality healthcare. Identifying the underlying cause and adjusting the treatment accordingly to address treatment failure is important. This is where approaches such as artificial intelligence, genetic screening, and whole genome sequencing can play a critical role. In this review, we suggest a set of particular clinical applications of these approaches, which might have the potential to influence decisions regarding the clinical management of tuberculosis patients.

Whole genome sequencing for tuberculosis in Victoria, Australia: A genomic implementation study from 2017 to 2020

BACKGROUND

Whole genome sequencing (WGS) is increasingly used by tuberculosis (TB) programs to monitor Mycobacterium tuberculosis (Mtb) transmission. We aimed to characterise the molecular epidemiology of TB and Mtb transmission in the low-incidence setting of Victoria, Australia, and assess the utility of WGS.

METHODS

WGS was performed on all first Mtb isolates from TB cases from 2017 to 2020. Potential clusters (≤12 single nucleotide polymorphisms [SNPs]) were investigated for epidemiological links. Transmission events in highly-related (≤5 SNPs) clusters were classified as likely or possible, based on the presence or absence of an epidemiological link, respectively. Case characteristics and transmission settings (as defined by case relationship) were summarised. Poisson regression was used to examine associations with secondary case number.

FINDINGS

Of 1844 TB cases, 1276 (69.2%) had sequenced isolates, with 182 (14.2%) in 54 highly-related clusters, 2-40 cases in size. Following investigation, 140 cases (11.0% of sequenced) were classified as resulting from likely/possible local-transmission, including 82 (6.4%) for which transmission was likely. Common identified transmission settings were social/religious (26.4%), household (22.9%) and family living in different households (7.1%), but many were uncertain (41.4%). While household transmission featured in many clusters (n = 24), clusters were generally smaller (median = 3 cases) than the fewer that included transmission in social/religious settings (n = 12, median = 7.5 cases). Sputum-smear-positivity was associated with higher secondary case numbers.

INTERPRETATION

WGS results suggest Mtb transmission commonly occurs outside the household in our low-incidence setting. Further work is required to optimise the use of WGS in public health management of TB.

FUNDING

The Victorian Tuberculosis Program receives block funding for activities including case management and contact tracing from the Victorian Department of Health. No specific funding for this report was received by manuscript authors or the Victorian Tuberculosis Program, and the funders had no role in the study design, data collection, data analysis, interpretation or report writing.

Predominance of the Mycobacterium tuberculosis Beijing strain amongst children from a high tuberculosis burden township in South Africa

The objective was to determine the molecular epidemiology and drug susceptibility patterns of Mycobacterium tuberculosis (MTB) of children and their household contacts (HHC) in Umlazi, a high TB-burden township in South Africa. Sixty eight MTBRifPLUS positive TB-infected children (TIC) (≤14 years) and 111 HHC were enrolled. Drug susceptibility testing (DST) was performed on sputum samples using the proportion method and GenoType® MTBDR. Genotyping of MTB was conducted using IS6110-restriction fragment length polymorphism (RFLP) and spoligotyping. Rifampicin (RIF) susceptibility was observed in 67/68 TIC. GenoType® MTBDRplus and phenotypic DST identified drug resistant strains in five of 16 culture-confirmed TIC. The Beijing strain was identified in six and the F15/LAM4/KZN strain in one of the 13 TIC respectively. Four patients with unknown RFLP strains belonged to spoligoclades S, T1, T3 variant and X2. The S-lineage and an unknown strain were identified in two HHC. MDR-TB and pre-XDR-TB were identified in one HHC each. Household transmission could not be determined as none of the culture-confirmed TIC resided with the six culture-confirmed contacts. The predominance of the hypervirulent Beijing strain and presence of drug-resistant strains must be considered in the implementation of effective TB control strategies and development of efficacious vaccines.

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.

Characterizing tuberculosis transmission dynamics in high-burden urban and rural settings

Mycobacterium tuberculosis transmission dynamics in high-burden settings are poorly understood. Growing evidence suggests transmission may be characterized by extensive individual heterogeneity in secondary cases (i.e., superspreading), yet the degree and influence of such heterogeneity is largely unknown and unmeasured in high burden-settings. We conducted a prospective, population-based molecular epidemiology study of TB transmission in both an urban and rural setting of Botswana, one of the highest TB burden countries in the world. We used these empirical data to fit two mathematical models (urban and rural) that jointly quantified both the effective reproductive number, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R$$\end{document}R, and the propensity for superspreading in each population. We found both urban and rural populations were characterized by a high degree of individual heterogeneity, however such heterogeneity disproportionately impacted the rural population: 99% of secondary transmission was attributed to only 19% of infectious cases in the rural population compared to 60% in the urban population and the median number of incident cases until the first outbreak of 30 cases was only 32 for the rural model compared to 791 in the urban model. These findings suggest individual heterogeneity plays a critical role shaping local TB epidemiology within subpopulations.

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.

Whole-Genome Sequencing Reveals Recent Transmission of Multidrug-Resistant Mycobacterium tuberculosis CAS1-Kili Strains in Lusaka, Zambia

Globally, tuberculosis (TB) is a major cause of death due to antimicrobial resistance. Mycobacterium tuberculosis CAS1-Kili strains that belong to lineage 3 (Central Asian Strain, CAS) were previously implicated in the spread of multidrug-resistant (MDR)-TB in Lusaka, Zambia. Thus, we investigated recent transmission of those strains by whole-genome sequencing (WGS) with Illumina MiSeq platform. Twelve MDR CAS1-Kili isolates clustered by traditional methods (MIRU-VNTR and spoligotyping) were used. A total of 92% (11/12) of isolates belonged to a cluster (≤12 SNPs) while 50% (6/12) were involved in recent transmission events, as they differed by ≤5 SNPs. All the isolates had KatG Ser315Thr (isoniazid resistance), EmbB Met306 substitutions (ethambutol resistance) and several kinds of rpoB mutations (rifampicin resistance). WGS also revealed compensatory mutations including a novel deletion in embA regulatory region (−35A > del). Several strains shared the same combinations of drug-resistance-associated mutations indicating transmission of MDR strains. Zambian strains belonged to the same clade as Tanzanian, Malawian and European strains, although most of those were pan-drug-susceptible. Hence, complimentary use of WGS to traditional epidemiological methods provides an in-depth insight on transmission and drug resistance patterns which can guide targeted control measures to stop the spread of MDR-TB.

Reporting practices for genomic epidemiology of tuberculosis: a systematic review of the literature using STROME-ID guidelines as a benchmark

BACKGROUND

Pathogen genomics have become increasingly important in infectious disease epidemiology and public health. The Strengthening the Reporting of Molecular Epidemiology for Infectious Diseases (STROME-ID) guidelines were developed to outline a minimum set of criteria that should be reported in genomic epidemiology studies to facilitate assessment of study quality. We evaluate such reporting practices, using tuberculosis as an example.

METHODS

For this systematic review, we initially searched MEDLINE, Embase Classic, and Embase on May 3, 2017, using the search terms "tuberculosis" and "genom* sequencing". We updated this initial search on April 23, 2019, and also included a search of bioRxiv at this time. We included studies in English, French, or Spanish that recruited patients with microbiologically confirmed tuberculosis and used whole genome sequencing for typing of strains. Non-human studies, conference abstracts, and literature reviews were excluded. For each included study, the number and proportion of fulfilled STROME-ID criteria were recorded by two reviewers. A comparison of the mean proportion of fulfilled STROME-ID criteria before and after publication of the STROME-ID guidelines (in 2014) was done using a two-tailed t test. Quasi-Poisson regression and tobit regression were used to examine associations between study characteristics and the number and proportion of fulfilled STROME-ID criteria. This study was registered with PROSPERO, CRD42017064395.

FINDINGS

976 titles and abstracts were identified by our primary search, with an additional 16 studies identified in bioRxiv. 114 full texts (published between 2009 and 2019) were eligible for inclusion. The mean proportion of STROME-ID criteria fulfilled was 50% (SD 12; range 16-75). The proportion of criteria fulfilled was similar before and after STROME-ID publication (51% [SD 11] vs 46% [14], p=0·26). The number of criteria reported (among those applicable to all studies) was not associated with impact factor, h-index, country of affiliation of senior author, or sample size of isolates. Similarly, the proportion of criteria fulfilled was not associated with these characteristics, with the exception of a sample size of isolates of 277 or more (the highest quartile). In terms of reproducibility, 100 (88%) studies reported which bioinformatic tools were used, but only 33 (33%) reported corresponding version numbers. Sequencing data were available for 86 (75%) studies.

INTERPRETATION

The reporting of STROME-ID criteria in genomic epidemiology studies of tuberculosis between 2009 and 2019 was low, with implications for assessment of study quality. The considerable proportion of studies without bioinformatics version numbers or sequencing data available highlights a key concern for reproducibility.

Probable long-term prevalence for a predominant Mycobacterium tuberculosis clone of a Beijing genotype in Colon, Panama

Beijing genotype Mycobacterium tuberculosis strains associate with increased virulence, resistance and/or higher transmission rates. This study describes a specific Beijing strain predominantly identified in the Panamanian province of Colon with one of the highest incidences of tuberculosis in the country. Retrospective mycobacterial interspersed repetitive unit/variable number of tandem repeats analysis of 42 isolates collected between January and August 2018 allowed to identify a cluster (Beijing A) with 17 (40.5%) Beijing isolates. Subsequent prospective strain-specific PCR-based surveillance from September 2019 to March 2020 confirmed the predominance of the Beijing A strain (44.1%) in this province. Whole-genome sequencing revealed higher-than-expected diversity within the cluster, suggesting long-term prevalence of this strain and low number of cases caused by recent transmission. The Beijing A strain belongs to the Asian African 3 (Bmyc13, L2.2.5) branch of the modern Beijing sublineage, with their closest isolates corresponding to cases from Vietnam, probably introduced in Panama between 2000 and 2012.

International tuberculosis contact-tracing notifications in Germany: analysis of national data from 2010 to 2018 and implications for efficiency

BACKGROUND

International contact-tracing (CT) following exposure during long-distance air travel is resource-intensive, whereas evidence for risk of tuberculosis (TB) transmission during international travel is weak. In this study, we systematically analyzed the information from international requests for CT received at the national level in Germany in order to evaluate the continued utility of the current approach and to identify areas for improvement.

METHODS

An anonymized archive of international CT notifications received by the Robert Koch Institute between 2010 and 2018 was searched for key parameters for data collection. A total of 31 parameters, such as characteristics of TB patients and their identified contacts, were extracted from each CT notification and collated into a dataset. Descriptive data analysis and trend analyses were performed to identify key characteristics of CT notifications, patients, and contacts over the years.

RESULTS

192 CT notifications, each corresponding to a single TB index case, were included in the study, increasing from 12 in 2010 to 41 in 2018. The majority of notifications (N = 130, 67.7%) concerned international air travel, followed by private contact (N = 39, 20.3%) and work exposure (N = 16, 8.3%). 159 (82.8%) patients had sputum smear results available, of which 147 (92.5%) were positive. Of 119 (62.0%) patients with drug susceptibility testing results, most (N = 92, 77.3%) had pan-sensitive TB, followed by 15 (12.6%) with multi-drug resistant TB. 115 (59.9%) patients had information on infectiousness, of whom 99 (86.1%) were considered infectious during the exposure period. 7 (5.3%) patients travelled on long-distance flights despite a prior diagnosis of active TB. Of the 771 contact persons, 34 (4.4%) could not be reached for CT measures due to lack of contact information.

CONCLUSION

The high variability in completeness of information contained within the international CT requests emphasizes the need for international standards for reporting of CT information. With the large proportion of TB patients reported to have travelled while being infectious in our study, we feel that raising awareness among patients and health professionals to detect TB early and prevent international long-distance travel during the infectious disease phase should be a cornerstone strategy to safeguard against possible transmission during international travel.

PathOGiST: A Novel Method for Clustering Pathogen Isolates by Combining Multiple Genotyping Signals

In this paper we study the problem of clustering bacterial isolates into epidemiologically related groups from next-generation sequencing data. Existing methods for this problem mainly use a single genotyping signal, and either use a distance-based method with a pre-specified number of clusters, or a phylogenetic tree-based method with a pre-specified threshold. We propose PathOGiST, an algorithmic framework for clustering bacterial isolates by leveraging multiple genotypic signals and calibrated thresholds. PathOGiST uses different genotypic signals, clusters the isolates based on these individual signals with correlation clustering, and combines the clusterings based on the individual signals through consensus clustering. We implemented and tested PathOGiST on three different bacterial pathogens - Escherichia coli, Yersinia pseudotuberculosis, and Mycobacterium tuberculosis - and we conclude by discussing further avenues to explore.