Patterns of antimicrobial resistance in Salmonella isolates from fattening pigs in Spain

Background

Swine are considered a major source of foodborne salmonellosis, a public health issue further complicated by the circulation of multidrug-resistant Salmonella strains that threaten the safety of the food chain. The current study aimed to identify patterns that can help to understand the epidemiology of antimicrobial resistance (AMR) in Salmonella in pigs in Spain through the application of several multivariate statistical methods to data from the AMR national surveillance programs from 2001 to 2017.

Results

A total of 1,318 pig Salmonella isolates belonging to 63 different serotypes were isolated and their AMR profiles were determined. Tetracycline resistance across provinces in Spain was the highest among all antimicrobials and ranged from 66.7% to 95.8%, followed by sulfamethoxazole resistance (range: 42.5% − 77.8%), streptomycin resistance (range: 45.7% − 76.7%), ampicillin resistance (range: 24.3% − 66.7%, with a lower percentage of resistance in the South-East of Spain), and chloramphenicol resistance (range: 8.5% − 41.1%). A significant increase in the percentage of resistant isolates to chloramphenicol, sulfamethoxazole, ampicillin and trimethoprim from 2013 to 2017 was observed. Bayesian network analysis showed the existence of dependencies between resistance to antimicrobials of the same but also different families, with chloramphenicol and sulfamethoxazole in the centre of the networks. In the networks, the conditional probability for an isolate susceptible to ciprofloxacin that was also susceptible to nalidixic acid was 0.999 but for an isolate resistant to ciprofloxacin that was also resistant to nalidixic acid was only 0.779. An isolate susceptible to florfenicol would be expected to be susceptible to chloramphenicol, whereas an isolate resistant to chloramphenicol had a conditional probability of being resistant to florfenicol at only 0.221. Hierarchical clustering further demonstrated the linkage between certain resistances (and serotypes). For example, a higher likelihood of multidrug-resistance in isolates belonging to 1,4,[5],12:i:- serotype was found, and in the cluster where all isolates were resistant to tetracycline, chloramphenicol and florfenicol, 86.9% (n = 53) of the isolates were Typhimurium.

Conclusion

Our study demonstrated the power of multivariate statistical methods in discovering trends and patterns of AMR and found the existence of serotype-specific AMR patterns for serotypes of public health concern in Salmonella isolates in pigs in Spain.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12917-022-03377-3.

Different Resolution Power of Multilocus Variable-Number Tandem Repeat Analysis and Whole-Genome Sequencing in the Characterization of S. 1,4,[5],12:i:- Isolates

Salmonella enterica serovar 1,4,[5],12:i:- has emerged over the last two decades as one of the most common serovars causing human salmonellosis in Europe. It is supposed to originate from Salmonella enterica serovar Typhimurium due to antigenic and genotypic similarities between the two serovars. Due to the high level of similarity, the multilocus variable-number tandem repeat analysis (MLVA) protocol designed for Salmonella Typhimurium routine typing is commonly used also for the characterization of S. 1,4,[5],12:i. Nevertheless, the Salmonella Typhimurium-based MLVA protocol often shows poor discriminatory power for S. 1,4,[5],12:i. Indeed, only a limited number of MLVA profiles have been described for S. 1,4,[5],12:i:-. Moreover, based on the MLVA clustering, S. 1,4,[5],12:i:- is supposed to display high clonality. The aim of the present work was to assess whether the five loci of Salmonella Typhimurium investigated by MLVA are sufficiently accurate to correctly assign S. 1,4,[5],12:i:- isolates. For this purpose, 38 epidemiologically unrelated S. 1,4,[5],12:i:- were subjected to whole-genome sequencing. Isolates were selected among a collection of monophasic strains isolated in Italy from different sources over the period 2014-2016 and belonging to the five most commonly detected MLVA profiles. Results confirmed the possible clonality for S. 1,4,[5],12:i:- serovar in the light of the scarce difference observed in terms of single-nucleotide polymorphisms (SNPs) among investigated isolates. Nevertheless, unrelated isolates on the basis of the difference of SNP number were characterized as indistinguishable by MLVA profile, thus suggesting an insufficient resolution of MLVA. Hence, we can conclude that MLVA-based approach does not seem a valuable proxy to deepen into the epidemiological relationship among S. 1,4,[5],12:i:- isolates. These evidences can be useful to avoid incorrect assignment especially when surveillance data are used for outbreak investigations.

A liquid bead array for the identification and characterization of fljB-positive and fljB-negative monophasic variants of Salmonella Typhimurium

Salmonella1,4,[5],12:i:- accounts currently for one of the most common serotypes observed worldwide. These isolates do not express the FljB flagellin and mostly derive from Salmonella Typhimurium. They are therefore termed Salmonella Typhimurium monophasic variants (STMV) and are considered of comparable public health risk. Since serological identification of the somatic and flagellar antigens of STMV is not sufficient to demonstrate relatedness with Salmonella Typhimurium, additional assays detecting genetic markers unique to Salmonella Typhimurium are required. In addition, identification of the mutations affecting expression of the flagellar gene fljB can be useful to support the monophasic character observed phenotypically. Finally, genetic subtyping of the various mono- and biphasic Salmonella Typhimurium clonal groups can facilitate their epidemiological follow-up. Here, we present a home-made liquid bead array able to fulfill these requirements. This array confirmed the monophasic character of 240 STMV isolates collected in Belgium during 2014-2015 and identified 10 genetic subtypes. Microevolution in and around the fljB locus linked to IS26 insertions is probably one of the driven force accounting for STMV population diversity. Thanks to its open design, other genetic signatures could later be merged to the assay to subtype additional STMV clonal groups and to detect rare mutations.