Scholarly article on topic 'Transcriptional and apoptotic responses of THP-1 cells to challenge with toxigenic, and non-toxigenic Bacillus anthracis'

Transcriptional and apoptotic responses of THP-1 cells to challenge with toxigenic, and non-toxigenic Bacillus anthracis Academic research paper on "Biological sciences"

Share paper
Academic journal
BMC Immunol
OECD Field of science

Academic research paper on topic "Transcriptional and apoptotic responses of THP-1 cells to challenge with toxigenic, and non-toxigenic Bacillus anthracis"

BMC Immunology

BioMed Central

Open Access

Research article

Transcriptional and apoptotic responses of THP-1 cells to challenge with toxigenic, and non-toxigenic Bacillus anthracis

Christopher Bradburne*1, Myung-Chul Chung2, Qin Zong3, Karen Schlauch4, Derong Liu3, Taissia Popova2, Anna Popova5, Charles Bailey2, Dan Soppet3 and Serguei Popov2

Address: Center for Bio/Molecular Science and Engineering Code 6900, US Naval Research Laboratory, Washington DC, USA, 2National Center for Biodefense and Infectious Disease, George Mason University, Manassas, VA, USA, 3Avalon Pharmaceuticals, Germantown, MD, USA, 4Department of Biochemistry and Molecular Biology, University of Nevada Reno, Reno, Nevada, USA and 5College of Arts and Sciences, University of Virginia, Charlottesville, VA, USA

Email: Christopher Bradburne* -; Myung-Chul Chung -;

Qin Zong -; Karen Schlauch -; Derong Liu -; Taissia Popova -;

Anna Popova -; Charles Bailey -; Dan Soppet -;

Serguei Popov -

* Corresponding author

Published: 13 November 2008 Received: 25 June 2007

, .,„„„,.,-, , ,„,,„,,, ,,„ Accepted: 13 November 2008

BMC Immunology 2008, 9:67 doi:l0.ll86/1471-2172-9-67 v

This article is available from: http://www.biomedcentral.eom/l47l-2l72/9/67

© 2008 Bradburne et al; licensee BioMed Central Ltd.

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 work is properly cited.


Background: Bacillus anthracis secretes several virulence factors targeting different host organs and cell types during inhalational anthrax infection. The bacterial expression of a key virulence factor, lethal toxin (LeTx) is closely tied to another factor, edema toxin (EdTx). Both are transcribed on the same virulence plasmid (pXOl) and both have been the subject of much individual study. Their combined effect during virulent anthrax likely modulates both the global transcriptional and the phenotypic response of macrophages and phagocytes. In fact, responses brought about by the toxins may be different than each of their individual effects.

Results: Here we report the transcriptional and apoptotic responses of the macrophage-like phagocytic cell line THP-l exposed to B. anthracis Sterne (pXOl+) spores, and B. anthracis A Sterne (pXOl-) spores. These cells are resistant to LeTx-induced cytolysis, a phenotype seen in macrophages from several mouse strains which are sensitive to toxigenic anthrax infection. Our results indicate that the pXOl-containing strain induces higher pro-inflammatory transcriptional responses during the first 4 hours of interaction with bacterium, evident in the upregulation of several genes relevant to Nf-KB, phosphatases, prostaglandins, and TNF-a, along with decreases in expression levels of genes for mitochondrial components. Both bacterial strains induce apoptosis, but in the toxigenic strain-challenged cells, apoptosis is delayed.

Conclusion: This delay in apoptosis occurs despite the much higher level of TNF-a secretion induced by the toxigenic-strain challenge. Interestingly, CFLAR, an important apoptotic inhibitor which blocks apoptosis induced by large amounts of extracellular TNF-a, is upregulated significantly during toxigenic-strain infection, but not at all during non-toxigenic-strain infection, indicating that it may play a role in blocking or delaying TNF-a-mediated apoptosis. The suppression of apoptosis by the toxigenic anthrax strain is consistent with the notion that apoptosis itself may represent a protective host cell response.


Highly pathogenic strains of Bacillus anthracis contain two plasmids, XO1 and XO2 encoding major virulence factors of this Gram-positive bacillus. Lethal toxin (LeTx) and edema toxin (EdTx) genes reside on the pXO1, while the anti-phagocytic capsule is encoded by pXO2. LeTx is necessary for pathogenicity, as deletion of its gene renders the microbe avirulent, while EdTx-knockout strains are only partially attenuated [1]. Capsule substantially contributes to the virulence of the microbe but unencapsulated strains, such as Sterne (34F2), are still capable of causing death in experimental animals [1,2] Therefore, the pXO1+, pXO2- Sterne strain serves as a convenient experimental toxigenic model of highly virulent strains.

During inhalational exposure to B. anthracis, the spores may enter alveoli where they become deposited on mucosal surfaces. Within hours of the initial interactions with the host, the spores can be engulfed by phagocytes, such as monocyte-derived macrophages or dendritic cells [3,4]. Both LeTx and EdTx are expressed upon germination within the macrophage phagosome [3] and seem to play important roles in suppressing the bactericidal innate immune mechanisms of the epithelium and the intra-phagocytic environment [3,5-9]. In a currently accepted model of anthrax, some of the phagocytosed spores en route to the mediastinal lymph nodes survive and multiply within the phagolysosome, kill the cell, and become released into the lymphatic system [5]. In the following process of hemorrhagic lymph node destruction, the bacteria gain access to the bloodstream and quickly become systemic by spreading to the spleen and other internal organs[10,11]. According to this mechanism, lung phagocytes such as macrophages and dendritic cells are critically involved in the initiation of the disease, and their response to anthrax spores, among other factors, determines whether the exposure to aerosolized spores results in the infectious process [12].

Macrophages were the first cell type discovered to die after exposure to LeTx [13]. LeTx consists of a heptameric protective antigen (PA) noncovalently associated with lethal factor (LeF). LeTx is a zinc metalloprotease, which cleaves and thus inhibits mitogen-activated protein kinase kinase (MAPKK) family members in vitro and in vivo, resulting in defective host cell signaling [14-16], with broad implications for the host innate and adaptive immune responses [17]. However the death of macrophages after exposure to LeTx in vitro does not correlate with the cleavage of the MAPKK substrates by LeTx, and the macrophages sensitive to LeTx reside in the strains of mice tending to be resistant to the lethal effect of LeTx [2,18]. This paradoxical effect did not get a satisfactory mechanistic explanation till it was discovered that LeTx was able to induce the process of programmed, apoptotic death in a murine macrophage

cell line RAW 264.7 [19]. It was suggested that for sensitive macrophages, and perhaps other cell types, undergoing apoptosis may serve as an in vivo sensor alerting the immune system and inducing the protective response [20]. Similar correlation between macrophage susceptibility to apoptosis and the outcome of infection in mice takes place, for example, in the case of M. tuberculosis [2123]. Generally, pneumonic macrophages appear to exhibit a broad apoptotic response during less virulent M. tuberculosis infection, while infection with more virulent strains correlates with a near complete loss of the same apoptotic indicators [23].

Several studies have found that resistance of cells to LeTx may vary depending on the conditions of growth or stimulation. Human peripheral blood monocytes with the LeTx-resistant phenotype become sensitive to LeTx upon growth media deprivation leading to a cellular stress [6]. This finding was further elaborated, when it was shown that the stress factors in vitro sensitizing murine macrophages to apoptosis include a range of the bacterial Tolllike receptor (TLR) agonists of different nature, including the pore-forming hemolysins, peptidoglycan or endo-toxin [24]. Human U-937, HL-60 and THP-1 cell lines of monocyte origin are resistant to LeTx but become sensitized upon phorbol myristate acetate (PMA)-induced differentiation into macrophages in culture [25], although relevance of this stimulation to the infection process in vivo is not clear. Additionally, PMA-induction may yield significant changes in the transcriptome which do not occur in vivo[26].

In contrast to the LeTx-exposed cells, human peripheral blood monocytes and monocyte-like undifferentiated THP-1 cells readily become apoptotic after exposure to B. anthracis spores [6]. These observations suggest that the process of infection may involve sensitization of cells to LeTx through the activity of unknown factors, such as the bacterial TLR agonists, similar to the effects of these agonists in vitro[24]. Another possibility is the existence of apoptotic factor(s) independent of LeTx or working in concert with it. In order to explore these hypotheses we decided to obtain insight into the transcriptional and apoptotic responses of THP-1 cells challenged with B. anthracis spores. To identify host signaling associated with the pathogenic activity of pXO1-encoded bacterial factors we compared the effects of the toxigenic (pXO1+, pXO2-) Sterne strain and the non-toxigenic (pXO1-, pXO2-) A Sterne strain. In experimental animals the toxigenic strain, even in the absence of capsule-encoding pXO2, is able to cause a lethal systemic infection, while a non-toxigenic strain is non-lethal and causes no clinical symptoms of disease. We therefore suggested that the use of the matched pair of isogenic strains would provide a means to

distinguish between the pathogenic and the protective host responses.

We demonstrate that the presence of toxin plasmid XO1 in the challenge strain results in higher pro-inflammatory responses. Both B. anthracis strains induce apoptosis, albeit by different mechanisms. Our data supports a hypothesis that in the presence of pXO1, apoptosis may proceed through the mitochondria-dependent pathway, while in the non-toxigenic infection it seems likely to be activated mainly through an extrinsic TNF-a receptor-mediated caspase-8 pathway. The onset of apoptosis induced by a toxigenic strain is markedly delayed, and key apoptotic inhibitors are clearly induced, indicating that the suppression of early apoptosis in phagocytic cells may be beneficial for the virulent microbe. The fact that pXO1 is not required for the induction of apoptosis suggests that in addition to LeTx, B. anthracis possesses chromosoma-lly-encoded pro-apoptotic factor(s).


Growth characteristics indicate increased bacterial replication in the presence of THP-1 during early growth

A Bacillus anthracis A Sterne strain was selected which was negative for the toxigenic XO1 plasmid, identical in genomic variable sequences checked, and which had identical growth characteristics to its parent Sterne strain (Additional files 1 and 2). THP-1 cells were challenged with spores of the Sterne strain and A Sterne strains, the bacterial growth monitored, and THP-1 mRNA was

extracted at early time points to be used for microarray analysis. Figure 1 shows the bacterial growth represented by simple absorbance at 2 to 8 h post-coincubation with THP-1. In the non-toxigenic challenge, there is no statistically significant difference between challenge with or without THP-1 other than the MOI-dependent difference. This is in contrast to the Sterne-challenge, which indicates more bacterial replication in the presence of THP-1 at both MOIs.

Principal component analysis indicates a suppressed response to toxigenic strains during early (2 h) infection

THP-1 mRNA was collected at 0, 2, and 4 h time points, and the transcriptional response determined by microar-rays. Principal component analysis of the microarray results shows time- and dose-dependent responses by THP-1 cells during challenge by both strains of B. anthracis (Fig. 2A). Variation attributed to each component is also shown (Fig 2B). Data points clustered in the PCA and confirmed that the principal components of the data set are the experimental variables such as MOI, post exposure time and the nature of the challenge strain. Therefore the mRNA expression levels reflect changes depending on the differences in treatment conditions rather than experimental artifacts. As expected, PCA also closely grouped the results with respect to chip replicates indicating their consistency, and allowed us to evaluate the overall dynamics of cell responses (Fig. 2). The results of the 2hour challenges by toxigenic Sterne strain spores at both MOIs cluster in the 3-dimensional parameter space close


2 3 4 5 6 7 8

Coincubation time with or without THP-1 (hours)

dS MOI 1 THP dS MOI 10 THP dS MOI 10 dS MOI 1

2 3 4 5 6 7 8

Coincubation time with or without THP-1 (hours)


Incubation of A) B. anthracis Sterne, or B) B. anthracis A Sterne, in serum-free media with or without ! * I06 THP-! cells, and at multiplicities of infection (MOI) of either !:! or 1:10 of THP-1 to bacterium.


A) Principle component analysis of microarrays run in this study confirms groupings of similar treatments, time points, and experimental replicates. Each color corresponds to the same experimental treatment, while each point represents a separate chip. Axes labeled PCl, PC2, and PC3 and represent Principle Component l, Principle Component 2, and Principle Component 3, respectively. B) Variance of principle components evaluated.

to the control samples within the same quadrant of our PCA, in contrast to the cell responses to 2-hour challenges by A Sterne spores. This demonstrates that THP-1 cells respond to the non-toxigenic infection faster, compared to the toxigenic one, in agreement with the hypothesis that the initial general response of macrophages to the virulent spores is suppressed due to the presence of the LeTx-expressing plasmid [17]. At the 4-hour time point, there are notable and discernable differences between Sterne and A Sterne strain-challenged cells indicating a clear distinction in cellular responses between the strains with regard to the presence of the toxigenic plasmid. The strain-specific separation at higher MOIs and time points is so distinct that it may be useful in delineating toxigenic vs. non-toxigenic strains. Further confidence in our gene chip data was obtained using Real-Time PCR (Fig. 3), which demonstrates the results overall correlating to the micro-array data within the experimental variations of both data sets.

Gene ontology comparisons yield known and unknown differences

Fig. 4 shows clustered heatmaps of important biological processes, molecular functions, and cellular component modules. In general, our results for both strains appear consistent with the macrophage and innate immune activation programs described in response to other bacterial pathogens [27,28]. Activation of innate and other immune defenses and inflammatory responses are typically observed for both spore strains, but in general the activation is greater in the case of the more virulent Sterne strain spore-challenge model (Fig. 4a) in spite of the overall initial delay noticed in the 2-hour challenges. Several of the NF-KB-controlled programs are upregulated similarly between both strains, indicating that they are pXO1-independent, while the others could be considered as characteristics of virulence. The most striking differences seen between strain challenges in our clustering analyses consist of a series of modules that are activated at the 4hour time points in the Sterne spore treatments, while no activation is seen in any of the A Sterne spore treatments. Among these modules, there are important signaling pathway processes, including at least 10 phosphatase modules, DNA modification modules, and transcrip-tional regulation modules. In addition, the pathways involved in eicosanoid and prostaglandin-type metabolism are similarly affected, which agrees with reports on drastic changes in prostaglandin levels during toxigenic-anthrax pathogenesis in mice [29]. Other cellular component functions encompassing catabolic and metabolic components do not show much variation between strains, similar to nucleotide/nucleic acid and ATP metabolism and catabolic processes (Fig 4b, c). Interestingly though, the cytochrome oxidase and oxidoreductase modules appear highly activated early in the avirulent A Sterne strain challenge, in comparison to the Sterne strain, mark-

ing the suppression of the transcriptional-level mitochondrial activation as an important virulence feature. Previous reports have indeed implicated LeTx in the dysfunction of mitochondria in murine macrophages and human peripheral blood monocytes [6,30].

Anti-apoptotic and stress-responses are favored during toxigenic challenge, while cell cycle/differentiation responses are favored during non-toxigenic challenge

To further examine specific differences in host response to each strain, an Analysis of Variance (ANOVA) was performed and genes found to be significant in the cell survival, stress response, and apoptosis were compared between different strain challenges under identical conditions (Tables 1, 2, 3). The use of ANOVA here allowed us to evaluate the differential expression of a particular human gene while varying 3 conditions: time, MOI, and the presence or absence of pXOl in the infecting B. anthracis strain. In our analysis, positive numbers indicate that the expression levels are higher in the non-toxigenic pXO 1 (-) spore treatment, while negative numbers reflect a shift in expression level favoring the more virulent toxigenic pXOl (+) spore treatment. A complete list of genes evaluated in the ANOVA can be found in Additional file 3 of the Additional Materials.

The ANOVA analysis shows a marked upregulation for several stress and inflammatory response genes relevant to the TNF-a biological module in the case of toxigenic Sterne spore challenge relative to the non-toxigenic A Sterne one (Table 1, 3). Contrasted against the apoptosis-related genes are the cell cycle genes, which seem to be affected in response to the non-toxigenic challenge (Table 2), indicating an initial shift towards proliferation and differentiation. CXCL2, the chemokine precursor to macrophage inflammatory protein-2a (MIP-2a) is strongly increased in Sterne spore-challenged cells, pointing to a higher innate response to the toxigenic strain. The apop-tosis inhibitor genes for cFLAR (CASP8/FADD-Like Apoptosis Regulator) and GADD45b (Growth Arrest and DNA-Damage-inducible, beta) also appear to be expressed higher in the toxigenic spore treatments (Table 1 and Fig. 5). The level of cFLAR protein determined by Western blot follows the same trend (Fig. 6a). The increased anti-apop-totic responses to the toxigenic strain challenge correlate with the induction of higher transcription of Nf-KBl (Table 1), a gene which may have some anti-apoptotic effects[ 31,32]. In addition the Nf-кВ inhibitor QNZ reduces cFLAR expression (Fig. 6b) in agreement with its higher level in the toxigenic strain infection (Fig. 6a).

TNF-a-dependent apoptosis during toxigenic infection may be countered by cFLAR-and NF-kB-dependent anti-apoptotic responses

The role of TNF-a, the primary mediator of the cellular inflammatory response, has attracted significant attention

Microarray Data Real-Time PCR

3.5 -,

<n 2.5

X HI 1.5

■o 1

u. 0.5

<d -0.5


4 h Sterne, 1:10 MOI


4 h & Sterne, 1:10 MOI

X „ m 2

4 h Sterne, 1:10 MOI


1.5 = 1

Ù5 -0.5

? -1 u_

-1.5 -2



<u 3 a.


2 h Sterne, 1:10 MOI

= 6 <2

"<n 5 -<n

£ 4 -ü

2 h Sterne, 1:10 MOI


c o m m a> 1.5

a. X HI 1

TD O U. 0.5

4 h Sterne, 1:10 MOI


c o 20

<n m <u 15

a. X UJ 10

o u. 5


4 h Sterne 1:10 MOI


Real-Time PCR of single gene treatments confirms trends in microarray expression levels. Y-axis in microarray data (left) represents the fold change in gene expression relative to untreated control samples (labeled control) of THP-1 cells. The Y-axis in Real-Time PCR data (right) represents the fold gene expression in the THP-1 cells within the sample relative to the level of endogenous 18-S gene expression, and was calculated using the 2-AACT method described by Livak and Schmitt-gen[69]. Labels within each panel indicate the time post- exposure, the challenge strain and the MOI.

A. Apoptosis and Immune Responses

B. Cell cycle

rteme Д^ете

^^ Нш„


BmiiSi""*" "ll-1 1L11'"-1 "*""*""

С. Stress response

Figure 4

Heatmaps of GO-type biological response modules show both similarities and differences in monocytic responses to Sterne and A Sterne B. anthracis strains. MOIs correspond to a 1:1 at the low end of a triangle, and a 1:10 at the high end of each triangle. Color ranges represent no activation (black), low activation (faint red), and extensive activation (full red).

Table 1: Log ratios of genes of interest in ANOVA analysis* for apoptotic genes



2 h 1:1 MOI 2 h 1:10 MOI 4 h 1:1 MOI 4 h 1:10 MOI

NM_ _00l 196 BID

NM_ _0l5675 GADD45B

XM_ 046543 MKPX

XM_ .010767 NCKAPl

NM_ _006107 OA48-I8

NM _004l80 TANK

NM _003879 CFLAR

XM_ .028204 NFKBI

XM_ .041847 TNF

XM_ .012894 ZNFI4

NM .014977 ACINUS

BH3 interacting domain death agonist growth arrest and DNA-damage-inducible, beta mitogen-activated protein kinase phosphatase x NCK-associated protein 1 acid-inducible phosphoprotein TRAF family member-associated NFKB activator CASP8 and FADD-like apoptosis regulator nuclear factor of kappa light polypeptide gene enhancer in B-cells 1 (p 105) tumor necrosis factor (TNF superfamily, member 2) zinc finger protein 14 (KOX 6) apoptotic chromatin condensation inducer in the nucleus

-0.0317 -0.0750 -0.6127 -0.57l7

-0.0800 -0.5800 -2.0967 -2.0000

0.1183 0.5667 0.8563 0.8250

0.1750 0.2333 0.4737 0.2750

0.0433 0.0333 0.5737 0.7933

0.2383 -0.2633 -0.7193 -0.6650

-0.0583 -0.1283 -0.7720 -1.1250

-0.1467 0.0283 -0.7677 -1.3933

-0.1683 -0.7250 -2.4107 -1.3967

0.3217 0.2067 0.2593 0.1017

0.2867 0.5983 0.0723 -0.5367

*All adjusted P-values < 0.05. Ratios represent comparisons of expression levels of non-toxigenic vs. toxigenic-challenge. (A more comprehensive list of genes is shown in the additional materials, Additional file 3).

in anthrax studies. Recent data indicate that its systemic levels during infection, or upon LeTx challenge in animals and cell cultures, are either undetectable or transient [12,20,33] and therefore cannot be considered as a primary cause of death, as it was previously suggested [34,35]. Nevertheless, some data suggest that the TNF-a response may play an important function in the local host antibacterial defense [36,37]. In our gene chip and RealTime PCR experiments, TNF-a gene expression is significantly induced relative to uninfected cells in a time- and concentration-dependent manner during infection by both strains (Fig. 5), however the induction by the toxi-genic strain is approximately 10 times higher, compared to the non-toxigenic one. This induction correlates with the increased protein levels of TNF-a determined by ELISA in the culture medium of cells after spore challenge (Fig. 7). The level of TNF-a is almost 3-fold higher (p < 0.01) in the cells challenged with the toxigenic spores after 4 h.

This large cellular TNF-a response to toxigenic spores compared to a lower one in the nontoxigenic spore-chal-

lenge is an interesting observation, suggesting that the issue of cytokine suppression during anthrax is not straightforward. It is well known that pXO1 modulates the expression of many chromosomal proteins in addition to its own, resulting in larger numbers of unique secretome proteins in pXO1 (+) strains compared to pXO1(-) strains[38]. This may contribute to the increased response to the plasmid (+) strain seen in our study. However LeTx, a well-known protease of the mitogen-activated protein kinase kinases (MEKs), would be expected to inhibit expression of cytokines, and in THP-1 cells, it has been shown to inhibit the release of TNF-a following induction by an agonist[39]. In our experiments, the expression of toxins takes place approximately 2 to 4 hours post-spore challenge, allowing time for a general innate response before cytokine suppression. Our study also agrees with a previous report from Pickering et. al., who found no inhibition of cytokine response by a toxi-genic strain, and an even lower response from a non-toxi-genic strain [12], while other studies, such as those by Tournier et. al., and Drysdale et. al. have demonstrated the

Table 2: Log ratios of genes of interest in ANOVA analysis* for cell cycle genes

Accession Name Description 2 h 1:1 MOI 2 h 1:10 MOI 4 h 1:1 MOI 4 h 1:10 MOI

XM_ 034567 CCND2 cyclin D2 0.I833 0.3083 0.5300 0.2267

NM 033379 CDC2 cell division cycle 2, Gl to S and G2 to M 0.0250 0.I783 0.4483 0.5400

NM 002357 MAD MAX dimerization protein 1 0.0I50 -0.0433 -0.8050 -0.7250

XM_ 037657 MYC v-myc myelocytomatosis viral oncogene homolog (avian) 0.I233 0.I783 0.6903 0.5983

NM 004725 BUB3 BUB3 budding uninhib. by benzimidazoles 3 homolog (yeast) -0.0500 0.I467 0.49I0 0.5767

XM_ 034725 NRPl neuropilin l -0.03I7 0.I583 0.2967 0.5467

XM_ 008679 TP53 tumor protein p53 0.6I83 0.2I67 0.68I0 -0.I967

NM 005620 SI00AII Sl00 calcium binding protein Al 1 (calgizzarin) 0.I667 0.2750 0.4893 0.2567

NM 004705 PRKRIR protein-kinase, interferon-inducible double stranded RNA dependent inhibitor, (P58 repressor) 0.0033 -0.0567 0.0303 0.08I7

*All adjusted P-values < 0.05 Ratios represent comparisons of expression levels of non-toxigenic vs. toxigenic-challenge.

Table 3: Log ratios of genes of interest in ANOVA analysis* for stress-response genes

Accession Name Description 2 h 1:1 MOI 2 h 1:10 MOI 4 h 1:1 MOI 4 h 1:10 MOI

NM_ _004040 ARHB ras homolog gene family, B -0.I783 0.0233 -I.7523 -I.5083

NM_ _002089 CXCL2 chemokine (C-X-C motif) ligand 2 -0.2850 -0.9700 -2.3083 -I.75I7

NM_ _00I4I6 EIF4AI eukaryotic translation initiation factor 4A, isoform I 0.I067 0.3233 0.5II3 0.5733

NM_ _005239 ETS2 v-ets erythroblastosis virus E26 oncogene homolog 2 (avian) -0.0417 -0.I3I7 -0.7237 -0.6I67

NM_ _006597 HSPA8 heat shock 70 kDa protein 8 -0.3I67 -0.3550 -0.5897 -0.2633

XM_ .058230 JUND jun D proto-oncogene -0.0883 -0.0I33 -0.9497 -0.7983

NM _005566 LDHA lactate dehydrogenase A -0.0I50 0.6I50 0.8707 I.I967

XM_ _03I8I8 LYN v-yes-I Yamaguchi sarcoma viral related oncogene homolog 0.0417 -0.2733 -0.6830 -I.06I7

XM_ 008855 NR2F6 nuclear receptor subfamily 2, group F, member 6 0.0I83 0.I483 0.3897 0.4650

XM_ .007258 TNFAIP2 tumor necrosis factor, alpha-induced protein 2 -0.0383 -0.2517 -I.0440 -I.5067

XM_ .002762 TNFAIP6 tumor necrosis factor, alpha-induced protein 6 -0.0350 -0.I5I7 -I.5390 -I.0600

NM 003295 TPTI tumor protein, translationally-controlled I -0.4950 -0.5I00 -0.3893 -0.4267

* All adjusted P-values < 0.05. Ratios represent comparisons of expression levels of non-toxigenic vs. toxigenic-challenge.

more expected suppression of TNF-a and cytokine responses upon toxigenic B. anthracis challenge [40,41]. These discrepancies may in part be explained by differences in the timing mentioned above, but also by the experimental conditions. Each of these studies attempts to remove extracellular bacteria and/or excreted proteases/ factors from their models in some way, therefore excluding their effect. In our conditions, we were interested in leaving these factors to allow for the evaluation of these extracellular components in the host cell response of this LeTx-cytolysis-resistant cell line.

Elevated TNF-a inflammatory response of the THP-1 cells to the pXO1 gene products is expected to enhance the protective, bactericidal capacity of macrophages and sensitize them to apoptosis [42,43]. It is therefore possible that B. anthracis has developed a mechanism to counter these macrophage defenses during the initial infectious activity. In sensitive macrophages TNF-a is expected to substantially contribute to the apoptotic death program through a TNF Receptor 1 (TNFRl)-dependent pathway including activation of caspase-8. An interesting feature of the TNF signaling network is the existence of extensive cross talk between the apoptosis and NF-kB signaling pathways that emanate from TNFR1. While in the absence of NF-kB activity the cellular susceptibility to TNF-a - induced apoptosis increases, activation of NF-kB protects against apoptosis[44,45]. Therefore, the increased TNF-a and NF-kB expression induced by the toxigenic strain may result either in the anti-apoptotic effect or in the delay of apop-tosis onset depending on the balance between the above stimuli in the TNF-related apoptotic pathway. In agreement with this hypothesis, our data show that apoptosis detected by TUNEL assay is delayed after 2 h post exposure to the toxigenic strain while after 6 h it takes place to a similar extent in the case of both strains (Fig. 8).


Our results indicate that the THP-1 monocytes present significant innate immune responses evident in activation of several biological modules during both toxigenic, and non-toxigenic challenge. In general, the early response is suppressed by the toxigenic strain but overall this strain demonstrates stronger activation for several stress, inflammatory and apoptosis-related genes relevant to the TNF-a biological module. Observed genotypic responses specific for the pXO1-encoded pathogenic factors are consistent with mitochondrial damage, prostaglandin level disruption, and phosphatase induction likely inhibiting phos-phoprotein cell signaling. These results agree with the transcriptional analysis of the murine macrophage RAW 264.7 cells in response to B. anthracis Sterne spore challenge [46]. The authors reported transcriptional activation of immune modules, induction of Nf-KB, along with the apoptotic inhibitor GADD45b. However in the case of un-stimulated macrophages treated with LeTx there were no significant changes observed. Analyzing the effect of LeTx on Salmonella lipopolysaccharide-stimulated RAW 264.7 cells, Tucker et al. detected suppression of the host immune response [47] however the antigen (LPS) used for stimulation is irrelevant to anthrax infection. According to our results, the plasmid-negative strain provides the background information, which is important for precise evaluation of the pXO1 contribution to the host cell response, such as silencing of protective signaling in response to infection and/or changes in signaling mechanisms.

We present evidence that both toxigenic and non-toxi-genic B. anthracis strains induce apoptosis, which can be considered as a form of innate immunity against some pathogens. Our data indicate that in the absence of toxi-genic pXO1, the infection with the avirulent strain pro-


Fold change of TNF-a and cFLAR expression in THP-1 cells as measured by Real-Time PCR. Left panels, B. anthracis Sterne infection. Right panels, B. anthracis A Sterne infection. Notice the difference in the Y scale between TNF-a panels.

250 200 150 100 50 0

TNF- a / sterne

■ 1:1 MOI □ 1:10 MOI

0 Min 120 min 240 min

25 20 15 10 5 0 -5

CFLAR / sterne

■ 1:1 MOI □ 1:10 MOI

0 Min 120 min 240 min

TNF- a / A sterne

■ 1:1 MOI □ 1:10 MOI

0 Min 120 min 240 min

CFLAR / A sterne

20 15 10 5 0 -5

■ 1:1 MOI □ 1:10 MOI

0 Min 120 min 240 min

ceeds without a substantial intrinsic insult to the cell. Therefore in the absence of cFLAR induction, and in conjunction with the observed TNF-a production, we hypothesize that infection in the absence of pXO1 probably induces apoptosis through a classical extrinsic TNF-a receptor-mediated caspase-8 pathway. In contrast, the pXO1-encoded complex of secreted toxins acquired by B. anthracis as an instrument of intracellular immune evasion causes a strong intracellular insult. The large induction of CFLAR likely plays some role in inhibiting TNFR-mediated apoptosis, and may be involved in its observed phenotypical delay, ultimately pushing the system toward

the intrinsic, mitochondria-dependent pathway of apoptosis some time later [48]. This pathway proceeds through mitochondrial dysfunction, and is a primary differential observation of our gene ontology modules between the strains. (Fig. 4), as indicated by the perturbation of modules containing mitochondrial components in the Sterne-infected cells. Future experimentation such as the TNFR1-dependent processing of caspase-8 vs. caspase-9, mito-chondrial membrane polarization/permeability, and the formation of apoptotic signaling complexes during challenge by each of these strains will help to determine if this is indeed the case.

Figure 6

CFLAR protein induction in THP-1 cells. A. Western blots of CFLAR induction at MOIs 1:1 and 1:10 at 4 h post infection. B. Inhibition of NF-kB with the inhibitor QNZ during challenge suppresses the production of cFLAR. Samples were taken at 8 h post infection at MOI 1:1. C denotes unchallenged control; S, challenge with Sterne strain; and AS, challenge with A Sterne strain. Actin is shown as a control.

Interestingly, the onset of apoptosis induced by a toxigenic strain is markedly delayed, and bacterial growth is significantly increased, indicating existence of a mechanism modulating apoptosis in favor of the microbe through the increased production of anti-apoptotic cFLAR and Nf-KB. This delay of apoptosis during anthrax infection could allow for more successful bacterial survival and dissemination, resulting in systemic disease. This supports a hypothesis that in anthrax the suppression of apoptosis in macrophage-like cells may contribute to disease virulence during the initial stages of infection, and that apop-tosis itself may represent a protective host cell response until it reaches a pathologic proportion.

Several lines of evidence in the literature support this hypothesis. Previous analyses of anthrax-infected RAW 264.7 cells [46] identified overexpression of ornithine decarboxylase (ODC), a normal biosynthetic enzyme involved in the conversion of putrescine to the polyamines spermine and spermidine. Since the overproduction of polyamines by ODC has been implicated in preventing apoptosis [49], this led the authors to speculate that ODC overexpression was involved in suppression of macrophage apoptosis. In the case of other bacterial

infections, apoptosis has been generally considered as a protective host response [50]. There is additional evidence that the delay of apoptosis may play roles in some viral and fungal infections as well [51-53].

Macrophages infected with viable virulent strains of Mycobacterium tuberculosis undergo apoptosis with far less frequency than macrophages infected with attenuated strains [21,22]. Strain virulence in mice correlates directly to the extent to which apoptosis is inhibited [23]. Thus virulent strains of M. tuberculosis apparently suppress apoptosis in host cells. Conversely, macrophage apoptosis appears to provide the host with a means for killing intracellular microorganisms, with host cells likely initiating apoptosis following initial infection [54,55]. Consistent with this, Mycobacterium Avium, a less virulent pathogen that is less adapted to the intra-macrophage growth species compared to M. tuberculosis does not suppress macrophage apoptosis [56].

Mycobacterium infections have been indicated to subvert pro-apoptotic TNF-a signaling by either reduction of TNF-a expression or downregulation of its receptor [57], as well as differential NF-kB activation [58]. It was found that both virulent and avirulent M. tuberculosis strains activated NF-kB after 4 h in THP-1 cells however after 48 h only the virulent strain maintained NF-kB activation, which lead to up-regulation of a bcl-2 family anti-apop-totic member, bfl-1/A1. These results indicate that NF-kB activation may be a determinant factor for the success of virulent mycobacteria within macrophages. L. pneu-mophila induces a robust activation of caspase-3 in alveolar macrophages however the apoptotic cell death is not executed until late stages of the infection, concomitant with the termination of intracellular replication [59].

Recent studies of hepatitis C infection directly implicate the sustained upregulation of the cFLAR gene product in preventing TNF-a-mediated apoptosis and contributing to disease initiation by the virus core protein [52]. The induction of cFLAR during toxigenic anthrax infection likely contributes to prevent apoptosis in a similar fashion, particularly when combined with the suppression of TNF-a production, as it has been shown that LeTx inhibited the TNF-a production in both endotoxin-stimulated macrophages [60,61] and anthrax bacterial cell wall-stimulated peripheral blood mononuclear cells [6]. A mechanism consistent with our observations has been suggested for THP-1 cells stimulated through TLR-4 with the pore-forming hemolysin, anthrolysin O (AnlO) [54]. In this report the NF-KB-dependent activation of Bfl-1/A1 gene by AnlO contributed to anti-apoptotic signaling, although overall the TLR-4 stimulation resulted in apoptosis. However the role of AnlO in vivo remains to be elucidated. Another possible pathway of apoptosis inhibition in mac-

11800 9800 E 7800 5800

та аз

0 2 4 6 Infection time (h)


Sterne(MOI, 1:1) Sterne(MOI, 1:10) dSterne(MOI, 1:1) dSterne(MOI, 1:10)


Production of TNF-a protein mirrors transcriptional results. ELISA of TNF-a produced by THP-1 cells challenged by Sterne and A Sterne over 8 h.

rophages includes EdTx-mediated activation of the anti-apoptotic Platelet Activator Inhibitor type 2 gene (PAI-2), although THP-1 cells carry an inactivating mutation in this gene [46].

The previous findings that anthrax LeTx inhibits the p38 signaling pathway by cleavage of MAPKKs in macrophages and dendritic cells [60,62] are consistent with the protective role of this stress-activated pathway in elimination of intracellular pathogens and activation of host immunity [63]. From this standpoint, the proteolytic activity of LeTx on the MAPK activation pathway may contribute to establishing infection by preventing cell activation in response to infection and thereby increasing the bacterial survival within phagocytes [19,60].


Using a matched pair of the pXO1 (+) and pXO1 (-) B. anthracis strains, we demonstrate that the monocyte-type

THP cells in vitro generate a distinct set of transcriptional responses to infection characteristic of the virulence factors contributed by the XO1 plasmid. Among these responses, the delay of apoptosis associated with the elevated expression of anti-apoptotic cFLAR represents a likely pathogenic strategy directed toward bacterial survival in the conditions of the increased host inflammatory response. This may be particularly relevant for hosts carrying mononuclear cells which, like THP-1, are resistant to LeTx-induced cytolysis. In the absence of pXO1, the non-toxigenic strain induces early apoptosis consistent with its protective role in several infectious diseases as a form of the host innate immune response. This finding implies that in addition to LeTx, virulent B. anthracis possesses a chromosome-encoded factor(s) triggering a pro-apoptotic host response, as well as the plasmid factors counterbalancing apoptosis.




! ". 1 ** *. . f/s ' • • • ■1 • • s * ' .. * •• * . *> * rJ* •*■ A m* * V * * ■ • » » .î • . v -, . . • " r • j 1.0

/1 ' "* • .* * 11» *. * • •- • 1 "¡." • ' • Äs, - • * * fc V* / , * * V • * 5.0 t # • •

*" • 4 f<>, * :*.V-- J'' *- . t\ ' ' ' ' • ■ • ' *# • Î . . • V i- • ' . • * »' • V- • ■ • ■ -t. ■ 0.8

' . \ i • / 9 ' ' " " * . V- * V • . • , • *.*'"•* " * • 1.0

• • t • ** " % h i t • : X'' '••':. ¿ v , * • •• * » j ***** • : * * M 5.1 •

». ' \ • ' . m , ! *s> • • ♦ • ^ fV » • . V 'i. • ' v . » • '. î **. • V ''j • » ' « » * " v. ' A' * 5.9 • • •

Figure 8

Apoptosis of THP-1 cells is delayed at 2 h post challenge by the toxigenic Sterne strain relative to the non-tox-igenic A Sterne strain. Apoptosis was detected with TUNEL staining. The numbers on the images show counts of objects with the intensities above background relative to corresponding untreated controls and indicate the extent of apoptosis. DAPI staining shows that equal number of cells was used in all experiments.


Bacterial strains and THP-1 cell cultures

Bacillus anthracis strain 34F2 Sterne (pXO1+, pXO2-) was obtained from Colorado Serum Company (CO). A Sterne strain, from the Collection of the National Center for Bio-defense and Infectious Diseases (George Mason University, Manassas, VA), is a plasmidless derivative of Sterne strain generated by curing of the pXO 1 plasmid by growth for 6 days in LB media at 42.5°C [64]. In our experiments, the absence of pXO1 in the cured strain was confirmed by PCR with the protective antigen-specific primers. B. anthracis is one of the most monophyletic bacterial strains known, however the introduction of genomic changes by the curing process cannot be discounted. Despite the uniformity of its genome, B. anthracis strains do contain sequences which exhibit variability in response to selective pressures and allow differentiation. As it is not possible to test all genomic sequence without considerable time and expense, we therefore applied a high-resolution multi locus sequence typing using the technique and primers of Helgeson [65] to ensure that both strains were at least identical for all variable sequences tested (Additional file 2) Growth curves for both strains were also

identical in 1% fetal bovine serum (FBS) cell culture medium, with 4 h corresponding to the mid point of the exponential phase of growth (Additional file). Spores were generated by growth for 4 days on LB agar plates followed by re-suspension in water and pelleting by centrif-ugation two to four times, then storage in water at 4°C until use.

THP-1 was obtained from ATCC (Manassas, VA) and grown in complete medium (RPMI 1640 with glutamine, 25 mM HEPES, 4.5 g/l glucose, and 10% FBS) for propagation, until 24 h before challenge with spores, at which point THP-1 cells were washed and re-suspended for 24 h in serum-free (low protein) Cellgro media (Media Tech, Herndon, Va.) at a concentration of 1.0 x 106 cells/ml. Cell cultures were challenged for 2 h or 4 h at a cell to spore multiplicity of infection (MOI) of either 1:1 or 1:10 for microarray experiments. The term MOI used here simply refers to the number of spores added per cell. Each timepoint/MOI combination was replicated at least 3 times for OD readings and microarray analysis. Viability was measured by a standard Trypan Blue permeability assay by counting dead cells with a hemocytometer. Infec-

tion of THP-1 monocytes by both strains for 4 h in the conditions of our experiments resulted in a viability of greater than 90% relative to untreated cells.

Microarrays and ANOVA analysis

RNA was collected and purfied immediately at the appo-priate timepoints by Trizol (Invitrogen Corporation, Carlsbad, CA), and quality determined by checking on a Lab-on-a-Chip Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA). Samples were similarly quality checked after both probe labelling and hybridizations steps. The microarray platform was developed, printed, hybridized and analyzed at Avalon Pharmaceuticals, Inc. (German-town, MD). Probe design and printing were all done using automated, proprietary processes (Avalon Pharmacueti-cals, Germaintown, Md) and contained quality control steps after each. Arrays consist of proprietary 80-mer long, spotted oligonucleotides representing 2,000 human genes. Each chip contains four replicates of each gene (8,000 total/chip) with 64 perfect match and 64 mismatch controls. Arrays were hybridized using standard Cy3/Cy5 sample/reference labeling and analyzed using a ScanArray reader (Perkin-Elmer, Wellesly, MA). Several sets of genes related to apoptosis, stress response, and cell differentiation were chosen for further evaluation. To determine whether these genes were differentially expressed between two experimental states, an ANOVA was performed on the log-transformed expression ratios as described in [66]. The following model was used for this analysis: yijkl = Ei + Mj + Tk + (EMT)ijkl + eijkl where yijkl denotes the log2 expression ratio measured for experimental state i, MOI j, time state k, replicate l, with 1 < i < 3, 1 < j < 2, 1 < k < 2, and 1 < l < 6. The term Ei measures the effect of the infection; Mj measures the effect of the MOI; Tk measures the effect of the time state, and the term (EMT)jkl measures the interaction effect of all three variables. An ANOVA was performed on each gene using the linear model above. Twelve contrasts were based on pair-wise differences between each experimental condition and the control state, as well as the experimental state which differed in the MOIs only. The R package limma was used for ANOVA methods [67]. A multiple testing correction [68] controlled the false discovery rate, which is the expected proportion of false positives in the rejected hypotheses. Genes with adjusted F-statistic ^-values < 0.05 were collected for further inspection.

Real-Time PCR

To confirm microarray data, Real-Time PCR was performed using the iQ SYBR Green kit from BIO-RAD (Hercules, CA) using the manufacturers specifications. Real time data for each gene represents the average of 3 biological replicates, each replicated 3 times for a total of 9 reactions per gene/treatment. Fold expression was calculated for each gene relative to an endogenous 18-S control gene using the 2-aact method as described[69]. Primers chosen

for Real-Time analysis were as follows: The human 18-S ribosomal gene (18S), 5'AGGAATTCCCAGTAAGTGCG and 5' GCCTCACTAAACCATCCAA; Tumor necrosis factor-a (TNF superfamily, member 2; TNF-a), 5'TGTCT-CAACCCCGCATCG and 5'AGGAACAGCCACCAATAAG C; Ras-related C3 botulinum toxin substrate 2 (Rho family, small GTP-binding protein Rac2; RAC2), 5'CCTCAT-CATCCAGTCCAACG and 5' GGCGAACTCCTGCTCCTC; Chemokine (C-X-C motif) ligand 2 (CXCL2), 5' TGTCT-CAACCCCGCATCG and 5' AGGAACAGCCACCAATAAG C; CASP8/FADD-like apoptosis regulator (cFLAR), 5'AGA TCACCAGCCTTACCG and 5' GATGAGCCACAGGAAAT GC. Tms used in Real-Time cyclers were 54.3°C for TNF-a, cFLAR, and 18-S genes; 52.9°C for the RAC gene; and 63.2°C for the CXCL2 gene.

Representation of biological response modules

Expression values of genes from each of three replicate chips for each treatment were averaged, and genes corresponding to specific gene ontology (GO) categorizations (modules) http://www. godataba se. org. were clustered using 'Cluster'

ware.htm[34,70]. Clustering was performed without additional filtering or adjustment, using the uncentered correlation coefficient, followed by complete linkage clustering.

Western Blot Analysis and ELISA

For Western blot analysis, spore-treated THP-1 cells were lysed in a buffer containing 25 mM Tris-HCl, pH 7.6, 150 mM NaCl, 1% NP-40, 1% sodium deoxycholate and 0.1% SDS (RIPA buffer, Pierce) and cell debris was removed by centrifugation. Resulting proteins were separated in a 420% gradient gel and electrophoretically transferred to a PVDF membrane. The membrane was probed with mouse monoclonal antibodies against cFLAR (Abnova), and actin (Chemicon), and then the horseradish peroxidisase-conjugated anti-mouse secondary antibody. To quantify TNF-a proteins secreted into culture medium, human TNF-a ELISA was performed according to the manufacturer's recommendation (BD Biosciences). The culture supernatants were collected by centrifugation of THP-1 cells cultured in a serum-free medium at the indicated time point post spore challenge. To investigate whether NF-kB is involved in cFLAR expression, 10 |M NF-kB inhibitor 6-amino-4-(4-phenoxyphenylethylamino) quinazoline, a cell-permeable quinazoline (QNZ, Sigma), was added to THP-1 cell culture 30 min prior to treatment with spores (MOI 1). After incubation for 8 h, cells were lysed in RIPA buffer and analyzed by Western blotting as described above.

Apoptosis Detection

For detection of apoptosis, terminal dUTP nick-end labeling (TUNEL) assay was performed according to manufacturer's recommendation (Roche Applied Science). Briefly,

the spore-challenged THP-1 cells were washed 3 times with phosphate-buffered saline (PBS), and fixed with 2% paraformaldehyde for 60 min at room temperature. The cells were washed with PBS and permeabilized with 0.1% Triton X-100 in 0.1% sodium citrate for 2 min. TUNEL reaction was carried out for 60 min at 37°C in a humidified atmosphere by adding 50 |L of reaction mixture containing terminal deoxynucleotidyl transferase. The TUNEL-positive cells were analyzed under a fluorescence microscope. The number of objects with the intensities above background on the images was determined using the Elements software (NikonNicon).

Authors' contributions

CB conceived this study, drafted the manuscript, generated the pXO1(-) strain, and performed the bacterial challenges of cell culture, qPCR, and GO clustering. MCC performed the inhibitor experiments/Western Blots, ELISAs, and TUNEL assays. QZ and DL performed all micro-array experiments and generated PCA, Variance, gene lists and GO profiles of results. KS performed ANOVA and contributed portions of the manuscript. TP grew and prepared anthrax strains and performed bacterial challenge experiments. AP performed qPCR experiments. CB and DS provided material support and encouragement for this work. SP provided material support and direction, and drafted significant portions of the manuscript.

Additional material

Additional file 1

Growth curves of B. anthracis Sterne and B. anthracis A Sterne in the conditions of THP-1 cells cultivation in serum-free medium supplemented with 1% FBS at 37°C. The data represent a comparison of growth curves for pXO1 (+), and pXO1 (-) strains. Click here for file


Additional file 2

Results of multi locus sequence typing (MLST) by method of[58]. The data represent a comparison of sequence homology between pXO1 (+), and pXO1 (-) strains, at certain variable sequence regions on the chromosome.

Click here for file


Additional file 3

Log ratios of genes evaluated in ANOVA (all adjusted P-values < 0.05). The data represent a comprehensive list of genes evaluated in ANOVA.

Click here for file



We would like to thank Meena Augustus at Avalon Pharmacueticals, for

supporting and encouraging the microarray portions of this work.


1. Pezard C, Berche P, Mock M: Contribution of individual toxin components to virulence of Bacillus anthracis. Infection and Immunity 1991, 59:3472-3477.

2. Welkos SL, Keener TJ, Gibbs PH: Differences in susceptibility of inbred mice to Bacillus anthracis. Infection and Immunity 1986, 51:795-800.

3. Guidi-Rontani C, Levy M, Ohayon H, Mock M: Fate of germinated Bacillus anthracis spores in primary murine macrophages. Molecular Microbiology 2001, 42:931 -938.

4. Brittingham KC, Ruthel G, Panchal RG, Fuller CL, Ribot WJ, Hoover TA, Young HA, Anderson AO, Bavari S: Dendritic Cells Endocy-tose Bacillus anthracis Spores: Implications for Anthrax Pathogenesis. Journal of Immunology 2005, 174:5545-5552.

5. Dixon TC, Fadl AA, Koehler TM, Swanson JA, Hanna PC: Early Bacillus anthracis-macrophage interactions: intracellular survival and escape. Cellular Microbiology 2000, 2:453-463.

6. Popov SG, Villasmil R, Bernardi J, Grene E, Cardwell J, Popova T, Wu A, Alibek D, Bailey C, Alibek K: Effect of Bacillus anthracis lethal toxin on human peripheral blood mononuclear cells. FEBS Letters 2002, 527:211-215.

7. Hu H, Sa Q, Koehler TM, Aronson AI, Zhou D: Inactivation of Bacillus anthracis spores in murine primary macrophages. Cell Microbiology 2006, 8:1634-1642.

8. Cote CK, Van Rooijen N, Welkos SL: Roles of Macrophages and Neutrophils in the Early Host Response to Bacillus anthracis Spores in a Mouse Model of Infection. Infection and Immunity 2006, 74:469-480.

9. Steffen Backert WK: Interplay of bacterial toxins with host defence: Molecular mechanisms of immunomodulatory signalling. International Journal of Medical Microbiology 2005, 295:519-530.

10. Lyons CR, Lovchik J, Hutt J, Lipscomb MF, Wang E, Heninger S, Ber-liba L, Garrison K: Murine Model of Pulmonary Anthrax: Kinetics of Dissemination, Histopathology, and Mouse Strain Susceptibility. Infection and Immunity 2004, 72:4801-4809.

11. Dixon TC, Meselson M, Guillemin J, Hanna PC: Anthrax. New England Journal of Medicine 1999, 341:815-826.

12. Pickering AK, Osorio M, Lee GM, Grippe VK, Bray M, Merkel TJ: Cytokine Response to Infection with Bacillus anthracis Spores. Infection and Immunity 2004, 72:6382-6389.

13. Friedlander AM: Macrophages are sensitive to anthrax lethal toxin through an acid-dependent process. Journal of Biological Chemistry 1986, 261:7123-7126.

14. Vitale G, Pellizzari R, Recchi C, Napolitani G, Mock M, Montecucco C: Anthrax lethal factor cleaves the N-terminus of MAPKKS and induces tyrosine/threonine phosphorylation of MAPKS in cultured macrophages. Biochemical Journal 1999, 87:288-288.

15. Duesbery NS, Webb CP, Leppla SH, Gordon VM, Klimpel KR, Cope-land TD, Ahn NG, Oskarsson MK, Fukasawa K, Paull KD, et al.: Proteolytic Inactivation of MAP-Kinase-Kinase by Anthrax Lethal Factor. Science 1998, 280:734-737.

16. Duesbery NS, Resau J, Webb CP, Koochekpour S, Koo HM, Leppla SH, Woude GF Vande: Suppression of ras-mediated transformation and inhibition of tumor growth and angiogenesis by anthrax lethal factor, a proteolytic inhibitor of multiple MEK pathways. Proceedings of the National Acadamies of Sciences, USA 2001, 98:4089-4094.

17. Guidi-Rontani C: The alveolar macrophage: the Trojan horse of Bacillus anthracis. Trends in Microbiology 2002, 10:405-409.

18. Tournier J-N, Quesnel-Hellmann A, Cleret A, Vidal DR: Contribution of toxins to the pathogenesis of inhalational anthrax. Cellular Microbiology 2007, 9:555-565.

19. VR Popov S, Bernardi J, Grene E, Cardwell J, Wu A, Alibek D, Bailey C, Alibek K: Lethal toxin of Bacillus anthracis causes apoptosis of macrophages. Biochemical and Biophysical Research Communications 2002, 293:349-355.

20. Popov SG, Popova TG, Grene E, Klotz F, Cardwell J, Bradburne C, Jama Y, Maland M, Wells J, Nalca A, et al.: Systemic cytokine response in murine anthrax. Cell Microbiology 2004, 6:225-233.

21. Danelishvili L, McGarvey J, Li Y-j, Bermudez LE: Mycobacterium tuberculosis infection causes different levels of apoptosis and necrosis in human macrophages and alveolar epithelial cells. Cellular Microbiology 2003, 5:649-660.

22. Lee J, Remold HG, leong MH, Kornfeld H: Macrophage Apoptosis in Response to High Intracellular Burden of Mycobacterium tuberculosis Is Mediated by a Novel Caspase-Independent Pathway. Journal of Immunology 2006, 176:4267-4274.

23. Ríos-Barrera Víctor A, Campos-Peña V, Aguilar-León D, Lascurain Ledesma R, Meraz-Ríos Marco A, Moreno J, Figueroa-Granados V, Hernández-Pando R: Macrophage and T lymphocyte apoptosis during experimental pulmonary tuberculosis: their relationship to mycobacterial virulence. European Journal of Immunology 2006, 36:345-353.

24. Park JM, Ng VH, Maeda S, Rest RF, Karin M: Anthrolysin O and Other Gram-positive Cytolysins Are Toll-like Receptor 4 Agonists. Journal of Experimental Medicine 2004, 200:1647-1655.

25. Kassam A, Der SD, Mogridge J: Differentiation of human mono-cytic cell lines confers susceptibility to Bacillus anthracis lethal toxin. Cell Microbiology 2005, 7:281-292.

26. Kohro T, Tanaka T, Murakami T, Wada Y, Aburatani H, Hamakubo T, Kodama T: A Comparison of Differences in the Gene Expression Profiles of Phorbol 12-myristate 13-acetate Differentiated THP-1 Cells and Human Monocyte-derived Macrophage. Journal of Atherosclerosis and Thrombosis 2004, 11:88-97.

27. Nau GJ, Richmond JFL, Schlesinger A, Jennings EG, Lander ES, Young RA: Human macrophage activation programs induced by bacterial pathogens. Proceedings of the National Acadamies of Sciences, USA 2002, 99:1503-1508.

28. Boldrick JC, Alizadeh AA, Diehn M, Dudoit S, Liu CL, Belcher CE, Botstein D, Staudt LM, Brown PO, Relman DA: Stereotyped and specific gene expression programs in human innate immune responses to bacteria. Proceedings of the National Acadamies of Sciences, USA 2002, 99:972-977.

29. SG Cherkasova TD, Shalygina NB, Proskurina VA, Burovtseva NP, Neliapin NM: The prostaglandin and cyclic nucleotide levels in the organs of mice with experimental anthrax intoxication. Zh Mikrobiol Epidemiol Immunobiol 1991, 3:70-73.

30. Alileche A, Squires RC, Muehlbauer SM, Lisanti MP, Brojatsch J: Mitochondrial Impairment Mediates Cytolysis in Anthrax Lethal Toxin-Treated Murine Macrophages. Cell Cycle 2006, 5:100-106.

31. Travert M, Ame-Thomas P, Pangault C, Morizot A, Micheau O, Semana G, Lamy T, Fest T, Tarte K, Guillaudeux T: CD40 Ligand Protects from TRAIL-Induced Apoptosis in Follicular Lymphomas through NF-{kappa}B Activation and Up-Regulation of c-FLIP and Bcl-xL. Journal of Immunology 2008, 181:1001-1011.

32. Allen C, Saigal K, Nottingham L, Arun P, Chen Z, Van Waes C: Bort-ezomib-Induced Apoptosis with Limited Clinical Response Is Accompanied by Inhibition of Canonical but not Alternative Nuclear Factor-{kappa}B Subunits in Head and Neck Cancer. Clinical Cancer Research 2008, 14:4175-4185.

33. Moayeri M, Haines D, Young HA, Leppla SH: Bacillus anthracis lethal toxin induces TNF-{alpha}-independent hypoxia-mediated toxicity in mice. Journal of Clinical Investigation 2003, 112:670-682.

34. Hanna P: Lethal toxin actions and their consequences. Journal of Applied Microbiology 1999, 87:285-287.

35. Hanna PC, Acosta D, Collier RJ: On the Role of Macrophages in

Anthrax. Proceedings of the National Acadamies of Sciences, USA 1993, 90:10198-10201.

36. Brieland JK, Remick DG, LeGendre ML, Engleberg NC, Fantone JC: In Vivo Regulation of Replicative Legionella pneumophila Lung Infection by Endogenous Interleukin-12. Infection and Immunity 1998, 66:65-69.

37. Poll T van der, Keogh CV, Buurman WA, Lowry SF: Passive immunization against tumor necrosis factor-alpha impairs host defense during pneumococcal pneumonia in mice. American Journal of Respiratory and Critical Care Medicine 1997, 155:603-608.

38. Lamonica JM, Wagner M, Eschenbrenner M, Williams LE, Miller TL, Patra G, DelVecchio VG: Comparative Secretome Analyses of Three Bacillus anthracis Strains with Variant Plasmid Contents. Infection and Immunity 2005, 73:3646-3658.

39. Ha S-D, Ng D, Pelech SL, Kim SO: Critical Role of the Phosphati-dylinositol 3-Kinase/Akt/Glycogen Synthase Kinase-3 Signaling Pathway in Recovery from Anthrax Lethal Toxin-induced Cell Cycle Arrest and MEK Cleavage in Macrophages. 2007, 282:36230-36239.

40. Tournier J-N, Quesnel-Hellmann A, Mathieu J, Montecucco C, Tang W-J, Mock M, Vidal DR, Goossens PL: Anthrax Edema Toxin Cooperates with Lethal Toxin to Impair Cytokine Secretion during Infection of Dendritic Cells. Journal of Immunology 2005, 174:4934-4941.

41. Drysdale M, Olson G, Koehler TM, Lipscomb MF, Lyons CR: Murine Innate Immune Response to Virulent Toxigenic and Nontox-igenic Bacillus anthracis Strains. Infection and Immunity 2007, 75:1757-1764.

42. Kim SO, Jing Q, Hoebe K, Beutler B, Duesbery NS, Han J: Sensitizing Anthrax Lethal Toxin-resistant Macrophages to Lethal Toxin-induced Killing by Tumor Necrosis Factor-alpha. Journal of Biological Chemistry 2003, 278:7413-7421.

43. McAllister RD, Singh Y, du Bois WD, Potter M, Boehm T, Meeker ND, Fillmore PD, Anderson LM, Poynter ME, Teuscher C: Susceptibility to Anthrax Lethal Toxin Is Controlled by Three Linked Quantitative Trait Loci. American Journal of Pathology 2003, 163:1735-1741.

44. Karin M, Lin A: NF-[kappa]B at the crossroads of life and death. Nat Immunol 2002, 3:221-227.

45. Chen G, Goeddel DV: TNF-RI Signaling: A Beautiful Pathway. Science 2002, 296:1634-1635.

46. Bergman NH, Passalacqua KD, Gaspard R, Shetron-Rama LM, Quack-enbush J, Hanna PC: Murine Macrophage Transcriptional Responses to Bacillus anthracis Infection and Intoxication. Infection and Immunity 2005, 73:1069-1080.

47. Tucker AE, Salles II, Voth DE, Ortiz-Leduc W, Wang H, Dozmorov I, Centola M, Ballard JD: Decreased glycogen synthase kinase 3-beta levels and related physiological changes in Bacillus anthracis lethal toxin-treated macrophages. Cellular Microbiology 2003, 5:523-532.

48. Hengartner MO: The biochemistry of apoptosis. Nature 2000, 407:770-776.

49. Park JK, Chung YM, Kang S, Kim J-U, Kim Y-T, Kim HJ, Kim YH, Kim JS, Yoo YD: c-Myc Exerts a Protective Function through Ornithine Decarboxylase against Cellular Insults. Molecular Pharmo-cology 2002, 62:1400-1408.

50. Hestvik ALK, Hmama Z, Av-Gay Y: Kinome Analysis of Host Response to Mycobacterial Infection: a Novel Technique in Proteomics. Infection and Immunity 2003, 71:5514-5522.

51. Saez-Lopez C, Ngambe-Tourere E, Rosenzwajg M, Petit J-C, Nicolas J-C, Gozlan J: Immediate-early antigen expression and modulation of apoptosis after in vitro infection of polymorphonu-clear leukocytes by human cytomegalovirus. Microbes and Infection 2005, 7:1139-1149.

52. Saito K, Meyer K, Warner R, Basu A, Ray RB, Ray R: Hepatitis C Virus Core Protein Inhibits Tumor Necrosis Factor AlphaMediated Apoptosis by a Protective Effect Involving Cellular FLICE Inhibitory Protein. Journal of Virology 2006, 80:4372-4379.

53. Allen HL, Deepe GS Jr: Apoptosis modulates protective immunity to the pathogenic fungus Histoplasma capsulatum. Journal of Clinical Investigation 2005, I 15:2875-2885.

54. Park JM, Greten FR, Wong A, Westrick RJ, Simon J, Arthur C, Otsu K, Hoffmann A, Montminy M, Karin M: Signaling Pathways and Genes that Inhibit Pathogen-Induced Macrophage Apoptosis - CREB and NF-kB as Key Regulators. Immunity 2005, 23:319-329.

55. Zhang J, Jiang R, Takayama H, Tanaka Y: Survival of Virulent Myco-bacterium tuberculosis Involves Preventing Apoptosis Induced by Bcl-2 Upregulation and Release Resulting from Necrosis in J774 Macrophages. MICROBIOLOGY and IMMUNOLOGY 2005, 49:845-852.

56. Roger P-M, Bermudez LE: Infection of Mice with Mycobacterium avium Primes CD8+ Lymphocytes for Apoptosis upon Exposure to Macrophages. Clinical Immunology 2001, 99:378-386.

57. Balcewicz-Sablinska MK, Keane J, Kornfeld H, Remold HG: Pathogenic Mycobacterium tuberculosis Evades Apoptosis of Host Macrophages by Release of TNF-R2, Resulting in Inactivation of TNF-{alpha}. Journal of Immunology 1998, 161:2636-2641.

58. RM Dhiman R, Majumdar S: Differential expression of NF-kB in mycobacteria infected THP-1 affects apoptosis. Biochimica et Biophysica Acta (BBA) - General Subjects 2007, 1770:649-658.

59. Abu-Zant A, Santic M, Molmeret M, Jones S, Helbig J, Abu Kwaik Y: Incomplete Activation of Macrophage Apoptosis during Intracellular Replication of Legionella pneumophila. Infection and Immunity 2005, 73:5339-5349.

60. Park JM, Greten FR, Li Z-W, Karin M: Macrophage Apoptosis by Anthrax Lethal Factor Through p38 MAP Kinase Inhibition. Science 2002, 297:2048-2051.

61. Ezzell JW, Welkos SL: The capsule of Bacillus anthracis, a review. Journal of Applied Microbiology 1999, 87:250-250.

62. Agrawal A, Lingappa J, Leppla SH, Agrawal S, Jabbar A, Quinn C, Pulendran B: Impairment of dendritic cells and adaptive immunity by anthrax lethal toxin. Nature 2003, 424:329-334.

63. Obata T, Brown GE, Yaffe MB: MAP kinase pathways activated by stress: The p38 MAPK pathway. Critical Care Medicine 2000, 28:N67-N77.

64. Louis Pasteur CC, Roux E: Nouvelle communication sur la rage. CR Acad Sci, Paris 1884, 98:457-463.

65. Helgason E, Tourasse NJ, Meisal R, Caugant DA, Kolsto A-B: Multi-locus Sequence Typing Scheme for Bacteria of the Bacillus cereus Group. Applied and Environmental Microbiology 2004, 70:191-201.

66. MM Kerr MK, Churchill GA: Analysis of Variance for Gene

Expression Microarray Data. Journal of Computational Biology 2000, 7:819-837.

67. Smyth GK: Limma: linear models for microarray data. Bioin-formatics and Computational Biology Solutions using R and Bioconductor. New York: Springer; 2005.

68. Benjamini Y, Hochberg Y: Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society 1995, 57:289-300.

69. Livak KJ, Schmittgen TD: Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2-[Delta][Delta]CT Method. Methods 2001, 25:402-408.

70. Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proceedings of the National Acadamies of Sciences, USA 1998, 95:14863-14868.

Publish with BioMecl Central and every scientist can read your work free of charge

"BioMed Central will be the most significant development for disseminating the results of biomedical research in our lifetime." Sir Paul Nurse, Cancer Research UK

Your research papers will be:

• available free of charge to the entire biomedical community

• peer reviewed and published immediately upon acceptance

• cited in PubMed and archived on PubMed Central

• yours — you keep the copyright

Submit your manuscript here: i } BioMedcentral ^