A hierarchical Bayesian latent class mixture model with censorship for detection of linear temporal changes in antibiotic resistance

Identifying and controlling the emergence of antimicrobial resistance (AMR) is a high priority for researchers and public health officials. One critical component of this control effort is timely detection of emerging or increasing resistance using surveillance programs. Currently, detection of temporal changes in AMR relies mainly on analysis of the proportion of resistant isolates based on the dichotomization of minimum inhibitory concentration (MIC) values. In our work, we developed a hierarchical Bayesian latent class mixture model that incorporates a linear trend for the mean log2MIC of the non-resistant population. By introducing latent variables, our model addressed the challenges associated with the AMR MIC values, compensating for the censored nature of the MIC observations as well as the mixed components indicated by the censored MIC distributions. Inclusion of linear regression with time as a covariate in the hierarchical structure allowed modelling of the linear creep of the mean log2MIC in the non-resistant population. The hierarchical Bayesian model was accurate and robust as assessed in simulation studies. The proposed approach was illustrated using Salmonella enterica I,4,[5],12:i:- treated with chloramphenicol and ceftiofur in human and veterinary samples, revealing some significant linearly increasing patterns from the applications. Implementation of our approach to the analysis of an AMR MIC dataset would provide surveillance programs with a more complete picture of the changes in AMR over years by exploring the patterns of the mean resistance level in the non-resistant population. Our model could therefore serve as a timely indicator of a need for antibiotic intervention before an outbreak of resistance, highlighting the relevance of this work for public health. Currently, however, due to extreme right censoring on the MIC data, this approach has limited utility for tracking changes in the resistant population.

[1]  Philippe Lambert,et al.  A Bayesian approach to the semiparametric estimation of a minimum inhibitory concentration distribution , 2016 .

[2]  Shaohua Zhao,et al.  Using machine learning to predict antimicrobial minimum inhibitory concentrations and associated genomic features for nontyphoidal Salmonella , 2018, bioRxiv.

[3]  Johan W Mouton,et al.  Breakpoints: current practice and future perspectives. , 2002, International journal of antimicrobial agents.

[4]  Geert Verbeke,et al.  A new semi-parametric mixture model for interval censored data, with applications in the field of antimicrobial resistance , 2014, Comput. Stat. Data Anal..

[5]  Marc Aerts,et al.  A trend analysis of antimicrobial resistance in commensal Escherichia coli from several livestock species in Belgium (2011-2014). , 2015, Preventive veterinary medicine.

[6]  Geert Verbeke,et al.  Application of the Vertex Exchange Method to estimate a semi-parametric mixture model for the MIC density of Escherichia coli isolates tested for susceptibility against ampicillin. , 2015, Biostatistics.

[7]  B A Craig.  Modeling approach to diameter breakpoint determination. , 2000, Diagnostic microbiology and infectious disease.

[8]  Shaohua Zhao,et al.  Antimicrobial Drug Resistance in Escherichia coli from Humans and Food Animals, United States, 1950–2002 , 2012, Emerging infectious diseases.

[9]  L. Warnick,et al.  Antimicrobial resistance trends among Salmonella isolates obtained from horses in the northeastern United States (2001-2013). , 2016, American journal of veterinary research.

[10]  Reza Mazloom,et al.  Potential Information Loss Due to Categorization of Minimum Inhibitory Concentration Frequency Distributions. , 2018, Foodborne pathogens and disease.

[11]  G Kahlmeter,et al.  Escherichia coli and Staphylococcus aureus: bad news and good news from the European Antimicrobial Resistance Surveillance Network (EARS-Net, formerly EARSS), 2002 to 2009. , 2011, Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin.

[12]  B. Mallick,et al.  Generalized Linear Models : A Bayesian Perspective , 2000 .

[13]  M A Hamilton,et al.  Descriptive statistical analyses of serial dilution data. , 1988, Statistics in medicine.

[14]  Arnošt Komárek,et al.  Bayesian estimation of multivariate normal mixtures with covariate‐dependent mixing weights, with an application in antimicrobial resistance monitoring , 2018, Biometrical journal. Biometrische Zeitschrift.

[15]  Jan van de Kassteele,et al.  New Statistical Technique for Analyzing MIC-Based Susceptibility Data , 2012, Antimicrobial Agents and Chemotherapy.

[16]  L. Dutil,et al.  Canadian Integrated Program for Antimicrobial Resistance Surveillance (CIPARS) Farm Program: Results from Finisher Pig Surveillance , 2010, Zoonoses and public health.

[17]  Arnost Komárek,et al.  A new R package for Bayesian estimation of multivariate normal mixtures allowing for selection of the number of components and interval-censored data , 2009, Comput. Stat. Data Anal..

[18]  Christel Faes,et al.  Development of statistical methods for the evaluation of data on antimicrobial resistance in bacterial isolates from animals and food , 2011 .

[19]  Bruce A Craig,et al.  Statistical properties and inference of the antimicrobial MIC test , 2005, Statistics in medicine.

[20]  J. Ruiz,et al.  From MIC creep to MIC decline: Staphylococcus aureus antibiotic susceptibility evolution over the last 4 years. , 2016, Clinical microbiology and infection : the official publication of the European Society of Clinical Microbiology and Infectious Diseases.