Remarkable genomic diversity among Escherichia isolates recovered from healthy chickens

Comparative genomic and phenotypic description of Escherichia ruysiae: a newly identified member of the gut microbiome of the domestic dog

Introduction

Escherichia ruysiae is a newly identified species within the Escherichia genus, yet its presence in domestic animals remains largely unexamined. This study characterizes four isolates detected for the first time in the domestic dog (Canis lupus familiaris), focusing on their phenotypic and genomic features.

Methods

We used culturomic methods to isolate four E. ruysiae isolates that were initially identified as Escherichia coli using MALDI-TOF mass spectrometry. Whole-genome sequencing confirmed that the isolates belonged to E. ruysiae, not E. coli. Phenotypic characterization included enzymatic activity assays and antimicrobial susceptibility testing. Comparative genomic analyses were performed on these four isolates, along with 14 additional E. ruysiae and representative genomes from the five other Escherichia species in order to assess genetic diversity and functional gene distribution.

Results and discussion

All strains exhibited similar enzymatic activities and resistance to clindamycin, erythromycin, and metronidazole. The pangenome analysis revealed that most missing gene orthologs are related to motility followed by metabolism, including synthetases, reductases, phosphatases, permeases, transferases, and epimerases, as well as structural genes like efflux pumps and transporters. Phylogroup typing using the ClermonTyping method identified two main groups within the E. ruysiae species, Clade III and IV. Typical virulence genes associated with E. coli are absent in these strains. The multiple approaches used in this study expand our understanding of the diverse aspects of the recently described species, E. ruysiae.

Impact of land-use and fecal contamination on Escherichia populations in environmental samples

Understanding the composition of complex Escherichia coli populations from the environment is necessary for identifying strategies to reduce the impacts of fecal contamination and protect public health. Metabarcoding targeting the hypervariable gene gnd was used to reveal the complex population diversity of E. coli and phenotypically indistinct Escherichia species in water, soil, sediment, aquatic biofilm, and fecal samples from native forest and pastoral sites. The resulting amplicons were cross-referenced against a database containing over 700 different partial gnd sequences from E. coli/non-E. coli Escherichia species. Alpha and beta measures of diversity of Escherichia populations were lowest in feces, soil and sediment compared to water and aquatic biofilm samples. Escherichia populations recovered from extensive freshwater catchments dominated by sheep, beef and dairy farming were extremely diverse but well-separated from a wetland dairy site. In contrast, Escherichia populations from the low-impact native forest site with fewer fecal sources were less diverse. Metabarcoding distinguished E. coli populations important to fecal contamination monitoring from non-E. coli Escherichia environmental populations. These data represent in-depth analysis and geographic stability of Escherichia populations from environmental samples with extensive heterogeneity, and reveal links with diverse fecal sources, land-use and the overall burden of fecal contamination at sample sites.

Host species shapes genotype, antimicrobial resistance, and virulence profiles of enterotoxigenic Escherichia coli (ETEC) from livestock in the United States

Enterotoxigenic Escherichia coli (ETEC) are significant pathogen in both cattle and pigs, causing diarrhea in these animals and leading to economic losses in the livestock industry. Understanding the dissimilarity in genotype, antimicrobial resistance (AMR), and virulence between bovine and swine ETEC is crucial for development of targeted preventive and therapeutic approaches for livestock. However, a comprehensive study on this area remains lacking. Here, we performed whole-genome sequencing-based analyses of bovine (n = 554) and swine (n = 623) ETEC collected in the United States over a 53-year period. We identified distinct ETEC genotypes (fimH type, O antigen, H antigen, sequence type) in cattle and pigs. Furthermore, specific AMR and virulence profiles were associated with bovine and swine ETEC. Compared to swine ETEC, bovine ETEC were less diverse in genotypes and had a significantly (P < 0.001) lower number of AMR genes per isolate but higher co-occurrence of Shiga toxin and enterotoxin genes. Our results provide an overview of the key genomic differences between bovine and swine ETEC in the United States, which might be attributed to host adaptation and antibiotic usage practice. Ongoing surveillance and research are essential to monitor the genetic diversity and AMR patterns of ETEC in different host species.

IMPORTANCE

Enterotoxigenic Escherichia coli (ETEC)-associated diarrhea represent one of the most economically important diseases in the livestock industry. By analyzing over a thousand livestock-derived ETEC samples in the United States, our study unveiled a clear distinction in ETEC's genetic traits (i.e., genotypes, antimicrobial resistance [AMR], and virulence profiles) that might be tied to the different use of antibiotics in cattle and pigs, and the bacteria's adaptation to their specific animal hosts. This understanding is crucial for tailoring preventive and therapeutic strategies. It also highlights the significance of ongoing surveillance and research into the evolution of bacterial pathogens like ETEC in livestock by using advanced techniques such as whole-genome sequencing.

Pathogenic profile and antimicrobial resistance of Escherichia coli, Escherichia marmotae and Escherichia ruysiae detected from hunted wild boars in Sardinia (Italy)

The objective of this study was to evaluate the occurrence of E. coli in hunted wild boars in Sardinia (Italy) and to further characterize the isolates with Whole Genome Sequencing to assess the genetic relatedness and the presence of virulence and antimicrobial resistance (AMR) genes. Samples were taken from 66 wild boars between 2020 and 2022 slaughtered in five hunting houses. A total of 181 samples were tested, including 66 samples from mesenteric lymph nodes, 66 samples from colon content and 49 samples from carcass surface. Isolates referable to Escherichia species were detected in all of the wild boars sampled. On a selection of 61 isolates, sequencing was conducted and antimicrobial susceptibility was tested. Among these, three isolates were confirmed to be two Escherichia marmotae (cryptic clade V) and one Escherichia ruysiae (cryptic clade III). E. coli pathotypes identified were UPEC (13 %), ExPEC-UPEC (5.6 %) and ETEC (3.7 %). Moreover, 3/6 E. marmotae isolates had typical ExPEC genes. Genetic similarity was observed in isolates collected from animals slaughtered in the same hunting house; this suggests epidemiological links deriving from the presence of animals infected with closely related strains or the result of cross-contamination. Antimicrobial resistance genes were detected in three non-pathogenic E. coli isolates: one isolate had sul2, tet(B), aph(6)-ld and aph(3″)-lb resistance genes and two had the fosA7 gene. This study confirmed that wild boars can act as reservoirs and spreaders of pathogenic Escherichia species and it provides information for future comparative genomic analysis in wildlife. Although isolates showed a limited resistome, the detection of resistance in non-pathogenic isolates underlines the need to monitor antimicrobial resistance in the wild boar population. To the best of our knowledge, this is the first detection of E. mamotae and E. ruysiae isolates in wild boars in Italy and the presence of this pathogen in wildlife and livestock need to be investigated further.

A systematic review and meta-analysis on the efficacy of vaccination against colibacillosis in broiler production

Colibacillosis, a disease caused by Escherichia coli in broiler chickens has serious implications on food safety, security, and economic sustainability. Antibiotics are required for treating the disease, while vaccination and biosecurity are used for its prevention. This systematic review and meta-analysis, conducted under the COST Action CA18217-European Network for Optimization of Veterinary Antimicrobial Treatment (ENOVAT), aimed to assess the efficacy of E. coli vaccination in broiler production and provide evidence-based recommendations. A comprehensive search of bibliographic databases, including, PubMed, CAB Abstracts, Web of Science and Agricola, yielded 2,722 articles. Following a defined protocol, 39 studies were selected for data extraction. Most of the studies were experimental infection trials, with only three field studies identified, underscoring the need for more field-based research. The selected studies reported various types of vaccines, including killed (n = 5), subunit (n = 8), outer membrane vesicles/protein-based (n = 4), live/live-attenuated (n = 16), and CpG oligodeoxynucleotides (ODN) (n = 6) vaccines. The risk of bias assessment revealed that a significant proportion of studies reporting mortality (92.3%) or feed conversion ratio (94.8%) as outcomes, had "unclear" regarding bias. The meta-analysis, focused on live-attenuated and CpG ODN vaccines, demonstrated a significant trend favoring both vaccination types in reducing mortality. However, the review also highlighted the challenges in reproducing colibacillosis in experimental setups, due to considerable variation in challenge models involving different routes of infection, predisposing factors, and challenge doses. This highlights the need for standardizing the challenge model to facilitate comparisons between studies and ensure consistent evaluation of vaccine candidates. While progress has been made in the development of E. coli vaccines for broilers, further research is needed to address concerns such as limited heterologous protection, practicability for application, evaluation of efficacy in field conditions and adoption of novel approaches.

High-resolution genomic analysis to investigate the impact of the invasive brushtail possum (Trichosurus vulpecula) and other wildlife on microbial water quality assessments

Escherichia coli are routine indicators of fecal contamination in water quality assessments. Contrary to livestock and human activities, brushtail possums (Trichosurus vulpecula), common invasive marsupials in Aotearoa/New Zealand, have not been thoroughly studied as a source of fecal contamination in freshwater. To investigate their potential role, Escherichia spp. isolates (n = 420) were recovered from possum gut contents and feces and were compared to those from water, soil, sediment, and periphyton samples, and from birds and other introduced mammals collected within the Mākirikiri Reserve, Dannevirke. Isolates were characterized using E. coli-specific real-time PCR targeting the uidA gene, Sanger sequencing of a partial gnd PCR product to generate a gnd sequence type (gST), and for 101 isolates, whole genome sequencing. Escherichia populations from 106 animal and environmental sample enrichments were analyzed using gnd metabarcoding. The alpha diversity of Escherichia gSTs was significantly lower in possums and animals compared with aquatic environmental samples, and some gSTs were shared between sample types, e.g., gST535 (in 85% of samples) and gST258 (71%). Forty percent of isolates gnd-typed and 75% of reads obtained by metabarcoding had gSTs shared between possums, other animals, and the environment. Core-genome single nucleotide polymorphism (SNP) analysis showed limited variation between several animal and environmental isolates (<10 SNPs). Our data show at an unprecedented scale that Escherichia clones are shared between possums, other wildlife, water, and the wider environment. These findings support the potential role of possums as contributors to fecal contamination in Aotearoa/New Zealand freshwater. Our study deepens the current knowledge of Escherichia populations in under-sampled wildlife. It presents a successful application of high-resolution genomic methods for fecal source tracking, thereby broadening the analytical toolbox available to water quality managers. Phylogenetic analysis of isolates and profiling of Escherichia populations provided useful information on the source(s) of fecal contamination and suggest that comprehensive invasive species management strategies may assist in restoring not only ecosystem health but also water health where microbial water quality is compromised.

Genomic dissection of Escherichia marmotae provides insights into diversity and pathogenic potential

Anthropogenic activities enhance the interconnection of human, animal, and environmental habitats and drive the evolution and inter-niche transmission of bacteria. Clear identification of emerging bacteria and pathogen control is therefore a public health priority. In 2015, the novel Escherichia species Escherichia marmotae was assigned, but due to the lack of appropriate detection and typing technologies, the One Health impact of this species is still being unraveled. E. marmotae represents a missing link in the impact of Escherichia spp. Here, we report 25 E. marmotae identified by next-generation sequencing that were previously phenotypically characterized as Escherichia coli during national zoonosis monitoring of food-producing animals. Applying fastANI to 153 738 published Escherichia spp. genome assemblies, we identified further 124 E. marmotae, originally classified as E. coli. Phylogenomics of all 149 isolates reveals an undefined population structure that is independent of the ecological niche. We highlight the phenotypic, genomic, and plasmid diversity of E. marmotae and provide evidence for gene flow across the species. The latter is illustrated by the acquisition of antibiotic resistance plasmids and pathogenicity islands, such as the type III secretion system. Thus, our comprehensive genomic overview of an emerging potential opportunistic pathogen underlines the importance of improved detection and characterization.

Zoonotic potential and antimicrobial resistance of Escherichia spp. in urban crows in Japan-first detection of E. marmotae and E. ruysiae

Little is known about the prevalence of antimicrobial-resistant bacteria and pathogenic Escherichia coli in crows (carrion and jungle crows). We studied the phylogeny, virulence and antimicrobial resistance gene profiles of crow E. coli isolates to investigate their zoonotic potential and molecular epidemiology. During the winter of 2021-2022, 34 putative E. coli isolates were recovered from 27 of the 65 fresh fecal samples collected in urban areas. Three strains of the B1-O88:H8-ST446-fimH54 lineage, classified as extraintestinal pathogenic E. coli (ExPEC) and necrotoxigenic E. coli type 2, were colistin-resistant and harbored mcr-1.1-carrying IncI2 plasmids. The blaCTX-M-55 was identified in a multidrug-resistant B1-O non-typeable:H23-ST224-fimH39 strain. In phylogroup B2, two lineages of O6:H1-ST73-fimH30 and O6:H5-ST83-fimH21 were classified as ExPEC, uropathogenic E. coli, and necrotoxigenic E. coli type 1 (O6:H5-ST83-fimH21), and contained several virulence genes associated with avian pathogenic E. coli. Noteworthy is that three isolates, identified as E. coli by MALDI-TOF MS, were confirmed to be two Escherichia marmotae (cryptic clade V) and one Escherichia ruysiae (cryptic clade III) based on ANI and dDDH analyses. Our results provide the first evidence of these new species in crows. E. marmotae and E. ruysiae isolates in this study were classified as ExPEC and contained the enteroaggregative E. coli heat-stable toxin 1 gene. In addition, these two E. marmotae isolates displayed a close genetic relationship with human isolates associated with septicemia. This study provides the first insights into the prevalence and zoonotic significance of Escherichia spp. in urban crows in Japan, posing a significant risk for their transmission to humans.

Characterization and comparative analysis of the Escherichia marmotae M-12 isolate from bank vole (Myodes glareolus)

The Escherichia marmotae is a bacterium of the Enterobacterales order, which was first isolated from the Himalayan marmot (Marmota himalayana). Recently E. marmotae has been shown to cause severe infections in humans. Wild animals were suggested to be a natural reservoir of this bacterium. The present study describes the first case of E. marmotae isolation from an apparently healthy wild bank vole (Myodes glareolus). Phenotype, as well as genotype-based techniques, were applied to characterize E. marmotae M-12 isolate. E. marmotae M-12 had the capsule-positive phenotype, high adhesion to human erythrocytes and HEp-2 cells as well as a low invasion into HEp-2 cells. E. marmotae M-12 was avirulent in mice. The phylogenomic analyses of E. marmotae showed dispersed phylogenetic structure among isolates of different origins. Virulome analysis of M-12 isolate revealed the presence of the following factors: siderophores, heme uptake systems, capsule synthesis, curli and type I fimbriae, flagella proteins, OmpA porin, etc. Comparative virulome analysis among available E. marmotae genomes revealed the presence of capsule K1 genes mostly in pathogenic isolates and OmpA porin presence among all strains. We assume that the K1 capsule and OmpA porin play a key role in the virulence of E. marmotae. Pathogenesis of the latter might be similar to extraintestinal pathogenic E. coli.

Escherichia ruysiae May Serve as a Reservoir of Antibiotic Resistance Genes across Multiple Settings and Regions

Gut colonization with multidrug-resistant Enterobacterales (MDR-Ent) has reached worrisome levels worldwide. In this context, Escherichia ruysiae is a recently described species mostly found in animals. However, its spread and impact on humans is poorly understood. A stool sample from a healthy individual living in India was screened for the presence of MDR-Ent using culture-based methods. Colonies were routinely identified using matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) and phenotypically characterized by broth microdilution. Illumina and Nanopore whole-genome sequencing (WGS) platforms were implemented to generate a complete assembly. E. ruysiae genomes deposited in international databases were used for a core genome phylogenetic analysis. An extended-spectrum β-lactamase (ESBL)-producing E. coli strain (S1-IND-07-A) was isolated from the stool. WGS confirmed that S1-IND-07-A was indeed E. ruysiae, belonged to sequence type 5792 (ST5792), core genome (cg) ST89059, serotype O13/O129-:H56-like, clade IV phylogroup, and possessed five virulence factors. A copy of blaCTX-M-15 and five other antimicrobial resistance genes (ARGs) were detected in a conjugative IncB/O/K/Z plasmid. A database search identified 70 further E. ruysiae strains from 16 countries (44, 15, and 11 strains isolated from animals, the environment, and humans, respectively). The core genome phylogeny revealed five major STs: ST6467, ST8084, ST2371, ST9287, and ST5792. Three out of the seventy strains possessed important ARGs: OTP1704 (blaCTX-M-14; ST6467), SN1013-18 (blaCTX-M-15; ST5792), and CE1758 (blaCMY-2; ST7531). These strains were of human, environmental, and wild animal origin, respectively. E. ruysiae may acquire clinically important ARGs and transmit them to other species. Due to its zoonotic potential, further efforts are needed to improve routine detection and surveillance across One Health settings. IMPORTANCE Escherichia ruysiae is a recently described species of the cryptic clades III and IV of the genus Escherichia and is commonly found in animals and the environment. This work highlights the zoonotic potential of E. ruysiae, as it has been shown to colonize the human intestinal tract. Importantly, E. ruysiae may be associated with conjugative plasmids carrying clinically relevant antibiotic resistance genes. Therefore, it is important to closely monitor this species. Overall, this study highlights the need for improved identification of Escherichia species and continued surveillance of zoonotic pathogens in One Health settings.