Next Article in Journal
Co-Fermentation of Edible Mushroom By-Products with Soybeans Enhances Nutritional Values, Isoflavone Aglycones, and Antioxidant Capacity of Douchi Koji
Next Article in Special Issue
Accurate Detection of Salmonella Based on Microfluidic Chip to Avoid Aerosol Contamination
Previous Article in Journal
High Levels of Policosanols and Phytosterols from Sugar Mill Waste by Subcritical Liquefied Dimethyl Ether
Previous Article in Special Issue
Use of Large-Scale Genomics to Identify the Role of Animals and Foods as Potential Sources of Extraintestinal Pathogenic Escherichia coli That Cause Human Illness
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantitative Risk Assessment of Susceptible and Ciprofloxacin-Resistant Salmonella from Retail Pork in Chiang Mai Province in Northern Thailand

by
Chaiwat Pulsrikarn
1,
Anusak Kedsin
2,
Parichart Boueroy
2,
Peechanika Chopjitt
2,
Rujirat Hatrongjit
3,
Piyarat Chansiripornchai
4,
Nipattra Suanpairintr
4,5 and
Suphachai Nuanualsuwan
5,6,*
1
National Institute of Health, Department of Medical Science, Ministry of Public Health, Nonthaburi 11000, Thailand
2
Faculty of Public Health, Kasetsart University, Chalermphrakiat Sakon Nakhon Province Campus, Sakon Nakhon 47000, Thailand
3
Faculty of Science and Engineering, Chalermphrakiat Sakon Nakhon Province Campus, Kasetsart University, Sakon Nakhon 47000, Thailand
4
Department of Pharmacology, Faculty of Veterinary Science, Chulalongkorn University, Bangkok 10330, Thailand
5
Center of Excellence for Food and Water Risk Analysis (FAWRA), Department of Veterinary Public Health, Faculty of Veterinary Science, Chulalongkorn University, Bangkok 10330, Thailand
6
Department of Veterinary Public Health, Faculty of Veterinary Sciences, Chulalongkorn University, Bangkok 10330, Thailand
*
Author to whom correspondence should be addressed.
Foods 2022, 11(19), 2942; https://doi.org/10.3390/foods11192942
Submission received: 14 July 2022 / Revised: 26 August 2022 / Accepted: 16 September 2022 / Published: 20 September 2022

Abstract

:
The adverse human health effects as a result of antimicrobial resistance have been recognized worldwide. Salmonella is a leading cause of foodborne illnesses while antimicrobial resistant (AMR) Salmonella has been isolated from foods of animal origin. The quantitative risk assessment (RA) as part of the guidelines for the risk analysis of foodborne antimicrobial resistance was issued by the Codex Alimentarius Commission more than a decade ago. However, only two risk assessments reported the human health effects of AMR Salmonella in dry-cured pork sausage and pork mince. Therefore, the objective of this study was to quantitatively evaluate the adverse health effects attributable to consuming retail pork contaminated with Salmonella using risk assessment models. The sampling frame covered pork at the fresh market (n = 100) and modern trade where pork is refrigerated (n = 50) in Chiang Mai province in northern Thailand. The predictive microbiology models were used in the steps where data were lacking. Susceptible and quinolone-resistant (QR) Salmonella were determined by antimicrobial susceptibility testing and the presence of AMR genes. The probability of mortality conditional to foodborne illness by susceptible Salmonella was modeled as the hazard characterization of susceptible and QR Salmonella. For QR Salmonella, the probabilistic prevalences from the fresh market and modern trade were 28.4 and 1.9%, respectively; the mean concentrations from the fresh market and modern trade were 346 and 0.02 colony forming units/g, respectively. The probability of illness (PI) and probability of mortality given illness (PMI) from QR Salmonella-contaminated pork at retails in Chiang Mai province were in the range of 2.2 × 10−8–3.1 × 10−4 and 3.9 × 10−10–5.4 × 10−6, respectively, while those from susceptible Salmonella contaminated-pork at retails were in the range 1.8 × 10−4–3.2 × 10−4 and 2.3 × 10−7–4.2 × 10−7, respectively. After 1000 iterations of Monte Carlo simulations of the risk assessment models, the annual mortality rates for QR salmonellosis simulated by the risk assessment models were in the range of 0–32, which is in line with the AMR adverse health effects previously reported. Therefore, the risk assessment models used in both exposure assessment and hazard characterization were applicable to evaluate the adverse health effects of AMR Salmonella spp. in Thailand.

1. Introduction

The adverse health effects posed by antimicrobial-resistant (AMR) bacteria have been increasing at an alarming rate and have recently been recognized worldwide [1]. Antimicrobial agents are beneficial to a wide variety of sectors, from human and veterinary medicine to animal and plant production. The use of antimicrobial agents in these areas is inevitably connected and this renders a circulating pool of both resistant bacteria and bacteria-borne resistant genes that are eventually delivered to humans [2]. Regardless of the environmental, genetic, or spatial boundary, mobile genetic elements containing resistance determinants can, directly and indirectly, propagate through horizontal transfer among bacteria from foods of animal origin and their environment to humans. Therefore, AMR risk management measures in terms of prevention and control strategy rely heavily on source attribution and risk assessment to evaluate the likelihood and severity of the consequences of AMR bacteria-contaminated foods [3].
The Codex Alimentarius Commission (CAC) has endorsed a systematic framework for foodborne AMR risk analysis. This framework is composed of preliminary risk management activities, risk assessment, and risk management. These three components are connected via risk communication, including the surveillance of AMR and other sources of information. The underlying rationale of the principle of risk analysis is to evaluate the risk to human health from foodborne AMR microorganisms and AMR determinants so that practical risk management measures can be implemented to prevent and control such human health risks [4].
The microbial risk assessment is a scientific process to evaluate the risk of consuming food contaminated with hazards. Hazard identification, hazard characterization, exposure assessment, and risk characterization constitute the four-step risk assessment. Hazard identification is the initial step to examining the risk of the hazards such as foodborne disease viruses, bacteria, protozoa, and parasites; in this study, the hazard is Salmonella. Hazard characterization determines the probability of illness upon getting a hazard into the host by a specific dose–response model. The exposure assessment determines the probability of getting hazards through consuming food. The last step is risk characterization, where the risk estimate is derived from the product of probabilities of exposure and illness from the preceding two steps [5,6]. Foodborne AMR risk assessment (AMR RA) is slightly different from the traditional methodology of microbiological risk assessment [5,6] in that hazard characterization is necessary to additionally include the adverse effects of AMR, e.g., antimicrobial treatment failure, prolonged treatment period, more illness severity or virulence, and higher mortality rate [4].
The prevalence of susceptible Salmonella spp. from swine manure was in the range of 2–61% and from swine farm swabs it was 95%, whereas that of AMR Salmonella spp. isolated from antimicrobial-use swine farms was lowest at 33% against florfenicol and highest at 66% against tetracycline [7]. Likewise, the prevalence of tetracycline-resistant Salmonella spp. was even higher at 90% in two independent studies [8,9]. However, isolated Salmonella spp. was sensitive to ceftiofur, ceftriaxone, and ciprofloxacin. In addition, fluoroquinolone-resistant Salmonella spp. warrant further surveillance by the World Health Organization as a tier group 2 [10].
AMR Salmonella spp. in retail pork could be derived either from a farm or abattoir. The prevalence of AMR Salmonella spp. from the environment of the abattoir was lowest at 4% against ceftiofur and highest at 86–89% against tetracycline [7,8]. Recently, we investigated a total of 387 non-typhoidal Salmonella enterica (NTS) isolated from abattoirs. Approximately 24% of NTS isolates were AMR, while only 6% of NTS isolates were susceptible to all antimicrobial agents tested. However, non-AMR NTS isolates carry extended-spectrum beta-lactamase (blaCTX-M) genes or narrow-spectrum beta-lactamase genes (blaTEM or blaSHV). The rest of the NTS isolates (70%) were susceptible to all fluoroquinolones as well as carbapenems and third-generation cephalosporins [11]. At retail, Salmonella spp. isolated from pork was susceptible to ampicillin, norfloxacin, and ciprofloxacin [8,12]. However, the prevalence of AMR Salmonella spp. isolated from retail pork were 100% against streptomycin and sulfamethoxazole [12] and 60% against tetracycline [8]. The source attribution of salmonellosis in children from pork was 11% [8].
Scientific evidence demonstrated that AMR Salmonella spp. is foodborne-transmitted. On the other hand, AMR Salmonella infection is seldom traced from patients in hospitals back through contaminated foods and even further back to animals along the food chain. Some of the implicated commodities in such reports were beef, pork, and milk, where the authors suggested that AMR Salmonella spp. in patients was attributable to farm animals [13,14,15]. While commonly found Salmonella serovars with either resistance or multiple resistance to Salmonella spp. in foods are Derby, Enteritidis, Hadar, Newport, Paratyphi, Typhimurium, and Virchow [16,17,18], Salmonella Typhimurium is the most prevalent serovar contaminating foods across continents [13,14,15,17,19,20]. Recently, both cephalosporin-resistant and extended-spectrum beta-lactamase-resistant Salmonella spp. have been frequently reported [21,22]. These reports implied that AMR Salmonella spp. has been widely circulated regardless of geographical borders, food commodities, serovars, resistance patterns, antimicrobial classes, and host ranges.
Even though 34 AMR RA relating to retail foods have been reported up to 2018, only eight articles investigated the adverse health effects of AMR Salmonella spp. Only half of these reports are related to the pork supply chain [23]. Two risk assessments reported the adverse health effects of AMR Salmonella in dry-cured pork sausage and pork mince [24,25]. Recently, a farm-to-fork quantitative risk assessment of Salmonella Heidelberg resistant to third-generation cephalosporins in broiler chickens was reported [26] while the AMR RA model was developed for anti-E. coli drugs [27]. However, a quantitative risk assessment using Monte Carlo simulation of QR Salmonella in retail pork has never been reported. In this study, QR Salmonella-contaminated pork was the hazard of interest. The sampling frame covered pork at retailers in Chiang Mai province in northern Thailand. The predictive microbiology models were used in the steps where data were lacking. The objective of this study was to comparatively evaluate the adverse health effects attributable to consuming pork contaminated with Salmonella susceptible and resistant to quinolone.

2. Materials and Methods

2.1. Pork Samples

The pork samples were collected in Chiang Mai province from both the fresh market and modern trade where pork is refrigerated. Ten pork samples were collected from each retailer. Eleven pork retailers from the fresh market and five butcher shops in the modern trade participated in this study. The sampling unit of pork was at least 100 g. Samples were collected using an aseptic technique to avoid undesirable cross-contamination from environmental fomite and then kept in a leak-proof container between 2 and 8 °C during transportation. The samples arrived at the laboratory and were analyzed within 8–10 h after being collected.

2.2. Enumeration of Salmonella

The ten-fold serial dilution of pork samples was achieved using buffered peptone water. For individual dilution, 1 mL of suspension was repeatedly transferred 3 times into 3 separate 9 mL tubes of Rappaport Vassiliades with soya (RVS) broth. Nine tubes of RVS broth for each sample were incubated at 42 °C for 24 h. Only RVS tubes with a turbid appearance and confirmed by xylose lysine desoxycholate agar and then triple sugar iron slant were counted as positive [28]. The concentrations of Salmonella in the Most Probable Number unit (MPN) were converted to colony-forming units (cfu) by multiplying by 0.8 since the MPN technique is more sensitive than a standard plate count by 25% [29]. The unit conversion of concentration is necessary to apply for a dose–response model using the dose unit as the cfu [30].

2.3. Antimicrobial Susceptibility Testing (AST)

Susceptibility testing for ampicillin, cefepime, cefotaxime, cefoxitin, chloramphenicol, ciprofloxacin, colistin, gentamicin, imipenem, meropenem, nalidixic acid, streptomycin, sulphamethoxazole, tetracycline, and trimethoprim was performed using a broth microdilution assay to determine the minimum inhibitory concentration (MIC) according to the M07 Clinical and Laboratory Standards Institute (CLSI) guidelines [31]. The results were interpreted according to the 2020 Clinical and Laboratory Standards Institute guidelines for the susceptibility testing of Salmonella isolates [32]. Escherichia coli ATCC 25922 was used as the control. The broth microdilution assay was performed using two-fold dilution at a concentration in a range of 0.03–64 μg/mL depending on the antimicrobial agents, which are suggested based on the 2020 CLSI.

2.4. Determination of Antimicrobial Resistance Genes

Antimicrobial resistance genes including quinolone, colistin, and carbapenem were conducted using a polymerase chain reaction (PCR). The PCR was carried out to determine the quinolone resistance determining region of gyrA and parC, and the plasmid-mediated quinolone resistance genes following are described elsewhere [33,34]. The PCR products of the quinolone resistance determining the region from the four genes were purified and subjected to Sanger sequencing (performed by Apical Scientific Sdn Bhd, Selangor, Malaysia) to determine their substitution by comparing with those of wild-type S. Typhimurium LT2 [11]. The presence of antibiotic resistance-conferring genes of colistin, including mcr-1 through mcr-9, and carbapenem consisting of blaNDM, blaOXA-48-like, blaIMP, and blaKPC was investigated using the PCR method described elsewhere [11]. All PCRs performed in this study are described in the Supplementary Materials (Tables S1–S4).

2.5. Risk Assessment Models

2.5.1. Exposure Assessment

  • Probabilistic prevalence variable
The range of prevalence is between zero (0%) and one (100%), inclusively applicable to the range of Beta distribution. The Beta distribution is characterized by 2 parameters, alpha and beta, as shown in Equation (1).
PPROB = Beta (α, β)
To describe the variability of prevalence, the alpha parameter is substituted by s + α, and the beta parameter is substituted by ns + β where s is the number of the successful trial (s) in the identical n trials of a binomial process, as shown in Equation (2). In this study, the successful trials were the QR Salmonella-contaminated (positive) samples where the identical n trials were the sample size.
PPROB = Beta (s + α, ns + β)
This study assumes that no prior prevalence of QR Salmonella was reported. The uniform probability distribution was assumed, which is equivalent to Beta (1, 1). Therefore, two parameters in Equation (2) were replaced with 1, as shown in Equation (3) [6].
PPROB = Beta (s + 1, ns + 1)
2.
Thermal inactivation model
The raw pork from retail was subjected to heat treatment before consumption. The cooking temperature and time were 64 °C for 2 min while the decimal reduction time at 64 °C (D64) is 0.48 min [30]. The log reduction of Salmonella is shown in Equation (4).
L R = t D 64
where LR = log reduction (LR) of susceptible or QR Salmonella in pork; D64 = decimal reduction time of Salmonella at 64 °C (min); t = cooking time (min).
3.
Concentration variable
If pork samples were all negative, the Salmonella concentration was determined by the maximum likelihood estimator (MLE) technique [29,35,36], as shown in Equation (5).
log   reduction   LR   C S = i = 1 k N i i = 1 k V i 10 L R
where CS = concentration of susceptible or QR Salmonella (g−1); Ni = no. of Salmonella detected in retail pork i to k; Vi = analytical unit of pork i to k (g); k = no. of pork retailers; LR = log reduction of Salmonella from heat treatment.
4.
Consumption variable (CP)
Food consumption data for Thailand in 2016 from the Agricultural Commodity and Food Standard report showed that the mean and 97.5th percentile consumption of pork among eaters more than 3 years old was 14.12 and 58.28 g/person/day, respectively. The triangular distribution was used to describe the variability of the consumption variable. The three parameters of triangular distribution (minimum, most likely, and maximum) were 0, 14.12, and 58.28 g/person/day, respectively.
5.
Dose of Salmonella ingested
The dose of Salmonella ingested was the product of Salmonella concentration after cooking and pork consumption per day. The equation for the dose of Salmonella ingested is shown in Equation (6) [6].
D = CS × CP
where D = dose of susceptible or QR Salmonella ingested per day (cfu); CS = concentration of susceptible or QR Salmonella (cfu/g); CP = pork consumption per day (g).
6.
Probability of exposure (PE)
PE is the likelihood of experiencing at least one cell of Salmonella from pork. Therefore, the input variables to model the PE are the concentration (CS) and prevalence (PPROB) of Salmonella, including pork consumption (6), as shown in Equation (7).
PE = PPROB (1 − exp − D)

2.5.2. Hazard Characterization

  • Probability of illness (PI)
The dose–response model was used to characterize the probability of illness caused by either residual susceptible or QR Salmonella-contaminated pork after cooking, as shown in Equation (8).
PI = 1 − (1+ (D/51.45))−0.1324
where PI = the probability of illness caused by an ingested dose of Salmonella; D = dose of susceptible or QR Salmonella ingested per day (cfu).
2.
Probability of mortality (PM)
Additional to the conventional hazard characterization of the microbial risk assessment, the adverse effects of AMR such as a higher mortality rate were included [4]. A previous study reported that the mortality rates caused by drug-susceptible and multidrug-resistant non-typhoid Salmonella were 0.2 and 3.4%, respectively [13]. Likewise, another study reported that the mortality rates caused by pan-susceptible and AMR Salmonella were 0.06 and 0.1%, respectively [18]. Therefore, in this study, the mean mortality rates as PM caused by susceptible and AMR Salmonella were averaged from these two previous reports as 0.13 and 1.75%, respectively.
3.
Probability of mortality given illness (PMI)
The integration of adverse health effects as the mortality conditional to the foodborne illness is the product of PI and PM, as shown in Equation (9).
PMI = PM × PI

2.5.3. Risk Characterization

In this study, the risk characterization is a two-step linked process of exposure assessment and hazard characterization. The probability of mortality given illness (PMI) is conditional on PE. Assuming that adverse health effects and hazard exposure are independent, the model for risk estimates in terms of the probability of foodborne mortality (PFM) is the product of PMI and PE, as shown in Equation (10).
PFM = PMI × PE
The probability of foodborne mortality from at least one day was calculated based on the binomial theorem [36]. The number of annual foodborne mortality cases per 100,000 population is calculated from Equation (11).
MAFM = (1 − (1 − PFM)365) × 100,000
where MAFM = annual foodborne mortality cases per 100,000 population; PFM = probability of foodborne mortality per day.
Simulations of MAFM were run for 10,000 iterations. The Simulación 4.0 freeware (developed by José Ricardo Varela) was used to run the Monte Carlo simulations.

2.6. Statistical Analysis

The MAFM of susceptible and QR Salmonella in pork from the fresh market and modern trade was determined for the statistical difference by one-way analysis of variance (ANOVA) [37]. Tukey’s multiple comparison test was followed to determine the pair-wise differences of MAFM. The IBM® SPSS® Statistics version 22 software (SPSS Inc., Chicago, IL, USA) was used to perform statistical analyses.

3. Results

3.1. Exposure Assessment

A total of 150 pork samples collected from pork retailers (fresh market (n = 100) and modern trade (n = 50)) in Chiang Mai province were analyzed for Salmonella contamination. The number of Salmonella-positive samples is shown in Table 1. All Salmonella isolates from positive samples were subject to the AST. We determined antimicrobial-resistant genes in QR isolates for colistin (mcr-1 through mcr-9), carbapenem, and fluoroquinolone including mcr, blaNDM, blaOXA-48-like, blaIMP, blaKPC, plasmid-mediated quinolone resistance, and the quinolone resistance-determining region of gyrA and parC. No isolates carried the mobile colistin resistance gene (mcr) and common carbapenemase genes (blaNDM, blaOXA-48-like, blaIMP, blaKPC). In the case of fluoroquinolone-resistant genes, among the QR isolates, five isolates carried qnrS, there were two substitutions in parC, and one isolate carried both qnrS and parC substitutions. No substitution occurred in gyrA in all isolates. PPROB and mean concentrations corresponding to susceptible and QR Salmonella contaminated in the pork samples are shown in Table 2. The PE to susceptible and QR Salmonella-contaminated pork at retail in Chiang Mai province was in the range of 2 × 10−7–0.03 (Table 3).

3.2. Hazard Characterization

The PI and PMI from QR Salmonella-contaminated pork at retails in Chiang Mai province were in the range of 2.2 × 10−8–3.1 × 10−4 and 3.9 ×10−10–5.4 × 10−6, respectively, while those from susceptible Salmonella-contaminated pork at retails were in the range of 1.8 × 10−4–3.2 ×10−4 and 2.3 × 10−7–4.2 × 10−7, respectively (Table 3).

3.3. Risk Characterization

The descriptive statistics and probability distributions of risk estimates in terms of PFM and MAFM from consuming retail pork contaminated with susceptible and QR Salmonella in Chiang Mai province, after performing a Monte Carlo simulation, are shown in Table 4 and Figure 1, Figure 2 and Figure 3. The mean PFM of susceptible Salmonella was lower than that of QR Salmonella from the fresh market. On the other hand, in the modern trade, the mean PFM of susceptible Salmonella became higher than that of QR Salmonella, essentially because the mean concentration of susceptible Salmonella was much higher than that of QR Salmonella.

4. Discussion

Two major approaches to AMR RA were determined by the data characteristics. The qualitative approach requires only a few calculations. The data variable is measured by the ordinal scale, e.g., low, moderate, and high. This could avoid complicated mathematical models and statistics, thus rendering risk assessment more straightforward, prolific, and time-saving. Nevertheless, the major drawback of qualitative AMR RA is the inherent subjectivity. One recommended solution to this dilemma is to transparently state or match the numerical values corresponding to individual descriptive terms for a qualitative variable [38,39]. Even though CAC encourages the quantitative technique to be performed as much as possible, the qualitative technique could not be discounted [4]. For the “quantitative technique”, the variables are measured by either interval or ratio scale. Two subcategories of quantitative AMR RA are deterministic and stochastic methods. Variables in the deterministic method possess only one single value, while those in the stochastic method encompass probability density corresponding to all possible values of a variable in the form of probability distribution [40,41,42]. This technique is more objective than the former technique, while complicated mathematical models are involved in almost every step of AMR RA (from hazard characterization to risk characterization) since the data in this study were allowed to quantitatively evaluate the mortality risk using Monte Carlo simulations. Therefore, the outputs from the mathematical models such as PE, PMI, and risk estimate are comparable whether between susceptible and QR Salmonella or fresh market and modern trade.
To better quantify the risk of exposure to the hazard, the types of hazard should be defined. Hazard, in the context of AMR RA, is either AMR pathogenic bacteria or an AMR determinant. The former hazard or sometimes so-called direct hazard in food is the AMR pathogenic microorganism being capable of colonizing and then infecting a human host. Furthermore, the direct hazard is also derived from handling contaminated food [43], while AMR bacteria harboring resistance genes directly transfer resistance genes to pathogenic bacteria or indirectly transfer to the commensal bacteria. The AMR determinant or resistant genes transferred through the last two mechanisms is a so-called indirect hazard [4]. This study determines the AMR hazard by both phenotypic and genotypic analyses; therefore, the AMR PPROB is more conservative and prevalent than taking into account only the AMR hazard from the genotypic analysis [44].
This study collected pork samples in Chiang Mai province in northern Thailand to investigate the risk of consuming pork contaminated with susceptible and QR Salmonella. The PPROB of susceptible and QR Salmonella isolated from the fresh market were in the narrow range of 28–30% (Table 2), while the PPROB of susceptible Salmonella was about 10 times higher than the PPROB of QR Salmonella isolated from the modern trade. The overall PPROB of (both susceptible and QR) Salmonella from the fresh market is eight times more than the PPROB of Salmonella from the modern trade. Likewise, QR Salmonella from the fresh market is almost 15 times more prevalent than susceptible Salmonella from the modern trade. In 2014, a similar study collected pork samples to compare the prevalence of susceptible and AMR Salmonella from the fresh market and the modern trade in Chiang Mai [45]. Even though 73% of fresh-market pork contaminated with Salmonella was more prevalent than only 10% of modern-trade pork contaminated with Salmonella, Salmonella prevalence from the fresh market in this previous study was slightly higher than the PPROB of Salmonella from the fresh market in our study. These compatible findings suggest that the sanitation along the pork supply chain of the fresh market in Chiang Mai province should have been improved.
Even though several Salmonella contaminations along the pork supply chain from farms and slaughterhouses to retail were reported in Chiang Mai province in northern Thailand [11,46,47,48,49,50], the magnitude of the contamination of Salmonella was reported as a percentage by the detection technique, since the risk assessment approach recommended by the Codex Alimentarius requires both the prevalence and concentration of Salmonella, particularly at the point of consumption. Only one previous study in Chiang Mai reported that Salmonella prevalence and concentration in pork from the fresh market were 39% (27/70) and 1.31 ± 0.25 log MPN/g, respectively [46]. The mean concentration of Salmonella from the previous study was lower than that of Salmonella from the fresh market in our study at 1.8 ± 0.8 log cfu/g (Table 2). We assume that MPN/g and cfu/g are compatible units and take into account the standard deviations from these two studies; so far Salmonella concentration in pork from the fresh market has never been changed. Note that the Commission Regulation on the microbiological criteria for foodstuffs indicated that Salmonella was not detected in the area tested per pig carcass after dressing but before chilling by the EN/ISO 6579 analytical reference method [50].
In this study, PE as a result of the exposure assessment step was derived from PPROB and the concentration of either susceptible or QR Salmonella, including the pork consumption of the Thai population, as shown in Equation (7) [6]. An alternative model to determine human exposure to AMR hazards per person per day requires additional parameters such as cross-contamination, which is dependent upon transfer rates between the food product and the environment [51]. The PE of QR Salmonella from fresh-market pork is considered low at 3 × 10−2, while the PE of QR Salmonella from modern-trade pork at 2 × 10−7 is considered negligible [52]. These results indicate that the PE of QR Salmonella from fresh market and modern trade followed the magnitude of both PPROB and the concentration of QR Salmonella.
In terms of hazard characterization, the consequence of hazard was determined by the dose–response model while AMR RA additionally includes the consequence of AMR [4] as the probability of mortality given illness (PMI) in this study. The PMI of QR Salmonella in the fresh market is much higher than PMI in the modern trade (Table 3), primarily because the probability of exposure (PE) of QR Salmonella in the fresh market is higher than the PE in the modern trade. In general, the PE model is determined by PPROB and the concentration (CS) of Salmonella (Equation (7)). This indicates that the adverse health effect of QR Salmonella from consuming fresh-market pork was higher than that from consuming modern-trade pork in Chiang Mai province.
So far, there have been very few risk assessments evaluating human health effects due to AMR Salmonella. One of these studies was the risk assessment of AMR Salmonella related to cattle [53,54]. A qualitative approach evaluated the additional risk of QR Salmonella recovered from minced pork as high [25]. Another qualitative risk assessment of human health effects from QR Salmonella Typhimurium in the EU upon using a (fluoro)quinolone in livestock (not necessarily swine) suggested the risk was low [55]. However, a quantitative risk assessment evaluated the human health effects of multi-resistant Salmonella Typhimurium DT104-contaminated Danish pork sausage [24]. The risk of salmonellosis from consuming such dry-cured pork sausages was in the range of 2.5 × 10−8–1.9 × 10−6, whereas in our study the mean mortality risks of QR Salmonella from modern-trade and fresh-market pork were as low as 7.4 × 10−17 and 2.0 × 10−7, respectively.
A previous study in Thailand reported that the annual mean mortality rate in 2009 (calculated from an average of the annual mortality cases of four major AMR bacteria (Acinetobacter baumannii, Staphylococcus aureus, Klebsiella pneumoniae, and E. coli)) was about 14.8 per 100,000 Thai population and was assumed to be the annual mean mortality rate for AMR salmonellosis [3]. In this study, the annual mortality rates for QR salmonellosis simulated by the risk assessment models were in the range of 0–32, which is in line with a previous study. The risk assessment models used in both exposure assessment and hazard characterization were applicable to evaluate the adverse health effects of AMR Salmonella in Thailand.

5. Conclusions

As far as we are aware, this is the first study of the quantitative microbial risk assessment of QR Salmonella in retail pork using a Monte Carlo simulation to comparatively report the human health adverse effects of susceptible and QR Salmonella from consuming retail pork from fresh market and modern trade, particularly in Thailand. The PPROB of both susceptible and QR Salmonella from the retail market are higher than the PPROB from modern trade. Likewise, the risk estimate in terms of the annual mortality rate of QR Salmonella from the fresh market is higher than that of QR Salmonella from modern trade and is also in line with a previous study reporting the mortality rate of AMR pathogens. The risk assessment models used in this study fit for evaluating the adverse health effects of QR Salmonella in Thailand and that of other foodborne AMR pathogens.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/foods11192942/s1, Table S1. Primers for carbapenemase genes; Table S2. Primers for plasmid-mediated colistin resistance genes (mcr-1-mcr-9); Table S3. Primers for plasmid-mediated quinolone resistance (PMQR) genes; Table S4. Primers for quinolone resistance-determining region (QRDR) genes. References [56, 57, 58] are cited in the Supplementary Materials.

Author Contributions

Conceptualization, S.N.; methodology, C.P., A.K., P.B., P.C. (Peechanika Chopjitt), R.H., P.C. (Piyarat Chansiripornchai) and N.S.; investigation, C.P., A.K., P.B., P.C. (Peechanika Chopjitt), R.H., P.C. (Piyarat Chansiripornchai) and N.S.; formal analysis, S.N.; writing—original draft preparation, C.P. and S.N.; writing—review and editing, C.P., A.K., P.B., P.C. (Peechanika Chopjitt), R.H., P.C. (Piyarat Chansiripornchai) and N.S.; supervision and review, A.K. and S.N.; project administration, S.N.; funding acquisition, A.K., P.B., P.C. (Peechanika Chopjitt) and R.H. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the National Research Council of Thailand (NRCT).

Data Availability Statement

All the data presented within the article is available upon request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Snary, E.; McEwen, S. Antimicrobial Resistance Risk Assessment. In Guide to Antimicrobial Use in Animals; Guardabassi, L., Jensen, L.B., Kruse, H., Eds.; Blackwell Publishing, Ltd.: Oxford, UK, 2008; pp. 27–43. [Google Scholar] [CrossRef]
  2. Booton, R.D.; Meeyai, A.; Alhusein, N.; Buller, H.; Feil, E.; Lambert, H.; Mongkolsuk, S.; Pitchforth, E.; Reyher, K.K.; Sakcamduang, W.; et al. One Health drivers of antibacterial resistance: Quantifying the relative impacts of human, animal and environmental use and transmission. One Health 2021, 12, 8. [Google Scholar] [CrossRef]
  3. Phumart, P.; Phodha, T.; Thamlikitkul, V.; Riewpaiboon, A.; Prakongsai, P.; Limwattananon, S. Health and economic impacts of antimicrobial resistant infections in Thailand: A Preliminary Study. J. Health Syst. Res. 2012, 6, 360. [Google Scholar]
  4. Codex Alimentarius Commission (CAC). Guidelines for Risk Analysis of Foodborne Antimicrobial Resistance (CXG 77-2011); FAO: Rome, Italy, 2011; p. 28. [Google Scholar]
  5. FAO; WHO. Principles and Guidelines for the Conduct of Microbiological Risk Assessment (CAC/GL-30), 2nd ed.; FAO: Rome, Italy; WHO: Geneva, Switzerland, 1999; p. 10. [Google Scholar]
  6. Khantasup, K.; Tungwongjulaniam, C.; Theerawat, R.; Lamaisri, T.; Piyalikit, K.; Nuengjamnong, C.; Nuanualsuwan, S. Cross-sectional risk assessment of zoonotic Streptococcus suis in pork and swine blood in Nakhon Sawan Province in northern Thailand. Zoonoses Public Health 2022, 69, 625–634. [Google Scholar] [CrossRef] [PubMed]
  7. Hanson, R.; Kaneene, J.B.; Padungtod, P.; Hirokawa, K.; Zeno, C. Prevalence of Salmonella and E. coli, and their resistance to antimicrobial agents, in farming communities in northern Thailand. Southeast Asian J. Trop. Med. Public Health 2002, 33 (Suppl. S3), 126. [Google Scholar]
  8. Padungtod, P.; Kaneene, J.B. Salmonella in food animals and humans in northern Thailand. Int. J. Food Microbiol. 2006, 108, 354. [Google Scholar] [CrossRef]
  9. Sanpong, P.; Theeragool, G.; Wajjwalku, W.; Amavisit, P. Characterization of multiple-antimicrobial resistant Salmonella isolated from pig farms in Thailand. Agric. Nat. Resour. 2010, 44, 651. [Google Scholar]
  10. Tacconelli, E.; Carrara, E.; Savoldi, A.; Harbarth, S.; Mendelson, M.; Monnet, D.L.; Pulcini, C.; Kahlmeter, G.; Kluytmans, J.; Carmeli, Y.; et al. Discovery, research, and development of new antibiotics: The WHO priority list of antibiotic-resistant bacteria and tuberculosis. Lancet Infect. Dis. 2018, 18, 327. [Google Scholar] [CrossRef]
  11. Poomchuchit, S.; Kerdsin, A.; Chopjitt, P.; Boueroy, P.; Hatrongjit, R.; Akeda, Y.; Tomono, K.; Nuanualsuwan, S.; Hamada, S. Fluoroquinolone resistance in non-typhoidal Salmonella enterica isolated from slaughtered pigs in Thailand. J. Med. Microbiol. 2021, 70, 9. [Google Scholar] [CrossRef] [PubMed]
  12. Angkititrakul, S.; Chomvarin, C.; Chaita, T.; Kanistanon, K.; Waethewutajarn, S. Epidemiology of antimicrobial resistance in Salmonella isolated from pork, chicken meat and humans in Thailand. Southeast Asian J. Trop. Med. Public Health 2005, 36, 1515. [Google Scholar]
  13. Holmberg, S.D.; Solomon, S.L.; Blake, P.A. Health and economic impacts of antimicrobial resistance. Rev. Infect. Dis. 1987, 9, 1078. [Google Scholar] [CrossRef]
  14. Molbak, K.; Baggesen, D.L.; Aarestrup, F.M.; Ebbesen, J.M.; Engberg, J.; Frydendahl, K.; Gerner-Smidt, P.; Petersen, A.M.; Wegener, H.C. An outbreak of multidrug-resistant, quinolone-resistant Salmonella enterica serotype typhimurium DT104. N. Engl. J. Med. 1999, 341, 1425. [Google Scholar] [CrossRef] [PubMed]
  15. Walker, R.A.; Lawson, A.J.; Lindsay, E.A.; Ward, L.R.; Wright, P.A.; Bolton, F.J.; Wareing, D.R.; Corkish, J.D.; Davies, R.H.; Threlfall, E.J. Decreased susceptibility to ciprofloxacin in outbreak-associated multiresistant Salmonella typhimurium DT104. Vet. Rec. 2000, 147, 395–396. [Google Scholar] [CrossRef] [PubMed]
  16. Meakins, S.; Fisher, I.S.; Berghold, C.; Gerner-Smidt, P.; Tschape, H.; Cormican, M.; Luzzi, I.; Schneider, F.; Wannett, W.; Coia, J.; et al. Antimicrobial drug resistance in human nontyphoidal Salmonella isolates in Europe 2000–2004: A report from the Enter-net International Surveillance Network. Microb. Drug Resist. 2008, 14, 35. [Google Scholar] [CrossRef] [PubMed]
  17. Threlfall, E.J.; Ward, L.R.; Skinner, J.A.; Graham, A. Antimicrobial drug resistance in non-typhoidal Salmonellas from humans in England and Wales in 1999: Decrease in multiple resistance in Salmonella enterica serotypes Typhimurium, Virchow, and Hadar. Microb. Drug Resist. 2000, 6, 325. [Google Scholar] [CrossRef] [PubMed]
  18. Varma, J.K.; Greene, K.D.; Ovitt, J.; Barrett, T.J.; Medalla, F.; Angulo, F.J. Hospitalization and antimicrobial resistance in Salmonella outbreaks, 1984–2002. Emerg. Infect. Dis. 2005, 11, 946. [Google Scholar] [CrossRef] [PubMed]
  19. Crook, P.D.; Aguilera, J.F.; Threlfall, E.J.; O’Brien, S.J.; Sigmundsdottir, G.; Wilson, D.; Fisher, I.S.; Ammon, A.; Briem, H.; Cowden, J.M.; et al. A European outbreak of Salmonella enterica serotype Typhimurium definitive phage type 204b in 2000. Clin. Microbiol. Infect. 2003, 9, 845. [Google Scholar] [CrossRef]
  20. Horby, P.W.; O’Brien, S.J.; Adak, G.K.; Graham, C.; Hawker, J.I.; Hunter, P.; Lane, C.; Lawson, A.J.; Mitchell, R.T.; Reacher, M.H.; et al. A national outbreak of multi-resistant Salmonella enterica serovar Typhimurium definitive phage type DT 104 associated with consumption of lettuce. Epidemiol. Infect. 2003, 130, 178. [Google Scholar] [CrossRef]
  21. Bertrand, S.; Weill, F.X.; Cloeckaert, A.; Vrints, M.; Mairiaux, E.; Praud, K.; Dierick, K.; Wildemauve, C.; Godard, C.; Butaye, P.; et al. Clonal emergence of extended-spectrum beta-lactamase (CTX-M-2)-producing Salmonella enterica serovar Virchow isolates with reduced susceptibilities to ciprofloxacin among poultry and humans in Belgium and France (2000 to 2003). J. Clin. Microbiol. 2006, 44, 2903. [Google Scholar] [CrossRef]
  22. Cloeckaert, A.; Praud, K.; Doublet, B.; Bertini, A.; Carattoli, A.; Butaye, P.; Imberechts, H.; Bertrand, S.; Collard, J.M.; Arlet, G.; et al. Dissemination of an extended-spectrum-beta-lactamase blaTEM-52 gene-carrying IncI1 plasmid in various Salmonella enterica serovars isolated from poultry and humans in Belgium and France between 2001 and 2005. Antimicrob. Agents Chemother. 2007, 51, 1875. [Google Scholar] [CrossRef]
  23. Caffrey, N.; Invik, J.; Waldner, C.; Ramsay, D.; Checkley, S. Risk assessments evaluating foodborne antimicrobial resistance in humans: A scoping review. J. Microb. Risk Anal. 2019, 11, 46. [Google Scholar] [CrossRef]
  24. Alban, L.; Olsen, A.M.; Nielsen, B.; Sorensen, R.; Jessen, B. Qualitative and quantitative risk assessment for human salmonellosis due to multi-resistant Salmonella Typhimurium DT104 from consumption of Danish dry-cured pork sausages. Prev. Vet. Med. 2002, 52, 265. [Google Scholar] [CrossRef]
  25. Doménech, E.; Jiménez-Belenguer, A.; Pérez, R.; Ferrús, M.A.; Escriche, I. Risk characterization of antimicrobial resistance of Salmonella in meat products. J. Food Control 2015, 57, 23. [Google Scholar] [CrossRef]
  26. Collineau, L.; Chapman, B.; Bao, X.; Sivapathasundaram, B.; Carson, C.A.; Fazil, A.; Reid-Smith, R.J.; Smith, B.A. A farm-to-fork quantitative risk assessment model for Salmonella Heidelberg resistant to third-generation cephalosporins in broiler chickens in Canada. Int. J. Food Microbiol. 2020, 330, 21. [Google Scholar] [CrossRef] [PubMed]
  27. Li, X.; Liang, B.; Xu, D.; Wu, C.; Li, J.; Zheng, Y. Antimicrobial Resistance Risk Assessment Models and Database System for Animal-Derived Pathogens. Antibiotics 2020, 9, 829. [Google Scholar] [CrossRef] [PubMed]
  28. Andrews, W.H.; Jacobson, A.; Hammack, T. Chapter 5: Salmonella. In Bacteriological Analytical Manual (BAM); U.S. Food and Drug Administration: Silver Spring, MD, USA, 2018. [Google Scholar]
  29. Haas, C.N.; Rose, J.B.; Gerba, C.P. Quantitative Microbial Risk Assessment; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 1999; p. 449. [Google Scholar]
  30. FAO; WHO. Risk Assessments of Salmonella in Eggs and Broiler Chickens; FAO: Rome, Italy; WHO: Geneva, Switzerland, 2002; p. 302. [Google Scholar]
  31. M07-A11; Methods for Dilution Antimicrobial Susceptibility Tests for Bacteria that Grow Aerobically. 10th ed. Clinical and Laboratory Standards Institute (CLSI): Wayne, PA, USA, 2018.
  32. M100-S30; Performance Standards for Antimicrobial Susceptibility Testing. 30th ed. Clinical and Laboratory Standards Institute (CLSI): Wayne, PA, USA, 2020.
  33. Ciesielczuk, H.; Hornsey, M.; Choi, V.; Woodford, N.; Wareham, D.W. Development and evaluation of a multiplex PCR for eight plasmid-mediated quinolone-resistance determinants. J. Med. Microbiol. 2013, 62, 1827. [Google Scholar] [CrossRef] [PubMed]
  34. Lu, Y.; Zhao, H.; Liu, Y.; Zhou, X.; Wang, J.; Liu, T.; Beier, R.C.; Hou, X. Characterization of quinolone resistance in Salmonella enterica serovar Indiana from chickens in China. Poult. Sci. 2015, 94, 460. [Google Scholar] [CrossRef]
  35. Alali, W.Q.; Mann, D.A.; Beuchat, L.R. Viability of Salmonella and Listeria monocytogenes in delicatessen salads and hummus as affected by sodium content and storage temperature. J. Food Prot. 2012, 75, 1056. [Google Scholar] [CrossRef]
  36. Ruchusatsawat, K.; Nuengjamnong, C.; Tawatsin, A.; Thiemsing, L.; Kawidam, C.; Somboonna, N.; Nuanualsuwan, S. Quantitative Risk Assessments of Hepatitis A Virus and Hepatitis E Virus from Raw Oyster Consumption. Risk Anal. 2022, 42, 965. [Google Scholar] [CrossRef]
  37. Nuanualsuwan, S.; Songkasupa, T.; Boonpornprasert, P.; Suwankitwat, N.; Lohlamoh, W.; Nuengjamnong, C. Thermal Inactivation of African Swine Fever Virus in Swill. Front. Vet. Sci. 2022, 9, 8. [Google Scholar] [CrossRef]
  38. Bywater, R.J.; Casewell, M.W. An assessment of the impact of antibiotic resistance in different bacterial species and of the contribution of animal sources to resistance in human infections. J. Antimicrob. Chemother. 2000, 46, 645. [Google Scholar] [CrossRef]
  39. Presi, P.; Stark, K.D.; Stephan, R.; Breidenbach, E.; Frey, J.; Regula, G. Risk scoring for setting priorities in a monitoring of antimicrobial resistance in meat and meat products. Int. J. Food Microbiol. 2009, 130, 100. [Google Scholar] [CrossRef] [PubMed]
  40. Berends, B.R.; van den Bogaard, A.E.; Van Knapen, F.; Snijders, J.M. Human health hazards associated with the administration of antimicrobials to slaughter animals. Part II. An assessment of the risks of resistant bacteria in pigs and pork. Vet. Q. 2001, 23, 21. [Google Scholar] [CrossRef]
  41. Hald, T.; Lo Fo Wong, D.M.; Aarestrup, F.M. The attribution of human infections with antimicrobial resistant Salmonella bacteria in Denmark to sources of animal origin. Foodborne Pathog. Dis. 2007, 4, 326. [Google Scholar] [CrossRef]
  42. Sommer, H.M.; Aabo, S.; Christensen, B.B.; Saadby, P.; Nielsen, N.; Nørrung, B.; Wong, D.L.F. Risk Assessment of the Impact on Human Health Related to Multiresistant Salmonella Typhimurium DT104 from Slaughter Pigs; Institute of Food Safety and Nutrition: Mørkhøj, Denmark, 2003; p. 104.
  43. European Food Safety Authority (EFSA). Scientific Opinion of the Panel on Biological Hazards on a request from the European Food Safety Authority on foodborne antimicrobial resistance as a biological hazard. EFSA 2008, 765, 87. [Google Scholar]
  44. Aslam, M.; Checkley, S.; Avery, B.; Chalmers, G.; Bohaychuk, V.; Gensler, G.; Reid-Smith, R.; Boerlin, P. Phenotypic and genetic characterization of antimicrobial resistance in Salmonella serovars isolated from retail meats in Alberta, Canada. Food Microbiol. 2012, 32, 117. [Google Scholar] [CrossRef]
  45. Patchanee, P.; Tansiricharoenkul, K.; Buawiratlert, T.; Wiratsudakul, A.; Angchokchatchawal, K.; Yamsakul, P.; Yano, T.; Boonkhot, P.; Rojanasatien, S.; Tadee, P. Salmonella in pork retail outlets and dissemination of its pulsotypes through pig production chain in Chiang Mai and surrounding areas, Thailand. Prev. Vet. Med. 2016, 130, 105. [Google Scholar] [CrossRef]
  46. Prasertsee, T.; Chokesajjawatee, N.; Santiyanont, P.; Chuammitri, P.; Deeudom, M.; Tadee, P.; Patchanee, P. Quantification and rep-PCR characterization of Salmonella spp. in retail meats and hospital patients in Northern Thailand. Zoonoses Public Health 2019, 66, 309. [Google Scholar] [CrossRef]
  47. Sanguankiat, A.; Pinthong, R.; Padungtod, P.; Baumann, M.P.; Zessin, K.H.; Srikitjakarn, L.; Fries, R. A cross-sectional study of Salmonella in pork products in Chiang Mai, Thailand. Foodborne Pathog. Dis. 2010, 7, 878. [Google Scholar] [CrossRef]
  48. Tadee, P.; Boonkhot, P.; Pornruangwong, S.; Patchanee, P. Comparative phenotypic and genotypic characterization of Salmonella spp. in pig farms and slaughterhouses in two provinces in northern Thailand. PLoS ONE 2015, 10, 11. [Google Scholar] [CrossRef]
  49. Tadee, P.; Kumpapong, K.; Sinthuya, D.; Yamsakul, P.; Chokesajjawatee, N.; Nuanualsuwan, S.; Pornsukarom, S.; Molla, B.Z.; Gebreyes, W.A.; Patchanee, P. Distribution, quantitative load and characterization of Salmonella associated with swine farms in upper-northern Thailand. J. Vet. Sci. 2014, 15, 334. [Google Scholar] [CrossRef]
  50. European Commission. Microbiological Criteria for Foodstuffs: Commission Regulation (EC) No 2073/2005 on 2005; European Commission: Brussels, Belgium, 2005; pp. 1–26. [Google Scholar]
  51. Hald, T.; Duarte, A.R.; Stärk, K. Quantifying human exposure to antimicrobial resistance from animals and food. In Proceedings of the International Conference on Animal Health Surveillance, Wellington, New Zealand, 30 April–4 May 2017; p. 2. [Google Scholar]
  52. Biosecurity Australia. Generic Import Risk Analysis Report for Chicken Meat: Final Report; Biosecurity Australia: Canberra, Australia, 2008.
  53. Hurd, H.S.; Vaughn, M.B.; Holtkamp, D.; Dickson, J.; Warnick, L. Quantitative risk from fluoroquinolone-resistant Salmonella and Campylobacter due to treatment of dairy heifers with enrofloxacin for bovine respiratory disease. J. Foodborne Pathog. 2010, 7, 1322. [Google Scholar] [CrossRef] [PubMed]
  54. Singer, R.S.; Ruegg, P.L.; Bauman, D.E. Quantitative Risk Assessment of Antimicrobial-Resistant Foodborne Infections in Humans Due to Recombinant Bovine Somatotropin Usage in Dairy Cows. J. Food Prot. 2017, 80, 1116. [Google Scholar] [CrossRef] [PubMed]
  55. European Agency for the Evaluation of Medicinal Products. Antibiotic Resistance in the European Union Associated with Therapeutic Use of Veterinary Medicines: Report and Qualitative Risk Assessment by the Committee for Veterinary Medicinal Products; European Agency for the Evaluation of Medicinal Products: Amsterdam, The Netherlands, 1999.
  56. Poirel, L.; Walsh, T.R.; Cuvillier, V.; Nordmann, P. Multiplex PCR for detection of acquired carbapenemase genes. Diagn. Microbiol. Infect Dis 2011, 70, 119–123. [Google Scholar] [CrossRef]
  57. Hatrongjit, R.; Kerdsin, A.; Akeda, Y.; Hamada, S. Detection of plasmid-mediated colistin-resistant and carbapenem-resistant genes by multiplex PCR. MethodsX 2018, 5, 532–536. [Google Scholar] [CrossRef]
  58. Khanawapee, A.; Kerdsin, A.; Chopjitt, P.; Boueroy, P.; Hatrongjit, R.; Akeda, Y.; Tomono, K.; Nuanualsuwan, S.; Hamada, S. Distribution and Molecular Characterization of Escherichia coli Harboring mcr Genes Isolated from Slaughtered Pigs in Thailand. Microb. Drug Resist. 2021, 27, 971–979. [Google Scholar] [CrossRef]
Figure 1. Annual mortality cases from susceptible Salmonella-contaminated pork from the fresh market in Chiang Mai.
Figure 1. Annual mortality cases from susceptible Salmonella-contaminated pork from the fresh market in Chiang Mai.
Foods 11 02942 g001
Figure 2. Annual mortality cases from QR Salmonella-contaminated pork from the fresh market in Chiang Mai.
Figure 2. Annual mortality cases from QR Salmonella-contaminated pork from the fresh market in Chiang Mai.
Foods 11 02942 g002
Figure 3. Annual mortality cases from susceptible Salmonella-contaminated pork from the modern trade in Chiang Mai.
Figure 3. Annual mortality cases from susceptible Salmonella-contaminated pork from the modern trade in Chiang Mai.
Foods 11 02942 g003
Table 1. No. of Salmonella positive samples collected from retailers in Chiang Mai province.
Table 1. No. of Salmonella positive samples collected from retailers in Chiang Mai province.
RetailNo. of SalmonellaTotal
SusceptibleQR
Fresh market302858 (n =100)
Modern trade606 (n = 50)
Table 2. PPROB and mean concentration of contaminants in the pork samples.
Table 2. PPROB and mean concentration of contaminants in the pork samples.
RetailPPROB (%)Mean Concentration ± SD (log cfu/g)
Salmonella spp.TotalSalmonella spp.Total *
SusceptibleQRSusceptibleQR
Fresh market30.428.457.81.5 ± 0.82.1 ± 0.71.8 ± 0.8
Modern trade13.51.96.91.9 ± 0.901.9 ± 0.9
* Accounted for only positive samples.
Table 3. Probabilities of exposure (PE), illness (PI), and mortality given illness (PMI) from susceptible and QR Salmonella spp.
Table 3. Probabilities of exposure (PE), illness (PI), and mortality given illness (PMI) from susceptible and QR Salmonella spp.
RetailPEPIPMI
Salmonella spp.Salmonella spp.Salmonella
SusceptibleQRSusceptibleQRSusceptibleQR
Fresh market0.0200.0301.8 × 10−43.1 × 10−42.3 × 10−75.4 × 10−6
Modern trade0.0162 × 10−73.2 × 10−42.2 × 10−84.2 × 10−73.9 × 10−10
Table 4. Descriptive statistics of risk estimate (PFM) and annual mortality rate (MAFM) from consuming pork contaminated with susceptible and AMR Salmonella spp.
Table 4. Descriptive statistics of risk estimate (PFM) and annual mortality rate (MAFM) from consuming pork contaminated with susceptible and AMR Salmonella spp.
Retail Risk EstimateAnnual Cases *
Salmonella spp.Salmonella spp.
SusceptibleAMRSusceptibleAMR
Fresh marketmin5.3 × 10−138.8 × 10−11<1<1
mean5.7 × 10−92.0 × 10−7<1 a7 b
max2.7 × 10−88.8 × 10−7132
Modern trademin1.5 × 10−124.2 × 10−21<1<1
mean7.9 × 10−97.4 × 10−17<1 c<1 d
max4.0 × 10−87.6 × 10−162<1
* Mean annual cases per 100,000 population (PAFM) with different letters implies that there are statistically significant differences (p < 0.05) (letters a through d).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pulsrikarn, C.; Kedsin, A.; Boueroy, P.; Chopjitt, P.; Hatrongjit, R.; Chansiripornchai, P.; Suanpairintr, N.; Nuanualsuwan, S. Quantitative Risk Assessment of Susceptible and Ciprofloxacin-Resistant Salmonella from Retail Pork in Chiang Mai Province in Northern Thailand. Foods 2022, 11, 2942. https://doi.org/10.3390/foods11192942

AMA Style

Pulsrikarn C, Kedsin A, Boueroy P, Chopjitt P, Hatrongjit R, Chansiripornchai P, Suanpairintr N, Nuanualsuwan S. Quantitative Risk Assessment of Susceptible and Ciprofloxacin-Resistant Salmonella from Retail Pork in Chiang Mai Province in Northern Thailand. Foods. 2022; 11(19):2942. https://doi.org/10.3390/foods11192942

Chicago/Turabian Style

Pulsrikarn, Chaiwat, Anusak Kedsin, Parichart Boueroy, Peechanika Chopjitt, Rujirat Hatrongjit, Piyarat Chansiripornchai, Nipattra Suanpairintr, and Suphachai Nuanualsuwan. 2022. "Quantitative Risk Assessment of Susceptible and Ciprofloxacin-Resistant Salmonella from Retail Pork in Chiang Mai Province in Northern Thailand" Foods 11, no. 19: 2942. https://doi.org/10.3390/foods11192942

APA Style

Pulsrikarn, C., Kedsin, A., Boueroy, P., Chopjitt, P., Hatrongjit, R., Chansiripornchai, P., Suanpairintr, N., & Nuanualsuwan, S. (2022). Quantitative Risk Assessment of Susceptible and Ciprofloxacin-Resistant Salmonella from Retail Pork in Chiang Mai Province in Northern Thailand. Foods, 11(19), 2942. https://doi.org/10.3390/foods11192942

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop