The in vivo transcriptome approach
We used the DNA macroarray technology to profile the transcriptome of
Listeria during mouse infection. We used the previously described
L. monocytogenes whole-genome arrays containing 500-bp-long PCR products specific for each gene
[6]. Ninety-nine per cent of the 2853 predicted ORFs of the
L. monocytogenes EGDe genome are represented on the arrays. They were used to analyze
Listeria transcription profiles under
in vitro growth in BHI in exponential phase at 37°C under aerobic conditions with shaking (pH 7) (
Figure S1), and under
in vivo growth conditions (mouse spleen) at 1, 2 and 3 days post intravenous infection (p.i.).
Listeria present in spleen were analyzed because this organ is with the liver one of the major sites of
L. monocytogenes infection. For unknown reasons, we never succeeded to prepare good quality bacterial RNAs from infected mouse livers. The time points chosen reflect key steps in the
Listeria infectious process.
Culture in BHI in exponential phase at 37°C with shaking was chosen as reference conditions because BHI is the
Listeria reference growth medium where bacteria divide in exponential growth phase at rates that are comparable to those observed for intracellular growth
[17]. In addition, these are the
in vitro reference conditions used in all previous studies analyzing the genome expression of
L. monocytogenes in vitro or intracellularly
[6]–
[15]. However, in order to analyze the potential impact of the
in vitro culture conditions used as reference on the relative gene expression
in vivo, we first analyzed the results obtained comparing transcriptome from
in vivo grown bacteria to transcriptomes from bacteria grown
in vitro in exponential or stationary phase (
Table S2). In addition, expression of known and potential virulence genes was analyzed by quantitative real time-PCR (qPCR) on RNAs extracted from bacteria cultured in BHI at 37°C in exponential or stationary growth phase, or in defined minimal medium
[18], and compared to
in vivo expression (). Results indicated that culture in exponential growth phase are closer conditions to those met by
Listeria in vivo (
Table S2). In addition, even if the expression of tested genes behaved differently in function of the
in vitro conditions, expression of all the genes was always lower
in vitro as compared to
in vivo, independently of the
in vitro growth conditions (). These experiments supported the choice of exponential growth phase in BHI as reference conditions and minimized the impact of the
in vitro growth conditions on the identification of genes differentially expressed
in vivo.
The reliability of the macroarray expression data was further assessed by qPCR. We selected a subset of 10 genes and performed qPCR on cDNA from bacteria grown in either standard medium or extracted from mouse spleens 48 h p.i.. qPCR results and array data exhibited a high correlation coefficient (0.7) (). This strong correlation was also observed for other infection time points (
Figure S2). However, the differences in gene expression, as measured by qPCR, were generally higher, indicating that
in vivo transcriptome data rather underestimate changes in gene expression.
The procedure used for bacterial RNA extraction from infected mouse spleens is an adaptation of the standard procedure originally used for transcriptional analysis of RNA extracted from pure culture. In order to test the effect of the RNA extraction method on gene expression, RNAs from bacteria grown in pure culture were extracted using the two methods. The relative expression of known virulence genes, cold shock genes and potential virulence genes was analyzed by qPCR in the two RNA pools. The results showed that the relative expression of the genes tested is not significantly affected by the RNA extraction procedure ().
For bacteria cultured in BHI at 37°C in exponential phase or extracted from infected mouse spleen at the different times p.i., two different RNA preparations from independent cultures (or infections) were used for cDNA synthesis and subsequent hybridization to two sets of arrays. To identify statistically significant differences in gene expression, we used the Statistical Analysis for Microarrays (SAM) program
[19]. Subsequently, all the genes showing statistically significant changes in the expression level and an at least two-fold change in their level of expression were considered in our analysis.
Important global changes in L. monocytogenes gene expression occur during in vivo growth
Overall, a total of 568 genes representing ≈20% of the total genome exhibited a differential expression during infection as compared to growth in BHI at 37°C in exponential phase. Among these 568 genes, 457 were up regulated (≈80%) and 111 (≈20%) were down regulated during mouse infection as compared to exponential growth in BHI medium (
Table S3).
In order to identify genes potentially implicated in virulence, all the genes differentially regulated
in vivo were analyzed for the presence of an ortholog in the nonpathogenic close relative
Listeria innocua strain CLIP11262
[20]. This analysis revealed that only 30 of the
in vivo regulated genes (25 up and 5 down regulated) were absent from
L. innocua (). Of these 30 genes, 20 were
L. monocytogenes “specific” (i.e. also present in
L. monocytogenes 1/2a F6854,
L. monocytogenes 4b F2365 and H7858
[21], and absent from
L. innocua). Interestingly, of these 20 genes, 16 were up regulated
in vivo. Among these 16 genes, 11 have been previously implicated in
Listeria virulence. The remaining 10
in vivo regulated genes, among which 9 up- and 1 down-regulated
in vivo, appeared lineage specific, i.e. present only in the sequenced serovar 1/2a strains ().
| Table 1L. monocytogenes EGDe genes absent from L. innocua and differentially regulated in the host. |
To identify genes regulated during different stages of listeriosis, gene expression levels of spleen-recovered bacteria at different time points p.i. were compared. This analysis revealed a core regulon of 106 genes (68 up and 38 down regulated) whose expression was significantly differentially regulated at all the time points of the infection as compared to bacteria grown in pure culture (). No gene appeared specifically differentially regulated at 24 h p.i. At two days p.i., a large proportion (245/457) of genes was up regulated. The largest number of down regulated genes was observed 72 h p.i..
Major virulence regulators and their downstream target genes are highly up regulated in vivo
As compared to
Listeria grown in BHI at exponential phase, bacteria extracted from mouse spleens showed a differential expression of genes belonging to various functional categories (). In particular, analysis of the expression profile of the 50 genes previously implicated in
Listeria virulence in the mouse model revealed that 29 were up regulated during infection, and two (
stp and
fbpA) down regulated
in vivo ( and
Table S1).
We observed that the entire virulence gene cluster of
L. monocytogenes comprising the genes
prfA,
plcA,
hly,
mpl,
actA and
plcB was highly activated during the 3 first days of infection (). In addition to the virulence gene cluster, genes encoding the two major
L. monocytogenes factors implicated in entry into eukaryotic cells (
inlA and
inlB)
[22], and
uhpT, a gene encoding a sugar phosphate transporter that mediates rapid intracellular proliferation
[23] were also activated during infection. PrfA is the principal regulator of the expression of not only these key virulence genes, but also of most other
L. monocytogenes genes involved in intracellular survival and virulence
[6]. The 12 genes previously reported to be preceded by a PrfA box and positively regulated by PrfA in a transcriptional analysis of the PrfA regulon
[6], were all highly up regulated in mouse spleens (). From the 53 other genes already shown as positively regulated by PrfA
[6], 20 were more expressed
in vivo. As previously shown
[8], 19 of these 20 genes are also controlled by SigB, including the LPXTG internalin-like protein
inlH known to be involved in
Listeria virulence
[24]. Two genes,
lmo0206 and
lmo0207, recently shown as regulated by PrfA and implicated in
L. monocytogenes intracellular survival
[14] were also activated in infected mice. Importantly, no gene previously shown under the PrfA positive regulation appeared down regulated during mouse infection.
| Table 2L. monocytogenes EGDe genes positively controlled by PrfA and up regulated in the host. |
VirR, another key
Listeria virulence regulator that mainly controls genes involved in the modification of bacterial surface components, is the response regulator of a two-component system (TCS) implicated in cell invasion and virulence
[13]. Using a transcriptomic approach, 17 genes were previously identified as regulated by VirR
in vitro [13]. In our study, 13 of the 17 VirR regulated genes, including the
dlt operon and
mprF, were up regulated
in vivo (). The
dlt operon is necessary for D-alanylation of lipoteichoic acid (LTA) and was reported to be important for
L. monocytogenes virulence
[13]. The VirR regulated
mprF encodes a protein shown to be required for lysinylation of phospholipids in the
Listeria cytoplasmic membrane and to confer
Listeria resistance to cationic antimicrobial peptides (CAMPs)
[25]. The
virR and
virS genes were themselves up regulated, constituting the only TCS whose expression of both components was induced in mouse spleens.
| Table 3L. monocytogenes EGDe genes regulated by VirR and up regulated in the host. |
In addition to VirRS, the
L. monocytogenes genome contains 15 additional predicted TCS systems
[26]. Genes encoding one component of three of these TCS (
degU,
resD and
phoR) were also up regulated
in vivo. DegU is an orphan response regulator (absence of the sensor kinase DegS in the
L. monocytogenes genome) and a pleiotropic regulatory system previously involved in
Listeria virulence
[27],
[28]. In particular, DegU has been implicated in the regulation of some
Listeria secreted proteins (
gap,
tsf,
sod,
lmo0644)
[26]. Interestingly, the expression of these four genes was also increased in mouse spleens.
Finally, OhrR a transcriptional regulator controlling OhrA, a hydroxyperoxidase implicated in intracellular survival of
Listeria [14], as well as several predicted transcriptional regulators were up regulated in infected mouse spleens.
Strong activation of genes encoding cell wall metabolism proteins during infection
In addition to genes already mentioned and involved in LTA modification (dltABCD), we observed that several genes implicated in peptidoglycan (PG) biosynthesis (lmo0516, lmo0540, lmo1438, lmo1521, lmo1855, lmo2522, lmo2526 and pbpB), cell shape determination (mreBC, lmo1713), cell wall peptide synthesis (murC) were up regulated in bacteria growing in mouse spleens ().
| Table 4L. monocytogenes EGDe genes implicated in cell wall metabolism and differentially regulated in the host. |
The expression of 3 genes encoding virulence factors involved in bacterial cell wall modifications (
murA,
iap, and
pgdA)
[29]–
[31] was also increased
in vivo. MurA and P60, the
iap gene product, are two SecA2-secreted autolysins required for
Listeria full virulence
[29],
[30].
pgdA encodes for the PG N-deacetylase of
L. monocytogenes that was demonstrated as playing an important role in virulence and evasion from host defenses
[31]. In addition,
spl [32] and
lmo2203 are two other autolysins encoding genes up regulated
in vivo, but until now never implicated in virulence.
Moreover,
prsA2, a gene encoding a surface protein involved in protein folding and previously shown as implicated in
Listeria intracellular survival and virulence
[14],
[33] was up regulated
in vivo. Interestingly, the gene encoding the sortase SrtB that covalently links proteins to the
Listeria peptidoglycan, and two genes encoding SrtB substrates (
svpA and
lmo2186)
[34], were also over expressed
in vivo ().
Differential expression of genes encoding specific surface and secreted proteins during infection
Whereas a total of 44 genes encoding potential surface proteins were up regulated
in vivo, only three were observed as down regulated during infection (
lspA,
lmo1851 and
lmo2642) (). In addition, among the 55 proteins previously identified in the cell wall subproteome of
L. monocytogenes [35], we found that 23 were up regulated
in vivo (
Table S4). The
L. monocytogenes genome encodes 41 LPXTG surface proteins
[20],
[36],
[37]. This class includes proteins containing leucine rich repeats (LRRs) and belonging to the internalin family. Four LPXTG-protein encoding genes were up regulated
in vivo. In addition to InlA and InlH,
lmo1290 and lmo2714 are the two other LPXTG encoding genes activated during infection (). Four genes encoding proteins associated to the cell wall via GW modules were also more expressed
in vivo:
inlB, the known invasion protein
[38], and
lmo1521,
lmo2203 and lmo2713.
actA [39] was the only gene encoding a protein with a carboxyl-terminal hydrophobic tail up regulated
in vivo. Genes encoding lipoproteins previously implicated in
Listeria virulence (TcsA and OppA)
[33],
[40] or in cell invasion (LpeA)
[41], were over expressed in mouse spleens. In addition, 10 genes predicted to encode other lipoproteins were activated
in vivo ().
| Table 5L. monocytogenes EGDe cell surface encoding genes differentially regulated in the host. |
Protein secretion is of key importance in both the colonization process and virulence of
Listeria [42]. Besides
L. monocytogenes virulence factors with a signal peptide (ActA, LLO, InlA, InlB, InlC, InlH, Mpl, MurA, PlcA, PlcB, P60 and SvpA), three other virulence proteins (Fri, TcsA and Sod) were also found secreted in the
Listeria culture supernatant
[43]. All the genes encoding these secreted virulence factors appeared activated in our
in vivo approach (
Table S5). The analysis of the products present in the
Listeria culture supernatant after growth
in vitro allowed the identification of 89 additional proteins
[43]. 29 of the genes encoding these secreted proteins were up regulated
in vivo (
Table S5). Most of the
Listeria secreted proteins are presumed to be secreted through the Sec translocation system. A gene encoding one component of the predicted Sec system,
secE, was observed up regulated
in vivo. SecA2 is an auxiliary secretory protein required for persistent colonization of host tissues, and responsible for the secretion of several
Listeria virulence factors (MurA, P60, Sod, OppA and TcsA)
[29],
[30],
[44]. We observed an
in vivo up regulation of the majority of the genes encoding SecA2-secreted proteins, including all the SecA2-secreted virulence factors (
Table S5).
In vivo high expression of genes involved in DNA metabolism, RNA and protein synthesis, cell division and multiplication
We observed an
in vivo up regulation of several genes involved in DNA synthesis (
dnaX and
lmo0162), DNA restriction/modifications and repair (
mutL,
uvrB,
lmo1639 and
lmo1782), DNA recombination (
recFRX,
codV and
lmo2267), and DNA packaging and segregation (
gyrA,
hup,
lmo1606 and
lmo2794) (
Table S6). In addition, the expression of 25 genes encoding ribosomal proteins, as well as genes involved in protein synthesis initiation (
infAC), elongation (
fus,
tsf,
lmo1067) and termination (
frr) was up regulated during infection. Genes encoding proteins implicated in chromosomal replication and segregation (
dnaABC,
ssb and
divIVA), and cell elongation and division (
mreBC,
ftsHX and
lmo0196) were also up regulated in mouse spleens (
Table S6).
Induction of genes implicated in stress responses during infection
In our study, genes belonging to the three principal classes of stress genes were up regulated in the host. Class I genes encode classical chaperones and are controlled by the HrcA repressor. Nine of the 25 genes previously shown as HrcA repressed
[11] were activated
in vivo, including genes encoding the molecular chaperones DnaK and GroEL respectively also shown as induced in macrophages and required for survival following phagocytosis
[45],
[46] (
Table S7). Inversely, 17 of the 36 genes shown to be indirectly positively regulated by HrcA
[11], were up regulated in mouse spleens. This list includes genes encoding ribosomal proteins, as well as a number of DNA replication, transcription or translation related genes.
The class II stress response is mediated by sigma B (SigB). A total of 30 genes that have been recently classified as SigB activated
[8] appeared here up regulated
in vivo (
Table S7). In particular we detected the up regulation of
inlH [24],
ltrC implicated in response to cold shock
[47], and
lmo1601 similar to general stress proteins. Interestingly, 40 genes previously classified as down regulated by SigB during the stationary growth phase
[8] were detected as activated
in vivo (
Table S7). These include
kat, a catalase involved in the oxidative stress response
[48], a large proportion of genes encoding ribosomal proteins or implicated in translation, cell division and cell wall biogenesis. Furthermore,
iap, the P60 gene
[29], is part of this group. Finally,
rsbU and
rsbX, two components of the complex regulation system of SigB
[49] were also up regulated in mouse spleens (
Table S7).
CtsR is a transcriptional repressor involved in the control of class III stress proteins and previously shown to be responsible for the repression of 42 genes
[12], 15 of which appeared up regulated in the host (
Table S7). In particular, CtsR regulates the expression of Clp proteases required for the degradation of abnormal proteins and implicated in bacterial escape from macrophage vacuoles and virulence in mice
[50]. Expression of
clpBCE was activated during infection, as well as
mcsA and
mcsB the modulators of the CtsR regulon.
In some host cells, bacteria are confronted with severe oxidative stress due to the release of reactive oxygen intermediates. We observed the
in vivo activation of an important number of oxidative stress resistance mechanisms. The
qoxABCD operon that encodes a quinol oxidase important for oxidative stress response, and two major proteins implicated in protection against superoxides and reactive oxygen species (ROS), Kat and Sod, were highly up regulated
in vivo (
Table S7). Sod was previously shown as required for
Listeria full virulence and is a target of Stp, a serine-threonine phosphatase also involved in
L. monocytogenes virulence
[44],
[51], and detected down regulated in the host. A decrease in the level of Stp was previously associated to an increase in phosphorylated Sod, accompanied by the secretion of active non-phosphorylated Sod by the SecA2 system
[44],
[51]. Furthermore, genes encoding a thioredoxin and two thioredoxin reductases involved in the response to oxidative stress (
lmo2152,
trxB and
lmo2390) were up regulated in our study (
Table S7). The ferritin protein Fri, that also provides protection against reactive oxygen species, is essential for virulence and is required for efficient bacterial growth at early stages of the infection process
[52],
[53]. Fri transcription is directly regulated by Fur, the ferric uptake regulator. The expression of
fri and
fur was activated during infection. In addition,
ohrA and
gap were up regulated
in vivo and encode two proteins respectively involved in hydroperoxide resistance
[54] and in resistance against reactive oxygen species produced by host phagocytic cells in
Leishmania [55] (
Table S7).
L. monocytogenes metabolism adaptation to in vivo conditions
Remarkably, 30% of the
in vivo regulated genes are involved in
L. monocytogenes metabolism (99 metabolism-related genes were up and 72 were down regulated) (
Table S8). As described above,
uhpT is an
in vivo highly up regulated virulence gene, regulated by PrfA and that promotes the uptake of phosphorylated hexoses (glucose-1-phosphate and glucose-6-phosphate)
[23],
[56]. Phosphorylated glucose is the product of glycogen hydrolysis in eukaryotic cells and there is experimental evidence that the PrfA-dependent utilization of this compound has a role in
L. monocytogenes virulence
[23],
[56].
We observed an
in vivo up regulation of several genes encoding enzymes involved in the glycolysis, like
gap,
pgi,
fbaA, and
pgm. Inversely, we found a down regulation of the expression of four genes involved in the non-oxidative phase of the pentose phosphate pathway (
lmo2660,
lmo2661,
lmo2662 and
lmo2674). The final step of glycolysis leads to pyruvate, which is then converted to acetyl-CoA by the pyruvate dehydrogenase complex. We found this complex partly up regulated
in vivo, as well as one of its activator, the lipoate ligase protein LplA2
[57],
[58]. The citric acid cycle is continuously supplied with acetyl-CoA during aerobic respiration. We observed an up regulation of three citric acid cycle genes (
citBCG) (
Table S8). The citric acid cycle is followed by oxidative phosphorylation. In this study, we found the up regulation of several genes implicated in biosynthesis and assembly of components of the respiratory chain (
menD,
lmo1677,
qoxABD,
ctaA,
cydA,
cydD,
atpD). In addition, genes encoding
resD, a regulator of aerobic and anaerobic respiration
[59] and
rex, a redox-sensing transcriptional repressor
[60], were also up regulated
in vivo. Genes encoding the pyruvate-formate lyase (
pfl) and pyruvate-formate lyase activating enzymes (
pflCB) are required for the anaerobic metabolism of pyruvate and were activated in the host (
Table S8).
Genes implicated in amino acid biosynthesis were also induced
in vivo, in particular
aroA and
pheA, two genes responsible for aromatic amino acid biosynthesis. Mutations in
aroA and
pheA were previously shown to induce an attenuation of virulence in the mouse model
[61],
[62]. Furthermore, the expression of genes implicated in the biosynthetic pathways of branched amino acids (
alsS,
ilvN and
lmo0978), and amino acids of the aspartate and glutamate families (
ansB,
lmo0594,
lmo1006,
lmo1011,
lmo2413 and
glnA,
lmo2770, respectively), was also increased
in vivo (
Table S8).
Significantly, mannose (
lmo0781–lmo0784), maltose (
lmo0278) and cellobiose (
lmo0301 and
lmo0915) -specific PTS encoding genes
[63] were up regulated
in vivo. Inversely, fructose (
lmo2733), galactitol (
lmo2665) and mannitol (
lmo2649) -specific PTS encoding genes appeared down regulated.
Among the genes involved in bacterial ion uptake systems, a potassium-transporting ATPase encoding gene (
kdpB) was down regulated
in vivo. Cobalt (
lmo1207), manganese (
lmo1424) and calcium (
lmo0841) transporter systems were, inversely, up regulated. As indicated above, the ferritin and ferric uptake protein encoding genes,
fri and
fur, shown to be activated under low iron concentration
[64],
[65], appeared highly up regulated
in vivo (
Table S8).
Detection of potential virulence genes by in vivo transcriptomics
A major goal of this work was the identification of genes that encode proteins that may play a role in the infectious process. To identify such virulence genes and in order to establish a short list, we arbitrarily used several criteria. The gene should be preferentially 1) highly activated during infection; 2) absent in the non pathogenic strain L. innocua and present in other L. monocytogenes strains from diverse serotypes; 3) a member of a specific protein family encoding gene (surface, secreted, stress) possibly involved in virulence; 4) controlled by virulence regulators (PrfA, VirR, CtsR, HrcA, SigB). Several candidates emerged, matching, at least, some of the above criteria ().
| Table 6L. monocytogenes EGDe genes differentially regulated in the host and potential virulence factors. |
lmo0206,
lmo0257,
lmo0915,
lmo1290 and
lmo2157 are genes that, as eleven already known virulence factors, are
L. monocytogenes species-specific and induced
in vivo.
lmo0206 and
lmo2157 are the only two genes activated
in vivo, controlled by PrfA, absent from
L. innocua and whose role in virulence was never investigated.
lmo0206,
orfX [66], is in addition located at the end of the
Listeria virulence cluster and was recently implicated in intracellular survival
[14]. The expression of
lmo2157 was shown to be controlled by PrfA and SigB
[6],
[8].
lmo1081, lmo1082, lmo1099 and lmo1102 are L. monocytogenes EGDe species-specific genes highly up regulated in vivo over the three time points of the infection (). Interestingly, these genes encode proteins potentially involved in cell wall metabolism and heavy metal detoxification.
Only two uncharacterized genes encoding LPXTG surface proteins (
lmo1290 and
lmo2714) and three encoding GW surface proteins (
lmo1521,
lmo2203 and
lmo2713) were up regulated within the host ().
lmo1521 and
lmo2203 are in addition predicted autolysins.
lmo2713 and
lmo2714 seem to be part of a genomic region over expressed at all time points of the infection and Lmo2714 was found in the
Listeria culture supernatant
[43]. Four genes (
lm0540,
lmo1438,
lmo1855 and
lmo2522) predicted to be involved in cell wall metabolism were up regulated
in vivo, and similar to
pgdA,
iap, and
murA [29]–
[31], could participate in
Listeria infection.
Twenty-five uncharacterized genes activated
in vivo encode secreted proteins that may interact with the host cells, including Lmo2201, a Tat-secreted protein
[42], and GAPDH. GAPDH was previously shown to be part of the
Listeria cell wall subproteome
[35], and to impair
Listeria phagosome maturation
[67]. GAPDH seems, in addition, to be implicated in the virulence of several other pathogens
[68]–
[70].
lmo0788 is highly activated in mouse spleens during infection and is the only gene of the group I PrfA-regulated genes (i.e. preceded by a PrfA-box and positively regulated by PrfA)
[6] whose role during infection has never been addressed ( and ).
lmo0788 encodes a protein similar to subunits (BadFG) of the benzoyl-CoA reductase used by facultative aerobes in absence of oxygen for reductive aromatic metabolism
[71].
VirR appears as the second main regulator of virulence genes and controls
lmo0604,
lmo1742,
lmo2114,
lmo2115,
lmo2177 and
lmo2439, whose expression was activated in the host ( and ).
lmo2114 and
lmo2115 are in addition part of a transcriptional unit co-regulated by CtsR and SigL
[7].
Several stress protein encoding genes that are under the control of different stress regulators were up regulated
in vivo. In particular,
lmo2048 is a stress protein-encoding gene that is co-controlled by CtsR and HrcA (). The 19 genes up regulated
in vivo and co-controlled by PrfA and SigB () could also be important for the infectious process. Among these,
lmo1601 and
lmo1602 are furthermore regulated by SigL
[7].
The use of such arbitrary criteria obviously not guaranteed that a selected gene was a virulence factor, and conversely probably excluded many virulence genes. In particular, it is worth mentioning that 91 genes encoding proteins similar to unknown proteins, and 31 encoding putative proteins with no similarity in public databases were differentially expressed
in vivo (
Table S3), representing a large reservoir of potential new virulence factors. Of these genes, those highly regulated all over the infectious process could be of special relevance for virulence.
Identification of new L. monocytogenes virulence factors
In order to validate our transcriptomics approach and identify new L. monocytogenes virulence factors, 6 genes (lmo1081, lmo1082, lmo1102, lmo2713, lmo2714 and gap) were selected for mutagenesis using the criteria presented above. As we were unable to produce a gap deletion mutant (probably because GAPDH is an essential protein), we constructed a GAPDH secretion mutant.
To analyze the potential role of the selected genes in virulence, we performed intravenous inoculations of BALB/c mice with wild type (wt) and mutant strains, and the number of bacteria in the mouse liver and spleen was determined 72 h after infection (). Mutants can be classified with respect to their virulence potential. Bacterial counts for lmo1081 and lmo2713 mutants were not significantly changed as compared to the wt strain, suggesting the non-implication of these genes in Listeria virulence in mice. For the lmo1082 mutant, bacterial counts were significantly affected (≈1 log) in mouse livers and at a lesser extent in the spleens. Interestingly, for lmo1102, lmo2714 and gap mutants we observed a remarkable decrease of bacterial counts in both mouse organs as compared to the wt. In particular, the number of bacteria was dramatically impaired in the liver 72 h after inoculation (≈2,5 to 4,5 log). The gap mutant appeared as the most attenuated mutant of our analysis with a considerable virulence decrease in both organs reaching 3,5 log in the spleen and 4,5 log in the liver as compared to the wild type ().
In order to better characterize virulence attenuated strains, mutants for
lmo1082,
lmo1102,
lmo2714 and
gap were complemented. The corresponding wild-type gene was inserted as a single copy under the control of its own promoter on the chromosome of the mutant strain, at the PSA bacteriophage attachment site using the pPL2 integration vector
[72]. Wild type, mutant and complemented strains were tested for growth in BHI at 37°C and for intracellular behavior after internalization in the murine macrophage cell line J774 ().
The growth rate observed in BHI at 37°C for the majority of the strains tested was comparable to that of the wild type (). However, the
gap secretion mutant exhibited an important
in vitro reduced growth rate and reduced density at the stationary phase. The growth defect observed for the
gap mutant was even accentuated in the complemented strain (). This is most probably the result of an over expression of intracellular GAPDH, expressed at the same time from the bacterial genome and from the plasmid harbored by this strain. Surprisingly, the
prsA2 mutant presented also a notable growth delay. This growth defect was not mentioned in previous studies implicating PrsA2 in intracellular behavior and virulence
[14],
[33].
Wild type, mutant and complemented strains were also tested for intracellular behavior. As shown in , all the strains grew with similar multiplication rates after internalization in J774 cells, indicating that the slight growth delay observed in BHI at 37°C for some strains has no consequences on intracellular multiplication.
In addition, complemented strains were analyzed after intravenous inoculations of BALB/c mice as compared to wt and mutant strains, and the number of bacteria in the mouse liver and spleen was determined 72 h after infection (). The virulence phenotype was restored, albeit partially in the case of lmo2714, in complemented strains, except for the gap mutant. The virulence defect of the gap complemented strain was even more severe in the spleen as compare to the corresponding mutant (). This was in correlation with the increased growth defect observed in BHI at 37°C for the gap complemented strain.
These results revealed a role for lmo1082 and lmo1102, and at less extent for lmo2714 and gap in Listeria virulence, validating our in vivo transcriptomics approach.