Keeping up with the pathogens: improved antimicrobial resistance detection and prediction from Pseudomonas aeruginosa genomes.

BACKGROUND: Antimicrobial resistance (AMR) is an intensifying threat that requires urgent mitigation to avoid a post-antibiotic era. Pseudomonas aeruginosa represents one of the greatest AMR concerns due to increasing multi- and pan-drug resistance rates. Shotgun sequencing is 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 AMR detection and prediction tool, ARDaP, 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 728 chromosomal variants, including 75 previously unreported chromosomal AMR variants, 10 variants associated with unusual antimicrobial susceptibility, and 281 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. ARDaP's superior performance was predominantly due to the inclusion of chromosomal AMR variants, which are generally not identified with most AMR identification tools. CONCLUSIONS: Our ARDaP software and associated AMR variant database provides an accurate tool for predicting AMR phenotypes in P. aeruginosa, far surpassing the performance of current tools. Implementation of ARDaP 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.

Antimicrobial susceptibility prediction from genomes: a dream come true?

Genome-based diagnostics provides relevant information to guide patient treatment and support pathogen and resistance surveillance. Recently, Coll et al. introduced a curated database for predicting antimicrobial resistance (AMR) from Enterococcus faecium genomics data, offering excellent predictive values for susceptibility to important antimicrobials. Challenges to predict resistance to last-resort antimicrobials remain.

Three Distinct Annotation Platforms Differ in Detection of Antimicrobial Resistance Genes in Long-Read, Short-Read, and Hybrid Sequences Derived from Total Genomic DNA or from Purified Plasmid DNA.

Recent advances and lower costs in rapid high-throughput sequencing have engendered hope that whole genome sequencing (WGS) might afford complete resistome characterization in bacterial isolates. WGS is particularly useful for the clinical characterization of fastidious and slow-growing bacteria. Despite its potential, several challenges should be addressed before adopting WGS to detect antimicrobial resistance (AMR) genes in the clinical laboratory. Here, with three distinct ESKAPE bacteria (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter spp.), different approaches were compared to identify best practices for detecting AMR genes, including: total genomic DNA and plasmid DNA extractions, the solo assembly of Illumina short-reads and of Oxford Nanopore Technologies (ONT) long-reads, two hybrid assembly pipelines, and three in silico AMR databases. We also determined the susceptibility of each strain to 21 antimicrobials. We found that all AMR genes detected in pure plasmid DNA were also detectable in total genomic DNA, indicating that, at least in these three enterobacterial genera, the purification of plasmid DNA was not necessary to detect plasmid-borne AMR genes. Illumina short-reads used with ONT long-reads in either hybrid or polished assemblies of total genomic DNA enhanced the sensitivity and accuracy of AMR gene detection. Phenotypic susceptibility closely corresponded with genotypes identified by sequencing; however, the three AMR databases differed significantly in distinguishing mobile dedicated AMR genes from non-mobile chromosomal housekeeping genes in which rare spontaneous resistance mutations might occur. This study indicates that each method employed in a WGS workflow has an impact on the detection of AMR genes. A combination of short- and long-reads, followed by at least three different AMR databases, should be used for the consistent detection of such genes. Further, an additional step for plasmid DNA purification and sequencing may not be necessary. This study reveals the need for standardized biochemical and informatic procedures and database resources for consistent, reliable AMR genotyping to take full advantage of WGS in order to expedite patient treatment and track AMR genes within the hospital and community.

Correlation between Phenotypic and In Silico Detection of Antimicrobial Resistance in Salmonella enterica in Canada Using Staramr.

Whole genome sequencing (WGS) of Salmonella supports both molecular typing and detection of antimicrobial resistance (AMR). Here, we evaluated the correlation between phenotypic antimicrobial susceptibility testing (AST) and in silico prediction of AMR from WGS in Salmonella enterica (n = 1321) isolated from human infections in Canada. Phenotypic AMR results from broth microdilution testing were used as the gold standard. To facilitate high-throughput prediction of AMR from genome assemblies, we created a tool called Staramr, which incorporates the ResFinder and PointFinder databases and a custom gene-drug key for antibiogram prediction. Overall, there was 99% concordance between phenotypic and genotypic detection of categorical resistance for 14 antimicrobials in 1321 isolates (18,305 of 18,494 results in agreement). We observed an average sensitivity of 91.2% (range 80.5-100%), a specificity of 99.7% (98.6-100%), a positive predictive value of 95.4% (68.2-100%), and a negative predictive value of 99.1% (95.6-100%). The positive predictive value of gentamicin was 68%, due to seven isolates that carried aac(3)-IVa, which conferred MICs just below the breakpoint of resistance. Genetic mechanisms of resistance in these 1321 isolates included 64 unique acquired alleles and mutations in three chromosomal genes. In general, in silico prediction of AMR in Salmonella was reliable compared to the gold standard of broth microdilution. WGS can provide higher-resolution data on the epidemiology of resistance mechanisms and the emergence of new resistance alleles.

Antimicrobial Resistance Profiling and Phylogenetic Analysis of Neisseria gonorrhoeae Clinical Isolates From Kenya in a Resource-Limited Setting.

BACKGROUND: Africa has one of the highest incidences of gonorrhea. Neisseria gonorrhoeae is gaining resistance to most of the available antibiotics, compromising treatment across the world. Whole-genome sequencing (WGS) is an efficient way of predicting AMR determinants and their spread in the population. Recent advances in next-generation sequencing technologies like Oxford Nanopore Technology (ONT) have helped in the generation of longer reads of DNA in a shorter duration with lower cost. Increasing accuracy of base-calling algorithms, high throughput, error-correction strategies, and ease of using the mobile sequencer MinION in remote areas lead to its adoption for routine microbial genome sequencing. To investigate whether MinION-only sequencing is sufficient for WGS and downstream analysis in resource-limited settings, we sequenced the genomes of 14 suspected N. gonorrhoeae isolates from Nairobi, Kenya. METHODS: Using WGS, the isolates were confirmed to be cases of N. gonorrhoeae (n = 9), and there were three co-occurrences of N. gonorrhoeae with Moraxella osloensis and N. meningitidis (n = 2). N. meningitidis has been implicated in sexually transmitted infections in recent years. The near-complete N. gonorrhoeae genomes (n = 10) were analyzed further for mutations/factors causing AMR using an in-house database of mutations curated from the literature. RESULTS: We observe that ciprofloxacin resistance is associated with multiple mutations in both gyrA and parC. Mutations conferring tetracycline (rpsJ) and sulfonamide (folP) resistance and plasmids encoding beta-lactamase were seen in all the strains, and tet(M)-containing plasmids were identified in nine strains. Phylogenetic analysis clustered the 10 isolates into clades containing previously sequenced genomes from Kenya and countries across the world. Based on homology modeling of AMR targets, we see that the mutations in GyrA and ParC disrupt the hydrogen bonding with quinolone drugs and mutations in FolP may affect interaction with the antibiotic. CONCLUSION: Here, we demonstrate the utility of mobile DNA sequencing technology in producing a consensus genome for sequence typing and detection of genetic determinants of AMR. The workflow followed in the study, including AMR mutation dataset creation and the genome identification, assembly, and analysis, can be used for any clinical isolate. Further studies are required to determine the utility of real-time sequencing in outbreak investigations, diagnosis, and management of infections, especially in resource-limited settings.

PARMAP: A Pan-Genome-Based Computational Framework for Predicting Antimicrobial Resistance.

Antimicrobial resistance (AMR) has emerged as one of the most urgent global threats to public health. Accurate detection of AMR phenotypes is critical for reducing the spread of AMR strains. Here, we developed PARMAP (Prediction of Antimicrobial Resistance by MAPping genetic alterations in pan-genome) to predict AMR phenotypes and to identify AMR-associated genetic alterations based on the pan-genome of bacteria by utilizing machine learning algorithms. When we applied PARMAP to 1,597 Neisseria gonorrhoeae strains, it successfully predicted their AMR phenotypes based on a pan-genome analysis. Furthermore, it identified 328 genetic alterations in 23 known AMR genes and discovered many new AMR-associated genetic alterations in ciprofloxacin-resistant N. gonorrhoeae, and it clearly indicated the genetic heterogeneity of AMR genes in different subtypes of resistant N. gonorrhoeae. Additionally, PARMAP performed well in predicting the AMR phenotypes of Mycobacterium tuberculosis and Escherichia coli, indicating the robustness of the PARMAP framework. In conclusion, PARMAP not only precisely predicts the AMR of a population of strains of a given species but also uses whole-genome sequencing data to prioritize candidate AMR-associated genetic alterations based on their likelihood of contributing to AMR. Thus, we believe that PARMAP will accelerate investigations into AMR mechanisms in other human pathogens.

Integrative Analysis of Whole Genome Sequencing and Phenotypic Resistance Toward Prediction of Trimethoprim-Sulfamethoxazole Resistance in Staphylococcus aureus.

As whole genome sequencing is becoming more accessible and affordable for clinical microbiological diagnostics, the reliability of genotypic antimicrobial resistance (AMR) prediction from sequencing data is an important issue to address. Computational AMR prediction can be performed at multiple levels. The first-level approach, such as simple AMR search relies heavily on the quality of the information fed into the database. However, AMR due to mutations are often undetected, since this is not included in the database or poorly documented. Using co-trimoxazole (trimethoprim-sulfamethoxazole) resistance in Staphylococcus aureus, we compared single-level and multi-level analysis to investigate the strengths and weaknesses of both approaches. The results revealed that a single mutation in the AMR gene on the nucleotide level may produce false positive results, which could have been detected if protein sequence analysis would have been performed. For AMR predictions based on chromosomal mutations, such as the folP gene of S. aureus, natural genetic variations should be taken into account to differentiate between variants linked to genetic lineage (MLST) and not over-estimate the potential resistant variants. Our study showed that careful analysis of the whole genome data and additional criterion such as lineage-independent mutations may be useful for identification of mutations leading to phenotypic resistance. Furthermore, the creation of reliable database for point mutations is needed to fully automatized AMR prediction.

Identification of Primary Antimicrobial Resistance Drivers in Agricultural Nontyphoidal Salmonella enterica Serovars by Using Machine Learning.

Nontyphoidal Salmonella (NTS) is a leading global cause of bacterial foodborne morbidity and mortality. Our ability to treat severe NTS infections has been impaired by increasing antimicrobial resistance (AMR). To understand and mitigate the global health crisis AMR represents, we need to link the observed resistance phenotypes with their underlying genomic mechanisms. Broiler chickens represent a key reservoir and vector for NTS infections, but isolates from this setting have been characterized in only very low numbers relative to clinical isolates. In this study, we sequenced and assembled 97 genomes encompassing 7 serotypes isolated from broiler chicken in farms in British Columbia between 2005 and 2008. Through application of machine learning (ML) models to predict the observed AMR phenotype from this genomic data, we were able to generate highly (0.92 to 0.99) precise logistic regression models using known AMR gene annotations as features for 7 antibiotics (amoxicillin-clavulanic acid, ampicillin, cefoxitin, ceftiofur, ceftriaxone, streptomycin, and tetracycline). Similarly, we also trained "reference-free" k-mer-based set-covering machine phenotypic prediction models (0.91 to 1.0 precision) for these antibiotics. By combining the inferred k-mers and logistic regression weights, we identified the primary drivers of AMR for the 7 studied antibiotics in these isolates. With our research representing one of the largest studies of a diverse set of NTS isolates from broiler chicken, we can thus confirm that the AmpC-like CMY-2 ß-lactamase is a primary driver of ß-lactam resistance and that the phosphotransferases APH(6)-Id and APH(3″-Ib) are the principal drivers of streptomycin resistance in this important ecosystem.IMPORTANCE Antimicrobial resistance (AMR) represents an existential threat to the function of modern medicine. Genomics and machine learning methods are being increasingly used to analyze and predict AMR. This type of surveillance is very important to try to reduce the impact of AMR. Machine learning models are typically trained using genomic data, but the aspects of the genomes that they use to make predictions are rarely analyzed. In this work, we showed how, by using different types of machine learning models and performing this analysis, it is possible to identify the key genes underlying AMR in nontyphoidal Salmonella (NTS). NTS is among the leading cause of foodborne illness globally; however, AMR in NTS has not been heavily studied within the food chain itself. Therefore, in this work we performed a broad-scale analysis of the AMR in NTS isolates from commercial chicken farms and identified some priority AMR genes for surveillance.

Pan-Genomic and Polymorphic Driven Prediction of Antibiotic Resistance in Elizabethkingia.

The Elizabethkingia are a genetically diverse genus of emerging pathogens that exhibit multidrug resistance to a range of common antibiotics. Two representative species, Elizabethkingia bruuniana and E. meningoseptica, were phenotypically tested to determine minimum inhibitory concentrations (MICs) for five antibiotics. Ultra-long read sequencing with Oxford Nanopore Technologies (ONT) and subsequent de novo assembly produced complete, gapless circular genomes for each strain. Alignment based annotation with Prokka identified 5,480 features in E. bruuniana and 5,203 features in E. meningoseptica, where none of these identified genes or gene combinations corresponded to observed phenotypic resistance values. Pan-genomic analysis, performed with an additional 19 Elizabethkingia strains, identified a core-genome size of 2,658,537 bp, 32 uniquely identifiable intrinsic chromosomal antibiotic resistance core-genes and 77 antibiotic resistance pan-genes. Using core-SNPs and pan-genes in combination with six machine learning (ML) algorithms, binary classification of clindamycin and vancomycin resistance achieved f1 scores of 0.94 and 0.84, respectively. Performance on the more challenging multiclass problem for fusidic acid, rifampin and ciprofloxacin resulted in f1 scores of 0.70, 0.75, and 0.54, respectively. By producing two sets of quality biological predictors, pan-genome genes and core-genome SNPs, from long-read sequence data and applying an ensemble of ML techniques, our results demonstrated that accurate phenotypic inference, at multiple AMR resolutions, can be achieved.

Classification of Mycobacterium tuberculosis DR, MDR,XDR Isolates and Identification of Signature MutationPattern of Drug Resistance.

Mycobacterium tuberculosis - a global threat, the recent breakout in MDR-TB and XDR-TB has challenged researchers in diagnosis to provide effective treatment. The main objective to combat drug resistance is to provide rapid, reliable and sensitive diagnostic methods in health care centres. This study focuses on development of an effective pipeline to identify drug resistance mutations in whole genome data of Mycobacterium tuberculosis utilizing the Next Generation Sequencing approach and classification of drug resistance strains based on genetic markers obtained from TGS-TB, tbvar and TBDReamDB. 74 isolates are characterized into 20 DR-TB, 16 MDR-TB, 16 XDR-TB and 6 nonresistant strains based on known drug resistance genetic markers. Results provide mutation pattern for each of the classified strains and profiling of drug resistance to the group of anti-TB drugs. The presence of specific mutation causing resistance to a drug will help set the dosage levels which play an important role in the treatment. Findings on amino acid changes and its respective codon positions in candidate genes will provide insights in drug sensitivity and a way for discovery of potent drugs. The implementation of these approaches in clinical setting provides rapid and sensitive diagnostics to combat the emerging drug resistance.

Characterization of antimicrobial resistance genes from Neisseria gonorrhoeae positive remnant Aptima urine specimens.

Aim: To ascertain the antimicrobial resistance and strain types (STs) of Neisseria gonorrhoeae from 50 remnant Aptima urine specimens using molecular methods. Methods: Mutations predictive of resistance to six antibiotics were identified in eight genes. STs were determined using NG-MAST and NG-STAR. Results: All eight antimicrobial resistance genes could be characterized in 36 specimens. A total of 17 specimens were predicted to be susceptible to all antibiotics, including ceftriaxone. Decreased susceptibility to cefixime and ciprofloxacin resistance was predicted in 11 specimens (PBP2 type 34.001). Overall, 38/50 specimens were predicted to be ciprofloxacin susceptible; three were azithromycin resistant. Nineteen NG-MAST and 21 NG-STAR STs were noted. Conclusion: Molecular analysis of remnant Aptima specimens enabled the prediction of emerging gonococcal cefixime and azithromycin resistance which would otherwise have been undetected.