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

AMRomics: a scalable workflow to analyze large microbial genome collections

Whole genome analysis for microbial genomics is critical to studying and monitoring antimicrobial resistance strains. The exponential growth of microbial sequencing data necessitates a fast and scalable computational pipeline to generate the desired outputs in a timely and cost-effective manner. Recent methods have been implemented to integrate individual genomes into large collections of specific bacterial populations and are widely employed for systematic genomic surveillance. However, they do not scale well when the population expands and turnaround time remains the main issue for this type of analysis. Here, we introduce AMRomics, an optimized microbial genomics pipeline that can work efficiently with big datasets. We use different bacterial data collections to compare AMRomics against competitive tools and show that our pipeline can generate similar results of interest but with better performance. The software is open source and is publicly available at https://github.com/amromics/amromics under an MIT license.

Persistent enrichment of multidrug-resistant Klebsiella in oral and nasal communities during long-term starvation

BACKGROUND

The human oral and nasal cavities can act as reservoirs for opportunistic pathogens capable of causing acute infection. These microbes asymptomatically colonize the human oral and nasal cavities which facilitates transmission within human populations via the environment, and they routinely possess clinically significant antibiotic resistance genes. Among these opportunistic pathogens, the Klebsiella genus stands out as a notable example, with its members frequently linked to nosocomial infections and multidrug resistance. As with many colonizing opportunistic pathogens, the essential transmission factors influencing the spread of Klebsiella species among both healthy and diseased individuals remain unclear.

RESULTS

Here, we explored a possible explanation by investigating the ability of oral and nasal Klebsiella species to outcompete their native microbial community members under in vitro starvation conditions, which could be analogous to external hospital environments or the microenvironment of mechanical ventilators. When K. pneumoniae and K. aerogenes were present within a healthy human oral or nasal sample, the bacterial community composition shifted dramatically under starvation conditions and typically became enriched in Klebsiella species. Furthermore, introducing K. pneumoniae exogenously into a native microbial community lacking K. pneumoniae, even at low inoculum, led to repeated enrichment under starvation. Precise monitoring of K. pneumoniae within these communities undergoing starvation indicated rapid initial growth and prolonged viability compared to other members of the microbiome. K. pneumoniae strains isolated from healthy individuals' oral and nasal cavities also exhibited resistance to multiple classes of antibiotics and were genetically similar to clinical and gut isolates. In addition, we found that in the absence of Klebsiella species, other understudied opportunistic pathogens, such as Peptostreptococcus, increased in relative abundance under starvation conditions.

CONCLUSIONS

Our findings establish an environmental and microbiome community circumstance that allows for the enrichment of Klebsiella species and other opportunistic pathogens. Klebsiella's enrichment may hinge on its ability to quickly outgrow other members of the microbiome. The ability to outcompete other commensal bacteria and to persist under harsh environmental conditions could be an important factor that contributes to enhanced transmission in both commensal and pathogenic contexts. Video Abstract.

Draft genome of a human-derived pks+ E. coli that caused spontaneous disseminated infection in a mouse

We present the draft genome of a novel human-derived Escherichia coli strain isolated from a healthy control human microbiota that, when put into a mouse, spontaneously disseminated from the gut to the kidneys.

Complete genome sequence of three Staphylococcus epidermidis strains isolated from patients with skin diseases in Japan

Staphylococcus epidermidis is a member of the human skin microbiota as a commensal organism but could be an important opportunistic pathogen for immunocompromised individuals. Here, we report the complete genome sequence of three S. epidermidis strains isolated from patients with skin diseases.

A multicenter genomic epidemiological investigation in Brazil, Chile, and Mexico reveals the diversity and persistence of Salmonella populations in surface waters

This study examined the diversity and persistence of Salmonella in the surface waters of agricultural regions of Brazil, Chile, and Mexico. Research groups (three in 2019-2020 and five in 2021-2022) conducted a long-term survey of surface water across 5-8 months annually (n = 30 monthly). On-site, each team filtered 10-L water samples with modified Moore Swabs to capture Salmonella, which were then isolated and identified using conventional microbiological techniques. Salmonella isolates were sequenced on Illumina platforms. Salmonella was present in 1,493/3,291 water samples (45.8%), with varying isolation rates across countries and years. Newport, Infantis, and Typhimurium were the most frequent among the 128 different serovars. Notably, 22 serovars were found in all three countries, representing almost half of the 1,911 different isolates collected. The resistome comprised 72 antimicrobial resistance (AMR) genes and six point mutations in three genes. At least one AMR determinant was observed in 33.8% (646/1,911) of the isolates, of which 47.4% (306/646) were potentially multidrug resistant. Phylogeny based on core genome multilocus sequence typing (cgMLST) showed that most isolates clustered according to sequence type and country of origin. Only 14 cgMLST multi-country clusters were detected among the 275 clusters. However, further analysis confirmed that close genetic relatedness occurred mostly among isolates from the same country, with three exceptions. Interestingly, isolates closely related phylogenetically were recovered over multiple years within the same country, indicating the persistence of certain Salmonella in those areas. In conclusion, surface waters in these regions are consistently contaminated with diverse Salmonella, including strains that persist over time.IMPORTANCESalmonella is a leading foodborne pathogen responsible for millions of illnesses, hospitalizations, and deaths annually. Although Salmonella-contaminated water has now been recognized as an important contamination source in the agrifood chain, there is a lack of knowledge on the global occurrence and diversity of Salmonella in surface water. Moreover, there has been insufficient research on Salmonella in surface waters from Latin American countries that are major producers and exporters of agricultural products. Incorporating genetic profiling of Salmonella isolates from underrepresented regions, such as Latin America, enhances our understanding of the pathogen's ecology, evolution, antimicrobial resistance, and pathogenicity. Moreover, leveraging genomic data derived from pathogens isolated from diverse geographical areas is critical for assessing the potential public health risk posed by the pathogen and expediting investigations of foodborne outbreaks. Ultimately, global efforts contribute significantly to reducing the incidence of foodborne infections.

ResFungi: A Novel Protein Database of Antifungal Drug Resistance Genes Using a Hidden Markov Model Profile

Fungal infections vary from superficial to invasive and can be life-threatening in immunocompromised and healthy individuals. Antifungal resistance is one of the main reasons for an increasing concern about fungal infections as they become more complex and harder to treat. The fungal "omics" databases help us find drug resistance genes, which is of great importance and extremely necessary. With that in mind, we built a new platform for drug resistance genes. We added seven drug classes of resistance genes to our database: azoles (without specifying which drug), fluconazole, voriconazole, itraconazole, flucytosine, micafungin, and caspofungin. Species with known resistance genes were used to validate the results from our database. This study describes a list of 261 candidate genes related to antifungal resistance, with several genes displaying transport functions involved in azole resistance. Over 65% of the candidate genes found were related to at least one type of azole. Overall, the candidate genes found have functional annotations consistent with genes or enzymes that have been linked to antifungal resistance in previous studies. Also, candidate antifungal resistance genes found exhibit functional annotations consistent with previously described resistance mechanisms. The existence of an HMM profile focusing on antifungal resistance genes allows in silico searches for candidate genes, helping future wet lab experiments, and hence, reducing costs when studying candidate antifungal genes without prior knowledge of the species or genes. Finally, ResFungi has proven to be a powerful tool to narrow down candidate antifungal-related genes and unravel mechanisms related to resistance to help in the design of experiments focusing on the genetic basis of antifungal resistance.

Multidrug-resistant Shigella flexneri outbreak associated with a high-mortality spillover event into nonhuman primates

Shigellosis is a gastrointestinal infection caused by species of Shigella. A large outbreak of Shigella flexneri serotype 2a occurred in Albuquerque, New Mexico (NM) between May 2021 and November 2023 that involved humans and nonhuman primates (NHP) from a local zoo. We analyzed the genomes of 202 New Mexico isolates as well as 15 closely related isolates from other states, and four from NHP. The outbreak was initially detected within men who have sex with men (MSM) but then predominantly affected people experiencing homelessness (PEH). Nearly 70% of cases were hospitalized and there was one human death. The outbreak extended into Albuquerque's BioPark Zoo, causing high morbidity and six deaths in NHPs. The NHP isolates were identical to those in the human outbreak. All isolates were multidrug-resistant, including towards fluoroquinolones, a first line treatment option which led to treatment failures in human and NHP populations. We demonstrate the transmission of this S. flexneri strain between humans and NHPs, causing fatalities in both populations. This study demonstrates the threat of antimicrobial resistant organisms to vulnerable human and primate populations and emphasizes the value of vigilant genomic surveillance within a One Health framework.

Clinical and genomic features of Mycobacterium avium complex: a multi-national European study

BACKGROUND

The Mycobacterium avium complex (MAC) comprises the most frequent non-tuberculous mycobacteria (NTM) in Central Europe and currently includes twelve species. M. avium (MAV), M. intracellulare subsp. intracellulare (MINT), and M. intracellulare subsp. chimaera (MCH) are clinically most relevant. However, the population structure and genomic landscape of MAC linked with potential pathobiological differences remain little investigated.

METHODS

Whole genome sequencing (WGS) was performed on a multi-national set of MAC isolates from Germany, France, and Switzerland. Phylogenetic analysis was conducted, as well as plasmids, resistance, and virulence genes predicted from WGS data. Data was set into a global context with publicly available sequences. Finally, detailed clinical characteristics were associated with genomic data in a subset of the cohort.

RESULTS

Overall, 610 isolates from 465 patients were included. The majority could be assigned to MAV (n = 386), MCH (n = 111), and MINT (n = 77). We demonstrate clustering with less than 12 SNPs distance of isolates obtained from different patients in all major MAC species and the identification of trans-European or even trans-continental clusters when set into relation with 1307 public sequences. However, none of our MCH isolates clustered closely with the heater-cooler unit outbreak strain Zuerich-1. Known plasmids were detected in MAV (325/1076, 30.2%), MINT (62/327, 19.0%), and almost all MCH-isolates (457/463, 98.7%). Predicted resistance to aminoglycosides or macrolides was rare. Overall, there was no direct link between phylogenomic grouping and clinical manifestations, but MCH and MINT were rarely found in patients with extra-pulmonary disease (OR 0.12 95% CI 0.04-0.28, p < 0.001 and OR 0.11 95% CI 0.02-0.4, p = 0.004, respectively) and MCH was negatively associated with fulfillment of the ATS criteria when isolated from respiratory samples (OR 0.28 95% CI 0.09-0.7, p = 0.011). With 14 out of 43 patients with available serial isolates, co-infections or co-colonizations with different strains or even species of the MAC were frequent (32.6%).

CONCLUSIONS

This study demonstrates clustering and the presence of plasmids in a large proportion of MAC isolates in Europe and in a global context. Future studies need to urgently define potential ways of transmission of MAC isolates and the potential involvement of plasmids in virulence.

The usefulness of nanopore sequencing in whole-genome sequencing-based genotyping of Listeria monocytogenes and Salmonella enterica serovar Enteritidis

Bacterial genotyping through whole-genome sequencing plays a crucial role in disease surveillance and outbreak investigations in public health laboratories. This study assessed the effectiveness of Oxford Nanopore Technologies (ONT) sequencing in the genotyping of Listeria monocytogenes and Salmonella enterica serovar Enteritidis. Our results indicated that ONT sequences, generated with the R10.4.1 flow cell and basecalled using the Dorado 0.5.0 Super Accurate 4.3 model, exhibited comparable accuracy to Illumina sequences, effectively discriminating among bacterial strains from outbreaks. These findings suggest that ONT sequencing has the potential to be a promising tool for rapid whole-genome sequencing of bacterial pathogens in public health laboratories for epidemiological investigations.

IMPORTANCE

This study unveils that Oxford Nanopore Technologies sequencing, by itself, holds the potential to serve as a whole-genome sequencing-based genotyping tool in public health laboratories, enabling routine subtyping of bacterial isolates for disease surveillance and outbreak investigations.

Different sampling strategies for optimal detection of the overall genetic diversity of methicillin-resistant Staphylococcus aureus

Surveillance schemes for methicillin-resistant Staphylococcus aureus (MRSA) are widely established at the national and international levels. Due to the simple standardization of the protocol, mainly isolates from bloodstream infections are used. However, the limitations of this simple surveillance system are well described. We conducted a comprehensive analysis of MRSA isolates in a large Slovenian region over 5 years to identify the optimal sample group for assessing the overall MRSA diversity. At the same time, this study provides to date non-available molecular characterization of Slovenian MRSA isolates. A total of 306 MRSA isolates from various sources were sequenced and phenotypically tested for resistance. The isolates exhibited significant molecular diversity, encompassing 30 multi locus sequence type (MLST) sequence types (STs), 39 ST-SCCmec genetic lineages, 49 spa types, and 29 antibiotic resistance profiles. Furthermore, the isolate pool comprised 57 resistance genes, representing 22 resistance mechanisms, and 96 virulence genes. While bloodstream isolates, commonly used in surveillance, provided insights into frequently detected clones, they overlooked majority of clones and important virulence and resistance genes. Blood culture isolates detected 21.3% spa types, 24.1% resistance phenotypes, and 28.2% MLST-SCCmec profiles. In contrast, strains from soft tissues demonstrated superior genomic diversity capture, with 65.3% spa types, 58.6% resistance phenotypes, and 71.8% MLST-SCCmec profiles. These strains also encompassed 100.0% of virulence and 82.5% of resistance genes, making them better candidates for inclusion in surveillance programs. This study highlights the limitations of relying solely on bloodstream isolates in MRSA surveillance and suggests incorporating strains from soft tissues to obtain a more comprehensive understanding of the epidemiology of MRSA.IMPORTANCEIn this study, we investigated the diversity of methicillin-resistant Staphylococcus aureus (MRSA), a bacterium that can cause infections that are difficult to treat due to its resistance to antimicrobial agents. Currently, surveillance programs for MRSA mainly rely on isolates from bloodstream infections, employing a standardized protocol. However, this study highlights the limitations of this approach and introduces a more comprehensive method. The main goal was to determine which group of samples is best suited to understand the overall diversity of MRSA and to provide, for the first time, molecular characterization of Slovenian MRSA isolates. Our results suggest that including MRSA strains from soft tissue infections rather than just blood infections provides a more accurate and comprehensive view of bacterial diversity and characteristics. This insight is valuable for improving the effectiveness of surveillance programs and for developing strategies to better manage MRSA infections.

Variation in the response to antibiotics and life-history across the major Pseudomonas aeruginosa clone type (mPact) panel

Pseudomonas aeruginosa is a ubiquitous, opportunistic human pathogen. Since it often expresses multidrug resistance, new treatment options are urgently required. Such new treatments are usually assessed with one of the canonical laboratory strains, PAO1 or PA14. However, these two strains are unlikely representative of the strains infecting patients, because they have adapted to laboratory conditions and do not capture the enormous genomic diversity of the species. Here, we characterized the major P. aeruginosa clone type (mPact) panel. This panel consists of 20 strains, which reflect the species' genomic diversity, cover all major clone types, and have both patient and environmental origins. We found significant strain variation in distinct responses toward antibiotics and general growth characteristics. Only few of the measured traits are related, suggesting independent trait optimization across strains. High resistance levels were only identified for clinical mPact isolates and could be linked to known antimicrobial resistance (AMR) genes. One strain, H01, produced highly unstable AMR combined with reduced growth under drug-free conditions, indicating an evolutionary cost to resistance. The expression of microcolonies was common among strains, especially for strain H15, which also showed reduced growth, possibly indicating another type of evolutionary trade-off. By linking isolation source, growth, and virulence to life history traits, we further identified specific adaptive strategies for individual mPact strains toward either host processes or degradation pathways. Overall, the mPact panel provides a reasonably sized set of distinct strains, enabling in-depth analysis of new treatment designs or evolutionary dynamics in consideration of the species' genomic diversity.

IMPORTANCE

New treatment strategies are urgently needed for high-risk pathogens such as the opportunistic and often multidrug-resistant pathogen Pseudomonas aeruginosa. Here, we characterize the major P. aeruginosa clone type (mPact) panel. It consists of 20 strains with different origins that cover the major clone types of the species as well as its genomic diversity. This mPact panel shows significant variation in (i) resistance against distinct antibiotics, including several last resort antibiotics; (ii) related traits associated with the response to antibiotics; and (iii) general growth characteristics. We further developed a novel approach that integrates information on resistance, growth, virulence, and life-history characteristics, allowing us to demonstrate the presence of distinct adaptive strategies of the strains that focus either on host interaction or resource processing. In conclusion, the mPact panel provides a manageable number of representative strains for this important pathogen for further in-depth analyses of treatment options and evolutionary dynamics.

A challenging STEC strain isolation from patients' stools: an O166:H15 STEC strain with the stx2 gene

Two patients with acute gastroenteritis tested positive for Shiga toxin-producing Escherichia coli (STEC) by polymerase chain reaction (PCR), and both strains carried the Shiga toxin 2 encoding gene. Since routine culture using CHROMagar STEC failed to recover these isolates, immunomagnetic separation (IMS) targeting the top six non-O157:H7 serotypes was used for isolate recovery. After two subsequent IMS runs, the STEC strains were isolated from trypticase soy broth with and without overnight enrichment for runs 1 and 2, respectively. Serotyping based on whole-genome sequencing revealed that both patients carried the strain O166:H15 STEC with the stx2 gene. Hence, the magnetic beads used in IMS appeared to have cross-reactivity with other E. coli serotypes. When the STEC isolates from both stools were cultured on CHROMagar STEC and sheep blood agar (BAP), two distinct colony sizes were apparent after overnight incubation. The small and large colonies were picked and separately cultured on both media, and colony growth was observed for 2 weeks at room temperature after an initial overnight incubation at 37°C. After 1 week, the colonies showed concentric ring structures with a darker center and a lighter surrounding on CHROMagar STEC and a "fried egg"-resembling structure with a raised circular center and a flat surrounding on BAP. Both colony types remained morphologically different on CHROMagar STEC throughout the 15 days. However, on BAP, their appearance was comparable by day 7.

IMPORTANCE

Shiga toxin-producing E. coli (STEC) infections can lead to severe complications such as bloody diarrhea and hemolytic uremic syndrome (HUS), especially in young children and the elderly. Strains that carry the shiga toxin 2 gene (stx2), such as O157:H7, have been mostly linked with severe disease outcomes. In recent years, outbreaks caused by non-O157:H7 strains have increased. E. coli O166:H15 has been previously reported causing a gastroenteritis outbreak in 1996 as a non-STEC strain, however the O166:H15 serotype we recovered carried the stx2 gene. It was particularly challenging to isolate this strain from stools by culture. Consequently, we tested immunomagnetic separation for the STEC recovery, which was a novel approach on clinical stools. Virulence genes were included for the characterization of these isolates.

Overcoming Microbial Inhibition of S. Sonnei Through the Exploitation of Genomically Predicted Antibiotic Resistance Profiles for the Development of Food Enrichment Media

Linking outbreaks of Shigella spp. to specific foods is challenging due to poor selectivity of current enrichment media. We have previously shown that enrichment media, tailored to the genomically-predicted antimicrobial resistance (AMR) of Shiga toxigenic E. coli strains, enhances their isolation from foods. This study investigates the application of this approach for Shigella isolation. The AMR gene profiles of 21,908 published S. sonnei genomes indicated a high prevalence of genes conferring resistance to streptomycin (aadA, aph(3″)-Ib, aph(6)-Id, 92.8%), sulfonamides (sul1, sul2, 74.8%), and/or trimethoprim (dfrA, 96.2%). Genomic analysis and antibiotic susceptibility testing conducted with a panel of 17 outbreak-associated S. sonnei strains confirmed the correlation of AMR gene detection with resistance phenotypes. Supplementation of Shigella Broth (SB) with up to 400 µg/mL of trimethoprim or sulfadiazine did not suppress the growth of sensitive strains, whereas 100 µg/mL of streptomycin increased the selectivity of this broth. All three antibiotics increased the selectivity of modified Tryptone Soya Broth (mTSB). Based on these results, supplemented media formulations were developed and assessed by measuring the relative growth of S. sonnei in cultures coinoculated with a strain of bacteriocin-producing E. coli that is inhibitory to Shigella growth. S. sonnei was not recovered from cocultures grown in SB or mTSB without antibiotics. In contrast, media supplemented with streptomycin at 50 and 100 µg/mL, trimethoprim at 25 and 50 µg/mL, and sulfadiazine at 100 µg/mL increased the relative proportion of S. sonnei in postenrichment cultures. The enhanced recovery of resistant S. sonnei strains achieved in this study indicates that, in cases where genomic data are available for clinical S. sonnei isolates, customization of selective enrichment media based on AMR gene detection could be a valuable tool for supporting the investigation of foodborne shigellosis outbreaks.

Five centuries of genome evolution and multi-host adaptation of Campylobacter jejuni in Brazil

Consumption of raw, undercooked or contaminated animal food products is a frequent cause of Campylobacter jejuni infection. Brazil is the world's third largest producer and a major exporter of chicken meat, yet population-level genomic investigations of C. jejuni in the country remain scarce. Analysis of 221 C. jejuni genomes from Brazil shows that the overall core and accessory genomic features of C. jejuni are influenced by the identity of the human or animal source. Of the 60 sequence types detected, ST353 is the most prevalent and consists of samples from chicken and human sources. Notably, we identified the presence of diverse bla genes from the OXA-61 and OXA-184 families that confer beta-lactam resistance as well as the operon cmeABCR related to multidrug efflux pump, which contributes to resistance against tetracyclines, macrolides and quinolones. Based on limited data, we estimated the most recent common ancestor of ST353 to the late 1500s, coinciding with the time the Portuguese first arrived in Brazil and introduced domesticated chickens into the country. We identified at least two instances of ancestral chicken-to-human infections in ST353. The evolution of C. jejuni in Brazil was driven by the confluence of clinically relevant genetic elements, multi-host adaptation and clonal population growth that coincided with major socio-economic changes in poultry farming.

Genomic study of European Clostridioides difficile ribotype 002/sequence type 8

Clostridioides difficile has significant clinical importance as a leading cause of healthcare-associated infections, with symptoms ranging from mild diarrhoea to severe colitis, and possible life-threatening complications. C. difficile ribotype (RT) 002, mainly associated with MLST sequence type (ST) 8, is one of the most common RTs found in humans. This study aimed at investigating the genetic characteristics of 537 C. difficile genomes of ST8/RT002. To this end, we sequenced 298 C. difficile strains representing a new European genome collection, with strains from Germany, Denmark, France and Portugal. These sequences were analysed against a global dataset consisting of 1,437 ST8 genomes available through Enterobase. Our results showed close genetic relatedness among the studied ST8 genomes, a diverse array of antimicrobial resistance (AMR) genes and the presence of multiple mobile elements. Notably, the pangenome analysis revealed an open genomic structure. ST8 shows relatively low overall variation. Thus, clonal isolates were found across different One Health sectors (humans, animals, environment and food), time periods, and geographical locations, suggesting the lineage's stability and a universal environmental source. Importantly, this stability did not hinder the acquisition of AMR genes, emphasizing the adaptability of this bacterium to different selective pressures. Although only 2.4 % (41/1,735) of the studied genomes originated from non-human sources, such as animals, food, or the environment, we identified 9 cross-sectoral core genome multilocus sequence typing (cgMLST) clusters. Our study highlights the importance of ST8 as a prominent lineage of C. difficile with critical implications in the context of One Health. In addition, these findings strongly support the need for continued surveillance and investigation of non-human samples to gain a more comprehensive understanding of the epidemiology of C. difficile.

Comparative genome analyses of clinical and non-clinical Clostridioides difficile strains

The pathogenic bacterium Clostridioides difficile is a worldwide health burden with increasing morbidity, mortality and antibiotic resistances. Therefore, extensive research efforts are made to unravel its virulence and dissemination. One crucial aspect for C. difficile is its mobilome, which for instance allows the spread of antibiotic resistance genes (ARG) or influence strain virulence. As a nosocomial pathogen, the majority of strains analyzed originated from clinical environments and infected individuals. Nevertheless, C. difficile can also be present in human intestines without disease development or occur in diverse environmental habitats such as puddle water and soil, from which several strains could already be isolated. We therefore performed comprehensive genome comparisons of closely related clinical and non-clinical strains to identify the effects of the clinical background. Analyses included the prediction of virulence factors, ARGs, mobile genetic elements (MGEs), and detailed examinations of the pan genome. Clinical-related trends were thereby observed. While no significant differences were identified in fundamental C. difficile virulence factors, the clinical strains carried more ARGs and MGEs, and possessed a larger accessory genome. Detailed inspection of accessory genes revealed higher abundance of genes with unknown function, transcription-associated, or recombination-related activity. Accessory genes of these functions were already highlighted in other studies in association with higher strain virulence. This specific trend might allow the strains to react more efficiently on changing environmental conditions in the human host such as emerging stress factors, and potentially increase strain survival, colonization, and strain virulence. These findings indicated an adaptation of the strains to the clinical environment. Further, implementation of the analysis results in pairwise genome comparisons revealed that the majority of these accessory genes were encoded on predicted MGEs, shedding further light on the mobile genome of C. difficile. We therefore encourage the inclusion of non-clinical strains in comparative analyses.

RATA: A novel class A carbapenemase with broad geographic distribution and potential for global spread

Carbapenem resistance's global proliferation poses a significant public health challenge. The primary resistance mechanism is carbapenemase production. In this study, we discovered a novel carbapenemase, RATA, located on the chromosome of Riemerella anatipestifer isolates. This enzyme shares ≤52 % amino acid sequence identity with other known β-lactamases. Antimicrobial susceptibility tests and kinetic assays demonstrated that RATA could hydrolyze not only penicillins and extended-spectrum cephalosporins but also monobactams, cephamycins, and carbapenems. Furthermore, its activity was readily inhibited by β-lactamase inhibitors. Bioinformatic analysis revealed 46 blaRATA-like genes encoding 27 variants in the NCBI database, involving 21 different species, including pathogens, host-associated bacteria, and environmental isolates. Notably, blaRATA-positive strains were globally distributed and primarily collected from marine environments. Concurrently, taxonomic analysis and GC content analysis indicated that blaRATA orthologue genes were predominantly located on the chromosomes of Flavobacteriaceae and shared a similar GC content as Flavobacteriaceae. Although no explicit mobile genetic elements were identified by genetic environment analysis, blaRATA-2 possessed the ability of horizontal transfer in R. anatipestifer via natural transformation. This work's data suggest that RATA is a new chromosome-encoded class A carbapenemase, and Flavobacteriaceae from marine environments could be the primary reservoir of the blaRATA gene.

Unravelling the mechanisms of antibiotic and heavy metal resistance co-selection in environmental bacteria

The co-selective pressure of heavy metals is a contributor to the dissemination and persistence of antibiotic resistance genes in environmental reservoirs. The overlapping range of antibiotic and metal contamination and similarities in their resistance mechanisms point to an intertwined evolutionary history. Metal resistance genes are known to be genetically linked to antibiotic resistance genes, with plasmids, transposons, and integrons involved in the assembly and horizontal transfer of the resistance elements. Models of co-selection between metals and antibiotics have been proposed, however, the molecular aspects of these phenomena are in many cases not defined or quantified and the importance of specific metals, environments, bacterial taxa, mobile genetic elements, and other abiotic or biotic conditions are not clear. Co-resistance is often suggested as a dominant mechanism, but interpretations are beset with correlational bias. Proof of principle examples of cross-resistance and co-regulation has been described but more in-depth characterizations are needed, using methodologies that confirm the functional expression of resistance genes and that connect genes with specific bacterial hosts. Here, we comprehensively evaluate the recent evidence for different models of co-selection from pure culture and metagenomic studies in environmental contexts and we highlight outstanding questions.

Metagenomic Investigation of the Short-Term Temporal and Spatial Dynamics of the Bacterial Microbiome and the Resistome Downstream of a Wastewater Treatment Plant in the Iskar River in Bulgaria

Waste Water Treatment Plants (WWTP) aim to reduce contamination in effluent water; however, studies indicate antimicrobial resistance genes (ARGs) persist post-treatment, potentially leading to their spread from human populated areas into the environment. This study evaluated the impact of a large WWTP serving 125,000 people on the Iskar River in Bulgaria, by characterizing the spatial and short-term temporal dynamics in bacterial community dynamics and resistance profiles of the surface water. Pairs of samples were collected biweekly on four dates from two different locations, one about 800 m after the WWTP effluents and the other 10 km downstream. Taxonomic classification revealed the dominance of Pseudomonodota and Bacteriodota, notably the genera Flavobacterium, Aquirufa, Acidovorax, Polynucleobacter, and Limnohabitans. The taxonomic structure corresponded with both lentic and lotic freshwater habitats, with Flavobacterium exhibiting a significant decrease over the study period. Principal Coordinate Analysis revealed statistically significant differences in bacterial community composition between samples collected on different dates. Differential abundance analysis identified notable enrichment of Polynucleobacter and Limnohabitans. There were shifts within the enriched or depleted bacterial taxa between early and late sampling dates. High relative abundance of the genes erm(B), erm(F), mph(E), msr(E) (macrolides); tet(C), tet(O), tet(W), tet(Q) and tet(X) (tetracyclines); sul1 and sul2 (sulphonamides); and cfxA3, cfxA6 (beta-lactams) were detected, with trends of increased presence in the latest sampling dates and in the location closer to the WWTP. Of note, genes conferring resistance to carbapenems blaOXA-58 and blaIMP-33-like were identified. Co-occurrence analysis of ARGs and mobile genetic elements on putative plasmids showed few instances, and the estimated human health risk score (0.19) according to MetaCompare2.0 was low. In total, 29 metagenome-assembled genomes were recovered, with only a few harbouring ARGs. This study enhances our understanding of freshwater microbial community dynamics and antibiotic resistance profiles, highlighting the need for continued ARGs monitoring.

Pathogenic Bacteroides fragilis strains can emerge from gut-resident commensals

Bacteroides fragilis is a prominent member of the human gut microbiota, playing crucial roles in maintaining gut homeostasis and host health. Although it primarily functions as a beneficial commensal, B. fragilis can become pathogenic. To determine the genetic basis of its duality, we conducted a comparative genomic analysis of 813 B. fragilis strains, representing both commensal and pathogenic origins. Our findings reveal that pathogenic strains emerge across diverse phylogenetic lineages, due in part to rapid gene exchange and the adaptability of the accessory genome. We identified 16 phylogenetic groups, differentiated by genes associated with capsule composition, interspecies competition, and host interactions. A microbial genome-wide association study identified 44 genes linked to extra-intestinal survival and pathogenicity. These findings reveal how genomic diversity within commensal species can lead to the emergence of pathogenic traits, broadening our understanding of microbial evolution in the gut.

Deep sequencing of Escherichia coli exposes colonisation diversity and impact of antibiotics in Punjab, Pakistan

Multi-drug resistant (MDR) E. coli constitute a major public health burden globally, reaching the highest prevalence in the global south yet frequently flowing with travellers to other regions. However, our comprehension of the entire genetic diversity of E. coli colonising local populations remains limited. We quantified this diversity, its associated antimicrobial resistance (AMR), and assessed the impact of antibiotic use by recruiting 494 outpatients and 423 community dwellers in the Punjab province, Pakistan. Rectal swab and stool samples were cultured on CLED agar and DNA extracted from plate sweeps was sequenced en masse to capture both the genetic and AMR diversity of E. coli. We assembled 5,247 E. coli genomes from 1,411 samples, displaying marked genetic diversity in gut colonisation. Compared with high income countries, the Punjabi population generally showed a markedly different distribution of genetic lineages and AMR determinants, while use of antibiotics elevated the prevalence of well-known globally circulating MDR clinical strains. These findings implicate that longitudinal multi-regional genomics-based surveillance of both colonisation and infections is a prerequisite for developing mechanistic understanding of the interplay between ecology and evolution in the maintenance and dissemination of (MDR) E. coli.

The Description and Analysis of the Complete Genome of Dermacoccus barathri FBCC-B549 Strain

Dermacoccus barathri is the first reported pathogen within the Dermacoccus genus to cause a catheter-related bloodstream infection, which occurred in 2015. In this study, the complete genome assembly of Dermacoccus barathri was constructed, and the complete genome of Dermacoccus barathri FBCC-B549 consists of a single chromosome (3,137,745 bp) without plasmids. The constructed genome of D. barathri was compared with those of two closely related species within the Dermacoccus genus. D. barathri exhibited a pattern similar to Dermacoccus abyssi in terms of gene clusters and synteny analysis. Contrary to previous studies, biosynthetic gene cluster (BGC) analysis for predicting secondary metabolites revealed the presence of the LAP biosynthesis pathway in the complete genome of D. barathri, predicting the potential synthesis of the secondary metabolite plantazolicin. Furthermore, an analysis to investigate the potential pathogenicity of D. barathri did not reveal any antibiotic resistance genes; however, nine virulence factors were identified in the Virulence Factor Database (VFDB). According to these matching results in the VFDB, despite identifying a few factors involved in biofilm formation, further research is required to determine the actual impact of D. barathri on pathogenicity. The complete genome of D. barathri is expected to serve as a valuable resource for future studies on D. barathri, which currently lack sufficient genomic sequence information.

Antimicrobial resistance and phylogenetic relatedness of extended-spectrum ß-lactamase (ESBL)-producing Escherichia coli in peridomestic rats (Rattus norvegicus and Rattus tanezumi) linked to city areas and animal farms in Hong Kong

Extended-spectrum β-lactamase-producing Escherichia (E.) coli (ESBL-EC) in the clinical setting have emerged as a major threat to public and animal health. Wildlife, including Rattus spp. may serve as reservoirs and spreaders of ESBL-EC in the environment. Peridomestic rats are well adapted to living in proximity to humans and animals in a variety of urban and agricultural environments and may serve as sentinels to identify variations of ESBL-EC within their different habitats. In this study, a set of 221 rats (Rattus norvegicus, R. tanezumi, R. andamanensis, and Niviventer huang) consisting of 104 rats from city areas, 44 from chicken farms, 52 from pig farms, and 21 from stables of horse-riding schools were screened for ESBL-EC. Overall, a total of 134 ESBL-EC were isolated from the caecal samples of 130 (59%) rats. The predominant blaESBL genes were blaCTX-M-14, blaCTX-M-15, blaCTX-M-55, and blaCTX-M-65. Phylogenetic analysis revealed a total of 62 sequence types (STs) and 17 SNP clusters. E. coli ST10 and ST155 were common to ESBL-EC from city areas and chicken farms, and ST44 were found among ESBL-EC from city areas and pig farms. Extra-intestinal pathogenic E. coli (ExPEC) ST69, ST131 and ST1193 were found exclusively among rats from city areas, and avian pathogenic E. coli (APEC) ST177 was restricted to ESBL-EC originating from chicken farms. Phylogenetic analysis showed that the populations of rodent ESBL-EC from city areas, chicken farms and pig farms were genetically different, suggesting a certain degree of partitioning between the human and animal locations. This study contributes to current understanding of ESBL-EC occurring in rats in ecologically diverse locations.

Type IV-A3 CRISPR-Cas systems drive inter-plasmid conflicts by acquiring spacers in trans

Plasmid-encoded type IV-A CRISPR-Cas systems lack an acquisition module, feature a DinG helicase instead of a nuclease, and form ribonucleoprotein complexes of unknown biological functions. Type IV-A3 systems are carried by conjugative plasmids that often harbor antibiotic-resistance genes and their CRISPR array contents suggest a role in mediating inter-plasmid conflicts, but this function remains unexplored. Here, we demonstrate that a plasmid-encoded type IV-A3 system co-opts the type I-E adaptation machinery from its host, Klebsiella pneumoniae (K. pneumoniae), to update its CRISPR array. Furthermore, we reveal that robust interference of conjugative plasmids and phages is elicited through CRISPR RNA-dependent transcriptional repression. By silencing plasmid core functions, type IV-A3 impacts the horizontal transfer and stability of targeted plasmids, supporting its role in plasmid competition. Our findings shed light on the mechanisms and ecological function of type IV-A3 systems and demonstrate their practical efficacy for countering antibiotic resistance in clinically relevant strains.

Complete genome sequence of Pseudomonas canadensis strain Pcan-CK-23 isolated from Zophobas morio larvae

We present the complete genome sequence of Pseudomonas canadensis. The strain (Pcan-CK-23) was isolated from Zophobas morio (superworm) larvae. The genome consisted of a 6,424,469 bp chromosome with a GC content of 60.3% and 5,973 genes. Pcan-CK-23 can be used as a reference genome for further studies with P. canadensis.

Complete genome sequence of the linezolid-resistant clinical Enterococcus faecalis N23-3408 linked to a livestock lineage in Switzerland

Here, we report the complete genome of human clinical linezolid-resistant Enterococcus faecalis N23-3408. N23-3408 harbored a 59.5 kb plasmid with antimicrobial resistance genes cat, erm(B), fexA, optrA, tet(L), and tet(M). Closely related E. faecalis harboring this plasmid was previously obtained from livestock animals and pet food in Switzerland.

Exploration of mobile genetic elements in the ruminal microbiome of Nellore cattle

Metagenomics has made it feasible to elucidate the intricacies of the ruminal microbiome and its role in the differentiation of animal production phenotypes of significance. The search for mobile genetic elements (MGEs) has taken on great importance, as they play a critical role in the transfer of genetic material between organisms. Furthermore, these elements serve a dual purpose by controlling populations through lytic bacteriophages, thereby maintaining ecological equilibrium and driving the evolutionary progress of host microorganisms. In this study, we aimed to identify the association between ruminal bacteria and their MGEs in Nellore cattle using physical chromosomal links through the Hi-C method. Shotgun metagenomic sequencing and the proximity ligation method ProxiMeta were used to analyze DNA, getting 1,713,111,307 bp, which gave rise to 107 metagenome-assembled genomes from rumen samples of four Nellore cows maintained on pasture. Taxonomic analysis revealed that most of the bacterial genomes belonged to the families Lachnospiraceae, Bacteroidaceae, Ruminococcaceae, Saccharofermentanaceae, and Treponemataceae and mostly encoded pathways for central carbon and other carbohydrate metabolisms. A total of 31 associations between host bacteria and MGE were identified, including 17 links to viruses and 14 links to plasmids. Additionally, we found 12 antibiotic resistance genes. To our knowledge, this is the first study in Brazilian cattle that connect MGEs with their microbial hosts. It identifies MGEs present in the rumen of pasture-raised Nellore cattle, offering insights that could advance biotechnology for food digestion and improve ruminant performance in production systems.

Genomic Characterization of Multidrug-Resistant Enterobacteriaceae Clinical Isolates from Southern Thailand Hospitals: Unraveling Antimicrobial Resistance and Virulence Mechanisms

The emergence and spread of antimicrobial resistance (AMR) among Enterobacteriaceae pose significant threats to global public health. In this study, we conducted a short-term surveillance effort in Southern Thailand hospitals to characterize the genomic diversity, AMR profiles, and virulence factors of Enterobacteriaceae strains. We identified 241 carbapenem-resistant Enterobacteriaceae, of which 12 were selected for whole-genome sequencing (WGS) and genome analysis. The strains included Proteus mirabilis, Serratia nevei, Klebsiella variicola, Klebsiella aerogenes, Klebsiella indica, Klebsiella grimontii, Phytobacter ursingii, Phytobacter palmae, Kosakonia spp., and Citrobacter freundii. The strains exhibited high levels of multidrug resistance, including resistance to carbapenem antibiotics. Whole-genome sequencing revealed a diverse array of antimicrobial resistance genes (ARGs), with strains carrying genes for ß-lactamase, efflux pumps, and resistance to other antibiotic classes. Additionally, stress response, metal tolerance, and virulence-associated genes were identified, highlighting the adaptability and pathogenic potential of these strains. A plasmid analysis identified several plasmid replicons, including IncA/C2, IncFIB(K), and Col440I, as well as several plasmids identical to those found globally, indicating the potential for the horizontal gene transfer of ARGs. Importantly, this study also identified a novel species of Kosakonia spp. PSU27, adding to the understanding of the genetic diversity and resistance mechanisms of Enterobacteriaceae in Southern Thailand. The results reported in this study highlight the critical importance of implementing effective antimicrobial management programs and developing innovative treatment approaches to urgently tackle AMR.

Assessment of the pathogen genomics landscape highlights disparities and challenges for effective AMR Surveillance and outbreak response in the East African community

The East African Community (EAC) grapples with many challenges in tackling infectious disease threats and antimicrobial resistance (AMR), underscoring the importance of regional and robust pathogen genomics capacities. However, a significant disparity exists among EAC Partner States in harnessing bacterial pathogen sequencing and data analysis capabilities for effective AMR surveillance and outbreak response. This study assesses the current landscape and challenges associated with pathogen next-generation sequencing (NGS) within EAC, explicitly focusing on World Health Organization (WHO) AMR-priority pathogens. The assessment adopts a comprehensive approach, integrating a questionnaire-based survey amongst National Public Health Laboratories (NPHLs) with an analysis of publicly available metadata on bacterial pathogens isolated in the EAC countries. In addition to the heavy reliance on third-party organizations for bacterial NGS, the findings reveal a significant disparity among EAC member States in leveraging bacterial pathogen sequencing and data analysis. Approximately 97% (n = 4,462) of publicly available high-quality bacterial genome assemblies of samples collected in the EAC were processed and analyzed by external organizations, mainly in Europe and North America. Tanzania led in-country sequencing efforts, followed by Kenya and Uganda. The other EAC countries had no publicly available samples or had all their samples sequenced and analyzed outside the region. Insufficient local NGS sequencing facilities, limited bioinformatics expertise, lack of adequate computing resources, and inadequate data-sharing mechanisms are among the most pressing challenges that hinder the EAC's NPHLs from effectively leveraging pathogen genomics data. These insights emphasized the need to strengthen microbial pathogen sequencing and data analysis capabilities within the EAC to empower these laboratories to conduct pathogen sequencing and data analysis independently. Substantial investments in equipment, technology, and capacity-building initiatives are crucial for supporting regional preparedness against infectious disease outbreaks and mitigating the impact of AMR burden. In addition, collaborative efforts should be developed to narrow the gap, remedy regional imbalances, and harmonize NGS data standards. Supporting regional collaboration, strengthening in-country genomics capabilities, and investing in long-term training programs will ultimately improve pathogen data generation and foster a robust NGS-driven AMR surveillance and outbreak response in the EAC, thereby supporting global health initiatives.

Salmonella Typhi Haplotype 58 (H58) Biofilm Formation and Genetic Variation in Typhoid Fever Patients with Gallstones in an Endemic Setting in Kenya

The causative agent of typhoid fever, Salmonella enterica serovar Typhi, is a human restricted pathogen. Human carriers, 90% of whom have gallstones in their gallbladder, continue to shed the pathogen after treatment. The genetic mechanisms involved in establishing the carrier state are poorly understood, but S. Typhi is thought to undergo specific genetic changes within the gallbladder as an adaptive mechanism. In the current study, we aimed to identify biofilm forming ability and the genetic differences in longitudinal clinical S. Typhi isolates from asymptomatic carriers with gallstones in Nairobi, Kenya. Whole genome sequences were analyzed from 22 S. Typhi isolates, 20 from stool and 2 from blood samples, all genotype 4.3.1 (H58). Nineteen strains were from four patients also diagnosed with gallstones, of whom, three had typhoid symptoms and continued to shed S. Typhi after treatment. All isolates had point mutations in the quinolone resistance determining region (QRDR) and only sub-lineage 4.3.1.2EA3 encoded multidrug resistance genes. There was no variation in antimicrobial resistance patterns among strains from the same patient/household. Non-multidrug resistant (MDR), isolates formed significantly stronger biofilms in vitro than the MDR isolates, p<0.001. A point mutation within the treB gene (treB A383T) was observed in strains isolated after clinical resolution from patients living in 75% of the households. Missense mutations in Vi capsular polysaccharide genes, tviE P263S was also observed in 18% of the isolates. This study provides insights into the role of typhoid carriage, biofilm formation, AMR genes and genetic variations in S. Typhi from asymptomatic carriers.

Advancing metagenome-assembled genome-based pathogen identification: unraveling the power of long-read assembly algorithms in Oxford Nanopore sequencing

Oxford Nanopore sequencing is one of the high-throughput sequencing technologies that facilitates the reconstruction of metagenome-assembled genomes (MAGs). This study aimed to assess the potential of long-read assembly algorithms in Oxford Nanopore sequencing to enhance the MAG-based identification of bacterial pathogens using both simulated and mock communities. Simulated communities were generated to mimic those on fresh spinach and in surface water. Long reads were produced using R9.4.1+SQK-LSK109 and R10.4 + SQK-LSK112, with 0.5, 1, and 2 million reads. The simulated bacterial communities included multidrug-resistant Salmonella enterica serotypes Heidelberg, Montevideo, and Typhimurium in the fresh spinach community individually or in combination, as well as multidrug-resistant Pseudomonas aeruginosa in the surface water community. Real data sets of the ZymoBIOMICS HMW DNA Standard were also studied. A bioinformatic pipeline (MAGenie, freely available at https://github.com/jackchen129/MAGenie) that combines metagenome assembly, taxonomic classification, and sequence extraction was developed to reconstruct draft MAGs from metagenome assemblies. Five assemblers were evaluated based on a series of genomic analyses. Overall, Flye outperformed the other assemblers, followed by Shasta, Raven, and Unicycler, while Canu performed least effectively. In some instances, the extracted sequences resulted in draft MAGs and provided the locations and structures of antimicrobial resistance genes and mobile genetic elements. Our study showcases the viability of utilizing the extracted sequences for precise phylogenetic inference, as demonstrated by the consistent alignment of phylogenetic topology between the reference genome and the extracted sequences. R9.4.1+SQK-LSK109 was more effective in most cases than R10.4+SQK-LSK112, and greater sequencing depths generally led to more accurate results.IMPORTANCEBy examining diverse bacterial communities, particularly those housing multiple Salmonella enterica serotypes, this study holds significance in uncovering the potential of long-read assembly algorithms to improve metagenome-assembled genome (MAG)-based pathogen identification through Oxford Nanopore sequencing. Our research demonstrates that long-read assembly stands out as a promising avenue for boosting precision in MAG-based pathogen identification, thus advancing the development of more robust surveillance measures. The findings also support ongoing endeavors to fine-tune a bioinformatic pipeline for accurate pathogen identification within complex metagenomic samples.

Genetic determinants of resistance to antimicrobial therapeutics are rare in publicly available Clostridioides difficile genome sequences

OBJECTIVES

To determine the frequencies and clonal distributions of putative genetic determinants of resistance to antimicrobials applied for treatment of Clostridioides difficile infection (CDI), as documented in the genomic record.

METHODS

We scanned 26 557 C. difficile genome sequences publicly available from the EnteroBase platform for plasmids, point mutations and gene truncations previously reported to reduce susceptibility to vancomycin, fidaxomicin or metronidazole, respectively. We measured the antimicrobial susceptibility of 143 selected C. difficile isolates.

RESULTS

The frequency of mutations causing reduced susceptibility to vancomycin and metronidazole, respectively, increased strongly after 2000, peaking at up to 52% of all sequenced C. difficile genomes. However, both mutations declined sharply more recently, reflecting major changes in CDI epidemiology. We detected mutations associated with fidaxomicin resistance in several major genotypes, but found no evidence of international spread of resistant clones. The pCD-METRO plasmid, conferring metronidazole resistance, was detected in a single previously unreported C. difficile isolate, recovered from a hospital patient in Germany in 2008. The pX18-498 plasmid, putatively associated with decreased vancomycin susceptibility, was confined to related, recent isolates from the USA. Phenotype measurements confirmed that most of those genetic features were useful predictors of antibiotic susceptibility, even though ranges of MICs typically overlapped among isolates with and without specific mutations.

CONCLUSIONS

Genomic data suggested that resistance to therapeutic antimicrobial drugs is rare in C. difficile. Public antimicrobial resistance marker databases were not equipped to detect most of the genetic determinants relevant to antibiotic therapy of CDI.

Genomic insights of Fictibacillus terranigra sp. nov., a versatile metabolic bacterium from Amazonian Dark Earths

The Amazon rainforest, a hotspot for biodiversity, is a crucial research area for scientists seeking novel microorganisms with ecological and biotechnological significance. A key region within the Amazon rainforest is the Amazonian Dark Earths (ADE), noted for supporting diverse plant and microbial communities, and its potential as a blueprint for sustainable agriculture. This study delineates the isolation, morphological traits, carbon source utilization, and genomic features of Fictibacillus terranigra CENA-BCM004, a candidate novel species of the Fictibacillus genus isolated from ADE. The genome of Fictibacillus terranigra was sequenced, resulting in 16 assembled contigs, a total length of 4,967,627 bp, and a GC content of 43.65%. Genome annotation uncovered 3315 predicted genes, encompassing a wide range of genes linked to various metabolic pathways. Phylogenetic analysis indicated that CENA-BCM004 is a putative new species, closely affiliated with other unidentified Fictibacillus species and Bacillus sp. WQ 8-8. Moreover, this strain showcased a multifaceted metabolic profile, revealing its potential for diverse biotechnological applications. It exhibited capabilities to antagonize pathogens, metabolize multiple sugars, mineralize organic matter compounds, and solubilize several minerals. These insights substantially augment our comprehension of microbial diversity in ADE and underscore the potential of Fictibacillus terranigra as a precious resource for biotechnological endeavors. The genomic data generated from this study will serve as a foundational resource for subsequent research and exploration of the biotechnological capabilities of this newly identified species.

Clinical, phenotypic, and genotypic characteristics of ESBL-producing Salmonella enterica bloodstream infections from Qatar

Background

Resistant Salmonella infections are a major global public health challenge particularly for multidrug-resistant (MDR) isolates manifesting as bloodstream infections (BSIs).

Objectives

To evaluate clinical, phenotypic, and genotypic characteristics of extended-spectrum beta-lactamase (ESBL) producing Salmonella enterica BSIs from Qatar.

Methods

Phenotypic ESBL Salmonella enterica from adult patients presenting with positive BSIs were collected between January 2019 to May 2020. Microbiological identification and characterization were performed using standard methods while genetic characteristics were examined through whole genome sequencing studies.

Results

Of 151 episodes of Salmonella enterica BSI, 15 (10%) phenotypic ESBL isolates were collected. Recent travel was recorded in most cases (80%) with recent exposure to antimicrobials (27%). High-level resistance to quinolines, aminoglycosides, and cephalosporins was recorded (80-100%) while meropenem, tigecycline and colistin demonstrated universal susceptibility. Genomic evaluation demonstrated dominance of serotype Salmonella Typhi sequence type 1 (93%) while antimicrobial resistance genes revealed dominance of aminoglycoside resistance (100%), qnrS1 quinolones resistance (80%), blaCTX-M-15 ESBLs (86.7%), and paucity of AmpC resistance genes (6.7%).

Conclusions

Invasive MDR Salmonella enterica is mainly imported, connected to patients from high prevalent regions with recent travel and antimicrobial use caused by specific resistant clones. In suspected cases of multidrug resistance, carbapenem therapy is recommended.

Approaching toxigenic Clostridia from a One Health perspective

Spore-forming pathogens have a unique capacity to thrive in diverse environments, and with temporal persistence afforded through their ability to sporulate. Their prevalence in diverse ecosystems requires a One Health approach to identify critical reservoirs and outbreak-associated transmission chains, given their capacity to freely move across soils, waterways, foodstuffs and as commensals or infecting pathogens in human and animal populations. Among anaerobic spore-formers, genomic resources for pathogens including C. botulinum, C. difficile, and C. perfringens enable our capacity to identify common and unique factors that support their persistence in diverse reservoirs and capacity to cause disease. Publicly available genomic resources for spore-forming pathogens at NCBI's Pathogen Detection program aid outbreak investigations and longitudinal monitoring in national and international programs in public health and food safety, as well as for local healthcare systems. These tools also enable research to derive new knowledge regarding disease pathogenesis, and to inform strategies in disease prevention and treatment. As global community resources, the continued sharing of strain genomic data and phenotypes further enhances international resources and means to develop impactful applications. We present examples showing use of these resources in surveillance, including capacity to assess linkages among clinical, environmental, and foodborne reservoirs and to further research investigations into factors promoting their persistence and virulence in different settings.

Paenalcaligenes faecalis sp. nov., a novel species of the family Alcaligenaceae isolated from chicken faeces

A Gram-negative, motile, rod-shaped aerobic and alkalogenic bacterium, designated as strain YLCF04T, was isolated from chicken faeces. Its growth was optimal at 28 °C (range, 10-40 °C), pH 8 (range, pH 6-9) and in 1 % (w/v) NaCl (range, 0-10 %). It was classified to the genus Paenalcaligenes and was most closely related to Paenalcaligenes hominis CCUG 53761AT (97.5 % similarity) based on 16S rRNA gene sequence analysis. Average nucleotide identity and digital DNA-DNA hybridization values between YLCF04T and P. hominis CCUG 53761AT were 76.3 and 18.2 %, respectively. Strain YLCF04T has a genome size of 2.7 Mb with DNA G+C content of 46.3 mol%. Based on its phylogenetic, genomic, phenotypic and biochemical characteristics, strain YLCF04T represents a novel species of the genus Paenalcaligenes, for which the name Paenalcaligenes faecalis sp. nov. is proposed. The type strain is YLCF04T (=CCTCC AB 2022359T= KCTC 92789T).

A panoramic view of the virosphere in three wastewater treatment plants by integrating viral-like particle-concentrated and traditional non-concentrated metagenomic approaches

Wastewater biotreatment systems harbor a rich diversity of microorganisms, and the effectiveness of biotreatment systems largely depends on the activity of these microorganisms. Specifically, viruses play a crucial role in altering microbial behavior and metabolic processes throughout their infection phases, an aspect that has recently attracted considerable interest. Two metagenomic approaches, viral-like particle-concentrated (VPC, representing free viral-like particles) and non-concentrated (NC, representing the cellular fraction), were employed to assess their efficacy in revealing virome characteristics, including taxonomy, diversity, host interactions, lifestyle, dynamics, and functional genes across processing units of three wastewater treatment plants (WWTPs). Our findings indicate that each approach offers unique insights into the viral community and functional composition. Their combined use proved effective in elucidating WWTP viromes. We identified nearly 50,000 viral contigs, with Cressdnaviricota and Uroviricota being the predominant phyla in the VPC and NC fractions, respectively. Notably, two pathogenic viral families, Asfarviridae and Adenoviridae, were commonly found in these WWTPs. We also observed significant differences in the viromes of WWTPs processing different types of wastewater. Additionally, various phage-derived auxiliary metabolic genes (AMGs) were active at the RNA level, contributing to the metabolism of the microbial community, particularly in carbon, sulfur, and phosphorus cycling. Moreover, we identified 29 virus-carried antibiotic resistance genes (ARGs) with potential for host transfer, highlighting the role of viruses in spreading ARGs in the environment. Overall, this study provides a detailed and integrated view of the virosphere in three WWTPs through the application of VPC and NC metagenomic approaches. Our findings enhance the understanding of viral communities, offering valuable insights for optimizing the operation and regulation of wastewater treatment systems.

Dataset of 569 metagenome-assembled genomes from the caeca of multiple chicken breeds from commercial and backyard farming setups of Pakistan

This article focuses the recovery of prokaryotic organisms including bacteria and archaea from 9 different groups of chicken raised in different farm setups in Pakistan. The groups comprise of three different breeds (Broilers, White Layers, and Black Australorp) of chicken raised in different farming setups that include antibiotic-free control, commercial (open and controlled shed), and backyard farms. We have recovered 569 Metagenomics-Assembled Genomes (MAGs) with a completeness of ≥50 % and contamination of ≤10 %. For each MAG, functional annotations were obtained that include KEGG modules, carbohydrate active enzymes (CAZymes), peptidases, geochemical cycles, antibiotic resistance genes, stress genes, and virulence genes. Furthermore, two different sets of Single Copy Genes (SCGs) were used to construct the phylogenetic trees. Based on the reconstructed phylogeny, phylogenetic gain of each MAG is calculated to give an account of novelty.

Faecalibacterium taiwanense sp. nov., isolated from human faeces

Two Gram-stain-negative, straight rods, non-motile, asporogenous, catalase-negative and obligately anaerobic butyrate-producing strains, HLW78T and CYL33, were isolated from faecal samples of two healthy Taiwanese adults. Phylogenetic analyses of 16S rRNA and DNA mismatch repair protein MutL (mutL) gene sequences revealed that these two novel strains belonged to the genus Faecalibacterium. On the basis of 16S rRNA and mutL gene sequence similarities, the type strains Faecalibacterium butyricigenerans AF52-21T(98.3-98.1 % and 79.0-79.5 % similarity), Faecalibacterium duncaniae A2-165T(97.8-97.9 % and 70.9-80.1 %), Faecalibacterium hattorii APC922/41-1T(97.1-97.3 % and 80.3-80.5 %), Faecalibacterium longum CM04-06T(97.8-98.0% and 78.3 %) and Faecalibacterium prausnitzii ATCC 27768T(97.3-97.4 % and 82.7-82.9 %) were the closest neighbours to the novel strains HLW78T and CYL33. Strains HLW78T and CYL33 had 99.4 % both the 16S rRNA and mutL gene sequence similarities, 97.9 % average nucleotide identity (ANI), 96.3 % average amino acid identity (AAI), and 80.5 % digital DNA-DNA hybridization (dDDH) values, indicating that these two strains are members of the same species. Phylogenomic tree analysis indicated that strains HLW78T and CYL33 formed an independent robust cluster together with F. prausnitzii ATCC 27768T. The ANI, AAI and dDDH values between strain HLW78T and its closest neighbours were below the species delineation thresholds of 77.6-85.1 %, 71.4-85.2 % and 28.3-30.9 %, respectively. The two novel strains could be differentiated from the type strains of their closest Faecalibacterium species based on their cellular fatty acid compositions, which contained C18 : 1  ω7c and lacked C15 : 0 and C17 : 1  ω6c, respectively. Phenotypic, chemotaxonomic and genotypic test results demonstrated that the two novel strains HLW78T and CYL33 represented a single, novel species within the genus Faecalibacterium, for which the name Faecalibacterium taiwanense sp. nov. is proposed. The type strain is HLW78T (=BCRC 81397T=NBRC 116372T).

Genetic diversity and antimicrobial susceptibility of Streptococcus suis from diseased Swiss pigs collected between 2019 - 2022

Streptococcus suis is an important pathogen causing severe disease in pigs and humans, giving rise to economic losses in the pig production industry. Out of 65 S. suis isolates collected from diseased pigs in Switzerland between 2019 and 2022, 57 isolates were thoroughly examined by phenotypic and whole genome sequence (WGS) based characterization. The isolates' genomes were sequenced allowing for a comprehensive analysis of their distribution in terms of serovar, sequence type (ST), clonal complex (CC), and classical virulence markers. Antimicrobial resistance (AMR) genes were screened, and phenotypic susceptibility to eight classes of antimicrobial agents was examined. Serovar 6, devoid of any resistance genes, was found to be most prevalent, followed by serovars 1, 3, 1/2, and 9. Thirty STs were identified, with ST1104 being the most prevalent. Serovar 2 and serovar 1/2 were associated with CC1, potentially containing the most virulent variants. Based on single nucleotide polymorphism (SNP) analyses, fifteen isolates belonged to one of seven putative transmission clusters each consisting of two or three isolates. High phenotypic AMR rates were detected for tetracyclines (80%) and macrolides (35%) and associated with the resistance genes tet(O) and erm(B), respectively. In contrast, susceptibility to β-lactam antibiotics and phenicols was high. Determination of phenotypic AMR profiling, including the minimum inhibitory concentrations (MICs) of the tested antimicrobial agents, sets a baseline for future studies. The study provides valuable insights into the genetic diversity and antimicrobial susceptibility of Swiss S. suis isolates, facilitating the identification of emerging clones relevant to public health concerns.

Genomic characterisation of Escherichia coli isolated from poultry at retail through Sink Surveillance in Dhaka, Bangladesh reveals high levels of multi-drug resistance

The surveillance of antimicrobial resistance (AMR) in commensal Escherichia coli from livestock at slaughter is widely employed to assess the potential for risk to humans. There is currently a limited understanding of AMR in Bangladesh poultry at retail in live bird markets, with studies focussing solely on phenotypic characterisation of resistance. To address this evidence gap we performed antimicrobial susceptibility testing and whole genome sequencing on E. coli obtained from chickens from live bird markets in Dhaka in 2018 (n = 38) and 2020 (n = 45). E. coli were isolated from caeca samples following ISO guidelines and sequenced using short and long read methods. Multidrug resistance was extremely common (n = 77) and there was excellent concordance between AMR phenotype and the presence of corresponding AMR genes or mutations. There was considerable genomic diversity, with 43 different sequence types detected. Public health considerations included the high occurrence of resistance to ciprofloxacin (n = 75) associated with plasmid-residing qnrS or mutations in the gyrA and parC chromosomal genes; and the detection of a tigecycline resistant isolate harbouring tet(X4) on an IncHI1A/B-IncFIA mosaic plasmid. Thirty-nine isolates were resistant to azithromycin and harboured mphA, with a significant increase in the incidence of resistance between 2018 and 2020. Although azithromycin is banned for veterinary use in Bangladesh it remains an important treatment option for humans. Interestingly, mphA confers high-level resistance to azithromycin and erythromycin, and the latter is commonly used on poultry farms in Bangladesh. Seven isolates were colistin resistant and carried mcr1. For two isolates hybrid assemblies revealed that mcr1 resided on a highly conserved IncHI2 plasmid that had 93% nucleotide identity to a plasmid from the published genome of an E. coli isolate of Bangladeshi human origin. Six isolates had resistance to third generation cephalosporins, associated with plasmid-residing bla CTX-M-55, bla CTX-M-65, or bla DHA-1. By employing phenotypic and genomic approaches for AMR surveillance we have provided new insights into the potential for One Health AMR linkages in Bangladesh. Employing similar approaches in human and environmental sectors will help inform the One Health approach to addressing AMR, and generate evidence to support mitigation measures such as improved antimicrobial stewardship.

Aeromonas trota Is Highly Refractory to Acquire Exogenous Genetic Material

Aeromonas trota is sensitive to most antibiotics and the sole species of this genus susceptible to ampicillin. This susceptibility profile could be related to its inability to acquire exogenous DNA. In this study, A. trota isolates were analyzed to establish their capacity to incorporate foreign DNA. Fourteen strains were identified as A. trota by multilocus phylogenetic analysis (MLPA). Minimal inhibitory concentrations of antibiotics (MIC) were assessed, confirming the susceptibility to most antibiotics tested. To explore their capacity to be transformed, A. trota strains were used as recipients in different horizontal transfer assays. Results showed that around fifty percent of A. trota strains were able to incorporate pBAMD1-2 and pBBR1MCS-3 plasmids after conjugal transfer. In all instances, conjugation frequencies were very low. Interestingly, several isoforms of plasmid pBBR1MCS-3 were observed in transconjugants. Strains could not receive pAr-32, a native plasmid from A. salmonicida. A. trota strains were unable to receive DNA by means of electroporation, natural transformation or vesiduction. These results confirm that A. trota species are extremely refractory to horizontal gene transfer, which could be associated to plasmid instability resulting from oligomerization or to the presence of defense systems against exogenous genetic material in their genomes. To explain the poor results of horizontal gene transfer (HGT), selected genomes were sequenced and analyzed, revealing the presence of defense systems, which could prevent the stable incorporation of exogenous DNA in A. trota.

Pseudomonas aeruginosa Activates Quorum Sensing, Antioxidant Enzymes and Type VI Secretion in Response to Oxidative Stress to Initiate Biofilm Formation and Wound Chronicity

Pseudomonas aeruginosa (PA) is an opportunistic pathogen frequently isolated from cutaneous chronic wounds. How PA, in the presence of oxidative stress (OS), colonizes chronic wounds and forms a biofilm is still unknown. The purpose of this study is to investigate the changes in gene expression seen when PA is challenged with the high levels of OS present in chronic wounds. We used a biofilm-forming PA strain isolated from the chronic wounds of our murine model (RPA) and performed a qPCR to obtain gene expression patterns as RPA developed a biofilm in vitro in the presence of high levels of OS, and then compared the findings in vivo, in our mouse model of chronic wounds. We found that the planktonic bacteria under OS conditions overexpressed quorum sensing genes that are important for the bacteria to communicate with each other, antioxidant stress genes important to reduce OS in the microenvironment for survival, biofilm formation genes and virulence genes. Additionally, we performed RNAseq in vivo and identified the activation of novel genes/pathways of the Type VI Secretion System (T6SS) involved in RPA pathogenicity. In conclusion, RPA appears to survive the high OS microenvironment in chronic wounds and colonizes these wounds by turning on virulence, biofilm-forming and survival genes. These findings reveal pathways that may be promising targets for new therapies aimed at disrupting PA-containing biofilms immediately after debridement to facilitate the treatment of chronic human wounds.

Multidrug-resistant Escherichia coli causing canine pyometra and urinary tract infections are genetically related but distinct from those causing prostatic abscesses

Despite extensive characterisation of uropathogenic Escherichia coli (UPEC) causing urinary tract infections (UTIs), the genetic background of non-urinary extraintestinal pathogenic E. coli (ExPEC) in companion animals remains inadequately understood. In this study, we characterised virulence traits of 104 E. coli isolated from canine pyometra (n = 61) and prostatic abscesses (PAs) (n = 38), and bloodstream infections (BSIs) in dogs (n = 2), and cats (n = 3). A stronger association with UPEC of pyometra strains in comparison to PA strains was revealed. Notably, 44 isolates exhibited resistance to third-generation cephalosporins and/or fluoroquinolones, 15 were extended-spectrum ß-lactamase-producers. Twelve multidrug-resistant (MDR) strains, isolated from pyometra (n = 4), PAs (n = 5), and BSIs (n = 3), along with 7 previously characterised UPEC strains from dogs and cats, were sequenced. Genomic characteristics revealed that MDR E. coli associated with UTIs, pyometra, and BSIs belonged to international high-risk E. coli clones, including sequence type (ST) 38, ST131, ST617, ST648, and ST1193. However, PA strains belonged to distinct lineages, including ST12, ST44, ST457, ST744, and ST13037. The coreSNPs, cgMLST, and pan-genome illustrated intra-clonal variations within the same ST from different sources. The high-risk ST131 and ST1193 (phylogroup B2) contained high numbers of ExPEC virulence genes on pathogenicity islands, predominating in pyometra and UTI. Hybrid MDR/virulence IncF multi-replicon plasmids, containing aerobactin genes, were commonly found in non-B2 phylogroups from all sources. These findings offer genomic insights into non-urinary ExPEC, highlighting its potential for invasive infections in pets beyond UTIs, particularly with regards to high-risk global clones.

A one health approach for monitoring antimicrobial resistance: developing a national freshwater pilot effort

Antimicrobial resistance (AMR) is a world-wide public health threat that is projected to lead to 10 million annual deaths globally by 2050. The AMR public health issue has led to the development of action plans to combat AMR, including improved antimicrobial stewardship, development of new antimicrobials, and advanced monitoring. The National Antimicrobial Resistance Monitoring System (NARMS) led by the United States (U.S) Food and Drug Administration along with the U.S. Centers for Disease Control and U.S. Department of Agriculture has monitored antimicrobial resistant bacteria in retail meats, humans, and food animals since the mid 1990's. NARMS is currently exploring an integrated One Health monitoring model recognizing that human, animal, plant, and environmental systems are linked to public health. Since 2020, the U.S. Environmental Protection Agency has led an interagency NARMS environmental working group (EWG) to implement a surface water AMR monitoring program (SWAM) at watershed and national scales. The NARMS EWG divided the development of the environmental monitoring effort into five areas: (i) defining objectives and questions, (ii) designing study/sampling design, (iii) selecting AMR indicators, (iv) establishing analytical methods, and (v) developing data management/analytics/metadata plans. For each of these areas, the consensus among the scientific community and literature was reviewed and carefully considered prior to the development of this environmental monitoring program. The data produced from the SWAM effort will help develop robust surface water monitoring programs with the goal of assessing public health risks associated with AMR pathogens in surface water (e.g., recreational water exposures), provide a comprehensive picture of how resistant strains are related spatially and temporally within a watershed, and help assess how anthropogenic drivers and intervention strategies impact the transmission of AMR within human, animal, and environmental systems.

Genomic Characterization of Carbapenemase-Producing Enterobacter hormaechei, Serratia marcescens, Citrobacter freundii, Providencia stuartii, and Morganella morganii Clinical Isolates from Bulgaria

Carbapenemase-producing Enterobacter spp. Serratia marcescens, Citrobacter freundii, Providencia spp., and Morganella morganii (CP-ESCPM) are increasingly identified as causative agents of nosocomial infections but are still not under systematic genomic surveillance. In this study, using a combination of whole-genome sequencing and conjugation experiments, we sought to elucidate the genomic characteristics and transferability of resistance genes in clinical CP-ESCPM isolates from Bulgaria. Among the 36 sequenced isolates, NDM-1 (12/36), VIM-4 (11/36), VIM-86 (8/36), and OXA-48 (7/36) carbapenemases were identified; two isolates carried both NDM-1 and VIM-86. The majority of carbapenemase genes were found on self-conjugative plasmids. IncL plasmids were responsible for the spread of OXA-48 among E. hormaechei, C. freundii, and S. marcescens. IncM2 plasmids were generally associated with the spread of NDM-1 in C. freundii and S. marcescens, and also of VIM-4 in C. freundii. IncC plasmids were involved in the spread of the recently described VIM-86 in P. stuartii isolates. IncC plasmids carrying blaNDM-1 and blaVIM-86 were observed too. blaNDM-1 was also detected on IncX3 in S. marcescens and on IncT plasmid in M. morganii. The significant resistance transfer rates we observed highlight the role of the ESCPM group as a reservoir of resistance determinants and stress the need for strengthening infection control measures.

AMRViz enables seamless genomics analysis and visualization of antimicrobial resistance

We have developed AMRViz, a toolkit for analyzing, visualizing, and managing bacterial genomics samples. The toolkit is bundled with the current best practice analysis pipeline allowing researchers to perform comprehensive analysis of a collection of samples directly from raw sequencing data with a single command line. The analysis results in a report showing the genome structure, genome annotations, antibiotic resistance and virulence profile for each sample. The pan-genome of all samples of the collection is analyzed to identify core- and accessory-genes. Phylogenies of the whole genome as well as all gene clusters are also generated. The toolkit provides a web-based visualization dashboard allowing researchers to interactively examine various aspects of the analysis results. Availability: AMRViz is implemented in Python and NodeJS, and is publicly available under open source MIT license at https://github.com/amromics/amrviz .

Demographic fluctuations in bloodstream Staphylococcus aureus lineages configure the mobile gene pool and antimicrobial resistance

Staphylococcus aureus in the bloodstream causes high morbidity and mortality, exacerbated by the spread of multidrug-resistant and methicillin-resistant S. aureus (MRSA). We aimed to characterize the circulating lineages of S. aureus from bloodstream infections and the contribution of individual lineages to resistance over time. Here, we generated 852 high-quality short-read draft genome sequences of S. aureus isolates from patient blood cultures in a single hospital from 2010 to 2022. A total of 80 previously recognized sequence types (ST) and five major clonal complexes are present in the population. Two frequently detected lineages, ST5 and ST8 exhibited fluctuating demographic structures throughout their histories. The rise and fall in their population growth coincided with the acquisition of antimicrobial resistance, mobile genetic elements, and superantigen genes, thus shaping the accessory genome structure across the entire population. These results reflect undetected selective events and changing ecology of multidrug-resistant S. aureus in the bloodstream.

Bloodstream Infection Caused by Erysipelothrix rhusiopathiae in an Immunocompetent Patient

Erysipelothrix rhusiopathiae is a facultative anaerobe Gram-positive bacillus, which is considered a zoonotic pathogen. E. rhusiopathiae causes erysipeloid, mainly in occupational groups such as veterinarians, slaughterhouse workers, farmers, and fishermen. Two cutaneous forms (localised and generalised) and a septicaemic form have been described. Here, we report the isolation of a strain of E. rhusiopathiae from a 56-year-old immunocompetent obese male admitted to Fondazione IRCCS Policlinico San Matteo Pavia (Italy). Blood cultures were collected and Gram-positive bacilli were observed. E. rhusiopathiae grew and was identified. Antimicrobial susceptibility tests were performed and interpreted with EUCAST breakpoints (PK-PD). The strain was susceptible to all the antibiotics tested, while it was intrinsically resistant to vancomycin. The clinical diagnosis of E. rhusiopathiae can be challenging, due to the broad spectrum of symptoms and potential side effects, including serious systemic infections such as heart diseases. In the case described, bacteraemia caused by E. rhusiopathiae was detected in a immunocompetent patient. Bacteraemia caused by E. rhusiopathiae is rare in immunocompetent people and blood cultures were proven to be essential for the diagnosis and underdiagnosis of this pathogen, which is possible due to its resemblance to other clinical manifestations.