MOB-suite: software tools for clustering, reconstruction and typing of plasmids from draft assemblies

Transmission of Klebsiella strains and plasmids within and between grey-headed flying fox colonies

The grey-headed flying fox (Pteropus poliocephalus) is an endemic Australian fruit bat, known to carry zoonotic pathogens. We recently showed they harbour bacterial pathogen Klebsiella pneumoniae and closely related species in the K. pneumoniae species complex (KpSC); however, the dynamics of KpSC transmission and gene flow within flying fox colonies are poorly understood. High-resolution genome comparisons of 39 KpSC isolates from grey-headed flying foxes identified five putative strain transmission clusters (four intra- and one inter-colony). The instance of inter-colony strain transmission of K. africana was found between two flying fox populations within flying distance, indicating either direct or indirect transmission through a common food/water source. All 11 plasmids identified within the KpSC isolates showed 73% coverage (mean) and ≥95% identity to human-associated KpSC plasmids, indicating gene flow between human clinical and grey-headed flying fox isolates. Along with strain transmission, inter-species horizontal plasmid transmission between K. pneumoniae and Klebsiella africana was also identified within a flying fox colony. Finally, genome-scale metabolic models were generated to predict and compare substrate usage to previously published KpSC models, from human and environmental sources. These models indicated no distinction on the basis of metabolic capabilities. Instead, metabolic capabilities were consistent with population structure and ST/lineage.

WGS-Based Lineage and Antimicrobial Resistance Pattern of Salmonella Typhimurium Isolated during 2000-2017 in Peru

Salmonella Typhimurium is associated with foodborne diseases worldwide, including in Peru, and its emerging antibiotic resistance (AMR) is now a global public health problem. Therefore, country-specific monitoring of the AMR emergence is vital to control this pathogen, and in these aspects, whole genome sequence (WGS)—based approaches are better than gene-based analyses. Here, we performed the antimicrobial susceptibility test for ten widely used antibiotics and WGS-based various analyses of 90 S. Typhimurium isolates (human, animal, and environment) from 14 cities of Peru isolated from 2000 to 2017 to understand the lineage and antimicrobial resistance pattern of this pathogen in Peru. Our results suggest that the Peruvian isolates are of Typhimurium serovar and predominantly belong to sequence type ST19. Genomic diversity analyses indicate an open pan-genome, and at least ten lineages are circulating in Peru. A total of 48.8% and 31.0% of isolates are phenotypically and genotypically resistant to at least one antibiotic, while 12.0% are multi-drug resistant (MDR). Genotype−phenotype correlations for ten tested drugs show >80% accuracy, and >90% specificity. Sensitivity above 90% was only achieved for ciprofloxacin and ceftazidime. Two lineages exhibit the majority of the MDR isolates. A total of 63 different AMR genes are detected, of which 30 are found in 17 different plasmids. Transmissible plasmids such as lncI-gamma/k, IncI1-I(Alpha), Col(pHAD28), IncFIB, IncHI2, and lncI2 that carry AMR genes associated with third-generation antibiotics are also identified. Finally, three new non-synonymous single nucleotide variations (SNVs) for nalidixic acid and eight new SNVs for nitrofurantoin resistance are predicted using genome-wide association studies, comparative genomics, and functional annotation. Our analysis provides for the first time the WGS-based details of the circulating S. Typhimurium lineages and their antimicrobial resistance pattern in Peru.

Capturing the antibiotic resistome of preterm infants reveals new benefits of probiotic supplementation

BACKGROUND

Probiotic use in preterm infants can mitigate the impact of antibiotic exposure and reduce rates of certain illnesses; however, the benefit on the gut resistome, the collection of antibiotic resistance genes, requires further investigation. We hypothesized that probiotic supplementation of early preterm infants (born < 32-week gestation) while in hospital reduces the prevalence of antibiotic resistance genes associated with pathogenic bacteria in the gut. We used a targeted capture approach to compare the resistome from stool samples collected at the term corrected age of 40 weeks for two groups of preterm infants (those that routinely received a multi-strain probiotic during hospitalization and those that did not) with samples from full-term infants at 10 days of age to identify if preterm birth or probiotic supplementation impacted the resistome. We also compared the two groups of preterm infants up to 5 months of age to identify persistent antibiotic resistance genes.

RESULTS

At the term corrected age, or 10 days of age for the full-term infants, we found over 80 antibiotic resistance genes in the preterm infants that did not receive probiotics that were not identified in either the full-term or probiotic-supplemented preterm infants. More genes associated with antibiotic inactivation mechanisms were identified in preterm infants unexposed to probiotics at this collection time-point compared to the other infants. We further linked these genes to mobile genetic elements and Enterobacteriaceae, which were also abundant in their gut microbiomes. Various genes associated with aminoglycoside and beta-lactam resistance, commonly found in pathogenic bacteria, were retained for up to 5 months in the preterm infants that did not receive probiotics.

CONCLUSIONS

This pilot survey of preterm infants shows that probiotics administered after preterm birth during hospitalization reduced the diversity and prevented persistence of antibiotic resistance genes in the gut microbiome. The benefits of probiotic use on the microbiome and the resistome should be further explored in larger groups of infants. Due to its high sensitivity and lower sequencing cost, our targeted capture approach can facilitate these surveys to further address the implications of resistance genes persisting into infancy without the need for large-scale metagenomic sequencing. Video Abstract.

ESBL plasmids in Klebsiella pneumoniae: diversity, transmission and contribution to infection burden in the hospital setting

Background

Resistance to third-generation cephalosporins, often mediated by extended-spectrum beta-lactamases (ESBLs), is a considerable issue in hospital-associated infections as few drugs remain for treatment. ESBL genes are often located on large plasmids that transfer horizontally between strains and species of Enterobacteriaceae and frequently confer resistance to additional drug classes. Whilst plasmid transmission is recognised to occur in the hospital setting, the frequency and impact of plasmid transmission on infection burden, compared to ESBL + strain transmission, is not well understood.

Methods

We sequenced the genomes of clinical and carriage isolates of Klebsiella pneumoniae species complex from a year-long hospital surveillance study to investigate ESBL burden and plasmid transmission in an Australian hospital. Long-term persistence of a key transmitted ESBL + plasmid was investigated via sequencing of ceftriaxone-resistant isolates during 4 years of follow-up, beginning 3 years after the initial study.

Results

We found 25 distinct ESBL plasmids. We identified one plasmid, which we called Plasmid A, that carried bla_CTX-M-15 in an IncF backbone similar to pKPN-307. Plasmid A was transmitted at least four times into different Klebsiella species/lineages and was responsible for half of all ESBL episodes during the initial 1-year study period. Three of the Plasmid A-positive strains persisted locally 3–6 years later, and Plasmid A was detected in two additional strain backgrounds. Overall Plasmid A accounted for 21% of ESBL + infections in the follow-up period.

Conclusions

Here, we systematically surveyed ESBL strain and plasmid transmission over 1 year in a single hospital network. Whilst ESBL plasmid transmission events were rare in this setting, they had a significant and sustained impact on the burden of ceftriaxone-resistant and multidrug-resistant infections. If onward transmission of Plasmid A-carrying strains could have been prevented, this may have reduced the number of opportunities for Plasmid A to transmit and create novel ESBL + strains, as well as reducing overall ESBL infection burden.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01103-0.

Comparative genomics of dairy-associated Staphylococcus aureus from selected sub-Saharan African regions reveals milk as reservoir for human-and animal-derived strains and identifies a putative animal-related clade with presumptive novel siderophore

Staphylococcus aureus infection is considered to be a neglected tropical disease with huge impact on human and animal health alike. Dairy production in sub-Saharan Africa (SSA) relies heavily on various animals such as cows, goats, and camels, depending on the region. S. aureus causes mastitis and exhibits high prevalence in raw milk. The population structure including genotypic and phenotypic traits of dairy S. aureus in relation to animal and human isolates is, however, unknown for SSA. In this work, 20 S. aureus dairy isolates from East and West Africa were selected for comparative genomics and phenotypic analysis. Comparing their population structure revealed a large diversity of different origins suggesting milk to be a reservoir for human and animal strains alike. Furthermore, a novel putative siderophore was detected in multiple strains in a distinct animal-clade with strains of global origin. This putative siderophore shares a high genetic identity with that from Streptococcus equi suggesting possible horizontal gene transfer. These findings combined with the virulence genes harbored by these dairy-derived strains such as pvl, human evasion factor scn, various enterotoxin, leucocidin and antibiotic resistance genes, stresses the need for an integrative One Health approach to tackle the problem of S. aureus infections in animals and humans in sub-Saharan Africa.

Prevalence and distribution of extended-spectrum β-lactamase and AmpC-producing Escherichia coli in two New Zealand dairy farm environments

Antimicrobial resistance (AMR) is a global threat to human and animal health, with the misuse and overuse of antimicrobials being suggested as the main driver of resistance. In a global context, New Zealand (NZ) is a relatively low user of antimicrobials in animal production. However, the role antimicrobial usage on pasture-based dairy farms, such as those in NZ, plays in driving the spread of AMR within the dairy farm environment remains equivocal. Culture-based methods were used to determine the prevalence and distribution of extended-spectrum β-lactamase (ESBL)- and AmpC-producing Escherichia coli from farm environmental samples collected over a 15-month period from two NZ dairy farms with contrasting management practices. Whole genome sequencing was utilised to understand the genomic epidemiology and antimicrobial resistance gene repertoire of a subset of third-generation cephalosporin resistant E. coli isolated in this study. There was a low sample level prevalence of ESBL-producing E. coli (faeces 1.7%; farm dairy effluent, 6.7% from Dairy 4 and none from Dairy 1) but AmpC-producing E. coli were more frequently isolated across both farms (faeces 3.3% and 8.3%; farm dairy effluent 38.4%, 6.7% from Dairy 1 and Dairy 4, respectively). ESBL- and AmpC-producing E. coli were isolated from faeces and farm dairy effluent in spring and summer, during months with varying levels of antimicrobial use, but no ESBL- or AmpC-producing E. coli were isolated from bulk tank milk or soil from recently grazed paddocks. Hybrid assemblies using short- and long-read sequence data from a subset of ESBL- and AmpC-producing E. coli enabled the assembly and annotation of nine plasmids from six E. coli, including one plasmid co-harbouring 12 antimicrobial resistance genes. ESBL-producing E. coli were infrequently identified from faeces and farm dairy effluent on the two NZ dairy farms, suggesting they are present at a low prevalence on these farms. Plasmids harbouring several antimicrobial resistance genes were identified, and bacteria carrying such plasmids are a concern for both animal and public health. AMR is a burden for human, animal and environmental health and requires a holistic “One Health” approach to address.

Staphylococcus aureus lineages associated with a free-ranging population of the fruit bat Pteropus livingstonii retained over 25 years in captivity

Conservation of endangered species has become increasingly complex, and costly interventions to protect wildlife require a robust scientific evidence base. This includes consideration of the role of the microbiome in preserving animal health. Captivity introduces stressors not encountered in the wild including environmental factors and exposure to exotic species, humans and antimicrobial drugs. These stressors may perturb the microbiomes of wild animals, with negative consequences for their health and welfare and hence the success of the conservation project, and ultimately the risk of release of non-native organisms into native ecosystems. We compared the genomes of Staphylococcus aureus colonising critically endangered Livingstone’s fruit bats (Pteropus livingstonii) which have been in a captive breeding programme for 25 years, with those from bats in the endemic founder population free ranging in the Comoros Republic. Using whole genome sequencing, we compared 47 isolates from captive bats with 37 isolates from those free ranging in the Comoros Republic. Our findings demonstrate unexpected resilience in the bacteria carried, with the captive bats largely retaining the same two distinctive lineages carried at the time of capture. In addition, we found evidence of genomic changes which suggest specific adaptations to the bat host.

Dominant Carbapenemase-Encoding Plasmids in Clinical Enterobacterales Isolates and Hypervirulent Klebsiella pneumoniae, Singapore

Dissemination of carbapenemase-encoding plasmids by horizontal gene transfer in multidrug-resistant bacteria is the major driver of rising carbapenem-resistance, but the conjugative mechanics and evolution of clinically relevant plasmids are not yet clear. We performed whole-genome sequencing on 1,215 clinical Enterobacterales isolates collected in Singapore during 2010–2015. We identified 1,126 carbapenemase-encoding plasmids and discovered pKPC2 is becoming the dominant plasmid in Singapore, overtaking an earlier dominant plasmid, pNDM1. pKPC2 frequently conjugates with many Enterobacterales species, including hypervirulent Klebsiella pneumoniae, and maintains stability in vitro without selection pressure and minimal adaptive sequence changes. Furthermore, capsule and decreasing taxonomic relatedness between donor and recipient pairs are greater conjugation barriers for pNDM1 than pKPC2. The low fitness costs pKPC2 exerts in Enterobacterales species indicate previously undetected carriage selection in other ecological settings. The ease of conjugation and stability of pKPC2 in hypervirulent K. pneumoniae could fuel spread into the community.

Differential Overlap in Human and Animal Fecal Microbiomes and Resistomes in Rural versus Urban Bangladesh

While the development of antibiotic resistance in animal gut microbiomes and subsequent transmission to humans has been demonstrated in intensive farming environments and high-income countries, evidence of zoonotic exchange of antibiotic resistance in LMIC communities is lacking. This research provides genomic evidence of overlap of antibiotic resistance genes between humans and animals, especially in urban communities, and highlights chickens as important reservoirs of antibiotic resistance.

Genomic surveillance for multidrug-resistant or hypervirulent Klebsiella pneumoniae among United States bloodstream isolates

BACKGROUND

Klebsiella pneumoniae strains have been divided into two major categories: classical K. pneumoniae, which are frequently multidrug-resistant and cause hospital-acquired infections in patients with impaired defenses, and hypervirulent K. pneumoniae, which cause severe community-acquired and disseminated infections in normal hosts. Both types of infections may lead to bacteremia and are associated with significant morbidity and mortality. The relative burden of these two types of K. pneumoniae among bloodstream isolates within the United States is not well understood.

METHODS

We evaluated consecutive K. pneumoniae isolates cultured from the blood of hospitalized patients at Northwestern Memorial Hospital (NMH) in Chicago, Illinois between April 2015 and April 2017. Bloodstream isolates underwent whole genome sequencing, and sequence types (STs), capsule loci (KLs), virulence genes, and antimicrobial resistance genes were identified in the genomes using the bioinformatic tools Kleborate and Kaptive. Patient demographic, comorbidity, and infection information, as well as the phenotypic antimicrobial resistance of the isolates were extracted from the electronic health record. Candidate hypervirulent isolates were tested in a murine model of pneumonia, and their plasmids were characterized using long-read sequencing. We also extracted STs, KLs, and virulence and antimicrobial resistance genes from the genomes of bloodstream isolates submitted from 33 United States institutions between 2007 and 2021 to the National Center for Biotechnology Information (NCBI) database.

RESULTS

Consecutive K. pneumoniae bloodstream isolates (n = 104, one per patient) from NMH consisted of 75 distinct STs and 51 unique capsule loci. The majority of these isolates (n = 58, 55.8%) were susceptible to all tested antibiotics except ampicillin, but 17 (16.3%) were multidrug-resistant. A total of 32 (30.8%) of these isolates were STs of known high-risk clones, including ST258 and ST45. In particular, 18 (17.3%) were resistant to ceftriaxone (of which 17 harbored extended-spectrum beta-lactamase genes) and 9 (8.7%) were resistant to meropenem (all of which harbored a carbapenemase genes). Four (3.8%) of the 104 isolates were hypervirulent K. pneumoniae, as evidenced by hypermucoviscous phenotypes, high levels of virulence in a murine model of pneumonia, and the presence of large plasmids similar to characterized hypervirulence plasmids. These isolates were cultured from patients who had not recently traveled to Asia. Two of these hypervirulent isolates belonged to the well characterized ST23 lineage and one to the re-emerging ST66 lineage. Of particular concern, two of these isolates contained plasmids with tra conjugation loci suggesting the potential for transmission. We also analyzed 963 publicly available genomes of K. pneumoniae bloodstream isolates from locations within the United States. Of these, 465 (48.3%) and 760 (78.9%) contained extended-spectrum beta-lactamase genes or carbapenemase genes, respectively, suggesting a bias towards submission of antibiotic-resistant isolates. The known multidrug-resistant high-risk clones ST258 and ST307 were the predominant sequence types. A total of 32 (3.3%) of these isolates contained aerobactin biosynthesis genes and 26 (2.7%) contained at least two genetic features of hvKP strains, suggesting elevated levels of virulence. We identified 6 (0.6%) isolates that were STs associated with hvKP: ST23 (n = 4), ST380 (n = 1), and ST65 (n = 1).

CONCLUSIONS

Examination of consecutive isolates from a single center demonstrated that multidrug-resistant high-risk clones are indeed common, but a small number of hypervirulent K. pneumoniae isolates were also observed in patients with no recent travel history to Asia, suggesting that these isolates are undergoing community spread in the United States. A larger collection of publicly available bloodstream isolate genomes also suggested that hypervirulent K. pneumoniae strains are present but rare in the USA; however, this collection appears to be heavily biased towards highly antibiotic-resistant isolates (and correspondingly away from hypervirulent isolates).

Population structure of blaKPC-harbouring IncN plasmids at a New York City medical centre and evidence for multi-species horizontal transmission

BACKGROUND

Carbapenem-resistant Enterobacterales (CRE) are highly concerning MDR pathogens. Horizontal transfer of broad-host-range IncN plasmids may contribute to the dissemination of the Klebsiella pneumoniae carbapenemase (KPC), spreading carbapenem resistance among unrelated bacteria. However, the population structure and genetic diversity of IncN plasmids has not been fully elucidated.

OBJECTIVES

We reconstructed blaKPC-harbouring IncN plasmid genomes to characterize shared gene content, structural variability, and putative horizontal transfer within and across patients and diverse bacterial clones.

METHODS

We performed short- and long-read sequencing and hybrid assembly on 45 CRE isolates with blaKPC-harbouring IncN plasmids. Eight serial isolates from two patients were included to assess intra-patient plasmid dynamics. Comparative genomic analysis was performed to assess structural and sequence similarity across plasmids. Within IncN sublineages defined by plasmid MLST and kmer-based clustering, phylogenetic analysis was used to identify closely related plasmids.

RESULTS

Comparative analysis of IncN plasmid genomes revealed substantial heterogeneity including large rearrangements in serial patient plasmids and differences in structure and content across plasmid clusters. Within plasmid sublineages, core genome content and resistance gene regions were largely conserved. Closely related plasmids (≤1 SNP) were found in highly diverse isolates, including ten pST6 plasmids found in eight bacterial clones from three different species.

CONCLUSIONS

Genomic analysis of blaKPC-harbouring IncN plasmids revealed the presence of several distinct sublineages as well as substantial host diversity within plasmid clusters suggestive of frequent mobilization. This study reveals complex plasmid dynamics within a single plasmid family, highlighting the challenge of tracking plasmid-mediated transmission of blaKPC in clinical settings.

Listeria ilorinensis sp. nov., isolated from cow milk cheese in Nigeria

During microbial assessment of cow milk cheese products in the city of Ilorin, Nigeria, a Listeria -like isolate was detected that could not be assigned to any known species. Whole-genome sequence analyses against all currently known 26 Listeria species confirmed that this isolate constitutes a new taxon within the genus Listeria , with highest similarity to Listeria costaricensis (average nucleotide identity blast of 82.66%, in silico DNA–DNA hybridization of 28.3%). Phenotypically, it differs from L. costaricensis by the inability to ferment sucrose, l-fucose and starch. The absence of haemolysis and Listeria pathogenic islands suggest that this novel species is not pathogenic for humans and animals. The name Listeria ilorinensis sp. nov. is proposed, with the type strain CLIP 2019/01311T (=CIP 111875T=DSM 111566T).

Molecular Evolution and Genomic Insights into Community-Acquired Methicillin-Resistant Staphylococcus aureus Sequence Type 88

Sequence type 88 (ST88) methicillin-resistant Staphylococcus aureus (MRSA) has been recognized as an important pathogen that causes infections in humans, especially when it has strong biofilm production and multidrug resistance (MDR). However, knowledge of the determinants of resistance or virulence and genomic characteristics of ST88 MRSA from China is still limited. In this study, we employed the antimicrobial resistance (AMR), biofilm formation, and genomic characteristics of ST88 MRSA collected from various foods in China and estimated the worldwide divergence of ST88 MRSA with publicly available ST88 genomes. All ST88 isolates studied were identified as having resistance genes, while 50% (41/82) harbored MDR genes. All isolates carried core virulence genes related to immune modulation, adherence, secreted enzymes, and hemolysin. In addition, all 20 Chinese ST88 isolates showed biofilm production capacity, three strongly so. Bayesian phylogenetic analysis showed that Chinese ST88 clones formed an independent MRSA lineage, with two subclades associated with acquisition of type IVc staphylococcal cassette chromosome mec (SCCmec) elements. In contrast, all African ST88 strains were subtyped as SCCmecIVa, where the African clades were mixed with a few European and American isolates, suggesting potential transmission from Africa to these regions. In summary, our results revealed the evolution of ST88 MRSA in humans, animals, and foods in Africa and Asia. The food-associated ST88 genomes in this study will remedy the lack of food-associated ST88 isolates, and the study in general will extend the discussion of the potential exchanges of ST88 between humans and foods or food animals.

IMPORTANCE ST88 MRSA has frequently been detected in humans, animals, and foods mainly in Africa and Asia. It can colonize and cause mild to severe infections in humans, especially children. Several studies from African countries have described its genotypic characteristics but, limited information is available on the evolution and characterization of ST88 MRSA in Asia, especially China. Meanwhile, the molecular history of its global spread remains largely unclear. In this study, we analyzed the genomic evolution of global ST88 MRSA strains in detail and identified key genetic changes associated with specific hosts or regions. Our results suggested geographical differentiation between ST88 MRSA’s evolution in Africa and its evolution in Asia, with a more recent clonal evolution in China. The introduction of ST88 MRSA in China was aligned with the acquisition of SCCmecIVc elements, specific virulent prophages, and AMR genes.

High diversity of pathogenic Escherichia coli clones carrying mcr-1 among gulls underlines the need for strategies at the environment-livestock-human interface

The expansion of mcr-carrying bacteria is a well-recognized public health problem. Measures to contain mcr spread have mainly been focused on the food-animal production sector. Nevertheless, the spread of MCR-producers at the environmental interface particularly driven by the increasing population of gulls in coastal cities has been less explored. Occurrence of mcr-carrying Escherichia coli in gull's colonies faeces on a Portuguese beach was screened over 7-months. Cultural, molecular, and genomic approaches were used to characterize their diversity, mcr plasmids and adaptive features. Multidrug-resistant mcr-1-carrying E. coli were detected for three consecutive months. Over time, multiple strains were recovered, including zoonotic-related pathogenic E. coli clones (e.g., B2-ST131-H22, A-ST10, and B1-ST162). Diverse mcr-1.1 genetic environments were mainly associated with ST2/ST4-HI2 (ST10, ST131, ST162, ST354 and ST4204) but also IncI2 (ST12990) plasmids or in the chromosome (ST656). Whole-genome sequencing revealed enrichment of these strains on antibiotic resistance, virulence, and metal tolerance genes. Our results underscore gulls as important spreaders of high priority bacteria and genes that may affect the environment, food-animals and/or humans, potentially undermining One-Health strategies to reduce colistin resistance. This article is protected by copyright. All rights reserved.

A 16th century Escherichia coli draft genome associated with an opportunistic bile infection

Escherichia coli – one of the most characterized bacteria and a major public health concern – remains invisible across the temporal landscape. Here, we present the meticulous reconstruction of the first ancient E. coli genome from a 16^th century gallstone from an Italian mummy with chronic cholecystitis. We isolated ancient DNA and reconstructed the ancient E. coli genome. It consisted of one chromosome of 4446 genes and two putative plasmids with 52 genes. The E. coli strain belonged to the phylogroup A and an exceptionally rare sequence type 4995. The type VI secretion system component genes appears to be horizontally acquired from Klebsiella aerogenes, however we could not identify any pathovar specific genes nor any acquired antibiotic resistances. A sepsis mouse assay showed that a closely related contemporary E. coli strain was avirulent. Our reconstruction of this ancient E. coli helps paint a more complete picture of the burden of opportunistic infections of the past.

Ancient DNA from an Italian mummy’s gallstone provides insight into opportunistic E. coli infection.

Draft Genome Sequences of Two Clostridium botulinum Group II Strains Carrying Phage-Like Plasmids

Clostridium botulinum is responsible for botulism, a potentially lethal foodborne intoxication. Here, we report the draft genome sequences of C. botulinum group II strains 202F (serotype F) and Hazen (serotype E). The genomes share many similarities, including multiple mobile genetic elements.

Comparative in silico characterization of Klebsiella pneumoniae hypervirulent plasmids and their antimicrobial resistance genes

Background

The hypervirulent pathotype of Klebsiella pneumoniae (hvKp) is mainly mediated by large virulent plasmids. It seems that these hypervirulent plasmids (HVPs) are accumulating antimicrobial resistance genes (ARGs) and are turning quickly into drug-resistant hypervirulent hybrids. Therefore, molecular mechanisms involved in this convergence needs to be investigated to control their global spread.

Methods

In this study, the complete sequence of 79 non-redundant hypervirulent plasmids were retrieved from GenBank and their genetic features, hypervirulence and antimicrobial resistance patterns (AMR) as well as their putative transmission capability were compared using bioinformatics tools.

Results

The majority of HVPs belonged to clonal complex (CC)23, and sequence type (ST)11. IncFIB and IncHI1B were the most prevalent plasmid replicon types. Out of 79 plasmids, 78 were positive for iutA and iucA. The iucC, iucB and iucD genes were found in 77 plasmids. Almost 26% of the HVPs were potentially conjugative of which 71% carried AGRs. ARGs against beta-lactams, carbapenems, quinolones, aminoglycosides, chloramphenicols, tetracyclines and macrolides were detected in 30% of HVPs. Class 1 integron and prophage structures harboring multiple ARGs were found in eight plasmids. Insertion sequences (IS)6, IS110 and IS1380 appeared to be important genetic elements in transmission of ARGs.

Conclusions

The high prevalence of iucA and iutA suggests their strong capability for rapid and accurate genetic markers for discrimination of hvKp in the laboratory. This study indicated the important role of mobile genetic elements (MGEs) in the emergence of drug-resistance in hypervirulent strains. The high prevalence of putative conjugative hybrids implies higher incidence of multidrug-resistant (MDR)-hvKp strains in near future.

Comprehensive pan-genomic, resistome and virulome analysis of clinical OXA-48 producing carbapenem-resistant Serratia marcescens strains

BACKGROUND

Carbapenem-resistant Enterobacteriaceae (CRE) have been thoroughly studied as the pathogens associated with hospital acquired infections. However, data on Serratia marcescens are not enough. S. marcescens is now becoming a propensity for its highly antimicrobial-resistant clinical infections.

METHODS

Four carbapenem-resistant S. marcescens (CR-SM) isolates were obtained from hospitalized patients through routine microbiological experiments. We assembled the isolates genomes using whole genome sequencing (WGS) and compared their resistome and virulome patterns.

RESULTS

The average length and CG content of chromosomes was 5.33 Mbp and 59.8%, respectively. The number of coding sequences (CDSs) ranged from 4,959 to 4,989. All strains had one single putative conjugative plasmid with IncL incompatibility (Inc) group. The strains harbored bla_CTX-M-15, bla_TEM-1 and bla_SHV-134. All plamsids were positive for bla_OXA-48. No bla_NDM-1, bla_KPC, bla_VIM and bla_IMP were identified. The bla_SRT-2 and aac(6')-Ic genes were chromosomally-encoded. Class 1 integron was detected in strains P8, P11 and P14. The Escher_RCS47 and Salmon_SJ46 prophages played major role in plasmid-mediated carraige of extended spectrum β-lactamases (ESBLs). The CR-SM strains were equipt with typical virulence factors of oppotunistic pathogens including biofilm formation, adhesins, secretory systems and siderophores. The strains did not have ability to produce prodigiosin but were positive for chitinase and EstA.

CONCLUSION

The presence of conjugative plasmids harboring major β-lactamases within prophage and class 1 integron structures highlights the role of different mobile genetic elements (MGEs) in distribution of AMR factors and more specifically carbapenemases. More molecular studies are required to determine the status of carbapenem resistance in clinical starins. However, appropriate strategies to control the global dissemination of CR-SM are urgent.

CRISPR-Cas systems are widespread accessory elements across bacterial and archaeal plasmids

Many prokaryotes encode CRISPR-Cas systems as immune protection against mobile genetic elements (MGEs), yet a number of MGEs also harbor CRISPR-Cas components. With a few exceptions, CRISPR-Cas loci encoded on MGEs are uncharted and a comprehensive analysis of their distribution, prevalence, diversity, and function is lacking. Here, we systematically investigated CRISPR-Cas loci across the largest curated collection of natural bacterial and archaeal plasmids. CRISPR-Cas loci are widely but heterogeneously distributed across plasmids and, in comparison to host chromosomes, their mean prevalence per Mbp is higher and their distribution is distinct. Furthermore, the spacer content of plasmid CRISPRs exhibits a strong targeting bias towards other plasmids, while chromosomal arrays are enriched with virus-targeting spacers. These contrasting targeting preferences highlight the genetic independence of plasmids and suggest a major role for mediating plasmid-plasmid conflicts. Altogether, CRISPR-Cas are frequent accessory components of many plasmids, which is an overlooked phenomenon that possibly facilitates their dissemination across microbiomes.

Mobility of antimicrobial resistance across serovars and disease presentations in non-typhoidal Salmonella from animals and humans in Vietnam

Non-typhoidal Salmonella (NTS) is a major cause of bacterial enterocolitis globally but also causes invasive bloodstream infections. Antimicrobial resistance (AMR) hampers the treatment of these infections and understanding how AMR spreads between NTS may help in developing effective strategies. We investigated NTS isolates associated with invasive disease, diarrhoeal disease and asymptomatic carriage in animals and humans from Vietnam. Isolates included multiple serovars and both common and rare phenotypic AMR profiles; long- and short-read sequencing was used to investigate the genetic mechanisms and genomic backgrounds associated with phenotypic AMR profiles. We demonstrate concordance between most AMR genotypes and phenotypes but identified large genotypic diversity in clinically relevant phenotypes and the high mobility potential of AMR genes (ARGs) in this setting. We found that 84 % of ARGs identified were located on plasmids, most commonly those containing IncHI1A_1 and IncHI1B(R27)_1_R27 replicons (33%), and those containing IncHI2_1 and IncHI2A_1 replicons (31%). The vast majority (95%) of ARGS were found within 10 kbp of IS6/IS26 elements, which provide plasmids with a mechanism to exchange ARGs between plasmids and other parts of the genome. Whole genome sequencing with targeted long-read sequencing applied in a One Health context identified a comparatively limited number of insertion sequences and plasmid replicons associated with AMR. Therefore, in the context of NTS from Vietnam and likely for other settings as well, the mechanisms by which ARGs move contribute to a more successful AMR profile than the specific ARGs, facilitating the adaptation of bacteria to different environments or selection pressures.

The Current Landscape of Antibiotic Resistance of Salmonella Infantis in Italy: The Expansion of Extended-Spectrum Beta-Lactamase Producers on a Local Scale

Salmonella enterica serovar Infantis is one of the five main causes of human salmonellosis in the European Union (EU) and in recent years, has been increasingly reported to carry multiple antimicrobial resistance determinants, including extended-spectrum beta-lactamase (ESBL) genes. In our study, we used WGS-based tools to characterize S. Infantis strains circulating in the Abruzzo and Molise regions of Italy between 2017 and 2020 and compared this local dataset to the S. Infantis population present in Italy over the last two decades. Phylogenetic analyses demonstrated that the majority of strains isolated from poultry and turkeys from Abruzzo and Molise were closely related and belonged to one of the two main genetic clusters present in Italy, which were grouped predominantly as ESBL-producing strains that harbored pESI-like plasmid. We showed that 60% of the local strains carried multiple antibiotic resistance genes, including ESBL gene bla_CTX–M–1 as well as aadA1, dfrA1, dfrA14, sul1, and tet(A) genes present on the pESI-like megaplasmid. The analysis of strains from Abruzzo and Molise and the publicly available Italian S. Infantis sequences revealed a dramatic increase in the number of identified AMR genes in the strains isolated after 2011. Moreover, the number of strains resistant to five or more antibiotic classes increased from 20–80% in the last decade likely due to the acquisition of the megaplasmid. The persistence of the ESBL-producing and the multidrug-resistant (MDR) clone of S. Infantis in poultry populations in Italy and in Europe requires rapid and efficient intervention strategies to prevent further expansion of the clone.

Precise genotyping of circular mobile elements from metagenomic data uncovers human-associated plasmids with recent common ancestors

Mobile genetic elements with circular genomes play a key role in the evolution of microbial communities. Their circular genomes correspond to circular walks in metagenome graphs, and yet, assemblies derived from natural microbial communities produce graphs riddled with spurious cycles, complicating the accurate reconstruction of circular genomes. We present DomCycle, an algorithm that reconstructs likely circular genomes based on the identification of so-called 'dominant' graph cycles. In the implementation we leverage paired reads to bridge assembly gaps and scrutinize cycles through a nucleotide-level analysis, making the approach robust to misassembly artifacts. We validated the approach using simulated and real sequencing data. Application of DomCycle to 32 publicly available DNA shotgun sequence data sets from diverse natural environments led to the reconstruction of hundreds of circular mobile genomes. Clustering revealed 20 highly prevalent and cryptic plasmids that have clonal population structures with recent common ancestors. This method facilitates the study of microbial communities that evolve through horizontal gene transfer.

Pathogenomes and variations in Shiga toxin production among geographically distinct clones of Escherichia coli O113:H21

Infections with globally disseminated Shiga toxin-producing Escherichia coli (STEC) of the O113:H21 serotype can progress to severe clinical complications, such as hemolytic uremic syndrome (HUS). Two phylogeographically distinct clonal complexes have been established by multi locus sequence typing (MLST). Infections with ST-820 isolates circulating exclusively in Australia have caused severe human disease, such as HUS. Conversely, ST-223 isolates prevalent in the US and outside Australia seem to rarely cause severe human disease but are frequent contaminants. Following a genomic epidemiology approach, we wanted to gain insights into the underlying cause for this disparity. We examined the plasticity in the genome make-up and Shiga toxin production in a collection of 20 ST-820 and ST-223 strains isolated from produce, the bovine reservoir, and clinical cases. STEC are notorious for assembly into fragmented draft sequences when using short-read sequencing technologies due to the extensive and partly homologous phage complement. The application of long-read technology (LRT) sequencing yielded closed reference chromosomes and plasmids for two representative ST-820 and ST-223 strains. The established high-resolution framework, based on whole genome alignments, single nucleotide polymorphism (SNP)-typing and MLST, includes the chromosomes and plasmids of other publicly available O113:H21 sequences and allowed us to refine the phylogeographical boundaries of ST-820 and ST-223 complex isolates and to further identify a historic non-shigatoxigenic strain from Mexico as a quasi-intermediate. Plasmid comparison revealed strong correlations between the strains' featured pO113 plasmid genotypes and chromosomally inferred ST, which suggests coevolution of the chromosome and virulence plasmids. Our pathogenicity assessment revealed statistically significant differences in the Stx_2a-production capabilities of ST-820 as compared to ST-223 strains under RecA-induced Stx phage mobilization, a condition that mimics Stx-phage induction. These observations suggest that ST-820 strains may confer an increased pathogenic potential in line with the strain-associated epidemiological metadata. Still, some of the tested ST-223 cultures sourced from contaminated produce or the bovine reservoir also produced Stx at levels comparable to those of ST-820 isolates, which calls for awareness and for continued surveillance of this lineage.

Inter-species geographic signatures for tracing horizontal gene transfer and long-term persistence of carbapenem resistance

BACKGROUND

Carbapenem-resistant Enterobacterales (CRE) are an urgent global health threat. Inferring the dynamics of local CRE dissemination is currently limited by our inability to confidently trace the spread of resistance determinants to unrelated bacterial hosts. Whole-genome sequence comparison is useful for identifying CRE clonal transmission and outbreaks, but high-frequency horizontal gene transfer (HGT) of carbapenem resistance genes and subsequent genome rearrangement complicate tracing the local persistence and mobilization of these genes across organisms.

METHODS

To overcome this limitation, we developed a new approach to identify recent HGT of large, near-identical plasmid segments across species boundaries, which also allowed us to overcome technical challenges with genome assembly. We applied this to complete and near-complete genome assemblies to examine the local spread of CRE in a systematic, prospective collection of all CRE, as well as time- and species-matched carbapenem-susceptible Enterobacterales, isolated from patients from four US hospitals over nearly 5 years.

RESULTS

Our CRE collection comprised a diverse range of species, lineages, and carbapenem resistance mechanisms, many of which were encoded on a variety of promiscuous plasmid types. We found and quantified rearrangement, persistence, and repeated transfer of plasmid segments, including those harboring carbapenemases, between organisms over multiple years. Some plasmid segments were found to be strongly associated with specific locales, thus representing geographic signatures that make it possible to trace recent and localized HGT events. Functional analysis of these signatures revealed genes commonly found in plasmids of nosocomial pathogens, such as functions required for plasmid retention and spread, as well survival against a variety of antibiotic and antiseptics common to the hospital environment.

CONCLUSIONS

Collectively, the framework we developed provides a clearer, high-resolution picture of the epidemiology of antibiotic resistance importation, spread, and persistence in patients and healthcare networks.

Nationwide surveillance in Thailand revealed genotype-dependent dissemination of carbapenem-resistant Enterobacterales

Carbapenem-resistant Enterobacterales (CRE) are a serious public health threat because of their rapid dissemination. To determine the epidemiological and genetic characteristics of CRE infections in Thailand, we performed whole-genome sequencing of 577 carbapenem-resistant Klebsiella pneumoniae isolates and 170 carbapenem-resistant Escherichia coli isolates from hospitals across the nation. The four most prevalent carbapenemase genes harboured by these bacteria were bla NDM-1, bla NDM-5, bla OXA-181 and bla OXA-232. The gene bla NDM-1 was identified in diverse sequence types. The gene bla NDM-5 was identified almost exclusively in E. coli . The genes bla OXA-181, bla OXA-232, and co-carriage of bla NDM-1 and bla OXA-232 were found in specific sequence types from certain provinces. Replicon typing revealed the diverse backbones of bla NDM-1- and bla NDM-5-harbouring plasmids and successful expansion of bla NDM-1-harbouring IncN2-type plasmids. Core-genome single-nucleotide polymorphism analysis suggested that bla OXA-181-, bla OXA-232-, bla NDM-5-, and co-carriage of bla NDM-1 and bla OXA-232-associated sub-clonal lineages have recently predominated in the provinces from where these isolates were isolated. Thus, we demonstrate genotype-dependent dissemination of CRE in Thailand, which is helpful for establishing infection-control strategies in CRE-endemic areas.

The Space-Exposed Kombucha Microbial Community Member Komagataeibacter oboediens Showed Only Minor Changes in Its Genome After Reactivation on Earth

Komagataeibacter is the dominant taxon and cellulose-producing bacteria in the Kombucha Microbial Community (KMC). This is the first study to isolate the K. oboediens genome from a reactivated space-exposed KMC sample and comprehensively characterize it. The space-exposed genome was compared with the Earth-based reference genome to understand the genome stability of K. oboediens under extraterrestrial conditions during a long time. Our results suggest that the genomes of K. oboediens IMBG180 (ground sample) and K. oboediens IMBG185 (space-exposed) are remarkably similar in topology, genomic islands, transposases, prion-like proteins, and number of plasmids and CRISPR-Cas cassettes. Nonetheless, there was a difference in the length of plasmids and the location of cas genes. A small difference was observed in the number of protein coding genes. Despite these differences, they do not affect any genetic metabolic profile of the cellulose synthesis, nitrogen-fixation, hopanoid lipids biosynthesis, and stress-related pathways. Minor changes are only observed in central carbohydrate and energy metabolism pathways gene numbers or sequence completeness. Altogether, these findings suggest that K. oboediens maintains its genome stability and functionality in KMC exposed to the space environment most probably due to the protective role of the KMC biofilm. Furthermore, due to its unaffected metabolic pathways, this bacterial species may also retain some promising potential for space applications.

Genomic characterization of multidrug-resistant Salmonella serovar Kentucky ST198 isolated in poultry flocks in Spain (2011-2017)

Salmonella Kentucky is commonly found in poultry and rarely associated with human disease. However, a multidrug-resistant (MDR) S. Kentucky clone [sequence type (ST)198] has been increasingly reported globally in humans and animals. Our aim here was to assess if the recently reported increase of S. Kentucky in poultry in Spain was associated with the ST198 clone and to characterize this MDR clone and its distribution in Spain. Sixty-six isolates retrieved from turkey, laying hen and broiler in 2011–2017 were subjected to whole-genome sequencing to assess their sequence type, genetic relatedness, and presence of antimicrobial resistance genes (ARGs), plasmid replicons and virulence factors. Thirteen strains were further analysed using long-read sequencing technologies to characterize the genetic background associated with ARGs. All isolates belonged to the ST198 clone and were grouped in three clades associated with the presence of a specific point mutation in the gyrA gene, their geographical origin and isolation year. All strains carried between one and 16 ARGs whose presence correlated with the resistance phenotype to between two and eight antimicrobials. The ARGs were located in the Salmonella genomic island (SGI-1) and in some cases (bla_SHV-12, catA1, cmlA1, dfrA and multiple aminoglycoside-resistance genes) in IncHI2/IncI1 plasmids, some of which were consistently detected in different years/farms in certain regions, suggesting they could persist over time. Our results indicate that the MDR S. Kentucky ST198 is present in all investigated poultry hosts in Spain, and that certain strains also carry additional plasmid-mediated ARGs, thus increasing its potential public health significance.

Whole-Genome Sequences of Two Listeria monocytogenes Biofilm Formers

Listeria monocytogenes strains from different lineages show different biofilm-forming abilities. In this study, two strains of L. monocytogenes were whole genome sequenced using single-molecule real-time (SMRT) technology and characterized.

The Bacillus cereus Strain EC9 Primes the Plant Immune System for Superior Biocontrol of Fusarium oxysporum

Antibiosis is a key feature widely exploited to develop biofungicides based on the ability of biological control agents (BCAs) to produce fungitoxic compounds. A less recognised attribute of plant-associated beneficial microorganisms is their ability to stimulate the plant immune system, which may provide long-term, systemic self-protection against different types of pathogens. By using conventional antifungal in vitro screening coupled with in planta assays, we found antifungal and non-antifungal Bacillus strains that protected the ornamental plant Kalanchoe against the soil-borne pathogen Fusarium oxysporum in experimental and commercial production settings. Further examination of one antifungal and one non-antifungal strain indicated that high protection efficacy in planta did not correlate with antifungal activity in vitro. Whole-genome sequencing showed that the non-antifungal strain EC9 lacked the biosynthetic gene clusters associated with typical antimicrobial compounds. Instead, this bacterium triggers the expression of marker genes for the jasmonic and salicylic acid defence pathways, but only after pathogen challenge, indicating that this strain may protect Kalanchoe plants by priming immunity. We suggest that the stimulation of the plant immune system is a promising mode of action of BCAs for the development of novel biological crop protection products.

Antimicrobial Resistance in Rivers: A Review of the Genes Detected and New Challenges

River ecosystems are very important parts of the water cycle and an excellent habitat, food, and drinking water source for many organisms, including humans. Antibiotics are emerging contaminants which can enter rivers from various sources. Several antibiotics and their related antibiotic resistance genes (ARGs) have been detected in these ecosystems by various research programs and could constitute a substantial problem. The presence of antibiotics and other resistance cofactors can boost the development of ARGs in the chromosomes or mobile genetic elements of natural bacteria in rivers. The ARGs in environmental bacteria can also be transferred to clinically important pathogens. However, antibiotics and their resistance genes are both not currently monitored by national or international authorities responsible for controlling the quality of water bodies. For example, they are not included in the contaminant list in the European Water Framework Directive or in the US list of Water-Quality Benchmarks for Contaminants. Although ARGs are naturally present in the environment, very few studies have focused on non-impacted rivers to assess the background ARG levels in rivers, which could provide some useful indications for future environmental regulation and legislation. The present study reviews the antibiotics and associated ARGs most commonly measured and detected in rivers, including the primary analysis tools used for their assessment. In addition, other factors that could enhance antibiotic resistance, such as the effects of chemical mixtures, the effects of climate change, and the potential effects of the coronavirus disease 2019 pandemic, are discussed. Environ Toxicol Chem 2022;41:687-714. © 2022 SETAC.

Whole-genome sequencing and gene sharing network analysis powered by machine learning identifies antibiotic resistance sharing between animals, humans and environment in livestock farming

Anthropogenic environments such as those created by intensive farming of livestock, have been proposed to provide ideal selection pressure for the emergence of antimicrobial-resistant Escherichia coli bacteria and antimicrobial resistance genes (ARGs) and spread to humans. Here, we performed a longitudinal study in a large-scale commercial poultry farm in China, collecting E. coli isolates from both farm and slaughterhouse; targeting animals, carcasses, workers and their households and environment. By using whole-genome phylogenetic analysis and network analysis based on single nucleotide polymorphisms (SNPs), we found highly interrelated non-pathogenic and pathogenic E. coli strains with phylogenetic intermixing, and a high prevalence of shared multidrug resistance profiles amongst livestock, human and environment. Through an original data processing pipeline which combines omics, machine learning, gene sharing network and mobile genetic elements analysis, we investigated the resistance to 26 different antimicrobials and identified 361 genes associated to antimicrobial resistance (AMR) phenotypes; 58 of these were known AMR-associated genes and 35 were associated to multidrug resistance. We uncovered an extensive network of genes, correlated to AMR phenotypes, shared among livestock, humans, farm and slaughterhouse environments. We also found several human, livestock and environmental isolates sharing closely related mobile genetic elements carrying ARGs across host species and environments. In a scenario where no consensus exists on how antibiotic use in the livestock may affect antibiotic resistance in the human population, our findings provide novel insights into the broader epidemiology of antimicrobial resistance in livestock farming. Moreover, our original data analysis method has the potential to uncover AMR transmission pathways when applied to the study of other pathogens active in other anthropogenic environments characterised by complex interconnections between host species.

Livestock have been suggested as an important source of antimicrobial-resistant (AMR) Escherichia coli, capable of infecting humans and carrying resistance to drugs used in human medicine. China has a large intensive livestock farming industry, poultry being the second most important source of meat in the country, and is the largest user of antibiotics for food production in the world. Here we studied antimicrobial resistance gene overlap between E. coli isolates collected from humans, livestock and their shared environments in a large-scale Chinese poultry farm and associated slaughterhouse. By using a computational approach that integrates machine learning, whole-genome sequencing, gene sharing network and mobile genetic elements analysis we characterized the E. coli community structure, antimicrobial resistance phenotypes and the genetic relatedness of non-pathogenic and pathogenic E. coli strains. We uncovered the network of genes, associated with AMR, shared across host species (animals and workers) and environments (farm and slaughterhouse). Our approach opens up new avenues for the development of a fast, affordable and effective computational solutions that provide novel insights into the broader epidemiology of antimicrobial resistance in livestock farming.

Use of genomics to explore AMR persistence in an outdoor pig farm with low antimicrobial usage

Food animals may be reservoirs of antimicrobial resistance (AMR) passing through the food chain, but little is known about AMR prevalence in bacteria when selective pressure from antimicrobials is low or absent. We monitored antimicrobial-resistant Escherichia coli over 1 year in a UK outdoor pig farm with low antimicrobial usage (AMU) compared to conventional pig farms in the United Kingdom. Short and selected long-read whole-genome sequencing (WGS) was performed to identify AMR genes, phylogeny and mobile elements in 385 E. coli isolates purified mainly from pig and some seagull faeces. Generally, low levels of antimicrobial-resistant E. coli were present, probably due to low AMU. Those present were likely to be multi-drug resistant (MDR) and belonging to particular Sequence Types (STs) such as ST744, ST88 or ST44, with shared clones (<14 Single Nucleotide Polymorphisms (SNPs) apart) isolated from different time points indicating epidemiological linkage within pigs of different ages, and between pig and the wild bird faeces. Although importance of horizontal transmission of AMR is well established, there was limited evidence of plasmid-mediated dissemination between different STs. Non-conjugable MDR plasmids or large AMR gene-bearing transposons were stably integrated within the chromosome and remained associated with particular STs/clones over the time period sampled. Heavy metal resistance genes were also detected within some genetic elements. This study highlights that although low levels of antimicrobial-resistant E. coli correlates with low AMU, a basal level of MDR E. coli can still persist on farm potentially due to transmission and recycling of particular clones within different pig groups. Environmental factors such as wild birds and heavy metal contaminants may also play important roles in the recycling and dissemination, and hence enabling persistence of MDR E. coli. All such factors need to be considered as any rise in AMU on low usage farms, could in future, result in a significant increase in their AMR burden.

Host ecology regulates interspecies recombination in bacteria of the genus Campylobacter

Horizontal gene transfer (HGT) can allow traits that have evolved in one bacterial species to transfer to another. This has potential to rapidly promote new adaptive trajectories such as zoonotic transfer or antimicrobial resistance. However, for this to occur requires gaps to align in barriers to recombination within a given time frame. Chief among these barriers is the physical separation of species with distinct ecologies in separate niches. Within the genus Campylobacter, there are species with divergent ecologies, from rarely isolated single-host specialists to multihost generalist species that are among the most common global causes of human bacterial gastroenteritis. Here, by characterizing these contrasting ecologies, we can quantify HGT among sympatric and allopatric species in natural populations. Analyzing recipient and donor population ancestry among genomes from 30 Campylobacter species, we show that cohabitation in the same host can lead to a six-fold increase in HGT between species. This accounts for up to 30% of all SNPs within a given species and identifies highly recombinogenic genes with functions including host adaptation and antimicrobial resistance. As described in some animal and plant species, ecological factors are a major evolutionary force for speciation in bacteria and changes to the host landscape can promote partial convergence of distinct species through HGT.

Whole-Genome Sequences of Vibrio Species from Warm-Water Shrimps Imported into Canada: Detection of Genetic Elements Associated with Antimicrobial Resistance and Potential Mobilizing Capacities

We present draft genome sequences of Vibrio species (Vibrio alginolyticus, Vibrio cholerae, and two Vibrio parahaemolyticus strains) that were isolated from warm-water shrimps imported into Canada. All four isolates harbor genetic elements associated with antimicrobial resistance (AMR), including mobile genetic elements that can promote horizontal transfer of AMR genes.

Closed Genome Sequences of Providencia alcalifaciens Isolates from Dogs

Eight Providencia alcalifaciens isolates from eight different dogs in Norway with acute hemorrhagic diarrhea were sequenced. Based on Illumina and Oxford Nanopore Technologies sequencing, all of the genomes were complete and closed after hybrid assembly.

Draft Genome Sequences of Two Lactobacillus johnsonii and Three Ligilactobacillus salivarius Strains Isolated from Intestinal Microbiomes of Chickens

This report describes the genome sequences of two Lactobacillus johnsonii strains (AER105 and AER25) and three Ligilactobacillus salivarius strains (AER35, AER36, and AER04) recovered from broiler chicken gastrointestinal tracts in the southeastern United States. These genome sequences will enhance our understanding of the ecology of lactobacilli in the chicken gut microbiome.

Virulence Genotype and Phenotype of Multiple Antimicrobial-Resistant Escherichia coli Isolates from Broilers Assessed from a "One-Health" Perspective

ABSTRACT

Extraintestinal pathogenic Escherichia coli (ExPEC) include several serotypes that have been associated with colibacillosis in poultry and with urinary tract infections (UTIs) and newborn meningitis in humans. In this study, 57 antimicrobial-resistant E. coli from apparently healthy broiler chickens were characterized for their health and safety risks. These isolates belonged to 12 serotypes, and isolates of the same serotype were clonal based on single nucleotide variant analysis. Most of the isolates harbored plasmids; IncC and IncFIA were frequently detected. The majority of the resistant isolates harbored plasmid-mediated resistance genes, including aph(3″)-Ib, aph(6)-Id, blaCMY-2, floR, sul1, sul2, tet(A), and tet(B), in agreement with their resistant phenotypes. The class 1 integron was detected in all E. coli serotypes except O124:H25 and O7:H6. Of the 57 broiler E. coli isolates, 27 were avian pathogenic, among which 18 were also uropathogenic E. coli and the remainder were other ExPEC. The two isolates of serotype O161:H4 (ST117) were genetically related to the control avian pathogenic strains and a clinical isolate associated with UTIs. A strain of serotype O159:H45 (ST101) also was closely related to a UTI isolate. The detected virulence factors included adhesins, invasins, siderophores, type III secretion systems, and toxins in combination with other virulence determinants. A broiler isolate of serotype O7:H18 (ST38) carried the ibeA gene encoding a protein involved in invasion of brain endothelium on a 102-kbp genetic island. This isolate moderately adhered and invaded Caco-2 cells and induced mortality (42.5%) in a day-old-chick infection model. The results of this study suggest that multiple antimicrobial-resistant E. coli isolates recovered from apparent healthy broilers can be pathogenic and act as reservoirs for antimicrobial resistance genes, highlighting the necessity of their assessment in a "One-Heath" context.

HIGHLIGHTS

Complete Genome Sequence and Benzophenone-3 Mineralisation Potential of Rhodococcus sp. USK10, A Bacterium Isolated from Riverbank Sediment

Benzophenone-3 (BP3) is an organic UV filter whose presence in the aquatic environment has been linked to detrimental developmental impacts in aquatic organisms such as coral and fish. The genus Rhodococcus has been extensively studied and is known for possessing large genomes housing genes for biodegradation of a wide range of compounds, including aromatic carbons. Here, we present the genome sequence of Rhodococcus sp. USK10, which was isolated from Chinese riverbank sediment and is capable of utilising BP3 as the sole carbon source, resulting in full BP3 mineralisation. The genome consisted of 9,870,030 bp in 3 replicons, a G+C content of 67.2%, and 9722 coding DNA sequences (CDSs). Annotation of the genome revealed that 179 of these CDSs are involved in the metabolism of aromatic carbons. The complete genome of Rhodococcus sp. USK10 is the first complete, annotated genome sequence of a Benzophenone-3-degrading bacterium. Through radiolabelling, it is also the first bacterium proven to mineralise Benzophenone-3. Due to the widespread environmental prevalence of Benzophenone-3, coupled with its adverse impact on aquatic organisms, this characterisation provides an integral first step in better understanding the environmentally relevant degradation pathway of the commonly used UV filter. Given USK10′s ability to completely mineralise Benzophenone-3, it could prove to be a suitable candidate for bioremediation application.

Development of a Genomics-Based Approach To Identify Putative Hypervirulent Nontyphoidal Salmonella Isolates: Salmonella enterica Serovar Saintpaul as a Model

While differences in human virulence have been reported across nontyphoidal Salmonella (NTS) serovars and associated subtypes, a rational and scalable approach to identify Salmonella subtypes with differential ability to cause human diseases is not available. Here, we used NTS serovar Saintpaul (S. Saintpaul) as a model to determine if metadata and associated whole-genome sequence (WGS) data in the NCBI Pathogen Detection (PD) database can be used to identify (i) subtypes with differential likelihoods of causing human diseases and (ii) genes and single nucleotide polymorphisms (SNPs) potentially responsible for such differences. S. Saintpaul SNP clusters (n = 211) were assigned different epidemiology types (epi-types) based on statistically significant over- or underrepresentation of human clinical isolates, including human associated (HA; n = 29), non-human associated (NHA; n = 23), and other (n = 159). Comparative genomic analyses identified 384 and 619 genes overrepresented among isolates in 5 HA and 4 NHA SNP clusters most significantly associated with the respective isolation source. These genes included 5 HA-associated virulence genes previously reported to be present on Gifsy-1/Gifsy-2 prophages. Additionally, premature stop codons in 3 and 7 genes were overrepresented among the selected HA and NHA SNP clusters, respectively. Tissue culture experiments with strains representing 4 HA and 3 NHA SNP clusters did not reveal evidence for enhanced invasion or intracellular survival for HA strains. However, the presence of sodCI (encoding a superoxide dismutase), found in 4 HA and 1 NHA SNP clusters, was positively correlated with intracellular survival in macrophage-like cells. Post hoc analyses also suggested a possible difference in intracellular survival among S. Saintpaul lineages. IMPORTANCE Not all Salmonella isolates are equally likely to cause human disease, and Salmonella control strategies may unintentionally focus on serovars and subtypes with high prevalence in source populations but are rarely associated with human clinical illness. We describe a framework leveraging WGS data in the NCBI PD database to identify Salmonella subtypes over- and underrepresented among human clinical cases. While we identified genomic signatures associated with HA/NHA SNP clusters, tissue culture experiments failed to identify consistent phenotypic characteristics indicative of enhanced human virulence of HA strains. Our findings illustrate the challenges of defining hypo- and hypervirulent S. Saintpaul and potential limitations of phenotypic assays when evaluating human virulence, for which in vivo experiments are essential. Identification of sodCI, an HA-associated virulence gene associated with enhanced intracellular survival, however, illustrates the potential of the framework and is consistent with prior work identifying specific genomic features responsible for enhanced or reduced virulence of nontyphoidal Salmonella.

The prevalence and characterization of ESBL-producing strains of Salmonella enterica circulating in the territory of the Russian Federation (2016–2020)

Objective. To analyze frequency and identify genetic determinants of resistance of non-typhoid Salmonella (NTS) producing extended-spectrum β-lactamase (ESBL) isolated in the Russian Federation over the period 2016 to 2020. Materials and Methods. Salmonella isolates, suspected to ESBL production, were collected by the All-Russia Reference Center of Salmonellosis during the national Salmonellosis surveillance program. Phenotypic resistance was determined by the broth microdilution method using G-I and G-II Mikrolatest®SensiLaTest MIC plates and by the double-disk synergy test. Whole genome sequencing was performed on the NextSeq platform (Illumina, USA), with subsequent de novo genome assembly (SPAdes 3.15.4), identification of plasmid types (MOB-suite v3.0.0), and identification of resistance genes (AMRFinderPlus v3.10.40). Results. Out of 1792 NTS isolates, 22 strains contained bla-genes of molecular classes A and D (blaTEM, blaCTX-M, blaSHV, blaOXA), one strain – AmpC (blaCMY-2) and three strains – combination ESBL of class A and AmpC (blaTEM, blaCMY-2, blaDHA). The frequency of occurrence of ESBL-producing Salmonella is 1.3%, AmpC – 0.2%. Additionally, strains were resistant to other non-β-lactam antibiotics. Six different types of plasmids were identified (IncI, IncFIB, IncC, IncHI2A, IncL/M and IncX1) in studied isolates. It was possible for 17 strains to identify location of resistance genes in plasmids of a certain type. Conclusions. The frequency of occurrence of Salmonella strains producing ESBL and AmpC was 1.45%, which were found in sporadic cases of human diseases, as well as food and environmental objects were sources of isolation. The fact of detection of such strains among various NTC serotypes and a wide range of sources of isolation confirms the relevance of monitoring antimicrobial resistance of Salmonella strains in the future.

A large-scale genomic snapshot of Klebsiella spp. isolates in Northern Italy reveals limited transmission between clinical and non-clinical settings

The Klebsiella group, found in humans, livestock, plants, soil, water and wild animals, is genetically and ecologically diverse. Many species are opportunistic pathogens and can harbour diverse classes of antimicrobial resistance genes. Healthcare-associated Klebsiella pneumoniae clones that are non-susceptible to carbapenems can spread rapidly, representing a high public health burden. Here we report an analysis of 3,482 genome sequences representing 15 Klebsiella species sampled over a 17-month period from a wide range of clinical, community, animal and environmental settings in and around the Italian city of Pavia. Northern Italy is a hotspot for hospital-acquired carbapenem non-susceptible Klebsiella and thus a pertinent setting to examine the overlap between isolates in clinical and non-clinical settings. We found no genotypic or phenotypic evidence for non-susceptibility to carbapenems outside the clinical environment. Although we noted occasional transmission between clinical and non-clinical settings, our data point to a limited role of animal and environmental reservoirs in the human acquisition of Klebsiella spp. We also provide a detailed genus-wide view of genomic diversity and population structure, including the identification of new groups.

Longitudinal Sampling Reveals Persistence of and Genetic Diversity in Extended-Spectrum Cephalosporin-Resistant Escherichia coli From Norwegian Broiler Production

There are knowledge gaps concerning dynamics of extended-spectrum cephalosporin (ESC)-resistant Escherichia coli and their plasmids in broiler production and the persistence of strains on broiler farms. Thus, we aimed at characterising ESC-resistant Escherichia coli collected from all flocks reared on 10 different farms during a six-months sampling period. All isolates (n = 43) were subjected to whole-genome sequencing, and a subset of isolates (n = 7) were also sequenced using oxford nanopore technology and subsequent hybrid assembly in order to do in-depth characterisation of the ESC resistance plasmids. The 43 isolates belonged to 11 different sequence types, and three different ESC resistance gene/plasmid combinations were present, namely, IncK2/bla_CMY-2 (n = 29), IncI1/bla_CMY-2 (n = 6) and IncI1/bla_CTX-M-1 (n = 8). ESC-resistant E. coli of different STs and with different ESC resistance gene/plasmid combinations could be present on the same farm, while a single ST and ESC resistance gene/plasmid displaying zero or few SNP differences were present on other farms. In-depth characterisation of IncK2/bla_CMY-2 plasmids revealed that at least two distinct variants circulate in the broiler production. These plasmids showed close homology to previously published plasmids from other countries. Our longitudinal study show that ESC-resistant E. coli belong to a multitude of different STs and that different ESC resistance genes and plasmids occur. However, there is also indication of persistence of both ESC-resistant E. coli strains and IncK2/bla_CMY-2 plasmids on farms. Further studies are warranted to determine the dynamics of strains, plasmids and ESC resistance genes within single broiler flocks.

Whole-Genome Sequencing-Based Antimicrobial Resistance Characterization and Phylogenomic Investigation of 19 Multidrug-Resistant and Extended-Spectrum Beta-Lactamase-Positive Escherichia coli Strains Collected From Hospital Patients in Benin in 2019

The increasing worldwide prevalence of extended-spectrum beta-lactamase (ESBL) producing Escherichia coli constitutes a serious threat to global public health. Surgical site infections are associated with high morbidity and mortality rates in developing countries, fueled by the limited availability of effective antibiotics. We used whole-genome sequencing (WGS) to evaluate antimicrobial resistance and the phylogenomic relationships of 19 ESBL-positive E. coli isolates collected from surgical site infections in patients across public hospitals in Benin in 2019. Isolates were identified by MALDI-TOF mass spectrometry and phenotypically tested for susceptibility to 16 antibiotics. Core-genome multi-locus sequence typing and single-nucleotide polymorphism-based phylogenomic methods were used to investigate the relatedness between samples. The broader phylogenetic context was characterized through the inclusion of publicly available genome data. Among the 19 isolates, 13 different sequence types (STs) were observed, including ST131 (n = 2), ST38 (n = 2), ST410 (n = 2), ST405 (n = 2), ST617 (n = 2), and ST1193 (n = 2). The bla _CTX-M-15 gene encoding ESBL resistance was found in 15 isolates (78.9%), as well as other genes associated with ESBL, such as bla _OXA-1 (n = 14) and bla _TEM-1 (n = 9). Additionally, we frequently observed genes encoding resistance against aminoglycosides [aac-(6')-Ib-cr, n = 14], quinolones (qnrS₁ , n = 4), tetracyclines [tet(B), n = 14], sulfonamides (sul2, n = 14), and trimethoprim (dfrA17, n = 13). Nonsynonymous chromosomal mutations in the housekeeping genes parC and gyrA associated with resistance to fluoroquinolones were also detected in multiple isolates. Although the phylogenomic investigation did not reveal evidence of hospital-acquired transmissions, we observed two very similar strains collected from patients in different hospitals. By characterizing a set of multidrug-resistant isolates collected from a largely unexplored environment, this study highlights the added value for WGS as an effective early warning system for emerging pathogens and antimicrobial resistance.

Molecular Analysis of blaKPC-2-Harboring Plasmids: Tn4401a Interplasmid Transposition and Tn4401a-Carrying ColRNAI Plasmid Mobilization from Klebsiella pneumoniae to Citrobacter europaeus and Morganella morganii in a Single Patient

The spread of Klebsiella pneumoniae carbapenemase (KPC)-producing Enterobacterales is a public health concern. KPC-encoding blaKPC is predominantly spread by strains of a particular phylogenetic lineage, clonal group 258, but can also be spread by horizontal transfer of blaKPC-carrying plasmids. Here, we report the transfer of a blaKPC-2-harboring plasmid via mobilization from K. pneumoniae to Citrobacter freundii complex and Morganella morganii strains in a single patient. We performed draft whole-genome sequencing to analyze 20 carbapenemase-producing Enterobacterales strains (15 of K. pneumoniae, two of C. freundii complex, and three of M. morganii) and all K. pneumoniae strains using MiSeq and/or MinION isolated from a patient who was hospitalized in New York and Montreal before returning to Japan. All strains harbored blaKPC-2-containing Tn4401a. The 15 K. pneumoniae strains each belonged to sequence type 258 and harbored a Tn4401a-carrying multireplicon-type plasmid, IncN and IncR (IncN+R). Three of these K. pneumoniae strains also possessed a Tn4401a-carrying ColRNAI plasmid, suggesting that Tn4401a underwent interplasmid transposition. Of these three ColRNAI plasmids, two and one were identical to plasmids harbored by two Citrobacter europaeus and three M. morganii strains, respectively. The Tn4401a-carrying ColRNAI plasmids were each 23,753 bp long and incapable of conjugal transfer via their own genes alone, but they mobilized during the conjugal transfer of Tn4401a-carrying IncN+R plasmids in K. pneumoniae. Interplasmid transposition of Tn4401a from an IncN+R plasmid to a ColRNAI plasmid in K. pneumoniae and mobilization of Tn4401a-carrying ColRNAI plasmids contributed to the acquisition of blaKPC-2 in C. europaeus and M. morganii. IMPORTANCE Plasmid transfer plays an important role in the interspecies spread of carbapenemase genes, including the Klebsiella pneumoniae carbapenemase (KPC)-coding gene, blaKPC. We conducted whole-genome sequencing (WGS) analysis and transmission experiments to analyze blaKPC-2-carrying mobile genetic elements (MGEs) between the blaKPC-2-harboring K. pneumoniae, Citrobacter europaeus, and Morganella morganii strains isolated from a single patient. blaKPC-2 was contained within an MGE, Tn4401a. WGS of blaKPC-2-carrying K. pneumoniae, C. europaeus, and M. morganii strains isolated from one patient revealed that Tn4401a-carrying ColRNAI plasmids were generated by plasmid-to-plasmid transfer of Tn4401a from a multireplicon-type IncN and IncR (IncN+R) plasmid in K. pneumoniae strains. Tn4401a-carrying ColRNAI plasmids were incapable of conjugal transfer in C. europaeus and M. morganii but mobilized from K. pneumoniae to a recipient Escherichia coli strain during the conjugal transfer of Tn4401a-carrying IncN+R plasmid. Therefore, Tn4401a-carrying ColRNAI plasmids contributed to the acquisition of blaKPC-2 in C. europaeus and M. morganii.

Phenolic compounds promote the horizontal transfer of antibiotic resistance genes in activated sludge

Phenolic compounds are common organic pollutants in wastewater. During the wastewater treatment process, these compounds may influence the microbial community structure and functions. However, the impact of the phenolic compounds in the wastewater treatment plants on the horizontal transfer of antibiotic resistance genes (ARGs) has not been well assessed. In this study, we investigated the horizontal transfer of ARGs under the stress of phenolic compounds. The results showed that in pure culture bacteria system, p-nitrophenol (PNP), p-aminophenol (PAP) and phenol (PhOH) (10-100 mg/L) can significantly increase the horizontal transfer frequency of ARGs by 2.2-4.6, 3.6-9.4 and 1.9-9.0 fold, respectively. And, the RP4 plasmid transfer from Escherichia coli HB101 (E. coli HB101) to the bacteria in activated sludge increased obviously under the stress of phenolic compounds. Further investigation revealed that the PNP and PhOH at the concentration of 10-100 mg/L increased the production of reactive oxygen species and the permeability of cell membrane in the donor and recipient, which could be the causes of horizontal transfer of RP4 plasmid. In addition, it was also found that PNP, PAP and PhOH stress inhibit the expression of the global regulatory genes korB and trbA in the RP4 plasmid, and increase the expression level of the traF gene, thereby promoting the conjugative transfer of the RP4 plasmid. Taken together, these results improved our understanding of the horizontal transfer of ARGs under the stress of phenolic compounds and provided basic information for management of the systems that treat wastewater containing phenolic compounds.

Genomic adaptations of Campylobacter jejuni to long-term human colonization

BACKGROUND

Campylobacter is a genus of bacteria that has been isolated from the gastrointestinal tract of humans and animals, and the environments they inhabit around the world. Campylobacter adapt to new environments by changes in their gene content and expression, but little is known about how they adapt to long-term human colonization. In this study, the genomes of 31 isolates from a New Zealand patient and 22 isolates from a United Kingdom patient belonging to Campylobacter jejuni sequence type 45 (ST45) were compared with 209 ST45 genomes from other sources to identify the mechanisms by which Campylobacter adapts to long-term human colonization. In addition, the New Zealand patient had their microbiota investigated using 16S rRNA metabarcoding, and their level of inflammation and immunosuppression analyzed using biochemical tests, to determine how Campylobacter adapts to a changing gastrointestinal tract.

RESULTS

There was some evidence that long-term colonization led to genome degradation, but more evidence that Campylobacter adapted through the accumulation of non-synonymous single nucleotide polymorphisms (SNPs) and frameshifts in genes involved in cell motility, signal transduction and the major outer membrane protein (MOMP). The New Zealand patient also displayed considerable variation in their microbiome, inflammation and immunosuppression over five months, and the Campylobacter collected from this patient could be divided into two subpopulations, the proportion of which correlated with the amount of gastrointestinal inflammation.

CONCLUSIONS

This study demonstrates how genomics, phylogenetics, 16S rRNA metabarcoding and biochemical markers can provide insight into how Campylobacter adapts to changing environments within human hosts. This study also demonstrates that long-term human colonization selects for changes in Campylobacter genes involved in cell motility, signal transduction and the MOMP; and that genetically distinct subpopulations of Campylobacter evolve to adapt to the changing gastrointestinal environment.

Comparative genomics of Chinese and international isolates of Escherichia albertii: population structure and evolution of virulence and antimicrobial resistance

Escherichia albertii is a recently recognized species in the genus Escherichia that causes diarrhoea. The population structure, genetic diversity and genomic features have not been fully examined. Here, 169 E. albertii isolates from different sources and regions in China were sequenced and combined with 312 publicly available genomes (from additional 14 countries) for genomic analyses. The E. albertii population was divided into two clades and eight lineages, with lineage 3 (L3), L5 and L8 more common in China. Clinical isolates were observed in all clades/lineages. Virulence genes were found to be distributed differently among lineages: subtypes of the intimin encoding gene eae and the cytolethal distending toxin gene cdtB were lineage associated, and the second type three secretion system (ETT2) island was truncated in L3 and L6. Seven new eae subtypes and one new cdtB subtype (cdtB-VI) were identified. Alarmingly, 85.9 % of the Chinese E. albertii isolates were predicted to be multidrug-resistant (MDR) with 35.9 % harbouring genes capable of conferring resistance to 10 to 14 different drug classes. The majority of the MDR isolates were of poultry source from China and belonged to four sequence types (STs) [ST4638, ST4479, ST4633 and ST4488]. Thirty-four plasmids with some carrying MDR and virulence genes, and 130 prophages were identified from 17 complete E. albertii genomes. The 130 intact prophages were clustered into five groups, with group five prophages harbouring more virulence genes. We further identified three E. albertii specific genes as markers for the identification of this species. Our findings provided fundamental insights into the population structure, virulence variation and drug resistance of E. albertii.

Mobile Colistin Resistance Genetic Determinants of Non-Typhoid Salmonella enterica Isolates from Russia

Polymyxin resistance, determined by mcr genes located on plasmid DNA, currently poses a high epidemiological threat. Non-typhoid Salmonella (NTS) are one of the key pathogens causing diarrheal diseases. Here, we report the isolation and whole genome sequencing of multidrug colistin-resistant/susceptible isolates of non-typhoid Salmonella enterica serovars carrying mcr genes. Non-typhoid strains of Salmonella enterica subsp. enterica were isolated during microbiological monitoring of the environment, food, and diarrheal disease patients between 2018 and 2020 in Russia (n = 586). mcr-1 genes were detected using a previously developed qPCR assay, and whole genome sequencing of mcr positive isolates was performed by both short-read (Illumina) and long-read (Oxford Nanopore) approaches. Three colistin-resistant isolates, including two isolates of S. Enteritidis and one isolate of S. Bovismorbificans, carried the mcr-1.1 gene located on IncX4 and IncI2 conjugative plasmids, respectively. The phenotypically colistin-susceptible isolate of S. Typhimurium carried a mcr-9 gene on plasmid IncHI2. In conclusion, we present the first three cases of mcr gene-carrying NTS isolates detected in Russia with both outbreak and sporadic epidemiological backgrounds.

Ruminant-associated Listeria monocytogenes isolates belong preferentially to dairy-associated hypervirulent clones: a longitudinal study in 19 farms

Studies have shown that ruminants constitute reservoirs of Listeria monocytogenes, but little is known about the epidemiology and genetic diversity of this pathogen within farms. Here we conducted a large-scale longitudinal study to monitor Listeria spp. in 19 dairy farms during three consecutive seasons (N = 3251 samples). L. innocua was the most prevalent species, followed by L. monocytogenes. Listeria monocytogenes was detected in 52.6% of farms and more frequently in cattle (4.1%) and sheep (4.5%) than in goat farms (0.2%). Lineage I accounted for 69% of L. monocytogenes isolates. Among animal samples, the most prevalent sublineages (SL) and clonal complexes (CC) were SL1/CC1, SL219/CC4, SL26/CC26 and SL87/CC87, whereas SL666/CC666 was most prevalent in environmental samples. Sixty-one different L. monocytogenes cgMLST types were found, 28% common to different animals and/or surfaces within the same farm and 21% previously reported elsewhere in the context of food and human surveillance. Listeria monocytogenes prevalence was not affected by farm hygiene but by season: higher prevalence was observed during winter in cattle, and during winter and spring in sheep farms. Cows in their second lactation had a higher probability of L. monocytogenes faecal shedding. This study highlights dairy farms as a reservoir for hypervirulent L. monocytogenes.