Assessment of Mycobacterium tuberculosis transmission in Oxfordshire, UK, 2007-12, with whole pathogen genome sequences: an observational study

Update on the diagnosis of tuberculosis

BACKGROUND

Tuberculosis (TB) remains a global public health threat, and the development of rapid and precise diagnostic tools is the key to enabling the early start of treatment, monitoring response to treatment, and preventing the spread of the disease.

OBJECTIVES

An overview of recent progress in host- and pathogen-based TB diagnostics.

SOURCES

We conducted a PubMed search of recent relevant articles and guidelines on TB screening and diagnosis.

CONTENT

An overview of currently used methods and perspectives in the following areas of TB diagnostics is provided: immune-based diagnostics, X-ray, clinical symptoms and scores, cough detection, culture of Mycobacterium tuberculosis and identifying its resistance profile using phenotypic and genotypic methods, including next-generation sequencing, sputum- and non-sputum-based molecular diagnosis of TB and monitoring of response to treatment.

IMPLICATIONS

A brief overview of the most relevant advances and changes in international guidelines regarding screening and diagnosing TB is provided in this review. It aims at reviewing all relevant areas of diagnostics, including both pathogen- and host-based methods.

Drug resistance and epidemiological success of modern Mycobacterium tuberculosis lineages in western India

BACKGROUND

Transmission is contributing to the slow decline of tuberculosis (TB) incidence globally. Drivers of TB transmission in India, the country estimated to carry a quarter of the World's burden, are not well studied. We conducted a genomic epidemiology study to compare epidemiological success, host factors and drug resistance (DR) among the four major Mycobacterium tuberculosis (Mtb) lineages (L1-4) circulating in Pune, India.

METHODS

We performed whole-genome sequencing (WGS) of Mtb sputum culture-positive isolates from participants in two prospective cohort studies and predicted genotypic susceptibility using a validated random forest model. We used maximum likelihood estimation to build phylogenies. We compared lineage specific phylogenetic and time-scaled metrics to assess epidemiological success.

RESULTS

Of the 642 isolates that underwent WGS, 612 met sequence quality criteria. Most isolates belonged to L3 (44.6%). The majority (61.1%) of multidrug-resistant isolates belonged to L2 (P < 0.001). In molecular dating, L2 demonstrated a higher rate and more recent resistance acquisition. We measured higher clustering, and time-scaled haplotypic density (THD) for L4 and L2 compared to L3 and/or L1 suggesting higher epidemiological success. L4 demonstrated higher THD and clustering (OR 5.1 (95% CI 2.3-12.3) in multivariate models controlling for host factors and DR.

CONCLUSION

L2 shows a higher frequency of DR and both L2 and L4 demonstrate evidence of higher epidemiological success than L3 or L1 in the study setting. Our findings highlight the need for contact tracing around TB cases, and heightened surveillance of TB DR in India.

Genetic Characterization and Population Structure of Drug-Resistant Mycobacterium tuberculosis Isolated from Brazilian Patients Using Whole-Genome Sequencing

The present study aimed to determine the genetic diversity of isolates of Mycobacterium tuberculosis (Mtb) from presumed drug-resistant tuberculosis patients from several states of Brazil. The isolates had been submitted to conventional drug susceptibility testing for first- and second-line drugs. Multidrug-resistant (MDR-TB) (54.8%) was the most frequent phenotypic resistance profile, in addition to an important high frequency of pre-extensive resistance (p-XDR-TB) (9.2%). Using whole-genome sequencing (WGS), we characterized 298 Mtb isolates from Brazil. Besides the analysis of genotype distribution and possible correlations between molecular and clinical data, we determined the performance of an in-house WGS pipeline with other online pipelines for Mtb lineages and drug resistance profile definitions. Sub-lineage 4.3 (52%) was the most frequent genotype, and the genomic approach revealed a p-XDR-TB level of 22.5%. We detected twenty novel mutations in three resistance genes, and six of these were observed in eight phenotypically resistant isolates. A cluster analysis of 170 isolates showed that 43.5% of the TB patients belonged to 24 genomic clusters, suggesting considerable ongoing transmission of DR-TB, including two interstate transmissions. The in-house WGS pipeline showed the best overall performance in drug resistance prediction, presenting the best accuracy values for five of the nine drugs tested. Significant associations were observed between suffering from fatal disease and genotypic p-XDR-TB (p = 0.03) and either phenotypic (p = 0.006) or genotypic (p = 0.0007) ethambutol resistance. The use of WGS analysis improved our understanding of the population structure of MTBC in Brazil and the genetic and clinical data correlations and demonstrated its utility for surveillance efforts regarding the spread of DR-TB, hopefully helping to avoid the emergence of even more resistant strains and to reduce TB incidence and mortality rates.

Prison as a driver of recent transmissions of multidrug-resistant tuberculosis in Callao, Peru: a cross-sectional study

Background

We sought to identify resistance patterns and key drivers of recent multidrug-resistant tuberculosis (MDR-TB) transmission in a TB-prevalent area in Peru.

Methods

Cross-sectional study including MDR Mycobacterium tuberculosis complex (Mtbc) strains identified in Callao-Peru between April 2017 and February 2019. Mtbc DNA was extracted for whole genome sequencing which was used for phylogenetic inference, clustering, and resistance mutation analyses. Clusters indicative of recent transmission were defined based on a strain-to-strain distance of ≤5 (D5) single nucleotide polymorphisms (SNPs). Epidemiologic factors linked to MDR-TB clustering were analyzed using Poisson regression.

Findings

171 unique MDR-Mtbc strains were included; 22 (13%) had additional fluoroquinolone resistance and were classified as pre-XDR. Six strains (3.5%) harboured bedaquiline (BDQ) resistance mutations and were classified as MDR + BDQ. 158 (92%) Mtbc strains belonged to lineage 4 and 13 (8%) to lineage 2. Using a cluster threshold of ≤5 SNPs, 98 (57%) strains were grouped in one of the 17 D5 clusters indicative of recent transmission, ranging in size from 2 to the largest cluster formed by 53 4.3.3 strains (group_1). Lineage 4.3.3 strains showed the overall highest cluster rate (43%). In multivariate analyses, current or previous imprisonment was independently associated with being part of any MDR-TB transmission clusters (adjusted prevalence ratio [aPR], 1.45; 95% CI, 1.09-1.92).

Interpretation

Pre-XDR-TB emerged in more than 10% of the MDR-TB strains investigated. Transmission of 4.3.3 Mtbc strains especially of the dominant group_1 clone is a major driver of the MDR-TB epidemic in Callao. Current or previous imprisonment was linked to recent MDR-TB transmissions, indicating an important role of prisons in driving the MDR-TB epidemic.

Funding

This work was supported in part by the ERANet-LAC Network of the European Union, Latin America and the Caribbean Countries on Joint Innovation and Research Activities, and FONDECYT. Additional support was received from Leibniz Science Campus Evolutionary Medicine of the Lung, the Deutsche Forschungsgemeinschaft (German Research Foundation, under Germany's Excellence Strategy-EXC 2167 Precision Medicine in Inflammation), and the Research Training Group 2501 TransEvo.

Polyketide synthases mutation in tuberculosis transmission revealed by whole genomic sequence, China, 2011-2019

Introduction: Tuberculosis (TB) is an infectious disease caused by a bacterium called Mycobacterium tuberculosis (Mtb). Previous studies have primarily focused on the transmissibility of multidrug-resistant (MDR) or extensively drug-resistant (XDR) Mtb. However, variations in virulence across Mtb lineages may also account for differences in transmissibility. In Mtb, polyketide synthase (PKS) genes encode large multifunctional proteins which have been shown to be major mycobacterial virulence factors. Therefore, this study aimed to identify the role of PKS mutations in TB transmission and assess its risk and characteristics. Methods: Whole genome sequences (WGSs) data from 3,204 Mtb isolates was collected from 2011 to 2019 in China. Whole genome single nucleotide polymorphism (SNP) profiles were used for phylogenetic tree analysis. Putative transmission clusters (≤10 SNPs) were identified. To identify the role of PKS mutations in TB transmission, we compared SNPs in the PKS gene region between "clustered isolates" and "non-clustered isolates" in different lineages. Results: Cluster-associated mutations in ppsA, pks12, and pks13 were identified among different lineage isolates. They were statistically significant among clustered strains, indicating that they may enhance the transmissibility of Mtb. Conclusion: Overall, this study provides new insights into the function of PKS and its localization in M. tuberculosis. The study found that ppsA, pks12, and pks13 may contribute to disease progression and higher transmission of certain strains. We also discussed the prospective use of mutant ppsA, pks12, and pks13 genes as drug targets.

Pangenome databases improve host removal and mycobacteria classification from clinical metagenomic data

BACKGROUND

Culture-free real-time sequencing of clinical metagenomic samples promises both rapid pathogen detection and antimicrobial resistance profiling. However, this approach introduces the risk of patient DNA leakage. To mitigate this risk, we need near-comprehensive removal of human DNA sequences at the point of sequencing, typically involving the use of resource-constrained devices. Existing benchmarks have largely focused on the use of standardized databases and largely ignored the computational requirements of depletion pipelines as well as the impact of human genome diversity.

RESULTS

We benchmarked host removal pipelines on simulated and artificial real Illumina and Nanopore metagenomic samples. We found that construction of a custom kraken database containing diverse human genomes results in the best balance of accuracy and computational resource usage. In addition, we benchmarked pipelines using kraken and minimap2 for taxonomic classification of Mycobacterium reads using standard and custom databases. With a database representative of the Mycobacterium genus, both tools obtained improved specificity and sensitivity, compared to the standard databases for classification of Mycobacterium tuberculosis. Computational efficiency of these custom databases was superior to most standard approaches, allowing them to be executed on a laptop device.

CONCLUSIONS

Customized pangenome databases provide the best balance of accuracy and computational efficiency when compared to standard databases for the task of human read removal and M. tuberculosis read classification from metagenomic samples. Such databases allow for execution on a laptop, without sacrificing accuracy, an especially important consideration in low-resource settings. We make all customized databases and pipelines freely available.

Signatures of transmission in within-host M. tuberculosis variation

Background

Because M. tuberculosis evolves slowly, transmission clusters often contain multiple individuals with identical consensus genomes, making it difficult to reconstruct transmission chains. Finding additional sources of shared M. tuberculosis variation could help overcome this problem. Previous studies have reported M. tuberculosis diversity within infected individuals; however, whether within-host variation improves transmission inferences remains unclear.

Methods

To evaluate the transmission information present in within-host M. tuberculosis variation, we re-analyzed publicly available sequence data from three household transmission studies, using household membership as a proxy for transmission linkage between donor-recipient pairs.

Findings

We found moderate levels of minority variation present in M. tuberculosis sequence data from cultured isolates that varied significantly across studies (mean: 6, 7, and 170 minority variants above a 1% minor allele frequency threshold, outside of PE/PPE genes). Isolates from household members shared more minority variants than did isolates from unlinked individuals in the three studies (mean 98 shared minority variants vs. 10; 0.8 vs. 0.2, and 0.7 vs. 0.2, respectively). Shared within-host variation was significantly associated with household membership (OR: 1.51 [1.30,1.71], for one standard deviation increase in shared minority variants). Models that included shared within-host variation improved the accuracy of predicting household membership in all three studies as compared to models without within-host variation (AUC: 0.95 versus 0.92, 0.99 versus 0.95, and 0.93 versus 0.91).

Interpretation

Within-host M. tuberculosis variation persists through culture and could enhance the resolution of transmission inferences. The substantial differences in minority variation recovered across studies highlights the need to optimize approaches to recover and incorporate within-host variation into automated phylogenetic and transmission inference.

Funding

NIAID: 5K01AI173385.

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.

Implications for TB control among migrants in large cities in China: A prospective population-based genomic epidemiology study in Shenzhen

Internal migrants are a challenge for TB control in large Chinese cities and understanding this epidemiology is crucial for designing effective control and prevention strategies. We conducted a prospective genomic epidemiological study of culture-positive TB patients diagnosed between June 1, 2018 and May 31, 2021 in the Longhua District of Shenzhen. Treatment status was obtained from local and national TB registries and all isolates were sequenced. Genomic clusters were defined as strains differing by ≤12 SNPs. Risk factors for clustering were identified with multivariable analysis and then Bayesian models and TransPhylo were used to infer the timing of transmission within clusters. Of the 2277 culture-positive patients, 70.1% (1596/2277) were migrants: 72.1% (1043/1446) of the migrants patients developed TB within two years of arriving in Longhua; 38.8% within 6 months of arriving; and 12.3% (104/843) had TB symptoms when they arrived. Only 15.4% of Longhua strains were in genomic clusters. More than one third (33.6%) of patients were not treated in Shenzhen but were involved in nearly one third of the recent transmission events. Clustering was associated with migrants not treated in Shenzhen, males, and teachers/trainers. TB in Longhua is prinicipally due to reactivation of infections in migrants, but a proportion may have had clinical or incipient TB upon arrival in the district. Patients diagnosed but not treated in Longhua were involved in recent local TB transmission. Controlling TB in Shenzhen will require strategies to comprehensively diagnose and treat active TB in the internal migrant population.

Genotypic and spatial analysis of transmission dynamics of tuberculosis in Shanghai, China: a 10-year prospective population-based surveillance study

Background

With improved tuberculosis (TB) control programs, the incidence of TB in China declined dramatically over the past few decades, but recently the rate of decrease has slowed, especially in large cities such as Shanghai. To help formulate strategies to further reduce TB incidence, we performed a 10-year study in Songjiang, a district of Shanghai, to delineate the characteristics, transmission patterns, and dynamic changes of the local TB burden.

Methods

We conducted a population-based study of culture-positive pulmonary TB patients diagnosed in Songjiang during 2011-2020. Genomic clusters were defined with a threshold distance of 12-single-nucleotide-polymorphisms based on whole-genome sequencing, and risk factors for clustering were identified by logistic regression. Transmission inference was performed using phybreak. The distances between the residences of patients were compared to the genomic distances of their isolates. Spatial patient hotspots were defined with kernel density estimation.

Findings

Of 2212 enrolled patients, 74.7% (1652/2212) were internal migrants. The clustering rate (25.2%, 558/2212) and spatial concentrations of clustered and unclustered patients were unchanged over the study period. Migrants had significantly higher TB rates but less clustering than residents. Clustering was highest in male migrants, younger patients and both residents and migrants employed in physical labor. Only 22.1% of transmission events occurred between residents and migrants, with residents more likely to transmit to migrants. The clustering risk decreased rapidly with increasing distances between patient residences, but more than half of clustered patient pairs lived ≥5 km apart. Epidemiologic links were identified for only 15.6% of clustered patients, mostly in close contacts.

Interpretation

Although some of the TB in Songjiang's migrant population is caused by strains brought by infected migrants, local, recent transmission is an important driver of the TB burden. These results suggest that further reductions in TB incidence require novel strategies to detect TB early and interrupt urban transmission.

Funding

Shanghai Municipal Science and Technology Major Project (ZD2021CY001), National Natural Science Foundation of China (82272376), National Research Council of Science and Technology Major Project of China (2017ZX10201302-006).

Ending tuberculosis: the cost of missing the World Health Organization target in a low-incidence country

BACKGROUND

Ending tuberculosis (TB) is a global priority and targets for doing so are outlined in the World Health Organization (WHO) End TB Strategy. For low-incidence countries, eliminating TB requires high levels of wealth, low levels of income inequality and effective TB programmes and services that can meet the needs of people who have not benefited from these and are still at risk of TB. In Ireland, numerous reports have noted a need for more funding for TB prevention and control.

AIM

The aim of this research was to estimate the cost of not meeting the WHO End TB target of a 90% reduction in TB incidence in Ireland between 2015 and 2035.

METHODS

The cost of projected TB cases between 2022 and 2035 is estimated based on trends in surveillance data for the period 2015 to 2019 and outcomes reported in the literature.

RESULTS

Between 2022 and 2035, it is projected that a failure to meet the WHO End TB Strategy target will result in an additional 989 cases of TB, 577.3 disability-adjusted life years and 35 deaths with TB in Ireland. The cost of this is estimated to be €70.779 million.

CONCLUSION

Given the estimated cost, Ireland's current prospects of eliminating TB and the tendency towards programmatic funding internationally, greater investment in TB prevention and control in Ireland is justifiable. A national elimination strategy with actions at the levels of the social determinants of health, the health system and the TB programme should be funded.

Trends in pulmonary tuberculosis mortality between 1985 and 2018: an observational analysis

BACKGROUND

Pulmonary tuberculosis (TB) is a major source of global morbidity and mortality. Latent infection has enabled it to spread to a quarter of the world's population. The late 1980s and early 1990s saw an increase in the number of TB cases related to the HIV epidemic, and the spread of multidrug-resistant TB. Few studies have reported pulmonary TB mortality trends. Our study reports and compares trends in pulmonary TB mortality.

METHODS

We utilized the World Health Organization (WHO) mortality database from 1985 through 2018 to analyze TB mortality using the International Classification of Diseases-10 codes. Based on the availability and quality of data, we investigated 33 countries including two countries from the Americas; 28 countries from Europe; and 3 countries from the Western Pacific region. Mortality rates were dichotomized by sex. We computed age-standardized death rates per 100,000 population using the world standard population. Time trends were investigated using joinpoint regression analysis.

RESULTS

We observed a uniform decrease in mortality in all countries across the study period except the Republic of Moldova, which showed an increase in female mortality (+ 0.12 per 100,000 population). Among all countries, Lithuania had the greatest reduction in male mortality (-12) between 1993-2018, and Hungary had the greatest reduction in female mortality (-1.57) between 1985-2017. For males, Slovenia had the most rapid recent declining trend with an estimated annual percentage change (EAPC) of -47% (2003-2016), whereas Croatia showed the fastest increase (EAPC, + 25.0% [2015-2017]). For females, New Zealand had the most rapid declining trend (EAPC, -47.2% [1985-2015]), whereas Croatia showed a rapid increase (EAPC, + 24.9% [2014-2017]).

CONCLUSIONS

Pulmonary TB mortality is disproportionately higher among Central and Eastern European countries. This communicable disease cannot be eliminated from any one region without a global approach. Priority action areas include ensuring early diagnosis and successful treatment to the most vulnerable groups such as people of foreign origin from countries with a high burden of TB and incarcerated population. Incomplete reporting of TB-related epidemiological data to WHO excluded high-burden countries and limited our study to 33 countries only. Improvement in reporting is crucial to accurately identify changes in epidemiology, the effect of new treatments, and management approaches.

Transmission and resistome of extremely drug-resistant tuberculosis in Beijing, China: A retrospective population-based epidemiological study

BACKGROUND

In this study, we utilized whole genome sequencing (WGS) of clinical extremely drug-resistant tuberculosis (EDR-TB) strains collected during 2014-2020 in Beijing to detect clustered strains.

METHODS

A retrospective cohort study was conducted by inclusion of EDR-TB patients with positive cultures in Beijing between 2014 and 2020.

RESULTS

A total of 95 EDR-TB patients were included in our analysis. Up on the WGS based genotyping, 94 (94/95, 98.9%) out of 95 were identified as lineage 2 (East Asia). The pairwise genomic distance analysis identified 7 clusters, ranging in size from 2 to 5 isolates. The clustering rate of EDR-TB was 21.1%; while no patients had significantly higher odds of clustering. All isolates harbor rpoB RRDR mutations that confer RIF resistance and katG or inhA promoter mutations that confer INH resistance. Of 95 EDR-TB isolates, a total of 15 mutation types were recorded in the transcriptional regulator mmpR5. In vitro susceptibility testing results revealed that 14 (14/15, 93.3%) out of 15 mutation types were resistant to CFZ; whereas only 3 (3/15, 20.0%) showed resistance to BDQ. Interestingly, 12 isolates harbored mutations within rrl locus, whereas only mutations at positions 2294 and 2296 conferred CLA resistance. Favorable outcomes of EDR-TB patients were positively associated with more effective drugs in the regimes.

CONCLUSION

WGS data demonstrate limited transmission of EDR-TB in this metropolis city. WGS-based drug susceptibility predictions will bring benefits to EDR-TB patients to formulate optimal therapeutic regimens.

Contact investigation of tuberculosis in Shanghai, China: A 13-year cohort

OBJECTIVES

Tuberculosis (TB) contact investigation is an effective strategy for TB control. We established a close-contact cohort to evaluate the effect of contact investigation in Shanghai, China.

METHODS

Close contacts of bacteria-positive TB cases diagnosed from 2009 to 2018 in the Songjiang District of Shanghai were screened and followed up until 2021 through symptom examination or other alternatives, and the incidence and risk factors for developing active TB among close contacts were assessed by survival analysis and WGS.

RESULTS

7018 close contacts of 2861 bacteria-positive TB cases were investigated. The median follow-up time was 7.4 years. 97 close contacts (185/100,000, 95%CI: 151-226) developed TB. Survival analysis showed that index cases who were in younger age groups, living in urban settings, smear-positive, and reported a healthcare-seeking delay for TB had higher risks of generating contact cases. Close contacts with matched strains developed TB at an average of 26.5 months based on WGS.

CONCLUSION

A 2-year ideally longer follow-up for close contacts would be beneficial. To improve early case detection, more attention needs to be paid to those contacts whose index cases were living in urban settings, were smear-positive, had a healthcare-seeking delay for TB, and/or were aged ≤18 years.

Risk for Prison-to-Community Tuberculosis Transmission, Thailand, 2017-2020

To determine contributions of previously incarcerated persons to tuberculosis (TB) transmission in the community, we performed a healthcare facility-based cohort study of TB patients in Thailand during 2017-2020. We used whole-genome sequencing of Mycobacterium tuberculosis isolates from patients to identify genotypic clusters and assess the association between previous incarceration and TB transmission in the community. We identified 4 large genotype clusters (>10 TB patients/cluster); 28% (14/50) of the patients in those clusters were formerly incarcerated. Formerly incarcerated TB patients were more likely than nonincarcerated patients to be included in large clusters. TB patients within the large genotype clusters were geographically dispersed throughout Chiang Rai Province. Community TB transmission in the community was associated with the presence of formerly incarcerated individuals in Thailand. To reduce the risk for prison-to-community transmission, we recommend TB screening at the time of entry and exit from prisons and follow-up screening in the community.

Advances in the Microbiological Diagnosis of Prosthetic Joint Infections

A significant number of prosthetic joint infections (PJI) are culture-negative and/or misinterpreted as aseptic failures in spite of the correct implementation of diagnostic culture techniques, such as tissue sample processing in a bead mill, prolonged incubation time, or sonication of removed implants. Misinterpretation may lead to unnecessary surgery and needless antimicrobial treatment. The diagnostic value of non-culture techniques has been investigated in synovial fluid, periprosthetic tissues, and sonication fluid. Different feasible improvements, such as real-time technology, automated systems and commercial kits are now available to support microbiologists. In this review, we describe non-culture techniques based on nucleic acid amplification and sequencing methods. Polymerase chain reaction (PCR) is a frequently used technique in most microbiology laboratories which allows the detection of a nucleic acid fragment by sequence amplification. Different PCR types can be used to diagnose PJI, each one requiring the selection of appropriate primers. Henceforward, thanks to the reduced cost of sequencing and the availability of next-generation sequencing (NGS), it will be possible to identify the whole pathogen genome sequence and, additionally, to detect all the pathogen sequences present in the joint. Although these new techniques have proved helpful, strict conditions need to be observed in order to detect fastidious microorganisms and rule out contaminants. Specialized microbiologists should assist clinicians in interpreting the result of the analyses at interdisciplinary meetings. New technologies will gradually be made available to improve the etiologic diagnoses of PJI, which will remain an important cornerstone of treatment. Strong collaboration among all specialists involved is essential for the correct diagnosis of PJI.

Evaluation of Nanopore sequencing for Mycobacterium tuberculosis drug susceptibility testing and outbreak investigation: a genomic analysis

BACKGROUND

Mycobacterium tuberculosis whole-genome sequencing (WGS) has been widely used for genotypic drug susceptibility testing (DST) and outbreak investigation. For both applications, Illumina technology is used by most public health laboratories; however, Nanopore technology developed by Oxford Nanopore Technologies has not been thoroughly evaluated. The aim of this study was to determine whether Nanopore sequencing data can provide equivalent information to Illumina for transmission clustering and genotypic DST for M tuberculosis.

METHODS

In this genomic analysis, we analysed 151 M tuberculosis isolates from Madagascar, South Africa, and England, which were collected between 2011 and 2018, using phenotypic DST and matched Illumina and Nanopore data. Illumina sequencing was done with the MiSeq, HiSeq 2500, or NextSeq500 platforms and Nanopore sequencing was done on the MinION or GridION platforms. Using highly reliable PacBio sequencing assemblies and pairwise distance correlation between Nanopore and Illumina data, we optimise Nanopore variant filters for detecting single-nucleotide polymorphisms (SNPs; using BCFtools software). We then used those SNPs to compare transmission clusters identified by Nanopore with the currently used UK Health Security Agency Illumina pipeline (COMPASS). We compared Illumina and Nanopore WGS-based DST predictions using the Mykrobe software and mutation catalogue.

FINDINGS

The Nanopore BCFtools pipeline identified SNPs with a median precision of 99·3% (IQR 99·1-99·6) and recall of 90·2% (88·1-94·2) compared with a precision of 99·6% (99·4-99·7) and recall of 91·9% (87·6-98·6) using the Illumina COMPASS pipeline. Using a threshold of 12 SNPs for putative transmission clusters, Illumina identified 98 isolates as unrelated and 53 as belonging to 19 distinct clusters (size range 2-7). Nanopore reproduced 15 out of 19 clusters perfectly; two clusters were merged into one cluster, one cluster had a single sample missing, and one cluster had an additional sample adjoined. Illumina-based clusters were also closely replicated using a five SNP threshold and clustering accuracy was maintained using mixed Illumina and Nanopore datasets. Genotyping resistance variants with Nanopore was highly concordant with Illumina, having zero discordant SNPs across more than 3000 SNPs and four insertions or deletions (indels), across 60 000 indels.

INTERPRETATION

Illumina and Nanopore technologies can be used independently or together by public health laboratories performing M tuberculosis genotypic DST and outbreak investigations. As a result, clinical and public health institutions making decisions on which sequencing technology to adopt for tuberculosis can base the choice on cost (which varies by country), batching, and turnaround time.

FUNDING

Academy for Medical Sciences, Oxford Wellcome Institutional Strategic Support Fund, and the Swiss South Africa Joint Research Award (Swiss National Science Foundation and South African National Research Foundation).

Distribution and origins of Mycobacterium tuberculosis L4 in Southeast Asia

Molecular and genomic studies have revealed that Mycobacterium tuberculosis Lineage 4 (L4, Euro-American lineage) emerged in Europe before becoming distributed around the globe by trade routes, colonial migration and other historical connections. Although L4 accounts for tens or hundreds of thousands of tuberculosis (TB) cases in multiple Southeast Asian countries, phylogeographical studies have either focused on a single country or just included Southeast Asia as part of a global analysis. Therefore, we interrogated public genomic data to investigate the historical patterns underlying the distribution of L4 in Southeast Asia and surrounding countries. We downloaded 6037 genomes associated with 29 published studies, focusing on global analyses of L4 and Asian studies of M. tuberculosis. We identified 2256 L4 genomes including 968 from Asia. We show that 81 % of L4 in Thailand, 51 % of L4 in Vietnam and 9 % of L4 in Indonesia belong to sub-lineages of L4 that are rarely seen outside East and Southeast Asia (L4.2.2, L4.4.2 and L4.5). These sub-lineages have spread between East and Southeast Asian countries, with no recent European ancestor. Although there is considerable uncertainty about the exact direction and order of intra-Asian M. tuberculosis dispersal, due to differing sampling frames between countries, our analysis suggests that China may be the intermediate location between Europe and Southeast Asia for two of the three predominantly East and Southeast Asian L4 sub-lineages (L4.2.2 and L4.5). This new perspective on L4 in Southeast Asia raises the possibility of investigating host population-specific evolution and highlights the need for more structured sampling from Southeast Asian countries to provide more certainty of the historical and current routes of dispersal.

Genome-wide host-pathogen analyses reveal genetic interaction points in tuberculosis disease

The genetics underlying tuberculosis (TB) pathophysiology are poorly understood. Human genome-wide association studies have failed so far to reveal reproducible susceptibility loci, attributed in part to the influence of the underlying Mycobacterium tuberculosis (Mtb) bacterial genotype on the outcome of the infection. Several studies have found associations of human genetic polymorphisms with Mtb phylo-lineages, but studies analysing genome-genome interactions are needed. By implementing a phylogenetic tree-based Mtb-to-human analysis for 714 TB patients from Thailand, we identify eight putative genetic interaction points (P < 5 × 10-8) including human loci DAP and RIMS3, both linked to the IFNγ cytokine and host immune system, as well as FSTL5, previously associated with susceptibility to TB. Many of the corresponding Mtb markers are lineage specific. The genome-to-genome analysis reveals a complex interactome picture, supports host-pathogen adaptation and co-evolution in TB, and has potential applications to large-scale studies across many TB endemic populations matched for host-pathogen genomic diversity.

Phylogeography and transmission of Mycobacterium tuberculosis spanning prisons and surrounding communities in Paraguay

Recent rises in incident tuberculosis (TB) cases in Paraguay and the increasing concentration of TB within prisons highlight the urgency of targeting strategies to interrupt transmission and prevent new infections. However, whether specific cities or carceral institutions play a disproportionate role in transmission remains unknown. We conducted prospective genomic surveillance, sequencing 471 Mycobacterium tuberculosis complex genomes, from inside and outside prisons in Paraguay's two largest urban areas, Asunción and Ciudad del Este, from 2016 to 2021. We found genomic evidence of frequent recent transmission within prisons and transmission linkages spanning prisons and surrounding populations. We identified a signal of frequent M. tuberculosis spread between urban areas and marked recent population size expansion of the three largest genomic transmission clusters. Together, our findings highlight the urgency of strengthening TB control programs to reduce transmission risk within prisons in Paraguay, where incidence was 70 times that outside prisons in 2021.

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.

Analysis of the twenty-six largest outbreaks of tuberculosis in Aragon using whole-genome sequencing for surveillance purposes

The incidence of tuberculosis in Aragon, Spain, is around ten cases per 100,000 inhabitants. Since 2004, a molecular surveillance protocol has been carried out; therefore, all M. tuberculosis strains are genotyped. Recently, whole-genome sequencing has been implemented for relevant isolates. The aim of this work is to characterise at the molecular level the causative strains of the 26 largest outbreaks of the community (including ten or more cases), genotyped by IS6110-RFLP and causing 26% of tuberculosis cases. To achieve this objective, two or three isolates of each IS6110-cluster belonging to different years were selected for sequencing. We found that strains of lineages L4.8, L4.3 and L4.1.2 were the most frequent. The threshold of 12 SNPs as the maximum distance for confirming the belonging to an outbreak was met for 18 of the 26 IS6110-clusters. Four pairs of isolates with more than 90 SNPs were identified as not belonging to the same strain, and four other pairs were kept in doubt as the number of SNPs was close to 12, between 14 and 35. The study of Regions of Difference revealed that they are lineage conserved. Moreover, we could analyse the IS6110 locations for all genome-sequenced isolates, finding some frequent locations in isolates belonging to the same lineage and certain IS6110 movements between the paired isolates. In the vast majority, these movements were not captured by the IS6110-RFLP pattern. After classifying the genes containing SNP by their functional category, we could confirm that the number of SNPs detected in genes considered as virulence factors and the number of cases the strain produced were not related, suggesting that a particular SNP is more relevant than the number. The characteristics found in the most successful strains in our community could be useful for other researchers in epidemiology, virulence and pathogenesis.

Transmission of multidrug-resistant tuberculosis in Beijing, China: An epidemiological and genomic analysis

Background

Understanding multidrug-resistant tuberculosis (MDR-TB) transmission patterns is crucial for controlling the disease. We aimed to identify high-risk populations and geographic settings of MDR-TB transmission.

Methods

We conducted a population-based retrospective study of MDR-TB patients in Beijing from 2018 to 2020, and assessed MDR-TB recent transmission using whole-genome sequencing of isolates. Geospatial analysis was conducted with kernel density estimation. We combined TransPhylo software with epidemiological investigation data to construct transmission networks. Logistic regression analysis was utilized to identify risk factors for recent transmission.

Results

We included 241 MDR-TB patients, of which 146 (60.58%) were available for genomic analysis. Drug resistance prediction showed that resistance to fluoroquinolones (FQs) was as high as 39.74% among new cases. 36 (24.66%) of the 146 MDR strains were grouped into 12 genome clusters, suggesting recent transmission of MDR strains. 44.82% (13/29) of the clustered patients lived in the same residential community, adjacent residential community or the same street as other cases. The inferred transmission chain found a total of 6 transmission events in 3 clusters; of these, 4 transmission events occurred in residential areas and nearby public places. Logistic regression analysis revealed that being aged 25-34 years-old was a risk factor for recent transmission.

Conclusions

The recent transmission of MDR-TB in Beijing is severe, and residential areas are common sites of transmission; high levels of FQs drug resistance suggest that FQs should be used with caution unless resistance can be ruled out by laboratory testing.

Connecting the dots: understanding how human mobility shapes TB epidemics

Tuberculosis (TB) remains a leading infectious cause of death worldwide. Reducing TB infections and TB-related deaths rests ultimately on stopping forward transmission from infectious to susceptible individuals. Critical to this effort is understanding how human host mobility shapes the transmission and dispersal of new or existing strains of Mycobacterium tuberculosis (Mtb). Important questions remain unanswered. What kinds of mobility, over what temporal and spatial scales, facilitate TB transmission? How do human mobility patterns influence the dispersal of novel Mtb strains, including emergent drug-resistant strains? This review summarizes the current state of knowledge on mobility and TB epidemic dynamics, using examples from three topic areas, including inference of genetic and spatial clustering of infections, delineating source-sink dynamics, and mapping the dispersal of novel TB strains, to examine scientific questions and methodological issues within this topic. We also review new data sources for measuring human mobility, including mobile phone-associated movement data, and discuss important limitations on their use in TB epidemiology.

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.

Epidemiological cut-off values for a 96-well broth microdilution plate for high-throughput research antibiotic susceptibility testing of M. tuberculosis

Drug susceptibility testing of M. tuberculosis is rooted in a binary susceptible/resistant paradigm. While there are considerable advantages in measuring the minimum inhibitory concentrations (MICs) of a panel of drugs for an isolate, it is necessary to measure the epidemiological cut-off values (ECOFF/ECVs) to permit comparison with qualitative data. Here we present ECOFF/ECVs for 13 anti-tuberculosis compounds, including bedaquiline and delamanid, derived from 20 637 clinical isolates collected by 14 laboratories based in 11 countries on five continents. Each isolate was incubated for 14 days on a dry 96-well broth microdilution plate and then read. Resistance to most of the drugs due to prior exposure is expected and the MIC distributions for many of the compounds are complex, and therefore a phenotypically wild-type population could not be defined. Since a majority of samples also underwent genetic sequencing, we defined a genotypically wild-type population and measured the MIC of the 99th percentile by direct measurement and via fitting a Gaussian using interval regression. The proposed ECOFF/ECVs were then validated by comparing with the MIC distributions of high-confidence genetic variants that confer resistance and with qualitative drug susceptibility tests obtained via the Mycobacterial Growth Indicator Tube (MGIT) system or Microscopic-Observation Drug Susceptibility (MODS) assay. These ECOFF/ECVs will inform and encourage the more widespread adoption of broth microdilution: this is a cheap culture-based method that tests the susceptibility of 12-14 antibiotics on a single 96-well plate and so could help personalise the treatment of tuberculosis.

Using Genomics to Understand the Epidemiology of Infectious Diseases in the Northern Territory of Australia

The Northern Territory (NT) is a geographically remote region of northern and central Australia. Approximately a third of the population are First Nations Australians, many of whom live in remote regions. Due to the physical environment and climate, and scale of social inequity, the rates of many infectious diseases are the highest nationally. Molecular typing and genomic sequencing in research and public health have provided considerable new knowledge on the epidemiology of infectious diseases in the NT. We review the applications of genomic sequencing technology for molecular typing, identification of transmission clusters, phylogenomics, antimicrobial resistance prediction, and pathogen detection. We provide examples where these methodologies have been applied to infectious diseases in the NT and discuss the next steps in public health implementation of this technology.

A data compendium associating the genomes of 12,289 Mycobacterium tuberculosis isolates with quantitative resistance phenotypes to 13 antibiotics

The Comprehensive Resistance Prediction for Tuberculosis: an International Consortium (CRyPTIC) presents here a data compendium of 12,289 Mycobacterium tuberculosis global clinical isolates, all of which have undergone whole-genome sequencing and have had their minimum inhibitory concentrations to 13 antitubercular drugs measured in a single assay. It is the largest matched phenotypic and genotypic dataset for M. tuberculosis to date. Here, we provide a summary detailing the breadth of data collected, along with a description of how the isolates were selected, collected, and uniformly processed in CRyPTIC partner laboratories across 23 countries. The compendium contains 6,814 isolates resistant to at least 1 drug, including 2,129 samples that fully satisfy the clinical definitions of rifampicin resistant (RR), multidrug resistant (MDR), pre-extensively drug resistant (pre-XDR), or extensively drug resistant (XDR). The data are enriched for rare resistance-associated variants, and the current limits of genotypic prediction of resistance status (sensitive/resistant) are presented by using a genetic mutation catalogue, along with the presence of suspected resistance-conferring mutations for isolates resistant to the newly introduced drugs bedaquiline, clofazimine, delamanid, and linezolid. Finally, a case study of rifampicin monoresistance demonstrates how this compendium could be used to advance our genetic understanding of rare resistance phenotypes. The data compendium is fully open source and it is hoped that it will facilitate and inspire future research for years to come.

Population-based sequencing of Mycobacterium tuberculosis reveals how current population dynamics are shaped by past epidemics

Transmission is a driver of tuberculosis (TB) epidemics in high-burden regions, with assumed negligible impact in low-burden areas. However, we still lack a full characterization of transmission dynamics in settings with similar and different burdens. Genomic epidemiology can greatly help to quantify transmission, but the lack of whole genome sequencing population-based studies has hampered its application. Here, we generate a population-based dataset from Valencia region and compare it with available datasets from different TB-burden settings to reveal transmission dynamics heterogeneity and its public health implications. We sequenced the whole genome of 785 Mycobacterium tuberculosis strains and linked genomes to patient epidemiological data. We use a pairwise distance clustering approach and phylodynamic methods to characterize transmission events over the last 150 years, in different TB-burden regions. Our results underscore significant differences in transmission between low-burden TB settings, i.e., clustering in Valencia region is higher (47.4%) than in Oxfordshire (27%), and similar to a high-burden area as Malawi (49.8%). By modeling times of the transmission links, we observed that settings with high transmission rate are associated with decades of uninterrupted transmission, irrespective of burden. Together, our results reveal that burden and transmission are not necessarily linked due to the role of past epidemics in the ongoing TB incidence, and highlight the need for in-depth characterization of transmission dynamics and specifically tailored TB control strategies.

Fine-grain population structure and transmission patterns of Mycobacterium tuberculosis in southern Mozambique, a high TB/HIV burden area

Genomic studies of the Mycobacterium tuberculosis complex (MTBC) might shed light on the dynamics of its transmission, especially in high-burden settings, where recent outbreaks are embedded in the complex natural history of the disease. To this end, we conducted a 1 year prospective surveillance-based study in Mozambique. We applied whole-genome sequencing (WGS) to 295 positive cultures. We fully characterized MTBC isolates by phylogenetics and dating evaluation, and carried out a molecular epidemiology analysis to investigate further associations with pre-defined transmission risk factors. The majority of strains (49.5%, 136/275) belonged to lineage (L) 4; 57.8 % of them (159/275) were in genomic transmission clusters (cut-off 5 SNPs), and a strikingly high proportion (45.5%) shared an identical genotype (0 SNP pairwise distance). We found two ‘likely endemic’ clades, comprising 67 strains, belonging to L1.2, which dated back to the late 19th century and were associated with recent spread among people living with human immunodeficiency virus (PLHIV). We describe for the first time the population structure of MTBC in our region, a high tuberculosis (TB)/HIV burden area. Clustering analysis revealed an unforeseen pattern of spread and high rates of progression to active TB, suggesting weaknesses in TB control activities. The long-term presence of local strains in Mozambique, which were responsible for large transmission among HIV/TB-coinfected patients, calls into question the role of HIV in TB transmission.

Whole-genome sequencing of Mycobacterium tuberculosis from Cambodia

Cambodia has one of the highest tuberculosis (TB) incidence rates in the WHO Western Pacific region. Remarkably though, the prevalence of multidrug-resistant TB (MDR-TB) remains low. We explored the genetic diversity of Mycobacterium tuberculosis (MTB) circulating in this unique setting using whole-genome sequencing (WGS). From October 2017 until January 2018, we collected one hundred sputum specimens from consenting adults older than 21 years of age, newly diagnosed with bacteriologically confirmed TB in 3 districts of Phnom Penh and Takeo provinces of Cambodia before they commence on their TB treatment, where eighty MTB isolates were successfully cultured and sequenced. Majority of the isolates belonged to Lineage 1 (Indo-Oceanic) (69/80, 86.25%), followed by Lineage 2 (East Asian) (10/80, 12.5%) and Lineage 4 (Euro-American) (1/80, 1.25%). Phenotypic resistance to both streptomycin and isoniazid was found in 3 isolates (3/80, 3.75%), while mono-resistance to streptomycin and isoniazid was identical at 2.5% (N = 2 each). None of the isolates tested was resistant to either rifampicin or ethambutol. The specificities of genotypic prediction for resistance to all drugs tested were 100%, while the sensitivities of genotypic resistance predictions to isoniazid and streptomycin were lower at 40% (2/5) and 80% (4/5) respectively. We identified 8 clusters each comprising of two to five individuals all residing in the Takeo province, making up half (28/56, 50%) of all individuals sampled in the province, indicating the presence of multiple ongoing transmission events. All clustered isolates were of Lineage 1 and none are resistant to any of the drugs tested. This study while demonstrating the relevance and utility of WGS in predicting drug resistance and inference of disease transmission, highlights the need to increase the representation of genotype-phenotype TB data from low and middle income countries in Asia and Africa to improve the accuracies for prediction of drug resistance.

Mutation of Mycobacterium tuberculosis and Implications for Using Whole-Genome Sequencing for Investigating Recent Tuberculosis Transmission

Tuberculosis (TB) control programs use whole-genome sequencing (WGS) of Mycobacterium tuberculosis (Mtb) for detecting and investigating TB case clusters. Existence of few genomic differences between Mtb isolates might indicate TB cases are the result of recent transmission. However, the variable and sometimes long duration of latent infection, combined with uncertainty in the Mtb mutation rate during latency, can complicate interpretation of WGS results. To estimate the association between infection duration and single nucleotide polymorphism (SNP) accumulation in the Mtb genome, we first analyzed pairwise SNP differences among TB cases from Los Angeles County, California, with strong epidemiologic links. We found that SNP distance alone was insufficient for concluding that cases are linked through recent transmission. Second, we describe a well-characterized cluster of TB cases in California to illustrate the role of genomic data in conclusions regarding recent transmission. Longer presumed latent periods were inconsistently associated with larger SNP differences. Our analyses suggest that WGS alone cannot be used to definitively determine that a case is attributable to recent transmission. Methods for integrating clinical, epidemiologic, and genomic data can guide conclusions regarding the likelihood of recent transmission, providing local public health practitioners with better tools for monitoring and investigating TB transmission.

Application of Next Generation Sequencing for Diagnosis and Clinical Management of Drug-Resistant Tuberculosis: Updates on Recent Developments in the Field

The World Health Organization’s End TB Strategy prioritizes universal access to an early diagnosis and comprehensive drug susceptibility testing (DST) for all individuals with tuberculosis (TB) as a key component of integrated, patient-centered TB care. Next generation whole genome sequencing (WGS) and its associated technology has demonstrated exceptional potential for reliable and comprehensive resistance prediction for Mycobacterium tuberculosis isolates, allowing for accurate clinical decisions. This review presents a descriptive analysis of research describing the potential of WGS to accelerate delivery of individualized care, recent advances in sputum-based WGS technology and the role of targeted sequencing for resistance detection. We provide an update on recent research describing the mechanisms of resistance to new and repurposed drugs and the dynamics of mixed infections and its potential implication on TB diagnosis and treatment. Whilst the studies reviewed here have greatly improved our understanding of recent advances in this arena, it highlights significant challenges that remain. The wide-spread introduction of new drugs in the absence of standardized DST has led to rapid emergence of drug resistance. This review highlights apparent gaps in our knowledge of the mechanisms contributing to resistance for these new drugs and challenges that limit the clinical utility of next generation sequencing techniques. It is recommended that a combination of genotypic and phenotypic techniques is warranted to monitor treatment response, curb emerging resistance and further dissemination of drug resistance.

A Whole-Genome Sequencing-Based Study to Delineate the Risk and Characteristics of Tuberculosis Transmission in an Insular Population Over 10 Years in Shanghai

Background

Tuberculosis (TB) has remained a tough problem in China. This study aims to identify the risk of tuberculosis transmission and to assess its characteristics.

Methods

We performed a molecular epidemiological study for patients with culture-positive Mycobacterium tuberculosis (M. tuberculosis) in Shanghai, from 2009 to 2018. Demographic information was obtained from the Tuberculosis Information Management System. Whole-genome sequencing (WGS) was conducted with a threshold of 12 single-nucleotide polymorphisms (SNPs) to distinguish the genomic cluster. To analyze the characteristics of TB transmission, the contact investigation for clustered cases was performed.

Results

In total, 94 (27.25%) of the 345 enrolled patients were grouped into 42 genomic clusters, indicating local transmission of M. tuberculosis strains. Compared to a health system delay <14 days, patients with a health system delay ≥14 days [adjusted odds ratios (AOR) = 2.57, 95% confidence interval (CI): 1.34–4.95] were more likely to be clustered. Patients under 65 years old (AOR = 3.11, 95% CI: 1.76–5.49), residents (AOR = 2.43, 95% CI: 1.18–4.99), and Beijing genotype strains (AOR = 3.35, 95% CI: 1.32–8.53) were associated with increased risk of clustering. Interestingly, patients with resistance to isoniazid (AOR = 2.36, 95% CI: 1.15–4.88) had a higher risk of transmission. Sixteen confirmed/probable epidemiological links were identified. Local transmission of imported cases and household transmission were prominent.

Conclusion

Health system delay is a crucial factor for TB transmission. Patients with resistance to isoniazid should be priority targets for contact investigation to reduce transmission.

SAM-TB: a whole genome sequencing data analysis website for detection of Mycobacterium tuberculosis drug resistance and transmission

Whole genome sequencing (WGS) can provide insight into drug-resistance, transmission chains and the identification of outbreaks, but data analysis remains an obstacle to its routine clinical use. Although several drug-resistance prediction tools have appeared, until now no website integrates drug-resistance prediction with strain genetic relationships and species identification of nontuberculous mycobacteria (NTM). We have established a free, function-rich, user-friendly online platform for MTB WGS data analysis (SAM-TB, http://samtb.szmbzx.com) that integrates drug-resistance prediction for 17 antituberculosis drugs, detection of variants, analysis of genetic relationships and NTM species identification. The accuracy of SAM-TB in predicting drug-resistance was assessed using 3177 sequenced clinical isolates with results of phenotypic drug-susceptibility tests (pDST). Compared to pDST, the sensitivity of SAM-TB for detecting multidrug-resistant tuberculosis was 93.9% [95% confidence interval (CI) 92.6–95.1%] with specificity of 96.2% (95% CI 95.2–97.1%). SAM-TB also analyzes the genetic relationships between multiple strains by reconstructing phylogenetic trees and calculating pairwise single nucleotide polymorphism (SNP) distances to identify genomic clusters. The incorporated mlstverse software identifies NTM species with an accuracy of 98.2% and Kraken2 software can detect mixed MTB and NTM samples. SAM-TB also has the capacity to share both sequence data and analysis between users. SAM-TB is a multifunctional integrated website that uses WGS raw data to accurately predict antituberculosis drug-resistance profiles, analyze genetic relationships between multiple strains and identify NTM species and mixed samples containing both NTM and MTB. SAM-TB is a useful tool for guiding both treatment and epidemiological investigation.

Exploiting genomics to mitigate the public health impact of antimicrobial resistance

Antimicrobial resistance (AMR) is a major global public health threat, which has been largely driven by the excessive use of antimicrobials. Control measures are urgently needed to slow the trajectory of AMR but are hampered by an incomplete understanding of the interplay between pathogens, AMR encoding genes, and mobile genetic elements at a microbial level. These factors, combined with the human, animal, and environmental interactions that underlie AMR dissemination at a population level, make for a highly complex landscape. Whole-genome sequencing (WGS) and, more recently, metagenomic analyses have greatly enhanced our understanding of these processes, and these approaches are informing mitigation strategies for how we better understand and control AMR. This review explores how WGS techniques have advanced global, national, and local AMR surveillance, and how this improved understanding is being applied to inform solutions, such as novel diagnostic methods that allow antimicrobial use to be optimised and vaccination strategies for better controlling AMR. We highlight some future opportunities for AMR control informed by genomic sequencing, along with the remaining challenges that must be overcome to fully realise the potential of WGS approaches for international AMR control.

Evaluating the Quality of Latent Tuberculosis Infection Screening in Ireland: A Single-Centre Retrospective Cohort Study

Ireland is a country with a low incidence of tuberculosis (TB) (5.6 cases per 100,000 population in 2019) that should be aiming for TB elimination (fewer than 1 case per million of population). To achieve TB elimination in low-incidence countries, programmatic latent tuberculosis infection (LTBI) management is important. This requires high-quality latent tuberculosis infection (LTBI) screening. Aim: To assess the quality of LTBI screening in a tertiary centre in Ireland using a framework. Methods: A retrospective review of the health care records of patients screened for TB in a tertiary centre in Ireland using an interferon-gamma release assay (IGRA) between 2016 and 2018 was performed. Three domains from the Institute of Medicine framework for health care quality, effectiveness, efficiency, and equity, were applied to measure the quality of LTBI screening. Results: Forty patients had LTBI and an indication for treatment, of whom 20% (8/40) were not offered treatment by the health care provider, 2.5% (1/40) did not accept treatment, and 10% (4/40) did not complete treatment. Seventy-five percent (6/8) of patients not offered treatment were non-Irish. The cost of screening per LTBI case identified was EUR 2048. Conclusions: This study evaluated the quality of LTBI screening using a framework and identified that LTBI screening in this tertiary centre needs to be scaled and expanded, and that treatment initiation needs to be improved, particularly among non-Irish nationals.

What Can Genetic Relatedness Tell Us About Risk Factors for Tuberculosis Transmission?

BACKGROUND

To stop tuberculosis (TB), the leading infectious cause of death globally, we need to better understand transmission risk factors. Although many studies have identified associations between individual-level covariates and pathogen genetic relatedness, few have identified characteristics of transmission pairs or explored how closely covariates associated with genetic relatedness mirror those associated with transmission.

METHODS

We simulated a TB-like outbreak with pathogen genetic data and estimated odds ratios (ORs) to correlate each covariate and genetic relatedness. We used a naive Bayes approach to modify the genetic links and nonlinks to resemble the true links and nonlinks more closely and estimated modified ORs with this approach. We compared these two sets of ORs with the true ORs for transmission. Finally, we applied this method to TB data in Hamburg, Germany, and Massachusetts, USA, to find pair-level covariates associated with transmission.

RESULTS

Using simulations, we found that associations between covariates and genetic relatedness had the same relative magnitudes and directions as the true associations with transmission, but biased absolute magnitudes. Modifying the genetic links and nonlinks reduced the bias and increased the confidence interval widths, more accurately capturing error. In Hamburg and Massachusetts, pairs were more likely to be probable transmission links if they lived in closer proximity, had a shorter time between observations, or had shared ethnicity, social risk factors, drug resistance, or genotypes.

CONCLUSIONS

We developed a method to improve the use of genetic relatedness as a proxy for transmission, and aid in understanding TB transmission dynamics in low-burden settings.

The Relevance of Genomic Epidemiology for Control of Tuberculosis in West Africa

Tuberculosis (TB), an airborne infectious disease caused by Mycobacterium tuberculosis complex (MTBC), remains a global health problem. West Africa has a unique epidemiology of TB that is characterized by medium- to high-prevalence. Moreover, the geographical restriction of M. africanum to the sub-region makes West Africa have an extra burden to deal with a two-in-one pathogen. The region is also burdened with low case detection, late reporting, poor treatment adherence leading to development of drug resistance and relapse. Sporadic studies conducted within the subregion report higher burden of drug resistant TB (DRTB) than previously thought. The need for more sensitive and robust tools for routine surveillance as well as to understand the mechanisms of DRTB and transmission dynamics for the design of effective control tools, cannot be overemphasized. The advancement in molecular biology tools including traditional fingerprinting and next generation sequencing (NGS) technologies offer reliable tools for genomic epidemiology. Genomic epidemiology provides in-depth insight of the nature of pathogens, circulating strains and their spread as well as prompt detection of the emergence of new strains. It also offers the opportunity to monitor treatment and evaluate interventions. Furthermore, genomic epidemiology can be used to understand potential emergence and spread of drug resistant strains and resistance mechanisms allowing the design of simple but rapid tools. In this review, we will describe the local epidemiology of MTBC, highlight past and current investigations toward understanding their biology and spread as well as discuss the relevance of genomic epidemiology studies to TB control in West Africa.

Import and transmission of Mycobacterium orygis and Mycobacterium africanum, Norway

Background

The aim of the current study was to improve our understanding of the origins and transmission of Mycobacterium africanum (MAF) in Norway.

Methods

Whole-genome sequences (WGS) were generated for all (n = 29) available clinical isolates received at the Norwegian National Reference Laboratory for Mycobacteria (NRL) and identified as MAF in Norway, in the period 2010–2020. Phylogenetic analyses were performed.

Results

The analyses indicated several imports of MAF lineage 6 from both East and West African countries, whereas MAF lineage 5 was restricted to patients with West African connections. We also find evidence for transmission of MAF in Norway. Finally, our analyses revealed that a group of isolates from patients originating in South Asia, identified as MAF by means of a commercial line-probe assay, in fact belonged to Mycobacterium orygis.

Conclusions

Most MAF cases in Norway are the result of import, but transmission is occurring within Norway.

Integrative transnational analysis to dissect tuberculosis transmission events along the migratory route from Africa to Europe

BACKGROUND

Growing international migration has increased the complexity of tuberculosis transmission patterns. Italy's decision to close its borders in 2018 made of Spain the new European porte entrée for migration from the Horn of Africa (HA). In one of the first rescues of migrants from this region at the end of 2018, tuberculosis was diagnosed in eight subjects, mainly unaccompanied minors.

METHODS

Mycobacterium tuberculosis isolates from these recently arrived migrants were analysed by Mycobacterial Interspersed Repetitive-Unit/Variable-Number of Tandem Repeat (MIRU-VNTR) and subsequent whole genome sequencing (WGS) analysis. Data were compared with those from collections from other European countries receiving migrants from the HA and a strain-specific PCR was applied for a fast searching of common strains. Infections in a cellular model were performed to assess strain virulence.

RESULTS

MIRU-VNTR analysis allowed identifying an epidemiological cluster involving three of the eight cases from Somalia (0 single-nucleotide polymorphisms between isolates, HA cluster). Following detailed interviews revealed that two of these cases had shared the same migratory route in most of the trip and had spent a long time at a detention camp in Libya. To confirm potential en route transmission for the three cases, we searched the same strain in collections from other European countries receiving migrants from the HA. MIRU-VNTR, WGS and a strain-specific PCR for the HA strain were applied. The same strain was identified in 12 cases from Eritrea diagnosed soon after their arrival in 2018 to the Netherlands, Belgium and Italy. Intracellular replication rate of the strain did not reveal abnormal virulence.

CONCLUSIONS

Our study suggests a potential en route transmission of a pan-susceptible strain, which caused at least 15 tuberculosis cases in Somalian and Eritrean migrants diagnosed in four different European countries.

Assessment of closely related Mycobacterium tuberculosis variants with different transmission success and in vitro infection dynamics

Whole genome sequencing (WGS) is able to differentiate closely related Mycobacterium tuberculosis variants within the same transmission cluster. Our aim was to evaluate if this higher discriminatory power may help identify and characterize more actively transmitted variants and understand the factors behind their success. We selected a robust MIRU-VNTR-defined cluster from Almería, Spain (22 cases throughout 2003–2019). WGS allowed discriminating, within the same epidemiological setting, between a successfully transmitted variant and seven closely related variants that did not lead to secondary cases, or were involved in self-limiting transmission (one single secondary case). Intramacrophagic growth of representative variants was evaluated in an in vitro infection model using U937 cells. Intramacrophage multiplication ratios (CFUs at Day 4/CFUs at Day 0) were higher for the actively transmitted variant (range 5.3–10.7) than for the unsuccessfully transmitted closely related variants (1.5–3.95). Two SNPs, mapping at the DNA binding domain of DnaA and at kdpD, were found to be specific of the successful variant.

Estimating the relative probability of direct transmission between infectious disease patients

BACKGROUND

Estimating infectious disease parameters such as the serial interval (time between symptom onset in primary and secondary cases) and reproductive number (average number of secondary cases produced by a primary case) are important in understanding infectious disease dynamics. Many estimation methods require linking cases by direct transmission, a difficult task for most diseases.

METHODS

Using a subset of cases with detailed genetic and/or contact investigation data to develop a training set of probable transmission events, we build a model to estimate the relative transmission probability for all case-pairs from demographic, spatial and clinical data. Our method is based on naive Bayes, a machine learning classification algorithm which uses the observed frequencies in the training dataset to estimate the probability that a pair is linked given a set of covariates.

RESULTS

In simulations, we find that the probabilities estimated using genetic distance between cases to define training transmission events are able to distinguish between truly linked and unlinked pairs with high accuracy (area under the receiver operating curve value of 95%). Additionally, only a subset of the cases, 10-50% depending on sample size, need to have detailed genetic data for our method to perform well. We show how these probabilities can be used to estimate the average effective reproductive number and apply our method to a tuberculosis outbreak in Hamburg, Germany.

CONCLUSIONS

Our method is a novel way to infer transmission dynamics in any dataset when only a subset of cases has rich contact investigation and/or genetic data.

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.

Genomic Epidemiology Analysis of Infectious Disease Outbreaks Using TransPhylo

Comparing the pathogen genomes from several cases of an infectious disease has the potential to help us understand and control outbreaks. Many methods exist to reconstruct a phylogeny from such genomes, which represents how the genomes are related to one another. However, such a phylogeny is not directly informative about transmission events between individuals. TransPhylo is a software tool implemented as an R package designed to bridge the gap between pathogen phylogenies and transmission trees. TransPhylo is based on a combined model of transmission between hosts and pathogen evolution within each host. It can simulate both phylogenies and transmission trees jointly under this combined model. TransPhylo can also reconstruct a transmission tree based on a dated phylogeny, by exploring the space of transmission trees compatible with the phylogeny. A transmission tree can be represented as a coloring of a phylogeny where each color represents a different host of the pathogen, and TransPhylo provides convenient ways to plot these colorings and explore the results. This article presents the basic protocols that can be used to make the most of TransPhylo. © 2021 The Authors. Basic Protocol 1: First steps with TransPhylo Basic Protocol 2: Simulation of outbreak data Basic Protocol 3: Inference of transmission Basic Protocol 4: Exploring the results of inference.

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

There were 290 multidrug-resistant (MDR)-TB cases diagnosed in Singapore from 2006 to 2018. Eighty-one percent were foreign-born. Spoligotyping and MIRU-VNTR methods identified 108 patients in 24 clusters. The Beijing spoligotype accounted for 22 clusters. Whole genome sequencing (WGS) analysis reduced the number of clustered patients and clusters to 43 and nine respectively. One MIRU cluster was redefined into three WGS clusters. All the clusters had foreign-born source cases. Forty percent of local-born, versus 9% of foreign-born, MDR-TB cases belonged to WGS clusters. WGS more accurately elucidated potential MDR-TB transmission which was overestimated by conventional genotyping methods in Singapore.

Case Report: Disseminated, rifampicin resistant Mycobacterium bovis (BCG) infection in an immunocompromised child

Background: Bacillus Calmette-Guérin (BCG) is a live-attenuated vaccine used world-wide for prevention of tuberculosis disease. In some immunocompromised hosts it has the potential to cause disease. As with other members of the M. tuberculosis complex it has the potential for acquiring drug resistance. Methods: We reviewed 10 years of paediatric clinical BCG strains referred to our clinical microbiology laboratory in Oxford where they underwent whole genome sequencing. We present a case series comparing clinical, pathogen genetic and pathogen phenotypic data, and consider the clinical implications. Results: We identified 15 BCG isolates from 8 children under 16 years old. Only one child had clinical disease with the other seven reported as local inoculation-site reactions. Case 1 suffered disseminated disease secondary to an undiagnosed IL-12/IFNγ receptor defect and the BCG isolates evolved two different rifampicin resistance mutations. Across all 15 isolates, phenotypic resistance to each first line drug was seen.  Conclusions: BCG is a safe and effective vaccine in children. Most clinical specimens in our series were not related to disease. However, in the context of rare pathogen-specific immunocompromise, BCG can cause pathology and acquire drug resistance under selection from therapy.

DNA Thermo-Protection Facilitates Whole-Genome Sequencing of Mycobacteria Direct from Clinical Samples

Mycobacterium tuberculosis is the leading cause of death from bacterial infection. Improved rapid diagnosis and antimicrobial resistance determination, such as by whole-genome sequencing, are required. Our aim was to develop a simple, low-cost method of preparing DNA for sequencing direct from M. tuberculosis-positive clinical samples (without culture). Simultaneous sputum liquefaction, bacteria heat inactivation (99°C/30 min), and enrichment for mycobacteria DNA were achieved using an equal volume of thermo-protection buffer (4 M KCl, 0.

Mycobacterium tuberculosis is the leading cause of death from bacterial infection. Improved rapid diagnosis and antimicrobial resistance determination, such as by whole-genome sequencing, are required. Our aim was to develop a simple, low-cost method of preparing DNA for sequencing direct from M. tuberculosis-positive clinical samples (without culture). Simultaneous sputum liquefaction, bacteria heat inactivation (99°C/30 min), and enrichment for mycobacteria DNA were achieved using an equal volume of thermo-protection buffer (4 M KCl, 0.05 M HEPES buffer, pH 7.5, 0.1% dithiothreitol [DTT]). The buffer emulated intracellular conditions found in hyperthermophiles, thus protecting DNA from rapid thermodegradation, which renders it a poor template for sequencing. Initial validation experiments employed mycobacteria DNA, either extracted or intracellular. Next, mock clinical samples (infection-negative human sputum spiked with 0 to 10⁵Mycobacterium bovis BCG cells/ml) underwent liquefaction in thermo-protection buffer and heat inactivation. DNA was extracted and sequenced. Human DNA degraded faster than mycobacteria DNA, resulting in target enrichment. Four replicate experiments achieved M. tuberculosis detection at 10¹ BCG cells/ml, with 31 to 59 M. tuberculosis complex reads. Maximal genome coverage (>97% at 5× depth) occurred at 10⁴ BCG cells/ml; >91% coverage (1× depth) occurred at 10³ BCG cells/ml. Final validation employed M. tuberculosis-positive clinical samples (n = 20), revealing that initial sample volumes of ≥1 ml typically yielded higher mean depths of M. tuberculosis genome coverage, with an overall range of 0.55 to 81.02. A mean depth of 3 gave >96% 1-fold tuberculosis (TB) genome coverage (in 15/20 clinical samples). A mean depth of 15 achieved >99% 5-fold genome coverage (in 9/20 clinical samples). In summary, direct-from-sample sequencing of M. tuberculosis genomes was facilitated by a low-cost thermo-protection buffer.

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

Background

Whole genome sequencing (WGS) is a reliable tool for studying tuberculosis (TB) transmission. WGS data are usually processed by custom-built analysis pipelines with little standardisation between them.

Aim

To compare the impact of variability of several WGS analysis pipelines used internationally to detect epidemiologically linked TB cases.

Methods

From the Netherlands, 535 Mycobacterium tuberculosis complex (MTBC) strains from 2016 were included. Epidemiological information obtained from municipal health services was available for all mycobacterial interspersed repeat unit-variable number of tandem repeat (MIRU-VNTR) clustered cases. WGS data was analysed using five different pipelines: one core genome multilocus sequence typing (cgMLST) approach and four single nucleotide polymorphism (SNP)-based pipelines developed in Oxford, United Kingdom; Borstel, Germany; Bilthoven, the Netherlands and Copenhagen, Denmark. WGS clusters were defined using a maximum pairwise distance of 12 SNPs/alleles.

Results

The cgMLST approach and Oxford pipeline clustered all epidemiologically linked cases, however, in the other three SNP-based pipelines one epidemiological link was missed due to insufficient coverage. In general, the genetic distances varied between pipelines, reflecting different clustering rates: the cgMLST approach clustered 92 cases, followed by 84, 83, 83 and 82 cases in the SNP-based pipelines from Copenhagen, Oxford, Borstel and Bilthoven respectively.

Conclusion

Concordance in ruling out epidemiological links was high between pipelines, which is an important step in the international validation of WGS data analysis. To increase accuracy in identifying TB transmission clusters, standardisation of crucial WGS criteria and creation of a reference database of representative MTBC sequences would be advisable.

Genomic variant-identification methods may alter Mycobacterium tuberculosis transmission inferences

Pathogen genomic data are increasingly used to characterize global and local transmission patterns of important human pathogens and to inform public health interventions. Yet, there is no current consensus on how to measure genomic variation. To test the effect of the variant-identification approach on transmission inferences for Mycobacterium tuberculosis, we conducted an experiment in which five genomic epidemiology groups applied variant-identification pipelines to the same outbreak sequence data. We compared the variants identified by each group in addition to transmission and phylogenetic inferences made with each variant set. To measure the performance of commonly used variant-identification tools, we simulated an outbreak. We compared the performance of three mapping algorithms, five variant callers and two variant filters in recovering true outbreak variants. Finally, we investigated the effect of applying increasingly stringent filters on transmission inferences and phylogenies. We found that variant-calling approaches used by different groups do not recover consistent sets of variants, which can lead to conflicting transmission inferences. Further, performance in recovering true variation varied widely across approaches. While no single variant-identification approach outperforms others in both recovering true genome-wide and outbreak-level variation, variant-identification algorithms calibrated upon real sequence data or that incorporate local reassembly outperform others in recovering true pairwise differences between isolates. The choice of variant filters contributed to extensive differences across pipelines, and applying increasingly stringent filters rapidly eroded the accuracy of transmission inferences and quality of phylogenies reconstructed from outbreak variation. Commonly used approaches to identify M. tuberculosis genomic variation have variable performance, particularly when predicting potential transmission links from pairwise genetic distances. Phylogenetic reconstruction may be improved by less stringent variant filtering. Approaches that improve variant identification in repetitive, hypervariable regions, such as long-read assemblies, may improve transmission inference.