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

Abstract

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.

Prevalence and trend analysis of antimicrobial resistance in clinical Escherichia coli isolates collected from diseased pigs in the USA between 2006 and 2016

“…For the network analysis approach, MIC data were analysed using Gaussian graphical models. This is a type of undirected graphical probabilistic model which has been used to estimate partial correlations between AMRs using MIC values directly and plotting them for easy graphical inferences, as previously described using 'Rnets' package in R version 3.5.3 (Love, Zawack, Booth, Grӧhn, & Lanzas, 2016). Briefly, MIC values for each year were first converted into empirical correlation matrices using Spearman's rank correlation method.…”

Monte Carlo Simulations Suggest Current Chlortetracycline Drug-Residue Based Withdrawal Periods Would Not Control Antimicrobial Resistance Dissemination from Feedlot to Slaughterhouse

“…Since our model was limited to modeling E. coli within the large intestine, we were unable to directly account for environmental factors that may impact resistance dissemination and persistence in beef feedlots. In addition, this model did not account for co-selection of CTC resistance as a result of administering other antimicrobials, which may be an important mechanism for resistance dissemination and persistence in food animals (Love et al, 2016 ). The model parameter value distributions were assigned based on the available published literature.…”

Establishing Statistical Equivalence of Data from Different Sampling Approaches for Assessment of Bacterial Phenotypic Antimicrobial Resistance

“…The most commonly used measurement of susceptibility of a bacterial isolate to an antimicrobial is the drug's MIC. When the MIC is measured using the broth microdilution assay based on serial 2-fold dilutions of the drug, the measurement is transformed to log 2 (MIC) for statistical analyses ( 12 , 15 ). The measurements for all the target bacterial species' isolates obtained via a given sampling or sample processing approach from the target stratum yield the antimicrobial's MIC distribution for the species in the stratum.…”