AMRFinderPlus and the Reference Gene Catalog facilitate examination of the genomic links among antimicrobial resistance, stress response, and virulence

Whole-genome sequencing of carbapenem-resistant Enterobacterales isolates in southeast Louisiana reveals persistent genetic clusters spanning multiple locations

BACKGROUND

We investigated 51 g-negative carbapenem-resistant Enterobacterales (CRE) isolates collected from 22 patients over a five-year period from six health care institutions in the Ochsner Health network in southeast Louisiana.

METHODS

Short genomic reads were generated using Illumina sequencing and assembled for each isolate. Isolates were classified as Enterobacter spp. (n = 20), Klebsiella spp. (n = 30), and Escherichia coli (n = 1) and grouped into 19 different multi-locus sequence types (MLST). Species and patient-specific core genomes were constructed representing ∼50% of the chromosomal genome.

RESULTS

We identified two sets of patients with genetically related infections; in both cases, the related isolates were collected > 6 months apart, and in one case, the isolates were collected in different locations. On the other hand, we identified four sets of patients with isolates of the same species collected within 21 days from the same location; however, none had genetically related infections. Genes associated with resistance to carbapenem drugs (blaKPC and/or blaCTX-M-15) were found in 76% of the isolates. We found three blaKPC variants (blaKPC-2, blaKPC-3, and blaKPC-4) associated with four different Enterobacter MLST variants, and two blaKPC variants (blaKPC-2, blaKPC-3) associated with seven different Klebsiella MLST variants.

CONCLUSIONS

Molecular surveillance is increasingly becoming a powerful tool to understand bacterial spread in both community and clinical settings. This study provides evidence that genetically related infections in clinical settings do not necessarily reflect temporal associations, and vice versa. Our results also highlight the regional genomic and resistance diversity within related bacterial lineages.

Genomic Insights into the Microbial Agent Streptomyces albidoflavus MGMM6 for Various Biotechnology Applications

Microbial biotechnology plays a crucial role in improving industrial processes, particularly in the production of compounds with diverse applications. In this study, we used bioinformatic approaches to analyze the genomic architecture of Streptomyces albidoflavus MGMM6 and identify genes involved in various metabolic pathways that have significant biotechnological potential. Genome mining revealed that MGMM6 consists of a linear chromosome of 6,932,303 bp, with a high G+C content of 73.5%, lacking any plasmid contigs. Among the annotated genes, several are predicted to encode enzymes such as dye peroxidase, aromatic ring-opening dioxygenase, multicopper oxidase, cytochrome P450 monooxygenase, and aromatic ring hydroxylating dioxygenases which are responsible for the biodegradation of numerous endogenous and xenobiotic pollutants. In addition, we identified genes associated with heavy metal resistance, such as arsenic, cadmium, mercury, chromium, tellurium, antimony, and bismuth, suggesting the potential of MGMM6 for environmental remediation purposes. The analysis of secondary metabolites revealed the presence of multiple biosynthesis gene clusters responsible for producing compounds with potent antimicrobial and metal-chelating activities. Furthermore, laboratory tests conducted under controlled conditions demonstrated the effectiveness of MGMM6 in inhibiting phytopathogenic microbes, decolorizing and degrading aromatic triphenylmethane dyes, particularly Blue Brilliant G250, from wastewater by up to 98 ± 0.15%. Overall, the results of our study highlight the promising biotechnological potential of S. albidoflavus MGMM6.

Genomic Diversity, Antimicrobial Resistance, Plasmidome, and Virulence Profiles of Salmonella Isolated from Small Specialty Crop Farms Revealed by Whole-Genome Sequencing

Salmonella is the leading cause of death associated with foodborne illnesses in the USA. Difficulty in treating human salmonellosis is attributed to the development of antimicrobial resistance and the pathogenicity of Salmonella strains. Therefore, it is important to study the genetic landscape of Salmonella, such as the diversity, plasmids, and presence antimicrobial resistance genes (AMRs) and virulence genes. To this end, we isolated Salmonella from environmental samples from small specialty crop farms (SSCFs) in Northeast Ohio from 2016 to 2021; 80 Salmonella isolates from 29 Salmonella-positive samples were subjected to whole-genome sequencing (WGS). In silico serotyping revealed the presence of 15 serotypes. AMR genes were detected in 15% of the samples, with 75% exhibiting phenotypic and genotypic multidrug resistance (MDR). Plasmid analysis demonstrated the presence of nine different types of plasmids, and 75% of AMR genes were located on plasmids. Interestingly, five Salmonella Newport isolates and one Salmonella Dublin isolate carried the ACSSuT gene cassette on a plasmid, which confers resistance to ampicillin, chloramphenicol, streptomycin, sulfonamide, and tetracycline. Overall, our results show that SSCFs are a potential reservoir of Salmonella with MDR genes. Thus, regular monitoring is needed to prevent the transmission of MDR Salmonella from SSCFs to humans.

Population structure of community-acquired extended-spectrum beta-lactamase producing Escherichia coli and methicillin resistant Staphylococcus aureus in a French region showed no difference between urban and rural areas

Antimicrobial resistance is a global health issue and extended-spectrum β-lactamase producing Escherichia coli (ESBL-Ec) and methicillin-resistant Staphylococcus aureus (MRSA) are of particular concern. Whole genome sequencing analysis of isolates from the community is essential to understand the circulation of those multidrug-resistant bacteria. Our main objective was to determine the population structure of clinical ESBL-Ec and MRSA isolated in the community setting of a French region. For this purpose, isolates were collected from 23 sites belonging to 6 private medical biology laboratories in the Bourgogne-Franche-Comté region. One hundred ninety ESBL-Ec and 67 MRSA were sequenced using the Illumina technology. Genomic analyses were performed to determine the bacterial typing, presence of antibiotic resistance genes, metal resistance genes as well as virulence genes. Analysis showed that ST131 was the major ESBL-Ec clone circulating in the region, representing 42.1% of the ESBL-Ec isolates. The blaCTX-M genes represented 98% of blaESBL with the majority being blaCTX-M-15 (53.9%). MRSA population consisted of mainly of CC8 (50.7%) and CC5 (38.8%) clonal complexes. Interestingly, we found a prevalence of 40% of the zinc resistance gene czrC in our MRSA population. We observed no differences in our ESBL-Ec or MRSA populations between urban and rural areas in our French region, suggesting no impact of population density or rural environment.

Genetic characterization of a multidrug-resistant Salmonella enterica serovar Agona isolated from a dietary supplement in Germany

Salmonella enterica subsp. enterica serovar Agona has a history of causing food-borne outbreaks and any emergence of multidrug-resistant (MDR) isolates in novel food products is of concern. Particularly, in food products frequently consumed without sufficient heating prior to consumption. Here, we report about the MDR isolate, 18-SA00377, which had been isolated from a dietary supplement in Germany in 2018 and submitted to the German National Reference Laboratory for Salmonella. WGS-based comparative genetic analyses were conducted to find a potential reservoir of the isolate itself or mobile genetic elements associated with MDR. As a phylogenetic analysis did not yield any closely related S. Agona isolates, either globally or from Germany, a detailed analysis of the largest plasmid (295,499 bp) was performed as it is the main carrier of resistances. A combined approach of long-read and short-read sequencing enabled the assembly of the isolate's chromosome and its four plasmids. Their characterization revealed the presence of 23 different antibiotic resistance genes (ARGs), conferring resistance to 12 different antibiotic drug classes, as well as genes conferring resistance to six different heavy metals. The largest plasmid, pSE18-SA00377-1, belongs to the IncHI2 plasmid family and carries 16 ARGs, that are organized as two distinct clusters, with each ARG associated with putative composite transposons. Through a two-pronged approach, highly similar plasmids to pSE18-SA00377-1 were identified in the NCBI database and a search for Salmonella isolates with a highly similar ARG resistance profile was conducted. Mapping and structural comparisons between pSE18-SA00377-1 and these plasmids and Salmonella isolates showed that both the plasmid backbone and identical or similar ARG clusters can be found not only in Salmonella isolates, originating mostly from a wide variety of livestock, but also in a diverse range of bacterial genera of varying geographical origins and isolation sources. Thus, it can be speculated that the host range of pSE18-SA00377-1 is not restricted to Salmonella and its spread already occurred in different bacterial populations. Overall, this hints at a complex history for pSE18-SA00377-1 and highlights the importance of surveilling multidrug-resistant S. enterica isolates, especially in novel food items that are not yet heavily regulated.

Exploring clonality and virulence gene associations in bloodstream infections using whole-genome sequencing and clinical data

Background

Bloodstream infections (BSIs) remain a significant cause of mortality worldwide. Causative pathogens are routinely identified and susceptibility tested but only very rarely investigated for their resistance genes, virulence factors, and clonality. Our aim was to gain insight into the clonality patterns of different species causing BSI and the clinical relevance of distinct virulence genes.

Methods

For this study, we whole-genome-sequenced over 400 randomly selected important pathogens isolated from blood cultures in our diagnostic department between 2016 and 2021. Genomic data on virulence factors, resistance genes, and clonality were cross-linked with in-vitro data and demographic and clinical information.

Results

The investigation yielded extensive and informative data on the distribution of genes implicated in BSI as well as on the clonality of isolates across various species.

Conclusion

Associations between survival outcomes and the presence of specific genes must be interpreted with caution, and conducting replication studies with larger sample sizes for each species appears mandatory. Likewise, a deeper knowledge of virulence and host factors will aid in the interpretation of results and might lead to more targeted therapeutic and preventive measures. Monitoring transmission dynamics more efficiently holds promise to serve as a valuable tool in preventing in particular BSI caused by nosocomial pathogens.

The Parkinson's drug entacapone disrupts gut microbiome homeostasis via iron sequestration

Increasing evidence shows that many human-targeted drugs alter the gut microbiome, leading to implications for host health. However, much less is known about the mechanisms by which drugs target the microbiome and how drugs affect microbial function. Here we combined quantitative microbiome profiling, long-read metagenomics, stable isotope probing and single cell chemical imaging to investigate the impact of two widely prescribed nervous system targeted drugs on the gut microbiome. Ex vivo supplementation of physiologically relevant concentrations of entacapone or loxapine succinate to faecal samples significantly impacted the abundance of up to one third of the microbial species present. Importantly, we demonstrate that the impact of these drugs on microbial metabolism is much more pronounced than their impact on abundances, with low concentrations of drugs reducing the activity, but not the abundance of key microbiome members like Bacteroides, Ruminococcus or Clostridium species. We further demonstrate that entacapone impacts the microbiome due to its ability to complex and deplete available iron, and that microbial growth can be rescued by replenishing levels of microbiota-accessible iron. Remarkably, entacapone-induced iron starvation selected for iron-scavenging organisms carrying antimicrobial resistance and virulence genes. Collectively, our study unveils the impact of two under-investigated drugs on whole microbiomes and identifies metal sequestration as a mechanism of drug-induced microbiome disturbance.

Machine Learning Suggests That Small Size Helps Broaden Plasmid Host Range

Plasmids mediate gene exchange across taxonomic barriers through conjugation, shaping bacterial evolution for billions of years. While plasmid mobility can be harnessed for genetic engineering and drug-delivery applications, rapid plasmid-mediated spread of resistance genes has rendered most clinical antibiotics useless. To solve this urgent and growing problem, we must understand how plasmids spread across bacterial communities. Here, we applied machine-learning models to identify features that are important for extending the plasmid host range. We assembled an up-to-date dataset of more than thirty thousand bacterial plasmids, separated them into 1125 clusters, and assigned each cluster a distribution possibility score, taking into account the host distribution of each taxonomic rank and the sampling bias of the existing sequencing data. Using this score and an optimized plasmid feature pool, we built a model stack consisting of DecisionTreeRegressor, EvoTreeRegressor, and LGBMRegressor as base models and LinearRegressor as a meta-learner. Our mathematical modeling revealed that sequence brevity is the most important determinant for plasmid spread, followed by P-loop NTPases, mobility factors, and β-lactamases. Ours and other recent results suggest that small plasmids may broaden their range by evading host defenses and using alternative modes of transfer instead of autonomous conjugation.

Knowing and Naming: Phage Annotation and Nomenclature for Phage Therapy

Bacteriophages, or phages, are viruses that infect bacteria shaping microbial communities and ecosystems. They have gained attention as potential agents against antibiotic resistance. In phage therapy, lytic phages are preferred for their bacteria killing ability, while temperate phages, which can transfer antibiotic resistance or toxin genes, are avoided. Selection relies on plaque morphology and genome sequencing. This review outlines annotating genomes, identifying critical genomic features, and assigning functional labels to protein-coding sequences. These annotations prevent the transfer of unwanted genes, such as antimicrobial resistance or toxin genes, during phage therapy. Additionally, it covers International Committee on Taxonomy of Viruses (ICTV)-an established phage nomenclature system for simplified classification and communication. Accurate phage genome annotation and nomenclature provide insights into phage-host interactions, replication strategies, and evolution, accelerating our understanding of the diversity and evolution of phages and facilitating the development of phage-based therapies.

Comparison of genomic diversity between single and pooled Staphylococcus aureus colonies isolated from human colonization cultures

The most common approach to sampling the bacterial populations within an infected or colonized host is to sequence genomes from a single colony obtained from a culture plate. However, it is recognized that this method does not capture the genetic diversity in the population. Sequencing a mixture of several colonies (pool-seq) is a better approach to detect population heterogeneity, but it is more complex to analyse due to different types of heterogeneity, such as within-clone polymorphisms, multi-strain mixtures, multi-species mixtures and contamination. Here, we compared 8 single-colony isolates (singles) and pool-seq on a set of 2286 Staphylococcus aureus culture samples to identify features that can distinguish pure samples, samples undergoing intraclonal variation and mixed strain samples. The samples were obtained by swabbing 3 body sites on 85 human participants quarterly for a year, who initially presented with a methicillin-resistant S. aureus skin and soft-tissue infection (SSTI). We compared parameters such as sequence quality, contamination, allele frequency, nucleotide diversity and pangenome diversity in each pool to those for the corresponding singles. Comparing singles from the same culture plate, we found that 18% of sample collections contained mixtures of multiple multilocus sequence types (MLSTs or STs). We showed that pool-seq data alone could predict the presence of multi-ST populations with 95% accuracy. We also showed that pool-seq could be used to estimate the number of intra-clonal polymorphic sites in the population. Additionally, we found that the pool may contain clinically relevant genes such as antimicrobial resistance markers that may be missed when only examining singles. These results highlight the potential advantage of analysing genome sequences of total populations obtained from clinical cultures rather than single colonies.

Systematic analysis of plasmids of the Serratia marcescens complex using 142 closed genomes

Plasmids play important roles in bacterial genome diversification. In the Serratia marcescens complex (SMC), a notable contribution of plasmids to genome diversification was also suggested by our recent analysis of >600 draft genomes. As accurate analyses of plasmids in draft genomes are difficult, in this study we analysed 142 closed genomes covering the entire complex, 67 of which were obtained in this study, and identified 132 plasmids (1.9-244.4 kb in length) in 77 strains. While the average numbers of plasmids in clinical and non-clinical strains showed no significant difference, strains belonging to clade 2 (one of the two hospital-adapted lineages) contained more plasmids than the others. Pangenome analysis revealed that of the 28 954 genes identified, 12.8 % were plasmid-specific, and 1.4 % were present in plasmids or chromosomes depending on the strain. In the latter group, while transposon-related genes were most prevalent (31.4 % of the function-predicted genes), genes related to antimicrobial resistance and heavy metal resistance accounted for a notable proportion (22.7 %). Mash distance-based clustering separated the 132 plasmids into 23 clusters and 50 singletons. Most clusters/singletons showed notably different GC contents compared to those of host chromosomes, suggesting their recent or relatively recent appearance in the SMC. Among the 23 clusters, 17 were found in only clinical or only non-clinical strains, suggesting the possible preference of their distribution on the environmental niches of host strains. Regarding the host strain phylogeny, 16 clusters were distributed in two or more clades, suggesting their interclade transmission. Moreover, for many plasmids, highly homologous plasmids were found in other species, indicating the broadness of their potential host ranges, beyond the genus, family, order, class or even phylum level. Importantly, highly homologous plasmids were most frequently found in Klebsiella pneumoniae and other species in the family Enterobacteriaceae, suggesting that this family, particularly K. pneumoniae, is the main source for plasmid exchanges with the SMC. These results highlight the power of closed genome-based analysis in the investigation of plasmids and provide important insights into the nature of plasmids distributed in the SMC.

Anti-persister efficacy of colistin and meropenem against uropathogenic Escherichia coli is dependent on environmental conditions

Antibiotic persistence is a phenomenon observed when genetically susceptible cells survive long-term exposure to antibiotics. These 'persisters' are an intrinsic component of bacterial populations and stem from phenotypic heterogeneity. Persistence to antibiotics is a concern for public health globally, as it increases treatment duration and can contribute to treatment failure. Furthermore, there is a growing array of evidence that persistence is a 'stepping-stone' for the development of genetic antimicrobial resistance. Urinary tract infections (UTIs) are a major contributor to antibiotic consumption worldwide, and are known to be both persistent (i.e. affecting the host for a prolonged period) and recurring. Currently, in clinical settings, routine laboratory screening of pathogenic isolates does not determine the presence or the frequency of persister cells. Furthermore, the majority of research undertaken on antibiotic persistence has been done on lab-adapted bacterial strains. In the study presented here, we characterized antibiotic persisters in a panel of clinical uropathogenic Escherichia coli isolates collected from hospitals in the UK and Australia. We found that a urine-pH mimicking environment not only induces higher levels of antibiotic persistence to meropenem and colistin than standard laboratory growth conditions, but also results in rapid development of transient colistin resistance, regardless of the genetic resistance profile of the isolate. Furthermore, we provide evidence for the presence of multiple virulence factors involved in stress resistance and biofilm formation in the genomes of these isolates, whose activities have been previously shown to contribute to the formation of persister cells.

Bacterial genome-wide association study substantiates papGII of Escherichia coli as a major risk factor for urosepsis

BACKGROUND

Urinary tract infections (UTIs) are among the most common bacterial infections worldwide, often caused by uropathogenic Escherichia coli. Multiple bacterial virulence factors or patient characteristics have been linked separately to progressive, more invasive infections. In this study, we aim to identify pathogen- and patient-specific factors that drive the progression to urosepsis by jointly analysing bacterial and host characteristics.

METHODS

We analysed 1076 E. coli strains isolated from 825 clinical cases with UTI and/or bacteraemia by whole-genome sequencing (Illumina). Sequence types (STs) were determined via srst2 and capsule loci via fastKaptive. We compared the isolates from urine and blood to confirm clonality. Furthermore, we performed a bacterial genome-wide association study (bGWAS) (pyseer) using bacteraemia as the primary clinical outcome. Clinical data were collected by an electronic patient chart review. We concurrently analysed the association of the most significant bGWAS hit and important patient characteristics with the clinical endpoint bacteraemia using a generalised linear model (GLM). Finally, we designed qPCR primers and probes to detect papGII-positive E. coli strains and prospectively screened E. coli from urine samples (n = 1657) at two healthcare centres.

RESULTS

Our patient cohort had a median age of 75.3 years (range: 18.00-103.1) and was predominantly female (574/825, 69.6%). The bacterial phylogroups B2 (60.6%; 500/825) and D (16.6%; 137/825), which are associated with extraintestinal infections, represent the majority of the strains in our collection, many of which encode a polysaccharide capsule (63.4%; 525/825). The most frequently observed STs were ST131 (12.7%; 105/825), ST69 (11.0%; 91/825), and ST73 (10.2%; 84/825). Of interest, in 12.3% (13/106) of cases, the E. coli pairs in urine and blood were only distantly related. In line with previous bGWAS studies, we identified the gene papGII (p-value < 0.001), which encodes the adhesin subunit of the E. coli P-pilus, to be associated with 'bacteraemia' in our bGWAS. In our GLM, correcting for patient characteristics, papGII remained highly significant (odds ratio = 5.27, 95% confidence interval = [3.48, 7.97], p-value < 0.001). An independent cohort of cases which we screened for papGII-carrying E. coli at two healthcare centres further confirmed the increased relative frequency of papGII-positive strains causing invasive infection, compared to papGII-negative strains (p-value = 0.033, chi-squared test).

CONCLUSIONS

This study builds on previous work linking papGII with invasive infection by showing that it is a major risk factor for progression from UTI to bacteraemia that has diagnostic potential.

Complete genome sequence of an avian pathogenic Escherichia coli strain isolated from poultry

Avian pathogenic Escherichia coli found in the avian intestinal tract can cause systemic disease in birds and act as a foodborne zoonotic pathogen associated with human disease. Here, we report the complete genome sequence of E. coli strain H1998 isolated from a chicken with colisepticemia.

Recovery of 52 bacterial genomes from the fecal microbiome of the domestic cat (Felis catus) using Hi-C proximity ligation and shotgun metagenomics

We used Hi-C proximity ligation with shotgun sequencing to retrieve metagenome-assembled genomes (MAGs) from the fecal microbiomes of two domestic cats (Felis catus). The genomes were assessed for completeness and contamination, classified taxonomically, and annotated for putative antimicrobial resistance (AMR) genes.

Draft genome sequence of biofilm-forming methicillin-resistant Staphylococcus aureus MTR_V1 strain isolated from a ready-to-eat food in Bangladesh

This announcement provides the genome sequence of the biofilm-forming methicillin-resistant Staphylococcus aureus MTR_V1 strain isolated from a ready-to-eat food sample in Bangladesh. Our assembled genome had a length of 2.8 Mb, 27 contigs, two CRISPR arrays, 38 predicted antibiotic resistance genes, and 66 predicted virulence factor genes.

Dispersal history and bidirectional human-fish host switching of invasive, hypervirulent Streptococcus agalactiae sequence type 283

Human group B Streptococcus (GBS) infections attributable to an invasive, hypervirulent sequence type (ST) 283 have been associated with freshwater fish consumption in Asia. The origin, geographic dispersion pathways and host transitions of GBS ST283 remain unresolved. We gather 328 ST283 isolate whole-genome sequences collected from humans and fish between 1998 and 2021, representing eleven countries across four continents. We apply Bayesian phylogeographic analyses to reconstruct the dispersal history of ST283 and combine ST283 phylogenies with genetic markers and host association to investigate host switching and the gain and loss of antimicrobial resistance and virulence factor genes. Initial dispersal within Asia followed ST283 emergence in the early 1980s, with Singapore, Thailand and Hong Kong observed as early transmission hubs. Subsequent intercontinental dispersal originating from Vietnam began in the decade commencing 2001, demonstrating ST283 holds potential to expand geographically. Furthermore, we observe bidirectional host switching, with the detection of more frequent human-to-fish than fish-to-human transitions, suggesting that sound wastewater management, hygiene and sanitation may help to interrupt chains of transmission between hosts. We also show that antimicrobial resistance and virulence factor genes were lost more frequently than gained across the evolutionary history of ST283. Our findings highlight the need for enhanced surveillance, clinical awareness, and targeted risk mitigation to limit transmission and reduce the impact of an emerging pathogen associated with a high-growth aquaculture industry.

ProbioMinServer: an integrated platform for assessing the safety and functional properties of potential probiotic strains

Motivation

ProbioMinServer is a platform designed to help researchers access information on probiotics regarding a wide variety of characteristics, such as safety (e.g. antimicrobial resistance, virulence, pathogenic, plasmid, and prophage genes) and functionality (e.g. functional classes, carbohydrate-active enzyme, and metabolite gene cluster profile). Because probiotics are functional foods, their safety and functionality are a crucial part of health care. Genomics has become a crucial methodology for investigating the safety and functionality of probiotics in food and feed. This shift is primarily attributed to the growing affordability of next-generation sequencing technologies. However, no integrated platform is available for simultaneously evaluating probiotic strain safety, investigating probiotic functionality, and identifying known phylogenetically related strains.

Results

Thus, we constructed a new platform, ProbioMinServer, which incorporates these functions. ProbioMinServer accepts whole-genome sequence files in the FASTA format. If the query genome belongs to the 25 common probiotic species collected in our database, the server performs a database search and analyzes the core-genome multilocus sequence typing. Front-end applications were implemented in JavaScript with a bootstrap framework, and back-end programs were implemented using PHP, Perl, and Python. ProbioMinServer can help researchers quickly and easily retrieve information on the safety and functionality of various probiotics.

Availability and implementation

The platform is available at https://probiomindb.imst.nsysu.edu.tw.

Identification and Characterization of ten Escherichia coli Strains Encoding Novel Shiga Toxin 2 Subtypes, Stx2n as Well as Stx2j, Stx2m, and Stx2o, in the United States

The sharing of genome sequences in online data repositories allows for large scale analyses of specific genes or gene families. This can result in the detection of novel gene subtypes as well as the development of improved detection methods. Here, we used publicly available WGS data to detect a novel Stx subtype, Stx2n in two clinical E. coli strains isolated in the USA. During this process, additional Stx2 subtypes were detected; six Stx2j, one Stx2m strain, and one Stx2o, were all analyzed for variability from the originally described subtypes. Complete genome sequences were assembled from short- or long-read sequencing and analyzed for serotype, and ST types. The WGS data from Stx2n- and Stx2o-producing STEC strains were further analyzed for virulence genes pro-phage analysis and phage insertion sites. Nucleotide and amino acid maximum parsimony trees showed expected clustering of the previously described subtypes and a clear separation of the novel Stx2n subtype. WGS data were used to design OMNI PCR primers for the detection of all known stx1 (283 bp amplicon), stx2 (400 bp amplicon), intimin encoded by eae (221 bp amplicon), and stx2f (438 bp amplicon) subtypes. These primers were tested in three different laboratories, using standard reference strains. An analysis of the complete genome sequence showed variability in serogroup, virulence genes, and ST type, and Stx2 pro-phages showed variability in size, gene composition, and phage insertion sites. The strains with Stx2j, Stx2m, Stx2n, and Stx2o showed toxicity to Vero cells. Stx2j carrying strain, 2012C-4221, was induced when grown with sub-inhibitory concentrations of ciprofloxacin, and toxicity was detected. Taken together, these data highlight the need to reinforce genomic surveillance to identify the emergence of potential new Stx2 or Stx1 variants. The importance of this surveillance has a paramount impact on public health. Per our description in this study, we suggest that 2017C-4317 be designated as the Stx2n type-strain.

The hologenome of Daphnia magna reveals possible DNA methylation and microbiome-mediated evolution of the host genome

Properties that make organisms ideal laboratory models in developmental and medical research are often the ones that also make them less representative of wild relatives. The waterflea Daphnia magna is an exception, by both sharing many properties with established laboratory models and being a keystone species, a sentinel species for assessing water quality, an indicator of environmental change and an established ecotoxicology model. Yet, Daphnia's full potential has not been fully exploited because of the challenges associated with assembling and annotating its gene-rich genome. Here, we present the first hologenome of Daphnia magna, consisting of a chromosomal-level assembly of the D. magna genome and the draft assembly of its metagenome. By sequencing and mapping transcriptomes from exposures to environmental conditions and from developmental morphological landmarks, we expand the previously annotates gene set for this species. We also provide evidence for the potential role of gene-body DNA-methylation as a mutagen mediating genome evolution. For the first time, our study shows that the gut microbes provide resistance to commonly used antibiotics and virulence factors, potentially mediating Daphnia's environmental-driven rapid evolution. Key findings in this study improve our understanding of the contribution of DNA methylation and gut microbiota to genome evolution in response to rapidly changing environments.

MGS2AMR: a gene-centric mining of metagenomic sequencing data for pathogens and their antimicrobial resistance profile

BACKGROUND

Identification of pathogenic bacteria from clinical specimens and evaluating their antimicrobial resistance (AMR) are laborious tasks that involve in vitro cultivation, isolation, and susceptibility testing. Recently, a number of methods have been developed that use machine learning algorithms applied to the whole-genome sequencing data of isolates to approach this problem. However, making AMR assessments from more easily available metagenomic sequencing data remains a big challenge.

RESULTS

We present the Metagenomic Sequencing to Antimicrobial Resistance (MGS2AMR) pipeline, which detects antibiotic resistance genes (ARG) and their possible organism of origin within a sequenced metagenomics sample. This in silico method allows for the evaluation of bacterial AMR directly from clinical specimens, such as stool samples. We have developed two new algorithms to optimize and annotate the genomic assembly paths within the raw Graphical Fragment Assembly (GFA): the GFA Linear Optimal Path through seed segments (GLOPS) algorithm and the Adapted Dijkstra Algorithm for GFA (ADAG). These novel algorithms improve the sensitivity of ARG detection and aid in species annotation. Tests based on 1200 microbiome samples show a high ARG recall rate and correct assignment of the ARG origin. The MGS2AMR output can further be used in many downstream applications, such as evaluating AMR to specific antibiotics in samples from emerging intestinal infections. We demonstrate that the MGS2AMR-derived data is as informative for the entailing prediction models as the whole-genome sequencing (WGS) data. The performance of these models is on par with our previously published method (WGS2AMR), which is based on the sequencing data of bacterial isolates.

CONCLUSIONS

MGS2AMR can provide researchers with valuable insights into the AMR content of microbiome environments and may potentially improve patient care by providing faster quantification of resistance against specific antibiotics, thereby reducing the use of broad-spectrum antibiotics. The presented pipeline also has potential applications in other metagenome analyses focused on the defined sets of genes. Video Abstract.

Chance Favors the Prepared Genomes: Horizontal Transfer Shapes the Emergence of Antibiotic Resistance Mutations in Core Genes

Bacterial lineages acquire novel traits at diverse rates in part because the genetic background impacts the successful acquisition of novel genes by horizontal transfer. Yet, how horizontal transfer affects the subsequent evolution of core genes remains poorly understood. Here, we studied the evolution of resistance to quinolones in Escherichia coli accounting for population structure. We found 60 groups of genes whose gain or loss induced an increase in the probability of subsequently becoming resistant to quinolones by point mutations in the gyrase and topoisomerase genes. These groups include functions known to be associated with direct mitigation of the effect of quinolones, with metal uptake, cell growth inhibition, biofilm formation, and sugar metabolism. Many of them are encoded in phages or plasmids. Although some of the chronologies may reflect epidemiological trends, many of these groups encoded functions providing latent phenotypes of antibiotic low-level resistance, tolerance, or persistence under quinolone treatment. The mutations providing resistance were frequent and accumulated very quickly. Their emergence was found to increase the rate of acquisition of other antibiotic resistances setting the path for multidrug resistance. Hence, our findings show that horizontal gene transfer shapes the subsequent emergence of adaptive mutations in core genes. In turn, these mutations further affect the subsequent evolution of resistance by horizontal gene transfer. Given the substantial gene flow within bacterial genomes, interactions between horizontal transfer and point mutations in core genes may be a key to the success of adaptation processes.

High-level ceftazidime/avibactam resistance in Escherichia coli conferred by the novel plasmid-mediated β-lactamase CMY-185 variant

OBJECTIVES

To characterize a blaCMY variant associated with ceftazidime/avibactam resistance from a serially collected Escherichia coli isolate.

METHODS

A patient with an intra-abdominal infection due to recurrent E. coli was treated with ceftazidime/avibactam. On Day 48 of ceftazidime/avibactam therapy, E. coli with a ceftazidime/avibactam MIC of >256 mg/L was identified from abdominal drainage. Illumina and Oxford Nanopore Technologies WGS was performed on serial isolates to identify potential resistance mechanisms. Site-directed mutants of CMY β-lactamase were constructed to identify amino acid residues responsible for ceftazidime/avibactam resistance.

RESULTS

WGS revealed that all three isolates were E. coli ST410. The ceftazidime/avibactam-resistant strain uniquely acquired a novel CMY β-lactamase gene, herein called blaCMY-185, harboured on an IncI-γ/K1 conjugative plasmid. The CMY-185 enzyme possessed four amino acid substitutions relative to CMY-2, including A114E, Q120K, V211S and N346Y, and conferred high-level ceftazidime/avibactam resistance with an MIC of 32 mg/L. Single CMY-2 mutants did not confer reduced ceftazidime/avibactam susceptibility. However, double and triple mutants containing N346Y previously associated with ceftazidime/avibactam resistance in other AmpC enzymes, conferred ceftazidime/avibactam MICs ranging between 4 and 32 mg/L as well as reduced susceptibility to the newly developed cephalosporin, cefiderocol. Molecular modelling suggested that the N346Y substitution confers the reduction of avibactam inhibition due to steric hindrance between the side chain of Y346 and the sulphate group of avibactam.

CONCLUSIONS

We identified ceftazidime/avibactam resistance in E. coli associated with a novel CMY variant. Unlike other AmpC enzymes, CMY-185 appears to require an additional substitution on top of N346Y to confer ceftazidime/avibactam resistance.

SHIP: identifying antimicrobial resistance gene transfer between plasmids

MOTIVATION

Plasmids are carriers for antimicrobial resistance (AMR) genes and can exchange genetic material with other structures, contributing to the spread of AMR. There is no reliable approach to identify the transfer of AMR genes across plasmids. This is mainly due to the absence of a method to assess the phylogenetic distance of plasmids, as they show large DNA sequence variability. Identifying and quantifying such transfer can provide novel insight into the role of small mobile elements and resistant plasmid regions in the spread of AMR.

RESULTS

We developed SHIP, a novel method to quantify plasmid similarity based on the dynamics of plasmid evolution. This allowed us to find conserved fragments containing AMR genes in structurally different and phylogenetically distant plasmids, which is evidence for lateral transfer. Our results show that regions carrying AMR genes are highly mobilizable between plasmids through transposons, integrons, and recombination events, and contribute to the spread of AMR. Identified transferred fragments include a multi-resistant complex class 1 integron in Escherichia coli and Klebsiella pneumoniae, and a region encoding tetracycline resistance transferred through recombination in Enterococcus faecalis.

AVAILABILITY AND IMPLEMENTATION

The code developed in this work is available at https://github.com/AbeelLab/plasmidHGT.

Evidence of potential Campylobacter jejuni zooanthroponosis in captive macaque populations

Non-human primates share recent common ancestry with humans and exhibit comparable disease symptoms. Here, we explored the transmission potential of enteric bacterial pathogens in monkeys exhibiting symptoms of recurrent diarrhoea in a biomedical research facility in China. The common zoonotic bacterium Campylobacter jejuni was isolated from macaques (Macaca mulatta and Macaca fascicularis) and compared to isolates from humans and agricultural animals in Asia. Among the monkeys sampled, 5 % (44/973) tested positive for C. jejuni, 11 % (5/44) of which displayed diarrhoeal symptoms. Genomic analysis of monkey isolates, and 1254 genomes from various sources in Asia, were used to identify the most likely source of human infection. Monkey and human isolates shared high average nucleotide identity, common MLST clonal complexes and clustered together on a phylogeny. Furthermore, the profiles of putative antimicrobial resistance genes were similar between monkeys and humans. Taken together these findings suggest that housed macaques became infected with C. jejuni either directly from humans or via a common contamination source.

Flavobacterium psychrotrophum sp. nov. and Flavobacterium panacagri sp. nov., Isolated from Freshwater and Soil

Two novel bacterial strains CJ74T and CJ75T belonging to the genus Flavobacterium were isolated from freshwater of Han River and ginseng soil, South Korea, respectively. Strain CJ74T was Gram-stain-negative, aerobic, rod-shaped, non-motile, and non-flagellated, and did not produce flexirubin-type pigments. Strain CJ75T was Gram-stain-negative, aerobic, rod-shaped, motile by gliding, and non-flagellated, and produced flexirubin-type pigments. Both strains were shown to grow optimally at 30 °C in the absence of NaCl on R2A medium. Phylogenetic analysis based on 16S rRNA gene sequences showed that strains CJ74T and CJ75T belonged to the genus Flavobacterium and were most closely related to Flavobacterium niveum TAPW14T and Flavobacterium foetidum CJ42T with 96.17% and 97.29% 16S rRNA sequence similarities, respectively. Genomic analyses including the reconstruction of phylogenomic tree, average nucleotide identity, and digital DNA-DNA hybridization suggested that they were novel species of the genus Flavobacterium. Both strains contained menaquinone 6 (MK-6) as the primary respiratory quinone and phosphatidylethanolamine as a major polar lipid. The predominant fatty acids of both strains were iso-C15:0 and summed feature 3 (C16:1 ω7c and/or C16:1 ω6c). Based on the polyphasic taxonomic study, strains CJ74T and CJ75T represent novel species of the genus Flavobacterium, for which names Flavobacterium psychrotrophum sp. nov. and Flavobacterium panacagri sp. nov. are proposed, respectively. The type strains are CJ74T (=KACC 19819T =JCM 32889T) and CJ75T (=KACC 23149T =JCM 36132T).

Prevalence and Characteristics of Invasive Staphylococcus argenteus among Patients with Bacteremia in Hong Kong

Staphylococcus argenteus is a novel Staphylococcus species derived from Staphylococcus aureus. Information on the prevalence and genetic characteristics of invasive S. argenteus in Asia is limited. In this study, 275 invasive S. aureus complex strains were retrieved from blood culture specimens in Hong Kong and re-analyzed using MALDI-TOF mass spectrometry and an in-house multiplex real-time PCR for S. argenteus. The prevalence of invasive S. argenteus in Hong Kong was found to be 4.0% (11/275). These strains were primarily susceptible to commonly used antibiotics, except penicillin. Whole-genome sequencing revealed the circulation of three S. argenteus genotypes (ST-2250, ST-1223, and ST-2854) in Hong Kong, with ST-2250 and ST-1223 being the predominant genotypes. The local ST-2250 and ST-1223 strains showed close phylogenetic relationships with isolates from mainland China. Antimicrobial-resistant genes (fosB, tet-38, mepA, blaI, blaZ) could be found in nearly all local S. argenteus strains. The ST-1223 and ST-2250 genotypes carried multiple staphylococcal enterotoxin genes that could cause food poisoning and toxic shock syndrome. The CRISPR/Cas locus was observed only in the ST-2250 strains. This study provides the first report on the molecular epidemiology of invasive S. argenteus in Hong Kong, and further analysis is needed to understand its transmission reservoir.

Towards interoperability in infection control: a standard data model for microbiology

The COVID-19 pandemic has made it clear: sharing and exchanging data among research institutions is crucial in order to efficiently respond to global health threats. This can be facilitated by defining health data models based on interoperability standards. In Germany, a national effort is in progress to create common data models using international healthcare IT standards. In this context, collaborative work on a data set module for microbiology is of particular importance as the WHO has declared antimicrobial resistance one of the top global public health threats that humanity is facing. In this article, we describe how we developed a common model for microbiology data in an interdisciplinary collaborative effort and how we make use of the standard HL7 FHIR and terminologies such as SNOMED CT or LOINC to ensure syntactic and semantic interoperability. The use of international healthcare standards qualifies our data model to be adopted beyond the environment where it was first developed and used at an international level.

Complete genome sequences of bacteria isolated from cockroaches collected in a Bangladeshi hospital

We report the complete genome sequences of four bacterial strains that were isolated from Blattella germanica (German cockroaches) that were found in three wards of the Rajshahi Medical College Hospital. Multiple antibiotic resistance genes were identified in each genome, with one genome containing multiple plasmid-encoded resistance genes.

Draft genome sequence of biofilm-forming Enterococcus faecalis BAU_Ef01 strain isolated from shrimp (Penaeus indicus) in Bangladesh

Here, we sequence and analyze a biofilm-forming strain of Enterococcus faecalis BAU_Ef01 isolated from a shrimp in Bangladesh. The whole genome of the strain had a length of 2,862,301 bp, 38 contigs, an average G+C content of 37.36%, 80.0× genome coverage, and 35 predicted antibiotic resistance and virulence genes each.

Enterobacterales draft genome sequences: 15 historical (1998-2004) and 30 contemporary (2015-2016) clinical isolates from Pakistan

The continued emergence and spread of antimicrobial resistance among pathogenic bacteria are ever-growing threats to health and economy. Here, we report the draft genomes for 45 Enterobacterales clinical isolates, including historical and contemporary drug-resistant organisms, obtained in Pakistan between 1998 and 2016: 5 Serratia, 3 Salmonella, 3 Enterobacter, and 34 Klebsiella.

Antimicrobial resistance, pathogenic potential, and genomic features of carbapenem-resistant Klebsiella pneumoniae isolated in Chile: high-risk ST25 clones and novel mobile elements

Multidrug- and carbapenem-resistant Klebsiella pneumoniae (CR-Kp) are critical threats to global health and key traffickers of resistance genes to other pathogens. Despite the sustained increase in CR-Kp infections in Chile, few strains have been described at the genomic level, lacking details of their resistance and virulence determinants and the mobile elements mediating their dissemination. In this work, we studied the antimicrobial susceptibility and performed a comparative genomic analysis of 10 CR-Kp isolates from the Chilean surveillance of carbapenem-resistant Enterobacteriaceae. High resistance was observed among the isolates (five ST25, three ST11, one ST45, and one ST505), which harbored 44 plasmids, most carrying genes for conjugation and resistance to several antibiotics and biocides. Ten plasmids encoding carbapenemases were characterized, including novel plasmids or variants with additional resistance genes, a novel genetic environment for blaKPC-2, and plasmids widely disseminated in South America. ST25 K2 isolates belonging to CG10224, a clone traced back to 2012 in Chile, which recently acquired blaNDM-1, blaNDM-7, or blaKPC-2 plasmids stood out as high-risk clones. Moreover, this corresponds to the first report of ST25 and ST45 Kp producing NDM-7 in South America and ST505 CR-Kp producing both NDM-7 and KPC-2 worldwide. Also, we characterized a variety of genomic islands carrying virulence and fitness factors. These results provide baseline knowledge for a detailed understanding of molecular and genetic determinants behind antibiotic resistance and virulence of CR-Kp in Chile and South America. IMPORTANCE In the ongoing antimicrobial resistance crisis, carbapenem-resistant strains of Klebsiella pneumoniae are critical threats to public health. Besides globally disseminated clones, the burden of local problem clones remains substantial. Although genomic analysis is a powerful tool for improving pathogen and antimicrobial resistance surveillance, it is still restricted in low- to middle-income countries, including Chile, causing them to be underrepresented in genomic databases and epidemiology surveys. This study provided the first 10 complete genomes of the Chilean surveillance for carbapenem-resistant K. pneumoniae in healthcare settings, unveiling their resistance and virulence determinants and the mobile genetic elements mediating their dissemination, placed in the South American and global K. pneumoniae epidemiological context. We found ST25 with K2 capsule as an emerging high-risk clone, along with other lineages producing two carbapenemases and several other resistance and virulence genes encoded in novel plasmids and genomic islands.

Acinetobacter baumannii IC2 and IC5 Isolates with Co-Existing blaOXA-143-like and blaOXA-72 and Exhibiting Strong Biofilm Formation in a Mexican Hospital

Acinetobacter baumannii is an opportunistic pathogen responsible for healthcare-associated infections (HAIs) and outbreaks. Antimicrobial resistance mechanisms and virulence factors allow it to survive and spread in the hospital environment. However, the molecular mechanisms of these traits and their association with international clones are frequently unknown in low- and middle-income countries. Here, we analyze the phenotype and genotype of seventy-six HAIs and outbreak-causing A. baumannii isolates from a Mexican hospital over ten years, with special attention to the carbapenem resistome and biofilm formation. The isolates belonged to the global international clone (IC) 2 and the Latin America endemic IC5 and were predominantly extensively drug-resistant (XDR). Oxacillinases were identified as a common source of carbapenem resistance. We noted the presence of the blaOXA-143-like family (not previously described in Mexico), the blaOXA-72 and the blaOXA-398 found in both ICs. A low prevalence of efflux pump overexpression activity associated with carbapenem resistance was observed. Finally, strong biofilm formation was found, and significant biofilm-related genes were identified, including bfmRS, csuA/BABCDE, pgaABCD and ompA. This study provides a comprehensive profile of the carbapenem resistome of A. baumannii isolates belonging to the same pulse type, along with their significant biofilm formation capacity. Furthermore, it contributes to a better understanding of their role in the recurrence of infection and the endemicity of these isolates in a Mexican hospital.

Genomic Characterization of IMP-Producing Pseudomonas aeruginosa in Bulgaria Reveals the Emergence of IMP-100, a Novel Plasmid-Mediated Variant Coexisting with a Chromosomal VIM-4

Multidrug-resistant (MDR) Pseudomonas aeruginosa infections represent a major public health concern and require comprehensive understanding of their genetic makeup. This study investigated the first occurrence of imipenemase (IMP)-carrying P. aeruginosa strains from Bulgaria. Whole genome sequencing identified a novel plasmid-mediated IMP-100 allele located in a a novel In4886 integron embedded in a putative Tn7700 transposon. Two other closely related chromosomal IMP variants, IMP-13 and IMP-84, were also detected. The IMP-producers were resistant to last-line drugs including cefiderocol (CFDC) (two out of three) and susceptible to colistin. The IMP-13/84 cassettes were situated in a In320 integron inserted in a Tn5051-like transposon as previously reported. Lastly, the p4782-IMP plasmid rendered the PA01 transformant resistant to CFDC, suggesting a transferable CFDC resistance. A variety of virulence factors associated with adhesion, antiphagocytosis, iron uptake, and quorum sensing, as well as secretion systems, toxins, and proteases, were confirmed, suggesting significant pathogenic potential consistent with the observed strong biofilm formation. The emergence of IMP-producing MDR P. aeruginosa is alarming as it remains unsusceptible even to last-generation drugs like CFDC. Newly detected IMP-100 was even located in a CFDC-resistant XDR strain.

Shotgun-metagenomics reveals a highly diverse and communal microbial network present in the drains of three beef-processing plants

Background

Multi-species biofilms pose a problem in various environments, especially food-processing environments. The diversity of microorganisms in these biofilms plays a critical role in their integrity and protection against external biotic and abiotic factors. Compared to single-species biofilms, mixed-species biofilms are more resistant to various stresses, including antimicrobials like sanitizers. Therefore, understanding the microbiome composition and diversity in biofilms and their metabolic potential is a priority when developing intervention techniques to combat foodborne pathogens in food processing environments.

Methods

This study aimed to describe and compare the microbiome profile of 75 drain biofilm samples obtained from five different locations (Hotscale, Hotbox, Cooler, Processing, & Grind room) of three beef-processing plants (Plant A, B & C) taken over two timepoints 2017-18 (T1) and 2021 (T2) by shotgun sequencing.

Results

Core microbiome analysis found Pseudomonas, Psychrobacter, and Acinetobacter to be the top three prevalent genera among the plants and locations. Alpha diversity analysis demonstrated a high diversity of microbiome present in all the plants and locations across the time points. Functional analysis showed the high metabolic potential of the microbial community with abundance of genes in metabolism, cell-adhesion, motility, and quorum sensing. Moreover, Quaternary Ammonium Compound (QAC) resistance genes were also observed, this is significant as QAC sanitizers are commonly used in many food processing facilities. Multi-functional genes such as transposases, polymerases, permeases, flagellar proteins, and Mobile Genetic Elements (MGEs) were found suggesting these are dynamic microbial communities that work together to protect themselves against environmental stresses through multiple defense mechanisms.

Conclusion

This study provides a framework for understanding the collective microbial network spanning a beef processing system. The results can be used to develop intervention strategies to best control these highly communicative microbial networks.

Phylogenetic lineages and antimicrobial resistance determinants of clinical Klebsiella oxytoca spanning local to global scales

Klebsiella oxytoca is an opportunistic pathogen causing serious nosocomial infections. Knowledge about the population structure and diversity of healthcare-associated K. oxytoca from a genomic standpoint remains limited. Here, we characterized the phylogenetic relationships and genomic characteristics of 20 K. oxytoca sensu stricto isolates recovered from bloodstream infections at the Dartmouth-Hitchcock Medical Center, New Hampshire, USA from 2017 to 2021. Results revealed a diverse population consisting of 15 sequence types (STs) that together harbored 10 variants of the intrinsic beta-lactamase gene bla OXY-2, conferring resistance to penicillins. Similar sets of antimicrobial resistance (AMR) determinants reside in multiple distinct lineages, with no one lineage dominating the local population. To place the New Hampshire K. oxytoca in a broader context, we compared them to 304 publicly available genomes of clinical isolates from 18 countries. This global clinical K. oxytoca sensu stricto population is represented by over 65 STs that together harbored resistance genes against 14 antimicrobial classes, including eight bla OXY-2 variants. Three dominant STs in the global population (ST2, ST176, ST199) circulate across multiple countries and were also present in the New Hampshire population. The global K. oxytoca population is genetically diverse, but there is evidence for broad dissemination of a few lineages carrying distinct set of AMR determinants. Our findings reveal the clinical diversity of K. oxytoca sensu stricto and its importance in surveillance efforts aimed at monitoring the evolution of this drug-resistant nosocomial pathogen. IMPORTANCE The opportunistic pathogen Klebsiella oxytoca has been increasingly implicated in patient morbidity and mortality worldwide, including several outbreaks in healthcare settings. The emergence and spread of antimicrobial resistant strains exacerbate the disease burden caused by this species. Our study showed that clinical K. oxytoca sensu stricto is phylogenetically diverse, harboring various antimicrobial resistance determinants and bla OXY-2 variants. Understanding the genomic and population structure of K. oxytoca is important for international initiatives and local epidemiological efforts for surveillance and control of drug-resistant K. oxytoca.

Serratia marcescens Outbreak at a Correctional Facility: Environmental Sampling, Laboratory Analyses and Genomic Characterization to Assess Sources and Persistence

Serratia marcescens is an environmental bacterium and clinical pathogen that can cause an array of infections. We describe an environmental sampling and comparative genomics approach used to investigate a multi-year outbreak of S. marcescens at a correctional facility. Whole genome sequencing analysis revealed a predominant cluster of clonally related S. marcescens from nine patient cases and items associated with illicit drug use. Closely related strains found among items associated with case-patient cells and diluted Cell Block 64 (CB64), a quaternary ammonium disinfectant, and Break Out (BO), a multipurpose cleaner, highlighted their role as environmental reservoirs for S. marcescens in this outbreak. Comparative genomic analysis suggested outbreak strains were both persistent (identical strains found over long periods and in multiple locations of the correctional facility) and diverse (strains clustered with multiple global samples from NCBI database). No correlation was found between antimicrobial resistance (AMR) genes of outbreak strains; NCBI strains have more AMR genes. Principal component analysis (PCA) of virulence factors associated with persistence and infectivity indicated variation based on phylogroups, including the predominant cluster; identifiable variations among environmental versus clinical strains were not observed. Identification of multiple distinct genetic groups highlights the importance of putting epidemiological genomic studies in a proper genetic context.

P-DOR, an easy-to-use pipeline to reconstruct bacterial outbreaks using genomics

SUMMARY

Bacterial Healthcare-Associated Infections (HAIs) are a major threat worldwide, which can be counteracted by establishing effective infection control measures, guided by constant surveillance and timely epidemiological investigations. Genomics is crucial in modern epidemiology but lacks standard methods and user-friendly software, accessible to users without a strong bioinformatics proficiency. To overcome these issues we developed P-DOR, a novel tool for rapid bacterial outbreak characterization. P-DOR accepts genome assemblies as input, it automatically selects a background of publicly available genomes using k-mer distances and adds it to the analysis dataset before inferring a Single-Nucleotide Polymorphism (SNP)-based phylogeny. Epidemiological clusters are identified considering the phylogenetic tree topology and SNP distances. By analyzing the SNP-distance distribution, the user can gauge the correct threshold. Patient metadata can be inputted as well, to provide a spatio-temporal representation of the outbreak. The entire pipeline is fast and scalable and can be also run on low-end computers.

AVAILABILITY AND IMPLEMENTATION

P-DOR is implemented in Python3 and R and can be installed using conda environments. It is available from GitHub https://github.com/SteMIDIfactory/P-DOR under the GPL-3.0 license.

Stress Resistance and Virulence Gene Profiles Associated with Phylogeny and Phenotypes of Escherichia coli from Cattle

Seven serogroups of E. coli (Top seven E. coli) are frequently implicated in foodborne outbreaks in North America, largely due to their carriage of Shiga toxin genes (stx). This study aimed to profile resistance genes and virulence factors (VF), and their potential association with phylogeny and phenotypes of Top seven E. coli originating from cattle in Canada. 155 Top seven E. coli isolates previously characterized for heat and acid resistance and biofilm-forming ability were whole-genome sequenced and analyzed for phylogeny, VF, and stress resistance genes. The 155 E. coli strains belonged to six phylogroups: A (n = 32), B1 (n = 93), C (n = 3), D (n = 11), E (n = 15), and G (n = 1). Different phylogroups were clearly separated on the core genome tree, with strains of the same serotype closely clustered. The carriage of stx and the transmissible locus of stress tolerance (tLST), the extreme heat resistance marker, was mutually exclusive, in 33 and 15 genomes, respectively. A novel O84:H2 strain carrying stx1a was also identified. In total, 70, 41, and 32 VF, stress resistance genes and antibiotic resistance genes were identified. The stress resistance genes included those for metal (n = 29), biocides/acid (n = 4), and heat (n = 8) resistance. All heat resistance genes and most metal-resistance genes that were differentially distributed among the phylogroups were exclusively in phylogroup A. VF were least and most present in phylogroups A and D, respectively. No specific genes associated with acid resistance or biofilm formation phenotypes were identified. VF were more abundant (P < 0.05) in the non-biofilm-forming population and acid-resistant population.

Identification of combinatorial mutations associated with colistin resistance in Shewanella algae

Colistin is a last-resort antibiotic used to treat infections caused by drug-resistant gram-negative bacteria. However, the genetic mechanisms underlying colistin resistance in Shewanella algae are not well understood. In this study, we sequenced and compared the genomes of 23 mcr-negative colistin-resistant and sensitive S. algae samples from various sources. We applied a computational approach to identify combinatorial mutations associated with colistin resistance. Our analysis revealed a combination of three mutations (PmrB 451, PmrE168, PmrH292) that were strongly associated with colistin resistance in S. algae. This study provides insights into the genetic mechanisms of colistin resistance in S. algae and demonstrates the utility of a computational approach for identifying epistatic interactions among mutations. Identifying the genetic mutations responsible for colistin resistance in S. algae can inform the development of new treatments or strategies to combat infections caused by this emerging pathogen.

Mge-cluster: a reference-free approach for typing bacterial plasmids

Extrachromosomal elements of bacterial cells such as plasmids are notorious for their importance in evolution and adaptation to changing ecology. However, high-resolution population-wide analysis of plasmids has only become accessible recently with the advent of scalable long-read sequencing technology. Current typing methods for the classification of plasmids remain limited in their scope which motivated us to develop a computationally efficient approach to simultaneously recognize novel types and classify plasmids into previously identified groups. Here, we introduce mge-cluster that can easily handle thousands of input sequences which are compressed using a unitig representation in a de Bruijn graph. Our approach offers a faster runtime than existing algorithms, with moderate memory usage, and enables an intuitive visualization, classification and clustering scheme that users can explore interactively within a single framework. Mge-cluster platform for plasmid analysis can be easily distributed and replicated, enabling a consistent labelling of plasmids across past, present, and future sequence collections. We underscore the advantages of our approach by analysing a population-wide plasmid data set obtained from the opportunistic pathogen Escherichia coli, studying the prevalence of the colistin resistance gene mcr-1.1 within the plasmid population, and describing an instance of resistance plasmid transmission within a hospital environment.

Evidence of host specificity in Lactobacillus johnsonii genomes and its influence on probiotic potential in poultry

To date, the selection of candidate strains for probiotic development in production animals has been largely based upon screens for desired phenotypic traits. However, increasing evidence indicates that the use of host-specific strains may be important, because coevolution with the animal host better prepares a bacterial strain to colonize and succeed in its respective host animal species. This concept was applied to Lactobacillus johnsonii in commercial poultry production because of its previous correlation with enhanced bird performance. Using 204 naturally isolated chicken- and turkey-source L. johnsonii, we demonstrate that there is a strong phylogenetic signal for coevolution with the animal host. These isolates differ phenotypically, even within host source, and these differences can be correlated with certain L. johnsonii phylogenetic clades. In commercial turkey poults, turkey-specific strains with strong in vitro phenotypes performed better early in life than strains lacking those phenotypes. A follow-up performance trial in broiler chickens demonstrated that chicken-specific strains result in better overall bird performance than nonchicken-specific strains. Collectively, this work provides evidence for the impact of host adaptation on a probiotic strain's potential. Furthermore, this top-down approach is useful for screening larger numbers of isolates for probiotic candidates.

Current concepts, advances, and challenges in deciphering the human microbiota with metatranscriptomics

Metatranscriptomics refers to the analysis of the collective microbial transcriptome of a sample. Its increased utilization for the characterization of human-associated microbial communities has enabled the discovery of many disease-state related microbial activities. Here, we review the principles of metatranscriptomics-based analysis of human-associated microbial samples. We describe strengths and weaknesses of popular sample preparation, sequencing, and bioinformatics approaches and summarize strategies for their use. We then discuss how human-associated microbial communities have recently been examined and how their characterization may change. We conclude that metatranscriptomics insights into human microbiotas under health and disease have not only expanded our knowledge on human health, but also opened avenues for rational antimicrobial drug use and disease management.

An Overview of Antimicrobial Resistance Profiles of Publicly Available Salmonella Genomes with Sufficient Quality and Metadata

Salmonella enterica (S. enterica) is a commensal organism or pathogen causing diseases in animals and humans, as well as widespread in the environment. Antimicrobial resistance (AMR) has increasingly affected both animal and human health and continues to raise public health concerns. A decade ago, it was estimated that the increased use of whole genome sequencing (WGS) combined with sharing of public data would drastically change and improve the surveillance and understanding of Salmonella epidemiology and AMR. This study aimed to evaluate the current usefulness of public WGS data for Salmonella surveillance and to investigate the associations between serovars, antibiotic resistance genes (ARGs), and metadata. Out of 191,306 Salmonella genomes deposited in European Nucleotide Archive and NCBI databases, 47,452 WGS with sufficient minimum metadata (country, year, and source) of S. enterica were retrieved from 116 countries and isolated between 1905 and 2020. For in silico analysis of the WGS data, KmerFinder, SISTR, and ResFinder were used for species, serovars, and AMR identification, respectively. The results showed that the five common isolation sources of S. enterica are human (29.10%), avian (22.50%), environment (11.89%), water (9.33%), and swine (6.62%). The most common ARG profiles for each class of antimicrobials are β-lactam (blaTEM-1B; 6.78%), fluoroquinolone [(parC[T57S], qnrB19); 0.87%], folate pathway antagonist (sul2; 8.35%), macrolide [mph(A); 0.39%], phenicol (floR; 5.94%), polymyxin B (mcr-1.1; 0.09%), and tetracycline [tet(A); 12.95%]. Our study reports the first overview of ARG profiles in publicly available Salmonella genomes from online databases. All data sets from this study can be searched at Microreact.

Intensive care unit sinks are persistently colonized with multidrug resistant bacteria and mobilizable, resistance-conferring plasmids

Contamination of hospital sinks with microbial pathogens presents a serious potential threat to patients, but our understanding of sink colonization dynamics is largely based on infection outbreaks. Here, we investigate the colonization patterns of multidrug-resistant organisms (MDROs) in intensive care unit sinks and water from two hospitals in the USA and Pakistan collected over 27 months of prospective sampling. Using culture-based methods, we recovered 822 bacterial isolates representing 104 unique species and genomospecies. Genomic analyses revealed long-term colonization by Pseudomonas spp. and Serratia marcescens strains across multiple rooms. Nanopore sequencing uncovered examples of long-term persistence of resistance-conferring plasmids in unrelated hosts. These data indicate that antibiotic resistance (AR) in Pseudomonas spp. is maintained both by strain colonization and horizontal gene transfer (HGT), while HGT maintains AR within Acinetobacter spp. and Enterobacterales, independent of colonization. These results emphasize the importance of proactive, genomic-focused surveillance of built environments to mitigate MDRO spread. IMPORTANCE Hospital sinks are frequently linked to outbreaks of antibiotic-resistant bacteria. Here, we used whole-genome sequencing to track the long-term colonization patterns in intensive care unit (ICU) sinks and water from two hospitals in the USA and Pakistan collected over 27 months of prospective sampling. We analyzed 822 bacterial genomes, representing over 100 different species. We identified long-term contamination by opportunistic pathogens, as well as transient appearance of other common pathogens. We found that bacteria recovered from the ICU had more antibiotic resistance genes (ARGs) in their genomes compared to matched community spaces. We also found that many of these ARGs are harbored on mobilizable plasmids, which were found shared in the genomes of unrelated bacteria. Overall, this study provides an in-depth view of contamination patterns for common nosocomial pathogens and identifies specific targets for surveillance.

Genomic Characterization of a Tetracycline-Resistant Strain of Brochothrix thermosphacta Highlights Plasmids Partially Shared between Various Strains

The Gram-positive bacterium Brochothrix thermosphacta is a spoilage agent commonly found on meat products. While the tet(L) gene, which confers resistance to tetracycline, has been identified in certain strains of B. thermosphacta, only a limited number of studies have investigated this gene and its potential presence on mobile DNA elements. This study aims to analyze the tetracycline-resistant strain B. thermosphacta BT469 at the genomic level to gain insight into the molecular determinants responsible for this resistance. Three plasmids have been identified in the strain: pBT469-1, which contains a tetR gene; pBT469-2, which harbours the tet(L) gene responsible for tetracycline resistance; and pBT469-3, which carries genes encoding for a thioredoxin and a phospholipase A2. Homology searches among sequences in public databases have revealed that the plasmid pBT469-2 is currently unique to the BT469 strain. However, the pBT469-1 plasmid is also found in three other strains of B. thermosphacta. Notably, sequences similar to pBT469-1 and pBT469-2 were also found in other bacterial genera, suggesting that these plasmids may be part of a diverse family present in several bacterial genera. Interestingly, sequences of various strains of B. thermosphacta show a high level of similarity with pBT469-3, suggesting that variants of this plasmid could be frequently found in this bacterium.

First Description of Mergibacter septicus Isolated from a Common Tern (Sterna hirundo) in Germany

Mergibacter septicus (M. septicus), previously known as Bisgaard Taxon 40, is a recently described species within the Pasteurellaceae family. In this study, we present a M. septicus strain isolated from a common tern (Sterna hirundo) chick that died just after fledging from the Banter See in Wilhelmshaven, Germany. The recovered M. septicus strain underwent microbiological phenotypic characterization, followed by whole genome sequencing on Illumina and Nanopore platforms. Phenotypically, M. septicus 19Y0039 demonstrated resistance to colistin, cephalexin, clindamycin, oxacillin, and penicillin G. The genome analysis revealed a circular 1.8 Mbp chromosome without any extrachromosomal elements, containing 1690 coding DNA sequences. The majority of these coding genes were associated with translation, ribosomal structure and biogenesis, followed by RNA processing and modification, and transcription. Genetic analyses revealed that the German M. septicus strain 19Y0039 is related to the American strain M. septicus A25201T. Through BLAST alignment, twelve putative virulence genes previously identified in the M. septicus type strain A25201T were also found in the German strain. Additionally, 84 putative virulence genes distributed across nine categories, including immune modulation, effector delivery system, nutrition/metabolic factors, regulation, stress survival, adherence, biofilm, exotoxin, and motility, were also identified.

Comparative Genomics Identified PenR E151V Substitution Associated with Carbapenem-Resistance Burkholderia cepacia Complex and a Novel Burkholderia cepacia Complex Specific OXA-1043 Subgroup

Purpose

Burkholderia cepacia complex (Bcc) is a known significant opportunistic pathogen causing morbidity and mortality, particularly in those with cystic fibrosis, chronic granulomatous disease, or immunocompromising host. Mortality of Bcc bloodstream infections among non-cystic fibrosis patients remained high. The antibiotic treatment for Bcc infection is quite challenging due to its intrinsic resistance to most antibiotics, and the resistance to carbapenems was the biggest concern among them. We aimed to realize the mechanism of carbapenem resistance in Bcc.

Patients and Methods

Ten strains of Bcc were identified by the MALDI-TOF MS, and the drug susceptibility test was using VITEK 2 system. The Burkholderia cepacia complex genomes were sequenced via Nanopore GridIon. We also downloaded another ninety-five strains of Bcc from the National Center for Biotechnology Information database to evaluate the divergence between carbapenem-resistance and carbapenem-sensitive strains.

Results

The genetic organization between carbapenem-sensitive and carbapenem-resistant strains of Bcc showed no difference. However, in the carbapenem-sensitive strain, E151V substitution in PenR was detected. In addition, a novel specific OXA family subgroup, blaOXA-1043 in Burkholderia cenocepacia was discovered.

Conclusion

The E151V substitution in PenR may be associated with carbapenem-sensitive in Bcc. Moreover, the V151E mutation in PenR may be related to the activation of PenB, leading to Bcc resistance to carbapenems. Besides, a novel OXA family subgroup, blaOXA-1043, was found in Burkholderia cenocepacia, which differs from the previous OXA family.

Community carriage of ESBL-producing Escherichia coli and Klebsiella pneumoniae: a cross-sectional study of risk factors and comparative genomics of carriage and clinical isolates

The global prevalence of infections caused by extended-spectrum β-lactamase-producing Enterobacterales (ESBL-E) is increasing, and for Escherichia coli, observations indicate that this is partly driven by community-onset cases. The ESBL-E population structure in the community is scarcely described, and data on risk factors for carriage are conflicting. Here, we report the prevalence and population structure of fecal ESBL-producing E. coli and Klebsiella pneumoniae (ESBL-Ec/Kp) in a general adult population, examine risk factors, and compare carriage isolates with contemporary clinical isolates. Fecal samples obtained from 4,999 participants (54% women) ≥40 years in the seventh survey of the population-based Tromsø Study, Norway (2015, 2016), were screened for ESBL-Ec/Kp. In addition, we included 118 ESBL-Ec clinical isolates from the Norwegian surveillance program in 2014. All isolates were whole-genome sequenced. Risk factors associated with carriage were analyzed using multivariable logistic regression. ESBL-Ec gastrointestinal carriage prevalence was 3.3% [95% confidence interval (CI) 2.8%-3.9%, no sex difference] and 0.08% (0.02%-0.20%) for ESBL-Kp. For ESBL-Ec, travel to Asia was the only independent risk factor (adjusted odds ratio 3.46, 95% CI 2.18-5.49). E. coli ST131 was most prevalent in both collections. However, the ST131 proportion was significantly lower in carriage (24%) versus clinical isolates (58%, P < 0.001). Carriage isolates were genetically more diverse with a higher proportion of phylogroup A (26%) than clinical isolates (5%, P < 0.001), indicating that ESBL gene acquisition occurs in a variety of E. coli lineages colonizing the gut. STs commonly related to extraintestinal infections were more frequent in clinical isolates also carrying a higher prevalence of antimicrobial resistance, which could indicate clone-associated pathogenicity.IMPORTANCEESBL-Ec and ESBL-Kp are major pathogens in the global burden of antimicrobial resistance. However, there is a gap in knowledge concerning the bacterial population structure of human ESBL-Ec/Kp carriage isolates in the community. We have examined ESBL-Ec/Kp isolates from a population-based study and compared these to contemporary clinical isolates. The large genetic diversity of carriage isolates indicates frequent ESBL gene acquisition, while those causing invasive infections are more clone dependent and associated with a higher prevalence of antibiotic resistance. The knowledge of factors associated with ESBL carriage helps to identify patients at risk to combat the spread of resistant bacteria within the healthcare system. Particularly, previous travel to Asia stands out as a major risk factor for carriage and should be considered in selecting empirical antibiotic treatment in critically ill patients.

Isolation and characterization of multidrug resistant Gallibacterium anatis biovar haemolytica strains from Polish geese and hens

Gallibacterium anatis biovar haemolytica is a bacterium that is frequently associated with infections of the reproductive tract and respiratory system in poultry. To assess the current prevalence and resistance profile of these bacteria in Poland, we collected and investigated 63 strains of Gallibacterium from diseased domestic poultry flocks including geese, laying hens, breeding hens and an ornamental hen. Detailed characterization of the isolates included the analysis of phenotypic antimicrobial resistance profiles and biofilm formation ability. Furthermore, the genetic background of 40 selected isolates regarding the presence of virulence and antimicrobial resistance genes and mobile genetic elements was determined. All investigated isolates were multidrug resistant, most prominently to β-lactams, fluoroquinolones, sulfonamides and macrolides. A total of 48 different resistance profiles were detected. Of all isolates, 50.8% formed a strong biofilm, where strains isolated from geese appeared to be better at biofilm formation than strains isolated from laying and breeding hens. Single-nucleotide polymorphism genotyping revealed that G. anatis bv. haemolytica strains are restricted in host and geographical distribution, and the geese isolates showed greater phylogenetic similarity. Whole genome sequencing enabled identification of 25 different antimicrobial resistance determinants. The most common resistance genes were tetB, blaROB-1, and blaTEM-1 which may be located on mobile genetic elements. All isolates possessed the toxin gene gtxA, and the fimbrial gene flfA was identified in 95% of strains. Our results indicated that all G. anatis bv. haemolytica isolates showed multidrug resistant phenotypes. Strains isolated from geese were characterized by the highest percentage of isolates resistant to selected antimicrobials, probably reflecting host-related adaptations.