Historic methicillin-resistant Staphylococcus aureus: expanding current knowledge using molecular epidemiological characterization of a Swiss legacy collection.

BACKGROUND: Few methicillin-resistant Staphylococcus aureus (MRSA) from the early years of its global emergence have been sequenced. Knowledge about evolutionary factors promoting the success of specific MRSA multi-locus sequence types (MLSTs) remains scarce. We aimed to characterize a legacy MRSA collection isolated from 1965 to 1987 and compare it against publicly available international and local genomes. METHODS: We accessed 451 historic (1965-1987) MRSA isolates stored in the Culture Collection of Switzerland, mostly collected from the Zurich region. We determined phenotypic antimicrobial resistance (AMR) and performed whole genome sequencing (WGS) using Illumina short-read sequencing on all isolates and long-read sequencing on a selection with Oxford Nanopore Technology. For context, we included 103 publicly available international assemblies from 1960 to 1992 and sequenced 1207 modern Swiss MRSA isolates from 2007 to 2022. We analyzed the core genome (cg)MLST and predicted SCCmec cassette types, AMR, and virulence genes. RESULTS: Among the 451 historic Swiss MRSA isolates, we found 17 sequence types (STs) of which 11 have been previously described. Two STs were novel combinations of known loci and six isolates carried previously unsubmitted MLST alleles, representing five new STs (ST7843, ST7844, ST7837, ST7839, and ST7842). Most isolates (83% 376/451) represented ST247-MRSA-I isolated in the 1960s, followed by ST7844 (6% 25/451), a novel single locus variant (SLV) of ST239. Analysis by cgMLST indicated that isolates belonging to ST7844-MRSA-III cluster within the diversity of ST239-MRSA-III. Early MRSA were predominantly from clonal complex (CC)8. From 1980 to the end of the twentieth century, we observed that CC22 and CC5 as well as CC8 were present, both locally and internationally. CONCLUSIONS: The combined analysis of 1761 historic and contemporary MRSA isolates across more than 50 years uncovered novel STs and allowed us a glimpse into the lineage flux between Swiss-German and international MRSA across time.

First Case of Meningoencephalitis and Bacteraemia with Flavobacterium Lindanitolerans.

Purpose: Flavobacterium lindanitolerans is an environmental Gram-negative, non-spore-forming rod which is usually not considered to be a human pathogen. Isolation from human clinical samples has been described only once. We report the first case of meningoencephalitis and bacteraemia with Flavobacterium lindanitolerans. Case description: A 76-year-old female presented with fever, headache, alteration of mental status, marked meningism and dysarthria. A lumbar puncture demonstrated cerebrospinal fluid findings consistent with bacterial meningitis, and a broad-spectrum antibiotic therapy was initiated. Blood and cerebrospinal fluid cultures revealed a growth of Flavobacterium lindanitolerans. Based on antimicrobial susceptibilities testing, antibiotic treatment was changed to levofloxacin, resulting in a remission of the clinical symptoms after 21 days of treatment. Conclusion: Flavobacterium species are extremely rare human pathogens. However, some of them have been reported to cause opportunistic infections. We describe the first case of meningoencephalitis and bacteraemia caused by Flavobacterium lindanitolerans which was effectively treated with levofloxacin for 21 days. LEARNING POINTS: This is the first case of a meningoencephalitis and bacteraemia with Flavobacterium lindanitolerans.Levofloxacin can be given for 21 days as a treatment of F. lindanitolerans.

Quality of MALDI-TOF mass spectra in routine diagnostics: results from an international external quality assessment including 36 laboratories from 12 countries using 47 challenging bacterial strains.

OBJECTIVES: Matrix assisted laser desorption ionization-time of flight (MALDI-TOF) mass spectrometry (MS) is a widely used method for bacterial species identification. Incomplete databases and mass spectral quality (MSQ) still represent major challenges. Important proxies for MSQ are the number of detected marker masses, reproducibility, and measurement precision. We aimed to assess MSQs across diagnostic laboratories and the potential of simple workflow adaptations to improve it. METHODS: For baseline MSQ assessment, 47 diverse bacterial strains, which are challenging to identify by MALDI-TOF MS, were routinely measured in 36 laboratories from 12 countries, and well-defined MSQ features were used. After an intervention consisting of detailed reported feedback and instructions on how to acquire MALDI-TOF mass spectra, measurements were repeated and MSQs were compared. RESULTS: At baseline, we observed heterogeneous MSQ between the devices, considering the median number of marker masses detected (range = [2-25]), reproducibility between technical replicates (range = [55%-86%]), and measurement error (range = [147 parts per million (ppm)-588 ppm]). As a general trend, the spectral quality was improved after the intervention for devices, which yielded low MSQs in the baseline assessment as follows: for four out of five devices with a high measurement error, the measurement precision was improved (p-values <0.001, paired Wilcoxon test); for six out of ten devices, which detected a low number of marker masses, the number of detected marker masses increased (p-values <0.001, paired Wilcoxon test). DISCUSSION: We have identified simple workflow adaptations, which, to some extent, improve MSQ of poorly performing devices and should be considered by laboratories yielding a low MSQ. Improving MALDI-TOF MSQ in routine diagnostics is essential for increasing the resolution of bacterial identification by MALDI-TOF MS, which is dependent on the reproducible detection of marker masses. The heterogeneity identified in this external quality assessment (EQA) requires further study.

Determinants of SARS-CoV-2 transmission to guide vaccination strategy in an urban area.

Transmission chains within small urban areas (accommodating ∼30 per cent of the European population) greatly contribute to case burden and economic impact during the ongoing coronavirus pandemic and should be a focus for preventive measures to achieve containment. Here, at very high spatio-temporal resolution, we analysed determinants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission in a European urban area, Basel-City (Switzerland). We combined detailed epidemiological, intra-city mobility and socio-economic data sets with whole-genome sequencing during the first SARS-CoV-2 wave. For this, we succeeded in sequencing 44 per cent of all reported cases from Basel-City and performed phylogenetic clustering and compartmental modelling based on the dominating viral variant (B.1-C15324T; 60 per cent of cases) to identify drivers and patterns of transmission. Based on these results we simulated vaccination scenarios and corresponding healthcare system burden (intensive care unit (ICU) occupancy). Transmissions were driven by socio-economically weaker and highly mobile population groups with mostly cryptic transmissions which lacked genetic and identifiable epidemiological links. Amongst more senior population transmission was clustered. Simulated vaccination scenarios assuming 60-90 per cent transmission reduction and 70-90 per cent reduction of severe cases showed that prioritising mobile, socio-economically weaker populations for vaccination would effectively reduce case numbers. However, long-term ICU occupation would also be effectively reduced if senior population groups were prioritised, provided there were no changes in testing and prevention strategies. Reducing SARS-CoV-2 transmission through vaccination strongly depends on the efficacy of the deployed vaccine. A combined strategy of protecting risk groups by extensive testing coupled with vaccination of the drivers of transmission (i.e. highly mobile groups) would be most effective at reducing the spread of SARS-CoV-2 within an urban area.

Direct antimicrobial resistance prediction from clinical MALDI-TOF mass spectra using machine learning.

Early use of effective antimicrobial treatments is critical for the outcome of infections and the prevention of treatment resistance. Antimicrobial resistance testing enables the selection of optimal antibiotic treatments, but current culture-based techniques can take up to 72 hours to generate results. We have developed a novel machine learning approach to predict antimicrobial resistance directly from matrix-assisted laser desorption/ionization-time of flight (MALDI-TOF) mass spectra profiles of clinical isolates. We trained calibrated classifiers on a newly created publicly available database of mass spectra profiles from the clinically most relevant isolates with linked antimicrobial susceptibility phenotypes. This dataset combines more than 300,000 mass spectra with more than 750,000 antimicrobial resistance phenotypes from four medical institutions. Validation on a panel of clinically important pathogens, including Staphylococcus aureus, Escherichia coli and Klebsiella pneumoniae, resulting in areas under the receiver operating characteristic curve of 0.80, 0.74 and 0.74, respectively, demonstrated the potential of using machine learning to substantially accelerate antimicrobial resistance determination and change of clinical management. Furthermore, a retrospective clinical case study of 63 patients found that implementing this approach would have changed the clinical treatment in nine cases, which would have been beneficial in eight cases (89%). MALDI-TOF mass spectra-based machine learning may thus be an important new tool for treatment optimization and antibiotic stewardship.

Colistin resistance in Gram-negative bacteria analysed by five phenotypic assays and inference of the underlying genomic mechanisms.

BACKGROUND: Colistin is used against multi-drug resistant pathogens, yet resistance emerges through dissemination of plasmid-mediated genes (mcr) or chromosomal mutation of genes involved in lipopolysaccharide synthesis (i.e. mgrB, phoPQ, pmrCAB). Phenotypic susceptibility testing is challenging due to poor diffusion of colistin in agar media, leading to an underestimation of resistance. Performance of five phenotypic approaches was compared in the context of different molecular mechanisms of resistance. We evaluated Vitek 2® (bioMérieux, AST N242), Colistin MIC Test Strip (Liofilchem Diagnostici), UMIC (Biocentric), and Rapid Polymyxin™ NP test (ELITechGroup) against the standard broth microdilution (BMD) method. We used whole genome sequencing (WGS) to infer molecular resistance mechanisms. We analysed 97 Enterobacterales and non-fermenting bacterial isolates, largely clinical isolates collected up to 2018. Data was analysed by comparing susceptibility categories (susceptible or resistant) and minimal inhibitory concentrations (MIC). Susceptibility category concordance is the percentage of test results sharing the same category to BMD. MIC concordance was calculated similarly but considering ±1 MIC titre error range. We determined genomic diversity by core genome multi locus sequencing typing (cgMLST) and identified putative antimicrobial resistance genes using NCBI and CARD databases, and manual annotation. RESULTS: Of 97 isolates, 54 (56%) were resistant with standard BMD. Highest susceptibility category concordance was achieved by Rapid Polymyxin™ NP (98.8%) followed by UMIC (97.9%), Colistin E-test MIC strip (96.9%) and Vitek 2® (95.6%). Highest MIC concordance was achieved by UMIC (80.4%), followed by Vitek 2® (72.5%) and Colistin E-test MIC strip (62.9%). Among resistant isolates, 23/54 (43%) were intrinsically resistant to colistin, whereas 31/54 (57%) isolates had acquired colistin resistance. Of these, mcr-1 was detected in four isolates and mcr-2 in one isolate. Non-synonymous mutations in mgrB, phoQ, pmrA, pmrB, and pmrC genes were encountered in Klebsiella pneumoniae, Escherichia coli, and Acinetobacter bereziniae resistant isolates. Mutations found in mgrB and pmrB were only identified in isolates exhibiting MICs of ≥16 mg/L. CONCLUSIONS: The Rapid Polymyxin™ NP test showed highest categorical concordance and the UMIC test provided MIC values with high concordance to BMD. We found colistin resistance in diverse species occurred predominantly through spontaneous chromosomal mutation rather than plasmid-mediated resistance.

Characterising the epidemic spread of influenza A/H3N2 within a city through phylogenetics.

Infecting large portions of the global population, seasonal influenza is a major burden on societies around the globe. While the global source sink dynamics of the different seasonal influenza viruses have been studied intensively, its local spread remains less clear. In order to improve our understanding of how influenza is transmitted on a city scale, we collected an extremely densely sampled set of influenza sequences alongside patient metadata. To do so, we sequenced influenza viruses isolated from patients of two different hospitals, as well as private practitioners in Basel, Switzerland during the 2016/2017 influenza season. The genetic sequences reveal that repeated introductions into the city drove the influenza season. We then reconstruct how the effective reproduction number changed over the course of the season. While we did not find that transmission dynamics in Basel correlate with humidity or school closures, we did find some evidence that it may positively correlated with temperature. Alongside the genetic sequence data that allows us to see how individual cases are connected, we gathered patient information, such as the age or household status. Zooming into the local transmission outbreaks suggests that the elderly were to a large extent infected within their own transmission network. In the remaining transmission network, our analyses suggest that school-aged children likely play a more central role than pre-school aged children. These patterns will be valuable to plan interventions combating the spread of respiratory diseases within cities given that similar patterns are observed for other influenza seasons and cities.

Identification of influenza urban transmission patterns by geographical, epidemiological and whole genome sequencing data: protocol for an observational study.

INTRODUCTION: Urban transmission patterns of influenza viruses are complex and poorly understood, and multiple factors may play a critical role in modifying transmission. Whole genome sequencing (WGS) allows the description of patient-to-patient transmissions at highest resolution. The aim of this study is to explore urban transmission patterns of influenza viruses in high detail by combining geographical, epidemiological and immunological data with WGS data. METHODS AND ANALYSIS: The study is performed at the University Hospital Basel, University Children's Hospital Basel and a network of paediatricians and family doctors in the Canton of Basel-City, Switzerland. The retrospective study part includes an analysis of PCR-confirmed influenza cases from 2013 to 2018. The prospective study parts include (1) a household survey regarding influenza-like illness (ILI) and vaccination against influenza during the 2015/2016 season; (2) an analysis of influenza viruses collected during the 2016/2017 season using WGS-viral genomic sequences are compared with determine genetic relatedness and transmissions; and (3) measurement of influenza-specific antibody titres against all vaccinated and circulated strains during the 2016/2017 season from healthy individuals, allowing to monitor herd immunity across urban quarters. Survey data and PCR-confirmed cases are linked to data from the Statistics Office of the Canton Basel-City and visualised using geo-information system mapping. WGS data will be analysed in the context of patient epidemiological data using phylodynamic analyses, and the obtained herd immunity for each quarter. Profound knowledge on the key geographical, epidemiological and immunological factors influencing urban influenza transmission will help to develop effective counter measurements. ETHICS AND DISSEMINATION: The study is registered and approved by the regional ethics committee as an observational study (EKNZ project ID 2015-363 and 2016-01735). It is planned to present the results at conferences and publish the data in scientific journals. TRIAL REGISTRATION NUMBER: NCT03010007.

Phenotypic and Genomic Analyses of Burkholderia stabilis Clinical Contamination, Switzerland.

A recent hospital outbreak related to premoistened gloves used to wash patients exposed the difficulties of defining Burkholderia species in clinical settings. The outbreak strain displayed key B. stabilis phenotypes, including the inability to grow at 42°C; we used whole-genome sequencing to confirm the pathogen was B. stabilis. The outbreak strain genome comprises 3 chromosomes and a plasmid, sharing an average nucleotide identity of 98.4% with B. stabilis ATCC27515 BAA-67, but with 13% novel coding sequences. The genome lacks identifiable virulence factors and has no apparent increase in encoded antimicrobial drug resistance, few insertion sequences, and few pseudogenes, suggesting this outbreak was an opportunistic infection by an environmental strain not adapted to human pathogenicity. The diversity among outbreak isolates (22 from patients and 16 from washing gloves) is only 6 single-nucleotide polymorphisms, although the genome remains plastic, with large elements stochastically lost from outbreak isolates.

Air-conditioner cooling towers as complex reservoirs and continuous source of Legionella pneumophila infection evidenced by a genomic analysis study in 2017, Switzerland.

IntroductionWater supply and air-conditioner cooling towers (ACCT) are potential sources of Legionella pneumophila infection in people. During outbreaks, traditional typing methods cannot sufficiently segregate L. pneumophila strains to reliably trace back transmissions to these artificial water systems. Moreover, because multiple L. pneumophila strains may be present within these systems, methods to adequately distinguish strains are needed. Whole genome sequencing (WGS) and core genome multilocus sequence typing (cgMLST), with their higher resolution are helpful in this respect. In summer 2017, the health administration of the city of Basel detected an increase of L. pneumophila infections compared with previous months, signalling an outbreak.AimWe aimed to identify L. pneumophila strains populating suspected environmental sources of the outbreak, and to assess the relations between these strains and clinical outbreak strains.MethodsAn epidemiological and WGS-based microbiological investigation was performed, involving isolates from the local water supply and two ACCTs (n = 60), clinical outbreak and non-outbreak related isolates from 2017 (n = 8) and historic isolates from 2003-2016 (n = 26).ResultsIn both ACCTs, multiple strains were found. Phylogenetic analysis of the ACCT isolates showed a diversity of a few hundred allelic differences in cgMLST. Furthermore, two isolates from one ACCT showed no allelic differences to three clinical isolates from 2017. Five clinical isolates collected in the Basel area in the last decade were also identical in cgMLST to recent isolates from the two ACCTs.ConclusionCurrent outbreak-related and historic isolates were linked to ACCTs, which form a complex environmental habitat where strains are conserved over years.

Effect of mutation and genetic background on drug resistance in Mycobacterium tuberculosis.

Bacterial factors may contribute to the global emergence and spread of drug-resistant tuberculosis (TB). Only a few studies have reported on the interactions between different bacterial factors. We studied drug-resistant Mycobacterium tuberculosis isolates from a nationwide study conducted from 2000 to 2008 in Switzerland. We determined quantitative drug resistance levels of first-line drugs by using Bactec MGIT-960 and drug resistance genotypes by sequencing the hot-spot regions of the relevant genes. We determined recent transmission by molecular methods and collected clinical data. Overall, we analyzed 158 isolates that were resistant to isoniazid, rifampin, or ethambutol, 48 (30.4%) of which were multidrug resistant. Among 154 isoniazid-resistant strains, katG mutations were associated with high-level and inhA promoter mutations with low-level drug resistance. Only katG(S315T) (65.6% of all isoniazid-resistant strains) and inhA promoter -15C/T (22.7%) were found in molecular clusters. M. tuberculosis lineage 2 (includes Beijing genotype) was associated with any drug resistance (adjusted odds ratio [OR], 3.0; 95% confidence interval [CI], 1.7 to 5.6; P < 0.0001). Lineage 1 was associated with inhA promoter -15C/T mutations (OR, 6.4; 95% CI, 2.0 to 20.7; P = 0.002). We found that the genetic strain background influences the level of isoniazid resistance conveyed by particular mutations (interaction tests of drug resistance mutations across all lineages; P < 0.0001). In conclusion, M. tuberculosis drug resistance mutations were associated with various levels of drug resistance and transmission, and M. tuberculosis lineages were associated with particular drug resistance-conferring mutations and phenotypic drug resistance. Our study also supports a role for epistatic interactions between different drug resistance mutations and strain genetic backgrounds in M. tuberculosis drug resistance.

"Pseudo-Beijing": evidence for convergent evolution in the direct repeat region of Mycobacterium tuberculosis.

BACKGROUND: Mycobacterium tuberculosis has a global population structure consisting of six main phylogenetic lineages associated with specific geographic regions and human populations. One particular M. tuberculosis genotype known as "Beijing" has repeatedly been associated with drug resistance and has been emerging in some parts of the world. "Beijing" strains are traditionally defined based on a characteristic spoligotyping pattern. We used three alternative genotyping techniques to revisit the phylogenetic classification of M. tuberculosis complex (MTBC) strains exhibiting the typical "Beijing" spoligotyping pattern. METHODS AND FINDINGS: MTBC strains were obtained from an ongoing molecular epidemiological study in Switzerland and Nepal. MTBC genotyping was performed based on SNPs, genomic deletions, and 24-loci MIRU-VNTR. We identified three MTBC strains from patients originating from Tibet, Portugal and Nepal which exhibited a spoligotyping patterns identical to the classical Beijing signature. However, based on three alternative molecular markers, these strains were assigned to Lineage 3 (also known as Delhi/CAS) rather than to Lineage 2 (also known as East-Asian lineage). Sequencing of the RD207 in one of these strains showed that the deletion responsible for this "Pseudo-Beijing" spoligotype was about 1,000 base pairs smaller than the usual deletion of RD207 in classical "Beijing" strains, which is consistent with an evolutionarily independent deletion event in the direct repeat (DR) region of MTBC. CONCLUSIONS: We provide an example of convergent evolution in the DR locus of MTBC, and highlight the limitation of using spoligotypes for strain classification. Our results indicate that a proportion of "Beijing" strains may have been misclassified in the past. Markers that are more phylogenetically robust should be used when exploring strain-specific differences in experimental or clinical phenotypes.