In this study, 117 S. Typhimurium isolates from humans (n = 43), food (n = 48) and swines (n = 26) in Brazil were compared using genomic analyses, such as phylogenetic neighbor joining, orthologous protein clusters detection, MLST analysis, and blast identification of prophages and resistance genes related to efflux pump.
The constructed neighbor joining trees grouped the 117 S. Typhimurium isolates into two and three groups, respectively. The majority of the 117 S. Typhimurium strains studied were grouped in a single cluster (≅ 90%) by the genome core (cgMLST) and (≅ 77%) by single copy marker genes (ggTree) suggesting that there is one prevalent subtype that has been successful in contaminating human, food and animal sources for 30 years in Brazil (Fig. 1 and Fig. 2). It is important to mention that the present study provided additional information about S. Typhimurium strains isolated from humans, food and swines in Brazil because such strains have never been studied together.
According to Jensen, homologous genes can be divided into orthologous and paralogs [17]. Orthologous genes are originated from a common ancestor during speciation events and kept the same function, on the other hand, paralogs genes are originated from duplication events and do not kept the same function [18].
Therefore, the OrthoVenn2 is a web server capable to annotate and compare orthologous protein clusters from the whole genome among different species [18]. In the present study, S. Typhimurium genomes more related to LT2, 14028S and D23580 references were compared and had their unique protein orthologous clusters revealed (Fig. 2). All S. Typhimurium isolates used for LT2 – comparison 1, 2, 3 and 4 were composed by a subcluster containing only strains of ST19 isolated from humans in the São Paulo State before the 1990s. The S. Typhimurium isolates used for 14028S – comparison 5 was formed by a subcluster containing only strains of ST19 isolated from food in the Rio Grande do Sul, Santa Catarina and Bahia States between 2006–2012. The S. Typhimurium isolates used for D23580 – comparisons 6 and 7 were formed by a subcluster containing strains of ST313 and ST19 isolated from humans and food in the São Paulo and Paraná States between 1995–2010.
Moreover, the different orthologous protein clusters found are involved in metabolic and regulatory processes, such as transposition, DNA replication, cell adhesion, formate oxidation, trehalose transport, lyase activity and response to mercury ion. These results show that despite being of the same serovar there are unique orthologous protein clusters in the strains studied in comparison to the reference strains and that may be important for their pathogenicity and were kept in these S. Typhimurium strains during natural selection and adaptation (Table 2).
In this study, MLST was performed only for swine isolates, because the STs for humans and food isolates were previously described in [19]. Of the 26 S. Typhimurium strains isolated from swine studied, 16 (61.5%) belonged to the ST19, three (11.5%) to the ST1921, two (7.6%) to the ST14, two (7.6%) to the ST64, one (3.8%) to the ST516, one (3.8%) to the ST639 and one did not have its ST detected.
According to Almeida and collaborators, the ST19 was the most common ST found, the ST313 was the second most prevalent. The STs 1649, 3343 and 1921 were also detected among the S. Typhimurium strains isolated from humans and different foods in Brazil [19].
S. Typhimurium isolates from swine showed greater diversity in the seven housekeeping genes studied despite having a lower number of strains (n = 26) in comparison to the number of S. Typhimurium strains isolated from humans (n = 43) and food (n = 48). In strains from swine the ST19 was the most common, the ST1921 was the second most prevalent and the STs 14, 64, 516 and 639 were also detected.
Accessing the Enterobase for Salmonella, it was observed that 29,572 samples isolated from human, reptile, ovine, swine, poultry, food and bovine in France, Mexico, China, Germany, Scotland, Portugal, Qatar, Korea, Ireland, United States (US), United Kingdom (UK) and Denmark deposited belonged to the ST19. For the ST313 3,049 samples isolated predominantly from human in Kenya, Ethiopia, Zimbabwe, Malawi, Mali and Nigeria were deposited [10].
Moreover, the STs 1649, 3343 and 1921 were found for 16, 4 and 7 strains respectively, isolated from humans, livestocks, food and swines in Venezuela, Ireland, US, UK, Colombia, Ecuador, Vietnam and Brazil [10]. Finally, the ST516, ST64, ST639 and ST14 were linked to 370, 3,850, 237 and 2,149 strains respectively, isolated from humans, poultry, food, aquatic animals, reptiles and the environment in the US, Mexico, Senegal, Germany, Portugal, Qatar, Canada, UK, India, Ghana, Thailand, Malaysia, Malta, Vietnam, Pakistan, Greece, France, Germany, China, Denmark, Scotland, Norway and South Korea [10].
In the present study, the Gifsy-2 prophage was detected in 79 (67.5%) S. Typhimurium isolates, Gifsy-1 in 63 (54%), Salmon 118970_sal3 in 46 (39%) and Haemop - HP1 in 21 (18%). Specifically, Gifsy-1 prophage was detected in 10 (23%) S. Typhimurium strains isolated from humans, 34 (71%) strains isolated from different foods and in 19 (73%) strains isolated from swine (Table 3). Gifsy-2 prophage was detected in 37 (86%) S. Typhimurium strains isolated from humans, 33 (68.7%) strains isolated from foods and 9 (34.6%) strains isolated from swine (Table 3).
It is important to be mentioned that Gifsy prophages carry genes that are related to virulence of S. Typhimurium in the host [20, 21]. The Gifsy-1 prophage encodes three genes involved in the intracellular survival of Salmonella spp. in the host, denominated gogB (leucine-rich repeat protein), sarA (anti-inflammatory response activator) and pagK2. In the same way, the Gifsy-2 prophage encodes a superoxide dismutase (sodC1) that contributes to the survival of Salmonella spp. destroying the toxic radicals of the host macrophages [12, 22]. It is important to emphasize that Gifsy prophages have been found only in S. Typhimurium strains, as well as the Fels-1 and Fels-2 prophages [21]. In the present study, Fels prophages were detected in four S. Typhimurium strains isolated from humans, two strains isolated from food and four strains isolated from swine (Table 3).
According to Brussow et al., the Fels-1 prophage encodes the sodCIII and nanH genes related to the production of superoxide dismutase and neuraminidase in S. Typhimurium, respectively [23]. Furthermore, the Fels-2 prophage carries genes that are apparently related to regulation and adhesion of S. Typhimurium to host cells [12].
The Gifsy and Fels prophages have already been described in S. Typhimurium isolated in various parts of the world, such as Australia, Europe, China, among others [24–26]. It is important to emphasize that other prophages were also found in the S. Typhimurium strains studied including Salmon 118970_sal3 and Haemop - HP1 (Table 3). Moreover, two dozen other prophages were detected in the S. Typhimurium strains studied, but there is less information about them related to pathogenicity and/or virulence of this serovar (Table 3).
In addition, S. Typhimurium isolates from swine showed 6 (23.1%) unique prophages despite having a lower number of strains analysed (n = 26) in comparison to S. Typhimurium strains isolated humans (n = 43) and food (n = 48) that presented 7 (16.3%) and 3 (6.25%) unique prophages, respectively, suggesting the greater diversity in these mobile genetic element for S. Typhimurium strains isolated from swine in Brazil (Table 3).
Resistance to multiple drugs in bacteria has been a serious public health problem worldwide [27]. It is known that there are four main mechanisms that can cause this resistance, such as target alteration, drug inactivation, decreased permeability and drug expulsion through the production of efflux pump [28].
In the present study, the acrA, acrB, macA, macB, mdtK, emrA, emrB, emrR, tolC, mdsA, mdsB, mdfA and cmlA1 genes were detected among the S. Typhimurium strains isolated from humans, food and swines. All of the isolates contained the acrA, acrB, macA, macB, mdtK, emrA, emrB, emrR and tolC genes (Table 4). Other genes related to production efflux pump, such as oqxAB and floR were previously reported in [15].
The AcrAB efflux system has been described as responsible for the intrinsic resistance to many antibiotics that can be used in medical practice for the treatment of S. Typhimurium, such as fluoroquinolones and beta-lactams [16]. According to the World Health Organization (WHO), Salmonella spp. was described as a high priority category pathogen in fluoroquinolones resistance of the Global Priority Pathogens List [27].
The macA and macB genes encode proteins that characterize an efflux pump related to macrolides resistance [29, 30]. According to the Universal Protein Resource (UniProt), the mdtK, emrA, emrB, emrR, mdsA, mdsB, mdfA and cmlA1 genes encode mainly proteins involved in multidrug efflux transporter and confers resistance to different antibiotics, such as aminoglycosides, tetracyclines, novobiocin, nalidixic acid, chloramphenicol and norfloxacin [30, 31]. Furthermore, the tolC gene has been described as important for the formation of some multidrug efflux systems (AcrAB, MacAB, EmrAB and MdsAB) in S. Typhimurium [31].