Genome evolution and the emergence of pathogenicity in avian Escherichia coli

Two high-quality Prototheca zopfii genomes provide new insights into their evolution as obligate algal heterotrophs and their pathogenicity

The majority of the nearly 10,000 described species of green algae are photoautotrophs; however, some species have lost their ability to photosynthesize and become obligate heterotrophs that rely on parasitism for survival. Two high-quality genomes of the heterotrophic algae Prototheca zopfii Pz20 and Pz23 were obtained using short- and long-read genomic as well as transcriptomic data. The genome sizes were 31.2 Mb and 31.3 Mb, respectively, and contig N50 values of 1.99 Mb and 1.26 Mb. Although P. zopfii maintained its plastid genome, the transition to heterotrophy led to a reduction in both plastid and nuclear genome size, including the loss of photosynthesis-related genes from both the nuclear and plastid genomes and the elimination of genes encoding for carotenoid oxygenase and pheophorbide an oxygenase. The loss of genes, including basic leucine-zipper (bZIP) transcription factors, flavin adenine dinucleotide-linked oxidase, and helicase, could have played a role in the transmission of autotrophy to heterotrophs and in the processes of abiotic stress resistance and pathogenicity. A total of 66 (1.37%) and 73 (1.49%) genes were identified as potential horizontal gene transfer events in the two P. zopfii genomes, respectively. Genes for malate synthase and isocitrate lyase, which are horizontally transferred from bacteria, may play a pivotal role in carbon and nitrogen metabolism as well as the pathogenicity of Prototheca and non-photosynthetic organisms. The two high-quality P. zopfii genomes provide new insights into their evolution as obligate heterotrophs and pathogenicity.

IMPORTANCE

The genus Prototheca, characterized by its heterotrophic nature and pathogenicity, serves as an exemplary model for investigating pathobiology. The limited understanding of the protothecosis infectious disease is attributed to the lack of genomic resources. Using HiFi long-read sequencing, both nuclear and plastid genomes were generated for two strains of P. zopfii. The findings revealed a concurrent reduction in both plastid and nuclear genome size, accompanied by the loss of genes associated with photosynthesis, carotenoid oxygenase, basic leucine-zipper (bZIP) transcription factors, and others. The analysis of horizontal gene transfer revealed the presence of 1.37% and 1.49% bacterial genes, including malate synthase and isocitrate lyase, which play crucial roles in carbon and nitrogen metabolism, as well as pathogenicity and obligate heterotrophy. The two high-quality P. zopfii genomes represent valuable resources for investigating their adaptation and evolution as obligate heterotrophs, as well as for developing future prevention and treatment strategies against protothecosis.

BGCFlow: systematic pangenome workflow for the analysis of biosynthetic gene clusters across large genomic datasets

Genome mining is revolutionizing natural products discovery efforts. The rapid increase in available genomes demands comprehensive computational platforms to effectively extract biosynthetic knowledge encoded across bacterial pangenomes. Here, we present BGCFlow, a novel systematic workflow integrating analytics for large-scale genome mining of bacterial pangenomes. BGCFlow incorporates several genome analytics and mining tools grouped into five common stages of analysis such as: (i) data selection, (ii) functional annotation, (iii) phylogenetic analysis, (iv) genome mining, and (v) comparative analysis. Furthermore, BGCFlow provides easy configuration of different projects, parallel distribution, scheduled job monitoring, an interactive database to visualize tables, exploratory Jupyter Notebooks, and customized reports. Here, we demonstrate the application of BGCFlow by investigating the phylogenetic distribution of various biosynthetic gene clusters detected across 42 genomes of the Saccharopolyspora genus, known to produce industrially important secondary/specialized metabolites. The BGCFlow-guided analysis predicted more accurate dereplication of BGCs and guided the targeted comparative analysis of selected RiPPs. The scalable, interoperable, adaptable, re-entrant, and reproducible nature of the BGCFlow will provide an effective novel way to extract the biosynthetic knowledge from the ever-growing genomic datasets of biotechnologically relevant bacterial species.

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.

Impressive pan-genomic diversity of E. coli from a wild animal community near urban development reflects human impacts

Human and domesticated animal waste infiltrates global freshwater, terrestrial, and marine environments, widely disseminating fecal microbes, antibiotics, and other chemical pollutants. Emerging evidence suggests that guts of wild animals are being invaded by our microbes, including Escherichia coli, which face anthropogenic selective pressures to gain antimicrobial resistance (AMR) and increase virulence. However, wild animal sources remain starkly under-represented among genomic sequence repositories. We sequenced whole genomes of 145 E. coli isolates from 55 wild and 13 domestic animal fecal samples, averaging 2 (ranging 1-7) isolates per sample, on a preserve imbedded in a human-dominated landscape in California Bay Area, USA, to assess AMR, virulence, and pan-genomic diversity. With single nucleotide polymorphism analyses we predict potential transmission routes. We illustrate the usefulness of E. coli to aid our understanding of and ability to surveil the emergence of zoonotic pathogens created by the mixing of human and wild bacteria in the environment.

Genome-wide association study reveals serovar-associated genetic loci in Riemerella anatipestifer

BACKGROUND

The disease caused by Riemerella anatipestifer (R. anatipestifer, RA) results in large economic losses to the global duck industry every year. Serovar-related genomic variation, such as the O-antigen and capsular polysaccharide (CPS) gene clusters, has been widely used for serotyping in many gram-negative bacteria. RA has been classified into at least 21 serovars based on slide agglutination, but the molecular basis of serotyping is unknown. In this study, we performed a pan-genome-wide association study (Pan-GWAS) to identify the genetic loci associated with RA serovars.

RESULTS

The results revealed a significant association between the putative CPS synthesis gene locus and the serological phenotype. Further characterization of the CPS gene clusters in 11 representative serovar strains indicated that they were highly diverse and serovar-specific. The CPS gene cluster contained the key genes wzx and wzy, which are involved in the Wzx/Wzy-dependent pathway of CPS synthesis. Similar CPS loci have been found in some other species within the family Weeksellaceae. We have also shown that deletion of the wzy gene in RA results in capsular defects and cross-agglutination.

CONCLUSIONS

This study indicates that the CPS synthesis gene cluster of R. anatipestifer is a serotype-specific genetic locus. Importantly, our finding provides a new perspective for the systematic analysis of the genetic basis of the R anatipestifer serovars and a potential target for establishing a complete molecular serotyping scheme.

Characteristics of antimicrobial resistance in Escherichia coli isolated from retail meat products in North Carolina

BACKGROUND

Escherichia coli is commonly used as an indicator for antimicrobial resistance (AMR) in food, animal, environment, and human surveillance systems. Our study aimed to characterize AMR in E. coli isolated from retail meat purchased from grocery stores in North Carolina, USA as part of the National Antimicrobial Resistance Monitoring System (NARMS).

MATERIALS AND METHODS

Retail chicken (breast, n = 96; giblets, n = 24), turkey (n = 96), and pork (n = 96) products were purchased monthly from different counties in North Carolina during 2022. Label claims on packages regarding antibiotic use were recorded at collection. E. coli was isolated from meat samples using culture-based methods and isolates were characterized for antimicrobial resistance using whole genome sequencing. Multi-locus sequence typing, phylogroups, and a single nucleotide polymorphism (SNP)-based maximum-likelihood phylogenic tree was generated. Data were analyzed statistically to determine differences between antibiotic use claims and meat type.

RESULTS

Of 312 retail meat samples, 138 (44.2%) were positive for E. coli, with turkey (78/138; 56.5%) demonstrating the highest prevalence. Prevalence was lower in chicken (41/138; 29.7%) and pork (19/138;13.8%). Quality sequence data was available from 84.8% (117/138) of the E. coli isolates, which included 72 (61.5%) from turkey, 27 (23.1%) from chicken breast, and 18 (15.4%) from pork. Genes associated with AMR were detected in 77.8% (91/117) of the isolates and 35.9% (42/117) were defined as multidrug resistant (MDR: being resistant to ≥3 distinct classes of antimicrobials). Commonly observed AMR genes included tetB (35%), tetA (24.8%), aph(3'')-lb (24.8%), and blaTEM-1 (20.5%), the majority of which originated from turkey isolates. Antibiotics use claims had no statistical effect on MDR E. coli isolates from the different meat types (X2 = 2.21, p = 0.33). MDR was observed in isolates from meat products with labels indicating "no claims" (n = 29; 69%), "no antibiotics ever" (n = 9; 21.4%), and "organic" (n = 4; 9.5%). Thirty-four different replicon types were observed. AMR genes were carried on plasmids in 17 E. coli isolates, of which 15 (88.2%) were from turkey and two (11.8%) from chicken. Known sequence types (STs) were described for 81 E. coli isolates, with ST117 (8.5%), ST297 (5.1%), and ST58 (3.4%) being the most prevalent across retail meat types. The most prevalent phylogroups were B1 (29.1%) and A (28.2%). Five clonal patterns were detected among isolates.

CONCLUSIONS

E. coli prevalence and the presence of AMR and MDR were highest in turkey retail meat. The lack of an association between MDR E. coli in retail meat and antibiotic use claim, including those with no indication of antimicrobial use, suggests that additional research is required to understand the origin of resistance. The presence of ST117, an emerging human pathogen, warrants further surveillance. The isolates were distinctly diverse suggesting an instability in population dynamics.

Disease-Associated Streptococcus pneumoniae Genetic Variation

Streptococcus pneumoniae is an opportunistic pathogen that causes substantial illness and death among children worldwide. The genetic backgrounds of pneumococci that cause infection versus asymptomatic carriage vary substantially. To determine the evolutionary mechanisms of opportunistic pathogenicity, we conducted a genomic surveillance study in China. We collected 783 S. pneumoniae isolates from infected and asymptomatic children. By using a 2-stage genomewide association study process, we compared genomic differences between infection and carriage isolates to address genomic variation associated with pathogenicity. We identified 8 consensus k-mers associated with adherence, antimicrobial resistance, and immune modulation, which were unevenly distributed in the infection isolates. Classification accuracy of the best k-mer predictor for S. pneumoniae infection was good, giving a simple target for predicting pathogenic isolates. Our findings suggest that S. pneumoniae pathogenicity is complex and multifactorial, and we provide genetic evidence for precise targeted interventions.

Step-by-Step Bacterial Genome Comparison

Thanks to advancements in genome sequencing and bioinformatics, thousands of bacterial genome sequences are available in public databases. This presents an opportunity to study bacterial diversity in unprecedented detail. This chapter describes a complete bioinformatics workflow for comparative genomics of bacterial genomes, including genome annotation, pangenome reconstruction and visualization, phylogenetic analysis, and identification of sequences of interest such as antimicrobial-resistance genes, virulence factors, and phage sequences. The workflow uses state-of-the-art, open-source tools. The workflow is presented by means of a comparative analysis of Salmonella enterica serovar Typhimurium genomes. The workflow is based on Linux commands and scripts, and result visualization relies on the R environment. The chapter provides a step-by-step protocol that researchers with basic expertise in bioinformatics can easily follow to conduct investigations on their own genome datasets.

Association Between Escherichia coli Load in the Gut and Body Weight Gain in Broiler Chickens: A Systematic Review and Meta-Analysis

Escherichia coli is one of the first commensal bacteria to colonize the chicken gut, where it predominates at an early stage of broiler chick life. Escherichia coli can potentially cause colibacillosis in chickens, spreading to extraintestinal systemic organs, which results in high economic losses in poultry industry, as well as a potential risk to public health. Many studies conducted to investigate the effectiveness of natural products as alternatives to antibiotics and to enhance the production performance in broiler chickens have assessed E. coli load in the chicken gut, but it is still unknown how the E. coli count is linked to broiler growth performance. A systematic search of published research articles, including key terms of interest such as broiler chickens, growth performance, and E. coli count, was conducted using two main databases (PubMed and the Web of Science). A random effects metaregression model was built to evaluate the association between E. coli count and weight gain in untreated groups of broilers (negative controls) from eligible studies. Of 2108 articles in the initial screening, 60 were included in the final meta-analysis. After data extraction, records from the ileum and cecum at 21, 35, and 42 days of age were considered for the meta-analysis. The meta-analysis showed that the average E. coli count in both the ileum and cecum at 21 days of age was positively associated with the average weight gain in the studied broiler chickens, while no statistically significant associations were found at 35 and 42 days of age. In conclusion, the positive association between E. coli load and body weight gain in young broiler chickens may be attributed to the relative dominance of E. coli in the gut of this age group when the microbial population is less diverse. The dynamic association between the production performance and the load of E. coli that has dubious pathogenic potential suggests the importance of careful assessment of commensal E. coli to develop strategies to enhance production, particularly in young broiler chickens.

Identification of plasmids in avian-associated Escherichia coli using nanopore and illumina sequencing

BACKGROUND

Avian pathogenic Escherichia coli (APEC) are the causative agents of colibacillosis in chickens, a disease which has significant economic impact on the poultry industry. Large plasmids detected in APEC are known to contribute to strain diversity for pathogenicity and antimicrobial resistance, but there could be other plasmids that are missed in standard analysis. In this study, we determined the impact of sequencing and assembly factors for the detection of plasmids in an E. coli whole genome sequencing project.

RESULTS

Hybrid assembly (Illumina and Nanopore) combined with plasmid DNA extractions allowed for detection of the greatest number of plasmids in E. coli, as detected by MOB-suite software. In total, 79 plasmids were identified in 19 E. coli isolates. Hybrid assemblies were robust and consistent in quality regardless of sequencing kit used or if long reads were filtered or not. In contrast, long read only assemblies were more variable and influenced by sequencing and assembly parameters. Plasmid DNA extractions allowed for the detection of physically smaller plasmids, but when averaged over 19 isolates did not significantly change the overall number of plasmids detected.

CONCLUSIONS

Hybrid assembly can be reliably used to detect plasmids in E. coli, especially if researchers are focused on large plasmids containing antimicrobial resistance genes and virulence factors. If the goal is comprehensive detection of all plasmids, particularly if smaller sized vectors are desired for biotechnology applications, the addition of plasmid DNA extractions to hybrid assemblies is prudent. Long read sequencing is sufficient to detect many plasmids in E. coli, however, it is more prone to errors when expanded to analyze a large number of isolates.

Plasmid-mediated colistin resistance from fresh meat and slaughtered animals in the Czech Republic: nation-wide surveillance 2020-2021

The aim of this study was to determine the occurrence of plasmid-mediated colistin resistance in domestic and imported meat and slaughter animals in the Czech Republic during 2020-2021 by using selective cultivation and direct PCR testing. A total of 111 colistin-resistant Escherichia coli isolates with mcr-1 gene were obtained from 65 (9.9%, n = 659) samples and subjected to whole-genome sequencing. Isolates with mcr were frequently found in fresh meat from domestic production (14.2%) as well as from import (28.8%). The mcr-1-positive E. coli isolates predominantly originated from meat samples (16.6%), mainly poultry (27.1%), and only minor part of the isolates came from the cecum (1.7%). In contrast to selective cultivation, 205 (31.1%) samples of whole-community DNA were positive for at least one mcr variant, and other genes besides mcr-1 were detected. Analysis of whole-genome data of sequenced E. coli isolates revealed diverse sequence types (STs) including pathogenic lineages and dominance of ST1011 (15.6%) and ST162 (12.8%). Most isolates showed multidrug-resistant profile, and 9% of isolates produced clinically important beta-lactamases. The mcr-1 gene was predominantly located on one of three conjugative plasmids of IncX4 (83.5%), IncI2 (7.3%), and IncHI2 (7.3%) groups. Seventy-two percent isolates of several STs carried ColV plasmids. The study revealed high prevalence of mcr genes in fresh meat of slaughter animals. Our results confirmed previous assumptions that the livestock, especially poultry production, is an important source of colistin-resistant E. coli with the potential of transfer to humans via the food chain. IMPORTANCE We present the first data on nation-wide surveillance of plasmid-mediated colistin resistance in the Czech Republic. High occurrence of plasmid-mediated colistin resistance was found in meat samples, especially in poultry from both domestic production and import, while the presence of mcr genes was lower in the gut of slaughter animals. In contrast to culture-based approach, testing of whole-community DNA showed higher prevalence of mcr and presence of various mcr variants. Our results support the importance of combining cultivation methods with direct culture-independent techniques and highlight the need for harmonized surveillance of plasmid-mediated colistin resistance. Our study confirmed the importance of livestock as a major reservoir of plasmid-mediated colistin resistance and pointed out the risks of poultry meat for the transmission of mcr genes toward humans. We identified several mcr-associated prevalent STs, especially ST1011, which should be monitored further as they represent zoonotic bacteria circulating between different environments.

Distribution of virulence-associated genes, antibiotic resistance and phylogenetic groups in Escherichia coli isolated from domestic and racing pigeons

Despite a lot of information about virulence and resistance of Escherichia coli (E. coli) in poultry, very limited data are currently available on its occurrence in pigeon isolates, although this poses a threat to human and animal health. Therefore, this study was conducted to explore the phylogenetic classification, antibiotic sensitivity, and virulence factors in E. coli isolated from cloacal swabs of domestic pigeons bred for meat (n = 47) and racing pigeons (n = 44). The most frequent phylogroup in racing pigeons was E (36, 82.00%), unlike domestic pigeons (B2- 19, 40.00%). The most abundant iron uptake system in both groups of bird was feoB (racing = 40, 90.90%; domestic = 44, 93.61%). The presence of ibeA (52, 57.10%) and kpsMTII (46, 50.50%) genes was detected in more than half of all strains belonging exclusively to phylogroups B2, D, E, F, clade I. Antibiotic resistance was higher in racing pigeons. All racing pigeon isolates were resistant to tetracycline and trimethoprim + sulphonamide. Resistance to ciprofloxacin was determined in three isolates (6.38%) of domestic and 33 isolates (75%) of racing pigeons. Aminoglycosides and β-lactamases resistance were also recorded. One of the important detected phenotypic mechanisms of resistance occurring in isolates from racing pigeons was AGL AAC(6´)I. Our study confirms that healthy pigeons are a reservoir of antibiotic-resistant E. coli containing an arsenal of virulence factors, thus capable of potentially causing infection. Pigeons with the option to fly to multiple places can transfer virulent and resistant bacteria. Direct contact with pigeons and their faeces and the contamination of water and food pose a threat of infection to humans and other animal species.

Comparative Genomic Analysis of Hypervirulence Carbapenem-Resistant Klebsiella pneumoniae from Inpatients with Infection and Gut Colonization, China

Background

The emergence and spread of hypervirulent carbapenem-resistant Klebsiella pneumoniae (hv-CRKP) is a potential epidemiological threat that needs to be monitored. However, the transmission and pathogenic characteristics of hv-CRKP in China remain unclear. We investigated the epidemiological characteristics of gut colonized hv-CRKP in a hospital in Guangdong Province, China.

Methods

A total of 46 gut colonized hv-CRKP isolates were collected from Sun Yat-Sen Memorial Hospital (Guangzhou, China) from August 31st to December 31st, 2021. Minimum inhibitory concentrations (MICs) were obtained for 15 antibiotics for 46 hv-CRKP isolates. BALB/C mice infection model and mucoviscosity assay was used to evaluate the virulence of the isolates. The characteristics of genome, phylogenetic relationship and the structure of the plasmid of 46 gut colonized hv-CRKP isolates were compared with pathogenic isolates from GeneBank based on whole-genome data.

Results

The hv-CRKP isolation rate of all gut colonized carbapenem-resistant Klebsiella pneumoniae was 17% (46/270), and the intestinal colonization rate of hv-CRKP was irrelevant to the sex, age, department of hospitalization, and history of antibiotic use of the host. The gut colonized hv-CRKP showed pandrug resistance and hypervirulence. The gut colonized hv-CRKP and pathogenic hv-CRKP prevalent in China were mainly ST11 hv-CRKP and had two major epidemic clades. The similarities in genomic characteristics between gut colonized hv-CRKP and pathogenic hv-CRKP were consistent. The gut colonized hv-CRKP carried an incomplete structure pK2044 virulence plasmid from hypervirulent K. pneumoniae NTUH-K2044 by analyzing the virulence plasmid structure.

Conclusion

Our results suggest that the gut colonized ST11 hv-CRKP may serve as a reservoir for the clinical pathogenic ST11 HV-CRKP. It is necessary to further strengthen the monitoring of gut colonized hv-CRKP and research the potential mechanism of infection caused by gut colonized hv-CRKP.

The genetic network underlying the evolution of pathogenicity in avian Escherichia coli

Introduction

Colibacillosis is a worldwide prevalent disease in poultry production linked to Escherichia coli strains that belong to the avian pathogenic E. coli (APEC) pathotype. While many virulence factors have been linked to APEC isolates, no single gene or set of genes has been found to be exclusively associated with the pathotype. Moreover, a comprehensive description of the biological processes linked to APEC pathogenicity is currently lacking.

Methods

In this study, we compiled a dataset of 2015 high-quality avian E. coli genomes from pathogenic and commensal isolates, based on publications from 2000 to 2021. We then conducted a genome-wide association study (GWAS) and integrated candidate gene identification with available protein-protein interaction data to decipher the genetic network underlying the biological processes connected to APEC pathogenicity.

Results

Our GWAS identified variations in gene content for 13 genes and SNPs in 3 different genes associated with APEC isolates, suggesting both gene-level and SNP-level variations contribute to APEC pathogenicity. Integrating protein-protein interaction data, we found that 15 of these genes clustered in the same genetic network, suggesting the pathogenicity of APEC might be due to the interplay of different regulated pathways. We also found novel candidate genes including an uncharacterized multi-pass membrane protein (yciC) and the outer membrane porin (ompD) as linked to APEC isolates.

Discussion

Our findings suggest that convergent pathways related to nutrient uptake from host cells and defense from host immune system play a major role in APEC pathogenicity. In addition, the dataset curated in this study represents a comprehensive historical genomic collection of avian E. coli isolates and constitutes a valuable resource for their comparative genomics investigations.

Whole Genome Sequencing of Avian Pathogenic Escherichia coli Causing Bacterial Chondronecrosis and Osteomyelitis in Australian Poultry

Bacterial chondronecrosis with osteomyelitis (BCO) impacts animal welfare and productivity in the poultry industry worldwide, yet it has an understudied pathogenesis. While Avian Pathogenic Escherichia coli (APEC) are known to be one of the main causes, there is a lack of whole genome sequence data, with only a few BCO-associated APEC (APEC_BCO) genomes available in public databases. In this study, we conducted an analysis of 205 APEC_BCO genome sequences to generate new baseline phylogenomic knowledge regarding the diversity of E. coli sequence types and the presence of virulence associated genes (VAGs). Our findings revealed the following: (i) APEC_BCO are phylogenetically and genotypically similar to APEC that cause colibacillosis (APEC_colibac), with globally disseminated APEC sequence types ST117, ST57, ST69, and ST95 being predominate; (ii) APEC_BCO are frequent carriers of ColV-like plasmids that carry a similar set of VAGs as those found in APEC_colibac. Additionally, we performed genomic comparisons, including a genome-wide association study, with a complementary collection of geotemporally-matched genomes of APEC from multiple cases of colibacillosis (APEC_colibac). Our genome-wide association study found no evidence of novel virulence loci unique to APEC_BCO. Overall, our data indicate that APEC_BCO and APEC_colibac are not distinct subpopulations of APEC. Our publication of these genomes substantially increases the available collection of APEC_BCO genomes and provides insights for the management and treatment strategies of lameness in poultry.

Deciphering risks of resistomes and pathogens in intensive laying hen production chain

Antimicrobial resistance (AMR) and pathogens derived from food animals and their associated environments have emerged as challenging threats to humans from a health perspective, but our understanding of these risks and their key prevention and control points in the current intensive breeding industry remains poor. By creating an integral composition and risk profile of the resistome and microbiome through metagenomics in feces, flies, dust, sewage, and soil along the four-stage laying hen production chain, we found that the whole production chain is a hotspot for antimicrobial resistance genes (ARGs) with 374 known subtypes and pathogens, including 157 human pathogenic bacteria (HPB). Feces and flies were identified as major risk sources for these contaminations. Also, we confirmed a twin-risk of AMR and pathogenicity prevailing throughout the chain, but with different frequencies in each stage; thus, high-risk ARGs in the young chicken stage and highly prioritized HPB in the chick stage contributed 37.33 % to the total AMR risk and 36.36 % to the pathogenic risks, respectively, thus rendering the two stages to be the key prevention points. Moreover, the prevalence of 112 binned ARG supercarriers (for example, Klebsiella pneumoniae harboring 20 ARGs) was unraveled along the production chain, especially in feces, flies, and dust, and 87 potential hosts exhibited high pathogenic risk, high-risk AMR, or both, with 262 ARGs and 816 virulence factor genes. Overall, this study provides first-hand comprehensive data on high-risk ARGs and their pathogenic hosts in the intensive laying hen production chain, and thus is fundamentally important for developing new measures to help control the global AMR crisis induced through the animal-environment-human pathway.

Survey of clinical and commensal Escherichia coli from commercial broilers and turkeys, with emphasis on high-risk clones using APECTyper

Molecular characterization of avian pathogenic Escherichia coli (APEC) is challenging due to the complex nature of its associated disease, colibacillosis, in poultry. Numerous efforts have been made toward defining APEC, and it is becoming clear that certain clonal backgrounds are predictive of an avian E. coli isolate's virulence potential. Thus, APEC can be further differentiated as high-risk APEC based upon their clonal background's virulence potential. However, less clear is the degree of overlap between clinical isolates of differing bird type, and between clinical and gastrointestinal isolates. This study aimed to determine genomic similarities and differences between such populations, comparing commercial broiler vs. turkey isolates, and clinical vs. gastrointestinal isolates. Differences were observed in Clermont phylogenetic groups between isolate populations, with B2 as the dominant group in turkey clinical isolates and G as the dominant group in broiler clinical isolates. Nearly all clinical isolates were classified as APEC using a traditional gene-based typing scheme, whereas 53.4% and 44.1% of broiler and turkey gastrointestinal isolates were classified as APEC, respectively. High-risk APEC were identified among 31.0% and 46.9% of broiler and turkey clinical isolates, compared with 5.7% and 2.9% of broiler and turkey gastrointestinal isolates. As found in previous studies, no specific known virulence or fitness gene sets were identified which universally differentiate between clinical and gastrointestinal isolates. This study further demonstrates the utility of a hybrid APEC typing approach, considering both plasmid content and clonal background, for the identification of dominant and highly virulent APEC clones in poultry production.

Genome-wide association reveals host-specific genomic traits in Escherichia coli

BACKGROUND

Escherichia coli is an opportunistic pathogen which colonizes various host species. However, to what extent genetic lineages of E. coli are adapted or restricted to specific hosts and the genomic determinants of such adaptation or restriction is poorly understood.

RESULTS

We randomly sampled E. coli isolates from four countries (Germany, UK, Spain, and Vietnam), obtained from five host species (human, pig, cattle, chicken, and wild boar) over 16 years, from both healthy and diseased hosts, to construct a collection of 1198 whole-genome sequenced E. coli isolates. We identified associations between specific E. coli lineages and the host from which they were isolated. A genome-wide association study (GWAS) identified several E. coli genes that were associated with human, cattle, or chicken hosts, whereas no genes associated with the pig host could be found. In silico characterization of nine contiguous genes (collectively designated as nan-9) associated with the human host indicated that these genes are involved in the metabolism of sialic acids (Sia). In contrast, the previously described sialic acid regulon known as sialoregulon (i.e. nanRATEK-yhcH, nanXY, and nanCMS) was not associated with any host species. In vitro growth experiments with a Δnan-9 E. coli mutant strain, using the sialic acids 5-N-acetylneuraminic acid (Neu5Ac) and N-glycolylneuraminic acid (Neu5Gc) as sole carbon source, showed impaired growth behaviour compared to the wild-type.

CONCLUSIONS

This study provides an extensive analysis of genetic determinants which may contribute to host specificity in E. coli. Our findings should inform risk analysis and epidemiological monitoring of (antimicrobial resistant) E. coli.

Migration Rates on Swim Plates Vary between Escherichia coli Soil Isolates: Differences Are Associated with Variants in Metabolic Genes

This study investigates migration phenotypes of 265 Escherichia coli soil isolates from the Buffalo River basin in Minnesota, USA. Migration rates on semisolid tryptone swim plates ranged from nonmotile to 190% of the migration rate of a highly motile E. coli K-12 strain. The nonmotile isolate, LGE0550, had mutations in flagellar and chemotaxis genes, including two IS3 elements in the flagellin-encoding gene fliC. A genome-wide association study (GWAS), associating the migration rates with genetic variants in specific genes, yielded two metabolic variants (rygD-serA and metR-metE) with previous implications in chemotaxis. As a novel way of confirming GWAS results, we used minimal medium swim plates to confirm the associations. Other variants in metabolic genes and genes that are associated with biofilm were positively or negatively associated with migration rates. A determination of growth phenotypes on Biolog EcoPlates yielded differential growth for the 10 tested isolates on d-malic acid, putrescine, and d-xylose, all of which are important in the soil environment. IMPORTANCE E. coli is a Gram-negative, facultative anaerobic bacterium whose life cycle includes extra host environments in addition to human, animal, and plant hosts. The bacterium has the genomic capability of being motile. In this context, the significance of this study is severalfold: (i) the great diversity of migration phenotypes that we observed within our isolate collection supports previous (G. NandaKafle, A. A. Christie, S. Vilain, and V. S. Brözel, Front Microbiol 9:762, 2018, https://doi.org/10.3389/fmicb.2018.00762; Y. Somorin, F. Abram, F. Brennan, and C. O'Byrne, Appl Environ Microbiol 82:4628-4640, 2016, https://doi.org/10.1128/AEM.01175-16) ideas of soil promoting phenotypic heterogeneity, (ii) such heterogeneity may facilitate bacterial growth in the many different soil niches, and (iii) such heterogeneity may enable the bacteria to interact with human, animal, and plant hosts.

Identification and characterization of microRNAs, especially gga-miR-181b-5p, in chicken macrophages associated with avian pathogenic E. coli infection

AbstractAvian pathogenic E. coli (APEC) is a common pathogen in the poultry industry, which can cause substantial economic losses. Recently, emerging evidence showed that the miRNAs were involved in various viral and bacterial infection. To elucidate the role of miRNAs in chicken macrophages in response to APEC infection, we attempted to investigate the miRNAs expression pattern upon APEC infection via miRNA-seq, and to identify the molecular mechanism of the important miRNAs by using RT-qPCR, Western blotting, dual-luciferase reporter assay, and CCK-8. Results showed that a total of 80 differentially expressed (DE) miRNAs were identified in the comparison of APEC vs. wild type group, which corresponded to 724 target genes. Moreover, the target genes of the identified DE miRNAs were mainly significantly enriched in MAPK signaling pathway, Autophagy-animal, mTOR signaling pathway, ErbB signaling pathway, Wnt signaling pathway, TGF-beta signaling pathway. Remarkably, gga-miR-181b-5p is capable to participate in host immune and inflammatory response against APEC infection via targeting of TGFBR1 to modulate the activation of TGF-beta signaling pathway. Collectively, this study provides a perspective of miRNA expression pattern in chicken macrophages upon APEC infection. These findings provide the insight into miRNAs against APEC infection and gga-miR-181b-5p might be a potential target for treating APEC infection.

A molecular epidemiological study on Escherichia coli in young chicks with colibacillosis identified two possible outbreaks across farms

Avian pathogenic Escherichia coli (APEC) is the cause of colibacillosis outbreaks in young poultry chicks, resulting in acute to peracute death. The high morbidity and mortality caused by colibacillosis results in poor animal welfare, reduced sustainability and economical loss worldwide. To advance the understanding of the molecular epidemiology, genomic relatedness and virulence traits of APEC, we performed systematic sampling from 45 confirmed colibacillosis broiler flocks with high first week mortality (FWM) during 2018-2021. From these flocks, 219 APEC isolates were whole genome sequenced (WGS) and bioinformatic analyses were performed. The bioinformatic analyses included sequence typing (ST), serotyping, detection of virulence-associated genes (VAGs) and phylogenetic analysis. Our results showed a high prevalence of ST23, ST429 and ST95 among APEC isolates from Norwegian broiler flocks, and identified ST23, ST429, ST117 and ST371 to cause disease more often alone, compared to ST95, ST69 and ST10. Phylogenetic analyses, together with associated metadata, identified two distinct outbreaks of colibacillosis across farms caused by ST429 and ST23 and gave insight into expected SNP distances within and between flocks identified with the same ST. Further, our results highlighted the need for combining two typing methods, such as serotyping and sequence typing, to better discriminate strains of APEC. Ultimately, systematic sampling of APEC from multiple birds in a flock, together with WGS as a diagnostic tool is important to identify the disease-causing APEC within a flock and to detect outbreaks of colibacillosis across farms.

Comparative Genomics of the Genus Pseudomonas Reveals Host- and Environment-Specific Evolution

Each Earth ecosystem has unique microbial communities. Pseudomonas bacteria have evolved to occupy a plethora of different ecological niches, including living hosts, such as animals and plants. Many genes necessary for the Pseudomonas-niche interaction and their encoded functions remain unknown. Here, we describe a comparative genomic study of 3,274 genomes with 19,056,667 protein-coding sequences from Pseudomonas strains isolated from diverse environments. We detected functional divergence of Pseudomonas that depends on the niche. Each group of strains from a certain environment harbored a distinctive set of metabolic pathways or functions. The horizontal transfer of genes, which mainly proceeded between closely related taxa, was dependent on the isolation source. Finally, we detected thousands of undescribed proteins and functions associated with each Pseudomonas lifestyle. This research represents an effort to reveal the mechanisms underlying the ecology, pathogenicity, and evolution of Pseudomonas, and it will enable clinical, ecological, and biotechnological advances. IMPORTANCE Microbes play important roles in the health of living beings and in the environment. The knowledge of these functions may be useful for the development of new clinical and biotechnological applications and the restoration and preservation of natural ecosystems. However, most mechanisms implicated in the interaction of microbes with the environment remain poorly understood; thus, this field of research is very important. Here, we try to understand the mechanisms that facilitate the differential adaptation of Pseudomonas-a large and ubiquitous bacterial genus-to the environment. We analyzed more than 3,000 Pseudomonas genomes and searched for genetic patterns that can be related with their coevolution with different hosts (animals, plants, or fungi) and environments. Our results revealed that thousands of genes and genetic features are associated with each niche. Our data may be useful to develop new technical and theoretical advances in the fields of ecology, health, and industry.

Genetic characterization of third- or fourth-generation cephalosporin-resistant avian pathogenic Escherichia coli isolated from broilers

The third- or fourth-generation cephalosporins (3GC or 4 GC) are classified as "critically important antimicrobials for human medicine" by WHO, but resistance to these drugs is increasing rapidly in avian pathogenic E. coli (APEC). This study investigated the distribution and genetic characteristics of 3GC- or 4 GC-resistant APEC isolates from five major integrated broiler operations in Korea. The prevalence of 3GC- or 4GC-resistant APEC isolates in 1-week-old broilers was the highest in farms of operation C (53.3%); however, the highest prevalence of these isolates in 4-week-old broilers was the highest on the farms of operation A (60.0%), followed by operations E (50.0%) and C (35.7%). All 49 3GC- or 4GC-resistant APEC isolates had at least one β-lactamase-encoding gene. The most common β-lactamase-encoding genes was extended-spectrum β-lactamase gene, bla _CTX-M-15, detected in 24 isolates (49.0%), followed by bla _TEM-1 (32.7%). Sixteen isolates (32.7%) harbored class 1 integrons, and four isolates (8.2%) showed different gene cassette-arrangements. However, only 1 of 26 isolates harboring class 2 integrons carried a gene cassette. Furthermore, both CRISPR 1 and 2 arrays were detected in most isolates (36 isolates; 73.5%), followed by CRISPR 2 (18.4%) and CRISPR 1 (4.1%). Interestingly, CRISPR 2 was significantly more prevalent in multidrug resistant (MDR)-APEC isolates than in non-MDR APEC isolates, whereas CRISPR 3 and 4 were significantly more prevalent in non-MDR APEC isolates (each 11.1%; p < 0.05). None of the protospacers of CRISPR arrays were directly associated with antimicrobial resistance. Our findings indicate that the distribution and characteristics of 3GC or 4GC-resistant APEC isolates differed among the integrated broiler operations; moreover, improved management protocols are needed to control the horizontal transmission of 3GC or 4GC-resistant APEC isolates.

Genomics and pathotypes of the many faces of Escherichia coli

Escherichia coli is the most researched microbial organism in the world. Its varied impact on human health, consisting of commensalism, gastrointestinal disease, or extraintestinal pathologies, has generated a separation of the species into at least eleven pathotypes (also known as pathovars). These are broadly split into two groups, intestinal pathogenic E. coli (InPEC) and extraintestinal pathogenic E. coli (ExPEC). However, components of E. coli's infinite open accessory genome are horizontally transferred with substantial frequency, creating pathogenic hybrid strains that defy a clear pathotype designation. Here, we take a birds-eye view of the E. coli species, characterizing it from historical, clinical, and genetic perspectives. We examine the wide spectrum of human disease caused by E. coli, the genome content of the bacterium, and its propensity to acquire, exchange, and maintain antibiotic resistance genes and virulence traits. Our portrayal of the species also discusses elements that have shaped its overall population structure and summarizes the current state of vaccine development targeted at the most frequent E. coli pathovars. In our conclusions, we advocate streamlining efforts for clinical reporting of ExPEC, and emphasize the pathogenic potential that exists throughout the entire species.

Genome-wide association study of Klebsiella pneumoniae identifies variations linked to carbapenems resistance

Klebsiella pneumoniae (KP) is one of the microorganisms that can acquire carbapenem-resistance (CR), and few antimicrobial therapy options exist for infections caused by Carbapenem-Resistant KP (CRKP). In recent years, with the increase of carbapenem resistance rates, treating CRKP has become a serious public health threat in clinical practice. We have collected 2,035 clinical KP isolates from a tertiary hospital in China. Whole genome sequencing data coupled with their binary antimicrobial susceptibility testing data were obtained to conduct the genome-wide association study using a bayesian-based method, including single nucleotide polymorphisms (SNPs) and genes. We identified 28 and 37 potential maker genes associated with imipenem and meropenem resistance, respectively. Among which 19 of them were selected in both drugs by genome-wide association study (GWAS), 11 genes among them were simultaneously validated in independent datasets. These genes were likely related to biofilm formation, efflux pump, and DNA repairing. Moreover, we identified 13 significant CR related SNPs in imipenem or meropenem, with one SNP located in the non-coding region and validated in the independent datasets. Our study indicates complex mechanisms of carbapenems resistance and further investigation of CRKP-related factors are warranted to better understand their contributions to carbapenems resistance. These identified biomarkers may provide targets for future drug interventions or treatments.

Autogenous Escherichia coli Vaccine Application as an Innovative Antimicrobial Therapy in Poultry Farming-A Case Report

Escherichia coli (E. coli) is one of the most common bacterial causes of infection in poultry farming. Whether the infection is localized or systemic, a primary or secondary disease, it is most frequently treated through the application of wide-spectrum antimicrobials. Excessive use of antimicrobials in agriculture is significantly contributing to the worldwide rise of antimicrobial resistance, but is also very expensive and often ineffective in the long term. Here, we present a case where a colibacillosis outbreak on a family farm of laying hens was treated using an autogenous vaccine. The birds had septicemia, cellulitis, and severe skin wounds. They were not vaccinated against E. coli, and did not receive any antimicrobials previously. E. coli strains were isolated from the daily mortalities on the farm and used for preparation of the vaccine. Each bird was given an intramuscular injection of the autogenous vaccine. The immunogenicity of the vaccine was tested by the determination of specific antibody levels in the sera of the birds using the in-house ELISA. Shortly after vaccination, the morbidity and mortality rates significantly decreased, and egg production was improved. The application of the autogenous vaccine served as a curative and preventive measure, and has proven to be a very efficient method of antimicrobial therapy.

Distribution of papG alleles among uropathogenic Escherichia coli from reproductive age women

BACKGROUND

Extraintestinal Escherichia coli (E. coli) causing urinary tract infections (UTIs), and often referred to as uropathogenic E. coli (UPEC), are a major contributor to the morbidity of UTIs and associated healthcare costs. UPEC possess several virulence factors (VFs) for infecting and injuring the host. We studied the papG allele distribution, and its association with other VF genes and phylogenetic groups, amongst 836 UPEC and fecal isolates from reproductive age women.

RESULTS

The papGII gene was highly prevalent amongst pyelonephritis isolates (68%), whilst the majority, albeit smaller proportion, of cystitis isolates (31%) harboured the papGIII gene. Among the pyelonephritis and cystitis isolates, papG positive isolates on average had higher VF gene scores, and were more likely to belong to phylogenetic group B2, than their negative counterparts. This was mostly due to the contribution of papGII isolates, which on average contained more VF genes than their papGIII counterparts, irrespective of the uro-clinical syndrome. However, the papGII isolates from the pyelonephritis cohort had higher VF gene scores than the cystitis ones, suggesting presence of possible papGII clones with differing inferred virulence potential. Furthermore, papGII isolates were more likely to possess an intact pap gene operon than their papGIII counterparts. Also of note was the high proportion of isolates with the papGI allele which was not associated with other pap operon genes; and this finding has not been described before.

CONCLUSIONS

The association of the papGII gene with several VF genes compared to the papGIII gene, appears to explain the abundance of these genes in pyelonephritis and cystitis isolates, respectively.

A comparative genomics approach for identifying genetic factors in Escherichia coli isolates associated with bovine diseases

AIMS

Escherichia coli are ubiquitously present bacterial pathogens that cause septicaemia, diarrhoea and other clinical illness in farm animals. Many pathogen factors can be associated with disease conditions. Currently, studies inferring E. coli genetic factors associated with infection in bovines are limited. Hence, the present study envisaged to determine the pathogen genetic factors associated with bovine disease conditions.

METHOD AND RESULTS

The comparative genomic analysis involved genome sequence data of 135 diseased and 145 healthy bovine origin E. coli strains. Phylogroups A and C, as well as pathotypes ExPEC and EPEC, were found to have a strong connection with bovine disease strains. STEC strains, including EHEC, seem to play a less important role in bovine disease. Sequence types (STs) predominant among strains from diarrhoeal origin were ST 301 (CC 165) and ST 342. Correlation of core genome phylogeny with accessory gene based clustering, phylogroups and pathotypes indicated lineage specific virulence factors mostly associated with disease conditions.

CONCLUSIONS

Comparative genomic analysis was applied to infer genetic factors significant in bovine disease origin E. coli strains. Isolates from bovine disease origin were enriched for the phylogroups A and C, and for the pathotypes ExPEC and EPEC. However, there was minimal evidence of STEC involvement. The study also indicated predominant genetic lineages and virulence genes (pap, sfa and afa) associated with disease origin strains. SIGNIFICANCE AND IMPACT OF STUDY;: The study revealed significant pathotypes, phylgroups, serotypes and sequence types associated with bovine disease conditions. These identified genetic factors can be applied for disease diagnosis, implementing vaccine and therapeutic measures. In addition, E. coli isolates from the bovine species revealed a complex pattern of disease epidemiology.

Bacteriocin-Producing Escherichia coli Isolated from the Gastrointestinal Tract of Farm Animals: Prevalence, Molecular Characterization and Potential for Application

Due to the spread of antibiotic-resistant bacteria, new alternatives to antibiotics and ways to prevent infections are being sought. Bacteriocin-producing bacteria are therefore attracting attention due to their probiotic potential as a safe alternative to antimicrobial drugs. The aim of this work was to determine the prevalence of bacteriocin-encoded genes among Escherichia coli strains from healthy farm animals and to characterize the presence of virulence-associated genes, the possibility of prophage induction, and hemolytic and bacterial antagonistic activity of the bacteriocin-producing E. coli in order to reveal their potential for application. It was found that 17 of 72 E. coli strains (23.6%) produced bacteriocins. Among them, 18 out of 30 bacteriocin genes were detected: the most prevalent genes were those for microcin M (58.8%), colicin E1 (52.9%), and colicin M (35.3%). Colicin Ia (29.4%), colicin E9, colicin Ib, colicin B (23.5%), and colicin E9 (17.7%) genes were also frequent, while the prevalence of genes encoding microcins V, B17, and H47 and colicins E3, K, N, U, Y, 5, and 10 did not exceed 11.8%. At least two different bacteriocin genes were detected in all 17 bacteriocinogenic strains; the highest number of different bacteriocin genes detected in one strain was seven genes. E. coli strains with combinations of colicin E1 and E or microcin M and colicin E1 genes were more prevalent than others (17.7%). Among the 17 bacteriocin-producing E. coli strains, 5.9% were hemolytic, 47.1% contained prophages, and 58.8% carried genes encoding toxins. Cell-free supernatants of bacteriocin-producing strains were shown to inhibit the growth of pathogenic E. coli strains belonging to the APEC, STEC, and ETEC pathotypes. Thus, among the studied bacteriocin-producing E. coli isolated from the gastrointestinal tract of farm animals, three strains with high antagonistic bacterial activity and the absence of pathogenicity genes, prophages, and hemolytic activity were identified and therefore have potential for application.

Combined comparative genomics and clinical modeling reveals plasmid-encoded genes are independently associated with Klebsiella infection

Members of the Klebsiella pneumoniae species complex frequently colonize the gut and colonization is associated with subsequent infection. To identify genes associated with progression from colonization to infection, we undertook a case-control comparative genomics study. Concordant cases (N = 85), where colonizing and invasive isolates were identical strain types, were matched to asymptomatically colonizing controls (N = 160). Thirty-seven genes are associated with infection, 27 of which remain significant following adjustment for patient variables and bacterial phylogeny. Infection-associated genes are not previously characterized virulence factors, but instead a diverse group of stress resistance, regulatory and antibiotic resistance genes, despite careful adjustment for antibiotic exposure. Many genes are plasmid borne, and for some, the relationship with infection is mediated by gut dominance. Five genes were validated in a geographically-independent cohort of colonized patients. This study identifies several genes reproducibly associated with progression to infection in patients colonized by diverse Klebsiella.

Probing into the unique relationship between a soil bacterium, Pseudomonas putida AKMP7 and Arabidopsis thaliana: A case of "conditional pathogenesis"

Plant growth-promoting rhizobacteria (PGPR) are beneficial soil bacteria that colonise the rhizosphere and help plants in growth, development, and stress tolerance. While there is a significant body of research elucidating their benefits to plants, studies on the "abnormal" or "unexpected" behavior of these bacteria are almost non-existent. One such study from our laboratory has previously reported a unique situation in which a certain strain of drought and thermo-tolerant PGPR, namely, Pseudomonas putida AKMP7, becomes pathogenic towards Arabidopsis thaliana under drought conditions, but not under normal (well-watered) conditions. In this study, we have probed deeper into this phenomenon of "conditional pathogenesis". We found that, AKMP7 imparts an enhancement in plant growth under well-watered conditions, while, causing a deterioration in plant health under drought conditions. In an attempt to understand the underlying reasons for this phenomenon, we analysed the phytohormones released by Pseudomonas putida AKMP7 using LC-ESI-MS/MS technique. We identified that AKMP7 releases zeatin (a cytokinin), the auxin derivative -indole acetamide and amino acid-conjugates of auxin (indole-3-acetyl-L-alanine, indole-3-acetyl-L-phenylalanine and indole-3-acetyl-L-aspartate) in the growth medium. By treating the plants with commercially obtained forms of these phytohormones, individually or in combination with AKMP7, we identified that zeatin and auxin derivative indole acetamide can play a crucial role in the conditional pathogenesis exhibited by this bacterium on A. thaliana under drought conditions. Our work lays a foundation for further understanding the precise molecular mechanisms involved in this unique phenomenon of conditional/opportunistic pathogenesis.

Longitudinal study on background lesions in broiler breeder flocks and their progeny, and genomic characterisation of Escherichia coli

In broiler breeders, background mortality is rarely addressed, however, it represents the death of a vast number of birds, a constant productivity loss, welfare concerns and it might affect chick quality. The study aimed to unveil lesions leading to mortality in a study population perceived as healthy, combined with whole-genome sequencing (WGS) of Escherichia coli, a well-known contributor to disease problems in poultry. Broiler breeders (n = 340) originating from three distinct, putative healthy flocks and their progeny (n = 154) were subjected to a comprehensive post-mortem examination, bacteriological sampling, and sequencing of 77 E. coli isolates. Productivity data confirmed an exemplary health status of the enrolled flocks, and post-mortem examination further verified the absence of general disease problems. Among the submitted broiler breeders, exudative peritonitis (31.2%) was the most frequent lesion linked to infectious disease, whereas airsacculitis, pericarditis, perihepatitis, and salpingitis occurred in 18.5%, 3.5%, 3.8% and 17%, respectively. Yolksacculitis occurred in 15.6% of the broilers, whilst pericarditis, perihepatitis and peritonitis were diagnosed in 9.7%, 7.1% and 9.1%, respectively. WGS revealed a diverse population where ST95 dominated the population retrieved from broiler breeders, whereas ST10 was highly prevalent among broilers. Both lineages could be isolated from extraintestinal sites of birds without lesions indicative of infection. In general, the genetic diversity within flocks was comparable to the diversity between farms, and the overall occurrence of resistance markers was low. In conclusion, a comprehensive insight into lesions associated with background mortality is presented, together with a vast diversity of E. coli isolated from extraintestinal sites during a non-outbreak situation.

Comparative Genomics of Thaumarchaeota From Deep-Sea Sponges Reveal Their Niche Adaptation

Thaumarchaeota account for a large portion of microbial symbionts in deep-sea sponges and are even dominant in some cases. In this study, we investigated three new sponge-associated Thaumarchaeota from the deep West Pacific Ocean. Thaumarchaeota were found to be the most dominant phylum in this sponge by both prokaryotic 16S rRNA amplicons and metagenomic sequencing. Fifty-seven published Thaumarchaeota genomes from sponges and other habitats were included for genomic comparison. Similar to shallow sponge-associated Thaumarchaeota, those Thaumarchaeota in deep-sea sponges have extended genome sizes and lower coding density compared with their free-living lineages. Thaumarchaeota in deep-sea sponges were specifically enriched in genes related to stress adapting, symbiotic adhesion and stability, host–microbe interaction and protein transportation. The genes involved in defense mechanisms, such as the restriction-modification system, clustered regularly interspaced short palindromic repeats (CRISPR)/Cas system, and toxin-antitoxin system were commonly enriched in both shallow and deep sponge-associated Thaumarchaeota. Our study demonstrates the significant effects of both depth and symbiosis on forming genomic characteristics of Thaumarchaeota, and provides novel insights into their niche adaptation in deep-sea sponges.

Genomic and Temporal Trends in Canine ExPEC Reflect Those of Human ExPEC

Companion animals and humans are known to share extraintestinal pathogenic Escherichia coli (ExPEC), but the extent of E. coli sequence types (STs) that cause extraintestinal diseases in dogs is not well understood. Here, we generated whole-genome sequences of 377 ExPEC collected by the University of Melbourne Veterinary Hospital from dogs over an 11-year period from 2007 to 2017. Isolates were predominantly from urogenital tract infections (219, 58.1%), but isolates from gastrointestinal specimens (51, 13.5%), general infections (72, 19.1%), and soft tissue infections (34, 9%) were also represented. A diverse collection of 53 STs were identified, with 18 of these including at least five sequences. The five most prevalent STs were ST372 (69, 18.3%), ST73 (31, 8.2%), ST127 (22, 5.8%), ST80 (19, 5.0%), and ST58 (14, 3.7%). Apart from ST372, all of these are prominent human ExPEC STs. Other common ExPEC STs identified included ST12, ST131, ST95, ST141, ST963, ST1193, ST88, and ST38. Virulence gene profiles, antimicrobial resistance carriage, and trends in plasmid carriage for specific STs were generally reflective of those seen in humans. Many of the prominent STs were observed repetitively over an 11-year time span, indicating their persistence in the dogs in the community, which is most likely driven by household sharing of E. coli between humans and their pets. The case of ST372 as a dominant canine lineage observed sporadically in humans is flagged for further investigation. IMPORTANCE Pathogenic E. coli that causes extraintestinal infections (ExPEC) in humans and canines represents a significant burden in hospital and veterinary settings. Despite the obvious interrelationship between dogs and humans favoring both zoonotic and anthropozoonotic infections, whole-genome sequencing projects examining large numbers of canine-origin ExPEC are lacking. In support of anthropozoonosis, we found that most STs from canine infections are dominant human ExPEC STs (e.g., ST73, ST127, ST131) with similar genomic traits, such as plasmid carriage and virulence gene burden. In contrast, we identified ST372 as the dominant canine ST and a sporadic cause of infection in humans, supporting zoonotic transfer. Furthermore, we highlight that, as is the case in humans, STs in canine disease are consistent over time, implicating the gastrointestinal tract as the major community reservoir, which is likely augmented by exposure to human E. coli via shared diet and proximity.

To kill or to be killed: pangenome analysis of Escherichia coli strains reveals a tailocin specific for pandemic ST131

BACKGROUND

Escherichia coli (E. coli) has been one of the most studied model organisms in the history of life sciences. Initially thought just to be commensal bacteria, E. coli has shown wide phenotypic diversity including pathogenic isolates with great relevance to public health. Though pangenome analysis has been attempted several times, there is no systematic functional characterization of the E. coli subgroups according to the gene profile.

RESULTS

Systematically scanning for optimal parametrization, we have built the E. coli pangenome from 1324 complete genomes. The pangenome size is estimated to be ~25,000 gene families (GFs). Whereas the core genome diminishes as more genomes are added, the softcore genome (≥95% of strains) is stable with ~3000 GFs regardless of the total number of genomes. Apparently, the softcore genome (with a 92% or 95% generation threshold) can define the genome of a bacterial species listing the critically relevant, evolutionarily most conserved or important classes of GFs. Unsupervised clustering of common E. coli sequence types using the presence/absence GF matrix reveals distinct characteristics of E. coli phylogroups B1, B2, and E. We highlight the bi-lineage nature of B1, the variation of the secretion and of the iron acquisition systems in ST11 (E), and the incorporation of a highly conserved prophage into the genome of ST131 (B2). The tail structure of the prophage is evolutionarily related to R2-pyocin (a tailocin) from Pseudomonas aeruginosa PAO1. We hypothesize that this molecular machinery is highly likely to play an important role in protecting its own colonies; thus, contributing towards the rapid rise of pandemic E. coli ST131.

CONCLUSIONS

This study has explored the optimized pangenome development in E. coli. We provide complete GF lists and the pangenome matrix as supplementary data for further studies. We identified biological characteristics of different E. coli subtypes, specifically for phylogroups B1, B2, and E. We found an operon-like genome region coding for a tailocin specific for ST131 strains. The latter is a potential killer weapon providing pandemic E. coli ST131 with an advantage in inter-bacterial competition and, suggestively, explains their dominance as human pathogen among E. coli strains.

Genomic Profiling of Antibiotic-Resistant Escherichia coli Isolates from Surface Water of Agricultural Drainage in North-Western Mexico: Detection of the International High-Risk Lineages ST410 and ST617

Aquatic environments are recognized as one of the main reservoirs for the emergence and dissemination of high-risk lineages of multidrug-resistant (MDR) bacteria of public health concern. However, the genomic characteristics of antibiotic-resistant Escherichia coli isolates from aquatic origins remain limited. Herein, we examined the antibiotic resistance and virulence genomic profiles of three E. coli recovered from surface water in northwest Mexico. Antimicrobial susceptibility testing, whole-genome sequencing (WGS), and in-depth in silico analysis were performed. Two E. coli exhibited MDR phenotypes. WGS-based typing revealed genetic diversity, and phylogenetic analysis corroborated a notable divergent relationship among the studied E. coli. One E. coli strain, harboring enterotoxigenic and extraintestinal pathogenic-associated virulence genes, was assigned to the ST4 lineage. MDR E. coli, belonging to the international high-risk clones ST410 and ST617, carried genes and mutations conferring resistance to aminoglycosides, β-lactams, quinolones, sulfonamides, tetracyclines, and trimethoprim. This study describes, for the first time, the detection and genomic profiling of high-risk lineages of E. coli ST410 and ST617 from surface water in Mexico. Additionally, our results underscore the role of surface water as a reservoir for critical pathogenic and MDR E. coli clones and the need for the surveillance and monitoring of aquatic environments via WGS from the One Health perspective.

Population genomics of Escherichia coli in livestock-keeping households across a rapidly developing urban landscape

Quantitative evidence for the risk of zoonoses and the spread of antimicrobial resistance remains lacking. Here, as part of the UrbanZoo project, we sampled Escherichia coli from humans, livestock and peri-domestic wildlife in 99 households across Nairobi, Kenya, to investigate its distribution among host species in this rapidly developing urban landscape. We performed whole-genome sequencing of 1,338 E. coli isolates and found that the diversity and sharing patterns of E. coli were heavily structured by household and strongly shaped by host type. We also found evidence for inter-household and inter-host sharing and, importantly, between humans and animals, although this occurs much less frequently. Resistome similarity was differently distributed across host and household, consistent with being driven by shared exposure to antimicrobials. Our results indicate that a large, epidemiologically structured sampling framework combined with WGS is needed to uncover strain-sharing events among different host populations in complex environments and the major contributing pathways that could ultimately drive the emergence of zoonoses and the spread of antimicrobial resistance.

Phylogeography and phylogenomic analyses of E. coli isolates collected from humans and domesticated and wild animals across 99 households in Nairobi reveal strong intra-household, and lower but detectable inter-household and inter-host, strain-sharing patterns.

Genome wide association study of Escherichia coli bloodstream infection isolates identifies genetic determinants for the portal of entry but not fatal outcome

Escherichia coli is an important cause of bloodstream infections (BSI), which is of concern given its high mortality and increasing worldwide prevalence. Finding bacterial genetic variants that might contribute to patient death is of interest to better understand infection progression and implement diagnostic methods that specifically look for those factors. E. coli samples isolated from patients with BSI are an ideal dataset to systematically search for those variants, as long as the influence of host factors such as comorbidities are taken into account. Here we performed a genome-wide association study (GWAS) using data from 912 patients with E. coli BSI from hospitals in Paris, France. We looked for associations between bacterial genetic variants and three patient outcomes (death at 28 days, septic shock and admission to intensive care unit), as well as two portals of entry (urinary and digestive tract), using various clinical variables from each patient to account for host factors. We did not find any association between genetic variants and patient outcomes, potentially confirming the strong influence of host factors in influencing the course of BSI; we however found a strong association between the papGII operon and entrance of E. coli through the urinary tract, which demonstrates the power of bacterial GWAS when applied to actual clinical data. Despite the lack of associations between E. coli genetic variants and patient outcomes, we estimate that increasing the sample size by one order of magnitude could lead to the discovery of some putative causal variants. Given the wide adoption of bacterial genome sequencing of clinical isolates, such sample sizes may be soon available.

Reduced chlorhexidine susceptibility is associated with tetracycline resistance tet genes in clinical isolates of Escherichia coli

Chlorhexidine is a widely used antiseptic in hospital and community healthcare. Decreased susceptibility to this compound has been recently described in Klebsiella pneumoniae and Pseudomonas aeruginosa, together with cross-resistance to colistin. Surprisingly, few data are available for Escherichia coli, the main species responsible for community and healthcare-associated infections. In order to decipher chlorhexidine resistance mechanisms in E. coli, we studied both in vitro derived and clinical isolates through whole-genome sequence analysis. Comparison of strains grown in vitro under chlorhexidine pressure identified mutations in the gene mlaA coding for a phospholipid transport system. Phenotypic analyses of single-gene mutant from the Keio collection confirmed the role of this mutation in the decreased susceptibility to chlorhexidine. However, mutations in mlaA were not found in isolates from large clinical collections. In contrast, genome wide association studies (GWAS) showed that, in clinical strains, chlorhexidine reduced susceptibility was associated with the presence of tetA genes of class B coding for efflux pumps and located in a Tn10 transposon. Construction of recombinant strains in E. coli K-12 confirmed the role of tetA determinant in acquired resistance to both chlorhexidine and tetracycline. Our results reveal two different evolutionary paths leading to chlorhexidine decreased susceptibility: one restricted to in vitro evolution conditions and involving a retrograde phospholipid transport system; the other observed in clinical isolates associated with efflux pump TetA. None of these mechanisms provides cross-resistance to colistin. This work demonstrates the GWAS power to identify new resistance mechanisms in bacterial species.

Genomic Analysis of Escherichia coli Longitudinally Isolated from Broiler Breeder Flocks after the Application of an Autogenous Vaccine

Escherichia coli is the main bacterial cause of major economic losses and animal welfare issues in poultry production. In this study, we investigate the effect of an autogenous vaccine on E. coli strains longitudinally isolated from broiler breeder flocks on two farms. In total, 115 E. coli isolates were sequenced using Illumina technologies, and compared based on a single-nucleotide polymorphism (SNP) analysis of the core-genome and antimicrobial resistance (AMR) genes they carried. The results showed that SNP-based phylogeny corresponds to a previous multilocus-sequence typing (MLST)-based phylogeny. Highly virulent sequence types (STs), including ST117-F, ST95-B2, ST131-B2 and ST390-B2, showed a higher level of homogeneity. On the other hand, less frequent STs, such as ST1485, ST3232, ST7013 and ST8573, were phylogenetically more distant and carried a higher number of antimicrobial resistance genes in most cases. In total, 25 antimicrobial genes were detected, of which the most prevalent were mdf(A) (100%), sitABCD (71.3%) and tet(A) (13.91%). The frequency of AMR genes showed a decreasing trend over time in both farms. The highest prevalence was detected in strains belonging to the B1 phylogenetic group, confirming the previous notion that commensal strains act as reservoirs and carry more resistance genes than pathogenic strains that are mostly associated with virulence genes.

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

F Plasmid Lineages in Escherichia coli ST95: Implications for Host Range, Antibiotic Resistance, and Zoonoses

Escherichia coli sequence type 95 (ST95) is an extraintestinal pathogenic E. coli (ExPEC) renowned for its ability to cause significant morbidity and mortality in humans and poultry. A core genome analysis of 668 ST95 isolates generated 10 clades (A to J), 5 of which are reported here for the first time. F plasmid replicon sequence typing showed that almost a third (178/668 [27%]) of the collection carry pUTI89 (F29:B10) and were restricted to clade A and a sublineage of clade B. In contrast, almost half (328/668 [49%]) of the collection across multiple clades harbor ColV plasmids (multiple F types). Strikingly, ST95 lineages with pUTI89 were almost exclusively from humans, while ColV⁺ ST95 lineages were sourced from poultry and humans. Clade I was notable because it comprises temporally and geographically matched ColV⁺ isolates sourced from human and retail poultry meat, suggesting interspecies transmission via food. Clade F contained ST95 isolates of bovine origin, none of which carried ColV or pUTI89 plasmids. Remarkably, an analysis of a cohort of 34,176 E. coli isolates comprising 2,570 sequence types mirrored what was observed in ST95: (i) pUTI89 was overwhelmingly linked to E. coli sourced from humans but almost entirely absent from 13,027 E. coli isolates recovered from poultry, pigs, and cattle, and (ii) E. coli isolates harboring ColV plasmids were from multiple sources, including humans, poultry, and swine. Overall, our data suggest that F plasmids influence E. coli host range, clade structure, and zoonotic potential in ST95 and ExPEC more broadly.

IMPORTANCEE. coli ST95 is one of five dominant ExPEC lineages globally and noted for causing urinary tract and bloodstream infections and neonatal meningitis in humans and colibacillosis in poultry. Using high-resolution phylogenomics, we show that F replicon sequence type is linked to ST95 clade structure and zoonotic potential. Specifically, human centric ST95 clades overwhelmingly harbor F29:B10 (pUTI89) plasmids, while clades carrying both human- and poultry-sourced isolates are typically ColV⁺ with multiple replicon types. Importantly, several clades identified clonal ColV⁺ ST95 isolates from human and poultry sources, but clade I, which housed temporally and spatially matched isolates, provided the most robust evidence. Notably, patterns of association of F replicon types with E. coli host were mirrored within a diverse collection of 34,176 E. coli genomes. Our studies indicate that the role of food animals as a source of human ExPEC disease is complex and warrants further investigation.

Increased microbial diversity and decreased prevalence of common pathogens in the gut microbiomes of wild turkeys compared to domestic turkeys

Turkeys (Meleagris gallopavo) provide a globally important source of protein and constitute the second most important source of poultry meat in the world. Bacterial diseases are common in commercial poultry production causing significant production losses for farmers. Due to the increasingly recognized problems associated with large-scale/indiscriminant antibiotic use in agricultural settings, poultry producers need alternative methods to control common bacterial pathogens. In this study we compared the cecal microbiota of wild and domestic turkeys, hypothesizing that environmental pressures faced by wild birds may select for a disease-resistant microbial community. Sequence analysis of 16S rRNA genes amplified from cecal samples indicate that free-roaming wild turkeys carry a rich and variable microbiota compared to domestic turkeys raised on large-scale poultry farms. Wild turkeys also had very low levels of Staphylococcus, Salmonella and E. coli when compared to domestic turkeys. E. coli strains isolated from wild or domestic turkey cecal samples also belong to distinct phylogenetic backgrounds and differ in their propensity to carry virulence genes. E. coli strains isolated from factory-raised turkeys were far more likely to carry genes for capsule (kpsII, kpsIII) or siderophore (iroN, fyuA) synthesis than those isolated from wild turkeys. These results suggest that the microbiota of wild turkeys may provide colonization resistance against common poultry pathogens. Importance Due to the increasingly recognized problems associated with antibiotic use in agricultural settings, poultry producers need alternative methods to control common bacterial pathogens. In this study we compare the microbiota of wild and domestic turkeys. Results suggest that free ranging wild turkeys carry a distinct microbiome when compared to farm raised turkeys. The microbiome of wild birds contains very low levels of poultry pathogens compared to farm raised birds. The microbiomes of wild turkeys may be used to guide development of new ways to control disease in large scale poultry production.

Functions predict horizontal gene transfer and the emergence of antibiotic resistance

Phylogenetic distance, shared ecology, and genomic constraints are often cited as key drivers governing horizontal gene transfer (HGT), although their relative contributions are unclear. Here, we apply machine learning algorithms to a curated set of diverse bacterial genomes to tease apart the importance of specific functional traits on recent HGT events. We find that functional content accurately predicts the HGT network [area under the receiver operating characteristic curve (AUROC) = 0.983], and performance improves further (AUROC = 0.990) for transfers involving antibiotic resistance genes (ARGs), highlighting the importance of HGT machinery, niche-specific, and metabolic functions. We find that high-probability not-yet detected ARG transfer events are almost exclusive to human-associated bacteria. Our approach is robust at predicting the HGT networks of pathogens, including Acinetobacter baumannii and Escherichia coli, as well as within localized environments, such as an individual’s gut microbiome.

In Vitro Characterization and In Vivo Efficacy Assessment in Galleria mellonella Larvae of Newly Isolated Bacteriophages against Escherichia coli K1

Extra-intestinal Escherichia coli express several virulence factors that increase their ability to colonize and survive in different localizations. The K1 capsular type is involved in several infections, including meningitis, urinary tract, and bloodstream infections. The aims of this work were to isolate, characterize, and assess the in vivo efficacy of phages targeting avian pathogenic E. coli (APEC) O18:K1, which shares many similarities with the human strains responsible for neonatal meningitis. Eleven phages were isolated against APEC O18:K1, and four of them presenting a narrow spectrum targeting E. coli K1 strains were further studied. The newly isolated phages vB_EcoS_K1-ULINTec2 were similar to the Siphoviridae family, and vB_EcoP_K1-ULINTec4, vB_EcoP_K1-ULINTec6, and vB_EcoP_K1-ULINTec7 to the Autographiviridae family. They are capsular type (K1) dependent and present several advantages characteristic of lytic phages, such as a short adsorption time and latent period. vB_EcoP_K1-ULINTec7 is able to target both K1 and K5 strains. This study shows that these phages replicate efficiently, both in vitro and in vivo in the Galleria mellonella model. Phage treatment increases the larvae survival rates, even though none of the phages were able to eliminate the bacterial load.

Flagellar rotor protein FliG is involved in the virulence of avian pathogenic Escherichia coli

Avian pathogenic Escherichia coli (APEC), a type of extraintestinal pathogenic E. coli, causes avian colibacillosis, a disease of significant economic importance to poultry producers worldwide, which is characterized by systemic infection. However, the pathogenesis of avian pathogenic E. coli strains is not well defined. Here, the role of a flagellar rotor protein encoded by the fliG gene of avian pathogenic E. coli strain AE17 was investigated. To study the role of FliG in the pathogenicity of APEC, fliG mutant and complemented strains were constructed and characterized. The inactivation of fliG had no effect on APEC growth, but significantly reduced bacterial motility. Compared with the wild type, the fliG mutant was highly attenuated in a chick infection model and showed severe defects in its adherence to and invasion of chicken embryo fibroblast DF-1 cells. The fliG mutant also showed reduced resistance to serum in chicks. The expression of the inflammatory cytokines interleukin 1β (IL1β), IL6, and IL8 was reduced in HD-11 macrophages infected with the fliG mutant strain compared with their expression in the wild-type strain. These results demonstrate that the FliG contributes to the virulence of APEC.

The E phylogroup of Escherichia coli is highly diverse and mimics the whole E. coli species population structure

To get a global picture of the population structure of the Escherichia coli phylogroup E, encompassing the O157:H7 EHEC lineage, we analysed the whole genome of 144 strains isolated from various continents, hosts and lifestyles and representative of the phylogroup diversity. The strains possess 4331 to 5440 genes with a core genome of 2771 genes and a pangenome of 33 722 genes. The distribution of these genes among the strains shows an asymmetric U-shaped distribution. E phylogenetic strains have the largest genomes of the species, partly explained by the presence of mobile genetic elements. Sixty-eight lineages were delineated, some of them exhibiting extra-intestinal virulence genes and being virulent in the mouse sepsis model. Except for the EHEC lineages and the reference EPEC, EIEC and ETEC strains, very few strains possess intestinal virulence genes. Most of the strains were devoid of acquired resistance genes, but eight strains possessed extended-spectrum beta-lactamase genes. Human strains belong to specific lineages, some of them being virulent and antibiotic-resistant [sequence type complexes (STcs) 350 and 2064]. The E phylogroup mimics all the features of the species as a whole, a phenomenon already observed at the STc level, arguing for a fractal population structure of E. coli.

From farm to fork: Colistin voluntary withdrawal in Portuguese farms reflected in decreasing occurrence of mcr-1-carrying Enterobacteriaceae from chicken meat

Expansion of mcr-carrying Enterobacteriaceae (MCR-E) is a well-recognized problem affecting animals, humans and the environment. Ongoing global control actions involve colistin restrictions among food-animal production, but their impact on poultry-derived products is largely unknown, justifying comprehensive farm-to-fork studies. Occurrence of MCR-E among 53 chicken-meat batches supplied from 29 Portuguese farms shortly after colistin withdrawal was evaluated. Strains (FT-IR/MLST/WGS), mcr plasmids and their adaptive features were characterized by cultural, molecular and genomic approaches. We found high rates of chicken-meat batches (80%-100% - 4 months; 12% - the last month) with multiple MDR + mcr-1-carrying Escherichia coli (Ec-including ST117 and ST648-Cplx) and Klebsiella pneumoniae (Kp-ST147-O5:K35) clones, some of them persisting over time. The mcr-1 was located in the chromosome (Ec-ST297/16-farms) or dispersed IncX4 (Ec-ST602/ST6469/5-farms), IncHI2-ST2/ST4 (Ec-ST533/ST6469/5 farms and Kp-ST147/6-farms) or IncI2 (Ec-ST117/1-farm) plasmids. WGS revealed high load and diversity in virulence, antibiotic resistance and metal tolerance genes. This study supports colistin withdrawal potential efficacy in poultry production and highlights both poultry-production chain as a source of mcr-1 and the risk of foodborne transmission to poultry-meat consumers. Finally, in the antibiotic reduction/replacement context, other potential co-selective pressures (e.g., metals-Cu as feed additives) need to be further understood to guide concerted, effective and durable actions under 'One Health' perspective.