Abstract
Antimicrobial resistance prediction from whole genome sequencing data (WGS) is an emerging application of machine learning, promising to improve antimicrobial resistance surveillance and outbreak monitoring. Despite significant reductions in sequencing cost, the availability and sampling diversity of WGS data with matched antimicrobial susceptibility testing (AST) profiles required for training of WGS-AST prediction models remains limited. Best practice machine learning techniques are required to ensure trained models generalize to independent data for optimal predictive performance. Limited data restricts the choice of machine learning training and evaluation methods and can result in overestimation of model performance. We demonstrate that the widely used random k-fold cross-validation method is ill-suited for application to small bacterial genomics datasets and offer an alternative cross-validation method based on genomic distance. We benchmarked three machine learning architectures previously applied to the WGS-AST problem on a set of 8,704 genome assemblies from five clinically relevant pathogens across 77 species-compound combinations collated from public databases. We show that individual models can be effectively ensembled to improve model performance. By combining models via stacked generalization with cross-validation, a model ensembling technique suitable for small datasets, we improved average sensitivity and specificity of individual models by 1.77% and 3.20%, respectively. Furthermore, stacked models exhibited improved robustness and were thus less prone to outlier performance drops than individual component models. In this study, we highlight best practice techniques for antimicrobial resistance prediction from WGS data and introduce the combination of genome distance aware cross-validation and stacked generalization for robust and accurate WGS-AST.
Highlights
Antimicrobial resistance (AMR) is a rising global threat to human health
Other representations of genomic data, such as amino acid k-mers or protein variants have been used for WGSAST model training as well (Kim et al, 2020; Valizadehaslani et al, 2020)
Genome assemblies coupled with antimicrobial susceptibility testing (AST) information were obtained from public databases for five human pathogens (A. baumannii, E. coli, P. aeruginosa, K. pneumoniae and S. aureus) and a total set of 77 organism/compound combinations
Summary
Antimicrobial resistance (AMR) is a rising global threat to human health. To ensure the continued efficacy of antimicrobial compounds, prudent use of this resource is crucial (O’Neill, 2016). Described WGS-AST techniques use nucleotide k-mer representations of genome assemblies or raw sequencing data, attempting to learn differences in k-mer counts or presence/ absence patterns that correlate with shifts in susceptibility to a target antibiotic (Drouin et al, 2016; Aun et al, 2018; Nguyen et al, 2018a; Drouin et al, 2019). This data-driven approach does not require expert knowledge of AMR mechanisms or prior information on AMR genes, and can be applied towards learning of models for novel antibiotics and unknown resistance mechanisms. Other representations of genomic data, such as amino acid k-mers or protein variants have been used for WGSAST model training as well (Kim et al, 2020; Valizadehaslani et al, 2020)
Published Version (Free)
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have