Use of High-Resolution Geospatial and Genomic Data to Characterize Recent Tuberculosis Transmission, Botswana

Combining genomic and geospatial data can be useful for understanding Mycobacterium tuberculosis transmission in high-burden tuberculosis (TB) settings. We performed whole-genome sequencing on M. tuberculosis DNA extracted from sputum cultures from a population-based TB study conducted in Gaborone, Botswana, during 2012-2016. We determined spatial distribution of cases on the basis of shared genotypes among isolates. We considered clusters of isolates with ≤5 single-nucleotide polymorphisms identified by whole-genome sequencing to indicate recent transmission and clusters of ≥10 persons to be outbreaks. We obtained both molecular and geospatial data for 946/1,449 (65%) participants with culture-confirmed TB; 62 persons belonged to 5 outbreaks of 10-19 persons each. We detected geospatial clustering in just 2 of those 5 outbreaks, suggesting heterogeneous spatial patterns. Our findings indicate that targeted interventions applied in smaller geographic areas of high-burden TB identified using integrated genomic and geospatial data might help interrupt TB transmission during outbreaks.

Back-to-Africa introductions of Mycobacterium tuberculosis as the main cause of tuberculosis in Dar es Salaam, Tanzania

In settings with high tuberculosis (TB) endemicity, distinct genotypes of the Mycobacterium tuberculosis complex (MTBC) often differ in prevalence. However, the factors leading to these differences remain poorly understood. Here we studied the MTBC population in Dar es Salaam, Tanzania over a six-year period, using 1,082 unique patient-derived MTBC whole-genome sequences (WGS) and associated clinical data. We show that the TB epidemic in Dar es Salaam is dominated by multiple MTBC genotypes introduced to Tanzania from different parts of the world during the last 300 years. The most common MTBC genotypes deriving from these introductions exhibited differences in transmission rates and in the duration of the infectious period, but little differences in overall fitness, as measured by the effective reproductive number. Moreover, measures of disease severity and bacterial load indicated no differences in virulence between these genotypes during active TB. Instead, the combination of an early introduction and a high transmission rate accounted for the high prevalence of L3.1.1, the most dominant MTBC genotype in this setting. Yet, a longer co-existence with the host population did not always result in a higher transmission rate, suggesting that distinct life-history traits have evolved in the different MTBC genotypes. Taken together, our results point to bacterial factors as important determinants of the TB epidemic in Dar es Salaam.

Risk of tuberculosis transmission by children in Hamburg, Germany

BACKGROUND

Data from a prospective molecular-epidemiologic study (1997-2021) in Hamburg, Germany, were evaluated to assess the transmission risk of Mycobacterium tuberculosis complex (Mtbc) by children <15 years in a low-incidence setting.

METHODS

Isolates of Mtbc were genotyped whole genome sequencing, applying a core genome multilocus sequence typing scheme. Close contacts of culture-confirmed children were examined for latent Mtbc infections (LTBI) with particular focus on IGRA testing.

RESULTS

Out of 3154 culture-confirmed tuberculosis (TB) cases, 79 (2.5%) were children <15 years. Of those, 52 (58%) had pulmonary TB. Genotyping revealed that 35 of the 52 children (67%) were epidemiologically confirmed secondary cluster members; all of their source cases were adults. Six immigrant children presented without a presumed source case; their TB diagnoses came on average 48 weeks (interquartile range [IQR] 71) after their arrival in Germany. Three German-born children were determined to have been infected by adult relatives while visiting their parents' home country. Of the 317 children's close contacts tested with QuantiFERON-TB Gold-In Tube for LTBI, only 21 (6.6%) were positive. Absent a history of prior exposure or immigration from a high-incidence country, none of the contacts of younger (<10 years) TB-afflicted children was latently infected, whereas 2 older children infected 12 of their contacts, children and adults. During a mean observational period of 551 weeks (IQR 735) on average, no secondary TB cases appeared.

CONCLUSIONS

Children with pulmonary TB disease, especially those aged below 10 years, rarely transmit Mtbc to their close contacts in a low-incidence setting.

Clinical implications of molecular drug resistance testing for Mycobacterium tuberculosis: a 2023 TBnet/RESIST-TB consensus statement

Drug-resistant tuberculosis is a substantial health-care concern worldwide. Despite culture-based methods being considered the gold standard for drug susceptibility testing, molecular methods provide rapid information about the Mycobacterium tuberculosis mutations associated with resistance to anti-tuberculosis drugs. This consensus document was developed on the basis of a comprehensive literature search, by the TBnet and RESIST-TB networks, about reporting standards for the clinical use of molecular drug susceptibility testing. Review and the search for evidence included hand-searching journals and searching electronic databases. The panel identified studies that linked mutations in genomic regions of M tuberculosis with treatment outcome data. Implementation of molecular testing for the prediction of drug resistance in M tuberculosis is key. Detection of mutations in clinical isolates has implications for the clinical management of patients with multidrug-resistant or rifampicin-resistant tuberculosis, especially in situations when phenotypic drug susceptibility testing is not available. A multidisciplinary team including clinicians, microbiologists, and laboratory scientists reached a consensus on key questions relevant to molecular prediction of drug susceptibility or resistance to M tuberculosis, and their implications for clinical practice. This consensus document should help clinicians in the management of patients with tuberculosis, providing guidance for the design of treatment regimens and optimising outcomes.

Phylogroup-specific variation shapes the clustering of antimicrobial resistance genes and defence systems across regions of genome plasticity in Pseudomonas aeruginosa

BACKGROUND

Pseudomonas aeruginosa is an opportunistic pathogen consisting of three phylogroups (hereafter named A, B, and C). Here, we assessed phylogroup-specific evolutionary dynamics across available and also new P. aeruginosa genomes.

METHODS

In this genomic analysis, we first generated new genome assemblies for 18 strains of the major P. aeruginosa clone type (mPact) panel, comprising a phylogenetically diverse collection of clinical and environmental isolates for this species. Thereafter, we combined these new genomes with 1991 publicly available P. aeruginosa genomes for a phylogenomic and comparative analysis. We specifically explored to what extent antimicrobial resistance (AMR) genes, defence systems, and virulence genes vary in their distribution across regions of genome plasticity (RGPs) and "masked" (RGP-free) genomes, and to what extent this variation differs among the phylogroups.

FINDINGS

We found that members of phylogroup B possess larger genomes, contribute a comparatively larger number of pangenome families, and show lower abundance of CRISPR-Cas systems. Furthermore, AMR and defence systems are pervasive in RGPs and integrative and conjugative/mobilizable elements (ICEs/IMEs) from phylogroups A and B, and the abundance of these cargo genes is often significantly correlated. Moreover, inter- and intra-phylogroup interactions occur at the accessory genome level, suggesting frequent recombination events. Finally, we provide here the mPact panel of diverse P. aeruginosa strains that may serve as a valuable reference for functional analyses.

INTERPRETATION

Altogether, our results highlight distinct pangenome characteristics of the P. aeruginosa phylogroups, which are possibly influenced by variation in the abundance of CRISPR-Cas systems and are shaped by the differential distribution of other defence systems and AMR genes.

FUNDING

German Science Foundation, Max-Planck Society, Leibniz ScienceCampus Evolutionary Medicine of the Lung, BMBF program Medical Infection Genomics, Kiel Life Science Postdoc Award.

The effect of M. tuberculosis lineage on clinical phenotype

Eight lineages of Mycobacterium tuberculosis sensu stricto are described. Single-country or small observational data suggest differences in clinical phenotype between lineages. We present strain lineage and clinical phenotype data from 12,246 patients from 3 low-incidence and 5 high-incidence countries. We used multivariable logistic regression to explore the effect of lineage on site of disease and on cavities on chest radiography, given pulmonary TB; multivariable multinomial logistic regression to investigate types of extra-pulmonary TB, given lineage; and accelerated failure time and Cox proportional-hazards models to explore the effect of lineage on time to smear and culture-conversion. Mediation analyses quantified the direct effects of lineage on outcomes. Pulmonary disease was more likely among patients with lineage(L) 2, L3 or L4, than L1 (adjusted odds ratio (aOR) 1.79, (95% confidence interval 1.49-2.15), p<0.001; aOR=1.40(1.09-1.79), p=0.007; aOR=2.04(1.65-2.53), p<0.001, respectively). Among patients with pulmonary TB, those with L1 had greater risk of cavities on chest radiography versus those with L2 (aOR=0.69(0.57-0.83), p<0.001) and L4 strains (aOR=0.73(0.59-0.90), p=0.002). L1 strains were more likely to cause osteomyelitis among patients with extra-pulmonary TB, versus L2-4 (p=0.033, p=0.008 and p=0.049 respectively). Patients with L1 strains showed shorter time-to-sputum smear conversion than for L2. Causal mediation analysis showed the effect of lineage in each case was largely direct. The pattern of clinical phenotypes seen with L1 strains differed from modern lineages (L2-4). This has implications for clinical management and could influence clinical trial selection strategies.

Genomic epidemiology of SARS-CoV-2 during the first four waves in Mozambique

Mozambique reported the first case of coronavirus disease 2019 (COVID-19) in March 2020 and it has since spread to all provinces in the country. To investigate the introductions and spread of SARS-CoV-2 in Mozambique, 1 142 whole genome sequences sampled within Mozambique were phylogenetically analyzed against a globally representative set, reflecting the first 25 months of the epidemic. The epidemic in the country was marked by four waves of infection, the first associated with B.1 ancestral lineages, while the Beta, Delta, and Omicron Variants of Concern (VOCs) were responsible for most infections and deaths during the second, third, and fourth waves. Large-scale viral exchanges occurred during the latter three waves and were largely attributed to southern African origins. Not only did the country remain vulnerable to the introductions of new variants but these variants continued to evolve within the borders of the country. Due to the Mozambican health system already under constraint, and paucity of data in Mozambique, there is a need to continue to strengthen and support genomic surveillance in the country as VOCs and Variants of interests (VOIs) are often reported from the southern African region.

2355. Prison as a Driver of Recent Transmissions of Multidrug-Resistant Tuberculosis in Callao, Peru

Background

Peru has one of the highest rates of multidrug-resistant tuberculosis (MDR-TB) in the Latin America region, with Callao being one of the hot spots. We sought to identify resistance patterns and key drivers of recent MDR-TB transmission in Callao, Peru.

Methods

Cross-sectional study including clinical specimens identified as MDR-TB strains in Callao, Peru between April 2017 and December 2019. DNA was extracted for whole genome sequencing and data used for phylogenetic classification, clustering, and resistance causing mutation analyses. Recent transmission was defined based on an isolate-to-isolate distance of ≤5 (D5) single nucleotide polymorphisms (SNPs). We used logistic regression models to analyze the relationship between MDR-TB clustering and epidemiologic factors including age, sex, history of diabetes mellitus, HIV infection, illicit drug use, known TB contact, prior TB disease, and imprisonment.

Results

171 unique MDR-TB strains were included; 93% were assigned to lineage 4 and 7% to lineage 2. The most prevalent sublineage was 4.3.3 LAM (57%), followed by 4.3.4.2 LAM (10%) and 4.1.2.1 Haarlem (9%). In the dominant 4.3.3 LAM sublineage, concomitant resistance was common, including resistance mutations to pyrazinamide (92%), ethambutol (22%), ethionamide (23%), and quinolones (10%); 4 isolates harbored bedaquiline resistance mutations. Seventy-four percent of 4.3.3 LAM isolates were D5-clustered, with 53 (30%) isolates within a single dominant cluster. Male sex (odds ratio [OR], 3.6; 95% confidence interval [CI], 1.5 – 8.9), drug use (OR, 4.9; 95% CI, 1.9 – 12.7), and current or prior history of imprisonment (OR, 12; 95% CI, 3.3 – 43.5) were associated with the dominant D5 cluster in unadjusted analyses. History of imprisonment remained independently associated with the dominant D5 cluster in multivariate analyses (adjusted OR, 8.9; 95% CI, 1.6 – 50.6). History of imprisonment was also independently associated with being part of any D5 clusters in multivariate analyses (adjusted OR, 4.8; 95% CI, 1.2 – 20.3).

Conclusion

History of imprisonment was linked to current MDR-TB transmissions, indicating an important role of prisons in driving the MDR-TB epidemic.

Disclosures

Moises A. Huaman, MD, MSc, Gilead: Grant/Research Support|Insmed: Grant/Research Support.

Case-Finding Strategies for Drug-Resistant Tuberculosis: Protocol for a Scoping Review

BACKGROUND

Transmission of drug-resistant tuberculosis (DR-TB) is ongoing. Finding individuals with DR-TB and initiating treatment as early as possible is important to improve patient clinical outcomes and to break the chain of transmission to control the pandemic. To our knowledge systematic reviews assessing effectiveness, cost-effectiveness, acceptability, and feasibility of different case-finding strategies for DR-TB to inform research, policy, and practice have not been conducted, and it is unknown whether enough research exists to conduct such reviews. It is unknown whether case-finding strategies are similar for DR-TB and drug-susceptible TB and whether we can draw on findings from drug-susceptible reviews to inform decisions on case-finding strategies for DR-TB.

OBJECTIVE

This protocol aims to describe the available literature on case-finding for DR-TB and to describe case-finding strategies.

METHODS

We will screen systematic reviews, trials, qualitative studies, diagnostic test accuracy studies, and other primary research that specifically sought to improve DR-TB case detection. We will exclude studies that invited individuals seeking care for TB symptoms, those including individuals already diagnosed with TB, or laboratory-based studies. We will search the academic databases including MEDLINE, Embase, The Cochrane Library, Africa-Wide Information, CINAHL, Epistemonikos, and PROSPERO with no language or date restrictions. We will screen titles, abstracts, and full-text articles in duplicate. Data extraction and analyses will be performed using Excel (Microsoft Corp).

RESULTS

We will provide a narrative report with supporting figures or tables to summarize the data. A systems-based logic model, developed from a synthesis of case-finding strategies for drug-susceptible TB, will be used as a framework to describe different strategies, resulting pathways, and enhancements of pathways. The search will be conducted at the end of 2021. Title and abstract screening, full text screening, and data extraction will be undertaken from January to June 2022. Thereafter, analysis will be conducted, and results compiled.

CONCLUSIONS

This scoping review will chart existing literature on case-finding for DR-TB-this will help determine whether primary studies on effectiveness, cost-effectiveness, acceptability, and feasibility of different case-finding strategies for DR-TB exist and will help formulate potential questions for a systematic review. We will also describe case-finding strategies for DR-TB and how they fit into a model of case-finding pathways for drug-susceptible TB. This review has some limitations. One limitation is the diverse, inconsistent use of intervention terminology within the literature, which may result in missing relevant studies. Poor reporting of intervention strategies may also cause misunderstanding and misclassification of interventions. Lastly, case-finding strategies for DR-TB may not fit into a model developed from strategies for drug-susceptible TB. Nevertheless, such a situation will provide an opportunity to refine the model for future research. The review will guide further research to inform decisions on case-finding policies and practices for DR-TB.

INTERNATIONAL REGISTERED REPORT IDENTIFIER (IRRID)

DERR1-10.2196/40009.

Using genetic data to identify transmission risk factors: Statistical assessment and application to tuberculosis transmission

Identifying host factors that influence infectious disease transmission is an important step toward developing interventions to reduce disease incidence. Recent advances in methods for reconstructing infectious disease transmission events using pathogen genomic and epidemiological data open the door for investigation of host factors that affect onward transmission. While most transmission reconstruction methods are designed to work with densely sampled outbreaks, these methods are making their way into surveillance studies, where the fraction of sampled cases with sequenced pathogens could be relatively low. Surveillance studies that use transmission event reconstruction then use the reconstructed events as response variables (i.e., infection source status of each sampled case) and use host characteristics as predictors (e.g., presence of HIV infection) in regression models. We use simulations to study estimation of the effect of a host factor on probability of being an infection source via this multi-step inferential procedure. Using TransPhylo-a widely-used method for Bayesian estimation of infectious disease transmission events-and logistic regression, we find that low sensitivity of identifying infection sources leads to dilution of the signal, biasing logistic regression coefficients toward zero. We show that increasing the proportion of sampled cases improves sensitivity and some, but not all properties of the logistic regression inference. Application of these approaches to real world data from a population-based TB study in Botswana fails to detect an association between HIV infection and probability of being a TB infection source. We conclude that application of a pipeline, where one first uses TransPhylo and sparsely sampled surveillance data to infer transmission events and then estimates effects of host characteristics on probabilities of these events, should be accompanied by a realistic simulation study to better understand biases stemming from imprecise transmission event inference.

Investigating resistance in clinical Mycobacterium tuberculosis complex isolates with genomic and phenotypic antimicrobial susceptibility testing: a multicentre observational study

BACKGROUND

Whole-genome sequencing (WGS) of Mycobacterium tuberculosis complex has become an important tool in diagnosis and management of drug-resistant tuberculosis. However, data correlating resistance genotype with quantitative phenotypic antimicrobial susceptibility testing (AST) are scarce.

METHODS

In a prospective multicentre observational study, 900 clinical M tuberculosis complex isolates were collected from adults with drug-resistant tuberculosis in five high-endemic tuberculosis settings around the world (Georgia, Moldova, Peru, South Africa, and Viet Nam) between Dec 5, 2014, and Dec 12, 2017. Minimum inhibitory concentrations (MICs) and resulting binary phenotypic AST results for up to nine antituberculosis drugs were determined and correlated with resistance-conferring mutations identified by WGS.

FINDINGS

Considering WHO-endorsed critical concentrations as reference, WGS had high accuracy for prediction of resistance to isoniazid (sensitivity 98·8% [95% CI 98·5-99·0]; specificity 96·6% [95% CI 95·2-97·9]), levofloxacin (sensitivity 94·8% [93·3-97·6]; specificity 97·1% [96·7-97·6]), kanamycin (sensitivity 96·1% [95·4-96·8]; specificity 95·0% [94·4-95·7]), amikacin (sensitivity 97·2% [96·4-98·1]; specificity 98·6% [98·3-98·9]), and capreomycin (sensitivity 93·1% [90·0-96·3]; specificity 98·3% [98·0-98·7]). For rifampicin, pyrazinamide, and ethambutol, the specificity of resistance prediction was suboptimal (64·0% [61·0-67·1], 83·8% [81·0-86·5], and 40·1% [37·4-42·9], respectively). Specificity for rifampicin increased to 83·9% when borderline mutations with MICs overlapping with the critical concentration were excluded. Consequently, we highlighted mutations in M tuberculosis complex isolates that are often falsely identified as susceptible by phenotypic AST, and we identified potential novel resistance-conferring mutations.

INTERPRETATION

The combined analysis of mutations and quantitative phenotypes shows the potential of WGS to produce a refined interpretation of resistance, which is needed for individualised therapy, and eventually could allow differential drug dosing. However, variability of MIC data for some M tuberculosis complex isolates carrying identical mutations also reveals limitations of our understanding of the genotype and phenotype relationships (eg, including epistasis and strain genetic background).

FUNDING

Bill & Melinda Gates Foundation, German Centre for Infection Research, German Research Foundation, Excellence Cluster Precision Medicine of Inflammation (EXC 2167), and Leibniz ScienceCampus EvoLUNG.

Molecular Epidemiology of Mycobacterium abscessus Isolates Recovered from German Cystic Fibrosis Patients

Infections due to Mycobacterium abscessus are a major cause of mortality and morbidity in cystic fibrosis (CF) patients. Furthermore, M. abscessus has been suspected to be involved in person-to-person transmissions. In 2016, dominant global clonal complexes (DCCs) that occur worldwide among CF patients have been described. To elucidate the epidemiological situation of M. abscessus among CF patients in Germany and to put these data into a global context, we performed whole-genome sequencing of a set of 154 M. abscessus isolates from 123 German patients treated in 14 CF centers. We used MTBseq pipeline to identify clusters of closely related isolates and correlate those with global findings. Genotypic drug susceptibility for macrolides and aminoglycosides was assessed by characterization of the erm(41), rrl, and rrs genes. By this approach, we could identify representatives of all major DCCs (Absc 1, Absc 2, and Mass 1) in our cohort. Intrapersonal isolates showed higher genetic relatedness than interpersonal isolates (median 3 SNPs versus 16 SNPs; P < 0.001). We further identified four clusters with German patients from same centers clustering with less than 25 SNPs distance (range 3 to 18 SNPs) but did not find any hint for in-hospital person-to-person transmission. This is the largest study investigating phylogenetic relations of M. abscessus isolates in Germany. We identified representatives of all reported DCCs but evidence for nosocomial transmission remained inconclusive. Thus, the occurrence of genetically closely related isolates of M. abscessus has to be interpreted with care, as a direct interhuman transmission cannot be directly deduced. IMPORTANCE Mycobacterium abscessus is a major respiratory pathogen in cystic fibrosis (CF) patients. Recently it has been shown that dominant global clonal complexes (DCCs) have spread worldwide among CF patients. This study investigated the epidemiological situation of M. abscessus among CF patients in Germany by performing whole-genome sequencing (WGS) of a set of 154 M. abscessus from 123 German patients treated in 14 CF centers. This is the largest study investigating the phylogenetic relationship of M. abscessus CF isolates in Germany.

Transcontinental spread and evolution of Mycobacterium tuberculosis W148 European/Russian clade toward extensively drug resistant tuberculosis

Transmission-driven multi-/extensively drug resistant (M/XDR) tuberculosis (TB) is the largest single contributor to human mortality due to antimicrobial resistance. A few major clades of the Mycobacterium tuberculosis complex belonging to lineage 2, responsible for high prevalence of MDR-TB in Eurasia, show outstanding transnational distributions. Here, we determined factors underlying the emergence and epidemic spread of the W148 clade by genome sequencing and Bayesian demogenetic analyses of 720 isolates from 23 countries. We dated a common ancestor around 1963 and identified two successive epidemic expansions in the late 1980s and late 1990s, coinciding with major socio-economic changes in the post-Soviet Era. These population expansions favored accumulation of resistance mutations to up to 11 anti-TB drugs, with MDR evolving toward additional resistances to fluoroquinolones and second-line injectable drugs within 20 years on average. Timescaled haplotypic density analysis revealed that widespread acquisition of compensatory mutations was associated with transmission success of XDR strains. Virtually all W148 strains harbored a hypervirulence-associated ppe38 gene locus, and incipient recurrent emergence of prpR mutation-mediated drug tolerance was detected. The outstanding genetic arsenal of this geographically widespread M/XDR strain clade represents a “perfect storm” that jeopardizes the successful introduction of new anti-M/XDR-TB antibiotic regimens.

An outbreak of the multidrug-resistant Mycobacterium tuberculosis lineage W148 has spread widely across Russia, Central Asia and Europe. Here, the authors use whole genome sequences of ~700 isolates of this lineage collected over ~20 years to analyze its spread, evolution of drug resistance, and impact of compensatory mutations.

High fluoroquinolone resistance proportions among multidrug-resistant tuberculosis driven by dominant L2 Mycobacterium tuberculosis clones in the Mumbai Metropolitan Region

Background

Multidrug-resistant (MDR) Mycobacterium tuberculosis complex (MTBC) strains are a serious health problem in India, also contributing to one-fourth of the global MDR tuberculosis (TB) burden. About 36% of the MDR MTBC strains are reported fluoroquinolone (FQ) resistant leading to high pre-extensively drug-resistant (pre-XDR) and XDR-TB (further resistance against bedaquiline and/or linezolid) rates. Still, factors driving the MDR/pre-XDR epidemic in India are not well defined.

Methods

In a retrospective study, we analyzed 1852 consecutive MTBC strains obtained from patients from a tertiary care hospital laboratory in Mumbai by whole genome sequencing (WGS). Univariate and multivariate statistics was used to investigate factors associated with pre-XDR. Core genome multi locus sequence typing, time scaled haplotypic density (THD) method and homoplasy analysis were used to analyze epidemiological success, and positive selection in different strain groups, respectively.

Results

In total, 1016 MTBC strains were MDR, out of which 703 (69.2%) were pre-XDR and 45 (4.4%) were XDR. Cluster rates were high among MDR (57.8%) and pre-XDR/XDR (79%) strains with three dominant L2 (Beijing) strain clusters (Cl 1–3) representing half of the pre-XDR and 40% of the XDR-TB cases. L2 strains were associated with pre-XDR/XDR-TB (P < 0.001) and, particularly Cl 1–3 strains, had high first-line and FQ resistance rates (81.6–90.6%). Epidemic success analysis using THD showed that L2 strains outperformed L1, L3, and L4 strains in short- and long-term time scales. More importantly, L2 MDR and MDR + strains had higher THD success indices than their not-MDR counterparts. Overall, compensatory mutation rates were highest in L2 strains and positive selection was detected in genes of L2 strains associated with drug tolerance (prpB and ppsA) and virulence (Rv2828c). Compensatory mutations in L2 strains were associated with a threefold increase of THD indices, suggesting improved transmissibility.

Conclusions

Our data indicate a drastic increase of FQ resistance, as well as emerging bedaquiline resistance which endangers the success of newly endorsed MDR-TB treatment regimens. Rapid changes in treatment and control strategies are required to contain transmission of highly successful pre-XDR L2 strains in the Mumbai Metropolitan region but presumably also India-wide.

Prediction of drug resistance by Sanger sequencing of Mycobacterium tuberculosis complex strains isolated from multidrug resistant tuberculosis suspect patients in Ethiopia

Background

Ethiopia is one of the high multidrug-resistant tuberculosis (MDR-TB) burden countries. However, phenotypic drug susceptibility testing can take several weeks due to the slow growth of Mycobacterium tuberculosis complex (MTBC) strains. In this study, we assessed the performance of a Sanger sequencing approach to predict resistance against five anti-tuberculosis drugs and the pattern of resistance mediating mutations.

Methods

We enrolled 226 MTBC culture-positive MDR-TB suspects and collected sputum specimens and socio-demographic and TB related data from each suspect between June 2015 and December 2016 in Addis Ababa, Ethiopia. Phenotypic drug susceptibility testing (pDST) for rifampicin, isoniazid, pyrazinamide, ethambutol, and streptomycin using BACTEC MGIT 960 was compared with the results of a Sanger sequencing analysis of seven resistance determining regions in the genes rpoB, katG, fabG-inhA, pncA, embB, rpsL, and rrs.

Result

DNA isolation for Sanger sequencing was successfully extracted from 92.5% (209/226) of the MTBC positive cultures, and the remaining 7.5% (17/226) strains were excluded from the final analysis. Based on pDST results, drug resistance proportions were as follows: isoniazid: 109/209 (52.2%), streptomycin: 93/209 (44.5%), rifampicin: 88/209 (42.1%), ethambutol: 74/209 (35.4%), and pyrazinamide: 69/209 (33.0%). Resistance against isoniazid was mainly mediated by the mutation katG S315T (97/209, 46.4%) and resistance against rifampicin by rpoB S531L (58/209, 27.8%). The dominating resistance-conferring mutations for ethambutol, streptomycin, and pyrazinamide affected codon 306 in embB (48/209, 21.1%), codon 88 in rpsL (43/209, 20.6%), and codon 65 in pncA (19/209, 9.1%), respectively. We observed a high agreement between phenotypic and genotypic DST, such as 89.9% (at 95% confidence interval [CI], 84.2%–95.8%) for isoniazid, 95.5% (95% CI, 91.2%–99.8%) for rifampicin, 98.6% (95% CI, 95.9–100%) for ethambutol, 91.3% (95% CI, 84.6–98.1%) for pyrazinamide and 57.0% (95% CI, 46.9%–67.1%) for streptomycin.

Conclusion

We detected canonical mutations implicated in resistance to rifampicin, isoniazid, pyrazinamide, ethambutol, and streptomycin. High agreement with phenotypic DST results for all drugs renders Sanger sequencing promising to be performed as a complementary measure to routine phenotypic DST in Ethiopia. Sanger sequencing directly from sputum may accelerate accurate clinical decision-making in the future.

One-year surveillance of SARS-CoV-2 transmission of the ELISA cohort: A model for population-based monitoring of infection risk

With newly rising coronavirus disease 2019 (COVID-19) cases, important data gaps remain on (i) long-term dynamics of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection rates in fixed cohorts (ii) identification of risk factors, and (iii) establishment of effective surveillance strategies. By polymerase chain reaction and antibody testing of 1% of the local population and >90,000 app-based datasets, the present study surveilled a catchment area of 300,000 inhabitants from March 2020 to February 2021. Cohort (56% female; mean age, 45.6 years) retention was 75 to 98%. Increased risk for seropositivity was detected in several high-exposure groups, especially nurses. Unreported infections dropped from 92 to 29% during the study. "Contact to COVID-19-affected" was the strongest risk factor, whereas public transportation, having children in school, or tourism did not affect infection rates. With the first SARS-CoV-2 cohort study, we provide a transferable model for effective surveillance, enabling monitoring of reinfection rates and increased preparedness for future pandemics.

Updating the approaches to define susceptibility and resistance to anti-tuberculosis agents: implications for diagnosis and treatment

Approximately 85 000 deaths globally in 2019 were due to drug-resistant tuberculosis (TB), which corresponds to 7% of global deaths attributable to bacterial antimicrobial resistance [1]. Yet concerns have been mounting that drug-resistant TB was being underestimated because the approaches to define susceptibility and resistance to anti-TB agents had not kept up with those used for other major bacterial pathogens [2–9]. Here, we outline the recent, evidence-based initiatives spearheaded by the World Health Organization (WHO) and others to update breakpoints (traditionally referred to as critical concentrations (CCs)) that are used for phenotypic antimicrobial susceptibility testing (AST), also called drug susceptibility testing in the TB literature.

Inappropriately high breakpoints have resulted in systematic false-susceptible AST results to anti-TB drugs. MIC, PK/PD and clinical outcome data should be combined when setting breakpoints to minimise the emergence and spread of antimicrobial resistance. https://bit.ly/3i43wb6

Sub-Lineage Specific Phenolic Glycolipid Patterns in the Mycobacterium tuberculosis Complex Lineage 1

“Ancestral” Mycobacterium tuberculosis complex (MTBC) strains of Lineage 1 (L1, East African Indian) are a prominent tuberculosis (TB) cause in countries around the Indian Ocean. However, the pathobiology of L1 strains is insufficiently characterized. Here, we used whole genome sequencing (WGS) of 312 L1 strains from 43 countries to perform a characterization of the global L1 population structure and correlate this to the analysis of the synthesis of phenolic glycolipids (PGL) – known MTBC polyketide-derived virulence factors. Our results reveal the presence of eight major L1 sub-lineages, whose members have specific mutation signatures in PGL biosynthesis genes, e.g., pks15/1 or glycosyltransferases Rv2962c and/or Rv2958c. Sub-lineage specific PGL production was studied by NMR-based lipid profiling and strains with a completely abolished phenolphthiocerol dimycoserosate biosynthesis showed in average a more prominent growth in human macrophages. In conclusion, our results show a diverse population structure of L1 strains that is associated with the presence of specific PGL types. This includes the occurrence of mycoside B in one sub-lineage, representing the first description of a PGL in an M. tuberculosis lineage other than L2. Such differences may be important for the evolution of L1 strains, e.g., allowing adaption to different human populations.

Emergence of bedaquiline resistance in a high tuberculosis burden country

RATIONALE

Bedaquiline has been classified as a group A drug for the treatment of multidrug-resistant tuberculosis (MDR-TB) by the World Health Organization; however, globally emerging resistance threatens the effectivity of novel MDR-TB treatment regimens.

OBJECTIVES

We analysed pre-existing and emerging bedaquiline resistance in bedaquiline-based MDR-TB therapies, and risk factors associated with treatment failure and death.

METHODS

In a cross-sectional cohort study, we employed patient data, whole-genome sequencing (WGS) and phenotyping of Mycobacterium tuberculosis complex (MTBC) isolates. We could retrieve baseline isolates from 30.5% (62 out of 203) of all MDR-TB patients who received bedaquiline between 2016 and 2018 in the Republic of Moldova. This includes 26 patients for whom we could also retrieve a follow-up isolate.

MEASUREMENTS AND MAIN RESULTS

At baseline, all MTBC isolates were susceptible to bedaquiline. Among 26 patients with available baseline and follow-up isolates, four (15.3%) patients harboured strains which acquired bedaquiline resistance under therapy, while one (3.8%) patient was re-infected with a second bedaquiline-resistant strain. Treatment failure and death were associated with cavitary disease (p=0.011), and any additional drug prescribed in the bedaquiline-containing regimen with WGS-predicted resistance at baseline (OR 1.92 per unit increase, 95% CI 1.15-3.21; p=0.012).

CONCLUSIONS

MDR-TB treatments based on bedaquiline require a functional background regimen to achieve high cure rates and to prevent the evolution of bedaquiline resistance. Novel MDR-TB therapies with bedaquiline require timely and comprehensive drug resistance monitoring.

Whole genome sequencing-based classification of human-related Haemophilus species and detection of antimicrobial resistance genes

Background

Bacteria belonging to the genus Haemophilus cause a wide range of diseases in humans. Recently, H. influenzae was classified by the WHO as priority pathogen due to the wide spread of ampicillin resistant strains. However, other Haemophilus spp. are often misclassified as H. influenzae. Therefore, we established an accurate and rapid whole genome sequencing (WGS) based classification and serotyping algorithm and combined it with the detection of resistance genes.

Methods

A gene presence/absence-based classification algorithm was developed, which employs the open-source gene-detection tool SRST2 and a new classification database comprising 36 genes, including capsule loci for serotyping. These genes were identified using a comparative genome analysis of 215 strains belonging to ten human-related Haemophilus (sub)species (training dataset). The algorithm was evaluated on 1329 public short read datasets (evaluation dataset) and used to reclassify 262 clinical Haemophilus spp. isolates from 250 patients (German cohort). In addition, the presence of antibiotic resistance genes within the German dataset was evaluated with SRST2 and correlated with results of traditional phenotyping assays.

Results

The newly developed algorithm can differentiate between clinically relevant Haemophilus species including, but not limited to, H. influenzae, H. haemolyticus, and H. parainfluenzae. It can also identify putative haemin-independent H. haemolyticus strains and determine the serotype of typeable Haemophilus strains. The algorithm performed excellently in the evaluation dataset (99.6% concordance with reported species classification and 99.5% with reported serotype) and revealed several misclassifications. Additionally, 83 out of 262 (31.7%) suspected H. influenzae strains from the German cohort were in fact H. haemolyticus strains, some of which associated with mouth abscesses and lower respiratory tract infections.

Resistance genes were detected in 16 out of 262 datasets from the German cohort. Prediction of ampicillin resistance, associated with bla_TEM-1D, and tetracycline resistance, associated with tetB, correlated well with available phenotypic data.

Conclusions

Our new classification database and algorithm have the potential to improve diagnosis and surveillance of Haemophilus spp. and can easily be coupled with other public genotyping and antimicrobial resistance databases. Our data also point towards a possible pathogenic role of H. haemolyticus strains, which needs to be further investigated.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01017-x.