Skip to main content

Two homologous Salmonella serogroup C1-specific genes are required for flagellar motility and cell invasion



Salmonella is a major bacterial pathogen associated with a large number of outbreaks of foodborne diseases. Many highly virulent serovars that cause human illness belong to Salmonella serogroup C1, and Salmonella ser. Choleraesuis is a prominent cause of invasive infections in Asia. Comparative genomic analysis in our previous study showed that two homologous genes, SC0368 and SC0595 in Salmonella ser. Choleraesuis were unique to serogroup C1. In this study, two single-deletion mutants (Δ0368 and Δ0595) and one double-deletion mutant (Δ0368Δ0595) were constructed based on the genome. All these mutants and the wild-type strain were subjected to RNA-Seq analysis to reveal functional relationships of the two serogroup C1-specific genes.


Data from RNA-Seq indicated that deletion of SC0368 resulted in defects in motility through repression of σ28 in flagellar regulation Class 3. Consistent with RNA-Seq data, results from transmission electron microcopy (TEM) showed that flagella were not present in 0368 and 03680595 mutants resulting in both swimming and swarming defects. Interestingly, the growth rates of two non-motile mutants 0368 and 03680595 were significantly greater than the wild-type, which may be associated with up-regulation of genes encoding cytochromes, enhancing bacterial proliferation. Moreover, the 0595 mutant was significantly more invasive in Caco-2 cells as shown by bacterial enumeration assays, and the expression of lipopolysaccharide (LPS) core synthesis-related genes (rfaB, rfaI, rfaQ, rfaY, rfaK, rfaZ) was down-regulated only in the 03680595 mutant. In addition, this study also speculated that these two genes might be contributing to serotype conversion for Salmonella C1 serogroup based on their apparent roles in biosynthesis of LPS and the flagella.


A combination of biological and transcriptomic (RNA-Seq) analyses has shown that the SC0368 and SC0595 genes are involved in biosynthesis of flagella and complete LPS, as well as in bacterial growth and virulence. Such information will aid to revealing the role of these specific genes in bacterial physiology and evolution within the serogroup C1.

Peer Review reports


Salmonella infections in humans and animals are a global public health problem [1]. The genus Salmonella is comprised of two species, Salmonella bongori and Salmonella enterica. Of the six S. enterica subspecies, subspecies I accounts for more than 99.5% of isolated Salmonella strains. The serogroup classification after the subspecies level relies on differences in the surface O-antigens, and the individual serovar is distinguished by flagellar H-antigens and biochemical tests [2]. There are > 2600 Salmonella serotypes, and some of the more virulent serovars are within Salmonella serogroup C1 [3]. All serovars in this serogroup share the same O7 antigen. For example, Salmonella ser. Choleraesuis (O6,7:Hc,1,5) was originally isolated from pig intestines and causes extra-intestinal or focal infections in humans. This serovar is associated with a higher mortality rate compared to other Salmonella serovars [4]. Salmonella ser. Paratyphi C (O6,7[Vi]:Hc,1,5) is a human-restricted pathogen that causes paratyphoid fever, a serious and potentially fatal systemic infection [5]. Salmonella ser. Infantis (O6,7,14:Hr,1,5) typically causes gastroenteritis and is the fourth most prevalent serovar causing human infections in Europe [6].

In our previous study, seven conserved and specific genes in Salmonella serogroup C1 were identified by comparative genomic analysis [7]. These C1-specific genes mainly encoded membrane proteins with high numbers of transmembrane segments (TMS). Five of the C1-specific genes are located in the C1 rfb gene (O-antigen gene) cluster that includes a putative flippase (SC2092) and polymerase (SC2098) for the O7-antigen [8, 9]. The other two genes (SC0368 and SC0595) that located outside of the O-antigen gene cluster are highly homologous (83%) to each other, which encode hypothetical proteins with 10 TMS [7]. Currently, no related research on these two homologous genes is available. Notably, genes downstream of SC0368 are gtrB (encoding a glycosyl transferase) and gtrA (encoding a bactoprenol-linked glucose translocase) [10, 11]. Downstream genes of SC0595 are yfdH and rfbI, which are putative orthologs of gtrA and gtrB, respectively [12, 13]. These gtrABC genes are located outside the main O-antigen chromosomal gene clusters and are associated with a bacteriophage involved in serotype conversion via O-antigen glucosylation [10, 14]. The yfdH and rfbI genes are also found upstream of SC0594 but are oriented in the opposite direction. These observations suggested that the two homologous C1-specific genes might have evolutionary relationships associated with serotype conversion.

One of the C1-specific genes, SC2092, is associated with biosynthesis of the O-antigen of lipopolysaccharide (LPS) [8]. In addition, an SC2092 null mutant displayed a NaCl-dependent deficiency in motility and produced fewer flagella than the wild-type strain [8]. Since SC0368 and SC0595 are also C1-specific genes, we wished to determine if the two orthologous genes SC0368 and SC0595 affect LPS biosynthesis and motility in Salmonella. To test this hypothesis, in-frame single (0368 and 0595) and double (05950368) deletion mutants were constructed using Salmonella ser. Choleraesuis ATCC 10708 as the parental strain. We used a combination of biological and transcriptomic (RNA-Seq) analyses to determine whether these genes are necessary for bacterial growth, motility, cell invasion and LPS biosynthesis. Information from this study will aid in evaluation of the effects of these genes on bacterial structure and physiology, and will assist in identifying their role in serotype conversion in serogroup C1.


Differentially expressed genes (DEGs) in three mutants

To reveal functional relationships of the two orthologous Salmonella serogroup C1-specific genes SC0368 and SC0595, we constructed single (Δ0368 and Δ0595) and double (Δ0368Δ0595) deletion mutants (Figure S1) and used these strains for RNA-Seq analysis. The sequencing yielded more than 28 million reads for 12 samples and 98.83–99.58% of the high-quality reads were mapped to the reference genome sequence consisting of 4500 Salmonella genes (Table S1). According to analysis with the DEGseq software, 75, 80, and 211 genes were considered differentially expressed in the 0368, 0595, and 03680595 mutants compared to the wild-type strain, respectively (> 2-fold changes, FDR < 0.05). There were 12 genes in strain 0368 and 35 in strain 03680595 that showed > 10-fold down-regulation. On the other hand, 2 genes in strain 0595 and 1 in strain 03680595 were > 5-fold up-regulated (Table 1).

Table 1 Statistics of differentially expressed genes in three mutants compared to wild type

To validate the RNA-Seq results, a total of 18 up- and down-regulated genes derived from the three different mutants were selected for RT-qPCR analysis (Fig. 1 and Table S2). As shown in Fig. 1A, RT-qPCR data correlated well with RNA-Seq results (R2 = 0.89) with a p-value < 0.001, indicating that the RNA-Seq data were very reliable.

Fig. 1

Differential gene expression of SC0368 and SC0595 Salmonella ser. Choleraesuis mutants. A Correlation of RT-qPCR and RNA-Seq results for 18 differentially expressed genes between the 03680595, 0595, 0368 mutants and the parental wild-type strain ATCC 10708 (R2 = 0.89, p < 0.001). Distributions of down-regulated genes (B) and up-regulated genes (C) with at least a 2-fold change in the mutants relative to the wild-type strain. D Down-regulated KEGG pathways. E Up-regulated KEGG pathways. Blue, mutant 0368; Yellow, mutant 0595 and Green, mutant 03680595. F Repression of DEGs encoding for LPS core synthesis (marked with red arrow in D) in strain 03680595

A Venn diagram was used to visualize the distribution of DEGs in each deletion mutant. The fliA gene encoding a flagellar biosynthesis sigma factor (σ28) was down-regulated in all three mutants (Fig. 1B and Table S3). There were 4 genes that were up-regulated in all three mutants, including SCTRNA84 and SCTRNA85 (both encoding tRNAs), SC3193 (encoding a hypothetical protein), and glgS (encoding a glycogen synthesis protein) (Fig. 1C and Table S4).

Down-regulated genes in the 03680595 mutant (Table S3) were classified into 10 functional KEGG categories (Fig. 1D, p < 0.05). The only genes shared between strains 0368 and 03680595 fell within 5 functional classes, including chemotaxis, flagellar assembly, microbial metabolism, biosynthesis of secondary metabolites, and carbon metabolism (Fig. 1D, Table S3). There were 4 enriched pathways (propanoate metabolism, butanoate metabolism, fatty acid degradation, and fatty acid metabolism) in up-regulated genes only found in the 0595 mutant (Fig. 1E, red arrow). In addition, 6 genes involved in synthesis of the LPS core component were down-regulated only in the double deletion mutant Δ0368Δ0595 (Fig. 1D red arrow, and Table S3) and RT-qPCR assays confirmed that genes rfaB, rfaI, rfaQ, rfaY, rfaK and rfaZ were > 2-fold down-regulated (from 2.6 to 6.5) only in the 03680595 mutant (Fig. 1F, Table S2).

Deletion of SC0368 resulted in defects in flagellar motility

The wild-type strain and 0595 mutant both showed a high degree of motility by both swimming and swarming assays (Fig. 2A). However, this was not the case for 03680595 and 0368 mutants, which had not spread far from the inoculation point (Fig. 2A), indicating that an in-frame deletion of the SC0368 gene significantly reduced motility relative to the wild-type strain. Total proteins of the wild-type and the three mutant strains were separated by sodium dodecylsulphate polyacrylamide gel electrophoresis (SDS-PAGE) (Fig. 2B). A band with size at ~ 50–52 kDa, likely representing a flagellin protein [15], was absent in the 03680595 and 0368 mutants (Fig. 2B). Furthermore, TEM analysis indicated that the non-motile mutants 03680595 and 0368 lacked flagella compared with the wild-type and 0595 strains (Fig. 2C). Moreover, these two flagella defect mutants appeared to have outer membrane blebs that were not seen in wild-type or 0595 mutant (Fig. 2C, marked with black rows), indicating that the outer membrane is less stable in the absence of SC0368.

Fig. 2

Bacterial motility and flagellar biosynthesis of deletion mutants compared to wild-type (WT). A The photos of one out of six independent experiments was randomly selected to show the motility behaviors of the indicated strains (B) SDS-PAGE analysis of bacterial total proteins (C) TEM photomicrographs of flagella and bacterial morphological features. Some spherical structures seem like outer membrane blebs were marked with black rows in these two flagella defect mutants

Bacterial growth rates were increased in SC0368 deletion mutants

The growth rates of the wild-type strain and three mutants were investigated to determine whether the differences in growth were caused by deletion of the SC0368 and SC0595 genes. Interestingly, the two non-motile mutants 0368 and 03680595 grew significantly faster and to a higher yield than the motile strains by cell counts (p < 0.05, n = 6) (Fig. 3A). The average size of bacterial colonies for 0368 or 03680595 was larger than that of 0595 strain at 12 h (Fig. 3B). Specifically, 03680595 strains formed the largest colonies at 3.78 ± 0.04 mm, while the smallest were for 0595 strains at 2.26 ± 0.03 mm.

Fig. 3

The growth of the wild-type strain and mutants in culture. A Growth curves, error bars indicate that each experiment was repeated six times and B colony cultures of three deletion mutants and the wild-type strain (WT) at 12 h. The average size of the bacterial colonies for each strains was marked in panel B

SC0595 and SC0368 are related with invasion capacity in Caco-2 cells

Both plate counting and flow cytometry were used to test the invasion capacity in Caco-2 cells of the three mutants and the wild-type strain. The results from plate counting indicated that the 0595 mutant exhibited the greatest invasion capacity with ~ 15.0% recovery rate, and the invasion capacity of the 0595 mutant was significantly higher than that of the wild-type strain (p < 0.01, Fig. 4A). On the other hand, a significantly lower invasion capacity with < 2.5% recovery rate was observed of the SC0368 deletion mutants (~ 1.8% for 0368 and ~ 2.4% for 03680595, respectively) compared with the wild-type strain (p < 0.05, Fig. 4A). Similar results were confirmed by bacterial cell enumeration using flow cytometry (Fig. 4B). The fluorescence value of 103 was chosen as the cut-off for dye-bacteria binding because 100% unstained bacteria were covered at this point in the negative control group. In 100 μL of cell lysates, the percentage of bacterial cells strained by SYTO 9 was 41.0, 22.6, 51.9, and 25.8% for wild-type, 0368, 0595, and 03680595, respectively, which was considered as the cell invasive capacity of each stain. In other words, the cell invasion capacity increased with the deletion of SC0595 (0595), weakened with deletion of SC0368 (0368 and 03680595).

Fig. 4

Caco2 cell invasion analysis by plate counting and flow cytometry. A Experiments were performed in duplicate and repeated three times. Statistical analyses of the invasive capacity compared to wild-type (WT) were performed using the SPSS 13.0 for Windows software; a p < 0.05 was considered to be significant, marked with * and p < 0.01 was considered to be extremely significant, marked with **; B Negative control is the same volume of bacterial suspension without simultaneous staining treatment by SYTO 9, and the value of 103 was chosen as the cut-off for cell invasion because 100% unstained bacteria were covered. Red cycle: the percentage of live invasion bacteria cells labelled by SYTO 9 fluorescence at above 103

Recovery of motility, growth, and cell invasion with complementation of the three mutants

To demonstrate if the SC0368 or SC0595 mutation alone was responsible for the observed phenotypes, three complemented strains (0368-C, 0595-C, and 03680595-C) were generated. Compared with the wild-type strain, the motility defect and increased growth did not occur in the 0368-C and 03680595-C strains (Fig. 5A and B). Meanwhile, the invasion capacity of these three complemented strains (0368-C, 0595-C, and 03680595-C) were at a similar lever with the wild-type strain at 4.6 ± 0.4% ~ 5.1 ± 0.2% recovery rate (Fig. 5C).

Fig. 5

Data on flagellar motility, bacterial proliferation and cell invasion of complementary strains relative to wild-type strain


During the past 10 years, comparative genomics has been used by our team to mine conserved and specific markers for Salmonella enterica [16] at the species [17, 18], serogroup [7] and serotype [19] levels. Based on these markers, rapid identification methods including multiplex PCR [7, 19], real-time PCR [17], and real-time reverse-transcriptase PCR [20] methods have been consecutively developed. Since all these specific markers have been verified with a large number of Salmonella strains, these PCR detection systems showed high specificity, which has been recognized by experts in other research teams [21, 22].

Recently, we investigated the specific biological functions and roles in evolution of these verified detection markers by using bioinformatics and molecular biological methods [8, 18]. GO (Gene Ontology) clusters analysis has shown that species-specific detection markers are significantly enriched in pathogenesis and the type III secretion system [18]. Most of the specific markers for Salmonella ser. Enteritidis, an important serotype frequently found in contaminated eggs, are predicted phage-related proteins [19], some of which have been proven to be necessary in bacterial survival in egg white [23]. The majority of serogroup-specific markers are located in the rfb gene cluster that is associated with LPS biosynthesis by sugar metabolism or glycosyl and O-actetyl transfer [7]. Furthermore, we have discovered that the C1-specific gene, SC2092, was involved in flagellar motility [8]. In the current study, the two homologous C1-specific genes (SC0368 and SC0595) that are located outside the rfb gene cluster, were observed to have important roles in flagella motility, bacterial growth, and cell invasion.

The motility defect was due to the reduced expression of flagellar class 3 genes regulated by σ28

Previous researches have shown that more than 60 genes are involved in flagellar regulation in Salmonella, and they are divided into three classes: class 1 (early), class 2 (middle), and class 3 (late) according to their temporal expression after induction of the flagellar regulon [24, 25]. The class 3 promoters transcribe the late genes for filament assembly, flagella rotation, and chemotaxis [26]. Genes whose expression is required late in flagellar assembly are primarily transcribed by FliA [24]. In the current study, a total of 20 motility-associated genes that fell within 2 functional classes (chemotaxis and flagellar assembly) were significantly repressed in the non-motile strains (Fig. 1D, Table S3). Significantly, this group of 20 down-regulated genes in the current study (Fig. 6A, Tables S2 and S3) that included fliA encoding sigma factor σ28 and flgM encoding its cognate anti-σ-factor are required for class 3 transcription [26, 27]. Therefore, both the transcriptomic (Fig. 6A, Tables S2 and S3) and biological data (Fig. 2) demonstrated that deletion of SC0368 resulted in defects in motility through repression of σ28 in flagellar regulation Class 3 (Fig. 6C). On the other hand, fliA gene had a less amount of reduction (~ 3 fold in both RNAseq and RT-qPCR assays) in 0595 mutant and no other flagella related genes had differential expression that there was no difference with wild-type in flagella and motility. It is hard to explain from the current data how these two membrane proteins with 10 TMS encoded by SC0368 and SC0595 could affect the expression of σ28. There is a possibility that these signaling linkages might be through a regulatory system sensing envelope stress just like Cpx and σE [28], and the differences of the structure and the expression level between SC0368 and SC0595 would diversify stress-sensing response. Notably, some outer membrane blebs might appear in the SC0368 deletion mutants (Fig. 2C), suggesting that the homeostasis of the bacterial cell membrane has been damaged or reformed by some unknow stress-sensing mechanisms.

Fig. 6

Heat map of differentially expressed genes and schematic diagram of SC0368 in bacterial motility regulation. A Down-regulated genes. B Up-regulated genes. Each column represents the different mutants as indicated. Each row represents individual gene. Scale is log2fold changes ranging from − 9 (dark blue) to − 1 (cyan) and from + 1 (red) to + 3 (dark red) relative to wild-type, respectively. * Deleted gene. # DEGs shared in three mutants. C The regulation of SC0368 in bacterial motility

The RNA-Seq data give clues to explain the increase in growth with deletion of SC0368

The RNA-Seq data showed that all 4 genes in the heme ABC exporter operon (ccmA, ccmB, ccmC, and ccmD), encoding cytochrome C [29] and 1 gene (nrfA) encoding cytochrome c552 nitrite reductase [30] were up-regulated in both the 0368 and 03680595 mutants (Fig. 6B and Table S4). These genes participate in a variety of biological processes, including oxygen transport, oxygen binding, electron transfer, and transcriptional regulation [31, 32]. Enhanced proliferation resulting from the up-regulation of these genes has been reported for Campylobacter jejuni [33] and Salmonella ser. Typhimurium [34]. Furthermore, the rpsV gene that showed increased expression in both the 0368 and 03680595 mutants (Fig. 6B and Table S4) is a 30S ribosomal subunit that is a growth regulator during stationary phase [35]. Likewise, there is another possibility that the growth advantage observed in nonmotile mutants may be simply due to the lack of energy consumption to make flagella, because flagella biosynthesis is an energy intensive process. Additionally, a previous report has shown that the mutant of FliA in E. coil elevated fitness by increasing growth rates and reducing flagellar transcription, which is similar with these flagella defect mutants (0368 and 03680595) in this current study that the gene expression of fliA was significantly inhibited.

The double deletion of SC0368 and SC0595 genes may affect LPS core biosynthesis

Results from protein tertiary structure predictions ( based on the deduced amino acid sequence of SC0368 or SC0595 indicated that these putative proteins shared similarity with LPS biosynthesis related proteins. SC0368 shares 11.90% identity to the LPS core biosynthesis protein RfaG, and SC0595 shares 17.16% identity to the oligosaccharide transferase that is essential for side branch synthesis of the core oligosaccharide [36]. Furthermore, the known or putative functions of these 6 down-regulated genes (i.e. rfaB, rfaI, rfaQ, rfaY, rfaK and rfaZ) only in the double deletion mutant Δ0368Δ0595 are involved in LPS core biosynthesis (Fig. 1D red arrow, Fig. 6A and Table S3, marked with a). This suggests that SC0368 and SC0595 probably have synergistic effects on LPS core synthesis. In other words, the existence of these two highly homologous genes in the Salmonella ser. Choleraesuis genome may ensure the complete synthesis of LPS.

LPS in most Gram-negative bacteria consists of the lipid A moiety linked to the short core oligosaccharide and the distal O-antigen polysaccharide chain [37]. Several research groups have described diverse effects on motility in those LPS biosynthesis defect mutant strains such as in Proteus mirabilis [38], Salmonella [39] and E. coli [40]. Our previous results have shown that mutation of a Salmonella serogroup-C1-specific gene (SC2092) abrogated O-antigen biosynthesis and displayed motility deficiency due to lower flagella synthesis [8]. Results of the present study also showed that the two serogroup-C1-specific genes probably affected LPS core oligosaccharide synthesis (Fig. 1D red arrow, Fig. 1E and Table S3, marked with a) and flagellar motility. Therefore, these serogroup-specific genes may participate in maintaining bacterial cell integrity by controlling the synthesis of complete LPS and flagella.

Notably, these 6 down-regulated genes encoded proteins associated with synthesis of the side branch of the core oligosaccharide [9]. These included RfaY that adds phosphate groups to the second Hep residues and RfaB that adds a Gal residue to the first Glc of the second Hep residue, as well as RfaI that adds the third Glc residues to the second Glc [12, 41]. However, agglutination tests with O7 antibody in the 03680595 mutant were positive and no apparent changes in the LPS profiles was found in the wild-type and three mutants (data not shown); probably because the proteins encoded by the 6 down-regulated genes were not the binding sites of the O-antigen [42], or this repressed effect on the LPS core was not sufficient for detectable phenotypic changes. Therefore, more detailed studies are needed to prove if the quantity and integrity of synthesized LPS are affected in these deletion mutants.

The evolution significances of two homologous genes (SC0368 and SC0595) with different phenotypes

Previous studies have provided evidences that the coexistence of multiple homologous genes or proteins is of great significance in the evolution of bacteria [43,44,45]. Specifically, homologous recombination, mainly mediated by prophage, drives the evolution in some serotype-associated genes, suggesting a critical role of these evolutionary mechanisms in serotype diversification of Salmonella [46, 47]. Coincidentally in this current study, these two homologous genes are surrounded by some bacteriophage-related genes (gtrA and gtrB), which involved in serotype conversion via O-antigen glucosylation [10,11,12,13,14]. The deletion of SC0368 resulted in the loss of flagella and the double deletion of two homologous genes (SC0368 and SC0595) was involved in LPS core biosynthesis. We boldly speculate that these two serogroup-C1-specific genes may have potential relevance to serotype conversion because flagella and LPS are important factors of serotyping for Salmonella.

Otherwise, the differential regulation of orthologous genes or allelic differences in orthologous sequences may have ecological consequences, determining the range of niches a microorganism can occupy [48, 49], such as colonization of the intestine that related with host adaptability. In this current study, Salmonella ser. Choleraesuis was chosen as a representative serotype for C1-group mainly because only Choleraesuis contains both these two genes (SC0368 and SC0595). As a primitive transitional serotype, Salmonella ser. Choleraesuis have a narrow host range and occasionally infect humans and included a high percentage of pseudogenes [50]. The replacement and deletion of these redundant genes potentially allowing more effective interactions with and invasion of host cells to occur or enhancing virulence and motility, which provide a simple evolutionary pathway [50]. Recently, gene loss via deletion, insertional inactivation or truncation has been considered important in the evolution of highly pathogenic Salmonella [51, 52]. This may be the reason that the 0595 mutant, alike host-restricted Salmonella ser. Paratyphi C, only possessing gene SC0368, displayed highest cell invasion capacity (Fig. 4). Regrettably, more than half (52.6%) of these differentially expressed genes in the 0595 mutant encoded hypothetical proteins with unknown functions; moreover, the growth and motility of the 0595 mutant were similar to those of the wild-type strain, which make it difficult to find the exactly role of gene SC0595 played in increasing the cell invasion capacity. Previous studies have shown that flagella-based motility is an important component in the initiation of contact with epithelial cells for Salmonella [53], implying that the lower invasion capacity caused by deletion of SC0368 is due to the lack of flagella. In all, we speculated that these two homologous genes were non-redundant for LPS core biosynthesis and may play different additional roles in motility and virulence.


Both the transcriptomic and biological data showed that the absence of flagella due to the repressed expression of σ28 involved in flagella regulation Class 3 was responsible for the defects in motility of the SC0368 gene deletion mutants. The increased growth rate observed in the two SC0368 deletion mutants was correlated with up-regulated expression of genes in cytochrome biogenesis as identified by RNA-Seq data. Moreover, these two C1-specific genes were related to virulence in epithelial cells and potentially involved in serotype diversification. Only the double deletion of the two homologous serogroup C1-specific genes (03680595) resulted in the reduced expression of genes related to side branching in LPS core oligosaccharide biosynthesis. Information from this study can be used to evaluate effects of these specific genes (SC0368 and SC0595) in bacterial structure, physiology and serotype conversion of Salmonella serogroup C1.


Bacterial strains and growth conditions

Strains and plasmids used in this study are listed in Table S5. A representative strain of Salmonella serogroup C1, Salmonella ser. Choleraesuis ATCC 10708, was obtained from the Shanghai Entry-Exit Inspection and Quarantine Bureau of China. The SC0368 deletion mutant (0368), the SC0595 deletion mutant (0595), the double mutant (03680595) and the corresponding complemented strains (0368-C, 0595-C and 03680595-C) were constructed in this study (Table S5). All of the Salmonella and Escherichia coli strains were cultured in Luria-Bertani (LB) broth (Oxoid, Cambridge, United Kingdom). When necessary, chloramphenicol and ampicillin were added at concentrations of 35 mg mL− 1 and 100 mg mL− 1, respectively. The incubation temperature was 37 °C, and all broth cultures were aerated by agitation at 180 rpm.

Generation of Salmonella ser. Choleraesuis ATCC 10708 0368, 0595, 03680595 in-frame deletion mutants and complementation

Deletion mutants 0368, 0595 and 03680595 were constructed using a homologous recombination knockout procedure as described previously [54]. Primers were designed according to the principle of overlap extension PCR (Table S6). The fragments of homologous arms were amplified by two rounds of PCR and then were cloned into the pMD18-T vector (Takara, Dalian, China) to generate pMD180368 and pMD180595. DNA sequencing was carried out to confirm the correct deletion in the plasmids. Both pMD180368 and pMD180595 plasmids were digested with Xba I and Sac I, respectively, and then ligated into suicide vector pRE112. pRE112 vector carries a chloramphenicol resistance gene and a sucrose-sensitivity gene sacB. The marker-free in-frame deletion mutants 0368, 0595, and 03680595, were screened by double selection as described previously [8]. The resulting plasmids pRE0368 and pRE0595 were transformed into E. coli SM10λpir by electroporation (2500 V, 5 ms). E. coli SM10λpir cells (1 mL ~ 106 CFU mL− 1) containing plasmid pRE0368 or pRE0595 and 3 mL of Salmonella wild-type cells or 0368 cells (~ 106 CFU mL− 1) were mixed in a culture flask (quantities to 10 mL) for 8 h at 37 °C to accomplish the conjugation process. Recipient cells were plated onto LA (LB with 1.5% agar) supplemented with chloramphenicol (35 mg mL− 1) to select the trans-conjugant strains that contained the plasmid integrated into the Salmonella ser. Choleraesuis or Salmonella ser. Choleraesuis 0368 genome as a single crossover. After overnight growth of these single crossover strains induced by LB containing 8% (w/v) sucrose, colonies that lost the pRE112 vector were selected.

For the construction of complemented mutant strains, the constructed plasmid pRE0368-C was transferred into the 0368 mutant, and pRE0595-C was transferred into 0595 mutant, and pRE0368-C and pRE0595-C were successively transferred into the 03680595 mutant by electroporation. The resulting strains, Salmonella ser. Choleraesuis 0368, 0595, 03680595 and complemented strains were confirmed by PCR using primers listed in Table S6.

RNA isolation, RNA-sequencing, and analysis of differentially expressed genes

For all experiments, a single colony from the overnight cultures was inoculated into LB broth and was grown at 37 °C for 16 h with agitation at 180 rpm. The cultures were diluted 1:100 in fresh LB broth and incubated at 37 °C for 8 h (late exponential phase) with agitation at 180 rpm. Approximately 109 CFU mL−1of each strain was pelleted at 5000 g at 4 °C for 10 min. The cell pellets were used for total RNA extraction using a RNeasy Micro Kit (Qiagen, Hilden, Germany) according to the protocol supplied by the manufacturer. Total RNA was extracted from three replicates of each strain. RNA integrity was assessed using a 2100 Bio-analyzer (Agilent, Foster City, CA). Total RNA was further purified by the RNase-Free DNase Set (Qiagen) after quantification. Total RNA samples were stored at − 80 °C until used.

Total RNA samples were submitted to the Shanghai Biotechnology Corporation for mRNA enrichment and subsequent RNA-Seq experiments. Removal of 16S and 23S rRNAs from total RNA was performed using a Ribo-Zero Magnetic Kit (Gram-Negative Bacteria) (Epicentre, Illumina, Madison, WI). The mRNA was used to prepare individually bar-coded (indexed) RNA-Seq libraries with a TruSeq RNA Sample Prep Kit (Illumina, San Diego, CA). Library QC and quantitation were performed on all individual libraries using the Qubit assay (Thermo Fisher Scientific, Waltham, MA) and the Agilent 2100 Bioanalyzer. RNA-Seq libraries were sequenced on a HiSeq2500 platform at the Shanghai Biotechnology Corporation using Version 3 chemistry. Reads were base called and quality filtered with the CASAVA (Consensus Assessment of Sequence and Variance) version 1.8 pipeline (Illumina) to generate single-ended 50-bp reads. The reference genome sequence and functional annotation information of Salmonella ser. Choleraesuis were obtained from the NCBI database (NC_006905.1). Mapping was based on a minimal length of 50 bp with an allowance of up to two mismatches, and > 90% of the each read length had to map to the reference sequence for it to be considered a mapped read. After reads were mapped with the CLC Workbench, the total numbers of mapped reads for each gene were counted. These read counts were used for further identification of differentially expressed genes.

The gene expression levels were measured in fragments per kilobase of transcript per million mapped reads (FPKM). The differentially expressed genes (DEGs) were identified using the DEGseq package [55] with the false discovery rate < 0.05 and at least two-fold changes. For the pathway enrichment analysis, all of the identified DEGs were mapped to terms in the KEGG database (p < 0.05) (

RT-qPCR assays

To measure gene expression levels in different strains, RT-qPCR was carried out using gene-specific primers (Table S2). Total RNA (1.0 μg) was reverse transcribed to generate cDNA as the template for RT-qPCR following the manufacturer’s instructions (Takara, Dalian, China). The RT-qPCR conditions were as follows: 10 μL SYBR mix (Takara, Dalian, China), 1 μL each primer (10 μM), 1 μL cDNA, and 7 μL ddH2O. All data in the samples were normalized to the level of the 16S rRNA or recA (DNA recombination and repair protein, recombinase A) genes. Three independent technical replicates were carried out for each reaction.

Growth and motility assays

The wild-type strain and mutants were grown overnight at 37 °C in 5 mL LB broth. Then, 0.1 mL of the cell suspensions (~ 108 CFU mL− 1) were transferred into 30 mL flasks containing 10 mL of LB broth and grown at 37 °C with agitation at 180 rpm. Cultures were measured at an optical density of 600 nm (OD600) every 2–3 h over the course of 48 h using a spectrophotometer (Tecan, Mannedorf, Switzerland). Three replicate samples were included for each time point. Each experiment was repeated six times.

Bacterial swimming and swarming motility were assayed on LB agar plates containing different concentrations of agar as described previously [56]. Using sterile toothpicks, single colonies from streak plates were stabbed onto the swimming plates (containing 0.3% agar, used to test for motility) and incubated for 24 h. At least six independent colonies were examined for each strain. For swarming motility assays, 5 μL of an overnight culture grown in LB were spotted directly on each plate (containing 0.6% agar, used to test for motility) and allowed to dry for 10 min without the lids. Plates were covered and incubated at 37 °C for 24 h before observation.

SDS-PAGE analysis of bacterial total protein

Bacterial total proteins were isolated and analyzed by SDS-PAGE with minor modifications [57]. All of the mutant and wild-type strains were grown overnight in LB broth at 37 °C to reach the densities of 109 CFU mL− 1. To characterize the total cellular proteins, the bacterial culture (1 mL) was pelleted by centrifugation (12,000×g, 5 min, 4 °C), suspended in SDS-PAGE loading buffer (20% SDS, 25% glycerol, 0.5% β-mercaptoethanol, 0.06 M Tris–HCl, pH 6.8, 0.15% bromophenol blue) and heated at 100 °C for 10 min. All samples were then separated on Ready Gel Precast Tris-HCl polyacrylamide gels with 15 and 5% acrylamide in the separating and stacking gels, respectively (Bio-Rad, Hercules, CA), and then fixed overnight in buffer with 10% acetic acid and 40% methanol. Gels were stained with 2.5 g L− 1Coomassie brilliant blue for 2 h. All of the solutions were prepared fresh before use. The experiment was repeated three times.

Transmission electron microscopy

The flagella and morphological features of the deletion mutants and the wild-type strain were examined by transmission electron microscopy (TEM) as descry [58]. A suspension of each Salmonella strain was placed on copper grids, allowed to form a film on the grid for 2 min before the excess solution was removed using absorbent paper, and the grids were dried at room temperature. All samples were examined using a Tecnai G2 spirit Biotwin microscope (FEI, Japan) operated at an accelerating voltage of 120 kV.

Cell invasion and bacterial cell enumeration using flow cytometry

The human intestinal epithelial cell line Caco-2 (ATCC HTB-37) was cultured in Dulbecco’s Modified Eagle media (DMEM) containing 4 mM glutamine, glucose, 10% (v/v) fetal bovine serum (HyClone, Beijing, China) containing 100 U mL− 1penicillin and 100 μg mL−1 streptomycin (Life Technologies, Beijing, China) and incubated at 37 °C in the presence of a 5% CO2 in a humidified atmosphere. Cells were used between passages 8 and 10. Bacterial invasion was characterized using the gentamicin protection assay as described previously [59]. Data were represented as mean percent recovery of individual Salmonella strains in relation to the original inoculum. Experiments were performed in duplicate and repeated three times.

A total of 100 μL each cell lysate that contained bacteria that invaded was re-suspended in 900 μL saline to make a 1-mL bacterial suspension in a 1.5-mL-microcentrifuge tube. For staining, 1 μl of SYTO 9 (ThermoFisher Scientific) that labels live bacterial cells was aliquoted into the microcentrifuge tube. Each sample was incubated in the dark for 15 min at room temperature to allow dye-bacteria binding. Meanwhile, bacterial suspension from the same volume cell lysate without staining by SYTO 9 was as negative control. All samples were evaluated using a LSR II Flow Cytometer (BD biosciences, San Jose, CA, USA). SYTO 9 fluorescence was collected using a 505 nm longpass filter and bandpass filter with transmission at 530/30 nm. The cut-off of SYTO 9 fluorescence value for dye-bacteria binding was chosen based on negative control group at the point that 100% unstained bacteria were covered, and the percentage of bacterial cells above this cut-off was used to characterize the cell invasion ability of each strain.

Availability of data and materials

All data generated or analyzed during this study are included in this published article and its supplementary information files. All of the raw RNA-Seq data have been deposited in the Gene Expression Omnibus (GEO) with accession number GSE114577.



Transmission electron microcopy




Transmembrane segments


Sodium dodecylsulphate polyacrylamide gel electrophoresis




Fragments per kilobase of transcript per million mapped reads


Differentially expressed genes


Dulbecco’s Modified Eagle media


Gene Ontology


  1. 1.

    Haselbeck AH, Panzner U, Im J, Baker S, Meyer CG, Marks F. Current perspectives on invasive nontyphoidal Salmonella disease. Curr Opin Infect Dis. 2017;30(5):498–503.

    Article  PubMed  PubMed Central  Google Scholar 

  2. 2.

    Grimont PAD, Weill FX. Antigenic formulae of the Salmonella serovars. In: WHO Collaborating Centre for Reference and Research on Salmonella. 9th ed. Paris: Institut Pasteur; 2007.

    Google Scholar 

  3. 3.

    Fuche FJ, Sow O, Simon R, Tennant SM. Salmonella serogroup C: current status of vaccines and why they are needed. Clin Vaccine Immunol. 2016;23(9):737–45.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Chiu CH, Chuang CH, Chiu S, Su LH, Lin TY. Salmonella enterica serotype Choleraesuis infections in pediatric patients. Pediatrics. 2006;117:1193–6.

    Article  Google Scholar 

  5. 5.

    Zou QH, Li QH, Zhu HY, Feng Y, Li YG, Johnston RN, et al. SPC-P1: a pathogenicity-associated prophage of Salmonella paratyphi C. BMC Genomics. 2010;11(1):729.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    EFSA and ECDC (European Food Safety Authority and European Centre for Disease Prevention and Control). The European union summary report on trends and sources of zoonoses, zoonotic agents and food-borne outbreaks in 2013. EFSA J. 2015;13:3991.

    Article  Google Scholar 

  7. 7.

    Liu B, Zhang L, Zhu X, Shi C, Chen J, Liu W, et al. PCR identification of Salmonella serogroups based on specific targets obtained by comparative genomics. Int J Food Microbiol. 2011;144(3):511–8.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Zhou X, Liu B, Shi C, Shi X. Mutation of a Salmonella serogroup-C1-specific gene abrogates O7-antigen biosynthesis and triggers NaCl-dependent motility deficiency. PLoS One. 2014;9(9):e106708.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Yoshida C, Franklin K, Konczy P, McQuiston JR, Fields PI, Nash JH, et al. Methodologies towards the development of an oligonucleotide microarray for determination of Salmonella serotypes. J Microbiol Methods. 2007;70(2):261–71.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Korres H, Mavris M, Morona R, Manning PA, Verma NK. Topological analysis of GtrA and GtrB proteins encoded by the serotype-converting cassette of Shigella flexneri. Biochem Biophys Res Commun. 2005;328(4):1252–60.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Wang Y, Chen AM, Yu AY, Luo L, Yu GQ, Zhu JB, et al. The GntR-type regulators gtrA and gtrB affect cell growth and nodulation of Sinorhizobium meliloti. J Microbiol. 2008;46(2):137–45.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    McClelland M, Sanderson KE, Spieth J, Clifton SW, Latreille P, Courtney L, et al. Complete genome sequence of Salmonella enterica serovar Typhimurium LT2. Nature. 2001;413(6858):852–6.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Jiang XM, Neal B, Santiago F, Lee SJ, Romana LK, Reeves PR. Structure and sequence of the rfb (O antigen) gene cluster of Salmonella serovar typhimurium (strain LT2). Mol Microbiol. 1991;5(3):695–713.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Wang L, Andrianopoulos K, Liu D, Popoff MY, Reeves PR. Extensive variation in the O-antigen gene cluster within one Salmonella enterica serogroup reveals an unexpected complex history. J Bacteriol. 2002;184(6):1669–77.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Simon R, Curtis B, Deumic V, Nicki J, Tennant SM, Pasetti MF, et al. A scalable method for biochemical purification of Salmonella flagellin. Protein Expr Purif. 2014;102:1–7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Yu S, Liu W, Shi C, Wang D, Dan X, Li X, et al. SMM-system: a mining tool to identify specific markers in Salmonella enterica. J Microbiol Methods. 2011;84(3):423–9.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Chen J, Zhang L, Paoli GC, Shi C, Tu SI, Shi X. A real-time PCR method for the detection of Salmonella enterica from food using a target sequence identified by comparative genomic analysis. Int J Food Microbiol. 2010;137(2–3):168–74.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Zhou X, Zhang L, Shi C, Fratamico PM, Liu B, Paoli GC, et al. Genome-scale screening and validation of targets for identification of Salmonella enterica and serovar prediction. J Food Prot. 2016;79(3):376–83.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Liu B, Zhou XJ, Zhang LD, Liu WB, Dan XL, Shi CL, et al. Development of a novel multiplex PCR assay for the identification of Salmonella enterica typhimurium and Enteritidis. Food Control. 2012;27(1):87–93.

    CAS  Article  Google Scholar 

  20. 20.

    Zhou M, Yang JL, Zhou XJ, Liu B, Liu DX, Yuan CG, et al. Development of a sigDE-based real-time reverse-transcriptase PCR for the detection of viable Salmonella enterica. Foodborne Pathog Dis. 2014;11(7):537–44.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Navaneethan U, Giannella RA. Infectious colitis. Curr Opin Gastroenterol. 2011;27(1):66–71.

    Article  PubMed  Google Scholar 

  22. 22.

    Zhou Y, Wylie KM, El Feghaly RE, Mihindukulasuriya KA, Elward A, Haslam DB, et al. Metagenomic approach for identification of the pathogens associated with diarrhea in stool specimens. J Clin Microbiol. 2016;54(2):368–75.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Raspoet R, Appia-Ayme C, Shearer N, Martel A, Pasmans F, Haesebrouck F, et al. Microarray-based detection of Salmonella enterica serovar Enteritidis genes involved in chicken reproductive tract colonization. Appl Environ Microbiol. 2014;80(24):7710–6.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  24. 24.

    Fitzgerald DM, Bonocora RP, Wade JT. Comprehensive mapping of the Escherichia coli flagellar regulatory network. PLoS Genet. 2014;10(10):e1004649.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Chilcott GS, Hughes KT. Coupling of flagellar gene expression to flagellar assembly in Salmonella enterica serovar typhimurium and Escherichia coli. Microbiol Mol Biol Rev. 2000;64(4):694–708.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Barembruch C, Hengge R. Cellular levels and activity of the flagellar sigma factor FliA of Escherichia coli are controlled by FlgM-modulated proteolysis. Mol Microbiol. 2007;65(1):76–89.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Karlinsey JE, Tanaka S, Bettenworth V, Yamaguchi S, Boos W, Aizawa SI, et al. Completion of the hook-basal body complex of the Salmonella typhimurium flagellum is coupled to FlgM secretion and fliC transcription. Mol Microbiol. 2000;37(5):1220–31.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Grabowicz M, Silhavy TJ. Envelope stress responses: an interconnected safety net. Trends Biochem Sci. 2017;42(3):232–42.

    CAS  Article  PubMed  Google Scholar 

  29. 29.

    Sanders C, Turkarslan S, Lee DW, Daldal F. Cytochrome c biogenesis: the Ccm system. Trends Microbiol. 2010;18(6):266–74.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Einsle O. Structure and function of formate-dependent cytochrome c nitrite reductase, NrfA. Methods Enzymol. 2011;496:399–422.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Schulz H, Hennecke H, Thöny-Meyer L. Prototype of a heme chaperone essential for cytochrome c maturation. Science. 1998;281(5380):1197–200.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Clarke TA, Cole JA, Richardson DJ, Hemmings AM. The crystal structure of the pentahaem c-type cytochrome NrfB and characterization of its solution-state interaction with the pentahaem nitrite reductase NrfA. Biochem J. 2007;406(1):19–30.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Pittman MS, Elvers KT, Lee L, Jones MA, Poole RK, Park SF, et al. Growth of Campylobacter jejuni on nitrate and nitrite: electron transport to NapA and NrfA via NrfH and distinct roles for NrfA and the globin Cgb in protection against nitrosative stress. Mol Microbiol. 2007;63(2):575–90.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    Mills PC, Rowley G, Spiro S, Hinton JC, Richardson DJ. A combination of cytochrome c nitrite reductase (NrfA) and flavorubredoxin (NorV) protects Salmonella enterica serovar typhimurium against killing by NO in anoxic environments. Microbiology. 2008;154(Pt 4):1218–28.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    Izutsu K, Wada C, Komine Y, Sako T, Ueguchi C, Nakura S, et al. Escherichia coli ribosome-associated protein SRA, whose copy number increases during stationary phase. J Bacteriol. 2001;183(9):2765–73.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Musumeci MA, Faridmoayer A, Watanabe Y, Feldman MF. Evaluating the role of conserved amino acids in bacterial O-oligosaccharyl transferases by in vivo, in vitro and limited proteolysis assays. Glycobiology. 2014;24(1):39–50.

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    Wang X, Quinn PJ. Lipopolysaccharide: biosynthetic pathway and structure modification. Prog Lipid Res. 2010;49(2):97–107.

    CAS  Article  PubMed  Google Scholar 

  38. 38.

    Morgenstein RM, Clemmer KM, Rather PN. Loss of the waaL O-antigen ligase prevents surface activation of the flagellar gene cascade in Proteus mirabilis. J Bacteriol. 2010;192(12):3213–21.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Toguchi A, Siano M, Burkart M, Harshey RM. Genetics of swarming motility in Salmonella enterica serovar typhimurium: critical role for lipopolysaccharide. J Bacteriol. 2000;182(22):6308–21.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Girgis HS, Liu Y, Ryu WS, Tavazoie S. A comprehensive genetic characterization of bacterial motility. PLoS Genet. 2007;3(9):1644–60.

    CAS  Article  PubMed  Google Scholar 

  41. 41.

    Bociek K, Ferluga S, Mardirossian M, Benincasa M, Tossi A, Gennaro R, et al. Lipopolysaccharide phosphorylation by the WaaY kinase affects the susceptibility of Escherichia coli to the human antimicrobial peptide LL-37. J Biol Chem. 2015;290(32):19933–41.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  42. 42.

    Wang Z, Wang J, Ren G, Li Y, Wang X. Influence of core oligosaccharide of lipopolysaccharide to outer membrane behavior of Escherichia coli. Mar Drugs. 2015;13(6):3325–39.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  43. 43.

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

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Cheng RA, Eade CR, Wiedmann M. Embracing diversity: differences in virulence mechanisms, disease severity, and host adaptations contribute to the success of nontyphoidal Salmonella as a foodborne pathogen. Front Microbiol. 2019;10:1368.

    Article  PubMed  PubMed Central  Google Scholar 

  45. 45.

    Soyer Y, Orsi RH, Rodriguez-Rivera LD, Sun Q, Wiedmann M. Genome wide evolutionary analyses reveal serotype specific patterns of positive selection in selected Salmonella serotypes. BMC Evol Biol. 2009;9(1):264.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  46. 46.

    Zhang J, Wang H, Huang Q, Zhang Y, Zhao L, Liu F, et al. Four superoxide dismutases of Bacillus cereus 0-9 are non-redundant and perform different functions in diverse living conditions. World J Microbiol Biotechnol. 2020;36(1):12.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Liao J, Orsi RH, Carroll LM, Kovac J, Ou H, Zhang H, et al. Serotype-specific evolutionary patterns of antimicrobial-resistant Salmonella enterica. BMC Evol Biol. 2019;19(1):132.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  48. 48.

    Winfield MD, Groisman EA. Phenotypic differences between Salmonella and Escherichia coli resulting from the disparate regulation of homologous genes. Proc Natl Acad Sci U S A. 2004;101(49):17162–7.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  49. 49.

    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.

    CAS  Article  PubMed  Google Scholar 

  50. 50.

    Chiu CH, Tang P, Chu C, Hu S, Bao Q, Yu J, et al. The genome sequence of Salmonella enterica serovar Choleraesuis, a highly invasive and resistant zoonotic pathogen. Nucleic Acids Res. 2005;33(5):1690–8.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  51. 51.

    Holt KE, Thomson NR, Wain J, Langridge GC, Hasan R, Bhutta ZA, et al. Pseudogene accumulation in the evolutionary histories of Salmonella enterica serovars Paratyphi A and Typhi. BMC Genomics. 2009;10(1):36.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  52. 52.

    Kisiela DI, Chattopadhyay S, Libby SJ, Karlinsey JE, Fang FC, Tchesnokova V, et al. Evolution of Salmonella enterica virulence via point mutations in the fimbrial adhesin. PLoS Pathog. 2012;8(6):e1002733.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  53. 53.

    Ibarra JA, Knodler LA, Sturdevant DE, Virtaneva K, Carmody AB, Fischer ER, et al. Induction of Salmonella pathogenicity island 1 under different growth conditions can affect Salmonella-host cell interactions in vitro. Microbiology. 2010;156(Pt 4):1120–33.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  54. 54.

    Ho SN, Hunt HD, Horton RM, Pullen JK, Pease LR. Site-directed mutagenesis by overlap extension using the polymerase chain reaction. Gene. 1989;77(1):51–9.

    CAS  Article  PubMed  Google Scholar 

  55. 55.

    Wang L, Feng Z, Wang X, Wang X, Zhang X. DEGseq: An R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010;26(1):136–8.

    CAS  Article  PubMed  Google Scholar 

  56. 56.

    Legendre G, Faÿ F, Linossier I, Vallée-Réhel K. Evaluation of antibacterial activity against Salmonella Enteritidis. J Microbiol. 2011;49(3):349–54.

    CAS  Article  PubMed  Google Scholar 

  57. 57.

    Elhadad D, Desai P, Rahav G, McClelland M, Gal-Mor O. Flagellin is required for host cell invasion and normal Salmonella pathogenicity island 1 expression by Salmonella enterica serovar Paratyphi A. Infect Immun. 2015;83(9):3355–68.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  58. 58.

    Burdett ID, Murray RG. Septum formation in Escherichia coli: characterization of septal structure and the effects of antibiotics on cell division. J Bacteriol. 1974;119(1):303–24.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  59. 59.

    Gokulan K, Khare S, Rooney AW, Han J, Lynne AM, Foley SL. Impact of plasmids, including those encodingVirB4/D4 type IV secretion systems, on Salmonella enterica serovar Heidelberg virulence in macrophages and epithelial cells. PLoS One. 2013;8(10):e77866.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

Download references


We would like to thank Dr. James Smith and Dr. Xiaoyuan Wang for critical reading of the manuscript.


This work was supported by the National Key R&D Program of China (grant number 2017YFC1601200) and the National Natural Science Foundation of China (grant No. 32072320).

Author information




Z.X., S.C., W.D. and S. X. conceived and designed the research. Z.X., L.B., Z.L., Z.J. and C.Y. performed the experiments. Z.X. and L.Y. wrote the manuscript. F.P., X. P. and S. X. revised the manuscript. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Xianming Shi.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no conflicts of interest in connection with the content of this article.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Table S1

. Summary of RNA-Seq data.

Additional file 2: Table S2

. Primers used in RT-qPCR assays.

Additional file 3: Table S3

. Down-regulated differentially expressed genes in strain 03680595.

Additional file 4: Table S4

. Up-regulated differentially expressed genes in strain 03680595.

Additional file 5: Table S5

. Strains and plasmids used in this study.

Additional file 6: Table S6

. Primers used for mutant construction and complementation.

Additional file 7: Figure S1

. PCR verification of gene deletion mutants. M1:1 Kb ladder; M2: 200 bp ladder; Lanes 1, 5, 9, 13, PCR products using primers SC0368-for and SC0368-rev; Lanes 2, 6, 10, 14 PCR products using primers SC0595-for and SC0595-rev; Lanes 3, 7, 11, 15 PCR products using primers qSC0368-for and qSC0368-rev; Lanes 4, 8, 12, 16 PCR products using primers qSC0595-for and qSC0595-rev. Primer sequences are listed in Table S6.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhou, X., Liu, B., Liu, Y. et al. Two homologous Salmonella serogroup C1-specific genes are required for flagellar motility and cell invasion. BMC Genomics 22, 507 (2021).

Download citation


  • Salmonella
  • Serogroup C1
  • RNA-Seq
  • Motility
  • Virulence
  • Growth
  • Choleraesuis