The relationship between transmission time and clustering methods in Mycobacterium tuberculosis epidemiology

YesBackground: Tracking recent transmission is a vital part of controlling widespread pathogens such as Mycobacterium tuberculosis. Multiple methods with specific performance characteristics exist for detecting recent transmission chains, usually by clustering strains based on genotype similarities. With such a large variety of methods available, informed selection of an appropriate approach for determining transmissions within a given setting/time period is difficult. Methods: This study combines whole genome sequence (WGS) data derived from 324 isolates collected 2005–2010 in Kinshasa, Democratic Republic of Congo (DRC), a high endemic setting, with phylodynamics to unveil the timing of transmission events posited by a variety of standard genotyping methods. Clustering data based on Spoligotyping, 24-loci MIRU-VNTR typing, WGS based SNP (Single Nucleotide Polymorphism) and core genome multi locus sequence typing (cgMLST) typing were evaluated. Findings: Our results suggest that clusters based on Spoligotyping could encompass transmission events that occurred almost 200 years prior to sampling while 24-loci-MIRU-VNTR often represented three decades of transmission. Instead, WGS based genotyping applying low SNP or cgMLST allele thresholds allows for determination of recent transmission events, e.g. in timespans of up to 10 years for a 5 SNP/allele cut-off. Interpretation: With the rapid uptake of WGS methods in surveillance and outbreak tracking, the findings obtained in this study can guide the selection of appropriate clustering methods for uncovering relevant transmission chains within a given time-period. For high resolution cluster analyses, WGS-SNP and cgMLST based analyses have similar clustering/timing characteristics even for data obtained from a high incidence setting.ERC grant [INTERRUPTB; no. 311725] to BdJ, FG and CJM; an ERC grant to TS [PhyPD; no. 335529]; an FWO PhD fellowship to PM [grant number 1141217N]; the Leibniz Science Campus EvolLUNG for MM and SN; the German Centre for Infection Research (DZIF) for TAK, MM, CU, PB and SN; a SNF SystemsX grant (TBX) to JP and TS and a Marie Heim-Vögtlin fellowship granted to DK by the Swiss National Science Foundation. The computational resources and services used in this work were provided by the VSC (Flemish Supercomputer Center), funded by the Research Foundation - Flanders (FWO) and the Flemish Government – department EWI

Influence of hepatic fibrosis and inflammation: Correlation between histopathological changes and Gd-EOB-DTPA-enhanced MR imaging

Objective To evaluate the influence of an active inflammatory process in the liver on Gd-EOB-DTPA-enhanced MR imaging in patients with different degrees of fibrosis/cirrhosis. Material and methods Overall, a number of 91 patients (61 men and 30 women; mean age 58 years) were included in this retrospective study. The inclusion criteria for this study were Gd-EOB-DTPA-enhanced MRI of the liver and histopathological evaluation of fibrotic and inflammatory changes. T1-weighted VIBE sequences of the liver with fat suppression were evaluated to determine the relative signal change (RE) between native and hepatobiliary phase (20 min). In simple and multiple linear regression analyses, the influence of liver fibrosis/cirrhosis (Ishak score) and the histopathological degree of hepatitis (Modified Hepatic Activity Index, mHAI) on RE were evaluated. Results RE decreased significantly with increasing liver fibrosis/cirrhosis (p < 0.001) and inflammation (mHAI, p = 0.004). In particular, a correlation between RE and periportal or periseptal boundary zone hepatitis (moth feeding necrosis, mHAI A, p = 0.001) and portal inflammation (mHAI D, p < 0.001) was observed. In multiple linear regression analysis, both the degree of inflammation and the degree of fibrosis were significant predictors for RE (p < 0.01). Conclusion The results of this study suggest that the MR-based hepatic enhancement index RE is not only influenced by the degree of fibrosis, but also by the degree of inflammation

Ancient and recent differences in the intrinsic susceptibility of Mycobacterium tuberculosis complex to pretomanid

OBJECTIVES: To develop a robust phenotypic antimicrobial susceptibility testing (AST) method with a correctly set breakpoint for pretomanid (Pa), the most recently approved anti-tuberculosis drug. METHODS: The Becton Dickinson Mycobacterial Growth Indicator Tube™ (MGIT) system was used at six laboratories to determine the MICs of a phylogenetically diverse collection of 356 Mycobacterium tuberculosis complex (MTBC) strains to establish the epidemiological cut-off value for pretomanid. MICs were correlated with WGS data to study the genetic basis of differences in the susceptibility to pretomanid. RESULTS: We observed ancient differences in the susceptibility to pretomanid among various members of MTBC. Most notably, lineage 1 of M. tuberculosis, which is estimated to account for 28% of tuberculosis cases globally, was less susceptible than lineages 2, 3, 4 and 7 of M. tuberculosis, resulting in a 99th percentile of 2 mg/L for lineage 1 compared with 0.5 mg/L for the remaining M. tuberculosis lineages. Moreover, we observed that higher MICs (≥8 mg/L), which probably confer resistance, had recently evolved independently in six different M. tuberculosis strains. Unlike the aforementioned ancient differences in susceptibility, these recent differences were likely caused by mutations in the known pretomanid resistance genes. CONCLUSIONS: In light of these findings, the provisional critical concentration of 1 mg/L for MGIT set by EMA must be re-evaluated. More broadly, these findings underline the importance of considering the global diversity of MTBC during clinical development of drugs and when defining breakpoints for AST

A Novel Metagenomic Short-Chain Dehydrogenase/Reductase Attenuates Pseudomonas aeruginosa Biofilm Formation and Virulence on Caenorhabditis elegans

In Pseudomonas aeruginosa, the expression of a number of virulence factors, as well as biofilm formation, are controlled by quorum sensing (QS). N-Acylhomoserine lactones (AHLs) are an important class of signaling molecules involved in bacterial QS and in many pathogenic bacteria infection and host colonization are AHL-dependent. The AHL signaling molecules are subject to inactivation mainly by hydrolases (Enzyme Commission class number EC 3) (i.e. N-acyl-homoserine lactonases and N-acyl-homoserine-lactone acylases). Only little is known on quorum quenching mechanisms of oxidoreductases (EC 1). Here we report on the identification and structural characterization of the first NADP-dependent short-chain dehydrogenase/reductase (SDR) involved in inactivation of N-(3-oxo-dodecanoyl)-L-homoserine lactone (3-oxo-C12-HSL) and derived from a metagenome library. The corresponding gene was isolated from a soil metagenome and designated bpiB09. Heterologous expression and crystallographic studies established BpiB09 as an NADP-dependent reductase. Although AHLs are probably not the native substrate of this metagenome-derived enzyme, its expression in P. aeruginosa PAO1 resulted in significantly reduced pyocyanin production, decreased motility, poor biofilm formation and absent paralysis of Caenorhabditis elegans. Furthermore, a genome-wide transcriptome study suggested that the level of lasI and rhlI transcription together with 36 well known QS regulated genes was significantly (≥10-fold) affected in P. aeruginosa strains expressing the bpiB09 gene in pBBR1MCS-5. Thus AHL oxidoreductases could be considered as potent tools for the development of quorum quenching strategies

Ancient and recent differences in the intrinsic susceptibility of Mycobacterium tuberculosis complex to pretomanid

Objectives To develop a robust phenotypic antimicrobial susceptibility testing (AST) method with a correctly set breakpoint for pretomanid (Pa), the most recently approved anti-tuberculosis drug. Methods The Becton Dickinson Mycobacterial Growth Indicator Tube™ (MGIT) system was used at six laboratories to determine the MICs of a phylogenetically diverse collection of 356 Mycobacterium tuberculosis complex (MTBC) strains to establish the epidemiological cut-off value for pretomanid. MICs were correlated with WGS data to study the genetic basis of differences in the susceptibility to pretomanid. Results We observed ancient differences in the susceptibility to pretomanid among various members of MTBC. Most notably, lineage 1 of M. tuberculosis, which is estimated to account for 28% of tuberculosis cases globally, was less susceptible than lineages 2, 3, 4 and 7 of M. tuberculosis, resulting in a 99th percentile of 2 mg/L for lineage 1 compared with 0.5 mg/L for the remaining M. tuberculosis lineages. Moreover, we observed that higher MICs (≥8 mg/L), which probably confer resistance, had recently evolved independently in six different M. tuberculosis strains. Unlike the aforementioned ancient differences in susceptibility, these recent differences were likely caused by mutations in the known pretomanid resistance genes. Conclusions In light of these findings, the provisional critical concentration of 1 mg/L for MGIT set by EMA must be re-evaluated. More broadly, these findings underline the importance of considering the global diversity of MTBC during clinical development of drugs and when defining breakpoints for AST

MDR M. tuberculosis outbreak clone in Eswatini missed by Xpert has elevated bedaquiline resistance dated to the pre-treatment era.

BACKGROUND: Multidrug-resistant (MDR) Mycobacterium tuberculosis complex strains not detected by commercial molecular drug susceptibility testing (mDST) assays due to the RpoB I491F resistance mutation are threatening the control of MDR tuberculosis (MDR-TB) in Eswatini. METHODS: We investigate the evolution and spread of MDR strains in Eswatini with a focus on bedaquiline (BDQ) and clofazimine (CFZ) resistance using whole-genome sequencing in two collections ((1) national drug resistance survey, 2009-2010; (2) MDR strains from the Nhlangano region, 2014-2017). RESULTS: MDR strains in collection 1 had a high cluster rate (95%, 117/123 MDR strains) with 55% grouped into the two largest clusters (gCL3, n = 28; gCL10, n = 40). All gCL10 isolates, which likely emerged around 1993 (95% highest posterior density 1987-1998), carried the mutation RpoB I491F that is missed by commercial mDST assays. In addition, 21 (53%) gCL10 isolates shared a Rv0678 M146T mutation that correlated with elevated minimum inhibitory concentrations (MICs) to BDQ and CFZ compared to wild type isolates. gCL10 isolates with the Rv0678 M146T mutation were also detected in collection 2. CONCLUSION: The high clustering rate suggests that transmission has been driving the MDR-TB epidemic in Eswatini for three decades. The presence of MDR strains in Eswatini that are not detected by commercial mDST assays and have elevated MICs to BDQ and CFZ potentially jeopardizes the successful implementation of new MDR-TB treatment guidelines. Measures to limit the spread of these outbreak isolates need to be implemented urgently

Ancient and recent differences in the intrinsic susceptibility of Mycobacterium tuberculosis complex to pretomanid.

OBJECTIVES: To develop a robust phenotypic antimicrobial susceptibility testing (AST) method with a correctly set breakpoint for pretomanid (Pa), the most recently approved anti-tuberculosis drug. METHODS: The Becton Dickinson Mycobacterial Growth Indicator Tube™ (MGIT) system was used at six laboratories to determine the MICs of a phylogenetically diverse collection of 356 Mycobacterium tuberculosis complex (MTBC) strains to establish the epidemiological cut-off value for pretomanid. MICs were correlated with WGS data to study the genetic basis of differences in the susceptibility to pretomanid. RESULTS: We observed ancient differences in the susceptibility to pretomanid among various members of MTBC. Most notably, lineage 1 of M. tuberculosis, which is estimated to account for 28% of tuberculosis cases globally, was less susceptible than lineages 2, 3, 4 and 7 of M. tuberculosis, resulting in a 99th percentile of 2 mg/L for lineage 1 compared with 0.5 mg/L for the remaining M. tuberculosis lineages. Moreover, we observed that higher MICs (≥8 mg/L), which probably confer resistance, had recently evolved independently in six different M. tuberculosis strains. Unlike the aforementioned ancient differences in susceptibility, these recent differences were likely caused by mutations in the known pretomanid resistance genes. CONCLUSIONS: In light of these findings, the provisional critical concentration of 1 mg/L for MGIT set by EMA must be re-evaluated. More broadly, these findings underline the importance of considering the global diversity of MTBC during clinical development of drugs and when defining breakpoints for AST

A cluster of multidrug-resistant Mycobacterium tuberculosis among patients arriving in Europe from the Horn of Africa: a molecular epidemiological study

SummaryBackground The risk of tuberculosis outbreaks among people fleeing hardship for refuge in Europe is heightened. We describe the cross-border European response to an outbreak of multidrug-resistant tuberculosis among patients from the Horn of Africa and Sudan. Methods On April 29 and May 30, 2016, the Swiss and German National Mycobacterial Reference Laboratories independently triggered an outbreak investigation after four patients were diagnosed with multidrug-resistant tuberculosis. In this molecular epidemiological study, we prospectively defined outbreak cases with 24-locus mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) profiles; phenotypic resistance to isoniazid, rifampicin, ethambutol, pyrazinamide, and capreomycin; and corresponding drug resistance mutations. We whole-genome sequenced all Mycobacterium tuberculosis isolates and clustered them using a threshold of five single nucleotide polymorphisms (SNPs). We collated epidemiological data from host countries from the European Centre for Disease Prevention and Control. Findings Between Feb 12, 2016, and April 19, 2017, 29 patients were diagnosed with multidrug-resistant tuberculosis in seven European countries. All originated from the Horn of Africa or Sudan, with all isolates two SNPs or fewer apart. 22 (76%) patients reported their travel routes, with clear spatiotemporal overlap between routes. We identified a further 29 MIRU-VNTR-linked cases from the Horn of Africa that predated the outbreak, but all were more than five SNPs from the outbreak. However all 58 isolates shared a capreomycin resistance-associated tlyA mutation. Interpretation Our data suggest that source cases are linked to an M tuberculosis clone circulating in northern Somalia or Djibouti and that transmission probably occurred en route before arrival in Europe. We hypothesise that the shared mutation of tlyA is a drug resistance mutation and phylogenetic marker, the first of its kind in M tuberculosis sensu stricto. Funding The Swiss Federal Office of Public Health, the University of Zurich, the Wellcome Trust, National Institute for Health Research (NIHR) Oxford Biomedical Research Centre (BRC), the Medical Research Council, BELTA-TBnet, the European Union, the German Center for Infection Research, and Leibniz Science Campus Evolutionary Medicine of the Lung (EvoLUNG)

The 2021 WHO catalogue of Mycobacterium tuberculosis complex mutations associated with drug resistance: a genotypic analysis.

Background: Molecular diagnostics are considered the most promising route to achievement of rapid, universal drug susceptibility testing for Mycobacterium tuberculosis complex (MTBC). We aimed to generate a WHO-endorsed catalogue of mutations to serve as a global standard for interpreting molecular information for drug resistance prediction. Methods: In this systematic analysis, we used a candidate gene approach to identify mutations associated with resistance or consistent with susceptibility for 13 WHO-endorsed antituberculosis drugs. We collected existing worldwide MTBC whole-genome sequencing data and phenotypic data from academic groups and consortia, reference laboratories, public health organisations, and published literature. We categorised phenotypes as follows: methods and critical concentrations currently endorsed by WHO (category 1); critical concentrations previously endorsed by WHO for those methods (category 2); methods or critical concentrations not currently endorsed by WHO (category 3). For each mutation, we used a contingency table of binary phenotypes and presence or absence of the mutation to compute positive predictive value, and we used Fisher's exact tests to generate odds ratios and Benjamini-Hochberg corrected p values. Mutations were graded as associated with resistance if present in at least five isolates, if the odds ratio was more than 1 with a statistically significant corrected p value, and if the lower bound of the 95% CI on the positive predictive value for phenotypic resistance was greater than 25%. A series of expert rules were applied for final confidence grading of each mutation. Findings: We analysed 41 137 MTBC isolates with phenotypic and whole-genome sequencing data from 45 countries. 38 215 MTBC isolates passed quality control steps and were included in the final analysis. 15 667 associations were computed for 13 211 unique mutations linked to one or more drugs. 1149 (7·3%) of 15 667 mutations were classified as associated with phenotypic resistance and 107 (0·7%) were deemed consistent with susceptibility. For rifampicin, isoniazid, ethambutol, fluoroquinolones, and streptomycin, the mutations' pooled sensitivity was more than 80%. Specificity was over 95% for all drugs except ethionamide (91·4%), moxifloxacin (91·6%) and ethambutol (93·3%). Only two resistance mutations were identified for bedaquiline, delamanid, clofazimine, and linezolid as prevalence of phenotypic resistance was low for these drugs. Interpretation: We present the first WHO-endorsed catalogue of molecular targets for MTBC drug susceptibility testing, which is intended to provide a global standard for resistance interpretation. The existence of this catalogue should encourage the implementation of molecular diagnostics by national tuberculosis programmes. Funding: Unitaid, Wellcome Trust, UK Medical Research Council, and Bill and Melinda Gates Foundation

The evolving SARS-CoV-2 epidemic in Africa: Insights from rapidly expanding genomic surveillance

INTRODUCTION Investment in Africa over the past year with regard to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) sequencing has led to a massive increase in the number of sequences, which, to date, exceeds 100,000 sequences generated to track the pandemic on the continent. These sequences have profoundly affected how public health officials in Africa have navigated the COVID-19 pandemic. RATIONALE We demonstrate how the first 100,000 SARS-CoV-2 sequences from Africa have helped monitor the epidemic on the continent, how genomic surveillance expanded over the course of the pandemic, and how we adapted our sequencing methods to deal with an evolving virus. Finally, we also examine how viral lineages have spread across the continent in a phylogeographic framework to gain insights into the underlying temporal and spatial transmission dynamics for several variants of concern (VOCs). RESULTS Our results indicate that the number of countries in Africa that can sequence the virus within their own borders is growing and that this is coupled with a shorter turnaround time from the time of sampling to sequence submission. Ongoing evolution necessitated the continual updating of primer sets, and, as a result, eight primer sets were designed in tandem with viral evolution and used to ensure effective sequencing of the virus. The pandemic unfolded through multiple waves of infection that were each driven by distinct genetic lineages, with B.1-like ancestral strains associated with the first pandemic wave of infections in 2020. Successive waves on the continent were fueled by different VOCs, with Alpha and Beta cocirculating in distinct spatial patterns during the second wave and Delta and Omicron affecting the whole continent during the third and fourth waves, respectively. Phylogeographic reconstruction points toward distinct differences in viral importation and exportation patterns associated with the Alpha, Beta, Delta, and Omicron variants and subvariants, when considering both Africa versus the rest of the world and viral dissemination within the continent. Our epidemiological and phylogenetic inferences therefore underscore the heterogeneous nature of the pandemic on the continent and highlight key insights and challenges, for instance, recognizing the limitations of low testing proportions. We also highlight the early warning capacity that genomic surveillance in Africa has had for the rest of the world with the detection of new lineages and variants, the most recent being the characterization of various Omicron subvariants. CONCLUSION Sustained investment for diagnostics and genomic surveillance in Africa is needed as the virus continues to evolve. This is important not only to help combat SARS-CoV-2 on the continent but also because it can be used as a platform to help address the many emerging and reemerging infectious disease threats in Africa. In particular, capacity building for local sequencing within countries or within the continent should be prioritized because this is generally associated with shorter turnaround times, providing the most benefit to local public health authorities tasked with pandemic response and mitigation and allowing for the fastest reaction to localized outbreaks. These investments are crucial for pandemic preparedness and response and will serve the health of the continent well into the 21st century