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.

Rapid molecular diagnostics of tuberculosis resistance by targeted stool sequencing

BACKGROUND

Stool is an important diagnostic specimen for tuberculosis in populations who struggle to provide sputum, such as children or people living with HIV. However, the culture of Mycobacterium tuberculosis (M. tuberculosis) complex strains from stool perform poorly. This limits the opportunity for phenotypic drug resistance testing with this specimen. Therefore, reliable molecular methods are urgently needed for comprehensive drug resistance testing on stool specimens.

METHODS

We evaluated the performance of targeted next-generation sequencing (tNGS, Deeplex® Myc-TB) for the detection of mutations associated with M. tuberculosis complex drug resistance on DNA isolated from stool specimens provided by participants from a prospective cohort of patients treated for tuberculosis in Eswatini (n = 66; 56 with and 10 participants without M. tuberculosis complex DNA detected in stool by real-time quantitative PCR), and an independent German validation cohort of participants with culture-confirmed tuberculosis (n = 21).

RESULTS

The tNGS assay detected M. tuberculosis complex DNA in 38 of 56 (68%) samples; for 28 of 38 (74%) samples, a full M. tuberculosis complex drug resistance prediction report was obtained. There was a high degree of concordance with sputum phenotypic drug susceptibility results (κ = 0.82). The ability to predict resistance was concentration-dependent and successful in 7/10 (70%), 18/25 (72%), and 3/21 (14%) of samples with stool PCR concentration thresholds of > 100 femtogram per microliter (fg/μl), 1 to 100 fg/μl, and < 1 fg/μl, respectively (p = 0.0004). The German cohort confirmed these results and demonstrated a similarly high concordance between stool tNGS and sputum phenotypic drug susceptibility results (κ = 0.84).

CONCLUSIONS

tNGS can identify drug resistance from stool provided by tuberculosis patients. This affords the opportunity to obtain critical diagnostic information for tuberculosis patients who struggle to provide respiratory specimens.

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.

OUP accepted manuscript

Background

Rifampicin- or multidrug-resistant (RR/MDR) Mycobacterium tuberculosis complex (MTBC) strains account for considerable morbidity and mortality globally. WGS-based prediction of drug resistance may guide clinical decisions, especially for the design of RR/MDR-TB therapies.

Methods

We compared WGS-based drug resistance-predictive mutations for 42 MTBC isolates from MDR-TB patients in Tanzania with the MICs of 14 antibiotics measured in the Sensititre™ MycoTB assay. An isolate was phenotypically categorized as resistant if it had an MIC above the epidemiological-cut-off (ECOFF) value, or as susceptible if it had an MIC below or equal to the ECOFF.

Results

Overall, genotypically non-wild-type MTBC isolates with high-level resistance mutations (gNWT-R) correlated with isolates with MIC values above the ECOFF. For instance, the median MIC value (mg/L) for rifampicin-gNWT-R strains was >4.0 (IQR 4.0–4.0) compared with 0.5 (IQR 0.38–0.50) in genotypically wild-type (gWT-S, P < 0.001); isoniazid-gNWT-R >4.0 (IQR 2.0–4.0) compared with 0.25 (IQR 0.12–1.00) among gWT-S (P = 0.001); ethionamide-gNWT-R 15.0 (IQR 10.0–20.0) compared with 2.50 (IQR; 2.50–5.00) among gWT-S (P < 0.001). WGS correctly predicted resistance in 95% (36/38) and 100% (38/38) of the rifampicin-resistant isolates with ECOFFs >0.5 and >0.125 mg/L, respectively. No known resistance-conferring mutations were present in genes associated with resistance to fluoroquinolones, aminoglycosides, capreomycin, bedaquiline, delamanid, linezolid, clofazimine, cycloserine, or p-amino salicylic acid.

Conclusions

WGS-based drug resistance prediction worked well to rule-in phenotypic drug resistance and the absence of second-line drug resistance-mediating mutations has the potential to guide the design of RR/MDR-TB regimens in the future.

OUP accepted manuscript

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.

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

Whole-genome sequencing (WGS) of Mycobacterium tuberculosis (MTB) isolates can be used to get an accurate diagnosis, to guide clinical decision making, to control tuberculosis (TB) and for outbreak investigations. We evaluated the performance of long-read (LR) and/or short-read (SR) sequencing for anti-TB drug-resistance prediction using the TBProfiler and Mykrobe tools, the fraction of genome recovery, assembly accuracies and the robustness of two typing approaches based on core-genome SNP (cgSNP) typing and core-genome multi-locus sequence typing (cgMLST). Most of the discrepancies between phenotypic drug-susceptibility testing (DST) and drug-resistance prediction were observed for the first-line drugs rifampicin, isoniazid, pyrazinamide and ethambutol, mainly with LR sequence data. Resistance prediction to second-line drugs made by both TBProfiler and Mykrobe tools with SR- and LR-sequence data were in complete agreement with phenotypic DST except for one isolate. The SR assemblies were more accurate than the LR assemblies, having significantly (P <0.05) fewer indels and mismatches per 100 kbp. However, the hybrid and LR assemblies had slightly higher genome fractions. For LR assemblies, Canu followed by Racon, and Medaka polishing was the most accurate approach. The cgSNP approach, based on either reads or assemblies, was more robust than the cgMLST approach, especially for LR sequence data. In conclusion, anti-TB drug-resistance prediction, particularly with only LR sequence data, remains challenging, especially for first-line drugs. In addition, SR assemblies appear more accurate than LR ones, and reproducible phylogeny can be achieved using cgSNP approaches.

Population structure, biogeography and transmissibility of Mycobacterium tuberculosis

Mycobacterium tuberculosis is a clonal pathogen proposed to have co-evolved with its human host for millennia, yet our understanding of its genomic diversity and biogeography remains incomplete. Here we use a combination of phylogenetics and dimensionality reduction to reevaluate the population structure of M. tuberculosis, providing an in-depth analysis of the ancient Indo-Oceanic Lineage 1 and the modern Central Asian Lineage 3, and expanding our understanding of Lineages 2 and 4. We assess sub-lineages using genomic sequences from 4939 pan-susceptible strains, and find 30 new genetically distinct clades that we validate in a dataset of 4645 independent isolates. We find a consistent geographically restricted or unrestricted pattern for 20 groups, including three groups of Lineage 1. The distribution of terminal branch lengths across the M. tuberculosis phylogeny supports the hypothesis of a higher transmissibility of Lineages 2 and 4, in comparison with Lineages 3 and 1, on a global scale. We define an expanded barcode of 95 single nucleotide substitutions that allows rapid identification of 69 M. tuberculosis sub-lineages and 26 additional internal groups. Our results paint a higher resolution picture of the M. tuberculosis phylogeny and biogeography.

Role of Epistasis in Amikacin, Kanamycin, Bedaquiline, and Clofazimine Resistance in Mycobacterium tuberculosis Complex

Antibiotic resistance among bacterial pathogens poses a major global health threat. Mycobacterium tuberculosis complex (MTBC) is estimated to have the highest resistance rates of any pathogen globally. Given the low growth rate and the need for a biosafety level 3 laboratory, the only realistic avenue to scale up drug susceptibility testing (DST) for this pathogen is to rely on genotypic techniques. This raises the fundamental question of whether a mutation is a reliable surrogate for phenotypic resistance or whether the presence of a second mutation can completely counteract its effect, resulting in major diagnostic errors (i.e., systematic false resistance results). To date, such epistatic interactions have only been reported for streptomycin that is now rarely used. By analyzing more than 31,000 MTBC genomes, we demonstrated that the eis C-14T promoter mutation, which is interrogated by several genotypic DST assays endorsed by the World Health Organization, cannot confer resistance to amikacin and kanamycin if it coincides with loss-of-function (LoF) mutations in the coding region of eis. To our knowledge, this represents the first definitive example of antibiotic reversion in MTBC. Moreover, we raise the possibility that mmpR (Rv0678) mutations are not valid markers of resistance to bedaquiline and clofazimine if these coincide with an LoF mutation in the efflux pump encoded by mmpS5 (Rv0677c) and mmpL5 (Rv0676c).

Transmission patterns of rifampicin resistant Mycobacterium tuberculosis complex strains in Cameroon: a genomic epidemiological study

Background

Determining factors affecting the transmission of rifampicin (RR) and multidrug-resistant (MDR) Mycobacterium tuberculosis complex strains under standardized tuberculosis (TB) treatment is key to control TB and prevent the evolution of drug resistance.

Methods

We combined bacterial whole genome sequencing (WGS) and epidemiological investigations for 37% (n = 195) of all RR/MDR-TB patients in Cameroon (2012–2015) to identify factors associated with recent transmission.

Results

Patients infected with a strain resistant to high-dose isoniazid, and ethambutol had 7.4 (95% CI 2.6–21.4), and 2.4 (95% CI 1.2–4.8) times increased odds of being in a WGS-cluster, a surrogate for recent transmission. Furthermore, age between 30 and 50 was positively correlated with recent transmission (adjusted OR 3.8, 95% CI 1.3–11.4). We found high drug-resistance proportions against three drugs used in the short standardized MDR-TB regimen in Cameroon, i.e. high-dose isoniazid (77.4%), ethambutol (56.9%), and pyrazinamide (43.1%). Virtually all strains were susceptible to fluoroquinolones, kanamycin, and clofazimine, and treatment outcomes were mostly favourable (87.5%).

Conclusion

Pre-existing resistance to high-dose isoniazid, and ethambutol is associated with recent transmission of RR/MDR strains in our study. A possible contributing factor for this observation is the absence of universal drug susceptibility testing in Cameroon, likely resulting in prolonged exposure of new RR/MDR-TB patients to sub-optimal or failing first-line drug regimens.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-021-06593-8.

WNT6/ACC2-induced storage of triacylglycerols in macrophages is exploited by Mycobacterium tuberculosis

In view of emerging drug-resistant tuberculosis (TB), host-directed adjunct therapies are urgently needed to improve treatment outcomes with currently available anti-TB therapies. One approach is to interfere with the formation of lipid-laden "foamy" macrophages in the host, as they provide a nutrient-rich host cell environment for Mycobacterium tuberculosis (Mtb). Here, we provide evidence that Wnt family member 6 (WNT6), a ligand of the evolutionarily conserved Wingless/Integrase 1 (WNT) signaling pathway, promotes foam cell formation by regulating key lipid metabolic genes including acetyl-CoA carboxylase 2 (ACC2) during pulmonary TB. Using genetic and pharmacological approaches, we demonstrated that lack of functional WNT6 or ACC2 significantly reduced intracellular triacylglycerol (TAG) levels and Mtb survival in macrophages. Moreover, treatment of Mtb-infected mice with a combination of a pharmacological ACC2 inhibitor and the anti-TB drug isoniazid (INH) reduced lung TAG and cytokine levels, as well as lung weights, compared with treatment with INH alone. This combination also reduced Mtb bacterial numbers and the size of mononuclear cell infiltrates in livers of infected mice. In summary, our findings demonstrate that Mtb exploits WNT6/ACC2-induced storage of TAGs in macrophages to facilitate its intracellular survival, a finding that opens new perspectives for host-directed adjunctive treatment of pulmonary TB.

High potency of sequential therapy with only β-lactam antibiotics

Evolutionary adaptation is a major source of antibiotic resistance in bacterial pathogens. Evolution-informed therapy aims to constrain resistance by accounting for bacterial evolvability. Sequential treatments with antibiotics that target different bacterial processes were previously shown to limit adaptation through genetic resistance trade-offs and negative hysteresis. Treatment with homogeneous sets of antibiotics is generally viewed to be disadvantageous as it should rapidly lead to cross-resistance. We here challenged this assumption by determining the evolutionary response of Pseudomonas aeruginosa to experimental sequential treatments involving both heterogenous and homogeneous antibiotic sets. To our surprise, we found that fast switching between only β-lactam antibiotics resulted in increased extinction of bacterial populations. We demonstrate that extinction is favored by low rates of spontaneous resistance emergence and low levels of spontaneous cross-resistance among the antibiotics in sequence. The uncovered principles may help to guide the optimized use of available antibiotics in highly potent, evolution-informed treatment designs.

Overuse of antibiotic drugs is leading to the appearance of antibiotic-resistant bacteria; this is, bacteria with mutations that allow them to survive treatment with specific antibiotics. This has made some bacterial infections difficult or impossible to treat. Learning more about how bacteria evolve resistance to antibiotics could help scientists find ways to prevent it and develop more effective treatments.

Changing antibiotics frequently may be one way to prevent bacteria from evolving resistance. That way if a bacterium acquires mutations that allow it to escape one antibiotic, another antibiotic will kill it, stopping it from dividing and preventing the appearance of descendants with resistance to several antibiotics. In order to use this approach, testing is needed to find the best sequences of antibiotics to apply and the optimal timings of treatment.

To find out more, Batra, Roemhild et al. grew bacteria in the laboratory and exposed them to different sequences of antibiotics, switching antibiotics at different time intervals. This showed that sequential treatments with different antibiotics can limit bacterial evolution, especially when antibiotics are switched quickly. Unexpectedly, one of the most effective sequences used very similar antibiotics. This was surprising because using similar antibiotics should lead to the evolution of cross-resistance, which is when a drug causes changes that make the bacterium less sensitive to other treatments. However, in the tested case, cross-resistance did not evolve when antibiotics were switched quickly, thereby ensuring efficiency of treatment.

Batra et al. show that alternating sequences of antibiotics may be an effective strategy to prevent drug resistance. Because the experiments were done in a laboratory setting it will be important to verify the results in studies in animals and humans before the approach can be used in medical or veterinary settings. If the results are confirmed, it could reduce the need to develop new antibiotics, which is expensive and time consuming.

Implementation of whole genome sequencing for tuberculosis diagnostics in a low-middle income, high MDR-TB burden country

Whole genome sequencing (WGS) is revolutionary for diagnostics of TB and its mutations associated with drug-resistances, but its uptake in low- and middle-income countries is hindered by concerns of implementation feasibility. Here, we provide a proof of concept for its successful implementation in such a setting. WGS was implemented in the Kyrgyz Republic. We estimated needs of up to 55 TB-WGS per week and chose the MiSeq platform (Illumina, USA) because of its capacity of up to 60 TB-WGS per week. The project's timeline was completed in 93-weeks. Costs of large equipment and accompanying costs were 222,065 USD and 8462 USD, respectively. The first 174 WGS costed 277 USD per sequence, but this was skewed by training inefficiencies. Based on real prices and presuming optimal utilization of WGS capacities, WGS costs could drop to 167 and 141 USD per WGS using MiSeq Reagent Kits v2 (500-cycles) and v3 (600-cycles), respectively. Five trainings were required to prepare the staff for autonomous WGS which cost 48,250 USD. External assessment confirmed excellent performance of WGS by the Kyrgyz laboratory in an interlaboratory comparison of 30 M. tuberculosis genomes showing complete agreeance of results.

Mycobacterium tuberculosis complex lineage 5 exhibits high levels of within-lineage genomic diversity and differing gene content compared to the type strain H37Rv

Pathogens of the Mycobacterium tuberculosis complex (MTBC) are considered to be monomorphic, with little gene content variation between strains. Nevertheless, several genotypic and phenotypic factors separate strains of the different MTBC lineages (L), especially L5 and L6 (traditionally termed Mycobacterium africanum) strains, from each other. However, this genome variability and gene content, especially of L5 strains, has not been fully explored and may be important for pathobiology and current approaches for genomic analysis of MTBC strains, including transmission studies. By comparing the genomes of 355 L5 clinical strains (including 3 complete genomes and 352 Illumina whole-genome sequenced isolates) to each other and to H37Rv, we identified multiple genes that were differentially present or absent between H37Rv and L5 strains. Additionally, considerable gene content variability was found across L5 strains, including a split in the L5.3 sub-lineage into L5.3.1 and L5.3.2. These gene content differences had a small knock-on effect on transmission cluster estimation, with clustering rates influenced by the selected reference genome, and with potential overestimation of recent transmission when using H37Rv as the reference genome. We conclude that full capture of the gene diversity, especially high-resolution outbreak analysis, requires a variation of the single H37Rv-centric reference genome mapping approach currently used in most whole-genome sequencing data analysis pipelines. Moreover, the high within-lineage gene content variability suggests that the pan-genome of M. tuberculosis is at least several kilobases larger than previously thought, implying that a concatenated or reference-free genome assembly (de novo) approach may be needed for particular questions.

Microevolution of Mycobacterium tuberculosis Subpopulations and Heteroresistance in a Patient Receiving 27 Years of Tuberculosis Treatment in Germany

Preexisting and newly emerging resistant pathogen subpopulations (heteroresistance) are potential risk factors for treatment failure of multi/extensively drug resistant (MDR/XDR) tuberculosis (TB). Intrapatient evolutionary dynamics of Mycobacterium tuberculosis complex (Mtbc) strains and their implications on treatment outcomes are still not completely understood. To elucidate how Mtbc strains escape therapy, we analyzed 13 serial isolates from a German patient by whole-genome sequencing. Sequencing data were compared with phenotypic drug susceptibility profiles and the patient's collective 27-year treatment history to further elucidate factors fostering intrapatient resistance evolution. The patient endured five distinct TB episodes, ending in resistance to 16 drugs and a nearly untreatable XDR-TB infection. The first isolate obtained, during the patient's 5th TB episode, presented fixed resistance mutations to 7 anti-TB drugs, including isoniazid, rifampin, streptomycin, pyrazinamide, prothionamide, para-aminosalicylic acid, and cycloserine-terizidone. Over the next 13 years, a dynamic evolution with coexisting, heterogeneous subpopulations was observed in 6 out of 13 sequential bacterial isolates. The emergence of drug-resistant subpopulations coincided with frequent changes in treatment regimens, which often included two or fewer active compounds. This evolutionary arms race between competing subpopulations ultimately resulted in the fixation of a single XDR variant. Our data demonstrate the complex intrapatient microevolution of Mtbc subpopulations during failing MDR/XDR-TB treatment. Designing effective treatment regimens based on rapid detection of (hetero) resistance is key to avoid resistance development and treatment failure.

Improving tuberculosis surveillance by detecting international transmission using publicly available whole genome sequencing data

IntroductionImproving the surveillance of tuberculosis (TB) is especially important for multidrug-resistant (MDR) and extensively drug-resistant (XDR) TB. The large amount of publicly available whole genome sequencing (WGS) data for TB gives us the chance to re-use data and to perform additional analyses at a large scale.AimWe assessed the usefulness of raw WGS data of global MDR/XDR Mycobacterium tuberculosis isolates available from public repositories to improve TB surveillance.MethodsWe extracted raw WGS data and the related metadata of M. tuberculosis isolates available from the Sequence Read Archive. We compared this public dataset with WGS data and metadata of 131 MDR- and XDR M. tuberculosis isolates from Germany in 2012 and 2013.ResultsWe aggregated a dataset that included 1,081 MDR and 250 XDR isolates among which we identified 133 molecular clusters. In 16 clusters, the isolates were from at least two different countries. For example, Cluster 2 included 56 MDR/XDR isolates from Moldova, Georgia and Germany. When comparing the WGS data from Germany with the public dataset, we found that 11 clusters contained at least one isolate from Germany and at least one isolate from another country. We could, therefore, connect TB cases despite missing epidemiological information.ConclusionWe demonstrated the added value of using WGS raw data from public repositories to contribute to TB surveillance. Comparing the German with the public dataset, we identified potential international transmission events. Thus, using this approach might support the interpretation of national surveillance results in an international context.