Markov Networks of Collateral Resistance: National Antimicrobial Resistance Monitoring System Surveillance Results from Escherichia coli Isolates, 2004-2012

Surveillance of antimicrobial resistance (AMR) is an important component of public health. Antimicrobial drug use generates selective pressure that may lead to resistance against to the administered drug, and may also select for collateral resistances to other drugs. Analysis of AMR surveillance data has focused on resistance to individual drugs but joint distributions of resistance in bacterial populations are infrequently analyzed and reported. New methods are needed to characterize and communicate joint resistance distributions. Markov networks are a class of graphical models that define connections, or edges, between pairs of variables with non-zero partial correlations and are used here to describe AMR resistance relationships. The graphical least absolute shrinkage and selection operator is used to estimate sparse Markov networks from AMR surveillance data. The method is demonstrated using a subset of Escherichia coli isolates collected by the National Antimicrobial Resistance Monitoring System between 2004 and 2012 which included AMR results for 16 drugs from 14418 isolates. Of the 119 possible unique edges, 33 unique edges were identified at least once during the study period and graphical density ranged from 16.2% to 24.8%. Two frequent dense subgraphs were noted, one containing the five β-lactam drugs and the other containing both sulfonamides, three aminoglycosides, and tetracycline. Density did not appear to change over time (p = 0.71). Unweighted modularity did not appear to change over time (p = 0.18), but a significant decreasing trend was noted in the modularity of the weighted networks (p < 0.005) indicating relationships between drugs of different classes tended to increase in strength and frequency over time compared to relationships between drugs of the same class. The current method provides a novel method to study the joint resistance distribution, but additional work is required to unite the underlying biological and genetic characteristics of the isolates with the current results derived from phenotypic data.

[1]  R. Cantón,et al.  Co-resistance: an opportunity for the bacteria and resistance genes. , 2011, Current opinion in pharmacology.

[2]  M. Salman,et al.  Factor analysis of minimum-inhibitory concentrations for Escherichia coli isolated from feedlot cattle to model relationships among antimicrobial-resistance outcomes. , 2003, Preventive veterinary medicine.

[3]  T. Nakashima,et al.  Effects of environmental novelty on fear-related behavior and stress responses of rats to emotionally relevant odors , 2009, Behavioural Brain Research.

[4]  D. Hooper,et al.  Quinolone resistance locus nfxD of Escherichia coli is a mutant allele of the parE gene encoding a subunit of topoisomerase IV , 1997, Antimicrobial agents and chemotherapy.

[5]  A. Khodursky,et al.  Topoisomerase IV is a target of quinolones in Escherichia coli. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

[6]  R. Tibshirani,et al.  Sparse inverse covariance estimation with the graphical lasso. , 2008, Biostatistics.

[7]  Mark Newman,et al.  Networks: An Introduction , 2010 .

[8]  N. Meinshausen,et al.  High-dimensional graphs and variable selection with the Lasso , 2006, math/0608017.

[9]  J. Laurie Snell,et al.  Markov Random Fields and Their Applications , 1980 .

[10]  T. Tsutsui,et al.  Unraveling Antimicrobial Resistance Genes and Phenotype Patterns among Enterococcus faecalis Isolated from Retail Chicken Products in Japan , 2015, PloS one.

[11]  M. Ferraro.  Performance standards for antimicrobial susceptibility testing , 2001 .

[12]  Morten Otto Alexander.  Use of collateral sensitivity networks to design drug cycling protocols that avoid resistance development. , 2016 .

[13]  W. Velicer,et al.  Relation of sample size to the stability of component patterns. , 1988, Psychological bulletin.

[14]  F. Harary,et al.  STRUCTURAL BALANCE: A GENERALIZATION OF HEIDER'S THEORY1 , 1977 .

[15]  R. Jennrich,et al.  Quartic rotation criteria and algorithms , 1988 .

[16]  Larry A. Wasserman,et al.  High Dimensional Semiparametric Gaussian Copula Graphical Models. , 2012, ICML 2012.

[17]  M E J Newman,et al.  Finding and evaluating community structure in networks. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[18]  P. Hawkey,et al.  The changing epidemiology of resistance. , 2009, The Journal of antimicrobial chemotherapy.

[19]  Antoinette Ludwig,et al.  Identifying associations in Escherichia coli antimicrobial resistance patterns using additive Bayesian networks. , 2013, Preventive veterinary medicine.

[20]  G. Jacoby,et al.  A functional classification scheme for beta-lactamases and its correlation with molecular structure , 1995, Antimicrobial agents and chemotherapy.

[21]  G. Bloemberg,et al.  Consequences of revised CLSI and EUCAST guidelines for antibiotic susceptibility patterns of ESBL- and AmpC β-lactamase-producing clinical Enterobacteriaceae isolates. , 2013, The Journal of antimicrobial chemotherapy.

[22]  Ronald N. Jones,et al.  Background and rationale for revised clinical and laboratory standards institute interpretive criteria (Breakpoints) for Enterobacteriaceae and Pseudomonas aeruginosa: I. Cephalosporins and Aztreonam. , 2013, Clinical infectious diseases : an official publication of the Infectious Diseases Society of America.

[23]  A. Wyke,et al.  Cloning and characterization of a DNA gyrase A gene from Escherichia coli that confers clinical resistance to 4-quinolones , 1989, Antimicrobial Agents and Chemotherapy.

[24]  R. Cattell.  The Scree Test For The Number Of Factors. , 1966, Multivariate behavioral research.

[25]  Trevor Hastie,et al.  Statistical Learning with Sparsity: The Lasso and Generalizations , 2015 .

[26]  S. Hopkins,et al.  Antimicrobial stewardship: English Surveillance Programme for Antimicrobial Utilization and Resistance (ESPAUR). , 2013, The Journal of antimicrobial chemotherapy.

[27]  J. Martínez,et al.  Contribution of a New Mutation in parE to Quinolone Resistance in Extended-Spectrum-β-Lactamase-Producing Escherichia coli Isolates , 2007, Journal of Clinical Microbiology.

[28]  Jianzhong Shen,et al.  Emergence of plasmid-mediated colistin resistance mechanism MCR-1 in animals and human beings in China: a microbiological and molecular biological study. , 2015, The Lancet. Infectious diseases.

[29]  J. Yamagishi,et al.  Nalidixic acid-resistant mutations of the gyrB gene of Escherichia coli , 1986, Molecular and General Genetics MGG.

[30]  Walker Rd,et al.  Standards for antimicrobial susceptibility testing. , 1999 .

[31]  M. Valcárcel,et al.  Unsegmented flow approach for on-line monitoring of pH, conductivity, dissolved oxygen and determination of nitrite and ammonia in aquaculture , 1994, The Journal of automatic chemistry.

[32]  A. M. George,et al.  Amplifiable resistance to tetracycline, chloramphenicol, and other antibiotics in Escherichia coli: involvement of a non-plasmid-determined efflux of tetracycline , 1983, Journal of bacteriology.

[33]  J. Franklin,et al.  The elements of statistical learning: data mining, inference and prediction , 2005 .

[34]  D. Paterson,et al.  Colistin resistance: a major breach in our last line of defence. , 2016, The Lancet. Infectious diseases.

[35]  Robert J. Clifford,et al.  Escherichia coli Harboring mcr-1 and blaCTX-M on a Novel IncF Plasmid: First Report of mcr-1 in the United States , 2016, Antimicrobial Agents and Chemotherapy.

[36]  P. M. Terry,et al.  Rapid development of ciprofloxacin resistance in methicillin-susceptible and -resistant Staphylococcus aureus. , 1991, The Journal of infectious diseases.

[37]  D. Hooper,et al.  Cross-resistance to fluoroquinolones in multiple-antibiotic-resistant (Mar) Escherichia coli selected by tetracycline or chloramphenicol: decreased drug accumulation associated with membrane changes in addition to OmpF reduction , 1989, Antimicrobial Agents and Chemotherapy.

[38]  H. Kaiser.  The Application of Electronic Computers to Factor Analysis , 1960 .

[39]  B. Spratt.  Properties of the penicillin-binding proteins of Escherichia coli K12,. , 1977, European journal of biochemistry.

[40]  Gábor Csárdi,et al.  The igraph software package for complex network research , 2006 .

[41]  Pablo Jensen,et al.  Analysis of community structure in networks of correlated data. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[42]  C. Spearman.  The proof and measurement of association between two things. , 2015, International journal of epidemiology.

[43]  D. Altman,et al.  Analysis by Categorizing or Dichotomizing Continuous Variables Is Inadvisable: An Example from the Natural History of Unruptured Aneurysms , 2011, American Journal of Neuroradiology.

[44]  R. Kishony,et al.  Multidrug evolutionary strategies to reverse antibiotic resistance , 2016, Science.

[45]  P. Phillips,et al.  Comparative quantitative genetics : evolution of the G matrix , 2002 .

[46]  F. Baquero.  From pieces to patterns: evolutionary engineering in bacterial pathogens , 2004, Nature Reviews Microbiology.

[47]  Pei Wang,et al.  Partial Correlation Estimation by Joint Sparse Regression Models , 2008, Journal of the American Statistical Association.

[48]  Valerii Fedorov,et al.  Consequences of dichotomization , 2009, Pharmaceutical statistics.