PT  - JOURNAL ARTICLE
AU  - Madden, Danielle E.
AU  - Baird, Timothy
AU  - Bell, Scott C.
AU  - McCarthy, Kate L.
AU  - Price, Erin P.
AU  - Sarovich, Derek S.
TI  - Keeping up with the pathogens: Improved antimicrobial resistance detection and prediction from <em>Pseudomonas</em> aeruginosa genomes
AID  - 10.1101/2022.08.11.22278689
DP  - 2023 Jan 01
TA  - medRxiv
PG  - 2022.08.11.22278689
4099  - http://medrxiv.org/content/early/2023/10/31/2022.08.11.22278689.short
4100  - http://medrxiv.org/content/early/2023/10/31/2022.08.11.22278689.full
AB  - Background Antimicrobial resistance (AMR) is an intensifying threat that requires urgent mitigation to avoid a post-antibiotic era. The ESKAPE pathogen, Pseudomonas aeruginosa, represents one of the greatest AMR concerns due to increasing multi- and pan-drug resistance rates. Shotgun sequencing is quickly gaining traction for in silico AMR profiling due to its unambiguity and transferability; however, accurate and comprehensive AMR prediction from P. aeruginosa genomes remains an unsolved problem.Methods We first curated the most comprehensive database yet of known P. aeruginosa AMR variants. Next, we performed comparative genomics and microbial genome-wide association study analysis across a Global isolate Dataset (n=1877) with paired antimicrobial phenotype and genomic data to identify novel AMR variants. Finally, the performance of our P. aeruginosa AMR database, implemented in our ARDaP software, was compared with three previously published in silico AMR gene detection or phenotype prediction tools – abritAMR, AMRFinderPlus, ResFinder – across both the Global Dataset and an analysis-naïve Validation Dataset (n=102).Results Our AMR database comprises 3639 mobile AMR genes and 733 AMR-conferring chromosomal variants, including 75 chromosomal variants not previously reported, and 284 chromosomal variants that we show are unlikely to confer AMR. Our pipeline achieved a genotype-phenotype balanced accuracy (bACC) of 85% and 81% across 10 clinically relevant antibiotics when tested against the Global and Validation Datasets, respectively, vs. just 56% and 54% with abritAMR, 58% and 54% with AMRFinderPlus, and 60% and 53% with ResFinder.Conclusions Our ARDaP software and associated AMR variant database provides the most accurate tool yet for predicting AMR phenotypes in P. aeruginosa, far surpassing the performance of current tools. Implementation of our ARDaP-compatible database for routine AMR prediction from P. aeruginosa genomes and metagenomes will improve AMR identification, addressing a critical facet in combatting this treatment-refractory pathogen. However, knowledge gaps remain in our understanding of the P. aeruginosa resistome, particularly the basis of colistin AMR.Competing Interest StatementThe authors have declared no competing interest.Funding StatementThis work was supported by Advance Queensland (AQIRF0362018 and AQRF13016-17RD2), the Wishlist Sunshine Coast Health Foundation (2019-14), and the National Health and Medical Research Council (GNT455919). The funders had no role in study design, data acquisition, analysis, interpretation, writing, or manuscript submission.Author DeclarationsI confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.YesThe details of the IRB/oversight body that provided approval or exemption for the research described are given below:Clinical specimen collection was approved by The Prince Charles Hospital (HREC/13/QPCH/127 and HREC/2019/QPCH/48013) and Royal Brisbane and Womens Hospital (DA: jl) Human Research Ethics Committees. Participants with cystic fibrosis or chronic obstructive pulmonary disease provided written consent. A waiver of informed consent was approved for other specimens due to a low-to-negligible risk assessment.I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.YesI understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).Yes I have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.YesThe datasets supporting the conclusions of this article are available in the NCBI GenBank (references NSXK00000000.1 and NSZK00000000.1) and Sequence Read Archive database via BioProject accession PRJNA761496 (Validation Dataset), and BioProject accessions PRJEB15036, PRJEB21341, PRJEB29539, PRJEB14771, PRJNA297679, PRJNA317143, PRJNA526797, PRJNA264310, PRJNA388357, PRJNA793523, PRJEB40140 and PRJNA532924 (Global Dataset). The P. aeruginosa AMR database constructed for this study is freely available at: https://github.com/dsarov/ARDaP/tree/master/Databases/Pseudomonas_aeruginosa_pao1. The ARDaP software is freely available at: https://github.com/dsarov/ARDaP. All code used in this manuscript is available at https://github.com/dsarov/P_aeruginosa_ARDaP_manuscript.The datasets supporting the conclusions of this article are available in the NCBI GenBank (references NSXK00000000.1 and NSZK00000000.1) and Sequence Read Archive database via BioProject accession PRJNA761496 (Validation Dataset), and BioProject accessions PRJEB15036, PRJEB21341, PRJEB29539, PRJEB14771, PRJNA297679, PRJNA317143, PRJNA526797, PRJNA264310, PRJNA388357, PRJNA793523, PRJEB40140 and PRJNA532924 (Global Dataset). The P. aeruginosa AMR database constructed for this study is freely available at: https://github.com/dsarov/ARDaP/tree/master/Databases/Pseudomonas_aeruginosa_pao1. The ARDaP software is freely available at: https://github.com/dsarov/ARDaP. All code used in this manuscript is available at https://github.com/dsarov/P_aeruginosa_ARDaP_manuscript.