Abstract
Background:
Salmonellosis, a disease caused by Salmonella, is a significant public health concern and economic burden worldwide. The ability of various Salmonella serovars to cause disease is closely linked to the virulence genes they possess.Objectives:
The aim of this study was to develop a multiplex fluorescence PCR for detecting ten major virulence genes (ssaR, spvC, pefA, sipA, fimA, sifA, sopE2, sopB, prgH, and stn) in Salmonella.Methods:
Primer pairs specific to ten target virulence genes were designed using Primer Premier 5.0 and distributed across two reaction tubes. The multiplex fluorescence PCR was optimized by adjusting one factor at a time.Results:
A total of sixty Salmonella strains were analyzed using the newly developed multiplex fluorescence PCR. All strains contained seven or more of the tested virulence genes. The positive rates of virulence genes ssaR, sipA, sopE2, sopB, prgH, and stn were high, each at 100%. The positive rate of sifA was also relatively high at 81.67%. However, the positive rates of spvC at 5% and pefA at 3.33% were relatively low.Conclusions:
The newly developed multiplex fluorescence PCR provides a straightforward, cost-effective, and high-throughput solution for detecting virulence genes in Salmonella. It has the potential to become a routine method for analyzing Salmonella virulence genes.Keywords
1. Background
Salmonellosis, a disease caused by Salmonella, represents a significant public health problem and economic burden worldwide (1-3). The transmission routes of Salmonella to humans are diverse, including consumption of contaminated food or water and contact with infected animals and humans. According to estimates from the World Health Organization (WHO), more than 90 million cases of human salmonellosis occur globally each year (4). Therefore, the control and prevention of salmonellosis are crucial.
The pathogenicity of Salmonella serovars depends on the presence of a variety of virulence genes, such as ssaR, spvC, pefA, sipA, fimA, sifA, sopE2, sopB, prgH, and stn. The ssaR and sifA genes are essential for the intracellular survival and subsequent replication of Salmonella (5-7), while the pefA, sipA, fimA, sopE2, sopB, and prgH genes are involved in adhesion and invasion of host cells (7-10). The spvC gene can suppress the innate immunity of infected hosts (11), and the stn gene is implicated in the complex mechanisms of diarrhea (12). Additionally, virulence genes can be transferred from one strain to another, enhancing the pathogenicity of Salmonella strains and potentially leading to disease outbreaks. For example, during the mid to late 1980s, there was an outbreak of typhoid fever in China characterized by severe symptoms and high mortality rates.
Subsequent studies identified the main causative agent of the epidemic as Salmonellaenterica serovar typhi (S. typhi) strains carrying the Salmonella plasmid virulence (spv) gene, which is typically absent in S. typhi (13). Investigating the presence of virulence genes in Salmonella strains is crucial for identifying pathogenicity, understanding potential transfer mechanisms, and assessing the risks posed by these strains. Thus, developing an efficient method for detecting Salmonella virulence genes is necessary for disease prevention and control.
PCR, including conventional and fluorescence PCR, is one of the most commonly used techniques for gene analysis. Fluorescence PCR introduces a fluorescent reporter (dsDNA binding dye or fluorophore-labeled probe) that binds to the amplification product and indicates its presence by fluorescence, allowing for real-time detection of the amplification product during the PCR reaction. The entire fluorescence PCR process can be completed in a closed-tube setting without any manual operation. Compared to conventional PCR, fluorescence PCR is simpler, faster, and easily automated for high-throughput testing (14). However, reports on the application of fluorescence PCR for analyzing Salmonella virulence genes are limited.
2. Objectives
The aim of this study is to develop a multiplex fluorescence PCR for detecting ten major virulence genes (ssaR, spvC, pefA, sipA, fimA, sifA, sopE2, sopB, prgH and stn) in Salmonella.
3. Methods
3.1. Strains
A total of 60 Salmonella strains, 2 Vibrio parahaemolyticus strains, and 3 Escherichia coli strains were utilized in this study. One S. typhimurium strain and two Salmonella London strains were isolated from meat, while other strains were isolated from patients with diarrhea. All strains were isolated and confirmed by biochemical characterization according to GB 4789 (national standards for microbiological examination of food in the People's Republic of China) and WS 271 (health standards for the diagnosis of infectious diarrhea in the People's Republic of China). Five Salmonella strains, comprising 3 strains of S. enteritidisand 2 strains of S. typhimurium (Table 1), were chosen as positive controls. Additionally, two V. parahaemolyticus and three E. coli strains served as negative controls. The virulence gene profiles of these control strains were detected by previously reported conventional PCR (15-20). Genomic DNA from all strains was extracted using the KAPA Express extract kit, following the manufacturer’s instructions.
Virulence Genes Profiles of Control Strains
| Isolates | Isolates No. | ssaR | spvC | pefA | sipA | fimA | sifA | sopE2 | sopB | prgH | stn |
|---|---|---|---|---|---|---|---|---|---|---|---|
| Salmonella enteritidis | 1 | + | - | - | + | + | + | + | + | + | + |
| 1 | + | - | - | + | + | - | + | + | + | + | |
| 1 | + | + | - | + | + | - | + | + | + | + | |
| Salmonella typhimurium | 1 | + | + | + | + | + | + | + | + | + | + |
| 1 | + | - | - | + | + | + | + | + | + | + | |
| Vibrio parahaemolyticus | 2 | - | - | - | - | - | - | - | - | - | - |
| Escherichia coli | 3 | - | - | - | - | - | - | - | - | - | - |
3.2. Oligonucleotide Primers
The DNA sequences of ten major virulence genes in Salmonella were downloaded from GenBank and aligned using Clustal X 1.8. Based on this sequence alignment, relatively conserved regions were selected for primer design using Primer Premier 5.0. The details of the ten primer pairs targeting these genes are shown in Table 2.
Primer Sequences Used in This Study
| Variables | Target Gene | Primer | Sequence |
|---|---|---|---|
| Tube1 | ssaR | ssaR-F | 5’- CAAAAAAACTCTGAAGAGAAGGAAGCC-3’ |
| ssaR-R | 5’- GCCTGCGTTAACTGACTCACCGT-3’ | ||
| spvC | spvC-F | 5’- TGGGGCGGAAATACCATCTACAAAT-3’ | |
| spvC-R | 5’- CTGAATAGTCAGGCACATCCAGCG-3’ | ||
| pefA | pefA-F | 5’- ACCAGCCGGGTAATTTTGTTGACT-3’ | |
| pefA-R | 5’- CAGCAGAAGCCCAGGAAACAGTG-3’ | ||
| sipA | sipA-F | 5’- CACGTCAGAAAAGGGCACTACGG-3’ | |
| sipA-R | 5’- CTCGACGCCAGCCAGTTGACTG-3’ | ||
| fimA | fimA-F | 5’-TAATGACCTCTACTATTGCGAGTCTGATGTT-3’ | |
| fimA-R | 5’- AGAAAGCCACGGCAGCGTTG-3’ | ||
| Tube2 | sifA | sifA-F | 5’- TTAAAAAGTGAAATCCTTACCAACTCCCC-3’ |
| sifA-R | 5’- TCCGCTTTTGCTTTGCCAGTAGAA-3’ | ||
| sopE2 | sopE2-F | 5’- TAACATAACACTATCCACCCAGCACTACAG-3’ | |
| sopE2-R | 5’- TAAAACGATCTGACAGACTCGTCGACA-3’ | ||
| sopB | sopB-F | 5’- ATAAAGATGGCGATCTACAGACGGTAAA-3’ | |
| sopB-R | 5’- CCCACATTAAATGCGGCGACG-3’ | ||
| prgH | prgH-F | 5’- CGGCTGTGAGTTTCCATTGCTG-3’ | |
| prgH-R | 5’- TTCTAACTTCTCAGGCTGCTCGGG-3’ | ||
| stn | stn-F | 5’- CGCCAGCCCTTTGATGGACG-3’ | |
| stn-R | 5’- GGCATCAGCGTTATCAGCGCTC-3’ |
3.3. Optimization of Multiplex Fluorescence PCR
Primer pairs for ten virulence genes were distributed across two reaction tubes, with the distribution of these pairs detailed in Table 2. The multiplex fluorescence PCR was optimized by adjusting one parameter at a time. The parameters evaluated included the target primer concentration, which ranged from 0.05 µM to 1 µM, and the annealing temperature, which varied from 55 °C to 65 °C.
The multiplex fluorescence PCR was conducted using a 20 µL reaction mixture for each tube, consisting of 1 µL of DNA template, 1 µL of EvaGreen dye, 10 µL of KAPA 2G Fast Multiplex PCR Mix, appropriate amounts of target primers, and the remaining volume filled with distilled water. The PCR amplification took place in a Bio-Rad CFX 96TM real-time PCR system with the following program: An initial denaturation at 95 °C for 3 minutes, followed by 30 cycles of denaturation at 95 °C for 15 seconds, annealing for 30 seconds at the specified annealing temperature, and extension at 72 °C for 10 seconds.
3.4. Melting Temperature Curve Analysis
Fluorescence melting temperature curve analysis was conducted following PCR amplification. The PCR products were initially cooled to 78 °C and subsequently heated to 91 °C at a rate of 0.1 °C per second. Fluorescence signals were continuously monitored throughout the melting temperature curve analysis process.
4. Results
4.1. Optimization of Multiplex Fluorescence PCR
To develop a multiplex fluorescence PCR capable of detecting ten primary virulence genes in Salmonella, a systematic study was conducted to optimize the conditions. The optimal primer concentrations were determined as follows: In tube 1, 0.4 µM for ssaR, 0.4 µM for spvC, 0.3 µM for pefA, 0.1 µM for sipA, and 0.1 µM for fimA; in tube 2, 0.4 µM for sifA, 0.2 µM for sopE2, 0.2 µM for sopB, 0.05 µM for prgH, and 0.1 µM for stn. The optimal annealing temperature was established at 60°C.
4.2. Melting Temperature Curve Analysis of PCR Products
After PCR amplification and melting temperature curve analysis, ten distinct melting peaks were observed. The melting peak temperatures for the PCR products were as follows: ssaR primer pairs at 80.7 ± 0.3°C, spvC primer pairs at 82.8 ± 0.2°C, pefA primer pairs at 85.5 ± 0.1°C, sipA primer pairs at 87.4 ± 0.2°C, and fimA primer pairs at 89.0 ± 0.2°C (Figure 1). Additionally, the melting peak temperatures for sifA primer pairs were at 79.6 ± 0.2°C, sopE2 primer pairs at 81.8 ± 0.3°C, sopB primer pairs at 84.3 ± 0.2°C, prgH primer pairs at 86.3 ± 0.1°C, and stn primer pairs at 88.1 ± 0.3 °C (Figure 2).
Melting peaks of virulence genes (ssaR, spvC, pefA, sipA and fimA) in tube 1
Melting peaks of virulence genes (sifA, sopE2, sopB, prgH and stn) in tube 2
4.3. Evaluation of Specificity of Multiplex Fluorescence PCR
Analysis of five positive controls and five negative controls using the newly developed multiplex fluorescence PCR showed that the five positive controls produced specific melting peaks, whereas no obvious melting peaks were observed in the five negative controls. The virulence genes profiles in the five positive controls detected by the newly developed multiplex fluorescence PCR were consistent with our previous results achieved by conventional PCR.
4.4. Analysis of Salmonella Strains
The distribution of 10 virulence genes in 60 Salmonella strains is summarized in Table 3. All strains contained seven or more of the tested virulence genes. The positive rates for virulence genes ssaR, sipA, sopE2, sopB, prgH, and stn were 100%, while sifA was 81.67%. However, the positive rates for spvC and pefA were relatively low, at 5% and 3.33% respectively. These results were confirmed by previously reported conventional PCR (15-20).
Virulence Genes Profiles of Tested Salmonella Strains
| Isolates | Isolates No. | ssaR | spvC | pefA | sipA | fimA | sifA | sopE2 | sopB | prgH | stn |
|---|---|---|---|---|---|---|---|---|---|---|---|
| Salmonella enteritidis | 11 | + | - | - | + | + | + | + | + | + | + |
| Salmonella enteritidis | 5 | + | - | - | + | + | - | + | + | + | + |
| Salmonella enteritidis | 1 | + | + | - | + | + | - | + | + | + | + |
| Salmonella typhimurium | 1 | + | + | + | + | + | + | + | + | + | + |
| Salmonella typhimurium | 8 | + | - | _ | + | + | + | + | + | + | + |
| Salmonella stanley | 7 | + | - | _ | + | + | + | + | + | + | + |
| Salmonella agona | 3 | + | - | _ | + | + | - | + | + | + | + |
| Salmonella london | 5 | + | - | _ | + | + | + | + | + | + | + |
| Salmonella london | 1 | + | - | + | + | - | + | + | + | + | |
| Salmonella essen | 3 | + | - | - | + | + | + | + | + | + | + |
| Salmonella give | 1 | + | + | + | + | + | + | + | + | + | + |
| Salmonella tshiongwe | 1 | + | - | - | + | + | - | + | + | + | + |
| Salmonella montevideo | 2 | + | - | - | + | + | + | + | + | + | + |
| Salmonella abony | 3 | + | - | - | + | + | + | + | + | + | + |
| Salmonella weltevreden | 5 | + | - | - | + | + | + | + | + | + | + |
| Salmonella dublin | 3 | + | - | - | + | + | + | + | + | + | + |
5. Discussion
Infection with Salmonella is still one of the most important public health problems in the world. It can cause diseases ranging from gastroenteritis to typhoid fever (21). The pathogenicity of Salmonella strains is determined by their abilities to adhere, invade and replicate inside host cells, as well as to escape from or neutralize the host’s defenses (7, 22). These abilities of Salmonella strains depend on the virulence genes that they possess. In many cases, when a Salmonella serotype acquires new virulence genes that were previously absent, it can enhance pathogenicity of Salmonella strains and cause disease outbreak (13). Therefore, analysis of virulence genes profiles of Salmonella strains is necessary for preventing outbreak of disease caused by Salmonella infection. As one of the most widely used techniques in genes analysis, several PCR methods have been developed for Salmonella virulence genes analysis (15-20). However, these reported PCR methods for Salmonella virulence genes analysis are conventional PCR.
The drawbacks of conventional PCR are obvious. First, post-PCR analysis depend on gel electrophoresis is necessary for conventional PCR. The procedure of post-PCR analysis is complex and labor-intensive, which results in low throughput of the conventional PCR. Second, conventional PCR has the relative high risk of cross contamination as the essential step of opening PCR tube for post-PCR analysis. Therefore, other method have been developed for Salmonella virulence genes analysis. Though these methods have their advantages, they can’t completely replace PCR in the field of Salmonella virulence genes analysis. For example, microarray has a higher throughput than conventional PCR in Salmonella virulence genes analysis, PCR is still needed to validate the microarray data demonstrating weaker signals.
In this study, we developed a multiplex fluorescence PCR for Salmonella virulence genes analysis. Fluorescence PCR differs from the traditional PCR. The whole detection process of fluorescence PCR could be completed in a closed-tube setting and the post-PCR analysis depend on gel electrophoresis is not required.Without post-PCR analysis, fluorescence PCR increases assay throughput and reduces the risk of cross contamination (23). Thus, our newly developed multiplex fluorescence PCR for Salmonella virulence genes analysis overcomes the drawbacks of conventional PCR. Our multiplex fluorescence PCR encompasses ten main virulence genes (ssaR, spvC, pefA, sipA, fimA, sifA, sopE2, sopB, prgH, and stn), which play crucial roles in the diseases caused by Salmonella. These genes are significant indicators of pathogenicity, facilitating the ability of Salmonella to adhere to, invade, and replicate inside host cells, and to escape from or neutralize the host's defenses (5, 12).
In this study, all analyzed strains contained seven or more of the tested virulence genes. The positive rates for virulence genes such as ssaR, sipA, sifA, sopE2, sopB, prgH, and stn were relatively high, whereas those for spvC and pefA were relatively low. The distribution of most virulence genes in our study aligns with results from previous reports. For instance, Yue et al. (15) analyzed 61 Salmonella isolates from the feces of children with acute diarrhea and found that all isolates (100%) tested positive for prgH, 51 isolates (83.61%) for sopB, and 11 isolates (18.03%) for pefA. Another study by Qiao et al. (24) reported prevalence rates for the virulence genes sipA, pefA, and spvC at 77.6%, 10.3%, and 1.9%, respectively.
5.1. Conclusions
In summary, the newly developed multiplex fluorescence PCR provides a straightforward, cost-effective, and high-throughput method for detecting virulence genes in Salmonella. Thus, it has the potential to become a routine method for analyzing Salmonella virulence genes.
References
-
1.
Li W, Ren Q, Ni T, Zhao Y, Sang Z, Luo R, et al. Strategies adopted by Salmonella to survive in host: a review. Arch Microbiol. 2023;205(12):362. [PubMed ID: 37904066]. https://doi.org/10.1007/s00203-023-03702-w.
-
2.
Popa GL, Papa MI. Salmonella spp. infection - a continuous threat worldwide. Germs. 2021;11(1):88-96. [PubMed ID: 33898345]. [PubMed Central ID: PMC8057844]. https://doi.org/10.18683/germs.2021.1244.
-
3.
Balasubramanian R, Im J, Lee JS, Jeon HJ, Mogeni OD, Kim JH, et al. The global burden and epidemiology of invasive non-typhoidal Salmonella infections. Hum Vaccin Immunother. 2019;15(6):1421-6. [PubMed ID: 30081708]. [PubMed Central ID: PMC6663144]. https://doi.org/10.1080/21645515.2018.1504717.
-
4.
Xu L, He Q, Tang Y, Wen W, Chen L, Li Y, et al. Multi‑locus sequence and drug resistance analysis of Salmonella infection in children with diarrhea in Guangdong to identify the dominant ST and cause of antibiotic‑resistance. Exp Ther Med. 2022;24(5):678. [PubMed ID: 36337293]. [PubMed Central ID: PMC9623445]. https://doi.org/10.3892/etm.2022.11614.
-
5.
Morasi RM, da Silva AZ, Nunez KVM, Dantas STA, Faganello C, Juliano LCB, et al. Overview of antimicrobial resistance and virulence factors in Salmonella spp. isolated in the last two decades from chicken in Brazil. Food Res Int. 2022;162(Pt A):111955. [PubMed ID: 36461285]. https://doi.org/10.1016/j.foodres.2022.111955.
-
6.
Dantas STA, Camargo CH, Tiba-Casas MR, Vivian RC, Pinto J, Pantoja JCF, et al. Environmental persistence and virulence of Salmonella spp. Isolated from a poultry slaughterhouse. Food Res Int. 2020;129:108835. [PubMed ID: 32036904]. https://doi.org/10.1016/j.foodres.2019.108835.
-
7.
Karacan Sever N, Akan M. Molecular analysis of virulence genes of Salmonella Infantis isolated from chickens and turkeys. Microb Pathog. 2019;126:199-204. [PubMed ID: 30403968]. https://doi.org/10.1016/j.micpath.2018.11.006.
-
8.
Turner M, Van Hulzen L, Pietri JE. The gut microbiota induces melanin deposits that act as substrates for fimA-mediated aggregation of Salmonella Typhimurium and enhance infection of the German cockroach vector. Microbiol Spectr. 2023;11(5). e0211923. [PubMed ID: 37606369]. [PubMed Central ID: PMC10580948]. https://doi.org/10.1128/spectrum.02119-23.
-
9.
De Azevedo EC, Martins BTF, Tiba Casas MR, Possebon FS, Araujo Junior JP, Nero LA, et al. Multidrug resistance and virulence profiles of Salmonella isolated from swine lymph nodes. Microb Drug Resist. 2021;27(4):562-70. [PubMed ID: 32865485]. https://doi.org/10.1089/mdr.2020.0120.
-
10.
Yin J, Cheng Z, Xu Z, Zhi L, Zhang Y, Yuan X, et al. Contribution of prgH gene for Salmonella Pullorum to virulence and the expression of NLRP3, Caspase-1 and IL-1beta in chickens. Microb Pathog. 2022;171:105744. [PubMed ID: 36049651]. https://doi.org/10.1016/j.micpath.2022.105744.
-
11.
Hawwas HAE, Aboueisha AM, Fadel HM, El-Mahallawy HS. Salmonella serovars in sheep and goats and their probable zoonotic potential to humans in Suez Canal Area, Egypt. Acta Vet Scand. 2022;64(1):17. [PubMed ID: 35906669]. [PubMed Central ID: PMC9336019]. https://doi.org/10.1186/s13028-022-00637-y.
-
12.
Sabry MA, Abdel-Moein KA, Abdel-Kader F, Hamza E. Extended-spectrum beta-lactamase-producing Salmonella serovars among healthy and diseased chickens and their public health implication. J Glob Antimicrob Resist. 2020;22:742-8. [PubMed ID: 32623001]. https://doi.org/10.1016/j.jgar.2020.06.019.
-
13.
Wei L, Jin Q, Chu Y, Wu S, Huang R. Suppression of dendritic cell and T-cell activation by the pRST(9)(8) Salmonella plasmid. Mol Med Rep. 2015;11(3):2306-14. [PubMed ID: 25385146]. https://doi.org/10.3892/mmr.2014.2919.
-
14.
Koo SB, Chi HG, Kim JD, Kim YS, Park JS, Park CY, et al. Multiple compact camera fluorescence detector for real-time PCR devices. Sensors (Basel). 2021;21(21). [PubMed ID: 34770319]. [PubMed Central ID: PMC8587052]. https://doi.org/10.3390/s21217013.
-
15.
Yue M, Li X, Liu D, Hu X. Serotypes, antibiotic resistance, and virulence genes of Salmonella in children with diarrhea. J Clin Lab Anal. 2020;34(12). e23525. [PubMed ID: 32797660]. [PubMed Central ID: PMC7755775]. https://doi.org/10.1002/jcla.23525.
-
16.
Obe T, Nannapaneni R, Schilling W, Zhang L, McDaniel C, Kiess A. Prevalence of Salmonella enterica on poultry processing equipment after completion of sanitization procedures. Poult Sci. 2020;99(9):4539-48. [PubMed ID: 32867998]. [PubMed Central ID: PMC7598133]. https://doi.org/10.1016/j.psj.2020.05.043.
-
17.
Fardsanei F, Soltan Dallal MM, Douraghi M, Memariani H, Bakhshi B, Zahraei Salehi T, et al. Antimicrobial resistance, virulence genes and genetic relatedness of Salmonella enterica serotype Enteritidis isolates recovered from human gastroenteritis in Tehran, Iran. J Glob Antimicrob Resist. 2018;12:220-6. [PubMed ID: 29045813]. https://doi.org/10.1016/j.jgar.2017.10.005.
-
18.
Krawiec M, Kuczkowski M, Kruszewicz AG, Wieliczko A. Prevalence and genetic characteristics of Salmonella in free-living birds in Poland. BMC Vet Res. 2015;11:15. [PubMed ID: 25636375]. [PubMed Central ID: PMC4316766]. https://doi.org/10.1186/s12917-015-0332-x.
-
19.
Campioni F, Moratto Bergamini AM, Falcao JP. Genetic diversity, virulence genes and antimicrobial resistance of Salmonella Enteritidis isolated from food and humans over a 24-year period in Brazil. Food Microbiol. 2012;32(2):254-64. [PubMed ID: 22986188]. https://doi.org/10.1016/j.fm.2012.06.008.
-
20.
Bhowmick PP, Devegowda D, Ruwandeepika HAD, Karunasagar I, Karunasagar I. Presence of Salmonella pathogenicity island 2 genes in seafood-associated Salmonella serovars and the role of the sseC gene in survival of Salmonella enterica serovar Weltevreden in epithelial cells. Microbiol (Reading). 2011;157(Pt 1):160-8. [PubMed ID: 20884689]. https://doi.org/10.1099/mic.0.043596-0.
-
21.
Takaya A, Yamamoto T, Tokoyoda K. Humoral immunity vs. Salmonella. Front Immunol. 2019;10:3155. [PubMed ID: 32038650]. [PubMed Central ID: PMC6985548]. https://doi.org/10.3389/fimmu.2019.03155.
-
22.
Lan TTQ, Gaucher ML, Nhan NTM, Letellier A, Quessy S. Distribution of virulence genes among Salmonella serotypes isolated from pigs in Southern Vietnam. J Food Prot. 2018;81(9):1459-66. [PubMed ID: 30084656]. https://doi.org/10.4315/0362-028X.JFP-17-408.
-
23.
Zuo X, Peng G, Xia Y, Xu L, Zhao Q, Zhu Y, et al. A quadruple fluorescence quantitative PCR method for the identification of wild strains of african swine fever and gene-deficient strains. Virol J. 2023;20(1):150. [PubMed ID: 37452402]. [PubMed Central ID: PMC10347796]. https://doi.org/10.1186/s12985-023-02111-1.
-
24.
Qiao J, Alali WQ, Liu J, Wang Y, Chen S, Cui S, et al. Prevalence of virulence genes in extended-spectrum β-lactamases (ESBLs)- producing Salmonella in retail raw chicken in China. J Food Sci. 2018;83(4):1048-52. [PubMed ID: 29577300]. https://doi.org/10.1111/1750-3841.14111.