Novel Organism Verification and Analysis (NOVA) study: identification of 35 clinical isolates representing potentially novel bacterial taxa using a pipeline based on whole genome sequencing.

BACKGROUND: Reliable species identification of cultured isolates is essential in clinical bacteriology. We established a new study algorithm named NOVA - Novel Organism Verification and Analysis to systematically analyze bacterial isolates that cannot be characterized by conventional identification procedures MALDI-TOF MS and partial 16 S rRNA gene sequencing using Whole Genome Sequencing (WGS). RESULTS: We identified a total of 35 bacterial strains that represent potentially novel species. Corynebacterium sp. (n = 6) and Schaalia sp. (n = 5) were the predominant genera. Two strains each were identified within the genera Anaerococcus, Clostridium, Desulfovibrio, and Peptoniphilus, and one new species was detected within Citrobacter, Dermabacter, Helcococcus, Lancefieldella, Neisseria, Ochrobactrum (Brucella), Paenibacillus, Pantoea, Porphyromonas, Pseudoclavibacter, Pseudomonas, Psychrobacter, Pusillimonas, Rothia, Sneathia, and Tessaracoccus. Twenty-seven of 35 strains were isolated from deep tissue specimens or blood cultures. Seven out of 35 isolated strains identified were clinically relevant. In addition, 26 bacterial strains that could only be identified at the species level using WGS analysis, were mainly organisms that have been identified/classified very recently. CONCLUSION: Our new algorithm proved to be a powerful tool for detection and identification of novel bacterial organisms. Publicly available clinical and genomic data may help to better understand their clinical and ecological role. Our identification of 35 novel strains, 7 of which appear to be clinically relevant, shows the wide range of undescribed pathogens yet to define.

External Quality Assessment of SARS-CoV-2 Sequencing: an ESGMD-SSM Pilot Trial across 15 European Laboratories.

This first pilot trial on external quality assessment (EQA) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) whole-genome sequencing, initiated by the European Society of Clinical Microbiology and Infectious Diseases (ESCMID) Study Group for Genomic and Molecular Diagnostics (ESGMD) and the Swiss Society for Microbiology (SSM), aims to build a framework between laboratories in order to improve pathogen surveillance sequencing. Ten samples with various viral loads were sent out to 15 clinical laboratories that had free choice of sequencing methods and bioinformatic analyses. The key aspects on which the individual centers were compared were the identification of (i) single nucleotide polymorphisms (SNPs) and indels, (ii) Pango lineages, and (iii) clusters between samples. The participating laboratories used a wide array of methods and analysis pipelines. Most were able to generate whole genomes for all samples. Genomes were sequenced to various depths (up to a 100-fold difference across centers). There was a very good consensus regarding the majority of reporting criteria, but there were a few discrepancies in lineage and cluster assignments. Additionally, there were inconsistencies in variant calling. The main reasons for discrepancies were missing data, bioinformatic choices, and interpretation of data. The pilot EQA was overall a success. It was able to show the high quality of participating laboratories and provide valuable feedback in cases where problems occurred, thereby improving the sequencing setup of laboratories. A larger follow-up EQA should, however, improve on defining the variables and format of the report. Additionally, contamination and/or minority variants should be a further aspect of assessment.

Description of Pseudoclavibacter triregionum sp. nov. from human blood and Pseudoclavibacter albus comb. nov., and revised classification of the genus Pseudoclavibacter: proposal of Caespitibacter gen. nov., with Caespitibacter soli comb. nov. and Caespitibacter caeni comb. nov.

We present polyphasic taxonomic data to demonstrate that strain 125703-2019T, a human blood isolate, represents a novel species within the genus Pseudoclavibacter, and to reclassify the illegitimate Zimmermannella alba Lin et al., 2004 as Pseudoclavibacter albus comb. nov. Upon primary isolation, strain 125703-2019T could not be identified reliably using MALDI-TOF mass spectrometry during routine diagnostic work, but partial 16S rRNA gene sequence analysis revealed that it belonged to the genus Pseudoclavibacter. Average nucleotide identity and digital DNA-DNA hybridisation analyses confirmed that it represented a novel species within this genus. A detailed physiological characterisation yielded differential tests between the novel species and its nearest neighbor taxa, which could also be differentiated using MALDI-TOF mass spectrometry. We propose to formally classify this strain into the novel species Pseudoclavibacter triregionum sp. nov., with strain 125703-2019T (= R-76471T, LMG 31777T, CCUG 74796T) as the type strain. The whole-genome assembly of strain 125703-2019T has a size of 2.4 Mb and a G + C content of 72.74%. A Pseudoclavibacter pangenome analysis revealed that 667 gene clusters were exclusively present in strain 125703-2019T. While these gene clusters were enriched in several COG functional categories, this analysis did not reveal functions that explained the occurrence of this species in human infection. Finally, several phylogenetic and phylogenomic analyses demonstrated that the genus Pseudoclavibacter is polyphyletic with Pseudoclavibacter soli and Pseudoclavibacter caeni representing a unique and deeply branching line of descent within the family Microbacteriaceae. We therefore also propose to reclassify both species into the novel genus Caespitibacter gen. nov. as Caespitibacter soli comb. nov. and Caespitibacter caeni comb. nov., respectively, and with C. soli comb. nov. as the type species.

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.

Gulosibacter hominis sp. nov.: a novel human microbiome bacterium that may cause opportunistic infections.

We present genomic, phylogenomic, and phenotypic taxonomic data to demonstrate that three human ear isolates represent a novel species within the genus Gulosibacter. These isolates could not be identified reliably using MALDI-TOF mass spectrometry during routine diagnostic work, but partial 16S rRNA gene sequence analysis revealed that they belonged to the genus Gulosibacter. Overall genomic relatedness indices between the draft genome sequences of the three isolates and of the type strains of established Gulosibacter species confirmed that the three isolates represented a single novel Gulosibacter species. A biochemical characterisation yielded differential tests between the novel and established Gulosibacter species, which could also be differentiated using MALDI-TOF mass spectrometry. We propose to formally classify these three isolates into Gulosibacter hominis sp. nov., with 401352-2018 T (= LMG 31778 T, CCUG 74795 T) as the type strain. The whole-genome sequence of strain 401352-2018 T has a size of 2,340,181 bp and a G+C content of 62.04 mol%. A Gulosibacter pangenome analysis revealed 467 gene clusters that were exclusively present in G. hominis genomes. While these G. hominis specific gene clusters were enriched in several COG functional categories, this analysis did not reveal functions that suggested a role in the human microbiome, nor did it explain the occurrence of G. hominis in ear infections. The absence of acquired antimicrobial resistance determinants and virulence factors in the G. hominis genomes, and an analysis of publicly available 16S rRNA gene sequences and 16S rRNA amplicon sequencing data sets suggested that G. hominis is a member of the human skin microbiota that may occasionally be involved in opportunistic infections.

A Novel Lineage of Ceftriaxone-resistant Salmonella Typhi From India That Is Closely Related to XDR S. Typhi Found in Pakistan.

Two MDR Salmonella Typhi isolates from India were found by whole genome sequencing to be closely related to the 2016 XDR S. Typhi outbreak strain from Pakistan. The Indian isolates have no chromosomal antimicrobial resistance cassette but carry the IncY plasmid p60006. Both isolates are susceptible to chloramphenicol, azithromycin, and carbapenems.

Tick-Borne Relapsing Fever Caused by Borrelia persica in Traveler to Central Asia, 2019.

We report a case of tick-borne relapsing fever caused by Borrelia persica in a traveler returning to Switzerland from central Asia. After the disease was diagnosed by blood smear microscopy, the causative Borrelia species was confirmed by shotgun metagenomics sequencing. PCR and sequencing techniques provide highly sensitive diagnostic tools superior to microscopy.

Comprehensive global genome dynamics of Chlamydia trachomatis show ancient diversification followed by contemporary mixing and recent lineage expansion.

Chlamydia trachomatis is the world's most prevalent bacterial sexually transmitted infection and leading infectious cause of blindness, yet it is one of the least understood human pathogens, in part due to the difficulties of in vitro culturing and the lack of available tools for genetic manipulation. Genome sequencing has reinvigorated this field, shedding light on the contemporary history of this pathogen. Here, we analyze 563 full genomes, 455 of which are novel, to show that the history of the species comprises two phases, and conclude that the currently circulating lineages are the result of evolution in different genomic ecotypes. Temporal analysis indicates these lineages have recently expanded in the space of thousands of years, rather than the millions of years as previously thought, a finding that dramatically changes our understanding of this pathogen's history. Finally, at a time when almost every pathogen is becoming increasingly resistant to antimicrobials, we show that there is no evidence of circulating genomic resistance in C. trachomatis.

Prevalence and phylogeny of Chlamydiae and hemotropic mycoplasma species in captive and free-living bats.

BACKGROUND: Bats are hosts for a variety of microorganisms, however, little is known about the presence of Chlamydiales and hemotropic mycoplasmas. This study investigated 475 captive and free-living bats from Switzerland, Germany, and Costa Rica for Chlamydiales and hemotropic mycoplasmas by PCR to determine the prevalence and phylogeny of these organisms. RESULTS: Screening for Chlamydiales resulted in a total prevalence of 31.4%. Positive samples originated from captive and free-living bats from all three countries. Sequencing of 15 samples allowed the detection of two phylogenetically distinct groups. These groups share sequence identities to Chlamydiaceae, and to Chlamydia-like organisms including Rhabdochlamydiaceae and unclassified Chlamydiales from environmental samples, respectively. PCR analysis for the presence of hemotropic mycoplasmas resulted in a total prevalence of 0.7%, comprising free-living bats from Germany and Costa Rica. Phylogenetic analysis revealed three sequences related to other unidentified mycoplasmas found in vampire bats and Chilean bats. CONCLUSIONS: Bats can harbor Chlamydiales and hemotropic mycoplasmas and the newly described sequences in this study indicate that the diversity of these bacteria in bats is much larger than previously thought. Both, Chlamydiales and hemotropic mycoplasmas are not restricted to certain bat species or countries and captive and free-living bats can be colonized. In conclusion, bats represent another potential host or vector for novel, previously unidentified, Chlamydiales and hemotropic mycoplasmas.

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.

Clinical impact of the type VI secretion system on virulence of Campylobacter species during infection.

BACKGROUND: The clinical course of Campylobacter infection varies in symptoms and severity depending on host factors, virulence of the pathogen and initiated therapy. The type VI secretion system (T6SS) has been identified as a novel virulence factor, which mediates contact-dependent injection of enzymes and toxins into competing bacteria or host cells and facilitates the colonisation of a host organism. We aimed to compare the clinical course of Campylobacter infection caused by strains with and without the T6SS and identify possible associations between this putative virulence factor and the clinical manifestations of disease. METHODS: From April 2015 to January 2017, patients with detection of Campylobacter spp. were identified at the University Hospital of Basel and the University Children's Hospital of Basel and included in this case-control study. Presence of the T6SS gene cluster was assayed by PCR targeting the hcp gene, confirmed with whole genome sequencing. Pertinent clinical data was collected by medical record review. Differences in disease- and host-characteristics between T6SS-positive (case) and -negative (control) were compared in a uni- and multi-variable analysis. Hospital admission, antibiotic therapy, admission to intensive care unit, development of bacteraemia and in-hospital mortality were considered as clinical endpoints. RESULTS: We identified 138 cases of Campylobacter jejuni infections and 18 cases of Campylobacter coli infections from a paediatric and adult population. Analyses were focused on adult patients with C. jejuni (n = 119) of which 16.8% were T6SS-positive. Comparisons between T6SS-positive and -negative C. jejuni isolates did not reveal significant differences regarding clinical manifestations or course of disease. All clinical endpoints showed a similar distribution in both groups. A higher score in the Charlson Comorbidity Index was associated with T6SS-positive C. jejuni isolates (p < 0.001) and patients were more likely to have a solid organ transplant and to be under immunosuppressive therapy. CONCLUSIONS: Our study does not provide evidence that T6SS is associated with a more severe clinical course. Interestingly, T6SS-positive isolates are more commonly found in immunocompromised patients: an observation which merits further investigation.

Patterns of genome evolution that have accompanied host adaptation in Salmonella.

Many bacterial pathogens are specialized, infecting one or few hosts, and this is often associated with more acute disease presentation. Specific genomes show markers of this specialization, which often reflect a balance between gene acquisition and functional gene loss. Within Salmonella enterica subspecies enterica, a single lineage exists that includes human and animal pathogens adapted to cause infection in different hosts, including S. enterica serovar Enteritidis (multiple hosts), S. Gallinarum (birds), and S. Dublin (cattle). This provides an excellent evolutionary context in which differences between these pathogen genomes can be related to host range. Genome sequences were obtained from ∼ 60 isolates selected to represent the known diversity of this lineage. Examination and comparison of the clades within the phylogeny of this lineage revealed signs of host restriction as well as evolutionary events that mark a path to host generalism. We have identified the nature and order of events for both evolutionary trajectories. The impact of functional gene loss was predicted based upon position within metabolic pathways and confirmed with phenotyping assays. The structure of S. Enteritidis is more complex than previously known, as a second clade of S. Enteritidis was revealed that is distinct from those commonly seen to cause disease in humans or animals, and that is more closely related to S. Gallinarum. Isolates from this second clade were tested in a chick model of infection and exhibited a reduced colonization phenotype, which we postulate represents an intermediate stage in pathogen-host adaptation.

Culture-independent genomics of a novel chlamydial pathogen of fish provides new insight into host-specific adaptations utilized by these intracellular bacteria.

Several Chlamydiales families are associated with epitheliocystis, a common condition of the fish gill epithelium. These families share common ancestors with the Chlamydiaceae and environmental Chlamydiae. Due to the lack of culture systems, little is known about the biology of these chlamydial fish pathogens. We investigated epitheliocystis in cultured Orange-spotted grouper (Epinephelus coioides) from North Queensland, Australia. Basophilic inclusions were present in the gills of 22/31 fish and the presence of the chlamydial pathogen in the cysts was confirmed by in situ hybridization. Giant grouper (Epinephelus lanceolatus) cultured in the same systems were epitheliocystis free. 16S rRNA gene sequencing revealed a novel member of the Candidatus Parilichlamydiaceae: Ca. Similichlamydia epinephelii. Using metagenomic approaches, we obtained an estimated 68% of the chlamydial genome, revealing that this novel chlamydial pathogen shares a number of key pathogenic hallmarks with the Chlamydiaceae, including an intact Type III Secretion system and several chlamydial virulence factors. This provides additional evidence that these pathogenic mechanisms were acquired early in the evolution of this unique bacterial phylum. The identification and genomic characterization of Ca. S. epinephelii provides new opportunities to study the biology of distantly-related chlamydial pathogens while shining a new light on the evolution of pathogenicity of the Chlamydiaceae.

Whole-genome sequences of Chlamydia trachomatis directly from clinical samples without culture.

The use of whole-genome sequencing as a tool for the study of infectious bacteria is of growing clinical interest. Chlamydia trachomatis is responsible for sexually transmitted infections and the blinding disease trachoma, which affect hundreds of millions of people worldwide. Recombination is widespread within the genome of C. trachomatis, thus whole-genome sequencing is necessary to understand the evolution, diversity, and epidemiology of this pathogen. Culture of C. trachomatis has, until now, been a prerequisite to obtain DNA for whole-genome sequencing; however, as C. trachomatis is an obligate intracellular pathogen, this procedure is technically demanding and time consuming. Discarded clinical samples represent a large resource for sequencing the genomes of pathogens, yet clinical swabs frequently contain very low levels of C. trachomatis DNA and large amounts of contaminating microbial and human DNA. To determine whether it is possible to obtain whole-genome sequences from bacteria without the need for culture, we have devised an approach that combines immunomagnetic separation (IMS) for targeted bacterial enrichment with multiple displacement amplification (MDA) for whole-genome amplification. Using IMS-MDA in conjunction with high-throughput multiplexed Illumina sequencing, we have produced the first whole bacterial genome sequences direct from clinical samples. We also show that this method can be used to generate genome data from nonviable archived samples. This method will prove a useful tool in answering questions relating to the biology of many difficult-to-culture or fastidious bacteria of clinical concern.

The emergence of epitheliocystis in the upper Rhone region: evidence for Chlamydiae in wild and farmed salmonid populations.

We present the first study comparing epitheliocystis in a wild and farmed salmonid in Europe. Sampling three tributaries to the Lake Geneva, including one from headwaters to river mouth, revealed an unequal distribution of epitheliocystis in brown trout (Salmo trutta). When evaluated histologically and comparing sites grouped as wild versus farm, the probability of finding infected trout is higher on farms. In contrast, the infection intensities, as estimated by the number of cysts per gill arch, were higher on average and showed maximum values in the wild trout. Sequence analysis showed the most common epitheliocystis agents were Candidatus Piscichlamydia salmonis, all clustering into a single clade, whereas Candidatus Clavichlamydia salmonicola sequences cluster in two closely related subspecies, of which one was mostly found in farmed fish and the other exclusively in wild brown trout, indicating that farms are unlikely to be the source of infections in wild trout. A detailed morphological analysis of cysts using transmission electron microscopy revealed unique features illustrating the wide divergence existing between Ca. P. salmonis and Ca. C. salmonicola within the phylum Chlamydiae.

Genomic Investigations unmask Mycoplasma amphoriforme, a new respiratory pathogen.

BACKGROUND: Mycoplasma amphoriforme has been associated with infection in patients with primary antibody deficiency (PAD). Little is known about the natural history of infection with this organism and its ability to be transmitted in the community. METHODS: The bacterial load was estimated in sequential sputum samples from 9 patients by quantitative polymerase chain reaction. The genomes of all available isolates, originating from patients in the United Kingdom, France, and Tunisia, were sequenced along with the type strain. Genomic data were assembled and annotated, and a high-resolution phylogenetic tree was constructed. RESULTS: By using high-resolution whole-genome sequencing (WGS) data, we show that patients can be chronically infected with M. amphoriforme manifesting as a relapsing-remitting bacterial load, interspersed by periods when the organism is undetectable. Importantly, we demonstrate transmission of strains within a clinical environment. Antibiotic resistance mutations accumulate in isolates taken from patients who received multiple courses of antibiotics. CONCLUSIONS: Mycoplasma amphoriforme isolates form a closely related species responsible for a chronic relapsing and remitting infection in PAD patients in the United Kingdom and from immunocompetent patients in other countries. We provide strong evidence of transmission between patients attending the same clinic, suggesting that screening and isolation may be necessary for susceptible patients. This work demonstrates the critical role that WGS can play in rapidly unraveling the biology of a novel pathogen.

Rapid detection of plasmid-mediated AmpC-producers by eazyplex® SuperBug AmpC assay compared to whole-genome sequencing.

Current methods for plasmid-mediated AmpC ß-lactamase (pAmpC) detection in routine microbiological laboratories are based on various phenotypic tests. Eazyplex®SuperBug AmpC assay is a molecular assay based on isothermal amplification for rapid detection of the most common pAmpC types from bacterial culture: CMY-2 group, DHA, ACC and MOX. Our aim was to evaluate the diagnostic performance of this assay. The assay was evaluated on 64 clinical isolates of Enterobacterales without chromosomal inducible AmpC, and with phenotypically confirmed AmpC production. The results were confirmed, and isolates further characterized by whole-genome sequencing (WGS). eazyplex®SuperBug AmpC assay correctly detected the two most common pAmpC types CMY-2 group (16/16) and DHA (19/19). Detection of ACC and MOX could not be evaluated on our set of isolates since there was only one isolate harbouring ACC and none with MOX. pAmpC encoding genes could be detected in only eight of 36 investigated Escherichia coli isolates. The remaining 28 E. coli isolates harboured previously described mutations in the blaEC promoter, leading to the overexpression of chromosomally encoded E. coli specific AmpC ß-lactamase. All results were 100% concordant with the results of WGS. eazyplex®SuperBug AmpC assay enabled rapid and reliable detection of pAmpC-encoding genes in Enterobacterales like Klebsiella spp. and Proteus spp. and the distinction between plasmid-mediated and chromosomally encoded AmpC in E. coli.

S-Gene Target Failure as an Effective Tool for Tracking the Emergence of Dominant SARS-CoV-2 Variants in Switzerland and Liechtenstein, Including Alpha, Delta, and Omicron BA.1, BA.2, and BA.4/BA.5.

During the SARS-CoV-2 pandemic, the Dr. Risch medical group employed the multiplex TaqPathTM COVID-19 CE-IVD RT-PCR Kit for large-scale routine diagnostic testing in Switzerland and the principality of Liechtenstein. The TaqPath Kit is a widely used multiplex assay targeting three genes (i.e., ORF1AB, N, S). With emergence of the B.1.1.7 (Alpha) variant, a diagnostic flaw became apparent as the amplification of the S-gene target was absent in these samples due to a deletion (ΔH69/V70) in the Alpha variant genome. This S-gene target failure (SGTF) was the earliest indication of a new variant emerging and was also observed in subsequent variants such as Omicron BA.1 and BA4/BA.5. The Delta variant and Omicron BA.2 did not present with SGTF. From September 2020 to November 2022, we investigated the applicability of the SGTF as a surrogate marker for emerging variants such as B.1.1.7, B.1.617.2 (Delta), and Omicron BA.1, BA.2, and BA.4/BA.5 in samples with cycle threshold (Ct) values < 30. Next to true SGTF-positive and SGTF-negative samples, there were also samples presenting with delayed-type S-gene amplification (higher Ct value for S-gene than ORF1ab gene). Among these, a difference of 3.8 Ct values between the S- and ORF1ab genes was found to best distinguish between "true" SGTF and the cycle threshold variability of the assay. Samples above the cutoff were subsequently termed partial SGTF (pSGTF). Variant confirmation was performed by whole-genome sequencing (Oxford Nanopore Technology, Oxford, UK) or mutation-specific PCR (TIB MOLBIOL). In total, 17,724 (7.4%) samples among 240,896 positives were variant-confirmed, resulting in an overall sensitivity and specificity of 93.2% [92.7%, 93.7%] and 99.3% [99.2%, 99.5%], respectively. Sensitivity was increased to 98.2% [97.9% to 98.4%] and specificity lowered to 98.9% [98.6% to 99.1%] when samples with pSGTF were included. Furthermore, weekly logistic growth rates (α) and sigmoid's midpoint (t0) were calculated based on SGTF data and did not significantly differ from calculations based on comprehensive data from GISAID. The SGTF therefore allowed for a valid real-time estimate for the introduction of all dominant variants in Switzerland and Liechtenstein.

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.