- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Comparative and Functional Genomics
Volume 2009 (2009), Article ID 914762, 16 pages
Comparative Analyses of Transcriptional Profiles in Mouse Organs Using a Pneumonic Plague Model after Infection with Wild-Type Yersinia pestis CO92 and Its Braun Lipoprotein Mutant
1Department of Microbiology and Immunology, The University of Texas Medical Branch, Galveston, TX 77555-1070, USA
2Virginia Bioinformatics Institute, Virginia Polytechnic and State University, Blacksburg, VA 24061, USA
Received 26 July 2009; Revised 28 September 2009; Accepted 18 October 2009
Academic Editor: Antoine Danchin
Copyright © 2009 Cristi L. Galindo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
We employed Murine GeneChips to delineate the global transcriptional profiles of the livers, lungs, and spleens in a mouse pneumonic plague infection model with wild-type (WT) Y. pestis CO92 and its Braun lipoprotein () mutant with reduced virulence. These organs showed differential transcriptional responses to infection with WT Y. pestis, but the overall host functional processes affected were similar across all three tissues. Gene expression alterations were found in inflammation, cytokine signaling, and apoptotic cell death-associated genes. Comparison of WT and mutant-infected mice indicated significant overlap in lipopolysaccharide- (LPS-) associated gene expression, but the absence of Lpp perturbed host cell signaling at critical regulatory junctions resulting in altered immune response and possibly host cell apoptosis. We generated a putative signaling pathway including major inflammatory components that could account for the synergistic action of LPS and Lpp and provided the mechanistic basis of attenuation caused by deletion of the lpp gene from Y. pestis in a mouse model of pneumonic plague.
The gram-negative bacterium Yersinia pestis is the etiological agent of plague. Y. pestis is transmitted to humans through the bite of an infected flea or inhalation of the organisms, resulting in bubonic, pneumonic, or septicemic forms of plague . Y. pestis has attracted much interest recently because of its potential as a weapon of bioterrorism. Following entry within a host, Y. pestis evades the host immune system and replicates in the lymph nodes, ultimately leading to lymph node necrosis and death if untreated [2–4]. Histological evidence indicates that bacteria within neutrophils are killed, while bacteria within macrophages and dendritic cells survive and go on to express various virulence determinants, which allow bacterial growth and their eventual release from the macrophages [5–7]. For example, F1 (capsular) antigen  and type III secretion system (T3SS) effectors  are expressed only at 37°C and have been shown to modulate the host response so that Y. pestis becomes resistant to subsequent phagocytosis. The use of these antiphagocytic mechanisms has led researchers to suggest that Y. pestis is predominantly an extracellular pathogen in the mammalian host [9, 10]. However, a strong cell-mediated immune response to Y. pestis infection is seen in immunized mice, suggesting that immune cells are also needed to clear either intracellular bacteria or extracellular Y. pestis that have been opsonized. A T-cell component of protection against Y. pestis, in the absence of antibody, has been established [11, 12]. In unvaccinated individuals, low doses of Y. pestis can be resolved following combined treatment with the T helper1- (Th1-) associated cytokines interferon (IFN)- and tumor necrosis factor (TNF)- . These studies suggest that cell-mediated immune responses are important for protection against Y. pestis.
The ability of Yersinia species to infect and replicate within a host is primarily due to the bacterial expression and implementation of the T3SS . T3SS is comprised of a molecular syringe-like complex that injects effector molecules into the target host cell enabling the bacteria to inhibit innate and acquired immune functions as well as to induce apoptosis. There are specific Yersinia outer membrane proteins (Yops) that have been studied extensively and characterized as inhibitors of specific biological processes that promote the survival of Yersinia species within the host. Specifically, the proteins YopE, -H, -J, -M, -O, -P, and -T disrupt cytoskeletal dynamics, inhibit innate and acquired immune functions, and promote apoptosis [15, 16].
The outer membrane of gram-negative bacteria is comprised of many different proteins that help maintain the structural integrity of the bacterial cell envelope. One particularly abundant lipoprotein, designated murein (or Braun) lipoprotein (Lpp), is associated with the outer membranes of bacteria within the family Enterobacteriaceae . Earlier studies indicated that Lpp (6.3 kDa) from enteropathogenic bacteria not only synergized with lipopolysaccharide (LPS) to induce septic shock but also evoked the production of TNF- and interleukin 6 (IL-6) in both LPS-responsive and LPS-nonresponsive mice and in mouse peritoneal exudate macrophages, suggesting an alternative signaling mechanism for Lpp . In fact, a subsequent study showed that Lpp signals through Toll-like receptor-2 (TLR-2) and not TLR-4, which LPS utilizes for cell signaling . Our more recent data provided evidence that lpp mutants of Y. pseudotuberculosis and Y. pestis KIM/D27 were attenuated in mice, an effect that could be complemented . In the latter strain of Y. pestis KIM/D27, a 102-kb pigmentation locus (pgm) was deleted, resulting in the attenuation of the WT bacteria . Importantly, immunization of mice with this mutant provided protection to animals against pneumonic plague invoked by intranasal inoculation of Y. pestis CO92 .
Most bacterial virulence genes are regulated upon entering the host. Global regulators have the ability to modulate multiple operons that belong to different metabolic pathways and are important for bacteria to adapt to new conditions. Recently, several laboratories have established an intranasal mouse model of the pneumonic plague infection and have investigated the host-pathogen interaction by pathological survey and bacterial gene expression microarrays . Liu et al.  examined the transcriptional profile of mice infected with Y. pestis strain 201, which is avirulent in humans, and reported upregulation of host cytokines that might mimic what would be observed during human infection . In our study, we investigated the transcriptional profiles of mice challenged by the intranasal route with Y. pestis CO92, a clinical isolate that is virulent in both mice and humans that would presumably better model human disease. We also examined the transcriptional profile of a Y. pestis CO92 lpp mutant and compared the results to mice infected with WT bacteria and found that Y. pestis CO92 lpp mutant infection of the lung caused upregulation of many genes encoding major proteins of the host immune system. Interestingly, we found a number of unique genes which were expressed differently in all three tissues of mice infected with the lpp mutant but were not altered by WT Y. pestis CO92 infection. This study provided new information on the dynamic of the liver, lung, and spleen host transcriptional responses to infection with WT Y. pestis CO92 and its lpp mutant.
2. Materials and Methods
2.1. Bacterial Strains
WT Y. pestis CO92 was obtained from the Centers for Disease Control and Prevention (CDC, Atlanta, GA) and maintained in our restricted access biosafety level- (BSL-) 2 laboratory. The construction and characterization of the strain deficient in the expression of the lpp gene were previously described in detail . All bacteria were grown in Brain Heart Infusion broth (BHI, Difco, Voigt Global Distribution Inc, Lawrence, KS) at 28°C prior to infection of mice.
2.2. Animal Studies
Swiss-Webster female mice (Charles River Laboratories, Wilmington, MA) 5-6 weeks of age were infected intranasally with 5 LD50 of either WT or lpp mutant of Y. pestis CO92 . Uninfected mice were used as controls. At either 12 or 48 hours post infection (p.i.), 3 mice per group were euthanized and the lungs, livers, and spleens were harvested and homogenized in 1 mL of RNALater (Ambion/Applied Biosystems, Austin, TX) using 50 mL tissue homogenizers (Kendell, Mansfield, MA). RNA was isolated from the tissue homogenates and purified using RNAqueous (Ambion). After an overnight precipitation, the RNA was resuspended in 20 L of diethylpyrocarbonate- (DEPC-) treated water and hybridized to Affymetrix GeneChip Mouse Genome 430 2.0 arrays, performed by the Molecular Genomics Core at UTMB Galveston, Texas, per manufacture protocols. The arrays had 45,000 probe sets representing more than 39,000 transcripts derived from 34,000 well-substantiated mouse genes. The experiments were performed in triplicate (biological replicates), generating a total of 45 arrays.
2.3. Normalization and Initial Characterization of Arrays
Data were Robust Multichip Average (RMA) normalized and log transformed using GeneSifter software (VizX Labs, Seattle, WA). Based on regression analysis of experimental replicates, there was an acceptable level of variation between each array (). Raw and processed data (a total of 45 arrays) were deposited in the Gene Expression Omnibus (GEO) online (http://www.ncbi.nlm.nih.gov/geo/) database (Accession GSE18293).
2.4. Analysis Methods
After normalization, data were considered separately for the three tissue types of mice. Data were further separated, based on time post infection and gene expression alterations that occurred in response to the WT Y. pestis CO92 and its lpp mutant. This resulted in four analyses per tissue type (uninfected versus WT-infected animals and WT-infected versus mutant-infected animals at 12 hours and 48 hours p.i.). ANOVA was performed for each comparison, and only genes with values of were considered for further analyses. Subsequent filtering was performed dependent upon group comparison types, as detailed below. Hierarchical clustering was employed on normalized and log transformed signals using GenSpring GX 10.0 (Agilent Technologies, Santa Clara, CA).
2.5. Data Analysis of Uninfected Controls versus WT-Infected Animals
For each time point, normalized signal values were averaged and pairwise comparisons were performed using GeneSifter. Only alterations (control versus WT-infected) of at least 2.0-fold were considered for further studies. Student’s -test with Benjamini and Hoshberg correction was also performed using GeneSifter. However, only the value without correction was used to filter data (), because natural biological variation was greater for some tissues than for others. All possible individual pairwise comparisons were performed using Spotfire DecisionSite 9.0 software (Spotfire, Inc., Sommerville, MA). An alteration of at least 1.5-fold was expected for each of the 9 possible comparisons between controls versus WT-infected samples (for each time point). Any alteration observed between uninfected and WT-infected animals was expected to be at least 50% greater than the fold change calculated for each uninfected control (C1 versus C2, C2 versus C3, and C1 versus C3).
2.6. Data Analysis of WT versus lpp Mutant-Infected Animals
For each time point, normalized signal values were averaged and pairwise comparisons were performed using GeneSifter. Only alterations (WT- versus lpp-infected) of at least 1.5-fold were considered for further analysis. Student’s test was performed using GeneSifter, with the expectation of a . All possible individual pairwise comparisons were performed using Spotfire DecisionSite 9.0 software (Spotfire, Inc.). An alteration of at least 1.5-fold was expected for each of the 9 possible comparisons between WT-infected and lpp-infected samples for each time point. Any alterations observed between WT-infected and lpp mutant-challenged animals were expected to be at least 50% greater than the fold change calculated for each uninfected control (C1 versus C2, C2 versus C3, and C1 versus C3). An alteration of at least 2.0-fold (on average) was expected between either uninfected versus WT-infected or uninfected versus lpp mutant-infected samples. This step was intended to eliminate any presumably spurious alterations observed between WT-infected and lpp mutant-challenged animals that was not normally affected by infection or altered in response to the lpp mutant as compared to healthy animals.
3.1. General Gene Expression Changes in All Tissues from WT Y. pestis CO92- and Its lpp Mutant-Infected Mice
The host transcriptional responses to infection with WT Y. pestis and its lpp mutant in an inhalational mouse model of pneumonic plague were studied. Mice were infected for 12 or 48 hours, and RNA was isolated from livers, lungs, and spleens for processing and hybridization to Affymetrix GeneChip Mouse Genome 430 2.0 arrays. Uninfected animals served as controls, and experiments were performed in triplicate, which generated a total of 45 arrays (see summary in Supplementary Table 1 available online at http://dx.doi.org/10.1155/2009/914762). A stringent data analysis method was employed, including analysis of replicate samples and subsequent elimination of naturally variable transcripts, to increase the reliability of the results and to greatly minimize false positives. Each experimental group was separately compared to appropriate uninfected control tissues, and the overall results of the analysis are shown in Table 1. In our initial GeneChip experiments, we also used times points of 24, 36, and 60 hours (data not shown); however, optimal transcriptional profiling changes were observed at 12 and 48 hours p.i.. Consequently, we focused only on these two time points.
Three general trends were apparent regarding the overall response of the host to Y. pestis infection: (1) host transcriptional responses increased dramatically between 12 hours and 48 hours p.i., (2) the liver transcriptome was more profoundly perturbed, compared to spleen or lung tissues, in WT Y. pestis-infected mice, and (3) there was a drastic difference in host transcriptional responses of mice infected with the lpp mutant, dependent on time course and tissue. For example, many genes were differentially expressed in lpp mutant-infected lungs (109 transcripts) and livers (256 transcripts) at 48 hours p.i., compared to WT Y. pestis-infected animals, whereas only modest differences were observed at the earlier time point (12 hours p.i.) in the spleens (25 transcripts) and livers (11 transcripts) of mice infected with the lpp mutant compared to WT bacteria (Table 1).
Hierarchical clustering of normalized and log transformed signal values for genes that were differentially expressed between the various tissue and infection types likewise indicated that the majority of gene expression differences between uninfected animals and mice infected with WT Y. pestis CO92 occurred at 48 hours post infection (Figure 2(a)). Moreover, there were some transcripts that represented a “generalized’’ host response at this latter time point, which is demonstrated in Figure 2(a) by the clustering of all 9 samples representing liver, lung, and spleen replicate samples from mice infected with the WT bacteria (bright red, right-hand side of Figure 2(a)). These higher expressed transcripts were further separated based on tissue type, as expected (Figure 2(a)), indicating that there was a high correlation between replicate samples for these differentially expressed transcripts. A similar concordance was obtained when the signal intensity values from Y. pestis CO92 lpp mutant-infected tissues were clustered. There were distinct transcriptional changes that characterized the livers (Figure 2(b)), lungs (Figure 2(c)), and spleens (Figure 2(d)) of mice infected with the lpp mutant, compared to uninfected animals and animals infected with the WT Y. pestis CO92.
The complete list of gene expression alterations in response to infection with WT Y. pestis is provided as Supplementary Table 2. Altered genes were mainly associated with immune responses and inflammation. For instance, CD14, several chemokines, INF-, interleukin 1 receptor, serine peptidase inhibitor, members 3G and 3N of clade A, and genes encoding guanylate binding proteins were observed to increase across the liver, lung, and spleen of WT infected mice at 12 or 48 hours p.i. (Supplementary Table 2). There were 30 genes whose expression was increased in all three tissues (liver, lung, and spleen) in response to WT Y. pestis infection, compared to uninfected control animals at 48 hours p.i. (Table 2). These genes also represented mainly immune and stress response functions. Importantly, no genes were altered across all tissue types at 12 hours p.i. with WT Y. pestis CO92. Although there were substantial differences in individual gene expression changes observed in the three different tissue types (liver, lung, and spleen) in response to WT Y. pestis infection (Figure 2(a)), overall functional processes, as determined using Ingenuity pathway analysis software, were remarkably similar (Figure 3).
3.2. Gene Expression Profiling of the Liver from WT Y. pestis CO92-Infected Mice
A total of 72 genes were altered in expression (33 upregulated and 39 downregulated) in the livers of mice infected with WT Y. pestis at 12 hours p.i. (Table 1). Upregulated genes were mainly those involved in stress and acute-phase responses, signal transduction, and regulation of various metabolic processes, while downregulated genes included those involved in the regulation of cell proliferation and differentiation, apoptosis, and immune cell activation. Contrary to what was observed at the earlier time point, there were a substantial number (1,407) of genes altered by WT Y. pestis infection in the liver of mice at 48 hours p.i. (966 upregulated and 441 downregulated, Table 1). Based on the KEGG report obtained using GeneSifter (Supplementary Table 3), the signaling pathways with which upregulated genes were significantly associated included those important for immune response signaling, cell adhesion, apoptosis, and stress responses. Downregulated genes were mainly those involved in various metabolic processes.
3.3. Gene Expression Profiling of the Lung from WT Y. pestis CO92-Infected Mice
A total of 37 different genes were upregulated in response to WT Y. pestis in the lungs at 12 hours p.i. compared to uninfected mice. These genes included those that code for several chemokines (e.g., Ccl20, Ccl9, Cxcl1, Cxcl2, Cxcl5, and IL6), stress/acute-phase molecules (e.g., Orm2, Serpina3n, Saa1, Saa3, Gclm, Hspa1a, and Srxn1), and regulators of cell cycle progression and apoptosis (e.g., Cdkn1a, Nupr1, and MafF) (Supplementary Table 2). We noted 11 genes, including cysteine rich protein 61 and gene encoded D site albumin promoter binding protein (Dbp), whose expression was downregulated in response to WT Y. pestis in the lung at 12 hours p.i.. At 48 hours p.i., 192 genes were altered in the lungs of mice in response to infection with WT Y. pestis CO92 (Table 1). Similar to what was observed at 12 hours p.i., the vast majority of altered genes were upregulated (162 genes), and comparatively fewer genes were downregulated (30 genes). Upregulated genes were mainly those involved in immune and acute-phase responses, inflammation, cell cycle regulation, and apoptosis (Supplementary Table 2).
3.4. Gene Expression Profiling of the Spleen from WT Y. pestis CO92-Infected Mice
A total of 48 genes (44 upregulated and 4 downregulated) were significantly altered in the spleens of mice 12 hours p.i. with WT Y. pestis, compared to uninfected control animals. These genes were involved in the regulation of transcription, cell growth and differentiation, and immune-specific functions (e.g., CD69, Cxcr4, Igh-6, and Src-like adaptor protein, Supplementary Table 2). At 48 hours, 77 genes were significantly upregulated and only 1 downregulated in the spleens of mice infected with the WT Y. pestis CO92 compared to uninfected control animals (Supplementary Table 2). Many of the upregulated genes seen at 48 hours, compared to uninfected control, were due to the host INF- response, as evident by Ingenuity pathway analysis of upregulated genes (Figure 3). Specifically, there were 18 transcripts (e.g., CXCL6, IL-1, and IL1R2), in addition to IFN- itself, that directly participate in IFN- signaling and were statistically upregulated in WT Y. pestis-infected mouse spleens. The mitogen-activated protein kinases ERK1/2 and JNK, transcription factors NF-B, CREB, and Akt, and apoptosis-associated caspases, all of which are typically regulated via nontranscriptional mechanism (e.g., phosphorylation), are integral components of IFN- signaling and were thus likely activated (Figure 3).
3.5. Comparison of WT Y. pestis CO92-Infected Mice to a Previous Study Utilizing a Strain (Y. pestis Strain 201) That Is Avirulent in Humans
Our study is the first to examine host global transcriptional responses to Y. pestis CO92 using an inhalation mouse model. However, Liu et al.  performed a similar study in Balb/c mice using a strain that is highly virulent in mice but not humans (Y. pestis strain 201). Because the entire gene expression data sets were not made publically available, we performed a comparison of the published results, which consisted primarily of cytokines and cytokine-related signaling molecules in order to gain insight into the potential differences and similarities in most responses to these two different strains. Despite the fact that different mice, array platforms, and analysis methods were employed, the majority of genes reported by Liu et al. as altered in response to infection with Y. pestis stain 201 were also identified as altered in our experiments (Table 3). However, there were some fundamental differences that could be contributed by the differences in virulence of the two strains. Most notably, all three tissue types (liver, lung, and spleen) responded similarly to infection with WT Y. pestis CO92 (Figure 2), particularly when major functional processes were considered (Figure 1), rather than individual gene alterations. However, some genes were altered differently depending on tissue type in mice infected with Y. pestis strain 201. Cd9, for instance, which was upregulated (2.8-fold) in our study in the livers of WT Y. pestis CO92-infected mice was upregulated in the liver (2.4-fold) and downregulated (−2.6-fold) in the lung of mice infected with Y. pestis strain 201 . Although Cd9 was not deemed significantly altered in our study in the lung or spleen, it was upregulated on average in both tissues (1.5-fold, and 2.1-fold, resp., data not shown). Likewise, Icam2 was reported as upregulated (1.9-fold) in the liver and downregulated (−3.1-fold) in the lung of mice infected with Y. pestis strain 201, whereas we found Icam2 to be upregulated (2.8-fold) in the liver (Table 3) and unaffected in the lung in response to infection with WT Y. pestis CO92 (data not shown). Most notably, IFN-, which we found to be upregulated in the liver, spleen, and lung of WT Y. pestis CO92-infected mice (Supplementary Tables 2 and 3) and also identified as a critical signaling pathway based on Ingenuity analysis of the entire gene expression dataset (Figure 3), was reported as downregulated in Balb/c mouse lungs in response to infection with the 201 strain .
While our study is the first to examine the entire host transcriptome in response to WT Y. pestis CO92, Lathem et al. previously performed a cytokine analysis of lung homogenates from C57BL/6 mice that were infected with Y. pestis CO92 via the intranasal route . They found that WT Y. pestis infection induced upregulation of IL12p70, TNF, IFN-, MCP-1 (also called CCL2), and IL-6. As shown in Supplementary Table 2, we also detected statistically significant upregulation of IFN- (4.3-fold), MCP-1 (also called CCL2, 3.4-fold), and IL-6 (34.3-fold) in the lung of WT Y. pestis-infected mice. Likewise, we detected an increase in IL-12a (2.6-fold), IL12b (1.9-fold), and TNF- (1.9-fold), although these differences were not deemed statistically significant (data not shown).
3.6. Gene Expression Profiling of the Liver, Lung, and Spleen of Mice Infected for 12 Hours with the lpp Mutant of Y. pestis CO92
For each experimental infection with WT Y. pestis, an experiment was also performed using the Y. pestis lpp mutant to determine the contribution of bacterial Lpp to host transcriptional responses. Based on a direct comparison of WT Y. pestis-infected mice and animals challenged with the lpp mutant, very few gene expression differences were observed at 12 hours p.i. (Table 1). In the liver of infected mice 12 hours p.i., the lpp mutant induced upregulation of 5 transcripts and downregulation of 6 transcripts, most of which were bacterial metabolic genes included on the array as controls (i.e., 28 probe sets representing 6 different genes, Supplementary Table 4). These alterations most likely represented differences in bacterial load at 12 hours in the livers of mice infected with the lpp mutant compared to WT Y. pestis, which is consistent with histopathological analysis of liver tissue . On the host side, 4 genes were upregulated in response to infection with the lpp mutant but not WT Y. pestis, including apoptosis inhibitor 5 (Api5), which suggests that Lpp might influence the host apoptotic response to infection (Table 4) and is consistent with our recently published data .
No genes were detected as differentially expressed in the lungs of mice infected for 12 hours with the lpp mutant, compared to WT Y. pestis (Table 1). In contrast, transcriptional differences in the spleen between WT Y. pestis-infected and lpp-challenged mice were limited to the earlier time point (i.e., differences were observed at 12 hours only). Although these alterations were few (only 25 genes were differential between WT Y. pestis-infected and lpp-challenged mice), the differences were profound. For instance, 18 probe sets (Affymetrix repeated transcripts) representing 16 different genes that were upregulated in WT-infected mouse spleens were not recapitulated by infection with the lpp mutant (Supplementary Table 5). Most of these genes are involved in the regulation of cell growth, including stress-associated cell proliferation (e.g., cyclin D3 and ERBB receptor feedback inhibitor 1). Only one gene was uniquely altered in the absence of Lpp (transthyretin that encodes a prealbumin carrier protein associated with acute phase response ), which was upregulated in lpp mutant-infected mice, compared to control animals. Of greater interest, immune-specific transcriptional responses (e.g., antihuman CD37 antibody, immunoglobulin kappa chain variable 28, and nemo-like kinase [Nik]) were downregulated in WT Y. pestis-infected mice and upregulated in lpp mutant-challenged animals, compared to uninfected controls (Table 4). Conversely, apoptosis-associated transcripts (e.g., Bak1 and BCL2L1) were downregulated in WT Y. pestis-infected mice and upregulated in lpp mutant-challenged mice, compared to uninfected controls (Table 4).
3.7. Gene Expression Profiling of the Liver, Lung, and Spleen of Mice Infected for 48 Hours with the lpp Mutant of Y. pestis CO92
The majority of transcriptional differences in host response to WT Y. pestis and the lpp mutant occurred at 48 hours p.i. in the liver and lung of infected animals. Most of these represented genes were commonly upregulated or downregulated in both WT Y. pestis and lpp mutant-infected mice, but by differing magnitudes. For example, CD96 antigen, which is important for macrophage activation and phagocytosis, was upregulated 4.4-fold in the livers of mice infected for 48 hours with WT Y. pestis and 11.1-fold in lpp mutant-infected mouse livers (Supplementary Table 6). However, there were also more profound differences, in which genes were altered uniquely by either WT Y. pestis- or lpp mutant-infected mice (Figure 1).
In mouse livers at 48 hours p.i., there were 27 genes that were specifically upregulated in lpp mutant-infected mice but not in animals challenged with WT bacteria. Induction of these genes, which included those involved in immune-specific signaling, inflammation, and the regulation of apoptosis, is therefore presumably repressed in the presence of Lpp. There were also 41 genes that were downregulated in the livers of lpp mutant-infected mice but not in WT Y. pestis-infected animals, compared to control animals (Supplementary Table 6). Two of these genes (i.e., complement factor H-related 1 and V-set immunoglobulin domain containing 4) (Table 4) are involved in the regulation of host immune responses, but the majority are associated with various metabolic processes (e.g., amino acid metabolism and gluconeogenesis).
There were 109 genes that were differentially expressed in the lung between WT Y. pestis-infected mice and animals challenged with lpp mutant bacteria (Table 1). Seventy of these genes were modestly upregulated (1.5-fold to 34.3-fold, mean = ) in response to WT Y. pestis infection but more profoundly upregulated (3-fold to 172-fold, mean = ) in response to infection with the lpp mutant (Supplementary Table 7). In contrast to what was observed at 12 hours or at 48 hours in liver tissue, the majority of differentially expressed genes in lungs were those critical for immune and stress responses, inflammation, and apoptosis. For instance, IL-6 and CXCL2 (also called Mip-2) were upregulated in the lungs of WT Y. pestis-infected mice 34.3-fold and 18.2-fold, respectively, compared to uninfected mice. In lpp mutant-infected mice, on the other hand, IL-6 and CXCL2 were upregulated 172-fold and 169.1-fold, respectively (i.e., 5-fold and 9.3-fold larger inductions), compared to uninfected control mice.
A total of 39 genes were upregulated exclusively in the lungs of mutant-infected mice after 48 hours of infection (i.e., not altered in the lungs of the WT Y. pestis group of infected mice) (Supplementary Table 7). Most of these genes were associated with apoptosis, inflammation, immune responses, and signaling pathways critical for immune cell activation, including apoptosis regulators Birc3 and Bcl2a1a, CD53, CXCL14 (also called Mip-2), coagulation factors III and X, IL-22, early growth response 1, leukemia inhibitory factor, and prostaglandin E synthase (Supplementary Table 7).
Based on the transcriptional profiles of liver, lung, and spleen of mice, the most profound differences between animals infected with WT Y. pestis versus the lpp mutant, literature searches, and known signaling pathways available in various online databases (e.g., NCBI, Biocarta, and The Protein Lounge), we created a putative Lpp-associated signaling pathway (Figure 4). For instance, we inferred that the most likely pathway for the production of the multiple cytokines that were identified as increased based on microarray results is phosphorylation and activation of NF-B and JNK via TLR-2 and TLR-4 induced activation of mitogen-activated protein kinases (MAPKs) (Figure 4). IFN- signaling, which is peripherally associated with this same pathway, was also inferred to be activated in response to WT Y. pestis and would possibly complement LPS signaling via TLR-4 to explain the induction of cytokines, albeit somewhat blunted, in the absence of the lpp gene (Figure 4). Reduction in some host responses, identified as expression alterations in lpp mutant-infected animals, could be partially explained by Lpp-mediated inhibition of leukemia inhibitory factor (Lif) and Dusp16 (Table 4), which downregulate activation of NF-B and JNK, respectively. The most directly affected process, based on WT versus lpp-infected animals, was apoptosis, possibly via inhibition of prostaglandin E synthase (Ptges) (Table 4) [24, 25] and perturbation of relative ratios of mitochondrial factors (e.g., Bcl-2 family members) (Figure 4). Overall, the three main signaling pathways induced by WT Y. pestis were TLR-4, TLR-2, and INF- signaling, which culminated in the production of multiple inflammatory cytokines, also detected as upregulated in infected mice in all three tissues examined.
In the present study, a transcriptional ontological assessment of significantly modulated genes in the liver, lung, and spleen from WT Y. pestis CO92-infected mice revealed a number of up- and downregulated transcripts that were associated with immune mechanisms. For example, an increase in CD14 transcript was observed across liver, lung, and spleen of mice infected with WT Y. pestis CO92 at 48 hours p.i., but the gene encoding IL-10 was not upregulated. CD14 exists as a membrane-bound or soluble form and serves as a coreceptor with TLRs or LPS-binding protein to associate with LPS from Gram-negative bacteria . During Y. enterocolitica infection, CD14 complexes with TLR-2 on macrophages and subsequently binds low calcium response antigen V (LcrV), which leads to a reduction in TNF- and an increase in IL-10 . This IL-10 induction by LcrV through binding to TLR-2/CD14 plays a key role in Y. enterocolitica immune evasion and pathogenicity . However, previous studies on Y. pestis indicated that IL-10 was not produced in the lungs of mice infected intranasally, and TLR-dependent IL-10 induction by LcrV did not contribute to the virulence of Y. pestis . Our results are consistent with these findings and suggest that IL-10 suppression might be an important virulence mechanism for enteropathogenic yersiniae.
The majority of transcriptional alterations identified in the liver, lung, and spleen of WT Y. pestis mice were those important for host immune responses, as expected. In addition to CD14, TLR-4 and TLR-2 were upregulated p.i. (2.6- and 9.2-fold, resp., Supplementary Table 2), as were several downstream targets of these two TLR signaling pathways (Figure 4). We also identified the IFN- signaling pathway as a central player in the host response to WT Y. pestis infection (Figure 3). INF-, which was induced at 48 hours p.i. in all of the tissues, is produced by activated natural killer cells and T cells and is critical for a successful immune response to intracellular pathogens [29–31]. Also upregulated in WT Y. pestis-infected mice were the IFN--regulated serine proteases Serpina3g and Serpina3n (Supplementary, Table 2), which can inhibit caspase-independent death  and assist in the development of memory CD8 T cells . Likewise, we noted WT Y. pestis-induced upregulation of suppressor of cytokine signaling 1 (socs1, 8-fold in the liver) and socs3 (4-15-fold in all three tissues), which regulate JAK-STAT signaling, and TNF--induced protein 3 (tnfaip3b, 2.7-15.4-fold in the liver and lung), which is essential for negative regulation of I-B kinase/NF-B cascade (Supplementary Table 2).
Other IFN--regulated molecules that were induced in response to WT Y. pestis infection included several guanylate binding proteins (GBP2, 4, 7, and 6/10), which were upregulated in all of the tissues collected from WT-infected mice. This IFN--induced family of proteins has been poorly characterized, but they have been shown to regulate endothelial cell proliferation during infection, possibly by slowing cell-to-cell spreading . IGTP (IRGM3) and TGTP (IRGB6), members of the p47 GTPases family, were also increased in all of the tissues of WT bacteria-infected mice (Supplementary Table 2). These molecules are similar to the GBPs but do not require de novo synthesis of transcription factors . Functionally, they have been shown to localize to infected vacuoles in a Toxoplasma gondii infection , which is followed by vesicle formation, disintegration of the vacuole, and the subsequent demise of the parasites . Consequently, these guanylate binding proteins could perform a similar function during Y. pestis intracellular infection.
Consistent with a strong host inflammatory response to infection, multiple cytokines and chemokines were upregulated in WT Y. pestis-infected animals in all three tissues examined (Supplementary Tables 2 and 3). For instance, CXCL10 and CCL2, which were profoundly upregulated (6-30.5-fold) in response to WT Y. pestis infection, are chemottractants for monocytes, T cells, and dendritic cells. Likewise, neutrophils, important to the amelioration of early bacteremia, are attracted by CXCL6 , which was upregulated in the liver (23.1-fold), lung (6.4-17.5-fold), and spleen (14.6-fold) post infection (Supplementary Tables 2 and 3). Induction of some of these inflammatory chemokines (e.g., CCL3) would specifically attract monocytes, which may benefit Yersinia by providing a safe haven for replication [38, 39]. The compendium of host responses identified in this study supports a strong host inflammatory response that culminates in the activation of immune effectors downstream of TLR-2 and TLR-4 and subsequent amplification of the inflammatory responses via production of IFN-.
We noted an upregulation of Lipocalin 2 (Lcn2) from 6.5-fold at 12 hours to 67.3-fold at 48 hours and downregulation of the HFE2 gene (8.1-fold) in the livers of WT-infected mice at 48 hours p.i. (Supplementary Table 2). Both Lcn2 and HFE2 are associated with iron regulation, and mutation in the HFE2 gene is causative for hematochormasis, which is characterized by iron overload . The increase in Lcn2 by WT bacteria in the liver might check bacterial growth by binding to siderophores and could be a mechanism of mediating innate immune response. No change in its level, as observed in the lpp mutant, would cause normal bacterial growth in the liver. Based on our recent results (19), the lpp mutant grew normally in liver but not in the spleen or blood.
In conjunction with our assessment of host transcriptional responses in WT Y. pestis-infected mice, we also investigated the effects of an lpp mutant on gene expression. Recently, we demonstrated that in Y. pestis, deletion of the lpp gene from the pgm-locus KIM/D27 background strain further attenuated its virulence. However, minimal differences were noted in pathogenicity between the WT- and the lpp-mutant strain of CO92 in a pneumonic plague mouse model, probably because Y. pestis CO92 strain is highly virulent and deletion of one gene causes only increases in mean time to death [20, 41]. Interestingly, when groups of mice infected with either the WT CO92 or its lpp mutant were given a subinhibitory dose of levofloxacin, we observed a significantly higher survival rate, less severe histopathological changes, and reduced cytokine/chemokine levels in the lpp mutant-infected group compared to WT-infected mice . These data indicated that Lpp contributed to virulence of Y. pestis CO92 and was dependent on bacterial load. We used an intranasal mouse model of infection to study host gene expression alterations in the liver, lung, and spleen at 12 hours and 48 hours p.i. that demonstrates the distinctions of virulence and pathogenic mechanism(s) between WT and lpp mutant strains of Y. pestis CO92 in a pneumonic plague model.
Our first observation of mice infected with the lpp-mutant strain of Y. pestis CO92, compared to WT-infected animals, was that transcriptional responses that could be due to TLR-4 activation via LPS (e.g., chemokines, JAK-STAT signaling molecules, etc.) were blunted in the absence of lpp gene expression (Supplementary Tables 4–7), which supports a synergistic role for Lpp and LPS to induce septic shock as well as the LPS-like signaling previously observed in an LPS-nonresponsive background strain of mice . More interesting were transcriptional responses that were completely perturbed in the absence of lpp, such as activation in WT Y. pestis-infected animals but not in those infected with the lpp-mutant. These results provided much greater insight into Lpp-specific host signaling in the context of Y. pestis infection and allowed us to propose a putative signaling pathway (Figure 4) that could explain the intertwined roles of LPS and Lpp and also how Y. pestis might survive inside host cells.
As shown in Figure 4, WT Y. pestis induces the upregulation of TLR-4, TLR-2, and CD14 independently of Lpp (i.e., these molecules were also upregulated in mice infected with the lpp-mutant). However, the LPS and Lpp share a common downstream signaling pathway, and even in the absence of Lpp, these intermediate inflammatory effectors (e.g., Myd88, IRAK, mitogen activated kinases, STATs, NF-B, c-Jun and Fos, and various proinflammatory cytokines) were increased during Y. pestis infection (Supplementary Table 2 and Figure 4). Nontranscriptional events (e.g., Nik-mediated phosphorylation of IKK and subsequent degradation of IBs and nuclear location of NF-B) that are likely to have occurred based on the transcriptional profiles of Y. pestis-infected mice and classical signaling pathways are included for clarity. In the context of this WT model of infection, Lpp-specific signaling events were also apparent. Nik, for instance, is a crucial regulatory point downstream of TLR and cytokine receptor engagement, and its upregulation in WT Y. pestis infected mice was not recapitulated when the lpp-mutant was used. Other mechanisms of IB phosporylation and degradation would presumably occur in the absence of Lpp, since proinflammatory cytokines are still produced in the absence of the lpp gene.
Cell death was a major process identified as statistically overrepresented in all three tissue types, based on Ingenuity pathway analysis of altered genes (Figure 3). The balance of proapoptotic and antiapoptotic factors often determines cell fate, and apoptosis regulators can also function differently depending on cell type. Its regulatory complexity makes apoptosis-related transcriptional responses difficult to interpret. However, the absence of the lpp gene clearly perturbed the effects of the WT Y. pestis infection by subtly altering some apoptotic related transcription responses and specifically inducing or depressing others. For instance, the expressions of two genes encoding for Bcl2 family proteins (Bak1 and Bcl2l1) that function to induce apoptosis  were suppressed in the spleen of WT-infected mice but not in animals infected with the lpp-mutant (Table 4). Likewise, Hk1 was uniquely downregulated in only WT-infected mice, suggesting that its suppression requires the presence of bacterial Lpp. Whereas suppression of Bak1 and Bcl2l1 would likely be cytoprotective, cytochrome c release is inhibited by Hk1 , and therefore its decrease could lead to increased apoptosis .
This study provided the first comprehensive assessment of the host transcriptional profile in the lung, liver, and, spleen of mice intranasally infected with a highly virulent strain of Y. pestis CO92. We further investigated the contributions of bacterial Lpp to host transcriptional responses and presented a putative host signaling pathway that plausibly explained the synergistic actions of LPS and Lpp in the context of Y. pestis infection. Our results supported a model in which Y. pestis induced a strong inflammatory response, mediated by both LPS and Lpp, but evaded immune clearance, possibly by Lpp-induced inhibition of host cell apoptosis.
C. L. Galindo, S. T. Moen, E. V. Kozlova, and J. Sha contributed equally to the manuscript. This research was supported by NIH/NIAID Grants AI06438 and N01-AI-30065. S. L. Agar was a Predoctoral Fellow supported by the NIAID T32 Emerging and Tropical Infectious Diseases (AI07526) and Biodefense (AI060549) training grants. C. L. Galindo was supported by the NIH/NIAID Western Regional Center of Excellence in Biodefense.
- T. V. Inglesby, D. T. Dennis, D. A. Henderson, et al., “Plague as a biological weapon: medical and public health management,” Journal of the American Medical Association, vol. 283, no. 17, pp. 2281–2290, 2000.
- F. Sebbane, D. Gardner, D. Long, B. B. Gowen, and B. J. Hinnebusch, “Kinetics of disease progression and host response in a rat model of bubonic plague,” American Journal of Pathology, vol. 166, no. 5, pp. 1427–1439, 2005.
- W. W. Lathem, S. D. Crosby, V. L. Miller, and W. E. Goldman, “Progression of primary pneumonic plague: a mouse model of infection, pathology, and bacterial transcriptional activity,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 49, pp. 17786–17791, 2005.
- S. S. Bubeck, A. M. Cantwell, and P. H. Dube, “Delayed inflammatory response to primary pneumonic plague occurs in both outbred and inbred mice,” Infection and Immunity, vol. 75, no. 2, pp. 697–705, 2007.
- I. Adkins, M. Köberle, S. Gröbner, E. Bohn, I. B. Autenrieth, and S. Borgmann, “Yersinia outer proteins E, H, P, and T differentially target the cytoskeleton and inhibit phagocytic capacity of dendritic cells,” International Journal of Medical Microbiology, vol. 297, no. 4, pp. 235–244, 2007.
- D. C. Cavanaugh and R. Randall, “The role of multiplication of Pasteurella pestis in mononuclear phagocytes in the pathogenesis of flea-borne plague,” Journal of Immunology, vol. 83, pp. 348–363, 1959.
- R. T. Dean and W. Jessup, Eds., Mononuclear Phagocytes: Physiology and Pathology, Elsevier, New York, NY, USA, 1985.
- Y. Du, R. Rosqvist, and A. Forsberg, “Role of fraction 1 antigen of Yersinia pestis in inhibition of phagocytosis,” Infection and Immunity, vol. 70, no. 3, pp. 1453–1460, 2002.
- G. R. Cornelis, “The Yersinia Yop virulon, a bacterial system to subvert cells of the primary host defense,” Folia Microbiologica, vol. 43, no. 3, pp. 253–261, 1998.
- R. R. Brubaker, “Factors promoting acute and chronic diseases caused by yersiniae,” Clinical Microbiology Reviews, vol. 4, no. 3, pp. 309–324, 1991.
- J. F. Wong and S. S. Elberg, “Cellular immune response to Yersinia pestis modulated by product(s) from thymus derived lymphocytes,” Journal of Infectious Diseases, vol. 135, no. 1, pp. 67–78, 1977.
- S. T. Smiley, “Cell-mediated defense against Yersinia pestis infection,” Advances in Experimental Medicine and Biology, vol. 603, pp. 376–386, 2007.
- R. Nakajima and R. R. Brubaker, “Association between virulence of Yersinia pestis and suppression of gamma interferon and tumor necrosis factor alpha,” Infection and Immunity, vol. 61, no. 1, pp. 23–31, 1993.
- G. R. Cornelis, “Yersinia type III secretion: send in the effectors,” Journal of Cell Biology, vol. 158, no. 3, pp. 401–408, 2002.
- G. R. Cornelis, “The Yersinia YSC-YOP ‘type III’ weaponry,” Nature Reviews Molecular Cell Biology, vol. 3, no. 10, pp. 742–752, 2002.
- M. Aepfelbacher, R. Zumbihl, and J. Heesemann, “Modulation of Rho GTPases and the actin cytoskeleton by YopT of Yersinia,” Current Topics in Microbiology and Immunology, vol. 291, pp. 167–175, 2005.
- A. A. Fadl, C. L. Galindo, J. Sha, G. R. Klimpel, V. L. Popov, and A. K. Chopra, “Global gene expression of a murein (Braun) lipoprotein mutant of Salmonella enterica serovar Typhimurium by microarray analysis,” Gene, vol. 374, no. 1-2, pp. 121–127, 2006.
- H. Zhang, I. Kaur, D. W. Niesel, et al., “Lipoprotein from Yersinia enterocolitica contains epitopes that cross-react with the human thyrotropin receptor,” Journal of Immunology, vol. 158, no. 4, pp. 1976–1983, 1997.
- A. O. Aliprantis, R.-B. Yang, M. R. Mark, et al., “Cell activation and apoptosis by bacterial lipoproteins through Toll-like receptor-2,” Science, vol. 285, no. 5428, pp. 736–739, 1999.
- J. Sha, S. L. Agar, W. B. Baze, et al., “Braun lipoprotein (Lpp) contributes to virulence of yersiniae: potential role of Lpp in inducing bubonic and pneumonic plague,” Infection and Immunity, vol. 76, no. 4, pp. 1390–1409, 2008.
- H. Liu, H. Wang, J. Qiu, et al., “Transcriptional profiling of a mice plague model: insights into interaction between Yersinia pestis and its host,” Journal of Basic Microbiology, vol. 49, no. 1, pp. 92–99, 2009.
- T. Liu, S. L. Agar, J. Sha, and A. K. Chopra, “Deletion of Braun lipoprotein gene (lpp) attenuates Yersinia pestis KIM/D27 strain: role of Lpp in modulating host immune response, NF-B activation and cell death,” Microbial Pathogenesis, vol. 48, no. 1, pp. 42–52, 2010.
- F. M. Campbell, M. Waterston, L. O. Andresen, N. S. Sorensen, P. M. H. Heegaard, and P. D. Eckersall, “The negative acute phase response of serum transthyretin following Streptococcus suis infection in the pig,” Veterinary Research, vol. 36, no. 4, pp. 657–664, 2005.
- L. Lalier, P.-F. Cartron, F. Pedelaborde, et al., “Increase in PGE2 biosynthesis induces a Bax dependent apoptosis correlated to patients' survival in glioblastoma multiforme,” Oncogene, vol. 26, no. 34, pp. 4999–5009, 2007.
- S. J. Myung and I. H. Kim, “Role of prostaglandins in colon cancer,” The Korean Journal of Gastroenterology, vol. 51, no. 5, pp. 274–279, 2008.
- E. S. Van Amersfoort, T. J. C. Van Berkel, and J. Kuiper, “Receptors, mediators, and mechanisms involved in bacterial sepsis and septic shock,” Clinical Microbiology Reviews, vol. 16, no. 3, pp. 379–414, 2003.
- A. Sing, D. Reithmeier-Rost, K. Granfors, J. Hill, A. Roggenkamp, and J. Heesemann, “A hypervariable N-terminal region of Yersinia LcrV determines Toll-like receptor 2-mediated IL-10 induction and mouse virulence,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 44, pp. 16049–16054, 2005.
- D. Reithmeier-Rost, J. Hill, S. J. Elvin, et al., “The weak interaction of LcrV and TLR2 does not contribute to the virulence of Yersinia pestis,” Microbes and Infection, vol. 9, no. 8, pp. 997–1002, 2007.
- J. A. Hamerman, F. Hayashi, L. A. Schroeder, et al., “Serpin 2a is induced in activated macrophages and conjugates to a ubiquitin homolog,” Journal of Immunology, vol. 168, no. 5, pp. 2415–2423, 2002.
- G. Trinchieri, “Cytokines acting on or secreted by macrophages during intracellular infection (IL-10, IL-12, IFN-),” Current Opinion in Immunology, vol. 9, no. 1, pp. 17–23, 1997.
- E. Jouanguy, R. Doffinger, S. Dupuis, A. Pallier, F. Altare, and J. L. Casanova, “IL-12 and IFN-gamma in host defense against mycobacteria and salmonella in mice and men,” Current Opinion in Immunology, vol. 11, no. 3, pp. 346–351, 1999.
- N. Liu, Y. Wang, and P. G. Ashton-Rickardt, “Serine protease inhibitor 2A inhibits caspase-independent cell death,” FEBS Letters, vol. 569, no. 1–3, pp. 49–53, 2004.
- N. Liu, T. Phillips, M. Zhang, et al., “Serine protease inhibitor 2A is a protective factor for memory T cell development,” Nature Immunology, vol. 5, no. 9, pp. 919–926, 2004.
- E. Guenzi, K. Töpolt, E. Cornali, et al., “The helical domain of GBP-1 mediates the inhibition of endothelial cell proliferation by inflammatory cytokines,” EMBO Journal, vol. 20, no. 20, pp. 5568–5577, 2001.
- F. Aberger, A. P. Costa-Pereira, J. F. Schlaak, et al., “Analysis of gene expression using high-density and IFN--specific low-density cDNA arrays,” Genomics, vol. 77, no. 1-2, pp. 50–57, 2001.
- S. Martens, I. Parvanova, J. Zerrahn, et al., “Disruption of Toxoplasma gondii parasitophorous vacuoles by the mouse p47-resistance GTPases,” PLoS Pathogens, vol. 1, no. 3, article e24, 2005.
- Y. M. Ling, M. H. Shaw, C. Ayala, et al., “Vacuolar and plasma membrane stripping and autophagic elimination of Toxoplasma gondii in primed effector macrophages,” Journal of Experimental Medicine, vol. 203, no. 9, pp. 2063–2071, 2006.
- R. A. Lukaszewski, D. J. Kenny, R. Taylor, D. G. C. Rees, M. G. Hartley, and P. C. F. Oyston, “Pathogenesis of Yersinia pestis infection in BALB/c mice: effects on host macrophages and neutrophils,” Infection and Immunity, vol. 73, no. 11, pp. 7142–7150, 2005.
- B. J. Hinnebusch, “Bubonic plague: a molecular genetic case history of the emergence of an infectious disease,” Journal of Molecular Medicine, vol. 75, no. 9, pp. 645–652, 1997.
- G. Papanikolaou, M. E. Samuels, E. H. Ludwig, et al., “Mutations in HFE2 cause iron overload in chromosome 1q-linked juvenile hemochromatosis,” Nature Genetics, vol. 36, no. 1, pp. 77–82, 2004.
- S. L. Agar, J. Sha, W. B. Baze, et al., “Deletion of Braun lipoprotein gene (lpp) and curing of pPCP1 dramatically alter the virulence of Yersinia pestis CO92 in a mouse model of pneumonic plague,” Microbiology, vol. 155, no. 10, pp. 3247–3259, 2009.
- E. C. Pietsch, E. Perchiniak, A. A. Canutescu, G. Wang, R. L. Dunbrack, and M. E. Murphy, “Oligomerization of BAK by p53 utilizes conserved residues of the p53 DNA binding domain,” Journal of Biological Chemistry, vol. 283, no. 30, pp. 21294–21304, 2008.
- S. Abu-Hamad, H. Zaid, A. Israelson, E. Nahon, and V. Shoshan-Barmatz, “Hexokinase-I protection against apoptotic cell death is mediated via interaction with the voltage-dependent anion channel-1: mapping the site of binding,” Journal of Biological Chemistry, vol. 283, no. 19, pp. 13482–13490, 2008.