Changing of the Genomic Pattern of Salmonella Enteritidis Strains Isolated in Brazil Over a 48 year-period revealed by Whole Genome SNP Analyses

Salmonella in eggs and egg-laying chickens: pathways to effective control

Eggs contaminated with Salmonella have been internationally significant sources of human illness for several decades. Most egg-associated illness has been attributed to Salmonella serovar Enteritidis, but a few other serovars (notably S. Heidelberg and S. Typhimurium) are also sometimes implicated. The edible interior contents of eggs typically become contaminated with S. Enteritidis because the pathogen's unique virulence attributes enable it to colonize reproductive tissues in systemically infected laying hens. Other serovars are more commonly associated with surface contamination of eggshells. Both research and field experience have demonstrated that the most effective overall Salmonella control strategy in commercial laying flocks is the application of multiple interventions throughout the egg production cycle. At the preharvest (egg production) level, intervention options of demonstrated efficacy include vaccination and gastrointestinal colonization control via treatments such as prebiotics, probiotics, and bacteriophages, Effective environmental management of housing systems used for commercial laying flocks is also essential for minimizing opportunities for the introduction, transmission, and persistence of Salmonella in laying flocks. At the postharvest (egg processing and handling) level, careful regulation of egg storage temperatures is critical for limiting Salmonella multiplication inside the interior contents.

A machine learning approach to identifying Salmonella stress response genes in isolates from poultry processing

We explored the potential of machine learning to identify significant genes associated with Salmonella stress response during poultry processing using whole genome sequencing (WGS) data. The Salmonella isolates (n = 177) used in this study were obtained from various chicken sources (skin before chiller, chicken carcass before chiller, frozen chicken, and post-chill chicken carcass). Six machine learning algorithms (random forest, neural network, cost-sensitive learning, logit boost, and support vector machine linear and radial kernels) were trained on Salmonella WGS data, and model fit was assessed using standard evaluation metrics such as the area under the receiver operating characteristic (AUROC) curve and confusion matrix statistics. All models achieved high performances based on the AUROC metric, with logit boost showing the best performance with an AUROC score of 0.904, sensitivity of 0.889, and specificity of 0.920. The significant genes identified included ybtX, which encodes a Yersiniabactin-associated zinc transporter, and the transferase-encoding genes yccK and thiS. Additionally, genes coding for cold (cspA, cspD, and cspE) and heat shock (rpoH and rpoE) responses were identified. Other significant genes included those involved in lipopolysaccharide biosynthesis (irp1, waaD, rfc, and rfbX), DNA repair and replication (traI), biofilm formation (ccdA and fyuA), and cellular metabolism (irtA).

Genome-Wide Searching Single Nucleotide-Polymorphisms (SNPs) and SNPs-Targeting a Multiplex Primer for Identification of Common Salmonella Serotypes

A rapid and high-quality single-nucleotide polymorphisms (SNPs)-based method was developed to improve detection and reduce salmonellosis burden. In this study, whole-genome sequence (WGS) was used to investigate SNPs, the most common genetic marker for identifying bacteria. SNP-sites encompassing 15 sets of primers (666–863 bp) were selected and used to amplify the target Salmonella serovar strains, and the amplified products were sequenced. The prevalent Salmonella enterica subspecies enterica serovars, including Typhimurium; Enteritidis, Agona, enterica, Typhi, and Abony, were amplified and sequenced. The amplified sequences of six Salmonella serovars with 15 sets of SNP-sites encompassing primers were aligned, explored SNPs, and SNPs-carrying primers (23 sets) were designed to develop a multiplex PCR marker (m-PCR). Each primer exists in at least two SNPs bases at the 3′ end of each primer, such as one was wild, and another was a mismatched base by transition or transversion mutation. Thus, twenty-three sets of SNP primers (242–670 bp), including 13 genes (SBG, dedA, yacG, mrcB, mesJ, metN, rihA/B, modA, hutG, yehX, ybiY, moeB, and sopA), were developed for PCR confirmation of target Salmonella serovar strains. Finally, the SNPs in four genes, including fliA gene (S. Enteritidis), modA (S. Agona and S. enterica), sopA (S. Abony), and mrcB (S. Typhimurium and S. Typhi), were used for detection markers of six target Salmonella serotypes. We developed an m-PCR primer set in which Salmonella serovars were detected in a single reaction. Nevertheless, m-PCR was validated with 21 Salmonella isolates (at least one isolate was taken from one positive animal fecal, and n = 6 reference Salmonella strains) and non-Salmonella bacteria isolates. The SNP-based m-PCR method would identify prevalent Salmonella serotypes, minimize the infection, and control outbreaks.

Characterisation of Salmonella Enteritidis ST11 and ST1925 Associated with Human Intestinal and Extra-Intestinal Infections in Singapore

Salmonella Enteritidis is a major foodborne pathogen worldwide. In this study, a total of 276 S. enteritidis isolates, collected between 2016 and 2017 from human, food and farm/slaughterhouse samples, were studied to enhance the understanding of the epidemiology of human salmonellosis in Singapore. Results showed all 276 isolates belonged either to ST1925 (70.3%) or ST11 (29.7%), with ST11 being significantly more frequent in extra-intestinal isolates and chicken isolates. Food isolates, most of which were from poultry, showed the highest prevalence of resistance (33–37%) against beta-lactams or beta-lactams/beta-lactamase inhibitor combination (ampicillin, piperacillin and ampicillin/sulbactam). The analysis showed the detection of genes associated with resistance to aminoglycoside genes (99.6%), tetracycline (55.1%), and beta-lactams (14.9%) of all isolates. Nine types of plasmids were found in 266 isolates; the most common incompatibility group profiles were IncFIB(S)-IncFII(S)-IncX1 (72.2%) and IncFIB(S)-IncFII(S) (15.8%). Most plasmid harbouring isolates from chicken (63.6%, 14/22) and from human (73.8%, 175/237) shared the same plasmid profile (IncFIB(S)-IncFII(S)-IncX1). SNP analysis showed clustering of several isolates from poultry food products and human isolates, suggesting phylogenetic relatedness among these isolates. Lastly, this study provides important epidemiological insights on the application of phenotypic and next-generation sequencing (NGS) tools for improved food safety and public health surveillance and outbreak investigation of S.enteritidis.

Characterization of Salmonella enterica Contamination in Pork and Poultry Meat from São Paulo/Brazil: Serotypes, Genotypes and Antimicrobial Resistance Profiles

Salmonellosis is a zoonosis of major relevance to global public health. Here we present the assessment of Salmonella enterica contamination in pork and poultry meat sold at retail markets in São Paulo, Brazil. A total of 780 meat samples (386 poultry meat and 394 pork samples) were collected from 132 markets. From these, 57 samples (7.3%) were positive for S. enterica isolation, including 32 (8.3%) poultry meat and 25 (6.3%) pork samples. S. enterica isolates were further characterized for serotyping, antimicrobial resistance and genotyping by amplified fragment length polymorphism and pulsed field gel electrophoresis. Antimicrobial resistance analysis demonstrated two main profiles: pork isolates were more resistant to macrolides, β-lactams, tetracycline, phenicols, and fluoroquinolones, and poultry meat isolates presented higher resistance to fluoroquinolones, sulfonamides, tetracycline, and β-lactams. A total of 72.4% of poultry meat isolates were identified as S. Heidelberg, while most of pork isolates were S. Typhimurium (31.7%) and S. Give (16.7%). Genotyping resulted in most clusters consisting exclusively of pork or poultry meat, no cross-contamination was detected, and a tendency to differentiate isolates according to their serotypes and markets of origin. High resistance rates to critically important antimicrobials reinforce the importance of controlling Salmonella contamination in meat production chains.

Phenotypic and Genotypic Characterization of Salmonella Isolates Recovered from Foods Linked to Human Salmonellosis Outbreaks in Minas Gerais State, Brazil

ABSTRACT

Salmonella is one of the primary pathogens that causes foodborne diseases worldwide. In the present study, we characterized Salmonella isolates recovered from foods linked to human salmonellosis outbreaks in Minas Gerais, Brazil, from 2003 to 2017. Serotype, antimicrobial susceptibility, presence of virulence genes, and genetic polymorphism as determined by repetitive element sequence-based PCR were determined for 70 Salmonella isolates. Thirteen Salmonella serotypes were identified, and the most prevalent were Enteritidis and Typhimurium, comprising 52 (74.3%) of the 70 isolates. Sixty-five (92.8%) of the isolates were resistant to at least 1 of the 15 antimicrobials tested. Ten isolates (14.2%) had a multidrug resistance phenotype. Isolates were screened for 16 virulence genes, which were found in 75.7 to 100% of the isolates. A statistical difference was found among Salmonella serotypes in the presence of the sipB, sopE, lfpA, sefA, and spvC genes. Based on their DNA fingerprints, 40 isolates of Salmonella Enteritidis from 16 outbreaks were separated into 14 groups and 12 isolates of Salmonella Typhimurium were separated into 6 groups. These serological patterns were similar to those reported by public health centers worldwide. Of concern is the high prevalence among the isolates in this study of both virulence genes and resistance to antimicrobials, especially to critically important drugs. Special attention should be given to Salmonella Enteritidis. Although the genomes of these Salmonella isolates were relatively variable, high genetic similarity was observed among them, and some had identical fingerprints. These results support the hypothesis of clonal circulation of Salmonella isolates causing human infections in Minas Gerais.

HIGHLIGHTS

Whole genome sequencing and protein structure analyses of target genes for the detection of Salmonella

Rapid and sensitive detection of Salmonella is a critical step in routine food quality control, outbreak investigation, and food recalls. Although various genes have been the targets in the design of rapid molecular detection methods for Salmonella, there is limited information on the diversity of these target genes at the level of DNA sequence and the encoded protein structures. In this study, we investigated the diversity of ten target genes (invA, fimA, phoP, spvC, and agfA; ttrRSBCA operon including 5 genes) commonly used in the detection and identification of Salmonella. To this end, we performed whole genome sequencing of 143 isolates of Salmonella serotypes (Enteritidis, Typhimurium, and Heidelberg) obtained from poultry (eggs and chicken). Phylogenetic analysis showed that Salmonella ser. Typhimurium was more diverse than either Enteritidis or Heidelberg. Forty-five non-synonymous mutations were identified in the target genes from the 143 isolates, with the two most common mutations as T ↔ C (15 times) and A ↔ G (13 times). The gene spvC was primarily present in Salmonella ser. Enteritidis isolates and absent from Heidelberg isolates, whereas ttrR was more conserved (0 non-synonymous mutations) than ttrS, ttrB, ttrC, and ttrA (7, 2, 2, and 7 non-synonymous mutations, respectively). Notably, we found one non-synonymous mutation (fimA-Mut.6) across all Salmonella ser. Enteritidis and Salmonella ser. Heidelberg, C → T (496 nt postion), resulting in the change at AA 166 position, Glutamine (Q) → Stop condon (TAG), suggesting that the fimA gene has questionable sites as a target for detection. Using Phyre² and SWISS-MODEL software, we predicted the structures of the proteins encoded by some of the target genes, illustrating the positions of these non-synonymous mutations that mainly located on the α-helix and β-sheet which are key elements for maintaining the conformation of proteins. These results will facilitate the development of sensitive molecular detection methods for Salmonella.

Discerning the Antimicrobial Resistance, Virulence, and Phylogenetic Relatedness of Salmonella Isolates Across the Human, Poultry, and Food Materials Sources in Malaysia

Salmonella enterica subspecies enterica serovar Enteritidis is one of the major foodborne zoonotic pathogens globally. It has significantly impacted human health and global trade. In this investigation, whole-genome sequencing was employed to determine the antimicrobial resistance (AMR) pattern of a collection of Salmonella Enteritidis isolated from humans, poultry, and food sources. The study also investigated the virulence genes profile of the isolates as well as the phylogenetic relationships among strains. Illumina NextSeq technology was used to sequence the genome of 82 Salmonella Enteritidis strains isolated over 3 years (2016-2018) in Peninsular Malaysia. The pattern of resistance showed that tetracycline had the highest frequency (37/82, 45.12%), and isolates from food samples showed the highest rate of 9/18 (50.00%), followed by human 17/35 (48.57%) and then poultry 11/29 (37.93%). The second drug with the highest resistance rate is ampicillin with 5/29 (17.24%) for poultry, 4/35 (11.43%) for human, and 0/18 (0.00%) for food isolates respectively. Similarly, a total of 19 antimicrobial resistance (AMR) genes corresponding to the nine drugs used in the disc diffusion assay were evaluated from the whole genome sequence data. The aminoglycoside resistance gene aac(6')-ly was detected in 79 of the 82 isolates (96.34%). While the phylogenetic analysis revealed distinct lineages isolated, the three sources indicating possible cross-contamination. In conclusion, the results showed that the genomic profile of Salmonella Enteritidis isolated from humans, poultry, and food samples share genetic traits, hence the need to institute measures at controlling the continuous spread of these resistant pathogens.

Retrospective whole-genome comparison of Salmonella enterica serovar Enteritidis from foodborne outbreaks in Southern Brazil

Salmonella enterica serovar Enteritidis is frequently isolated from animal-source foods associated with human salmonellosis outbreaks. This serovar was spread to animal (mainly poultry) farms worldwide in the 1980s, and it is still detected in foods produced in many countries, including Brazil. The present study reports a retrospective genome-wide comparison of S. Enteritidis from foodborne outbreaks in Southern Brazil in the last two decades. Fifty-two S. Enteritidis isolates were obtained from foodborne outbreaks occurring in different cities of the Brazilian southernmost State, Rio Grande do Sul (RS), from 2003 to 2015. Whole-genome sequences (WGS) from these isolates were obtained and comparatively analyzed with 65 additional genomes from NCBI. Phylogenetic and Bayesian analyses were performed to study temporal evolution. Genes related to antibiotic resistance and virulence were also evaluated. The results demonstrated that all S. Enteritidis isolates from Southern Brazil clustered in the global epidemic clade disseminated worldwide originally in the 1980s. Temporal analysis demonstrated that all Brazilian isolates had a tMRCA (time to most recent common ancestor) in 1986 with an effective population size (N_e) increase soon after until 1992, then becoming constant up to now. In Southern Brazil, there was a significant decrease in the spreading of S. Enteritidis in the last decade. In addition, three antibiotic resistance genes were detected in all isolates: aac(6')-Iaa, mdfA, and tet(34). These results demonstrate the high frequency of one only specific S. Enteritidis lineage (global epidemic clade) in foodborne outbreaks from Southern Brazil in the last two decades.

Genomic investigation of antimicrobial resistance determinants and virulence factors in Salmonella enterica serovars isolated from contaminated food and human stool samples in Brazil

This study investigated the antimicrobial resistance determinants, virulence factors and identified serovars in 37 Salmonella enterica strains isolated from human stool and contaminated foods linked to outbreaks that occurred in Brazil over 7 years using whole genome sequencing (WGS). Phylogenetic analysis of selected serovars (S. Typhimurium, S. Infantis, S. London, and S. Johannesburg) was performed. Ten distinct serovars were identified and, 51% of the tested strains (n = 19) showed disagreement with the previous conventional serotyping. The antimicrobial resistance (AMR) determinants or plasmids varied among the strains. Resistome analysis revealed the presence of resistance genes to aminoglycosides [aac (6')-laa, aph (3″)-lb, aph (6)-ld, aadA1 and aadA2], sulfonamides (sul1), trimethoprin (dfrA8), fosfomycin (fosA7) and tetracyclines (tetA, tetB, tetC), as well as point mutations in parC (T57S) and gyrA (S83F). Plasmidome showed the presence of IncHI2, IncHI2A, IncFIB (S), IncFII (S), IncI1 and p0111 plasmids. Eight Salmonella pathogenicity islands and up to 102 stress and/or virulence genes were identified in the evaluated genomes. Virulence genes of K88 fimbrial adhesin were first reported in S. enterica (S. Pomona, S. Bredeney and S. Mbandaka strains). pilW gene was first identified in S. Pomona. Phylogenetic analysis showed that some serovars circulated in Brazil for decades, primarily within the poultry production chain. Findings highlighted the virulence and AMR determinants in strains that may lead to recurring food outbreaks.

Phylogenetic structure of Salmonella Enteritidis provides context for a foodborne outbreak in Peru

Salmonella Enteritidis, an important foodborne zoonosis, has a dramatically increased number of cases around the world. To explore the phylogenetic structure of Peruvian Salmonella Enteritidis strains and their relationship with an outbreak occurred in 2018, we analyzed a comprehensive strains of S. Enteritidis received by the National Institute of Health during the period 2000-2018. A total of 180 strains were characterized by microbiological procedures, serotyping and whole genome sequencing. Based on genome sequences annotated, virulence factors and accessory genes were identified. Phylogenetic and population structure analysis were also analyzed based on SNPs. The phylogenetic analysis grouped the genomes into two well-supported clades that were consistent with population structure analysis. The clinical and food strains corresponding to the outbreak were included in the same cluster, which presented the sdhA gene, related to the increase of the virulence of this pathogen. The phylogenetic relationship of Peruvian S. Enteritidis suggests the presence of four S. enteritidis population with high epidemiological importance.

Phenotypic analyses of Salmonella enterica serovar Enteritidis strains isolated in the pre- and post-epidemic period in Brazil

Salmonella Enteritidis has caused, since the 1980s, a sustained epidemic of human infections in many countries. This study analyzed S. Enteritidis strains isolated before and after the epidemic period in Brazil regarding their capacities to survive to acid, oxidative, and high-temperature stresses, and capacity to grow in egg albumen. Moreover, the ability to invade human epithelial cells (Caco-2) and to survive inside human (U937) and chicken (HD11) macrophages was checked. Post-epidemic strains showed a better ability to survive after 10 min under acid stress at 37 °C (P ≤ 0.05). However, both groups of strains showed similar ability to survive after 1 h under acid stress at 37 °C and at 42 °C independently of the time of exposure. Similar ability was verified in both groups of strains regarding oxidative stress, growth in egg albumen, high-temperature stress, invasion to Caco-2 cells, and invasion and survival in macrophages. In conclusion, post-epidemic S. Enteritidis strains showed a better ability to survive under the acid stress found in the stomach, which might be an advantage to reach the intestine and colonize chickens and humans. However, both groups of strains did not differ significantly in the majority of the phenotypic tests analyzed in this study.

Draft Genome Sequences of 30 Salmonella enterica Serovar Enteritidis Isolates Associated with Multiple Outbreaks in Brazil

Salmonella enterica subsp. enterica serovar Enteritidis has been the prevalent serovar isolated from gastroenteritis cases in Brazil since the 1990s. Here, we report the draft genomes of 30 S. Enteritidis isolates originating from a variety of patients and implicated foods during outbreaks between 1999 and 2006 in Brazil.

Salmonella enterica subsp. enterica serovar Enteritidis has been the prevalent serovar isolated from gastroenteritis cases in Brazil since the 1990s. Here, we report the draft genomes of 30 S. Enteritidis isolates originating from a variety of patients and implicated foods during outbreaks between 1999 and 2006 in Brazil.

Comparative genomics of Salmonella enterica serovar Enteritidis ST-11 isolated in Uruguay reveals lineages associated with particular epidemiological traits

Salmonella enterica serovar Enteritidis is a major cause of foodborne disease in Uruguay since 1995. We used a genomic approach to study a set of isolates from different sources and years. Whole genome phylogeny showed that most of the strains are distributed in two major lineages (E1 and E2), both belonging to MLST sequence type 11 the major ST among serovar Enteritidis. Strikingly, E2 isolates are over-represented in periods of outbreak abundance in Uruguay, while E1 span all epidemic periods. Both lineages circulate in neighbor countries at the same timescale as in Uruguay, and are present in minor numbers in distant countries. We identified allelic variants associated with each lineage. Three genes, ycdX, pduD and hsdM, have distinctive variants in E1 that may result in defective products. Another four genes (ybiO, yiaN, aas, aceA) present variants specific for the E2 lineage. Overall this work shows that S. enterica serovar Enteritidis strains circulating in Uruguay have the same phylogenetic profile than strains circulating in the region, as well as in more distant countries. Based on these results we hypothesize that the E2 lineage, which is more prevalent during epidemics, exhibits a combination of allelic variants that could be associated with its epidemic ability.

Single-stranded DNA (ssDNA) Aptamer targeting SipA protein inhibits Salmonella Enteritidis invasion of intestinal epithelial cells

Salmonella Enteritidis is an important pathogen that can invade the intestinal cells of its host causing salmonellosis. SipA protein, an effector protein secreted by T3SS, maintains invasion of host cells more efficient. Thus, inhibitory aptamers against SipA protein were developed using magnetic bead-based Systematic Evolution of Ligands by Exponential Enrichment (SELEX) method. The enriched sequences were obtained after 9 SELEX rounds. Among which, an aptamer namely Apt17 displayed K_d values equivalent to 114.9 and 63.4 nM at 27 °C and 37 °C, respectively. The effect of Apt17 on adhesion and invasion of Caco-2 cells by the tested strains was determined. While the adhesion and invasion of Salmonella Enteritidis TM 6 were inhibited by 70% and 37.7%, those of Salmonella Enteritidis TM 68 were inhibited by 45.71% and 39.5% respectively. These results represent a corner stone for future studies that could aim to develop putative inhibitors against Salmonellosis.

Salmonella control in poultry flocks and its public health impact

An increase in confirmed human salmonellosis cases in the EU after 2014 triggered investigation of contributory factors and control options in poultry production. Reconsideration of the five current target serovars for breeding hens showed that there is justification for retaining Salmonella Enteritidis, Salmonella Typhimurium (including monophasic variants) and Salmonella Infantis, while Salmonella Virchow and Salmonella Hadar could be replaced by Salmonella Kentucky and either Salmonella Heidelberg, Salmonella Thompson or a variable serovar in national prevalence targets. However, a target that incorporates all serovars is expected to be more effective as the most relevant serovars in breeding flocks vary between Member State (MS) and over time. Achievement of a 1% target for the current target serovars in laying hen flocks is estimated to be reduced by 254,400 CrI95[98,540; 602,700] compared to the situation in 2016. This translates to a reduction of 53.4% CrI95[39.1; 65.7] considering the layer-associated human salmonellosis true cases and 6.2% considering the overall human salmonellosis true cases in the 23 MSs included in attribution modelling. A review of risk factors for Salmonella in laying hens revealed that overall evidence points to a lower occurrence in non-cage compared to cage systems. A conclusion on the effect of outdoor access or impact of the shift from conventional to enriched cages could not be reached. A similar review for broiler chickens concluded that the evidence that outdoor access affects the occurrence of Salmonella is inconclusive. There is conclusive evidence that an increased stocking density, larger farms and stress result in increased occurrence, persistence and spread of Salmonella in laying hen flocks. Based on scientific evidence, an impact of Salmonella control programmes, apart from general hygiene procedures, on the prevalence of Campylobacter in broiler flocks at the holding and on broiler meat at the end of the slaughter process is not expected.