Metagenomics-Based Approach to Source-Attribution of Antimicrobial Resistance Determinants – Identification of Reservoir Resistome Signatures

Metagenomics can unveil the genetic content of the total microbiota in different environments, such as food products and the guts of humans and livestock. It is therefore considered of great potential to investigate the transmission of foodborne hazards as part of source-attribution studies. Source-attribution of antimicrobial resistance (AMR) has traditionally relied on pathogen isolation, while metagenomics allows investigating the full span of AMR determinants. In this study, we hypothesized that the relative abundance of fecal resistome components can be associated with specific reservoirs, and that resistomes can be used for AMR source-attribution. We used shotgun-sequences from fecal samples of pigs, broilers, turkeys- and veal calves collected across Europe, and fecal samples from humans occupationally exposed to livestock in one country (pig slaughterhouse workers, pig and broiler farmers). We applied both hierarchical and flat forms of the supervised classification ensemble algorithm Random Forests to classify resistomes into corresponding reservoir classes. We identified country-specific and -independent AMR determinants, and assessed the impact of country-specific determinants when attributing AMR resistance in humans. Additionally, we performed a similarity percentage analysis with the full spectrum of AMR determinants to identify resistome signatures for the different reservoirs. We showed that the number of AMR determinants necessary to attribute a resistome into the correct reservoir increases with a larger reservoir heterogeneity, and that the impact of country-specific resistome signatures on prediction varies between countries. We predicted a higher occupational exposure to AMR determinants among workers exposed to pigs than among those exposed to broilers. Additionally, results suggested that AMR exposure on pig farms was higher than in pig slaughterhouses. Human resistomes were more similar to pig and veal calves’ resistomes than to those of broilers and turkeys, and the majority of these resistome dissimilarities can be explained by a small set of AMR determinants. We identified resistome signatures for each individual reservoir, which include AMR determinants significantly associated with on-farm antimicrobial use. We attributed human resistomes to different livestock reservoirs using Random Forests, which allowed identifying pigs as a potential source of AMR in humans. This study thus demonstrates that it is possible to apply metagenomics in AMR source-attribution.

[1]  F. Aarestrup,et al.  Description and determinants of the faecal resistome and microbiome of farmers and slaughterhouse workers: A metagenome-wide cross-sectional study. , 2020, Environment international.

[2]  P. Njage,et al.  Application of Whole‐Genome Sequences and Machine Learning in Source Attribution of Salmonella Typhimurium , 2020, Risk analysis : an official publication of the Society for Risk Analysis.

[3]  Liqing Zhang,et al.  Identification of discriminatory antibiotic resistance genes among environmental resistomes using extremely randomized tree algorithm , 2019, Microbiome.

[4]  T. Hald,et al.  Attributable sources of community-acquired carriage of Escherichia coli containing β-lactam antibiotic resistance genes: a population-based modelling study. , 2019, The Lancet. Planetary health.

[5]  F. Aarestrup,et al.  Associations between antimicrobial use and the faecal resistome on broiler farms from nine European countries , 2019, The Journal of antimicrobial chemotherapy.

[6]  F. Aarestrup,et al.  The antimicrobial resistome in relation to antimicrobial use and biosecurity in pig farming, a metagenome-wide association study in nine European countries , 2019, The Journal of antimicrobial chemotherapy.

[7]  Jukka Corander,et al.  Fast and flexible bacterial genomic epidemiology with PopPUNK , 2018, bioRxiv.

[8]  Beau B. Bruce,et al.  Zoonotic Source Attribution of Salmonella enterica Serotype Typhimurium Using Genomic Surveillance Data, United States , 2019, Emerging infectious diseases.

[9]  M. Dec,et al.  Identification and antibiotic susceptibility of lactobacilli isolated from turkeys , 2018, BMC Microbiology.

[10]  J. Corander,et al.  Fast hierarchical Bayesian analysis of population structure , 2018, bioRxiv.

[11]  S. Khorram,et al.  Supervised Classification , 2018, Image Processing and Data Analysis with ERDAS IMAGINE®.

[12]  O. Lund,et al.  Abundance and diversity of the faecal resistome in slaughter pigs and broilers in nine European countries , 2018, Nature Microbiology.

[13]  Xu Li,et al.  Characterization of Natural and Affected Environments Tracking the Sources of Antibiotic Resistance Genes in an Urban Stream during Wet Weather using Shotgun Metagenomic Analyses , 2018 .

[14]  S. Sheppard,et al.  Ruminant and chicken: important sources of campylobacteriosis in France despite a variation of source attribution in 2009 and 2015 , 2018, Scientific Reports.

[15]  T. Hald,et al.  Source Attribution and Risk Assessment of Antimicrobial Resistance. , 2018, Microbiology spectrum.

[16]  Tong Zhang,et al.  Tracking antibiotic resistance gene pollution from different sources using machine-learning classification , 2018, Microbiome.

[17]  Palma Blonda,et al.  Comparing the performance of flat and hierarchical Habitat/Land-Cover classification models in a NATURA 2000 site , 2018 .

[18]  W. van Pelt,et al.  New paradigms for Salmonella source attribution based on microbial subtyping. , 2017, Food microbiology.

[19]  M. Dec,et al.  Assessment of antibiotic susceptibility in Lactobacillus isolates from chickens , 2017, Gut Pathogens.

[20]  Ole Lund,et al.  MGmapper: Reference based mapping and taxonomy annotation of metagenomics sequence reads , 2017, PloS one.

[21]  M. Hipsey,et al.  Source tracking using microbial community fingerprints: Method comparison with hydrodynamic modelling. , 2017, Water research.

[22]  Jurgen E. Chardon,et al.  Comparative Exposure Assessment of ESBL-Producing Escherichia coli through Meat Consumption , 2017, PloS one.

[23]  D. Mccarthy,et al.  Into the deep: Evaluation of SourceTracker for assessment of faecal contamination of coastal waters. , 2016, Water research.

[24]  Mrudula Patel,et al.  Detection of cfxA2, cfxA3, and cfxA6 genes in beta-lactamase producing oral anaerobes , 2016, Journal of applied oral science : revista FOB.

[25]  T. Hald,et al.  Application of Molecular Typing Results in Source Attribution Models: The Case of Multiple Locus Variable Number Tandem Repeat Analysis (MLVA) of Salmonella Isolates Obtained from Integrated Surveillance in Denmark , 2016, Risk analysis : an official publication of the Society for Risk Analysis.

[26]  W. van Pelt,et al.  Salmonella source attribution based on microbial subtyping: does including data on food consumption matter? , 2014, International journal of food microbiology.

[27]  Tine Hald,et al.  Source attribution of human salmonellosis: an overview of methods and estimates. , 2014, Foodborne pathogens and disease.

[28]  M. Stephens,et al.  fastSTRUCTURE: Variational Inference of Population Structure in Large SNP Data Sets , 2014, Genetics.

[29]  Jukka Corander,et al.  Hierarchical and Spatially Explicit Clustering of DNA Sequences with BAPS Software , 2013, Molecular biology and evolution.

[30]  Zhengwei Zhu,et al.  CD-HIT: accelerated for clustering the next-generation sequencing data , 2012, Bioinform..

[31]  S. Rasmussen,et al.  Identification of acquired antimicrobial resistance genes , 2012, The Journal of antimicrobial chemotherapy.

[32]  Rob Knight,et al.  Bayesian community-wide culture-independent microbial source tracking , 2011, Nature Methods.

[33]  R. Knight,et al.  Supervised classification of human microbiota. , 2011, FEMS microbiology reviews.

[34]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[35]  F. Aarestrup,et al.  The attribution of human infections with antimicrobial resistant Salmonella bacteria in Denmark to sources of animal origin. , 2007, Foodborne pathogens and disease.

[36]  Adam Godzik,et al.  Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences , 2006, Bioinform..

[37]  Tine Hald,et al.  A Bayesian Approach to Quantify the Contribution of Animal‐Food Sources to Human Salmonellosis , 2004, Risk analysis : an official publication of the Society for Risk Analysis.

[38]  Leo Breiman,et al.  Random Forests , 2001, Machine Learning.

[39]  A. Gupta,et al.  A Bayesian Approach to , 1997 .

[40]  K. R. Clarke,et al.  Non‐parametric multivariate analyses of changes in community structure , 1993 .

[41]  J. H. Ward Hierarchical Grouping to Optimize an Objective Function , 1963 .

[42]  Jacob Cohen A Coefficient of Agreement for Nominal Scales , 1960 .