Core Genome Multilocus Sequence Research Articles

Genomic Diversity, Antimicrobial Susceptibility, and Biofilm Formation of Clinical Acinetobacter baumannii Isolates from Horses.

Acinetobacter (A.) baumannii is an opportunistic pathogen that causes severe infections in humans and animals, including horses. The occurrence of dominant international clones (ICs), frequent multidrug resistance, and the capability to form biofilms are considered major factors in the successful spread of A. baumannii in human and veterinary clinical environments. Since little is known about A. baumannii isolates from horses, we studied 78 equine A. baumannii isolates obtained from clinical samples between 2008 and 2020 for their antimicrobial resistance (AMR), clonal distribution, biofilm-associated genes (BAGs), and biofilm-forming capability. Based on whole-genome sequence analyses, ICs, multilocus (ML) and core-genome ML sequence types (STs), and AMR genes were determined. Antimicrobial susceptibility testing was performed by microbroth dilution. A crystal violet assay was used for biofilm quantification. Almost 37.2% of the isolates were assigned to IC1 (10.3%), IC2 (20.5%), and IC3 (6.4%). Overall, the isolates revealed high genomic diversity. We identified 51 different STs, including 22 novel STs (ST1723-ST1744), and 34 variants of the intrinsic oxacillinase (OXA), including 8 novel variants (OXA-970 to OXA-977). All isolates were resistant to ampicillin, amoxicillin/clavulanic acid, cephalexin, cefpodoxime, and nitrofurantoin. IC1-IC3 isolates were also resistant to gentamicin, enrofloxacin, marbofloxacin, tetracycline, and trimethoprim/sulfamethoxazole. All isolates were susceptible to imipenem. Thirty-one multidrug-resistant (MDR) isolates mainly accumulated in the IC1-IC3 groups. In general, these isolates showed less biofilm formation (IC1 = 25.0%, IC2 = 18.4%, IC3 = 15.0%) than the group of non-IC1-IC3 isolates (58.4%). Isolates belonging to the same ICs/STs revealed identical BAG patterns. BAG blp1 was absent in all isolates, whereas bfmR and pgaA were present in all isolates. At the level of the IC groups, the AMR status was negatively correlated with the isolates' ability to form a biofilm. A considerable portion of equine A. baumannii isolates revealed ICs/STs that are globally present in humans. Both an MDR phenotype and the capability to form biofilms might lead to therapeutic failures in equine medicine, particularly due to the limited availability of licensed drugs.

Clostridioides difficile in South American Camelids in Germany: First Insights into Molecular and Genetic Characteristics and Antimicrobial Resistance.

Little is known about zoonotic pathogens and their antimicrobial resistance in South American camelids (SAC) in Germany including Clostridioides (C.) difficile. The aim of this study was to investigate prevalence, molecular characteristics and antimicrobial resistance of C. difficile in SAC. Composite SAC faecal samples were collected in 43 husbandries in Central Germany and cultured for C. difficile. Toxinotyping and ribotyping was done by PCR. Whole genome sequencing was performed with Illumina® Miseq™. The genomes were screened for antimicrobial resistance determinants. Genetic relatedness of the isolates was investigated using core genome multi locus sequence typing (cgMLST) and single nucleotide polymorphism analysis. Antimicrobial susceptibility testing was done using the Etest® method. Eight C. difficile isolates were recovered from seven farms. The isolates belonged to different PCR ribotypes. All isolates were toxinogenic. cgMLST revealed a cluster containing isolates recovered from different farms. Seven isolates showed similar resistance gene patterns. Different phenotypic resistance patterns were found. Agreement between phenotypic and genotypic resistance was identified only in some cases. Consequently, SAC may act as a reservoir for C. difficile. Thus, SAC may pose a risk regarding zoonotic transmission of toxinogenic, potentially human-pathogenic and resistant C. difficile isolates.

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

BackgroundMultidrug-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.MethodsIn 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.ResultsIn 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.ConclusionsOur 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.

Evaluation of within-host evolution of methicillin-resistant Staphylococcus aureus (MRSA) by comparing cgMLST and SNP analysis approaches

Whole genome sequencing (WGS) of methicillin-resistant Staphylococcus aureus (MRSA) provides high-resolution typing, facilitating surveillance and outbreak investigations. The aim of this study was to evaluate the genomic variation rate in MRSA, by comparing commonly used core genome multilocus sequencing (cgMLST) against single nucleotide polymorphism (SNP) analyses. WGS was performed on 95 MRSA isolates, collected from 20 carriers during years 2003–2019. To assess variation and methodological-related differences, two different cgMLST schemes were obtained using Ridom SeqSphere+ and the cloud-based 1928 platform. In addition, two SNP methods, 1928 platform and Northern Arizona SNP Pipeline (NASP) were used. The cgMLST using Ridom SeqSphere+ and 1928 showed a median of 5.0 and 2.0 allele variants/year, respectively. In the SNP analysis, performed with two reference genomes COL and Newman, 1928 showed a median of 13 and 24 SNPs (including presumed recombination) and 3.8 respectively 4.0 SNPs (without recombination) per individual/year. Accordantly, NASP showed a median of 5.5 and 5.8 SNPs per individual/year. In conclusion, an estimated genomic variation rate of 2.0–5.8 genetic events per year (without recombination), is suggested as a general guideline to be used at clinical laboratories for surveillance and outbreak investigations independently of analysis approach used.

Guidelines for reporting Whole Genome Sequencing‐based typing data through the EFSA One Health WGS System

In 2019, the European Food Safety Authority (EFSA) and the European Centre for Disease Prevention and Control (ECDC) received a mandate from the European Commission (EC) to implement and manage two interoperable systems for the collection and analysis of Whole Genome Sequence (WGS) data of Salmonella enterica, Listeria monocytogenes and Escherichia coli, including Shiga toxin-producing E. coli (STEC) from humans and from food, feed, animals and the related environment, to contribute to the epidemiological investigations of food-borne outbreaks and to the identification of emerging health threats. As response to this EC mandate, EFSA has developed the EFSA One Health WGS System that interoperates with the ECDC Molecular Typing system exchanging core genome Multi Locus Sequence Typing (cgMLST) profiles and minimum metadata. Interested European Union (EU) Member States (MS) and European Economic Area (EEA) countries can submit their typing data from food, feed and animal samples to the EFSA One Health WGS system either through an user-friendly interface or programmatically. This technical report describes the EFSA One Health WGS system, including his architecture, the users and their relative data visibility, as well as data ownership and intellectual property. This report provides also detailed information on the types and formats of data that can be submitted to the EFSA One Health WGS system. Moreover, these guidelines describe the individual steps of the reporting procedure in order to support the EU MSs and EEA countries in data submission and data management.

Sequencing-Based Genotyping of Pakistani Burkholderia mallei Strains: A Useful Way for Investigating Glanders Outbreaks

Burkholderia (B.) mallei is a host-adapted equine pathogen that causes glanders, a re-emerging zoonotic disease, which is endemic in Pakistan and other developing countries and seriously impacts the global equine movement. Due to globalization, the geographical restriction of diseases vanishes and the lack of awareness of and experience with eradicated diseases in industrialized countries also promotes the re-introduction of infections in these regions. Owing to the high equine population, the Pakistani province Punjab is a potential hotspot where several glanders outbreaks have been seen over last two decades. For determining the genomic diversity of B. mallei in this and other equine-populated prefectures, the genomes of 19 B. mallei strains isolated between 1999 and 2020 in different locations were sequenced and their genotypes were determined. Particularly, for genetically highly homogenous pathogens like B. mallei genotyping techniques require a high discriminatory power for enabling differentiation on the strain level. Thus, core-genome single nucleotide polymorphism (cgSNP) analysis was applied for distinguishing the highly similar strains. Furthermore, a whole-genome sequence-based core genome multi locus sequence typing (cgMLST) scheme, specific to B. mallei, was developed and additionally applied to the data. It was found that B. mallei genotypes in Pakistan persisted over time and space and genotype clusters preferred connection with a time point rather than the place of isolation, probably due to frequent equine movement, which promotes the spread of glanders. The cgMLST approach proved to work in accord with SNP typing and may help to investigate future glanders outbreaks.

Invasive listeriosis outbreaks and salmon products: a genomic, epidemiological study

ABSTRACT Invasive listeriosis, caused by Listeria (L.) monocytogenes, is a severe foodborne infection, especially for immunocompromised individuals. The aim of our investigation was the identification and analysis of listeriosis outbreaks in Germany with smoked and graved salmon products as the most likely source of infection using whole-genome sequencing (WGS) and patient interviews. In a national surveillance programme, WGS was used for subtyping and core genome multi locus sequence typing (cgMLST) for cluster detection of L. monocytogenes isolates from listeriosis cases as well as food and environmental samples in Germany. Patient interviews were conducted to complement the molecular typing. We identified 22 independent listeriosis outbreaks occurring between 2010 and 2021 that were most likely associated with the consumption of smoked and graved salmon products. In Germany, 228 cases were identified, of 50 deaths (22%) reported 17 were confirmed to have died from listeriosis. Many of these 22 outbreaks were cross-border outbreaks with further cases in other countries. This report shows that smoked and graved salmon products contaminated with L. monocytogenes pose a serious risk for listeriosis infection in Germany. Interdisciplinary efforts including WGS and epidemiological investigations were essential to identifying the source of infection. Uncooked salmon products are high-risk foods frequently contaminated with L. monocytogenes. In order to minimize the risk of infection for consumers, food producers need to improve hygiene measures and reduce the entry of pathogens into food processing. Furthermore, susceptible individuals should be better informed of the risk of acquiring listeriosis from consuming smoked and graved salmon products.

Population Analysis of O26 Shiga Toxin-Producing Escherichia coli Causing Hemolytic Uremic Syndrome in Italy, 1989-2020, Through Whole Genome Sequencing.

Shiga toxin-producing Escherichia coli (STEC) belonging to the O26 serogroup represent an important cause of Hemolitic Uremic Syndrome (HUS) in children worldwide. The localization of STEC virulence genes on mobile genetic elements allowed the emergence of clones showing different assets of this accessory genomic fraction. A novel O26 STEC clone belonging to Sequence Type (ST) 29 and harboring stx2a, ehxA and etpD plasmid-borne genes has emerged and spread in Europe since the mid-1990s, while another ST29 clone positive for stx2d and lacking plasmid-borne virulence genes was recently described as emerging in France. In Italy, O26 has been the most frequently detected STEC serogroup from HUS cases since the late 1990s. In this study we describe the genomic characterization and population structure of 144 O26 STEC strains isolated from human sources in Italy in the period 1989-2020. A total of 89 strains belonged to ST21, 52 to ST29, two to ST396 and one to ST4944. ST29 strains started to be isolated from 1999. 24 strains were shown to harbour stx1a, alone (n=20) or in combination with stx2a (n=4). The majority of the strains (n=118) harbored stx2a genes only and the two ST396 strains harbored stx2d. A Hierarchical Clustering on Principal Components (HCPC) analysis, based on the detection of accessory virulence genes, antimicrobial resistance (AMR) genes and plasmid replicons, classified the strains in seven clusters identified with numbers from 1 to 7, containing two, 13, 39, 63, 16, 10 and one strain, respectively. The majority of the genetic features defining the clusters corresponded to plasmid-borne virulence genes, AMR genes and plasmid replicons, highlighting specific assets of plasmid-borne features associated with different clusters. Core genome Multi Locus Sequence Typing grouped ST21 and ST29 strains in three clades each, with each ST29 clade exactly corresponding to one HCPC cluster. Our results showed high conservation of either the core or the accessory genomic fraction in populations of ST29 O26 STEC, differently from what observed in ST21 strains, suggesting that a different selective pressure could drive the evolution of different populations of these pathogens possibly involving different ecological niches.

Population dynamics in colonizing vancomycin-resistant Enterococcus faecium isolated from immunosuppressed patients

Vancomycin-resistant Enterococcus faecium and Enterococcus faecalis (VRE) are a common cause of healthcare-associated infections. Whole genome sequencing-based typing methods yield the highest discriminatory power for outbreak surveillance in the hospital. We analysed the clonal composition of enteric VRE populations of at-risk patients over several weeks to characterise VRE population diversity and dynamics. Five bone marrow transplant recipients (three colonised with vanA-positive isolates, two colonised with vanB-positive isolates) contributed three rectal swabs over a course of several weeks. Fourteen VRE colonies per swab were analysed by core genome multi locus sequence typing (cgMLST) and typing of the van-element. VRE populations were clonally diverse in three of five patients, and population composition changed dynamically over the time of observation. Besides new acquisition of VRE isolates, shared van-elements localised on nearly identical plasmids between clonally different isolates indicate horizontal gene transfer as a mechanism behind VRE population diversity within single patients. Outbreak detection relies on typing of isolates, usually by analysing one isolate per patient. We here show that this approach is insufficient for outbreak surveillance of VRE in highly vulnerable patients, as it does not take into account VRE population heterogeneity and horizontal gene transfer of the resistance element.

Genetic diversity of Listeria monocytogenes strains contaminating food and food producing environment as single based sample in Italy (retrospective study)

Human listeriosis outbreaks are often associated with food products, which could be contaminated, at the same time, also by different clones of Listeria monocytogenes. This emphasize the need to type more than one L.monocytogenes isolate found in a single food or environmental sample.The purpose of the present study was to evaluate the presence of different L.monocytogenes strains in food and food production environment in order to understand if there is need to type more isolates from the same sample in case of presence of L.monocytogenes. Between 2011 and 2015, at the Italian National Reference Laboratory for L.monocytogenes, for each positive sample, from two to twenty-three isolates of L.monocytogenes were collected. All the isolates were characterized by conventional serotyping and pulsed field gel electrophoresis (PFGE). Moreover, isolates from the same sample, having indistinguishable PFGE profile, were subjected to whole genome sequencing in order to perform core genome Multi Locus Sequence Typing (cgMLST).Within each sample, more than one serotype and one pulsotype were found in 11.9% and 27.5%, respectively. For indistinguishable PFGE patterns the cgMLST analysis showed 96.2% of concordance demonstrating the added value of new sequencing technologies.This study has demonstrated the need to select and type more than one L.monocytogenes colony in one food or food environmental sample to detect the diversity of L.monocytogenes strains and facilitate downstream investigations and effective source attribution in foodborne outbreak inquiry.

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

BackgroundColistin 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.ResultsOf 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.ConclusionsThe 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.

A core genome multilocus sequence typing scheme for Mycoplasma hyorhinis

Mycoplasma (M.) hyorhinis is a commensal and pathobiont residing in the upper respiratory tract in swine and with the ability to spread systemically, mainly causing polyserositis and polyarthritis in nursery pigs. Since little is known on the epidemiology of M. hyorhinis infection, whole genome sequences of 73 strains isolated from pigs in Austria (n = 71) and Germany (n = 2), that have been isolated from clinically affected pigs during routine diagnostics, and publicly available genomes of eight M. hyorhinis strains were analyzed in the presented study. For this purpose, a core genome multi locus sequence typing (cgMLST) scheme encompassing 453 target genes was developed using the Ridom© SeqSphere + software. Results were compared to two previously described conventional MLST schemes and to a core genome single nucleotide polymorphism (cgSNP) analysis approach. Core genome MLST showed high diversity among the M. hyorhinis strains studied and while certain isolates from one farm or a single animal formed cgMLST clusters (≤ 8 allele differences), no isolates with identical allele profiles were identified. In addition, cgMLST had superior discriminatory power (Simpson’s ID = 0.995) over conventional MLST (Simpson’s ID = 0.952 and 0.985), while demonstrating a lack of congruence between conventional MLST and genome-wide relationship. Core genome SNP results were highly congruent with cgMLST results but lacked in resolution when comparing closely related isolates. Thus, cgMLST is the most suitable method for epidemiological investigations such as outbreak analysis, and to gain insights into M. hyorhinis population structure.

Molecular Epidemiology of Carbapenem-Resistant Acinetobacter baumannii Isolates from Northern Africa and the Middle East.

At the Bundeswehr Hospitals of Hamburg and Westerstede, patients repatriated from subtropical war and crisis zones of Northern Africa and the Middle East were medically treated, including microbiological assessment. Within a six-year interval, 16 Acinetobacter spp. strains, including 14 Acinetobacter baumannii (Ab) isolates with resistance against carbapenems and origins in Afghanistan (n = 4), Iraq (n = 2), Libya (n = 2), and Syria (n = 8) were collected. While clonal relationships of Libyan and Syrian strains had been assessed by superficial next generation sequencing (NGS) and “DiversiLab” repetitive elements sequence-based (rep-)PCR so far, this study provides core genome-based sequence typing and thus more detailed epidemiological information. In detail, sequencing allowed a definitive species identification and comparison with international outbreak-associated Ab strains by core genome multi locus sequence typing (cgMLST) and the identification of MLST lineages, as well as the identification of known resistance genes. The sequence analysis allowed for the confirmation of outbreak-associated clonal clusters among the Syrian and Afghan Ab isolates, indicating likely transmission events. The identified acquired carbapenem resistance genes comprised blaOXA-23, blaOXA-58, blaNDM-1, and blaGES-11, next to other intrinsic and acquired, partly mobile resistance-associated genes. Eleven out of 14 Ab isolates clustered with the previously described international clonal lineages IC1 (4 Afghan strains), IC2 (6 Syrian strains), and IC7 (1 Syrian strain). Identified Pasteur sequence types of the 14 Ab strains comprised ST2 (Syrian), ST25 (Libyan), ST32 (Iraqi), ST81 (Afghan), ST85 (Libyan), and ST1112 (Syrian), respectively. In conclusion, the study revealed a broad spectrum of resistance genes in Ab isolated from war-injured patients from Northern Africa and the Middle East, thereby broadening the scarcely available data on locally abundant clonal lineages and resistance mechanisms.

Comparison of Whole Genome (wg-) and Core Genome (cg-) MLST (BioNumericsTM) Versus SNP Variant Calling for Epidemiological Investigation of Pseudomonas aeruginosa.

Whole genome sequencing (WGS) is increasingly used for epidemiological investigations of pathogens. While SNP variant calling is currently considered as the most suitable method, the choice of a representative reference genome and the isolate dependency of results limit standardization and affect resolution in an unknown manner. Whole or core genome Multi Locus Sequence Typing (wg-, cg-MLST) represents an attractive alternative. Here, we assess the accuracy of wg- and cg-MLST by comparing results of four Pseudomonas aeruginosa datasets for which epidemiological and genomic data were previously described. Three datasets included 155 isolates from three different sequence types (ST) of P. aeruginosa collected in our ICUs over a 5-year period. The fourth dataset consisted of 10 isolates from an investigation of P. aeruginosa contaminated hand soap. All isolates were previously analyzed by a core SNP approach. In this study, wg- and cg-MLST were performed in BioNumericsTM using a scheme developed by Applied-Maths. Correlation between SNP calling and wg- or cg-MLST results were evaluated by calculating linear regressions and their coefficient of correlations (R2) between the number of SNPs and the number of allele differences in pairwise comparison of isolates. The number of SNPs and allele difference between isolates with close epidemiological linkage varies between 0–26 and 0–13, respectively. When compared to core-SNP calling, a higher coefficient of correlation was obtained with cgMLST (R2 of 0.92–0.99) than with wgMLST (0.78–0.99). In one dataset, a putative homologous recombination of a large DNA fragment (202 loci) was identified among these isolates, affecting its phylogeny, but with no impact on the epidemiological analysis of outbreak isolates. In conclusion, we showed that the P. aeruginosa wgMLST scheme in BioNumericsTM is as discriminatory as the core-SNP calling approach and apparently useful for outbreak investigations. We also showed that epidemiological linked isolates showed less than 26 SNPs or 13 allele differences. These are important figures for the distinction between outbreak and non-outbreak isolates when interpreting WGS results. However, as P. aeruginosa is highly recombinant, a cgMLST approach is preferable and caution should be addressed to possible recombination of large DNA fragments.

Identification of Antimicrobial Resistance-Associated Genes through Whole Genome Sequencing of Mycoplasma bovis Isolates with Different Antimicrobial Resistances.

Antimicrobial resistance (AMR) in Mycoplasma bovis has been previously associated with topoisomerase and ribosomal gene mutations rather than specific resistance-conferring genes. Using whole genome sequencing (WGS) to identify potential new AMR mechanisms for M. bovis, it was found that a 2019 clinical isolate with high MIC (2019-043682) for fluoroquinolones, macrolides, lincosamides, pleuromutilins and tetracyclines had a new core genome multilocus sequencing (cgMLST) type (ST10-like) and 91% sequence similarity to the published genome of M. bovis PG45. Closely related to PG45, a 1982 isolate (1982-M6152) shared the same cgMLST type (ST17), 97.2% sequence similarity and low MIC results. Known and potential AMR- associated genetic events were identified through multiple sequence alignment of the three genomes. Isolate 2019-043682 had 507 genes with non-synonymous mutations (NSMs) and 67 genes disrupted. Isolate 1982-M6152 had 81 NSMs and 20 disruptions. Using functional roles and known mechanisms of antimicrobials, a 55 gene subset was assessed for AMR potential. Seventeen were previously identified from other bacteria as sites of AMR mutation, 38 shared similar functions to them, and 11 contained gene-disrupting mutations. This study indicated that M. bovis may obtain high AMR characteristics by mutating or disrupting other functional genes, in addition to topoisomerases and ribosomal genes.

Population structure-guided profiling of antibiotic resistance patterns in clinical Listeria monocytogenes isolates from Germany identifies pbpB3 alleles associated with low levels of cephalosporin resistance

ABSTRACT Numbers of listeriosis illnesses have been increasing in Germany and the European Union during the last decade. In addition, reports on the occurrence of antibiotic resistance in Listeria monocytogenes in clinical and environmental isolates are accumulating. The susceptibility towards 14 antibiotics was tested in a selection of clinical L. monocytogenes isolates to get a more precise picture of the development and manifestation of antibiotic resistance in the L. monocytogenes population. Based on the population structure determined by core genome multi locus sequence typing (cgMLST) 544 out of 1220 sequenced strains collected in Germany between 2009 and 2019 were selected to cover the phylogenetic diversity observed in the clinical L. monocytogenes population. All isolates tested were susceptible towards ampicillin, penicillin and co-trimoxazole – the most relevant antibiotics in the treatment of listeriosis. Resistance to daptomycin and ciprofloxacin was observed in 493 (91%) and in 71 (13%) of 544 isolates, respectively. While all tested strains showed resistance towards ceftriaxone, their resistance levels varied widely between 4 mg/L and >128 mg/L. An allelic variation of the penicillin binding protein gene pbpB3 was identified as the cause of this difference in ceftriaxone resistance levels. This study is the first population structure-guided analysis of antimicrobial resistance in recent clinical isolates and confirms the importance of penicillin binding protein B3 (PBP B3) for the high level of intrinsic cephalosporin resistance of L. monocytogenes on a population-wide scale.

Whole genome sequencing reveals resemblance between ESBL-producing and carbapenem resistant Klebsiella pneumoniae isolates from Austrian rivers and clinical isolates from hospitals

In 2016, the Austrian Agency for Health and Food Safety started a pilot project to investigate antimicrobial resistance in surface water. Here we report on the characterisation of carbapenem resistant and ESBL-producing K. pneumoniae isolates from Austrian river water samples compared to 95 clinical isolates recently obtained in Austrian hospitals.Ten water samples were taken from four main rivers, collected upstream and downstream of major cities in 2016. For subtyping and comparison, public core genome multi locus sequence typing (cgMLST) schemes were used. The presence of AMR genes, virulence genes and plasmids was extracted from whole genome sequence (WGS) data.In total three ESBL-producing strains and two carbapenem resistant strains were isolated. WGS based comparison of these five water isolates to 95 clinical isolates identified three clusters. Cluster 1 (ST11) and cluster 2 (ST985) consisted of doublets of carbapenem resistant strains (one water and one clinical isolate each). Cluster 3 (ST405) consisted of three ESBL-producing strains isolated from one water sample and two clinical specimens. The cities, in which patient isolates of cluster 2 and 3 were collected, were in concordance with the water sampling locations downstream from these cities. The genetic concordance between isolates from river water samples and patient isolates raises concerns regarding the release of wastewater treatment plant effluents into surface water. From a public health perspective these findings demand attention and strategies are required to minimize the spread of multiresistant strains to the environment.

A Listeria monocytogenes ST2 clone lacking chitinase ChiB from an outbreak of non-invasive gastroenteritis

ABSTRACT An outbreak with a remarkable Listeria monocytogenes clone causing 163 cases of non-invasive listeriosis occurred in Germany in 2015. Core genome multi locus sequence typing grouped non-invasive outbreak isolates and isolates obtained from related food samples into a single cluster, but clearly separated genetically close isolates obtained from invasive listeriosis cases. A comparative genomic approach identified a premature stop codon in the chiB gene, encoding one of the two L. monocytogenes chitinases, which clustered with disease outcome. Correction of this premature stop codon in one representative gastroenteritis outbreak isolate restored chitinase production, but effects in infection experiments were not found. While the exact role of chitinases in virulence of L. monocytogenes is still not fully understood, our results now clearly show that ChiB-derived activity is not required to establish L. monocytogenes gastroenteritis in humans. This limits a possible role of ChiB in human listeriosis to later steps of the infection.

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

BackgroundTracking recent transmission is a vital part of controlling widespread pathogens such as Mycobacterium tuberculosis. Multiple methods with specific performance characteristics exist for detecting recent transmission chains, usually by clustering strains based on genotype similarities. With such a large variety of methods available, informed selection of an appropriate approach for determining transmissions within a given setting/time period is difficult.MethodsThis study combines whole genome sequence (WGS) data derived from 324 isolates collected 2005–2010 in Kinshasa, Democratic Republic of Congo (DRC), a high endemic setting, with phylodynamics to unveil the timing of transmission events posited by a variety of standard genotyping methods. Clustering data based on Spoligotyping, 24-loci MIRU-VNTR typing, WGS based SNP (Single Nucleotide Polymorphism) and core genome multi locus sequence typing (cgMLST) typing were evaluated.FindingsOur results suggest that clusters based on Spoligotyping could encompass transmission events that occurred almost 200 years prior to sampling while 24-loci-MIRU-VNTR often represented three decades of transmission. Instead, WGS based genotyping applying low SNP or cgMLST allele thresholds allows for determination of recent transmission events, e.g. in timespans of up to 10 years for a 5 SNP/allele cut-off.InterpretationWith the rapid uptake of WGS methods in surveillance and outbreak tracking, the findings obtained in this study can guide the selection of appropriate clustering methods for uncovering relevant transmission chains within a given time-period. For high resolution cluster analyses, WGS-SNP and cgMLST based analyses have similar clustering/timing characteristics even for data obtained from a high incidence setting.