Skip to main content

Advertisement

Association of Salmonella virulence factor alleles with intestinal and invasive serovars

Article metrics

  • 1120 Accesses

  • 1 Citations

Abstract

Background

The role of Salmonella virulence factor (VF) allelic variation in modulating pathogenesis or host specificity has only been demonstrated in a few cases, mostly through serendipitous findings. Virulence factor (VF) alleles from Salmonella enterica subsp. enterica genomes were compared to identify potential associations with the host-adapted invasive serovars Typhi, Dublin, Choleraesuis, and Gallinarum, and with the broad host-range intestinal serovars Typhimurium, Enteritidis, and Newport.

Results

Through a bioinformatics analysis of 500 Salmonella genomes, we have identified allelic variants of 70 VFs, many of which are associated with either one of the four host-adapted invasive Salmonella serovars or one of the three broad host-range intestinal serovars. In addition, associations between specific VF alleles and intra-serovar clusters, sequence types (STs) and/or host-adapted FimH adhesins were identified. Moreover, new allelic VF associations with non-typhoidal S. Enteritidis and S. Typhimurium (NTS) or invasive NTS (iNTS) were detected.

Conclusions

By analogy to the previously shown association of specific FimH adhesin alleles with optimal binding by host adapted Salmonella serovars, lineages or strains, we predict that some of the identified association of other VF alleles with host-adapted serovars, lineages or strains will reflect specific contributions to host adaptation and/or pathogenesis. The identification of these allelic associations will support investigations of the biological impact of VF alleles and better characterize the role of allelic variation in Salmonella pathogenesis. Most relevant functional experiments will test the potential causal contribution of the detected FimH-associated VF variants in host adapted virulence.

Background

In recent years, it has been shown that adaptation of pathogen to certain niches is accompanied by genome degradation [1]. One of the models used for this observation is Salmonella enterica, a versatile Gram-negative Enterobacteriaceae, that includes over 2,500 different serovars and accounts for a vast number of human and animal infections [2, 3]. S. enterica strains causing human disease are divided into human-restricted typhoidal serovars (Typhi and Paratyphi) and nontyphoidal (NTS) serovars (e.g. Enteritidis, Typhimurium, Newport and most others S. enterica serovars). The latter serovars are designated generalists, since they also infect a broad-range of mammals and birds, causing gastroenteritis and/or persistent infections with potential for long-term shedding. Animals have also more invasive serovars capable of extra-intestinal spread responsible for septicemia, such as S. Gallinarum, which is restricted to poultry and many wild birds [4], and S. Dublin and S. Choleraesuis, that are mostly adapted to cattle and pigs, respectively, with rarer human or other mammalian infections. Some Salmonella serovars can cause invasive NTS infections (iNTS) by septicemic spread [5], and are often not associated with diarrhea. In the last decade, such strains were isolated mostly from Africa, frequently in malaria or HIV-infected patients. Genomic studies have shown that they are mainly caused by specific lineages of S. Typhimurium [6] and S. Enteritidis [7].

From an evolutionary point of view, the most parsimonious explanation for the host-restricted serovar is that they emerged from an ancestor with a broad host range [8]. The Salmonella adaptation process to a reduced number of hosts and a more efficient escape from the microbial competitive intestinal milieu by invading deeper tissues resulted in the loss of genes by gradual pseudogenization and subsequent deletions. Nuccio et al. demonstrated that 469 genes of the central anaerobic metabolism network underwent some degree of degradation in extraintestinal Salmonella serovars [9]. During the re-annotation of 15 S. enterica genomes they found that 471 intact genes had been misidentified as pseudogenes. Despite progress in pseudogene annotation, the precise mechanisms associated with host adaptation remain mostly unexplained. One infrequently addressed caveat with pseudogenes is that some of them still express full proteins due to frameshifting or t-RNA-mediated suppression over premature stop codon, particularly relevant possibilities when only low levels of gene product is required for full function [10, 11].

In addition to pseudogenization and gene loss, host adaptation in Salmonella has also been influenced by gene gains via horizontal genetic transfer, such as the acquisition of virulence plasmids and pathogenicity islands that encode genes used for adherence, host cell invasion, induction of innate immune and/or pathophysiological responses [8, 12, 13]. More subtle effects can result from point mutations that modify DNA structure and recognition sites for regulatory binding proteins or cause amino acid substitutions in open reading frames that might contribute to host adaptation of S. enterica serovars [14]. Overall, adaptation to human or animal species is a complex phenotype [15] and the relative importance of the participation of gene loss, gain or modification for a preferential intestinal or systemic pathotype remains understudied.

Early investigations using multi locus enzyme electrophoresis and multilocus sequence typing (MLST) already revealed that in contrast to generalists, host-adapted Salmonella serovar showed reduced allelic variation in their house-keeping genes [16]. Partial core genome analysis of several generalist and host-adapted Salmonella serovars was consistent with the separation of clear lineages distinguishing generalists and specialists [17]. Core or whole genome analyses have been very useful for phylogeny and epidemiological typing [18,19,20], but their potential to identify critical allelic variants participating in host preference and virulence remains obscured by an excess of mutations that may not be adaptive or relate to other phenotypes, such as survival in different non-host environments. Relevant Salmonella properties linked to intestinal immunopathology, septicemic capacities and/or adaptation to certain host species are more likely enriched among allelic variants of virulence genes such as surface or exported proteins that interact with host cells and molecules [14]. Notably, we previously showed that a distinctive genomic property of S. Typhimurium’s association with diverse hosts related at least partially to allelic variants of surface or exported proteins, as established with some fimbrial adhesins [21, 22]. Allele-specific adhesive properties confirmed the biological significance of the detected associations. For example, the FimH allele most frequently present on bovine S. Typhimurium strains mediated best bacterial binding to bovine intestinal epithelial cells, whereas bacteria expressing a FimH allele associate most frequently with human isolates adhered preferentially to human intestinal epithelial cells. Comparable associations were determined for the host-adapted binding properties of the FimH alleles from other intestinal and invasive serovars, such as Typhi and Choleraesuis. Further studies showed that even though strain-specific combination of alleles for different adhesins of the broad host range serovar Newport mediated preferential binding to one host species, they kept at least one allele mediating binding to cells from another host. These findings suggested that allelic variation in multiple adhesins of a generalist serovar can contribute to bacterial adaptation to certain preferential hosts without losing the capacity to maintain a broad host range [22].

Since Salmonella evolution to increased levels of host-adaptation may accompany more aggressive and invasive pathotypes, we hypothesized that not only adhesins, but also other proteins involved in virulence might have allelic variants that are associated with either intestinal or septicemic Salmonella serovars or strains. To examine this issue, we collected the sequences of 70 virulence genes from 500 Salmonella enterica subsp. enterica strains. We deduced the protein sequences of the genes and compared the alleles. To evaluate the potential involvement of allelic variation in the more aggressive and invasive S. Typhimurium and Enteritidis iNTS strains, we analyzed them separately. In support of our hypothesis, we found that not only specific degraded/absent virulence gene, but also certain allelic virulence protein sequences were linked to the intestinal or invasive lifestyle of the various Salmonella serovars or lineages. These findings narrow down the list of genes and allelic variants that may relate to host adapted pathologies worth studying for host-determined effects. The confirmation of previously described virulence levels associated with specific allelic or absent virulence factors [21, 23,24,25,26] supports the potential biological relevance of the new host-virulence factor allele associations detected in this study.

Results

Collection of virulence factors from intestinal and systemic Salmonella serovars

For this study, 70 known VFs from 500 Salmonella genomes were collected and their protein sequences were deduced to analyze and compare their alleles. Selection of the VFs (Additional file 1: Table S1) was based on their presence in most of the serovars studied and their capacity to interact directly with host molecules such as known adhesins and effector proteins exported by the two Salmonella type 3 secretion systems (T3SS) [27,28,29]. We included in the list outer membrane proteins and LPS modifying proteins that are regulated by the PhoPQ two component system which modulates the surface of intracellular Salmonella [30, 31]. The analysis included proteins from four host-adapted or restricted invasive serovars and three typical broad host range intestinal serovars (Additional file 2: Table S2).

Host-related VF gene degradation

The most host-adapted invasive and pseudogenized S. Typhi had also the most disrupted or absent VF genes from a collection of 70 selected virulence genes (27.1%) (Additional file 3: Table S3 and Fig. 1). The VF genes of this human restricted serovar underwent even a greater selective pressure than the corresponding genome itself (5-10%)[9, 32]. Similarly, the predicted reduction in VFs for the other three invasive Salmonella serovars (Gallinarum, Choleraesuis and Dublin) was significantly higher than that for the three intestinal Salmonella serovars (9.5% versus 2.0%; p < 0.0001, two-tailed Fisher’s exact test). Thus, in addition to the inactivation of genes for unnecessary metabolic pathways used for intestinal colonization [9, 32], many VFs of intestinal Salmonella generalists are apparently superfluous and possibly deleterious for invasive-septicemic infections. In addition, the more host-restricted serovars Typhi and Gallinarum had less than half as many predicted virulence genes as the host-adapted serovars Choleraesuis and Dublin, including the SPI-2 effector proteins SseI/SrfH, GtgE, SopA, SseK2, and the outer membrane proteins RatB and SadA (Additional file 3: Table S3). Concurring findings were obtained with increased tendency for VF gene degradation of iNTS versus non-invasive NTS from serovars Typhimurium and Enteritidis (Fig. 2a and b), in agreement with earlier whole genome studies [7, 33]. By including only one sub-Saharan representative of each serovar to minimize same strain/lineage over-representation [6, 7] statistically significant differences in gene degradation between iNTS and non-invasive NTS was found for SseI, SspH2, ShdA, SadA and SiiE (p < 0.02) in S. Enteritidis, and SseI, ShdA and SiiE (p < 0.01) in S. Typhimurium (Fig. 2a and b, Additional file 3: Table S3). S. Typhimurium iNTS of the ST313 lineage showed also increased degradation for SspH2 and SadA as compared with non-invasive NTS in general (p < 0.0001). All sub-Saharan S. Enteritidis iNTS strains lacked MisL, unlike all the other iNTS strains (p < 0.0001) and all non-invasive NTS strains, suggesting that they may form a separate lineage. Although Salmonella pseudogenization has generally been recognized as a successful evolutionary adaptation to increased virulence for fewer hosts as well as reduced or lost virulence for generalists [9, 32], as shown for SseI in specific infection models [23, 34,35,36], the biological relevance of most pseudogenized or missing genes remains to be demonstrated.

Fig. 1
figure1

Virulence factor predicted to be missing due to disrupted or absent genes among the genomes of 7 Salmonella serovars and some of their lineages. Serovars and lineages are indicated on the top and and VFs on the right. The heatmap shows genes present in most (> 80 %, in black), some (≤ 80%, in grey), or no genomes (white) due to deletion or pseudogenization

Fig. 2
figure2

Distribution of the virulence factor alleles most frequently degraded among iNTS and non-iNTS. a S. Typhimurium and (b) S. Enteritidis genomes. VFs are shown on the top and strains on the right side. The heatmap shows the numbers of different alleles for each gene product for each strain represented in black (one allele), different gray gradations (2 to ≥ 5 alleles), or white (absent gene product)

Allelic variation among the 70 VFs

DNA sequence changes were evaluated in 70 genes encoding VFs and the corresponding deduced protein sequences were used for all the comparative analyses of this study. Protein alignments identified allelic protein variants that were assigned incremental numbers (Additional file 4: Table S4). Comparisons were based on protein and not gene sequences to detect potentially functional associations with different pathotypes of the Salmonella serovars and lineages. The number of amino acid substitutions in the deduced protein sequences was compared between each allele and the dominant allele of the corresponding 70 proteins. We calculated Hamming distances (number of amino acid changes between the most common allele and other alleles for each of the studied protein) and found the PhoPQ-regulated proteins were the least diverse (Additional file 5: Figure S1). As can be seen from the figure, the most variable VFs were StfH, BapA, SiiE, RatB, FliC, SteC, SlrP and ZirS. The most diverse proteins, as determined by the Simpson’s Diversity index, were StfH, SlrP, SteA, ZirS, SteC, FliC, SspH2, SseL, RatB, BapA and SiiE (Additional file 6: Table S5). The number of alleles varied from 4 to 98, with a mean value of 9.4 ± 2.1 per protein (confidence level 95%; median value of 16). Most proteins with the highest numbers of alleles were high molecular weight membrane proteins such as SiiE (595 kDa), RatB (206 kDa) and BapA (386 kDa). The ratio of VF alleles to the total residues per VF varied from 0.8% (SadA) to 15% (SssA), with a mean of 5.5% (Additional file 6: Table S5). On average, 11.9% of the residue positions varied for all VF, with a median value of 7.4%. Nearly a third (32.2%) of all the VF were represented by one dominant allele, defined as including at least 50% of the strains (Fig. 3). Notably, adhesins FimH, BcfD, StfH, SteD, SinH, BapA, SadA, MisL and SiiE were more variable than other VFs, as determined by genetic diversity h [37] with a mean of 35.6 for the adhesin alleles and 16.7 for the other VF alleles (p = 0.08, Mann–Whitney U test; Additional file 7: Table S6). A few proteins had two major alleles, such as the fimbrial adhesins BcfD, StfH and FimH, the large membrane proteins ShdA, RatB, BapA, CigR and SiiE, the secreted proteins SlrP, SseI, PipB, SifA, SpiC, SseC, SteA, SrfA, SteC, SopA, SseL, SopE and ZirS, and the outer membrane protein OmpD (Fig. 3). Some proteins (21%) had a most frequent allele for less than 20% of the strains, including the high molecular weight membrane-associated protein ShdA, SiiE, BapA and SadA. In summary, Salmonella VFs demonstrate a significant level of allelic variation.

Fig. 3
figure3

Allelic frequencies for 70 studied virulence factors among 500 Salmonella. VFs are organized in six groups as labeled on the bottom of the graph (PhoPQ-regulated, “independent”, SPI1 and/or SPI2 effectors, adhesins/outer membrane proteins). Each VF has several alleles whose number of strains per allele are each represented by column heights, as labeled on the right axis. The most frequent to the least frequent allele (from allele 1 to > 20) are shown from the back to the front of the graph, each in a different color (as best visualized with the first VF, EptA, on the left)

VF alleles in intestinal and invasive Salmonella serovars

Even though dominant alleles of most VFs were shared almost equally between the intestinal and invasive serovars of Salmonella (Additional file 8: Figure S2), some dominant protein alleles were only present in the former (StfH, SseI, SifB, SopA, SseK2, ShdA, BapA, SadA, SiiE and SopE) or latter serovars (FimH, SlrP, SteA, SrfA, ZirS, SteC, FliC, SspH2, SseL, RatB and MgtB). Of the 70 proteins studied, 58 (83%) had alleles that were shared between two or more serovars. In contrast, 16 proteins had serovar-restricted alleles, the adhesin FimH, the flagellin FliC, the membrane-associated proteins ShdA, RatB, SadA, SiiE and BapA, the PhoPQ-regulated protein EptA, and the secreted proteins SteA, SrfA, SteC, SopA, PipB2, SseK1, SseK2 and ZirS. The alleles of the phage-encoded SopE were shared only between intestinal serovars. Allelic variation of each protein within one serovar was evaluated by calculating h values (Additional file 7: Table S6). Comparisons of monoallelic proteins (h = 0) suggested that the intestinal serovars were twice as diverse as the invasive serovars (Additional file 7: Table S6, mean number of monoallelic VFs = 30.3 versus 13.2, respectively; p = 0.01, Mann–Whitney U test). Moreover, dominant protein alleles from either intestinal or invasive serovars with h values below 0.1 confirmed a corresponding trend (mean number of VFs with h < 0.1 = 43.2 versus 32.7 alleles; p = 0.25 Mann–Whitney U test). Even though host-restricted S. Typhi had over three alleles for 17 VFs (Additional file 7: Table S6), it was the least diverse serovar with the largest number of monoallelic proteins (65.7%), followed by the other invasive serovars (38.7%) and the intestinal serovars (18.9%). Taken together, these results support a Salmonella evolutionary model whereby the invasive serovars are younger than the intestinal serovars, with the accumulation of allelic variation being time dependent. In addition, purifying selection linked to host adaptation most likely participated in restricting allelic variation.

Among all the VFs, SteC and FimH had the highest Simpson’s diversity indices with values close to one (> 0.9, Additional file 6: Table S5). Linkage disequilibrium (LD) was evaluated by determining indices of associations (IAS) and mean gene product diversity (H) for each serovar, Newport lineages and Gallinarum biovars (Additional file 7: Table S6). LDs correlated with gene product diversities. Diversity increased from serovar Typhi to serovars Newport, consistent with the previously described polyphyletic nature of the latter serovar [38, 39]. Interestingly, biovar Pullorum showed greater VF allelic variability than biovar Gallinarum; this might be related to a broader host range, particularly in wild birds, of the former biovar [40].

To evaluate whether alleles of the 70 VFs can be used to identify subdivisions among all Salmonella studied irrespective of serovars, a range of assumed population groups (K = 2-13) was investigated, and K = 10 was found to be the best estimate to discriminate serovars and lineages with the Structure v. 2.3.4 program [41,42,43] (Additional files 9 and 10: Figures S3 and S4). In addition to recognizing the individual serovars, the program distinguished the two known Newport lineages and Gallinarum biovars at this K value, with a separate cluster for strains in lineage III of S. Newport. A minimum spanning tree obtained with goeBURST located the intestinal serovars in the central nodes and the invasive serovars on tip nodes (Additional file 11: Figure S5), and thus best illustrated how the invasive serovars represent dead-ends in pathogen evolution as they relate to host adaption. In summary, VF variants were found associated with unique sets of flagellar and O-antigens by separating the seven serovars and to further divide two serovars in their biovars or lineages, in agreement with their distinct evolutionary host adaptations.

Hierarchical clusters of Salmonella VF alleles

Hierarchical clustering analyses were performed to determine the distribution of VF alleles among the different serovars. The colored heat map showing each VF’s 20 most dominant alleles irrespective of serovars not only confirmed that in contrast to all other serovars Typhi has more deleted or pseudogenized VF genes, but also showed that Typhi has unique alleles for 66% of its VFs (Fig. 4). Significant variation was only detected in the repetitive bacterial immunoglobulin domains of the giant Ca-binding inner membrane adhesin SiiE [44]. The heat map also illustrated the highly clonal property of S. Typhi, in agreement with earlier whole genome studies [45].

Fig. 4.
figure4

Heatmap showing hierarchical clustering of 500 Salmonella strains with corresponding alleles across the 70 virulence factors studied. The seven serovars each with 52-75 strains are shown on the left of the heatmap, VFs on the top, arrayed in the order genes are located on the chromosome. Clade numbers, biovars (S. Gallinarum) and lineages (S. Newport) are shown on the left. Colors indicate different alleles, with black being the dominant allele for each VF among all 500 Salmonella; blue-green, the 2nd most frequent allele, etc. up to the 20th (and more) most frequent allele in purple, as shown on the top left (1-20); missing VFs are shown in white

S. Enteritidis was divided in two clusters of different alleles for the seemingly unrelated proteins SifB, ZirS, BapA, SptP and SopD (Fig. 4). The smaller cluster 2 consisted of half of the Malawian iNTS strains with a signature allele for the T3SS effectors SptP, SifB, SopD and ZirS, and the PhoPQ-regulators LpxO and CptA (p < 0.05). The other Malawian iNTS and the non-African iNTS isolates belonged to cluster 1, suggesting that separate and/or additional routes of genomic evolution can lead to an invasive phenotype with the likely participation of a weakened host immune system related to concomitant infections or malnutrion [46].

For S. Dublin, hierarchical clustering detected minor cluster 3 of human isolates from Uruguay that was characterized by unique VF alleles of SipD and LpxR (Fig. 4). MLST did not distinguish these strains in earlier studies [47, 48].

Two major S. Typhimurium groups could be distinguished by hierarchical clustering with one group subdivided in three smaller clusters (Fig. 4). The majority of S. Typhimurium iNTS lineage consisted of the African iNTS isolates (ST313 and ST2080), represented by cluster 1. All the ST313 strains had different alleles for SopD2, GtgA, SseC, SopA and SipB (p < 0.03). The other few iNTS were dispersed in cluster 2, 3 and 4.

Hierarchical clustering clearly separated S. Gallinarum biovar Gallinarum from biovar Pullorum with separate alleles for 23 VFs (StfH, FimH, OmpX, SopD2, GtgA, PipB, SopB, SifA, PagD, SseB, SseC, OmpD, SteA, SrfA, SifB, ZirS, SteC, PagL, PipB2, AvrA, SipA, SopD and CigR), as well as the presence or absence of STM3031and SseK1 (Fig. 4). Biovar Pullorum was further divided in two clusters that had different alleles for 11 VFs (BcfD, StfH, SpiC, SseB, SseC, SseG, SifB, Ais, PipB2, SipC and SopD), whereas SrfA was absent in one cluster. Cluster 2 of this biovar shared several alleles with biovar Gallinarum (BcfD, SpiC, SseG and Ais).

S. Choleraesuis was divided in two hierarchical clusters distinguished by FimH, StfH, SifA, LpxR, SifB, IagB, IacP, SipB, SpaR, SopD and MgtB (Fig. 4). The two clusters appeared to not overlap with the annotation for the two major biovars of Choleraesuis, senso stricto and Kunzendorf (Additional file 2: Table S2). However, since a recent study discovered that there were many errors for the biovar annotations and that all fimH- strains were biovar senso stricto [49], our data suggest that all biovar senso stricto strains belong to cluster 2.

For S. Newport, hierarchical clustering of VF alleles detected two major groups (Fig. 4) that corresponded to the known lineages II and III of this serovar [38]. Lineage II included an additional discrete cluster of Asian isolates recognizable by a novel StfH allele and separate alleles for SrfN, SteA, RatB, SinH, AvrA, MisL, CigR and SopE. Lineage III had a separate cluster of strains that were characterized by specific alleles for FimH, SlrP, PagC, LpxR, SteB, ZirS, SteC, Udg, SinH, BapA, SadA and PhoN, and by the presence of AvrA. As expected, none of the strains studied related to the rare Europe-associated S. Newport lineage I [38, 50].

Taken together these results showed that hierarchical clustering of Salmonella virulence protein sequences complements whole genomes studies by proposing the identification of new minor strain clusters of potential pathobiological relevance.

VF alleles associations with STs, adhesin alleles and hosts

Several Salmonella serovars and some intra-serovar strain lineages or sequence types are known to demonstrate different levels of host specificities [14, 51, 52]. We previously reported that Salmonella strains and serovars with particular adhesin alleles that are more frequently isolated from certain hosts than others typically mediated better bacterial binding to enterocytes from these hosts [21, 22]. Here we suggest parallel host-associations for VF alleles as determined by their genomic ties with alleles of the FimH adhesin. Since an excessive number of the randomly selected strains from this study lacked the required metadata, specifically the source of isolation (Additional file 2: Table S2), we relied on our previously described adhesin allele-host associations to propose potential roles of VF alleles in host-specific virulence.

Two serovars, Typhi and Enteritidis, had no significant FimH allelic variation and thus no potentially possible host associations with VF alleles. S. Typhi had only one FimH allele for all 75 genomes evaluated, and similarly, S. Enteritidis had one dominant FimH allele (92% of the strains) with only four new minor alleles mostly in Malawian strains that were not associated with STs or other VF alleles (Additional file 2: Table S2 and Additional file 12: Figure S6). This was consistent with the host restriction or preference of both serovars, Typhi as a human pathogen and Enteritidis as an effective bird colonizer and poultry-derived food product contaminant responsible for human infections.

S. Dublin isolates were mainly ST10 strains (93.3% of all isolates) (Additional file 2: Table S2). A newly observed second FimH allele (30.7% of 75 strains; Additional file 12: Figure S6) associated with an allele of SseC (Fig. 5) with no detectable host preference.

Fig. 5
figure5

Degree of association of FimH alleles with other virulence factor alleles for five Salmonella serovars. VF alleles are indicated on the top of each heatmap while on the left the major FimH alleles (Additional file 7: Table S6) to which they are associated in S. Newport, S. Gallinarum, S. Choleraesuis, S. Typhimurium, or S. Dublin. Different alleles are shown in different colors (allele 1 to 4, blue, yellow, green and brown, respectively). Missing VFs are shown in white. Numbers indicate the prevalence of the most dominant allele for the corresponding serovar/lineage

Most S. Typhimurium were either non-invasive NTS ST19 or iNTS ST313 strains (Additional file 2: Table S2) with mainly the human-associated FimH alleles [21]. A few ST19 strains formed cluster 4 (Additional file 12: Figure S6 and Additional file 2: Table S2) and had the FimH allele previously found to be preferentially associated with bovines and to bind best to bovine intestinal epithelial cells [21]. These strains also had unique alleles for the surface-exposed proteins BcfD, BapA, ShdA, RatB and MisL (Fig. 5), suggesting that like FimH, the five other proteins might participate in some level of adaptation to the bovine host. In addition to the previously described five FimH alleles one new minor allele was discovered (Fig. 5, Additional file 2: Table S2).

In addition to different sequence types (Additional file 2: Table S2), S. Gallinarum biovars Gallinarum and Pullorum had different alleles for 44% of the VFs. Strains of biovar Gallinarum had either one of three FimH alleles (Additional file 2: Table S2), even though most of their other VFs were represented by only one allele. In contrast, biovar Pullorum had a fourth FimH allele with mostly different VF alleles than biovar Gallinarum (Fig. 5, and Additional file 12: Figure S6), which likely relates to the different pathogenesis of these bird-restricted biovars [53]. Although all discovered FimH, StfH, and BcfD adhesin alleles were not shared across serovars, the BcfD allele 1 was an exception by being shared between both biovars of S. Gallinarum and S. Enteritidis, potentially participating in the concurring affinity of these serovars for avian hosts [54].

The majority of S. Choleraesuis of known origin were porcine isolates that included mainly ST145 and ST66 strains (Additional file 2: Table S2 and Additional file 13: Figure S7). All S. Choleraesuis strains of cluster 1 (Additional file 12: Figure S6) had either one of the two previously reported FimH alleles, the porcine-associated FimH104 (here Fim8) or the human-associated FimH105 (here FimH9) [21]. Even though there were very few confirmed human isolates in the current study to corroborate these associations [21], the two FimH alleles were each present with a different SopD allele (Fig. 5). Cluster 2 strains (21.6%) consisted of both of human and porcine isolates lacking fimH and typically associated with specific alleles, as shown in Fig. 4 (Additional file 2: Table S2).

MLST delineated S. Newport lineages II and III, in agreement with the previously described host [38, 55, 56] and allelic adhesin associations [22]. The two lineages were also differentiated by distinct alleles for many VFs (Additional files 11 and 12: Figures. S5 and S6, Additional file 2: Table S2). Consistent with previous findings, lineage II was mostly populated by strains with the three adhesin alleles A/A/A1 for FimH/BcfD/StfH whereas the lineage III strains had mainly alleles B/B/B1 for these adhesins [22]. Lineage II formed the largest cluster that mainly consisted of ST45 strains which were typically US isolates. A few ST31 and ST46 strain (together with its single locus variants ST157, ST2364) were from Asia (Additional file 12: Figure S6 and Additional file 2: Table S2), with unique alleles for ST31 (MisL, SseL, PipB2, CigR, SteA) and all these Asian strains (StfH, SinH, SiiE, RatB, AvrA, LpxO). Lineage III had two separate clusters of VFs mainly populated by either ST5 (21 strains in cluster 1) or ST118 strains (16 strains in cluster 2). The small cluster 2 of lineage III included strains with distinct set of alleles and the consistent presence of AvrA (Additional files 11 and 12: Figures S5 and S6, Additional file 2: Table S2).

Finally, to determine whether VF generally show more allelic variation than other gene products, we evaluated the protein sequences of 40 other gene products chosen to be only indirectly associated to virulence as a stringent comparative group (Additional file 14: Table S7). Most serovars, lineages or biovars showed significant greater allelic variations among the 70 VFs than their corresponding control group, in support of their potential role in host adaptation (Additional file 15: Figure S8). In addition, the VF alleles that formed different strain clusters and serovars that concomitantly share reported host-associated allelic FimH might represent further protein sequence adaptations of VFs to certain hosts and pathotypes, as suggested for adhesin alleles [14, 21, 22].

Discussion

Having previously confirmed the biological relevance of detected associations of Salmonella adhesin alleles with host species [21, 22, 27], we extended our studies by interrogating the extent of protein sequence variabilities of 70 VFs in 500 genomes from 7 Salmonella serovars. Strains from both intestinal generalists and invasive-septicemic host-adapted serovars were evaluated. Distinct new allelic protein variants and gene losses/pseudogenes were observed to be associated not just with specific serovars or a few STs, as previously described [9, 47, 57], but more broadly with many unique STs and with strain lineages or clusters within one serovar. For example, by screening the genes encoding VFs, we found that shdA and siiE were pseudogenized in all isolates belong to S. Typhimurium ST313 lineage, in addition to confirming the previously described inactivation of the genes for SseI, RatB, and SspH2 for ST313 strains [6]. Whether shdA and/or siiE pseudogenization plays a role in host adaptation or systemic infections remains to be determined. Most significantly, we found that in addition to sadA, sseI and siiE were also pseudogenized in the sub-Saharan African S. Enteritidis iNTS strains, and we confirmed the reported pseudogenes for SspH2, ShdA, and RatB in these strains [7]. More generally, gene degradation was found to be significantly higher for VFs than previously reported with whole genome studies on host restricted and adapted serovars. In addition to pseudogenization, a decreasing VF allelic variability from intestinal broad host range to invasive/septicemic host-adapted serovars highlighted the involvement of virulence factor alleles in evolutionary adaptation to fewer hosts with the concomitant change in virulence and disease severity.

The role of pseudogenization in changing Salmonella into a systemic pathogen has been well described. For example, pseudogenization of sopD2, sopA and sopE2 in S. Typhi increased pathogen invasiveness [24, 26]. Similarly, whereas SseI supported a generalist’s persistence in its host by inhibiting dendritic cell migration [35], its absence participated in the spreading of iNTS strain 313 in the host [23]. However, the latter example might be Salmonella strain and infection model-dependent, since SseI was previously shown to participate in invasiveness [34,35,36].

Although shdA which encodes a Peyer’s patch and cecum colonization factor was shown to be pseudogenized in the typical host adapted serovars [3, 48], a pseudogenized shdA in S. Typhi was shown to be functional [25]. Apparently, shorter proteins and/or bypassed stop codons by translational frameshifts with potential recoding can result in functional proteins [10]. Thus, even though the shdA gene was pseudogenized in most iNTS of serovars Enteritidis (92%) and Typhimurium (87%), and in all but one S. Newport strains, it remains possible that a functional ShdA protein is expressed in theses strains. Alternatively, loss of ShdA-mediated colonization in some strains might be compensated by the activity of one or more other bacterial factors [22].

Each VF showed an extensive range of allelic variants, most being shared between 2 or more separate serovars and only 17% being serovar-restricted. Analysis of the VF allelic variants by three methods separated the serovars from one another, including the Newport lineages and Gallinarum biotypes, suggesting a co-adaptative evolution between the Salmonella VF alleles and its O- and flagellar antigens that define its serovars. Results from the goeBURST program visualized the intestinal generalists as nodes on a central branch with branches radiating separately out of it for each host-adapted invasive serovar, in support of separate adaptation histories for their VF alleles. Hierarchical cluster analysis further separated VF alleles into different clusters. Interestingly, both serovar Typhimurium and Enteritidis iNTS formed separate clusters that included either only iNTS strains, or both iNTS and non-invasive NTS. Hence, iNTS comprise both successful outbreak strains such as ST313, as well as epidemically less efficient iNTS whose invasive classification is likely based on clinical symptoms of sepsis, as related to a weakened immune system of individual patients. Serovar Choleraesuis was clearly separated in two clusters by two main sets of VFs alleles, in contrast to serovar Dublin which had one main set of VFs, with a few strains forming 2 minor clusters. Allelic profiles of these host-adapted Salmonella might relate to host-determined variables such a pig and cattle breeds [58], and differential geographic distribution of livestock breeds. As expected, the two biovars of serovar Gallinarum formed separate clusters, albeit biovar Pulllorum showed more diversity with an additional sub-cluster. However, each biotype had relatively few alleles for each VF, in agreement with their adaptation to the restricted avian host, and their rare isolation from wild birds [59]. S. Newport has the most sets of diverse VF alleles, with lineages II and III each having two separate clusters. One clade and a subclade of lineage II were mainly occupied by Asian isolates from humans or edible non-mammalian vertebrates, suggesting parallel adaptation to different geographical niches with different dietary sources [60].

In addition to the described presence of specific VF alleles in ST313 of iNTS, several major STs of serovar Typhimurium, Gallinarum, Choleraesuis and Newport showed associations with sets of VF alleles. However, the separate evolution of intra-serovar lineages and within some STs was better revealed by the numerous sets of allelic VF variants that are potentially involved in successful adaptations to unique or multiple distinct niches. A physical association between the alleles of three fimbrial adhesin genes in 262 S. Newport isolates previously showed two major sets of allelic variants for most strains, one representative of bovine and porcine isolates, and the other of environmental and human isolates [22]. Each strain had at least two allelic adhesins binding significantly better than the other alleles to an intestinal epithelial cell line from the original host(s), leaving one adhesin allele binding best to corresponding cells from another host. Thus, allelic variation in multiple adhesins of S. Newport was suggested to contribute to bacterial adaptation to certain preferential hosts while still maintaining a large host range [22]. Based on these results and the established genomic link between FimH and other VF alleles for each of the 500 strains of this study, we inferred that some of the VF alleles were also host-associated, and possibly host-adapted. Accordingly, sets of VF alleles, each with an associated host-adapted FimH were found in serovar Typhimurium, Choleraesuis, Gallinarum and Newport.

Although this study hints at the possibility that some of the described serovar, lineage, biotype and clonal associations for VF alleles might play a role in host specificity and pathogenesis, the biological relevance of the detected associations need to be tested for a definite proof, as addressed previously for a few adhesins [21, 22, 61]. It is also expected that a significant number of amino acid substitutions in various alleles are essentially silent with no functional effect on host adaption. Other substitutions might have indirect effects. For example, different SptP alleles between S. Typhi and Typhimurium affect export efficiency of this SPI-1 effector protein [62]. Some of the variable amino acids of the SipD tip protein of S. Typhi and Typhimurium might not represent adaption to different pathogen-host interfaces, but provide recalibrated contacts with variable PrgI needle subunit proteins for structural assembly and function of the type 3 secretion system [63]. Lastly, amino acid replacements might also inactivate the function of the protein and represent a first steps in pseudogene formation.

To optimize the significance of our comparative analyses and minimize biases, 500 genomes were collected randomly, choosing similar numbers for each serovar, biotype or Newport lineages, and for Typhimurium and Enteritidis iNTS and non-invasive NTS, based genome data availability. Taxonomically misidentified genomes were rejected and several serovar Gallinarum biovars were classified by their genomic signatures [64]. Unfortunately, many Salmonella genomes lacked metadata on host origins for the corresponding isolates, with the possibility that some strains were epidemiological clones. However, when clones were suspected as for the African ST313 isolates, only one representative strain was used for comparisons. A main strength of this study was the inclusion of a relatively large number of Salmonella VFs and genomes from several intestinal and septicemic serovars. This allowed us not only to confirm a few previous observations on allelic VF variants, but also to bring attention to VFs whose sequences that associate with specific metadata might play a role in the different pathogenic properties of the various serovars and strain lineages.

Conclusions

In conclusion, the comparative analysis of the allelic variants of 70 VFs from three intestinal/broad range Salmonella serovars and four invasive/septicemic serovars is providing new information on amino acid substitutions that might modulate the function of Salmonella proteins known to play a role in host specificity and pathogenesis. These findings should help to better target wet lab and animal model experiments aimed at identifying allelic variants responsible for causal effects.

Methods

Bacterial genomes

A total of 500 Salmonella enterica subsp. enterica genomes were used for the study: 75 S. Enteritidis, 75 S. Typhimurium, 74 S. Newport, 75 S. Dublin, 74 S. Choleraesuis, 75 S. Typhi, and 52 S. Gallinarum biovars Gallinarum and Pullorum strains (Additional file 2: Table S2). Genomic data with corresponding metadata were downloaded from NCBI RefSeq database. Raw sequencing data with metadata were collected from NCBI SRA or EBI ENA repositories and were assembled using SPAdes genome assembler implemented through Galaxy/GVL server [65]. Sampling was performed randomly with higher assembly level being collected first. With the exception of S. Gallinarum that had only 52 available genomes 75 genomes were collected from each serovar. Incorrectly serotyped strains were determined for one S. Newport and one S. Choleraesuis genome, as determined by their sequences with SISTR [19], which found that 5% of all Salmonella enterica genomes or short reads deposited in publicly available databases were misidentified. The MLST v. 1.8 server was used for the detection of sequence types (ST) [66]. ST and metadata for some of the genomes were obtained from EnteroBase [67].

We collected and translated 70 genes for known and putative Salmonella VFs [29], including adhesins, SPI-1 and SPI-2 effectors and PhoPQ-regulated proteins (Additional file 1: Table S1). A control group of 40 genes, encoding proteins indirectly related to virulence were listed in the same table.

VF allele’s determination

Most protein sequences were identified with TBLASTN (e-value cut-off of 0.001, minimum query coverage >95%), and additional proteins were detected with MegaBLAST, using S. Typhimurium LT2 (Genbank # NC_003197) nucleotide sequence information (e-value cut-off of 0.01, minimum query coverage of >95% and identity >95%). Translation from DNA to protein sequences was performed in MEGA7 [68]. Multiple sequence alignment was performed in BioEdit v. 7.2.6.1 (Ibis Biosciences, an Abbott Co.). The SignalP v. 4.1 server was used for the prediction and location of signal peptide cleavage sites [69]. Variable amino acid residues were extracted from protein alignments using FaBox v. 1.41 online tool [70].

Alleles were defined as protein sequences without signal peptide, which have at least one amino acid change (nonsynonymous single nucleotide polymorphism; nsSNP) or indel enumerated incrementally in order of their frequency. A script was made to number protein allelic variants and calculate Hamming distances (number of amino acid substitutions between sequences) for each allelic protein of the 500 genomes to determine dominant allele(s) [71].

Association analysis and statistics

The standardized index of association, IAS, was used as a measure of multilocus linkage disequilibrium (LD) [37]. The null hypothesis of linkage equilibrium, IAS = 0, was tested with 10,000 Monte Carlo simulations using the LIAN v. 3.7 program [37]. Allelic diversity at each locus (h) was also estimated using LIAN. The mean of the two groups (“intestinal” vs. “invasive”) of allelic frequencies was compared using a Mann–Whitney U test performed with Prism v. 6.0 (GraphPad Software Inc.). P-values that compared allele numbers from iNTS vs. non-iNTS isolates were computed by 2x2 contingency table with the two-tailed Fisher’s exact test in QuickCalcs (GraphPad Software Inc.). P-values for the genetic diversity comparisons between the 70 virulence factors and 40 virulence-associated proteins within each serovar, biotype or lineage studied were done with an unpaired t-test using Prism v. 6.0.

Population and phylogenetic analysis

Global optimal eBURST (goeBURST) algorithm implemented by PHYLOViZ Online (N locus variant = 25) was used for phylogenetic inference and data visualization of minimum spanning tree [72]. To evaluate population structure, we used Structure v. 2.3.4 [43]; four independent runs were performed for each value of the number of populations K ranging from 2 to 13. Each run consisted of 100,000 Markov Chain Monte Carlo iterations (the first 50,000 being discarded as burn-in iterations). We used the Evanno K method [41] representing the highest median likelihood values with the online CLUMPAK server [42]. The Comparing Partitions online tool was used for calculating Simpson’s diversity indexes [73] with confidence interval (CI)[74]. CINA is non-approximated CI from the original Simpson study [75]. Other metrics were calculated from sequence information for each VF (number of aminoacid for the mature VF, number of alleles divided by the total number of amino acids, number of variable amino acids, number of variable amino acids divided by the total number of amino acids).

Visualization softwares

All heatmap images were produced by using the complete linkage method with Euclidian distance metric for hierarchical clustering of the Morpheus server (Broad Institute, MIT, https://software.broadinstitute.org/morpheus ). Population structure results were visualized with the CLUMPAK server [42]. The Venn diagram was made by eulerAPE v. 3 [76].

Availability of data and materials

Not Applicable.

Abbreviations

(i)NTS:

(invasive) non-typhoidal Salmonella

IAS :

Index of association

LD:

Linkage disequilibrium

MLST:

Multilocus sequence typing

S. enterica :

Salmonella enterica

SNP:

Single nucleotide polymorphism

ST(s):

Sequence type(s)

T3SS:

Type 3 secretion system

VF(s):

Virulence factor(s)

References

  1. 1.

    Sakharkar KR, Chow VT. Strategies for genome reduction in microbial genomes. Genome Inform. 2005;16(2):69–75.

  2. 2.

    Graziani C, Losasso C, Luzzi I, Ricci A, Scavia G, Pasquali P: Chapter 5 - Salmonella. In: Foodborne Diseases (Third edition). Academic Press; 2017: 133-169.

  3. 3.

    Matthews TD, Schmieder R, Silva GG, Busch J, Cassman N, Dutilh BE, Green D, Matlock B, Heffernan B, Olsen GJ, et al. Genomic Comparison of the Closely-Related Salmonella enterica Serovars Enteritidis, Dublin and Gallinarum. PLoS One. 2015;10(6):e0126883.

  4. 4.

    Feng Y, Johnston RN, Liu GR, Liu SL. Genomic comparison between Salmonella Gallinarum and Pullorum: differential pseudogene formation under common host restriction. PLoS One. 2013;8(3):e59427.

  5. 5.

    Gordon MA. Invasive nontyphoidal Salmonella disease: epidemiology, pathogenesis and diagnosis. Curr Opin Infect Dis. 2011;24(5):484–9.

  6. 6.

    Okoro CK, Barquist L, Connor TR, Harris SR, Clare S, Stevens MP, Arends MJ, Hale C, Kane L, Pickard DJ, et al. Signatures of Adaptation in Human Invasive Salmonella Typhimurium ST313 Populations from Sub-Saharan Africa. PLoS Negl Trop Dis. 2015;9(3):e0003611.

  7. 7.

    Feasey NA, Hadfield J, Keddy KH, Dallman TJ, Jacobs J, Deng X, Wigley P, Barquist Barquist L, Langridge GC, Feltwell T, et al. Distinct Salmonella Enteritidis lineages associated with enterocolitis in high-income settings and invasive disease in low-income settings. Nat Genet. 2016;48(10):1211–7.

  8. 8.

    Kingsley RA, Bäumler AJ. Host adaptation and the emergence of infectious disease: the Salmonella paradigm. Mol Microbiol. 2000;36(5):1006–14.

  9. 9.

    Nuccio SP, Baumler AJ. Comparative analysis of Salmonella genomes identifies a metabolic network for escalating growth in the inflamed gut. MBio. 2014;5(2):e00929–14.

  10. 10.

    Feng Y, Chien KY, Chen HL, Chiu CH. Pseudogene recoding revealed from proteomic analysis of Salmonella serovars. J Proteome Res. 2012;11(3):1715–9.

  11. 11.

    Jalajakumari MB, Thomas CJ, Halter R, Manning PA. Genes for biosynthesis and assembly of CS3 pili of CFA/II enterotoxigenic Escherichia coli: novel regulation of pilus production by bypassing an amber codon. Mol Microbiol. 1989;3(12):1685–95.

  12. 12.

    Keestra-Gounder AM, Tsolis RM, Baumler AJ. Now you see me, now you don't: the interaction of Salmonella with innate immune receptors. Nat Rev Microbiol. 2015;13(4):206–16.

  13. 13.

    Deng L, Song J, Gao X, Wang J, Yu H, Chen X, Varki N, Naito-Matsui Y, Galan JE, Varki A. Host adaptation of a bacterial toxin from the human pathogen Salmonella Typhi. Cell. 2014;159(6):1290–9.

  14. 14.

    Yue M, Schifferli DM. Allelic variation in Salmonella: an underappreciated driver of adaptation and virulence. Front Microbiol. 2014;4:419.

  15. 15.

    Bäumler AJ, Tsolis RM, Ficht TA, Adams LG. Evolution of host adaptation in Salmonella enterica. Infect Immun. 1998;66(10):4579–87.

  16. 16.

    Selander RK, Beltran P, Smith NH, Helmuth R, Rubin FA, Kopecko DJ, Ferris K, Tall BD, Cravioto A, Musser JM. Evolutionary genetic relationships of clones of Salmonella serovars that cause human typhoid and other enteric fevers. Infect Immun. 1990;58(7):2262–75.

  17. 17.

    Didelot X, Bowden R, Street T, Golubchik T, Spencer C, McVean G, Sangal V, Anjum MF, Achtman M, Falush D, et al. Recombination and population structure in Salmonella enterica. PLoS Genet. 2011;7(7):e1002191.

  18. 18.

    Leekitcharoenphon P, Lukjancenko O, Friis C, Aarestrup FM, Ussery DW. Genomic variation in Salmonella enterica core genes for epidemiological typing. BMC Genomics. 2012;13:88.

  19. 19.

    Yoshida CE, Kruczkiewicz P, Laing CR, Lingohr EJ, Gannon VP, Nash JH, Taboada EN. The Salmonella In Silico Typing Resource (SISTR): An Open Web-Accessible Tool for Rapidly Typing and Subtyping Draft Salmonella Genome Assemblies. PLoS One. 2016;11(1):e0147101.

  20. 20.

    Liu YY, Chen CC, Chiou CS. Construction of a Pan-Genome Allele Database of Salmonella enterica Serovar Enteritidis for Molecular Subtyping and Disease Cluster Identification. Front Microbiol. 2016;7:2010.

  21. 21.

    Yue M, Han X, De Masi L, Zhu C, Ma X, Zhang J, Wu R, Schmieder R, Kaushik RS, Fraser GP, et al. Allelic variation contributes to bacterial host specificity. Nat Commun. 2015;6:8754.

  22. 22.

    De Masi L, Yue M, Hu C, Rakov AV, Rankin SC, Schifferli DM. Cooperation of Adhesin Alleles in Salmonella-Host Tropism. mSphere. 2017;2(2).

  23. 23.

    Carden SE, Walker GT, Honeycutt J, Lugo K, Pham T, Jacobson A, Bouley D, Idoyaga J, Tsolis RM, Monack D. Pseudogenization of the Secreted Effector Gene sseI Confers Rapid Systemic Dissemination of S. Typhimurium ST313 within Migratory Dendritic Cells. Cell Host Microbe. 2017;21(2):182–94.

  24. 24.

    Trombert AN, Rodas PI, Mora GC. Reduced invasion to human epithelial cell lines of Salmonella enterica serovar Typhi carrying S. Typhimurium sopD2. FEMS Microbiol Lett. 2011;322(2):150–6.

  25. 25.

    Urrutia IM, Fuentes JA, Valenzuela LM, Ortega AP, Hidalgo AA, Mora GC. Salmonella Typhi shdA: pseudogene or allelic variant? Infect Genet Evol. 2014;26:146–52.

  26. 26.

    Valenzuela LM, Hidalgo AA, RodrÌguez L, Urrutia IM, Ortega AP, Villagra NA, Paredes-Sabja D, CalderÛn IL, Gil F, Saavedra CP, et al. Pseudogenization of sopA and sopE2 is functionally linked and contributes to virulence of Salmonella enterica serovar Typhi. Infect Genet Evol. 2015;33:131–42.

  27. 27.

    Yue M, Rankin SC, Blanchet RT, Nulton JD, Edwards RA, Schifferli DM. Diversification of the Salmonella fimbriae: a model of macro- and microevolution. PLoS One. 2012;7(6):e38596.

  28. 28.

    Chen L, Zheng D, Liu B, Yang J, Jin Q. VFDB 2016: hierarchical and refined dataset for big data analysis--10 years on. Nucleic Acids Res. 2016;44(D1):D694–7.

  29. 29.

    Heffron F, Niemann G, Yoon H, Kidwai A, Brown RNE, McDermott JD, Smith R, Adlkins JN: Salmonella-secreted virulence factors. In: Salmonella : from genome to function. Edited by Porwollik S. Norfolk, UK: Caister Academic Press; 2011: 187-223.

  30. 30.

    Chen HD, Groisman EA. The biology of the PmrA/PmrB two-component system: the major regulator of lipopolysaccharide modifications. Annu Rev Microbiol. 2013;67:83–112.

  31. 31.

    Dalebroux ZD, Miller SI. Salmonellae PhoPQ regulation of the outer membrane to resist innate immunity. Curr Opin Microbiol. 2014;17:106–13.

  32. 32.

    McClelland M, Sanderson KE, Clifton SW, Latreille P, Porwollik S, Sabo A, Meyer R, Bieri T, Ozersky P, McLellan M, et al. Comparison of genome degradation in Paratyphi A and Typhi, human-restricted serovars of Salmonella enterica that cause typhoid. Nat Genet. 2004;36(12):1268–74.

  33. 33.

    Kingsley RA, Msefula CL, Thomson NR, Kariuki S, Holt KE, Gordon MA, Harris D, Clarke L, Whitehead S, Sangal V, et al. Epidemic multiple drug resistant Salmonella Typhimurium causing invasive disease in sub-Saharan Africa have a distinct genotype. Genome Res. 2009;19(12):2279–87.

  34. 34.

    Lawley TD, Chan K, Thompson LJ, Kim CC, Govoni GR, Monack DM. Genome-wide screen for Salmonella genes required for long-term systemic infection of the mouse. PLoS Pathog. 2006;2(2):e11.

  35. 35.

    McLaughlin LM, Govoni GR, Gerke C, Gopinath S, Peng K, Laidlaw G, Chien YH, Jeong HW, Li Z, Brown MD, et al. The Salmonella SPI2 effector SseI mediates long-term systemic infection by modulating host cell migration. PLoS Pathog. 2009;5(11):e1000671.

  36. 36.

    Worley MJ, Nieman GS, Geddes K, Heffron F. Salmonella typhimurium disseminates within its host by manipulating the motility of infected cells. Proc Natl Acad Sci U S A. 2006;103(47):17915–20.

  37. 37.

    Haubold B, Hudson RR. LIAN 3.0: detecting linkage disequilibrium in multilocus data. Linkage Analysis. Bioinformatics. 2000;16(9):847–8.

  38. 38.

    Sangal V, Harbottle H, Mazzoni CJ, Helmuth R, Guerra B, Didelot X, Paglietti B, Rabsch W, Brisse S, Weill FX, et al. Evolution and population structure of Salmonella enterica serovar Newport. J Bacteriol. 2010;192(24):6465–76.

  39. 39.

    Fabre L, Zhang J, Guigon G, Le Hello S, Guibert V, Accou-Demartin M, de Romans S, Lim C, Roux C, Passet V, et al. CRISPR typing and subtyping for improved laboratory surveillance of Salmonella infections. PLoS One. 2012;7(5):e36995.

  40. 40.

    Fowl Typhoid and Pullorum Disease [http://www.cfsph.iastate.edu/Factsheets/pdfs/fowl_typhoid.pdf]

  41. 41.

    Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14(8):2611–20.

  42. 42.

    Kopelman NM, Mayzel J, Jakobsson M, Rosenberg NA, Mayrose I. Clumpak: a program for identifying clustering modes and packaging population structure inferences across K. Mol Ecol Resour. 2015;15(5):1179–91.

  43. 43.

    Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155(2):945–59.

  44. 44.

    Barlag B, Hensel M. The giant adhesin SiiE of Salmonella enterica. Molecules. 2015;20(1):1134–50.

  45. 45.

    Holt KE, Parkhill J, Mazzoni CJ, Roumagnac P, Weill FX, Goodhead I, Rance R, Baker S, Maskell DJ, Wain J, et al. High-throughput sequencing provides insights into genome variation and evolution in Salmonella Typhi. Nat Genet. 2008;40(8):987–93.

  46. 46.

    Uche IV, MacLennan CA, Saul A. A Systematic Review of the Incidence, Risk Factors and Case Fatality Rates of Invasive Nontyphoidal Salmonella (iNTS) Disease in Africa (1966 to 2014). PLoS Negl Trop Dis. 2017;11(1):e0005118.

  47. 47.

    Langridge GC, Fookes M, Connor TR, Feltwell T, Feasey N, Parsons BN, Seth-Smith HM, Barquist L, Stedman A, Humphrey T, et al. Patterns of genome evolution that have accompanied host adaptation in Salmonella. Proc Natl Acad Sci U S A. 2015;112(3):863–8.

  48. 48.

    Betancor L, Yim L, MartÌnez A, Fookes M, Sasias S, Schelotto F, Thomson N, Maskell D, Chabalgoity JA. Genomic Comparison of the Closely Related Salmonella enterica Serovars Enteritidis and Dublin. Open Microbiol J. 2012;6:5–13.

  49. 49.

    Zhou Z, Lundstrom I, Tran-Dien A, Duchene S, Alikhan NF, Sergeant MJ, Langridge G, Fotakis AK, Nair S, Stenoien HK, et al. Pan-genome Analysis of Ancient and Modern Salmonella enterica Demonstrates Genomic Stability of the Invasive Para C Lineage for Millennia. Curr Biol. 2018;28(15):2420–2428 e2410.

  50. 50.

    Cao G, Meng J, Strain E, Stones R, Pettengill J, Zhao S, McDermott P, Brown E, Allard M. Phylogenetics and Differentiation of Salmonella Newport Lineages by Whole Genome Sequencing. PLoS One. 2013;8(2):e55687.

  51. 51.

    Gal-Mor O, Boyle EC, Grassl GA. Same species, different diseases: how and why typhoidal and non-typhoidal Salmonella enterica serovars differ. Front Microbiol. 2014;5:391.

  52. 52.

    Baumler A, Fang FC. Host specificity of bacterial pathogens. Cold Spring Harb Perspect Med. 2013;3(12):a010041.

  53. 53.

    Shivaprasad HL, Methner U, Barrow PA: Salmonella infections in the domestic fowl. In: Salmonella in domestic animals. Edited by Barrow PA, Methner U, 2nd edn. Wallingford, UK: CABI; 2013: 162-192.

  54. 54.

    Rabsch W, Hargis BM, Tsolis RM, Kingsley RA, Hinz KH, Tschape H, Baumler AJ. Competitive exclusion of Salmonella Enteritidis by Salmonella Gallinarum in poultry. Emerg Infect Dis. 2000;6(5):443–8.

  55. 55.

    Cao G, Allard M, Strain E, Stones R, Zhao S, Brown E, Meng J. Genetic diversity of Salmonella pathogenicity islands SPI-5 and SPI-6 in Salmonella Newport. Foodborne Pathog Dis. 2014;11(10):798–807.

  56. 56.

    Alcaine SD, Soyer Y, Warnick LD, Su WL, Sukhnanand S, Richards J, Fortes ED, McDonough P, Root TP, Dumas NB, et al. Multilocus sequence typing supports the hypothesis that cow- and human-associated Salmonella isolates represent distinct and overlapping populations. Appl Environ Microbiol. 2006;72(12):7575–85.

  57. 57.

    Eswarappa SM, Janice J, Nagarajan AG, Balasundaram SV, Karnam G, Dixit NM, Chakravortty D. Differentially evolved genes of Salmonella pathogenicity islands: insights into the mechanism of host specificity in Salmonella. PLoS One. 2008;3(12):e3829.

  58. 58.

    Zhang C, Plastow G. Genomic Diversity in Pig (Sus scrofa) and its Comparison with Human and other Livestock. Curr Genomics. 2011;12(2):138–46.

  59. 59.

    Friend M, Franson JC: Salmonellosis. In: Field Manual of Wildlife Diseases General Field Procedures and Diseases of Birds. Edited by Ciganovich EA. Washington, DC: USGS; 1999: 99-109.

  60. 60.

    Pan H, Zhou X, Chai W, Paudyal N, Li S, Zhou X, Zhou K, Wu Q, Wu B, Li G, et al. Diversified sources for human infections by Salmonella enterica serovar Newport. Transbound Emerg Dis. 2018;66(2):1044–8.

  61. 61.

    Yue M, Schmieder R, Edwards RA, Rankin SC, Schifferli DM. Microfluidic PCR combined with pyrosequencing for identification of allelic variants with phenotypic associations among targeted Salmonella genes. Appl Environ Microbiol. 2012;78(20):7480–2.

  62. 62.

    Johnson R, Byrne A, Berger CN, Klemm E, Crepin VF, Dougan G, Frankel G. The Type III Secretion System Effector SptP of Salmonella enterica Serovar Typhi. J Bacteriol. 2017;199:4.

  63. 63.

    Lunelli M, Hurwitz R, Lambers J, Kolbe M. Crystal structure of PrgI-SipD: insight into a secretion competent state of the type three secretion system needle tip and its interaction with host ligands. PLoS Pathog. 2011;7(8):e1002163.

  64. 64.

    Federhen S. Type material in the NCBI Taxonomy Database. Nucleic Acids Res. 2015;43(Database issue):D1086–98.

  65. 65.

    Afgan E, Sloggett C, Goonasekera N, Makunin I, Benson D, Crowe M, Gladman S, Kowsar Y, Pheasant M, Horst R, et al. Genomics Virtual Laboratory: A Practical Bioinformatics Workbench for the Cloud. PLoS One. 2015;10(10):e0140829.

  66. 66.

    Larsen MV, Cosentino S, Rasmussen S, Friis C, Hasman H, Marvig RL, Jelsbak L, Sicheritz-PontÈn T, Ussery DW, Aarestrup FM, et al. Multilocus sequence typing of total-genome-sequenced bacteria. J Clin Microbiol. 2012;50(4):1355–61.

  67. 67.

    Alikhan NF, Zhou Z, Sergeant MJ, Achtman M. A genomic overview of the population structure of Salmonella. PLoS Genet. 2018;14(4):e1007261.

  68. 68.

    Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol Biol Evol. 2016;33(7):1870–4.

  69. 69.

    Nielsen H. Predicting Secretory Proteins with SignalP. Methods Mol Biol. 2017;1611:59–73.

  70. 70.

    Villesen P. FaBox: an online toolbox for fasta sequences. Mol Ecol Notes. 2007;7(6):965–8.

  71. 71.

    Ekmekci B, McAnany CE, Mura C. An Introduction to Programming for Bioscientists: A Python-Based Primer. PLoS Comput Biol. 2016;12(6):e1004867.

  72. 72.

    Nascimento M, Sousa A, Ramirez M, Francisco AP, CarriÁo JA, Vaz C. PHYLOViZ 2.0: providing scalable data integration and visualization for multiple phylogenetic inference methods. Bioinformatics. 2017;33(1):128–9.

  73. 73.

    Pinto FR, Melo-Cristino J, Ramirez M. A confidence interval for the wallace coefficient of concordance and its application to microbial typing methods. PLoS One. 2008;3(11):e3696.

  74. 74.

    Grundmann H, Hori S, Tanner G. Determining confidence intervals when measuring genetic diversity and the discriminatory abilities of typing methods for microorganisms. J Clin Microbiol. 2001;39(11):4190–2.

  75. 75.

    Simpson EH. Measurement of diversity. Nature. 1949;163:688.

  76. 76.

    Micallef L, Rodgers P. eulerAPE: drawing area-proportional 3-Venn diagrams using ellipses. PLoS One. 2014;9(7):e101717.

Download references

Acknowledgements

We thank all the contributors for publicly sharing their Salmonella genomic data through the National Center for Biotechnology Information, European Bioinformatics Institute and Wellcome Sanger Institute Salmonella genome projects. We would like to thank Min Yue for preliminary data and valuable advice for some software usage, and Leslie King for some text corrections.

The opinions expressed in this manuscript are solely the responsibility of the authors and do not necessarily represent the official views and policy of the National Institutes of Health.

Funding

This study was supported by funds from USDA National Institute of Food and Agriculture grant 2013-67015-21285, NIH/NIAID grant AI117135, and the PennVet Center for Host-Microbial Interactions to DMS. AVR was supported by NIH/NIAID grant AI117135. The funding bodies had no role in in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

AVR: designed and performed all the bioinformatics analyses, prepared figures and Tables and wrote parts of the manuscript; EM: assisted with the design of Hamming distance analyses; SLL: participated in the redaction of the manuscript; DMS: designed and coordinated the study and wrote most of the manuscript. All authors read and approved the final version of the manuscript.

Correspondence to Dieter M. Schifferli.

Ethics declarations

Ethics approval and consent to participate

Not Applicable.

Consent for publication

Not Applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. List of virulence factors and “control” proteins selected for this study. (XLSX 20 kb)

Additional file 2:

Table S2. List of the 500 Salmonella genomes (7 serovars) selected and metadata. (XLSX 58 kb)

Additional file 3:

Table S3. List of virulence factor genes absent in some Salmonella strains or serovars. Genes were either present (+) or absent (-) in all or some strains (+/-). (XLSX 14 kb)

Additional file 4:

Table S4. List of allele numbers for the 70 virulence factors in the 500 Salmonella strains. (XLSX 114 kb)

Additional file 5:

Figure S1. Allelic diversity for 70 studied virulence factors among 500 Salmonella assessed by Hamming distance. The heatmap shows VFs represented as columns and organized in six groups as labeled on the top of the heatmap. Distance corresponds to gradient colors varied from blue (main allele/no amino acid changes) to red (most amino acid changes). Missing VFs are shown in white. (PPT 264 kb)

Additional file 6:

Table S5. Genetic diversity for 70 virulence factors among all serovars. Simpson’s index with confidence interval (CI), non-approximated confidence interval (CINA), number of amino acids (aa) in the mature protein (without the signal peptide, SP), number of alleles to the number of total amino acids, number of variable amino acids, percent of variable amino acids per total amino acids. (XLSX 17 kb)

Additional file 7:

Table S6. Genetic diversity for each serovar, lineage or biotype. Mean genetic diversity, and index of association IAS, were determined by the LIAN v. 3.7 software. Number of monoallelic VFs, VFs with allelic diversity h < 0.1, total number of monoallelic VFs and VFs with allelic diversity h < 0.1, were calculated with corresponding 95% confidence intervals. Cells were highlighted with a color gradient (white = 0 to darkest blue = 0.8655). (XLSX 15 kb)

Additional file 8:

Figure S2. Distribution of dominant alleles versus other alleles for 70 studied virulence factors. VFs are represented as horizontal bars and organized in six groups as labeled on the left side of the graph. The numbers of dominant alleles are shown separately for the intestinal serovars in green (S. Typhimurium, S. Enteritidis and S. Newport) and the invasive serovars in red (S. Typhi, S. Dublin, S. Gallinarum, S. Choleraesuis); all the other alleles are grouped and shown in white. Bars that do not reach 500 (number of Salmonella studied) represent missing VFs. (PPT 100 kb)

Additional file 9:

Figure S3. Best stratified population based on 70 virulence factor sequences from 500 Salmonella. The best estimated K value was equal to 10, as calculated by the Evanno K plot that represents the highest median likelihood values for each K, using the CLUMPAK server. (PPT 64 kb)

Additional file 10:

Figure S4. Population stratification by 70 virulence factor sequences from 500 Salmonella. The range of numbers of assumed populations K was tested from 2 to 13 as shown on the left using the Structure 2.3.4 program. The 70 VF alleles of 500 Salmonella were grouped by serovars, and lineages (S. Newport) or biotypes (S. Gallinarum), as shown on the top, and represented by thin vertical lines. The coloring of each vertical line was proportional to the ancestry of each isolate for each K population. (PPT 202 kb)

Additional file 11:

Figure S5. Minimum spanning tree based on the alleles of the 70 virulence factors. The minimum spanning tree was built with the goeBURST algorithm, using an N locus variant level equal to 25 to link all nodes with distances equal or above this level. The tree created nine clusters represented as red circle for the generalists (S. Enteritidis, S. Typhimurium, S. Newport lineage II and S. Newport lineage III) and as black circles for the septicemic serovars (S. Dublin, S. Choleraesuis, S. Gallinarum biovar Gallinarum, S. Gallinarum biovar Pullorum, and S. Typhi). Circle sizes correspond to the number of strains for each cluster. (PPT 91 kb)

Additional file 12:

Figure S6. Reorganized heatmap of the 70 virulence factors and 40 control gene products. This heatmap was based on the function of the virulence factors and using the hierarchical clustering data of the 500 Salmonella strains as determined in Fig. 4. In addition to the VFs studied, 40 gene products for virulence-associated proteins were evaluated as a comparative group. The 7 serovars each with 52-75 strains are shown on the left of the heatmap. VFs on the top, arrayed in 7 functional groups, starting with FimH, BcfD and StfH in the adhesin group and ending with the comparative group. Clade numbers, biovars (S. Gallinarum) and lineages (S. Newport) are shown on the left. Colors indicate different alleles, with black being the dominant allele for each VF among all 500 Salmonella; blue-green, the 2nd most frequent allele, etc. up to the 20th (and more) most frequent allele in purple, as shown in Fig. 3 (1-20); missing VFs are shown in white. (PPT 1058 kb)

Additional file 13:

Figure S7. Proportional Venn diagram for STs and their FimH alleles in S. Choleraesuis. The diagram displays the number of strains with a specific ST and FimH allele. As an example, for ST145 strains, 23 have a predicted FimH8 allele (red), 16 have a predicted FimH9 allele (blue-purple) and 1 is predicted to lack FimH (yellow), whereas all 13 ST66 strains are predicted to have the FimH9 allele. (PPT 82 kb)

Additional file 14:

Table S7. Genetic diversity for 40 virulence factors-associated proteins among all serovars. Values as described in Additional file: Table S5. (XLSX 12 kb)

Additional file 15:

Figure S8. Mean genetic diversity (H) of 70 virulence factors (blue) and 40 virulence-associated proteins (red) for each serovar/lineage/biotype studied. The two groups H values were compared using an unpaired t-test. Statistically significant differences are marked by asterisks. (PPT 115 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Rakov, A.V., Mastriani, E., Liu, S. et al. Association of Salmonella virulence factor alleles with intestinal and invasive serovars. BMC Genomics 20, 429 (2019) doi:10.1186/s12864-019-5809-8

Download citation

Keywords

  • Salmonella enterica
  • virulence factors
  • allelic variation
  • host adaptation
  • FimH
  • STs
  • iNTS