Abstract:Background: Antimicrobial resistance (AMR) in Escherichia coli is a global problem associated with substantial morbidity and mortality. AMR-associated genes are typically annotated based on similarity to variants in a curated reference database, with the implicit assumption that uncatalogued genetic variation within these is phenotypically unimportant. In this study, we evaluated the performance of the AMRFinder tool and, subsequently, the potential for discovering new AMR-associated gene families and characterising variation within existing ones to improve genotype-to-susceptibility phenotype predictions in E coli. Methods: In this cross-sectional study of international genome sequence data, we assembled a global dataset of 9001 E coli sequences from five publicly available data collections predominantly deriving from human bloodstream infections from: Norway, Oxfordshire (UK), Thailand, the UK, and Sweden. 8555 of these sequences had linked antibiotic susceptibility data. Raw reads were assembled using Shovill and AMR genes (relevant to amoxicillin-clavulanic acid, ampicillin, ceftriaxone, ciprofloxacin, gentamicin, piperacillin-tazobactam, and trimethoprim) extracted using the National Center for Biotechnology Information AMRFinder tool (using both default and strict [100%] coverage and identity filters). We assessed the predictive value of the presence of these genes for predicting resistance or susceptibility against US Food and Drug Administration thresholds for major and very major errors. Mash was used to calculate the similarity between extracted genes using Jaccard distances. We empirically reclustered extracted gene sequences into AMR-associated gene families (≥70% match) and antibiotic-resistance genes (ARGs; 100% match) and categorised these according to their frequency in the dataset. Accumulation curves were simulated and correlations between gene frequency in the Oxfordshire and other datasets calculated using the Spearman coefficient. Firth regression was used to model the association between the presence of blaTEM-1 variants and amoxicillin-clavulanic acid or piperacillin-tazobactam resistance, adjusted for the presence of other relevant ARGs. Findings: The performance of the AMRFinder database for genotype-to-phenotype predictions using strict 100% identity and coverage thresholds did not meet US Food and Drug Administration thresholds for any of the seven antibiotics evaluated. Relaxing filters to default settings improved sensitivity with a specificity cost. For all antibiotics, most explainable resistance was associated with the presence of a small number of genes. There was a proportion of resistance that could not be explained by known ARGs; this ranged from 75·1% for amoxicillin-clavulanic acid to 3·4% for ciprofloxacin. Only 18 199 (51·5%) of the 35 343 ARGs detected had a 100% identity and coverage match in the AMRFinder database. After empirically reclassifying genes at 100% nucleotide sequence identity, we identified 1042 unique ARGs, of which 126 (12·1%) were present ten times or more, 313 (30·0%) were present between two and nine times, and 603 (57·9%) were present only once. Simulated accumulation curves revealed that discovery of new (100% match) ARGs present more than once in the dataset plateaued relatively quickly, whereas new singleton ARGs were discovered even after many thousands of isolates had been included. We identified a strong correlation (Spearman coefficient 0·76 [95% CI 0·73-0·80], p<0·0001) between the number of times an ARG was observed in Oxfordshire and the number of times it was seen internationally, with ARGs that were observed six times in Oxfordshire always being found elsewhere. Finally, using the example of blaTEM-1, we showed that uncatalogued variation, including synonymous variation, is associated with potentially important phenotypic differences; for example, two common, uncatalogued blaTEM-1 alleles with only synonymous mutations compared with the known reference were associated with reduced resistance to amoxicillin-clavulanic acid (adjusted odds ratio 0·58 [95% CI 0·35-0·95], p=0·031) and piperacillin-tazobactam (0·50 [95% CI 0·29-0·82], p=0·005). Interpretation: We highlight substantial uncatalogued genetic variation with respect to known ARGs, although a relatively small proportion of these alleles are repeatedly observed in a large international dataset suggesting strong selection pressures. The current approach of using fuzzy matching for ARG detection, ignoring the unknown effects of uncatalogued variation, is unlikely to be acceptable for future clinical deployment. The association of synonymous mutations with potentially important phenotypic differences suggests that relying solely on amino acid-based gene detection to predict resistance is unlikely to be sufficient. Finally, the inability to explain all resistance using existing knowledge highlights the importance of new target gene discovery. Funding: National Institute for Health and Care Research, Wellcome, and UK Medical Research Council.

Exploring antimicrobial resistance determinants in the Neanderthal microbiome

Identifying ancient antibiotic resistance genes in archaeological dental calculus

Profiling Antibiotic Resistance Determinants in Ancient Permafrost Microbiomes

Inferring diet, disease and antibiotic resistance from ancient human oral microbiomes

Unraveling the Microbiotas and Key Genetic Contexts Identified on Stone Heritage Using Illumina and Nanopore Sequencing Platforms

Neanderthal behaviour, diet, and disease inferred from ancient DNA in dental calculus

Revealing antibiotic resistance's ancient roots: insights from pristine ecosystems

Metagenome-wide Analysis of Antibiotic Resistance Genes in a Large Cohort of Human Gut Microbiota

Clinical considerations on antimicrobial resistance potential of complex microbiological samples

Ancient Genomes From Bronze Age Remains Reveal Deep Diversity and Recent Adaptive Episodes for Human Oral Pathobionts

The past, present and future of ancient bacterial DNA

Metagenomic Analysis Reveals the Distribution of Antibiotic Resistance Genes in a Large-Scale Population of Healthy Individuals and Patients With Varied Diseases

Appearance of β-lactam Resistance Genes in Agricultural Soils and Clinical Isolates over the 20th Century

Unravelling the complex interplay between antibiotic consumption and adaptive changes in methicillin-resistant Staphylococcus aureus

A Metagenomic Study of Antibiotic Resistance Across Diverse Soil Types and Geographical Locations

Temporal dynamics and persistence of resistance genes to broad spectrum antibiotics in an urban community

Within-host diversity of MRSA antimicrobial resistances

Metagenomic-based surveillance systems for antibiotic resistance in non-clinical settings

Rolling the evolutionary dice: Neisseria commensals as proxies for elucidating the underpinnings of antibiotic resistance mechanisms and evolution in human pathogens

Exploring uncatalogued genetic variation in antimicrobial resistance gene families in Escherichia coli: an observational analysis

Deeper Exploration of Gut Microbiome: Profile of Resistome, Virome and Viral Auxiliary Metabolic Genes of Three Ethnic Indian Groups