Open access
Environmental Microbiology
Research Article
26 June 2024

A multicenter genomic epidemiological investigation in Brazil, Chile, and Mexico reveals the diversity and persistence of Salmonella populations in surface waters

ABSTRACT

This study examined the diversity and persistence of Salmonella in the surface waters of agricultural regions of Brazil, Chile, and Mexico. Research groups (three in 2019–2020 and five in 2021–2022) conducted a long-term survey of surface water across 5–8 months annually (n = 30 monthly). On-site, each team filtered 10-L water samples with modified Moore Swabs to capture Salmonella, which were then isolated and identified using conventional microbiological techniques. Salmonella isolates were sequenced on Illumina platforms. Salmonella was present in 1,493/3,291 water samples (45.8%), with varying isolation rates across countries and years. Newport, Infantis, and Typhimurium were the most frequent among the 128 different serovars. Notably, 22 serovars were found in all three countries, representing almost half of the 1,911 different isolates collected. The resistome comprised 72 antimicrobial resistance (AMR) genes and six point mutations in three genes. At least one AMR determinant was observed in 33.8% (646/1,911) of the isolates, of which 47.4% (306/646) were potentially multidrug resistant. Phylogeny based on core genome multilocus sequence typing (cgMLST) showed that most isolates clustered according to sequence type and country of origin. Only 14 cgMLST multi-country clusters were detected among the 275 clusters. However, further analysis confirmed that close genetic relatedness occurred mostly among isolates from the same country, with three exceptions. Interestingly, isolates closely related phylogenetically were recovered over multiple years within the same country, indicating the persistence of certain Salmonella in those areas. In conclusion, surface waters in these regions are consistently contaminated with diverse Salmonella, including strains that persist over time.

IMPORTANCE

Salmonella is a leading foodborne pathogen responsible for millions of illnesses, hospitalizations, and deaths annually. Although Salmonella-contaminated water has now been recognized as an important contamination source in the agrifood chain, there is a lack of knowledge on the global occurrence and diversity of Salmonella in surface water. Moreover, there has been insufficient research on Salmonella in surface waters from Latin American countries that are major producers and exporters of agricultural products. Incorporating genetic profiling of Salmonella isolates from underrepresented regions, such as Latin America, enhances our understanding of the pathogen’s ecology, evolution, antimicrobial resistance, and pathogenicity. Moreover, leveraging genomic data derived from pathogens isolated from diverse geographical areas is critical for assessing the potential public health risk posed by the pathogen and expediting investigations of foodborne outbreaks. Ultimately, global efforts contribute significantly to reducing the incidence of foodborne infections.

INTRODUCTION

Salmonella is a major global health concern, responsible for over 95 million cases and more than 50,000 deaths annually, predominantly in developing countries (1). The pathogen causes disease in a wide range of hosts, survives in various environments for extended periods, and is mainly transmitted to humans through contaminated food and water (2), leading to multiple salmonellosis outbreaks each year (24). While traditionally associated with meat and poultry, recent outbreaks increasingly implicate produce, such as onions, alfalfa sprouts, peaches, and papayas (5). One of the primary sources of produce contamination is irrigation water (69).
In Latin America, as in other regions, surface waters such as rivers, canals, dams, and lakes are used extensively for crop irrigation. Moreover, studies have reported these water sources are vulnerable to contamination with pathogens and other hazards (9, 10). These water sources are critical in understanding the interconnected health of humans, animals, and the environment, as per the One Health approach. They serve as a reliable indicator of regional health and potential reservoirs of various pathogens, including Salmonella.
Global studies have revealed widespread Salmonella contamination in surface waters (1115). A meta-analysis examining published reports from several countries estimates the presence of Salmonella in surface water to be 31% (16). Recently, high rates of Salmonella contamination in the surface waters of California (56.4%) and Pennsylvania (49%) have also been reported (17). However, not all Salmonella serovars are epidemiologically relevant, and some seem to pose higher risks to humans (2). Studies have identified both high-risk serovars, such as Newport, Typhimurium, and Enteritidis, and less common types in surface water (11, 1720).
Research specific to surface waters in Latin America is limited, with existing studies often focusing on small areas or single countries and are of a cross-sectional nature (18, 19, 21, 22). Comprehensive and longitudinal research across the region could reveal the diversity and commonality of Salmonella serovars, their resistome, and their ability to persist and disseminate across different regions in Latin America. This investigation aimed to understand the broader picture of Salmonella’s population in agricultural waters in different regions of Latin America to obtain further insights into its role in food safety and its global diversity.

MATERIALS AND METHODS

Samples and isolates

Water samples were obtained in Mexico (n = 689) and Chile (n = 2,036) from 2019 to 2022. Additionally, through 2021 and 2022, 542 samples were obtained in Brazil (Table 1; Fig. 1). Reports for Salmonella isolation rates in 2019 and 2020 were published by Toro et al. for Chile and Ballesteros-Nova et al. for Mexico (21, 22). The sampling procedure for four of the teams involved on-site filtration of 10 L of surface water (rivers, canals, streams, ponds, and dams, among others) through a modified Moore swab as described by Sbodio et al. (23), using a peristaltic pump. Samples from Paraíba State in Brazil were taken in triplicates, filtering a total of 30 L per sample. Samples were transported on ice to one of the five participant local teams’ laboratories: Laboratorio de Microbiología y Probióticos, Universidad de Chile, Santiago, Chile; Laboratorio de Inocuidad Veterinaria, Pontificia Universidad Católica de Chile, Santiago, Chile; Laboratorio de Epidemiologia Veterinaria, Universidad Autónoma de México, Mexico City, Mexico; Laboratório de Avaliação de Produtos de Origem Animal, Universidade Federal da Paraíba, Areia, Paraíba, Brazil; Laboratório de Investigação em Microbiologia Médica, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil. Samples were stored in refrigeration until processing within 7 days. The samples were processed following the U.S. Food and Drug Administration (FDA)’s Bacteriological Analytical Manual methodology for Salmonella spp. with minor modifications (https://www.fda.gov/food/laboratory-methods-food/bam-chapter-5-salmonella), as described by Toro et al. (21). Presumptive colonies (1–10 per sample) were confirmed as Salmonella through invA PCR with primers F 5′ GAATCCTCAGTTTTTCAACGTTTC 3′ and R 5′ TAGCCGTAACAACCAATACAAATG 3′ and kept in brain heart infusion (BHI) broth with 20% glycerol at −80°C for further analysis (24).
TABLE 1
TABLE 1 Salmonella isolation from surface water by country and yeara
CountryYearTotal
2019202020212022
Chile20.0% (120/600)38.1% (183/480)34.4% (165/480)43.4% (208/479)33.1% (676/2,039)
Mexico62.6% (122/195)75.1% (139/185)57.8% (67/116)55.0% (110/200)62.9% (438/696)
Brazil0065.0% (206/317)72.4% (173/239)68.2% (379/556)
Total per year30.4% (242/795)48.4% (322/665)48.0% (438/913)53.5% (491/918)45.4% (1,493/3,291)
a
Sample sizes varied by country and year, and samples were collected in different seasons. Brazil joined the project in 2021, so there are no results for 2019 or 2020 for Brazil. Values in parentheses indicate samples in which Salmonella was isolated/total number of samples.
Fig 1
A map of Latin America highlights the regions of Central Mexico, Mexico; Paraíba State, Brazil; Rio de Janeiro State, Brazil; Santiago Metropolitan Region, Chile; and Region del Maule, Chile. Insets detail specific areas within each region.
Fig 1 Sampling sites of surface waters in Latin America from 2019 to 2022: (a) Central Mexico, Mexico, (b) Paraíba, Brazil, (c) Rio de Janeiro, Brazil, (d) Santiago Metropolitan Region, Chile, and (e) Region del Maule, Chile. Mapping performed on Microrreact Version 252 (Centre for Genomic Pathogen Surveillance) (25).

Whole-genome sequencing, assembly, and quality control

Over 4,255 Salmonella isolates from 1,169 samples were sequenced at the Center for Food Safety and Applied Nutrition (CFSAN/FDA). All isolates were confirmed as Salmonella with VitekMS (Biomerieux, Marcy-l'Étoile, France) before DNA extraction. Then, individual colonies were grown overnight in TSB broth (BD, Sparks, MD, USA) at 37°C. DNA was purified on a Maxwell RSC-48 Instrument using the Cultured Cells DNA Kit (Promega, Madison, WI, USA). Libraries were prepared with the Illumina DNA Prep kit (Illumina, San Diego, CA, USA) on the Sciclone G3 NGSx iQ Workstation (Perkin Elmer, MA, USA), and sequencing was performed on the Illumina NextSeq 2000 using the NextSeq 1000/2000 P2 reagents 300 cycles with the 150 paired-end chemistry (Illumina).
Read quality control procedures were performed following the CFSAN MicroRunQC pipeline in GalaxyTrakr at CFSAN, FDA (26). The pipeline includes Trimmomatic for read trimming (27), skesa (V2.4.0) for assembly (28), mlst (V2.23.0) to identify species and determine sequence type (29, 30), and fastq-scan (V1.0.1) for assembly statistics (31). Genomes that passed quality control (32) were uploaded as SRA files to the NCBI (Dataset S1) in batches and made publicly available through the NCBI Pathogen Detection browser (PD). NCBI-curated assembled genomes, public metadata, assembly metrics, and antimicrobial resistance (AMR) gene results from the PD browser were downloaded for the current data analysis (https://www.ncbi.nlm.nih.gov/pathogens/). Metadata were curated using each participant’s laboratory database. No genomes were available from 324 positive samples due to contamination, unavailability, low sequencing quality, and other unforeseen circumstances. Genomic data from isolates collected in Mexico during 2019 and 2020 were published in Ballesteros-Nova et al. (22).

Genomic characterization of Salmonella enterica isolates (MLST, in silico serotyping, and AMR detection)

Achtman multilocus sequence typing (MLST) types with genes aroC, dnaN, hemD, hisD, purE, sucA, and thrA were used for Salmonella MLST determination with the Ridom Seqsphere software (Ridom GmbH, Germany) (33). For serovar definition, we downloaded the serovar prediction from PD and analyzed each selected genome with SISTR (34) in Ridom. Disagreements between PD and SISTR’s serovar predictions were resolved with manual curation based on PD clustering. AMRFinderPlus results from the PD browser were used to identify the presence of AMR genes and assign AMR class to each AMR gene detected as NCBI definition for every gene (35). Genomes carrying AMR determinants against three or more antimicrobial classes were identified as multidrug-resistant. Kruskal-Wallis and Dunn’s post hoc (adjusted for multiple comparisons) tests were run in Past v4.03 (36) to investigate the potential differences in the number of AMR genes between isolates from different countries.

Phylogenomic relatedness among Salmonella isolates

In the first phase, we cleaned our database by removing clonal isolates from the same sample (i.e., isolates obtained from the same day and location). First, available genomes (n = 4,255) were retrieved from the PD browser. Then, we examined each single nucleotide polymorphism (SNP) cluster comprising genomes from the project and determined the SNP count for all genomes from the same sample. Clones were defined as genomes with 20 or fewer SNP and monophyletic relationships (37). The genome with the greatest quality (defined by N50 and the number of contigs) was chosen to represent the clonal group in each sample. A total of 1,911 genomes were selected for further studies.
Then, we analyzed the relatedness among selected genomes (n = 1,911) through a core genome multilocus sequence typing (cgMLST) scheme in Ridom SeqSphere v9.0.2 (Ridom GmbH, Germany), which comprises 3,002 core genes for Salmonella species (33). A Minimum Spanning Tree (MST) was crafted from the resulting distance matrix. We removed 1,887 columns with missing values in at least one sample and excluded one genome that was missing more than 10% of alleles in distance columns (CLA-161-1, Chile) for a total of 1,910 genomes used for this analysis. Additionally, we examined the phylogenetic relatedness among genomes from the same serovar and different countries that resulted in mixed country clusters in the previous MST using the same methodology.

Worldwide genomic comparison

We investigated the relatedness among the genomes in our collection and genomes from isolates collected elsewhere as per PD. For that, we searched all selected genomes in the PD cluster and downloaded the SNP cluster information into a spreadsheet (22 January 2024). Each PD cluster was examined, and clustering with genomes of other origins was examined. Information was added to each record, such as the country where the genome came from, serovar, the presence of genomes from clinical isolates in the cluster, and other notes associated with the origin of the isolates sequenced.

RESULTS

Salmonella in surface water

Our study revealed a notable occurrence of Salmonella in surface water across three Latin American countries, with an overall isolation rate of 45.4% (1,493/3,291) over the period 2019–2022 (Table 1). The isolation rates varied significantly by country, with Brazil reporting the highest at 68.2%, followed by Mexico at 62.9% and Chile at 33.1% (Table 1). In samples testing positive for Salmonella, 1–10 isolates per sample underwent genomic sequencing, characterization, and phylogenomic analysis.

Genomic diversity and analysis of Salmonella from surface water

The number of non-clonal isolates per sample ranged from 1–6, and two samples carried 10 (Dataset S1; Table S1). Most samples (58.4% or 683 out of 1,169 examined) contained only a single clonal group (Fig. 2; Table S2). Notably, the samples from Brazil carried an average of 2.01 non-clonal isolates per sample, and two samples carried 10 different isolates representing 10 different serovars (Fig. 2; Table S2). Meanwhile, the average for Mexican samples was 1.77, while the average for non-clonal isolates from Chile was 1.34. Statistical analysis revealed that Chile carried the least number of non-clonal isolates per sample (P < 0.001), but no differences were found between the number of isolates per sample collected in Brazil and Mexico. In general, non-clonal isolates from the same sample were of different serovars (Dataset S1). Mexico (n = 686) contributed the most genomes to the collection, followed by Chile (n = 684) and Brazil (n = 541) (Tables S1 and S2).
Fig 2
A bar graph plots the percentage of samples with the number of isolates per sample from Brazil, Chile, and Mexico. The highest percentages are for 1 isolate per sample, especially in Chile, followed by Mexico and Brazil.
Fig 2 Non-clonal Salmonella isolates per sample (%). Clonality was examined in the Pathogen Detection browser (NCBI). Clonal isolates had less than 20 SNP differences and a monophyletic clustering. Genomes available for analysis were from 269 samples from Brazil, 513 from Chile, and 387 from Mexico, for a total of 1,169 samples.

Sequence types and serovars of Salmonella from surface water

We identified 201 distinct sequence types (STs) among 1,911 representative isolates, but 12 were not attributed to any ST due to the incomplete sequence of one gene required for complete typing. The most prevalent STs were ST32 (7.4%; Salmonella Infantis), ST19 (4.8%; Salmonella Typhimurium), and ST13 (4.6%; Salmonella Agona). Chile exhibited the least diversity with 86 STs, while Mexico and Brazil had slightly higher diversity with 91 and 98 STs, respectively. A notable finding was that only 17 STs were present in all three countries with most STs detected in a single country (147/201; Table S3).
Similarly, the sequenced isolates were classified into 128 serovars based on the genomics analysis. The predominant serovars across the study were Newport (9.6%), Infantis (7.7%), and Typhimurium (6.7%) (Table 2). The same three serovars were the most common for the Chilean genomes, while Brazil and Mexico showed different top three serovars (Newport, Saintpaul, and Rubislaw for Brazil; Newport, Senftenberg, and Anatum for Mexico) (Table 2). Most serovars (77 out of 128) were unique to a single country, with only 22 serovars shared among all three nations (Fig. 3). Interestingly, Salmonella Enteritidis was identified predominantly in Chile (47/49 S. Enteritidis isolates) but was absent in samples from Brazil (Table 2).
TABLE 2
TABLE 2 Number of genomes of each Salmonella serovar in surface water samples in Latin America (n = 1,911)
SerovarBrazilChileMexicoTotalNumber of countries
Abaetetuba1 452
Adelaide1528343
Agona93939873
Albany4119243
Altona  111
Anatum5843563
Bareilly  111
Bovismorbificans256133
Braenderup19622473
Brandenburg2202243
Bredeney  16161
Bulbay1  11
Bullbay1  11
Businga1  11
Carrau19  191
Cerro2127213
Corvallis19711373
Derby 317202
Dublin 2 21
Edinburg 25 251
Enteritidis 472492
Freetown1  11
Fresno 1232
Gaminara5  51
Give 225272
Glostrup  331
Goldcoast 5 51
Grumpensis1  11
Hadar1012133
Havana  661
Heidelberg3  31
I 1,3,19:c:-1  11
I 1,4,[5],12:b:-3  31
I 1,4,[5],12:d:-27 92
I 1,4,[5],12:i:-286163
I 16:r:e,n,z151  11
I 18:d:-3  31
I 3,10:d:-1  11
I 4:-:1,51  11
I 7:-:1,52  21
I 7:k:-1  11
I 7:l,v:-3  31
I 7:z36,z38:-1  11
I K:y:1,5  111
II 42:r:-5  51
IIIb 16:z10:e,n,x,z152 21
IIIb 18:i:z 2 21
IIIb 18:k:z 5 51
IIIb 18:z10:e,n,x,z15 1 11
IIIb 35:r:z1  11
IIIb 38:(k):z351  11
IIIb 48:i:z 15 151
IIIb 50:r:z26 82
IIIb 58:k:z 15 151
IIIb 60:r:e,n,x,z15 331
IIIb 61:i:z 7 71
IIIb 61:l,v:z1  11
IIIb 65:(k):z 13 131
IIIb P:k:z351  11
Infantis3488251473
Inganda1  11
Isangi  331
IV [1],53:g,z51:-1  11
IV 16:z4,z32:-2  21
IV 18:m,t:-1  11
IV 21:z4,z23:-2  21
IV 38:g,z51:-1  11
IV 40:z4,z24:- 1 11
IV 43:z4,z23:-44 82
IV 43:z4,z24:-11  111
IV 45:g,z51:-1  11
IV 48:g,z51:-1  11
IV 50:z4,z23:-1  11
IV 6,7:z4,z24:-1  11
IV R:z4,z24:- 1 11
Javiana1912223
Johannesburg 4 41
Jos3  31
Kedougou 2 21
Kentucky11793
Kiambu1 11122
Litchfield  111
Liverpool  221
Livingstone 211222
Lomita2  21
London  21211
Madelia11  111
Manhattan 3142
Mbandaka5813263
Meleagridis1 22232
Miami1  11
Michigan2  21
Minnesota4 592
Mississippi  331
Molade3 252
Montevideo 1892
Muenchen13531493
Muenster3 11142
Newport5750761833
Ohio1 672
Oran13  131
Oranienburg5813263
Oslo4 6102
Othmarschen1  11
Panama341623733
Paratyphi B 7 71
Pomona1 452
Poona11 3142
Reading  331
Rhydyfelin4  41
Rissen 1452
Rubislaw36 5412
Saintpaul46 7532
Sandiego168 242
Santiago or Belem816 242
Saphra7  71
Schwarzengrund525123
Senftenberg 2452762
Soerenga17193
Stanley 34 341
Tennessee 3142
Thompson 181192
Tucson1  11
Typhimurium1879311283
Uganda  551
Urbana4 482
Weltevreden  111
Worthington 1 11
Total5416846861,911 
Fig 3
A Venn diagram depicts the overlap of isolates among four countries. Brazil has 44 unique isolates, Chile has 18, and Mexico has 15. Overlaps include 22 shared by all three countries: Brazil and Mexico: 13, Chile and Mexico: 11, Chile and Brazil: 5.
Fig 3 Number of serovars per country (n = 128 serovars).

AMR determinants in Salmonella from surface water

Our analysis detected a total of 78 AMR elements; 72 of those were acquired antimicrobial resistance genes (ARGs) and 6 were point mutations in genes gyrA (D87G, D87Y, S83F, and S83Y), parE (H462Y), and ramR (T18P) (Dataset S2). These AMR elements suggest resistance against a broad range of antimicrobial classes, including aminoglycosides, beta-lactams, fosfomycin, macrolides, phenicols, quinolones, sulfonamides, tetracyclines, and trimethoprim (Table S4). A small subset of genomes (2.5%; 47/1,911) carried ARGs associated with resistance to less common antimicrobial classes such as bleomycin, colistin, lincosamides, streptothricin, and rifamycin (Fig. 4). Overall, one-third of genomes (33.9%; 647/1,911) carried 1–17 AMR elements. However, almost half of them (47.8%; 309/647) contained a single AMR element (Dataset S2; Table S5). A similar percentage (47.3%; 306/647) was predicted to be multidrug-resistant, carrying AMR elements against three or more AMR classes (Table 3). Interestingly, 69 genomes carried 10 or more AMR determinants simultaneously (Table S5). For instance, one Salmonella Senftenberg (MPSPSA19321-1) carried 14 AMR elements and 4 Salmonella Derby isolates (MPSPSA2067-1, MPSPSA2068-2, MPSPSA2094-1, and MPSPSA2149-1) carried 13 of them; these were predicted to be resistant to nine AMR families (Dataset S2).
Fig 4
A bar graph compares antimicrobial classes versus genome percentage across Brazil, Chile, Mexico, and the total. Quinolones and tetracyclines have the highest percentages overall, with Mexico depicting the highest values in several categories.
Fig 4 Genomes (%) with predicted antimicrobial resistance to 14 antimicrobial classes by country. The number of available genomes per country: Brazil n = 541, Chile n = 685, and Mexico n = 686.
TABLE 3
TABLE 3 Number of antimicrobial resistances predicted by isolate and countrya
Number of predicted AMR resistances/isolateNo. (%) of isolates
BrazilChileMexicoTotal
90 (0.0)0 (0)5 (0.7)5 (0.3)
80 (0.0)7 (1.0)45 (6.6)52 (2.7)
71 (0.2)35 (5.1)20 (2.9)56 (2.9)
62 (0.4)18 (2.6)29 (4.2)49 (2.6)
59 (1.7)14 (2.0)34 (5.0)57 (3.0)
44 (0.7)18 (2.6)23 (3.4)45 (2.4)
314 (2.6)3 (0.4)25 (3.6)42 (2.2)
28 (1.5)8 (1.2)17 (2.5)33 (1.7)
168 (12.6)84 (12.3)155 (22.6)307 (16.1)
0435 (80.4)497 (72.7)333 (48.5)1,265 (66.2)
Total isolates541 (100)684 (100)686 (100)1,911 (100)
a
Prediction was based on the presence of AMR elements (genes and mutations) in each genome and phenotype associated with each element in the NCBI AMRFinderPlus database.
Quinolone resistance determinants were the most frequently detected in each country (22.8%) (Fig. 4). The gene qnrB19, linked to reduced susceptibility to quinolones and previously described in plasmids, was the most prevalent ARG (14.3%) (Table S6). Additionally, four different point mutations in the gyrA gene, often associated with quinolone resistance, were detected in 5.8% of the genomes (110/1,911; Table S6). Tetracycline resistance was predicted in 15.8% (301/1,911) of the genomes (Fig. 4), with tet(A) as the second most frequent ARG (13.7%; Table S6). The resistance to this antimicrobial class was among the three most frequent in the three countries. The floR gene conferring resistance to phenicols ranked third in the prevalence (10.5%), although phenicol resistance was only the fifth most common in the AMR resistance prediction (Fig. 4; Table S6). Moreover, 17 different genes linked to aminoglycoside-modifying enzymes (AME) were identified, and 129 genomes carried simultaneously three or more AME genes (Table S6). Resistance to aminoglycosides was among the four most important in each country. Notably, the mobile gene mcr-9.1, a newly described colistin resistance gene, was found in three genomes (two S. Agona from Mexico and one from S. 1,4,[5],12:d:- from Chile) (Dataset S2; Fig. 4). Finally, resistance to fosfomycin was the third most frequently detected in Salmonella from Brazil (5.0% of Brazilian genomes), but it only reached seventh place in Chile and eighth in Mexico (Fig. 4).
Data analysis revealed the genomes from Mexico generally carried the highest number of AMR determinants and had the most frequent occurrence of ARG per antimicrobial class (P < 0.001) (Fig. 4; Table S6). For instance, 35.3% of genomes from Mexico carried AMR determinants against quinolones, compared to only 16.1% from Chile and 15.4% from Brazil (Fig. 4). This trend was consistent across all analyzed antimicrobial classes. Moreover, Mexican genomes collectively carried AMR determinants against 14 antimicrobial classes, whereas the Chilean and Brazilian genomes had AMR determinants against 13 and 10 families, respectively (Dataset S2).

Genomic diversity of Salmonella from surface water through core genome MLST analysis

The cgMLST analysis grouped 1,910 Salmonella genomes into 275 clusters with ≤7 gene differences (Fig. 5). Out of the 3,002 genes in the cgMLST scheme for Salmonella, only 1,115 were ubiquitously present across all genomes and were used for distance calculations. The largest cluster (Cluster 1) consisted of 61 S. Infantis isolates from Chile, 6 from Mexico, and 1 from Brazil. The next two largest clusters comprised S. Newport (Cluster 2) and S. Enteritidis (Cluster 3) from Chile. The fourth largest cluster (Cluster 4) included 29 S. Newport from Mexico collected between 2019 and 2021 and 1 from Chile in 2021. The most numerous clusters formed by Brazilian genomes included 20 Salmonella Saintpaul (Cluster 11) collected in 2021 and 2022 (Fig. 5).
Fig 5
A network of various clusters connected by lines indicates the interactions. Specific clusters are labeled, including Cluster 1, Cluster 2, Cluster 3, Cluster 9, Cluster 11, Cluster 12, Cluster 18, Cluster 39, and Cluster 46. Nodes are color-coded.
Fig 5 Minimum spanning tree (MST) of 1,910 Salmonella isolates obtained from the surface waters in Latin America. The MST was crafted with 1,115 genes present in all isolates. Each circle represents one isolate or a group of clonal isolates. Colors represent isolation countries: Brazil in red, Chile in green, and Mexico in purple. Gray background linking two or more circles indicates clonal groups of seven or fewer core genome gene difference. The number next to the lines between circles indicates cg gene differences. Clusters 1, 2, 3, 4, 5, 11, 12, 38, and 65 mentioned in the Results section are highlighted with a larger cluster name and a red circle around the genomes.
Although most clusters with seven or fewer allele differences were composed of genomes of the same serovar and country, we detected 17 multi-country clusters (6.2%) (Fig. 5; Table S7). Specifically, only five clusters contained related genomes from all three countries, representing serovars Agona (n = 28; Cluster 5), Corvallis (n = 20; Cluster 12), Infantis (n = 68; Cluster 1), Mbandaka (n = 5; Cluster 65), and Soerenga (n = 9; Cluster 38), clusters that collectively accounted for 130 genomes. Secondary cgMLST-specific analyses including only genomes from individual serovars revealed only three cases of closely related isolates from different countries: Genomes of S. Enteritidis from Chile and Mexico (6-MAP-28-3 and MPSPSA2069-5; seven allele difference in Cluster 6), S. Infantis from Brazil and Chile (4C1011TX2 and 6-MAP-07B-2; 0 allele difference in Cluster 19), and S. Newport from Chile and Mexico (6-MAP-26-3 and MPSPSA2012-6/2014-4; two allele difference in Cluster 2) (Fig. 6a through c). We did not find close phylogenetic relatedness among the isolates from the same serovar and different countries in other clusters (Fig. 7a through d; Fig. S1a through e).
Fig 6
A network diagram of Enteritidis, Infantis, and Newport. Nodes are color-coded by country: Brazil, Chile, and Mexico. Each part depicts clusters of nodes indicating the interactions, with some clusters highlighted in dark boxes.
Fig 6 Minimum spanning tree (MST) of phylogenetically related isolates from a single serovar isolated in the three countries. MST generated from core genome MLST of individual serovars clustering together from different countries in the general analysis cgMLST. Each circle represents one genome or a group of clonal genomes. Each circle color represents isolation countries: Brazil in red, Chile in green, and Mexico in purple. Gray background linking two or more circles indicates clonal groups of seven or fewer core genome gene difference. The number next to the lines between circles indicates cg gene differences. The number of isolates and genes considered for distance calculation are as follows: (a) Salmonella Enteritidis genomes (n = 49) had 2,826 common genes. (b) The Salmonella Infantis collection (n = 147) had 2,656 common genes. (c) Salmonella Newport genomes (n = 183) had 2,509 common genes. Dark boxes highlight clusters comnprised of genomes from different countries.
Fig 7
A network diagram of Agona, Corvallis, Mbandaka, and Soerenga. Nodes are color-coded by country: Brazil, Chile, and Mexico. Each part depicts clusters of nodes indicating the interactions, with some clusters highlighted.
Fig 7 Minimum spanning tree (MST) of non-closed phylogenetically related genomes from a single serovar isolated in the three countries. (a) Agona, (b) Corvallis, (c) Mbandaka, and (d) Soerenga. MSTs generated from core genome MSLT of genomes of individual serovars clustering together from different countries in the general analysis cgMLST. Each circle represents one genome or a group of clonal genomes. Each circle color represents isolation countries: Brazil in red, Chile in green, and Mexico in purple. Gray background linking two or more circles indicates clonal groups of seven or fewer core genome gene difference. The number next to the lines between circles indicates cg gene differences. The number of isolates and genes considered for distance calculation are as follows: (a) Salmonella Agona (n = 87) had 2,786 common genes. (b) Salmonella Corvallis genomes (n = 37) had 2,872 common genes. (c) The Salmonella Mbandaka (n = 26) group had 2,616 genes in common. (d) Salmonella Soerenga (n = 9) had 2,896 common genes.
Interestingly, we observed that closely related genomes from the same country were isolated over long periods. For example, S. Braenderup from Mexico isolated from 2019 to 2021 formed a cluster of nine genomes and S. Enteritidis from Chile (n = 12) were isolated from 2019 to 2021 (Fig. 6a; Fig. S1b). S. Infantis from Chile (Cluster 1; Fig. 5) were collected throughout every year from 2019 to 2022 and formed a cluster of more than 40 highly related genomes (Fig. 6b). Moreover, subclusters were formed by undistinguishable S. Infantis from Chile collected over 3 (April–June) and 6 (July through December) months in 2019, and some of them were isolated in the same sampling site (Fig. 5 and 6b).

Genomic comparison of Salmonella from surface waters with global Salmonella collections

In a broader context, the 1,911 genomes were grouped into 613 clusters according to the Pathogen Detection Browser (as of 22 January 2024). These PD clusters comprised 2–22,924 genomes and represented 115 different serovars. Most PD clusters (57.8%; 354/613) included only isolates collected in this study and were country specific. In total, 240 PD clusters contained clinical isolates (39.0%) (Table S8). Among those, nine clusters showed clinical and water isolates with 0 SNP differences (PDS PDS000026867.62, PDS000007781.1003, PDS000031814.184, PDS000176795.35, PDS000001852.755, PDS000027266.25, PDS000042357.49, PDS000013845.430, and PDS000123320.1; cluster information updated on 12 April 2024). We identified 13 PD clusters (2.1%) of 11 serovars that contained genomes from two or three countries in the study. Most of these results coincide with the cgMLST results, and PD clusters of serovars Agona, Corvallis, Enteritidis, Infantis, Kiambu, Mbandaka, Newport, Senftenberg, Soerenga, and Typhimurium contained isolates from two or three countries in the study. For instance, PD Cluster PDS000007781.892 contained genomes from serovar Newport from Chile and Mexico (6-MAP-26-3 and MPSPSA2057-1) that were separated by three SNPs in PD and were in the same cluster in the S. Newport cgMLST tree (Fig. 6c; Fig. S2).
Furthermore, we identified multiple single-country clusters produced by a combination of genomes from water and other sources, such as poultry (PDS000121961.2, Salmonella Corvallis from Chile; PDS000113904.5, Salmonella Molade or Salmonella Wipra from Brazil; and PDS000020042.7 S. Infantis from Brazil), pork (PDS000123296.2, Salmonella Anatum from Mexico), bovine (PDS000051281.2, Salmonella Fresno from Mexico), horse (PDS000053294.12, S. Typhimurium from Chile), and vegetable origin foods (PDS000031203.7, Salmonella Meleagridis from Mexico; PDS000047156.4, S. Anatum from Mexico; and PDS000075580.4, Salmonella Sandiego from Chile).

DISCUSSION

In this study, we described the presence of a large diversity of Salmonella serovars in surface waters of Latin America, with genomes genetically closely related within each country but more distant or unrelated when comparing among countries.
First, Salmonella isolation rates across different regions varied. This was expected because of the regional differences in ecosystem characteristics and climatic conditions, which could affect Salmonella spread and survival in the environment. For instance, sampled regions in Central Chile have a temperate Mediterranean climate; Central Mexico has arid climates; and Brazil has tropical climates that are usually humid. This variability highlights the complexity of Salmonella’s ecological presence, affected by factors such as temperature and rainfall patterns (38). The research teams from Chile, using data from the first 2 years of the project, already explored environmental and anthropogenic factors that could influence the likelihood of isolating Salmonella locally. They found that the sampling month was the factor with more impact on Salmonella isolation and that the anthropogenic variables evaluated only marginally contributed to this result (21). The other teams in the project are currently exploring local factors influencing the phenomena.
Co-occurrence of multiple Salmonella serovars was observed in almost 40% of samples. Samples from Brazil contained the highest number of serovars per sample (n = 10). This can be explained by the fact that the group in Paraiba State, Brazil, used a triplicate-sampling approach, resulting in the filtering of 30 L of water per sample compared with the 10-L samples used by the other groups, which used uniplicate sampling (10 L). In the remaining groups, we observed up to six isolates per sample. This could also have affected the isolation rate comparability among groups as a larger sample volume resulted in higher detection rates and greater diversity across samples. Other studies have described the presence of more than one Salmonella serovar in a single sample. For instance, a study in the Susquehanna River watershed in Pennsylvania reported the presence of up to 10 different Salmonella serovars per sample, with an average of three (17), which is higher than the average of 1.6 from this study. However, Deaven et al. (17) used CRISPR-SeroSeq, an amplicon-based sequencing tool that allows for all Salmonella present in a sample to be identified. In our study, we picked colonies from plates and sequenced up to 10 isolates per sample. Since one serovar might be more abundant than others by several orders of magnitude, picking colonies representing different populations becomes challenging, and it might underestimate the actual Salmonella diversity in our samples. This has been observed in other matrices, such as chicken carcasses and droppings (39, 40). Consequently, our diversity per sample could be even more substantial than reported.
The frequent isolation of S. Newport, S. Typhimurium, and S. Infantis in our samples is noteworthy, given their high relevance in foodborne salmonellosis outbreaks and the emergence of antimicrobial-resistant strains belonging to these serovars. S. Newport has been frequently implicated in foodborne outbreaks linked to produce in the U.S.: alfalfa sprouts in 2010, cantaloupe in 2012, cucumbers in 2014, papayas in 2017, and onions in 2020 (4143). S. Newport has also been frequently isolated from tomatoes grown on the East Coast in the U.S. (44). This serovar has been reported among the top five serovars causing human clinical cases in Brazil and Mexico for several years (45, 46). Accordingly, this was the most common serovar found in surface water samples from both countries. Although S. Newport is not among the top five most prevalent human serovars in the last Chilean national report (47), it was among the most frequently identified serovars from Chile. S. Typhimurium has been reported as the first or second most frequent cause of salmonellosis in the three participant countries (4548), and it has been isolated from several other sources in Latin America, such as domestic animals (49), food of animal origin (46), and water (18, 19). This serovar has also caused foodborne outbreaks in the U.S. linked to tomatoes in 2006, cantaloupes in 2012, and alfalfa sprouts in 2022 (5052). S. Typhimurium high frequency in surface water indicates that this serovar is widespread in the environments surveyed. Interestingly, S. Typhimurium was detected in only one of the two participating groups in Brazil, probably indicating potential risk factors associated with this serovar. S. Infantis is an emerging serovar worldwide, characterized by the presence of a plasmid carrying multiple AMR genes (53). In Chile, it was the second leading cause of salmonellosis in 2018 (47), and multi-resistant isolates have been isolated from poultry facilities (32), chicken meat (33), pet foods (34), and water (19), among other sources. In Brazil, it has also been isolated from diverse sources, including the poultry industry (32), and it has been among the top 10 clinical human serovars since 2013 (45). In Mexico, the presence of this serovar has been associated with beef (46). Moreover, this serovar was linked to an outbreak involving papayas in the U.S. in 2017 (54). Altogether, the presence of these three serovars among the most frequently isolated across the whole study indicates that high-risk salmonellae contaminate surface waters used for irrigation in Latin America and are widespread in the surveyed areas, potentially becoming a risk to consumers.
Other clinically relevant serovars were found in this study. For instance, serovars that have caused outbreaks linked to produce, such as Agona 4.6% (n = 87), Anatum 2.9% (n = 56), Panama 3.8% (n = 73), Muenchen 2.6% (n = 49), Braenderup 2.5% (n = 42), Adelaide 1.8% (n = 34), Oranienburg 1.4% (n = 26), and Javiana 1.2% (n = 22), were present in all three countries in variable frequencies (Table 2). On the other hand, other relevant serovars were not collected from all three countries. For example, Salmonella Enteritidis has been the top serovar linked to human disease in Chile for several years, and it represents a top disease-causing Salmonella in Brazil and Mexico (4547). Interestingly, almost all S. Enteritidis isolates were collected in Chile, and none came from Brazil. Salmonella Heidelberg has been consistently reported as a top serovar causing salmonellosis in Brazil (45, 48), and it is linked with several foodborne outbreaks in the U.S. (5557); this serovar was only isolated in samples from Brazil. These variations were expected and might be due to natural and anthropogenic conditions such as climates, fauna, local policies, and sanitation conditions. For example, in Brazil, vaccination in poultry against S. Enteritidis is widespread (58). Since poultry is the main reservoir for this serovar, this practice could restrict the spread of S. Enteritidis into the environment.
It is estimated that by 2050, AMR will be responsible for the deaths of 10 million people annually (59). Enterobacteriaceae is among the bacterial families most significantly impacted by the rise in AMR rates, and Salmonella is no exception (60). In our study, almost one-third of the Salmonella isolates contained AMR determinants, indicating that AMR genetic elements are widespread in the surveyed areas. Moreover, some isolates carried over 10 AMR determinants and were predicted to be resistant to nine antimicrobial families simultaneously. These findings showed that water environments act as a reservoir for AMR determinants, facilitating their dissemination in the environment and possible transfer to other bacteria. Antimicrobial residues in surface waters might increase this problem, which has been described in at least one of the watersheds sampled in this study (61). These findings encourage us to further investigate the phenotypic antimicrobial resistance traits of Salmonella from this study in the future.
The phylogenetic analysis showed intra-serovar diversity. Isolates from the same serovar but collected in different countries were generally not clonal as shown by the multiple cgMLST analysis, including genomes of single serovars, incrementing the technique’s resolution. This highlights the relevance of performing additional analysis in front of apparently closely related genomes to detect whether the phylogenetic relatedness observed is actual or an artifact of the analysis. Some exceptions include some isolates of worldwide spread serovars Enteritidis, Infantis, and Typhimurium; large PD clusters of these serovars contained not only water isolates from multiple countries in the study but also hundreds of genomes of diverse origins, including clinical, animal, and environmental isolates. Numerous plausible hypotheses can be made about the international transmission of these isolates; however, in the absence of a trace-back investigation, these remain speculative at best. This makes it exceedingly challenging to draw definitive conclusions, as the causes may include international food trade, travel, wildlife migration (e.g., avian migration), biofilms in transoceanic boats, or other factors that remain undefined.
Furthermore, we observed that multiple isolates from the same country and serovar were clonal or closely related, even from samples collected 3 years apart, suggesting the long-term persistence of specific strains in water sources. This persistence could be attributed to various factors, including the ability of Salmonella to form biofilms and survive in diverse hosts, among others (62). This is interesting since surface waters can integrate the environment with both animal production and human activities, serving as a valuable sample in the One Health context.
Interestingly, over 50% of the PD clusters comprised genomes generated in this study did not include genomes from other countries, suggesting that other unique, undiscovered Salmonella may exist in unexplored environments and that diversity is widely more extensive than expected. However, some genomes originated from clinical isolates were phylogenetically close to water genomes of this study with differences as close as 0 SNPs. The closeness between some of the water isolates and isolates from clinical cases, especially in the U.S. and the UK, cannot be dismissed, but the epidemiological link between these genomes has yet to be determined. More importantly, this finding highlights the presence of high-risk Salmonella in surface waters in Latin America.
A significant portion of the study was conducted during the COVID-19 pandemic. The disruptions caused by this worldwide event in this study included changes in sampling schedules and locations, and it could have further impacted the representativeness of our findings. As a result, it is important to note that the Salmonella isolation rates obtained from this research are only applicable to the study design and sampled areas and should not be interpreted as representative of the pathogen’s prevalence in surface waters across Latin America or an individual country.
Foodborne outbreaks associated with produce have been linked to irrigation water and, specifically, to surface waters (9, 63). In Latin America, surface waters are one of the primary water sources for crop irrigation (64); therefore, the presence of high-risk foodborne pathogens in these waters could represent a risk to consumers—the three countries in this study, Brazil, Chile, and Mexico, are prominent food exporters worldwide. Mexico is a leading exporter of produce to its neighbor country, the United States, and Chile exports fruits in counter-season to the northern hemisphere. Brazil mainly exports meats, corn, and soybeans, and it is one of the leading producers of fruit juices worldwide. To date, some vegetable-origin foods from these countries have been linked to foodborne outbreaks in the U.S. For example, Mexican onions contaminated with Salmonella Oranienburg were attributed to causing an outbreak in 2021 (65). Accordingly, our study provides valuable information showing that Salmonella isolates present in surface waters in Latin America are highly genetically related to Salmonella found in different food products and to others that have caused clinical cases of salmonellosis. Therefore, it is crucial to investigate foodborne diseases using the One Health perspective to better understand the dynamics of these types of pathogens in the environment and prevent more cases of foodborne diseases.
Surface waters in agricultural regions of Brazil, Chile, and Mexico were investigated in this study. The genomic information gleaned from our sampling efforts sheds light on the persistence and diversity of Salmonella in these regions. The extensive sampling and sequencing employed in our study distinguish it from previous investigations. For a more comprehensive understanding of the ecological dynamics of Salmonella in agricultural surface waters, we integrated these diverse components. While we observed distinctions and similarities in the distribution of Salmonella serovar prevalence among different countries and time periods, we also examined the resistome and evolutionary connections of the isolates. Even though antibiotic resistance, serovar risks, and epidemics are well-known topics, this study provides current data and context for the agricultural regions surveyed. Furthermore, we highlight the persistence of certain Salmonella strains over time, which has important implications for public health interventions.

Conclusion

This study highlights an important public health issue regarding surface waters contaminated with clinically significant Salmonella serovars in Latin America, some carrying multiple antimicrobial-resistance determinants. The findings underline the increased risks the pathogen poses to human and animal health and reveal the remarkable diversity of Salmonella in the region. We observed long-term contamination in surface waters within countries, which could be explained either by the persistence of strains in the aquatic environments or arising from repeated contamination events from other constantly disseminating sources. These findings suggest that most strains are rather restricted to specific regions in Latin America. The presence of many high-risk Salmonella serovars in surface waters, a key source of irrigation for produce production, poses a significant risk not only to those in direct contact with contaminated waters but also to the broader population through the consumption of produce irrigated with these waters. Improved, coordinated surveillance, control, and prevention strategies are needed to mitigate the spread and impact of Salmonella on public health to ensure food safety and protect public health across the region and beyond.
This study represents a comprehensive analysis of Salmonella diversity in surface waters across agricultural regions of Brazil, Chile, and Mexico, offering significant insights into the genomic diversity and persistence of Salmonella in Latin America that has been relatively underrepresented in previous studies.

ACKNOWLEDGMENTS

This research is supported by the FDA of the U.S. Department of Health and Human Services (HHS) as part of financial assistance award U01FDU001418.
We acknowledge the work of Maria Balkey in data submission, the CFSAN strain curation team and the Salmonella HPOP team for strain management at the FDA. Also, we thank Leonela Díaz, Raul Guevara, and Sebastian Gutierrez from INTA, University of Chile; Francisca Álvarez and Constanza Díaz from Pontificia Universidad Católica de Chile; personnel, social service, undergraduate, and graduate students at the Department of Veterinary Preventive Medicine, Faculty of Veterinary Medicine, Universidad Autónoma de Mexico; Alan Douglas L. Rocha, Laiorayne A. Lima, Gustavo F. C. Sales, Elma L. Leite, Nadyra Jeronimo, and Juliana Alves from LAPOA/CCA/ Federal University of Paraiba; Ana Beatriz S. R. da Silva, Ana Paula S. da Silva, Arthur L.L.de Araújo, Dennys M. Girão, Esther H. R. B. Prado, Francisca E. S. Almeida, Luca O. Valdez, Rossiane M. Souza, and Vinicius C. Moura Medical Microbiology Research Laboratory, Paulo de Góes Institute of Microbiology, Federal University of Rio de Janeiro.

SUPPLEMENTAL MATERIAL

Figure S1 - mbio.00777-24-s0001.tif
MST serovars not closely related.
Figure S2 - mbio.00777-24-s0002.tif
Pathogen detection cluster.
Supplemental material - mbio.00777-24-s0003.xlsx
Data Sets S1 and S2 and Tables S1 to S8.
ASM does not own the copyrights to Supplemental Material that may be linked to, or accessed through, an article. The authors have granted ASM a non-exclusive, world-wide license to publish the Supplemental Material files. Please contact the corresponding author directly for reuse.

REFERENCES

1.
GBD 2017 Causes of Death Collaborators. 2018. Global, regional, and national age-sex-specific mortality for 282 causes of death in 195 countries and territories, 1980–2017: a systematic analysis for the global burden of disease study 2017. Lancet 392:1736–1788.
2.
Cheng RA, Eade CR, Wiedmann M. 2019. Embracing diversity: differences in virulence mechanisms, disease severity, and host adaptations contribute to the success of nontyphoidal Salmonella as a foodborne pathogen. Front Microbiol 10:1368.
3.
White AE, Tillman AR, Hedberg C, Bruce BB, Batz M, Seys SA, Dewey-Mattia D, Bazaco MC, Walter ES. 2022. Foodborne illness outbreaks reported to national surveillance, United States, 2009–2018. Emerg Infect Dis 28:1117–1127.
4.
Chanamé Pinedo L, Mughini-Gras L, Franz E, Hald T, Pires SM. 2022. Sources and trends of human salmonellosis in Europe, 2015–2019: an analysis of outbreak data. Int J Food Microbiol 379:109850.
5.
List of Multistate Foodborne outbreak notices. CDC. Accessed April 22, April 2024. https://www.cdc.gov/foodborne-outbreaks/active-investigations/all-foodborne-outbreak-notices.html.
6.
Walsh KA, Bennett SD, Mahovic M, Gould LH. 2014. Outbreaks associated with cantaloupe, watermelon, and honeydew in the United States, 1973–2011. Foodborne Pathog Dis 11:945–952.
7.
Hanning IB, Nutt JD, Ricke SC. 2009. Salmonellosis outbreaks in the United States due to fresh produce: sources and potential intervention measures. Foodborne Pathog Dis 6:635–648.
8.
Bell RL, Kase JA, Harrison LM, Balan KV, Babu U, Chen Y, Macarisin D, Kwon HJ, Zheng J, Stevens EL, Meng J, Brown EW. 2021. The persistence of bacterial pathogens in surface water and its impact on global food safety. Pathogens 10:1391.
9.
Liu H, Whitehouse CA, Li B. 2018. Presence and persistence of Salmonella in water: the impact on microbial quality of water and food safety. Front Public Health 6:159.
10.
Malakar A, Snow DD, Ray C. 2019. Irrigation water quality-a contemporary perspective. Water 11:1482.
11.
Bell RL, Zheng J, Burrows E, Allard S, Wang CY, Keys CE, Melka DC, Strain E, Luo Y, Allard MW, Rideout S, Brown EW. 2015. Ecological prevalence, genetic diversity, and epidemiological aspects of Salmonella isolated from tomato agricultural regions of the Virginia Eastern Shore. Front Microbiol 6:415.
12.
Song Q, Zhang D, Gao H, Wu J. 2018. Salmonella species’ persistence and their high level of antimicrobial resistance in flooded man-made rivers in China. Microb Drug Resist 24:1404–1411.
13.
Jokinen CC, Koot J, Cole L, Desruisseau A, Edge TA, Khan IUH, Koning W, Lapen DR, Pintar KDM, Reid-Smith R, Thomas JL, Topp E, Wang LY, Wilkes G, Ziebell K, van Bochove E, Gannon VPJ. 2015. The distribution of Salmonella enterica serovars and subtypes in surface water from five agricultural regions across Canada. Water Res 76:120–131.
14.
Ho Y-N, Tsai H-C, Hsu B-M, Chiou C-S. 2018. The association of Salmonella enterica from aquatic environmental and clinical samples in Taiwan. Sci Total Environ 624:106–113.
15.
Mahagamage M, Pathirage M, Manage PM. 2020. Contamination status of Salmonella spp., Shigella spp. and Campylobacter spp. in surface and groundwater of the Kelani River basin, Sri Lanka. Water 12:2187.
16.
Rocha AD de L, Ferrari RG, Pereira WE, de Lima LA, Givisiez PEN, Moreno-Switt AI, Toro M, Delgado-Suárez EJ, Meng J, de Oliveira CJB. 2022. Revisiting the biological behavior of Salmonella enterica in hydric resources: a meta-analysis study addressing the critical role of environmental water on food safety and public health. Front Microbiol 13:802625.
17.
Deaven AM, Ferreira CM, Reed EA, Chen See JR, Lee NA, Almaraz E, Rios PC, Marogi JG, Lamendella R, Zheng J, Bell RL, Shariat NW. 2021. Salmonella genomics and population analyses reveal high inter- and intraserovar diversity in freshwater. Appl Environ Microbiol 87:e02594-20.
18.
Martínez MC, Retamal P, Rojas-Aedo JF, Fernández J, Fernández A, Lapierre L. 2017. Multidrug-resistant outbreak-associated Salmonella strains in irrigation water from the Metropolitan Region, Chile. Zoonoses Public Health 64:299–304.
19.
González-López I, Medrano-Félix JA, Castro-del Campo N, López-Cuevas O, González-Gómez JP, Valdez-Torres JB, Aguirre-Sánchez JR, Martínez-Urtaza J, Gómez-Gil B, Lee BG, Quiñones B, Chaidez C. 2022. Prevalence and genomic diversity of Salmonella enterica recovered from river water in a major agricultural region in northwestern Mexico. Microorganisms 10:1214.
20.
McEgan R, Mootian G, Goodridge LD, Schaffner DW, Danyluk MD. 2013. Predicting Salmonella populations from biological, chemical, and physical indicators in Florida surface waters. Appl Environ Microbiol 79:4094–4105.
21.
Toro M, Weller D, Ramos R, Diaz L, Alvarez FP, Reyes-Jara A, Moreno-Switt AI, Meng J, Adell AD. 2022. Environmental and anthropogenic factors associated with the likelihood of detecting Salmonella in agricultural watersheds. Environ Pollut 306:119298.
22.
Ballesteros-Nova NE, Sánchez S, Steffani JL, Sierra LC, Chen Z, Ruíz-López FA, Bell RL, Reed EA, Balkey M, Rubio-Lozano MS, Soberanis-Ramos O, Barona-Gómez F, Brown EW, Allard MW, Meng J, Delgado-Suárez EJ, Dudley EG. 2022. Genomic epidemiology of Salmonella enterica circulating in surface waters used in agriculture and aquaculture in central Mexico. Appl Environ Microbiol 88:e0214921.
23.
Sbodio A, Maeda S, Lopez-Velasco G, Suslow TV. 2013. Modified Moore swab optimization and validation in capturing E. coli O157: H7 and Salmonella enterica in large volume field samples of irrigation water. Food Res Int 51:654–662.
24.
Kim JS, Lee GG, Park JS, Jung YH, Kwak HS, Kim SB, Nam YS, Kwon ST. 2007. A novel multiplex PCR assay for rapid and simultaneous detection of five pathogenic bacteria: Escherichia coli O157:H7, Salmonella, Staphylococcus aureus, Listeria monocytogenes, and Vibrio parahaemolyticus. J Food Prot 70:1656–1662.
25.
Argimón S, Abudahab K, Goater RJE, Fedosejev A, Bhai J, Glasner C, Feil EJ, Holden MTG, Yeats CA, Grundmann H, Spratt BG, Aanensen DM. 2016. Microreact: Visualizing and sharing data for genomic epidemiology and phylogeography. Microb Genom 2:e000093.
26.
Gangiredla J, Rand H, Benisatto D, Payne J, Strittmatter C, Sanders J, Wolfgang WJ, Libuit K, Herrick JB, Prarat M, Toro M, Farrell T, Strain E. 2021. GalaxyTrakr: a distributed analysis tool for public health whole genome sequence data accessible to non-bioinformaticians. BMC Genomics 22:114.
27.
Bolger AM, Lohse M, Usadel B. 2014. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30:2114–2120.
28.
Souvorov A, Agarwala R, Lipman DJ. 2018. SKESA: strategic k-mer extension for scrupulous assemblies. Genome Biol 19:153.
29.
Seemann T. 2020. MLST. Available from: https://github.com/tseemann/mlst
30.
Maiden MCJ, Jansen van Rensburg MJ, Bray JE, Earle SG, Ford SA, Jolley KA, McCarthy ND. 2013. MLST revisited: the gene-by-gene approach to bacterial genomics. Nat Rev Microbiol 11:728–736.
31.
Petit RA. 2022. fastq-scan. Available from: https://github.com/rpetit3/fastq-scan
32.
Timme RE, Wolfgang WJ, Balkey M, Venkata SLG, Randolph R, Allard M, Strain E. 2020. Optimizing open data to support one health: best practices to ensure Interoperability of genomic data from bacterial pathogens. One Health Outlook 2:20.
33.
Alikhan NF, Zhou Z, Sergeant MJ, Achtman M. 2018. A genomic overview of the population structure of Salmonella. PLoS Genet 14:e1007261.
34.
Yoshida CE, Kruczkiewicz P, Laing CR, Lingohr EJ, Gannon VPJ, Nash JHE, Taboada EN. 2016. The Salmonella in silico typing resource (SISTR): an open web-accessible tool for rapidly typing and subtyping draft Salmonella genome assemblies. PLoS One 11:e0147101.
35.
Feldgarden M, Brover V, Gonzalez-Escalona N, Frye JG, Haendiges J, Haft DH, Hoffmann M, Pettengill JB, Prasad AB, Tillman GE, Tyson GH, Klimke W. 2021. AMRFinderPlus and the reference gene catalog facilitate examination of the genomic links among antimicrobial resistance, stress response, and virulence. Sci Rep 11:12728.
36.
Hammer DAT, Ryan PD, Hammer Ø, Harper DAT. 2001. PAST: paleontological statistics software package for education and data analysis. Palaeont Electr 4:4–9.
37.
Pightling AW, Pettengill JB, Luo Y, Baugher JD, Rand H, Strain E. 2018. Interpreting whole-genome sequence analyses of foodborne bacteria for regulatory applications and outbreak investigations. Front Microbiol 9:1482.
38.
Cho S, Jackson CR, Frye JG. 2020. The prevalence and antimicrobial resistance phenotypes of Salmonella, Escherichia coli and Enterococcus sp. in surface water. Lett Appl Microbiol 71:3–25.
39.
Thompson CP, Doak AN, Amirani N, Schroeder EA, Wright J, Kariyawasam S, Lamendella R, Shariat NW. 2018. High-resolution identification of multiple Salmonella serovars in a single sample by using CRISPRSeroSeq. Appl Environ Microbiol 84:e01859-18.
40.
Rasamsetti S, Shariat NW. 2023. Biomapping Salmonella serovar complexity in broiler carcasses and parts during processing. Food Microbiol 110:104149.
41.
CDC. Multistate outbreak of Salmonella Newport - June 29, 2010 - Salmonella. Available from: https://www.cdc.gov/salmonella/2010/newport-alfalfa-sprout-6-29-10.html. Retrieved 14 Nov 2023.
42.
Outbreak of Salmonella Newport infections linked to cucumbers — United States,2014. Available from: https://www.cdc.gov/mmwr/preview/mmwrhtml/mm6406a3.htm?s_cid=mm6406a3_e. Retrieved 14 Nov 2023.
43.
CDC. Outbreak of Salmonella Newport infections linked to onions. Available from: https://www.cdc.gov/salmonella/newport-07-20/index.html. Retrieved 14 Nov 2023.
44.
Bennett SD, Littrell KW, Hill TA, Mahovic M, Behravesh CB. 2015. Multistate foodborne disease outbreaks associated with raw tomatoes, United States, 1990-2010: a recurring public health problem. Epidemiol Infect 143:1352–1359.
45.
Santos AF das M, Amparo LFV, Machado SCA, Dias TS, Berto LH, Abreu DL da C, Aquino MHC de, Rodrigues D dos P, Pereira VL de A. 2022. Salmonella serovars associated with human salmonellosis in Brazil (2011-2020). RSD 11:e28011830533.
46.
Godínez-Oviedo A, Tamplin ML, Bowman JP, Hernández-Iturriaga M. 2020. Salmonella enterica in Mexico 2000-2017: epidemiology, antimicrobial resistance, and prevalence in food. Foodborne Pathog Dis 17:98–118.
47.
ISP. 2019. Boletin de vigilancia de Laboratorio. Salmonella spp. 2014-2018 9:13. Retrieved Nov 26 Nov 2023. https://www.ispch.cl/sites/default/files/BoletínSalmonella-12052020A.pdf.
48.
Fernandes SA, Tavechio AT, Ghilardi ÂCR, Almeida EA de, Silva JML da, Camargo CH, Tiba-Casas MR. 2022. Salmonella enterica serotypes from human and nonhuman sources in Sao Paulo state. Rev Inst Med Trop Sao Paulo 64:e66.
49.
Toro M, Rivera D, Toledo V, Campos-Vargas R, Allard MW, Hamilton-West C, Moreno-Switt AI. 2018. Genomics of Salmonella contaminating backyard production systems reveals persistence and transmission of genetically related Salmonella on a farm basis. Zoonoses Public Health 65:1008–1014.
50.
CDC. Salmonella Typhimurium and Salmonella Newport infections linked to cantaloupe - Salmonella. Available from: https://www.cdc.gov/salmonella/typhimurium-cantaloupe-08-12/index.html. Retrieved 26 Nov 2023.
51.
Salmonella outbreak linked to alfalfa sprouts. CDC. Available from: https://www.cdc.gov/salmonella/typhimurium-12-22/index.html. Retrieved 26 Nov 2023.
52.
CDC. Salmonella Typhimurium infections linked to tomatoes - Salmonella. Available from: https://www.cdc.gov/salmonella/2006/tomatoes-11-2006.html. Retrieved 26 Nov 2023.
53.
Alvarez DM, Barrón-Montenegro R, Conejeros J, Rivera D, Undurraga EA, Moreno-Switt AI. 2023. A review of the global emergence of multidrug-resistant Salmonella enterica subsp. enterica Serovar infantis. Int J Food Microbiol 403:110297.
54.
Salmonella Newport and Salmonella Infantis infections linked to imported Maradol papayas | September 2017 | Salmonella. CDC. Available from: https://www.cdc.gov/salmonella/newport-09-17/index.html. Retrieved 26 Nov 2023.
55.
CDC. Salmonella Heidelberg infections linked to ground turkey. Available from: https://archive.cdc.gov/www_cdc_gov/salmonella/2011/ground-turkey-11-10-2011.html. Retrieved 14 Nov 2023.
56.
CDC. Multistate outbreak of multidrug-resistant Salmonella Heidelberg infections linked to foster farms brand chicken. Available from: https://www.cdc.gov/salmonella/heidelberg-10-13/index.html. Retrieved 14 Nov 2023.
57.
CDC. Salmonella infections linked to kosher broiled chicken livers. Available from: https://archive.cdc.gov/www_cdc_gov/salmonella/2011/chicken-liver-1-11-2012.html. Retrieved 14 Nov 2023.
58.
Alikhan N-F, Moreno LZ, Castellanos LR, Chattaway MA, McLauchlin J, Lodge M, O’Grady J, Zamudio R, Doughty E, Petrovska L, Cunha MPV, Knöbl T, Moreno AM, Mather AE. 2022. Dynamics of Salmonella enterica and antimicrobial resistance in the Brazilian poultry industry and global impacts on public health. PLoS Genet 18:e1010174.
59.
United Nations Environment Programme. Bracing for superbugs: strengthening environmental action in the one health response to antimicrobial resistance. Available from: https://www.unep.org/resources/superbugs/environmental-action. Retrieved 14 Nov 2023.
60.
Murray CJL, Ikuta KS, Sharara F, Swetschinski L, Robles Aguilar G, Gray A, Han C, Bisignano C, Rao P, Wool E, et al. 2022. Global burden of bacterial antimicrobial resistance in 2019: a systematic analysis. The Lancet 399:629–655.
61.
Wilkinson JL, Boxall ABA, Kolpin DW, Leung KMY, Lai RWS, Galbán-Malagón C, Adell AD, Mondon J, Metian M, Marchant RA, et al. 2022. Pharmaceutical pollution of the world’s rivers. Proc Natl Acad Sci USA 119:e2113947119.
62.
Douesnard-Malo F, Daigle F, Frédéric D-M, France D. 2011. Increased persistence of Salmonella enterica serovar Typhi in the presence of Acanthamoeba castellanii. Appl Environ Microbiol 77:7640–7646.
63.
Gelting RJ, Baloch MA, Zarate-Bermudez MA, Selman C. 2011. Irrigation water issues potentially related to the 2006 multistate E. coli O157:H7 outbreak associated with spinach. Agric Water Manag 98:1395–1402.
64.
Thenkabail PS, Biradar CM, Noojipady P, Dheeravath V, Li Y, Velpuri M, Gumma M, Gangalakunta ORP, Turral H, Cai X, Vithanage J, Schull MA, Dutta R. 2009. Global irrigated area map (GIAM), derived from remote sensing, for the end of the last millennium. Int J Remote Sens 30:3679–3733.
65.
Outbreak investigation of Salmonella Oranienburg: whole, fresh onions (October 2021). FDA. Available from: https://www.fda.gov/food/outbreaks-foodborne-illness/outbreak-investigation-salmonella-oranienburg-whole-fresh-onions-october-2021. Retrieved 14 Nov 2023.

Information & Contributors

Information

Published In

cover image mBio
mBio
Volume 15Number 717 July 2024
eLocator: e00777-24
Editor: Francisco Diez-Gonzalez, University of Georgia Center for Food Safety, Griffin, Georgia, USA
PubMed: 38920393

History

Received: 13 March 2024
Accepted: 21 May 2024
Published online: 26 June 2024

Keywords

  1. water microbiological contamination
  2. WGS
  3. Salmonella persistence
  4. Latin America

Data Availability

All reads were submitted as FASTQ files to the Sequence Read Archive (SRA), NCBI, under BioProject numbers PRJNA186035 and PRJNA560080 as part of FDA’s GenomeTrakr Project (Table S1).

Contributors

Authors

Joint Institute for Food Safety and Applied Nutrition, University of Maryland, College Park, Maryland, USA
Center for Food Safety and Security Systems, University of Maryland, College Park, Maryland, USA
Author Contributions: Conceptualization, Data curation, Formal analysis, Investigation, and Writing – review and editing.
Andrea I. Moreno-Switt
Escuela de Medicina Veterinaria, Pontificia Universidad Católica de Chile, Santiago, Chile
Author Contributions: Funding acquisition, Investigation, Project administration, Resources, Supervision, and Writing – review and editing.
Angelica Reyes-Jara
Institute of Nutrition and Food Technology, University of Chile, Santiago, Chile
Author Contributions: Funding acquisition, Investigation, Project administration, Resources, Supervision, and Writing – review and editing.
Enrique Delgado Suarez
Faculty of Veterinary Medicine and Zootechnics, National Autonomous University of Mexico, Mexico City, Mexico
Author Contributions: Funding acquisition, Investigation, Project administration, Resources, and Writing – review and editing.
Aiko D. Adell
Escuela de Medicina Veterinaria, Facultad de Ciencias de la Vida, Universidad Andrés Bello, Santiago, Chile
Author Contributions: Investigation, Methodology, Supervision, and Writing – review and editing.
Celso José Bruno Oliveira https://orcid.org/0000-0002-7761-0697
Department of Animal Science, Federal University of Paraiba, Areia, Brazil
Author Contributions: Funding acquisition, Investigation, Project administration, Resources, Supervision, and Writing – review and editing.
Medical Microbiology Research Laboratory, Paulo de Góes Institute of Microbiology, Federal University of Rio de Janeiro, Rio de Janeiro, Brazil
Author Contributions: Funding acquisition, Investigation, Project administration, Resources, Supervision, and Writing – review and editing.
Xinyang Huang
Joint Institute for Food Safety and Applied Nutrition, University of Maryland, College Park, Maryland, USA
Center for Food Safety and Security Systems, University of Maryland, College Park, Maryland, USA
Author Contributions: Investigation and Writing – review and editing.
Eric Brown
Center for Food Safety and Applied Nutrition, Food and Drug Administration, College Park, Maryland, USA
Author Contributions: Conceptualization, Funding acquisition, Resources, and Writing – review and editing.
Marc Allard
Center for Food Safety and Applied Nutrition, Food and Drug Administration, College Park, Maryland, USA
Author Contributions: Conceptualization, Resources, and Writing – review and editing.
Christopher Grim
Center for Food Safety and Applied Nutrition, Food and Drug Administration, College Park, Maryland, USA
Author Contributions: Data curation, Investigation, Resources, and Writing – review and editing.
Rebecca Bell
Center for Food Safety and Applied Nutrition, Food and Drug Administration, College Park, Maryland, USA
Author Contributions: Conceptualization, Investigation, Methodology, Resources, and Writing – review and editing.
Jianghong Meng
Joint Institute for Food Safety and Applied Nutrition, University of Maryland, College Park, Maryland, USA
Center for Food Safety and Security Systems, University of Maryland, College Park, Maryland, USA
Department of Nutrition and Food Science, University of Maryland, College Park, Maryland, USA
Author Contributions: Conceptualization, Funding acquisition, Investigation, Project administration, Resources, Supervision, and Writing – review and editing.
Joint Institute for Food Safety and Applied Nutrition, University of Maryland, College Park, Maryland, USA
Center for Food Safety and Security Systems, University of Maryland, College Park, Maryland, USA
Institute of Nutrition and Food Technology, University of Chile, Santiago, Chile
Author Contributions: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – original draft, and Writing – review and editing.

Editor

Francisco Diez-Gonzalez
Editor
University of Georgia Center for Food Safety, Griffin, Georgia, USA

Notes

The authors declare no conflict of interest.

Metrics & Citations

Metrics

Note:

  • For recently published articles, the TOTAL download count will appear as zero until a new month starts.
  • There is a 3- to 4-day delay in article usage, so article usage will not appear immediately after publication.
  • Citation counts come from the Crossref Cited by service.

Citations

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. For an editable text file, please select Medlars format which will download as a .txt file. Simply select your manager software from the list below and click Download.

View Options

Figures

Tables

Media

Share

Share

Share the article link

Share with email

Email a colleague

Share on social media

American Society for Microbiology ("ASM") is committed to maintaining your confidence and trust with respect to the information we collect from you on websites owned and operated by ASM ("ASM Web Sites") and other sources. This Privacy Policy sets forth the information we collect about you, how we use this information and the choices you have about how we use such information.
FIND OUT MORE about the privacy policy