Abstract
OBJECTIVE: Pseudomonas aeruginosa has strong drug resistance and can tolerate a variety of antibiotics, which is a major problem in the management of antibiotic-resistant infections. Direct prediction of multi-drug resistance (MDR) resistance phenotypes of P. aeruginosa isolates and clinical samples by genotype is helpful for timely antibiotic treatment. METHODS: In the study, whole genome sequencing (WGS) data of 494 P. aeruginosa isolates were used to screen key anti-microbial resistance (AMR)-associated genes related to imipenem (IPM), meropenem (MEM), piperacillin/tazobactam (TZP), and levofloxacin (LVFX) resistance in P. aeruginosa by comparing genes with copy number differences between resistance and sensitive strains. Subsequently, for the direct prediction of the resistance of P. aeruginosa to four antibiotics by the AMR-associated features screened, we collected 74 P. aeruginosa positive sputum samples to sequence by metagenomics next-generation sequencing (mNGS), of which 1 sample with low quality was eliminated. Then, we constructed the resistance prediction model. RESULTS: We identified 93, 88, 80, 140 AMR-associated features for IPM, MEM, TZP, and LVFX resistance in P. aeruginosa. The relative abundance of AMR-associated genes was obtained by matching mNGS and WGS data. The top 20 features with importance degree for IPM, MEM, TZP, and LVFX resistance were used to model, respectively. Then, we used the random forest algorithm to construct resistance prediction models of P. aeruginosa, in which the areas under the curves of the IPM, MEM, TZP, and LVFX resistance prediction models were all greater than 0.8, suggesting these resistance prediction models had good performance. CONCLUSION: In summary, mNGS can predict the resistance of P. aeruginosa by directly detecting AMR-associated genes, which provides a reference for rapid clinical detection of drug resistance of pathogenic bacteria.