# Intracellular Salmonella Paratyphi A is motile and differs in the expression of flagella-chemotaxis, SPI-1 and carbon utilization pathways in comparison to intracellular S. Typhimurium

Helit Cohen, Claire Hoede, Felix Scharte, Charles Coluzzi, Emiliano Cohen, Inna Shomer, Ludovic Mallet, Sébastien Holbert, Remy Felix Serre, Thomas Schiex, et al.

► To cite this version:

Helit Cohen, Claire Hoede, Felix Scharte, Charles Coluzzi, Emiliano Cohen, et al. Intracellular Salmonella Paratyphi A is motile and differs in the expression of flagella-chemotaxis, SPI-1 and carbon utilization pathways in comparison to intracellular S. Typhimurium. PLoS Pathogens, 2022, 18 (4), 35 p. 10.1371/journal.ppat.1010425. [https://doi.org/10.1371/journal.ppat.1010425](https://doi.org/10.1371/journal.ppat.1010425)

HAL Id: hal-03690615

https://hal.inrae.fr/hal-03690615

Submitted on 8 Jun 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

# Intracellular *Salmonella* Paratyphi A is motile and differs in the expression of flagella-chemotaxis, SPI-1 and carbon utilization pathways in comparison to intracellular *S. Typhimurium*

Helit Cohen<sup>1,°</sup>, Claire Hoede<sup>2,°</sup>, Felix Scharte<sup>3,°</sup>, Charles Coluzzi<sup>4</sup>, Emiliano Cohen<sup>1</sup>, Inna Shomer<sup>1</sup>, Ludovic Mallet<sup>2</sup>, Sébastien Holbert<sup>5</sup>, Remy Felix Serre<sup>6</sup>, Thomas Schiex<sup>7</sup>, Isabelle Virlogeux-Payant<sup>5</sup>, Guntram A. Grassl<sup>8</sup>, Michael Hensel<sup>3,9*</sup>, Hélène Chiapello<sup>2,4*</sup>, Ohad Gal-Mor<sup>1,10*</sup>

1 The Infectious Diseases Research Laboratory, Sheba Medical Center, Tel-Hashomer, Israel, 2 Université Fédérale de Toulouse, INRAE, BioinfOmics, UR MIAT, GenoToul Bioinformatics facility, 31326, Castanet-Tolosan, France, 3 Abt. Mikrobiologie, Universität Osnabrück, Osnabrück, Germany, 4 INRAE, Université Paris-Saclay, MalAGE, Jouy-en-Josas, France, 5 INRAE, Université de Tours, UMR ISP, Nouzilly, France, 6 INRAE, GeT-PlaGe, Genotoul, Castanet-Tolosan, France, 7 Université Fédérale de Toulouse, ANITI, INRAE, Toulouse, France, 8 Institute of Medical Microbiology and Hospital Epidemiology, Hannover Medical School and German Center for Infection Research (DZIF), Hanover, Germany, 9 CellNanOs-Center of Cellular Nanoanalistics Osnabrück, Universität Osnabrück, Germany, 10 Department of Clinical Microbiology and Immunology, Faculty of Medicine, Tel-Aviv University, Tel-Aviv, Israel

© These authors contributed equally to this work.

* michael.hensel@uni-osnabrueck.de (MH); helene.chiapello@inrae.fr (HC); Ohad.Gal-Mor@sheba.health.gov.il (OG-M)

<figure><img src="image_3.png" /></figure>

## OPEN ACCESS

**Citation:** Cohen H, Hoede C, Scharte F, Coluzzi C, Cohen E, Shomer I, et al. (2022) Intracellular *Salmonella* Paratyphi A is motile and differs in the expression of flagella-chemotaxis, SPI-1 and carbon utilization pathways in comparison to intracellular *S. Typhimurium*. PLoS Pathog 18(4): e1010425. https://doi.org/10.1371/journal.ppat.1010425

**Editor:** Denise M. Monack, Stanford University School of Medicine, UNITED STATES

**Received:** October 18, 2021

**Accepted:** March 9, 2022

**Published:** April 5, 2022

**Copyright:** © 2022 Cohen et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** The raw reads data for this study have been deposited in the European Nucleotide Archive (ENA) at EMBL-EBI under accession number PRJEB46495 and as specified in S7 Table.

**Funding:** This work was supported by the Infect-Era program, SalHostTrop project ("Understanding the Human-Restricted Host Tropism of Typhoidal Salmonella", 2016-2020) and by the France

## Abstract

Although *Salmonella* Typhimurium (STM) and *Salmonella* Paratyphi A (SPA) belong to the same phylogenetic species, share large portions of their genome and express many common virulence factors, they differ vastly in their host specificity, the immune response they elicit, and the clinical manifestations they cause. In this work, we compared their intracellular transcriptomic architecture and cellular phenotypes during human epithelial cell infection. While transcription induction of many metal transport systems, purines, biotin, PhoPQ and SPI-2 regulons was similar in both intracellular SPA and STM, we identified 234 differentially expressed genes that showed distinct expression patterns in intracellular SPA vs. STM. Surprisingly, clear expression differences were found in SPI-1, motility and chemotaxis, and carbon (mainly citrate, galactonate and ethanolamine) utilization pathways, indicating that these pathways are regulated differently during their intracellular phase. Concurring, on the cellular level, we show that while the majority of STM are non-motile and reside within *Salmonella*-Containing Vacuoles (SCV), a significant proportion of intracellular SPA cells are motile and compartmentalized in the cytosol. Moreover, we found that the elevated expression of SPI-1 and motility genes by intracellular SPA results in increased invasiveness of SPA, following exit from host cells. These findings demonstrate unexpected flagellum-dependent intracellular motility of a typhoidal *Salmonella* serovar and intriguing differences
Génomique National infrastructure, funded as part of "Investissement d'avenir" program managed by Agence Nationale pour la Recherche (ANR-10-INBS-09 contract). The work at the Gal-Mor laboratory was supported by grant numbers: 2616/18 from the joint ISF-Broad Institute program; 3-12435 from Infect-Era /Chief Scientist Ministry of Health; I-41-416.6-2018 from the German-Israeli Foundation for Scientific Research and Development (GIF, awarded to OGM and MH); and A128055 from the Research Cooperation Lower Saxony - Israel (The Volkswagen Foundation, awarded to OGM, MH and GG). The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.

**Competing interests:** The authors have declared that no competing interests exist.

in intracellular localization between typhoidal and non-typhoidal salmonellae. We propose that these differences facilitate new cycles of host cell infection by SPA and may contribute to the ability of SPA to disseminate beyond the intestinal lamina propria of the human host during enteric fever.

## Author summary

*Salmonella enterica* is a ubiquitous, facultative intracellular animal and human pathogen. Although non-typhoidal *Salmonella* (NTS) and typhoidal *Salmonella* serovars belong to the same phylogenetic species and share many virulence factors, the disease they cause in humans is very different. While the underlying mechanisms for these differences are not fully understood, one possible reason expected to contribute to their different pathogenicity is a distinct expression pattern of genes involved in host-pathogen interactions. Here, we compared the global gene expression and intracellular phenotypes, during human epithelial cell infection of S. Paratyphi A (SPA) and S. Typhimurium (STM), as prototypical serovars of typhoidal and NTS, respectively. Interestingly, we identified different expression patterns in key virulence and metabolic pathways, cytosolic motility and increased reinvasion of SPA, following exit from infected cells. We hypothesize that these differences contribute to the invasive and systemic disease developed following SPA infection in humans.

## Introduction

*Salmonella enterica* (*S. enterica*) is an abundant Gram-negative, facultative intracellular animal and human pathogen. This highly diverse bacterial species contains more than 2600 antigenically distinct serovars (biotypes) that are different in their host specificity and the disease they cause. The developed clinical manifestation is the result of multiple factors, but largely depends on the characteristics of the infecting serovar and the immunological status of the host [1]. Many of the *S. enterica* serovars, including the ubiquitous *S. enterica* serovar Typhimurium (*S. Typhimurium*) are generalist pathogens and are capable of infecting a broad range of host species. Infection of immunocompetent humans by non-typhoidal *Salmonellae* (NTS) typically leads to a self-limiting, acute inflammatory gastroenteritis, confined to the terminal ileum and colon. In contrast, a few serovars including Typhi, Paratyphi A and Sendai, collectively referred to as 'typhoidal salmonellae' are restricted to the human host and known as the causative agents of enteric (typhoid) fever. In most cases, enteric fever is a non-inflammatory, systemic life-threatening disease, presented as bacteremia and dissemination of the pathogen to systemic sites such as the spleen, liver and lymph nodes [2-4].

*S. enterica* infections are still considered a significant cause of mortality and morbidity with an annual incidence of over 27 million cases of enteric fever [5], and 78.7 million cases of gastroenteritis [6] worldwide. In recent years, the global prevalence of *S. Paratyphi A* (SPA) is increasing and in some countries (especially in eastern and southern Asia), SPA infections are accountable for up to 50% of all enteric fever cases [7,8]. The lack of a commercial SPA vaccine and the increased occurrence of antibiotic resistant strains illuminate SPA as a significant public health concern that is still an understudied pathogen [9].

Active invasion into eukaryotic non-phagocytic cells is one of the key virulence-associated phenotypes of all *S. enterica* serovars. This unique capability facilitates *Salmonella* intestinal
epithelium crossing of the small intestine [10] and is mediated by a designated type three
secretion system (T3SS) encoded within the *Salmonella* Pathogenicity Island (SPI)-1. This
sophisticated syringe-like nanosystem is evolutionary related to the flagellar apparatus [11]
and is used to translocate an array of effector proteins directly into the host cell cytoplasm.
These effectors trigger cytoskeletal rearrangements and *Salmonella* penetration of intestinal
barriers [12]. In addition, T3SS-1 and its associated effectors significantly contribute to intesti-
nal inflammation [13] that helps *Salmonella* to compete with the gut microbiota [14].

Within host cells, *Salmonella* are first compartmentalized into a modified intracellular phagosome, known as the *Salmonella* containing vacuole (SCV). Inside the SCV *Salmonella* manipulates the phagosome-lysosome membrane fusion and other cellular pathways, and starts to replicate intracellularly [13]. These activities require a second T3SS encoded by genes on SPI-2 and the translocation of a distinct set of effectors proteins [15]. Later, a subpopulation of this pathogen can escape into the cytosol where environmental conditions are more prone to bacterial growth. This step may lead to hyper-replication, epithelial cell lysis and bacterial release [16].

Typhoid fever disease and host-response to S. Typhi infection is largely attributed to the function of the Vi polysaccharide capsule and its associated regulator, TviA encoded on SPI-7 [17–20]. Nonetheless, since SPA does not harbor the SPI-7, nor expresses the Vi capsule, different mechanisms such as long O antigen chains in the LPS of SPA were suggested to play a role in evasion from the host-immune response and the clinically indistinguishable enteric fever disease caused by SPA and S. Typhi [21,22]. Previously, we showed that in *Salmonella* culture grown in LB to the late logarithmic phase, SPI-1 genes and T3SS-1 effectors are expressed and secreted at significantly lower levels by SPA compared to S. Typhimurium (STM) [23], and demonstrated differences in the regulatory setup of the flagella-chemotaxis pathway between these serovars [24,25].

Here, to further understand the different pathogenicity of SPA vs. STM as prototypic NTS,
we set out to compare their transcriptional landscape during non-phagocytic host cell infec-
tion. Intriguingly, we show that in contrast to STM, intracellular SPA is motile, largely resides
in the host cell cytosol and induces the expression of SPI-1 and flagellar genes. In addition, we
demonstrate that SPA and STM diverge in the intracellular expression of genes for citrate and
ethanolamine metabolism and in their ability to utilize these carbon sources *in vitro*. We
hypothesize that these differences contribute to the distinct diseases resulting from SPA vs.
STM infection, and their different interactions with the human host.

**Results**

Intracellular and extracellular STM and SPA exhibit distinct gene
expression profiles

Although STM and SPA belong to the same species, share about 89% of their genes and express
many common virulence factors [26], the disease they cause in immunocompetent humans is
very different. One possible reason, expected to contribute to their different pathogenicity is a
distinct expression of genes involved in host-pathogen interactions. To test this hypothesis, we
compared the transcriptomic architecture and gene activity of SPA or STM in a relevant con-
text of epithelial cell infection. We chose to sample *Salmonella* infection at 8 h post infection
(p.i.) when the pathogen is considered to be well engaged in replication and optimally adapted
to the intracellular environment [27]. Thus, HeLa cells were infected with fluorescent STM
and SPA for 8 h and sorted by flow cytometry to isolate *Salmonella*-infected cells. Overall,
RNA was extracted from six independent HeLa cell cultures infected with SPA (N = 3) and
STM (N = 3). In each experiment, 2–2.5×10⁶ GFP-positive HeLa cells, (infected with either
STM or SPA) were collected. In parallel, RNA was extracted from six independent extracellular STM (N = 3) and SPA (N = 3) cultures that were grown in LB to stationary phase microaerobically, prior to cells infection. Isolated RNA from these 12 samples was deep-sequenced (RNA-Seq), aligned to the corresponding *Salmonella* genome and compared between all conditions. Principal Component Analysis (PCA) of normalized gene read counts of these RNA-Seq data indicated a unique expression profile for each condition and good reproducibility between independent biological triplicates (Fig 1A).

## The intracellular transcriptomic landscape of S. Typhimurium

RNA extracted from three independent Fluorescence-Activated Cell Sorting (FACS) HeLa cell cultures infected with STM have generated 496–575 million RNA-Seq reads (human and bacterial reads) per experiment. Of which, 4.9–6.6 million sequence reads, assigned to non-rRNA bacterial genes were used for intracellular bacterial transcriptome analysis. The three extracellular STM cultures grown in LB to stationary phase under microaerobic conditions, have generated 4.3–13.8 million RNA-Seq informative reads that were assigned to non-rRNA bacterial genes. To analyze the STM gene expression, we calculated for each gene the number of transcripts per kilobase million (TPM) from its feature read counts. The expression level threshold was set to a TPM value of 10 or above [28]. Altogether, we identified 1,018 distinct genes that were exclusively expressed by STM during intracellular infection of HeLa cells and 191 distinct genes that were expressed only when STM was grown in LB. 3,127 transcripts were commonly expressed at both conditions (Fig 1B and S1 Table).

RNA-Seq successfully identified 365 upregulated and 465 downregulated genes that were changed by at least twofold inside epithelial cells at 8 h p.i. (adjusted p-value ≤ 0.05) by STM, relative to their expression in LB (Fig 2A and S2 Table). In agreement with previous studies [14,29,30], during STM infection, at least 27 genes belonging to the flagella-motility regulon and several genes from the type I fimbriae (*fimAWYZ*) were significantly repressed. Interestingly, the expression of genes encoding the Saf fimbriae (*safABCD*) was upregulated during HeLa cells infection by STM.

In contrast to the flagella and the *fim* genes, at least 31 genes from the PhoPQ regulon were highly induced intracellularly (Fig 2B). The two-component system PhoPQ, orchestrates *Salmonella* adaptation to the intracellular milieu and regulates a wide array of genes central for *Salmonella* virulence including those encoded on SPI-2 [31,32]. Correspondingly, we identified upregulation of the SPI-2 regulon including *ssaBCDEGHIJKLMVNOPRSTU*, *sseABCDEF-GIJL*, *sseK3*, *sspH2*, *pipB*, *pipB2*, and *steC* genes (Fig 2B). These results are consistent with the known notion that SPI-2 expression is induced during intracellular replication of *Salmonella* [15,33] and with previous reports that studied STM expression during macrophages infection [34,35]. Other PhoP-regulated genes including *pgtE* (involved in bacterial resistance to antimicrobial peptides), *pagN* (encoding an adhesion/invasion protein), *phoN* (acid phosphatase) and *pagCDJKO* were also induced (Fig 2B).

Iron is an essential micronutrient required by nearly all bacterial species, including *Salmonella*. Limiting the availability of the nutrient metals to intracellular pathogens is one of the main mechanisms of nutritional immunity. In return, *Salmonella* uses a variety of high-affinity iron uptake systems to compete with the host for essential transition metals required for its intracellular growth [36]. Congruently, we identified elevated expression of multiple iron (*entABEFH* and *iroNBCDE*), iron/manganese (*sitABCD*), zinc (*zntR*, *zitT* and *znuA*), and magnesium (*mgtAB*) acquisition system genes, highlighting the essentiality of these metal ions for STM intracellular replication (Fig 2B).
A

<figure><img src="image_4.png" /></figure>

B

<figure><img src="image_5.png" /></figure>

**S. Typhimurium**

C

<figure><img src="image_6.png" /></figure>

**S. Paratyphi A**

**Fig 1.** The number of STM and SPA genes expressed in LB medium and during intracellular growth. (A) Principal component analysis (PCA) of RNA-Seq data. The PCA was performed by using factomineR package on the normalized expression of 3851 shared genes of three independent SPA and STM cultures grown in LB medium (M1 to M3) and of three independent HeLa cell cultures (C1 to C3) infected with STM and SPA (as described in Materials and Methods). The percentages on each axis represent the percentages of variation explained by the principal components. Points that are closer together are more similar in gene expression patterns, and the barycenter of the triplicates used for each condition is shown
by a larger symbol. (B and C) *Salmonella* RNA-Seq reads were aligned to S. Typhimurium 14028S and S. Paratyphi A 45157 genomes and normalized by transcripts per million (TPM) transformation. Genes with TPM ≥ 10 were considered as expressed. The number of expressed genes in STM (B) and SPA (C) growing in LB medium to stationary phase under microaerobic conditions and within HeLa cells at 8 h p.i. is shown.

https://doi.org/10.1371/journal.ppat.1010425.g001

More than 100 suitable carbon substrates as well as various nitrogen, phosphorus and sulfur sources are available to invading pathogens, at different niches of the vertebrate hosts [37,38]. Interestingly, multiple genes belonging to the ethanolamine operon (*eutQTMGAC*) and genes involved in glycerol (*glpABDK*), and maltose (*malEKM*) metabolism were particularly induced during STM infection (Fig 2B). Nevertheless, the most highly induced metabolic gene was *uhpT*, encoding a hexose phosphate transporter that was upregulated by more than 40-fold at 8 h p.i. These results suggest that sugar phosphates are major substrates required for cytosolic intracellular proliferation of STM within HeLa cells, but also that ethanolamine, glycerol and maltose could be utilized as alternative carbon sources in epithelial cells.

Similarly, biotin (*bioABCDF*), purine (*nrdEFHI* and *purDEHKLT*), and sulfate/thiosulfate (*cysACHITW* and *sufABCDE*) biosynthetic metabolic pathways were highly transcribed intra-cellularly by STM. Moreover, genes involved in biosynthesis of the amino acids histidine (*hisACDFH* and *hutG*), tryptophan (*trpABDE*), and glutamine (*glnAGKL*), and genes that cata-lyze the cleavage of N-acetylneuraminic acid (sialic acid; Neu5Ac) to form pyruvate and N-acetylmannosamine (ManNAc) (*nanAEKT* and *nagB*) were also upregulated by STM during HeLa cell infections (Fig 2B).

A large number of genes involved in phosphate metabolism were identified to be induced within HeLa cells. For example, the transcriptional levels of *pstABCS* (encoding an ATP-dependent phosphate uptake system), which is responsible for inorganic phosphate uptake

<figure><img src="image_7.png" /><figcaption>**Fig 2.** DEGs of intracellular S. Typhimurium compared to growth in LB medium. (A) a Volcano plot showing the fold change (log2 ratio) in the expression of STM genes grown intracellularly in HeLa cells (8 h p.i.) vs. the growth in LB medium (X-axis) plotted against the -log₁₀ adjusted p-value (Y-axis). Each dot on the plot represents the mean value (from three independent cultures) of one gene. Genes that were changed by more than twofold are colored in red. (B) Heatmap based on RNA-Seq results showing the relative transcription of STM genes of interest during intracellular infection (C1 to C3) vs. their extracellular expression in LB medium (M1 to M3).</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1010425.g002
during phosphate starvation, *ugpABCEQ*, *phoR* and *apeE* were markedly increased. The *ugp*
operon is involved in hydrolysis of diesters during their transport at the cytoplasmic side of
the inner membrane, generating sn-glycerol-3-phosphate (G3P) which is used as a phosphate
source [39]. The *upg* operon is regulated by a two-component system, PhoBR, which also regu-
lates the outer membrane esterase, ApeE [40]. Additional genes belonging to the PhoB regu-
lon, such as *phnSVWU* encoding the phosphonatase pathway, which is responsible for
breaking down phosphonate and yielding cellular phosphate [41] were also significantly
induced.

Collectively, these results suggest that STM likely experiences iron, biotin, purine and phosphate deprivation in host epithelial cells and induces the expression of genes involved in their transport and utilization. In addition, these results imply that ethanolamine, glycerol and maltose are available sources of carbon for STM in epithelial cells.

The intracellular transcriptomic landscape of S. Paratyphi A

RNA-Seq analysis was similarly applied for intracellular SPA. 313.8 to 805.9 million RNA-Seq
reads were obtained from three independent FACS-sorted HeLa cell cultures infected with
SPA. Out of these, 14.6 to 36.6 million sequence reads were assigned to non-rRNA SPA genes
and used to determine the intracellular SPA transcriptome (S3 Table). As a reference control,
RNA extracted from extracellular SPA cultures grown in LB, generated 7.3 to 23 million infor-
mative reads (assigned to non-rRNA SPA genes). Overall, we identified 322 genes that were
exclusively expressed (TPM ≥10) by SPA during infection of HeLa cells compared to 576
genes that were expressed only by extracellular SPA in LB culture. 2,995 SPA genes were
expressed at both conditions (Fig 1C and S3 Table).

At 8 h p.i. of HeLa cells, 322 SPA genes were upregulated and 579 genes were down-regu-
lated (Fig 3A and S4 Table) by at least twofold (adjusted p-value ≤ 0.05). Like in intracellular
STM, the genes for type I (Fim) fimbriae (*fimCIWYZ*) and additional fimbria clusters includ-
ing Bcf (*bcfDE*), the Curli (*csgBCDE*), Stb (*stbBE*), Stc (*stcABC*), Stf (*stfAD*), Sth (*sthABCD*),
and Stk (*stkABC*) were all strongly downregulated intracellularly (Fig 3B). Noteworthy, the
downregulation of type I (*fimACIZ*), Curli (*csgBDE*) and Bcf (STM14_0037 and *bcfB*) fimbriae
was more pronounced in SPA than in STM (S1 Fig), indicating that the SPA fimbriome is
robustly repressed during intracellular replication, while some basal level of expression of
these fimbrial genes is still maintained during STM infection.

In contrast, but consistent with our findings for intracellular STM, genes belonging to the
PhoPQ (including polymyxins resistance) and SPI-2 regulons were significantly induced by
intracellular SPA. Moreover, various transition metals import systems including manganese
and iron transporter (*sitABCD*), enterobactin biosynthesis system (*entABCDEF* and *fepAG*)
and the iron ABC transporter IroC were considerably upregulated during HeLa cell infection
by SPA (Fig 3B).

In addition, the expression of several metabolic pathways that were significantly elevated
during STM infection were induced by intracellular SPA as well, including biotin (*bioABCDF*)
and purine (*purCDEHKLT*) biosynthesis genes. Further overlap in the intracellular transcrip-
tome of these pathogens includes several amino acid biosynthesis pathways such as glutamine
(*glnAGKLPQ*), tryptophan (*trpABCDE*) and histidine (*hisAFGHI*) that were upregulated in
both intracellular STM and SPA.

Nevertheless, comparison of the intracellular transcriptome between STM and SPA, identi-
fied 234 differentially expressed genes (DEGs; S5 Table), indicating that some key pathways
are distinctly expressed during the intracellular phase of STM vs. SPA. Such differences
included expression of several amino acids biosynthesis genes, which was more pronounced in
<figure><img src="image_8.png" /><figcaption>**Fig 3. DEGs of intracellular S. Paratyphi A compared to growth in LB medium.** (A) a volcano plot showing the fold change (log2 ratio) in the expression of SPA genes grown intracellularly in HeLa cells (8 h p.i.) vs. the growth in LB medium (X-axis) plotted against the -log₁₀ adjusted p-value (Y-axis). Each dot on the plot represents the mean value (from three independent cultures) of one gene. Genes that were changed by more than twofold are colored in red. (B) Heatmap based on RNA-Seq results showing the relative transcription of SPA genes of interest during intracellular infection (C1 to C3) vs. their extracellular expression in LB medium (M1 to M3).</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1010425.g003

intracellular SPA than in intracellular STM. These include biosynthesis genes of glycine (*glyASU*), isoleucine (*ilvCDEG*), leucine (*leuABCD*), lysine (*lysAC* and *sucBD*), methionine (*metCEFHKLR*), serine (*serAC*), and threonine (*thrBC*). In contrast, the threonine degradation pathway (*tdcABCDEG*) was found to be repressed in intracellular SPA. These results may suggest increased protein synthesis in SPA vs. STM during epithelial cell infection. Additional pathways that were upregulated in SPA, but not in STM are the oxidative phosphorylation (*nuoEFGHIJKLMN* and *atpCD*), and peptidoglycan cell wall formation (*murACDEG*) pathway genes (Fig 3B). Remarkably, further differences between the intracellular transcriptome of SPA vs. STM included different expression profiles of the SPI-1, motility and chemotaxis, and carbon utilization regulons as detailed in the following sections.

## Differences in carbon catabolism during intracellular infection of S. Paratyphi A and S. Typhimurium

In order to survive, persist and replicate in host cells, intracellular pathogens must adapt their metabolic pathways to the specific physical conditions and available nutrients found in the intracellular environment. Carbon catabolism provides bacteria with energy by means of reducing equivalents, ATP and essential biosynthetic precursors. Different studies have suggested that hexose monosaccharides such as glucose and glucose-6P are the main source of carbon during *Salmonella* infection [42–46]. Comparison of the metabolic gene expression between intracellular STM and SPA revealed significant differences in the expression profile of several carbon catabolic pathways. For example, the gene *uhpT* encoding an antiporter for external hexose 6-phosphate and internal inorganic phosphate [47], presented about 90-fold higher expression by intracellular STM than SPA (Fig 4A). This different expression might be
<figure><img src="image_9.png" /><figcaption>**Fig 4.** Differences in carbon utilization gene expression between STM and SPA. (A) Heatmap of carbon catabolism gene expression by intracellular STM vs. intracellular SPA. Pseudogenes in the SPA genome are marked with an asterisk. (B) The fold change in the transcription of the glucose-6P (*uhpT*), citrate (*citA* and *citX*), and ethanolamine (*eutP*, *eutQ* and *eutG*) utilization pathways in intracellular STM relative to their expression in intracellular SPA was determined using qRT-PCR. RNA was extracted from FACS-sorted intracellular salmonellae 8 h p.i. The mean and the standard error of the mean (SEM) of 3–8 independent qRT-PCR reactions are shown. Student's t-test was used to determine statistical significance. **, p < 0.01; ***, p < 0.001; ND, not detected. (C) Growth of STM and SPA in M9 in the presence (+) or absence (-) of casamino acids (AA) and ethanolamine (ETA), under anaerobic conditions. (D) Growth of STM and SPA in M9 supplemented with citrate as a sole carbon source.</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1010425.g004

because the sensor-regulator system encoded by *uhpABC* is defective by pseudogene formation in SPA [48]. Similarly, the gene encoding for galactonate permeases (*dgoT*) was upregulated in intracellular STM, but not in SPA, suggesting that STM utilizes more hexose 6-phosphate and galactonate than SPA as a carbon source during intracellular growth.

Similarly, RNA-Seq identified that at least seven genes involved in citrate metabolism including *citCDEFTX* and *citC2* exhibited higher transcription levels in intracellular STM than SPA (Fig 4A). Independent qRT-PCR further confirmed that the intracellular expression of *citA* and *citX* encoding citrate-proton symporter and apo-citrate lyase phosphoribosyl-dephospho-CoA transferase, respectively is lower by about 3- and 10-fold in intracellular SPA vs. STM (Fig 4B).

STM can utilize ethanolamine as a sole source of carbon, nitrogen, and energy in a cobalamin (vitamin B12)-dependent manner. Ethanolamine may be important to STM growth in the host, since it is derived from the membrane phospholipid phosphatidylethanolamine that is particularly prevalent in the gastrointestinal tract, in a process involving breakdown of the precursor molecule phosphatidylethanolamine by phosphodiesterases to glycerol and ethanolamine [49]. In STM, ethanolamine catabolism involves 17 genes organized in the *eut* operon, controlled by the transcriptional activator EutR [50]. At least six ethanolamine catabolism
genes (*eutSPQTGH*) were expressed at higher levels by intracellular STM compared to SPA (Fig 4A). These differences were independently verified by qRT-PCR that showed significantly reduced expression of *eutG*, *eutP* and *eutQ* genes in intracellular SPA vs. STM (Fig 4B).

Collectively, these observations suggest that STM is able to better utilize intracellular citrate, ethanolamine and galactonate as intracellular carbon sources than SPA, and that SPA may use an alternative carbon source such as amino acids during intracellular growth.

To further test this hypothesis, we compared the growth of SPA and STM in a defined M9 medium supplemented with ethanolamine or citrate as a carbon source. Growth under anaerobic conditions in M9 medium supplemented with ethanolamine, without amino acids, did not allow bacterial growth of neither STM nor SPA. Adding casamino acids, without additional carbon source, resulted in a minimal and comparable growth of both STM and SPA, likely due to the utilization of amino acids as a restricted carbon source. Nevertheless, adding ethanolamine as a *bona fide* carbon source to the medium resulted in increased and rapid growth of STM, but only a minor improvement of SPA growth (relative to its growth on amino acids only; Fig 4C), indicating the ability of STM, but not SPA, to efficiently ferment ethanolamine. Replacing the carbon source with citrate under aerobic growth conditions allowed delayed, but substantial growth of STM, in comparison to restricted and slower growth of SPA (Fig 4D). These results are concurring with the transcriptomic data and indicate a better ability of STM than SPA to utilize both ethanolamine and citrate as a sole carbon source. These differences may also contribute to a higher intracellular replication pace of STM than SPA found in HeLa cells (S2 Fig).

**Intracellular S. Paratyphi A expresses SPI-1 genes at higher levels than intracellular S. Typhimurium**

Previously we showed that in LB medium, under aerobic growth conditions to the late logarithmic phase, SPI-1 gene expression, as well as secretion of SPI-1-T3SS effector proteins, occur at significantly lower levels in SPA compared to STM [23]. Here, in sharp contrast, we found that intracellular SPA transcribed significantly higher levels of SPI-1 and T3SS-1-effector genes than intracellular STM (Fig 5A). This was clearly evident (log 2 change ≥ 2 and p-value ≤ 0.05) for at least 14 SPI-1 genes including the SPI-1 regulators *hilA* and *hilD*, the effector genes *sopBDF*, the T3SS-1 structural genes *invCHIJ* and *spaOPQ*, *sigE* encoding a chaperon (a.k.a. *pipC*) and *iagB* (S5 Table).

Independent qRT-PCR was used to verify these observations and demonstrated that intracellular SPA transcribed *hilA*, *hilD*, *invH*, *sopB*, *sopD* and *spaO* at 4- to 27-fold higher levels than STM (Fig 5B). Moreover, Western blotting against a 2HA-tagged SopB and SopE2 further indicated elevated translation of these two SPI-1 effectors in intracellular SPA vs. STM (Fig 5C). Collectively, these results indicate that SPA expresses SPI-1 genes at higher levels during infection of non-phagocytic cells, in comparison to intracellular STM.

**Intracellular S. Paratyphi A is motile and expresses flagella-chemotaxis genes at higher levels than S. Typhimurium**

Functional flagella are necessary for *Salmonella* motility and invasion of host cells [24,51,52]. In *S. enterica*, more than 50 genes are involved in the assembly and motion of flagella, organized in one regulon that composes three ordered regulatory phases corresponding to classes I (early), II (middle), and III (late) [53].

Previous studies have shown that the genes coding for the flagellar machinery, as well as many genes involved in chemotaxis, are strongly repressed during STM intracellular lifestyle [30,34,35]. Here, RNA-Seq analysis has demonstrated that at least 38 out of 50 chemotaxis and
<figure><img src="image_10.png" /><figcaption>**Fig 5.** *Intracellular SPA expresses higher levels of SPI-1 genes than intracellular STM.* (A) Heat map of RNA-Seq results showing the relative transcription of SPI-1 genes in three independent HeLa cell infections with STM and SPA. (B) The fold change in the transcription of six SPI-1 genes (*hilA*, *hilD*, *invH*, *sopB*, *sopD* and *spaO*) in intracellular SPA relative to their expression in intracellular STM was determined by qRT-PCR. RNA was extracted from FACS-sorted intracellular salmonellae 8 h p.i. The results show the mean of 3–6 qRT-PCR reactions and the error bars indicate the SEM. (C) SPA and STM strains expressing GFP and a 2HA-tagged version of the SPI-1 effectors SopE2 and SopB were used to infect HeLa cells. At 8 h p.i., cells were fixed and FACS-sorted for infected (GFP-positive cells). Equal number of infected HeLa cells (1×10⁶ for each sample) were collected and separated by 12% SDS-PAGE. Western blotting using anti-hemagglutinin (HA) tag antibodies was used to detect the intracellular expression of SopE2 and SopB. Anti-DnaK antibodies were used as a loading control.</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1010425.g005

flagella genes, belonging to class II and III of the flagella regulon are expressed at significantly higher levels in intracellular SPA compared to STM (Fig 6A and S5 Table). These include the flagellin gene *fliC* that was expressed at 11-fold higher levels by SPA vs. STM, and many other flagellar genes including *flgC* that demonstrated 64-fold difference in its expression between intracellular SPA and STM. In addition to the flagellar genes, SPA residing in the host cell presented much higher expression levels of chemotaxis genes (e.g. *cheB*, *cheR*, *cheM*, *cheW*, *cheZ*, and *cheY*) than intracellular STM.

Independent qRT-PCR analysis confirmed these results and showed that in sorted infected HeLa cells, the expression of *fliC*, *fliZ*, *flgM*, *flgF*, *cheA* and *cheB* was 8- to 130-fold higher in intracellular SPA than in intracellular STM (Fig 6B). In agreement with the RNA-Seq and
<figure><img src="image_11.png" /><figcaption>**Fig 6.** Intracellular SPA expresses higher levels of motility-chemotaxis genes than intracellular STM. (A) Heat map of RNA-Seq results showing the relative transcription of the motility-chemotaxis regulon in three independent HeLa cells infections with STM and SPA. (B) The fold change in the transcription of six flagellum genes (*fliC*, *fliZ*, *glgM*, *flgF*, *cheA* and *cheB*) in intracellular SPA relative to their expression in intracellular STM was analyzed by qRT-PCR. RNA was extracted from FACS-sorted intracellular salmonellae 8 h p.i. The results show the mean of 3–5 independent reactions and the error bars indicate the SEM. (C) SPA and STM strains expressing GFP were used to infect HeLa cells. At 8 h p.i., cells were fixed and FACS-sorted. Equal numbers of 1×10⁶ infected HeLa cells were collected, lysed, and proteins were separated on 12% SDS-PAGE. Western blotting using anti-FliC antibodies were used to detect synthesis of FliC by intracellular *Salmonella*. Anti-DnaK antibodies were used as a loading control.</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1010425.g006

qRT-PCR analyses, Western blotting using anti-FliC antibodies demonstrated significantly higher translation of FliC in intracellular SPA vs. STM (Fig 6C). Taken together, these findings uncover significantly higher expression levels of motility and chemotaxis genes by SPA vs. STM within infected HeLa cells.

To further analyze if the increased expression of motility genes by intracellular SPA has physiological consequences, we investigated the possible presence of flagella in intracellular SPA. Immunostaining of intracellular SPA in either HeLa (Fig 7A) or Caco-2 cells (Fig 7B) at 8 or 4 h p.i., respectively, indicated the presence of flagella filaments on intracellular SPA. Flagella staining was absent for a Δ*fliC* strain of SPA, and restored in the plasmid-complemented SPA Δ*fliC*. No immunolabeling of flagella was observed for intracellular STM WT or STM Δ*fliC* Δ*fljB* (Fig 7A), indicating that this phenotype is SPA-specific.
<figure><img src="image_12.png" /><figcaption>B Caco-2</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Fig 7.** Formation of flagella by SPA within the host cell cytoplasm. (A) HeLa cells stably expressing LAMP1-GFP (green) were infected with SPA WT, SPA *fliC*, SPA *fliC* [*fliC*], STM WT, or STM *fliC* *fljB* as indicated, each constitutively expressing mCherry (red). Infected cells were fixed 8 h p.i. and subsequently immunolabeled with serovar-specific antibody against FliC (blue). (B) Caco-2 well were infected with SPA WT constitutively expressing mCherry (red). Cells were fixed 4 h p.i. and immunolabeled for FliC (blue). Scale bars: 10 µm.</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1010425.g007
Differences in intracellular compartmentation and SIF-formation between STM and SPA

Phenotypic comparison by live-cell imaging between host cells infected with STM or SPA was
performed. To distinguish *Salmonella* in the SCV from *Salmonella* in other compartments or
in the cytosol, and to enable live-cell imaging of infected host cells, we used LAMP1-GFP
expressing HeLa host cells (Fig 8, Movies 1, 4). Lifeact-GFP expressing HeLa cells were used to
investigate potential interactions with the host cell actin cytoskeleton (Movie 1–3).

We quantified the intracellular phenotype of SPA in comparison to STM after infection of HeLa LAMP1-GFP cells using live-cell imaging. Representative recordings of live-cell imaging sessions are shown in Movies 5 and 6 for SPA WT and STM WT, respectively. Quantification revealed that more than 80% of STM WT were present in SCV with associated *Salmonella*-induced filaments (SIFs; Figs 8A and 9). STM WT was observed residing in SCV and the majority of infected host cells showed formation of extensive SIF networks. The STM *ssaV* strain was typically present in SCV, but tubular endosomal filaments were almost absent. A smaller proportion of host cells (ca. 10%) harbored cytosolic STM, and for these cells, massive cytosolic proliferation or hyper-replication was evident (Fig 8A). For intracellular SPA WT, more heterogeneous phenotypes were observed (Fig 8B–8E). SPA was located in the SCV, or present in the cytosol without association with LAMP1-positive membranes. Intracellular SPA also induced SIF in a T3SS-2-dependent manner, but the frequency of SIF-positive cells was much lower with ca. 9% for SPA WT compared to more than 80% for STM WT (Figs 8B, 8C, and 9). SPA-infected host cells often contained thin SIF, that are likely composed of single-membrane tubules, in contrast to SIF with more intense LAMP1-GFP signal, indicating double-membrane tubules (Fig 8B and 8C). For SPA WT, SIF formation was significantly less frequent, yet dependent on T3SS-2. The frequency of SPA associated with SCV marker LAMP1 was highly reduced compared to STM. As intracellular SPA was observed to reside in tightly enclosing SCV [54], we considered flagella expression by vacuolar SPA less likely than by cytosolic bacteria. The transcriptional analyses showing expression of flagellar genes by intracellular SPA, detection of flagella filaments on intracellular SPA, and cytosolic SPA in a major proportion of infected HeLa cells (~ 40%) all pointed to functional role of flagella expression. Thus, we next investigated potential function of flagella expressed by intracellular SPA by live-cell imaging of *Salmonella*-infected host cells (Movies 1–4). Indeed, we observed actively motile intracellular SPA in about 40% of the infected host cells (Fig 9).

As motility in intracellular microbes is typically mediated by polymerization of host cell
actin by microbial surface proteins [55], we investigated F-actin dynamics in SPA-infected
HeLa Lifeact-GFP cells (Movies 2 and 3). No association of F-actin with static or motile intra-
cellular SPA was evident, and F-actin trails as observed for intracellular motile *Listeria mono-*
*cytogenes* [56] were absent in SPA infected cells. No cytosolic motility was observed for
intracellular STM. To control if cytosolic motility is linked to the observed flagella expression,
we used a SPA *fliC* mutant strain. This non-motile mutant completely lacked cytosolic motil-
ity, but showed increased presence in SCV (Fig 9, Movie 1). Cytosolic motility of SPA was not
affected in the T3SS-deficient *ΔssaR* strain (Fig 9, Movie 1). Live-cell imaging at high frame
rates (Movies 2 and 3) demonstrated the dynamics of intracellular motility of SPA in real time,
with characteristic alternations between straight swimming and stops with tumbling. Analyses
of the vectors of individual SPA indicated collisions with the host cell plasma membrane, but
there were no indications for formations of protrusions, showing no evidence for lateral cell-
to-cell movement as observed for intracellular motile *L. monocytogenes* [56].

We next set out to analyze the broader relevance of intracellular motility of SPA and investi-
gated a set of four additional clinical SPA isolates [25], and two SPA reference strains of SARB
<figure><img src="image_14.png" /><figcaption>**Fig 8.** Distinct phenotypes of intracellular SPA. HeLa cells stably transfected for expression of LAMP1-GFP (green) were infected with STM WT (A, MOI 25) or SPA WT (B-E, MOI 75), each harboring pFPV-mCherry for constitutive expression of mCherry (red). Live-cell imaging was performed at 8 h p.i. Images are shown as maximum intensity projection. Arrowheads indicate subpopulations with distinct phenotypes in infected cells: SIF (green), thin SIF (green with asterisk), *Salmonella* residing in SCV (yellow), cytosolic *Salmonella* (red), hyper-replicating cytosolic *Salmonella* (red with asterisk). Scale bars: 10 µm (overview), 5 µm (inset).</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1010425.g008
<figure><img src="image_15.png" /><figcaption>**Fig 9.** Quantification of intracellular phenotypes. HeLa cells stably expressing LAMP1-GFP (green) were infected with *Salmonella* strains expressing mCherry (red). Quantification of intracellular phenotypes was performed with a Zeiss Cell Observer with 63x objective. From 6–8 h p.i. at least 100 infected cells per strain (for SPA *fliC* at least 50) were examined for intracellular phenotype. Data represent the means with standard deviations (SD) for three biological replicates for each strain and condition. Statistical analysis was performed with One-way ANOVA and is indicated as n.s., not significant; *, p < 0.05; **, p < 0.01; ***, p < 0.001.</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1010425.g009

collection [57]. We observed that all seven SPA strains tested showed intracellular motility in HeLa cells (Movie 4). Expression of flagella genes and presence of flagella filaments was previously demonstrated for STM in dead host cells extruded from polarized epithelial monolayers [58]. As these experiments were performed with STM strain SL1344, we also analyzed this strain in our experimental setting. As for STM ATCC 14028, also STM SL1344 showed cytosolic hyper-replication in infected HeLa cells, but cytosolic motility was absent (Movie 4).

To investigate the intracellular motility of SPA in other host cell types, infection of human colonic epithelial cell line Caco-2 and human monocyte cell line U937 by SPA WT was analyzed. We observed intracellular motility of SPA WT in Caco-2 (Movie 7), a result in line with the detection of flagella (Fig 7B). Motile SPA WT were also detected in U937 (Movie 8). Due to the localized distribution of motile SPA in U937 and the frequent association with endosomal cargo Dextran-Alexa647, we conclude that motility of SPA in U937 macrophages occurs within spacious phagosomes rather than in the cytosol of host cells.

These data confirm that the expression of flagella by intracellular SPA and its localization outside of the SCV leads to flagella-mediated motility, and that these phenotypes are conserved traits specific to SPA and distinct from the cellular phenotypes characterizing intracellular STM.

## Intracellular SPA are primed for reinfection of naïve host cells

We set out to investigate potential functional consequences of the increased expression of SPI-1 and motility genes by intracellular SPA. We hypothesized that increased expression of SPI-1, higher amounts of T3SS-1 effector proteins, and motility could affect the interaction of SPA released from infected cells with naïve host cells. As host cell decay and release of the intracellular *Salmonella* is asynchronous, and quantification of new infection events is difficult in genetic protection assays, we set up a reinvasion assay that provides experimental control.
<figure><img src="image_16.png" /><figcaption>**Fig 10.** SPA is more invasive after release from infected host cells. (A) HeLa cells were infected with STM WT or SPA WT at MOI 5, and invasion by each strain was determined by gentamicin protection assays. At 8 h p.i., one portion of infected cells was lysed by hypotonic conditions, and the released bacteria were used to reinfect naïve HeLa cells. (B) Intracellular CFU counts were determined 1 h p.i. after initial infection (invasion, black bars), or 1 h p.i. after infection with *Salmonella* released from host cells at 8 h p.i. (reinvasion, grey bars). The percentage of *Salmonella* in inoculum able to invade host cells was calculated. (C) Multiplicity of infection (MOI) calculated for invasion and reinvasion experiments shown in (B). For invasion, MOI was adjusted based on cell density of the inoculum and controlled by CFU determination, while MOI for reinvasion was calculated based on CFU determination of lysates of infected cells. Depicted are means and SD of three biological replicates, each performed in triplicate. For statistical analysis, Student's t-test was performed. Significance is indicated as n.s., not significant; *, p < 0.05; **, p < 0.01.</figcaption></figure>

https://doi.org/10.1371/journal.ppat.1010425.g010

HeLa cells were infected by invasive *Salmonella*, at a defined time point of 8 h p.i., host cells were osmotically lysed, and released bacteria were used to reinfect new epithelial cells. The number of viable bacteria in inoculum, released from initially infected host cells, and the number of bacteria in the second round of invasion were quantified (Fig 10A).

For this reinvasion analysis we have used SPA culture that was grown to stationary phase microaerobically and STM culture that was grown to the late logarithmic phase under aerobic conditions. These growth conditions have resulted in a similar initial infection rate for SPA and STM (3.25% and 3.41% of inoculum, respectively). Interestingly, reinfection assays indicated that 6.88% of the SPA released from HeLa cells at 8 h p.i. were capable of reinvasion to naïve HeLa cells. In contrast, only 0.91% of released STM were able to reinvade, and reinvasion of STM after microaerobic culture was lower with 0.75% (Fig 10B). These data indicate that higher expression of SPI-1 genes and bacterial motility by intracellular SPA results in increased invasiveness of SPA following exit from host cells.
Discussion

Despite high genetic similarity, typhoidal and non-typhoidal *Salmonella* serovars differ hugely in their host specificity, the immune response they elicit in humans, and the clinical manifestations they cause [59]. Previously, we showed differences in host cell invasion and SPI-1 expression between STM and SPA, following extracellular growth in rich LB medium under aerobic conditions [23] and demonstrated that elevated physiological temperature affects differently their motility and host cell entry [25,54].

In this work, we focused on their intracellular phase of infection and have taken a deep
RNA-Seq approach [60] to gain a systematic overview on the intracellular transcriptomic sig-
nature of STM and SPA, and to identify differences in the activity of their genes in the intracel-
lular milieu. Following transition from extracellular to the intravacuolar environment,
*Salmonella* undergoes extensive adaptation and global modulation of gene expression to
respond to the intracellular environment. However, the nutritional contents of the SCV and
the metabolic pathways required during intracellular growth are largely undefined. Amino
acids and purines seem to be deficient since previous studies have shown that many auxotro-
phic *Salmonella* mutants are unable to grow intracellularly and are severely attenuated for vir-
ulence in mice [61–64]. These observations suggest that intracellular *Salmonella* depend on de
*novo* synthesis of metabolic precursors.

Recently, a similar RNA-Seq approach was applied to determine the transcriptome of STM populations that are either free-living in the cytosol or residing within the SCV of HeLa cells 8 h p.i. [65]. Earlier studies have implemented a proteomic approach to determine the expres- sion of STM proteins within HeLa cells at 6 [66] and 18 h [30] p.i. Also, the STM transcriptome inside macrophages [34,35] and epithelial cells [29] has been reported previously. Nevertheless, a global transcriptomic comparison between SPA and STM during intracellular infection has not been published so far.

Previous efforts to identify changes in gene expression using quantitative proteomic approach (LC-MS/MS analyses) have revealed that among about 3,000 detected proteins, the expression of 100 *Salmonella* proteins significantly changed (61 were upregulated and 39 proteins were downregulated) during the transition of STM from the extracellular to the intracellular niche of HeLa cells, at 6 h p.i. [66]. A microarray hybridization-based technique has identified 605 upregulated and 616 downregulated STM genes at 6 h p.i. of HeLa cells [29]. Our current transcriptomic analysis identified 830 differentially expressed STM genes (365 upregulated and 465 downregulated) that were altered at 8 h p.i., relative to their extracellular expression. Similar numbers of DEGs were also detected for SPA, which presented 322 upregulated and 579 downregulated genes, during epithelial cell infection.

In agreement with these previous analyses, we observed substantial downregulation in the
expression of genes associated with SPI-1, Type I fimbriae, flagella and chemotaxis systems in
STM, but significant upregulation of STM genes involved in metal ions uptake, purines and
biotin synthesis. In fact, nearly all known *Salmonella* iron-, zinc-, or manganese-responsive
pathways were found to be intracellularly induced. Since the SCV is characterized by low levels
of magnesium, manganese, iron and zinc metals [67] and given the pivotal role of these metal
ions in *Salmonella* infection [68,69], upregulation of iron, zinc and manganese acquisition sys-
tem genes is not surprising. Nonetheless, these results highlight the importance of biotin and
purines to the intracellular phase of *Salmonella* during infection. This notion is consistent with
(i) SCVs of both macrophage and epithelial cells are limited for purines and pyrimidines [70];
(ii) the avirulent phenotype of purine auxotrophic mutants of S. Typhimurium and S. Dublin
in mice [63]; and (iii) the reported contribution of biotin sulfoxide reductase to oxidative stress
tolerance and virulence of *Salmonella* in mice [71].
Interestingly, in contrast to the proteomics [66] and the microarray-based transcriptome [29] analyses, we were able to identify significant upregulation of SPI-2 and PhoPQ regulon genes both in STM and SPA. These differences may reflect the time point differences (8 h p.i. in our analysis vs. 6 h p.i. in previous analyses), or the superior sensitivity of the current RNA-Seq techniques.

While the change in the expression of metal transport systems, purines, biotin, PhoPQ and SPI-2 regulons was common to both SPA and STM, intriguing differences in the expression profile of SPI-1, motility and chemotaxis, and carbon utilization pathways (mainly citrate, galactonate and ethanolamine) were identified between intracellular SPA and STM. Ethanolamine is a ubiquitous molecule within vertebrate hosts and serves as a carbon and nitrogen source for bacteria in the intestine, as well as within epithelial cells [72]. *Salmonella* can utilize ethanolamine as an exclusive source of carbon and nitrogen, generating acetaldehyde and ammonia. The ammonia can then provide a cellular supply of reduced nitrogen, while the acetaldehyde can be converted into the metabolically useful compound acetyl-CoA [73]. Our data indicated much higher expression of multiple ethanolamine utilization genes by STM vs. SPA during host cell infection and more efficient utilization *in-vitro*. Recently, it was shown that ethanolamine metabolism in the intestine enables STM to establish infection and that EutR, the transcription factor of the ethanolamine cluster in *Salmonella* directly activates the expression of SPI-2 in the intracellular environment [74,75].

The sugar acid galactonate is another example of an available carbon and energy source that is produced in the mammalian gut. The gut microbiota releases sugar acids including galactonate from ingested polysaccharides, or from simple sugars as catabolic intermediates of metabolism [76]. These results may suggest that STM utilizes galactonate, citrate and ethanolamine more effectively than SPA as available carbon sources, and that STM may better coordinate ethanolamine metabolism and virulence via EutR during infection. This notion is consistent with our *in-vitro* data showing that while STM is able to utilize ethanolamine and citrate as a sole carbon source, SPA cannot grow or grows much slower on these sources. Previous studies have shown that deficiency for citrate lyase activity (*citF*) leads to significant reduction in STM survival inside macrophages [77] and that an STM mutant that cannot utilize citrate is highly attenuated in the mouse [78]. Similarly, ethanolamine is released by host tissues during inflammation and a study by Thiennimitr and colleagues in mice has shown that ethanolamine metabolism provides a growth advantage to STM during intestinal colonization [75]. Since SPA provokes a noninflammatory disease, it is possible that SPA has lost the ability to utilize ethanolamine as it provides no advantage during noninflammatory infection.

It is well established that the expression of the flagella and chemotaxis are highly repressed by intracellular STM [30,34,66,79], supporting a nonmotile state of intracellular STM. In agreement with these reports, we have also observed extensive repression of the entire motility-chemotaxis regulon in intracellular STM, however we found significantly higher expression of this regulon in intracellular SPA. Similar observations were also obtained for the evolutionary related T3SS-1 genes. Most salmonellae are motile and express peritrichous flagella around its membrane. Like T3SS-1, *Salmonella* motility is playing an important role in host colonization and virulence [80,81]. Previous reports have shown a regulatory association between T3SS-1 and the flagella in STM [82]. Subsequent to STM entry into host cells, the intracellular expression of T3SS-1 [30,34,83] as well as flagella production and chemotaxis [30,34,66,84] genes were shown to be highly repressed, possibly to prevent inflammatory response by the NAIP/NLRC4 inflammasome [85]. Therefore, it was unexpected to find that the flagella regulon is readily expressed during SPA intracellular infection and that it maintains an intracellular motility during epithelial cell infection.
Intracellular motility is long known in bacterial pathogens such as *Listeria monocytogenes*, *Shigella flexneri*, *Rickettsia* spp., and *Burkholderia* spp. that employ intracellular actin-based motility to spread from cell to cell [86]. Nonetheless, data indicating a flagellum-dependent intracellular motility of a typhoidal *Salmonella* serovar in living host cells, were unexpected. While we did not find evidence for lateral intercellular spread, it is possible that this mechanism facilitates escape of SPA from the SCV to the cytoplasm. In agreement with this possibility, we identified that a significant portion of intracellular SPA were motile and resided in the cytosol, compared to STM that were mainly non-motile and compartmentalized in SIF-associated SCVs. Previous studies have shown that a small subpopulation of hyper-replicating STM resides within the cytosol of epithelial cells and serves as a reservoir for dissemination. These motile bacteria retain both SPI-1 and flagella expression, but present a low expression of SPI-2 genes, and therefore are primed for invasion [65,87,88], and fuels high luminal density that facilitates fecal shedding [89]. This expression profile is somewhat similar to the identified intracellular profile of SPA. Indeed, we showed that intracellular SPA possesses higher ability to reinfect non-phagocytic cells than intracellular STM. We propose that elevated expression of motility-chemotaxis and SPI-1 genes, in conjunction with flagella-mediated motility, prime SPA toward a new cycle of host cell infection, supporting systemic dissemination of this pathogen in the human body. Fig 11 summarizes the main differences found in this study between intracellular SPA and STM. These differences in gene expression and cellular phenotypes may result from the distinct intracellular localization of SPA vs. STM (i.e. cytosolic localization vs. SCV), but may also contribute to the invasive nature of paratyphoid fever and the ability of SPA to disseminate beyond the intestinal lamina propria. Future work is expected to reveal which cells or tissues of infected host organisms are affected by intracellular motility of SPA and what is the effect of elevated expression of FliC and T3SS-1 components by intracellular SPA on inflammasome activation and host cell death pathways. Nonetheless, such analyses are hampered by the lack of animal models for SPA infection. Infection models with human organoids [90], for examples based on gall bladder epithelial organoids [91] as a tissue type important for SPA persistence and dissemination, may offer new perspectives for analyses of cellular interactions.

## Material and methods

### Bacterial strains, growth conditions and data sets

Bacterial strains utilized in this study are listed in S6 Table. S. Typhimurium 14028 S [92] and S. Paratyphi A strain 45157, an epidemic strain responsible for a paratyphoid outbreak in Nepal [9] were used as the wild-type strains. Bacterial cultures were routinely maintained in Luria-Bertani (LB; BD Difco) liquid medium at 37°C. An M9 medium [33.7 mM Na₂HPO₄, 22 mM KH₂PO₄, 8.5 mM NaCl, 18.7 mM NH₄Cl, 12.1% Trizma base, pH = 7.4] supplemented with 2 mM MgSO₄, 0.3 mM CaCl₂, 3.7 µM thiamine, 4.1 µM biotin, 0.134 mM EDTA, 31 µM FeCl₃, 6.2 µM ZnCl₂, 7.6 nM CuCl₂, 4.2 nM CoCl₂, 1.62 µM H₃BO₃, 8.1 nM MnCl₂, with or without 0.1% Casamino acids, and with one of the following carbon sources glucose (11 mM), ethanolamine hydrochloride (51 mM) or sodium citrate (34 mM) was used as a minimal growth medium.

For HeLa cell infections, *Salmonella* cultures were grown to the stationary phase under microaerobic conditions (see below). The final dataset for RNA-Seq included 12 samples resulting from three independent cultures in LB or infection experiments that were used for RNA extraction: (i) S. Typhimurium 14028 S bacteria grown in LB medium (hereafter named STM_m) (ii) S. Paratyphi A strain 45157 grown in LB medium (hereafter named SPA_m) (iii) Typhimurium 14028 S bacteria-infected HeLa cells (hereafter named STM_i) and (iv) S. Paratyphi A strain 45157 bacterium in HeLa cells (hereafter named SPA_i).
<figure><img src="image_17.png" /></figure>

**Fig 11.** A model summarizing the identified differences between intracellular SPA and STM in epithelial cells.
After invasion, SPA is present in an SCV (A). SPA can either stays in a SIF-less SCV or becomes cytosolic (B).
Cytosolic SPA express high levels of flagella-chemotaxis, amino acids biosynthesis, and SPI-1 genes (C). This
expression profile results in pronounced intracellular motility and enhanced reinvasion capabilities. In contrast,
intracellular STM is almost exclusively present in SIF-associated SCV and express higher levels of citrate, galactonate
and ethanolamine metabolic genes relative to intracellular SPA (D).

https://doi.org/10.1371/journal.ppat.1010425.g011

Infection of HeLa cells

Human epithelial HeLa (ATCC CCL-2) cells were purchased from the American Type Culture
Collection and were cultured in a high-glucose (4.5 g/l) DMEM supplemented with 10% FBS,
1 mM pyruvate and 2 mM L-glutamine at 37°C in a humidified atmosphere with 5% CO2.
Cells were seeded at 6×10⁶ cells/ml in a onewell plate (Greiner) tissue culture dish 18 h prior to
bacterial infection and infected at multiplicity of infection (MOI) of ~200 (bacteria per cell).
For RNA extraction, five plates for S. Typhimurium and six plates for S. Paratyphi A were
used. *Salmonella* strains carrying the pBR-GFP2 plasmid were grown in 2 ml LB containing
20 µg/ml tetracycline at 37°C with shaking (250 rpm) for 6–8 h, and then diluted 1:75 into 10
ml LB supplemented with 20 µg/ml tetracycline and grown without shaking at 37°C for 16 h.
Bacterial cultures were centrifuged for 5 min at 4,629 × g, at room temperature (RT), and
resuspended in 3 ml prewarmed complete DMEM. 3 ml of DMEM suspended salmonellae
(~2.5×10⁹ CFU) were used to replace the medium of seeded cells and plates were centrifuged for 5 min at 142 × g at RT. After centrifugation, 7 ml of prewarmed DMEM were added and plates were incubated at 37°C under 5% CO₂ atmosphere. At 2 h p.i. the medium was aspirated, cells were washed three times with 10 ml PBS containing CaCl₂ and MgCl₂ (PBS +/-), and incubated with 10 ml prewarmed DMEM containing 100 µg/ml gentamicin. After 60 min, the media was replaced with DMEM containing 10 µg/ml gentamicin and cells were further incubated as before. At 8 h p.i., the cells were washed three times with 10 ml PBS +/- and trypsinized with 2 ml prewarmed trypsin at 37°C. When the cells were detached, prewarmed DMEM was added to quench the trypsin. The cells were transferred to a 50 ml conical tube, pelleted by 5 min centrifugation at 247 x g and washed in 15 ml PBS+/+. Finally, the cells were resuspended in 1 ml / plate fixation reagent (RNA later: RNA protect 350:1), incubated at RT for 20 min, and then stored at 4°C until sorting.

HeLa cell lines stably transfected for expression of LAMP1-GFP or Lifeact-GFP have been described before [93] and were maintained and infected as described above.

## Cell sorting and RNA extraction

RNA was extracted from three independent infection experiments of HeLa cells with STM strain ATCC 14028S or SPA 45157. As controls, RNA was also extracted from three independent cultures of uninfected cells and from bacterial cultures grown in LB to stationary phase under microaerobic conditions. HeLa cells (infected or uninfected) were fixed, centrifuged at 247 × g for 5 min and resuspended in 3 ml PBS+/+. The cells were filtered through 30 µm pre-separation filter into 5 ml FACS tubes (SARSTEDT) and sorted on a BD-FACSAria IIu (100 µm nozzle; ND filter 2) at a rate of approximately 300 cells/sec. 2–2.5×10⁶ GFP-positive (infected) HeLa cells were collected and centrifuged at 247 × g for 5 min. RNA extraction was performed using Hybrid-R RNA Purification System (GeneAll). RNA samples were treated with RNase-free DNase (Qiagen) according to manufacturer's protocol, followed by ethanol precipitation and Ribosomal RNA depletion that was performed by Epidemiology Ribo-Zero Gold rRNA Removal Kit (Illumina).

## RNA sequencing

RNAseq was performed at the GeT-PlaGe core facility, INRAE Toulouse. RNA-Seq libraries were prepared using the Illumina TruSeq Stranded mRNA Library Prep, without a poly-A mRNA selection, according to the manufacturer's protocol. Briefly, RNAs were fragmented to generate double-stranded cDNA and adaptors were ligated. 11 cycles of PCR were applied to amplify the libraries and their quality was assessed using a Fragment Analyser. Libraries were quantified by qPCR using the Kapa Library Quantification Kit (Roche). To reach 5 million bacterial reads for each sample, as recommended in [94], RNA-Seq reads were generated from three lanes of an Illumina HiSeq3000 and two lanes of an Illumina NOVAseq6000 using a paired-end read length of 2×150 bp. Ribosomal RNA (rRNA) reads were filtered out using sortmerna-2.1b [95] against Silva and Rfam databases. The raw reads data for this study have been deposited in the European Nucleotide Archive (ENA) at EMBL-EBI under accession number PRJEB46495 and as specified in S7 Table.

## Genome annotation and quality control

A flowchart summarizing the RNA-seq pipeline and bioinformatic analyses is shown in S3 Fig. Annotation of SPA 45157 and STM 14028S genomes was performed using EugenePP 1.2 [96] with general and targeted evidence to predict coding genes and noncoding features (S1 Fig).
Functional prediction was produced with InterProscan v5.15–54.0 with GO terms, COG version COG2014 [97], REBASE v708 [98], IslandViewer4 [99] and PHAST [100].

Structural and regulatory non-coding RNA were annotated using RNAspace software v1.2.1 [101] together with the following tools: RNAmmer 1.2 for rRNA [102], tRNAScan-SE 1.23 for tRNA [103] and blastn 2.2.25 against RFAM 10.0 for other ncRNA [104]. The gff3 files produced by Eugene and RNA annotation were merged by a semi-manually procedure using BEDTools-2.26.0 [105].

FastQC v0.11.7 was used for quality control of the raw sequence data. Cutadapt 1.8.3 was used to remove Illumina Truseq adaptors and to filter out all reads with length shorter than 36 nucleotides.

## Reads mapping and counting

The remaining, non rRNA reads were mapped with STAR-2.6.0c [106], with default parameters, against a “hybrid” genome composed by the Ensembl GRCh38 human primary assembly (release 92) and the SPA (accession number CP076727) or STM (accession number CP001362 –CP001363) genomes. In addition, in order to analyze bacterial samples, we indexed each bacterial genome separately by decreasing the parameter “genomeSAindexNbases” to 10 for STM and to 9 for SPA, while producing the alignment (.bam) files sorted by coordinates. Samtools flagstat v1.8 was used to obtain mapping statistics. Then we added read groups to trace the origin of each read with AddOrReplaceReadGroups of Picard tools v2.18.2. Next, we merged the bam files that originated from resequencing of the same sample using the Merge-SamFiles command of Picard tools. The featureCounts command from the Subread-1.6.0 tool [107] was used to compute fractional counting for all annotated genes.

## Clustering of orthologous genes

Orthologous genes were computed using Roary v.3.11.2 [108] with default parameters (minimum of 95% identity between coding sequences of the same cluster) except for the “-s” option that avoids paralog splitting in different clusters. Based on the orthologous groups, we produced raw count matrices that compare the expression level between orthologous genes of SPA and STM. We also added groups including paralog genes on either SPA or STM genome into the raw count matrix to sum up the raw reads count of each putative paralog gene.

## Differential gene expression analysis

Normalization and differential analysis were carried out according to the DESeq2 model and package [109] using adjusted p-value with FDR correction and the Wald test to infer the probability value. Genes that had too few reads (the sum of reads for all the samples was less than 10) were discarded. To compare DEGs between SPA and STM, we created a count matrix composed of the shared protein-coding and ncRNA genes between these two serovars. The heatmaps were created by applying DESeq2 variance stabilizing transformations (the vst command) on the counts followed by subtracting the mean of the gene across all samples. The heatmaps were drawn using the ComplexHeatmap R package [110].

## Reverse transcription and quantitative real-time PCR (qRT-PCR)

RNA was extracted from infected and FACS-sorted HeLa cells, as detailed above. 500 ng of RNA were reverse-transcribed using qScript cDNA synthesis kit (Quanta-bio). Real-time PCR was performed as previously described [23] on a StepOnePlus Real-Time PCR System
(Applied Biosystems). Relative quantity of transcripts was calculated as 2^(-ΔΔCt), using the rpoD gene as an endogenous normalization control.

Western blotting

HeLa cells were infected with *Salmonella* strains carrying pBR-GFP2 plasmid and pWSK29 carrying HA-tagged alleles of *sopE2* or *sopB*. The infection procedure was performed as above with *Salmonella* cultures that were grown overnight under microaerobic conditions. Cell fixation was performed by incubation with 5 ml of 4% formaldehyde for 15 min under gentle shaking. Five ml of PBS +/- were added and the cells were centrifuged for 5 min at 2,700 × g. Subsequently, the cells were washed three times with PBS +/- and finally resuspended in PBS +/- containing 1% BSA. The cells were stored at 4°C until sorting. 1×10⁶ GFP-positive cells were centrifuged for 5 min at 247 × g, and resuspended in 100 µl sample buffer. Western blot of the infected cells was performed as previously described [111] using the anti-2HA (Abcam ab18181; 1:1,000), anti-FliC (Abcam ab93713; 1:200,000) and anti-DnaK (Abcam ab69617; 1:10,000) antibodies.

Host cell infection and live-cell imaging

HeLa cells stably transfected with LAMP1-GFP were seeded in surface-treated 8-well slides (ibidi) 24 or 48 h prior to infection to reach ~80% confluency (~80,000 cells) and were then used for infection. STM strains for infection experiments were subcultured from an overnight culture (1:31) in fresh LB medium and grown for 3.5 h at 37°C. For infection, SPA strains were grown for 8 h under aerobic conditions, subcultured (1:100) in fresh LB medium and stationary-phase subcultures were grown for 16 h under microaerophilic conditions as described in [24]. Bacteria were adjusted to an optical density of 0.2 at 600 nm in PBS and used for infection with the respective MOI (between 25 and 100, dependent on strain). Bacteria were centrifuged onto the cells for 5 min at 500 × g to synchronize infection and the infection was allowed to proceed for 25 min. After washing thrice with PBS, cells were incubated in medium containing 100 µg × ml⁻¹ of gentamicin to kill extracellular bacteria. Afterwards, the cells were maintained in medium supplemented with 10 µg × ml⁻¹ gentamicin. Live-cell imaging (LCI) was performed with Cell Observer microscope (Zeiss) equipped with a Yokogawa Spinning Disc Unit CSU-X1a5000, an incubation chamber, 63× objective (α-Plan-Apochromat, NA 1.4), two ORCA Flash 4.0 V3 cameras (Hamamatsu) and appropriate filters for the respective fluorescence proteins. For overview sequences, HeLa LAMP1-GFP cells in 8-well slides were infected at MOI 30 as described above. LCI was performed at 6 to 8 h p.i. as described above using a 40x objective (Plan-Apochromat, NA 1.4). A tile area of 641 x 664 x 2 µm was imaged for 180 s and stitched.

Caco-2 cells were seeded in 8-well slides 48 h prior to infection to reach ~80% confluency (ca. 70,000 cells) and then infected. If indicated, cells were pulse-chased with Dextran Alexa Fluor 488/568/647 (Thermo Fisher Scientific) 16 h prior infection. Infection was conducted as described before at MOI 50. Imaging was performed from 3 to 5 h p.i. as described above. Infected Caco 2 cells in 8-well slides were fixed and used for flagella staining as described below.

Human monocyte cell line U937 cells were seeded in 8-well slides 72 h prior infection to reach ~80% confluency (ca. 200,000 cells) and treated with 50 ng x µl-1 phorbol 12 myristate 13-acetate (PMA) for differentiation and cell attachment until medium was exchanged after 48 h. Cells were pulse-chased with Dextran Alexa Fluor and infected with overnight cultures at MOI 30 as described above. LCI was performed from 3 to 5 h p.i. as described above.
## Immunostaining of infected cells

HeLa cells stably transfected with LAMP1-GFP were seeded in surface-treated 24-well plates (TPP) on coverslips 24 h or 48 h prior infection to reach ~80% confluency (~180,000 cells) and were used for infection. Cultivation of bacteria and infection of HeLa cells was carried out as described above. After the desired incubation time the cells were washed thrice with PBS and fixed with 3% PFA in PBS for 15 min at RT, followed by immunostaining. Fixed cells were washed thrice with PBS and incubated in blocking solution (2% goat serum, 2% bovine serum albumin, 0.1% saponin, 0.2% Triton X-100) for 30 min. Cells were stained 1 h at RT with primary antibodies against *Salmonella* (H:i, BD Difco 228241; H:a, Sifin TR1401; O:5, Sifin TR5303). After washing thrice with PBS, cells were incubated with the appropriate secondary antibodies for 1 h at RT. Coverslips were mounted with Fluoroshield (Sigma) and sealed with Entellan (Merck). Microscopy of fixed samples was either performed using a Cell Observer Spinning Disc microscope (Zeiss) as described above, or using a Leica SP5 confocal laser-scanning microscope using 100× objective (HCX PL APO CS, NA 1.4–0.7) and polychroic mirror TD 488/543/633.

## Reinvusion assays

HeLa cells were seeded in surface-treated 6-well plates (TPP) 24 h or 48 h prior to infection to reach ~100% confluency (1×10⁶ cells) on the day of infection. Infection cultures of STM and SPA strains were grown as either 3.5 h aerobic or 16 h microaerobic subcultures and used for infection as described above. At 1 h p.i., cells were washed thrice with PBS and lysed using 0.1% Triton X-100. CFU were determined by plating serial dilutions of lysates and inoculum on Mueller-Hinton II (MH) agar and incubated overnight at 37°C. The percentage of internalized bacteria of the inoculum (from subcultures) was calculated. 8 h p.i., cells were washed twice with PBS, once with ice-cold ddH₂O, and lysed using ddH₂O. Bacteria liberated from host cells by hypotonic lysis were used to infect naïve HeLa cells, and 1 h p.i. cells were washed and lysed and CFU were determined as described above. The percentage of internalized bacteria of the inoculum (from liberated bacteria) was calculated.

## Supporting information

**S1 Fig. Expression of fimbriae genes in intracellular STM vs. SPA.** (A) A Heat map of RNA-Seq results showing the relative transcription of genes encoding curli (Csg), Fim, or Bcf fimbriae in three independent HeLa cells infections with STM and SPA. (B) The fold change in the expression of *csgA* and *csgD* in intracellular SPA relative to their expression in intracellular STM was analyzed by qRT-PCR. RNA was extracted from FACS-sorted intracellular salmonellae at 8 h p.i. The results show the means of 4–8 independent reactions and the error bars indicate the SEM.
(TIF)

**S2 Fig. Intracellular replication of SPA and STM in HeLa cells.** Intracellular growth of STM vs. SPA in HeLa cells at 2, 4, 6, 8 and 24 h p.i. Bacterial replication was determined by the genetic protection assay. Replication was calculated as the ratio between the intracellular bacteria (CFU) recovered at each time point and the number of CFU at 2 p.i. The means and the SEM of four independent infections are shown.
(TIF)

**S3 Fig. A flowchart of RNA-Seq and bioinformatic analyses.** The different steps and main tools used for the bioinformatic analyses illustrated. RNA extracted from STM and SPA grown in culture and from infected cells was subjected to NGS on Illumina sequencing. RNA-Seq reads were cleaned and quality controlled using FastQC. Cutadapt was used to remove
Illumina Truseq adaptors and to filter out all reads shorter than 36 nucleotides. SortMeRNA was used to filter out ribosomal RNA reads from raw data. Human and bacterial genomes were merged into "Hybrid genomes" and indexed with STAR. Bacterial genomes were reannotated using EugenePP. STAR was used to align reads to reference genomes and sort resulting bam files by coordinates. To compute fractional counting for all annotated coding and non-coding genes the Subread tool with parameters "-s 2 (reversely stranded), -p (count fragments instead of reads) -M, -O,—fraction" was used. Clustering of *Salmonella* homologous genes was computed using Roary. Normalization and differential expressed genes analysis were carried out according to the DESeq2 model and package using adjusted P-value and the Wald test to infer the probability value. To obtain DEGs between STM and SPA a count matrix composed of the shared and coding genes and ncRNAs between these two serovars was created.
(TIF)

**S1 Table. The transcriptome of STM grown in LB and intracellularly.** The RNA-Seq reads obtained from STM grown on LB or extracted from HeLa cells infected with STM were mapped to the STM 12048S genome. This table shows the normalized (TPM) featurecount counts of reads assigned to the STM genes. This table also lists the genes that were found expressed in LB only, intracellularly only, and at both conditions.
(XLSX)

**S2 Table. DEGs of intracellular STM vs. LB.**
(XLSX)

**S3 Table. The transcriptome of SPA grown in LB or intracellular.** The RNA-Seq reads obtained from SPA grown on LB or extracted from HeLa cells infected with SPA were mapped to the SPA 45157 genome. This table shows the normalized (TPM) featurecount counts of reads assigned to the SPA genes. This table also lists the genes that were found expressed in LB only, intracellular only, and under both conditions.
(XLSX)

**S4 Table. DEGs of intracellular SPA vs. LB.**
(XLSX)

**S5 Table. DEGs between intracellular STM and intracellular SPA.**
(XLSX)

**S6 Table. Bacterial strains and primers utilized in this study.**
(DOCX)

**S7 Table. RNA-Seq metadata and accession numbers of the 12 samples analyzed in this study.**
(XLSX)

**S1 Movie. Live-cell imaging of STM and SPA strains reveals intracellular motility of SPA.**
To investigate intracellular motility, HeLa cells stably expressing LAMP1-GFP (green) or Life-Act-GFP (green) as indicated were infected with SPA WT, SPA *ssaR*, SPA *fliC*, or STM ATCC 14028s WT, each harboring pFPV-mCherry for constitutive expression of mCherry (red). Individual clips of about 30 seconds are shown successively. Time stamps in upper left corner indicate sec.ms. Scale bars: 10 µm.
(MP4)

**S2 Movie. Real time imaging of intracellular motility of SPA WT.** HeLa cells stably expressing LifeAct-GFP (green) were infected with SPA WT expression mCherry (red) and imaged at
a frame rate of ~40 frames per second. Time stamps in upper left corner indicate sec.ms. Scale bar: 10 µm.
(MP4)

**S3 Movie. Real time imaging of intracellular motility of SPA WT.** HeLa cells stably expressing LifeAct-GFP (green) were infected with SPA WT expression mCherry (red) and imaged at a frame rate of ~40 frames per second. Time stamps in upper left corner indicate sec.ms. Scale bar: 10 µm.
(MP4)

**S4 Movie. Live-cell imaging of further intracellular STM and SPA strains.** HeLa cells stably expressing LAMP1-GFP (green) were infected with STM SL1344 WT, or 7 isolates of SPA as indicated in the lower left corner of movie sections. All strains harbored pFPV-mCherry for constitutive expression of mCherry (red) Time stamps in upper left corner indicate sec.ms.
Scale bars: 10 µm.
(MP4)

**S5 Movie. Live-cell imaging for scoring of intracellular phenotypes of SPA.** HeLa cells stably expressing LAMP1-GFP (green) were infected with SPA WT harboring pFPV-mCherry for constitutive expression of mCherry (red). Live-cell imaging was performed by SD microscopy and infected host cell in several field of view were scored for SPA residing within LAMP1-positive compartments (yellow arrows), hyper-replicate in host cell cytosol (violet arrows), or showing motility in host cell cytosol (red arrows). A representative analyses session for data for SPA WT shown in Fig 8 was recorded, and overviews and details of 3 fields of view with 4 x 3 tile grid with 5 slices each (641 x 664 x µm total), 15 images in 180 seconds of view were stitched.
(MP4)

**S6 Movie. Live-cell imaging for scoring of intracellular phenotypes of STM.** HeLa cells stably expressing LAMP1-GFP (green) were infected with STM WT harboring pFPV-mCherry for constitutive expression of mCherry (red). Live-cell imaging was performed by SD microscopy and infected host cells in several field of view were scored for SPA residing within LAMP1-positive compartments (yellow arrows), hyper-replicate in host cell cytosol (violet arrows), or showing motility in host cell cytosol (red arrows). A representative analyses session for data for STM WT shown in Fig 9 was recorded, and overviews and details of 3 fields of view were stitched as described or Movie 6.
(MP4)

**S7 Movie. Live-cell imaging of intracellular motility of SPA WT in Caco-2 cells.** Caco-2 cells were infected with SPA WT harboring pFPV-mCherry for constitutive expression of mCherry (red). Caco-2 cells were pulsed-chased with Dextran Alexa-488 (green) prior infection to label endosomal compartments. Time stamps in upper left corner indicate sec.ms, Scale bar: 10 µm.
(MP4)

**S8 Movie. Live-cell imaging of intracellular motility of SPA WT in U937 cells.** Human macrophages cell line U937 was pulse-chased with Dextran Alexa-647 (blue) prior infection to label endosomal compartments. Cells were infected with SPA WT harboring pFPV-mCherry for constitutive expression of mCherry (red). Time stamps in upper left corner indicate sec.ms, Scale bar: 10 µm.
(MP4)
**Acknowledgments**

We thank C. Gaspin for providing annotation of ncRNA genes of the STM and SPA genomes.

**Author Contributions**

**Conceptualization:** Claire Hoede, Felix Scharte, Thomas Schiex, Isabelle Virlogeux-Payant, Guntram A. Grassl, Michael Hensel, Ohad Gal-Mor.

**Data curation:** Helit Cohen, Claire Hoede, Felix Scharte, Charles Coluzzi, Emiliano Cohen, Ludovic Mallet, Hélène Chiapello.

**Formal analysis:** Claire Hoede, Felix Scharte, Charles Coluzzi, Emiliano Cohen, Inna Shomer, Ludovic Mallet, Sébastien Holbert, Remy Felix Serre, Michael Hensel, Hélène Chiapello, Ohad Gal-Mor.

**Funding acquisition:** Thomas Schiex, Isabelle Virlogeux-Payant, Guntram A. Grassl, Michael Hensel, Ohad Gal-Mor.

**Investigation:** Helit Cohen, Claire Hoede, Felix Scharte, Charles Coluzzi, Emiliano Cohen, Inna Shomer, Ludovic Mallet, Sébastien Holbert, Isabelle Virlogeux-Payant, Guntram A. Grassl, Michael Hensel, Hélène Chiapello, Ohad Gal-Mor.

**Methodology:** Helit Cohen, Claire Hoede, Felix Scharte, Charles Coluzzi, Inna Shomer, Remy Felix Serre, Guntram A. Grassl, Michael Hensel, Ohad Gal-Mor.

**Project administration:** Thomas Schiex, Isabelle Virlogeux-Payant, Guntram A. Grassl, Michael Hensel, Hélène Chiapello, Ohad Gal-Mor.

**Resources:** Remy Felix Serre, Michael Hensel, Ohad Gal-Mor.

**Software:** Ludovic Mallet, Hélène Chiapello.

**Supervision:** Thomas Schiex, Isabelle Virlogeux-Payant, Guntram A. Grassl, Michael Hensel, Hélène Chiapello, Ohad Gal-Mor.

**Validation:** Helit Cohen, Emiliano Cohen, Sébastien Holbert, Ohad Gal-Mor.

**Visualization:** Felix Scharte, Emiliano Cohen, Michael Hensel, Ohad Gal-Mor.

**Writing – original draft:** Felix Scharte, Michael Hensel, Ohad Gal-Mor.

**Writing – review & editing:** Felix Scharte, Sébastien Holbert, Isabelle Virlogeux-Payant, Guntram A. Grassl, Michael Hensel, Ohad Gal-Mor.

**References**

1. Gal-Mor O. Persistent Infection and Long-Term Carriage of Typhoidal and Nontyphoidal Salmonellae. Clin Microbiol Rev. 2019; 32(1). https://doi.org/10.1128/CMR.00088-18 PMID: 30487167; PubMed Central PMCID: PMC6302356.

2. Kraus MD, Amatya B, Kimula Y. Histopathology of typhoid enteritis: morphologic and immunophenotypic findings. Mod Pathol. 1999; 12(10):949-55. Epub 1999/10/26. PMID: 10530559.

3. Nguyen QC, Everest P, Tran TK, House D, Murch S, Parry C, et al. A clinical, microbiological, and pathological study of intestinal perforation associated with typhoid fever. Clin Infect Dis. 2004; 39(1):61-7. Epub 2004/06/19. https://doi.org/10.1086/421555 PMID: 15206054.

4. Sprinz H, Gangarosa EJ, Williams M, Hornick RB, Woodward TE. Histopathology of the upper small intestines in typhoid fever. Biopsy study of experimental disease in man. Am J Dig Dis. 1966; 11(8):615-24. Epub 1966/08/01. https://doi.org/10.1007/BF02233509 PMID: 5943767.

5. Crump JA, Luby SP, Mintz ED. The global burden of typhoid fever. Bull World Health Organ. 2004; 82(5):346-53. Epub 2004/08/10. PMID: 15298225; PubMed Central PMCID: PMC2622843.
6. Havelaar AH, Kirk MD, Torgerson PR, Gibb HJ, Hald T, Lake RJ, et al. World Health Organization Global Estimates and Regional Comparisons of the Burden of Foodborne Disease in 2010. PLoS Med. 2015; 12(12):e1001923. https://doi.org/10.1371/journal.pmed.1001923 PMID: 26633896; PubMed Central PMCID: PMC4668832.

7. Meltzer E, Schwartz E. Enteric fever: a travel medicine oriented view. Curr Opin Infect Dis. 2010; 23 (5):432-7. Epub 2010/07/09. https://doi.org/10.1097/QCO.0b013e32833c7ca1 PMID: 20613510.

8. Ochiai RL, Wang X, von Seidlein L, Yang J, Bhutta ZA, Bhattacharya SK, et al. *Salmonella* paratyphi A rates, Asia. Emerg Infect Dis. 2005; 11(11):1764-6. Epub 2005/12/02. https://doi.org/10.3201/ eid1111.050168 PMID: 16318734.

9. Gal-Mor O, Suez J, Elhadad D, Porwollik S, Leshem E, Valinsky L, et al. Molecular and cellular characterization of a *Salmonella enterica* serovar Paratyphi a outbreak strain and the human immune response to infection. Clin Vaccine Immunol. 2012; 19(2):146-56. Epub 2011/12/23. https://doi.org/ 10.1128/CVI.05468-11 PMID: 22190395; PubMed Central PMCID: PMC3272918.

10. Liu SL, Ezaki T, Miura H, Matsui K, Yabuchi E. Intact motility as a *Salmonella typhi* invasion-related factor. Infect Immun. 1988; 56(8):1967-73. Epub 1988/08/01. https://doi.org/10.1128/iai.56.8.1967-1973.1988 PMID: 2840399; PubMed Central PMCID: PMC259509.

11. Abby SS, Rocha EP. The non-flagellar type III secretion system evolved from the bacterial flagellum and diversified into host-cell adapted systems. PLoS Genet. 2012; 8(9):e1002983. Epub 2012/10/03. https://doi.org/10.1371/journal.pgen.1002983 PMID: 23028376; PubMed Central PMCID: PMC3459982.

12. Zhou D, Galan J. *Salmonella* entry into host cells: the work in concert of type III secreted effector proteins. Microbes Infect. 2001; 3(14-15):1293-8. Epub 2002/01/05. https://doi.org/10.1016/s1286-4579 (01)01489-7 PMID: 11755417.

13. Haraga A, Ohlson MB, Miller SL. *Salmonellae* interplay with host cells. Nat Rev Microbiol. 2008; 6 (1):53-66. Epub 2007/11/21. nrmicro1788 [pii] https://doi.org/10.1038/nrmicro1788 PMID: 18026123.

14. Sekirov I, Gill N, Jogova M, Tam N, Robertson M, de Llanos R, et al. *Salmonella* SPI-1-mediated neutrophil recruitment during enteric colitis is associated with reduction and alteration in intestinal microbiota. Gut Microbes. 2010; 1(1):30-41. Epub 2011/02/18. https://doi.org/10.4161/gmic.1.1.10950 PMID: 21327114; PubMed Central PMCID: PMC3035138.

15. Jennings E, Thurston TLM, Holden DW. *Salmonella* SPI-2 Type III Secretion System Effectors: Molecular Mechanisms And Physiological Consequences. Cell Host Microbe. 2017; 22(2):217-31. Epub 2017/08/12. https://doi.org/10.1016/j.chom.2017.07.009 PMID: 28799907.

16. Knodler LA. *Salmonella enterica*: living a double life in epithelial cells. Curr Opin Microbiol. 2015; 23:23-31. Epub 2014/12/03. https://doi.org/10.1016/j.mib.2014.10.010 PMID: 25461569.

17. Raffatellu M, Chessa D, Wilson RP, Dusold R, Rubino S, Baumler AJ. The Vi capsular antigen of *Salmonella enterica* serotype Typhi reduces Toll-like receptor-dependent interleukin-8 expression in the intestinal mucosa. Infect Immun. 2005; 73(6):3367-74. Epub 2005/05/24. https://doi.org/10.1128/IAI. 73.6.3367-3374.2005 PMID: 15908363; PubMed Central PMCID: PMC1111811.

18. Wangdi T, Winter SE, Baumler AJ. Typhoid fever: "you can't hit what you can't see". Gut Microbes. 2012; 3(2):88-92. Epub 2011/12/14. https://doi.org/10.4161/gmic.18602 PMID: 22156762; PubMed Central PMCID: PMC3370952.

19. Winter SE, Winter MG, Thiennimitr P, Gerriets VA, Nuccio SP, Russmann H, et al. The TviA auxiliary protein renders the *Salmonella enterica* serotype Typhi RcsB regulon responsive to changes in osmolarity. Mol Microbiol. 2009; 74(1):175-93. Epub 2009/08/26. https://doi.org/10.1111/j.1365-2958.2009.06859.x PMID: 19703107; PubMed Central PMCID: PMC2763492.

20. Winter SE, Winter MG, Poon V, Keestra AM, Sterzenbach T, Faber F, et al. *Salmonella enterica* Serovar Typhi conceals the invasion-associated type three secretion system from the innate immune system by gene regulation. PLoS Pathog. 2014; 10(7):e1004207. Epub 2014/07/06. https://doi.org/10.1371/journal.ppat.1004207 PMID: 24992093; PubMed Central PMCID: PMC4081808.

21. Hiyoshi H, Wangdi T, Lock G, Saechao C, Raffatellu M, Cobb BA, et al. Mechanisms to Evade the Phagocyte Respiratory Burst Arose by Convergent Evolution in Typhoidal *Salmonella* Serovars. Cell Rep. 2018; 22(7):1787-97. Epub 2018/02/15. https://doi.org/10.1016/j.celrep.2018.01.016 PMID: 29444431; PubMed Central PMCID: PMC5826628.

22. Mylona E, Sanchez-Garrido J, Hoang Thu TN, Dongol S, Karkey A, Baker S, et al. Very long O-antigen chains of *Salmonella* Paratyphi A inhibit inflammasome activation and pyroptotic cell death. Cell Microbiol. 2021; 23(5):e13306. Epub 2020/12/24. https://doi.org/10.1111/cmi.13306 PMID: 33355403; PubMed Central PMCID: PMC8609438.

23. Elhadad D, Desai P, Grassl GA, McClelland M, Rahav G, Gal-Mor O. Differences in Host Cell Invasion and SPI-1 Expression between *Salmonella enterica* serovar Paratyphi A and the Non-Typhoidal
Serovar Typhimurium. Infect Immun. 2016. Epub 2016/02/10. https://doi.org/10.1128/IAI.01461-15
PMID: 26857569.

24. 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. Epub 2015/06/10. https://doi.org/10.1128/IAI.00468-15 PMID: 26056383.

25. Elhadad D, McClelland M, Rahav G, Gal-Mor O. Feverlike Temperature is a Virulence Regulatory Cue Controlling the Motility and Host Cell Entry of Typhoidal *Salmonella*. J Infect Dis. 2015; 212(1):147-56. Epub 2014/12/11. https://doi.org/10.1093/infdis/jiu663 PMID: 25492917.

26. 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. Epub 2001/10/26. https://doi.org/10.1038/35101614 [pii]. PMID: 11677609.

27. Knodler LA, Steele-Mortimer O. Taking possession: biogenesis of the *Salmonella*-containing vacuole. Traffic. 2003; 4(9):587-99. Epub 2003/08/13. https://doi.org/10.1034/j.1600-0854.2003.00118.x
PMID: 12911813.

28. Kroger C, Colgan A, Srikumar S, Handler K, Sivasankaran SK, Hammarlof DL, et al. An infection-relevant transcriptomic compendium for *Salmonella enterica* Serovar Typhimurium. Cell Host Microbe. 2013; 14(6):683-95. Epub 2013/12/18. https://doi.org/10.1016/j.chom.2013.11.010 PMID: 24331466.

29. Hautefort I, Thompson A, Eriksson-Ygberg S, Parker ML, Lucchini S, Danino V, et al. During infection of epithelial cells *Salmonella enterica* serovar Typhimurium undergoes a time-dependent transcriptional adaptation that results in simultaneous expression of three type 3 secretion systems. Cell Microbiol. 2008; 10(4):958-84. Epub 2007/11/23. https://doi.org/10.1111/j.1462-5822.2007.01099.x PMID: 18031307; PubMed Central PMCID: PMC2343689.

30. Liu Y, Yu K, Zhou F, Ding T, Yang Y, Hu M, et al. Quantitative Proteomics Charts the Landscape of *Salmonella* Carbon Metabolism within Host Epithelial Cells. J Proteome Res. 2017; 16(2):788-97. Epub 2017/02/06. https://doi.org/10.1021/acs.jproteome.6b00793 PMID: 28152601.

31. Deiwick J, Nikolaus T, Erdogan S, Hensel M. Environmental regulation of *Salmonella* pathogenicity island 2 gene expression. Mol Microbiol. 1999; 31(6):1759-73. Epub 1999/04/21. https://doi.org/10.1046/j.1365-2958.1999.01312.x PMID: 10209748.

32. Prost LR, Miller SI. The *Salmonellae* PhoQ sensor: mechanisms of detection of phagosome signals. Cell Microbiol. 2008; 10(3):576-82. Epub 2008/01/10. https://doi.org/10.1111/j.1462-5822.2007.01111.x PMID: 18182085.

33. Hensel M, Shea JE, Waterman SR, Mundy R, Nikolaus T, Banks G, et al. Genes encoding putative effector proteins of the type III secretion system of *Salmonella* pathogenicity island 2 are required for bacterial virulence and proliferation in macrophages. Mol Microbiol. 1998; 30(1):163-74. Epub 1998/10/24. https://doi.org/10.1046/j.1365-2958.1998.01047.x PMID: 9786193.

34. Eriksson S, Lucchini S, Thompson A, Rhen M, Hinton JC. Unravelling the biology of macrophage infection by gene expression profiling of intracellular *Salmonella enterica*. Mol Microbiol. 2003; 47(1):103-18. Epub 2002/12/21. https://doi.org/10.1046/j.1365-2958.2003.03313.x PMID: 12492857.

35. Srikumar S, Kroger C, Hebrard M, Colgan A, Owen SV, Sivasankaran SK, et al. RNA-seq Brings New Insights to the Intra-Macrophage Transcriptome of *Salmonella* Typhimurium. PLoS Pathog. 2015; 11(11):e1005262. Epub 2015/11/13. https://doi.org/10.1371/journal.ppat.1005262 PMID: 26561851; PubMed Central PMCID: PMC4643027.

36. Ratledge C. Iron metabolism and infection. Food Nutr Bull. 2007; 28(4 Suppl):S515-23. Epub 2008/02/27. https://doi.org/10.1177/15648265070284S405 PMID: 18297890.

37. Brown SA, Palmer KL, Whiteley M. Revisiting the host as a growth medium. Nat Rev Microbiol. 2008; 6(9):657-66. Epub 2008/08/06. https://doi.org/10.1038/nrmicro1955 PMID: 18679171; PubMed Central PMCID: PMC3115587.

38. Eisenreich W, Dandekar T, Heesemann J, Goebel W. Carbon metabolism of intracellular bacterial pathogens and possible links to virulence. Nat Rev Microbiol. 2010; 8(6):401-12. Epub 2010/05/11. https://doi.org/10.1038/nrmicro2351 PMID: 20453875.

39. Brzoska P, Boos W. Characteristics of a *Ugp*-encoded and *phoB*-dependent glycerophosphoryl diester phosphodiesterase which is physically dependent on the *Ugp* transport system of *Escherichia coli*. J Bacteriol. 1988; 170(9):4125-35. Epub 1988/09/01. https://doi.org/10.1128/jb.170.9.4125-4135. 1988 PMID: 2842304; PubMed Central PMCID: PMC211418.

40. Conlin CA, Tan SL, Hu H, Segar T. The *apeE* gene of *Salmonella enterica* serovar Typhimurium is induced by phosphate limitation and regulated by *phoBR*. J Bacteriol. 2001; 183(5):1784-6. Epub 2001/02/13. https://doi.org/10.1128/JB.183.5.1784-1786.2001 PMID: 11160112; PubMed Central PMCID: PMC95066.
41. Jiang W, Metcalf WW, Lee KS, Wanner BL. Molecular cloning, mapping, and regulation of Pho regulon genes for phosphonate breakdown by the phosphonatase pathway of *Salmonella* typhimurium LT2. J Bacteriol. 1995; 177(22):6411-21. Epub 1995/11/01. https://doi.org/10.1128/jb.177.22.6411-6421. 1995 PMID: 7592415; PubMed Central PMCID: PMC177490.

42. Bowden SD, Rowley G, Hinton JC, Thompson A. Glucose and glycolysis are required for the successful infection of macrophages and mice by *Salmonella enterica* serovar Typhimurium. Infect Immun. 2009; 77(7):3117-26. Epub 2009/04/22. https://doi.org/10.1128/IAI.00093-09 PMID: 19380470; PubMed Central PMCID: PMC2708584.

43. Bowden SD, Hopper-Chidlaw AC, Rice CJ, Ramachandran VK, Kelly DJ, Thompson A. Nutritional and metabolic requirements for the infection of HeLa cells by *Salmonella enterica* serovar Typhimurium. PLoS One. 2014; 9(5):e96266. Epub 2014/05/07. https://doi.org/10.1371/journal.pone.0096266 PMID: 24797930; PubMed Central PMCID: PMC4010460.

44. Lundberg BE, Wolf RE Jr., Dinauer MC, Xu Y, Fang FC. Glucose 6-phosphate dehydrogenase is required for *Salmonella* Typhimurium virulence and resistance to reactive oxygen and nitrogen intermediates. Infect Immun. 1999; 67(1):436-8. Epub 1998/12/24. https://doi.org/10.1128/IAI.67.1.436-438.1999 PMID: 9864251; PubMed Central PMCID: PMC96332.

45. Romero M, Laorden ML, Hernandez J, Serrano JS. Evidence for involvement of catecholamines in the effect of morphine on ventricular automaticity in the rat. J Auton Pharmacol. 1991; 11(2):93-9. Epub 1991/04/01. https://doi.org/10.1111/j.1474-8673.1991.tb00248.x PMID: 2045387.

46. Tchawa Yimga M, Leatham MP, Allen JH, Laux DC, Conway T, Cohen PS. Role of gluconeogenesis and the tricarboxylic acid cycle in the virulence of *Salmonella enterica* serovar Typhimurium in BALB/c mice. Infect Immun. 2006; 74(2):1130-40. Epub 2006/01/24. https://doi.org/10.1128/IAI.74.2.1130-1140.2006 PMID: 16428761; PubMed Central PMCID: PMC1360343.

47. Gotz A, Goebel W. Glucose and glucose 6-phosphate as carbon sources in extra- and intracellular growth of enteroinvasive *Escherichia coli* and *Salmonella enterica*. Microbiology (Reading). 2010; 156 (Pt 4):1176-87. Epub 2010/01/16. https://doi.org/10.1099/mic.0.034744-0 PMID: 20075042.

48. 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:36. Epub 2009/01/23. https://doi.org/10.1186/1471-2164-10-36 PMID: 19159446; PubMed Central PMCID: PMC2658671.

49. Kaval KG, Garsin DA. Ethanolamine Utilization in Bacteria. mBio. 2018; 9(1). Epub 2018/02/22. https://doi.org/10.1128/mBio.00066-18 PMID: 29463652; PubMed Central PMCID: PMC5821096.

50. Kofoid E, Rappleye C, Stojiljkovic I, Roth J. The 17-gene ethanolamine (*eut*) operon of *Salmonella* Typhimurium encodes five homologues of carboxysome shell proteins. J Bacteriol. 1999; 181 (17):5317-29. Epub 1999/08/28. https://doi.org/10.1128/JB.181.17.5317-5329.1999 PMID: 10464203; PubMed Central PMCID: PMC94038.

51. Horstmann JA, Zschieschang E, Truschel T, de Diego J, Lunelli M, Rohde M, et al. Flagellin phase-dependent swimming on epithelial cell surfaces contributes to productive *Salmonella* gut colonisation. Cell Microbiol. 2017; 19(8). Epub 2017/03/16. https://doi.org/10.1111/cmi.12739 PMID: 28295924.

52. Misselwitz B, Barrett N, Kreibich S, Vonaesch P, Andritschke D, Rout S, et al. Near surface swimming of *Salmonella* Typhimurium explains target-site selection and cooperative invasion. PLoS Pathog. 2012; 8(7):e1002810. Epub 2012/08/23. https://doi.org/10.1371/journal.ppat.1002810 PMID: 22911370; PubMed Central PMCID: PMC3406100.

53. 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. Epub 2000/12/06. https://doi.org/10.1128/MMBR.64.4.694-708.2000 PMID: 11104815; PubMed Central PMCID: PMC99010.

54. Reuter T, Scharte F, Franzkoch R, Liss V, Hensel M. Single cell analyses reveal distinct adaptation of typhoidal and non-typhoidal *Salmonella enterica* serovars to intracellular lifestyle. PLoS Pathog. 2021; 17(6):e1009319. Epub 2021/06/19. https://doi.org/10.1371/journal.ppat.1009319 PMID: 34143852; PubMed Central PMCID: PMC8244875.

55. Welch MD, Way M. Arp2/3-mediated actin-based motility: a tail of pathogen abuse. Cell Host Microbe. 2013; 14(3):242-55. Epub 2013/09/17. https://doi.org/10.1016/j.chom.2013.08.011 PMID: 24034611; PubMed Central PMCID: PMC3933244.

56. Pizarro-Cerda J, Cossart P. *Listeria monocytogenes*: cell biology of invasion and intracellular growth. Microbiol Spectr. 2018; 6(6). Epub 2018/12/14. https://doi.org/10.1128/microbiolspec.GPP3-0013-2018 PMID: 30523778.

57. Boyd EF, Wang FS, Beltran P, Plock SA, Nelson K, Selander RK. *Salmonella* reference collection B (SARB): strains of 37 serovars of subspecies I. J Gen Microbiol. 1993; 139 Pt 6:1125-32. Epub 1993/06/01. https://doi.org/10.1099/00221287-139-6-1125 PMID: 8360609.
58. 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. Epub 2009/12/26. https://doi.org/10.1099/mic.0.032896-0 PMID: 20035008; PubMed Central PMCID: PMC2848694.

59. 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. Epub 2014/08/20. https://doi.org/10.3389/fmicb.2014.00391 PMID: 25136336; PubMed Central PMCID: PMC4120697.

60. Westermann AJ, Forstner KU, Amman F, Barquist L, Chao Y, Schulte LN, et al. Dual RNA-seq unveils noncoding RNA functions in host-pathogen interactions. Nature. 2016; 529(7587):496–501. Epub 2016/01/21. https://doi.org/10.1038/nature16547 PMID: 26789254.

61. Hoiseth SK, Stocker BA. Aromatic-dependent *Salmonella typhimurium* are non-virulent and effective as live vaccines. Nature. 1981; 291(5812):238–9. Epub 1981/05/21. https://doi.org/10.1038/291238a0 PMID: 7015147.

62. Fields PI, Swanson RV, Haidaris CG, Heffron F. Mutants of *Salmonella typhimurium* that cannot survive within the macrophage are avirulent. Proc Natl Acad Sci U S A. 1986; 83(14):5189–93. Epub 1986/07/01. https://doi.org/10.1073/pnas.83.14.5189 PMID: 3523484; PubMed Central PMCID: PMC323916.

63. McFarland WC, Stocker BA. Effect of different purine auxotrophic mutations on mouse-virulence of a Vi-positive strain of Salmonella dublin and of two strains of *Salmonella* Typhimurium. Microb Pathog. 1987; 3(2):129–41. Epub 1987/08/01. https://doi.org/10.1016/0882-4010(87)90071-4 PMID: 2849016.

64. O'Callaghan D, Maskell D, Liew FY, Easmon CS, Dougan G. Characterization of aromatic- and purine-dependent *Salmonella* Typhimurium: attention, persistence, and ability to induce protective immunity in BALB/c mice. Infect Immun. 1988; 56(2):419–23. Epub 1988/02/01. https://doi.org/10.1128/iai.56.2.419-423.1988 PMID: 3276625; PubMed Central PMCID: PMC259298.

65. Powers TR, Haeberle AL, Predeus AV, Hammarlof DL, Cundiff JA, Saldana-Ahuactzi Z, et al. Intracellular niche-specific profiling reveals transcriptional adaptations required for the cytosolic lifestyle of *Salmonella enterica*. PLoS Pathog. 2021; 17(8):e1009280. Epub 2021/08/31. https://doi.org/10.1371/journal.ppat.1009280 PMID: 34460873; PubMed Central PMCID: PMC8432900.

66. Liu Y, Zhang Q, Hu M, Yu K, Fu J, Zhou F, et al. Proteomic Analyses of Intracellular *Salmonella enterica* Serovar Typhimurium Reveal Extensive Bacterial Adaptations to Infected Host Epithelial Cells. Infect Immun. 2015; 83(7):2897–906. Epub 2015/05/06. https://doi.org/10.1128/IAI.02882-14 PMID: 25939512; PubMed Central PMCID: PMC4468536.

67. Osman D, Cavet JS. Metal sensing in *Salmonella*: implications for pathogenesis. Adv Microb Physiol. 2011; 58:175–232. Epub 2011/07/05. https://doi.org/10.1016/B978-0-12-381043-4.00005-2 PMID: 21722794.

68. Ammendola S, Pasquali P, Pistoia C, Petrucci P, Petrarca P, Rotilio G, et al. High-affinity Zn²⁺ uptake system ZnuABC is required for bacterial zinc homeostasis in intracellular environments and contributes to the virulence of *Salmonella enterica*. Infect Immun. 2007; 75(12):5867–76. Epub 2007/10/10. https://doi.org/10.1128/IAI.00559-07 PMID: 17923515; PubMed Central PMCID: PMC2168356.

69. Schaible UE, Kaufmann SH. Iron and microbial infection. Nat Rev Microbiol. 2004; 2(12):946–53. Epub 2004/11/20. https://doi.org/10.1038/nrmicro1046 PMID: 15550940.

70. Finlay BB, Chatfield S, Leung KY, Dougan G, Falkow S. Characterization of a *Salmonella choleraesuis* mutant that cannot multiply within epithelial cells. Can J Microbiol. 1991; 37(7):568–72. Epub 1991/07/01. https://doi.org/10.1139/m91-095 PMID: 1913359.

71. Denkel LA, Rhen M, Bange FC. Biotin sulfoxide reductase contributes to oxidative stress tolerance and virulence in *Salmonella enterica* serovar Typhimurium. Microbiology (Reading). 2013; 159(Pt 7):1447–58. Epub 2013/05/10. https://doi.org/10.1099/mic.0.067256-0 PMID: 23657680.

72. Bakovic M, Fullerton MD, Michel V. Metabolic and molecular aspects of ethanolamine phospholipid biosynthesis: the role of CTP:phosphoethanolamine cytidylyltransferase (Pcyt2). Biochem Cell Biol. 2007; 85(3):283–300. Epub 2007/07/07. https://doi.org/10.1139/o07-006 PMID: 17612623.

73. Garsin DA. Ethanolamine utilization in bacterial pathogens: roles and regulation. Nat Rev Microbiol. 2010; 8(4):290–5. Epub 2010/03/18. https://doi.org/10.1038/nrmicro2334 PMID: 20234377; PubMed Central PMCID: PMC2950637.

74. Joseph B, Przybilla K, Stuhler C, Schauer K, Slaghuis J, Fuchs TM, et al. Identification of *Listeria monocytogenes* genes contributing to intracellular replication by expression profiling and mutant screening. J Bacteriol. 2006; 188(2):556–68. Epub 2005/12/31. https://doi.org/10.1128/JB.188.2.556-568.2006 PMID: 16385046; PubMed Central PMCID: PMC1347271.

75. Thiennimitr P, Winter SE, Winter MG, Xavier MN, Tolstikov V, Huseby DL, et al. Intestinal inflammation allows *Salmonella* to use ethanolamine to compete with the microbiota. Proc Natl Acad Sci U S A.
2011; 108(42):17480–5. Epub 2011/10/05. https://doi.org/10.1073/pnas.1107857108 PMID: 21969563; PubMed Central PMCID: PMC3198331.

76. Singh B, Arya G, Kundu N, Sangwan A, Nongthombam S, Chaba R. Molecular and Functional Insights into the Regulation of d-Galactonate Metabolism by the Transcriptional Regulator DgoR in *Escherichia coli*. J Bacteriol. 2019; 201(4). Epub 2018/11/21. https://doi.org/10.1128/JB.00281-18 PMID: 30455279; PubMed Central PMCID: PMC6351737.

77. Arunima A, Yelamanchi SD, Padhi C, Jaiswal S, Ryan D, Gupta B, et al. "Omics" of Food-Borne Gas- troenteritis: Global Proteomic and Mutagenic Analysis of *Salmonella enterica* Serovar Enteritidis. OMICS. 2017; 21(10):571–83. Epub 2017/10/20. https://doi.org/10.1089/omi.2017.0112 PMID: 29049011.

78. Utley M, Franklin DP, Krogfelt KA, Laux DC, Cohen PS. A *Salmonella* Typhimurium mutant unable to utilize fatty acids and citrate is avirulent and immunogenic in mice. FEMS Microbiol Lett. 1998; 163 (2):129–34. Epub 1998/07/22. https://doi.org/10.1111/j.1574-6968.1998.tb13036.x PMID: 9673014.

79. Nunez-Hernandez C, Tierrez A, Ortega AD, Pucciarelli MG, Godoy M, Eisman B, et al. Genome expression analysis of nonproliferating intracellular *Salmonella enterica* serovar Typhimurium unrav- els an acid pH-dependent PhoP-PhoQ response essential for dormancy. Infect Immun. 2013; 81 (1):154–65. Epub 2012/10/24. https://doi.org/10.1128/IAI.01080-12 PMID: 23090959; PubMed Cen- tral PMCID: PMC3536153.

80. Schmitt CK, Ikeda JS, Darnell SC, Watson PR, Bispham J, Wallis TS, et al. Absence of all components of the flagellar export and synthesis machinery differentially alters virulence of *Salmonella enterica* ser- ovar Typhimurium in models of typhoid fever, survival in macrophages, tissue culture invasiveness, and calf enterocolitis. Infect Immun. 2001; 69(9):5619–25. Epub 2001/08/14. https://doi.org/10.1128/ IAI.69.9.5619-5625.2001 PMID: 11500437; PubMed Central PMCID: PMC98677.

81. Stecher B, Hapfelmeier S, Muller C, Kremer M, Stallmach T, Hardt WD. Flagella and chemotaxis are required for efficient induction of *Salmonella enterica* serovar Typhimurium colitis in streptomycin-pre- treated mice. Infect Immun. 2004; 72(7):4138–50. Epub 2004/06/24. https://doi.org/10.1128/IAI.72.7. 4138-4150.2004 PMID: 15213159; PubMed Central PMCID: PMC427403.

82. Lin D, Rao CV, Slauch JM. The *Salmonella* SPI1 type three secretion system responds to periplasmic disulfide bond status via the flagellar apparatus and the RcsCDB system. J Bacteriol. 2008; 190 (1):87–97. Epub 2007/10/24. https://doi.org/10.1128/JB.01323-07 PMID: 17951383; PubMed Central PMCID: PMC2223759.

83. Boddicker JD, Jones BD. Lon protease activity causes down-regulation of *Salmonella* pathogenicity island 1 invasion gene expression after infection of epithelial cells. Infect Immun. 2004; 72(4):2002–13. Epub 2004/03/25. https://doi.org/10.1128/IAI.72.4.2002-2013.2004 PMID: 15039320; PubMed Central PMCID: PMC375200.

84. Nunez-Hernandez C, Alonso A, Pucciarelli MG, Casadesus J, Garcia-del Portillo F. Dormant intracel- lular *Salmonella enterica* serovar Typhimurium discriminates among *Salmonella* pathogenicity island 2 effectors to persist inside fibroblasts. Infect Immun. 2014; 82(1):221–32. Epub 2013/10/23. https:// doi.org/10.1128/IAI.01304-13 PMID: 24144726; PubMed Central PMCID: PMC3911833.

85. Fattinger SA, Sellin ME, Hardt WD. Epithelial inflammasomes in the defense against *Salmonella* gut infection. Curr Opin Microbiol. 2021; 59:86–94. Epub 2020/11/01. https://doi.org/10.1016/j.mib.2020. 09.014 PMID: 33128958.

86. Kuehl CJ, Dragoi AM, Talman A, Agaisse H. Bacterial spread from cell to cell: beyond actin-based motility. Trends Microbiol. 2015; 23(9):558–66. Epub 2015/05/30. https://doi.org/10.1016/j.tim.2015. 04.010 PMID: 26021574; PubMed Central PMCID: PMC4560970.

87. Knodler LA, Vallance BA, Celli J, Winfree S, Hansen B, Montero M, et al. Dissemination of invasive *Salmonella* via bacterial-induced extrusion of mucosal epithelia. Proc Natl Acad Sci U S A. 2010; 107 (41):17733–8. https://doi.org/10.1073/pnas.1006098107 PMID: 20876119; PubMed Central PMCID: PMC2955089.

88. Knodler LA, Nair V, Steele-Mortimer O. Quantitative assessment of cytosolic *Salmonella* in epithelial cells. PLoS One. 2014; 9(1):e84681. Epub 2014/01/09. https://doi.org/10.1371/journal.pone.0084681 PMID: 24400108; PubMed Central PMCID: PMC3882239.

89. Chong A, Cooper KG, Kari L, Nilsson OR, Hillman C, Fleming BA, et al. Cytosolic replication in epithe- lial cells fuels intestinal expansion and chronic fecal shedding of *Salmonella* Typhimurium. Cell Host Microbe. 2021; 29(7):1177–85 e6. Epub 2021/05/28. https://doi.org/10.1016/j.chom.2021.04.017 PMID: 34043959; PubMed Central PMCID: PMC8282707.

90. Dutta D, Clevers H. Organoid culture systems to study host-pathogen interactions. Curr Opin Immu- nol. 2017; 48:15–22. Epub 2017/08/02. https://doi.org/10.1016/j.coi.2017.07.012 PMID: 28756233; PubMed Central PMCID: PMC7126332.
91. Sepe LP, Hartl K, Iftekhar A, Berger H, Kumar N, Goosmann C, et al. Genotoxic Effect of *Salmonella* Paratyphi A Infection on Human Primary Gallbladder Cells. mBio. 2020; 11(5). Epub 2020/09/24.
https://doi.org/10.1128/mBio.01911-20 PMID: 32963006; PubMed Central PMCID: PMC7512552.

92. Jarvik T, Smillie C, Groisman EA, Ochman H. Short-term signatures of evolutionary change in the *Salmonella enterica* serovar Typhimurium 14028 genome. J Bacteriol. 2010; 192(2):560-7. Epub 2009/11/10.
https://doi.org/10.1128/JB.01233-09 PMID: 19897643; PubMed Central PMCID: PMC2805332.

93. Lorkowski M, Felipe-Lopez A, Danzer CA, Hansmeier N, Hensel M. *Salmonella enterica* invasion of polarized epithelial cells is a highly cooperative effort. Infect Immun. 2014; 82(6):2657-67. Epub 2014/04/09.
https://doi.org/10.1128/IAI.00023-14 PMID: 24711567; PubMed Central PMCID: PMC4019164.

94. Westermann AJ, Gorski SA, Vogel J. Dual RNA-seq of pathogen and host. Nat Rev Microbiol. 2012; 10(9):618-30. Epub 2012/08/15.
https://doi.org/10.1038/nrmicro2852 PMID: 22890146.

95. Kopylova E, Noe L, Touzet H. SortMeRNA: fast and accurate filtering of ribosomal RNAs in metatranscriptomic data. Bioinformatics. 2012; 28(24):3211-7. Epub 2012/10/17.
https://doi.org/10.1093/bioinformatics/bts611 PMID: 23071270.

96. Sallet E, Gouzy J, Schiex T. EuGene-PP: a next-generation automated annotation pipeline for polykaryotic genomes. Bioinformatics. 2014; 30(18):2659-61.
https://doi.org/10.1093/bioinformatics/btu366 PMID: 24880686.

97. Galperin MY, Makarova KS, Wolf YI, Koonin EV. Expanded microbial genome coverage and improved protein family annotation in the COG database. Nucleic Acids Res. 2015; 43 part 6. Epub 2015/09/09.
https://doi.org/10.1093/nar/gku1223 PMID: 25428365; PubMed Central PMCID: PMC4383993.

98. Roberts RJ, Vincze T, Posfaj J, Macelis D. REBASE—a database for DNA restriction and modification: enzymes, genes and genomes. Nucleic Acids Res. 2010; 38 part 10. Epub 2009/10/23.
https://doi.org/10.1093/nar/gkp874 PMID: 19846593; PubMed Central PMCID: PMC2808884.

99. Bertelli C, Laird MR, Williams KP, Simon Fraser University Research Computing G, Lau BY, Hoad G, et al. IslandViewer 4: expanded prediction of genomic islands for larger-scale datasets. Nucleic Acids Res. 2017; 45(W1):W30-W5. Epub 2017/05/05.
https://doi.org/10.1093/nar/gkx343 PMID: 28472413; PubMed Central PMCID: PMC5570257.

100. Zhou Y, Liang Y, Lynch KH, Dennis JJ, Wishart DS. PHAST: a fast phage search tool. Nucleic Acids Res. 2011; 39 (Web Server issue): W347-52. Epub 2011/06/16.
https://doi.org/10.1093/nar/gkr485 PMID: 21672955; PubMed Central PMCID: PMC3125810.

101. Cros MJ, de Monte A, Mariette J, Bardou P, Grenier-Boley B, Gautheret D, et al. RNAspace.org: An integrated environment for the prediction, annotation, and analysis of ncRNA. RNA. 2011; 17 (11):1947-56. Epub 2011/09/29.
https://doi.org/10.1261/rna.2844911 PMID: 21947200; PubMed Central PMCID: PMC3198588.

102. Lagesen K, Hallin P, Rodland EA, Staerfeldt HH, Rognes T, Ussery DW. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007; 35(9):3100-8. Epub 2007/04/25.
https://doi.org/10.1093/nar/gkm160 PMID: 17452365; PubMed Central PMCID: PMC1888812.

103. Chan PP, Lowe TM. tRNAscan-SE: Searching for tRNA Genes in Genomic Sequences. Methods Mol Biol. 2019; 1962:1-14. Epub 2019/04/26.
https://doi.org/10.1007/978-1-4939-9173-0_1 PMID: 31020551; PubMed Central PMCID: PMC6768409.

104. Gardner PP, Daub J, Tate JG, Nawrocki EP, Kolbe DL, Lindgreen S, et al. Rfam: updates to the RNA families database. Nucleic Acids Res. 2009; 37 (Database issue): D136-40. Epub 2008/10/28.
https://doi.org/10.1093/nar/gkn766 PMID: 18953034; PubMed Central PMCID: PMC2686503.

105. Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010; 26(6):841-2. Epub 2010/01/30.
https://doi.org/10.1093/bioinformatics/btq033 PMID: 20110278; PubMed Central PMCID: PMC2832824.

106. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013; 29(1):15-21. Epub 2012/10/30.
https://doi.org/10.1093/bioinformatics/bts635 PMID: 23104886; PubMed Central PMCID: PMC3530905.

107. Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014; 30(7):923-30. Epub 2013/11/15.
https://doi.org/10.1093/bioinformatics/btt656 PMID: 24227677.

108. Page AJ, Cummins CA, Hunt M, Wong VK, Reuter S, Holden MT, et al. Roary: rapid large-scale polykaryote pan genome analysis. Bioinformatics. 2015; 31(22):3691-3. Epub 2015/07/23.
https://doi.org/10.1093/bioinformatics/btv421 PMID: 26198102; PubMed Central PMCID: PMC4817141.
109. Boddu D, George R, Nair S, Bindra M, L GM. Hydroa Vacciniforme-Like Lymphoma: A Case Report From India. J Pediatr Hematol Oncol. 2014. Epub 2014/07/30. https://doi.org/10.1097/MPH. 0000000000000221 PMID: 25072366.

110. Gu Z, Eils R, Schlesner M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics. 2016; 32(18):2847–9. Epub 2016/05/22. https://doi.org/10.1093/ bioinformatics/btw313 PMID: 27207943.

111. Aviv G, Gal-Mor O. Western Blotting Against Tagged Virulence Determinants to Study Bacterial Pathogenicity. Methods Mol Biol. 2018; 1734:47–54. Epub 2017/12/31. https://doi.org/10.1007/978-1- 4939-7604-1_6 PMID: 29288446.